python等值面平滑_最佳拟合矩形网格到平滑的三维曲面

如果有定义曲面的节点，则意味着坐标和相应值的栅格不规则。因此，您可以从中生成一个三角剖分(您用来显示这些填充轮廓的工具很可能在屏幕后面使用相同的方法)。在

以下是同一个示例中的基本步骤，我对此进行了注释，但是matplotlib贡献者值得称赞：import matplotlib.pyplot as plt

import matplotlib.tri as mtri

import numpy as np

# Create triangulation.

coords, earthquake_fault = get_coordinate_data() # to be filled in by you

x = coords['x']

y = coords['y']

triang = mtri.Triangulation(x, y)

# Interpolate to regularly-spaced quad grid.

z = earthquake_fault # the "height" data

xi, yi = np.meshgrid(np.linspace(x.min(), x.max() 20), np.linspace(y.min(), y.max(), 20))

interp_lin = mtri.LinearTriInterpolator(triang, z)

zi_lin = interp_lin(xi, yi)

# Plot the triangulation.

plt.subplot(121)

plt.tricontourf(triang, z)

plt.triplot(triang, 'ko-')

plt.title('Triangular grid')

# Plot linear interpolation to quad grid.

plt.subplot(122)

plt.contourf(xi, yi, zi_lin)

plt.title('Rectangular grid')