def test_evaluate_fk5_model():
geom = WcsGeom.create(width=(5, 5), binsz=0.1, frame="icrs")
model = GaussianSpatialModel(
lon_0="0 deg", lat_0="0 deg", sigma="0.1 deg", frame="fk5"
)
data = model.evaluate_geom(geom)
assert data.sum() > 0
def test_evaluate_on_fk5_map():
# Check if spatial model can be evaluated on a map with FK5 frame
# Regression test for GH-2402
header = {}
header["CDELT1"] = 1.0
header["CDELT2"] = 1.0
header["CTYPE1"] = "RA---TAN"
header["CTYPE2"] = "DEC--TAN"
header["RADESYS"] = "FK5"
header["CRVAL1"] = 0
header["CRVAL2"] = 0
header["CRPIX1"] = 5
header["CRPIX2"] = 5
wcs = WCS(header)
geom = WcsGeom(wcs, npix=(10, 10))
model = GaussianSpatialModel(lon_0="0 deg", lat_0="0 deg", sigma="1 deg")
data = model.evaluate_geom(geom)
assert data.sum() > 0
def test_evaluate_fk5_model():
geom = WcsGeom.create(width=(5, 5), binsz=0.1, coordsys="CEL")
model = GaussianSpatialModel("0 deg", "0 deg", "0.1 deg", frame="fk5")
data = model.evaluate_geom(geom)
assert data.sum() > 0
