[Insight-users] check is point inside a closed surface or not

[shirani kannangara]

Hi itk users

I am attempting  to check whether  a point is inside a closed surface(poly
dat mesh) or not. For instance I created a sphere  and used the codes given
under vtk examples *VTK/Examples/Cxx/PolyData/PointInsideObject. *I have to
feed data as an array as my original program (based on vhp abdomen data
set) needs to feed a  large range of data.( eg. x= 675 , y= 401,  z=450 ).
The program was built with zero errors.But it does not check inside
correctly.(is inside=1,outside=0).I spent a long time to make it correct
but failed.Can itk do this.
.Can any body help me to get it corrected or propose a suitable way to do
this task.My codes are attached herewith  for you reference.

Many thanks
Shirani
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.itk.org/pipermail/insight-users/attachments/20130426/be0f1378/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: sphere1.cxx
Type: application/octet-stream
Size: 9119 bytes
Desc: not available
URL: <http://www.itk.org/pipermail/insight-users/attachments/20130426/be0f1378/attachment.obj>