I have been dabbling with SimpleITK in python for a while now, and just now I ran in to an issue I cannot find the answer to. From the documentation it is possible to create images with complex type (sitkComplexFloat32, sitkComplexFloat64) (https://simpleitk.readthedocs.io/en/master/Documentation/docs/source/fundamentalConcepts.html)
SimpleITK images are either 2D, 3D, or 4D and can have an arbitrary number of channels with a scalar or complex value in each channel.
However, I can’t go much further than to create an empty complex image. If I try to export that image as a numpy array I get the error message
RuntimeError: Images of Complex Pixel types currently are not supported.
I get similar error if I try to create a image from a complex numpy array, or if I try the slicing operators on my complex image.
Will this support be added? Otherwise I think it is a bit misleading to state that SimpleITK support complex values, when all you can do is create an empty complex image.
Thank you in advance!