import json import numpy as np import pysiaf from pysiaf.utils.rotations import attitude_matrix from astropy.table import Table from astropy.coordinates import SkyCoord import astropy.units as u ECSV = ( "/Users/adam/repos/roman_notebooks/notebooks/" "footprint_visualization/aux_data/roman_gps.sim.ecsv" ) print("Loading Roman SIAF …") rsiaf = pysiaf.Siaf("Roman") wfi_cen = rsiaf['WFI_CEN'] V2Ref = wfi_cen.V2Ref # +1546.38" (WFI focal-plane offset from V3 boresight) V3Ref = wfi_cen.V3Ref # -892.79" sensors = [rsiaf[f"WFI{j:02d}_FULL"] for j in range(1, 19)] print("Reading roman_gps.sim.ecsv …") t = Table.read(ECSV) def unique_pointings(target_name): sub = t[t["TARGET_NAME"] == target_name] pts = {} for r in sub: key = (round(float(r["RA"]), 6), round(float(r["DEC"]), 6), round(float(r["PA"]), 3)) pts[key] = True return sorted(pts.keys()) def all_sca_polygons(pointings): """Return flat list of all SCA corner polygons for a list of pointings.""" polys = [] for ra, dec, pa in pointings: attmat = attitude_matrix(V2Ref, V3Ref, ra, dec, pa) for sensor in sensors: sensor.set_attitude_matrix(attmat) c = sensor.corners("sky") # (ra_4pts, dec_4pts) polys.append([[round(float(c[0][i]), 5), round(float(c[1][i]), 5)] for i in range(4)]) return polys print("Computing Deep+Spec_W51 …") w51_pts = unique_pointings("Deep+Spec_W51") deep_polys = all_sca_polygons(w51_pts) print(f" {len(deep_polys)} SCA polygons") print("Computing Disk7 (may take ~30 s) …") d7_pts = unique_pointings("Disk7") disk7_polys = all_sca_polygons(d7_pts) print(f" {len(disk7_polys)} SCA polygons") deep_poly_js = json.dumps(deep_polys, separators=(',', ':')) disk7_poly_js = json.dumps(disk7_polys, separators=(',', ':')) html = """ Roman RGPS Footprint — W51 — Aladin Lite

Roman RGPS — W51 Footprint

Background survey
Image overlay (click to toggle)
Legend
Deep_W51 (l=49.4°, b=−0.2°)
Disk7 region (l=29–50.5°)
Image overlay (active)
Source: Full_Roman_Target_list.csv
(github.com/rachel3834/rgps)
Real WFI SCA footprints via pysiaf
Source: roman_gps.sim.ecsv (APT)
Roman RGPS · pysiaf WFI geometry · Deep+Spec_W51 (8 pointings) · Disk7 (632 pointings) · Source: roman_gps.sim.ecsv
""" with open("/Users/adam/Downloads/roman_footprint_rgps_w51.html", "w") as f: f.write(html) print("Written roman_footprint_rgps_w51.html")