;THIS PROGRAM ILLUSTRATES PHASE-I OF THE PROJECT.

pro dilationbse1

file1 = filepath('subtractionbse1.tif')               ;open the image file
imagesize =[1024,768]                                 ;dimensions of the original image
width=1024
height=768
img1 = read_image(file1)                              ;load the image file.
img1=Bytscl(img1,Top=!D.Table_Size-1)

; Threshold and make binary image:
B = img1

; Create the shape operator:
S = REPLICATE(1, 3, 3)

; "Opening" operator:
img1 = DILATE(B, S)

scaled=Bytscl(img1,Top=!D.Table_Size-1)

write_tiff,"dilationbse1(r=3).tif",scaled



END