Astrometry with the IPP

The following page reviews my process of extracting astrometry information from Megacam images

---

Image Injection

ipp_serial_inject.pl --telescope CFHT --camera MEGACAM --workdir /data/ipp035.0/beaumont/reduce.sa_98/ --label sa.98 -- dbname beaumont

Pantasks Macro File

I use the following macro file (named input) to control injection, camera, and chip reduction with Pantasks.

macro setup

  $MODULES:0 = /home/kiawe/eugene/src/panstarrs/ipp/ippTasks
   module pantasks.pro

  # $LABEL = prpics.20081216
  # $POLLLIMIT = 100

  add.database beaumont
  load.hosts
end

macro load.hosts
  controller host add ipp023 -threads 4
  <etc etc>
end

macro check.books
  book list
  echo ""
end

macro register
 module register.pro
end

macro registerOn
 register.on
end

macro registerOff
  register.off
end

macro detrend
 module detrend.process.pro
 module detrend.stack.pro
 module detrend.norm.pro
 module detrend.resid.pro
 module detrend.reject.pro
end

macro chip
 module chip.pro
end

macro chipOn
 chip.on
end

macro chipOff
 chip.off
end

macro camera
 module camera.pro
end

macro cameraOn
 camera.on
end

macro cameraOff
 camera.off
end

Then, in pantasks:

(ipp035:~/data) beaumont% pantasks
inputno file /home/panstarrs/beaumont/.pantasksrc

Welcome to pantasks - parallel task scheduler


pantasks: input input
pantasks: setup
starting controller connection
Connected
HostID: 1
HostID: 2
HostID: 3
HostID: 4
HostID: 5
HostID: 6
HostID: 7
HostID: 8
HostID: 9
HostID: 10
HostID: 11
HostID: 12
HostID: 13
HostID: 14
pantasks: register
pantasks: chip
pantasks: camera
pantasks: run
pantasks: add.label sa.98

DVO Creation

The above commands produce the .smf files which can be injected into DVO. Adding images to a DVO database proceeds as follows:

<cd into working directory>
set x = `ls */*.smf` # selects all smf files
foreach a ($x)
  addstar -D CATDIR catdir.98 -D CAMERA megacam -update $a
end

Running addstar with updates adds new entries to the end of the DVO database. At the end of this process, running addstar with -resort re-arranges the database so that measurements of common objects appear consecutively

addstar -D CATDIR catdir.98 -D CAMERA megacam -resort

Next, I ran relphot to populate the secfilt DVO table. Note that this probably does not produce very good photometry measurements; this is just a workaround so that relastro can access the information it needs from the secfilt tables.

relphot -D CATDIR catdir.98 -D CAMERA megacam -averages -update -D GRID_TOOFEW 10

Relastro

Relastro has several different modes. I iterated between two of them. Update-objects updates values in the AVERAGE DVO tables by taking the weighted mean of each object's measurements. Update-chips compares, on a chip-by-chip basis, the measured stellar positions to the average positions. It uses this information to update the astrometric solution (in essence, it resolves for the polynomial to map focal plane coordinates to sky coordinates):

relastro -D CATDIR catdir.98 -D CAMERA megacam -region 0 360 -2 2
    -D RELASTRO_SIGMA_LIM .1 -update-objects +photflags -update
relastro -D CATDIR catdir.98 -D CAMERA megacam -region 0 360 -2 2
    -D RELASTRO_SIGMA_LIM .05 -update-chips +photflags -update -instmag -12 -8

I found that repeating these commands 2-3 times converged on a stable solution.

Further discussion on configuration file options, command line options, etc coming soon.