if (mArgsInfo.mask_given) this->AddInputFilename(mArgsInfo.mask_arg);
}
//--------------------------------------------------------------------
if (mArgsInfo.mask_given) this->AddInputFilename(mArgsInfo.mask_arg);
}
//--------------------------------------------------------------------
double volume = 1000; //1 cc into mc
const double PI = 3.141592653589793238463;
double radius = std::pow(3*volume/(4*PI),1./3);
double volume = 1000; //1 cc into mc
const double PI = 3.141592653589793238463;
double radius = std::pow(3*volume/(4*PI),1./3);
typename ImageType::Pointer kernel = ComputeMeanFilterKernel<ImageType>(input->GetSpacing(), radius);
// Perform the convolution
typename ImageType::Pointer kernel = ComputeMeanFilterKernel<ImageType>(input->GetSpacing(), radius);
// Perform the convolution
typedef itk::ImageRegionConstIteratorWithIndex<ImageType> IteratorType;
typedef itk::ImageRegionConstIteratorWithIndex<MaskImageType> MIteratorType;
IteratorType iters(output, output->GetLargestPossibleRegion());
typedef itk::ImageRegionConstIteratorWithIndex<ImageType> IteratorType;
typedef itk::ImageRegionConstIteratorWithIndex<MaskImageType> MIteratorType;
IteratorType iters(output, output->GetLargestPossibleRegion());
//--------------------------------------------------------------------
template<class args_info_type>
template<class ImageType>
//--------------------------------------------------------------------
template<class args_info_type>
template<class ImageType>
SUVPeakGenericFilter<args_info_type>::ComputeMeanFilterKernel(const typename ImageType::SpacingType & spacing, double radius)
{
// Some kind of cache to speed up a bit
SUVPeakGenericFilter<args_info_type>::ComputeMeanFilterKernel(const typename ImageType::SpacingType & spacing, double radius)
{
// Some kind of cache to speed up a bit