Kriging z barierami w ArcGIS Desktop?

10

Próbuję dokonać interpolacji w strefie przybrzeżnej z punktami danych w kilku rzekach opróżniających się do zatoki i po przeciwnych stronach półwyspów.

Utworzyłem warstwę wybrzeża, aby wykorzystać ją jako barierę, ale nie mogę wymyślić, jak to zrobić w ArcGIS 10.

Nigdzie nie widzę obiektu barier w żadnym z narzędzi kriging.

Próbowałem po prostu wykonać kriging, a następnie przyciąć go do ziemi, ale otrzymuję nierealistyczne wyniki, w których punkty danych blisko siebie, ale oddzielone przez ziemię, stwarzają problemy.

Widzę dokumentację online dotyczącą robienia tego w skryptach geoprzetwarzania w 9.3, ale nic podobnego dla 10.

arcgis-desktop arcgis-10.0 interpolation kriging barriers elayne927
źródło
Aby dodać linki do podglądu zmian, porównaj dokumentację ArcGIS 9.3 i ArcGIS 10
Mike T

Odpowiedzi:

3

Czy rzeki płyną ogólnie w określonym kierunku? Jeśli tak, możesz użyć binowania kierunkowego wzdłuż linii środkowych rzek, aby ograniczyć empiryczne semiwariogramy tylko do par punktowych biegnących wzdłuż rzek. (Aby uzyskać więcej informacji, zobacz tutaj: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Exploring_directional_autocorrelation/0031000000m2000000/ )

Inną opcją byłoby podzielenie próbek punktów na osobne zestawy danych i uruchomienie kriging na każdym zestawie osobno. To i tak może być bardziej odpowiednie, ponieważ każda rzeka może mieć dla nich inne modele statystyczne; ale nie zadziałałoby, jeśli próbka jest mała dla każdej rzeki.

Czy potrzebujesz metody geostatystycznej, czy może zadziałałaby metoda deterministyczna? Myślę, że możesz osiągnąć to, co chcesz, używając interpolacji dyfuzyjnej z barierami. Możesz wykorzystać swoją ziemię jako niezwykle wysoką barierę kosztową. http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//00310000002s000000.htm

Również tutaj jest interpolacja dyfuzyjna z absolutnymi barierami jako narzędzie skryptowe do geoprzetwarzania. Twoje wybrzeże stanowiłoby idealną absolutną barierę. http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Diffusion_Interpolation_With_Barriers/003000000005000000/

blord-castillo
źródło
0

Czy próbowałeś narzędzia Topo to Raster? Możesz określić typ obiektów wejściowych (Wysokość punktu, Strumień, zlew, kontur, granica, jezioro). Możesz spróbować ustawić funkcje warstwy linii brzegowej na „TopoLake”

Ponadto wynikiem działania narzędzi Topo to Raster jest raster „poprawny hydrologicznie”

Cytując z pomocy: TopoLake Klasa cech wielokąta, która określa lokalizację jezior. Wszystkie wyjściowe komórki rastrowe w jeziorze zostaną przypisane do minimalnej wartości wysokości wszystkich komórek wzdłuż linii brzegowej.

Jakub Sisak GeoGraphics
źródło
Uważam, że zamiast budowania reprezentacji topografii, OP chce stworzyć powierzchnię pomiarów, takich jak jakość wody, dokonywanych w jednolitych częściach wód. Topo2Raster tego nie zrobi. (To też nie Krige.)
whuber
Tak, staram się to dokładnie zrobić.
elayne927
Masz całkowitą rację Bill, Topo do Raster nie pozwala na kriging. Miałem wrażenie, że Elayne odnosi się ogólnie do interpolacji. (Często słyszę, jak ludzie
zbijają
Elayne, myślę, że wiem, co rozumiesz przez nierealne wyniki. Możesz spróbować przekonwertować jeziora na siatkę regulare rozmieszczonych punktów z wartością elewacji jeziora przypisaną do każdego punktu i dodać go jako dane wejściowe do narzędzia Topo do Raster jako „Elewacja punktu”. To również nie zapewni idealnego wyniku, szczególnie wzdłuż linii brzegowej ale będą wyglądać bardziej realistycznie. Czy próbowałeś już numeru NIP i linii nieciągłości, a następnie przekształcałeś wynikowy NIP w raster? Czy możesz zamieścić link do skryptów geoprzetwarzania, o których mowa?
Jakub Sisak GeoGraphics
0

Jestem spóźnionym adoptującym, więc mogę tylko powiedzieć, co działa w wersji 9.3 z rozszerzeniem Spatial Analyst - publikowanie na wypadek, gdyby metoda nadal działała w wersji 10. Najpierw utwórz zamknięty znacznik liniowy reprezentujący obszar zainteresowania (brzmi jak odwrotny twojej linii brzegowej tutaj). Określ „kontener” w oknie dialogowym opcji Spatial Analyst na karcie zasięgu. Przeprowadź analizę, wybierając Kriging z Interpolacji do Rastra w menu Analizatora przestrzennego.

Powodzenia!

rec. thegeom
źródło