plot,xx,yy,psym=-2


ww=xx*0+100.
ww[61:87]=1.e-3

xxx=xx[61:87]
yyy=yy[61:87]
www=ww[61:87]

plot,xxx,yyy,psym=-2

par=[0,0,0, yy[(where(yy eq max(yy)))[0]], 88.65,0.1, yy[(where(yy eq max(yy)))[0]], 88.75,0.1]*1.
res=mgaussfit(xxx,yyy,www,par)
print,par

oplot,xxx,par[0]+par[1]*(xxx-par[4])/par[5]+par[2]*((xxx-par[4])/par[5])^2 $
+par[3]*exp(-0.5*((xxx-par[4])/par[5])^2) $
+par[6]*exp(-0.5*((xxx-par[7])/par[8])^2),color=254

.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

たまにスパムが入る。

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

テスト。