Estoy usando el python api de liblas para leer/escribir archivos de las. Me gustaría obtener la información de referencia espacial de la archivo las liblas, pero han sido capaces de descubrirlo. ¿Alguien puede ofrecer algunos consejos sobre cómo ir sobre cómo obtener esta información? Gracias de antemano.
Respuestas
¿Demasiados anuncios?Para ser honesto, creo que la manera más fácil sería para shell a lasinfo con el --xml argumento y, a continuación, utilizar ElementTree o algo para arrancar lo que usted necesita desde el XML.
Se puede llegar a todo esto de los enlaces Python demasiado, pero es un poco de un desastre. En definitiva, abrir el archivo, ir a buscar el encabezado, y buscar el srs, y luego pedir la wkt de los srs. Todo esto supone que usted tiene libLAS + GeoTIFF + GDAL habilitado, y degrada mal si usted no tiene todas esas tres cosas vinculadas entre sí.
Además, si todo lo que quiero hacer es interactuar con los archivos LAS, laspy es una alternativa mucho mejor a libLAS en este momento. laspy no permite hacer nada con el sistema de coordenadas de cosas en LAS otras que para ir a buscar el GeoTIFF bytes del archivo, pero es una completa biblioteca en todos los demás. Es puro Python + numpy, fácil de implementar, y se comporta de forma natural, a diferencia de libLAS (descargo de responsabilidad, yo soy el libLAS autor y me ayudó a bootstrap laspy como un esfuerzo a la obsolescencia de los libLAS).
Espero que esto ayude,
Howard
No he utilizado liblas antes. Pero estaba leyendo la documentación y se encontró que puede necesitar el uso de las2las para el sistema de coordenadas y descripción de la alteración.
Los argumentos que usted necesita son:
--a_srs arg
o --t_srs arg
lasinfo también podría ser útil para conseguir extensiones.