問題タブ [automap]

私は ggplot に比較的慣れていないので、私の問題のいくつかが本当に単純であるか、まったく解決できない場合は、ご容赦ください。

私がやろうとしているのは、形状の塗りつぶしが連続している国の「ヒート マップ」を生成することです。さらに、国の形をしてい.RDataます。Hadley wickham のスクリプトを使用して、SpatialPolygon データをデータ フレームに変換しました。データ フレームの long と lat のデータは次のようになります

この緯度経度データは、ドイツの輪郭を描きます。残りのデータ フレームは必要ないと思うのでここでは省略します。特定の経度/緯度点の値の 2 番目のデータ フレームもあります。これはこのように見えます

私が今やりたいことは、そのポイントから特定の距離内にあるすべての固定ポイントの平均値に従って、マップの各ポイントに色を付けることです。そうすれば、国の地図全体を (ほぼ) 連続して塗りつぶすことができます。私がこれまでに持っているのは、ggplot2 でプロットされた国の地図です。

私の最初のアイデアは、描画されたマップ内にあるポイントを生成し、生成されたすべてのポイントの値を次のmy_generated_pointように計算することでした

ただし、これらのポイントを生成する方法は見つかりませんでした。そして、問題全体と同様に、この方法で解決できるかどうかさえわかりません。私の質問は、これらのポイントを生成する方法が存在するかどうか、および/または解決策に到達する別の方法があるかどうかです。

解決

ポールのおかげで、私はほとんど欲しいものを手に入れることができました。オランダのサンプル データの例を次に示します。

米国チェサピーク湾のさまざまな場所で採取された種の豊富さの調査データがあり、そのデータを「ヒート マップ」としてグラフィカルに提示したいと考えています。

緯度/経度座標と豊かさの値のデータフレームがあり、これを に変換し、automap パッケージ SpatialPointsDataFrameの関数を使用して補間値を生成しました。autoKrige()

まず、autoKrige()関数を正しく実装しているかどうかについてコメントできますか?

次に、データをプロットして地域の地図を重ねるのに問題があります。または、ベイの境界を反映するように補間グリッドを指定できますか (ここで提案されているように)? どうすればそれができるか、どこでその情報を入手できるかについて何か考えはありますか? グリッドを提供するのautoKrige()は簡単に見えます。

編集: 非常に役立つ投稿をしてくれた Paul に感謝します! これが私が今持っているものです。ggplot が補間されたデータと地図投影の両方を受け入れるのに問題があります。

クリギング関数を使用して海洋トランセクトの化学パラメータの垂直マップを作成しようとしていますが、なかなか始められません。

私のデータは次のようになります。

パッケージ内のautoKrige関数は、automap私のために仕事をする必要があるように見えますが、 class のオブジェクトを取りますSpatialPointsDataFrame。私が知る限りspTransform、パッケージ内の関数はオブジェクトrgdalを作成しSpatialPointsDataFrameますが、2 つの問題があります。

• これの OSX バイナリはCRAN から入手できず、OXS 10.7 で実行されている RStudio の私のコピーはそれをインストールできないようです。

• この関数は、緯度/経度データと地球の曲率の正しい距離値で機能するようです。私は垂直面（および短い距離、数百メートルのスケール）を扱っているので、距離を修正したくありません。

R でのクリギングに関する優れた議論がここにありますが、上記の問題により、特定の問題にそれを適用する方法がよくわかりません。

化学パラメーターの補間値を含む点のグリッドを記述するマトリックスまたはデータフレームが必要です。これをプロットできます (理想的には ggplot2 を使用)。私の問題の解決策は、私が考えているよりもかなり簡単だと思います。

パッケージのR関数を使用していautokrigeますautomapが、エラーが発生し、解決方法がわかりません。ヒントはありますか？

ありがとうございました！

補間を作成するためにgstatライブラリを使用するコードをRで記述しようとしています。私はすでにgstatマニュアルを読み、インターネット上のいくつかの例に基づいて、このコードを書くことができました（これはほんの一部です）：

私の問題は、通常のクリギング補間で結果を取得する代わりに、最後のコマンド（predict）を実行すると、逆距離加重（IDW）で結果を取得することです。gstatのマニュアルで、「バリオグラムが指定されていない場合、逆距離加重補間がデフォルトのアクションです。バリオグラムが指定されている場合、デフォルトの予測方法は通常のクリギングです」と読みました。

しかし、私のコードでわかるように、経験的バリオグラムと理論的バリオグラムの両方を指定します。なぜ私が通常のクリギングの代わりにIDWを取得し続けるのか知っていますか？私が持っている座標のタイプと関連付けることができますか？たとえば、座標が互いに近い場合、または関心領域が大きすぎる場合はどうなりますか？どんな助けでも本当に役に立ちます。

よろしくお願いしますディミトリス

私は R の初心者です。autofitVariogramを 50 ステーションの毎日の降水量データに使用しています。サンプル データを以下に示します。一部のステーションには、"NaN" 値で表される欠損値があります。

私の質問は、バリオグラムフィットに関するものです。バリオグラムは 60,000m の距離しかカバーしていません。60Km を超えるビンのポイントがプロットされないのはなぜですか。空間相関プロットから、経度-緯度情報からの最大距離が 200Km を超えていることがわかりました。

緯度と経度の情報を以下にまとめます。summary(lonlat) lon lat
Min. ：74.78分 :15.77
1st Qu.:75.14 1st Qu.:16.04
中央値:75.56 中央値:16.33
平均値:75.54 平均値:16.37
3rd Qu.:75.94 3rd
Qu.:16.66 :76.31 最大 ：17.23

Q2) バリオグラム フィットの png ファイルをループで保存するにはどうすればよいですか。私が行った図を保存するたびに dev.off() を使用する必要があることは理解していますが、図を保存することはできません。どんな助けでも大歓迎です。

ありがとう、

何か提案をいただければ幸いです。

geoRパッケージのca20-Datasetの相互検証を行いたい。たとえば、meuse-datasetの場合、これは正常に機能しますが、このデータセットの場合、SpatialPointsDataFrameのディメンションで奇妙な問題が発生します。たぶん、これを自分で試して、autoKrige.cv関数が機能しない理由を説明することができます（いくつかのnfold-valuesを試しましたが、これはエラーメッセージのlocations-valueのみを変更します...）：

誰かが問題を再現できることを願っています。私のRバージョンは2.15で、すべてのパッケージが最新です（少なくとも1か月ほど前ではありません...）。

ご協力いただきありがとうございます！！

autokrige を使用して降雨データを補間したい automap ライブラリを使用しています。結果として得られる補間されたマップには、予測ポイント マップとは異なる次元 (行と列の数が異なる) があります。これはどのように可能ですか？以下の私のコードを見てください：

autoKrig 関数に問題があり、ここで再現可能な例を作成しようとしています。

これにより、次のようになります。

元のデータでも同じエラーが発生します。助けていただければ幸いです。