Mehr

Ausdehnungen so anpassen, dass sich CRS-Shapefiles in QGIS überlappen?

Ausdehnungen so anpassen, dass sich CRS-Shapefiles in QGIS überlappen?


Ich habe zwei Dateien mit identischen Projektionen/CRS. Wenn ich OTF ein- und ausschalte, kann ich sie nicht überlappen lassen: Ein Shapefile enthält Küstenpolygone von Massachusetts und das andere Shapefile enthält Küstenlinien für MA.

+proj=laea +lat_0=45 +lon_0=-100 +x_0=0 +y_0=0 +a=6370997 +b=6370997 +Einheiten=m +no_defs +proj=laea +lat_0=45 +lon_0=-100 +x_0 =0 +y_0=0 +a=6370997 +b=6370997 +Einheiten=m +no_defs

aber die Ausmaße sind

xMin,yMin 2321544.64,258328.06 : xMax,yMax 2333464.74,277553.12

vs.

xMin,yMin -176.685,17.9108 : xMax,yMax -65.2246,71.3413

Ich habe die Schritte in http://qgis.spatialthoughts.com/2012/04/tutorial-working-with-projections-in.html befolgt, mit ein- und ausgeschaltetem OTF versucht: Die Ebenen überlappen sich immer noch nicht.

BEARBEITEN:

Ich habe das SHP der MA-Küstenlinie (Ebene> speichern unter) mit jedem derLayer-CRS,Projekt CRSundAusgewählte CRSOptionen: Der Umfang bleibt gleich. Soweit ich das beurteilen kann, ist die gewählte Projektion (EPSG:2163) ein PCS, kein GCS.


Die Koordinaten des zweiten sinds NICHT in Laea-Projektion, aber in Grad.

Sie müssen also mit der rechten Maustaste klicken ->CRS für Layer festlegendie falschen CRS-Informationen zu WGS84 oder NAD 83 oder NAD27 zu korrigieren.

Ich weiß nicht, was Sie falsch gemacht haben, aber Sie haben möglicherweise den gleichen Fehler auch mit der ersten Schicht gemacht, und sie werden möglicherweise immer noch nicht ausgerichtet, nachdem nur die zweite Schicht desinfiziert wurde.

Am besten beginnen Sie bei Null mit einem frischen Download der Datenquellen in ein neues qgis-Projekt.


Ausdehnungen so anpassen, dass sich CRS-Shapefiles in QGIS überlappen? - Geografisches Informationssystem

Ausdehnungen des Grundwasserleiters und der Begrenzungseinheiten des Aquifersystems der atlantischen Küstenebene von North und South Carolina (pp1773_extents) Vector Digital Data Set (Polygon) Professional Paper 1773

https://water.usgs.gov/lookup/getspatial?pp1773_extents Campbell, B.G. Coes, A. L.

Grundwasserverfügbarkeit in der atlantischen Küstenebene von North und South Carolina Online-Publikation Professional Paper 1773

Die Ausdehnung jeder der 16 Modellschichten wurde verwendet, um das hydrogeologische Gerüst der Grundwasserleiter und eingrenzenden Einheiten zu beschreiben, die die atlantische Küstenebene von North und South Carolina bilden. Weitere Informationen zur Verwendung dieser Daten finden Sie in den Abschnitten zu den Untersuchungsmethoden von Professional Paper 1773: Gellici, JA, Lautier, JC, Kapitel B. Hydrogeological Framework of the Atlantic Coastal Plain, North and South Carolina of large work Campbell, BG, and Coes, AL, Hrsg., 2010, Grundwasserverfügbarkeit in der atlantischen Küstenebene von North Plain of North and South Carolina: US Geological Survey Professional Paper 1773, 241 S., 7pls.

Die Schichten, die die Ausdehnungen der Grundwasserleiter und begrenzenden Einheiten des Grundwasserleitersystems darstellen, wurden ursprünglich als CAD-Dateien formatiert und von Bruce Campbell in proprietäre Esri-Coverages umgewandelt. Die Coverages wurden dann für Archivierungszwecke in Shapefiles umgewandelt. Erscheinungsdatum 2010

Keine geplant -83.310243915 -75.429738739 37.588308642 31.218516725 USGS-Thesaurus Binnengewässer Grundwasser Grundwasserverfügbarkeit Grundwasserleiter Begrenzungseinheit Hydrogeologie hydrogeologische Modellierung Grundwassermodellierung hydrogeologisches Gerüst Oberflächenwasserleiter ISO 19115 Themenkategorie geowissenschaftlicheInformationen Binnengewässer Umgebung

Informationssystem für geografische Namen

US Geological Survey, Southeast Region Bruce G Campbell Hydrologist Postanschrift 720 Gracern Road Columbia SC

US 803-750-6161 803-750-6181 [email protected]

