仙台での連星・変光星・低温度星研究会で、Mhhさんに、適当なサーベイ画像を渡して、Mhhの動かされているシステムで、テストをしていただきました。天体検出、カタログマッチングによるゼロ点の決定と測光、既知の変光星の抽出ができます。
U Monの画像だったのですが、Z CMaとかも検出されていて、同じのがうちでも動けば、サーベイ用画像から副産物で、沢山の変光星の測光データが得られることになります。
で、その野望に一歩近づくために、まず、SExtratorのinstallから。
#え、以前入れたんではないかって、あれは、テスト用のPCだったし、installできたところで止まっているので。
メモ
AstrOmaticから、SExtractorのrpmファイルを落とす。
#とりあえず、現時点で、version 2.8.6
システムエンジニア奮闘記 〜ダメダメSEの備忘録〜のUbuntuでRPMをインストールする件を参考に、rpmをdebにする。
まず、alienは、synapticパッケージマネージャーでいれた。
変換は、
できたsextractor_2.8.6-2_i386.debからのinstallは、Ubuntuソフトウェアセンターでした。
早速、動かしてみる
$sex V1723Aql-S001-R001-C001-V.fts
> WARNING: default.sex not found, using internal defaults
- SExtractor 2.8.6 started on 2010-11-25 at 15:22:21 with 4 threads
> Setting catalog parameters
> *ERROR*: can't read default.param
default.paramが無いって怒られました。
マニュアルによると、default.paramは、localに、default.sexってファイルを置けばよいらしいです。
パラメータの確認は、
$ sex -d
# Default configuration file for SExtractor 2.8.6
# EB 2009-04-09
#
#-------------------------------- Catalog ------------------------------------
CATALOG_NAME test.cat # name of the output catalog
CATALOG_TYPE ASCII_HEAD # NONE,ASCII,ASCII_HEAD, ASCII_SKYCAT,
# ASCII_VOTABLE, FITS_1.0 or FITS_LDAC
PARAMETERS_NAME default.param # name of the file containing catalog contents
#------------------------------- Extraction ----------------------------------
DETECT_TYPE CCD # CCD (linear) or PHOTO (with gamma correction)
DETECT_MINAREA 5 # minimum number of pixels above threshold
DETECT_THRESH 1.5 #or , in mag.arcsec-2
ANALYSIS_THRESH 1.5 #or , in mag.arcsec-2
FILTER Y # apply filter for detection (Y or N)?
FILTER_NAME default.conv # name of the file containing the filter
DEBLEND_NTHRESH 32 # Number of deblending sub-thresholds
DEBLEND_MINCONT 0.005 # Minimum contrast parameter for deblending
CLEAN Y # Clean spurious detections? (Y or N)?
CLEAN_PARAM 1.0 # Cleaning efficiency
MASK_TYPE CORRECT # type of detection MASKing: can be one of
# NONE, BLANK or CORRECT
#------------------------------ Photometry -----------------------------------
PHOT_APERTURES 5 # MAG_APER aperture diameter(s) in pixels
PHOT_AUTOPARAMS 2.5, 3.5 # MAG_AUTO parameters:,
PHOT_PETROPARAMS 2.0, 3.5 # MAG_PETRO parameters:,
#
SATUR_LEVEL 50000.0 # level (in ADUs) at which arises saturation
SATUR_KEY SATURATE # keyword for saturation level (in ADUs)
MAG_ZEROPOINT 0.0 # magnitude zero-point
MAG_GAMMA 4.0 # gamma of emulsion (for photographic scans)
GAIN 0.0 # detector gain in e-/ADU
GAIN_KEY GAIN # keyword for detector gain in e-/ADU
PIXEL_SCALE 1.0 # size of pixel in arcsec (0=use FITS WCS info)
#------------------------- Star/Galaxy Separation ----------------------------
SEEING_FWHM 1.2 # stellar FWHM in arcsec
STARNNW_NAME default.nnw # Neural-Network_Weight table filename
#------------------------------ Background -----------------------------------
BACK_SIZE 64 # Background mesh:or ,
BACK_FILTERSIZE 3 # Background filter:or ,
BACKPHOTO_TYPE GLOBAL # can be GLOBAL or LOCAL
#------------------------------ Check Image ----------------------------------
CHECKIMAGE_TYPE NONE # can be NONE, BACKGROUND, BACKGROUND_RMS,
# MINIBACKGROUND, MINIBACK_RMS, -BACKGROUND,
# FILTERED, OBJECTS, -OBJECTS, SEGMENTATION,
# or APERTURES
CHECKIMAGE_NAME check.fits # Filename for the check-image
#--------------------- Memory (change with caution!) -------------------------
MEMORY_OBJSTACK 3000 # number of objects in stack
MEMORY_PIXSTACK 300000 # number of pixels in stack
MEMORY_BUFSIZE 1024 # number of lines in buffer
#----------------------------- Miscellaneous ---------------------------------
VERBOSE_TYPE NORMAL # can be QUIET, NORMAL or FULL
WRITE_XML N # Write XML file (Y/N)?
XML_NAME sex.xml # Filename for XML output
NTHREADS 0 # Number of simultaneous threads for
それぞれのパラメータは、自分の画像に合わせて決めて、パラメータファイルに書き込んでおく。
現在、ここまで。
なお、Mhhさんに教えていただいたところによると、パラメータにX_WORLDとY_WORLDを書いておけば、x、yでなく赤緯赤経を書き出してくれるそうです。
で、とりあえず今後の予定
画像の一時処理、WCSの書き込みは、MaxImDL+Pinpointsでやる。
Sextractorで、星の検出
検出された星のlistから変光星を抜き出すところは、Mhhさんに教えてもらう。
#場合によっては、カタログ自体は、Mhhからもらう。
MySQLで、変光星名で管理するデータベースをつくる。
抜き出した変光星の測光結果をstd形式にする所は、Rでできそう?