| | 53 | |
| | 54 | ------ |
| | 55 | |
| | 56 | Message from EAM regarding pswarp modifications: |
| | 57 | |
| | 58 | i'm about to check a set of changes into the trunk related to pswarp. the changes are driven by the goal of warping stack skycells back onto the exposures. to do this, I needed to make some significant changes to the pswarp structure. I hopefully have cleaned up the code organization somewhat to make it easier to read in the process. |
| | 59 | |
| | 60 | Two major changes relate to the input and output data sources. |
| | 61 | |
| | 62 | First, the output target may now be any of our camera formats : in theory, it should be able to warp to any chip-level multi-extension or even multi-file output. the output astrometry is defined by a cmf/smf style file, or any other image with WCS elements. If the config system does not recognize the format of the output header, it will define it to be a skycell for the input camera. |
| | 63 | |
| | 64 | Second, multiple input images may now be specified. It has always been possible to warp from multiple inputs of the same exposure, but now I've added the option to warp from multiple 'exposures'. This is necessary because all skycells are defined to be equivalent to an 'exposure', and in order to populate more than a small fraction of an input exposure, it is necessary to supply at least a few skycells. The multiple input exposures are specified by using an input mdc file (like ppStack or psphotStack) and identifying it to pswarp with "-input foo.mdc". The -file / -list options are still valid, but mutually exclusive with -input. |
| | 65 | |
| | 66 | Note that in multiple exposure mode, the background model and sources are not currently transformed as there was too much ambiguity for me to deal with at this time. Also, the covariances and averaged metadata may more be calculated correctly (though the output WCS astrometry *is* calculated and populated correctly). |
| | 67 | |
| | 68 | Also note that in mutliple exposure mode, there is no attempt to rationally combine information from different exposures into the same pixels. that starts to look like ppStack, and is beyond the scope of this work. Instead, I used the simple rule that the first non-NAN value generated for an output pixel is kept. |
| | 69 | |
| | 70 | To make these changes, and to be efficient, I changed the way the code determines the transformations. It now first loads the astrometric header information of both source and target, determines which of the output chips will be overlapped by any of the input chips, then generates only those output pixels. Only as it does the warping process does it read in the input pixels and input sources. After all inputs have been warped, it then generates the output PSF model. |
| | 71 | |
| | 72 | I've put together a set of tests and minimal test scripts at ~ipp/ippc18/warptest.20130407. |
| | 73 | |