songsenand
/
cgal
mirror of https://github.com/CGAL/cgal
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

duplicate image before application of filters

This commit is contained in:
Jane Tournois 2021-06-15 16:01:21 +02:00
parent 16db150f65
commit c07c71275d
1 changed files with 7 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
Mesh_3/include/CGAL/Mesh_3/generate_weights_from_labeled_image.h
View File

@ -18,6 +18,7 @@
#include <CGAL/Image_3.h>
#include <itkImage.h>
#include <itkImageDuplicator.h>
#include <itkThresholdImageFilter.h>
#include <itkRescaleIntensityImageFilter.h>
#include <itkSmoothingRecursiveGaussianImageFilter.h>
@ -85,6 +86,7 @@ void generate_weights(const CGAL::Image_3& image,
boost::container::flat_set<Image_word_type> labels;
internal::convert_image_3_to_itk(image, *itk_img, labels);
using DuplicatorType = itk::ImageDuplicator<ImageType>;
using IndicatorFilter = itk::ThresholdImageFilter<ImageType>;
using RescaleFilterType = itk::RescaleIntensityImageFilter<ImageType, ImageType>;
using GaussianFilterType = itk::SmoothingRecursiveGaussianImageFilter<ImageType, ImageType>;
@ -93,9 +95,13 @@ void generate_weights(const CGAL::Image_3& image,
for (Image_word_type label : labels)
{
DuplicatorType::Pointer duplicator = DuplicatorType::New();
duplicator->SetInputImage(itk_img);
duplicator->Update();
//compute "indicator image" for "label"
IndicatorFilter::Pointer indicator = IndicatorFilter::New();
indicator->SetInput(itk_img);
indicator->SetInput(duplicator->GetOutput());
indicator->SetOutsideValue(0);
indicator->ThresholdOutside(label, label);
indicator->Update();