Czy tworzysz wielokąty wokół zasięgu podobnych punktów za pomocą ArcGIS Desktop?


11

Mam miliony punktów zebranych za warstwę litologiczną.

Użyli kodowania dla różnych rodzajów skał.

Muszę stworzyć wielokąt wokół podobnych punktów.

Jaki jest najprostszy sposób na uzyskanie wielokąta zamiast ręcznej digitalizacji.

Szukałem narzędzi, które konwertują punkty na wielokąty, ale wygląda na to, że ich nie ma.

Widziałem narzędzia do przekształcania linii w wielokąty, wielokąty w linie i punkty, ale nie do punktów w wielokąty.


4
Zobacz komentarz moderatora do poprzedniego pytania . Gdy przestrzegasz zasad, ludzie doceniają to i częściej oferują szybkie odpowiedzi, gdy naprawdę ich potrzebujesz.
whuber

Czy twój problem „muszę stworzyć wielokąt w zakresie podobnych punktów”? Wypróbuj prostokąt (punkt (min (x), min (y)), punkt (maks. (X), maks. (Y))) ... lub bardziej konkretnie.
Remigijus Pankevičius

w jakim środowisku jesteś
Ragi Yaser Burhum

Wygląda na to, że chcesz wypukłej funkcji kadłuba. Jeśli tego właśnie szukasz, to pytanie zostało już rozwiązane.
grad

Spróbuj „zrobić jeden wielokąt z punktów” „XToolsPro” w Arcgis
Dr. ADPrasad

Odpowiedzi:


22

Należy pamiętać, że dane te są próbkami dyskretnych domen litologicznych. Często granica między dwiema takimi domenami nie może zostać zidentyfikowana w terenie, dlatego nie można oczekiwać, że wiele przykładowych lokalizacji będzie leżało dokładnie wzdłuż granic. Prawidłowym rozwiązaniem będzie podział obszaru badań, a każdy wielokąt w obrębie tego podziału może (i często będzie) wykraczać poza lokalizacje próbek, które go określają. Z wyjątkiem przybliżonych przybliżeń, wyklucza to wszelkie podejście, które wykorzystuje lokalizacje próbek jako wierzchołki powstałych wielokątów .

W przypadku prac o wysokiej jakości najlepszą metodą jest dopasowanie uogólnionego liniowego modelu przestrzennego do procesu wielomianowego. Ta procedura wymaga znacznej wiedzy i wysiłku. Jako substytut, możesz rozważyć rozszerzenie każdego punktu próbki do jego wielokąta wpływu (inaczej wielokąta Thiessena, wielokąta Voronoi lub komórki Dirichleta). Dobrym pomysłem jest ograniczenie ekspansji do obszarów lądowych; można to zrobić za pomocą siatki maski.

Aby to zilustrować, rozważ ten znacznie mniejszy zestaw danych (z 14 136 punktów) reprezentujący 12 klas litologicznych w odróżnieniu od koloru:

Próbki

Oto szczegół ze środka płata wschodniego, pokazujący nieregularne położenie punktów i stosunkowo szybkie zmiany litologii w tym miejscu. Ręczne śledzenie tego byłoby trudną i arbitralną procedurą:

Próbki - szczegóły

Dokonałem ekspansji, przekształcając te punkty w siatkę (około 800 wierszy i 1000 kolumn) i obliczając ich alokację euklidesową , używając maski, która ograniczyła obliczenia do gruntów nie zlodowaconych. (Schemat kolorów na następnych dwóch rysunkach różni się od schematu na poprzednim.)

Alokacja euklidesowa

Dla porównania, oto szczegółowa mapa litologiczna tego samego obszaru narysowana w tej samej skali z tą samą symboliką:

Oryginalna mapa

Dzięki naprawdę dużemu zestawowi danych lub zawiłemu obszarowi badań, może być szybkie ułożenie regionu i wykonanie tej procedury osobno na każdym kafelku, w razie potrzeby mozaikując wyniki w jeden wyjściowy raster. Aby to zadziałało, płytki muszą nieznacznie zachodzić na siebie, aby uniknąć efektów krawędziowych (a następnie powinny być równomiernie przycięte przed mozaikowaniem).

Głównymi powodami, dla których warto przejść do reprezentacji rastrowej, są (1) szybkie i łatwe do obliczenia oraz (2) trudne do znalezienia dokładne rozwiązania oparte na wektorze. Jeśli spróbujesz zderzaków, wypukłych kadłubów, wklęsłych kadłubów itp., Przekonasz się, że wszystkie one przecinają się wzajemnie i nadal pozostawiają luki: innymi słowy, nie wytworzą spójnego topologicznie podziału przestrzeni na odrębne domeny litologiczne.

Jedną z metod opartych na wektorze, która będzie działać, jest obliczenie ograniczonej teselacji punktów Voronoi ( dobre metody wymagają O (n * log (n)) czasu na n punktów), przestrzenne scalenie komórek Voronoi zgodnie z atrybutami litologicznymi ich powiązanych punkty, a następnie podziel powstałe wieloboki na połączone ze sobą komponenty (jeśli chcesz). Jeśli jednak wszystko, czego potrzebujesz, to wyjście wektorowe , łatwiej jest pogrupować wyniki rastra w region i przekonwertować je na format wektorowy.


Dzięki za miłe wyjaśnienie pozwól mi spróbować zrobić coś podobnego do tego.
Wielkie

jak stworzyłeś siatkę, szukam narzędzi, które konwertują punkty na siatkę. Dzięki
Ramakrishna Billakanti,

