[FOSSGIS-Talk] Problem mit r.lake in GRASS 7.6.1, eingebunden in QGIS Desktop 3.6.1

[Markus Neteler]

Hallo Dirk,

On Wed, Apr 3, 2019 at 3:36 PM Dirk Müller <dm at umweltplan.de> wrote:
>
> Hallo in die Runde,
>
> ich habe mit der GRASS-Funktion r.lake (GRASS 7.6.1, eingebunden in QGIS Desktop
> 3.6.1) folgendes Problem, welches mutmaßlich mit der vom Tool intern verwendeten
> Projektion zusammenhängt:
>
> Im QGIS-Projekt verwende ich als Koordinatenbezug EPSG 5650, auch das Höhenmodell
> und die als Quellpunkt angegebene Koordinate befinden sich in der gleichen Projektion. Bei
> der Berechnung werden einige Probleme protokolliert (anliegend). Das Ergebnis liefert eine
> plausible "Seefläche", allerdings mit einem Versatz von ca. 7.870 m (gemessen) nach
> Westen und geringer Abweichung von 6m nach Süden.
>
> Die Berechnung mit denselben Daten mit GRASS 7.4.1, eingebunden in QGIS Desktop
> 2.18.24 (hier ist die Funktion noch als r.lake.coords implementiert) führt dagegen zu einem
> korrekten Ergebnis (Protokoll und Eingabe-Dialog ebenfalls anliegend).
>
> Ich habe schon einiges probiert, leider alles ohne Erfolg. Zugegebenermaßen habe ich
> bisher nicht allzuviel Erfahrung mit GRASS. Für Tipps wäre ich äußerst dankbar!

Hier einige Anmerkungen:
- der Vergleich der Log-Ausgaben sehe ich keine großen Unterschiede
(außer, dass der Seed point leicht abweicht)
- was "[Freilassen, um min. Abdeckungsausdehnung zu nutzen]" tut,
weiss ich nicht, aber es sollte sicherlich so eingestellt werden, dass
die Berechnungsregion sich auf die Rasterkarte bezieht bzw. v.a. die
Pixelgeometrie korrekt ist (aktuell wurde die "lakes" Karte mit 4,99xx
m und nicht 5.0 m Pixelauflösung berechnet).
- ist deine QGIS 3 vielleicht so eingestellt, dass es "Reprojection on
the fly" durchführt, derweil die QGIS 2 Version das nicht tut?

Viele Grüße,
Markus


--
Markus Neteler, PhD
https://www.mundialis.de - free data with free software
https://grass.osgeo.org
https://courses.neteler.org/blog