Procesamiento de una carpeta de ficheros LAZ

mail@pastecode.io avatar
unknown
sh
a month ago
1.4 kB
1
Indexable
Never
#!/bin/bash

# Itera sobre cada archivo LAZ en el directorio ./laz/original que no tenga 'copc' en su nombre.
for laz in $(find ./laz/original -type f | grep -v copc); do
    echo $laz  # Imprime el nombre del archivo LAZ actual.
    step_value=1  # Define el valor de paso para las operaciones las2dem64.

    # Genera nombres de archivos derivados de $laz.
    tif_name=${laz//.laz/.tif}
    mds_file=${tif_name//original/mds}
    mdt_file=${tif_name//original/mdt}
    hillshade_file=${tif_name//original/hillshade}
    hillshade2_file=${tif_name//original/hillshade_2}
    rgb_file=${tif_name//original/rgb}

    # MDS - Modelo Digital de Superficie
    las2dem64 -step $step_value -i $laz -o $mds_file

    # MDT - Modelo Digital del Terreno
    las2dem64 -keep_class 2 -step $step_value -i $laz -o $mdt_file

    # RGB
    las2dem64 -step $step_value -i $laz -o $rgb_file -rgb

    # HILLSHADE
    las2dem64 -i $laz -hillshade -step $step_value -o $hillshade_file

    # HILLSHADE para la clase 2
    las2dem64 -keep_class 2 -i $laz -hillshade -step $step_value -o $hillshade2_file
done

# Construye archivos VRT (Virtual Raster Table) para cada tipo de producto raster.
gdalbuildvrt laz/mdt.vrt laz/mdt/*.tif
gdalbuildvrt laz/mds.vrt laz/mds/*.tif
gdalbuildvrt laz/hillshade.vrt laz/hillshade/*.tif
gdalbuildvrt laz/hillshade_2.vrt laz/hillshade_2/*.tif
gdalbuildvrt laz/rgb.vrt laz/rgb/*.tif
Leave a Comment