https://water.usgs.gov/GIS/browse/pp1773_extents.jpg
Darstellung des Datensatzes
jpg Diese Studie wurde vom U.S. Geological Survey (USGS) Groundwater Resources Program finanziert. Umgebung ab Metadatenerstellung: Microsoft Windows 7 Version 6.1 (Build 7601) Service Pack 1 ESRI ArcCatalog 10.2 (Build 3497) Service Pack [N/A] (Build [N/A]) Joseph A. Gellici Jeff C. Lautier

Bohrlochstandorte und Oberflächenhöhen jeder der 16 hydrostratigraphischen Schichten der atlantischen Küstenebene von North und South Carolina (pp1773_unit_alt_boreholes) Vector Digital Data Set (Point) Professional Paper 1773

https://water.usgs.gov/lookup/getspatial?pp1773_unit_alt_boreholes Bruce G. Campbell Alissa L. Coes Joseph A. Gellici Jeff C. Lautier

Oberflächenhöhe hydrogeologischer Schichten der atlantischen Küstenebene von North und South Carolina Vector Digital Dataset (Polygon)

Eine formale Genauigkeitsbewertung der horizontalen Positionsinformationen im Datensatz wurde nicht durchgeführt. Ein informeller visueller Vergleich der Linienführung zeigte eine korrekte horizontale Ausrichtung. Eine formale Genauigkeitsbewertung der vertikalen Positionsinformationen im Datensatz wurde entweder nicht durchgeführt oder ist nicht anwendbar.

Das Virginia Coastal Plain Hydrogeological Framework Digitale und/oder gedruckte Ressourcen

https://pubs.usgs.gov/pp/2006/1731/ Digital and/or Hardcopy Resources 2006 Veröffentlichungsdatum PP1731 Hydrogeologische Einheiten und Modellschichten für das Untersuchungsgebiet im südlichen Virginia Winner, M.D., und Coble, R.W.

Hydrogeologisches Framework der North Carolina Coastal Plain Digital- und/oder Hardcopy-Ressourcen

http://pubs.er.usgs.gov/publication/pp1404I Digital and/or Hardcopy Resources 1996 Veröffentlichungsdatum PP 1404-I Hydrogeologische Einheiten und Modellschichten für das Untersuchungsgebiet in North Carolina Coes, AL, Campbell, BG, Petkewich, MD und Fine, JM

Simulation der Grundwasserströmung in der atlantischen Küstenebene, North und South Carolina und Teilen von Eastern Georgia und Southern Virginia Vorentwicklung bis 2004 Digitale und/oder gedruckte Ressourcen

https://pubs.usgs.gov/pp/1773/ Digital and/or Hardcopy Resources 1900 2004 ground condition PP1773, Kapitel C Hydrogeologische Einheiten und Modellschichten für das Untersuchungsgebiet in South Carolina Clarke, JS, Brooks, R. und Faye, RE

Hydrogeologie der Dublin- und Midville-Aquifer-Systeme von Ost-Zentral-Georgia Digital- und/oder Hardcopy-Ressourcen

http://ga.water.usgs.gov/publications/ggs/ic-74 Digital and/or Hardcopy Resources 1985 Veröffentlichungsdatum IC-74 Hydrogeologische Einheiten und Modellschichten für das Untersuchungsgebiet im ost-zentralen Georgia Brooks, R., Clarke, JS und Faye, RE

Hydrogeologie des Gordon Aquifer Systems of East-Central Georgia Digital and/oder Hardcopy Resources

http://ga.water.usgs.gov/publications/ggs/ic-75/ Digital and/or Hardcopy Resources 1985 Veröffentlichungsdatum IC-75 Hydrogeologische Einheiten und Modellschichten für das Untersuchungsgebiet in Ost-Zentralgeorgien

Siehe PP 1773 für Methoden zur Generierung hydrogeologischer Einheitenausdehnungen.

Erhaltene Daten aus Projektarchiv, Projektion South Carolina State Plane, NAD83 Datum: 1. Polygonabdeckung der Ebenenausdehnungen (ursprünglich CAD-Ausgabe) 2. Alle in Shapefile-Format konvertiert 3. Überflüssige Attribute für die Archivierung entfernt 4. Attribut NAME gefüllt mit Layer with Nummer 5. Attribut STRAT hinzugefügt mit Name(n) des Grundwasserleiters oder der einschließenden Einheit, wie in North Carolina/South Carolina verwendet used

State-Plane-Koordinatensystem 1983 3900 32,5 34,8333333333 -81,0 31,83333333333 1999996,0 0,0


Georeferenzierung eines Rasters auf einen Vektor

