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:
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ą:
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.)
Dla porównania, oto szczegółowa mapa litologiczna tego samego obszaru narysowana w tej samej skali z tą samą symboliką:
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.