QGIS-Tipp: DotMaps – Punktdichtekarten

Animation: Einwohnerdaten im QGIS als Choropleten und DotMaps mit 10 und 100 Elementen pro Punkt

Für die speziell Symbolisierungsart der Punktdichtekarten (DotMaps) gibt es im QGIS ein Plugin “Dot Map” [1], in der Beschreibung des Plugins heißt es:

“Eine Punktdichtekarte ist eine Art von thematischer Karte, die Punkte oder andere Symbole auf der Karte verwendet, um die Werte eines oder mehrerer numerischer Datenfelder darzustellen. Jeder Punkt auf einer Punktdichtekarte repräsentiert eine bestimmte Menge an Daten. QGIS hat keine Symbologie für diese Art von Daten. Aus diesem Grund kann dieses Plugin einen Punktdichte-Layer aus einem Polygon-Layer mit einem Integer-Feld mit den Daten erstellen, die der Benutzer in Punkte umwandeln möchte.” [1].

Ich habe dieses Plugin mal kurz mit Einwohnerdaten aus dem Open Data Portal der Stadt Halle (Saale) [2] mit 10 und 100 Elementen pro Punkt getestet, hier die Ergebnisse:

Noch weiter gehen die DotGridMaps, Details dazu findet Ihr in [3] und [4].

Alternativ kann in QGIS ab 3.16 auch die eingebaute Funktion “Vektor / Forschungswerkzeuge / Zufällige Punkte in Polygonen…” genutzt werden, vgl. auch “Making a Dot Density Map in QGIS” [5]. Mir scheint das o. g. Plugin jedoch einfach in der Bedienung.

[1] … https://plugins.qgis.org/plugins/DotMap/
[2] … https://www.halle.de/de/Verwaltung/Online-Angebote/Offene-Verwaltungsdaten/Mit-Kartenbezug/index.aspx?ID=fb7c7826-c2e1-82f9-395c-bf3c65da4786
[3] … https://www.maproomblog.com/2023/02/in-praise-of-dot-grid-maps/
[4] … https://medium.com/earthrisemedia/we-need-more-dot-grid-maps-7e4374f79acd
[5] … https://www.youtube.com/watch?v=UlG9kwsbnQQ

QGIS-Tipp: Die Unfallschwerpunkte @Youtube

Screenshot: „Heatmap zu Verkehrsunfällen in Chemnitz | QGIS“ auf Youtube [3]

Am Freitag hatte ich Euch im QGIS-Tipp “Halle und die Unfallschwerpunkte” [1] vorgestellt und es gab auch gleich einige Rückfragen. Kaum wollte ich anfangen, diese zu beantworten, teilte mir PyQGIS alias Ivo Partschefeld via Twitter mit “Da musste ich heute noch ein eigenes Video machen.” [2]. So muss Netzwerken! Das Youtube-Video [3]: Interessant, kurz(weilig), treffend, anschaulich, empfehlenswert. Danke Ivo!

[1] … https://geoobserver.de/2023/02/24/qgis-tipp-halle-und-die-unfallschwerpunkte/
[2] … https://twitter.com/PyQgis/status/1629234155379691523?s=20
[3] … https://www.youtube.com/watch?v=kme_G718bwA&feature=youtu.be

QGIS-Tipp: Datenquelle reparieren

Immer wieder kommt es vor, dass QGIS in einem Projekt nicht mehr auf die Datenquelle zurückgreifen kann, vielleicht hat irgend jemand einen Pfad oder Dateinamen umbenannt. Wie Ihr das im QGIS recht einfach fixen könnt, seht Ihr im Youtube-Video “QGIS 3.16 Datenmanagement – Datenquelle reparieren” [1] von Marshal Mappers. “Auto-Finden” kann gut weiter helfen.

[1] … https://www.youtube.com/watch?v=jZZvqXW6nSI

QGIS-Tipp: Heatmaps oder die Apotheken von Halle

Also zuerst, ich mag diese Heatmaps! Heatmap steht für “Wärmebild”, nie war es einfacher, komplexe Inhalte grafisch so gut zu visualisieren, aber man sollte beim Einsatz unbedingt Einiges beachten:

Animation: Die Apotheken von Halle (Saale) mit 1000m, 2000m und 3000m Radius visualisiert

Pro Heatmaps:

  • gute, schnell erfassbare Aussage aus sehr großer Datenmengen bzgl. Dichte der Daten
  • lenkt Aufmerksamtkeit
  • für qualitative Interpretation OK

Contra Heatmaps:

  • Kein/wenig/nicht sichtbarer Bezug zu Originaldaten – mitunter fragt sich der Benutzer, wie es zu dieser Dichte kommt – quantitative Aussagen können nicht gemacht werden
  • Unscharfes Modell durch den (kritischen) Wert des Radius*. Wie kommt der zustande? Reale Werte durch Vorschrift (wenn ja, welche?), Ziel eine schönen Grafik, Zufall oder Defaulteinstellung, Screen- oder Realeinheiten?

Meine Empfehlung:

  • mind. Radius angeben und begründen
  • Originaldaten einblenden, wenigstens schaltbar

Weitere Details zu Heatmaps findet Ihr in “Was ist eine Heatmap?” [1]. Manche würden Heatmaps sogar eher abschaffen [3].
Der Radius* ist in der QGIS-Dokumentation [2] wie folgt erklärt:

* … Radius: Wird verwendet um den Heatmap Suchradius (oder Kernbandbreite) in Metern oder Karteneinheiten anzugeben. Der Radius gibt den Abstand um einen Punkt ab dem der Einfluss des Punktes spürbar wird an. Größere Werte haben eine stärkere Glättung zur Folge, kleinere Werte können aber feinere Details und eine Abwechslung in der Punktdichte zeigen.” [1]

Ich habe mal versucht, das Ganze mit den Apothekenstandorten von Halle (Saale) zu verdeutlichen. Die Daten stammen aus dem OSM-Projekt uns wurden mittels QuickOSM-Plugin ins QGIS übernommen. Dort dann mit der Symbolisierung “Heatmap” für mit 1000m, 2000m und 3000m Radius visualisiert. Welcher der drei Radien verdeutlicht nun die hallesche Apothekendichte am besten? In jedem Fall bemerkenswert die hohe Dichte am Markt und am Reileck – mindestens drei Apotheken auf engstem Raum.

Ihr könnt gern in den Kommentaren Eure Erfahrungen zu Heatmaps ergänzen …

[1] … https://www.omt.de/lexikon/heatmap/#worin-liegt-der-nutzen-einer-heatmap-analyse-von-webseiten-dank-w%C3%A4rmebild
[2] … https://docs.qgis.org/2.14/de/docs/user_manual/plugins/plugins_heatmap.html
[3] … https://geoobserver.de/2019/09/06/heatmaps-pro-contra/

QGIS-Tipp: “QGIS Output Resolution” by Jess Zimmerman

Screenshot: Youtube-Video “QGIS Output Resolution” [1] von Jess Zimmerman

Wie kommt man von 30m Bodenauflösung zu ausreichend guter Druckauflösung? Das interessante Youtube-Video “QGIS Output Resolution” [1] von Jess Zimmerman erklärt Euch …

“… den Unterschied zwischen der Abtastung und der Ausgabeauflösung von schattierten Reliefs. Es beschreibt zwei Methoden zur Steuerung der Ausgabeauflösung von aus QGIS exportierten Dateien, um die gewünschte Ausgabeauflösung für Druckprojekte zu erzeugen.” [1]

[1] … https://www.youtube.com/watch?v=nbzXsEkBvew

KomGIS+: Deutsche Gründlichkeit & die ALK-ALB-Flächendifferenz?

Dass die Amtliche Liegenschaftskarte und das Automatisierte Liegenschaftsbuch z. T. erhebliche Abweichungen in den Flächenangaben haben ist vermutlich vielen hier bekannt. Wir haben in unserem KomGIS+ [1] aus diesem Grund extra ein Thema eingeführt, welches bei jeder neue ALKIS-Datenübernahme diese Differenz für den Nutzer ausweist, vgl. folgende Screenshots.

Screenshot: ALK-ALB-Flächendifferenz im KomGIS+ (Identifikation)
Screenshot: ALK-ALB-Flächendifferenz im KomGIS+ (Sachdaten-Tabelle)

Wer mehr zu den Hintergründen wissen möchte, dem sei das Youtube-Video “Wie können Flächenfehler im Grundbuch/Kataster entstehen? | QGIS” von pyQGIS [2] empfohlen.

[1] … https://itc-halle.de/loesungen/geoinformationssysteme/KomGIS
[2] … https://www.youtube.com/watch?v=qJxk9bG3qXA

Topologie – die heilige Kuh?!

Screenshot: “TOPOLOGY: A BIG WORD WITH BIG IMPLICATIONS FOR YOUR DATA” von Adam Thomas (Quelle [1])

Früher in den 90ern, ich habe 1991 GIS gelernt, war Topologie in Geodaten Standard, ein hohes Gut, eben eine heilige Kuh! Auch beim Marktführer, aus guten Gründen! Clean & Build nach jeder Editiersitzung waren Pflicht und es war gut so. Irgendwann ging die Pflicht leider verloren, spätestens bei der Einführung des ESRI-Shape-Files, für mich der größte Architekturfehler diese Formates, welches durchaus seine Verdienste hatte [6]. Das Topologie und topologisch korrekte Daten immer noch unendlich wichtig sind – für Konsistenz und Präzision – könnt Ihr im Artikel “TOPOLOGY: A BIG WORD WITH BIG IMPLICATIONS FOR YOUR DATA” von Adam Thomas [1] im mapbox-Blog nachlesen.

