Erstellen von Polygonen um die Ausdehnung ähnlicher Punkte mit ArcGIS Desktop?

11

Ich habe Millionen von Punkten für die Lithologieschicht gesammelt.

Sie haben eine Codierung für verschiedene Gesteinsarten verwendet.

Ich muss ein Polygon um die Ausdehnung ähnlicher Punkte erstellen.

Was ist der einfachste Weg, um das Polygon zu erhalten, anstatt manuell zu digitalisieren?

Ich habe nach Werkzeugen gesucht, die Punkte in Polygone konvertieren, aber es sieht so aus, als gäbe es keine.

Ich habe Werkzeuge zum Konvertieren von Linien in Polygone, Polygone in Linien und Punkte gesehen, aber nicht für Punkte in Polygone.

Ramakrishna Billakanti
quelle
4
Bitte lesen Sie den Moderator-Kommentar zu Ihrer vorherigen Frage . Wenn Sie die Regeln befolgen, wissen die Leute das zu schätzen und geben mit größerer Wahrscheinlichkeit schnelle Antworten, wenn Sie sie wirklich brauchen.
whuber
Ist Ihr Problem "Ich muss ein Polygon um die Ausdehnung ähnlicher Punkte erstellen"? Versuchen Sie es mit einem Rechteck (Punkt (min (x), min (y)), Punkt (max (x), max (y))) oder genauer.
Remigijus Pankevičius
In welcher Umgebung bist du?
Ragi Yaser Burhum
Es hört sich so an, als ob Sie eine konvexe Rumpffunktion wünschen. Wenn Sie danach suchen, wurde diese Frage bereits beantwortet .
begrüßt
Versuchen Sie "ein Polygon aus Punkten machen" von "XToolsPro" in arcgis
Dr. ADPrasad

Antworten:

22

Wir müssen berücksichtigen, dass diese Daten Proben diskreter lithologischer Domänen sind. Oft kann die Grenze zwischen zwei solchen Domänen im Feld nicht identifiziert werden, und daher ist nicht zu erwarten, dass viele der Probenorte genau entlang der Grenzen liegen. Eine korrekte Lösung ist eine Partition des Untersuchungsgebiets, und jedes Polygon innerhalb dieser Partition kann (und wird häufig) über die Positionen der Proben hinausgehen , die es bestimmen. Mit Ausnahme von groben Näherungen schließt dies jeden Ansatz aus, bei dem die Stichprobenpositionen als Eckpunkte der resultierenden Polygone verwendet werden .

Für qualitativ hochwertige Arbeiten besteht die beste Methode darin, ein verallgemeinertes lineares räumliches Modell für einen multinomialen Prozess anzupassen. Dies ist ein Verfahren, das beträchtliches Fachwissen und Aufwand erfordert. Als Ersatz können Sie erwägen , jeden Stichprobenpunkt in sein Einflusspolygon (auch bekannt als Thiessen-Polygon, Voronoi-Polygon oder Dirichlet-Zelle) zu erweitern. Es ist eine gute Idee, die Expansion auf Landflächen zu beschränken. Dies kann mit einem Maskenraster erfolgen.

Betrachten Sie zur Veranschaulichung diesen viel kleineren Datensatz (von 14.136 Punkten), der 12 lithologische Klassen darstellt, die sich durch Farbe unterscheiden:

Proben

Hier ist ein Detail aus der Mitte des Ostlappens, das die unregelmäßigen Positionen der Punkte und die relativ schnellen Änderungen der Lithologie dort zeigt. Dies manuell zu verfolgen wäre ein schwieriges und willkürliches Verfahren:

Proben - Detail

Ich habe die Erweiterung durchgeführt, indem ich diese Punkte in ein Raster (ca. 800 Zeilen und 1000 Spalten) umgewandelt und ihre euklidische Zuordnung mithilfe einer Maske berechnet habe, die die Berechnung auf nicht vergletschertes Land beschränkte. (Das Farbschema in den nächsten beiden Abbildungen unterscheidet sich von dem der vorhergehenden.)

Euklidische Allokation

Zum Vergleich hier eine detaillierte lithologische Karte desselben Gebiets im gleichen Maßstab mit derselben Symbolisierung:

Originalkarte

Bei einem wirklich großen Datensatz oder einem verschlungenen Untersuchungsgebiet kann es schnell sein, die Region zu kacheln und dieses Verfahren für jede Kachel separat durchzuführen, wobei die Ergebnisse auf Wunsch in einem Ausgaberaster mosaikiert werden. Damit dies funktioniert, müssen sich die Fliesen leicht überlappen, um Randeffekte zu vermeiden (und sollten dann vor dem Mosaikieren gleichmäßig zugeschnitten werden).

Die Hauptgründe für den Besuch einer Rasterdarstellung sind: (1) die schnelle und einfache Berechnung und (2) genaue vektorbasierte Lösungen sind schwer zu finden. Wenn Sie Puffer, konvexe Hüllen, konkave Hüllen oder was auch immer ausprobieren, werden Sie feststellen, dass sich alle gegenseitig überschneiden und dennoch Lücken hinterlassen. Mit anderen Worten, sie erzeugen keine topologisch konsistente Aufteilung des Raums in verschiedene lithologische Domänen.

Eine vektorbasierte Methode, die funktioniert, besteht darin, eine eingeschränkte Voronoi-Tessellation der Punkte zu berechnen ( gute Methoden benötigen O (n * log (n)) Zeit für n Punkte) und die Voronoi-Zellen räumlich gemäß den lithologischen Attributen ihrer zugehörigen zu verschmelzen Punkte, und trennen Sie dann die resultierenden Multi-Polygone in ihre verbundenen Komponenten (wenn Sie möchten). Wenn alles , was Sie benötigen jedoch Vektor ausgegeben , dann ist es einfacher , das Raster Ergebnisses und Konvertit Region dass auf Vektor - Format.

whuber
quelle
Vielen Dank für eine nette Erklärung. Lassen Sie mich versuchen, etwas Ähnliches zu tun. Vielen Dank
Ramakrishna Billakanti
Wie haben Sie das Raster erstellt? Ich suche nach Werkzeugen, die Punkte in Raster konvertieren. Vielen Dank
Ramakrishna Billakanti
Ich habe die Punkte als Raster gespeichert . Diese Prozedur initialisiert jede Zelle mit NoData und aktualisiert dann die Zellen, die Punkte enthalten, mit den Punktattributen: Es ist keine Interpolation beteiligt.
whuber
Können Sie das Werkzeug angeben, mit dem Sie die Punkte im Raster gespeichert haben? Es tut mir leid, dass ich Sie gestört habe, aber ich bin wirklich schlecht in der Rasteranalyse. Danke nochmal.
Ramakrishna Billakanti
Bitte folgen Sie dem Link in meinem vorherigen Kommentar zur Hilfeseite.
whuber
3

Das ESRI Resouce Center verfügt über ein Tool zum Erstellen eines "konkaven Rumpfs". Dies kann zu einem Polygon führen, das sich im Vergleich zu einer konvexen Hülle besser an die Kante Ihrer Punkte anpasst. Die Eingabe ist eine Punkt-Feature-Class und erzeugt ein Polygon.
Konkaver Rumpfschätzer

klewis
quelle
3

Sie können versuchen, Thiessen-Polygone aus den Punkten zu erstellen und dann die resultierenden Polygone mithilfe des Gesteinstypattributs aufzulösen .

nmpeterson
quelle
Ich diskutiere diese Lösung am Ende meiner Antwort: Thiessen-Polygone sind gleichbedeutend mit einer Voronoi-Tessellation. Vielen Dank für die Bereitstellung eines Links zur Software!
whuber
1
Der aufgelistete Thread erfordert eine ArcInfo-Lizenz. Wenn Sie eine Arcview-Lizenz haben, können Sie resources.arcgis.com/gallery/file/geoprocessing/…
Ich wäre sehr interessiert, @Dan, zu wissen, wie gut Ihre Lösung mit "Millionen von Punkten" umgeht. :-) Es kann schwierig sein, einen so großen Datensatz nur zum Testen zu generieren. Also, Ramakrishna, wenn Sie eine der hier genannten Vektorlösungen ausprobieren, wären Sie so freundlich, uns später von ihrer Leistung zu erzählen?
whuber
1
_ @ Bill, ich habe mich nicht beleidigt, Begriff schließt eine Bewertung eines Geschenks aus, aber es ist jetzt eine Haftnotiz auf meinem Monitor :) Ich werde viel später berichten. In der Zwischenzeit konnte jeder mit einer ArcInfo-Lizenz die Implementierungszeiten von Arc für Punktdateien unterschiedlicher Größe angeben, damit ich die reine Python-Implementierung damit vergleichen konnte. Grüße
1
@whuber Die Tools aus ArcGIS haben nicht viel Zeit für die Verarbeitung von 28 Millionen Datensätzen benötigt. Es war schneller, die Punkte zu lesen und sie im Raster und dann in der euklidischen Zuordnung zu verarbeiten. Ich schätze es sehr, dass ich mir Zeit genommen und Ihre Antworten auf dem Blog veröffentlicht habe. Danke nochmal.
Ramakrishna Billakanti
2

