From: dsarrut Date: Thu, 3 Mar 2011 10:13:16 +0000 (+0000) Subject: replace old resize with new CropLikeImage X-Git-Tag: v1.2.0~216 X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=commitdiff_plain;h=9e5e119f00eb9051aa56e480aba3e86d053cdb80;p=clitk.git replace old resize with new CropLikeImage --- diff --git a/itk/clitkAddRelativePositionConstraintToLabelImageFilter.h b/itk/clitkAddRelativePositionConstraintToLabelImageFilter.h index c365571..c801dc0 100644 --- a/itk/clitkAddRelativePositionConstraintToLabelImageFilter.h +++ b/itk/clitkAddRelativePositionConstraintToLabelImageFilter.h @@ -78,7 +78,7 @@ namespace clitk { typedef itk::Image FloatImageType; /** Orientation types */ - typedef enum { RightTo = 0, LeftTo = 1, + typedef enum { AtRightTo = 0, AtLeftTo = 1, AntTo = 2, PostTo = 3, InfTo = 4, SupTo = 5, Angle = 6 } OrientationTypeEnumeration; @@ -96,9 +96,9 @@ namespace clitk { std::string GetOrientationTypeString(int i) { return m_OrientationTypeString[i]; } std::vector & GetOrientationTypeString() { return m_OrientationTypeString; } - itkGetConstMacro(ResampleBeforeRelativePositionFilter, bool); - itkSetMacro(ResampleBeforeRelativePositionFilter, bool); - itkBooleanMacro(ResampleBeforeRelativePositionFilter); + itkGetConstMacro(IntermediateSpacingFlag, bool); + itkSetMacro(IntermediateSpacingFlag, bool); + itkBooleanMacro(IntermediateSpacingFlag); itkGetConstMacro(IntermediateSpacing, double); itkSetMacro(IntermediateSpacing, double); @@ -140,7 +140,7 @@ namespace clitk { PixelType m_ObjectBackgroundValue; std::vector m_Angle1; std::vector m_Angle2; - bool m_ResampleBeforeRelativePositionFilter; + bool m_IntermediateSpacingFlag; bool m_AutoCropFlag; bool m_InverseOrientationFlag; bool m_RemoveObjectFlag; @@ -150,7 +150,7 @@ namespace clitk { virtual void GenerateInputRequestedRegion(); virtual void GenerateData(); - typedef itk::PasteImageFilter PadFilterType; + typedef itk::PasteImageFilter PasteFilterType; typename ImageType::Pointer working_image; typename ImageType::Pointer object_resampled; typename FloatImageType::Pointer relPos; diff --git a/itk/clitkAddRelativePositionConstraintToLabelImageFilter.txx b/itk/clitkAddRelativePositionConstraintToLabelImageFilter.txx index 8e12c22..f56c4f4 100644 --- a/itk/clitkAddRelativePositionConstraintToLabelImageFilter.txx +++ b/itk/clitkAddRelativePositionConstraintToLabelImageFilter.txx @@ -49,7 +49,7 @@ AddRelativePositionConstraintToLabelImageFilter(): SetBackgroundValue(0); SetObjectBackgroundValue(0); ClearOrientationType(); - ResampleBeforeRelativePositionFilterOn(); + IntermediateSpacingFlagOn(); SetIntermediateSpacing(10); AutoCropFlagOn(); InverseOrientationFlagOff(); @@ -116,8 +116,8 @@ AddOrientationTypeString(std::string t) { m_OrientationTypeString.push_back(t); switch (t[0]) { - case 'L' : AddOrientationType(LeftTo); break; - case 'R' : AddOrientationType(RightTo);break; + case 'L' : AddOrientationType(AtLeftTo); break; + case 'R' : AddOrientationType(AtRightTo);break; case 'A' : AddOrientationType(AntTo);break; case 'P' : AddOrientationType(PostTo);break; case 'S' : AddOrientationType(SupTo);break; @@ -179,11 +179,11 @@ AddOrientationType(OrientationTypeEnumeration orientation) { m_OrientationType.push_back(orientation); switch (orientation) { - case LeftTo: + case AtRightTo: m_Angle1.push_back(clitk::deg2rad(0)); m_Angle2.push_back(clitk::deg2rad(0)); break; - case RightTo: + case AtLeftTo: m_Angle1.push_back(clitk::deg2rad(180)); m_Angle2.push_back(clitk::deg2rad(0)); break; @@ -230,83 +230,86 @@ GenerateData() // Get input pointer input = dynamic_cast(itk::ProcessObject::GetInput(0)); object = dynamic_cast(itk::ProcessObject::GetInput(1)); - - //-------------------------------------------------------------------- - //-------------------------------------------------------------------- static const unsigned int dim = ImageType::ImageDimension; - StartNewStep("Initial resample"); - // Step 1 : resample - if (m_ResampleBeforeRelativePositionFilter) { - typedef clitk::ResampleImageWithOptionsFilter ResampleFilterType; - typename ResampleFilterType::Pointer resampleFilter = ResampleFilterType::New(); - resampleFilter->SetInput(object); - resampleFilter->SetOutputIsoSpacing(m_IntermediateSpacing); - resampleFilter->SetGaussianFilteringEnabled(false); - // resampleFilter->SetVerboseOptions(true); - resampleFilter->Update(); - working_image = resampleFilter->GetOutput(); - } - else { - working_image = object; - } - StopCurrentStep(working_image); // Step 2: object pad to input image -> we want to compute the // relative position for each point belonging to the input image // domain, so we have to extend (pad) the object image to fit the // domain size + working_image = object; if (!clitk::HaveSameSizeAndSpacing(input, working_image)) { - StartNewStep("Pad object to image size"); - typename ImageType::Pointer output = ImageType::New(); - SizeType size; - for(unsigned int i=0; iGetLargestPossibleRegion().GetSize()[i]* - input->GetSpacing()[i])/(double)working_image->GetSpacing()[i]); - } + StartNewStep("Pad (resize) object to input size"); + + if (0) { // OLD VERSION (TO REMOVE) + StartNewStep("Pad object to image size"); + typename ImageType::Pointer output = ImageType::New(); + SizeType size; + for(unsigned int i=0; iGetLargestPossibleRegion().GetSize()[i]* + input->GetSpacing()[i])/(double)working_image->GetSpacing()[i]); + } - // The index of the input is not necessarily zero, so we have to - // take it into account (not done) - RegionType region; - IndexType index = input->GetLargestPossibleRegion().GetIndex(); - region.SetSize(size); - for(unsigned int i=0; iGetLargestPossibleRegion().GetIndex(); + region.SetSize(size); + for(unsigned int i=0; iSetLargestPossibleRegion(region); + output->SetRegions(region); + output->SetSpacing(working_image->GetSpacing()); + PointType origin = input->GetOrigin(); + for(unsigned int i=0; iGetSpacing()[i] + input->GetOrigin()[i]; + } + output->SetOrigin(origin); + // output->SetOrigin(input->GetOrigin()); + + output->Allocate(); + output->FillBuffer(m_BackgroundValue); + typename PasteFilterType::Pointer padFilter = PasteFilterType::New(); + // typename PasteFilterType::InputImageIndexType index; + for(unsigned int i=0; iGetSpacing()[i]/(double)working_image->GetSpacing()[i] + + lrint((working_image->GetOrigin()[i] - input->GetOrigin()[i])/working_image->GetSpacing()[i]); + } + padFilter->SetSourceImage(working_image); + padFilter->SetDestinationImage(output); + padFilter->SetDestinationIndex(index); + padFilter->SetSourceRegion(working_image->GetLargestPossibleRegion()); + padFilter->Update(); + working_image = padFilter->GetOutput(); } - // output->SetLargestPossibleRegion(region); - output->SetRegions(region); - output->SetSpacing(working_image->GetSpacing()); - PointType origin = input->GetOrigin(); - for(unsigned int i=0; iGetSpacing()[i] + input->GetOrigin()[i]; - } - output->SetOrigin(origin); - // output->SetOrigin(input->GetOrigin()); - - output->Allocate(); - output->FillBuffer(m_BackgroundValue); - typename PadFilterType::Pointer padFilter = PadFilterType::New(); - // typename PadFilterType::InputImageIndexType index; - for(unsigned int i=0; iGetSpacing()[i]/(double)working_image->GetSpacing()[i] - + lrint((working_image->GetOrigin()[i] - input->GetOrigin()[i])/working_image->GetSpacing()[i]); - } - padFilter->SetSourceImage(working_image); - padFilter->SetDestinationImage(output); - padFilter->SetDestinationIndex(index); - padFilter->SetSourceRegion(working_image->GetLargestPossibleRegion()); - padFilter->Update(); - working_image = padFilter->GetOutput(); + + // Resize object like input + working_image = clitk::ResizeImageLike(working_image, input, GetBackgroundValue()); StopCurrentStep(working_image); } - else { - // DD("[debug] RelPos : same size and spacing : no padding"); + + //-------------------------------------------------------------------- + //-------------------------------------------------------------------- + // Step 1 : resample + if (m_IntermediateSpacingFlag) { + StartNewStep("Resample object to intermediate spacing"); + typedef clitk::ResampleImageWithOptionsFilter ResampleFilterType; + typename ResampleFilterType::Pointer resampleFilter = ResampleFilterType::New(); + resampleFilter->SetInput(working_image); + resampleFilter->SetDefaultPixelValue(0); + resampleFilter->SetOutputIsoSpacing(m_IntermediateSpacing); + resampleFilter->SetGaussianFilteringEnabled(false); + // resampleFilter->SetVerboseOptions(true); + resampleFilter->Update(); + working_image = resampleFilter->GetOutput(); + StopCurrentStep(working_image); } + // Keep object image (with resampline and pad) object_resampled = working_image; - // StopCurrentStep(working_image); // Step 3: compute rel pos in object StartNewStep("Relative Position Map"); @@ -394,8 +397,8 @@ GenerateData() //-------------------------------------------------------------------- //-------------------------------------------------------------------- // Step 5: resample to initial spacing - if (m_ResampleBeforeRelativePositionFilter) { - StartNewStep("Resample to get the same sampling than input"); + if (m_IntermediateSpacingFlag) { + StartNewStep("Resample to come back to the same sampling than input"); typedef clitk::ResampleImageWithOptionsFilter ResampleFilterType; typename ResampleFilterType::Pointer resampleFilter = ResampleFilterType::New(); resampleFilter->SetDefaultPixelValue(m_BackgroundValue); @@ -419,7 +422,7 @@ GenerateData() temp->SetRegions(input->GetLargestPossibleRegion()); // Do not forget !! temp->Allocate(); temp->FillBuffer(m_BackgroundValue); - typename PadFilterType::Pointer padFilter2 = PadFilterType::New(); + typename PasteFilterType::Pointer padFilter2 = PasteFilterType::New(); padFilter2->SetSourceImage(working_image); padFilter2->SetDestinationImage(temp); // DD(input->GetLargestPossibleRegion().GetIndex());