[1] … https://www.mapbox.com/blog/topology-data-implications
[2] … https://geoobserver.de/2018/07/19/shapefile-die-neverending-story/
[3] … https://geoobserver.de/2017/10/26/shapefile-eine-hassliebe/
[4] … http://switchfromshapefile.org/
[5] … https://geoobserver.de/2016/03/29/theshapefilechallenge-the-winner-is/
[6] … https://geoobserver.de/2017/05/29/25-jahre-esri-shapefile-herzlichen-glueckwunsch/
[7] … https://twitter.com/Mapbox/status/1600928950053703689

QGIS-Tipp: Objekte nach Abstand einfärben

Screenshot: MeinTestprojekt [2], das Zentrum bei der Pauluskirche (Einfärbung nach Abstand siehe Beschreibung)

Innerhalb der #30DayMapChallenge hat eurojam eine coole Visualisierung veröffentlicht, nämlich abhängig von der Entfernung zu einem Punkt die Objekte einzufärben. Ich habe versucht, es nachzuvollziehen und … es klappt. Für alle, die diesbezüglich ein bisschen Hilfe brauchen, hier das Ganze mal etwas untersetzter. Ich habe die Gebäude des OSM-Projektes, hier der Layer “osm_buildings_test_2398” genutzt, dazu einen Layer “center” mit der Spalte “name” und dem Eintrag “Test”. Dann lautet der adaptierte und entscheidende String für die Füllfarbe:

color_hsv(scale_linear(distance( geometry(get_feature('center','name','Test')),$geometry),0,1000,100,0), 100, 100)

Aus dem Layer “center” wird also in der Spalte “name” der Eintrag “Test” gesucht. Zu dessen Geometrie (get_feature)wird die Entfernung (distance) zur aktuellen Gebäude-Geometrie ($geometry) des Datensatzes der Gebäude ermittelt und entsprechend dem aktuellen Abstand eingefärbt (colore_hsv) über einen linearen Farbverlauf (scale_linear). So einfach … Alles klar? Schaut Euch die Syntax der Funktionen get_feature, distance, scale_linear und color_hsv an, es wird klarer, versprochen 😉
Um es einfacher nachvollziehen zu können, könnt Ihr gern mein Testprojekt [2] nutzen.

Animation: Vier Center-Standorte und das eingefärbte Ergebnis

Danke eurojam, hier der Original-Tweet [1]:

[1] … https://twitter.com/27aea577f17d4c4/status/1591081607108390913
[2] … https://www.geoObserver.de/Download/QGIS_Center_Building_Demo1.zip

QGIS-Tipp: Rasterauflösung erhöhen mit GWR

Screenshot (Quelle [1])

Mittels GWR – Geografisch gewichteter Regression kann unter bestimmten Bedingungen die Auflösung eines recht groben Rasters deutlich verfeinert werden. In “‘Increasing’ the resolution of a raster using Downscaling in QGIS“[1] hat Autor franzpc [2] die Vorgehensweise auf GIScrack beschrieben, genutzt wird im QGIS die Saga-Funktion “GWR for Raster Downscaling” [3]. Hochinteressant, ich werde es mal testen.

Hier auch als Video auf Youtube [4] (leider “nur” in Spanisch):

[1] … https://giscrack.com/increasing-the-resolution-of-a-raster-using-downscaling-in-qgis/
[2] … https://giscrack.com/author/franzpc/
[3] … https://saga-gis.sourceforge.io/saga_tool_doc/7.2.0/statistics_regression_14.html
[4] … https://www.youtube.com/watch?v=PBmwy4t0gEQ
[5] … https://acolita.com/aumentar-la-resolucion-de-un-raster-usando-la-tecnica-downscaling-en-qgis/

QGIS-Tipp: Spatial SQL

Screenshot: Tweet mit dem Spatial SQL Spickzettel von Mo Sarwat (Quelle [3])

Wer mit GIS (hier QGIS) arbeitet, kommt bei Datenabfragen irgendwann nicht mehr um SQL und räumliche SQL-Abfragen herum. Einfache Hilfsmittel zum Einstieg findet Ihr z. B. in den Vortragsfolien von Martin Werner [1] und in den Spatial SQL Spickzettel von Mo Sarwat (@MoSarwat) [2], [3]. Viel Spaß beim Ausprobieren.

[1] … https://www.martinwerner.de/teaching/geoinf/files/lecture_04.pdf
[2] … https://twitter.com/MoSarwat/status/1571910560471732225
[3] … https://twitter.com/MoSarwat/status/1571973081111105536