[FOSSGIS-Talk] QGIS Beginner mit Abschlussarbeit - Vektorgitter Analyse?

Maximilian Krambach Maximilian.Krambach at gmx.de
Mi Mär 16 15:26:59 CET 2016


   Hallo,
   Wenn du eh zellenweise rechnest, würde das vielleicht am einfachsten
   mit Rasterdaten gehen. (per QGIS/python/ mit GDAL) Erstmal die Daten
   mit den Gebäuden rasterisieren (das sollte irgendwie gehen, dass die
   Summe der Geschosszahlen pro Rasterzelle darin auftauchen. Kann aber
   gerade nix ausprobieren). Danach dann mit dem Raster rechnen:

   #python-Skript mit rasterisierter Datei; zu erstellen/ausführen in
   Erweiterungen-- Python-Konsole
   -----berechnung.py:-----
   import gdal, numpy
   #öffne Raster
   raster = gdal.Open("Rasterlayer.tif")
   band = raster.GetRasterBand(1)
   #lese Raster in ein Array
   array = band.ReadAsArray(0, 0, xgroesse, ygroesse).astype(numpy.float)
   array2 = numpy.zeros(??) # ein neues array mit dem output
   #Loope durch die Pixel des Arrays
   for zeile in range(0,zeilenanzahl):
        array2[zeile][0] = Anfangswert für Spalte 0

   for zeile in range(0, zeilenanzahl):
       for spalte in range(1,spaltenanzahl):
           #berechne die Operationen für array2[zeile][spalte].
           #dabei ist der wert "links" array2[zeile][spalte-1]
            array2[zeile][spalte] = array2[zeile][spalte-1] *
   array[zeile][spalte]
   #schreibe das array2 in ein neues Raster. #TODO kann ich nicht
   auswendig, steht aber in der angegebenen Website:
   https://pcjericks.github.io/py-gdalogr-cookbook/index.html

   Grüße, Maximilian


   Gesendet: Mittwoch, 16. März 2016 um 14:48 Uhr
   Von: "Keine Werbung" <dezentech at web.de>
   An: fossgis-talk-liste at fossgis.de
   Betreff: Re: [FOSSGIS-Talk] QGIS Beginner mit Abschlussarbeit -
   Vektorgitter Analyse?
   Rauhigkeiten muß Du gar nicht ermitteln - gibt doch die CORINE-Daten.
   Grüße sendet Lars
   Am 16.03.2016 um 13:59 schrieb Hayley Haydon:
   > Hi!
   >
   > Danke für die Antwort!
   >
   > Es geht auch nicht um eine exakte Analyse, was ja auf Grund der
   Bebauung
   > zusätzlich komplex ist, sondern eben um die quadratkilometerweise
   > Abschätzung.
   >
   > Ich kann Dir gern das Paper schicken, auf das ich mich beziehen würde
   > (Drew et al, Estimating the potential yield of small wind turbines in
   > urban areas - a case study for Greater London, UK 20013).
   > Das Vorgehen ist etwa wie von mir beschrieben: Rauhigkeit ermitteln
   > (darum geht es mir!), von einem Startpunkt (Wettermast) ausgehend
   > kalkulieren, Ergebnisse in Karte zusammenfassen.
   >
   > Grüße
   >
   > H.
   >
   > Am 16.03.2016 um 13:50 schrieb Keine Werbung:
   >> Das Windpotential zu ermitteln, egal wo, ist eine hochkomplexe
   Sache!
   >>
   >> Dafür gibt es kostenpflichtige Spezialsoftware, z.B. WindPRO von
   EMD.
   >>
   >> Da ich mit beidem arbeite, kann ich mir beim besten Willen nicht
   >> vorstellen, wie das mit QGIS funktionieren soll.
   >>
   >> Grüße sendet Lars
   >>
   >>
   >> Am 16.03.2016 um 11:04 schrieb Hayley Haydon:
   >>> Erst einmal Hallo in die Runde!
   >>>
   >>> Ich bin hier gelandet, weil ich als Student der Umwelttechnik
   gerade im
   >>> Rahmen meiner Abschlussarbeit zum ersten Mal in engeren Kontakt mit
   GIS
   >>> komme. Konkret geht es um die Abschätzung des Windenergiepotenzials
   über
   >>> bebautem Gebiet.
   >>> Die Firma nutzt QGIS in Verbindung mit kommunalen Gebäudedaten
   >>> (Attribute wie Gebäudefläche und Geschossanzahl). Leider gibt es
   >>> niemanden, der/die mir bei der Bearbeitung wirklich hilfreich zu
   Seite
   >>> stehen kann. Daher ist dies hier quasi mein letzter Versuch, die
   Aufgabe
   >>> doch noch so wie gedacht zum Abschluss zu bringen - no pressure ;-)
   .
   >>>
   >>> Die Hauptaufgabe besteht darin, mit o.g. Attributen (evtl.
   zuzüglich der
   >>> Geländehöhe aus einem zweiten, noch zu integrierenden Datensatz)
   >>> quadratkilometerweise die Geländerauhigkeit nach gegebenen Formel
   zu
   >>> bestimmen.
   >>>
   >>> Meine Vorstellung ist in etwa:
   >>>
   >>> - Vektorgitter erzeugen
   >>> - Abfrage pro Zelle: Sind Gebäude vorhanden?
   >>> - - höhe = geschosszahl * pauschale
   >>> - - Rauhigkeit = f(Summe_höhe; Summe_fläche)
   >>> - Startwert = Wert aus Messung
   >>> - Zellenweise Berechnung:
   >>> - - Windenergiepotential = f (Windenergiepotential_n-1; Rauhigkeit)
   >>>
   >>> Das Vektorgitter habe ich schon gefunden, mir ist nur überhaupt
   nicht
   >>> klar, wie ich zu einer automatischen, zellenweisen Analyse (z.B.
   mit
   >>> einem Skript) komme. Kann QGIS das überhaupt? Wonach muss ich genau
   >>> suchen? Mir fehlen natürlich die grundlegenden Begriffe - ich will
   mich
   >>> zwar einarbeiten, aber ungern sämtliche geografischen Grundlagen
   >>> erarbeiten müssen, falls es anders ginge. Mit einem Tipp, wonach
   ich
   >>> konkreter suchen muss wäre mir also wahrscheinlich schon geholfen,
   wie
   >>> auch mit Links oder Literatur.
   >>> Auch, falls mein Vorhaben offensichtlich zu groß für eine
   Bachelorarbeit
   >>> sein sollte wäre ich für einen Hinweis dankbar!
   >>>
   >>> Vielen Dank u. Gruß
   >>>
   >>> H. Haydon
   >>>
   >>> --
   >>>
   ....................................................................
   >>> FOSSGIS 2016, Die Konferenz für Open Source GIS mit OpenData und
   >>> OpenStreetMap in Zusammenarbeit mit der AGIT 2016!
   >>> 4.-6. Juli 2016 in Salzburg (3. Juli OpenStreetMap Workshoptag)
   >>> [1]http://www.fossgis.de/konferenz/2016/
   >>>
   >>> AGIT 2016 vom 6.-8. Juli 2016
   >>> [2]http://agit.at/
   >>>
   >>> FOSS4G 2016 Bonn - annual global event of the Open Source
   Geospatial
   >>> Foundation (OSGeo) - 24.-26. August 2016 in Bonn (zusätzlich noch
   >>> FOSS4G Hacking Event und Workshops)
   >>> [3]http://2016.foss4g.org
   >>>
   >>>
   >>> FOSSGIS e.V, der Verein zur Förderung von Freier Software aus dem
   >>> GIS-Bereich und Freier Geodaten!
   >>> [4]http://www.fossgis.de/ [5]https://twitter.com/fossgis_eV
   >>>
   >>>
   ____________________________________________________________________
   >>> FOSSGIS-Talk-Liste mailing list
   >>> FOSSGIS-Talk-Liste at fossgis.de
   >>> [6]https://lists.fossgis.de/mailman/listinfo/fossgis-talk-liste
   >>>
   >> --
   >> ....................................................................
   >> FOSSGIS 2016, Die Konferenz für Open Source GIS mit OpenData und
   >> OpenStreetMap in Zusammenarbeit mit der AGIT 2016!
   >> 4.-6. Juli 2016 in Salzburg (3. Juli OpenStreetMap Workshoptag)
   >> [7]http://www.fossgis.de/konferenz/2016/
   >>
   >> AGIT 2016 vom 6.-8. Juli 2016
   >> [8]http://agit.at/
   >>
   >> FOSS4G 2016 Bonn - annual global event of the Open Source Geospatial
   >> Foundation (OSGeo) - 24.-26. August 2016 in Bonn (zusätzlich noch
   >> FOSS4G Hacking Event und Workshops)
   >> [9]http://2016.foss4g.org
   >>
   >>
   >> FOSSGIS e.V, der Verein zur Förderung von Freier Software aus dem
   >> GIS-Bereich und Freier Geodaten!
   >> [10]http://www.fossgis.de/ [11]https://twitter.com/fossgis_eV
   >>
   >> ____________________________________________________________________
   >> FOSSGIS-Talk-Liste mailing list
   >> FOSSGIS-Talk-Liste at fossgis.de
   >> [12]https://lists.fossgis.de/mailman/listinfo/fossgis-talk-liste
   >>
   >
   > --
   > ....................................................................
   > FOSSGIS 2016, Die Konferenz für Open Source GIS mit OpenData und
   > OpenStreetMap in Zusammenarbeit mit der AGIT 2016!
   > 4.-6. Juli 2016 in Salzburg (3. Juli OpenStreetMap Workshoptag)
   > [13]http://www.fossgis.de/konferenz/2016/
   >
   > AGIT 2016 vom 6.-8. Juli 2016
   > [14]http://agit.at/
   >
   > FOSS4G 2016 Bonn - annual global event of the Open Source Geospatial
   > Foundation (OSGeo) - 24.-26. August 2016 in Bonn (zusätzlich noch
   > FOSS4G Hacking Event und Workshops)
   > [15]http://2016.foss4g.org
   >
   >
   > FOSSGIS e.V, der Verein zur Förderung von Freier Software aus dem
   > GIS-Bereich und Freier Geodaten!
   > [16]http://www.fossgis.de/ [17]https://twitter.com/fossgis_eV
   >
   > ____________________________________________________________________
   > FOSSGIS-Talk-Liste mailing list
   > FOSSGIS-Talk-Liste at fossgis.de
   > [18]https://lists.fossgis.de/mailman/listinfo/fossgis-talk-liste
   >
   --
   ....................................................................
   FOSSGIS 2016, Die Konferenz für Open Source GIS mit OpenData und
   OpenStreetMap in Zusammenarbeit mit der AGIT 2016!
   4.-6. Juli 2016 in Salzburg (3. Juli OpenStreetMap Workshoptag)
   [19]http://www.fossgis.de/konferenz/2016/
   AGIT 2016 vom 6.-8. Juli 2016
   [20]http://agit.at/
   FOSS4G 2016 Bonn - annual global event of the Open Source Geospatial
   Foundation (OSGeo) - 24.-26. August 2016 in Bonn (zusätzlich noch
   FOSS4G Hacking Event und Workshops)
   [21]http://2016.foss4g.org
   FOSSGIS e.V, der Verein zur Förderung von Freier Software aus dem
   GIS-Bereich und Freier Geodaten!
   [22]http://www.fossgis.de/ [23]https://twitter.com/fossgis_eV
   ____________________________________________________________________
   FOSSGIS-Talk-Liste mailing list
   FOSSGIS-Talk-Liste at fossgis.de
   [24]https://lists.fossgis.de/mailman/listinfo/fossgis-talk-liste

References

   1. http://www.fossgis.de/konferenz/2016/
   2. http://agit.at/
   3. http://2016.foss4g.org/
   4. http://www.fossgis.de/
   5. https://twitter.com/fossgis_eV
   6. https://lists.fossgis.de/mailman/listinfo/fossgis-talk-liste
   7. http://www.fossgis.de/konferenz/2016/
   8. http://agit.at/
   9. http://2016.foss4g.org/
  10. http://www.fossgis.de/
  11. https://twitter.com/fossgis_eV
  12. https://lists.fossgis.de/mailman/listinfo/fossgis-talk-liste
  13. http://www.fossgis.de/konferenz/2016/
  14. http://agit.at/
  15. http://2016.foss4g.org/
  16. http://www.fossgis.de/
  17. https://twitter.com/fossgis_eV
  18. https://lists.fossgis.de/mailman/listinfo/fossgis-talk-liste
  19. http://www.fossgis.de/konferenz/2016/
  20. http://agit.at/
  21. http://2016.foss4g.org/
  22. http://www.fossgis.de/
  23. https://twitter.com/fossgis_eV
  24. https://lists.fossgis.de/mailman/listinfo/fossgis-talk-liste