ArcGIS: Daten aus CF-NetCDF-Datei importieren
Aus BAWiki
Vorbemerkung
Diese Seite erläutert den Import von Daten aus einer Datei des Typs CF-NetCDF nach ArcGIS. Neben den geografischen Koordinaten kann dieser Dateityp auch die X- und Y-Koordinaten des Projektionssystems enthalten. Das zu Grunde liegende geografische Koordinatensystem und die verwendete Projektion sind im Attribut comment der den Daten zugeordneten Grid Mapping Variable benannt. Die geografischen Koordinaten beziehen sich auf European Terrestrial Reference System 1989 (ETRS_1989). Für die Projektion wird in diesem Beispiel Gauß-Krüger-Streifen_2 angenommen.
Daten importieren, bzw. Featurelayer erzeugen
- ArcMap öffnen
- Neuen Data Frame anlegen und aktivieren
Anmerkung: Koordinatensystem des verwendeten Kartenrahmens soll undefiniert sein. Prüfe in seinen Eigenschaften, ob current coordinate system auf no projection steht. - ArcToolbox öffnen
- Ordner Multidimension Tools der ArcToolbox öffnen
- Tool Make NetCDF Feature Layer per Doppelklick starten
- Dialog Make NetCDF Feature Layer ausfüllen:
- unter Input netCDF File die CF-NetCDF-Datei auswählen
- unter Variables die verschiedenen Daten für den Import angeben
Bspw. Mesh0_node_depth mit den Tiefen an den Knoten
Achtung: hier nicht die Koordinatenvariablen nennen! - unter X Variable die X-Koordinate der Projektion oder die Variable mit der geografische Länge (z.B. Mesh0_node_lon) einstellen.
- unter Y Variable die Y-Koordinate der Projektion oder die Variable mit der geografische Breite (z.B. Mesh0_node_lat) einstellen.
- unter Output Feature Layer dem neuen Layer einen Namen geben
- unter Row Dimensions die Dimension mit der Anzahl der Geo-Objekte (z.B. nMesh0_node) einstellen
- Dialog Make NetCDF Feature Layer mit OK-Button beenden
- Warten bis Importvorgang abgeschlossen (Completed)
- Tool Make NetCDF Feature Layer mit CLOSE-Button schließen
- Neu angelegten Feature Layer prüfen (z.B Open Attribute Table im Layer-Kontextmenü)
Featurelayer als Shapefile exportieren
- Ordner Conversion Tools der ArcToolbox öffnen
- Ordner To Shapefile in den Conversion Tools öffnen
- Tool Feature Class To Shapefile (multiple) per Doppelklick starten
- Dialog Feature Class To Shapefile (multiple) ausfüllen:
- unter Input Features den eben angelegten Feature Layer auswählen
- unter Output Folder ein Zielverzeichnis für den Shapefile einstellen
- Dialog Feature Class To Shapefile (multiple) mit OK-Button beenden
- Warten bis Exportvorgang abgeschlossen (Completed)
- Tool Feature Class To Shapefile (multiple) mit CLOSE-Button schließen
Koordinatensystem des Shapefile setzen
Leider ist ArcGIS nicht in der Lage aus den Angaben der in der Datei enthaltenen Grid Mapping Variablen das zu Grunde liegende Koordinatensystem exakt zu zu ordnen. Um den Shapefile für die ArcGIS-Welt eindeutig zu beschreiben,
muss das Koordinatensystem noch mal explizit gesetzt werden.
Der hier beschriebene Weg verwendet dazu das Programm ArcCatalog.
- ArcMap schließen
- ArcCatalog öffnen
- Den erstellten Shapefile in ArcCatalog ansteuern
- Über Kontextmenü der Shapedatei den Shapefile Properties-Dialog öffnen
- Registerblatt XY Coordinate System wählen
- Button Select... klicken um das entsprechende Koordinatensystem einzustellen
- Beispiel: Gauß-Krüger-Streifen_2
Browse to- Projected Coordinate Systems
- National Grids
- DHDN 3 Degree Gauss Zone 2.prj
- Beispiel: ETRS_1989
- Geographic Coordinate Systems
- Europe
- ETRS 1989.prj
- Beispiel: Gauß-Krüger-Streifen_2
- Den Dialog Browse for Coordinate System per Add-Button schließen.
- Button Übernehmen klicken um die Eigenschaft zu übernehmen
- Button OK klicken Eigenschaften-Dialog zu schließen
- ArcCatalog schließen
Tabelle: Zuordnung der ArcGIS-Koordinatensysteme
Eintrag in CF-NetCDF-Grid-Mapping-Variable | ArcGIS PRJ-File |
---|---|
Ellipsoid - European Terrestrial Reference System 1989 | ETRS 1989.prj |
3 degree Gauss-Kruger zone 2 - Germany west of 7.5 E - Bessel 1841 | DHDN 3 Degree Gauss Zone 2.prj |
3 degree Gauss-Kruger zone 3 - Germany 7.5 E to 10.5 E - Bessel 1841 | DHDN 3 Degree Gauss Zone 3.prj |
3 degree Gauss-Kruger zone 4 - Germany 10.5 E to 13.5 E - Bessel 1841 | DHDN 3 Degree Gauss Zone 4.prj |
3 degree Gauss-Kruger zone 5 - Germany 13.5 E to 16.5 E - Bessel 1841 | DHDN 3 Degree Gauss Zone 5.prj |
UTM zone 31N - Europe 0 E to 6 E - ETRS89 | ETRS 1989 UTM Zone 31N.prj - ( 6-, 7- oder 8-stellig ) |
UTM zone 32N - Europe 6 E to 12 E - ETRS89 | ETRS 1989 UTM Zone 32N.prj - ( 6-, 7- oder 8-stellig ) |
UTM zone 33N - Europe 12 E to 18 E - ETRS89 | ETRS 1989 UTM Zone 33N.prj - ( 6-, 7- oder 8-stellig ) |
Temporäre Koordinaten-Transformation zu Darstellungszwecken
ArcGIS kann Daten, die sich auf unterschiedliche Koordinatensysteme beziehen, gemeinsam in einem Kartenrahmen darstellen. Die Koordinaten müssen dafür transformiert werden, damit räumlich identische Objekte übereinander dargestellt werden. ArcGIS bietet eine Auswahl an Transformationsmethoden an, es erlaubt aber auch das Einrichten eigener Methoden.
Das Koordinatensystem eines Kartenrahmens lässt sich explizit über seine Eigenschaften setzen. Ist es nicht gesetzt wird das Koordinatensystem des ersten hinzugefügten Featurelayers übernommen.
Wird einem Kartenrahmen mit definiertem Koordinatensystem ein Shapefile mit Koordinaten eines
abweichenden Systems zugefügt wird der Geographic Coordinate Systems Warning Dialog angezeigt, der zur Angabe einer Transformation auffordert.
Hinweis: Dieser Mechanismus greift nur, wenn das Koordinatensystem des Shapefiles, wie im vorangehenden Abschnitt beschrieben, ArcGIS-verständlich gesetzt wurde!
Im Geographic Coordinate Systems Warning Dialog:
- Button Transformations... klicken, um den Geographic Coordinate System Transformations Dialog zu öffnen
- unter Convert from: das Koordinatensystem der hinzuzufügenden Daten einstellen
- unter Into: das Ziel-Koordinatensystem angeben
- unter Using: eine geeignete Transformations-Vorschrift bestimmen
- Dialog mit OK-Button schließen
- Dialog über Close-Button schließen
Der Shapefile bleibt von dieser Maßnahme unberührt!
Den Geographic Coordinate System Transformations-Dialog erreicht man auch direkt über das Registerblatt Coordinate System im Eigenschaften-Dialog des Kartenrahmens.
Tabelle: Transformationsmethoden
Von | Nach | Transformation |
---|---|---|
DHDN 3 Degree Gauss Zone 2.prj | ETRS 1989.prj | DHDN_To_ETRS_1989_8_NTv2 |
DHDN 3 Degree Gauss Zone 3.prj | ETRS 1989.prj | DHDN_To_ETRS_1989_8_NTv2 |
DHDN 3 Degree Gauss Zone 4.prj | ETRS 1989.prj | DHDN_To_ETRS_1989_8_NTv2 |
DHDN 3 Degree Gauss Zone 5.prj | ETRS 1989.prj | DHDN_To_ETRS_1989_8_NTv2 |
zurück zu ArcGIS-Anwendungen