問題タブ [maptools]

R の GIS は初めてで、緯度/経度または UTM 座標をシェープファイルに追加しようとしています。ここからシカゴ市の境界 (City_Boundaries.shp) をダウンロードしました: http://www.cityofchicago.org/city/en/depts/doit/supp_info/gis_data.html

maptools と rgeos ライブラリをロードしました:

データを R に取り込み、ゾーン 16T の UTM コードを追加しようとしました。

ただし、結果のデータは意味がありません。「city」内の「coords」スロットを見ると、座標値は、たとえば、X = 1092925、Y = 1944820 です。外部ツール ( http://home.hiwaay.net/~taylorc/toolbox/geography/geoutm.html ) を使用して GPS 座標を見つけました。結果は、lat=17.511、long=-81.42 でした。これは、ジャマイカとホンジュラスの海岸の間のどこかにあります。シェープファイルの座標は独自の宇宙に存在するようです。シェープファイルという名前が示すように、座標は都市の形状を正確に表していますが、これらの座標は地球上に自動的にマッピングされません。

私の最終的な目標は、ジオタグが付けられた多くのイベントがシカゴ内で行われたかどうかを判断することです。以下のように、緯​​度/経度でリストされているこれらのポイントを UTM に変換することに問題はありません。

また、それがより理にかなっている場合は、イベント データを元の形式で処理することも考えています。最終的には、シカゴの都市境界をポリゴンに変換して、イベント データと照合できるようにする必要があります。このようなもの：

ポイントが市内にあるかどうかを確認できるように、シェープファイルとイベント データを共通の (正確な) 参照フレームに入れるにはどうすればよいですか? あなたが提供できるどんな助けも大歓迎です。

R の等高線から標高プロットを作成できるようにしたいと考えています。シェープ ファイルを使用するのは初めてです。

現時点 では、英国全土の .shp ファイルを提供するデータをここからダウンロードしました。

また、英国のトポロジーを要約する等高線も提供します。

標高プロットでは、x、y、z 値を示すデータ出力を生成するために、等間隔の点 (互いに 100 m 離れている) の1data.frameつまたは1 つが必要です。data.tableここで、x と y は緯度と経度 (または東座標と北座標) を表し、z は高さ (海抜メートル) を表します。

おそらく、補間を自動的に実行するツールがいくつかあると思いますが、地理空間データでどのように機能するかはわかりません.

これが私の基本的なスタートです...

これはおそらく簡単な質問ですが、ここからボーダーをダウンロードして、各出力エリアを地図上に表示したいと思います。2011 OA を選択し、完全な解像度、領域の範囲を取得します。

私はそれをアップしました...そして、次を使用してshpファイルを読み込みました

各出力領域の座標をプロットまたは表示するにはどうすればよいですか?

私はRが初めてで、現在コロプレスマップを作成しようとしています。私の問題は、休憩が意味をなさないことです。私がやりたいことは、等しいブレークを作成することです:

59%未満 60-69% 70-79% 80-89% 90%以上

ただし、代わりに得られるのは:
60% 未満 60-70% 70-80% 80-90% 90% 以上

これは、values=80 が両方のカテゴリで取り上げられているということですか?

これらの休憩を作成するための私のコードは次のとおりです。

注: Edzer Pebesma の提案により、この質問は R-sig-geo (こちら) にクロスポストされ、いくつかの良い反応がありました。

を使用すると、次の予期しない結果が発生しましたcheckPolygonsHoles。

ここで厄介な副作用の 1 つは、元のオブジェクトwrld_simplも変更されることです。この結果は私にはバグのように見えますか、それとも何か見逃していますか?

shape_with_holePS: 以前に編集されたオブジェクトcheckPolygonsHolesは、引き続き奇妙な動作をします: