大陸のみにデータをプロットしようとしていました。以前にこれについて質問しましたが、ここに答えがあります。これをコードに実装すると、プログラムはエラーなしで実行されます。
ただし、それは永遠に実行されます。で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)