[Rtk-users] FDK reconstruction error

[Akshara P K]

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/fecce95f/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/fecce95f/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/fecce95f/attachment-0003.png>