;;;;;;;;;;;;;;;;;;; syori1.pro ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; dir ='E:\25CM\PMIDATA\2009\091016\' ; <== データのあるディレクトリ outdir ='E:\25CM\data\091016\' ; <== 結果を出力するディレクトリ fname ='pr00_*.pmi' ; <== シーケンスファイル名の*指定 dark ='dark00.pmi' ; <== ダークの入ったファイル名 flat ='flat00.pmi' ; <== フラットの入ったファイル名 fdark ='dark00.pmi' ; <== フラット用ダークの入ったファイル名 ;------------------------------------------------------------------------- cd,dir ; <== ディレクトリをファイルのある場所に nkrget,dark,hd,drk ; <== ダークを読み込みます nkrget,flat,hf,flt ; <== フラットを読み込みます nkrget,fdark,hfd,fdrk ; <== フラット用ダークを読み込みます sbtdrk,hf,flt,hfd,fdrk ; <== フラットの準備をします files=findfile(fname) ; <== シーケンスファイル名が読み込まれます n=n_elements(files) ; <== 処理するファイルの数 for i=0, n-1 do begin ; <== endfor と挟まれた部分をi=0からn-1まで繰り返す nkrget,files(i),h,img ; <== データ読み込み sbtdrk,h,img,hd,drk ; <== ダーク引き算 divflt,h,img,hf,flt ; <== フラット補正 tvscl,img ; <== 結果を表示します outfile=strmid(files(i),0, $ ; <== 出力ファイル名(outfile) を strlen(files(i))-3)+'dat' ;  'a_*.pmi'から 'a_*.dat'とします nkrsave,outdir+outfile,h,img ; <== outfile に保存します endfor print,'** syori1 finish! **' ; <== 終わったといいます end ; <== プログラムのおわりです ;;;;;;;;;;;;;;;;;;; end of program ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;