1

我目前能够在 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)
4

1 回答 1

0

我意识到这是一个老问题,但我最近遇到了完全相同的问题,也找不到任何答案。您需要做的就是添加Delete(view)afterDelete(reader)以使您的脚本正常工作。

于 2015-01-09T11:35:49.147 回答