Skip to content

GeoSentinel
GeoSentinel
Commit 26acb51d authored by Mario Chirinos Colunga's avatar Mario Chirinos Colunga 💬
Browse files
Options

ndvi

parent e163a1bc
Hide whitespace changes
Inline Side-by-side
Showing with 42 additions and 42 deletions
tools/productListToIndex.sh
View file @ 26acb51d
...@@ -2,67 +2,67 @@ ...@@ -2,67 +2,67 @@
# Get Normalized Index (B1-B2)/(B1+B2) for all products in PRODCUTSDIR # Get Normalized Index (B1-B2)/(B1+B2) for all products in PRODCUTSDIR
PRODCUTSDIR=$1 #Products directory PRODCUTSDIR=$1 #Products directory
JP2DIR=$2 #JP2 Directory JP2DIR=$2 #JP2 Directory
WINDOW=$3 #Crop Window POLYGON=$3 #WKT Polygon
OUTPUTDIR=$4 #Process output OUTPUTDIR=$4 #Process output
B1=$5 B1=$5
B2=$6 B2=$6
count=1 count=1
ext=.jpg echo $POLYGON
BOX=$(python3 "$HOME/GeoSentinel/geosentinel/polygonToBox.py" "$WINDOW") BOX=$(python3 "$HOME/GeoSentinel/geosentinel/polygonToBox.py" "$POLYGON")
echo $BOX echo $BOX
cd $PRODCUTSDIR #cd $PRODCUTSDIR
for f in $(find . -type f -name '*.zip') #for f in $(find . -type f -name '*.zip')
do #do
filepattern=$(echo $f | cut -d"_" -f6)"_"$(echo $f | cut -d"_" -f3)"_B" # filepattern=$(echo $f | cut -d"_" -f6)"_"$(echo $f | cut -d"_" -f3)"_B"
filename=$JP2DIR$filepattern # filename=$JP2DIR$filepattern
echo $filename # echo $filename
if [ ! -e $filename$B1".jp2" ]; then # if [ ! -e $filename$B1".jp2" ]; then
unzip" -n -j $f *"$B2".jp2 -d " $JP2DIR # unzip" -n -j $f *"$B2".jp2 -d " $JP2DIR
fi # fi
if [ ! -e $filename$B2".jp2" ]; then # if [ ! -e $filename$B2".jp2" ]; then
unzip" -n -j $f *"$B1".jp2 -d " $JP2DIR # unzip" -n -j $f *"$B1".jp2 -d " $JP2DIR
fi # fi
#
cd $JP2DIR # cd $JP2DIR
$B1FN=$(ls $filepattern$B1".jp2" -t |head -1) # $B1FN=$(ls $filepattern$B1".jp2" -t |head -1)
$B2FN=$(ls $filepattern$B1".jp2" -t |head -1) # $B2FN=$(ls $filepattern$B1".jp2" -t |head -1)
echo "B1-"$B1FN # echo "B1-"$B1FN
echo "B2-"$B2FN # echo "B2-"$B2FN
if [ ! -d $OUTPUTDIR ]; then # if [ ! -d $OUTPUTDIR ]; then
mkdir $OUTPUTDIR # mkdir $OUTPUTDIR
fi # fi
fileout1=$(echo $filepattern | cut -d"_" -f1) # fileout1=$(echo $filepattern | cut -d"_" -f1)
fileout2=$(echo $filepattern | cut -d"_" -f2) # fileout2=$(echo $filepattern | cut -d"_" -f2)
if [ ! -e $OUTPUTDIR"/"$fileout2"_"$fileout1".tiff" ]; then # if [ ! -e $OUTPUTDIR"/"$fileout2"_"$fileout1".tiff" ]; then
python3 $HOME/GeoSentinel/geosentinel/normalizedIndex.py $B1FN $B2FN $OUTPUTDIR"/"$fileout2"_"$fileout1".tiff" # python3 $HOME/GeoSentinel/geosentinel/normalizedIndex.py $B1FN $B2FN $OUTPUTDIR"/"$fileout2"_"$fileout1".tiff"
fi # fi
cd $OUTPUTDIR # cd $OUTPUTDIR
if [ ! -d $fileout1 ]; then # if [ ! -d $fileout1 ]; then
mkdir $fileout1 # mkdir $fileout1
fi # fi
cd .. # cd ..
if [ ! -e $OUTPUTDIR$fileout1/$fileout2"_"$fileout1".jpg" ]; then # if [ ! -e $OUTPUTDIR$fileout1/$fileout2"_"$fileout1".jpg" ]; then
gdal_translate -projwin $BOX -projwin_srs WGS84 -ot Byte -scale 0 255 0 255 -of JPEG $OUTPUTDIR$$fileout2"_"$fileout1".tiff" $OUTPUTDIR$fileout1/$fileout2"_"$fileout1".jpg" # gdal_translate -projwin $BOX -projwin_srs WGS84 -ot Byte -scale 0 255 0 255 -of JPEG $OUTPUTDIR$$fileout2"_"$fileout1".tiff" $OUTPUTDIR$fileout1/$fileout2"_"$fileout1".jpg"
fi # fi
#
cd $PRODCUTSDIR # cd $PRODCUTSDIR
done #done
cd $JP2DIR"/video_ndvi" #cd $JP2DIR"/video_ndvi"
ffmpeg -i %*.jpg -c:v libx264 -vf fps=10 -pix_fmt yuvj422p ndvi.mp4 #ffmpeg -i %*.jpg -c:v libx264 -vf fps=10 -pix_fmt yuvj422p ndvi.mp4
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