Suite

Comment remplir les trous/valeurs manquantes dans les données XYZ/un DEM avec QGIS ?

Comment remplir les trous/valeurs manquantes dans les données XYZ/un DEM avec QGIS ?


J'ai des données de point (texte) XYZ équidistantes que j'ai utilisées pour créer une grille de points. Cependant, ce n'est pas parfait - il y a quelques trous. J'ai ensuite créé un raster DSM à partir de celui-ci, qui a évidemment aussi les trous. Les trous vont d'un seul point à des zones plus grandes (~50 points ou plus manquants).

Comment colmater ces trous ? J'ai essayé d'interpoler le raster, mais cela réduit la qualité et ne corrige pas les "trous" zéro - ils sont juste légèrement surélevés par la suite. Idéalement, j'aimerais simplement remplir ces valeurs zéro avec les valeurs Z de leurs carrés voisins.

Un autre problème est que l'ensemble de données XYZ n'est pas rectangulaire mais comme dans la troisième image jointe (chemin rouge). Si j'utilise une méthode de remplissage, comment éviter de remplir toute la zone environnante ?


Je suggère de rompre le processus dans les étapes suivantes.

  1. Si les points sont également espacés (ce qui est évidemment le cas), déterminez l'espacement dans les directions x et y
  2. Interpolez un raster initial. Je suggérerais d'utiliser un simple algorithme de pondération de distance inverse.

QGIS Processing Toolbox SAGA Grid - Gridding Distance inverse pondérée

    - Donnez l'espacement comme argument pour le paramètre Cellsize, et
    - spécifiez l'étendue de sortie comme l'étendue de vos points moins la moitié de l'espacement pour xmin/ymin et plus la moitié de l'espacement pour xmax/ymax. Vous obtenez l'étendue de vos points à partir des propriétés de la couche - métadonnées
    - Laissez les autres paramètres inchangés
    - Cela devrait donner une grille avec chaque cellule égale à la valeur du point en plein milieu de la cellule, et les trous remplis selon l'algorithme IDW, et les valeurs dans la zone environnante
  1. Pour se débarrasser de ces cellules environnantes, vous devez masquer la grille.
    - Construisez la coque concave des points d'entrée comme polygone de délimitation minimum à l'aide de ce plugin

http://plugins.qgis.org/plugins/concavehull/

ou alors

Boîte à outils de traitement QGIS Géoalgorithmes QGIS Outils de géométrie vectorielle Coque concave

  1. Convertir la coque concave en grille

QGIS Processing Toolbox SAGA Grid - Gridding Shapes to grid

  1. et enfin masquer la grille interpolée avec la coque quadrillée

QGIS Processing Toolbox SAGA Grid - Tools Grid masking


Vous pouvez utiliser le plugin d'interpolation dans QGIS. Il prend les points en entrée, interpole et donne un raster en sortie.


Voir la vidéo: 16 GIS course Video QGIS DEM GRASS - from topo