Authors belong to:
- University of LYON http://www.universite-lyon.fr/
- - Léon Bérard cancer center http://oncora1.lyon.fnclcc.fr
+ - Léon Bérard cancer center http://www.centreleonberard.fr
- CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr
This software is distributed WITHOUT ANY WARRANTY; without even
- BSD See included LICENSE.txt file
- CeCILL-B http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
- ======================================================================-====*/
+ ===========================================================================**/
#ifndef CLITKAUTOCROPFILTER_TXX
#define CLITKAUTOCROPFILTER_TXX
// clitk
#include "clitkCommon.h"
-#include "clitkSegmentationUtils.h"
// itk
#include "itkAutoCropLabelMapFilter.h"
AutoCropFilter():itk::ImageToImageFilter<ImageType, ImageType>() {
this->SetNumberOfRequiredInputs(1);
m_BackgroundValue = 0;
+ UseBorderOff();
}
//--------------------------------------------------------------------
typedef itk::AutoCropLabelMapFilter<LabelMapType> AutoCropFilterType;
typename AutoCropFilterType::Pointer autoCropFilter = AutoCropFilterType::New();
autoCropFilter->SetInput(imageToLabelFilter->GetOutput());
- autoCropFilter->ReleaseDataFlagOff();
-
+ // autoCropFilter->ReleaseDataFlagOff();
+ if (GetUseBorder()) {
+ DD("Crop UseBorder : not correctly implemented do not use (use PadLabelMapFilter)");
+ // http://www.itk.org/Doxygen/html/classitk_1_1AutoCropLabelMapFilter.html#a54f49fdff8d9f2d2313134109d510285
+ exit(0);
+ typename ImageType::SizeType s;
+ for(uint i=0; i<ImageType::ImageDimension; i++) s[i] = 1;
+ autoCropFilter->SetCropBorder(s);
+ }
+ autoCropFilter->ReleaseDataFlagOn();
+
// Convert to LabelImage
typedef itk::LabelMapToLabelImageFilter<LabelMapType, ImageType> MapToImageFilterType;
typename MapToImageFilterType::Pointer labelToImageFilter = MapToImageFilterType::New();
m_Region = m_labeImage->GetLargestPossibleRegion();
// Sometimes the index is 9223372036854775807 ???
if (m_Region.GetIndex()[0] > 99999) {
+ std::cerr << "Warning !! BUG int clitkAutoCropFilter ?" << std::endl;
typename ImageType::IndexType index;
index.Fill(0);
m_Region.SetIndex(index);
}
+
+ // Set the region to output
output->SetLargestPossibleRegion(m_Region);
output->SetRequestedRegion(m_Region);
output->SetBufferedRegion(m_Region);
// Get input pointers
ImageConstPointer input = dynamic_cast<const ImageType*>(itk::ProcessObject::GetInput(0));
- // Extract the region
+ // Extract the region with RegionOfInterestImageFilter or ExtractImageFilter ?
+ // The second is when reducing the nb of dimension (index always zero)
+ // The first keep index.
+ // OLD : typedef itk::ExtractImageFilter<ImageType, ImageType> CropFilterType;
+ // OLD : cropFilter->SetExtractionRegion(m_Region);
- //TO CHANGE WITH EXTRACTIMAGEFILTER !
-
- // typedef itk::RegionOfInterestImageFilter<ImageType, ImageType> CropFilterType;
- typedef itk::ExtractImageFilter<ImageType, ImageType> CropFilterType;
+ typedef itk::RegionOfInterestImageFilter<ImageType, ImageType> CropFilterType;
m_labeImage->SetRequestedRegion(m_labeImage->GetLargestPossibleRegion());
typename CropFilterType::Pointer cropFilter = CropFilterType::New();
cropFilter->SetInput(m_labeImage);
-
- // cropFilter->SetRegionOfInterest(m_Region);
- cropFilter->SetExtractionRegion(m_Region);
-
- cropFilter->ReleaseDataFlagOff();
+ cropFilter->SetReleaseDataFlag(this->GetReleaseDataFlag());
+ cropFilter->SetRegionOfInterest(m_Region);
// Go !
cropFilter->Update();