如何解决用 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 举报,一经查实,本站将立刻删除。