タグ付けされた質問 「kriging」

理想的にはオープンソースソフトウェアパッケージを使用して、Krigeにしたいポイントデータセットがあります。可能であれば、プロセス中にセミバリオグラムモデルを選択して、推定を改善したいと思います。

42 open-source-gis  interpolation  spatial-statistics  kriging 

私は、約10m離れた非常に規則的な間隔のポイントデータを使用してDEMを作成しようとしています。私が内挿しているエリアは、多くの平坦な駐車場とサッカー場がある教育施設ですが、駐車場に停滞することが多いかなり急な丘がまだあります。これらの既知のプラトーのため、スプライン法は除外しました。ただし、IDWとKrigingのどちらの方法を使用するかはまだわかりません。両方を試してみても大きな違いは見られませんが、少し調査してもまだ決定していません。 誰かがこれを私のために解決するための知恵の言葉を持っていますか？

23 dem  interpolation  kriging  inverse-distance-weighted 

私は、クリギング法を使用して補間するためのリクエストとともに、いくつかの標本でいくつかのデータを取得します。 調査の結果、クリギングの結果（ArcGIS Geostatistical Analystでデフォルトのパラメーターを使用して実行）は満足のいくものではないようでした。補間された値は、測定値（特に一番上の値）とは大きく異なり、表面は信頼できるようには見えません。写真は次のとおり です。主な問題はサンプル数が不十分だと思います。 信頼できる結果を得るには、何ポイント使用する必要がありますか？ それとも、クリギング法はそのような多様な価値に適していませんか？

19 arcgis-10.1  interpolation  kriging  geostatistical-analyst  sampling 

私はGISデータの使用に非常に慣れており、Rをあまり経験していません。spatial-analyst.netPDFブックを使用して空間データを分析する方法について読んでいるので、完全に失われたわけではありませんが、私の問題と人々はアイデアを提案するかもしれません。 さまざまな緯度/経度座標で約2000の測定値を含むデータセットがありますが、データが3年間にわたって収集され、条件が時間とともに変化したため、おそらくこのデータセットを細分化します。測定対象の変数を「IP」と呼びましょう。 サンプルデータに対してKrigingまたは他の補間方法を使用して、問題のある領域全体にIPのマップを作成したい。次に、さまざまなIPバケットの土地の量を測定するヒストグラムを作成します。また、各バケット内のサンプル数を示すヒストグラムを作成する必要があります（サンプルは、クリギングがその土地に対して予測するものよりも実際のIPが大きい場合も小さい場合もあります）。 データをSpatialPointsDataFrameに読み込み、クリギング分析を実行する方法に従いますが、その問題を抱えているのは、ヒストグラム分析を行うためにそのデータをグリッドデータフレームに変換する方法です。 ポイントをグリッドに変換するための提案はありますか？

14 raster  r  kriging 

Rで通常のクリギングを使用して地下水標高マップを作成しようとしていますが、常にメモリ不足になります サイズ146.5 Mbのベクトルを割り当てることができません 割り当て（windows xp 32bit）。メモリサイズを増やすべきか、どのようにすればいいのか、グリッドサイズを減らすにはどうすればいいのか（Rの初心者のように）？ この例http://casoilresource.lawr.ucdavis.edu/drupal/node/438に従ってください

11 grass  r  error  kriging 

