[Rtk-users] FDK reconstruction error

[Akshara P K]

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/5a800333/attachment.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/5a800333/attachment.png>