[Paraview] Combining vtk and paraview python scripts

[Jeff Becker]

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/paraview/attachments/20150720/5dccf4cb/attachment.html>