World Coordinate Systems & Distortion

The simulator has two options for modeling distortion and world coordinate systems. The first is to use the routines in the galsim.roman package; see GalSim’s documentation for more information. The second is to use distortion reference files from the Calibration References Data System (CRDS).

The latter system works by grabbing reference distortion maps for the appropriate detector and filter combinations from the Roman CRDS server. These distortion maps are then wrapped into a galsim WCS object and fed to galsim’s rendering pipeline.

romanisim.wcs Module

Roman WCS interface for galsim.

galsim.roman has an implementation of Roman’s WCS based on some SIP coefficients for each SCA. This is presumably plenty good, but here we take the alternative approach of using the distortion functions provided in CRDS. These naturally are handled by the gWCS library, but galsim only naturally supports astropy WCSes via the ~legacy interface. So this module primarily makes the wrapper that interfaces gWCS and galsim.CelestialWCS together.

This presently gives rather different world coordinates given a specific telescope boresight. Partially this is not doing the same roll_ref determination that galsim.roman does, which could be fixed. But additionally the center of the SCA looks to be in a different place relative to the boresight for galsim.roman than for what I get from CRDS. This bears more investigation.

Functions

convert_wcs_to_gwcs(wcs)	Convert a GalSim WCS object into a GWCS object.
fill_in_parameters(parameters, coord[, ...])	Add WCS info to parameters dictionary.
get_wcs(image[, usecrds, distortion])	Get a WCS object for a given sca or set of CRDS parameters.
make_wcs(targ_pos, distortion[, roll_ref, ...])	Create a gWCS from a target position, a roll, and a distortion map.
wcs_from_fits_header(header)	Convert a FITS WCS to a GWCS.

Classes

GWCS(gwcs[, origin])

This WCS uses gWCS to implent a galsim CelestialWCS.

Class Inheritance Diagram