問題タブ [spatial-data-frame]

Spatial*DataFrame(* = ポイント、ライン、ポリゴン、ピクセル、グリッドなど) のデータをどうにかして簡単に操作することは可能ですか? 特に、値を割り当てて操作するのが困難です。

どちらの場合も、かなり醜い回避策を見つけました。

しかし、これらは単なる初心者の試みであり、あまりにも醜くて複雑です.Rでははるかに簡単に違いありません!

私は時空間データを持っています: 36 測定ステーションの 425 日間の毎日の測定。Rでこれらのデータを分析したいので、データからデータフレームを作成しました。これは次のようになります。

ステーション X10004、X10007 のそれぞれについて、緯度/経度の値がありますが、R の利用可能な分析ツールを使用するために、これらの情報をデータ フレームに正しく追加する方法がわかりません。

どうやってするの？または、R の他のデータ構造の可能性を使用する必要がありますか?

ポイント（約47K）と最も近いX国（世界のすべての国）との間の距離を計算しようとしています。ポイントの緯度/経度を SpatialPoints としてインポートし、世界地図を SpatialPolygons としてロードしました。ここで与えられたアドバイスから構築できると思います: SpatialLinesDataFrame: 最小値を計算する方法. 点と線の 間の距離 すべての国とすべての点の間の距離を計算し、最も近い X を抽出する必要があるようです。要するに、ポリゴン制限を課す方法はありますか? そうでない場合は、何を提案しますか-私の唯一の考えは、少数のポイントをインポートしてから、このコードをループすることです(私は新しいRユーザーです)。ありがとう！

さまざまな場所 (X1、X2、...) での温度データの毎日の平均値があり、それらをマップに補間したいと考えています。次のようなフォーマットされた Excel シートからロードして、長いフォーマットのデータ オブジェクトを作成しました。

場所 X1、X2、... の緯度経度についても同じです。

それらを長い形式にマージしました：

と

spplot を使用して、それらを正しい場所にプロットすることができます。また、国境も使用できます。

1 日間の観測に IDW を使用したいのですが、見つけたパッケージ (gstat など) の idw メソッドには、他の「グリッド化された」データ オブジェクトが必要です。そのようなメソッドで補間するために、そのようなデータ オブジェクトを作成するにはどうすればよいですか?

RでrandomForestsパッケージを使用する目的で、ポリゴンデータとラスターデータを1つのデータフレームにマージしたいと思います。
これには、最初にポリゴンごとの平均ラスター値を抽出することが含まれます。

これまでのところ、私は次のものを持っています：

出力は行列で、データフレームに書き込むことができます。しかし、空間座標や元のポリゴン ID がないため、出力を同じデータベースに結合する方法がわかりません。sp=TRUE 引数がこれを行うと思っていましたが、機能していないようです。

RandomForests の目的で実際にポリゴンをポイントに変換する必要があることに注意してください (重心法を使用しますか?)。そのため、ポリゴンではなくポイントに結合された平均ラスター値を結合することが本当に必要であると推測できます。

どんな提案でも大歓迎です。ありがとうございました！！

パッケージ 'maptools' のデータによって提供される世界地図の国境に従ってトリミングしたいラスター グリッドがあります。crop()いくつかの再検索を行うことで、関数を使用してから関数を使用する必要があることがわかりましたmask()が、エラーメッセージが表示されます。

これが私のコードです:

以下のデータフレームがあり、パーセンテージ列のパーセンテージに基づいて、val 列の合計値を Cat 列に戻したいと考えています。

合計行を取得し、グループごとのパーセンテージを入力することもでき、いくつかの列に基づいて2つのデータフレームをマージすることはうまくいきましたが、今は立ち往生しています...

望ましい出力:

または、マージを行うときに値を計算/分散する必要がありますか? 私は助けが必要です...

編集：入力データを提供しました。

ワークブックからこのコードを実行しようとしています

「ミューズ」はSpatialDataFrame

「meuse.grid」はSpatialGridDataFrame

この関数overlayは非推奨なのでover代わりに使用していますが、結果の "meuse.ov" ファイルは S4 クラス データではありません。SpatialDataFrameコードの 2 行目に進むには、"meuse.ov" が である必要があると思います。「meuse.ov」をに変換する方法を誰かが提案できますか、それを生成SpatialDataFrameする代替案を提案できoverますか?

Windows10でR（V3.2.2）を実行しているRStudio（V 0.99.467）で作業しています。rgdalを使用してループでRに読み込み、クラス「SpatialPolygonsDataFrame」のオブジェクトを作成する、ArcGISからのベクターマップのコレクションがあります。

同じループで実行したいデータ QAQC チェックと修正がいくつかあります。たとえば、データ フレームのエントリからスペースを削除したいと考えています。これはループ外で簡単に実行できますが、ループ内での参照と名前の割り当てに苦労しています。ループの外では、このコマンドはデータ内のすべてのスペースを削除します。

この関数はループ内でも機能しますが、スペースのない新しいデータを再割り当てして元のデータを置き換えることはできません。これは、私が何をしようとしているのか、どこで行き詰まっているのかを示すサンプルコードです。

しかし、この最後の行は次のエラーを返します。割り当ての左側に関数を配置できないと思いますか?「get(name)@data[]」と入力すると、R はエラーなしで正しい答えを返すからです。

また、一時データフレームを作成し、一時データフレームを変更し、割り当て関数を使用してマップオブジェクト内のデータフレームを置き換えるさまざまな方法も試しました。ただし、繰り返しますが、これはループ外 (assign(mymap@data[],d)) で機能し、ループ内 (assign(get(name)@data[],d)) では機能しません。

しかし、このメソッドは次のエラーを返します。おそらく、割り当ての最初の引数は、データフレームではなく、単一の要素 (変数名) であると想定されているためでしょうか?

また、変数名を反復処理して列ごとに修正するサブループも作成しましたが、メイン ループ内の 1 行ですべての修正を実行できる場合、これは非常に非効率的です。上に貼り付けた2つの例のいずれかを使用するための解決策を誰かが持っていることを願っています。

StackOverflow 内のクラス SpatialPolygonsDataFrame のオブジェクトの生データの例を添付する方法がわかりません。ここに十分な情報がない場合は、GitHub にパブリック フォルダーを作成してリンクを提供できますか? しかし、私はまだRとプログラミングに慣れていないので、これはかなり基本的な割り当ての問題だと思います.