Escribir un archivo LAZ con laspy

Question

Soy nuevo en el manejo de datos lidar. En realidad estoy usando laspy para leer y escribir archivos laz. Estoy usando la siguiente sintaxis para escribir el archivo LAZ. Se crea el archivo LAZ en la ruta de la carpeta especificada, pero tomó más tamaño en comparación con el LAZ de entrada (Por ejemplo: laz de entrada 50 mb, LAZ de salida 450 mb).

¿Cómo lo resuelvo?

        output_File="c:\test.LAZ"
        out_file = laspy.file.File(output_File, mode = "w", header=in_file.header)
        out_file.points = point_records
        out_file.close()

Answer 1

https://github.com/grantbrown/laspy dice

Laspy es una biblioteca pitónica para leer, modificar y escribir LAS de los archivos. El soporte para LAZ está limitado a la lectura de archivos LAS versión 1.0-1.3 de la versión 1.0-1.3.

laspy no tiene soporte para escribir archivos en formato LAZ.

Answer 2

Recomiendo encarecidamente PDAL para el manejo y procesamiento de estos archivos DTM. No estoy seguro de que laspy soporte directamente el formato LIDAR comprimido o requiera un plugin extra, pero con PDAL puedes convertir entre muchos formatos, aplicar filtros como elevación mínima/máxima, reproyectar, etc.

A continuación, un pequeño ejemplo.

import pdal
import json
pdal_pipeline = {
        "pipeline": [   
            *las_files,
            {
                "type":"filters.merge",
                "tag":"merged",              

            },
            {
                "type":"filters.reprojection",
                "in_srs":"EPSG:4326",
                "out_srs":"EPSG:32635"
            },
            {           
              "type":"writers.las",                     
              "filename": "merged_reprojected.laz"
            }
        ]
    }
pipeline = pdal.Pipeline(json.dumps(pdal_pipeline))
pipeline.execute()