問題タブ [rgeo-shapefile]

座標参照系のない数百のシェープファイルがあります。私の目標は、空間ポリゴンを WorldClim ラスター レイヤーにオーバーレイすることです。以前はこのアプローチを問題なく使用していました。ただし、今回はシェープファイルの座標が奇妙です。境界ボックスとポリゴン内の座標の各座標は、8 桁の数字で構成され、小数点以下を区切るためにコンマやドットは使用されません。これは、いずれかの形状のバウンディング ボックスです。

WorldClim ラスター レイヤーのバウンディング ボックスとは明らかに異なります。

プロット コマンドを使用してシェープファイルをラスター レイヤーに重ね​​ようとしても、何も起こりません。

これは投影の問題であることを理解しました。次に、CRS 関数を使用して、シェープファイルに WorldClim ラスター レイヤーの同じ座標系を割り当てようとしました。ただし、結果は同じままです (つまり、シェープファイルはラスターを超えません)。シーケンスでは、rgdal パッケージの spTransform 関数を使用して、シェープファイルの座標を再投影しようとしました。ただし、シェープファイルには参照システムがないため、関数は機能せず、ラスター レイヤーと一致させるためにシェープファイルを再投影する方法がわかりません。私はこの問題に対処する方法について数日間調査してきましたが、参照システムがないことが問題の重要なポイントであると考えています。しかし、私はこの問題を克服できず、誰かがこの状況に対処する方法を手伝ってくれるかどうか知りたい.

次のコードは、座標ペアが特定のポリゴンに属しているかどうかを確認します。ループなどを使用して、より多くの座標データを確認する必要があります。

コードを呼び出すと、すべての値がチェックされますが、国外の座標が 1 つしかない場合、リストのすべてのポイントのステータスは「false」になります。