import numpy as np import matplotlib.pyplot as plt from scipy.interpolate import griddata

随机生成南海区域风速数据

x = np.random.uniform(100, 130, 100) # 经度 y = np.random.uniform(0, 20, 100) # 纬度 z = np.random.uniform(0, 20, 100) # 风速

创建网格坐标

xi = np.linspace(100, 130, 100) yi = np.linspace(0, 20, 100) xi, yi = np.meshgrid(xi, yi)

内插函数

zi = griddata((x, y), z, (xi, yi), method='cubic')

绘制等高线图

plt.contourf(xi, yi, zi, cmap='rainbow') plt.colorbar() plt.scatter(x, y, c=z, cmap='rainbow', s=10) plt.title('南海区域风速数据插值') plt.xlabel('经度') plt.ylabel('纬度') plt.show()