;THIS PROGRAM ILLUSTRATES PHASE-I OF THE PROJECT.
;This program realized an image subtraction from threshold image bse1 to morphedclose image bse2.

pro subtractionbse1
file1 = filepath('thresholded(60) bse1.tif')
imagesize =[1046,790]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
img1 = read_image(file1)           ;THIS IS TO READ A IMAGE FILE.

file2 = filepath('morphedclose(r=5.5)bse1.tif')
imagesize =[1046,790]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
img2 = read_image(file2)           ;THIS IS TO READ A IMAGE FILE.

initialcolor=!d.table_size-1
tvlct,0,0,0,initialcolor
img3=replicate(byte(initialcolor),imagesize[0],imagesize[1])


for i=0,1045 do begin
    for j=0,789 do begin
       img3[i,j]=img1[i,j]-img2[i,j]
    end
end

write_tiff, 'subtractionbse1.tif' , img3        ;TO WRITE FINAL OUTPUT IMAGE.



END