問題タブ [rgdal]

特定の色のマトリックスに基づいて GTiff を作成するにはどうすればよいのでしょうか? これは、色が量的変数を表すのではなく、人口統計クラスターなどの質的空間現象を表す場合に発生する可能性があります。たとえば、色値の 1000x1000 マトリックス ('#FF0000' のような文字列) がある場合、これらを簡単にプロットできますgrid.raster(mtrx, interpolate=FALSE)が、これを他のシェープファイルと重ね合わせるのは困難です。また、フィールド名が「x」、「y」、「col」の SpatialPixelsDataFrame に変換することもできます。しかし、これらのオブジェクトのいずれかを GTiff に書き込む方法に困惑しています。例えば：

収量

ポインタに非常に感謝します。

例：

R でシェープファイルを開く 2 つの基本的な方法を見つけました - と を使用rgdalしmaptoolsます。

どちらの場合も、データ構造はまったく同じように見えます (クラス SpatialLinesDataFrame、パッケージ sp)。ただし、rgdalプロジェクションを適切に読み取りますが、そうでmaptoolsはありません(CRS を手動で割り当てる必要がある可能性があります)。

では、なぜmaptoolsシェイプ ファイルを開くために を使用するのでしょうか? CRS を手動で割り当てるのを間違えるだけかもしれません。

どちらの方法も同等であると結論付けてもrgdalよろしいですか?

私はこの質問で尋ねられたのと同じことをしようとしています. Cartogram + choropleth map in Rですが、 SpatialPolygonsDataFrame から始めて、同じタイプのオブジェクトになることを望んでいます。

オブジェクトをシェープファイルとして保存し、scapetoadを使用して再度開き、元に戻すこともできますが、手順を完全に再現できるように、また数十のバリエーションを自動的にコーディングできるように、すべてを R 内に置きたいと思います。

私は github で Rcartogram コードをフォークし、これまでの努力をここに追加しました。

基本的に、このデモが行うことは、マップ上に SpatialGrid を作成し、グリッドの各ポイントで人口密度を検索し、これをcartogram()が作業するために必要な形式の密度マトリックスに変換することです。ここまでは順調ですね。

しかし、の出力に基づいて元のマップ ポイントを補間する方法はcartogram()?

ここには 2 つの問題があります。1 つ目は、マップとグリッドを同じ単位にして、補間を可能にすることです。2 つ目は、すべてのポリゴンのすべてのポイントにアクセスして補間し、それらをすべて正しい順序に保つことです。

グリッドはグリッド単位で、マップは投影単位です (例の longlat の場合)。グリッドを longlat に投影するか、マップをグリッド単位に投影する必要があります。spTransform()私の考えでは、偽の CRS を作成し、これを の関数と一緒に使用します。これpackage(rgdal)は、最小限の手間でオブジェクト内のすべてのポイントを処理するためです。

SpPDFオブジェクトにいくつかのレイヤーがあるため、すべてのポイントにアクセスすることは困難です。オブジェクト>ポリゴン>ポリゴン>ライン>座標だと思います。マップ全体の構造をそのまま維持しながら、これらにアクセスする方法はありますか?

これを行う方法のかなり単純な例を見つけましたが、うまくいきません。私はRにかなり慣れていません

しかし、これは例の数字です。変換する必要がある何千もの座標があり、それらをテーブルからこのスクリプトに取得する方法がわかりません

私のデータセットには ID、X、Y の 3 つの列があります。この式を使用してそれらを変換するにはどうすればよいですか? 私は何週間もこれにこだわっています

EPSG 2790 ( http://spatialreference.org/ref/epsg/2790/ ) またはおそらく ESRI 102672 ( http://spatialreference.org/ref/esri/102672/ )。

これは間違いなく以前に尋ねられました。私のコードは、ここの回答に基づいています ( rgdal R パッケージの spTransform の「非有限変換が検出されました」およびhttp://r-sig-geo.2731867.n2.nabble.com/Converting-State-Plane-Coordinates-td5457204. html )。

しかし、何らかの理由で私はそれを動作させることができません:

与えます:

IPython Notebook を使用して、Python と R コードを 1 つのワークフローに統合し、いくつかの空間分析を生成しています。以前は、このタスクを実行するために必要なインフラストラクチャを両方の名前空間にインポートできていましたが、その後の更新により何かが変わったと思います。特に、 %%Rセルでrgdalライブラリを呼び出すことができなくなりました。

これが奇妙である理由は、ベース R と RStudio でライブラリを呼び出すことができるためです (出力は同じです)。

RStudio と Notebook のバージョン情報は同じです...:

...Rの場所は次のとおりです。

Notebook 内からrgdalをインストールしようとしたことに注意してください...:

...そして、RStudio に正常にロードされたものと同じバージョンです。

それで、私はここで何が起こっているのか途方に暮れています。RStudio で作業を進めることはできますが、作成しようとしている製品が断片化されています。この問題に関するアドバイスをいただければ幸いです。

補遺:これは %rmagic だけの問題ではないようです。Python で別の操作を実行しようとして、osgeoとgdalの両方をインポートしようとすると、次のエラーが発生しました。

文字通りGDALを再度インストールしたので、最新バージョンを使用しています。考え？

readOGR を使用してポリゴン シェープファイルを正常にインポートしましたが、.prj ファイルの情報がインポートされません。

テスト クラス: SpatialPolygonsDataFrame nfeatures: 19407 範囲: 35551.4、1585917、6318047、9408727 (xmin、xmax、ymin、ymax) 座標。参照。: なし

summary(test) クラス SpatialPolygonsDataFrame のオブジェクト 座標: 最小 最大 x 35551.4 1585917 y 6318047.3 9408727 射影: NA proj4string : [NA]

射影は utm21n 複合ゾーンである必要があります

R でラスター ファイルの上にシェープファイルをプロットしたいのですが、それらを完全にオーバーラップさせることはできません。ラスターが反時計回りに数度回転しているように見えます。投影の問題でしょうか。

次の MWE を考慮してください