from tvtk.api import tvtk
from tvtkfunc import ivtk_scene,event_loop
#读入plot3D数据
plot3d=tvtk.MultiBlockPLOT3DReader(xyz_file_name="combxyz.bin",#网格文件q_file_name="combq.bin",#空气动力学结果文件scalar_function_number=100,#设置标量数据数量vector_function_number=200#设置矢量数据数量
        )
plot3d.update()
grid = plot3d.output.get_block(0)#获取读入的数据集
con = tvtk.ContourFilter()#创建等值面对象
con.set_input_data(grid)#将网格与其进行绑定
#范围由数组Scalars.range决定,颜色也由scalars决定
con.generate_values(10,grid.point_data.scalars.range)#创建10个等值面
#默认映射表：最小值为红色，最大值为蓝色
#对映射器进行构造
m=tvtk.PolyDataMapper(scalar_range=grid.point_data.scalars.range,input_connection=con.output_port)
#对标量行为属性进行赋值，设置成新增加的数组的取值范围
a=tvtk.Actor(mapper=m)
#由于这十个等值面会相互嵌套，为观察等值面内部结构，修改Actor对象透明度为0.5
a.property.opacity=0.5
#绘制交互窗口
win=ivtk_scene(a)
win.scene.isometric_view()
event_loop()

from tvtk.api import tvtk
from tvtkfunc import ivtk_scene,event_loop
#读入plot3D数据
plot3d=tvtk.MultiBlockPLOT3DReader(xyz_file_name="combxyz.bin",#网格文件q_file_name="combq.bin",#空气动力学结果文件scalar_function_number=100,#设置标量数据数量vector_function_number=200#设置矢量数据数量
        )
plot3d.update()
grid = plot3d.output.get_block(0)#获取读入的数据集(StructureGrid)
#对数据集中的数据进行随机选取，每50个点选择一个点
mask = tvtk.MaskPoints(random_mode = True,on_ratio = 50)#为数据进行了降采样
#每50个点选择一个点，为数据进行了降采样。为了观察效果可以在shell中进行输出降采样的效果
mask.set_input_data(grid)#将grid与musk相连
#创建表示箭头的PolyData数据集
glyph_source = tvtk.ArrowSource()#ArrowSource可以改成ConeSource()出来的箭头变成圆锥,scale_factor=2设置防缩系数
#在Mask采样后的PolyData数据集每个点上放置一个箭头
#箭头的方向、长度和颜色由于点对应的矢量和标量数据决定
#在本例中箭头的方向表示速度的方向，大小和颜色表示密度
#箭头越大，该点标量值越大。箭头的颜色表表示标量值的大小
#红色对应的标量值越小，蓝色对应的标量值越大
glyph=tvtk.Glyph3D(input_connection=mask.output_port,scale_factor=4)#4表示符号的共同放缩系数
#tvtk的可视化技术，它输入数据的每个点都被拷贝一个符号。
#符号本身是通过Glyph3D的filter的第二个输入函数来接收vtkPolyData的类型数据
#然后通过ArrowSource在每个点上放置一个箭头
glyph.set_source_connection(glyph_source.output_port)
m=tvtk.PolyDataMapper(scalar_range=grid.point_data.scalars.range,input_connection=glyph.output_port)
a=tvtk.Actor(mapper=m)
#交互窗口绘制
win=ivtk_scene(a)
win.scene.isometric_view()
event_loop()

from tvtk.api import tvtk
from tvtk.common import configure_input
from tvtkfunc import ivtk_scene,event_loop
#读入plot3D数据
plot3d=tvtk.MultiBlockPLOT3DReader(xyz_file_name="combxyz.bin",#网格文件q_file_name="combq.bin",#空气动力学结果文件scalar_function_number=100,#设置标量数据数量vector_function_number=200#设置矢量数据数量
        )
plot3d.update()#让plot3D计算出其输出数据
grid = plot3d.output.get_block(0)#获取读入的数据集(StructureGrid)
outline = tvtk.StructuredGridOutlineFilter()#计算表示外边框（轮廓）的PolyData对象
configure_input(outline,grid)#调用将外框计算与数据集产生关联
#两个参数：一个是外框PolyData对象一个是流体数据grid
m = tvtk.PolyDataMapper(input_connection=outline.output_port)a = tvtk.Actor(mapper=m)#创建一个Action实体Actor
a.property.color= 0.3, 0.3, 0.3
#窗口绘制
win = ivtk_scene(a)