Verschneidung DGM mit Nutzungsshape

Gespeichert von Margret am
im Forum ArcGIS
Hallo,
ich möchte eine Sichtbarkeitsanalyse für Windkraftanlagen erstellen. Hierzu habe ich ein DGM50, das ich mit Hilfe von IDW zu einem Raster interpoliert habe und dann habe ich ein DLM25 als shape, aus dem ich die Wälder und die Siedlungen herausgefiltert habe. Wie kann ich daraus ein Oberflächenmodell erstellen? Und wie weise ich den Nutzungsshapes die jeweilige Höhe zu?

Gespeichert von Maxilla am Mi., 13.04.2011 - 12:22

Permalink
Erzeuge im DLM25 ein Zahlenfeld, selektiere nutzungsspezifisch die Datensätze und weise der jeweiligen Nutzung eine Strukturhöhe zu (z.B. 12 m für Wald usw).

Stelle Analyseumgebung und Zellgröße des Spatial Analysten entsprechend des DGM50-Rasters ein. Wandle dein DLM25 über das Feld mit der Strukturhöhe in ein Raster um (Featur to Raster)

Öffne den Raster Calculator und berechne [DGM50] + [Strukturhöhe]. Damit hast du dann ein vereinfachtes Oberflächenmodell (DOM) erzeugt.

Gespeichert von Margret am Mi., 13.04.2011 - 12:34

Permalink
Vielen Dank für die schnelle Antwort
Die Strukturhöhe habe ich schon mal zugewiesen.

Zu Punkt 2: Meinst du damit, dass ich die Analyseumgebung und die Zellgröße des DLM auf das DGM50-Raster einstellen muss? Wo und wie kann ich diese Einstellungen ändern und woher weiß ich, wie meine Einstellungen im DGM sind?

Mit welchem Befehl wandle ich dann mein DLM25 in ein Raster um?

Du merkst, ich bin leider nicht so versiert bei diesen Anwendungen.

Gespeichert von Maxilla am Mi., 13.04.2011 - 13:24

Permalink
Rasterweite des DGM steht in seinen Layer-Eigenschaften (Rechtsklick - Eigenschaften - 'Source'-Registerkarte)

Analyseumgebung und Zellgröße:

Kannst du in 9.x über das Spatial-Analyst-Menü unter Optionen einstellen: Dort Extent auf 'Same as Layer DGM' und Cellsize entsprechend der DGM-Auflösung einstellen. In ArcGIS 10 liegt das unter dem Eintrag 'Environments' im 'Geoprocessing'-Menü

Die Umwandlung machst du mit dem Tool 'Feature to Raster'. Dazu einen Tipp: Tools, deren Namen du kennst, findest du am Einfachsten über die Index-Registerkarte der ArcToolbox oder - in ArcGIS 10 - über die Search-Box.

Gespeichert von Margret am Mi., 13.04.2011 - 14:19

Permalink
Es hat nun alles geklappt bis auf die Berechnung DGM50 + Strukturhöhe (Beispielhaft habe ich zunächst die Siedlungen verwendet). Eigentlich müsste doch als Ergebnis die Fläche des DGM angezeigt werden, welches die Siedlungen mit der entsprechenden Höhe enthält aber auch alle anderen Flächen, oder?
Mein Ergebnis enthält jedoch nur die Siedlungen, die wahrscheinlich die Höhe des DGM enthalten, aber alle anderen Bereiche außerhalb der Siedlungen sind hier nicht enthalten.
Ich hoffe, das war verständlich beschrieben, wo das Problem liegt?
Vielen Dank für deine Mühe.

Gespeichert von Maxilla am Mi., 13.04.2011 - 18:08

Permalink
Liegt wahrscheinlich an den NoData-Werten deines Siedlungs-Rasters. Du solltest alle NoDatas in diesem Raster vor der Addition auf den Wert 0 setzen. Das geht z.B. mit der Raster-Calculator-Formel:

Con ( isNull ( [DeinGrid]) , 0 , [DeinGrid])

