[vtkusers] Grid intersections?

[David.Pont at ensisjv.com]

"Andrea Gavana" <andrea.gavana at gmail.com> wrote on 14/12/2006 06:59:12:

> Hi all,
>
>    I have a 3D unstructured grid (a reservoir representation) and I
> am faced with a challenge for which I don't even know if a solution
> exists. If it exists, I don't know where to look for it. I've googled
> around quite a bit but got zero clues.
> Based on 3 user inputs (width, height, depth), I would like to
> construct some kind of parallelepiped and then I would like to move it
> (in x, y, z directions sequentially) inside my 3D unstructred grid:
> for evey move, I would like to know which cells of my 3D grid are
> intersected by this parallelepiped. Actually, the parallelepiped will
> be distorted, as the 3D grid is not orthogonal (is a corner point
> grid), so the problem looks like a 3D unstructured grid that
> intersects a smaller 3D unstructured grid, and I would like to find
> out which cells of the bigger grid are intersected by the smaller one.
> It is like a moving 3D window inside the bigger grid.
> Does anyone have a suggestion or a comment to share?
>
> Thank you for every hint.
>
> --
> Andrea.
>
> "Imagination Is The Only Weapon In The War Against Reality."
> http://xoomer.virgilio.it/infinity77/
> _______________________________________________
> This is the private VTK discussion list.
> Please keep messages on-topic. Check the FAQ at: http://www.vtk.
> org/Wiki/VTK_FAQ
> Follow this link to subscribe/unsubscribe:
> http://www.vtk.org/mailman/listinfo/vtkusers

Hi Andrea,
   You could look at vtkCutter or vtkClipDataSet using an instance of
vtkPlanes (implicit function) for your parallelepiped. If you really just
want to know which cells are intersected (output cell Ids?) you may have to
code something yourself starting with one of them.
  regards
     Dave P
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.vtk.org/pipermail/vtkusers/attachments/20061214/4b8b1497/attachment.htm>