Diskussion:Reliefkarten erstellen: Unterschied zwischen den Versionen

Aktuelle Version vom 15. August 2014, 23:42 Uhr

Notizen:

Erstellung einer Reliefkarte ausschließlich mit Quantum GIS:

Testumgebung Quantum Gis Lisboa Windows 7 64bit

• SRTM Geotiff Kacheln laden (Layer → Rasterlayer hinzufügen: (GDAL) GeoTIF (*.tif) )
• Kacheln zusammenfügen (Raster → Sonstiges → Verschmelzen)
  • Eingabedateien auswählen
  • Ausgabedateinamen festlegen
  • OK klicken

• Erstellte Datei laden (Layer → Rasterlayer hinzufügen: (GDAL) GeoTIF (*.tif) )
• Quellkacheln aus Projekt entfernen (Vermeidung von Datenmüll)

Gebietsrahmen erstellen:

• MTB Schnitt wird als Gebietsrahmen gewünscht
• MTB-Grenzpolygone als Shapedatei laden (Wichtig WGS84 EPSG 4326 als Georeferenz, da sonst nicht für das Ausschneiden des Reliefs geeignet)
• Gewünschte MTB durch Rechteck auswählen (Ansicht → Auswahl → Objekte durch Rechteck wählen)
• Mit dem Objekte abfragen Tool (I) kontrollieren, ob Objekte korrekt abgefragt
• Ausgewähltes MTB-Raster als neues Shape abspeichern: Rechte Maustaste auf Layereintrag in Leiste links, dort "Auswahl speichern als" wählen.
• Datei Karte hinzufügen, wenn nicht bereits beim Speichern ausgewählt.
• Deutschland-MTB Raster kann jetzt entfernt werden.

• Verschmelzen der Raster zu einem Gebietsumriß (Projektrahmen): Vektor → Geoverarbeitungswerkzeuge → Auflösen:
  • Im Dialogfenster den richtigen Layer auswählen, Auflösungsfeld: "--- Alles auflösen ---" wählen. Ausgabeshapedatei (Namen und Pfad) festlegen. Mit [OK] starten.
• Datei laden, wenn nicht beim Speichern ausgewählt
• Höhen-Rasterdaten anhand des eben erstellen Projektrahmens ausschneiden: Raster → Extraktion → Clipper
  • Im Dialog "Clipper" Eingabedatei auswählen, Auschnittmodus auf Maskenlayer setzen, Projektrahmen als Maskenlayer auswählen
  • Nach Abschluß "zu Karte hinzufügen" nicht auswählen, führt gelegentlich zu Fehlern.

• Dialog schließen (Close), Raster hinzufügen, großes Höhenrsater kann entfernt werden.
• Schummerung erstellen: Raster → Geländeanalyse → Schummerung wählen.
  • Im Dialog Schummerung folgende Einträge: Geländehöhenlayer: auswählen, Ausgabelayer: (Pfad und Name) festlegen, Ausgabeformat: Geotiff, Z-Faktor: mindestens 11000, für Überhöung höher, z.B. 55560, Azimuth: 300, Vertikaler Winkel: 40. Für bessere Ergebnisse mit den Werten experimentieren.
• Erstellte Datei laden. Rechte Maustaste auf Layereintrag in Leiste links, dort "Eigenschaften" wählen.
• Im Dialogfenster Reiter "Transparenz" wählen. Kein Datum-Wert: auf -32768 setzten, ebneso transparente Pixel Liste: -32768, 100%. Reiter Stil: Kontrastverbesserung: Strecke auf MinMax wählen.

Erstellen eines Kartenrahmens

Erstellen eines Kartenrahmens für MTB-Rasterkarten

• Einbinden des Quadranten-Layers:
  • Layer → Vektorlayer hinzufügen, Layer d_mtbq_epsg3395 auswählen (aus Deutschland Resource)
  • umrahmende MTB-Reihe auswählen mittels Rechteck-Auswahl-Werkzeug
  • Auswahl als Layer abspeichern: Layer → Auswahl als Vektordatei speichern. Namen vergeben. ZB: Kartenrahmen_qu_epsg3395. Gespeicherte Datei als Layer hinzufügen.
  • Deutschland Layer kann jetzt entfernt werden (d_mtbq_epsg3395)
  • Umwandlung des Polygonlayers in einen Linienlayer:
    • Vektor → Geometriewerkzeuge → Polygon zu Linie. Im Dialogfenster Layer Kartenrahmen_qu_eps3395 wählen. Dateinamen für Ausgabedatei festlegen, z.B. Kartenrahmen_qu_li_epsg3395, bei Nachfrage dem Projekt als neuen Layer hinzufügen