3D-Interpolation zum Erstellen von Volumenkörpern und dann horizontaler Schnitt in der Tiefe, um die Lithologie zu durchschneiden und Polygon zu erhalten. Basierend auf Ihren neuesten Kommentaren scheint es sich um 3D-Bohrlochdaten zu handeln. Dies bedeutet, dass Sie zuerst 3D-Volumenkörper (triangulierte Netze) aus Ihren Daten erstellen müssen. Hierfür gibt es zwei Möglichkeiten: Digitalisieren der Kontakte im 3D-Raum, um die Lithologie-Volumenkörper oder 3D-Interpolation zu erstellen. Um dies manuell zu tun, benötigen Sie Software wie GEOMCOM GEMS oder ähnliches. Das einzige mir bekannte Paket, das dies dynamisch ausführen kann, ist Leapfrog Mining. (Was ich benutze) Die manuellen Modelle sind in der Regel einfacher und ermöglichen eine menschliche Interpretation der Geologie, aber die Aktualisierung ist schwierig. Die Einrichtung der dynamischen Modelle kann einige Zeit in Anspruch nehmen. Wenn sich Ihre Programme jedoch weiterentwickeln und neue Daten verfügbar werden, können Sie aktualisierte Lithologiemodelle einfach aktualisieren und neu generieren. Beide Techniken sind hier recht komplex zu erklären. Ich würde wahrscheinlich empfehlen, für einfachere Modelle manuell zu digitalisieren oder wenn Sie das Modell später nicht mit neuen Daten aktualisieren werden. Leapfrog Mining ist eine wirklich ausgefeilte Anwendung und verfügt über eine gesamte 3D-Domänenstruktur, in der Sie beispielsweise die Lithologie vom ältesten bis zum jüngsten korrekt definieren können. Um dieses Konzept zu verstehen, ist jedoch eine Schulung erforderlich.

Sobald Sie Ihr Lithologiemodell erstellt haben, müssen Sie lediglich einen horizontal ausgerichteten Schnitt in einer bestimmten Tiefe erstellen. Anschließend können Sie die Umrisse der Lithologie in Polygone exportieren, die die Grundlage Ihrer Lithologiekarte bilden. Sie können die Schicht in regelmäßigen Abständen mehrmals absenken, um dann zu vergleichen, wie sich die Lithologie in verschiedenen Tiefen ändert.

Dies kann auch in Leapfrog visualisiert werden, aber ich exportiere die 3D-Modelle häufig nach DXF und verwende andere Anwendungen wie Geosoft Target , um traditionellere Abschnitte mit diesen Volumenkörpern zu erstellen.

Ich habe nur meinen Workflow beschrieben, bin mir aber sicher, dass es andere Lösungen gibt. Es ist möglich, ein Lithologiemodell in Geosoft Target zu generieren, aber ich habe es nicht einfach versucht, weil ich nicht gerne mit dieser Software arbeite. Ich bin mir ziemlich sicher, dass Datamine Studio auch verwendet werden kann oder meinen Techniken sogar überlegen sein kann.

Jakub Sisak GeoGraphics
quelle
-1

Sie können versuchen, den Puffer-Assistenten in ArcMap zu verwenden.

Sie können versuchen, jede Art von Metrik oder Entfernung zu verwenden, aber Sie müssen sie wahrscheinlich herausfinden, um sie anzupassen, oder möglicherweise die Extent-Vorlage verwenden, falls es eine gibt.

Möglicherweise möchten Sie es in ArcScripts oder in der Geoverarbeitung auf der ESRI-Support-Website finden.

PROBERT
quelle