Bei der Georeferenzierung sollten Sie in Ihren Bildern nach gut definierten Objekten suchen, z. B. Straßenkreuzungen oder Landmerkmalen. Auf diese Weise können Sie sicher sein, dass Sie sowohl im Raster als auch im ausgerichteten Layer auf dieselbe Position verweisen. Sie können ein Raster-Dataset, einen Raster-Layer, der Raster-Funktionen enthält, einen Image-Service oder einen Mosaik-Layer georeferenzieren. Einige der Georeferenzierungsoptionen können hilfreich sein, um Ihre Georeferenzierungsaufgabe zu vereinfachen.

Die Layer-Liste der Werkzeugleiste Georeferenzierung zeigt Raster-Layer und CAD-Layer als gültige Datentypen an. Die Layer müssen sich entweder im selben Koordinatensystem wie der Datenrahmen befinden oder keinen definierten Raumbezug haben.

  1. Fügen Sie in ArcMap die Layer hinzu, die sich in Kartenkoordinaten befinden, und fügen Sie das Raster-Dataset hinzu, das Sie georeferenzieren möchten.

Es empfiehlt sich, die Daten zuerst mit dem Kartenkoordinatensystem hinzuzufügen, sodass Sie das Datenrahmen-Koordinatensystem nicht festlegen müssen.

Es kann hilfreich sein, den Befehl „Ausdehnung verwendet von voller Ausdehnung“ in den Datenrahmeneigenschaften auf Ihren Untersuchungsbereich einzustellen, damit das Werkzeug „Auf volle Ausdehnung zoomen“ automatisch auf die volle Ausdehnung Ihres Untersuchungsbereichs zoomt.

Dadurch wird das Raster-Dataset im selben Bereich wie die Ziel-Layer angezeigt. Sie können auch die Umschalttaste verwenden und drehen Werkzeuge zum Verschieben des Raster-Datasets nach Bedarf. Um alle Datensätze anzuzeigen, passen Sie deren Reihenfolge im Inhaltsverzeichnis an.

Sie können Ihre Links auch im Vergrößerungsfenster oder im Viewer-Fenster hinzufügen. Wenn Sie Polygone als referenzierte Ebene verwenden, können Sie die Symbolleiste Effekte öffnen, um die Transparenz beim Hinzufügen Ihrer Links anzupassen.

Sie benötigen mindestens 1 Verknüpfung für ein Polynom nullter Ordnung, 3 Verknüpfungen für ein Polynom erster Ordnung oder eine Anpassungstransformation, 4 Verknüpfungen für eine projektive Transformation, 6 Verknüpfungen für ein Polynom zweiter Ordnung und 10 Verknüpfungen für ein Polynom dritter Ordnung oder Spline-Transformation.

Sie können den Restfehler für jede Verbindung und den quadratischen Mittelwertfehler (RMS) untersuchen. Wenn Sie mit der Registrierung zufrieden sind, können Sie die Eingabe von Links beenden.

Sie können Ihr Raster-Dataset nach der Georeferenzierung dauerhaft transformieren, indem Sie den Befehl "Rektifizieren" verwenden (klicken Sie auf das Dropdown-Menü "Georeferenzierung" und dann auf "Rektifizieren" ), das Werkzeug "Verkrümmen" oder das Werkzeug "Aus Datei verzerren".

Beim Aktualisieren eines Raster-Layers, eines Image-Services oder eines Mosaik-Layers wird nur der Layer in Ihrem Kartendokument aktualisiert. Die Georeferenzierungsinformationen werden nicht in der Quelle gespeichert.


Bestimmen Sie unbekannte Koordinatensysteme: Basierend auf Werten der räumlichen Ausdehnung

Ich entwickle derzeit ein Python-Modul, um Shapefiles und CAD-Zeichnungen (DWG, DXF) mit bekannten und unbekannten Koordinatensystemen neu zu projizieren. Der folgende Code befasst sich mit den bekannten Koordinatensystemen der Shapefiles, von denen es 1500 Shapefiles und 3600 CAD-Zeichnungen gibt. Mir ist bewusst, dass es Einschränkungen gibt, das richtige Koordinatensystem von ShapefilesCAD-Zeichnungen mit unbekannten Koordinatensystemen zu identifizieren. Ich habe folgende zwei Artikel gefunden:

Projizierte Koordinatensystemtabellen: Datum Arc 1960

Projizierte Koordinatensystemtabellen: Datum WGS 84

Ich hatte gehofft, dass ich durch den Vergleich der räumlichen Ausdehnung der ShapefilesCAD-Zeichnungen mit der Ausdehnung der Koordinatensysteme (d. h. Arc_1960_UTM_Zone_37S, WGS_1984_UTM_Zone_37S) das richtige Koordinatensystem bestimmen kann. Das erste Problem besteht darin, dass der Umfang in Dezimalgrad angegeben ist und aus irgendeinem Grund in WGS 84.

Tabellen des projizierten Koordinatensystems: Werte basieren auf WGS84

