CFD Online Discussion Forums - [General] Generating isosurface from vtk file with python script

I have a vtk file that is generated via auxiliary code and would like to generate isosurfaces of this (3D) data. I can read this in correctly and create the necessary isosurfaces manually, and can even do the same with the python shell when the gui is up by doing the following:

Code:

import paraview.simple as ps



# reads in a vtk file

test = ps.OpenDataFile("/path/to/test.vtk")

c = ps.Contour(Input=test, Isosurfaces=[0.5])



ps.Show(test)

ps.Show(c)

ps.Render()

ps.WriteImage("/path/to/check.png")

If I enter these commands line-by-line into the python shell, the correct isosurface is generated; however, if I run the script by using either the "Run Script" button or pvbatch, the isosurface is not created. In both cases, the vtk data is (test) is generated correctly. Interestingly, when I click the "Apply" button (with no other input on my end) on the isosurface in the paraview window, the isosurface is also generated correctly, which means that I must be missing some scriptable "Apply" button.

The closest solution to this can be found here: https://colorsfordirectors.wordpress...ython-scripts/

Code:

# 1. import the paraview.simple module (here renammed 'pvs')

import paraview.simple as pvs

 

# 2. create the filter with desired input

my_filter = pvs.FilterName(Input=my_data)

 

# 3. set properties

my_filter.PropertyOne = [0, 0, 0]

 

# 4. update the filter (optionnal)

my_filter.UpdatePipeline()

I added the UpdatePipeline() function to my script and this seemed to do nothing.

I will test this out on another paraview installation / linux distribution tomorrow when I get back into the office. I am currently running Arch and paraview 5.4.1 .

If anyone has any advice, I would love to hear it! I'm new to paraview python scripting and am having a bit of trouble figuring everything out. Thanks for reading and any help you can provide!