Cómo Calcular Curvas de Nivel desde LiDAR en QGIS: Guía de Ingeniería 3D para CAD

1. Filtrado de la Nube de Puntos: El Algoritmo SMRF

La consolidación del núcleo de procesamiento de nubes de puntos en QGIS se fundamenta en la integración nativa de la biblioteca Point Data Abstraction Library (PDAL). A través del motor de ejecución binario pdal_wrench, QGIS proporciona un entorno de procesamiento multiproceso que realiza el análisis, filtrado y transformación de datos LiDAR directamente sobre archivos binarios .las y .laz. Esta arquitectura elimina por completo la dependencia histórica de herramientas propietarias o complementos externos con restricciones de licenciamiento, como LAStools.

Dentro de la Caja de Herramientas de Procesos de QGIS, el grupo de algoritmos de "Nubes de puntos" (Point clouds) se compone de tres familias fundamentales de herramientas nativas: Conversión de Nubes de Puntos (Point Cloud Conversion), Gestión de Datos de Nubes de Puntos (Point Cloud Data Management) y Extracción de Nubes de Puntos (Point Cloud Extraction). La equivalencia funcional directa con respecto a las herramientas clásicas se detalla en la siguiente matriz de correspondencia técnica:

Configuración de Parámetros en pdal:classifyground

Para aislar el terreno natural, el algoritmo nativo pdal:classifyground aplica el filtro SMRF, el cual requiere una parametrización geométrica precisa para evitar la omisión de relieve real o la inclusión de vegetación baja. Las variables críticas del algoritmo se definen a continuación:

• CELL_SIZE (Tamaño de celda de rejilla): Tipo de dato doble (m). Determina la resolución espacial de la cuadrícula de análisis morfológico inicial. Valores excesivamente bajos incrementan el ruido en el modelo, mientras que valores muy altos omiten variaciones topográficas locales.
• SCALAR (Escalar): Tipo de dato doble. Factor de multiplicación para pendientes pronunciadas. Actúa como un umbral para terrenos rugosos o montañosos, donde se requiere un valor más elevado para evitar que pendientes reales sean clasificadas erróneamente como objetos sobre el suelo.
• SLOPE (Pendiente): Tipo de dato doble (m/m). Umbral de pendiente máxima tolerable para clasificar puntos como suelo continuo. Debe configurarse con valores altos en relieves abruptos.
• THRESHOLD (Umbral): Tipo de dato doble (m). Altura máxima permitida para que un punto se considere parte de la superficie del suelo de referencia estimada.
• WINDOW_SIZE (Tamaño de ventana): Tipo de dato doble (m). Diámetro de búsqueda del filtro de apertura morfológica que discrimina objetos grandes (como naves industriales o coberturas forestales densas) del suelo continuo.

Aislamiento de la Clase 2 (Bare Earth)

Una vez clasificado el terreno, el algoritmo nativo para filtrar y exportar los puntos de suelo a un Modelo Digital del Terreno (MDT) en formato ráster es pdal:exportraster. Su sintaxis interna y la parametrización de sus variables clave se estructuran de la siguiente manera:

• INPUT: Capa de nube de puntos de entrada (archivo .las, .laz o índice virtual .vpc).
• ATTRIBUTE: Dimensión de la nube a extraer, definida como Z para procesar cotas de elevación.
• RESOLUTION: Dimensión lineal del píxel del ráster de salida (m).
• TILE_SIZE: Tamaño de la tesela para el procesamiento paralelo en la CPU (por defecto, 1000 metros).
• FILTER_EXPRESSION: Expresión lógica nativa de PDAL para filtrar los puntos de entrada. Para aislar de forma estricta los puntos clasificados como terreno, se define el siguiente parámetro sintáctico:
  Classification == 2

La carga y gestión de archivos masivos se realiza mediante la indexación espacial automática basada en el formato Cloud Optimized Point Cloud (COPC). Al cargar un archivo binario .las o .laz convencional, QGIS ejecuta un subproceso de lectura y reestructuración para generar una jerarquía de árbol octogonal (octree) que permite el acceso indexado bajo demanda. Si el archivo carece de índice, se genera una estructura Entwine Point Tile (ept_...) o un archivo .copc.laz adyacente que evita la lectura secuencial completa del archivo binario.

