基本

  • ds9を立ち上げておく。
  • irafのホームディレクトリで
    cl
  • M51
    disp dev$pix
    plot
    surface dev$pix

daofind

test01.jpeg

点源をガウシアンで近似し、検出する。x=194,y=288

 
  • 表示
    disp test01.fits
    daofind00.jpg
 
  • imexamine「pixel valueを見る。」
    imexamine
    で黒丸が現れる。これをds9の画像上に持っていき、'r'をクリックすると図のようなRadius-Pixel Valueの図が現れる。
    daofind01.jpg
     
    ソースの無い箇所で'r'をクリックし、標準偏差を調べる。今の場合'20'であることが分かる。
    daofind02.jpg
     
    'q'で終わる。
 
  • Task
    noao
    digiphot
    apphot
 
  • daofind
    ap> daofind
    Input image(s) (03.fits): test01.fits
    
    FWHM of features in scale units (2.5) (CR or value): 10
            New FWHM of features: 10. scale units  10. pixels
    Standard deviation of background in counts (INDEF) (CR or value): 20
            New standard deviation of background: 20. counts
    Detection threshold in sigma (4.) (CR or value): 
            New detection threshold: 4. sigma 80. counts
    Minimum good data value (INDEF) (CR or value): 
            New minimum good data value: INDEF counts
    Maximum good data value (INDEF) (CR or value): 
            New maximum good data value: INDEF counts
    確実に一点だけを検出するようにFWHMは大きめにとる。ここでは'10'。 標準偏差はさっきの'20'。
    #K IRAF       = NOAO/IRAFV2.12.2-EXPORT version    %-23s     
    #K USER       = chinone                 name       %-23s     
    #K HOST       = magi                    computer   %-23s     
    #K DATE       = 2007-01-05              yyy-mm-dd  %-23s     
    #K TIME       = 16:43:27                hh:mm:ss   %-23s     
    #K PACKAGE    = apphot                  name       %-23s     
    #K TASK       = daofind                 name       %-23s     
    #
    #K SCALE      = 1.                      units      %-23.7g   
    #K FWHMPSF    = 10.                     scaleunit  %-23.7g   
    #K EMISSION   = yes                     switch     %-23b     
    #K DATAMIN    = INDEF                   counts     %-23.7g   
    #K DATAMAX    = INDEF                   counts     %-23.7g   
    #K EXPOSURE   = ""                      keyword    %-23s     
    #K AIRMASS    = ""                      keyword    %-23s     
    #K FILTER     = ""                      keyword    %-23s     
    #K OBSTIME    = ""                      keyword    %-23s     
    #
    #K NOISE      = poisson                 model      %-23s     
    #K SIGMA      = 20.                     counts     %-23.7g   
    #K GAIN       = ""                      keyword    %-23s     
    #K EPADU      = 1.                      e-/adu     %-23.7g   
    #K CCDREAD    = ""                      keyword    %-23s     
    #K READNOISE  = 0.                      e-         %-23.7g   
    #
    #K IMAGE      = test01.fits             imagename  %-23s     
    #K FWHMPSF    = 10.                     scaleunit  %-23.7g   
    #K THRESHOLD  = 4.                      sigma      %-23.7g   
    #K NSIGMA     = 1.5                     sigma      %-23.7g   
    #K RATIO      = 1.                      number     %-23.7g   
    #K THETA      = 0.                      degrees    %-23.7g   
    #
    #K SHARPLO    = 0.2                     number     %-23.7g   
    #K SHARPHI    = 1.                      number     %-23.7g   
    #K ROUNDLO    = -1.                     number     %-23.7g   
    #K ROUNDHI    = 1.                      number     %-23.7g   
    #
    #N XCENTER   YCENTER   MAG      SHARPNESS   SROUND          GROUND      ID         \
    #U pixels    pixels        #        #           #           #           #          \
    #F %-13.3f   %-10.3f   %-9.3f   %-12.3f     %-12.3f     %-12.3f     %-6d       \
    #
       193.609   194.892   -0.826   0.743       -0.077      0.008       1     
    これよりx=193.609, y'=194.892であることが分かる。先の座標系に合わせると、 y=480-194.892=285.108
 
  • tvmark
    tvmark test01.fits.coo.1
    でピンク色の点が検出されたソースの中心に載ることを確認する。
    daofind03.jpg
 

jpeg -> fits

jpeg画像をfitsに変換する為のラッパー(参考

#!/usr/bin/perl

# USAGE: 
# This program convert jpeg files to fits files.
#   jpeg2fits [jpegfiles....]
#
# This program need "convert" in ImageMagick.
#   Shigetomo Shiki, 4th Dec 2001

foreach $jpegfile (@ARGV){
    $_=$jpegfile;
    s/(\.jpg|\.jpeg)/\.fits/;
    $fitsfile = $_;
    print stderr "convert $jpegfile $fitsfile\n";
    system "convert $jpegfile $fitsfile";
}

使い方

jpeg2fits *.jpeg

添付ファイル: filedaofind03.jpg 798件 [詳細] filedaofind02.jpg 873件 [詳細] filedaofind01.jpg 848件 [詳細] filedaofind00.jpg 834件 [詳細] filetest01.jpeg 817件 [詳細]