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

PRO finalpassbse1

device , decomposed = 0
loadct , 0
file = filepath('finaloutput.tif' , $
subdirectory = ['idldata1','bse1'])
imagesize =[1068,812]
image = read_image(file)

file = filepath('bse1.tif' , $
subdirectory = ['idldata1'])
imagesize1 =[1024,768]
original = read_image(file)

greencolor = !d.table_size-1
tvlct ,0,255,0,greencolor

image1 = bytscl(image,top=greencolor-1)



for i=0,811 do begin
    for j=0,1067 do begin
        IF image1[j,i] EQ 199 THEN BEGIN
           image1[j,i]=greencolor
        END
    END
END

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

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-3
tvlct ,0,0,255, bluecolor

for i=0,811 do begin
    for j=0,1067 do begin
        IF image1[j,i] EQ 99 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 greencolor)||(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, 'finaloutput.tif' , image1        ;TO WRITE FINAL OUTPUT IMAGE.

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

end