3D-Visualisierung von DGM-Daten

Gespeichert von Wibke Steer am
im Forum ArcGIS
Hallo ihr,
ich arbeite mich nach Jahren mit AV 3.2 jetzt in AV 9.2 ein, und es fällt mir doch schwerer als ich dachte. Folgendes Problem ist gerade aktuell:

Ich habe Daten vom LVA im ascii-Format mit der Endung *.xyz, die sehen so aus:

2546000.000 5472000.000 328.759
2546012.500 5472000.000 329.396
2546025.000 5472000.000 329.814
2546037.500 5472000.000 330.115
2546050.000 5472000.000 330.393
2546062.500 5472000.000 330.903
2546075.000 5472000.000 331.332
2546087.500 5472000.000 331.552
2546100.000 5472000.000 331.598
2546112.500 5472000.000 331.668
2546125.000 5472000.000 331.506

usw. (insgesamt 2,5 Mio. Zeilen). So, wie sie sind (mit Leerzeichen vorne und leerzeichengetrennt), kriege ich sie nicht in AV geladen. Ich habe schon mehrere Versuche mit Kommatrennung

(inkl. entsprechender Anpassung in schema.ini ->
[test3.txt]
Format=Delimited(,)
DecimalSymbol=.)

hinter mir, die Daten werden dann auch gelesen, aber die Performance läßt schwer zu wünschen übrig, ich frage mich, ob das nicht irgendwie anders geht. Inzwischen kommt auch die Fehlermeldung, dass die Daten nicht angezeigt werden können, obwohl er sie eben noch angezeigt hat. Puuh! Wisst ihr da was?

Was ich eigentlich will, ist ja dann in etwa: die Punkte in ein Grid umwandeln (habe ich schon versucht, da kommt nur eine Fehlermeldung, ich weiß auch nicht, ob man da ascii to raster oder punkt to raster nimmt, oder was ganz anderes), um aus diesem ein TIN zu erstellen, welches dann die 3D-Oberfläche für meine Visualisierung bilden soll. Ist das der richtige Weg? Wie kriege ich diese Daten gehandelt? Hilfe!

Danke für eure Hilfe
Wibke

Gespeichert von Obiwan am Di., 17.07.2007 - 11:13

Permalink
Ersetze die Dezimalpunkte durch Kommas mit einem Text Editor.
Der Import müsste dann klappen. s. auch hier

Gespeichert von Wibke Steer am Di., 17.07.2007 - 13:19

Permalink
hallo nochmal,

das einlesen klappt inzwischen einigermaßen, die Daten sind jetzt kommagetrennt, haben Spaltenüberschriften (x,y,z) und die schema.ini passt auch. Aber jetzt, wie weiter?

Gruß
Wibke

Gespeichert von Wibke Steer am Di., 17.07.2007 - 14:29

Permalink
Hallo,

danke für die Antworten. Das Raster habe ich inzwischen, leider kommt bei dem Versuch, dieses in ein TIN umzuwandeln, folgende Meldung:

Executing (RasterTin_19): RasterTin Natural_DGM52 "natural_dgm52_RasterTin1" 26.6848030090332 1500000 1
Start Time: Tue Jul 17 14:31:47 2007
System Error (0): Unexpected Error.
Failed to execute (RasterTin_19).
End Time: Tue Jul 17 14:31:48 2007 (Elapsed Time: 1.00 seconds)

Was heißt jetzt das? Ich bin für jeden Hinweis dankbar!

Gruß
Wibke

Gespeichert von Obiwan am Di., 17.07.2007 - 14:58

Permalink
Verkürze mal den Namen Deines Rasters.

Gespeichert von Wibke Steer am Di., 17.07.2007 - 15:20

Permalink
Hey, deine Antworten sind ja kurz und knackig, was? :-)

Den Rasternamen denkt sich das Programm selbst aus! Also, das TIN habe ich jetzt auch. Blödes neues Problem: (ich bin echt neu in AV9)

Wie erstellt man ein neues Punktthema (möchte damit Sichtbarkeitsanalyse machen)?

Gruß
Wibke

Gespeichert von Maxilla am Di., 17.07.2007 - 18:03

Permalink
... Über ArcCatalog - Rechtsklick auf Ordner - New - Shapefile
... Über ArcToolbox - CreateFeatureDataset (in eine GDB)
... Über den Enhanced Shapefile Creator
... Über eine Text- oder dBase-Datei (siehe oben)
... Über 'Export Data' eines bestehenden Punktlayers
...

Gespeichert von Wibke Steer am Mi., 18.07.2007 - 08:51

Permalink
Hey, danke,

aber leider weiß ich jetzt nicht, wie ich dem Shape Daten hinzufüge. Ich habe meinen Punkt als Grafik erstellt und will ihn jetzt einfach in das Shape bringen. Wie editiere ich???

Gruß
Wibke

Gespeichert von Wibke Steer am Mi., 18.07.2007 - 09:07

Permalink
Alles zurück, ich habs gefunden! Bin wohl auch schon blind! (peinlich)

Gespeichert von Maxilla am Mi., 18.07.2007 - 09:24

Permalink
Nu hab ich hier schon die ganze Batterie aufgefahren :)

Wenn's nur ein Punkt ist: Bearbeitung des Punktshapes starten, mit dem Bleistift den Punkt malen, Doppelklick drauf, Cursor auf dem Punkt lassen, Rechtsklick - Properties öffnen und die genauen Koordinaten eingeben.

Ansonsten kannst du zur Umwandlung von Grafiken in Features die entsprechende XToolsPro-Funktion, Dieses Skript oder Convert Graphics to Features benutzen.

Gespeichert von Wibke Steer am Mi., 18.07.2007 - 09:53

Permalink
Danke für deine Mühe! Leider komme ich mit dem Versuch, die Extensions herunterzuladen nicht weit, da im Browser nichts angezeigt wird. Aber erstmal egal.

Wie kriege ich denn eine Sichtbarkeitsanalyse erstellt, die einigermaßen realistisch ist, mit den folgenden Vorgaben: der Punkt ist ein Turm, der 160m über das Geländeniveau hinausragt, und von dieser Höhe aus soll die Sichtbarkeit berechnet werden. Das geht irgendwie mit Viewshed im Spatial Analyst, aber da habe ich bisher nichts brauchbares herausbekommen. Was besagt z.B. der Z-Faktor? Ich dachte, das ich damit die Höhe der "Augen", also meine 160m einstellen kann. Scheint aber nicht so...

Gruß
Wibke

Gespeichert von Maxilla am Mi., 18.07.2007 - 11:05

Permalink
Nö, das ist nicht so. Der Z-Faktor gibt die Überhöhung des DGM in der Berechnung an und ist mit dem Defaultwert 1 meist ganz gut belegt. Die Beobachterhöhe stellst du ein, indem du im Punktlayer ein Feld namens OFFSETA anlegst, und dort einen Wert einträgst (in der ArcGIS-Desktop Hilfe steht noch mehr dazu).

Gespeichert von Wibke Steer am Do., 19.07.2007 - 15:06

Permalink
Hallo nochmal,

ich habe inzwischen die Sichtbarkeit berechnet, indem ich aus der Koordinatenliste zunächst mal ein Shape erstellt habe (das hat vier Stunden gedauert!), in welchem ich dann einige Punkte um meinen Standort herum um 160m höher gemacht habe (Z-Feld), das Abspeichern der Änderung hat nochmal 2 Stunden gedauert. Dann habe ich aus diesem Shape mein Grid berechnet, mit derselben Zellgröße wie der Abstand zwischen den Punkten war (12,5m) und auf diesem habe ich dann meinen Punkt auf den "Buckel" gesetzt und die Sichtbarkeit berechnet. Funktioniert soweit.

Jetzt weiß ich aber nicht, wie das mit den Genauigkeitsanforderungen bei solchen Aufgaben ist. Welche Parameter muss man da wie einstellen, bzw. was besagen die? Und das Grid selbst kann man auch auf verschiedene Arten herstellen (Natural Neighbor, was macht das genau?, oder mit dem Conversion-Tool "to Raster" oder, oder).

Wer kann mir da was genaueres sagen? Maxilla, du kennst dich doch da anscheinend gut aus. Ich habe leider auch nur eine englische Oberfläche und kenne mich mit den mathematischen Methoden nicht so aus.

