我正在绘制一些 NETCDF 数据,由 Matplotlib 中的 xarray 使用 Cartopy 包装器进行地图投影处理。当我生成填充等值线图(contourf)时,我最终在经度零(我的数据经度开始和结束的地方)处出现了“接缝”。是否有一种内联方法可以在经度[-1]和经度[0]之间进行插值,也许可以通过将数据视为周期性数据或其他方式?
请参阅附图示例:
您应该使用
add_cyclic_point
函数向数据数组及其经度坐标添加循环点: http://scitools.org.uk/cartopy/docs/v0.16/cartopy/util/util.html#cartopy .util.add_circular_point.
还有更通用的 add_cyclo 函数。
使用示例:
import cartopy.util as cutil
lon2d, lat2d = np.meshgrid(lon,lat)
cdata, clon2d, clat2d = cutil.add_cyclic(data,lon2d,lat2d)
然后将其用于绘图。