我目前能够在 ParaView 中为模拟的每个时间步正确地可视化一个 .vtp 文件,并为每个时间步打印一个屏幕截图。我想批量执行此操作,但我想为每个状态(视点、应用的过滤器等)保持相同的状态。我已经将状态保存到一个 .psvm 文件中,并且我尝试编写一个 python 脚本,在 pvbatch 运行之后,它将(希望)打印屏幕截图。但是,不幸的是,它不起作用。我试图通过处理状态文件并进行搜索和替换来更改状态中的文件名,但它仍然无法正常工作。例如,它只绘制第一个数据输入,即使当前文件不同(尽管 GetSources() 显示一个总是增加的源列表)。我在雪豹中使用 ParaView 3.14.0。我相信这很容易,但是我对关于 python 和 paraview 的大量信息感到不知所措,而没有提到这个特定的问题。拜托,拜托,非常欢迎任何建议,如果之前已经回答过,我很抱歉(我查看了谷歌、paraview 邮件列表和这里)。下面是我的脚本,也可以在http://pastebin.com/4xiLNrS0。此外,您可以在http://goo.gl/XjPpE中找到一些示例文件和状态。
#!/bin/python
import glob, string, os, commands
from paraview.simple import *
#help(servermanager)
# vtp files are inside the local subdir DISPLAY
files = (commands.getoutput("ls DISPLAY/data-*.vtp | grep -v contacts")).split()
# process each file
for filename in files:
fullfn = commands.getoutput("ls $PWD/" + filename)
fn = filename.replace('DISPLAY/', '')
#os.system("cp ../dem_git/addons/paraview_state.pvsm tmp.pvsm")
os.system("cp ~/Desktop/state.pvsm tmp.pvsm")
os.system("sed -i.bck 's/DATA.vtp/" + fullfn.replace('/','\/') + "/1' tmp.pvsm") # replace first intance with full path
os.system("sed -i.bck 's/DATA.vtp/" + fullfn.replace('/','\/') + "/1' tmp.pvsm") # replace second intance with full path
os.system("sed -i.bck 's/DATA.vtp/" + fn + "/1' tmp.pvsm") # replace third with just the filename path
servermanager.LoadState("tmp.pvsm")
pm = servermanager.ProxyManager()
reader = pm.GetProxy('sources', fullfn)
reader.FileName = fullfn
reader.FileNameChanged()
reader.UpdatePipeline()
view = servermanager.GetRenderView()
SetActiveView(view)
view.StillRender()
WriteImage(filename.replace(".vtp", ".png"))
os.system("rm -f tmp.pvsm")
os.system("rm -f tmp.pvsm.bck")
Delete(reader)