Gruß
Wibke

Gespeichert von Maxilla am Do., 19.07.2007 - 15:44

Permalink
Falls dein Höhenmodell ein DGM5 ist, musst du daraus erst mal ein digitales Oberflächenmodell (DOM, DSM) erzeugen, indem du die Strukturhöhen draufaddierst, sonst ist's mit der Genauigkeit von vornherein nicht weit her... Ansonsten verweise ich zunächst mal auf diesen Thread: Sichtbarkeitsanalysen

Gespeichert von Wibke Steer am Fr., 20.07.2007 - 15:10

Permalink
Vielen Dank für Deine Hilfe, Maxilla, so langsam fange ich an, die Problematik des ganzen zu erkennen. Ich habe aber doch immer mal wieder konkrete Fragen, ich hoffe, es macht dir nichts aus. Zunächst mal:

- mit welcher Methode mache ich am besten aus meinen Höhenpunkten das Grid?
- wie berechne ich denn ein Oberflächenmodell, welche Daten brauche ich dafür, Landnutzungsdaten im Vektorformat, z.B. ATKIS?

Weiter geht´s am Montag! :-)

Gruß
Wibke

Gespeichert von Maxilla am Fr., 20.07.2007 - 16:40

Permalink
Nein, das macht mir nichts aus (brauche ja auch nicht zu antworten ;o)

Da du genau einen Punkt pro Rasterzelle hast, ist eine Interpolation erst mal überflüssig und du kannst einfach aus der SpatialAnalyt-Toolbar: Convert - Features to Raster benutzen.

Ein DSM kannst du (eventuell) genauso kaufen wie ein DGM - schön z.B. ein aus Laserscanbefliegungen generiertes Modell - du kannst dir aber auch selber eins generieren, wenn dir Strukturhöhen beliebiger Quelle vorliegen.

Wenn diese direkt nicht zu bekommen sind, kannst du sie aus Landnutzungsdaten erzeugen indem du z.B. mittlere klassifizierte Nutzungshöhen vergibst (so in der Art: Wald = 12,5m usw).

ATKIS-B-DLM oder ALK-Nutzungsdaten sind zwei Möglichkeiten, es ist ein wenig lästig sie aufzubereiten, weil sie sehr nutzungs- und wenig strukturbezogen sind (Attribute wie z.B. 'Elektrizität' und so was), aber es geht schon.

Wenn das Gebiet nicht zu groß ist, kann man's auch einfach von der TK25 abdigitalisieren. Außerdem könntest du noch automatisch Satelitenbilder auswerten oder die freien CORINE-Landnutzungsdaten verwenden. Und sicherlich gibt's noch zahlreiche andere Quellen, die mir gerade nicht einfallen...

Wie auch immer, letzendlich machst du ein Raster draus und addierst es auf das DGM.

Schönes Wochenende,

Mx.

Gespeichert von Wibke Steer am Mo., 23.07.2007 - 07:58

Permalink
Guten Morgen in der neuen Woche!

Danke, Maxilla, es hat sich herausgestellt, dass ich zunächst mal die Strukturhöhen nicht brauche. Ich melde mich wieder, wenn ich neue Fragen habe (was bestimmt bald der Fall ist!).

Gruß
Wibke

Gespeichert von alte Banane am Mi., 12.09.2007 - 23:24

Permalink
...das mit den Strukturhöhen ist doch wohl nen Selbstgänger.
Als wenn Du in Rottenburg einen abgebröselt hättest...

Gespeichert von Matze am Fr., 07.03.2008 - 09:28

Permalink
Hallo Wibke,

bin hier mal mit großem Interesse auf dieses Thema bei Dir gestoßen - das Problem türmt sich nun gerade bei mir auf ;-) Hoffe, Du bist inzwischen im Bereich 3D, Höhendaten, Sichtbarkeitsanalysen fit.
Bezüglich Leitungsmastenbau > möchte ich demnächst evtl. über Spatial Analyst ein ähnliches ergebnis wie Du realisieren. Gelingt sowas überhaupt über ArcViewGis 3.3 bzw. 9.1 (von dem ich noch kaum Ahnung habe)?

Mal beste Grüße an der stelle von mir!!!!

Matze