Erstellen passender Ausschnitte aus den Vektordaten

• Einbinden der Vektor-Layer, wie z.B. Flüsse:
  • Layer → Vektorlayer hinzufügen, z.B. Layer d_gw_fluesse1_linie_epsg3395 auswählen (aus Deutschland Ressource)
  • Vektor → Geoverarbeitungsverkzeuge → Clip
  • Im Dialogfenster wählen:
    • Eingabevektorlayer: d_gw_fluesse1_linie_epsg3395;
    • Layer zuschneiden: rahmen_epsg3395 (der Rahmen, der zur Abgrenzung des Gebietes erstellt wurde, d.h. die Maskendatei);
    • Ausgabeshapedatei: z.B. rheinland_gw_fluesse1_linie_epsg3395
  • ►Wichig: die Maskendatei muss die gleiche Georeferenzierung haben wie die Datei, aus der ausgeschnitten werden soll.
  • Nacheinander alle benötigten Layer einbinden und nach dem obigen Schema ausschneiden

Schön gestalteten Kartenrahmen erstellen

Um einen schön gestalteten Kartenrahmen zu erhalten, sind zahlreiche Einzelschritte notwendig. Ich zeige hier einen möglichen Weg auf, dieser kann jedoch je nach Designwünschen erheblich variieren.

Folgendes Design ist vorgegeben: Um eine Karte, die durch die Eckkoordinaten eines Messtischblattrahmens definiert wurde, soll ein Kartenrahmen erstellt werden, der oben die Spalten der Messtischblattnummern, links die Reihen der Messtischblattnummern, unten eine Auswahl an Längengraden und rechts entsprechende Breitengrade enthält. Der Rahmen soll gleichmäßig in einem in Metern definierten Abstand um das Kartenbild erstellt werden. MTB-Reihen/Spaltennummern mittig.

Pufferlayer erstellen

Zunächst wird ein Puffer um das Kartenbild erstellt. Die Maßzahl für den Pufferabstand ist vom Koordinatenbezugssystem (KBS/CRS) des Quelllayers abhängig. Wir wählen den bereits erstellten Rahmenlayer in Raumreferenzsystem "WGS 84/World Mercator", epsg 3395. Ein passender Pufferabstand für Karten mit ca. 16 MTB-Spalten liegt bei diesem RBZ bei 120000.

• Vektor → Geoverarbeitungswerkzeuge → Puffer
• Im Dialogfenster folgende Einstellung:
  • Eingabevektorlayer: [Name des Rahmenlayers]
  • Pufferabstand: 12000
  • beim Rest die Standardwerte belasssen, Ausgabelayer festlegen, z.B. mit puffer_12000 bezeichnen

QGIS erstellt nun einen Puffer mit abgrundeten Ecken. Dieses Verhalten lässt sich in der Version 2.0 Dufour nicht ändern. Deswegen ist ein weiterer Schritt notwendig:

• Vektor → Forschungswerkzeuge → Polygon aus Layergrenzen
• Im Dialogfeld folgende Einstellung:
  • Eingabelayer: [Name des soeben erstellten Layers, hier puffer_12000]
  • Ausgabepolygonlayer: z.B. kartenrahmen_epsg3395

Bei Nachfrage den soeben erstellte Polygonlayer dem Projekt hinzufügen

Der Pufferlayer puffer_12000 kann jetzt aus dem Projekt entfernt und dann auch von der Festplatte gelöscht werden.

Es wird noch ein zweiter Pufferlayer benötigt, der mit dem halben Abstand analog zu oben erstellt wird. Er kann z.B. als kartenhilfsrahmen_epsg3395 bezeichnet werden.

Rastergitter erzeugen

Im nächsten Schritt wird ein Rastergitter erzeugt, das jeweils die Mitte der MTB-Raster darstellt. Dazu folgende Überlegungen: Das deutsche MTB-Raster ist ein geographisches Raster, das auf Basis des Bessel-Ellipsoids und des RBS DHDN definiert ist. Der EPSG-Code dafür ist 4314. Jedes MTB ist 10 geographische Minuten breit und 6 hoch.

Für die nächsten Schritte muss das KBS des Projektes umgestellt werden:

