問題タブ [rgdal]

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 は高さ (海抜メートル) を表します。

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

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

ダウンロードしたラスターがありますが、なんらかの理由で読み込み時に「proj4string」がありません。

ドキュメントには地理参照されているとありますが、データを定義するためにカスタムの「proj4string」を作成できるかどうか疑問に思っています。ARCGIS にもアクセスできますが、よくわかりません。投影の中心をリセットする方法がわかりません。

提供されるジオリファレンス情報は次のとおりです。

geos、proj-4.8.0、および gdal-1.9.2 をインストールしましたが、R パッケージ rgeos をインストールした後、rgdal のインストールに失敗します。

これが私のホームディレクトリです。

. . . そして私が実行しているエラーinstall.packages("rgdal", type="source")：

** パッケージ インデックスの構築 ** ビネットのインストール ** インストールされたパッケージをロードできるかどうかのテストredhat-linux-gnu-library/3.1/rgdal/libs/rgdal.so': libgdal.so.1: 共有オブジェクト ファイルを開けません: そのようなファイルやディレクトリはありません エラー: 読み込みに失敗しました 実行が停止しました エラー: 読み込みに失敗しました * '/ を削除しますmainstorage/michaeld/R/x86_64-redhat-linux-gnu-library/3.1/rgdal'</p>

ダウンロードされたソース パッケージは '/tmp/Rtmpt51iFG/downloaded_pa​​ckages' にあります

このエラーは、依存ライブラリと関係があると推測されます。. . しかし、それは私が集めるすべてです。

ここでsudo yum install libgdal1 libgdal1-dev libgeos libgeos-dev提案されているように実行すると、次の結果が得られることに注意してください。

パッケージ libgdal1 がありません。パッケージ libgdal1-dev がありません。パッケージ geos-3.4.2-1.3.amzn1.x86_64 は既にインストールされており、最新バージョンのパッケージ libgeos-dev は利用できません。

devしたがって、おそらくより良い質問は、これらのパッケージをどのようにインストールするのですか?

方向性をありがとう。

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

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

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

上記の 2 つのパッケージ "rgdal" と "raster" を使用して、GeoTiff ラスター ファイルをトリミングしたいと思います。結果として得られる出力 tif の品質が非常に悪く、カラーではなくグレースケールであることを除いて、すべてが正常に機能します。元のデータは、スイス連邦地形局の高品質ラスター マップです。サンプル ファイルは、ここからダウンロードできます。

これは私のコードです：

この例を再現するには、サンプル データをダウンロードし、「c:/files/」フォルダに展開します。奇妙なことに、サンプル データを使用すると、トリミングされた画像の品質は問題ありませんが、グレースケールのままです。

オプション「datatype」、「format」を使用して遊んでいましたが、それでどこにも行きませんでした。誰でも解決策を指摘できますか？入力データにさらに情報を提供する必要がありますか?

編集:ジョシュの例は、サンプルデータ2でうまく機能します。残念ながら、私が持っているデータは古く、多少異なっているようです。次の GDALinfo を読んだ場合、どのオプションを選択するか教えてください。

rgdal ライブラリを使用してシェープファイルを読み込もうとしていますが、うまくいきません。

次の構文を使用してインポートしようとすると:

エラーが発生しました

問題を診断して修正するにはどうすればよいですか? どうもありがとう。

R: Crop GeoTiff Raster using packages "rgdal" and "raster" で提起された質問と同様に、元のカラー テーブルを維持しながら、パッケージ "rgdal" および "raster" を使用して、スイス連邦地形局から地図をトリミングしようとしています。 . 単一のバンド *.tif ファイルの場合、トリミングされたイメージはカラー テーブル情報を失うため、適切に表示されません (結果のイメージはほとんど黒になります)。

入力ファイルはここからダウンロードでき、「C:/files」フォルダーに展開する必要があります。これはコードです：

前述の投稿で提示されたソリューションは、3 バンドの *.tif に対してのみ機能し、1 バンドの *.tif (サンプル ファイルなど) に対しては機能しませんでした。

機能するはずの解決策は、カラー テーブルを含むシングル バンドの rasterLayer を 3 バンド RGB の rasterStack に変換することです (前述の投稿のコメントで概説されているように)。

ただし、カラー テーブルを保持しながら、単一バンド *.tif を 3 バンド RGB rasterStack に変換する方法がわかりません。この変換を行う方法を知っている人はいますか、または問題を解決するためのより良いアイデアを持っている人はいますか?