Más

¿Cómo llenar huecos / valores faltantes en datos XYZ / un DEM con QGIS?

¿Cómo llenar huecos / valores faltantes en datos XYZ / un DEM con QGIS?


Tengo datos equidistantes de puntos XYZ (texto) que utilicé para crear una cuadrícula de puntos. Sin embargo, no es perfecto, hay algunos agujeros. Luego creé un ráster DSM a partir de él, que obviamente también tiene los agujeros. Los agujeros van desde un solo punto hasta áreas más grandes (faltan ~ 50 o más puntos).

¿Cómo reparo esos agujeros? Traté de interpolar el ráster, pero eso reduce la calidad y no corrige los "agujeros" cero, solo están ligeramente elevados después. Idealmente, me gustaría completar esos valores cero con los valores Z de sus cuadrados vecinos.

Otro problema es que el conjunto de datos XYZ no es rectangular, sino como en la tercera imagen adjunta (ruta roja). Si utilizo un método de relleno, ¿cómo evito llenar todo el área circundante?


Sugiero romper el proceso en los siguientes pasos.

  1. Si los puntos están igualmente espaciados (lo que obviamente es el caso), determine el espaciado en las direcciones x e y
  2. Interpolar un ráster inicial. Sugeriría usar un algoritmo simple ponderado por distancia inversa.

Caja de herramientas de procesamiento QGIS SAGA Grid - Gridding Distancia inversa ponderada

    - Dar el espaciado como argumento para el parámetro Cellsize, y
    - especifique la Extensión de salida como la extensión de sus puntos menos la mitad del espaciado para xmin / ymin y más la mitad del espaciado para xmax / ymax. Obtienes la extensión de tus puntos de las propiedades de la capa - metadatos
    - Deje los otros parámetros sin cambios
    - Esto debería resultar en una cuadrícula con cada celda igual al valor del punto justo en el medio de la celda, y los agujeros llenados de acuerdo con el algoritmo IDW y los valores en el área circundante.
  1. Para deshacerse de estas celdas circundantes, debe enmascarar la cuadrícula.
    - Construya el casco cóncavo de los puntos de entrada como el polígono delimitador mínimo usando este complemento

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

o

Caja de herramientas de procesamiento QGIS Geoalgoritmos QGIS Herramientas de geometría vectorial Casco cóncavo

  1. Convierta el casco cóncavo en una cuadrícula

Caja de herramientas de procesamiento QGIS SAGA Grid - Gridding Shapes to grid

  1. y finalmente enmascarar la cuadrícula interpolada con el casco cuadriculado

Caja de herramientas de procesamiento QGIS SAGA Grid - Tools Grid masking


Puede usar el complemento de interpolación en QGIS. Toma los puntos como entrada, los interpola y da ráster como salida.


Ver el vídeo: Cortar Ráster en QGis - Extracción por Máscara con polígono