Raster "glätten"
Lisa Römer
lisa.roemer at fischer-teamplan.de
Mo Jun 10 12:19:17 CEST 2024
Hallo zusammen,
ich möchte gerne ein Raster "glätten" Dazu habe ich mit dem Modellbuilder ein Modell erstellt.
Schritt1: mit r.neighbors über "range" die Standartabweichung zur Nachbarzelle aus dem Input-Raster zu erfragen.
Schritt2: Mit dem Rastercalkulator nur die Zellen Ausgeben, die > meiner gewählten Standartabweichung sind.
Schritt3: mit r.Null Nullwerte dafür vergeben
Schritt3: mit r.neighbors den Werten aus r.Null den Median der Nachbarzellen auf das Input-Raster übertragen und so ein geglättetes Raster erstellen.
Ich weiß nicht, ob ich das zu kompliziert angehe. Vielleicht kennt jemand eine einfachere Lösung. Vor allem das Werkzeug r.neighbors funktioniert in meiner Version 3.34.4 nicht. Wenn ich das Tool einzeln ausführe, bekomme ich folgende Fehlermeldung, die ich leider nicht verstehe:
QGIS-Version: 3.34.4-Prizren
QGIS-Codeversion: decd43dc
Qt-Version: 5.15.3
Python-Version: 3.9.18
GDAL-Version: 3.8.4
GEOS-Version: 3.12.1-CAPI-1.18.1
PROJ-Version: Rel. 9.3.1, December 1st, 2023
PDAL-Version: 2.6.0 (git-version: 3fced5)
Algorithmus gestartet um: 2024-06-10T12:16:02
Algorithmus r.neighbors startet.
Eingabeparameter:
{ '-a' : False, '-c' : True, 'GRASS_RASTER_FORMAT_META' : '', 'GRASS_RASTER_FORMAT_OPT' : '', 'GRASS_REGION_CELLSIZE_PARAMETER' : 0, 'GRASS_REGION_PARAMETER' : None, 'gauss' : None, 'input' : 'C:/temp/processing_RYRiDP/9a74de468b554f01b3ed0f86885e4831/OUTPUT.tif', 'method' : 5, 'output' : 'TEMPORARY_OUTPUT', 'quantile' : '', 'selection' : None, 'size' : 3, 'weight' : '' }
g.region n=-1.7976931348623157e+308 s=1.7976931348623157e+308 e=-1.7976931348623157e+308 w=1.7976931348623157e+308 res=100.0
r.neighbors input=rast_6666c128b0da43 selection=rast_6666c2bd0e9855 method="range" size=3 -c output=output9022b35c361642198338ab041d524d1e --overwrite
g.region raster=output9022b35c361642198338ab041d524d1e
r.out.gdal -t -m input="output9022b35c361642198338ab041d524d1e" output="C:\temp\processing_RYRiDP\af6dfec3137c4a0fae17cfbc9375f401\output.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
WARNUNG: Sperren gleichzeitiger Zugriffe auf ein Mapset ist unter Windows nicht möglich.
C:\Users\Public\Desktop\QGIS 3.34.4>chcp 1252 1>NUL
C:\Users\Public\Desktop\QGIS 3.34.4>g.region n=-1.7976931348623157e+308 s=1.7976931348623157e+308 e=-1.7976931348623157e+308 w=1.7976931348623157e+308 res=100.0
FEHLER: North must be larger than South, but -1.79769e+308 (north) <= 1.79769e+308 (south
C:\Users\Public\Desktop\QGIS 3.34.4>r.neighbors input=rast_6666c128b0da43 selection=rast_6666c2bd0e9855 method="range" size=3 -c output=output9022b35c361642198338ab041d524d1e --overwrite
Opening selection map <rast_6666c2bd0e9855>
FEHLER: Rasterkarte <rast_6666c2bd0e9855> wurde nicht gefunden.
C:\Users\Public\Desktop\QGIS 3.34.4>g.region raster=output9022b35c361642198338ab041d524d1e
FEHLER: Rasterkarte <output9022b35c361642198338ab041d524d1e> konnte nicht gefunden werden.
C:\Users\Public\Desktop\QGIS 3.34.4>r.out.gdal -t -m input="output9022b35c361642198338ab041d524d1e" output="C:\temp\processing_RYRiDP\af6dfec3137c4a0fae17cfbc9375f401\output.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
FEHLER: Rasterkarte oder Gruppe <output9022b35c361642198338ab041d524d1e> nicht gefunden.
C:\Users\Public\Desktop\QGIS 3.34.4>exit
Drücken Sie eine beliebige Taste . . .
Execution completed in 1.18 Sekunden
Ergebnisse:
{'output': 'C:\\temp\\processing_RYRiDP\\af6dfec3137c4a0fae17cfbc9375f401\\output.tif'}
Lade Ergebnis Layer
Die folgenden Layer wurden nicht erzeugt.
. C:/temp/processing_RYRiDP/af6dfec3137c4a0fae17cfbc9375f401/output.tif
Im 'Protokoll-Fenster' im QGIS-Hauptfenster sind mehr Informationen zur Ausführung des Algorithmus zu finden.
Viele Grüße,
Lisa Römer
Mehr Informationen über die Mailingliste FOSSGIS-Talk-Liste