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

周囲のポイントの既知の値に基づく非サンプリングポイントの値の推定。

2
ArcGIS IDW空間内挿を使用する場合、ポイントデータは等距離で投影されますか?
ArcGIS IDW空間内挿を使用するときに、ポイントデータを等距離投影に投影する必要があるかどうかを誰かが明確にできるかどうか疑問に思っていますか? 私は、北アメリカ西部の緯度約30度にわたるデータセットで作業しています。データは現在、緯度/経度(NAD83)です。 ArcGISは「オンザフライで投影」するか、補間手順中にサンプルポイントまでの距離を計算するときに緯度を調整するのですか、それとも距離を維持する投影ですべてを提供する必要がありますか?

1
ナゲット効果の高いセミバリオグラムを解釈しますか?
パッケージgstat、variogram()関数を使用して、Rでセミバリオグラムを作成しました。モデルの残差に空間的自己相関があるかどうかを確認したい(glmmを使用して、数kmから900 km間隔で離れたサイト間で、生息地の関数としての種の存在量)。 私の単位はkmです。したがって、私の解釈では、空間的自己相関が「問題」でなくなるまで、範囲は100 kmを少し超えていると解釈されます。ナゲットがこんなに高いように見える理由を誰かが説明できるかどうか疑問に思っていますか?これは、同じような場所でも比較的大きな差があることを意味しますか?または、この波状のバリオグラムは、より一般的な形状になるまでラグ数とラグ距離を調整する必要があることを意味しますか? もう少し詳しく調べるために、私variog()はパッケージgeoR の関数も使用し、を使用breaks=seq(0,100,10)して、より近い距離を(同じ点と同じモデルの残差を使用して)調べようとしました。これは、最も近い点がより異なっていることを示していますが、これも意味がありません。多分これは空間的自己相関がないことを示しており、私のモデルはすでにこれを説明しています。 私はこの優れた情報源「涙のないジオスタット」を見つけました。51ページで、バリオグラムのフィッティングについていくつかの良いアドバイスがあります。このアドバイスにより、私の最初のものは正しい範囲を持っているようです。これは最初の質問に戻ります-これをどのように解釈しますか?

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

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

5
サンプルポイントから土壌深度を補間するための最良の技術は何ですか
これが私の最初の投稿です。 土壌深さ調査からいくつかのポイントがあり、領域全体の深さを補間する必要があります。これに最も適したテクニックは何ですか?また、このテクニックを使用するときに考慮する必要があることはありますか?チュートリアルへのポインタは素晴らしいでしょう! ありがとう

1
R gstat krige()-位置[5.88,47.4,0]で特異な共分散行列:スキップ
データテーブルで使用する値に応じて、クリギングを実行したい場合にのみ機能します。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 <- …

1
QGISを使用して次数データを補間しますか?
風マップを作成したい。ArcGISおよびQGISの紹介があります。Spatial Analystにアクセスできないため、QGISバージョンを試してみました。ただし、次数データを補間する方法がわかりません。 たとえば、風向は1つのポイントで359°、2番目のポイントで1°です。このポイント間の値は、おそらく359°、0°、または1°になるはずです。 補間はそれが次数データであることを認識していないため、補間された値は約180°のように見えます。 QGISを使用して学位データをどのように補間できますか?

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

1
QGISで等高線から標高を正確に補間しますか?
QGISの等高線(1m間隔で〜70000、〜20km ^ 2の範囲)から標高を補間しようとしています。 プラグインで補間を使用Inverse Distance Weighting (IDW)するInterpolationと、隣接する輪郭だけでなく、すべての輪郭線が使用されますか? 使用するGDAL Rasterプラグイン(パワーとは逆の距離)は、すべての輪郭を使用していないラインとしてではなくポイントとして等高線を扱うように見えますか? またTriangular interpolation (TIN)、等高線を構造線として扱ってみましたが、おそらく等高線の数と密度が原因でQGISがクラッシュしました。 5未満の累乗の場合、GDALラスタープラグインを使用すると、結果のDEMに非常にノイズが多く、輪郭点が極大または極小になる傾向があります。 つまり、等高線が走っているが等高線を定義するポイントが含まれていないグリッドセルの場合、補間された標高は、等高線の標高が割り当てられるのではなく、すべての等高線データの平均標高に向かってシフトされます。それらを通して; これは、補間がすべての単一のコンターポイントの標高を不必要に考慮しているためです(ただし、逆距離で重み付けされます)。 この問題は、より高い累乗でIDWを使用すると軽減されますが、比較的平坦な領域は、小さな傾斜のある大きな領域のようではなく、ステップのように最近傍補間のように見え始めます。 他の投稿から、GRASS r.surf.contourアルゴリズムにも同様の問題があるようです。 検索半径を使用することをお勧めしますが、検索半径よりも小さいフィーチャを平滑化し、グリッドセルに関連しない等高線が補間された標高に影響を与えることを許可します。 QGISで、全方向に直接隣接する等高線のみを使用して補間する方法はありますか? QGISの等高線データから補間するためのベストプラクティスは何ですか?等高線データは、実際には1つの大きな点または線のセットとしてではなく、等高線として扱われますか?

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


