問題タブ [sp]

R を使用して、poly_center (pracma パッケージ) を使用して、このファイルで読み込まれた図によって定義されたポリゴンの重心を計算しています。

奇妙なことに、私が得ている結果はポリゴンの境界から外れており、中心点が x (経度) と y (緯度) 軸/座標の最小値の両方よりも低くなっています (以下のコードを参照)。私が間違っていること、またはここで何が起こっているのか、誰かが私にヒントを与えることができますか?

アイデアを与えるために、マップの空間ポリゴンを定義するエントリを含むリスト (Stata dta ファイルからインポート) があります。

これを R でマップを作成するために使用できるデータ構造に変換したいと思います (より大きな困難や広範な専門知識を必要とせずに、これまでに行ったことがない)。 sp パッケージの SpacialPolygons タイプは最も簡単な選択。さらに、定義から、SpacialPolygons メソッド (このパッケージで定義されています。パッケージのドキュメント、79 ページを参照) は、リストからこのデータ型に変換する正しいメソッドであるように見えます。

残念ながら、この方法で作業するのはそれほど簡単ではなく、助けが必要です。私の（素朴な）試みは、私が理解できないエラーを引き起こし、Google検索で興味深い結果をもたらしません。

上記のリストを SpacialPolygon に変換できますか? もしそうなら、どのように？そうでない場合、代わりにどの形式を選択すればよいですか? ありがとう。

ポリゴン シェイプ ファイルと一連のデータ フレームがあります。私がやりたいことは、シェイプ ファイルとそれらすべてのデータ フレームの間で「空間結合」を行うループを作成することですが、エラーが発生し続けます。これに関するアイデアはありますか？

これが私のループのコードです

以下は、3 つのステップで再現可能な例です。

(1) ランダムな多角形から正格子状の多角形を作成する

(2) データフレームの作成

(3) ループしてみる

空間結合を 1 つずつ行うと正しい結果が得られますが、データ セットが多すぎるため、これでは問題は解決しません。

フィーチャ ジオメトリ属性を持つ MySQL データベースから取得したデータ フレームがあります。たとえば、geom 属性の 1 つのレコードの値は次のとおりです。

テーブル内のすべてのフィーチャはポリゴンです。

検索してきましたが、これらの座標文字列から空間データ フレームを作成するきれいな方法が見つかりません。

空間ポリゴンと空間ラインが交差するかどうかを迅速に判断する必要があります。現在、ポリゴンを空間線に変換して使用してgIntersection()います。より迅速な方法を提案できる人はいますか? おそらく、空間線などの代わりにラスターを使用しています。これを何千回も行う必要があります。

gDifference2 つのポリゴンの間で実行しようとしていますが、エラーが発生しました。

だから私はこれをしました：

しかし、今は別のエラーが発生しています:

これを修正する方法はありますか？ArcGIS の Erase 関数でまったく同じことを行うと、まったく問題なく動作します。

Rでspパッケージを使用しています。あまり違いがないように思えますか？

SpatialPointsDataframe でより多くの属性を格納できるという唯一の違いはありますか?! もしそうなら、spatialPointsを作成する迂回路を使わずに、既存のデータフレームからSpatialPointDataframeを作成できますか(このデータフレームに座標が存在する場合)?