我正在处理此script,需要查询我通过netCDF4
(第一个示例)获取的数据对象,以检索特定{{1}的值}和latitude
坐标。我不确定如何使用纬度/长度(以度为单位)索引示例中的数据对象,或者如何将坐标映射到meshgrid并从那里查询。想法有人吗?
答案 0 :(得分:0)
根据NOAA示例脚本并在ww3邮件列表的帮助下,我想到了这个:
为了从数据对象中获取值,您可以
from mpl_toolkits.basemap import Basemap,interp
lat = GPSLat
lon = (360 + GPSLon)
value = interp(data, lon, lat, np.asarray( [[ (360+ GPSLon) % 360 ]] ), np.asarray( [[ GPSLat ]] ), checkbounds=True, masked=True, order=1)