我正在使用网格上的海洋数据,其中地球的每个度数都对应一个数据点,即我可以使用data[lat_index, lon_index]
访问一个点,其中lat_index
和lon_index
对应到网格点的经纬度。我也有对应于网格点的数组:
lat_arr = np.arange(-90,91)
lon_arr = np.arange(0,360)
我正在寻找一种在子午线上切片数据的方法。例如,我希望能够检索与对应于-20S到10N和-30W(即360-30 E)到30E的矩形相对应的切片。显然data[-20:11, 330:31)
不能解决问题,因为330:31
不检索任何数据。有没有一种pythonic的方法可以在lon_index=0
的子午线上切出我的数组?我正在处理大型数据集,因此效率是关键,并且我不想创建许多新数组。