[Insight-users] Bug in itkImageRegion

[Koen Van Leemput]

Hi all,=20

I think there is a bug in itkImageRegion::IsInside(const=20
ContinuousIndex<TCoordRepType,ImageDimension> &): indices outside a regio=
n=20
are actually reported to fall inside. Here is some example code reproduci=
ng=20
the error:

  typedef itk::ImageRegion<1>        RegionType;
  RegionType::IndexType    index =3D {0};
  RegionType::SizeType     size =3D {100};
  RegionType region(index, size);
  itk::ContinuousIndex<double, 1>  continuousIndex;
  continuousIndex[0] =3D 99.9;
  std::cout << region.IsInside( continuousIndex ) << std::endl;

In itkImageRegion::IsInside(const=20
ContinuousIndex<TCoordRepType,ImageDimension> &), the test "if( index[i] =
>=3D=20
bound )" should probably be replaced with "if ( index[i] > bound-1 )"

Cheers,

Koen


--=20
*************************************************************************=
****=20
Koen Van Leemput, Ph.D.                email : koen.vanleemput@hus.fi=20
Department of Radiology                  phone: +358 9 471 71331       =20
Helsinki University Central Hospital   mobile: +358 9 471 62300       =20
P.O. Box 340                                  fax:      +358 9 471 71342 =
      =20
FIN-00029 HUS                                                            =
              =20
FINLAND                                                                  =
                 =20
*************************************************************************=
****=20