在我的应用程序中,我有一个从 meshgrid 命令获得的值矩阵及其坐标(经度、纬度)。
我想根据经度和纬度限制提取该矩阵的特定子区域。我已经尝试过这个解决方案,但它不起作用。我需要三个矩阵作为输出,一个用于数据,另外两个用于网格。
Lons, Lats = meshgrid(X, Y)
indexes = np.where((Lons < MLon) & (Lons > mLon) & (Lats < MLat) & (Lats > mLat))
newLons = Lons[indexes]
newLats = Lats[indexes]
newData = Data[indexes]
获得的新值是一维数组而不是矩阵。
我怎样才能解决这个问题?
没有任何保证np.where
的观点是,您将提取构成连续矩形子矩阵的值,因此它将返回平坦的值。你可以重塑它们,但为此你需要弄清楚它们的形状。更好、更通用的解决方案是找到边界框,然后提取它:
Xspan = np.where((X < MLon) & (X > mLon))[0][[0, -1]]
Yspan = np.where((Y < MLat) & (Y > mLat))[0][[0, -1]]
# Create a selection
sel = [slice(Xspan[0], Xspan[1] + 1), slice(Yspan[0], Yspan[1] + 1)]
# Extract
newLons = Lons[sel] # == Lons[Xspan[0]:Xspan[1]+1, Yspan[0]:Yspan[1]+1]
newLats = Lats[sel]
newData = Data[sel]
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)