Michael,
thanks for pointing me there.
I gave it a try, but I'm actually stuck.
Config: winXp(SP2) + msvc7.0 + vtk5.0.
So I downloaded everything required from bioengineering web site
(vtkBioeng.zip), ran cmake and eventually built vtkBioeng library as a
shared library (vtkBioeng.dll). So far so good.
Then I translated most of the testCollisionDetection.tcl sample given
from tcl to c++ (see code below)
When I run it, I end up with an error located at the following line:
vtkCollisionDetectionFilter *collide =
vtkCollisionDetectionFilter::New();
Trying to debug it, an access violation is reported in the constructor
of vtkCollisionDetectionFilter, at line #4 :
this->Inputs[0] = NULL;
Apparently, from what I understand, it looks like the array Inputs[]
(inherited from vtkPolyDataSource) is not initialized (value is 0x00000000).
In case you have you already tested this sample :
- did you experience similar problem?
- if no, was it in c++ or in tcl, with which vtk version 4.x or 5.x?
Thanks for helping
Jean-Do
File testCollisionDetection.cpp
#include "vtkProperty.h"
#include "vtkPolyDataMapper.h"
#include "vtkActor.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkRenderer.h"
#include "vtkCamera.h"
#include "vtkSphereSource.h"
#include "vtkMatrix4x4.h"
#include "vtkCollisionDetectionFilter.h"
#include "vtkGlyph3D.h"
#include "vtkInteractorStyleJoystickActor.h"
#include "vtkTextActor.h"
int _tmain(int argc, _TCHAR* argv[])
{
vtkSphereSource *sphere0 = vtkSphereSource::New();
sphere0->SetPhiResolution(30);
sphere0->SetThetaResolution(30);
sphere0->SetCenter(-0.7, 0, 0);
sphere0->Update();
vtkSphereSource *sphere1 = vtkSphereSource::New();
sphere1->SetPhiResolution(30);
sphere1->SetThetaResolution(30);
sphere1->Update();
vtkMatrix4x4 *matrix0 = vtkMatrix4x4::New();
vtkMatrix4x4 *matrix1 = vtkMatrix4x4::New();
vtkCollisionDetectionFilter *collide =
vtkCollisionDetectionFilter::New();
collide->SetInput(0, sphere0->GetOutput());
collide->SetMatrix(0, matrix0);
collide->SetInput(1, sphere1->GetOutput());
collide->SetMatrix(1, matrix1);
collide->SetBoxTolerance(0.0);
collide->SetCellTolerance(0.0);
collide->SetNumberOfCellsPerBucket(2);
//collide->SetCollisionModeToFirstContact();
collide->GeneratePolydataOutputOn();
// collide->AddObserver->EndEvent(cbCollision());
vtkTextActor *txt = vtkTextActor::New();
vtkSphereSource *point = vtkSphereSource::New();
point->SetRadius(0.01);
vtkPolyData *pointdata = vtkPolyData::New();
pointdata->SetPoints(collide->GetContactPoints());
vtkGlyph3D *points = vtkGlyph3D::New();
points->SetInput(pointdata);
points->SetSource(point->GetOutput());
vtkPolyDataMapper *mapper1 = vtkPolyDataMapper::New();
mapper1->SetInput(collide->GetOutput());
vtkActor *actor1 = vtkActor::New();
actor1->SetMapper(mapper1);
(actor1->GetProperty())->BackfaceCullingOn();
actor1->SetUserMatrix(matrix0);
vtkPolyDataMapper *mapper2 = vtkPolyDataMapper::New();
mapper2->SetInput(collide->GetOutput(1));
vtkActor *actor2 = vtkActor::New();
actor2->SetMapper(mapper2);
(actor2->GetProperty())->BackfaceCullingOn();
actor2->SetUserMatrix(matrix1);
vtkPolyDataMapper *mapper3 = vtkPolyDataMapper::New();
mapper3->SetInput(points->GetOutput());
vtkActor *actor3 = vtkActor::New();
actor3->SetMapper(mapper3);
(actor3->GetProperty())->SetColor(0,0,0);
vtkRenderer *ren = vtkRenderer::New();
ren->AddActor(actor1);
ren->AddActor(actor2);
//ren AddObserver EndEvent ChangeOrigin
//ren AddActor actor3
ren->AddActor(txt);
ren->SetBackground(0.5,0.5,0.5);
vtkRenderWindow *renWin = vtkRenderWindow::New();
renWin->AddRenderer(ren);
vtkInteractorStyleJoystickActor *istyle =
vtkInteractorStyleJoystickActor::New();
vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();
iren->SetRenderWindow(renWin);
iren->SetInteractorStyle(istyle);
iren->Initialize();
renWin->Render();
iren->Start();
}
http://www.bioengineering-research.com/vtk/vtkCollisionDetectionFilter.htm
Post by Jean-Dominique BarnichonOne of the surface is not a plane, it is a polyplane (i.e. a 2D (xy)
polyline extruded along the z direction.
Jean-Do
Post by Bill LorensenIf one surface is a plane, you can use ClipPolyData and define an
implicit plane as the clipper.
Bill
Post by Jean-Dominique BarnichonDear all,
I need to find the intersection between two surfaces, each one
being defined as a polydata (actually one surface is a vertical
polyplane, and the other one is a topographical surface).
Such intersection operation is already effective in other libraries
such as CGAL (http://www.cgal.org) and GTS
(http://gts.sourceforge.net).
However, it is still not clear to me which approach i should follow
to perform such an operation with vtk.