| | 1 | = Astrometry with the IPP = |
| | 2 | |
| | 3 | The following page reviews my process of extracting astrometry information from Megacam images |
| | 4 | |
| | 5 | |
| | 6 | ---- |
| | 7 | === Image Injection === |
| | 8 | |
| | 9 | |
| | 10 | {{{ |
| | 11 | ipp_serial_inject.pl --telescope CFHT --camera MEGACAM --workdir /data/ipp035.0/beaumont/reduce.sa_98/ --label sa.98 -- dbname beaumont |
| | 12 | }}} |
| | 13 | |
| | 14 | |
| | 15 | === Pantasks Macro File === |
| | 16 | I use the following macro file (named input) to control injection, camera, and chip reduction with Pantasks. |
| | 17 | |
| | 18 | {{{ |
| | 19 | macro setup |
| | 20 | |
| | 21 | $MODULES:0 = /home/kiawe/eugene/src/panstarrs/ipp/ippTasks |
| | 22 | module pantasks.pro |
| | 23 | |
| | 24 | # $LABEL = prpics.20081216 |
| | 25 | # $POLLLIMIT = 100 |
| | 26 | |
| | 27 | add.database beaumont |
| | 28 | load.hosts |
| | 29 | end |
| | 30 | |
| | 31 | macro load.hosts |
| | 32 | controller host add ipp023 -threads 4 |
| | 33 | <etc etc> |
| | 34 | end |
| | 35 | |
| | 36 | macro check.books |
| | 37 | book list |
| | 38 | echo "" |
| | 39 | end |
| | 40 | |
| | 41 | macro register |
| | 42 | module register.pro |
| | 43 | end |
| | 44 | |
| | 45 | macro registerOn |
| | 46 | register.on |
| | 47 | end |
| | 48 | |
| | 49 | macro registerOff |
| | 50 | register.off |
| | 51 | end |
| | 52 | |
| | 53 | macro detrend |
| | 54 | module detrend.process.pro |
| | 55 | module detrend.stack.pro |
| | 56 | module detrend.norm.pro |
| | 57 | module detrend.resid.pro |
| | 58 | module detrend.reject.pro |
| | 59 | end |
| | 60 | |
| | 61 | macro chip |
| | 62 | module chip.pro |
| | 63 | end |
| | 64 | |
| | 65 | macro chipOn |
| | 66 | chip.on |
| | 67 | end |
| | 68 | |
| | 69 | macro chipOff |
| | 70 | chip.off |
| | 71 | end |
| | 72 | |
| | 73 | macro camera |
| | 74 | module camera.pro |
| | 75 | end |
| | 76 | |
| | 77 | macro cameraOn |
| | 78 | camera.on |
| | 79 | end |
| | 80 | |
| | 81 | macro cameraOff |
| | 82 | camera.off |
| | 83 | end |
| | 84 | }}} |
| | 85 | |
| | 86 | Then, in pantasks: |
| | 87 | {{{ |
| | 88 | (ipp035:~/data) beaumont% pantasks |
| | 89 | inputno file /home/panstarrs/beaumont/.pantasksrc |
| | 90 | |
| | 91 | Welcome to pantasks - parallel task scheduler |
| | 92 | |
| | 93 | |
| | 94 | pantasks: input input |
| | 95 | pantasks: setup |
| | 96 | starting controller connection |
| | 97 | Connected |
| | 98 | HostID: 1 |
| | 99 | HostID: 2 |
| | 100 | HostID: 3 |
| | 101 | HostID: 4 |
| | 102 | HostID: 5 |
| | 103 | HostID: 6 |
| | 104 | HostID: 7 |
| | 105 | HostID: 8 |
| | 106 | HostID: 9 |
| | 107 | HostID: 10 |
| | 108 | HostID: 11 |
| | 109 | HostID: 12 |
| | 110 | HostID: 13 |
| | 111 | HostID: 14 |
| | 112 | pantasks: register |
| | 113 | pantasks: chip |
| | 114 | pantasks: camera |
| | 115 | pantasks: run |
| | 116 | pantasks: add.label sa.98 |
| | 117 | }}} |
| | 118 | |
| | 119 | === DVO Creation === |
| | 120 | The above commands produce the .smf files which can be injected into DVO. Adding images to a DVO database proceeds as follows: |
| | 121 | |
| | 122 | {{{ |
| | 123 | <cd into working directory> |
| | 124 | set x = `ls */*.smf` # selects all smf files |
| | 125 | foreach a ($x) |
| | 126 | addstar -D CATDIR catdir.98 -D CAMERA megacam -update $a |
| | 127 | end |
| | 128 | }}} |
| | 129 | |
| | 130 | 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 |
| | 131 | |
| | 132 | {{{ |
| | 133 | addstar -D CATDIR catdir.98 -D CAMERA megacam -resort |
| | 134 | }}} |
| | 135 | |
| | 136 | 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. |
| | 137 | |
| | 138 | {{{ |
| | 139 | relphot -D CATDIR catdir.98 -D CAMERA megacam -averages -update -D GRID_TOOFEW 10 |
| | 140 | }}} |
| | 141 | |
| | 142 | === Relastro === |
| | 143 | |
| | 144 | 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): |
| | 145 | |
| | 146 | |
| | 147 | {{{ |
| | 148 | relastro -D CATDIR catdir.98 -D CAMERA megacam -region 0 360 -2 2 |
| | 149 | -D RELASTRO_SIGMA_LIM .1 -update-objects +photflags -update |
| | 150 | relastro -D CATDIR catdir.98 -D CAMERA megacam -region 0 360 -2 2 |
| | 151 | -D RELASTRO_SIGMA_LIM .05 -update-chips +photflags -update -instmag -12 -8 |
| | 152 | }}} |
| | 153 | |
| | 154 | |
| | 155 | I found that repeating these commands 2-3 times converged on a stable solution. |
| | 156 | |
| | 157 | Further discussion on configuration file options, command line options, etc coming soon. |