Add moment map generation information to FITS header

[markccchiang]

Closes #1024. Adding some information to the FITS header history about moment map generation. As shown below:

[kswang1029]

@markccchiang quick questions, does this new info exist in the header of the in-memory moment images (the image auto appended to carta after the moment image is generated)? I am not sure where I can locate this new info. 🤔

I suggest we refine the extra moment generation info as the following:

1. add a string like CARTA MOMENT MAP GENERATOR LOG before add extra info from CARTA
2. the extra info should include
   a. the input image for moment map generation
   b. region info (pixel coordinate should be fine)
   c. spectral range (channel, freq, and velocity). This is included already. 👍
   d. mask info

[markccchiang]

@kswang1029 as for the moment region info about pixel coordinates. Does it mean the control points of regions (e.g., point, rectangle, polygon..)?

[kswang1029]

@kswang1029 as for the moment region info about pixel coordinates. Does it mean the control points of regions (e.g., point, rectangle, polygon..)?

yes.

[markccchiang]

Now the new moment map generation info from the header of moment images would look like this:

[kswang1029]

Now the new moment map generation info from the header of moment images would look like this:

this works nicely! Could you add region type too? Also we need rotation angle be part of the info.

HISTORY = / Region control points: [402.9424, 374.0267], [242.3408, 294.6644]

I suggest the format like this (please refer to the header of the exported tsv file from the spectral profiler)

rotbox[[402.942341pix, 374.026678pix], [242.340792pix, 294.664372pix], 0.000000deg]

[kswang1029]

👍

[kswang1029]

@markccchiang just realized that the polarization info is missing. Could you add it too?

[kswang1029]

thanks for adding polarization info. It works with usual Stokes. Now CARTA supports computed polarization quantities too such as polarized intensity, angle, etc, so when these are selected, the displayed text is unknown.

HISTORY = / Stokes: Unknown

I suggest we call show this as

Polarization: Stokes Q
Polarization: Fractional polarized intensity
Polarization: XY

[markccchiang]

I suggest we call show this as

Polarization: Stokes Q
Polarization: Fractional polarized intensity
Polarization: XY

Please confirm are the full name of the computed stokes correct?

Ptotal: Total linearly polarized intensity
Plinear: Linearly polarized intensity
PFtotal: Total fractional linear polarization
PFlinear: Fractional linear polarization
Pangle: Linearly polarized angle

[kswang1029]

@markccchiang please consider

Ptotal: total polarization intensity
Plinear: linear polarization intensity
PFtotal: fractional total polarization intensity
PFlinear: fractional linear polarization intensity
Pangle: polarization angle

[kswang1029]

looking good! 👍