¿Existe alguna herramienta para convertir todos los valores de los atributos Z (elipsoide h) en altura ortométrica (H)? La idea es, teniendo una nube de puntos (.las) y la ondulación del geoide EGM08 (.asc) obtener una nube de puntos (.las) con valores ortométricos Z.
Para obtener alturas ortométricas es necesario calcular N (altura del geoide) y para ello se utiliza una interpolación bilineal.
¿Cómo puede hacerse esta interpolación?
https://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.interpolate.interp2d.html
0 votos
Terrascan puede hacerlo, aunque no estoy seguro de cómo pulsar todos los botones para ello. ¿De dónde ha sacado el LAS? ¿Puedes pedirle a tu proveedor que te dé los archivos LAS en altura ortométrica? ¿No es EGM08 (Earth Gravitational Model) y no EMG08? Puede parecer pedante, pero podría suponer una gran diferencia a la hora de buscar el término correcto.
0 votos
Hay muchas formas de atacar este problema. Asumiendo que quieres permanecer en Python, si tienes tus propios archivos de cuadrícula, podrías leer los puntos con LASpy, leer el archivo de cuadrícula usando GDAL, y escribir un nuevo archivo con LASpy representando tus alturas ortométricas. Si quieres simplemente ejecutar un comando, PDAL puede convertir datos LAS de elipsoide (EPSG:4326) a ortométricos (EPSG:5773) usando las capacidades de transformación de datum vertical de Proj.4. La pregunta necesita más detalles sobre lo que tienes y lo que intentas conseguir.
0 votos
¿Has mirado en VDatum de NOAA? vdatum.noaa.gov
0 votos
Tengo un .las de LiDAR con más de 4M puntos. Compruebo la altura de estos puntos comparándolos con un modelo de superficie, y descubro que los valores Z son la altura del elipsoide. Necesito transformar los valores Z en ortométricos. @MichaelStimson estas en el camino, hice un script en python que lee un .asc (EGM08) y construye una funcion ((scipy.interpolate.interp2d) para obtener una altura geoide dando latitud y longitud. He comprobado los resultados con algunos puntos y está bien. Para evitar hacer un módulo para esta tarea, pedí una herramienta que pueda hacer esto. Pero creo que estoy muy cerca de la solución con LASpy.
0 votos
Me parece justo. Aunque Terrascan puede hacerlo Escribí mi propio programa de consola C++ para hacer lo que estás haciendo en python; 4 millones de puntos (¿o era pulsos?) es un muy pequeño número de registros cuando se trata de LAS, por lo que debe ser capaz de procesar que en python, pero si esto es algo que estás haciendo con regularidad o para áreas más grandes que usted debe considerar algún software de procesamiento LiDAR como Terrascan, sobre todo para la fijación cuando se encuentran anomalías en el LiDAR que afecta a su salida.
0 votos
@HowardButler Las celdas de la GRID son muy grandes, en consecuencia la diferencia de valores entre una celda y la adyacente es significativa, si simplemente lees el valor en la GRID acabarás con "escalones" en el LiDAR... tu debe interpolar Sin embargo, tu segundo método parece interesante, siempre que obtengas los datos verticales correctos, pero aún así dudo que se utilice el modelo gravitacional correcto, ya que pueden existir varios (he visto al menos 3 para la costa este de Australia, 2 para Nueva Zelanda, así como una docena de datos verticales MSL), así que ¿qué modelo utiliza PDAL?
0 votos
PDAL utiliza el modelo definido en vertcs.csv de Proj4 para el código EPSG dado. Proj.4 interpola, y no sólo de forma escalonada.
0 votos
¿Sólo para hacer eso? ¿Es Terrascan licencia gratuita sólo para hacer eso?. Sí, tengo que hacer esto como un lote, pero esto no es un problema. ¿no?
0 votos
@MichaelStimson, he utilizado un modelo generado por un instituto gubernamental y es muy preciso (está basado en EGM08), y los datos de prueba son de una estación geodésica y el error es inferior al centímetro.
0 votos
Terrascan no es gratuito (que yo sepa), se ejecuta en Microstation (que tampoco es gratuito). Si ya tienes Esri y 3d analyst (o LP360) puedes hacer alguna edición menor en el conjunto de datos LAS pero es muy lento LASTools tiene un editor... No quiero promocionar un software sobre los demás. Es posible que ya dispongas de las herramientas que necesitas. Es probable que el modelo aprobado por la agencia geodésica de su gobierno sea el mejor para su situación, pero tenga cuidado si intenta comparar con otro conjunto de datos de otro país, ya que las diferencias pueden ser significativas.