From 1c36198476c52684f55d4b6e4d30241d3a4c2e68 Mon Sep 17 00:00:00 2001
From: Romulo Pinho <pinho@lyon.fnclcc.fr>
Date: Mon, 30 May 2011 14:48:16 +0200
Subject: [PATCH] debugging messages
---
registration/clitkGenericMetric.txx | 17 +++++++++++------
1 file changed, 11 insertions(+), 6 deletions(-)
diff --git a/registration/clitkGenericMetric.txx b/registration/clitkGenericMetric.txx
index 6f10e72..b5df490 100644
--- a/registration/clitkGenericMetric.txx
+++ b/registration/clitkGenericMetric.txx
@@ -256,11 +256,15 @@ GenericMetric<args_info_type,FixedImageType, MovingImageType>::GetMetricPointer(
}
- //typedef itk::ImageMaskSpatialObject<itkGetStaticConstMacro(FixedImageDimension)> ImageMaskSpatialObjectType;
- //typename ImageMaskSpatialObjectType::ConstPointer mask = dynamic_cast<const ImageMaskSpatialObjectType*>(m_FixedImageMask.GetPointer());
+ typedef itk::ImageMaskSpatialObject<itkGetStaticConstMacro(FixedImageDimension)> ImageMaskSpatialObjectType;
+ typename ImageMaskSpatialObjectType::ConstPointer mask = NULL;
+ if (m_FixedImageMask.IsNotNull())
+ mask = dynamic_cast<const ImageMaskSpatialObjectType*>(m_FixedImageMask.GetPointer());
- //typedef typename ImageMaskSpatialObjectType::RegionType ImageMaskRegionType;
- //ImageMaskRegionType mask_region = mask->GetAxisAlignedBoundingBoxRegion();
+ typedef typename ImageMaskSpatialObjectType::RegionType ImageMaskRegionType;
+ ImageMaskRegionType mask_region;
+ if (mask.IsNotNull())
+ mask_region = mask->GetAxisAlignedBoundingBoxRegion();
// Common properties
if( m_FixedImageMask.IsNotNull() )
@@ -300,7 +304,7 @@ GenericMetric<args_info_type,FixedImageType, MovingImageType>::GetMetricPointer(
// Calculate the number
const unsigned int totalNumberOfPixels = m_FixedImageRegion.GetNumberOfPixels();
- //const unsigned int totalNumberOfPixels = mask_region.GetNumberOfPixels();
+ const unsigned int totalNumberOfMaskPixels = mask_region.GetNumberOfPixels();
const unsigned int numberOfDemandedPixels = static_cast< unsigned int >( (double) totalNumberOfPixels *m_ArgsInfo.samples_arg );
// --------------------------------------------------
@@ -359,7 +363,7 @@ GenericMetric<args_info_type,FixedImageType, MovingImageType>::GetMetricPointer(
//RandomIterator randIter( m_FixedImage, mask_region );
if (m_Verbose) std::cout << "Search region " << m_FixedImageRegion << std::endl;
- //if (m_Verbose) std::cout << "Mask search region " << mask_region << std::endl;
+ if (m_Verbose) std::cout << "Mask search region " << mask_region << std::endl;
// Randomly sample the image
short att = 1;
@@ -427,6 +431,7 @@ GenericMetric<args_info_type,FixedImageType, MovingImageType>::GetMetricPointer(
if (m_Verbose) std::cout<<"A fraction of "<<m_ArgsInfo.samples_arg<<" spatial samples was requested..."<<std::endl;
double fraction=(double)numberOfValidPixels/ (double) totalNumberOfPixels;
if (m_Verbose) std::cout<<"Found "<<numberOfValidPixels <<" valid pixels for a total of "<<totalNumberOfPixels<<" (a fraction of "<<fraction<<")..."<<std::endl;
+ if (m_Verbose) std::cout<<"number of mask pixels "<<totalNumberOfMaskPixels<<std::endl;
}
--
2.49.0