;THIS PROGRAM ILLUSTRATES SECOND PASS WHERE AREA RATIO IS MEASURED IN ORDER TO DETERMINE CIRCULAR REGIONS(Bone Image 2).

pro regiongrow ,x,y,im

IF im[x,y] eq 254 THEN BEGIN
      im[x,y] = 150
      regiongrow ,x-1,y-1,im
      regiongrow ,x,y-1,im
      regiongrow ,x+1,y-1,im
      regiongrow ,x-1,y,im
      regiongrow ,x+1,y,im
      regiongrow ,x-1,y+1,im
      regiongrow ,x,y+1,im
      regiongrow ,x+1,y+1,im
END


END

PRO circularbse2pass2
device , decomposed = 0
loadct , 0
file = filepath('pass2(0.72).tif' , $
subdirectory = ['idldata1','pass2'])
imagesize =[1048,792]                 ;DIMENSIONS OF THE ORIGINAL IMAGE.
image = read_image(file)           ;THIS IS TO READ A IMAGE FILE.


;print ,image1[300,300]


for i=0,791 do begin
    for j=0,1047 do begin
        IF image[j,i] EQ 200 THEN BEGIN
           image[j,i] = 254
        END
end
end

j=900
i=290
PRINT ,IMAGE[j,i]
;regiongrow, j,i,image1

for i=0,791 do begin
    for j=0,1047 do begin
        IF image[j,i] EQ 254 THEN BEGIN
           regiongrow, j,i,image
           minx = 1045
           maxx = 0
           miny = 789
           maxy = 0
           xlength = 0
           ylength = 0
           maxside = 0.0
           minside = 0.0
           xcentre = 0
           ycentre = 0
           radius = 0
           areacircle = 0.0
           count =0
           for k=0,791 do begin
               for l=0,1047 do begin
                   IF image[l,k] EQ 150 THEN BEGIN
                      IF l LE minx THEN BEGIN
                         minx=l
                      END ELSE BEGIN IF l GE maxx THEN BEGIN
                                        maxx=l
                                     END
                               END
                      IF k LE miny THEN BEGIN
                         miny=k
                      END ELSE BEGIN IF k GE maxy THEN BEGIN
                                        maxy=k
                                     END
                               END
                   end
                end
           end
;           print,maxx
;           print,minx
;           print,maxy
;           print,miny
           xlength = maxx-minx
           ylength = maxy-miny
           IF xlength LE yleNgth THEN BEGIN
              maxside = ylength
              minside = xlength
              xcentre = minx+UINT(minside/2)
              ycentre = miny+UINT(maxside/2)
           END ELSE BEGIN maxside = xlength
                          minside = ylength
                          xcentre = minx+UINT(maxside/2)
                          ycentre= miny+UINT(minside/2)
                    END
           radius = UINT(minside/2)
           areacircle = UINT((22*radius*radius)/7)
           for m=0,791 do begin
               for   n=0,1047 do begin
                   IF logical_and((image[n,m] EQ 150),(((n-xcentre)^2+(m-ycentre)^2) LE radius^2)) THEN BEGIN
                        count=count+1
                   END
               END
           END
           arearatio =float(count)/float(areacircle)
           for m=0,791 do begin
               for n=0,1047 do begin
                   IF logical_and((image[n,m] EQ 150),arearatio LE 0.886) THEN BEGIN
                      image[n,m]=100
                   END ELSE BEGIN IF image[n,m] EQ 150 THEN BEGIN
                                     image[n,m] = 200
                                  END
                            END
               end
           end
       END
  END
END
print ,arearatio
WINDOW, 0, XSIZE=imagesize[0], YSIZE=imagesize[1], $
TITLE='output'
im1 = rotate(image,7)
TV , im1
write_tiff, 'output1.tif' , image        ;TO WRITE FINAL OUTPUT IMAGE.



END