<!doctype html public "-//w3c//dtd html 4.0 transitional//en">
<html>
Hi, there,
<p>I am implementing moving sphere on the streamline, which need to access
the data in the vtkStreamArray. Unfortunately, vtk4.0 seems does not support
this. Below is a testing program to output the number of all points on
the streamline onto terminal screen, it is very simple, any gurus want
to take a look at it and help me to find out what's wrong ?????
<br>--------------------------------------------------------------------------------
<p>#include "vtkUnstructuredGridReader.h"
<br>#include "vtkStreamLine.h"
<br>#include "vtkRungeKutta4.h"
<br>#include "vtkTubeFilter.h"
<br>#include "vtkGeometryFilter.h"
<br>#include "vtkPolyDataMapper.h"
<br>#include "vtkActor.h"
<br>#include "vtkLookupTable.h"
<br>#include "vtkPolyDataMapper.h"
<br>#include "vtkRenderWindow.h"
<br>#include "vtkRenderWindowInteractor.h"
<br>#include <iostream.h>
<br>
<p>class strm : public vtkStreamLine
<br>{
<br> public:
<br> void init(vtkUnstructuredGrid* in);
<br> vtkPolyData* getoutput(){return this->GetOutput();}
<br>};
<p>void strm::init(vtkUnstructuredGrid* in)
<br>{
<br> vtkRungeKutta4* integ = vtkRungeKutta4::New();
<p> SetInput( in );
<br> SetStartPosition( 1.4f, 2.0f, 1.2f );
<br> SetMaximumPropagationTime( 5.67517 );
<br> SetIntegrationStepLength( 0.0348257);
<br> SetStepLength( 0.0348257 );
<br> SetIntegrationDirectionToIntegrateBothDirections( );
<br> SetNumberOfThreads( 1 );
<br> SetIntegrator( integ );
<br> Update();
<p> cout<<this->GetStreamers()->GetNumberOfPoints()<<endl;
//output number of points to terminal scrren
<br>}
<p>int main( int argc, char *argv[] )
<br>{
<br> vtkUnstructuredGridReader* reader = vtkUnstructuredGridReader::New(
);
<br> vtkTubeFilter* tubeFilter = vtkTubeFilter::New();
<br> vtkGeometryFilter* filter = vtkGeometryFilter::New( );
<br> vtkPolyDataMapper* mapper = vtkPolyDataMapper::New( );
<br> vtkActor* actor = vtkActor::New( );
<p> reader->SetFileName("flowdata.vtk"); //a dataset with the type
of "vtkUnstructuredGrid"
<br> reader->Update();
<p> strm* s = new strm;
<br> s->init(reader->GetOutput());
<p>//the visualization routines are not important for this problem
<br> vtkLookupTable* lut = vtkLookupTable::New();
<br> lut->SetHueRange( 0.66667, 0.0);
<br> lut->SetNumberOfColors(64);
<p> tubeFilter->SetInput( s->getoutput( ) );
<br> tubeFilter->SetRadius( 0.0174128 );
<p> filter->SetInput( tubeFilter->GetOutput( ) );
<br> mapper->SetInput( filter->GetOutput( ) );
<br> mapper->SetColorModeToMapScalars( );
<br> mapper->SetScalarRange( 0.0,200.0 );
<br> mapper->SetLookupTable( lut );
<br> actor->SetMapper( mapper );
<br> actor->GetProperty( )->SetSpecularPower( 20.0f );
<p> vtkCamera *camera = vtkCamera::New();
<br> camera->SetPosition(1,1,1);
<br> camera->SetFocalPoint(0,0,0);
<p> // Create the rendering objects.
<br> vtkRenderer* ren = vtkRenderer::New();
<br> ren->AddActor(actor);
<br> vtkRenderWindow* renWin = vtkRenderWindow::New();
<br> renWin->AddRenderer(ren);
<br> vtkRenderWindowInteractor* iren = vtkRenderWindowInteractor::New();
<br> iren->SetRenderWindow(renWin);
<br> ren->SetActiveCamera(camera);
<br> ren->ResetCamera();
<br> ren->SetBackground(1,1,1);
<br> renWin->Render();
<br> iren->Start();
<p> reader->Delete();
<br> filter->Delete();
<br> lut->Delete();
<br> delete s;
<br> tubeFilter->Delete();
<br> mapper->Delete();
<br> actor->Delete();
<br> ren->Delete();
<br> renWin->Delete();
<br> iren->Delete();
<br> return 0;
<br>}
<br>-----------------------------------------------------------------------------------
<p>I don't know why in vtk4.0 the method "vtkStreamer::GetStreamers()"
is changed to protected(in vtk3.2 it is public) which make a little more
trouble to write the program......
<p>I can see a streamline in the rendering window when running this program.
If compile it using vtk3.2, the terminal output is 72, which is very reasonable,
if compile the same program using vtk4.0, the output is 1 !!!!!
<p>Is this indicates a bug of vtk4.0 ?? I don't have any more ideas......I
am working on SGI's platform. Anyone can give me more informations ?? Thank
a lot!!!
<pre>--
Song Li
VRAC(Virtual Reality Application Center),Iowa State University
Email: lisong@vrac.iastate.edu
Homepage: <A HREF="http://www.vrac.iastate.edu/~lisong">http://www.vrac.iastate.edu/~lisong</A></pre>
</html>