• Projekt → Projekteinstellungen
  • Im Dialogfenster, linke Leiste, KBS (Koordinatenbezugysystem) wählen.
  • Dort DHDN, epsg 4314 einstellen und bestätigen.

Die Karte wird jetzt umprojeziert und sieht zunächst etwas merkwürdig verzerrt aus. Jetzt wird das Vektorgitter erstellt.

Zunächst die groben Eckkoordinaten der Karte in Grad ermitteln. Dazu die Darstellung der Koordinaten umstellen:

• Projekt → Projekteinstellungen
  • Im Dialogfenster, linke Leiste, "Allgemein" wählen.
  • Dort Karteneinheiten: Grad und Gradanzeige: Grad, Minuten, Sekunden wählen und bestätigen

Danach mittels der Maus die West-, Ost-, Nord- und Südgrenze des Rahmens ermitteln. Im Beispiel sind es die Koordinaten:

 W 5°50' E 8°30' N 52°00', S 49°30'

Um ein entsprechendes Rastergitter zu erhalten, wird das Raster um 15 Minuten nach Westen bzw. Osten und 9 Minuten nach Süden bzw. Norden ausgedehnt:

 W 5°35', E 8°45', N 52°09', S 49°21'

Diese Daten müssen in Dezimalnotation umgerechnet werden:

x min 5.5833333333333333
x max 8.75
y min 49.35
y max 52.15

• Vektor → Forschungswerkzeuge → Vektorgitter
  • Im Dialogfenster für das oben genannte Beispiel folgende Angaben:
  • Gittergrenze des Layer:
    • kartenrahmen_epsg3395 (???)
    • x-Min 5.5833333333333333
    • x-Max 8.75
    • y-Min 49.35
    • y-Max 52.15
  • Parameter
    • X: 0,1666666666666666
    • Y: 0,1

Danach muss das korrekte Koordinatenbezugssystem dem erstellen Layer zugeordnet werden:

• Vektorgitter-Layer aus der Layerleiste auswählen.
• mit rechter Maustaste Kontext-Menü öffnen, auswählen: KBS für Layer setzen.
• im Dialogfenster KBS DHDN, EPSG:4314 einstellen und bestätigen.

Spalten und Reihenlinien erzeugen

• Notation der Spalten errechen und zuordnen.
  • Vektorgitter-Layer aus der Layerleiste auswählen.
  • mit rechter Maustaste Kontext-Menü öffnen, auswählen: Bearbeitungsstatus umschalten
  • mit rechter Maustaste Kontext-Menü öffnen, auswählen: Attributtabelle öffnen
  • Im Fenster Attributtabelle das Werkzeug Feldrechner (ganz rechts) anklicken
  • im Dialogfenster Feldrechner Folgendes eingeben/auswählen:
    • nur ausgewählte Objekte aktualisieren: deselektiert
    • neues Feld anlegen: auswählen
    • Ausgabefeldname: Spalte
    • Ausgabefeldtyp: Text
    • Ausgabefeldbreite: 2
    • ins Textfeld Ausdruck folgenden Code kopieren, danach Berechnung mit [OK] starten:

CASE WHEN  "COORD" > 49 THEN '' ELSE  right('0' || tostring (floor( (("COORD" - (5.833333333333)) * 6) + 1 )) , 2) END

• Notation der Reihen errechen und zuordnen:
  • Im noch offenen Fenster Attributtabelle erneut das Werkzeug Feldrechner (ganz rechts) anklicken
  • im Dialogfenster Feldrechner Folgendes eingeben/auswählen
    • nur ausgewählte Objekte aktualisieren: deselektiert
    • neues Feld anlegen: auswählen
    • Ausgabefeldname: Reihe
    • Ausgabefeldtyp: Text
    • Ausgabefeldbreite: 2
    • ins Textfeld Ausdruck folgenden Code kopieren, danach Berechnung mit [OK] starten

CASE WHEN  "COORD" > 49 THEN right('0' || tostring( floor( ( (55.1 - "COORD")/0.1) +9)) ,2) ELSE  ' ' END

• nicht vergessen: Änderungen speichern!

ab hier nur mehr Stichpunkte

grad

floor("COORD")

min

floor(("COORD"-(floor("COORD")))*59.999999999999)+1

• Kartenhilfsrahmen in KBS DHDN EPSG:4312 abspeichern und laden
• Kartenhilfsrahmen EPSG:4312 in Linie umwandeln
• Schnittpunkte: