3

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

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

horiz=rep(1:5, 5)
depth=runif(25)
value = horiz+runif(25)/5
df <- data.frame(horiz, depth, value)

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

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

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

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

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

4

1 に答える 1

4

ですから、答えてほしい質問がいくつかあります。

  1. このspTransform関数は SPDF を作成しませんが、投影間で変換します。SPDF を作成するには、最初に単純なものを使用できますdata.frame。dfSPDFに変換するには:

    coordinates(df) = c("horiz", "depth")

  2. の OS X バイナリはrgdalhttp://www.kyngchaos.com にあります。しかし、rgdal が必要かどうかは疑問です。

  3. spTransformlatlong データだけでなく、投影データも操作できます。しかし、rgdal または は必要ないと思いますspTransform。ポイント 1 も参照してください。

  4. ポイント1を使用してSPDFを作成したら、言及した投稿の情報を使用して続行できます。

于 2012-06-13T15:47:33.323 に答える