我试图仅在大陆上绘制数据。我之前问过一个关于这个的问题,我在这里得到了答案。现在,当我在代码中实现它时,程序将运行而不会出现错误。
但是,它永远运行。当我使用ctrl + C
强制停止它时,它会追溯到"if not (map.is_land(X[i], Y[i])):"
,所以我认为这里有问题。但我无法弄清楚。
这是涉及挑选非大陆数据并删除它们的代码部分:
X, Y = map(Lon,Lat)
ocean = []
for i in range(len(X)):
if not (map.is_land(X[i], Y[i])):
ocean.append(i)
X_new = np.delete(X, ocean)
Y_new = np.delete(Y, ocean)
HDO_new = np.delete(HDO, ocean)