Bug #2866

RADECSYS keyword missing in skymaps/cubes created by GammaLib

Added by Tibaldo Luigi over 5 years ago. Updated over 5 years ago.

Status:ClosedStart date:04/08/2019
Priority:NormalDue date:
Assigned To:Knödlseder Jürgen% Done:

100%

Category:-
Target version:1.6.0
Duration:

Description

This can create problems when read using external software such as ds9.

gamma_cygni_maps.png (405 KB) Tibaldo Luigi, 04/08/2019 05:43 PM

Crab_50h_93600s_MAPCUBE_Gal.png (87.1 KB) Knödlseder Jürgen, 04/09/2019 09:29 AM

comparison.png (196 KB) Knödlseder Jürgen, 04/09/2019 09:31 AM

Gamma_cygni_maps Crab_50h_93600s_mapcube_gal Comparison

Recurrence

No recurrence.

History

#1 Updated by Tibaldo Luigi over 5 years ago

Issue found by Giovanni Piano when simulating gamma Cygni

Maps 1 and 3 in celestial coordinates show a shift w.r.t. the simulated disk (green). Maps in galactic coordinates (2 e 4) match expectations.

#2 Updated by Knödlseder Jürgen over 5 years ago

  • Status changed from New to Pull request
  • % Done changed from 0 to 100

I added some code to the GWcs class that sets the m_radesys and m_equinox members if they were not yet set. This leads now to the attachment of the RADESYS=ICRS keyword to the FITS header. Hope this solves the issue.

#3 Updated by Knödlseder Jürgen over 5 years ago

Here a comparison of the map generate with ctmapcube to the input map. Top is the original comparison by Giovanni Piano, bottom is a comparison I did after changing the code. Both images, displayed with ds9, now match.

#4 Updated by Knödlseder Jürgen over 5 years ago

  • Status changed from Feedback to Closed

No feedback probably means that issue was solved

Close issue now.

Also available in: Atom PDF