を使用してデジタル標高モデルからルートの標高データを抽出しています
my.elev <- over(new.points, mygrid)
- new.points は、mygrid の CRS で変換された約 7000 ポイントの座標 (経度/緯度) を持つ SpatialPoints オブジェクトです。
- mygrid は、800 万を超える要素を持つ SpatialGridDataFrame です
(詳細については、前の質問を参照してください) 。
my.elev にいくつかの NA 値があるので、コードをデバッグしたところ、new.points のほぼすべてのポイントが複数回繰り返されていることがわかりました (私のルートでは、いくつかのセグメントが 2 回交差しています)
。 .elev 値が正しい
- 2 番目の値は NA 値 (場合によってはまったく異なる値) です。
new.points で重複した値を削除する問題は簡単に解決できますが、over 関数が同じポイントに対して同じ値を返さないのはなぜでしょうか。