PRO findedge,name,c0,r0,radius,nosobel=nosobel

namelen=STRLEN(STRTRIM(name,2))
cam=FIX(STRMID(name,namelen-9,1))
exp=FIX(STRMID(name,namelen-7,3))
rdkimg,name,h,im
rdmeandk,cam,exp,dk
cim=im-dk
IF KEYWORD_SET(nosobel) THEN sim=cim ELSE sim=SOBEL(cim)
sh=HISTOGRAM(sim,MIN=0,BIN=1)
csh=sh & FOR i=1,N_ELEMENTS(csh)-1 DO csh(i)=csh(i)+csh(i-1)
ceiling=MAX(WHERE(csh LT N_ELEMENTS(sim)/2))*10
dd=dist(256)
IF NOT (KEYWORD_SET(c0) OR KEYWORD_SET(r0) OR $
        KEYWORD_SET(radius)) THEN BEGIN
    c0=131
    r0=126
    radius=143
ENDIF
resp=''
REPEAT BEGIN
    tvsi,sim<ceiling,0,0,255,255
    CONTOUR,SHIFT(dd,c0,r0),/FOLLOW,LEVELS=[radius],C_LABELS=[1] $
    ,/DEVICE,POSITION=[80,40,529,489],COLOR=!D.N_COLORS-1 $
    ,/OVERPLOT,XSTYLE=1,YSTYLE=1 ;$
    ;,CLIP=[80,40,529,489],/NOERASE
    PRINT,c0,r0,radius $
         ,FORMAT='("Center=(",I3,",",I3,"), radius=",I3)'
    PRINT,'OK? (Y/[N])'
    READ,resp
    again=STRUPCASE(resp) NE 'Y'
    IF again THEN BEGIN
        PRINT,'Enter new values for center column and row, and radius'
        READ,c0,r0,radius
    ENDIF
ENDREP UNTIL NOT again
PRINT,c0,r0,radius $
     ,FORMAT='("Final values: Center=(",I3,",",I3,"), radius=",I3)'
RETURN
END