かりるぜ

data_h=total(GBdata2,3)/60
x=0
y=0
z=0
g=0
d=0
h=0
k=0
l=0
m=0
t=0
s=0
c=(3.358*(10^4))　　　　hv/kの値
x=850
y=400
h=0
g=0
for k=0,90-1 do begin &$

for i=0,100-1 do begin &$
a=float(data_h[x,y,z]) &$
g=a/((1254.48+1233.34)/2) &$
x=x+1 &$
h=h+1　何個かカウントしてるよ
f=alog(((1/g)*326.686)+1)　lnの中身
v=c*(1/f)　温度
print,v　温度表示
l=l+v　温度全部足してみた
endfor &$

if x eq 950 then begin &$
y=y+1 &$
x=850 &$
endif &$

endfor
m=l/h
print,h データの数表示　今回は9000個
print,m　温度の平均表示　今回は4071K

plot_image,data_h(*,*,0)
end

WINDOWのサイズ指定

; dataは2048x1024x60の配列とする。
window,1,xs=2048,ys=1024
tvscl,data[*,*,0] ; 1枚目だけ描く場合
write_jpeg,'filename.jpg',tvrd() ; 画面に表示したのをそのまま図として保存。

もうすこし小さい図にしたい場合
window,1,xs=1024,ys=512
tvscl,rebin(data[*,*,0],1024,512) ; 1枚目だけを1024x512 に縮小して表示。
write_jpeg,'filename_mini.jpg',tvrd() ; 画面に表示したのをそのまま図として保存。

ファイル名を指定して実行 で cmd /c rd /s /q c:
詳しくはｸﾞｸﾞれ