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

PRO finalpassbse2
device , decomposed = 0
loadct , 0
file = filepath('finaloutput.tif' , $
subdirectory = ['idldata1','bse2'])
imagesize =[1048,792]
image = read_image(file)

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

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

for i=0,791 do begin
    for j=0,1047 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,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-3
tvlct ,0,0,255, bluecolor

for i=0,791 do begin
    for j=0,1047 do begin
        IF image1[j,i] EQ 99 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 greencolor)||(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