用 Cartopy 绘制 csv

问题描述

我对 Python 还很陌生，所以请原谅我的经验不足。我尝试用 cartopy 创建一个轮廓图。我以前用 netCDF4 数据做过，但这次我必须处理普通的文本数据格式。问题是 pcolormesh 需要一个二维数组作为 C 组件。

我的数据看起来像这样：

       lon    lat    Q_upt
0    28.18 -25.73 -999.990
1    27.98 -25.83    0.388
2    27.78 -25.97    0.212
3    27.61 -26.13    0.082
4    27.47 -26.31    0.000

我尝试使用 interpolate.griddata 对数据进行插值，但这根本不起作用。

header_list = ['time','lon','lat','Q_upt']
dt = pd.read_csv('~/1996062812_1996070812_PRCupt.txt',skiprows=4,sep='\s+',names=header_list)

#create DataFrame with required variable
dtt = dt[['lon','Q_upt']]

#transform to numpy.array
df = dtt.to_numpy()

#create coordinate array
x = np.arange(-180,180,0.1)
y = np.arange(-90,90,0.1)

X,Y = np.meshgrid(x,y)

#interpolate the variable into an 2d array dependent on chosen coordinates
Z = interpolate.griddata((df[:,0],df[:,1]),2],(X,Y),method='nearest')

我的下一个想法： 我创建了一个带有经度和纬度的 m x n 矩阵，现在我必须将它与每一行的 lon lat 值进行比较，并将其替换为匹配的 Q_upt。但我不知道如何编码。 如果有人可以帮助我，我会很高兴。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）