La estabilidad del sistema durante la visualización y análisis de datos masivos se regula mediante la gestión del "Presupuesto de puntos" (Point budget). El límite de estabilidad de renderizado en memoria RAM por defecto está parametrizado en un millón de puntos ($1.0 \times 10^6\text{ pts}$) para la visualización en el lienzo de trabajo. El motor de ejecución nativo evita el colapso del sistema administrando la densidad de visualización según el nivel de detalle (LoD) requerido por la escala del mapa, mientras que las tareas de procesamiento pesado mediante pdal_wrench se ejecutan por bloques espaciales (teselas) independientes con buffers de solape (collars) para mantener constante el consumo de recursos físicos.

2. Generación del MDT y Matriz de Rasterización Óptima

La determinación de la resolución espacial óptima del píxel (tamaño de celda en metros) para un Modelo Digital del Terreno (MDT) derivado de LiDAR se basa en directrices técnicas estandarizadas de agencias cartográficas líderes, como el Servicio Geológico de los Estados Unidos (USGS) en sus especificaciones para el programa 3DEP, y el Instituto Geográfico Nacional (IGN) de España en el marco del proyecto PNOA-LiDAR.

Criterio de Muestreo Nyquist-Shannon para LiDAR

El criterio de muestreo de Nyquist-Shannon adaptado a cartografía digital establece que para representar fielmente un terreno sin introducir distorsiones por submuestreo o redundancia estructural, el tamaño de celda teórica del ráster de salida ($G$, en metros) debe determinarse a partir del espaciamiento nominal de impulsos de la nube de puntos ($S$, en metros). El espaciamiento nominal de impulsos se relaciona directamente con la densidad nominal de puntos por metro cuadrado ($D$, en $\text{pt/m}^2$) mediante la expresión:

Para garantizar la estabilidad geométrica y la continuidad topográfica del MDT, evitando la generación de artefactos, el estándar de rasterización óptima se define mediante la aplicación de un factor multiplicador de escala sobre el espaciamiento nominal de impulsos:

La siguiente matriz detalla las equivalencias de resolución y densidad según las especificaciones técnicas del USGS y del IGN de España:

Ejecución de pdal:exportrastertin (Evitar el Efecto Ojo de Buey)

El proceso de clasificación de la nube de puntos genera áreas vacías de datos de suelo (clase 2) en zonas donde la vegetación densa absorbe o intercepta la señal del pulso láser, impidiendo retornos en la superficie terrestre. La elección del método de interpolación matemática para rellenar estas lagunas espaciales determina la calidad final del relieve:

