search

CARTAvis / carta-frontend

Source code repository for the frontend component of CARTA, a new visualization tool designed for the ALMA, the VLA and the SKA pathfinders.
https://cartavis.github.io/
GNU General Public License v3.0
20 stars 10 forks source link

In 4.0, spectral overlays don't work? #2253

Open keflavich opened 1 year ago

keflavich commented 1 year ago

In 4.0b, it is possible to overlay spectra from matched cubes:

image

in 4.0, it doesn't work:

image
keflavich commented 1 year ago

I guess the bug is actually that they all show (hidden), but I don't know why. They are successfully matched.

kswang1029 commented 1 year ago

What are their pixel units, respectively?

keflavich commented 1 year ago

Their brightness units are Jy/beam - is that what you meant? Or are you referring to something else?

kswang1029 commented 1 year ago

Could you share with me all the image headers please? (Via private email is fine if those info is sensitive) please also check the settings dialog and the conversion tab. There is a unit conversion dropdown. In v4.0 we do check intensity units from different cubes to have multiple profiles in one plot.

keflavich commented 1 year ago

Headers look like this (copy-pasted from CARTA, so no formatting...). BUNIT is unset.

BPA = 8.984436624563E+01BEAM = Beam: BMAJ=2.2139292858953508 arcsec BMIN=1.6701402389415263 &CONTINUE  'arcsec BPA = 8.984436624563E+01SLICE = [[(None, None, None), (None, None, 9), (None, None, 9)]]WCSAXES = 3 / Number of coordinate axes CRPIX1 = 1533.7777777778 / Pixel coordinate of reference point CRPIX2 = 624.77777777778 / Pixel coordinate of reference point CRPIX3 = 150 / Pixel coordinate of reference point CDELT1 = -4.999999998870E-04 / [deg] Coordinate increment at reference point CDELT2 = 4.999999998870E-04 / [deg] Coordinate increment at reference point CDELT3 = 1.470155020000E+00 / [m/s] Coordinate increment at reference point CUNIT1 = deg / Units of coordinate increment and value CUNIT2 = deg / Units of coordinate increment and value CUNIT3 = km.s-1 / Units of coordinate increment and value CTYPE1 = GLON-TAN / galactic longitude, gnomonic projection CTYPE2 = GLAT-TAN / galactic latitude, gnomonic projection CTYPE3 = VRAD / Radio velocity (linear) CRVAL1 = 1.426959252752E-01 / [deg] Coordinate value at reference point CRVAL2 = 0.000000000000E+00 / [deg] Coordinate value at reference point CRVAL3 = 0.000000000000E+00 / [m/s] Coordinate value at reference point LONPOLE = 1.800000000000E+02 / [deg] Native longitude of celestial pole LATPOLE = 0.000000000000E+00 / [deg] Native latitude of celestial pole RESTFRQ = 1.000763000000E+11 / [Hz] Line rest frequency MJDREF = 0.000000000000E+00 / [d] MJD of fiducial time SPECSYS = LSRK / Reference frame of spectral coordinates BUNIT =
keflavich commented 1 year ago

Setting BUNIT solved the issue. It would be nice if that could also be set within CARTA, though, and it would be helpful to indicate that's the reason for not showing the data.