IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 12138


Ignore:
Timestamp:
Mar 1, 2007, 10:53:16 AM (19 years ago)
Author:
Paul Price
Message:

Adding shutter correction. Changing to use SIMTEST camera.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/simtest/sim.pro

    r12102 r12138  
    1 
    2 #list MODULES -add /home/kiawe/eugene/src/panstarrs/ipp/ippTasks
     1### Generate simulated test data for testing the Image Processing Pipeline
    32
    43module images.pro
    54
    6 $flatlevel = 10000
     5$flatlevel = 1000
    76$biaslevel = 500
    87$skyrate  = 200
    98$darkrate = 100
     9$shutter = 0.5
    1010
    1111macro cleanup
    … …  
    6161    keyword bias CCDSUM   -w "1 1"
    6262    keyword bias TELESCOP -w simscope
    63     keyword bias INSTRUME -w SIMPLE
     63    keyword bias INSTRUME -w SIMTEST
    6464    keyword bias OBSTYPE  -w BIAS
    6565    keyword bias OBJECT   -w BIAS
    6666    sprintf name "$simdir/raw/bias.%02d.fits" $i
     67    keyword bias GAIN     -wf 1.0
     68    keyword bias RDNOISE  -wf $mkiRD
    6769    wd bias $name
    6870  end
    … …  
    8789    keyword dark CCDSUM   -w "1 1"
    8890    keyword dark TELESCOP -w simscope
    89     keyword dark INSTRUME -w SIMPLE
     91    keyword dark INSTRUME -w SIMTEST
    9092    keyword dark OBSTYPE  -w DARK
    9193    keyword dark OBJECT   -w DARK
     94    keyword dark GAIN     -wf 1.0
     95    keyword dark RDNOISE  -wf $mkiRD
    9296    sprintf name "$simdir/raw/dark.%02d.fits" $i
    9397    wd dark $name
    9498  end 
     99end
     100
     101macro mkshutter
     102  echo "generating 10 shutter images"
     103  for i 0 20
     104    $exptime = 3 * ($i + 1)
     105    $opentime = $exptime - $shutter
     106    mkimages.init 1024 1024 5 1000 0 2.0
     107    set shutter = blank + sqrt($mkiRD*$mkiRD + $darkrate*$exptime + $flatlevel*$opentime)*mki_noisei + $flatlevel*$opentime + $darkrate*$exptime + $biaslevel
     108    keyword shutter AIRMASS  -wf 1.0
     109    keyword shutter FILTER   -w  NONE
     110    keyword shutter POSANGLE -wf 0.0
     111    keyword shutter RA       -wf 0.0
     112    keyword shutter DEC      -wf 0.0
     113    keyword shutter EXPTIME  -wf $exptime
     114    keyword shutter DARKTIME -wf $exptime
     115    keyword shutter DATE-OBS -w "2006/10/01"
     116    sprintf time "00:01:%02d" $i
     117    keyword shutter UTC-OBS  -w $time
     118    keyword shutter CCDSUM   -w "1 1"
     119    keyword shutter TELESCOP -w simscope
     120    keyword shutter INSTRUME -w SIMTEST
     121    keyword shutter OBSTYPE  -w SHUTTER
     122    keyword shutter OBJECT   -w SHUTTER
     123    keyword shutter GAIN     -wf 1.0
     124    keyword shutter RDNOISE  -wf $mkiRD
     125    sprintf name "$simdir/raw/shutter.%02d.fits" $i
     126    wd dark $name
     127  end
    95128end
    96129
    … …  
    98131  echo "generating 20 flat images in 2 filters"
    99132  $exptime = 10.0
     133  $opentime = $exptime - $shutter
    100134  for i 0 10
    101135    mkimages.init 1024 1024 5 1000 0 2.0
    102     set flat = blank + sqrt($mkiRD*$mkiRD + $darkrate*$exptime + $flatlevel)*mki_noisei + $darkrate*$exptime + $flatlevel + $biaslevel
     136    set flat = blank + sqrt($mkiRD*$mkiRD + $darkrate*$exptime + $flatlevel*$opentime)*mki_noisei + $darkrate*$exptime + $flatlevel*$opentime + $biaslevel
    103137    keyword flat AIRMASS  -wf {1 + 0.1*$i}
    104138    keyword flat FILTER   -w  V
    … …  
    113147    keyword flat CCDSUM   -w "1 1"
    114148    keyword flat TELESCOP -w simscope
    115     keyword flat INSTRUME -w SIMPLE
     149    keyword flat INSTRUME -w SIMTEST
    116150    keyword flat OBSTYPE  -w FLAT
    117151    keyword flat OBJECT   -w FLAT
     152    keyword flat GAIN     -wf 1.0
     153    keyword flat RDNOISE  -wf $mkiRD
    118154    sprintf name "$simdir/raw/flat.%02d.fits" $i
    119155    wd flat $name
    … …  
    122158  for i 10 20
    123159    mkimages.init 1024 1024 5 1000 0 2.0
    124     set flat = blank + sqrt($mkiRD*$mkiRD + $darkrate*$exptime + $flatlevel)*mki_noisei + $darkrate*$exptime + $flatlevel + $biaslevel
     160    set flat = blank + sqrt($mkiRD*$mkiRD + $darkrate*$exptime + $flatlevel*$opentime)*mki_noisei + $darkrate*$exptime + $flatlevel*$opentime + $biaslevel
    125161    keyword flat AIRMASS  -wf {1 + 0.1*$i}
    126162    keyword flat FILTER   -w  B
    … …  
    135171    keyword flat CCDSUM   -w "1 1"
    136172    keyword flat TELESCOP -w simscope
    137     keyword flat INSTRUME -w SIMPLE
     173    keyword flat INSTRUME -w SIMTEST
    138174    keyword flat OBSTYPE  -w FLAT
    139175    keyword flat OBJECT   -w FLAT
     176    keyword flat GAIN     -wf 1.0
     177    keyword flat RDNOISE  -wf $mkiRD
    140178    sprintf name "$simdir/raw/flat.%02d.fits" $i
    141179    wd flat $name
    … …  
    147185  for i 0 5
    148186    $exptime = 3 * ($i + 1)
    149     mkimages.init 1024 1024 5 1000 0 2.0
    150     set obj = blank + sqrt($mkiRD*$mkiRD + $darkrate*$exptime + $skyrate*$exptime)*mki_noisei + $darkrate*$exptime + $skyrate*$exptime + $biaslevel
     187    $opentime = $exptime - $shutter
     188    mkimages.init 1024 1024 5 1000 0 2.0
     189    set obj = blank + sqrt($mkiRD*$mkiRD + $darkrate*$exptime + $skyrate*$opentime)*mki_noisei + $darkrate*$exptime + $skyrate*$opentime + $biaslevel
    151190    keyword obj AIRMASS  -wf {1 + 0.1*$i}
    152191    keyword obj FILTER   -w  V
    … …  
    161200    keyword obj CCDSUM   -w "1 1"
    162201    keyword obj TELESCOP -w simscope
    163     keyword obj INSTRUME -w SIMPLE
     202    keyword obj INSTRUME -w SIMTEST
    164203    keyword obj OBSTYPE  -w OBJECT
    165204    keyword obj OBJECT   -w OBJECT
     205    keyword obj GAIN     -wf 1.0
     206    keyword obj RDNOISE  -wf $mkiRD
    166207    sprintf name "$simdir/raw/obj.%02d.fits" $i
    167208    wd obj $name
    … …  
    170211  for i 5 10
    171212    $exptime = 3 * ($i - 4)
    172     mkimages.init 1024 1024 5 1000 0 2.0
    173     set flat = blank + sqrt($mkiRD*$mkiRD + $darkrate*$exptime + $skyrate*$exptime)*mki_noisei + $darkrate*$exptime + $skyrate*$exptime + $biaslevel
    174     keyword flat AIRMASS  -wf {1 + 0.1*$i}
    175     keyword flat FILTER   -w  B
    176     keyword flat POSANGLE -wf 0.0
    177     keyword flat RA       -wf 0.0
    178     keyword flat DEC      -wf 0.0
    179     keyword flat EXPTIME  -wf $exptime
    180     keyword flat DARKTIME -wf $exptime
    181     keyword flat DATE-OBS -w "2006/10/01"
     213    $opentime = $exptime - $shutter
     214    mkimages.init 1024 1024 5 1000 0 2.0
     215    set obj = blank + sqrt($mkiRD*$mkiRD + $darkrate*$exptime + $skyrate*$opentime)*mki_noisei + $darkrate*$exptime + $skyrate*$opentime + $biaslevel
     216    keyword obj AIRMASS  -wf {1 + 0.1*$i}
     217    keyword obj FILTER   -w  B
     218    keyword obj POSANGLE -wf 0.0
     219    keyword obj RA       -wf 0.0
     220    keyword obj DEC      -wf 0.0
     221    keyword obj EXPTIME  -wf $exptime
     222    keyword obj DARKTIME -wf $exptime
     223    keyword obj DATE-OBS -w "2006/10/01"
    182224    sprintf time "00:02:%02d" $i
    183     keyword flat UTC-OBS  -w $time
    184     keyword flat CCDSUM   -w "1 1"
    185     keyword flat TELESCOP -w simscope
    186     keyword flat INSTRUME -w SIMPLE
    187     keyword flat OBSTYPE  -w OBJECT
    188     keyword flat OBJECT   -w OBJECT
     225    keyword obj UTC-OBS  -w $time
     226    keyword obj CCDSUM   -w "1 1"
     227    keyword obj TELESCOP -w simscope
     228    keyword obj INSTRUME -w SIMTEST
     229    keyword obj OBSTYPE  -w OBJECT
     230    keyword obj OBJECT   -w OBJECT
     231    keyword obj GAIN     -wf 1.0
     232    keyword obj RDNOISE  -wf $mkiRD
    189233    sprintf name "$simdir/raw/obj.%02d.fits" $i
    190234    wd flat $name
Note: See TracChangeset for help on using the changeset viewer.