Wenn jemand einen Workflow zur Identifizierung unbekannter Koordinatensysteme basierend auf der räumlichen Ausdehnung der ShapefilesFeature-ClassesCAD-Zeichnungen entwickelt hat, wäre ich für jeden Rat dankbar.


Ausdehnungen so anpassen, dass sich CRS-Shapefiles in QGIS überlappen? - Geografisches Informationssystem

Das Ziel der Arbeit mit einer georeferenzierten Zeichnung in Vectorworks besteht darin, den internen Ursprung der Zeichnung und damit ihre Geometrie an einer bestimmten Position auf der Erde auszurichten. Georeferenzierungsoptionen sind komplex, was für Architekturworkflows funktioniert, möglicherweise nicht für einen Planungsworkflow geeignet. Darüber hinaus haben verschiedene Regionen der Erde unterschiedliche geografische Anforderungen. Da die Erde keine perfekte Kugel ist, ist es kein einfacher Prozess, die dreidimensionale Erdoberfläche mit einer zweidimensionalen Karte zu korrelieren. Koordinatensysteme verwenden mathematische Formeln, um die Form der Erde auf eine 2D-Seite zu übertragen. Da bei der Translation Daten verloren gehen, sind für verschiedene Regionen der Erde unterschiedliche Koordinatensysteme vorzuziehen.

Ein richtiges Koordinatensystem hat einen großen Einfluss auf Entfernungsmessungen. Eine Entfernung in der Zeichnung und dieselbe Entfernung auf der Karte stimmen je nach ausgewähltem Koordinatensystem mehr oder weniger überein. Je näher die Übereinstimmung mit der bearbeiteten Region ist, desto genauer und genauer sind die Messungen und Abmessungen.

Wenn die Georeferenzierung aktiviert ist, speichert jeder Design-Layer Informationen, die die Abbildung des kartesischen Koordinatensystems dieses Layers auf ein geografisches Koordinatensystem mit Breiten- und Längengrad oder Rechts- und Hochwert beschreiben. Der Ursprungsversatz bestimmt, wie sich das geografische Koordinatensystem zum internen Ursprung (der Mitte der Zeichnung in der Vectorworks-Datei) verhält. Dieser Versatz kann auf Dokumentebene oder auf Ebenenebene angepasst werden.

Die Parameter Georeferenzierung   richten die Georeferenzierung des Dokuments ein. Normalerweise verfügen alle oder die meisten Design-Layer über dieselben Georeferenzierungsinformationen wie die Dokument-Layer ohne aktivierte Georeferenzierung, und übernehmen automatisch die Einstellungen für die Dokument-Georeferenzierung. Ein georeferenzierter Layer bewegt sich mit wechselnder Geoposition (im Wesentlichen bleibt er auf dem Globus), während ein nicht georeferenzierter Layer immer um den internen Ursprung herum verbleibt. Nachdem die Georeferenzierung des Dokuments festgelegt wurde, können Sie die einzelnen Georeferenzierungseinstellungen eines Layers unterschiedlich aktivieren und anpassen, dies wird jedoch selten benötigt.

Beim Importieren eines Shapefiles kann die Georeferenzierung des Dokuments automatisch an das Koordinatensystem der importierten Datei angepasst werden.

So legen Sie die Georeferenzierung für das Dokument fest:

Das Dialogfeld Dokument-Georeferenzierung wird geöffnet.

Passen Sie den Ursprung und die Ausrichtung der georeferenzierten Daten an die Dokumentkoordinaten an

Diese Optionen beeinflussen die Art und Weise, wie sich die Vectorworks-Zeichnung (interner Ursprung) auf den geografischen Standort bezieht. Sie können die Position der georeferenzierten Daten so verschieben, dass sie der Beziehung des Koordinatensystems zum internen Ursprung der Zeichnung entspricht.

Normalerweise sollten diese Optionen aktiviert und auf einen Standort in der Nähe des Standorts eingestellt werden. Dadurch wird sichergestellt, dass Ihre Zeichenarbeit auf der Site innerhalb der Grenzen der doppelten Genauigkeit liegt, die für genaue mathematische Berechnungen und Operationen in Vectorworks erforderlich sind.

