[Insight-users] exception in itk::BSplineScatteredDataPointSetToImageFilter ?

[Nicholas Tustison]

That sounds reasonable.  To be comprehensive, though, 
you might also want to mention the origin since the points 
could also be outside the left boundar(y/ies).
What that error is telling you is that the "image" domain,
which defines the parametric domain, does not contain
all of your points.

Instead of "It may have to do with the size and spacing 
values set in this filter." perhaps something like "It is required
that the parametric domain encompass all the points.  The
parametric domain can be modified and enlarged by changing 
the origin, size, and spacing."

Nick

  


On Jun 6, 2012, at 10:50 AM, Arnaud Gelas wrote:

> Hi Nick,
> 
> What about
> 
>     << "The reparameterized point component " << p[i] " is outside the corresponding parametric domain of [0, " << totalNumberOfSpans << "]. It may have to do with size and spacing values set in this filter."
> 
> Or something similar (just to provide an indication to users)?
> 
> Arnaud
> 
> On 06/05/2012 03:23 PM, Nicholas Tustison wrote:
>> 
>> That's good.
>> 
>> What would you suggest?  
>> 
>> 
>> On Jun 5, 2012, at 7:49 AM, Arnaud Gelas wrote:
>> 
>>> By slightly increasing the spacing, or by adding 1 to the image size, it works!
>>> 
>>> Would it be possible either to update the exception message, either add a comment here to know what could be the reasons and possible fixes?
>>> Would it be possible to detect such an error at an earlier stage?
>>> 
>>> Thanks,
>>> Arnaud
>>> 
>>> On 06/05/2012 01:34 PM, Nicholas Tustison wrote:
>>>> 
>>>> Yeah, so it's a floating point thing where your last point is ending
>>>> up on the boundary.  Try making the spacing slightly larger.
>>>> 
>>>> 
>>>> On Jun 5, 2012, at 3:37 AM, Arnaud Gelas wrote:
>>>> 
>>>>> Hi Nick,
>>>>> 
>>>>> Here is the print out:
>>>>> 
>>>>> terminate called after throwing an instance of 'itk::ExceptionObject' 
>>>>> what(): /home/ajg23/DEVEL/ITK/Modules/Core/Common/src/itkMultiThreader.cxx:359: 
>>>>> itk::ERROR: MultiThreader(0xaacf610): Exception occurred during SingleMethodExecute 
>>>>> /home/ajg23/DEVEL/ITK/Modules/Filtering/ImageGrid/include/itkBSplineScatteredDataPointSetToImageFilter.hxx:625: 
>>>>> itk::ERROR: PointSetToImageFilter(0xaab0d40): The reparameterized point component 2.00002 is outside the corresponding parametric domain of [0, 2].
>>>>> 
>>>>> It seems to work with 1 or 2 level, but start throwing an exception from level 3 onward.
>>>>> 
>>>>> Arnaud
>>>>> 
>>>>> On Mon, Jun 4, 2012 at 5:24 PM, Nicholas Tustison <ntustison at gmail.com> wrote:
>>>>> THat looks right.  What is the print out you get
>>>>> when the exception is thrown?
>>>>> 
>>>>> 
>>>>> On Jun 4, 2012, at 11:02 AM, Arnaud Gelas wrote:
>>>>> 
>>>>> > Nick,
>>>>> >
>>>>> > Here is the code (see below), I have been using (adapted from the test)
>>>>> >
>>>>> > Any idea?
>>>>> >
>>>>> > Thanks,
>>>>> > Arnaud
>>>>> >
>>>>> >      const unsigned int ParametricDimension = 1;
>>>>> >      typedef itk::PointSet< Vector3DType, ParametricDimension >  PointSetType;
>>>>> >
>>>>> >      int length = m_Lateral - m_AX;
>>>>> >
>>>>> >      CoordinateType dt =
>>>>> >          static_cast< CoordinateType >( 1. ) / static_cast< CoordinateType >( length );
>>>>> >
>>>>> >      CoordinateType t = 0.;
>>>>> >
>>>>> >      typename PointSetType::Pointer pointSet = PointSetType::New();
>>>>> >
>>>>> >      IndexType idx;
>>>>> >      idx.Fill( 0 );
>>>>> >
>>>>> >      int i = 0, k = 0;
>>>>> >      for( int y = m_AX; y <= m_Lateral; y++ )
>>>>> >        {
>>>>> >        if( ( y == ax ) || ( y >= m_Bottom ) )
>>>>> >          {
>>>>> >          typename PointSetType::PointType pt;
>>>>> >          pt[0] = t;
>>>>> >
>>>>> >          pointSet->SetPoint( i, pt );
>>>>> >
>>>>> >          Vector3DType v;
>>>>> >          v[0] = static_cast< typename PointSetType::PointType::CoordRepType >( contour[y] );
>>>>> >          v[1] = static_cast< typename PointSetType::PointType::CoordRepType >( y );
>>>>> >          v[2] = 0.;
>>>>> >
>>>>> >          pointSet->SetPointData( i, v );
>>>>> >
>>>>> >          i++;
>>>>> >          }
>>>>> >
>>>>> >        t += dt;
>>>>> >        }
>>>>> >
>>>>> >      typedef itk::Image< Vector3DType, ParametricDimension > ParametricImageType;
>>>>> >
>>>>> >      typename ParametricImageType::SpacingType spacing;
>>>>> >      spacing.Fill( dt );
>>>>> >
>>>>> >      typename ParametricImageType::SizeType size;
>>>>> >      size.Fill( length + 1 );
>>>>> >
>>>>> >      typename ParametricImageType::PointType origin;
>>>>> >      origin.Fill( 0. );
>>>>> >
>>>>> >      typedef itk::BSplineScatteredDataPointSetToImageFilter< PointSetType, ParametricImageType > FilterType;
>>>>> >      typename FilterType::Pointer filter = FilterType::New();
>>>>> >      filter->SetInput( pointSet );
>>>>> >      filter->SetSize( size );
>>>>> >      filter->SetOrigin( origin );
>>>>> >      filter->SetSpacing( spacing );
>>>>> >      filter->SetSplineOrder( 3 );
>>>>> >      filter->SetNumberOfLevels( 1 );
>>>>> >
>>>>> >      typename FilterType::ArrayType ncps;
>>>>> >      ncps.Fill( 4 );
>>>>> >      filter->SetNumberOfControlPoints( ncps );
>>>>> >
>>>>> >      filter->SetNumberOfThreads( 1 );
>>>>> >      filter->SetGenerateOutputImage( true );
>>>>> >      filter->Update();
>>>>> >
>>>>> > On 06/04/2012 04:55 PM, Nicholas Tustison wrote:
>>>>> >> Hi Arnaud,
>>>>> >>
>>>>> >> Yes, so that means that one of your data points is outside the
>>>>> >> user-defined parametric domain.  Even if it is right on the edge,
>>>>> >> the parametric domain is open-ended on the right side which
>>>>> >> will cause this exception to be thrown.  How are you setting up
>>>>> >> your parametric domain?  In other words, what does your call
>>>>> >> look like when you're setting up the filter?
>>>>> >>
>>>>> >> Nick
>>>>> >>
>>>>> >>
>>>>> >> On Jun 4, 2012, at 10:45 AM, Arnaud Gelas wrote:
>>>>> >>
>>>>> >>> Hi Nick,
>>>>> >>>
>>>>> >>> When using itk::BSplineScatteredDataPointSetToImageFilter, I get the following exception
>>>>> >>>
>>>>> >>> itkExceptionMacro( "The reparameterized point component "<<  p[i]
>>>>> >>> <<  " is outside the corresponding parametric domain of [0,"
>>>>> >>> <<  totalNumberOfSpans<<  "]." );
>>>>> >>>
>>>>> >>>
>>>>> >>> I could not figure out if I made a mistake when setting my parameters? or anything else?
>>>>> >>>
>>>>> >>> Any idea? or suggestion?
>>>>> >>>
>>>>> >>> Thanks,
>>>>> >>> Arnaud
>>>>> >
>>>>> 
>>>>> 
>>>> 
>>> 
>> 
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.itk.org/pipermail/insight-users/attachments/20120606/b816df01/attachment.htm>