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

QGIS: Wie erklär ich’s meinem Kind?

Screenshot (Quelle [2])

Wahrscheinlich wissen wir alle, was QGIS (und GIS) ist, wenn wir aber dann versuchen, es Leuten zu erzählen, die wenig bis keine Ahnung davon haben, wird es meistens schwer, wir stehen vor dem typischen “Wie erklär ich’s meinem Kind?”-Problem. Alasdair Rae (@undertheraedar) hat es getan und es ist ein wunderbares, leicht verständliches PDF [1] entstanden; Nerdwörter werden vermieden! Eine überarbeitete Version findet Ihr in [2], die Weiterverwendung wird ausdrücklich erlaubt. Als QGIS-Projekt steht Euch auch alles unter [3] zur Verfügung, Downloads via Dropbox [4] und GoogleDrive [5].
Danke Alasdair, großartig!

Screenshot (Quelle [3])
Screenshot: Das Ganze als deutsche Version (Dank an Torsten Wolff), Download: PDF (353 KByte) [6]

Hier der initiale Tweet [1]:

[1] … https://twitter.com/undertheraedar/status/1511340898839408652
[2] … https://twitter.com/undertheraedar/status/1552313734592450565
[3] … https://twitter.com/undertheraedar/status/1552606985463398400
[4] … https://www.dropbox.com/sh/vag4yakoshmqna8/AAAGo-RonMpJxexwBukqo7-ya?dl=0
[5] … https://drive.google.com/drive/folders/11eYJ45vX147Y0NWq7Zg88GB_IrevHJ6C
[6] … https://geoobserver.de/wp-content/uploads/2022/08/qgis_what-is-qgis_german.pdf

QGIS-Tipp: Extraktion gemeinsamer Grenzen und einseitiger Puffer

Wer wissen möchte, wie man “eine gemeinsame Grenze zwischen benachbarten Verwaltungsregionen extrahiert und eine Pufferzone auf jeder Seite entlang der Grenze berechnet”, sollte sich das QGIS-Video “Extracting Shared Border and One-Sided Buffers in QGIS” [1] von Ujaval Gandhi (@spatialthoughts) anschauen.

Hier der Original-Tweet [2]:

[1] … https://www.youtube.com/watch?v=YRcrxB3zNhw
[2] … https://twitter.com/spatialthoughts/status/1544609803691266048

QGIS-Tipp: LineLabels korrigiert

Wer kennt das nicht? Linienbeschriftungen haben so Ihre Eigendynamik, das GIS plaziert den Text (irgendwo) auf der Linie, fachlich richtig, optisch eher wenig gefällig. Welche Möglichkeiten es gibt, die Beschriftung zu verbessern, zeigt der Artikel “Turn line labels in QGIS to a new point shape layer” [1]. Mein Favorit ist in den folgenden Screenshots (Bildquelle [1]).

Hier der Original-Tweet [1]:

[1] … https://gis.stackexchange.com/questions/256662/turn-line-labels-in-qgis-to-a-new-point-shape-layer
[2] … https://twitter.com/vanuatugis/status/1544109453036834816