Skip to content

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

ndvi

parent a590935c
Show whitespace changes
Inline Side-by-side
Showing with 106 additions and 106 deletions
tools/productListToIndex.sh
View file @ f186f7df
...@@ -11,58 +11,58 @@ echo $POLYGON ...@@ -11,58 +11,58 @@ echo $POLYGON
BOX=$(python3 $HOME/GeoSentinel/geosentinel/polygonToBox.py "$POLYGON") 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
# $B1FN=$(ls $filepattern$B1".jp2" -t |head -1)
# $B2FN=$(ls $filepattern$B1".jp2" -t |head -1)
# echo "B1-"$B1FN cd $JP2DIR
# echo "B2-"$B2FN $B1FN=$(ls $filepattern$B1".jp2" -t |head -1)
$B2FN=$(ls $filepattern$B1".jp2" -t |head -1)
# if [ ! -d $OUTPUTDIR ]; then echo "B1-"$B1FN
# mkdir $OUTPUTDIR echo "B2-"$B2FN
# fi
# fileout1=$(echo $filepattern | cut -d"_" -f1) if [ ! -d $OUTPUTDIR ]; then
# fileout2=$(echo $filepattern | cut -d"_" -f2) mkdir $OUTPUTDIR
fi
# if [ ! -e $OUTPUTDIR"/"$fileout2"_"$fileout1".tiff" ]; then fileout1=$(echo $filepattern | cut -d"_" -f1)
# python3 $HOME/GeoSentinel/geosentinel/normalizedIndex.py $B1FN $B2FN $OUTPUTDIR"/"$fileout2"_"$fileout1".tiff" fileout2=$(echo $filepattern | cut -d"_" -f2)
# fi
# cd $OUTPUTDIR if [ ! -e $OUTPUTDIR"/"$fileout2"_"$fileout1".tiff" ]; then
# if [ ! -d $fileout1 ]; then python3 $HOME/GeoSentinel/geosentinel/normalizedIndex.py $B1FN $B2FN $OUTPUTDIR"/"$fileout2"_"$fileout1".tiff"
# mkdir $fileout1 fi
# fi
# cd ..
# if [ ! -e $OUTPUTDIR$fileout1/$fileout2"_"$fileout1".jpg" ]; then cd $OUTPUTDIR
# 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" if [ ! -d $fileout1 ]; then
# fi mkdir $fileout1
fi
cd ..
# if [ ! -e $OUTPUTDIR$fileout1/$fileout2"_"$fileout1".jpg" ]; then
# cd $PRODCUTSDIR 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
#done
#cd $JP2DIR"/video_ndvi" cd $PRODCUTSDIR
#ffmpeg -i %*.jpg -c:v libx264 -vf fps=10 -pix_fmt yuvj422p ndvi.mp4
done
cd $JP2DIR"/video_ndvi"
ffmpeg -i %*.jpg -c:v libx264 -vf fps=10 -pix_fmt yuvj422p ndvi.mp4
tools/productListToRGB.sh
View file @ f186f7df
...@@ -2,73 +2,73 @@ ...@@ -2,73 +2,73 @@
PRODCUTSDIR=$1 #Products directory PRODCUTSDIR=$1 #Products directory
JP2DIR=$2 #JP2 Directory JP2DIR=$2 #JP2 Directory
POLYGON=$3 #Crop Window POLYGON=$3 #Crop Window
echo $POLYGON
BOX=$(python3 $HOME/GeoSentinel/geosentinel/polygonToBox.py "$POLYGON") BOX=$(python3 $HOME/GeoSentinel/geosentinel/polygonToBox.py "$POLYGON")
echo $BOX
#cd $PRODCUTSDIR
#for f in $(find . -type f -name '*.zip')
#do
# filepattern=$(echo $f | cut -d"_" -f6)"_"$(echo $f | cut -d"_" -f3)"_B"
# filename=$JP2DIR$filepattern
# echo $filename
# if [ ! -e $filename"02.jp2" -o ! -e $filename"03.jp2" -o ! -e $filename"04.jp2" ]; then
# unzip -n -j $f *B02.jp2 *B03.jp2 *B04.jp2 -d $JP2DIR
# else
# echo PASS $f
# fi
#
# cd $JP2DIR
# red=$(ls $filepattern"04.jp2" -t |head -1)
# green=$(ls $filepattern"03.jp2" -t |head -1)
# blue=$(ls $filepattern"02.jp2" -t |head -1)
# echo "red-"$red
# echo "gre-"$green
# echo "blu-"$blue
# if [ ! -d "rgb" ]; then
# mkdir rgb
# fi
# fileout1=$(echo $filepattern | cut -d"_" -f1) cd $PRODCUTSDIR
# fileout2=$(echo $filepattern | cut -d"_" -f2) for f in $(find . -type f -name '*.zip')
do
# filepattern=$(echo $f | cut -d"_" -f6)"_"$(echo $f | cut -d"_" -f3)"_B"
# if [ ! -e rgb/$fileout2"_"$fileout1".tiff" ]; then filename=$JP2DIR$filepattern
# gdal_merge.py -separate -co PHOTOMETRIC=RGB -o rgb/$fileout2"_"$fileout1".tiff" $filepattern"04.jp2" $filepattern"03.jp2" $filepattern"02.jp2" echo $filename
# fi if [ ! -e $filename"02.jp2" -o ! -e $filename"03.jp2" -o ! -e $filename"04.jp2" ]; then
unzip -n -j $f *B02.jp2 *B03.jp2 *B04.jp2 -d $JP2DIR
else
echo PASS $f
fi
# cd rgb cd $JP2DIR
# if [ ! -d $fileout1 ]; then red=$(ls $filepattern"04.jp2" -t |head -1)
# mkdir $fileout1 green=$(ls $filepattern"03.jp2" -t |head -1)
# fi blue=$(ls $filepattern"02.jp2" -t |head -1)
# cd ..
echo "red-"$red
echo "gre-"$green
echo "blu-"$blue
if [ ! -d "rgb" ]; then
mkdir rgb
fi
# if [ ! -e rgb/$fileout1/$fileout2"_"$fileout1".jpg" ]; then fileout1=$(echo $filepattern | cut -d"_" -f1)
fileout2=$(echo $filepattern | cut -d"_" -f2)
# gdal_translate -projwin $BOX -projwin_srs WGS84 -ot Byte -scale 0 4096 0 255 -of JPEG rgb/$fileout2"_"$fileout1".tiff" rgb/$fileout1/$fileout2"_"$fileout1".jpg"
## convert rgb/$fileout1/$fileout2"_"$fileout1"_t.jpg" -resize 640x480\! rgb/$fileout1/$fileout2"_"$fileout1".jpg" if [ ! -e rgb/$fileout2"_"$fileout1".tiff" ]; then
## rm rgb/$fileout1/$fileout2"_"$fileout1"_t.jpg" gdal_merge.py -separate -co PHOTOMETRIC=RGB -o rgb/$fileout2"_"$fileout1".tiff" $filepattern"04.jp2" $filepattern"03.jp2" $filepattern"02.jp2"
fi
cd rgb
if [ ! -d $fileout1 ]; then
mkdir $fileout1
fi
cd ..
if [ ! -e rgb/$fileout1/$fileout2"_"$fileout1".jpg" ]; then
gdal_translate -projwin $BOX -projwin_srs WGS84 -ot Byte -scale 0 4096 0 255 -of JPEG rgb/$fileout2"_"$fileout1".tiff" rgb/$fileout1/$fileout2"_"$fileout1".jpg"
# convert rgb/$fileout1/$fileout2"_"$fileout1"_t.jpg" -resize 640x480\! rgb/$fileout1/$fileout2"_"$fileout1".jpg"
# rm rgb/$fileout1/$fileout2"_"$fileout1"_t.jpg"
fi
cd $PRODCUTSDIR
done
# if [ ! -d "video" ]; then
# mkdir video
# fi # fi
# cd $PRODCUTSDIR cd $JP2DIR"rgb/"
#done for d in */; do
## if [ ! -d "video" ]; then # Will print */ if no directories are available
## mkdir video cd $d
## fi # VIDEONAME = $(echo $d/ | cut -d"/" -f1)
# ffmpeg -i %*.jpg -c:v libx264 -vf fps=10 -pix_fmt yuvj422p $VIDEONAME.mp4
#cd $JP2DIR"rgb/" rm *.xml
#for d in */; do cd ..
# # Will print */ if no directories are available done
# cd $d
## VIDEONAME = $(echo $d/ | cut -d"/" -f1) #ffmpeg -i %*.jpg -c:v libx264 -vf fps=10 -pix_fmt yuvj422p rgb.mp4
## ffmpeg -i %*.jpg -c:v libx264 -vf fps=10 -pix_fmt yuvj422p $VIDEONAME.mp4
# rm *.xml
# cd ..
#done
##ffmpeg -i %*.jpg -c:v libx264 -vf fps=10 -pix_fmt yuvj422p rgb.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