[Insight-developers] ImageMaskSpatialObject IsInside and ValueAt methods give inconsistent results

[Andras Lasso]

Hi all,

 

We worked on a single-slice-to-volume registration algorithm at the NA-MIC
week and found that the sampling often fails if we define thin
(single-slice) image mask. We identified the root cause of the problem
(ImageSpatialObject::ValueAt uses
m_Interpolator->EvaluateAtContinuousIndex() for continuous to integer index
conversion, while ImageMaskSpatialObject::IsInside uses static_cast<int>)
and we proposed a fix. It affects all registration algorithms that use a
mask. The impact is significant when many mask voxels are on the boundary
compared to the total number of mask voxels.

 

See more details and proposed fix at:
http://public.kitware.com/Bug/view.php?id=11699

 

Could someone have a look at the problem and commit a fix? (I don't have
write access to the ITK repository)

 

Thanks.

Andras

 

________________________________

Andras Lasso

Senior Medical Image Computing & Applications Engineer

Laboratory for Percutanous Surgery

School of Computing

Queen's University

Kingston, ON, CANADA

Email:  <mailto:lasso at cs.queensu.ca> lasso at cs.queensu.ca

Web:  <http://perk.cs.queensu.ca/> http://perk.cs.queensu.ca

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.itk.org/mailman/private/insight-developers/attachments/20110114/d7396ad8/attachment.htm>