[vtkusers] [Paraview] Combining vtk and paraview python scripts
Jeff Becker
jeffrey.c.becker at nasa.gov
Tue Jul 21 12:14:06 EDT 2015
On 07/20/2015 05:43 PM, Berk Geveci wrote:
> You probably need to update vtkImageExport and shallow copy its output
> to the programmable source's output. You also need to add something to
> the RequestInformation setting the whole extent. You need something like:
>
> from paraview import util
> util.SetOutputWholeExtent(self, [0,nx-1,0,ny-1,0,nz-1])
>
> By the way, there are better ways of getting numpy arrays into VTK
> than the image import stuff. See some of the earlier blogs here:
Right. I found the RequestInformation info above from searching the web.
Then I was able to follow the instructions in the ParaView Guide section
13.2.4 Reading binary 2D image. I basically replaced the image import
stuff with code analogous to the output assignment,
Thanks for your help.
-jeff
>
> http://www.kitware.com/blog/home/user/53/2
>
> Best,
> -berk
>
> On Mon, Jul 20, 2015 at 2:08 PM, Jeff Becker
> <jeffrey.c.becker at nasa.gov <mailto:jeffrey.c.becker at nasa.gov>> wrote:
>
> Hi,
>
> On 07/18/2015 05:17 PM, Berk Geveci wrote:
>> To elaborate on what Andy said, you can stick your entire VTK
>> script in side the programmable source. Which you can use in
>> batch mode as well. If you want to keep your VTK script as a
>> reusable code, I recommend making it a module that you import
>> from the programmable source.
>
> So if the last thing my VTK script does is set up a vtkImageImport
> object should that work as a programmable source (I tried but
> ParaView didn't get any bounds or scalar field name even though I
> set them up)? If not, what else do I need to do? Thanks.
>
> -jeff
>
>
>>
>> Best,
>> -berk
>>
>> On Tue, Jul 14, 2015 at 3:43 PM, Jeff Becker
>> <jeffrey.c.becker at nasa.gov <mailto:jeffrey.c.becker at nasa.gov>> wrote:
>>
>> On 07/14/2015 12:20 PM, Andy Bauer wrote:
>>> ParaView's Python Programmable Filter and Python
>>> Programmable Source use VTK Python wrapping to do the work
>>> while making it available in Paraview. Check out
>>> http://www.paraview.org/Wiki/Python_Programmable_Filter.
>>
>> Thanks. That's pretty neat, but what I'd really like is a
>> standalone script I can run in batch mode (no interaction) to
>> generate one png file per timestep for all the timesteps in
>> the data set. Essentially I'd like a way to construct a
>> pipeline between the vtk.vtkImageImport() object and the
>> ParaView renderView1.
>>
>> -jeff
>>
>>>
>>> On Tue, Jul 14, 2015 at 2:59 PM, Jeff Becker
>>> <jeffrey.c.becker at nasa.gov
>>> <mailto:jeffrey.c.becker at nasa.gov>> wrote:
>>>
>>> Hi.
>>>
>>> As a proof of concept, I have a vtk script that reads
>>> some binary data, and produces a vti file. I then have a
>>> second script that I generated using ParaView's tracing
>>> facility while viewing the data. Now I'd like to combine
>>> them, so as to eliminate the intermediate file, i.e., go
>>> from binary data directly to rendering. To be specific,
>>> my vtk script ends like this:
>>>
>>> steinbmag = steinbmag.reshape(nx,ny,nz).T
>>>
>>> vtkImporter = vtk.vtkImageImport()
>>> vtkImporter.CopyImportVoidPointer(steinbmag,
>>> steinbmag.nbytes)
>>>
>>> vtkImporter.SetDataScalarTypeToFloat()
>>> vtkImporter.SetNumberOfScalarComponents(1)
>>> vtkImporter.SetDataExtent(0, nx-1, 0, ny-1, 0, nz-1)
>>> vtkImporter.SetWholeExtent(0, nx-1, 0, ny-1, 0, nz-1)
>>> vtkImporter.SetScalarArrayName('B field magnitude')
>>>
>>> writer=vtk.vtkXMLImageDataWriter()
>>> writer.SetFileName(out_fname)
>>> writer.SetInputConnection(vtkImporter.GetOutputPort())
>>> writer.Write()
>>>
>>> and the ParaView script starts like this:
>>>
>>> from paraview.simple import *
>>>
>>> paraview.simple._DisableFirstRenderCameraReset()
>>>
>>> filepfx = '/nobackup/jcbecker/steinBmag_0'
>>> infile = filepfx+'.vti'
>>> # create a new 'XML Image Data Reader'
>>> steinBmag_ = XMLImageDataReader(FileName=infile)
>>> steinBmag_.CellArrayStatus = []
>>> steinBmag_.PointArrayStatus = ['B field magnitude']
>>>
>>> # get animation scene
>>> animationScene1 = GetAnimationScene()
>>>
>>> # update animation scene based on data timesteps
>>> animationScene1.UpdateAnimationUsingDataTimeSteps()
>>>
>>> # get active view
>>> renderView1 = GetActiveViewOrCreate('RenderView')
>>> # uncomment following to set a specific view size
>>> renderView1.ViewSize = [1090, 761]
>>>
>>> I'm looking at the VTK examples and pvpython
>>> documentation, but any help is appreciated.
>>>
>>> Thanks.
>>>
>>> -jeff
>>> _______________________________________________
>>> Powered by www.kitware.com <http://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 ParaView
>>> Wiki at: http://paraview.org/Wiki/ParaView
>>>
>>> Search the list archives at:
>>> http://markmail.org/search/?q=ParaView
>>>
>>> Follow this link to subscribe/unsubscribe:
>>> http://public.kitware.com/mailman/listinfo/paraview
>>>
>>>
>>
>>
>> _______________________________________________
>> Powered by www.kitware.com <http://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 ParaView Wiki at:
>> http://paraview.org/Wiki/ParaView
>>
>> Search the list archives at:
>> http://markmail.org/search/?q=ParaView
>>
>> Follow this link to subscribe/unsubscribe:
>> http://public.kitware.com/mailman/listinfo/paraview
>>
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/vtkusers/attachments/20150721/f7b5be0f/attachment.html>
More information about the vtkusers
mailing list