[vtkusers] Stream lines from custom data
Berk Geveci
berk.geveci at kitware.com
Fri Apr 16 20:09:08 EDT 2010
I don't see a SetDimensions on the structured grid. That is probably
part of the problem.
-berk
On Wed, Apr 14, 2010 at 11:09 AM, Gomez, Alberto
<alberto.gomez at kcl.ac.uk> wrote:
> Thanks for the quick answer,
>
> Unfortunately, it doesn't work for me: the vtkCellToPointData filter
> crashes. Would you have a sample code that I can compare to mine? That would
> be very useful :-).
>
> Anyway, I found out a tournaround that, I think, gets me closer to the
> visualisation of streamlines (but does not work yet):
> ...
> ...
> vtkSmartPointer<vtkPolyData> sGrid = vtkSmartPointer<vtkPolyData>::New();
> sGrid->GetPointData()->SetVectors(vecArr);
> sGrid->SetPoints(pts);
> vtkSmartPointer<vtkStreamLine> streamer =
> vtkSmartPointer<vtkStreamLine>::New();
> streamer->SetInput(sGrid);
>
> Then I still see black but I get
>
> Debug: In /data/workspace/vtk-5.4.2/src/Filtering/vtkAlgorithm.h, line 150
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): setting AbortExecute
> to 0
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.cxx, line 475
> vtkStreamLine (0x22fd980): Generating streamers
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 174
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning Integrator
> address 0x2300f50
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 174
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning Integrator
> address 0x2300f50
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 174
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning Integrator
> address 0x2300f50
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 174
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning Integrator
> address 0x2300f50
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 149
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> TerminalSpeed of 0
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 107
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> MaximumPropagationTime of 200
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 165
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> SavePointInterval of 0.1
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 290
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> NumberOfStreamers of 1
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 290
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> NumberOfStreamers of 1
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 174
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning Integrator
> address 0x2300f50
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 174
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning Integrator
> address 0x2300f50
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 149
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> TerminalSpeed of 0
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 107
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> MaximumPropagationTime of 200
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 165
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> SavePointInterval of 0.1
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 290
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> NumberOfStreamers of 1
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 290
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> NumberOfStreamers of 1
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 174
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning Integrator
> address 0x2300f50
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 174
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning Integrator
> address 0x2300f50
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 149
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> TerminalSpeed of 0
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 107
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> MaximumPropagationTime of 200
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 165
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> SavePointInterval of 0.1
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 290
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> NumberOfStreamers of 1
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 290
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> NumberOfStreamers of 1
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 149
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> TerminalSpeed of 0
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 107
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> MaximumPropagationTime of 200
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 165
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> SavePointInterval of 0.1
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 290
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> NumberOfStreamers of 1
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.h, line 290
> vtkStreamLine (0x22fd980): vtkStreamLine (0x22fd980): returning
> NumberOfStreamers of 1
>
> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamLine.cxx, line 178
> vtkStreamLine (0x22fd980): Created 0 points, 0 lines
>
>
> Any clue?
>
> Thanks
> Sebastian Gatzka wrote:
>>
>> I had the same problem.
>>
>> You will have to use the vtkCellDataToPointData filter first!
>> Input to the filter should be your structured grid.
>>
>> And then you can use the output from this filter as input to the stream
>> line.
>>
>> -----Ursprüngliche Nachricht-----
>> Von: vtkusers-bounces at vtk.org [mailto:vtkusers-bounces at vtk.org] Im Auftrag
>> von Gomez, Alberto
>> Gesendet: Mittwoch, 14. April 2010 16:05
>> Cc: vtkusers
>> Betreff: [vtkusers] Stream lines from custom data
>>
>> Hi all,
>>
>> I have been looking for a way to represent a vector field with stream
>> lines, but I have not found a solution so far. My case differs from the
>> example in http://www.itk.org/Wiki/VTK/Examples/StreamLines in two things:
>>
>> I would like to use vtkStreamTracer as vtkStreamLine may be deprecated
>> soon.
>> I am using custom data and not reading from a file.
>>
>> So I have my points in the variable positions[][] ( size n x 3) and my
>> vectors in the variable directions[][] (size n x 3 also).
>>
>> I have tried both StreamTracer and StreamLines but all I see is a black
>> background. When using vtkStreamLines I get the error :
>>
>> Debug: In /data/workspace/vtk-5.4.2/src/Filtering/vtkAlgorithm.h, line 150
>> vtkStreamLine (0x22a0880): vtkStreamLine (0x22a0880): setting AbortExecute
>> to 0
>>
>> Debug: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.cxx, line 475
>> vtkStreamLine (0x22a0880): Generating streamers
>>
>> ERROR: In /data/workspace/vtk-5.4.2/src/Graphics/vtkStreamer.cxx, line 488
>> vtkStreamLine (0x22a0880): No vector data defined!
>>
>> But I actually define the vectors, so I do not know what I am doing wrong.
>>
>> This is my source code (an extract of it)
>>
>>
>> vtkSmartPointer<vtkPoints> pts = vtkSmartPointer<vtkPoints>::New();
>>
>> vtkSmartPointer<vtkDoubleArray> vecArr =
>> vtkSmartPointer<vtkDoubleArray>::New();
>> vecArr->SetNumberOfComponents(3);
>>
>> for (int i = 0; i < n; i++) {
>> pts->InsertNextPoint(positions[i][0], positions[i][1],
>> positions[i][2]);
>> vecArr->InsertNextTuple3(directions[i][0] , directions[i][1] ,
>> directions[i][2] );
>> }
>>
>> // Generate the vector field
>> vtkSmartPointer<vtkStructuredGrid> sGrid =
>> vtkSmartPointer<vtkStructuredGrid>::New();
>> sGrid->GetCellData()->SetVectors(vecArr);
>> sGrid->SetPoints(pts);
>> //Generate start points
>> vtkSmartPointer<vtkPointSource> seeds =
>> vtkSmartPointer<vtkPointSource>::New();
>> seeds->SetRadius(20.0);
>> seeds->SetCenter(positions[0][0], positions[0][1],
>> positions[0][2]);
>> seeds->SetNumberOfPoints(100);
>>
>> //vtkSmartPointer<vtkStreamTracer> streamer =
>> vtkSmartPointer<vtkStreamTracer>::New();
>> vtkSmartPointer<vtkStreamLine> streamer =
>> vtkSmartPointer<vtkStreamLine>::New();
>>
>> streamer->SetInput(sGrid);
>> streamer->SetSource(seeds->GetOutput());
>> streamer->SetIntegrationDirection(0);
>> streamer->DebugOn();
>> streamer->Update();
>> vtkSmartPointer<vtkPolyDataMapper> gMapper = vtkSmartPointer<
>> vtkPolyDataMapper>::New();
>> gMapper->SetInputConnection(streamer->GetOutputPort());
>> ...
>>
>>
>> Many thanks!!!!
>>
>>
>
>
> --
>
> Alberto Gómez
>
> /Division of Imaging Sciences
> The Rayne Institute
> 4th Floor, Lambeth Wing
> St Thomas' Hospital
> London SE1 7EH /
>
> phone: +44 (0) 20 718 88364
> email: alberto.gomez at kcl.ac.uk <mailto:alberto.gomez at kcl.ac.uk>
>
> _______________________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Please keep messages on-topic and check the VTK FAQ at:
> http://www.vtk.org/Wiki/VTK_FAQ
>
> Follow this link to subscribe/unsubscribe:
> http://www.vtk.org/mailman/listinfo/vtkusers
>
More information about the vtkusers
mailing list