IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Dec 3, 2018, 6:23:32 AM (8 years ago)
Author:
eugene
Message:

updates to the tessellation descriptions

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/doc/release.2015/ps1.datasystem/datasystem.tex

    r40560 r40563  
    782782grid laid out on the chips in the camera to a system of pixels with
    783783consistent geometry for a location on the sky.  The new image
    784 coordinate system is defined by one of a number ``tessellations'' of
    785 the sky.  Each tessellation specifies the way in which the sky is
    786 broken into individual images by defining a collection of projection
    787 centers.  For each projection center, positions on the sky are
    788 transformed to image pixels via a projection with a specified pixel
    789 scale and rotation.  In general, the pixel grid within the projection
    790 is defined as a simplified grid with the y-axis aligned to the
    791 Declination lines and no distortion terms.  The projection centers are
    792 typically separated by several degrees on the sky; for pixel scales
    793 appropriate to GPC1, the resulting collection of pixels would be
    794 unwieldy in terms of memory in the processing computer.  The pixel
    795 grid is thus subdivided into smaller sub-images called 'skycells'. 
    796 
    797 moves the data from a given exposure beyond
    798 away from being camera specific and towards a uniform sky oriented
    799 arrangement.  There are a number of ``tessellations'' defined and used
    800 by the IPP to define the extent and scaling of images on this uniform
    801 arrangement.  A tessellation can be defined for a limited region, such
    802 as M31 or other fields of particular interest that can be well
    803 described by a single tangent plane projection, or for larger regions
    804 which have multiple projection centers.  For the $3\pi$ survey, the
    805 \ippmisc{RINGS.V3} tessellation was used that arrange projection centers
    806 spaced every four degrees in both RA and DEC, with $0\farcs{}25$
    807 pixels.  These projections are further broken down into ``skycells''
    808 that form a $10\times{}10$ grid within the projection, with an overlap
    809 region of 60\arcsec\ between adjacent skycells to ensure that objects are not
    810 split on all images.
    811 
    812 These tessellations are stored in the DVO format, with
    813 \ippdbtable{SkyTable} entries defining the projection centers and
    814 image boundaries for all the skycells.  The first step of the
    815 \ippstage{warp} stage is determining which skycells overlap with the
    816 input exposure.  These overlaps are determined by the
    817 \ippprog{dvoImageOverlaps} program, which compares the astrometrically
    818 calibrated catalog from the \ippstage{camera} stage to the
    819 \ippdbtable{SkyTable} entries.  The output of this command is used to
    820 populate the \ippdbtable{warpSkyCellMap} table in the database, which
    821 contains a row for each skycell and OTA that overlap.  This results in
    822 more rows than there are OTAs, as each skycell can contain
     784coordinate system is defined by one of a number of ``tessellations''
     785which specify how the sky is divided into individual images.  A single
     786tessellation starts with a collection of projection centers
     787distributed across the sky.  A grid of image pixels about each
     788projection center corresponds to sky positions via a projection with a
     789specified pixel scale and rotation.  In general, the pixel grid within
     790the projection is defined as a simplified grid with the y-axis aligned
     791to the Declination lines and no distortion terms.  The projection
     792centers are typically separated by several degrees on the sky; for
     793pixel scales appropriate to GPC1, the resulting collection of pixels
     794would be unwieldy in terms of memory in the processing computer.  The
     795pixel grid is thus subdivided into smaller sub-images called
     796'skycells'.
     797
     798A tessellation can be defined for a limited region, with only a small
     799number of projection centers (e.g., for processing the M31 region), or
     800even a single projection center (e.g., for the Medium Deep fields).
     801For the $3\pi$ survey, the tessellation contains projection centers
     802covering the entire sky.  The version used to for the PV3 analysis is
     803called the \ippmisc{RINGS.V3}.  In this tessellation, projection
     804centers are spaced every four degrees in DEC and the RA spacing is
     805approximately four degrees as well, adjusted to ensure an integer
     806number of equal-sized regions.  \ippmisc{RINGS.V3} uses a pixel scale
     807of $0\farcs{}25$ per pixel.  The projections subdivided into a
     808$10\times{}10$ grid of skycells, with an overlap region of
     80960\arcsec\ between adjacent skycells to ensure that objects of modest
     810size are not split on all images.  The coordinate system used for
     811these images matches the parity of the sky, with north in the positive
     812$y$ direction and east to the negative $x$ direction.  The
     813tessellations used by the IPP are stored in the DVO format (see
     814Section~\ref{sec:DVO}), with \ippdbtable{SkyTable} entries defining
     815the projection centers and image boundaries for all skycells.
     816
     817The first step of the \ippstage{warp} stage is to determine which
     818skycells overlap with the input exposure.  These overlaps are
     819determined by the \ippprog{dvoImageOverlaps} program, which compares
     820the astrometrically calibrated catalog from the \ippstage{camera}
     821stage to the DVO database defining the target tessellation.  The
     822output of this command is used to populate the
     823\ippdbtable{warpSkyCellMap} table in the database, which contains a
     824row for each skycell and OTA that overlap.  Each skycell may contain
    823825contributions from multiple OTAs.
    824826
    825 Once this mapping has been defined, jobs to construct each skycell are
    826 run, passing the \ippstage{camera} stage catalog and the
     827Once this mapping has been defined, jobs to warp the pixels onto each
     828skycell are run, passing the \ippstage{camera} stage catalog and the
    827829\ippstage{chip} stage images (including the variance images and the
    828830updated masks) to the \ippprog{pswarp} program.  For details on the
     
    830832are the geometrically transformed images containing all input pixels
    831833warped to the common skycell pixel grid, which can subsequently be
    832 used for stacking and difference image analysis.  The image, mask, and
    833 variance generated at this stage will be available from the image
    834 extraction tools at the MAST archive at STScI as part of the DR2 data
    835 release.  A catalog is also generated containing the locations of
    836 sources from the input catalog that fall within area of the
    837 \ippstage{warp}.
    838 
    839 When the jobs have completed, an entry for the skycell is added to the
    840 \ippdbtable{warpSkyfile} database table, linked to the
     834used for stacking and difference image analysis.  For the $3\pi$
     835survey data, the signal, mask, and variance images generated at this
     836stage are being made available from the image extraction tools at the
     837MAST archive at STScI as part of the DR2 data release. 
     838
     839%% A catalog is
     840%% also generated containing the locations of sources from the input
     841%% catalog that fall within area of the \ippstage{warp}.
     842
     843When the \ippstage{warp} jobs have completed, an entry for the skycell
     844is added to the \ippdbtable{warpSkyfile} database table, linked to the
    841845\ippdbtable{warpRun} entry by a common \ippdbcolumn{warp_id}.  An
    842846\ippmisc{advance} task again checks that all potential skycells have
    843847been generated.  At this point, the direct promotion of exposures from
    844848one stage to the next stops, as the logic for matching exposures for
    845 combination is more complicated than simply adding a single entry (as
    846 discussed above).
     849other combinations is more complicated than simply adding a single
     850entry.
    847851
    848852\subsection{Stack Combination}
Note: See TracChangeset for help on using the changeset viewer.