我只想在大陆上绘制数据。我以前问过一个问题,我得到了答案here。 现在,当我在代码中实现它时,程序将运行而不会给出错误。
然而,它永远运行。当我使用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)
由于使用了
for
循环,程序无论如何都应该停止,唯一的问题是X
有多大。我只能提出两个建议:1)在输入for cicle之前,尝试打印
len(X)
的值,以查看其大小2)如果使用python2.x,请尝试使用}
xrange
,而不是{还要记住,使用包含大约十万个元素的list似乎很慢。 如果可以的话,最终尝试使用较小的数据样本
相关问题 更多 >
编程相关推荐