Master Thesis
im Rahmen des
Universitätslehrganges „Geographical Information Science & Systems“
(UNIGIS MSc) am Zentrum für GeoInformatik (Z_GIS) der Paris Lodron-Universität Salzburg
zum Thema
„Kanalnetzüberrechnung mit Hilfe von GIS“
vorgelegt von
Bernhard Resch
UP10568, UNIGIS MSc Jahrgang 2012
Zur Erlangung des Grades
„Master of Science (Geographical Information Science & Systems) – MSc(GIS)”
Gutachter:
Ao. Univ. Prof. Dr. Josef Strobl
Neukirchen an der Vöckla, 18.02.2015
Erklärung
Ich versichere, diese Master Thesis ohne fremde Hilfe und ohne Verwendung anderer als der angeführten Quellen angefertigt zu haben, und dass die Arbeit in gleicher oder ähnlicher Form noch keiner anderen Prüfungsbehörde vorgelegen hat. Alle Ausführungen der Arbeit die wörtlich oder sinngemäß übernommen wurden sind entsprechend gekennzeichnet.
Neukirchen am 18.02.2015 Bernhard Resch
Kurzfassung
In den letzten Jahren wurde die Wichtigkeit der Kanalnetzberechnung immer deutlicher.
Die Versieglung von immer mehr Grünflächen führt zu einer erhöhten Wasserableitung in die Kanalisation, dies wiederum führt bei Starkregenereignissen zu einer Überlastung der Kanäle. Überflutungen sind die Folge. Kanalnetzberechnungen dienen zur Feststellung der Leistungsfähigkeit von bestehenden Kanälen sowie zur Dimensionierung von neuen Kanalsystemen.
In der Masterthesis wurde untersucht, ob und wie GIS in der Kanalnetzberechnung eingesetzt werden kann. Es sollten Möglichkeiten gefunden werden, wo GIS einfach und schnell für Zivilingenieure eingesetzt werden kann. Eine Voraussetzung war, dass keine neuen Investitionen in GIS-Software getätigt werden müssen. Aus diesem Grund wurden die Analysen mit Open Source Produkten durchgeführt.
Die Kanalnetzberechnungen erfolgten mit HydroCad, eine Berechnungssoftware, die in die Oberfläche des Barthauer Netzinformationssystems BASYS integriert ist. Für die Berechnung benötigt HydroCad zahlreiche Umweltparameter. Einige dieser Parameter können mit einem GIS sehr gut ermittelt werden. Dies sind der Anteil der befestigten Flächen, die Hangneigung und die Bodenart. Alle diese Parameter wurden den Einzugsflächen mittels Verschneidungen und Überlagerungen zugewiesen.
Die Ermittlung der befestigten Flächen erfolgte durch die Verwendung des NDVI. Für diesen wurden die Orthophotos in echt Farbe und als Infrarot Aufnahmen benötigt. Die Berechnung der Hangneigung erfolgte durch die Erstellung eines DGM, basierend auf einen ALS. Der Bodenaufbau wurde aus der digitalen Bodenkarte von Österreich übernommen. Anhand der ermittelten Werte konnten den Flächen die Oberflächentypen, welche in HydroCad vorab definiert wurden, zugewiesen werden.
Nachdem alle Werte in einer shp-Datei vorhanden waren, wurde eine csv-Datei erstellt.
Diese Datei wurde anschließend mittels ASCII-Import in die Kanaldatenbank importiert. Der Aufbau des ASCII-Imports ist in einer XML-Datei verspeichert, welche an die Anforderungen angepasst wurde.
Abstract
In recent years, the importance of the sewer simulation has become increasingly clear.
The sealing of more and more green spaces leads to increased water discharge into drains, as a result of this heavy rainfall lead to an overloading of the channels. Floods are the result.
Sewer simulations are used to determine the efficiency of existing channels and dimensioning of new sewer systems. There was a research in the master thesis, whether and how GIS can be used in the calculation of sewer networks. Ways should be found where GIS can be used easily and quickly for civil-engineers. One condition was that there would be no new investing in any GIS software. As a reason of this, the analyzes were carried out using open source products. The sewer simulation were performed with HydroCAD, a calculation software, which is integrated into the surface of Barthauer network information system BASYS. HydroCAD requires numerous environmental parameters for the calculation. Some of these parameters can be determined very well by using a GIS-system. These are the percentage of paved areas, the slope and soil type.
All these parameters were assigned to the catchment areas through intersections and overlaps. The determination of the paved areas was carried out by the use of the NDVI.
The requirements were the orthophotos in genuine color and infrared-admittance shots.
The calculation of the slope was carried out by the generation of a DGM, based on an ALS. The floor structure was taken from the digital soil map of Austria. The surfaces could be assigned the surface-types based on the determined values, which were defined in advance by using HydroCAD. After all values were present in a shp file, a csv file was created. This file was imported into the channel database by using ASCII-import.
The structure of the ASCII import is saved as a XML file, which was adapted to the requirements.
Abbildungsverzeichnis
Abbildung 1: Aufbau der Thesis ... 3
Abbildung 2: Ablaufprozess einer Kanalmodellierung aus (HOCHEDLINGER, et al., 2012) ... 5
Abbildung 3: Gegenüberstellung Messdaten und Simulation (aus dem Thesis Projekt) 6 Abbildung 4: Radiale Versetzung durch Höhenunterschiede aus (ALBERTZ, 2009) .. 10
Abbildung 5: Abtastender Laser im Flugzeug, aus (KRAUS, 2004) ... 13
Abbildung 6: First- und Last-Pulse Echo eines Laserstrahls (www.toposys.de, 2014) 13 Abbildung 7: Karte von Oberösterreich Datenquelle: (DORIS, 2013) ... 16
Abbildung 8: Einzugsflächen von Schwanenstadt ... 17
Abbildung 9: HydroCad Einzugsflächenformular ... 18
Abbildung 10: HydroCad Oberflächentypformular ... 20
Abbildung 11: Auswirkungen des Befestigungsgrades ... 21
Abbildung 12: Zusammenspiel von Boden und Gefälle ... 21
Abbildung 13: Landbedeckung eines LISA Testdatensatzes in Oberösterreich (GRILLMAYER, et al., 2010) ... 25
Abbildung 14: vom Orthophoto zum NDVI-Raster ... 26
Abbildung 15: links: digitalisierte befestigte Flächen für EZGB 102a; rechts: NDVI für EZGB102a ... 27
Abbildung 16: NDVI-Wert einer Ackerfläche ... 27
Abbildung 17: links: Gelände Relief des Untersuchungsgebietes; rechts: DGM des Untersuchungsgebietes erstellt in SAGA ... 30
Abbildung 18: Eingabeformular des "Attribute nach Position zusammenführen" in QGIS ... 32
Abbildung 19: Ausschnitt aus der digitalen Bodenkarte des BFW im Untersuchungsgebiet, bei den weißen Flächen sind keine Bodentypen vorhanden ... 33
Abbildung 20: Einteilung der Oberflächentypen ... 34
Abbildung 21: Auszug aus der geänderten XML-Vorlage ... 36
Abbildung 22: Gegenüberstellung der in der Thesis ermittelten Parameter, und für die Berechnung verwendeten Parameter ... 38
Inhalt
Erklärung ... II Kurzfassung... III Abstract ... IV Abbildungsverzeichnis ... V
1.Einleitung ... 1
1.1 Motivation ... 1
1.2 Ziel der Thesis ... 1
1.3 Nicht behandelte Themen ... 1
1.4 Lösungsansatz ... 2
1.5 Werkzeuge... 2
1.6 Erwartete Ergebnisse ... 3
1.7 Zielgruppe ... 3
1.8 Aufbau der Thesis ... 3
2. Theoretische Grundlagen ... 4
2.1 Kanalnetzberechnung ... 4
2.1.1 Ablauf einer Kanalmodellierung ... 4
2.1.2 Einzugsgebiete ... 7
2.1.3 Oberflächenabflussmodellierung ... 7
2.1.4 Berechnungssoftware ... 8
2.1.4.1 BASYS - HydroCad ... 8
2.1.4.2 HYSTEM-EXTRAN ... 8
2.1.5 Hydrodynamische Kanalnetzberechnung mit der GVM ... 9
2.2 Geodaten ... 9
2.2.1 Bodenkarte ... 9
2.2.2 Fernerkundungsdaten ... 10
2.2.2.1 Satellitenbilder ... 10
2.2.2.2 Orthophotos ... 10
2.2.2.3 Auswerten der digitalen Luftaufnahmen ... 11
2.2.2.4 NDVI ... 11
2.2.2.5 Laserscanmessungen ... 12
2.2.3 LISA ... 14
2.3.1 Open Source GIS... 14
2.3.2 XML ... 15
3. Grundlagen für die Analysen ... 16
3.1 Untersuchungsgebiet ... 16
3.2 Erforderliche Daten in der Berechnungssoftware ... 17
3.2.1 Einzugsflächen ... 17
3.2.2 Oberflächentypen ... 18
3.3 Möglichkeiten für den Einsatz von GIS ... 20
3.4 Auswirkungen der gewählten Parameter ... 20
3.5 Datensichtung ... 22
3.5.1 Kanalnetz... 22
3.6.2 Einzugsflächen ... 22
3.6.3 Luftbilder... 22
3.6.4 Digitales Geländemodell ... 23
3.6.5 Bodenkarte ... 23
4. Ermittlung der Parameter ... 24
4.1 Berechnung des Befestigungsgrades ... 24
4.1.1Ermittlung des NDVI ... 25
4.1.2 Grenzwertbestimmung ... 26
4.1.3 Reklassifizierung ... 28
4.1.4 Befestigte Flächen den Einzugsflächen zuordnen... 28
4.2 Berechnung der Hangneigung ... 29
4.2.1 Erstellung eines DGM ... 29
4.2.2 Berechnung des Gefälles ... 30
4.3 Zuordnen des Bodenaufbaus ... 31
4.4 Automatische Zuweisung des Oberflächentyps ... 33
4.5 Automatischer Import der Parameter ... 35
4.5.1 Erstellen einer csv-Datei ... 35
4.5.2 Erweitern einer BASYS-Vorlage ... 35
5. Analyse der Ergebnisse ... 37
5.1 In welchen Bereichen können GIS gestützte Analysen sinnvoll für die hydrodynamische Kanalnetzberechnung eingesetzt werden? ... 37
5.2 Welche Möglichkeiten gibt es, um befestigte Flächen aus Orthophotos zu ermitteln?... 37
5.3 Liefert die Berechnung des Befestigungsgrades anhand des NDVI ausreichend
genaue Ergebnisse für die Kanalnetzberechnung? ... 37
5.4 Wie funktioniert der Datentransfer zwischen GIS und Kanaldatenbank? ... 39
5.5 Kommt es durch den Einsatz der gewählten Analysen zu einer Zeit- und Kostenersparnis ... 39
6. Zusammenfassung und Ausblick ... 41
6.1 Zusammenfassung ... 41
6.2 Ausblick ... 42
Literaturverzeichnis... 43
Anhang ... 46
Erfahrungswerte für die Berechnung nach Horton ... 46
1.Einleitung
1.1 Motivation
Die Masterthesis entstand in enger Zusammenarbeit mit dem Zivilingenieurbüro DI Hitzfelder und DI Pillichshammer aus Vöcklabruck, welches den Auftrag für die hydrodynamische Kanalnetzberechnung von zwei Kanalabschnitten im Verbandsgebiet des Reinhaltungsverband (RHV) Schwanenstadt und Umgebung hatte.
Der RHV Schwanenstadt befindet sich im Bezirk Vöcklabruck, besteht aus 11 Mitgliedsgemeinden, und betreut ein Kanalnetz mit einer Gesamtlänge von ca. 90km.
Teile dieses Kanalnetzes sind bereits über 50 Jahre alt und weisen bereits erhebliche Schäden auf. Um diese Schäden zu erfassen und einen Überblick über den gesamten Zustand des Kanalnetzes zu bekommen, werden in Regelmäßigen Abständen TV- Untersuchungen durchgeführt. Die Ergebnisse der Untersuchungen dienen als Grundlage für ein Sanierungskonzept. Für einige Stränge sieht das Sanierungskonzept eine Sanierung durch Neubau vor. Damit nach der Sanierung, die Leistungsfähigkeit der Kanäle noch ausreichend ist, wird eine Kanalnetzberechnung durchgeführt.
Die Kanalnetzberechnung erfolgt EDV gestützt und erfordert die Eingabe von sämtlichen Kanalstammdaten (Sohlhöhen, Haltungslängen usw.), und zahlreichen Umweltparametern für die angeschlossenen Einzugsgebiete.
1.2 Ziel der Thesis
In der Thesis soll geprüft werden, ob Analysen mit Geoinformationssystemen (GIS), anhand von aufbereiteten Fernerkundungsdaten, die Kanalnetzberechnung unterstützen und somit Arbeitszeit verringern können.
Es soll dies mit einfach nachvollziehbaren und kostengünstigen Methoden geschehen, welche auch von nicht GIS-Experten jederzeit durchgeführt werden können.
1.3 Nicht behandelte Themen
Auf die technischen Details der Kanalnetzberechnung wird in der Thesis nicht genauer eingegangen. Ebenso wird die Erfassung der Fernerkundungsdaten nur kurz erläutert.
1.4 Lösungsansatz
Besonderes Augenmerk der Thesis liegt auf den Parametern der Einzugsflächen und deren Import in die Kanaldatenbank. Beides erfolgt derzeit manuell, und sollte zukünftig mit wenigen Mausklicks erledigt sein.
Folgende Parameter können, auf Grund ihrer geografischen Eigenschaften, mit GIS- Analysen ermittelt werden:
• Befestigungsgrad
• Gefälle der Einzugsflächen
• Bodenaufbau
Inspiriert durch eine Methode der Anstalt für Wasserversorgung, Abwasserreinigung und Gewässerschutz (eawag), in welcher der Befestigungsgrad der Einzugsflächen anfangs anhand des Farbtones der einzelnen Pixel von Orthophotos ermittelt wird (FANKHAUSER, 1999), soll in diesem Fall, eine Bestimmung des Befestigungsgrades nur durch den normalized vegetation index (NDVI) ermittelt werden.
Zur Berechnung des Oberflächenabflusses, sind unter anderem das Gefälle und der Bodenaufbau wichtige Parameter (ZILCH, et al., 2002). Diese werden aus einem Geländemodell, anhand eines Laserscanns bzw. aus der digitalen Bodenkarte abgeleitet.
Damit die Ergebnisse in die Kanaldatenbank übernommen werden können, wird ein ASCII-Import benötigt, welcher über eine XML-Vorlage definiert ist.
1.5 Werkzeuge
Da der Auftraggeber mit einer BASYS-Kanaldatenbank arbeitet, wird für die hydrodynamische Kanalnetzberechnung das Programm HydroCAD verwendet.
HydroCad ist ein in BASYS eingebettetes Berechnungsmodul und erspart dadurch die Verwendung eines externen Programmes.
Die Analysen der Geodaten erfolgt mit Open Source Produkten. Zur Verwendung kommen QGIS, gvSIG mit der Erweiterung Sextante und SAGA.
1.6 Erwartete Ergebnisse
Folgende Fragen sollen durch die Thesis beantwortet werden:
• In welchen Bereichen können GIS gestützte Analysen sinnvoll für die hydrodynamische Kanalnetzberechnung eingesetzt werden?
• Welche Möglichkeiten gibt es, um befestigte Flächen aus Orthophotos zu ermitteln?
• Liefert die Berechnung des Befestigungsgrades anhand des NDVI ausreichend genaue Ergebnisse für die Kanalnetzberechnung?
• Wie funktioniert der Datentransfer zwischen GIS und Kanaldatenbank?
• Kommt es durch den Einsatz der gewählten Analysen zu einer Zeit- und Kostenersparnis?
1.7 Zielgruppe
Die Thesis richtet sich an Ingenieurbüros, welche, mit einfachen Mitteln, GIS sinnvoll in der Kanalnetzberechnung einsetzen wollen.
1.8 Aufbau der Thesis
Abbildung 1: Aufbau der Thesis
1.
• Allgemeine Einleitung
2.
• Theoretische Grundlagen
3.
• Notwendige Grundlagen für die Berechnung
4.
• Berechnung der Parameter
5.
• Analyse der Ergebnisse
6.
• Zusammenfassung
2. Theoretische Grundlagen
2.1 Kanalnetzberechnung
In den letzten Jahren wurde die Wichtigkeit der Kanalnetzberechnung immer deutlicher.
Die Versieglung von immer mehr Grünflächen führt zu einer erhöhten Wasserableitung in die Kanalisation, dies wiederum führt bei Starkregenereignissen zu einer Überlastung der Kanäle. Überflutungen sind die Folge.
Kanalnetzberechnungen dienen zur Feststellung der Leistungsfähigkeit von bestehenden Kanälen sowie zur Dimensionierung von neuen Kanalsystemen. Das ÖWAV Regelblatt 11 unterscheidet zwischen einfachen Modellen, wie dem Fließzeitverfahren, und komplexen Modellen wie dem hydrologischen oder hydrodynamischen Modell. Es sollten jedoch einfache Modelle nur noch für die Neubemessung von kleinen Kanalnetzen verwendet werden. Hydrologische Modelle sollen zur Berechnung von Schmutzfrachtsimulationen, und hydrodynamische Modelle für den Nachweis größerer Kanalnetze verwendet werden (ATV-M 165, 1994).
Die Vorteile der hydrodynamischen Kanalnetzberechnung umfassen die Berücksichtigung von Druckabfluss, Fließrichtungsumkehr, Verzweigungen, Sonderbauwerken, gesteuerte Elemente (Schieber, Pumpen), sowie des Speichervolumens von Kanalschächten und die Miteinbeziehung des Oberflächenabflusses. Die hydrodynamischen Berechnungsmodelle werden auf der Grundlage der Saint-Venant-Differentialgleichungen durchgeführt, welche das Ablaufverhalten von Abwasserwellen nummerisch lösen (ZILCH, et al., 2002). Eine Erläuterung zu den Gleichungen findet man in der Fachliteratur wie z.B. bei Engel (2006) bzw. Hochedlinger et al (2012).
2.1.1 Ablauf einer Kanalmodellierung
Bei der Kanalnetzberechnung gibt es eine gewisse Reihenfolge zu beachten.
HOCHEDLINGER et al (2012) zum Beispiel, unterteilen die Kanalmodellierung in folgende Teilbereiche:
• Pre-Processing
• Kalibration - Verifikation
• Ingenieurmäßige Interpretation
• Post-Processing
Abbildung 2: Ablaufprozess einer Kanalmodellierung aus (HOCHEDLINGER, et al., 2012)
Der erste Schritt einer Kanalnetzberechnung kann als Pre-Processing bezeichnet werden. In diesem Prozess werden alle wesentlichen Stammdaten der Kanäle (zB.
Dimension, Material, Länge, Gefälle usw.) in das Berechnungsprogramm eingegeben bzw. importiert. Es wird auch überprüft, ob das Kanalmodell berechnet werden kann.
So ist es etwa möglich, dass komplexe Bauwerke nicht berechnet werden können. Hier muss man sich durch das Einfügen von fiktiven Anlageteilen ein sogenanntes Ersatzsystem erschaffen.
Damit man bei der Kanalnetzberechnung seriöse Ergebnisse erzielt, ist es unumgänglich eine Kalibrierung des Kanalnetzes vorzunehmen (HOCHEDLINGER, et al., 2012). Für die Kalibrierung des Kanalnetzes müssen Durchfluss- und Niederschlagsmessungen über den gleichen Zeitraum durchgeführt werden. Anhand des Niederschlages erfolgt eine Berechnung des Kanalnetzes. Im optimal Fall stimmen die berechneten Ergebnissen mit des Ergebnissen der Durchflussmessungen überein. Meistens muss das Kanalnetz noch adaptiert werden, wobei die Schieberöffnungszeiten oder der Befestigungsgrad so geändert werden, bis das Ergebnis der Simulation mit der Durchflussmessung annähernd übereinstimmt.
Abbildung 3: Gegenüberstellung Messdaten und Simulation (aus dem Thesis Projekt)
Ebenfalls in den Bereich der Kalibrierung-Verifikation fällt die Parametervariation. Hier werden den Einzugsgebieten die entsprechenden Umweltparameter zugeordnet. Die Parametervariation ist auch jener Bereich, für den in der Thesis versucht wird, mögliche Verbesserungen mit Hilfe der Geoinformatik zu finden. Dieser Aspekt wird nachstehend näher erläutert.
Nachdem die Kalibrierung abgeschlossen ist, kann mit der eigentlichen hydrodynamischen Kanalnetzberechnung begonnen werden, dem Processing. In diesem Prozess wird eine Kanalnetzberechnung anhand von statistischen Starkregenereignissen durchgeführt. Die Wiederkehrzeit des gewählten Bemessungsregen ist von der Gebietsnutzung abhängig, und wird in der ÖNORM EN 752 (2008) geregelt. Die gelieferten Ergebnisse werden ingenieursmäßig auf Plausibilität geprüft. Dies bedeutet, dass die berechneten Ergebnisse mit den tatsächlichen Ereignissen (Überflutungen usw.) verglichen werden (HOCHEDLINGER, et al., 2012).
Am Ende der hydrodynamischen Berechnung steht das Post-Processing. Es geht nun darum, die Ergebnisse visuell darzustellen. Dies kann durch Überstaupläne oder Belastungslängenschnitte erfolgen. Daraus werden mögliche Schwachstellen im Kanalnetz ersichtlich.
2.1.2 Einzugsgebiete
Ein Einzugsgebiet ist das Gebiet, welches einen Abfluss zu einer Abwasserhaltung, einem Abwasserkanal, oder einem Gewässer hat (ÖNORM EN 752, 2008). Für die genauere Darstellung der Realität unterteilt man ein Einzugsgebiet in mehrere kleine Einzugsflächen, welche bei der Kanalnetzberechnung den Haltungen zugewiesen werden, in die sie entwässern.
2.1.3 Oberflächenabflussmodellierung
Die auf den Boden auftretende Niederschlagsmenge teilt sich in Verdunstung, Versickerung und den Oberflächenabfluss (welcher ein wesentlicher Punkt in der Parametrisierung ist) auf. Unter Oberflächenabfluss versteht man den Teil des Niederschlages, der an der Oberfläche des Bodens abgeleitet wird. Er setzt dann ein, wenn die Niederschlagsmenge die Versickerungsrate des Bodens übersteigt. Die Höhe des Oberflächenabflusses hängt, neben der Niederschlagsmenge, noch von der Lage (Hangneigung, Hangform) und dem Aufbau des Bodens ab (Humusgehalt des Oberbodens, Porenverhältnisse usw.) (WOHLRAB, 1992). Ein weiterer wesentlicher Parameter für die Berechnung des Oberflächenabflusses ist der Befestigungsgrad (PRAXL, 2010). Er kann basierend auf Luftbildaufnahmen, Naturbestandsaufnahmen oder ähnlichen genauen Grundlagen bestimmt werden (ÖWAV Regelblatt 11, 2009).
Als befestigte Flächen gelten, undurchlässige (Asphaltflächen, Ziegel-, Metall-, Glasdächer) oder unterschiedlich durchlässige Oberflächen (Schotter) (ÖWAV Regelblatt 11, 2009) . Weitere benötigte Parameter sind die Hangneigung, der Muldenverlust, die Dachneigungen, der Benetzungsverlust und der Bodenaufbau. Alle erforderlichen Parameter, für die Berechnung des Oberflächenabflussmodells, können bei jedem Einzugsgebiet unterschiedlich sein, und müssen daher einzeln zugeordnet werden.
Für die Berechnung des Oberflächenabflusses gibt es zwei häufig verwendete Methoden. Dies sind das hydrologische und das hydrodynamische Berechnungsmodell.
Das hydrodynamischen Modell ist eine implizite Berechnung, welche eine längere Berechnungszeit mit sich trägt. Die Berechnungsergebnisse sind dafür realitätsgetreuer.
Bei dem hydrologischen Modell handelt es sich um eine explizites Berechnungsmodell, mit dem Vorteil von kürzeren Rechenzeiten (HOCHEDLINGER, et al., 2012). Ein weiterer wesentlicher Unterschied zwischen den Modellen ist, dass bei einem
hydrodynamischen Modell ein Rückstau (das Wasser kann nicht ungehindert von oben nach unten abfließen) berücksichtigt werden kann (VEIT, 2009),
In der verwendeten Software wird das hydrodynamische Modell verwendet, wobei der Abfluss nach den Gleichungen von Horton berechnet wird. Diese Gleichung gibt den Abfluss bei einem Starkregenereignis an. Die maximale Wasserkapazität der Bodensäule wird dabei noch nicht erreicht. Die, auf Grund des intensiven Regens, vorhandenen Wassermengen können aber trotzdem nicht mehr in den Boden infiltrieren.
Eine Erläuterung dieser Formeln ist in ILLGEN(2009) nachzulesen.
2.1.4 Berechnungssoftware 2.1.4.1 BASYS - HydroCad
BASYS ist ein serverbasiertes Netzinformationssystem der Firma Barthauer. Kern des Programmes ist eine Datenbank, die entweder auf Oracle oder Microsoft SQL Server bereitgestellt wird. In dieser Datenbank sind alle Angaben zu den Objekten des Netzes gespeichert. BASYS verfügt über zahlreich Schnittstellen (ISYBAU, ASCII usw.), welche alle über XML-Schemas definiert sind, und beliebig bearbeitet und erweitert werden können.
In BASYS integriert ist das HydroCAD-Modul der Firma Dorsch Consult. Mit HydroCAD wird die hydrodynamische Kanalnetzberechnung durchgeführt. Dies geschieht auf Basis der impliziten Ganglinien-Volumen-Methode (GVM) (BARTHAUER, 2013). Der Oberflächenabfluss wird hier ebenfalls hydrodynamisch gerechnet.
2.1.4.2 HYSTEM-EXTRAN
Als eine alternative zu HydroCad sei hier HYSTEM-EXTRAN erwähnt. Dies ist ein, in der Praxis oft verwendetes Programm zur hydrodynamischen Kanalnetzberechnung.
HYSTEM-EXTRAN wurde vom Institut für technisch-wissenschaftliche Hydrologie (itwh) in Deutschland entwickelt. Es setzt sich zusammen aus dem Programmteilen HYSTEM, zur hydrologischen Oberflächenabflussberechnung, und EXTRAN zur hydrodynamischen Kanalnetzberechnung (Abwasser, 2014).
2.1.5 Hydrodynamische Kanalnetzberechnung mit der GVM
Die GVM ist ein Niederschlag-Abfluss-Modell zur Berechnung des Abflusses in ein offenes oder geschlossenes Gerinne mit einem impliziten Differenzenschema. Bei dieser Methode wird der Weg des Niederschlages vom Auftreffen am Boden bis zur Ableitung durch die Kanalisation wirklichkeitsgetreu wiedergegeben. Bei der GVM erfolgt eine vollständige implizite Lösung der Saint-Venant-Gleichungen. Vorteile des impliziten Verfahren sind unter anderem, dass das Volumen der Kanalschächte in die Berechnungen mit einfließt, dass die Drucklinien über die Geländeoberkante ansteigen können, und dass das übergelaufene Wasser wieder in das Kanalnetz zurückfließen kann. Ein bereits erwähnter Nachteil sind die längeren Rechenzeiten.
2.2 Geodaten 2.2.1 Bodenkarte
In Österreich wurde im Auftrag des Bundesministerium für Land. und Forstwirtschaft eine digitale Bodenkarte erstellt. Diese Bodenkarte ist im Wesentlichen eine Übernahme von grafischen und textlichen Informationen der seit 1958 durchgeführten Bodenkartierung (Lebensminsiterium, 2013). Die Karte, im Maßstab 1:25.000, ist unterteilt in 219 Kartierungsbereiche, und enthält 10500 Bodenformen. Der große Maßstab bildet eine gute Übersicht über die Bodenformen, ist aber für Details nur bedingt geeignet (MURER, 2009). Die Kartierung erfolgt anhand von Begehungen vor Ort, Bohrproben und Schürfungen bis zu einer Tiefe von 1,5m.
Die detaillierte Beschreibung über den Aufbau der einzelnen Bodenformen findet man auf der Homepage über die Bodenkarte (Bundesforschungs- und Ausbildungszentrum für Wald, 2014). Die Bodenkarte kann zu einem Fixpreis pro km² erworben, oder über die kostenlose Web-GIS Applikation eBOD genutzt werden.
Die Bodenkarte enthält alle, für die Berechnung des Oberflächenabflusses, notwendigen Informationen über die Bodenbeschaffenheit.
2.2.2 Fernerkundungsdaten 2.2.2.1 Satellitenbilder
Rund um die Erde kreisen zahlreiche Satelliten, mit der Aufgabe die Erdoberfläche zu beobachten. Damit alle Punkte auf der Erde aufgenommen werden können, wurde für die Satelliten eine polnahe kreisförmige Umlaufbahn gewählt (ALBERTZ, 2009).
Vorreiter bei der Luftbildaufnahme durch Satelliten war das Landsat-Programm.
Aktuell steht mit der Pléiades Satellitengeneration eine hoch automatisierte Möglichkeit zur Landbeobachtung zur Verfügung. Diese Satelliten liefern ein täglich aktuelles Bild eines Punktes auf der Erdoberfläche, mit einer Auflösung von bis zu 50cm im panchromatischen Kanal (STEINOCHER, et al., 2014). Für die technische Ausrüstung der Satelliten wird auf die Fachliteratur verwiesen.
2.2.2.2 Orthophotos
Ein anderes Verfahren der Landbeobachtung sind digitale Laufaufnahmen anhand von Befliegungen. Hier nimmt eine, an einem Flugzeug befestigte, Kamera die Oberfläche auf. Durch die Zentralprojektion, die unterschiedlichen Entfernungen der Objekte zur Kamera, und den Schieflagen des Flugzeuges entstehen Verzerrungen des Luftbildes.
Insbesondere durch die Zentralprojektion kommt es, bei Objekten die über den Boden ragen, zum Beispiel Bäume und Häuser, zu starken Verszerrungen. Diese Verzerrungen werden in der Fachliteratur als radiale Versetzungen bezeichnet, und nehmen am Rand der Bilder stark zu (ALBERTZ, 2009).
Abbildung 4: Radiale Versetzung durch Höhenunterschiede aus (ALBERTZ, 2009)
Ein Orthophoto ist ein entzerrtes Luftbild. Die Verzerrungen werden durch den Vorgang
2.2.2.3 Auswerten der digitalen Luftaufnahmen
Eine Möglichkeit digitale Bilder auszuwerten, ist das Verfahren der multispektralen Klassifizierung (auch pixelbasierte Klassifizierung). Die multispektrale Klassifizierung basiert auf den spektralen Reflektionswerten der Objekte, die in verschiedenen Spektralkanälen unterschiedlicher spektraler Empfindlichkeit aufgezeichnet werden (ULRICH, 2011). Man unterscheidet zwischen der überwachten und der unüberwachten Klassifizierung. Eine unüberwachten Klassifizierung ist die Cluster-Analyse. Bei dieser Methode werden keine Testgebiete benötigt, daher sind die Ergebnisse meist unzufriedenstellend (ALBERTZ, 2009) und können daher nur für grobe Analysen verwendet werden. Als überwachte Klassifizierungen gelten das Minimum-Distance- Verfahren, das Parallelepiped-Verfahren das Maximum-Likelihood-Verfahren und die hierarchische Klassifikation. Ein Nachteil der pixelbasierten Verfahren ist, dass es bei hoher Auflösung der Luftbilder zu dem sogenannten "salt & pepper" Effekt kommt (BLASCHKE, 2000).
Dem gegenüber steht das Verfahren der objektbasierten Klassifizierung (OBIA).
Hierbei werden nicht die einzelnen Pixel analysiert, sondern ein Regionen mit ähnlichen Pixel. Bevor mit der Klassifizierung begonnen werden kann, muss ein Segmentierungsschritt vorgeschaltet werden. Dabei werden die Pixel aufgrund bestimmter Eigenschaften zu Segmenten mit ähnlichen Pixeln zusammengefasst (ALBERTZ, 2009). Diese Eigenschaften können Form, Farbe, Höhe, Nachbarschaftseigenschaften, oder auch der "noramlized difference vegetation index"
(NDVI) sein. Einen sehr guten Überblick über Projekte wo OBIA erfolgreich eingesetzt wird bietet BLASCHKE (2010).
2.2.2.4 NDVI
Der NDVI ist ein häufig verwendeter Vegetationsindex. Er beruht darauf, dass jedes Objekt auf der Erdoberfläche Lichtstrahlen unterschiedlich reflektiert bzw. absorbiert.
Dies ist abhängig von der Zusammensetzung des Objektes, sowie von der Wellenlänge der betroffenen Strahlung. Für die Vegetation bedeutet dies zum Beispiel, dass bei photosynthetisch aktiver Vegetation eine starke Absorption im sichtbaren rotem Bereich, und eine hohe Reflexion im nahen Infrarot vorhanden ist. Bei Böden hingegen gibt es diesen Unterschied zwischen den Kanälen nicht. Man kann daher den NDVI zur Bestimmung für die Vitalität der Vegetation heranziehen, und somit Klassifizierungen
von unterschiedlichen Vegetationstypen, oder zur Unterscheidung von Vegetation / keine Vegetation heranziehen. Des Weiteren kann der NDVI zur Beurteilung des Wassergehaltes von Pflanzen verwendet werden. Der NDVI basiert auf einer Ratio- Berechnung, bei der die Spektralwerte des nahen Infrarot (NIR), welche die Vitalität der Vegetation besonders widerspiegeln, mit den Werten des roten Spektralbereichs (RED), auf die sich der Zustand der Vegetation nicht auswirkt, in Beziehung setzt (ALBERTZ, 2009):
NDVI = (NIR – RED) / (NIR + RED)
Das Ergebnis führt zu Werten, die ausschließlich zwischen –1 und +1 liegen, wobei nahe bei +1 liegende Werte charakteristisch für gesunde, photosynthetisch aktive Vegetation sind, und nahe bei -1 liegende Werte kennzeichnend für Wasser sind (ULRICH, 2011). Abgestorbene Vegetation sowie vegetationslose Bereiche weisen niedrige Werte im Bereich von Null auf. Ein Vorteil des NDVI ist, dass Unterschiede in den Beleuchtungsverhältnissen, wie z.B. durch Schlagschatten, und Einflüsse der Geländeneigung durch die dargestellte Ratiobildung weitgehend kompensiert werden können (ALBERTZ, 2009).
2.2.2.5 Laserscanmessungen
Laserscanning ist ein aktives Messverfahren, bei dem die Entfernung zwischen Laser und Objekt gemessen wird. Airborne Laserscaning (ALS) ist ein flugzeuggestütztes Messverfahren und erlaubt eine präzise Aufnahme des Geländes. Bei dieser Methode ist ein Scanner am Flugzeug oder Helikopter montiert, und sendet in regelmäßigen Abständen Laserstrahlen zur Erdoberfläche. Aus der Zeit die vergeht vom Absenden des Signals, bis zum wieder eintreffen am Scanner wird die Entfernung berechnet (KRAUS, 2004).
Abbildung 5: Abtastender Laser im Flugzeug, aus (KRAUS, 2004)
Zur genauen Positionierung des Scanners, und somit auch zur koordinativen Bestimmung der gemessenen Punkte, wird ein Satellitenpositionierungssystem (GPS) und eine inertiale Messeinrichtung (IMU) benötigt. Beide sind ständig mit dem Bordcomputer des Flugzeuges verbunden (KRAUS, 2004). Laut WEVER (1999) hat der gesendete Laserimpuls an der Oberfläche einen Durchmesser von etwa 25cm, das führt dazu, dass der Impuls öfter reflektiert wird. Der als erstes reflektierte Punkt, First- Pulse, (zB.: Baumkronen oder Hausdächer) wird für das digitale Oberflächenmodell (DOM) verwendet. Der Impuls, der als letzter (Last-Pulse) aufgezeichnet wird, spiegelt den Boden wieder und wird für das DGM verwendet.
Abbildung 6: First- und Last-Pulse Echo eines Laserstrahls (www.toposys.de, 2014)
2.2.3 LISA
Das Österreichische Landinformationssystem (LISA) ist ein Projekt, in dem die oben erwähnten Analysemethoden erfolgreich kombiniert und angewendet werden.
LISA wurde 2009 gestartet und ist ein, mit öffentlichen Mitteln finanziertes, Projekt zur Darstellung der Landnutzung und Landbedeckung in Österreich. In der ersten Phase des Projektes wurde die Landbedeckung mit einer objektbasierten Klassifizierung durchgeführt. Durch die Hinzunahme von LIDAR-Daten konnte in einer zweiten Phase auch die Landbedeckung ermittelt werden (GRILLMAYER, et al., 2010). Eine genau technische Beschreibung des Projektes ist nachzulesen in BANKO, et al. (2010).
Derzeit sind die Daten für 16 Testgebiete vorhanden und können auf www.landinformationssystem.at heruntergeladen werden. Wann es für Österreich einen flächendeckenden Datensatz geben wird ist aus finanziellen Gründen noch unklar. Es sollen jedoch bis September 2014 ein Großteil der österreichischen Ballungszentren kartiert werden (STEMBERGER, 2013).
2.3 Technische Grundlagen 2.3.1 Open Source GIS
In der Masterthesis wird versucht, die anfallenden Analysen mit Open Source GIS Programmen zu lösen. Open Source Produkte sind Programme, die kostenlos zur Verfügung gestellt, und nur durch Sponsoren finanziert werden.
Es werden im Wesentlichen drei Programme zur Anwendung kommen: Quantum GIS (QGIS), SAGA und gvSIG. Alle drei Programme sind unter der GNU General Public License (GPL) verfügbar und sind Projekte der Open Source Geospatil Foundation (OSGeo) (MITCHELL & CHRISTL, 2008). Die Softwareprodukte sind von ihren Raster- und Vektorfunktionalitäten nahezu ident, unterscheiden sich aber von der Benutzeroberfläche doch einigermaßen. Die SEXTANTE Erweiterung erlaubt es, dass Befehle oberflächenunabhängig verwendet werden können. Das bedeutet, dass z.B.
SAGA Befehle in QGIS funktionieren.
2.3.2 XML
Die Extensible Markup Language (XML) ist eine Auszeichnungssprache. Sie wurde ursprünglich erfunden um Daten zu speichern, und nicht um Daten anzuzeigen.
Mittlerweile ist XML aber überall zu finden und dient als meist genütztes Tool zur Übertragung von Daten zwischen diversen Anwendungen (W3, 2013). XML ist eine Metasprache, mit der konkrete Sprachen definiert werden können. Damit sind alle XML-Sprachen aus den gleichen Konstrukt aufgebaut. Ein Programm, das die XML- Regeln kennt, kann jede XML-Sprache lesen und in ihre Bestandteile zerlegen (Parser).
Eine wesentliche Eigenschaft von XML-Sprachen ist, dass sie selbst beschreibend sind (FISCHER & HUBER, 2013).
Ein XML-Dokument ist hierarchisch aufgebaut. Elemente können geschachtelt sein, womit die hierarchische Abhängigkeit zwischen Elementen dargestellt wird. Wichtig ist, das Anfangs- und End-Tag richtig verschachtelt sind, wobei das End-Tag des ersten Elementes immer am Schluss kommt (FISCHER & HUBER, 2013). Die Tags werden in spitzen Klammern geschrieben.
In der Thesis wird XML unter anderem als Schnittstelle zwischen den einzelnen Software Produkten verwendet.
3. Grundlagen für die Analysen
3.1 UntersuchungsgebietWie bei Punkt 1.1 bereits erwähnt, handelt es sich bei der Kanalnetzüberrechnung um Kanalstränge des RHV Schwanenstadt. Die Stadt Schwanenstadt befindet sich im Nord- Osten des Bezirkes Vöcklabruck. Sie liegt auf einer Seehöhe von ca. 390 müA und hat eine Fläche von 2.6 km². Von der Einwohnerdichte mit 1588.8 EW/km² lässt sich für Schwanenstadt ein eher urbanes Siedlungsgebiet ableiten.
Abbildung 7: Karte von Oberösterreich Datenquelle: (DORIS, 2013)
Die zu berechnenden Kanalstränge (VBOB01 und VBNT01) befinden sich im Norden und im Süden des Stadtgebietes. In das Kanalnetz der zwei Stränge münden insgesamt 72 Einzugsflächen mit einer Gesamtfläche von 0,43km².
Bezirk Vöcklabruck
Schwanenstadt
Abbildung 8: Einzugsflächen von Schwanenstadt
3.2 Erforderliche Daten in der Berechnungssoftware 3.2.1 Einzugsflächen
Damit eine hydrodynamische Kanalnetzberechnung durchgeführt werden kann, müssen in der verwendeten Software (HydroCad) den Einzugsflächen folgende Parameter zwingend zugewiesen werden:
• Bezeichnung
• Oberflächentyp
• Befestigungsgrad
• Gesamtfläche
• Haltung
Die Bezeichnung der Einzugsflächen ist aus einer älteren Überrechnung vorgegeben und wird in weiterer Folge als Verknüpfungsfeld für den ASCII-Import verwendet. Der Oberflächentyp wird in einem eigenen Eingabeformular definiert. Die Zuweisung zu den jeweiligen Einzugsflächen erfolgt im GIS-Programm und wird mittels ASCII- Import importiert. Der Befestigungsgrad wird ebenfalls im GIS ermittelt und importiert. Die Größe der befestigten Flächen ergibt sich aus den Voreinstellungen des
Oberflächentyps. Der Oberflächentyp und der Befestigungsgrad sind essentielle Bestandteile der Oberflächenabflussberechnung. Die Gesamtfläche kann entweder manuell eingegeben, oder rechnerisch aus einem digitalen Umriss ermittelt werden. Da für die GIS Analysen ein Polygon der Einzugsflächen vorhanden sein muss, benötigt man eine "shp" bzw. eine "dxf" Datei. Diese Datei wird aus der Kanaldatenbank mit Hilfe des Barthauer Moduls BASYS-Plan A erzeugt. Die Zuweisung der Haltung zur Einzugsfläche erfolgt manuell in einem speziellen Einzugsflächenformular.
Abbildung 9: HydroCad Einzugsflächenformular
3.2.2 Oberflächentypen
Alle für den Oberflächenabfluss relevanten Parameter werden in dem Formular
"Oberflächentyp" eingegeben. Bevor mit den Analysen in einer GIS Anwendung begonnen wird, müssen hier die Oberflächentypen definiert werden. Für das Untersuchungsgebiet wurden 18 verschiedene Oberflächentypen definiert:
• STR NG I, NG II, NG III, NG IV
• STZ LS NG I/II, NG III/IV
• STZ L NG I/II, NG III/IV
• STZ KS NG I/II, NG III/IV
• WGB L NG I/II, NG III/IV
• WGB LS NG I/II, NG I/II Acker, NG III/IV
• WGB SK NG I/II, NG III/IV
Die Bezeichnung der Oberflächentypen ergibt sich aus der Lage, den Bodenverhältnissen und der Hangneigung der Einzugsflächen. STR steht für Straße, STZ für Stadtzentrum und WGB für Wohngebiet. An zweiter Stelle stehen die Bodenverhältnisse, hier bedeutet LS Lehmiger Sand, L steht für Lehm und SK für Sandiger Kies. Die Unterteilung nach den Bodenverhältnisse könnte noch mehrere Untergruppen aufweisen, in Schwanenstadt sind jedoch, lt. digitaler Bodenkarte, nicht mehr unterschiedliche Bodenverhältnisse vorhanden. An der dritten Stelle der Bezeichnung befindet sich die Neigungsklasse. Um eine Vielzahl an verschiedenen Oberflächentypen zu vermeiden, wurden immer zwei Neigungsgruppen zusammengefasst. Die Gruppen NG I/II beinhalten Neigungen kleiner 4% Gefälle, und die NG III/IV die Flächen mit steileren Gelände.
Bei jedem Oberflächentyp müssen für die befestigen Flächen (Verkehrsflächen und Dachflächen) und den nicht befestigten Flächen noch folgende Parameter eingegeben werden:
Rezession K, Start- und Endwert nach Horton sowie die Rauigkeit werden aus einer Tabelle mit Erfahrungsberichten (Anhang 1) übernommen. Diese 4 Parameter sind von den geologischen Bodenbeschaffenheit abhängig.
Im Feld Neigung wird der Mittelwert der betreffenden Neigungsgruppe eingegeben.
Die Streifenbreite wird in HydroCad zu Berechnung der flächenmäßigen Anteile von Dachflächen und Verkehrsflächen benötigt.
Der Muldenverlust gibt den Verlust des Niederschlages an, welcher in Mulden bzw.
Geländeunebenheiten zu Stande kommt.
Unter Benetzungsverlust versteht man wiederum einen Verlust des Niederschlages, welcher entsteht, bis eine Fläche vollkommen benässt ist.
Abbildung 10: HydroCad Oberflächentypformular
3.3 Möglichkeiten für den Einsatz von GIS
Anhand der von HydroCad geforderten Parameter für die hydrodynamische Kanalnetzberechnung, lässt sich erkennen, wo Geoinformationssysteme unterstützend eingesetzt werden können. Dies sind auf Grund ihrer geometrischen Eigenschaften
• der Befestigungsgrad
• die Zuweisung des Oberflächentyps anhand o der Bodenverhältnisse und
o der Hangneigung
Ein weiterer Ansatz wäre, das Abflussverhalten der Einzugsflächen mittels GIS Analysen (Flow accumulation) zu bestimmen, und dadurch den Einmündungspunkt der Oberflächenwässer in die Kanalisation zu definieren. Da aber bereits durch die Kamerabefahrung alle Anschlüsse erhoben wurden, wird diese Möglichkeit in der Thesis nicht behandelt.
3.4 Auswirkungen der gewählten Parameter
Wie in den folgenden Diagrammen Dargestellt, wirken sich die oben genannten
Für die Diagramme wurde ein 10 jährliches Regenereignis und eine Einzugsfläche von ca. 2ha angenommen.
Abbildung 11 soll verdeutlichen, wie sehr sich der Befestigungsgrad auf den Oberflächenabfluss auswirkt.
Abbildung 11: Auswirkungen des Befestigungsgrades
Werden bei 10% Befestigungsgrad nur ca. 40l/s weitergeleitet, so sind es bei 50%
bereits das Dreifache und bei 90% das Sechsfache des Anfangswertes.
Bei Abbildung 12 wird ersichtlich, wie sich Bodenart und Gefälle bei gleichem Befestigungsgrad (10%) auf den Oberflächenabfluss auswirken.
Abbildung 12: Zusammenspiel von Boden und Gefälle
Es wird ersichtlich, dass Kiesböden um ca. 1/5 weniger Wasser in die Kanalisation weiterleiten als Lehmböden. Einen deutlichen Unterschied des Oberflächenabflusses erkennt man auch bei der Neigung der Flächen. So geben etwa steile Hänge um ca.80%
mehr Wasser ins Kanalnetz ab als ebene Flächen.
3.5 Datensichtung 3.5.1 Kanalnetz
Das Kanalnetz der zu berechnenden Kanalabschnitte ist in einer BASYS-Datenbank vorhanden. Alle erforderlichen Stammdaten, wie Dimension, Sohlhöhen, Längen und der Gleichen sind darin eingegeben. Das Kanalnetz kann als shp-Datei oder im AutoCad dwg-Format ausgegeben werden.
3.6.2 Einzugsflächen
Die Einzugsflächen sind aus einer älteren Überrechnung übernommen worden, und nur geringfügig verändert bzw. ergänzt worden. Sie wurden mittels BASYS-Plan A Modul in die Datenbank eingearbeitet und können ebenfalls als dwg oder shp exportiert werden.
3.6.3 Luftbilder
Da vom Auftraggeber Orthophotos aus dem Jahr 2013 für dieses Projekt zur Verfügung gestellt wurden, wurde auf Analysen an Hand von Sattelitenbilder verzichtet.
In Österreich werden seit 1998 Orthophotos erstellt. Seit 2010 werden diese in einem regelmäßigen Abstand von drei Jahren, von der gesamten Landesfläche erzeugt. Durch die Einführung von Bildmessflügen wurde 2005 eine Verbesserung der geometrischen Bodenauflösung erreicht. Die Auflösung der Orthophotos beträgt nun 20cm. Seit 2009 werden die Orthophotos als Echtfarbbilder (RGB) und Farbinfrarotbilder (CIR) geliefert (DORIS, 2013), und können daher auch für die Berechnung des NDVI verwendet werden. Die Orthophotos von Oberösterreich können über das Online Portal des BEV (www.bev.gv.at), oder über die Abteilung für Geoinformation der oberösterreichischen Landesregierung bezogen werden. Die Infrarotbilder, welche für die Berechnung des NDVI notwendig sind, können nicht über das Online Portal bestellt werden, sondern müssen schriftlich bestellt werden.
Die für das Untersuchungsgebiet notwendigen Orthophotos sind: 4933-5002, -5003, -5200 und -5201. Die Orthophotos wurden vom Auftraggeber zur Verfügung gestellt, die Falschfarbbilder wurden vom Büro Hitzfelder & Pillichshammer erworben.
3.6.4 Digitales Geländemodell
In Oberösterreich ist für die Erstellung von Höhenmodellen das Landesverwaltung zuständig. Seit 2011 gibt es von ganz Oberösterreich Laserscandaten, welche eine Punktdichte von ca. 16Pkt./m² und eine Höhengenauigkeit von 5-10cm aufweisen. Auf Grund eines Kooperationsvertrages zwischen den Gemeinden und der OÖ- Landesregierung werden den Kommunen 1m Höhenschichtlinien bzw. 10m Raster Punkte kostenlos zur Verfügung gestellt. Ein 1m Punktraster kann gebührenpflichtig bei der Abteilung für Geoinformation des Landes Oberösterreich angefordert werden (GREIFENEDER, 2013).
3.6.5 Bodenkarte
Die digitale Bodenkarte im Untersuchungsgebiet wurde vom Auftraggeber erworben, und steht als shp-Datei zur Verfügung. Die Bodenkarte enthält im Stadtzentrum keine Angaben über den Bodenaufbau. Die notwendigen Angaben müssen daher von Erfahrungen aus Tiefbauprojekten abgeleitet werden.
4. Ermittlung der Parameter
4.1 Berechnung des Befestigungsgrades
Wie bereits in den vorherigen Kapiteln öfters erwähnt, ist der Befestigungsgrad eine der wichtigsten Voreinstellungen für die hydrodynamische Kanalnetzberechnung. Für die Bestimmung dieses Parameters stehen daher auch einige Möglichkeiten und Methoden zur Verfügung.
Zu Zeiten in denen Orthophotos und Satellitenbilder noch nicht allgegenwärtig waren, wurde der Befestigungsgrad durch Vor-Ort-Begehung und Vermessungen ermittelt.
Diese Methode ist wahrscheinlich die genaueste Methode, wird heute aber aus Zeit- und Kostengründen jedoch nur noch in seltenen Fällen verwendet. Mittlerweile stehen, von fast jedem Punkt der Erdoberfläche, Luftbilder zur Verfügung. Diese stellen eine Ideale Grundlage für die Bestimmung des Befestigungsgrades dar. Eine sehr einfache Methode dafür ist die manuelle Digitalisierung. Hierbei werden die befestigten Flächen manuell mittels CAD oder GIS Programm über die Luftbilder gezeichnet werden. Diese Methode ist für ein kleines Einzugsgebiet mit nur wenigen kleinen Einzugsflächen durchaus empfehlenswert. Für größere Einzugsgebiete ist diese Vorgehensweise jedoch eine sehr zeitaufwändige, und dadurch auch kostenintensive Lösung.
Eine automatische Ermittlung der befestigten Flächen bieten Analysen mittels OBIA.
Wie bereits in Kapitel 2.2.2.3 erwähnt, werden hier anhand mehrerer gleichen Eigenschaften Objekte gebildet. Mit dieser Methode wurden in zahlreichen Projekten schon sehr genaue Aussagen über die Landnutzung getätigt (BLASCHKE, 2010). Für dieses Verfahren sind jedoch sehr gute GIS-Kenntnisse und Erfahrungen in der Fernerkundung erforderlich. Da diese Kenntnisse nur in den wenigsten Zivilingenieurbüros anzutreffen sind, wird diese Methode nicht in der Thesis angewendet. LISA-Daten können als Ergebnis einer OBIA bezeichnet werden. Diese Daten werden zukünftig eine wesentliche Rolle in der Bestimmung des Befestigungsgrades spielen. Da aber noch nicht genau gesagt werden kann, ab wann eine flächendeckende Datenstruktur vorhanden ist, und wie hoch die Kosten für die Anschaffung sind, werden keine Analysen mit diesen Daten in der Thesis durchgeführt.
Abbildung 13: Landbedeckung eines LISA Testdatensatzes in Oberösterreich (GRILLMAYER, et al., 2010)
In einem Projekt der eawag von 1999 wurde versucht, den Befestigungsgrad der Einzugsflächen mittels Farbton der einzelnen Pixel zu ermitteln. Um die Ergebnisse zu verbessern, wurde zu den Farbtönen noch der NDVI hinzugefügt. Das Ergebnis brachte eine Differenz zwischen den von Hand digitalisierten Flächen und den automatisierten Flächen von kleiner 10% (FANKHAUSER, 1999). Aufgrund des hohen Aufwandes für die Klassifizierung durch die Farbtöne, wird in der Thesis versucht, den Befestigungsgrad nur anhand des NDVI zu bestimmen. Den Grundgedanken hierzu bildet die Annahme, dass Oberflächen mit Vegetation sickerfähig sind, und Oberflächen ohne Vegetation versiegelt sind.
4.1.1Ermittlung des NDVI
Wie bereits erläutert, wird bei der Thesis versucht die durchzuführenden Analysen mittels Open Source Software zu bewältigen. Aus diesem Grund werden die einzelnen Funktionen in den Applikationen etwas genauer beschrieben. Für die Ermittlung des NDVI wird die Erweiterung Sextante von der gvSig Oberfläche aus ausgeführt.
Da die Orthophotos eine sehr hohe Auflösung aufweisen (0,2x0,2m), würde ein Zusammenfügen der Bilder zu einem Objekt, die Rechenzeit und die Leistungsfähigkeit der Software und der herkömmlichen Hardware überfordern. Aus diesem Grund, wird der NDVI für jedes Falschfarbbild einzeln berechnet.
Der benötigte "Red Layer" ist in den herkömmlichen Orthophotos enthalten. Für den
"near infrared Layer" wird das dazugehörige Falschfarbfoto ausgewählt. Bei beiden Bildern ist der "Red Band" das Band 1.
Das Ergebnis des NDVI Befehls ist ein Raster mit Werten zwischen -1 und 1. Dieser Raster hat die selbe Auflösung wie die Eingangsparameter.
Abbildung 14: vom Orthophoto zum NDVI-Raster
4.1.2 Grenzwertbestimmung
Da von der Überlegung ausgegangen wird, dass Flächen mit einer gesunden Vegetation versickerungsfähig sind (NDVI nahe 1), und Flächen mit keiner Vegetation befestigt sind (NDVI nahe 0), muss nun ein Wert gefunden werden, der die Grenze am Besten wiederspiegelt. Hierfür wird die Testeinzugsfläche 102a ausgewählt. Bei dieser Einzugsfläche werden nun die befestigten und nicht befestigten Flächen per Hand eingezeichnet. Für diese Fläche wird auch eine, wie in Punkt 4.1.3 beschriebene, Reklassifizierung des NDVI Rasters solange durchgeführt, bis die digitalisierten und berechneten Flächen annähernd übereinstimmen.
Abbildung 15: links: digitalisierte befestigte Flächen für EZGB 102a; rechts: NDVI für EZGB102a
Die Einzugsfläche hat eine Gesamtfläche von 8049m². Die digitalisierten befestigten Flächen (gelb markiert) ergeben eine Fläche von 5300m², was einen Befestigungsgrad von 68,33% bedeutet. Reklassifiziert man den NDVI Raster mit dem Wert von 0.06, so erhält man eine befestigte Fläche von 65%. Die geringe Abweichung von 3% ist tolerierbar. Somit kann der Wert 0.06 für die weiteren Berechnungen verwendet werden.
Probleme zeichnen sich bei Ackerflächen und Straßen ab. Ackerflächen haben keine Vegetation, und daher einen niedrigen NDVI-Wert, sind aber sickerfähig (Abbildung 10). Da Ackerflächen jedoch keinen Kanalanschluss aufweisen, werden sie auch in der Berechnung nicht berücksichtigt. Straßen sind versiegelte Flächen, weisen aber teilweise einen Wert von über 0,1 auf. Es wird daher, bei reinen Straßeneinzugsflächen vorab ein Befestigungsgrad vergeben.
Abbildung 16: NDVI-Wert einer Ackerfläche
4.1.3 Reklassifizierung
Durch die Berechnung im Testgebiet hat sich ein Grenzwert, zwischen befestigten und nicht befestigten Flächen, von 0.06 ergeben. Mit einem "reclassify" Tool wird der NDVI-Raster neu eingeteilt. Um die Rechenzeit zu verringern, wird jedes NDVI-Bild extra klassifiziert. Bei der Reklassifizierung werden die Parameter so eingegeben, dass alle Werte die größer als 0.06 sind (versickerungsfähig), zu "no Data" werden, und alle Werte kleiner 0.06 (keine Versickerung) erhalten den Wert "1".
Bei der Reklassifizierung kann auch die Auflösung des Rasters verändert werden. Da die neu erstellten Raster zu einem gesamten Raster zusammengefügt werden, wird zur Verringerung der Rechenzeit, eine geringe Auflösung (1x1) verwendet. Nach dem alle NDVI-Rasterbilder reklassifiziert wurden, werden sie mittels "merge" Befehl zu einem Raster für das Einzugsgebiet zusammengefügt.
4.1.4 Befestigte Flächen den Einzugsflächen zuordnen
Um herauszufinden wie groß die befestigte Fläche einer Einzugsfläche ist, müssen der Raster und die Einzugsflächen übereinander gelegt werden. In ArcGIS funktioniert dies relativ einfach mit dem "Zonal Statistics as Table" Befehl des "Spatial Analyst" Moduls.
Das Ergebnis des Befehles ist eine Tabelle, mit den Eigenschaften des Rasters bezogen auf die einzelnen Polygone.
Damit man mit einem Open Source Produkt das gleich Ergebnis bekommt, muss man etwas in diversen Foren stöbern, um eine Lösung zu finden. Ein zufriedenstellendes Ergebnis liefert das experimentelle Plugin "ZonalStat" für QGIS. ZonalStat benötigt die Phyton-Erweiterung "scipy" welche im Internet (Scipy-developers, 2013) heruntergeladen werden kann. Das Ergebnis der ZonalStat-Funktion ist eine csv-Datei, in welcher die Anzahl der Zellen mit dem Wert "1" pro Einzugsfläche in der Spalte
"area" eingetragen wird.
Die Zuweisung der befestigten Flächen der jeweiligen Einzugsflächen in das Einzugsflächenshape erfolgt durch einen "join"-Befehl. Ein Join Befehl verbindet zwei verschiedene Datensätze miteinander. Voraussetzung dafür ist, dass Beide ein Tabelle besitzen, in der die gleiche Bezeichnung oder Id eingetragen ist.
Das Zusammenführen der csv Tabelle und des Shapes erfolgt ebenfalls in QGIS. Damit der join-Befehl ausgeführt werden kann, müssen beide Dateien in die Software geladen werden. Als Join-Field dient die Spalte "Bezeichnung". Das Shape enthält jetzt die Summe der befestigen Flächen pro Einzugsfläche, jedoch noch nicht die gesamte Größe der einzelnen Einzugsflächen. Die Ermittlung der Fläche, sowie die anschließende Berechnung des Anteils der befestigten Flächen erfolgt über den QGIS Feldrechner.
Das Feld für die Gesamtfläche erhält den Namen "gesamtFl". Den nötigen Befehl zur Berechnung findet man in der Funktionsliste unter Geometrie. Für die Berechnung des Anteils der befestigten Flächen wird das vorhandene Feld "Abflussbei" erneuert. Der Ausdruck für die Berechnung des Feldes lautet: "area" / "gesamtFl". Die Feldnamen können ebenfalls über die Funktionsliste, unter "Felder und Werte" ausgewählt werden.
4.2 Berechnung der Hangneigung
Die Berechnung der Hangneigung erfolgt auf der Grundlage eines Airbornlaserscans.
Dieser wurde bei der OÖ-Landesregierung erworben, und steht als xyz-Punktdatei zur Verfügung.
Eine Alternative zu dem hier beschriebenen Weg wäre eine Ermittlung mittels 1 Meter Höhenschichtlinien, welche den Kommunen von der OÖ-Landesregierung zur Verfügung gestellt werden. Hierbei kann die Hangneigung manuell oder automatisch ermittelt werden. Bei der manuellen Methode wird anhand des Höhenunterschiedes zwischen der höchsten und der die niedrigsten Höhenlinie, und der Länge der Einzugsfläche die Neigung berechnet. Um diesen Ablauf zu automatisieren, sind die selben Schritte notwendig wie ab Punkt 4.2.1 beschrieben.
4.2.1 Erstellung eines DGM
Bevor mit der Ermittlung der Hangneigung begonnen werden kann, muss zuerst ein DGM aus den Punktraster erstellt werden. Wie bereits erwähnt, hat der Laserscann eine Punktdichte von ca. 16Pkt./m², was zu einer sehr hohen Datengröße führt, und die Open Source Produkte an die Leistungsgrenze bringt. Um eine schnellere Berechnung zu ermöglichen, stünde auch eine Punktwolke mit einer geringeren Punktdichte zur Verfügung. Dieser könnte jedoch zu Verlusten in der Genauigkeit führen, und daher in diesem Fall nicht angewandt wird.
Für das Arbeiten mit einer hohen Datenmenge hat sich SAGA als leistungsstärkstes Programm behauptet.
Vorweg müssen in SAGA die ALS-Punkte geladen werden. Dies erfolgt über den Befehl "Import Point Cloud from Text File". Aus diesen Punkten wird nun das DGM Raster erstellt. Da keine Bruchkanten vom Untersuchungsgebiet vorliegen, und auch nicht unbedingt benötigt werden, wird direkt eine Interpolation mit der "Inverse Distance Weighting" (IDW) Methode durchgeführt. Für nähere Informationen zum Thema Interpolationsverfahren wird auf die Fachliteratur der Geostatistik verwiesen.
Das Ergebniss der Analyse ist ein Raster mit einer Auflösung von 1m. Für Geländeanalysen und Interpretationen stehen in SAGA, unter "Terrain Analyses", sehr viele Anwendungen zur Verfügung. Zur besseren Veranschaulichung des Untersuchungsgebietes wurde mit dem Tool "Analytical Hillshading" ein Geländerelief erzeugt. Anhand der Abbildungen 13 und 14 ist sehr gut zuerkennen, dass es sich um ein sehr ebenes Gelände handelt.
Abbildung 17: links: Gelände Relief des Untersuchungsgebietes; rechts: DGM des Untersuchungsgebietes erstellt in SAGA
4.2.2 Berechnung des Gefälles
Mit dem Erstellten DGM kann nun das Gefälle der Einzugsflächen ermittelt werden.
Die Berechnung des Gefälles der Einzelnen Rasterzellen, erfolgt mit dem Befehl
"Slope, Aspect, Curvature".
Das Ergebnis der Berechnung sind Raster für das Gefälle (Slope), und die Ausrichtung
Gefälleraster von Bedeutung. In SAGA wird als Einheiten für die Winkel das Bogenmaß verwendet.
Wie in QGIS gibt es auch in SAGA den Befehl zur statistischen Auswertung eines Rasters mit einem darüber liegenden Polygonshape. Die Anwendung heißt hier "Grid Statistics for Polygons". Man kann hier aus mehreren Statistikfunktionen auswählen.
Für die Berechnung des Oberflächenabflusses wird das durchschnittliche Gefälle eingegeben, dafür wird die Statistikoption "mean" ausgewählt. Im Unterschied zu QGIS wird hier keine csv-Datei erzeugt sonder ein Shape.
Das Gefälle, welches mit Bogenmaß im Shape eingetragen ist, muss noch mit Hilfe eines Feldrechners in Prozent umgerechnet werden. Die Umrechnung des Gefälles in Prozent erfolgt nach der Formel "tan()*100". Aufgrund dessen, dass auch QGIS in Bogenmaßeinheiten rechnet, müssen die Werte nicht nach Grad umgerechnet werden.
Der Ausdruck im QGIS Feldrechner für die neue Spalte "Prozent" lautet daher
"tan(mean)*100".
4.3 Zuordnen des Bodenaufbaus
Der Bodenaufbau wird, wie bereits erwähnt, aus der digitalen Bodenkarte abgeleitet.
Für den Oberflächenabfluss entscheidend ist die Tabelle "BODENART_G", welche die Bodenart unter der Humusschicht beschreibt.
Die Übernahme der Attribute aus der Bodenkarte zu den Einzugsflächen wird wieder in QGIS durchgeführt. Hier stehen mehrere GIS-Operation zur Auswahl. Eine Möglichkeit ist, die beiden Polygonshapes miteinander zu verschneiden (intersect). Das Problem bei dieser Methode ist jedoch, dass wenn sich eine Einzugsfläche über zwei verschiedene Bodenarten erstreckt, das Programm eine weitere Einzugsfläche erstellt. Da für jede Einzugsfläche nur ein Oberflächentyp vergeben werden kann, erscheint die nötige Nacharbeit des Befehles zu aufwändig. Eine andere Methode die zwei Shapes miteinander zu verknüpfen, stellt der Befehl "Attribute nach Position zusammenführen".
Hier werden die Attribute aufgrund ihrer räumlichen Überlagerung miteinander verknüpft. Diese Methode ist vor allem bei der Übernahme von Polygon-Attributen zu Punkt-Attributen sehr häufig in Gebrauch, funktioniert aber auch mit 2 Polygonshapes.
Der Vorteil gegenüber der Verschneidung ist, dass im Ziellayer die Anzahl der
Einzugsflächen nicht verändert wird. Der Befehl dafür ist in QGIS unter "Vektor" -
"Datenmanagement-Werkzeuge" zu finden. Mit den in Abbildung 16 getätigten Einstellungen wird das gewünschte Ergebnis geliefert.
Abbildung 18: Eingabeformular des "Attribute nach Position zusammenführen" in QGIS
Auch bei dieser Methode muss das Ergebnis noch etwas nachgebessert werden. Man kann nicht auswählen, welche Spalten übernommen werden sollen. Dies führt dazu, dass die Einzugsflächen alle Tabellen aus der Bodenkarte erhalten, und somit die Attributtabelle sehr unübersichtlich wird. Ein größeres Problem stellen die Einzugsflächen dar, welche verschiedene Bodenarten überlagern. Hier wird jede Einzugsfläche einzeln nachgebessert, und die Bodenart mit den größeren Anteil eingetragen.
Leider weist die digitale Bodenkarte im urbanen Bereich häufig Flächen auf, in denen keine Bodentypen eingetragen sind (Abbildung 17). Dies liegt daran, dass in versiegelten Stadtgebieten nur sehr wenige Bodenproben vorhanden sind. Um für diese Flächen einen Bodentyp zu erfassen, muss auf analoge Karten oder Erfahrungswerten zurückgegriffen werden. Diese Werte werden manuell in die Shapes eingetragen. Der Spaltenname im Einzugsflächenshape für die Bodenart lautet "Bodenart_G".
Abbildung 19: Ausschnitt aus der digitalen Bodenkarte des BFW im Untersuchungsgebiet, bei den weißen Flächen sind keine Bodentypen vorhanden
4.4 Automatische Zuweisung des Oberflächentyps
Da die Einzugsflächen bereits Befestigungsgrad, Gefälle und Bodenart enthalten, kann jetzt der Oberflächentyp zugeordnet werden. Die Zuweisung der Oberflächentypen erfolgt nach dem unter Punkt 3.1.2 erläuterten Kriterien. In der Abbildung 18 werden diese noch einmal vereinfacht dargestellt.
Abbildung 20: Einteilung der Oberflächentypen
Damit eine automatische Zuweisung der Oberflächentypen zu den Einzugsflächen erfolgen kann, wird eine, von QGIS im Feldrechner unterstützte, "if when then"- Funktion verwendet. Die Formel für die oben dargestellte Einteilung lautet wie folgt:
case when "Abflussbei" <= 0.5 AND Prozent < 4 AND BODENART_G ILIKE '%lehmiger sand%' then 'WGB LS NG I/II' when "Abflussbei" <= 0.5 AND Prozent >= 4 AND "Prozent"
<0.9 AND BODENART_G ILIKE '%lehmiger sand%' then 'WGB LS NG III/IV' when
"Abflussbei" <= 0.5 AND Prozent < 4 AND BODENART_G ILIKE '%Lehm%' then 'WGB L NG I/II' when "Abflussbei" <= 0.5 AND Prozent >= 4 AND "Prozent" <0.9 AND
BODENART_G ILIKE '%Lehm%' then 'WGB L NG III/IV' when "Abflussbei" <= 0.5 AND Prozent < 4 AND BODENART_G ILIKE '%sandiger Kies%' then 'WGB KS NG I/II' when
"Abflussbei" <= 0.5 AND Prozent >= 4 AND "Prozent" <0.9 AND BODENART_G ILIKE '%sandiger Kies%' then 'WGB KS NG III/IV' when "Abflussbei" > 0.5 AND Prozent < 4 AND BODENART_G ILIKE '%lehmiger sand%' then 'STZ LS NG I/II' when "Abflussbei" > 0.5 AND Prozent >= 4 AND "Prozent" <0.9 AND BODENART_G ILIKE '%lehmiger sand%' then 'STZ LS NG III/IV' when "Abflussbei" > 0.5 AND Prozent < 4 AND BODENART_G ILIKE '%Lehm%' then 'STZ L NG I/II' when "Abflussbei" > 0.5 AND Prozent >= 4 AND "Prozent"
<0.9 AND BODENART_G ILIKE '%Lehm%' then 'STZ L NG III/IV' when "Abflussbei" > 0.5 AND Prozent < 4 AND BODENART_G ILIKE '%sandiger Kies%' then 'STZ KS NG I/II' when
"Abflussbei" > 0.5 AND Prozent >= 4 AND "Prozent" <0.9 AND BODENART_G ILIKE '%sandiger Kies%' then 'STZ KS NG III/IV' when "Abflussbei" = 0.9 and "Prozent" <=2.5 then 'STR NG I'when "Abflussbei" = 0.9 and "Prozent" >2.5 and "Prozent" <=4 then 'STR NG II' when "Abflussbei" = 0.9 and "Prozent" >4 and "Prozent" <=7.5 then 'STR NG III' when
"Abflussbei" = 0.9 and "Prozent" >7.5 then 'STR NG IV' else '' end
Im Feldrechner können keine Ausdrücke gespeichert bzw. importiert werden. Daher muss man sich mit "copy" und "paste" abhelfen und den Befehl in einem externen Dokument für zukünftige Anwendungen abspeichern.
4.5 Automatischer Import der Parameter 4.5.1 Erstellen einer csv-Datei
Nachdem alle Einzugsflächen einen Oberflächentyp erhalten haben, können sie mittels ASCII-Import in die Kanaldatenbank importiert werden. Damit das Polygonshape mit den Einzugsflächen etwas übersichtlicher wird, werden in QGIS die überflüssigen Spalten gelöscht. Somit bleiben nur noch die Spalten Bezeichnung, Abflussbeiwert und Oberflächentyp übrig.
Für den ASCII-Import in BASYS wird eine Textdatei benötigt. Besonders gut geeignet sind hier csv-Datein (comma separated Value). Kommaseparierte Daten werden durch Trennzeichen (zb. Komma oder Strichpunkt) getrennt. Diese Dateien werden von den meisten Tabellenkalkulationsprogrammen und Datenbanken unterstützt und sind dadurch bestens für Import/Export Funktionen geeignet (Goyvaerts & Levithan, 2010).
QGIS bietet die Möglichkeit Vektordaten in zahlreiche Formate zu speichern wie zB.:
GML, Interlis, KML, DXF oder csv. Für die Einzugsflächen wird das Format "Komma- separierte Werte [csv]" ausgewählt. Der Dateiname wird unter "Speichern unter"
eingegeben.
4.5.2 Erweitern einer BASYS-Vorlage
In BASYS gibt es bereits eine Fülle an ASCII-Import Vorlagen, jedoch enthält keine dieser Vorlagen alle ermittelten Parameter der Einzugsflächen. Aus diesem Grund wird eine Vorlage so erweitert, dass alle Parameter auf einmal importiert werden können.
Die Vorlagen des ASCII-Imports sind in einem XML-Dokument definiert und können auch in diesem verändert werden. Für die Bearbeitung eines XML-Dokuments stehen mehrere Programme zur Auswahl, wie zB. auch der herkömmliche Texteditor von Windows. Da im Windows Editor nur eine einfache Textansicht der XML-Datei zur Verfügung steht, wird die Verwendung eines professionellen XML-Editor bevorzugt.
Ein sehr umfangreicher XML-Editor ist zB. die gebührenpflichtige Software XMLSpy der Firma altova. Dem gegenüber steht mit Notepad++, ein Open Source Programm, welches die gewünschten Aufgaben ebenfalls vollständig erfüllen kann. Da es ein Ziel der Thesis ist, alle Aufgaben mit OpenSource Produkten zu lösen, wird für die Bearbeitung der XML-Datei Notepad++ verwendet.
In Abbildung 19 ist die Vorlage eines ASCII-Imports zu sehen. In Zeile 1 ist, wie in jeder XML-Datei, die XML-Version eingetragen. Unter "Configuration", in Zeile 5, ist bei "designation" "Einzugsflaeche" eingetragen, dies bedeutet, dass es sich um einen Import für das Einzugsflächenformular handelt. Die Parameter welche importiert werden, sind unter "Destination" zu finden (Zeile 7-25). Die in den vorhergegangenen Schritten ermittelten Parameter befinden sich in den Zeilen 18 (Abflussbeiwert), und in Zeile 25 (Oberflächentyp). Die Zeile 25, für den Import des Oberflächentyps, wurde in notepad++ neu erstellt.
Zwei wichtige tags im XML-Dokument sind "designation" und "opath". Ersteres ist der Name des zu Importierenden Wertes, welcher im ASCII-Import Modul angezeigt wird, und "opath" ist die Bezeichnung des Datenbankfeldes in BASYS.
Abbildung 21: Auszug aus der geänderten XML-Vorlage
Das fertige XML-Dokument wird in BASYS, in der Schnittstelle "ASCII-Import", als Vorlage geladen. Bevor der Import gestartet werden kann, müssen die Feldpositionen der csv-Datei den richtigen Feldern des Imports zugewiesen werden. Nachdem alle Parameter erfolgreich eingegeben wurden, kann nun mit der hydrodynamischen Kanalnetzberechnung begonnen werden.
5. Analyse der Ergebnisse
Im folgenden Abschnitt werden die unter Punkt 1.6 gestellten Fragen beantwortet.
5.1 In welchen Bereichen können GIS gestützte Analysen sinnvoll für die hydrodynamische Kanalnetzberechnung eingesetzt werden?
Die oben erläuterten Analysen beschreiben den Ablauf für die Berechnung mit HydroCad. Es ist aber bei allen Kanalnetzberechnungsprogrammen notwendig, die befestigten Flächen, die Hangneigung und der Bodenaufbau zu ermitteln. Bei der Definition dieser Parameter macht es durchaus Sinn, wenn man auf die Unterstützung von GIS-Programmen setzt. Besonders bei der Ermittlung der befestigten Flächen ist der Einsatz von Geoinformationsprogrammen nicht mehr wegzudenken.
5.2 Welche Möglichkeiten gibt es, um befestigte Flächen aus Orthophotos zu ermitteln?
Bei Punkt 4.1 wurden einige Methoden für die Ermittlung der befestigten Flächen aufgezeigt. Stand der Technik, für die Bestimmung der Oberflächennutzung, ist aber auf jeden Fall eine Analyse mittels OBIA, für die jedoch viel Erfahrung nötigt ist (ALBERTZ, 2009). Man hat nun die Möglichkeit diese Analysen an einen GIS- Spezialisten zu vergeben, oder man hat das Glück, dass für das betroffenen Gebiet LISA-Daten zur Verfügung stehen.
Da jedoch Erfahrungen mit OBIA Analysen in den meisten Zivilingenieurbüros nicht vorhanden sind, und auch noch keine flächendeckenden LISA-Daten zur Verfügung stehen, wurde in der Thesis versucht, mit einer einfacheren Methode die benötigten Ergebnisse zu erzielen.
5.3 Liefert die Berechnung des Befestigungsgrades anhand des NDVI ausreichend genaue Ergebnisse für die Kanalnetzberechnung?
Mit den importierten Daten wird ein erster Berechnungsvorgang in BASYS-Hydrocad durchgeführt. Da es mit sehr hoher Wahrscheinlichkeit zu Abweichungen zwischen den Aufzeichnungen der Durchflussmessungen und den berechneten Werten kommt, dient der erste Berechnungslauf nur zur Kalibrierung des Kanalnetzes. Bei dieser Kalibrierung wird zuerst kontrolliert, ob die Funktion aller Schieber und Wehre richtig eingegeben wurden. Als nächster Schritt werden die Befestigungsanteile der
