usage: romanisim-gaia-catalog [-h] [--radius RADIUS] ra dec time filename

Generate a romanisim input catalog from Gaia.

positional arguments:
  ra               right ascension (degree)
  dec              declination (degree)
  time             time (ISO 8601)
  filename         output catalog

options:
  -h, --help       show this help message and exit
  --radius RADIUS  find sources within this radius of query (default: None)

EXAMPLE: romanisim-gaia-catalog 270.0 66.0 2026-01-01T00:00:00 gaia.ecsv

romanisim-gaia-catalog 6 -72 --radius 0.5 2016-01-01T00:00:00 gaia-catalog_ngc104.ecsv

usage: romanisim-make-image [-h] [--bandpass BANDPASS] [--boresight] [--catalog CATALOG] [--config CONFIG] [--date DATE] [--level LEVEL] [--ma_table_number MA_TABLE_NUMBER] [--nobj NOBJ] [--previous PREVIOUS] [--radec RADEC RADEC] [--rng_seed RNG_SEED] [--roll ROLL] [--sca SCA]
                            [--usecrds] [--psftype {epsf,galsim,stpsf}] [--webbpsf] [--stpsf] [--truncate TRUNCATE] [--pretend-spectral PRETEND_SPECTRAL] [--drop-extra-dq] [--scale-factor SCALE_FACTOR] [--extra-counts EXTRA_COUNTS [EXTRA_COUNTS ...]]
                            filename

Make a demo image.

positional arguments:
  filename              output image (asdf). {} and {bandpass} strings will be automatically populated with detector and bandpass information.

options:
  -h, --help            show this help message and exit
  --bandpass BANDPASS   bandpass to simulate (default: F087)
  --boresight           radec specifies location of boresight, not center of WFI. (default: False)
  --catalog CATALOG     input catalog (ecsv) (default: None)
  --config CONFIG       input parameter override file (yaml) (default: None)
  --date DATE           UTC Date and Time of observation to simulate in ISOT format. (default: 2026-01-01T00:00:00.000)
  --level LEVEL         1 or 2, for L1 or L2 output (default: 2)
  --ma_table_number MA_TABLE_NUMBER
  --nobj NOBJ
  --previous PREVIOUS   previous simulated file in chronological order used for persistence modeling. (default: None)
  --radec RADEC RADEC   ra and dec (deg) (default: None)
  --rng_seed RNG_SEED
  --roll ROLL           Position angle (North towards YIdl) measured at the V2Ref/V3Ref of the aperture used. (default: 0)
  --sca SCA             SCA to simulate. Use -1 to generate images for all SCAs; include {} in filename for this mode to indicate where the detector number should be filled, e.g. l1_{}.asdf (default: 7)
  --usecrds             Use CRDS references. (default: False)
  --psftype {epsf,galsim,stpsf}
                        Type of PSF generator to use. If None and --usecrds, then "epsf" will be used otherwise "galsim" is used. (default: galsim)
  --webbpsf             Use stpsf for PSF (deprecated, use `--psftype stpsf`) (default: False)
  --stpsf               Use stpsf for PSF (deprecated, use `--psftype stpsf`) (default: False)
  --truncate TRUNCATE   If set, truncate the MA table at given number of resultants. (default: None)
  --pretend-spectral PRETEND_SPECTRAL
                        Pretend the image is spectral. exposure.type and instrument.element are updated to be grism / prism. (default: None)
  --drop-extra-dq       Do not store the optional simulated dq array. (default: False)
  --scale-factor SCALE_FACTOR
                        Velocity aberration-induced scale factor. If negative, use given time to calculated based on orbit ephemeris. (default: -1.0)
  --extra-counts EXTRA_COUNTS [EXTRA_COUNTS ...]
                        An optional FITS file to read to get an array of counts to add into the simulated image. Useful for wrapping idealized images. If 2 arguments are sent in, then the second argument is assumed to be the HDU to use (default=0). (default: None)

EXAMPLE: romanisim-make-image output_image.asdf

romanisim-make-image --bandpass F158 --catalog gaia-catalog_ngc104.ecsv --level 2 --ma_table_number 1002 --radec 6 -72 --roll 0 --sca 1 --usecrds --psftype stpsf example_ngc104.asdf

import roman_datamodels as rdm

from astropy.io import fits
from astropy import visualization
from astropy import units as u

import matplotlib
import matplotlib.pyplot as plt
from matplotlib import colors
from matplotlib import colormaps as cm

import numpy as np

filename = 'example_ngc104.asdf'

with rdm.open(filename) as file:
    wfi = file.meta.instrument.detector # we also extract the string corresponding to the WFI detector that has been simulated
    image = np.array(file.data)
    wcs = file.meta.wcs

# font settings
matplotlib.rcParams['mathtext.fontset'] = 'stix'
matplotlib.rcParams['font.family'] = 'STIXGeneral'

# color-map settings
cmap = cm['YlOrBr_r']

# figure settings
fig, ax = plt.subplots(figsize=(8, 6), subplot_kw={'projection':wcs})
ax.set_title(f"Roman I-Sim {wfi}")
ax.imshow(image, cmap=cm["YlOrBr_r"], origin="lower", norm=colors.LogNorm(0.1, 200.0))
ax.coords.grid(linestyle='dashed', color='grey', alpha=0.5)

# axes settings
lon = ax.coords[0]
lat = ax.coords[1]
lon.set_axislabel('R.A. (degrees)')
lat.set_axislabel('Dec. (degrees)')
lon.set_major_formatter('d.dd')
lat.set_major_formatter('d.dd')
lon.set_ticks(spacing=2 * u.arcminute)
lat.set_ticks(spacing=2 * u.arcminute)

plt.show()