La interpolación por Inverso de la Distancia al Cuadrado (IDW) estima el valor de elevación de cada celda de origen desconocido basándose en un promedio ponderado de los valores de los puntos circundantes, donde el peso disminuye exponencialmente con la distancia ($w_i = d_i^{-p}$). En coberturas boscosas densas donde los puntos de suelo válidos están sumamente aislados y dispersos, IDW genera el efecto "ojo de buey" (bull's eye). Este artefacto topográfico se caracteriza por la aparición de depresiones o sobreelevaciones artificiales concéntricas y circulares alrededor de los puntos aislados, alterando la morfología real del terreno. Además, si la densidad local disminuye por debajo del umbral de puntos mínimos definidos dentro del radio de búsqueda del algoritmo, IDW genera celdas con valores NoData (píxeles vacíos), interrumpiendo la continuidad matemática del MDT.

Por el contrario, la interpolación basada en Triangulación (TIN) construye en primera instancia una superficie continua tridimensional compuesta por caras triangulares planas empleando la triangulación de Delaunay sobre todos los puntos clasificados como suelo. Al rasterizar esta malla continua para generar el MDT, la elevación de cada píxel intermedio se deduce mediante una interpolación lineal simple basada en la posición planar exacta del centro del píxel dentro de su faceta triangular correspondiente. Este procedimiento asegura que los vacíos de información geométrica generados por el filtrado de vegetación densa queden cubiertos por facetas inclinadas continuas, eliminando por completo la formación de "píxeles vacíos" (NoData) y el efecto "ojo de buey". El relieve resultante mantiene la coherencia geométrica de las pendientes, lo que resulta indispensable para un trazado limpio y homogéneo de las curvas de nivel vectoriales.

3. Extracción Vectorial de Contornos con gdal_contour

La extracción de curvas de nivel propiamente dicha a partir de nuestro Modelo Digital del Terreno (MDT) se efectúa de manera óptima utilizando el backend GDAL integrado. Este proceso lee los valores de elevación del archivo ráster de entrada y genera líneas que conectan píxeles con la misma cota de forma continua.

Configuración del Intervalo y Atributo de Elevación

Para configurar de forma correcta la extracción de curvas, se debe acceder al menú superior de QGIS y ejecutar el algoritmo de GDAL a través de la ruta: Ráster ➜ Extracción ➜ Curvas de nivel.... En la interfaz del algoritmo gdal_contour, se deben parametrizar las siguientes variables esenciales:

• Capa de entrada: El MDT en formato ráster .tif generado en el paso anterior mediante interpolación TIN (pdal:exportrastertin).
• Intervalo entre curvas de nivel: Define la equidistancia en metros (p. ej., 1.0 m para curvas finas y 5.0 m para curvas maestras).
• Nombre del atributo: Define de forma estricta ELEV o ELEVATION. Este campo de tipo flotante en la tabla de atributos de salida almacenará la cota numérica extraída.
• Curvas de nivel: Ruta y nombre del archivo vectorial intermedio (p. ej., formato GeoPackage o Shapefile).

4. Resolución del Error del Eje Z: Inyección de Geometría 3D

Al exportar un archivo de curvas de nivel desde un entorno SIG a un archivo DXF/DWG, es común que AutoCAD interprete las polilíneas como elementos bidimensionales planos en cota cero ($Z=0$), a pesar de que la tabla de atributos asociada contenga el valor correcto de altitud.

La causa raíz exacta de este comportamiento reside en el tipo de entidad vectorial empleada en la conversión. Por defecto, QGIS exporta las geometrías lineales planas (tipo de objeto LineString) mapeándolas en el formato DXF como entidades de tipo LWPOLYLINE (Lightweight Polyline). La especificación técnica de la entidad LWPOLYLINE está optimizada para el dibujo bidimensional eficiente y no admite valores de coordenada $Z$ independientes en cada vértice. En su lugar, almacena un único valor de elevación uniforme para toda la entidad en el código de grupo DXF 38. Si la geometría de la polilínea en QGIS carece de componentes tridimensionales reales en sus vértices (es decir, si la geometría es puramente 2D), el motor de exportación de vectores no genera el registro para el código de grupo 38 en el archivo DXF, provocando que AutoCAD sitúe la polilínea en su plano por defecto de elevación cero ($Z=0$).

Por el contrario, para que AutoCAD reconozca curvas de nivel con elevaciones variables en cada vértice utilizando entidades tridimensionales de tipo 3DPOLYLINE, los datos de elevación deben estar inyectados en los códigos de grupo de vértice DXF 10 (coordenada $X$), 20 (coordenada $Y$) y 30 (coordenada $Z$), lo cual requiere de forma estricta que la geometría origen en QGIS sea de tipo tridimensional LineStringZ.

Paso 1: Extracción de curvas de nivel desde el MDT

Se accede al menú superior de QGIS y se ejecuta el algoritmo de extracción de GDAL a través de la ruta: Ráster ➜ Extracción ➜ Curvas de nivel.... En la interfaz del algoritmo gdal_contour, se configuran las siguientes variables de trabajo:

• Capa de entrada: El Modelo Digital del Terreno en formato ráster .tif derivado de la triangulación de la nube de puntos (pdal:exportrastertin).
• Intervalo entre curvas de nivel: Se define la equidistancia deseada en metros (p. ej., 1.0 o 5.0 metros).
• Nombre del atributo: Se define de forma fija en ELEV o ELEVATION. Este campo de la tabla de atributos almacenará el valor numérico flotante de la cota extraída de cada píxel de origen.
• Curvas de nivel: Ruta y nombre del archivo vectorial intermedio (p. ej., formato GeoPackage o Shapefile).

Paso 2: Inyección de la coordenada Z en la geometría del vector

Para transformar las curvas de nivel vectoriales planas de 2D a polilíneas 3D reales, abre la Caja de Herramientas de Procesos de QGIS y localiza el algoritmo: Geometría vectorial ➜ Establecer valor Z (ID técnico: native:setzvalue). Ajusta las variables bajo las siguientes directrices técnicas:

• Capa de entrada: La capa de curvas de nivel vectoriales generada en el paso anterior.
• Valor Z: Presiona el botón del menú desplegable a la derecha del parámetro para activar la Suplantación definida por datos (Data defined override) ➜ Tipo de campo (Field type) y selecciona el atributo numérico "ELEV". Este proceso de forma dinámica vincula la coordenada geométrica vertical con la cota escalar del campo tabular.
• Z agregado: Ruta de salida para la nueva capa vectorial tridimensional. Al inspeccionar las propiedades de este archivo de salida, se confirmará que la definición geométrica ha cambiado de tipo de LineString a LineStringZ.

Paso 3: Conversión de proyección a coordenadas métricas

Bajo ninguna circunstancia se deben exportar curvas de nivel en un Sistema de Coordenadas Geográficas (grados decimales como el elipsoide WGS 84 / EPSG:4326). AutoCAD interpreta las unidades angulares de longitud y latitud del archivo de intercambio como unidades métricas planas, lo que compactará todo tu proyecto a una escala microscópica concentrada en las cercanías inmediatas del origen físico $(0,0)$.

Para solventar esto, se debe aplicar el algoritmo: Vectorial ➜ Herramientas de gestión de datos ➜ Reproyectar capa. Elige un Sistema de Referencia de Coordenadas Proyectado (PCS) métrico de tipo UTM adecuado para la cuadrícula geográfica del proyecto (p. ej., EPSG:25830 - ETRS89 / UTM zone 30N para España peninsular).

Puedes automatizar este proceso de conversión geométrica mediante la consola de comandos de Python de QGIS ejecutando la siguiente lógica de reproyección sistemática:

5. Exportación Limpia a Formato DXF para AutoCAD

Para consolidar el archivo final de intercambio sin pérdida de datos en la coordenada Z, utiliza la herramienta nativa del sistema accediendo a: Proyecto ➜ Importar/Exportar ➜ Exportar proyecto a DXF.... En el diálogo de configuración, se deben definir exactamente los siguientes parámetros:

• Guardar como: Ruta y nombre del archivo de dibujo con extensión .dxf.
• SRC: Selecciona de manera estricta el SRC proyectado métrico de trabajo (p. ej., EPSG:25830).
• Modo de simbología: Configurado como Sin simbología. Esto exportará polilíneas limpias y ligeras sin la adición de bloques complejos de renderizado que saturan la memoria del software CAD.
• Exportar capas vectoriales: Marca únicamente la casilla correspondiente a la capa de curvas de nivel reproyectada en formato LineStringZ.
• Incluir dimensión Z: Marca de forma obligatoria la casilla (✓). Esta instrucción le ordena al backend del driver OGR que fuerce la escritura tridimensional de cada vértice del contorno en el archivo DXF final.

Al abrir el archivo resultante en AutoCAD, Civil 3D o MicroStation, cada polilínea conservará su posición 3D real, permitiendo su uso inmediato en el modelado de superficies del terreno y cálculos de movimiento de tierras.

7. Guía de Solución de Errores Críticos (FAQ de Foros OSGeo)

Durante el procesamiento masivo de datos geoespaciales se reportan de forma repetida incidencias derivadas de los límites estructurales de formatos y hilos de procesamiento. Mapeamos las tres incidencias principales de los foros de soporte técnico en 2026: