Procesamiento de una carpeta de ficheros LAZ
unknown
sh
2 years ago
1.4 kB
6
Indexable
#!/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
Editor is loading...
Leave a Comment