Annonce
L'appel à participation est ouvert jusqu'au 19 janvier 2025!
Pages: 1
- Sujet précédent - QGIS 2.6.1 / MNT : Méthode plus proche voisin : lisser des valeurs - Sujet suivant
#1 Mon 31 August 2015 09:31
- Pivès001
- Participant actif
- Date d'inscription: 23 Jan 2008
- Messages: 138
QGIS 2.6.1 / MNT : Méthode plus proche voisin : lisser des valeurs
Bonjour,
j'utilise des MNT SRTM récupérés sur le site USGS (http://earthexplorer.usgs.gov/) et j'ai des valeurs négatives extrêmes. En regardant de plus près cela correspond à des "défauts" de numérisation (des zones à -32767 dans les vallées Alpines principalement).
Comment obtenir un mnt en lissant ces valeurs selon les valeurs voisines ?
J'ai cherché du coté des fonctions resample en essayant d'appliquer des méthodes de plus proche voisin mais sans succès. J'aimerais obtenir quelque chose de continu et plus cohérent du point de vue topo (en interpolant les valeurs par exemple).
Merci par avance.
Hors ligne
#2 Mon 31 August 2015 13:34
- dominique.lys
- Participant assidu
- Date d'inscription: 5 Oct 2006
- Messages: 473
- Site web
Re: QGIS 2.6.1 / MNT : Méthode plus proche voisin : lisser des valeurs
Bonjour,
Voir du coté des outils suivants :
-GRASS r.fillnulls
-GDAL fillnodata
-SAGA close gaps
Dans tous les cas il faut au préalable convertir les valeurs erronées en "nodata". Pour cela vous pouvez utiliser la calculatrice raster de SAGA avec une formule du type:
Code:
ifelse(lt(a,0),-99999,a)
-99999 étant considéré comme nodata par SAGA, cf doc http://www.saga-gis.org/saga_module_doc … lus_1.html et http://sourceforge.net/p/saga-gis/wiki/ … 20library/
Hors ligne
#3 Tue 01 September 2015 09:04
- Pivès001
- Participant actif
- Date d'inscription: 23 Jan 2008
- Messages: 138
Re: QGIS 2.6.1 / MNT : Méthode plus proche voisin : lisser des valeurs
Bonjour et merci de cette réponse.
Code:
ifelse(lt(a,0),-99999,a)
Ce 1er point a très bien fonctionné pour les valeurs inférieures à 0. J'ai bien du "no data" sur mes zones aux valeurs marginales.
Voir du coté des outils suivants :
-GRASS r.fillnulls
-GDAL fillnodata
-SAGA close gaps
La méthode SAGA close gaps s'est révélée être la plus rapide et efficace. J'ai traité par lot mes 129 dalles SRTM en une 30aine de minutes pour un résultat tout a fait satisfaisant visuellement.
Encore merci pour tout.
EDIT
Sur la question ce doc est très intéressant (p3 Void filling):
https://grass.osgeo.org/newsletter/GRASSNews_vol3.pdf
Dernière modification par Pivès001 (Wed 02 September 2015 14:23)
Hors ligne
Pages: 1
- Sujet précédent - QGIS 2.6.1 / MNT : Méthode plus proche voisin : lisser des valeurs - Sujet suivant
