def miri_rate_zero_crossing():
xsize = 1032
ysize = 1024
shape = (ysize, xsize)
im = ImageModel(shape)
im.var_rnoise = np.random.random(shape)
im.meta.wcsinfo = {
'dec_ref': 2.16444343946559e-05,
'ra_ref': -0.00026031780056776,
'roll_ref': 0.0,
'v2_ref': -415.0690466121227,
'v3_ref': -400.575920398547,
'v3yangle': 0.0,
'vparity': -1
}
im.meta.instrument = {
'detector': 'MIRIMAGE',
'filter': 'P750L',
'name': 'MIRI'
}
im.meta.observation = {'date': '2019-01-01', 'time': '17:00:00'}
im.meta.subarray = {
'fastaxis': 1,
'name': 'FULL',
'slowaxis': 2,
'xsize': xsize,
'xstart': 1,
'ysize': ysize,
'ystart': 1
}
im.meta.exposure = {
'duration': 11.805952,
'end_time': 58119.85416,
'exposure_time': 11.776,
'frame_time': 0.11776,
'group_time': 0.11776,
'groupgap': 0,
'integration_time': 11.776,
'nframes': 1,
'ngroups': 100,
'nints': 1,
'nresets_between_ints': 0,
'nsamples': 1,
'readpatt': 'FAST',
'sample_time': 10.0,
'start_time': 58119.8333,
'type': 'MIR_LRS-FIXEDSLIT',
'zero_frame': False
}
return im
im.meta.wcsinfo.roll_ref = header_slp['ROLL_REF']
references = {
"area": "crds/jwst_nircam_area_0017.fits",
"distortion": "crds/jwst_nircam_distortion_0093.asdf",
"drizpars": "crds/jwst_nircam_drizpars_0001.fits",
"flat": "crds/jwst_nircam_flat_0337.fits",
"photom": "crds/jwst_nircam_photom_0074.fits"
}
load_wcs(im, references)
im.data = hdul_rate['SCI'].data
im.err = hdul_rate['ERR'].data
im.dq = hdul_rate['DQ'].data
im.var_poisson = hdul_rate['VAR_POISSON'].data
im.var_rnoise = hdul_rate['VAR_RNOISE'].data
im.meta.target.proper_motion_epoch = "2000"
im.meta.dither.primary_type = "IMAGING"
ra_ref = header_slp['RA_REF']
dec_ref = header_slp['DEC_REF']
crval1_ra = header_slp['CRVAL1']
crval2_dec = header_slp['CRVAL2']
targ_ra = header_slp['TARG_RA']
targ_dec = header_slp['TARG_DEC']
c_ref = SkyCoord(ra=ra_ref * u.degree, dec=dec_ref * u.degree, frame='fk5')
c_crval = SkyCoord(ra=crval1_ra * u.degree,
dec=crval2_dec * u.degree,
frame='fk5')
