Details
-
Type:
Improvement
-
Status: Done
-
Resolution: Done
-
Fix Version/s: None
-
Component/s: ts_aos
-
Labels:
-
Story Points:2
-
Epic Link:
-
Team:Telescope and Site
Description
Update in `phosim_repackager.py` the repackaged FITS header `ROTANG` , so that the correct WCS is achieved. Currently, one needs to feed `90+/-ROTANG` to `lsst.obs.base.createInitialSkyWcsFromBoresight` to achieve the match. An illustration (attached) shows a simulated `lsstCam` defocal image, with red circles corresponding to locations of stars inferred from the input catalog assuming `90+ROTANG`
minimal code used to make the image is attached- wcs_code.py
Crucial part is
new_wcs = createInitialSkyWcsFromBoresight(pointing_loc,
|
(postIsr_md['ROTANG']+90)*lsst.geom.degrees, |
exposure.getDetector(),
|
flipX=False,)
|
since we desire to be able to use just {{postIsr_md['ROTANG'] }}, without adding 90 degrees to make the input catalog match the WCS information.
Attachments
Activity
This PR looks good to me. Please discuss with Bryce for the possible break in ts_wep for the assumption of 90 degree. Thanks!
Merged and released as v0.2.4, discussed with Bryce the changes (he suggested them originally)
Opened PR https://github.com/lsst-dm/phosim_utils/pull/16