Volume Process
Load the Data
# Load the a TensorMesh and some already processed model data
mesh = discretize.TensorMesh.readUBC('craig_chile.msh', directory='data/Craig-Chile')
models = {'lpout': mesh.readModelUBC(fileName='Lpout.mod', directory='data/Craig-Chile')}
空间网格
属性数据:
Volume
cell_sizes = mesh.toVTK(models=models).compute_cell_sizes()
cell_sizes.plot(scalars='lpout')
Slice
Sample
grid=pyvista.create_grid(slices_v[0], dimensions=(256, 256,1))
result = grid.sample(slices_v[0])
temp=result.get_array("lpout").reshape(256,256)*10000
data=temp.astype(np.int16)
data
写入图片:
from numpngw import write_png
import numpy as np
write_png("data/s0.png",data)
Extract