Commit d048126d authored by Simon Gascoin's avatar Simon Gascoin
Browse files

clean script to compute latest clear observation in hal

parent a1a2f9e4
......@@ -24,7 +24,7 @@
# example: ./computeLatestClearObservation.sh 31TCH /work/OT/siaa/Theia/Neige/Theia/ /tmp/ 20
# output: geotiff file with colormap cyan=snow, grey=no-snow, white=cloud, black=no-data
# tile ID (ex. 31TCH)
# Tile ID (ex. 31TCH)
tileId=$1
# Input directory where to look for the products
inputDir=$2
......@@ -32,39 +32,35 @@ inputDir=$2
outputDir=$3
# Number of days to go back
nDays=$4
# Input directory containing the snow products
#inputDir="/datalake/L2B-SNOW/${tileId}" # faster than /datalake/L2B-SNOW
#inputDir="/work/OT/siaa/Theia/Neige/Theia/"
# Today in yyyymmdd format
today=$(date +"%Y%m%d")
# name of the output composite
# Name of the output composite
outputFn="${outputDir}/LCO_${tileId}_${today}.tif"
# temporary file
# Temporary file
tmpFn="${outputDir}/tmp_${tileId}_${today}.tif"
# get date list of the n latest days
# Get date list of the n latest days
dateList=$(for i in $(seq $nDays); do date -d "$today -$i day" +"%Y%m%d"; done)
# search matching products in the tile folder and store in an array
# Search matching products in the tile folder and store in an array
productList=($(for d in ${dateList}; do find ${inputDir} -name "SENTINEL2[A,B]_${d}*${tileId}*SNW_R2.tif"; done))
echo $productList
# test if the product list is empty
# Test if the product list is empty
if [ -z "$productList" ]
then
# no products available: exit
# No products available: exit
echo "could not find products matching these dates: "$dateList
exit 1
else
# load OTB
# Load OTB
module load otb/7.0
# number of products to be composited
# Number of products to be composited
nProduct=${#productList[@]}
# init band math expression
# Init band math expression
cmd=im1b1
# build band math expression
# Build band math expression
for i in $(seq $((nProduct-1))); do cmd=$cmd" <= 100 ? im${i}b1 : im$((i+1))b1" ; done
# execute band math
# Execute band math
otbcli_BandMath -progress true -il ${productList[@]} -out "$tmpFn" uint8 -exp "$cmd"
# apply color table to composite
# Apply color table to composite
otbcli_ColorMapping -progress true -in "$tmpFn" -out "$outputFn""?&gdal:co:COMPRESS=DEFLATE" uint8 -method.custom.lut "/work/OT/siaa/Theia/hpc_scripts/LIS_SEB_style_OTB.txt"
# delete temporary file
# Delete temporary file
rm "$tmpFn"
fi
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment