| 1 |
25 |
gam |
\chapter{Handling of \index{image} image data}
|
| 2 |
40 |
bertin |
|
| 3 |
|
|
\section{Image format}
|
| 4 |
25 |
gam |
{\sc SExtractor} accepts \index{image} images stored in FITS\footnote{\it Flexible Image Transport System}
|
| 5 |
19 |
gam |
format (Wells \etal 1981, see also {\tt http://fits.gsfc.nasa.gov}).
|
| 6 |
|
|
Both ``Basic FITS'' (one single header and one single body) and ``Multi-Extension-FITS'' (MEF)
|
| 7 |
25 |
gam |
\index{image} images are recognized. Binary {\sc SExtractor} catalogues produced from MEF \index{image} images are MEF files
|
| 8 |
19 |
gam |
themselves. If catalogue output is in ASCII format, all catalogues from the individual extensions
|
| 9 |
|
|
are concatenated in one big file; the {\tt EXT\_NUMBER} catalogue parameter must be used to tell
|
| 10 |
|
|
which extension the source belongs to.
|
| 11 |
|
|
|
| 12 |
25 |
gam |
For \index{image} images with ${\tt NAXIS} > 2$, only the first data-plane is loaded.
|
| 13 |
|
|
If \index{WCS} WCS\footnote{\it World Coordinate System} information (Greisen \& Calabretta 1995,
|
| 14 |
19 |
gam |
{\tt http://www.cv.nrao.edu/fits/documents/wcs/wcs.all.ps}) is available in the
|
| 15 |
|
|
header, it is automatically used by {\sc SExtractor} to compute astrometric parameters. Other
|
| 16 |
|
|
astrometric descriptions
|
| 17 |
|
|
like AST ({\it Starlink} format) \gam{Provide URL in footnote} or the solution coefficients of the DSS
|
| 18 |
|
|
\footnote{{\it Digital Sky Survey} \gam{URL in footnote}} plates are not recognized by the software.
|
| 19 |
|
|
|
| 20 |
|
|
In {\sc SExtractor}, as in all similar programs, FITS axis ``1'' is traditionaly referred as the
|
| 21 |
|
|
{\tt X} axis, and FITS axis ``2'' as the {\tt Y} axis.
|
| 22 |
|
|
|
| 23 |
40 |
bertin |
\section{Double image mode}
|
| 24 |
|
|
If data are available for the same field in several photometric bands,
|
| 25 |
|
|
it is usually desirable to measure object characteristics such as
|
| 26 |
|
|
magnitudes exactly at the same positions and through the same apertures
|
| 27 |
|
|
for the different bands, to derive precise colour indices for example.
|
| 28 |
|
|
{\sc SExtractor} makes this possible by providing a special mode dubbed
|
| 29 |
|
|
``double image mode'' where detections are made on one image while
|
| 30 |
|
|
measurements are carried out on another (both images must have exactly the
|
| 31 |
|
|
same number of pixels). By repeatedly running {\sc SExtractor}
|
| 32 |
|
|
with various ``measurement images'' while keeping the same ``detection image'',
|
| 33 |
|
|
one ends up with a set of catalogues having the same sources measured through
|
| 34 |
42 |
bertin |
different channels. The detection image will generally be chosen in the band
|
| 35 |
|
|
where the data are the deepest. Alternatively, using a ``$\chi2$ image''
|
| 36 |
|
|
\citep{szalay:al:1999} (produced e.g., by
|
| 37 |
|
|
\href{http://astromatic.net/software/swarp}{\sc SWarp}) as a detection image,
|
| 38 |
|
|
will allow most of the sources present in at least one channel to be detected
|
| 39 |
|
|
and measured.
|
| 40 |
40 |
bertin |
|
| 41 |
42 |
bertin |
Double image mode is automatically engaged by providing {\sc SExtractor} with
|
| 42 |
|
|
two images instead of one:
|
| 43 |
40 |
bertin |
|
| 44 |
42 |
bertin |
\% {\tt sex detection.fits,measurement.fits}\\
|
| 45 |
40 |
bertin |
|
| 46 |
42 |
bertin |
A space may be used instead of a coma. In the example above,
|
| 47 |
|
|
{\tt sex detection.fits} is used as a detection image, while measurements are
|
| 48 |
|
|
carried out on {\tt measurement.fits}.
|
| 49 |
40 |
bertin |
|
| 50 |
|
|
You then obtain two different catalogues for the two images, which
|
| 51 |
|
|
have exactly the same numbers of lines and can be put together.
|
| 52 |
|
|
|
| 53 |
|
|
If you have images in other bands covering the same regions, you can
|
| 54 |
|
|
measure them in double image mode as well, always taking as a
|
| 55 |
|
|
reference image 1, and then combine all the catalogues obtained.
|
| 56 |
|
|
|