;THIS PROGRAM MERGES GREEN REGIONS SURVIVED FROM PASS 1 & 2 ON THE "cpl" IMAGES OF THE BONE(Bone Image 2).

pro mergecpl2
device , decomposed = 0
loadct , 0
file = filepath('finaloutput.tif' , $
subdirectory = ['idldata1','bse2'])
imagesize1 =[1048,792]
image1 = read_image(file)

file = filepath('cpl2.tif' , $
subdirectory = ['idldata1','bse2'])
imagesize2 =[1024,768]
image2 = read_image(file)

for i=0,767 do begin
   for j=0,1023 do begin
       IF image2[j,i] EQ 255 THEN BEGIN
          image2[j,i]=254
       END
    END
END

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

for i=0,767 do begin
   for j=0,1023 do begin
       IF image1[j+12,i+12] EQ 200 THEN BEGIN
          image2[j,i]=greencolor
       END
    END
END

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

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

end