Buenos días,
Estoy interesada en tratar datos lidar-PNOA para obtener un ráster con alturas. ¿Qué tipo de interpolación me recomiendan?
Muchas gracias.
Tratamiento de datos lidar
Re: Tratamiento de datos lidar
¡Hola!
Desde el menú “Fichero | Importación”, importa los datos LAS o LAZ (formato LAS comprimido) a formato PNT [módulo LASPNT]. Ello tiene la ventaja de disponer los datos y metadatos de manera más agradable (campos de tiempo de captación en dia-mes-año y en hora UTC separados para mejor filtro de consulta, tratamiento en tablas separadas (con lo que se consigue un tratamiento mas rápido) de datos secundarios como imágenes RGB que a menudo no son de la fecha de captación lidar, disposición ordenada y clara del número de pulsos retornados, el número concreto de retorno, las direcciones de captación, etc).
Si lo necesitas, puedes programar fácilmente una importación masiva de muchas hojas de distribución en un BAT o PS1 a través del botón disponible en la propia caja de diálogo contrastando la sintaxis con la que se expone en el botón Ayuda.
Abre el módulo InterPNT desde el menú “Herramientas | Interpolación y análisis del terreno | Interpolación a partir de puntos”. Observa que una vez seleccionado el fichero PNT de la importación puedes utilizar todos los puntos de la capa o aplicar una selección que te permitirá trabajar solamente con, por ejemplo, registros con intensidades más altas o más bajas, o acotadas entre determinadas Z para evitar cotas absurdas (retorno sobre pájaros, errores, etc), en determinadas fechas de vuelo, etc. Como método de interpolación te sugiero Estadísticas de vecindad. A través del botón Opciones podrás escoger la función estadística que más te convenga de un amplio abanico de posibilidades.
Entiendo que deseas obtener altitudes no “del terreno desnudo” (Modelo Digital de Elevaciones) sino de los objetos sobre el mismo (árboles, edificios, puentes, etc: Modelo Digital de Superficies, o MDS). En este caso lo mejor suele ser optar por la función “máximo” e indicar aplicarla en un radio (o en un cuadrado, puedes escoger el que mejor se adapte a tus necesidades) de tamaño típicamente algo mayor al tamaño de celda a obtener.
Observa que puedes obtener sendos rásters con el día, mes y año en cada celda, simultáneamente al MDS, con lo que al consultar cualquier punto del resultado también sabrás de qué fecha es el pulso utilizado para obtener la altura de la superficie.
De nuevo, si tienes que trabajar con muchas hojas, puede ser conveniente automatizar el tratamiento de todas en un fichero BAT o PS1 como he indicado antes. Una vez generado el MDS de cada hoja, puedes hacer un mosaico de todas ellas con el módulo Mosaic desde el menú “Herramientas | Organización espacial | Mosaico de capas”.
En el botón Avanzadas podrás escoger también una máscara, ráster o polígono, para evitar interpolar en zonas fuera de tu interés y acelerar también el cálculo.
Finalmente, si dispones también de un MDE (o lo has creado) puedes realizar la resta MDS-MDE para tener las alturas sobre el terreno. Para ello puedes utilizar el módulo CalcImg, accesible desde "Herramientas | Calculadora/Reclasificadora de rásters, vectores y tablas".
Espero que estas indicaciones te sean útiles.
Saludos.
Desde el menú “Fichero | Importación”, importa los datos LAS o LAZ (formato LAS comprimido) a formato PNT [módulo LASPNT]. Ello tiene la ventaja de disponer los datos y metadatos de manera más agradable (campos de tiempo de captación en dia-mes-año y en hora UTC separados para mejor filtro de consulta, tratamiento en tablas separadas (con lo que se consigue un tratamiento mas rápido) de datos secundarios como imágenes RGB que a menudo no son de la fecha de captación lidar, disposición ordenada y clara del número de pulsos retornados, el número concreto de retorno, las direcciones de captación, etc).
Si lo necesitas, puedes programar fácilmente una importación masiva de muchas hojas de distribución en un BAT o PS1 a través del botón disponible en la propia caja de diálogo contrastando la sintaxis con la que se expone en el botón Ayuda.
Abre el módulo InterPNT desde el menú “Herramientas | Interpolación y análisis del terreno | Interpolación a partir de puntos”. Observa que una vez seleccionado el fichero PNT de la importación puedes utilizar todos los puntos de la capa o aplicar una selección que te permitirá trabajar solamente con, por ejemplo, registros con intensidades más altas o más bajas, o acotadas entre determinadas Z para evitar cotas absurdas (retorno sobre pájaros, errores, etc), en determinadas fechas de vuelo, etc. Como método de interpolación te sugiero Estadísticas de vecindad. A través del botón Opciones podrás escoger la función estadística que más te convenga de un amplio abanico de posibilidades.
Entiendo que deseas obtener altitudes no “del terreno desnudo” (Modelo Digital de Elevaciones) sino de los objetos sobre el mismo (árboles, edificios, puentes, etc: Modelo Digital de Superficies, o MDS). En este caso lo mejor suele ser optar por la función “máximo” e indicar aplicarla en un radio (o en un cuadrado, puedes escoger el que mejor se adapte a tus necesidades) de tamaño típicamente algo mayor al tamaño de celda a obtener.
Observa que puedes obtener sendos rásters con el día, mes y año en cada celda, simultáneamente al MDS, con lo que al consultar cualquier punto del resultado también sabrás de qué fecha es el pulso utilizado para obtener la altura de la superficie.
De nuevo, si tienes que trabajar con muchas hojas, puede ser conveniente automatizar el tratamiento de todas en un fichero BAT o PS1 como he indicado antes. Una vez generado el MDS de cada hoja, puedes hacer un mosaico de todas ellas con el módulo Mosaic desde el menú “Herramientas | Organización espacial | Mosaico de capas”.
En el botón Avanzadas podrás escoger también una máscara, ráster o polígono, para evitar interpolar en zonas fuera de tu interés y acelerar también el cálculo.
Finalmente, si dispones también de un MDE (o lo has creado) puedes realizar la resta MDS-MDE para tener las alturas sobre el terreno. Para ello puedes utilizar el módulo CalcImg, accesible desde "Herramientas | Calculadora/Reclasificadora de rásters, vectores y tablas".
Espero que estas indicaciones te sean útiles.
Saludos.
Re: Tratamiento de datos lidar
Genial! Ha funcionado muy bien. ¡Gracias!