[Rtk-users] FDK reconstruction error

[Simon Rit]

I don't know what could be the issue. Can you try the following code alone
using FFTType = itk::RealToHalfHermitianForwardFFTImageFilter<class
itk::Image<float,3>>;
typename FFTType::Pointer fft = FFTType::New();
If it doesn't work, you can ask on discourse.itk.org, it's not an RTK issue.
Simon


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

> It's already OFF in Cmake. Here is the screenshot.
> [image: image.png]
>
> Regards,
> Akshara
>
>
> On Tue, 4 Jun 2024 at 11:46, Simon Rit <simon.rit at creatis.insa-lyon.fr>
> wrote:
>
>> 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/8aed1b95/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/8aed1b95/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 21424 bytes
Desc: not available
URL: <http://www.creatis.insa-lyon.fr/pipermail/rtk-users/attachments/20240604/8aed1b95/attachment-0003.png>