;   THIS PRGRAM IS TO INSERT RED AND BLUE COLOR IN THE REGIONS EXTRACTED FROM PASS 1(Bone Image 2).

PRO pass1bse2
device , decomposed = 0
loadct , 0
file = filepath('pass2(0.72).tif' , $
subdirectory = ['idldata1','pass2'])
imagesize =[1048,792]
image1 = read_image(file)

file = filepath('bse2.tif' , $
subdirectory = ['idldata1'])
imagesize1 =[1023,767]
original = read_image(file)

redcolor = !d.table_size-1
tvlct ,255,0,0,redcolor
image1 = bytscl(image1,top=redcolor-1)

for i=0,791 do begin
    for j=0,1047 do begin
        IF image1[j,i] EQ 49 THEN BEGIN
           image1[j,i]=redcolor
        END
    END
END

bluecolor = !d.table_size-2
tvlct ,0,0,255, bluecolor

for i=0,791 do begin
    for j=0,1047 do begin
        IF image1[j,i] EQ 199 THEN BEGIN
           image1[j,i]= bluecolor
        END
    END
END



for i=12,779 do begin
    for j=12,1035 do begin
        IF ((image1[j,i] EQ redcolor)||(image1[j,i] EQ bluecolor)) THEN BEGIN
           original[j-12,i-12]= image1[j,i]
        END
    END
END

WINDOW, 0, XSIZE=imagesize1[0], YSIZE=imagesize1[1], $
TITLE='ouput'
im1 = rotate(original,7)
TV , im1
;write_tiff, 'ouput.tif' , image1        ;TO WRITE FINAL OUTPUT IMAGE.

void = TVRead(/tiff, File='myimage')

end