Precisely how you extract world coordinate system (WCS) information from a dataset obviously depends on what type of dataset it is. Usually, however, you should be able to obtain a set of FITS header cards which contain the WCS information (and probably much more besides). Suppose that ``cards'' is a pointer to an array of strings containing a complete set of FITS header cards and ``ncard'' is the number of cards. Then proceed as follows:
AstFitsChan *fitschan; AstFrameSet *wcsinfo; char *cards[ MAXCARD ]; int icard, ncard; ... /* Create a FitsChan and fill it with FITS header cards. */ fitschan = astFitsChan( NULL, NULL, "" ); for ( icard = 0; icard < ncard; icard++ ) astPutFits( fitschan, cards[ icard ], 0 ); /* Rewind the FitsChan and read WCS information from it. */ astClear( fitschan, "Card" ); wcsinfo = astRead( fitschan );
The result should be a pointer, ``wcsinfo'', to a FrameSet which contains the WCS information. This pointer can now be used to perform many useful tasks, some of which are illustrated in the following recipes.
Some datasets which do not easily yield FITS header cards may require
a different approach, possibly involving use of a Channel
() rather than a FitsChan. In the case of the
Starlink NDF data format, for example, all the above may be replaced
by a single call to the function
ndfGtwcs--see SUN/33. The
whole process can probably be encapsulated in a similar way for
most data systems, whether they use FITS header cards or not.
For more details about reading WCS information from datasets, see
and
. For a more general description of
FitsChans and their use with FITS header cards, see
and
. For more details
about FrameSets, see
and
.
AST A Library for Handling World Coordinate Systems in Astronomy