rgb update

parent c83dcfa6
......@@ -3,55 +3,55 @@ PRODCUTSDIR=$1 #Products directory
JP2DIR=$2 #JP2 Directory
POLYGON=$3 #Crop Window
BOX=$(python3 "$HOME/GeoSentinel/geosentinel/polygonToBox.py" "$POLYGON")
cd $PRODCUTSDIR
for f in $(find . -type f -name '*.zip')
do
#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
# 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)
# 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
# echo "red-"$red
# echo "gre-"$green
# echo "blu-"$blue
if [ ! -d "rgb" ]; then
mkdir rgb
fi
# if [ ! -d "rgb" ]; then
# mkdir rgb
# fi
fileout1=$(echo $filepattern | cut -d"_" -f1)
fileout2=$(echo $filepattern | cut -d"_" -f2)
# fileout1=$(echo $filepattern | cut -d"_" -f1)
# fileout2=$(echo $filepattern | cut -d"_" -f2)
if [ ! -e rgb/$fileout2"_"$fileout1".jp2" ]; then
gdal_merge.py -separate -co PHOTOMETRIC=RGB -o rgb/$fileout2"_"$fileout1".jp2" $filepattern"04.jp2" $filepattern"03.jp2" $filepattern"02.jp2"
fi
# if [ ! -e rgb/$fileout2"_"$fileout1".jp2" ]; then
# gdal_merge.py -separate -co PHOTOMETRIC=RGB -o rgb/$fileout2"_"$fileout1".jp2" $filepattern"04.jp2" $filepattern"03.jp2" $filepattern"02.jp2"
# fi
if [ ! -d "video" ]; then
mkdir video
fi
# if [ ! -d "video" ]; then
# mkdir video
# fi
if [ ! -e video/$fileout2"_"$fileout1".jpg" ]; then
# if [ ! -e video/$fileout2"_"$fileout1".jpg" ]; then
gdal_translate -projwin $BOX -projwin_srs WGS84 -ot Byte -scale 0 4096 0 255 -of JPEG rgb/$fileout2"_"$fileout1".jp2" video/$fileout2"_"$fileout1".jpg"
convert video/$fileout2"_"$fileout1".jpg" -resize 640x480\! video/$fileout2"_"$fileout1".jpg"
fi
# gdal_translate -projwin $BOX -projwin_srs WGS84 -ot Byte -scale 0 4096 0 255 -of JPEG rgb/$fileout2"_"$fileout1".jp2" video/$fileout2"_"$fileout1".jpg"
# convert video/$fileout2"_"$fileout1".jpg" -resize 640x480\! video/$fileout2"_"$fileout1".jpg"
# fi
cd $PRODCUTSDIR
done
# cd $PRODCUTSDIR
#done
cd $JP2DIR"/rgb/video"
ffmpeg -i %*.jpg -c:v libx264 -vf fps=10 -pix_fmt yuvj422p rgb.mp4
#cd $JP2DIR"/rgb/video"
#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