I zapisywane punkty w postaci siatki . Ta procedura inicjuje każdą komórkę za pomocą NoData, a następnie aktualizuje komórki zawierające punkty z atrybutami punktów: nie jest wymagana interpolacja.
whuber

Czy możesz określić narzędzie używane do zapisywania punktów na siatce. Przepraszam, że przeszkadzam, ale jestem naprawdę słaby w analizie rastrowej. Jeszcze raz dziękuję.
Ramakrishna Billakanti

Proszę kliknąć link w poprzednim komentarzu do strony pomocy.
whuber

3

Centrum zasobów ESRI posiada narzędzie do tworzenia „wklęsłego kadłuba”. Może to spowodować powstanie wielokąta, który lepiej pasuje do krawędzi twoich punktów w porównaniu do wypukłego kadłuba. Dane wejściowe to klasa obiektów punktowych i tworzy wielokąt.
Estymator kadłuba wklęsłego


3

Możesz spróbować utworzyć wielokąty Thiessena z punktów, a następnie rozpuścić powstałe wielokąty za pomocą atrybutu typu skały.


Omawiam to rozwiązanie na końcu mojej odpowiedzi: wielokąty Thiessena są synonimem teselacji Voronoi. Dziękujemy za udostępnienie linku do oprogramowania!
whuber

1
Wymieniony wątek wymaga licencji ArcInfo, jeśli masz licencję Arcview, możesz wypróbować resources.arcgis.com/gallery/file/geoprocessing/…

Byłbym bardzo zainteresowany, @Dan, aby dowiedzieć się, jak dobrze twoje rozwiązanie radzi sobie z „milionami punktów”. :-) Wygenerowanie tak dużego zestawu danych do testów może być uciążliwe; więc, Ramakrishna, jeśli spróbujesz któregokolwiek z wymienionych tutaj rozwiązań wektorowych, czy byłbyś tak uprzejmy, aby powiedzieć nam o ich wydajności później?
whuber

1
_ @ Bill, nie obraziłem się, termin wyklucza ocenę prezentu, ale teraz jest to lepka notatka na moim monitorze :) Zgłosię się znacznie później. Tymczasem każdy, kto ma licencję ArcInfo, może zgłosić czasy implementacji Arc dla plików punktowych o różnych rozmiarach, dzięki czemu mógłbym porównać z nią czystą implementację Pythona. Pozdrawiam

1
@whuber Narzędzia ArcGIS nie zajęły dużo czasu na przetworzenie 28 milionów rekordów. Szybsze było odczytywanie punktów i przetwarzanie ich do siatki, a następnie do alokacji euklidesowej. Naprawdę doceniam to, że poświęciłem czas i opublikowałem swoje odpowiedzi na blogu. Jeszcze raz dziękuję.
Ramakrishna Billakanti

2

Interpolacja 3D w celu zbudowania brył, a następnie przekrój przekroju poziomego na głębokości w celu przecięcia litologii w celu uzyskania wielokąta - Na podstawie ostatnich komentarzy wydaje się, że mamy do czynienia z danymi wiercenia 3D. Oznacza to, że najpierw musisz zbudować bryły 3D (siatki triangulowane) z danych. Można to zrobić na 2 sposoby: Digitalizacja kontaktów w przestrzeni 3D w celu utworzenia brył litologii lub interpolacji 3D. Aby to zrobić ręcznie, będziesz potrzebować oprogramowania, takiego jak GEOMCOM GEMS lub podobny, a jedynym pakietem, który wiem, że może to robić dynamicznie, jest Leapfrog Mining. (Którego używam) Modele ręczne są zwykle prostsze i pozwalają na ludzką interpretację geologii, ale aktualizacja jest trudna. Konfiguracja modeli dynamicznych może trochę potrwać, ale wraz z postępem programów i udostępnieniem nowych danych możesz po prostu zaktualizować i zregenerować zaktualizowane modele litologiczne. Obie techniki są tutaj dość skomplikowane do wyjaśnienia. Prawdopodobnie zaleciłbym ręczną digitalizację w przypadku prostszych modeli lub tam, gdzie nie będziesz aktualizować modelu o nowe dane. Leapfrog Mining jest naprawdę dobrze dopracowaną aplikacją i ma całą strukturę domen 3D, w której możesz poprawnie zdefiniować litologię, na przykład od najstarszej do najmłodszej, ale wymaga szkolenia, aby zrozumieć tę koncepcję.

Po utworzeniu modelu litologii wystarczy po prostu utworzyć wycinek przekroju zorientowany poziomo na określonej głębokości. Następnie możesz wyeksportować kontury litologii do wielokątów, które stanowią podstawę mapy litologii. Możesz kilkakrotnie obniżyć plasterek w odstępach, aby następnie porównać zmiany litologii na różnych głębokościach.

Można to również wizualizować w Leapfrog, ale często eksportuję modele 3D do DXF i używam innych aplikacji, takich jak Geosoft Target, do tworzenia bardziej tradycyjnych przekrojów przy użyciu tych brył.

Opisałem tylko mój przepływ pracy, ale jestem pewien, że istnieją inne rozwiązania; możliwe jest wygenerowanie modelu litologii w Geosoft Target, ale nie próbowałem po prostu dlatego, że nie lubię pracować z tym oprogramowaniem. Jestem prawie pewien, że Datamine Studio może być również używane, a nawet przewyższyć moje techniki.


-1

Możesz spróbować użyć Kreatora bufora w ArcMap.

Możesz spróbować użyć dowolnego rodzaju wskaźników lub odległości, ale prawdopodobnie musisz je dopasować lub użyć szablonu zasięgu, jeśli taki istnieje.

Możesz go znaleźć w ArcScripts lub geoprzetwarzaniu ze strony wsparcia ESRI.

Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.