Como interpolar espacialmente fácilmente más de 2 millones de puntos con GDAL - Tutorial

Interpolar grandes volúmenes de datos espaciales a través de la interfaz gráfica de QGIS suele ser un dolor de cabeza: el proceso tarda demasiado o, peor aún, el programa se congela consume todos los recursos del sistema. En este tutorial te enseñaremos a interpolar casi 2 millones de puntos en menos de dos minutos utilizando OSGeo4W Shell y comandos de GDAL con el método IDW. A través de la creación de un archivo de capa virtual (.vrt), aprenderás un flujo de trabajo altamente eficiente, seguro y ligero que te permitirá trabajar con millones de datos sin depender de los límites de la interfaz gráfica.

Tutorial

Codigo

Estructura del layer virtual:

<OGRVRTDataSource>
    <OGRVRTLayer name="points">
        <SrcDataSource relativeToVRT="1">CSV:2022-08-19-UTM-25cm.xyz</SrcDataSource>
        <SrcLayer>2022-08-19-UTM-25cm</SrcLayer>
        <OpenOptions>
            <OOI key="SEPARATOR">TAB</OOI>
            <OOI key="HEADERS">NO</OOI>
        </OpenOptions>
        <GeometryType>wkbPoint</GeometryType>
        <GeometryField encoding="PointFromColumns" x="field_1" y="field_2" z="field_3"/>
    </OGRVRTLayer>
</OGRVRTDataSource>

Comando de GDAL en OSGeo4W Shell

gdal_grid -a invdistnn:power=2.0:max_points=12:radius=20:nodata=-9999 -txe 697086.0 697829.0 -tye 4345594.0 4346001.0 -tr 1.0 1.0 -zfield "field_3" -of GTiff --config GDAL_NUM_THREADS ALL_CPUS points.vrt output_1m_grid.tif

Datos de ingreso

Puedes descarga los datos desde este enlace.