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 <[email protected]> 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 > [email protected] > https://public.kitware.com/mailman/listinfo/rtk-users >
_______________________________________________ Rtk-users mailing list [email protected] https://public.kitware.com/mailman/listinfo/rtk-users
