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

PRO pass1bse1
device , decomposed = 0
loadct , 0
file = filepath('pass1(0.95).tif' , $
subdirectory = ['idldata1','pass1'])
imagesize =[1068,812]
image1 = read_image(file)

file = filepath('bse1.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,811 do begin
    for j=0,1067 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,811 do begin
    for j=0,1067 do begin
        IF image1[j,i] EQ 199 THEN BEGIN
           image1[j,i]= bluecolor
        END
    END
END



for i=22,789 do begin
    for j=22,1045 do begin
        IF ((image1[j,i] EQ redcolor)||(image1[j,i] EQ bluecolor)) THEN BEGIN
           original[j-22,i-22]= 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