From 068c559f9f0779848d75e581a2769a7b31a24efd Mon Sep 17 00:00:00 2001 From: Stream Date: Wed, 13 May 2026 22:10:04 +0300 Subject: [PATCH] auto-sync: 2026-05-13 22:10:01 --- tasks/enduro-trails/scripts/gen_tri_full.sh | 85 +++++++++++++++++++++ 1 file changed, 85 insertions(+) create mode 100644 tasks/enduro-trails/scripts/gen_tri_full.sh diff --git a/tasks/enduro-trails/scripts/gen_tri_full.sh b/tasks/enduro-trails/scripts/gen_tri_full.sh new file mode 100644 index 0000000..fad3b0b --- /dev/null +++ b/tasks/enduro-trails/scripts/gen_tri_full.sh @@ -0,0 +1,85 @@ +#!/bin/bash +set -e + +echo '=== FULL TRI REGENERATION (32E-50E) ===' +echo "Start: $(date)" + +SRTM_DIR='/home/slin/enduro-trails/data/srtm' +TERRAIN_DIR='/home/slin/enduro-trails/data/terrain' +TRI_DIR="${TERRAIN_DIR}/tri" +TMP='/tmp/tri_gen' +mkdir -p $TMP + +echo '' +echo '=== Step 1: Merge all 101 HGT files ===' +ls $SRTM_DIR/*.hgt > $TMP/hgt_list.txt +echo "Files: $(wc -l < $TMP/hgt_list.txt)" + +gdalbuildvrt -input_file_list $TMP/hgt_list.txt $TMP/merged_full.vrt +echo "VRT built" + +echo '' +echo '=== Step 2: Compute TRI ===' +gdaldem TRI $TMP/merged_full.vrt $TMP/tri_raw_full.tif \ + -of GTiff -co COMPRESS=LZW -co BIGTIFF=YES +echo "TRI computed" + +echo '' +echo '=== Step 3: Color relief (thresholds x1.3) ===' +cat > $TMP/tri_color_v3.txt << 'RAMP' +nv 0 0 0 0 +0 0 0 0 0 +3 0 0 0 0 +4 255 255 100 60 +5 255 230 50 100 +6.5 255 200 0 140 +9 255 150 0 180 +13 255 100 0 200 +20 230 50 0 220 +26 200 0 0 235 +39 180 0 50 245 +65 150 0 100 255 +RAMP + +gdaldem color-relief \ + $TMP/tri_raw_full.tif \ + $TMP/tri_color_v3.txt \ + $TMP/tri_colored_full.tif \ + -alpha \ + -of GTiff \ + -co COMPRESS=LZW \ + -co BIGTIFF=YES +echo "Color relief done" + +echo '' +echo '=== Step 4: Generate tiles ===' +rm -rf ${TRI_DIR} +mkdir -p ${TRI_DIR} + +gdal2tiles.py \ + --zoom=5-12 \ + --processes=4 \ + --tilesize=256 \ + --webviewer=none \ + $TMP/tri_colored_full.tif \ + ${TRI_DIR}/ +echo "Tiles generated" + +echo '' +echo '=== Step 5: Stats ===' +for z in 5 6 7 8 9 10 11 12; do + total=$(find ${TRI_DIR}/$z/ -name '*.png' 2>/dev/null | wc -l) + nonempty=$(find ${TRI_DIR}/$z/ -name '*.png' -size +400c 2>/dev/null | wc -l) + echo " zoom $z: $total total, $nonempty non-empty" +done + +echo '' +du -sh ${TRI_DIR} + +# Also update merged_srtm.tif for future use +cp $TMP/tri_raw_full.tif $TMP/tri_raw.tif +gdal_translate -of GTiff -co COMPRESS=LZW $TMP/merged_full.vrt ${TERRAIN_DIR}/merged_srtm.tif + +echo '' +echo '=== DONE ===' +echo "End: $(date)"