如何使用Metpy.interpaly从一组数据中插值温度

Question

我正在尝试使用Cartopy和Metpy.interpaly生成地图，我希望最终产品与下面的图片相似从ndawn网站拍摄的图像。我有一组具有坐标和温度的数据（LAT，LON，TEMP）：

lat：[49.8134，49.8134，49.8134，49.8134，49.8134，49.8134，49.8134，49.8134， 49.8134，49.8134，49.8134]

lon：[-100.3721，-100.3721，-100.3721，-100.3721，-100.3721， -100.3721，-100.3721，-100.3721，-100.3721，-100.3721]

tem：[-8.45，-4.026，-5.993，-3.68，-7.35，-7.421，-6.477，-8.03，-8.03， -3.834，-13.04]

我正在尝试插入温度并使用Contourf将其抽出。但是，我不知道如何使用Metpy的插值函数将它们插入。有人可以帮忙吗？

Answer 1

这是一个例子。我必须从上面更改LAT/LON值，因为这是一个仅在上面重复的点（您无法插入

import cartopy.crs as ccrs
import matplotlib.pyplot as plt
import metpy.interpolate as minterp
from metpy.plots import StationPlot
import numpy as np

# Generate some locations for the data--originals were a single point
tem = np.array([-8.45, -4.026, -5.993, -3.68, -7.35, -7.421,
                -6.477, -8.03, -3.834, -13.04])
lat = 49 + np.random.rand(len(tem))
lon = -100 - np.random.rand(len(tem))

# Interpolating directly on the sphere with lat/lon is problematic
# so interpolate to a grid on a LambertConformal projection
proj = ccrs.LambertConformal(central_latitude=49, central_longitude=-100)
x, y, _ = proj.transform_points(ccrs.PlateCarree(), lon, lat).T

# Interpolate to a grid with 5km spacing
x_grid, y_grid, tem_grid = minterp.interpolate_to_grid(x, y, tem, interp_type='barnes',
                                                       hres=5000)

# Create figure
fig = plt.figure(figsize=(8, 8))
ax = fig.add_subplot(1, 1, 1, projection=proj)

# Contourf plot
ax.contourf(x_grid, y_grid, tem_grid)

# Use MetPy's StationPlot to plot the text values
sp = StationPlot(ax, x, y, transform=proj, fontsize=16)
sp.plot_parameter('C', tem)

# Show grid locations
ax.scatter(x_grid, y_grid)

的插件：

有关更多信息，我建议您查找更多信息at 在METPY文档中此示例

的，因为间距不相等

M从LAT/LON中转换为LAT/LON，因为用度插值是有问题 在调用transform_points（）的调用中，使用platecarree（）是您与造型的传达方式，所给出的坐标为LON/LAT坐标。