跨子午线切片经度数组的Pythonic方法

时间:2019-03-11 15:31:53

标签: python numpy

我正在使用网格上的海洋数据,其中地球的每个度数都对应一个数据点,即我可以使用data[lat_index, lon_index]访问一个点,其中lat_indexlon_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的子午线上切出我的数组?我正在处理大型数据集,因此效率是关键,并且我不想创建许多新数组。

0 个答案:

没有答案