[Rtk-users] ConjugateGradientConeBeamReconstructionFilter

[Simon Rit]

Hi,
The conjugate gradient algorithm solves weighted least squares. If you want
to use unweighted least squares, simply pass an image filled with ones
which has the same size as the projections.
Simon

On Sat, Jul 13, 2019 at 8:01 AM Ali ghorbanian <e.a.gh72 at gmail.com> wrote:

> Hi,
> I read the projections of a fan Beam CT by below syntax:
>
> using ReaderType = itk::ImageSeriesReader< OutputImageType >;
> ReaderType::Pointer reader = ReaderType::New();
> const unsigned int first = std::stoi("10000");
> const unsigned int last = std::stoi("10403");
> using NameGeneratorType = itk::NumericSeriesFileNames;
> NameGeneratorType::Pointer nameGenerator = NameGeneratorType::New();
> nameGenerator->SetSeriesFormat("F:/Data Aqu/98.04.17/Tiff-1/%04d.tiff");
> nameGenerator->SetStartIndex(first);
> nameGenerator->SetEndIndex(last);
> nameGenerator->SetIncrementIndex(1);
> reader->SetImageIO(itk::TIFFImageIO::New());
> reader->SetFileNames(nameGenerator->GetFileNames());
>
> How <https://fastdic.com/word/How> to <https://fastdic.com/word/to> set
> <https://fastdic.com/word/set> uniformWeightsSource
> <https://fastdic.com/word/Ali>'s parameter
> <https://fastdic.com/word/parameter> in <https://fastdic.com/word/in>
> ConjugateGradientConeBeamReconstructionFilter
> <https://fastdic.com/word/Reza>?
>
>
> typedef rtk::ConjugateGradientConeBeamReconstructionFilter<
> OutputImageType > ConjugateGradientType;
> ConjugateGradientType::Pointer conjugategradient =
> ConjugateGradientType::New();
> conjugategradient->SetInput(tomographySource->GetOutput());
> conjugategradient->SetInput(1, reader->GetOutput());
> conjugategradient->SetInput(2, *uniformWeightsSource*->GetOutput());
> conjugategradient->SetGeometry(geometry);
> conjugategradient->SetNumberOfIterations(1);
> conjugategradient->SetDisableDisplacedDetectorFilter(true);
> conjugategradient->SetGamma(1);
>
> Best regards
> ذث
> _______________________________________________
> Rtk-users mailing list
> Rtk-users at public.kitware.com
> https://public.kitware.com/mailman/listinfo/rtk-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://public.kitware.com/pipermail/rtk-users/attachments/20190715/11ef7d4b/attachment.html>