問題タブ [rgdal]

シェープファイル オブジェクトの背後に「ラスター」オブジェクトをプロットするにはどうすればよいですか? 両方とも単独でうまくプロットされますが、ポイントはラスタ上にプロットされません:

Rを使用して、すでに色付きのポリゴンが含まれている別のシェープファイルの上にポリゴンのシェープファイルをプロットしようとしています。maptoolsとrgdalを使用しています。下部のシェープファイル（複数のポリゴンがある）を、4つのポリゴンがある上部のシェープファイル（私は思う）で区切られた領域にクリップしたいと思います。もちろん、plot（）内でcol =コマンドを使用して、一番上のポリゴンの内側を塗りつぶす方法を知っています。しかし、私がやりたいのは、上部のポリゴンの外側を白で塗りつぶし、下部の色付きのポリゴンを上部のポリゴンの内側から見えるようにすることです。これができたら、下のポリゴンの境界を再描画します。私はR2.13.2と最新版のmaptoolsとrgdal、およびそれらのすべての依存関係をWindows7マシンで使用しています。

これらのシェープファイルに依存するコードがあります：ftp： //ftpext.usgs.gov/pub/er/wi/la.crosse/McKann/bcr%20arcinfo%20files/

範囲マップの内側を色にすることができます

しかし、私がやりたいのは、外側を白にして（後で州と州の境界線を描画できるようにするため）、内側を透明にして、範囲内の州と州の色を確認できるようにすることです。

さらに詳しい情報が必要な場合はお知らせください。サポートが必要な場合は事前にお知らせください。たぶんもっと簡単な方法があります...

こんにちは、CRS Australian Geodetic Datum 66/84 (簡潔にするためにAGD66
) にある多数の経度/緯度座標があります。これらの座標を AGD66 から WGS84 に変更したいのは、それらの間に約 200m の違いがあり、WGS84 に他の座標とレイヤーがあるためです。私はこれをやろうとしました:

このエラーが発生する理由を知っている人、またはこれらの座標を AGD66 から WGS84 に変更する方法について提案がある人はいますか? 事前にご協力いただきありがとうございます。

乾杯、
アダム

座標があり、すべてDCに配置する必要がありますが、それらをNAD 83からRの緯度と経度に変換する方法がわかりませんspTransform()。rgdalパッケージの関数を使用しているため、不適合に関するエラーが発生します。データ。

proj4stringsの他の組み合わせでも、同じエラーが発生します。エラーは座標が大きすぎるためだと思いますが、なぜそうなるのかわかりません。座標のドキュメントは以下のとおりです。

値は、メリーランド州平面メーターNAD83地図投影法にあります。

私はマッピングとプロジェクションに非常に慣れていないので、どんな助けでもありがたいです。

https://github.com/hadley/ggplot2/wiki/plating-polygon-shapefiles http://www.nceas.ucsb.edu/scicomp/に従って、zipレベルのシェープファイルをロードしてプロットを実行しようとしてい ます。ユースケース/ReadWriteESRIShapeFilesなど

私のコード：

このエラーが発生します：

ドライバーを確認しましたが、正直なところ、出力を解釈できません。

file.infoは以下を提供します：

私が試したインターネット検索を考えると、この問題を抱えているのは私だけではないようですが、答えを見つけることができませんでした。問題が.zipフォルダー内にあるシェープファイルに関連しているかどうかはわかりません。仕事用のコンピューターなので、シェープファイルを抽出して.zipフォルダーの外で試すことができるように、ITがWinZipをインストールするのを待つ必要があります。うまくいけば、それは明日になるでしょう。

また、パッケージヘルプには、「データソースディレクトリ（* .dbfなど）内の漂遊ファイルにより、付随する*.shpが欠落しているという重大なエラーが発生する可能性があることに注意してください」と記載されています。.zipフォルダーには次のものがあります：a.dbf、b.prj、c.shp、d.shp.xml、e.shx。

あなたが提供できるどんな助けも素晴らしいでしょう！-アレックス

これらの 3 行は、netcdf ファイルから 1 つの変数を開いて読み取ります。

抽出されたデータを見ると、非常に小さい負の値があることがわかりました。次に、gdalinfo を使用したところ、最大値が 400 で最小値が 0 であることがわかりました (これは正常であり、値は次のようになります)。番号を間違って取得している理由はありますか?

以前にシェープファイルに情報を正常に追加しました（http://ruse​​rgroup.swansea.ac.uk/Healthmap.ashx?HL=mapの私の投稿を参照してください）。

ただし、わずかに異なるシェープファイル（ウェールズの新しいローカルヘルスボード）を使用してもう一度実行しようとすると、コードがspCbind「行名が同一ではないエラー」で失敗します。

以前は私の行の名前が異なっていたので、問題は発生しませんでした。wds.xtraの列のラベルを「NEW_LABEL」に一致するように変更しましたが、それは役に立ちません。

ラベルとラベルの順序は、wales.lonlatとwds.xtraの間で正確に一致します。

（私はR2.13.2上に構築されたRevolutionR 5.0を使用しています）

非常に大きなシェープファイルがいくつかあります。関数SpatialPolygonsDataFrameを使用してそれらを に読み込むことはできますが、ファイルごとに非常に長い時間がかかります。私が実際に興味があるのは、スロットに表示されるのみです。リソースを集中的に使用するポリゴンをスキップして、データだけを読み取る方法はありますか?rgdalreadOGRdata.frame@data

コード例:

SpatialPixelsDataFrame オブジェクトに一連の集水域ポリゴンを (シェープファイルから) オーバーレイするときに奇妙な問題が発生します。70 の流域がありますがNA、60 の流域で s を取得しています。

データは Dropbox で入手できます。次のスクリプトは、ファイルのダウンロード、解凍などを行い、再現可能な例を作成します。

問題はNA次のとおりです。

グリッドが英国地域全体をカバーし、集水域が小さいため、 sはないはずですNA。おそらく、何か間違ったことをしているでしょう! もしそうなら、私は何を間違っていますか?どうすれば修正できますか? の 70 の集水域ごとにDeposition(データ)の平均値を求めたいと思います。SpatialPixelsDataFrameshp

[グリッドはダミー データであることに注意してください。実際の堆積グリッドを配布することはできませんが、データとオブジェクトは基本的に (正確に) と同じdepo.spです。

追加情報

rgdal を使用して R にロードした形状ファイルを強化しようとしていますが、次のエラーが発生します。

maplibrary.org のアフリカ大陸のシェープ ファイルを使用しています。こちらのドロップボックスから入手できます: https://www.dropbox.com/s/etqdw3nky52czv4/Africa%20map.zip

次のコードを使用しています。

そして、前に述べたエラーが発生します。Rにはポリゴンに問題があると思います-これを回避する方法はありますか?