Von vornherein vermeiden könntest du die NoDatas, indem du im Ausgangs-Shape einfach alle Nutzungen, die nicht berücksichtigt werden sollen mit der Strukturhöhe 0 belegst, statt sie "herauszulöschen".

Gespeichert von Margret am Do., 14.04.2011 - 11:24

Permalink
Wenn ich das richtig verstehe, sind NoData-Werte, Zeilen, die keine Inhalte zur Siedlung enthalten? Ich glaube, dass es in meinem shape keine NoData-Werte gibt.
Mein DLM besteht aus ganz vielen verschiedenen shapes und ich habe zum Test erstmal nur das Shape Siedlungen verwendet. Ich habe also keine Werte herausgelöscht. Bedeutet das, dass ich für die Addition aus DGM und Nutzung ein vollständiges Nutzungsshape haben muss? Ich kann also nicht mein DGM mit dem shape Siedlungen addieren?

Gespeichert von Maxilla am Do., 14.04.2011 - 11:41

Permalink
Nein, die NoData-Werte sind wie gesagt nur im Siedlungs-Raster, das du aus dem Siedlungs-Shape erstellt hast - und zwar überall dort, wo keine Siedlung ist.

Und nein, du musst kein vollständiges Shape haben. Setz einfach, die Nutzung, die du addieren möchtest in ein Raster um. In diesem Fall entstehen aber dort, wo keine Siedlung ist, im Raster NoData-Werte, die du wie gesagt mit 0 ersetzen musst.

Gespeichert von Margret am Do., 14.04.2011 - 11:52

Permalink
Das kapiere ich noch nicht ganz. Denn mein shape besteht ja nur aus Siedlungen. Das shape enthält ja gar keine anderen Nutzungen. Ich habe also nur die Siedlungen in ein Raster umgewandelt. Ich habe jetzt auch einfach mal versucht diese Formel auf mein Siedlungsraster anzuwenden, aber da bekomme ich nur eine Fehlermeldung.

Gespeichert von Maxilla am Do., 14.04.2011 - 12:04

Permalink
Ein Raster ist immer flächendeckend und rechteckig. Wenn du dein Siedlungsshape in ein Raster umwandelst, hast du im Bereich der Siedlungen Strukturhöhen, und im Bereich zwischen den Siedlungen nicht "nichts" sondern eben NoData-Werte. Da muss 0 hin. Wegen der Formel: Hast du da an der Stelle [DeinGrid] den Namen deines Siedlungsrasters eingefügt? (In ArcGIS werden statt der eckigen Klammern übrigens Anführungsstriche um die Grids gesetzt).

Gespeichert von Margret am Do., 14.04.2011 - 12:19

Permalink
Ja, ich habe genau die Formel eingegeben, die du mir vorgegeben hast und habe mein Raster eingefügt, aber in der Fehlermeldung sagt er mir, dass er das input-grid nicht öffnen kann.

Gespeichert von Maxilla am Do., 14.04.2011 - 12:23

Permalink
Kannst du bitte mal deine Formel hier reinposten?

Gespeichert von Margret am Do., 14.04.2011 - 12:27

Permalink
Con(isNull([siedl_r],0,[siedl_r]))
Ich habe auch mal Leerzeichen eingefügt, oder mal hinten eine Klammer wegeglassen, aber es kommt immer ne Fehlermeldung

Gespeichert von Maxilla am Do., 14.04.2011 - 12:33

Permalink
Okay, das ist nicht die Formel, die ich oben gepostet habe (beachte die schließenden Klammern!).
Richtig ist:

Con ( isNull ( [siedl_r]) , 0 , [siedl_r])

Gespeichert von Margret am Do., 14.04.2011 - 13:05

Permalink
Entschuldigung, du hast Recht, ich habe die Klammer mittendrin vergessen. Vielen Dank

Gespeichert von Maxilla am Do., 14.04.2011 - 13:34

Permalink
Gut, hoffe, es klappt jetzt...

Gespeichert von Margret am Do., 14.04.2011 - 16:40

Permalink
Ja es hat nun alles, auch die Sichtbarkeitsanalyse sieht nun ganz brauchbar aus.
Víelen Dank für deine Hilfe