Cannot properly set vopt values of markups

[chukus]

Hello. Beautiful module. I am slowly learning but I really admire this work.

I have a CO JCMT datacube for a molecular cloud, as well as a list of sources with RA, DEC and radial velocity values.
I used the example listed as a python script (the sofia catalog markup example) and I succesfully loaded the image and the markups, which are correctly displayed in RA, DEC but not in vopt.
I tried to modify the restfreq value in the header from 1.42 GHz (HI) to 115.2 GHz (CO) but then I get the following error...

vtkFITSReader::AllocateWCS: WCSlib failed to create WCSstruct.
World coordinates will not be displayed.

[Punzo]

Hello. Beautiful module. I am slowly learning but I really admire this work.

thanks!

I have a CO JCMT datacube for a molecular cloud, as well as a list of sources with RA, DEC and radial velocity values. I used the example listed as a python script (the sofia catalog markup example) and I succesfully loaded the image and the markups, which are correctly displayed in RA, DEC but not in vopt. I tried to modify the restfreq value in the header from 1.42 GHz (HI) to 115.2 GHz (CO) but then I get the following error...

vtkFITSReader::AllocateWCS: WCSlib failed to create WCSstruct. World coordinates will not be displayed.

uhm if I remember well, in default for SlicerAstro I set to use radio velocities.

Have you tried to switch to optical velocity in the WCS for the volume? i.e. these methods:
https://github.com/Punzo/SlicerAstro/blob/master/AstroVolume/MRML/vtkMRMLAstroVolumeDisplayNode.cxx#L262-L365

from the UI:

If instead you have already set the optical velocity, then it would be a bug either in SlicerAstro or in the data.
For the coordinate transformation, SlicerAstro uses WCSLIB version 7.1 https://github.com/Punzo/wcslib. The latest is 8.3 (https://www.atnf.csiro.au/people/mcalabre/WCS/). You may try the latest version with a Python script and see if you can still reproduce the issue with your data.

[Punzo]

btw: if the WCSlib fails, there should be more logging printed by https://github.com/Punzo/SlicerAstro/blob/master/vtkFits/vtkFITSReader.cxx#L2017-L2800

looking at those would help.