SCI Seg3D Mailing List

Text archives Help


[Seg3D] Re: Re: Re: usage of the Intensity Correction filter


Chronological Thread 
  • From: Michael Callahan <callahan@sci.utah.edu>
  • To: seg3d@sci.utah.edu, ramon.casero@comlab.ox.ac.uk
  • Subject: [Seg3D] Re: Re: Re: usage of the Intensity Correction filter
  • Date: Tue, 05 May 2009 19:18:52 -0600

Here is the patch to fix this (for the trunk, in the Seg3D
subdirectory):

Set the itk::DiscreteGaussianImageFilter::SetUserImageSpacingOff to
reflect that we are using this filter in image space.  Otherwise for
volumes with spacings not close to 1.0 the Discrete Gausian Image Filter
is broken, Inhomogeneity Correction crashes, and the Speed functions are
computed improperly.

  Michael

On Mon, 2009-04-06 at 15:08 +0100, Ramón Casero Cañas wrote:
> Michael Callahan wrote:
> > If you want to send me your data file I can take a look at it when I get
> > a chance.
> 
> Hi Michael,
> 
> Thanks, and sorry for the late reply. I had to check and get 
> authorisation to share the data. The data is confidential for the 
> moment, so I cannot post it in the list, but I'll get in touch directly 
> with you and see whether we can arrange something.
> 
> Cheers,
> 
> .::r
> 
Index: InhomogeneityCorrectionFilter.cc
===================================================================
--- InhomogeneityCorrectionFilter.cc	(revision 44233)
+++ InhomogeneityCorrectionFilter.cc	(working copy)
@@ -372,6 +372,7 @@
   dgfilter->SetInput(itkdgimg);
   dgfilter->SetVariance(1.0);
   dgfilter->SetMaximumKernelWidth(8);
+  dgfilter->SetUseImageSpacingOff();
   dgfilter->Update();
   SCIRun::ITKDatatypeHandle output1_img = new SCIRun::ITKDatatype();
   output1_img->data_ = dgfilter->GetOutput();
Index: ITKDiscreteGaussianImageFilterTool.cc
===================================================================
--- ITKDiscreteGaussianImageFilterTool.cc	(revision 44233)
+++ ITKDiscreteGaussianImageFilterTool.cc	(working copy)
@@ -82,6 +82,7 @@
   const int max_kernel_width = Max_Error_;
 
   filter_->SetVariance(variance);
+  filter_->SetUseImageSpacingOff();
   if (max_kernel_width > 0)
   {
     filter_->SetMaximumKernelWidth(max_kernel_width);
Index: SpeedToPathTool.cc
===================================================================
--- SpeedToPathTool.cc	(revision 44233)
+++ SpeedToPathTool.cc	(working copy)
@@ -377,6 +377,7 @@
   const int max_kernel_width = 40;
 
   blur_filter->SetVariance(variance);
+  blur_filter->SetUseImageSpacingOff();
   if (max_kernel_width > 0)
   {
     blur_filter->SetMaximumKernelWidth(max_kernel_width);



Archive powered by MHonArc 2.6.16.

Top of page