カテゴリカルデータをグリッド化する方法を探しています。海図とフィールドシートから、海底の性質を特定するポイントのコレクションを抽出しました。これらのデータは、数値ではなくカテゴリであり、定期的またはランダムにサンプリングされることもありません。航海図は、ナビゲーションとアンカーを支援するために作成されます。それらは生息地をマッピングするために作成されていません。そのため、比較的浅い水深が航行に危険をもたらす可能性があり、船が停泊する傾向がある海岸近くで、より多くの探査が行われます。岸からさらに遠く、航行に十分な深さがあり、アンカーが実用的でない場合、探査はそれほど頻繁に行われません。 他の誰かが海図から格子状の基質マップを作成しようとしましたか？ 私はティーセン（ボルノイ）ポリゴンを調べましたが、海岸に沿った音響の集中は、海岸沿いの細かい「ハニカム」、オフショア、オフショアに伸びる長いパイ形ポリゴンの間にある大きなポリゴンにつながります。最近傍を使用したグリッドは、ほぼ同じ結果になります。 岸に近い浅いポイントの影響を制限する方法が必要です。長いパイ形のポリゴンを制限する方法です。深い海では、海底の性質が岸に近い底部の続きになるとは思わない。私は2つの線に沿って考え始めました-両方とも深さを使用しています。1つは、グリッドセルと隣接する点の間の深さの違いを使用して、「最も近い」隣接の選択に重みを付けることです。もう1つは、特定の許容差よりも深さが異なる隣接ポイントの選択を解除することです。または、事前に指定された許容値ではなく、深度範囲をビン化して、隣接するポイントの選択を同じ深度範囲またはビン内のポイントに制限することもできます。 これらの2つのオプションのいずれかを実装する方法についての考えはありますか？ 他のフォーラムで同僚と話し合って以来、私は他にもいくつかのアプローチを検討してきました。1つ目は、岸壁のデータの影響を制限するために、バリア（深度100mのコンター）を使用することです。このアプローチの課題は、バリアを使用できるすべてのESRI補間ルーチンが、不連続データではなく連続データを処理するように設計されていることです。ティーセンポリゴンを作成する前に、バリアを使用してポイントを岸の近くの浅いポイントとより深いポイントに分割できます。ただし、ArcGISは複雑な領域ではなく長方形の領域に対してティエセンポリゴンを作成するため、エッジ効果が蔓延すると予想しています。 2番目のアプローチ-複数の同僚によって提案された-はクリギングでした。私は、継続的なデータのためにそれを検討したことがないので、最初は手に負えないクリギングを却下しました。クリギングの課題は、クリギングもカテゴリカルデータ用に設計されていないことです。今、私は表面の深さと性質を使ったcokrigingを見ていますが、どのようなタイプのクリギングでも、表面の性質に整数の数値コードを使用する必要があります。その後、結果の浮動小数点数値コードは、元の整数コーディングに戻す必要があります。きれいじゃない。 他の方針に従うことを提案できますか？（おそらく、地形分析を使用できます。たとえば、安息角より急な斜面は堆積物にはなりません。もっと単純なものを探しています。とにかく、十分な空間解像度のデータがありません。） よろしく、

10 raster  interpolation  kriging  voronoi-thiessen 

私は、湾内に流れ込む半島の反対側のいくつかの川のデータポイントを使用して、沿岸地域で補間を試みています。 海岸線レイヤーを作成してバリアとして使用しましたが、ArcGIS 10でこれを行う方法がわかりません。 クリギングツールのどこにもバリアオブジェクトが表示されません。 クリギングを行い、それを土地にクリッピングするだけで試しましたが、データポイントが互いに接近しているが、土地によって分離されていると問題が発生するという非現実的な結果が得られます。 9.3のジオプロセシングスクリプトでこれを行うためのオンラインドキュメントを参照してください。

10 arcgis-desktop  arcgis-10.0  interpolation  kriging  barriers 

データテーブルで使用する値に応じて、クリギングを実行したい場合にのみ機能します。krige関数の結果として、私はのために取得var1.pred: NA NA NA ...し、var1.var: NA NA NA ...（私は私のDataTableの「間違った」値を使用する場合のみです。） 例えば： それが動作します常に（今のところ）私は10個の値を使用する場合 50個の値を使用すると機能しますが、特定の値でのみ機能します 50個の値と「間違った」値を使用すると機能しません 25個の値と前述の「間違った」値を使用すると機能します なぜ機能しているのか、機能していないのかわかりません。奇妙なことはZwiesel;49.02999878;13.22999954;2.2、データテーブルに追加するときに、20未満の値を使用すると機能するが、50を超える値を使用すると機能しないということです... 私の間違いはどこですか？ myWeatherTable.csv： Place;Latitude;Longitude;Temperature Aachen;50.77999878;6.09999990;3 Abbikenhausen;53.52999878;8.00000000;7.9 Adelbach;49.04000092;9.76000023;3.1 Adendorf;51.61999893;11.69999981;1.9 Alberzell;48.45999908;11.34000015;4.6 ... ... クリギング補間を実行するための私のコード WeatherData <- read.csv(file="myWeatherTable", header = TRUE, sep ";") coordinates(WeatherData) = ~Longitude + Latitude vario <- variogram(log(Temperature) ~1, WeatherData) vario.fit <- fit.variogram(vario, vgm("Sph")) min_lon <- …

10 r  interpolation  kriging  gstat 

のドキュメントで、関数spのcoords引数SpatialPolygonsは2列しか受け入れないことがわかりました。3次元ポリゴンを作成することはできませんか？のgstat3番目の空間次元として時間を追加することにより、パッケージを使用して時空間クリギングを行っていSpatialPointsDataFrameます。これを拡張して、ポリゴンをブロックとするクリギングをブロックします。それを行う方法はありますか？

9 r  interpolation  spatial-statistics  kriging  gstat 

エリアをカバーするために六角ビニンググリッドとして作成したテッセレーションされた六角形のテーブルがあります。 次に、いくつかの六角形に値を与えるためにポイントセットをオーバーレイしました。図を参照してください： 値が「0」のシェイプには値がありません。クリギングに従って、これらを隣接の値から補間したいと思います。私のテーブルが次のようになっているとしたら、どうすればpostgresでこれを達成できますか？ hex_grid_data ( gid SERIAL, wkb_geometry (geometry), value INTEGER )

8 postgis  interpolation  kriging  hexagon-binning