meineko’s blog

元つくばの某独立行政法人勤務の植物屋です。最近は、ほぼ、突発天体の話題です。

Sextractor

とりあえず、以下のようなスクリプトを作って、Icの3日分だけ測りました。

#!/bin/bash
sex -DETECT_THRESH 2.0 -PHOT_APERTURES 5.0 -FILTER N -CATALOG_NAME temp.out -PARAMETERS_NAME ./sextractor_param/wcstool.param -BACKPHOTO_THICK 15.0 -BACKPHOTO_TYPE LOCAL $1.fit
python scat2id.py temp.out > temp.id
python identifyIc.py temp.id > $1.mag
grep "GCVS" $1.mag > $1.gcvs

それだけでも、のべ1094星になりました。
#Vは、数がありすぎて、勘弁です(w


観測日時については、MIRAにlogを作らせて、そこから転記することに。
メニューからcatalogを選んで、画像ファイルをまとめて選択すると、簡単に一覧表ができます。


とりあえず、結果をぱらぱら見ていたら、U Pupが、Icで明るかったです。
#U Pupって、std形式ではPUPUって書くのですが、「ぷぷ」ってなんかいいな(w


(追記)
まとめて解析しやすいように改変

#!/bin/bash
sex -DETECT_THRESH 2.0 -PHOT_APERTURES 5.0 -FILTER N -CATALOG_NAME temp.out -PARAMETERS_NAME ./sextractor_param/wcstool.param -BACKPHOTO_THICK 15.0 -BACKPHOTO_TYPE LOCAL $1.fit
python scat2id.py temp.out > temp.id
python identify.py temp.id > $1.mag
echo $1 >> results.gcvs
grep "GCVS" $1.mag >> results.gcvs


(追記2)
2011年に入ってからの分を全て測り終えました。
Vが、1267星4342、Icが、678星1094でした。
でも、同定ミスが沢山ありそうです。
SCOTちゃう、それ、M80や。