[Rtk-users] FDK reconstruction error

[Simon Rit]

What about ITK_USE_CUFFTW? I would keep it OFF to see if it fixes the
problem, RTK has its own implementation of this.
Simon

On Tue, Jun 4, 2024 at 8:11 AM Akshara P K <akshara at advitech.in> wrote:

> Hi Simon,
> Yes, I already tried it. Tried enabling 'ITK_USE_FFTWD' and
> 'ITK_USE_FFTWF' and made use of those libs and includes in the project.
> Still the error is coming.
> [image: image.png]
>
> Regards,
> Akshara
>
> On Tue, 4 Jun 2024 at 10:29, Simon Rit <simon.rit at creatis.insa-lyon.fr>
> wrote:
>
>> Hi,
>> I have never encountered this issue. Have you set ITK_USE_FFTWF and/or
>> ITK_USE_CUFFTW to ON in your CMake options?
>> Simon
>>
>> On Tue, Jun 4, 2024 at 6:50 AM Akshara P K <akshara at advitech.in> wrote:
>>
>>> Hi all,
>>>
>>>   We're currently engaged in CBCT reconstruction and have successfully
>>> developed a Python-based code for the task, which functions without issue.
>>> However, encountering an error during the conversion process to C++,
>>> specifically at the stage of FDK reconstruction filter implementation, has
>>> prompted our inquiry. The error message reads as follows:
>>>
>>> itk::ExceptionObject (000000EBA52FDC68)
>>> Location: "unknown"
>>> File:
>>> C:\ITK\src\Modules\Filtering\FFT\include\itkRealToHalfHermitianForwardFFTImageFilter.h
>>> Line: 85
>>> Description: ITK ERROR: Object factory failed to instantiate class
>>> itk::RealToHalfHermitianForwardFFTImageFilter<class
>>> itk::Image<float,3>,class itk::Image<class std::complex<float>,3> >
>>>
>>>   Our FDK implementation in C++ is as follows:
>>>
>>> using FDKGPUType = rtk::CudaFDKConeBeamReconstructionFilter;
>>>
>>> // Creating the FDK reconstruction filter
>>> FDKGPUType::Pointer feldkamp = FDKGPUType::New();
>>>
>>> // Setting input images for the FDK reconstruction filter
>>> feldkamp->SetInput(0, constantImageSource2->GetOutput());
>>> feldkamp->SetInput(1, ParkerFilter->GetOutput()); // This is the
>>> projection stack from rtkProjectionsReader
>>>
>>> // Setting geometry for the FDK reconstruction filter
>>> feldkamp->SetGeometry(geometry);
>>>
>>> // Setting ramp filter parameters
>>> feldkamp->GetRampFilter()->SetTruncationCorrection(0.0);
>>> feldkamp->GetRampFilter()->SetHannCutFrequency(0.0);
>>>
>>> // Executing the reconstruction filter
>>> try {
>>>     feldkamp->Update();
>>> }
>>>
>>> ChatGPT
>>>
>>> We're currently engaged in CBCT reconstruction and have successfully
>>> developed a Python-based code for the task, which functions without issue.
>>> However, encountering an error during the conversion process to C++,
>>> specifically at the stage of FDK reconstruction filter implementation, has
>>> prompted our inquiry. The error message reads as follows:
>>>
>>> vbnet
>>> Copy code
>>> itk::ExceptionObject (000000EBA52FDC68)Location: "unknown"File: C:\ITK\src\Modules\Filtering\FFT\include\itkRealToHalfHermitianForwardFFTImageFilter.hLine: 85Description: ITK ERROR: Object factory failed to instantiate class itk::RealToHalfHermitianForwardFFTImageFilter<class itk::Image<float,3>,class itk::Image<class std::complex<float>,3> >
>>>
>>> Our FDK implementation in C++ is as follows:
>>>
>>> cpp
>>> Copy code
>>> using FDKGPUType = rtk::CudaFDKConeBeamReconstructionFilter;
>>> // Creating the FDK reconstruction filter
>>> FDKGPUType::Pointer feldkamp = FDKGPUType::New();
>>> // Setting input images for the FDK reconstruction filter
>>> feldkamp->SetInput(0, constantImageSource2->GetOutput());
>>> feldkamp->SetInput(1, ParkerFilter->GetOutput()); // This is the projection stack from rtkProjectionsReader
>>> // Setting geometry for the FDK reconstruction filter
>>> feldkamp->SetGeometry(geometry);
>>> // Setting ramp filter parameters
>>> feldkamp->GetRampFilter()->SetTruncationCorrection(0.0);
>>> feldkamp->GetRampFilter()->SetHannCutFrequency(0.0);
>>> // Executing the reconstruction filtertry {
>>>     feldkamp->Update();
>>> }
>>>
>>> We're employing ITK version 5.3 for this implementation. Any insights
>>> into the cause of this error would be greatly appreciated.
>>>
>>> Regards,
>>> Akshara
>>>
>>> _______________________________________________
>>> Rtk-users mailing list
>>> rtk-users at openrtk.org
>>> https://www.creatis.insa-lyon.fr/mailman/listinfo/rtk-users
>>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.creatis.insa-lyon.fr/pipermail/rtk-users/attachments/20240604/33f91824/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 2682 bytes
Desc: not available
URL: <http://www.creatis.insa-lyon.fr/pipermail/rtk-users/attachments/20240604/33f91824/attachment-0001.png>