CFD Online Logo CFD Online URL
Home > Forums > ParaView

Extract Paraview data (vtkUnstructured Grid) using pvpython

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   October 28, 2013, 14:02
Default Extract Paraview data (vtkUnstructured Grid) using pvpython
New Member
Join Date: Oct 2013
Posts: 1
Rep Power: 0
lu_siyah is on a distinguished road
I am working with large FE simulation data and I wish to extract information on only coordinates (at different timesteps, precisely a 3xN List of coordinates for every timestep considered or alternatively the displacements, which I know to be stored in a Point Array) and the basic connectivity of the grid into python arrays.
For now, I am working with Paraview in the pvpython-shell and vtk and struggling with finding and arranging the information. So far I used this code to get the coordinates in the beginning (data considered originates from LS-Dyna, hence 3dplot):
import numpy as np 
from paraview import simple  
reader.PartArrays = ['Part1 (Matl2000001)'] 
# reader: <paraview.servermanager.LSDynaReader object at 0x376b5d0>  
data = simple.servermanager.Fetch(reader) 
# data: (vtkMultiBlockDataSet) 
grid = data.GetBlock(0) 
# grid: (vtkUnstructuredGrid) 
npoints = grid.GetNumberOfPoints() 
coord = np.zeros((npoints,3)) 
for i in range(npoints):
  coord[i,:] = np.array(grid.GetPoint(i))
Any ideas on how to proceed from there? That is: Get a matrix of connectivities, and get acess to the data saved in the Point Arrays OR the coordinates at different time steps??
Any help will be appreciated as I am fairly new to this sort of programming.
lu_siyah is offline   Reply With Quote

Old   November 3, 2013, 07:10
Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,964
Blog Entries: 34
Rep Power: 92
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings lu_siyah and welcome to the forum!

Sorry, I'm a bit confused with your description. Is what you're trying to do similar to what's described on this thread: Mapping in ParaView

Because if I understand you correctly, you want to resample/probe the simulation volume you have from the FE simulation, with an array of point/probe locations.
The part that I don't understand is: why do you need a connectivity matrix?

Best regards,
wyldckat is online now   Reply With Quote


pvpython, python, vtk

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Animating unsteady data from Fluent in Paraview.. frodooon ParaView 10 May 9, 2014 10:39
Extracting ParaView Data into Python Arrays Jeffzda ParaView 7 December 6, 2012 09:10
structured and unstructured grids user Main CFD Forum 6 November 25, 2010 02:14
PARAVIEW Grid visualization simone Marras ParaView 2 January 23, 2007 09:05
Grid Independent Solution Chuck Leakeas Main CFD Forum 2 May 26, 2000 11:18

All times are GMT -4. The time now is 09:59.