Open chukus opened 2 months ago
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.
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.
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.