.r
for ii=0,29999 do begin
if ii eq 0 then num=nn[ii]
if num eq nn[ii] then plots,[0,0]+xx[ii],[0,0]+yy[ii],/dev,color=count*1.0/max(nn)*180+74 else begin
print,num
num=nn[ii]
count++
gk
endelse
endfor
end

DATAFILE ='TRACE_19980519.fits'
LOOPFILE ='TRACE_19980519.dat'
IMAGE1 =READFITS(datafile,header)
NSM1 =3 ;lowpass filter
RMIN =30 ;minimum curvature radius of loop (pixels)
LMIN =25 ;minimum loop length (in pixels)
NSTRUC =1000 ;maximum limit of traced structures used in array dimension
NLOOPMAX =1000 ;maximum number of detected loops
NGAP =0 ;number of pixels in loop below flux threshold (0,...3)
THRESH1 =0.0 ;ratio of image base flux to median flux
THRESH2 =3 ;threshold in detected structure (number of significance ;levels with respect to the median flux in loop profile
TEST =1001 ;option for display of traced structures if TEST < NSTRUC
PARA2 =[NSM1,RMIN,LMIN,NSTRUC,NLOOPMAX,NGAP,THRESH1,THRESH2]
autoloop_a4,IMAGE1,IMAGE2,LOOPFILE,PARA2,OUTPUT,TEST

たまにスパムが入る。

掲示板を置いておく意味はもはやほとんどないかも。

テスト。