問題タブ [rgdal]

UTM33 プロジェクションで作成されたラスター ファイルがあります。これを WGS84 プロジェクションに変更して、他のレイヤー (シェープファイルとラスター) と一致させたいと考えています。次に、他のファイルと一致するように、上記のラスター (投影を変更したばかりのもの) の範囲も変更する必要があると思います。Rでこれらの両方を行うにはどうすればよいですか? raster および rgdal パッケージを使用する必要があることはわかっており、ここで説明されている他のトピックを調べましたが、まだ解決できません。あらゆるアドバイスをありがとうございました！

R のプロジェクトで複数のシェープファイルを操作する必要があります。

したがって、これらのファイルをリストにロードするために readOGR 関数を繰り返したいと思います。後で、plyr を使用するか、ループを作成して、リスト内のすべてのシェープファイルに対して同じ操作を実行できます。

最小限の再現可能な例を次に示します。

エラーメッセージは

空間ポリゴン オブジェクトと空間ポイント オブジェクトがあります。後者は、データフレーム内の xy 緯度データ ベクトル (それぞれ緯度と経度と呼ばれます) から作成されましたが、前者は rgdal を使用して直接 R に読み込まれました。私のコードは次のとおりです。

この時点では、どちらの空間オブジェクトも投影されません。これらのオブジェクトを次のように投影すると:

両方のオブジェクトが投影され、両方とも次のように plot() 関数で正常にマップされます。

ただし、それらを一緒にプロットしようとすると:

私は最初のプロットのみを取得し、rats オブジェクトを boros に重ねません。ただし、これは大きな問題です。エラー メッセージが表示されないため、これら 2 つの空間オブジェクトが互いに通信できる間の切断がどこにあるのかを判断できませんでした。どちらのコマンドもエラーや警告なしでスムーズに実行されますが、プロットは 1 つしか残っていません。そして、proj4string() で各オブジェクトの射影をチェックすると、各オブジェクトに対して同じ射影が返されます。

plot() を使用して一緒にマッピングできるように、CRS と投影が一致する 2 つの空間オブジェクトを作成するさまざまな方法を試すために、数日間にわたって何時間も費やしてきました。ちなみに、私がとったアプローチの 1 つは、ArcGIS でラット オブジェクトのシェープファイルを作成することでした。これは、ファイルを作成するのにうまく機能しました。しかし、R で 2 つの空間オブジェクトを連携させることができないという問題がまだ残っています。どんな助けでも大歓迎です。上記の 2 つのデータ ファイルを含むドロップボックス リンクも含めました: https://www.dropbox.com/sh/x0urdo6guprnw8y/tQdfzSZ384

ArcMap を使用して、空間データ サーバー上の GIS データにアクセスしています。R内で同じことを行う方法を理解したい.

シェープファイルを R に読み込む方法を知っています。maptools と rgdal を使用して、ローカルに保存されたシェープファイルを開いてマップすることに成功しました (例: http://www.nceas.ucsb.edu/scicomp/usecases/ReadWriteESRIShapeFiles )

私の問題は、データがローカルに保存されておらず、アプリケーション サーバー上にある場合です。私はそれがOracleデータベースだと信じています。1. サーバー 2. インスタンス (数字) 3. データベース 4. ユーザーと 5. パスワードに関する情報を受け取りました。通常、例を含めますが、外部ユーザーがサーバーにアクセスできるかどうかは疑わしいです。

たとえば、Rでローカルファイルを読み取ってプロットする方法は次のとおりです

「。」ローカル ディレクトリを指します。サーバー上のデータにアクセスするには、これをどのように変更しますか? コードの実際の構文が役立ちます。

私は一連の地震の経度と緯度を持っています。陸と海に分けられるようにしたいです。

それを行うr関数はありますか？