Gibt den Breiten- und Längengrad am Dokumentursprung (0,0) an. Geben Sie Dezimalgrad oder Grad/Minuten/Sekunden (z. B. 39,18, 39° 10' 32” oder 39d 10m 32s) mit einem vorangestellten Minuszeichen (-) ein wenn angemessen

Gibt an, welche Richtung nördlich vom Dokumentursprung ist Verwenden Sie eine Azimut-Grad-Messung (0° ist entlang der Y-Achse, 90° ist entlang der positiven X-Achse)

Wählen Sie ein Koordinatensystem

Wählen Sie eines der Koordinatensysteme aus, basierend auf der geografischen Position der Zeichnung und den Ausdehnungen der Zeichnung zusätzliche Parameter, die erforderlich sind, um das ausgewählte Koordinatensystem im unteren Bereich abzubilden

Verwenden Sie WGS84/Pseudo-Mercator für den weltweiten Einsatz

Dieses globale Koordinatensystem ist für allgemeine Zwecke gut geeignet, aber da es den gesamten Globus anzeigt, kann eine Verzerrung die Genauigkeit beim Abgleich der Vectorworks-Entfernungen und der Kartenentfernungen verringern. Wählen Sie nach Möglichkeit ein regionalspezifischeres Koordinatensystem aus.

Verwenden Sie einen EPSG-Code oder CRS-Namen

Verwendet den EPSG-Code des geodätischen Parameterdatensatzes der International Association of Oil and Gas Producer oder den Namen des Koordinatenreferenzsystems (CRS), um das Koordinatensystem anzugeben. Klicken Sie auf EPSG/CRS Lookup  , um das gewünschte Referenzsystem aus dem Dialogfeld EPSG/CRS Lookup auszuwählen. Das ausgewählte System wird im Dialogfeld "Dokument-Georeferenzierung" angezeigt.

Geben Sie Text aus dem EPSG-Code oder CRS-Namen in das Feld Filter ein, um die Liste zu filtern.

Verwenden Sie eine Projektionsdatei (.prj) auf Ihrem Computer

Wenn Sie eine bestimmte Projektionsdatei für das Georeferenzierungs-Koordinatensystem haben, klicken Sie auf Durchsuchen, um sie auszuwählen

Verwenden Sie ein gemeinsames Koordinatensystem oder geben Sie Well Known Text (WKT) ein.

Wählen Sie eines der Vectorworks-Standardkoordinatensysteme aus, das von GDAL (Geospatial Data Abstraction Library) unterstützt wird, oder wählen Sie die Verwendung der WKT-Markup-Sprache aus. Geben Sie aus der Liste Koordinatensystem   einen oder mehrere Begriffe in das Suchfeld ein, um die Liste zu filtern.

Je nach Auswahl werden im darunter liegenden Bereich weitere Parameter angezeigt.

Zeigt eine Beschreibung der Koordinatensystemdaten an.

Sie können den Text des Koordinatensystems bearbeiten und zur späteren Verwendung speichern. Klicken Sie auf Speichern   oder Löschen  , um ein benutzerdefiniertes Element in der Liste Koordinatensystem hinzuzufügen/zu entfernen. Geben Sie beim Speichern einen eindeutigen CRS-Namen an, wenn Sie die vorhandenen Daten nicht ersetzen möchten.

Klicken Sie auf EPSG/CRS-Lookup, um das gewünschte Referenzsystem aus dem Dialogfeld EPSG/CRS-Lookup auszuwählen.

Wählen Sie das entsprechende Referenzellipsoid aus, das zum Erstellen des Koordinatensystems verwendet wird

Wählen Sie für das Koordinatensystem Universal Transverse Mercator (UTM) die UTM-Gitterzone und die UTM-Halbkugel aus, die für diese Zeichnung gelten

Wählen Sie für ein Zustandsebenen-Koordinatensystem (NAD83 oder NAD27) den Zustand und die Zone aus, die für diese Zeichnung gelten

Geben Sie den Breiten- und Längengrad des Mittelpunkts des ausgewählten Koordinatensystems an

Geben Sie für das transversale Mercator-Koordinatensystem den Skalierungsfaktor entlang der zentralen Länge an, für das stereografische Koordinatensystem den Skalierungsfaktor am Mittelpunkt at

Geben Sie für das Lambert Conformal Conic-Koordinatensystem die beiden Standardparallelen ein, die für das Koordinatensystem verwendet werden sollen

Wählen Sie die entsprechende Option zum Transformieren der Geometrie von Objekten auf dem Layer, wenn sich der Koordinatensystemtyp des Layers ändert:

● Geometrie dieses Layers nicht transformieren:  Layer-Objekte werden nicht transformiert, wenn das Koordinatensystem des Layers geändert wird.

● Transformieren Sie die gesamte Geometrie einzeln:  Am besten für die Projektion von geografischen Daten wie Gemeindegrenzen, Straßen, Flüssen und Standorten auf einer Karte.

Diese Option stellt sicher, dass die Zeichnungselemente mit der Georeferenzierung des Koordinatensystems des Layers übereinstimmen. Wenn Sie manuelle Anpassungen des Breiten-/Längen-Offsets oder des Nordwinkels vorgenommen haben, wählen Sie diese Option.

● Verschieben und drehen Sie die gesamte Geometrie als Gruppe:  Am besten zum Projizieren von nicht-geografischen Zeichnungsobjekten, wie z. B. Architekturplänen.

Transformieren Sie Bilder so, dass sie dem Koordinatensystem des Layers entsprechen

Stellen Sie bei importierten georeferenzierten Bildern sicher, dass sich die Bildpixel auf dem richtigen Breiten- und Längengrad befinden, um dem Koordinatensystem des Layers zu entsprechen, wenn ausgewählt, Bildobjekte werden auch transformiert, wenn der Transformationsmodus geändert wird.

Verwirft alle vorgenommenen Änderungen und setzt alle Felder auf ihre früheren Werte zurück

2. Aktivieren Sie die Georeferenzierung für jeden Design-Layer, der dies erfordert, wie unter Festlegen der Georeferenzierung für Design-Layer beschrieben.


Textdateien mit Trennzeichen können als Punktdaten importiert werden, sofern sie Koordinatenwerte enthalten. MAPublisher unterstützt die folgenden Trennzeichen zwischen Datenwerten: Komma, Return, Zeilenende und Tabulator. Geben Sie beim Importieren von Textdaten mit Trennzeichen die Zeichenkodierung und das Format an: Projizierte Einheiten, Dezimalgrad, Grad mit Trennzeichen Minuten Sekunden, Grad.Minuten, Grad.MinutesSekunden, Gepacktes DMS mit Dezimalstellen und Gepacktes DMS. Geben Sie im Abschnitt Koordinatenspalten die Felder in der Quelldatei an, die die entsprechenden Koordinatenwerte enthalten. Die Beschriftungen der Dropdown-Liste zeigen Längengrad und Breitengrad oder X und Y an, je nach ausgewähltem Format. Wenn die erste Zeile der Textdatei Spaltenüberschriften enthält, aktivieren Sie das Kontrollkästchen Erste Zeile als Kopfzeile verwenden (MAPublisher erkennt dies in der Regel automatisch). Wenn die Datei keine Spaltenüberschriften zu enthalten scheint, weist MAPublisher den ausgewählten Feldern beim Importieren Standardüberschriftennamen (Spalte1 und Spalte2) zu. Weisen Sie den Attributspalten im Rahmen Schema auswählen Datenschematypen zu.

Die meisten Formate verfügen über eine Einstellung für die Zeichencodierung. Als Attribute beim Import werden erweiterte und internationale Zeichensätze unterstützt. Um einen für Ihren ausgewählten Datensatz geeigneten Zeichencodec zuzuweisen, wählen Sie den entsprechenden Wert aus der Dropdown-Liste Kodierung (in der ersten Spalte wird der Name angezeigt und in der zweiten Spalte wird der Code oder der alternative Name angezeigt, falls verfügbar).


Um ein Projekt zu unterstützen, müssen wir häufig räumliche Datensätze aus verschiedenen Quellen und/oder Daten sammeln, die unterschiedliche räumliche Ausdehnungen abdecken . Geodaten aus unterschiedlichen Quellen, die unterschiedliche Ausdehnungen abdecken, befinden sich oft in unterschiedlichen Koordinatenreferenzsystemen (CRS).

Einige Gründe dafür, dass sich Daten in unterschiedlichen CRS befinden, sind:

  1. Die Daten werden in einer bestimmten CRS-Konvention gespeichert, die vom Datenlieferanten verwendet wird, der eine Bundesbehörde oder ein Landesplanungsamt sein kann.
  2. Die Daten werden in einem bestimmten CRS gespeichert, das auf eine Region zugeschnitten ist. Zum Beispiel ziehen es viele Staaten vor, a . zu verwenden Staatsebene Projektion, die für diesen Zustand angepasst ist.

Sehen Sie sich dieses kurze Video an, in dem gezeigt wird, wie Kartenprojektionen Kontinente proportional größer oder kleiner erscheinen lassen, als sie tatsächlich sind!

In diesem Tutorial lernen wir, wie Sie räumliche Daten in verschiedenen Projektionen identifizieren und verwalten. Wir werden lernen, wie Sie die Daten neu projizieren, sodass sie sich in derselben Projektion befinden, um das Plotten / Mapping zu unterstützen. Beachten Sie, dass diese Kenntnisse auch für jede Geoverarbeitung / räumliche Analyse erforderlich sind. Die Daten müssen sich im selben CRS befinden, um genaue Ergebnisse zu gewährleisten.


Unser Array hat einen nodata-Wert, -9999.0 , wodurch die Farbe unseres Diagramms über einen zu großen Bereich gestreckt wird. Wir möchten nur gültige Werte anzeigen, daher können wir vor dem Plotten die nodata-Werte mit der where()-Funktion und dem .rio.nodata-Attribut unseres DataArray herausfiltern.

Wenn wir Terrain_HARV_UTM18 in einer Datei gespeichert und dann mit dem masked=True Argument von open_rasterio eingelesen hätten, wäre der nodata-Wert des Rasters maskiert und wir müssten die Funktion where() nicht verwenden, um die Maskierung vor dem Plotten durchzuführen.


Zuschneiden eines Rasters mit Vektorausdehnung

Wir können die Funktion crop() verwenden, um ein Raster auf die Ausdehnung eines anderen räumlichen Objekts zuzuschneiden. Dazu müssen wir das zu beschneidende Raster und das zum Zuschneiden des Rasters verwendete räumliche Objekt angeben. R verwendet die Ausdehnung des räumlichen Objekts als Zuschneidegrenze.

Um dies zu veranschaulichen, beschneiden wir das Canopy Height Model (CHM), um nur den Bereich von Interesse (AOI) einzuschließen. Beginnen wir damit, den gesamten Umfang der CHM-Daten zu zeichnen und zu überlagern, wo der AOI darin liegt. Die Grenzen des AOI werden blau gefärbt, und wir verwenden fill = NA, um den Bereich transparent zu machen.

``` ggplot() + geom_raster(data = CHM_HARV_df, aes(x = x, y = y, fill = HARV_chmCrop)) + scale_fill_gradientn(name = „Canopy Height“, colors = Terrain.colors(10)) + geom_sf(data = aoi_boundary_HARV , Farbe = „blau“, Füllung = NA) + coord_sf()

Jetzt können wir die zugeschnittenen CHM-Daten zusammen mit einem Begrenzungsrahmen darstellen, der die vollständige CHM-Ausdehnung zeigt. Denken Sie jedoch daran, dass es sich um Rasterdaten handelt, die wir in einen Datenrahmen konvertieren müssen, um mit ggplot zu plotten. Um die Boundary Box von CHM zu erhalten, extrahiert st_bbox() die 4 Ecken des Rechtecks, die alle in diesem Objekt enthaltenen Features umfassen. st_as_sfc() wandelt diese 4 Koordinaten in ein Polygon um, das wir plotten können:

``` CHM_HARV_Cropped_df <- as.data.frame(CHM_HARV_Cropped, xy = TRUE)

ggplot() + geom_sf(data = st_as_sfc(st_bbox(CHM_HARV)), fill = „grün“, color = „grün“, alpha = .2) +
geom_raster(data = CHM_HARV_Cropped_df, aes(x = x, y = y, fill = HARV_chmCrop)) + scale_fill_gradientn(name = „Canopy Height“, colors = Terrain.colors(10)) + coord_sf()

Wir können uns den Umfang aller unserer anderen Objekte für diese Feldstelle ansehen.

``` st_bbox(CHM_HARV) st_bbox(CHM_HARV_Cropped) st_bbox(aoi_boundary_HARV) st_bbox(plot_locations_sp_HARV)


Symbologie der Layer ändern

Nehmen Sie sich einen Moment Zeit, um die beiden Datenschichten zu erkunden. Öffnen Sie ihre Attributtabellen, um herauszufinden, welche Art von Daten sie enthalten – und welche Felder Sie möglicherweise zur Visualisierung der Daten verwenden möchten.

Denken Sie daran, dass kategoriale Daten (Codes, Gruppen) mit der Funktion Kategorisiert symbolisiert werden können, während Messdaten (Zahlen, Entfernungen, Bevölkerung) mit der Funktion Abgestuft visualisiert werden sollten.

Visualisieren Sie die HOLC-Grade

Schauen Sie sich die HOLC-Karten auf der Mapping Inequality-Website noch einmal an. Beachten Sie das Farbschema der Kartographen:

Klicken Sie im Bedienfeld „Layer“ mit der rechten Maustaste auf den Layer „cartodb-query“ und klicken Sie auf „Eigenschaften“ und navigieren Sie dann zur Registerkarte „Symbologie“. Im Moment ist der Symbologietyp auf set eingestellt Einzelsymbol. Ändere dies in Kategorisiert. Wählen Sie für das Feld Wert das Feld aus, das Sie symbolisieren möchten, und klicken Sie dann auf die Schaltfläche Klassifizieren in der unteren linken Ecke, um alle möglichen Werte zu laden.

Ändern Sie die Farbe für jeden Grad, damit er den HOLC-Karten entspricht. Dazu gibt es mehrere Möglichkeiten:

  1. Markieren Sie das Element im Symbologie-Fenster und klicken Sie dann auf den kleinen Dropdown-Pfeil rechts neben dem Farbstreifen darüber.
  2. Doppelklicken Sie auf jedes kleine Farbfeld, um die Symbolauswahl aufzurufen. Stellen Sie sicher, dass Füllen ausgewählt ist, und klicken Sie dann auf den kleinen Dropdown-Pfeil rechts neben dem Farbstreifen.
  3. Wählen Sie alternativ Einfache Füllung, wenn Sie den Fülltyp anpassen, die Farbe und Breite des Umrisses (Strichs) ändern möchten usw.

Sie können die Transparenz des Layers anpassen, indem Sie Layer-Rendering im Symbologie-Fenster erweitern und dann den Schieberegler Deckkraft anpassen. 0 % ist völlig transparent und 100 % ist völlig undurchsichtig.

Entfernen Sie schließlich alle anderen Werte, indem Sie sie markieren und auf das rote Minussymbol klicken, dann klicken Sie auf OK .

Visualisieren Sie die Volkszählungsdaten

Die andere Datenschicht hat viele, viele weitere Felder, die zur Visualisierung verwendet werden könnten. Tatsächlich ist die Attributtabelle etwas überwältigend. Es gibt Dutzende von Feldern, aber alle haben codierte Titel, die für sich genommen nicht viel Sinn machen.

Um die Tabelle zu entschlüsseln, werfen Sie einen Blick auf die Metadaten zur Datei im GitHub-Repository der MGGG. (Beachten Sie, dass die erste Gruppe von Metadatenfeldern für ein anderes Shapefile ist. Scrollen Sie weiter, bis Sie die Informationen zum Shapefile 2012-2016 sehen). Welches der Felder wäre Ihrer Meinung nach interessant zu kartieren?

Lassen Sie uns zunächst die demografischen Informationen untersuchen, die diese Datei enthält. Nehmen wir an, wir sind daran interessiert, den Prozentsatz der Bevölkerung zu kartieren, der sich als nicht-hispanisches weiß. Eine schnelle Möglichkeit wäre, das Feld NH_WHITE: Weiße, nicht-hispanische Bevölkerung in der Volkszählung 2010 durch das Feld TOTPOP: Gesamtbevölkerung in der Volkszählung 2010 zu teilen.

Führen Sie diese Schritte aus, um ein neues Feld in der Attributtabelle zu berechnen:

  • Öffnen Sie die Attributtabelle für den Layer MA_prencts12_16
  • Klicken Sie auf das Stiftsymbol, um mit der Bearbeitung zu beginnen
  • Klicken Sie auf die Schaltfläche Feldrechner öffnen
  • Stellen Sie sicher, dass Neues Feld erstellen aktiviert ist
  • Ändern Sie den Namen des Ausgabefelds in Per_Whi
  • Ändern Sie den Ausgabefeldtyp in Dezimalzahl (real)
  • Geben Sie im Feld Ausdruck diese Formel ein: „NH_WEISS“ / „TOTPOP“
  • OK klicken

Stellen Sie sicher, dass Sie Ihre Änderungen speichern und deaktivieren Sie das Stiftsymbol, um die Bearbeitung zu beenden.

Jetzt können Sie dieses neue Feld verwenden, um die Symbologie für den Layer zu ändern. Öffnen Sie das Symbologie-Fenster und ändern Sie den Symbologietyp in Absolvent. Wählen Sie für Wert das gerade erstellte neue Feld aus und klicken Sie dann auf die Schaltfläche Klassifizieren in der unteren linken Ecke, um alle möglichen Werte zu laden. Wählen Sie einen Farbverlauf, der Ihnen am besten gefällt.

Hinweis: Sie können auch mit verschiedenen Modi experimentieren, die die Gruppierung der Daten ändern. Versuchen Sie, den Modus zu ändern, und klicken Sie dann auf Übernehmen, um zu sehen, wie sich dies auf die Datenvisualisierung auswirkt.

Wenn Sie fertig sind, klicken Sie auf OK .

Visualisierung der HOLC-Kategorien

Visualisierung demografischer Daten

Erkunde die zwei Schichten

Schließlich! Dies ist der Teil, in dem die Magie beginnt. Verbringen Sie einige Zeit damit, die räumliche Beziehung zwischen diesen beiden Schichten zu untersuchen.

  • Welche Muster scheinen aufzufallen?
  • Gibt es offensichtliche Cluster in Bezug auf HOLC-Grade oder demografische Zusammensetzung? Irgendwelche Ausreißer?
  • Wie würden Sie den räumlichen Zusammenhang zwischen der historischen HOLC-Bezeichnung der Quartiere in diesen Städten und der demografischen Zusammensetzung ihrer heutigen Bewohner beschreiben?
  • Welche Hypothesen könnten Sie weiter untersuchen?

Um mehr über das Mapping-Programm der Home Owners’ Loan Corporation und die langfristigen Auswirkungen der Rassentrennung zu erfahren, die in Städten in den USA resultierte, lesen Sie die Einführung des Mapping Inequality-Projekts für einen Überblick und Links zu weiterführenden Informationen.

Vergessen Sie nicht, Ihre Karte zu speichern!


Schau das Video: How to change Coordinate System CRS of Shapefile in QGIS #Projections and CRS