diff --git a/changelog/112.doc.rst b/changelog/112.doc.rst new file mode 100644 index 0000000..de914cc --- /dev/null +++ b/changelog/112.doc.rst @@ -0,0 +1 @@ +Correct calculation of the HPC coordinate of the center of the STIX pointing in the imaging demo. diff --git a/examples/imaging_demo.py b/examples/imaging_demo.py index ac120df..c5b0ed1 100644 --- a/examples/imaging_demo.py +++ b/examples/imaging_demo.py @@ -129,7 +129,7 @@ header = make_fitswcs_header(bp_image, coord, telescope='STIX', observatory='Solar Orbiter', scale=[10,10]*u.arcsec/u.pix) fd_bp_map = Map((bp_image, header)) -hpc_ref = Helioprojective(pointing[0], pointing[1], observer=solo, obstime=fd_bp_map.date) +hpc_ref = coord.transform_to(Helioprojective(observer=solo, obstime=fd_bp_map.date)) # Center of STIX pointing in HPC header_hp = make_fitswcs_header(bp_image, hpc_ref, scale=[10, 10]*u.arcsec/u.pix, rotation_angle=90*u.deg+roll) hp_map = Map((bp_image, header_hp)) hp_map_rotated = hp_map.rotate()