I want to change the code to CUDA version in inline.py, but it seems to be different from the cpu version. my code:
extractor = itk.ExtractImageFilter[image_type, image_type].New(Input=reader.GetOutput()) cuda_input = cuda_type.New() cuda_input.SetPixelContainer(extractor.GetOutput().GetPixelContainer()) cuda_input.CopyInformation(extractor.GetOutput()) cuda_input.SetBufferedRegion(extractor.GetOutput().GetBufferedRegion()) cuda_input.SetRequestedRegion(extractor.GetOutput().GetRequestedRegion()) parker = rtk.ParkerShortScanImageFilter[image_type].New(Input=extractor.GetOutput(), Geometry=geometry_rec) fdk = rtk.CudaFDKConeBeamReconstructionFilter.New(Geometry=geometry_rec)
However, in fdk.GetOutput().UpdateOutputInformation(),
ITK ERROR: ExtractImageFilter(0x64d27b9d57e0) appears: The number of zero sized dimensions in the input image Extraction Region is not consistent with the dimensionality of the output image.
Expected the extraction region size ([0, 0, 0]) to contain 0 zero sized dimensions to collapse. Mistake.
I carefully compared the cpu version, the shape of each variable is the same, why the CUDA version will appear such a mistake? It is worth mentioning that I changed the code for reading and saving the projection to:
projection = itk.GetImageFromArray(projections[i], is_vector=False) projection.SetOrigin([0, 0]) projection.SetSpacing([0.198 * 3, 0.198 * 3]) projection.Allocate()