これは難しい

宇宙的な色らしい。
投稿者
メッセージ  (タグ使用可)   ひらがながないと投稿できません。リンクは書き込めません。

削除キー 

 投稿者: joten    投稿日: 2018年04月25日 (水) 20時33分34秒
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

 投稿者: joten    投稿日: 2018年03月26日 (月) 22時06分06秒
.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

 投稿者: joten    投稿日: 2018年03月23日 (金) 16時16分00秒
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

 投稿者: じょうてん    投稿日: 2016年01月21日 (木) 04時53分37秒
たまにスパムが入る。

 投稿者: じょうてん    投稿日: 2015年02月24日 (火) 13時31分21秒
掲示板を置いておく意味はもはやほとんどないかも。

 投稿者: じょうてん    投稿日: 2014年07月03日 (木) 14時44分45秒
テスト。

    過去ログ  削除
MASTER
エース