1
非サンプリングポイントの値の推定
沿岸干潟のバイオマスを測定したい。ポリゴン内のポイントにしかアクセスできません。ポリゴン内部のポイントの値に基づいて、ポリゴン外部のポイントの値を推定できる方法はありますか? set.seed(5) x <- rnorm(50, -1.841, 0.01) set.seed(50) y <- rnorm(50, 55.663, 0.01) xy <- data.frame(x,y, values=rnorm(50)) coordinates(xy) <- c("x", "y") proj4string(xy) <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84") plot(xy) makePolygons <- function(coordsx, coordsy){ coords <- matrix(c(c(coordsx, coordsy)), ncol=2) p <- Polygon(coords) p <- Polygons(list(p), ID = "p") myPoly <- SpatialPolygons(list(p)) spdf = …

1
ラスター内挿法の選択方法は?
私は、約130万の住宅建物のデータセットを使用しています。それらのそれぞれは、0から100の範囲の正規分布変数によって記述されます。データセットはスイス全体をカバーしているため、非常に密な点密度と非常に疎な点密度の両方の領域があります。 データプロバイダーによる制限は、5つの建物よりも詳細なデータセットをリリースできないことです。 このポイントデータをラスターサーフェスに変換することを考えていました。この演習の主な目的は、誰かが対象のポイント位置のデータを抽出することにより、このラスターからポイント値を取得しようとしたときに発生するエラーを最小限に抑えるサーフェスを作成することです(元のデータセットからデータを抽出するのと比較して)ポイントの場所)。 そのためにはどのような方法が最も適しているのか知りたいのですが。Spatial Analystツールボックスから選択するものはほとんどなく、さらにGeostatisticalツールボックスから選択するものもあるので、どこから始めればよいかを非常に感謝します。

3
補間されたOSMポイントを使用した最短経路
OSMネットワークにスナップした一連のGPSポイントがあります。下のスクリーンショットでは、GPSポイントは赤、スナップされたポイントは緑です。 これらすべての緑のウェイポイントを含む最短経路を計算したいと思います。私の解決策は、ポイントの各ペア間の最短経路を計算し、最終的に結果を連結することです。 私の問題は、dijkstra_spがOSMネットワーク上の任意のポイントを受け入れないことです。スナップされたポイントは、次のロジックを使用して計算されたため、必ずしもウェイテーブルにあるわけではありません。 特定のGPSポイントに最も近い方法を見つけます。 補間を使用して、この方法でGPSポイントに最も近いポイントを見つけます。 スナップされたポイントは補間によって導出されたため、ウェイテーブルにはありません。 だから私の質問は次のとおりです:ウェイテーブルに必ずしもないOSMネットワーク上の2点間の最短経路を計算するにはどうすればよいですか?

1
ラスター(DEM)を部分的に更新していますか?
ArcView 9.3で.img DEMを更新する方法 部品が少ない場合にのみ更新する必要があるDEMがあります。更新される領域は、バリア付きのスプライン補間を使用した別のDEMで作成されています。 更新される領域自体はS字型ですが、より大きな長方形のDEMを受け取りました(ところで、バリアの外側で-9999の値を受け取りたいと思っていましたが、外側の領域も補間されているようです)。 私の目的は、ポリライン/バリア内の領域のみを更新することです。どうやってするか?
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.