前回の実験から計算するとReadout Noiseが45e位になり、大きすぎたので再チェック。
2.) 直接ライトを当てると明るすぎるようなので、一面のみ開けたダンボール箱内面に反射させた光を当てる(下の写真) ![]() 試し撮りの結果、
3.) 以前にダンボール箱の中で実験したときのように、カメラのシャッター前方を封筒で遮蔽した ![]() 試し撮りの結果、
以下、このセットアップで実験を進める ●撮像、解析(2009/09/15 :小西、利川) 以下のように進めた #!/bin/sh #100msec-4000msecで各2枚ずつバイアス、ダーク、オブジェクト画像の撮像 awk -v n=2 'BEGIN{for(i=100;i<=4000;i+=300){for (j=1;j<=n;j++){print "./fliexp 1 1 e 0"};for (j=1;j<=n;j++){print "./fliexp 1 1 d",i/1000};for (j=1;j<=n;j++){print "./fliexp 1 1 e", i/1000};printf("echo \"echo sleep 1\"\n");print "sleep 1"}}' | sh #画像リストファイルの作成 cd ~/data/20090612 rm conversion.090915* awk 'BEGIN{ns=2013;for(i=0;i<84;i++)printf("ANOA%08d.fits\n",ns+i)}' > conversion.090915 #ヘッダの書き出し cat <<EOF|cl hselect @conversion.090915 "\$I,DATA-TYP,EXPTIME" yes > conversion.090915.head logout EOF #各リストファイルの作成 awk '$3==0{print $1}' conversion.090915.head > conversion.090915.bias awk '$2=="DARK"{print $1}' conversion.090915.head > conversion.090915.dark awk '$2=="TEST" && $3>0{print $1}' conversion.090915.head > conversion.090915.obj awk '{print "sub"$1"[384:640,384:640]"}' conversion.090915.obj > conversion.090915.sub.sectバイアスの撮像からReadout Noiseをチェック #バイアスを撮像したものの標準偏差からReadout Noise(ADU)がわかる cat <<EOF|cl imstat @conversion.090915.bias nclip=10 fields="image,mean,midpt,stddev" # IMAGE MEAN MIDPT STDDEV # ANOA00002013.fits 1982. 1981. 8.003 # ANOA00002014.fits 1982. 1981. 8.016 # ANOA00002019.fits 1982. 1981. 8.009 # ANOA00002020.fits 1982. 1982. 8.031 # ANOA00002025.fits 1982. 1981. 8.005 # ANOA00002026.fits 1982. 1982. 8.036 # ANOA00002031.fits 1982. 1982. 8.005 # ANOA00002032.fits 1982. 1982. 8.031 # ANOA00002037.fits 1982. 1982. 8.032 # ANOA00002038.fits 1982. 1982. 8.007 # ANOA00002043.fits 1982. 1981. 8.004 # ANOA00002044.fits 1982. 1982. 8.04 # ANOA00002049.fits 1982. 1982. 8.026 # ANOA00002050.fits 1982. 1982. 8.028 # ANOA00002055.fits 1982. 1982. 8.035 # ANOA00002056.fits 1982. 1982. 8.032 # ANOA00002061.fits 1982. 1982. 8.032 # ANOA00002062.fits 1982. 1982. 8.005 # ANOA00002067.fits 1982. 1981. 8.004 # ANOA00002068.fits 1982. 1982. 8.035 # ANOA00002073.fits 1982. 1982. 8.021 # ANOA00002074.fits 1982. 1982. 8.033 # ANOA00002079.fits 1982. 1981. 8.013 # ANOA00002080.fits 1982. 1982. 8.027 # ANOA00002085.fits 1982. 1982. 8.019 # ANOA00002086.fits 1982. 1982. 8.031 # ANOA00002091.fits 1982. 1981. 8.019 # ANOA00002092.fits 1982. 1981. 8.01Readout Noiseはだいたい8.01-8.03[ADU]であるとわかる カウントの、露出時間に対するリニアリティをチェックする #積分時間に対する、カウントのリニアリティをチェックする #オブジェクト画像からダークを引く !awk '{print "subdark"$1}' conversion.090915.dark > conversion.090915.sub.dark imarith @conversion.090915.obj - @conversion.090915.dark @conversion.090915.sub.dark #ダークを引いた画像のexptimeを書き出す hselect @conversion.090915.obj "\$EXPTIME" yes > conversion.090915.head.sub.dark.exptime #ダークを引いた画像の、光量の一様な領域の、カウントの平均値をとる !awk '{print $1"[384:640,384:640]"}' conversion.090915.sub.dark > conversion.090915.sub.dark.sect imstat @conversion.090915.sub.dark.sect nclip=10 fields="mean" format- > conversion.090915.sub.dark.mean #exptime-meanのテーブルファイルの作成 !paste conversion.090915.head.sub.dark.exptime conversion.090915.sub.dark.mean > conversion.090915.sub.dark.exp-mean logout EOF #gnuplotでフィット cat <<EOF|gnuplot set term png set output "conversion.090915.liner.png" f(x)=a*x+b fit f(x) " |