BOX

tools/productListToNDVI.sh

@@ -8,56 +8,55 @@ ext=.jpg
 BOX=$(python3 "$HOME/git/GeoSentinel/geosentinel/polygonToBox.py" "$3")
 echo $BOX
 
-#cd $BASEDIR
+cd $BASEDIR
 
-#for f in $(find . -type f -name '*.zip')
-#do 
+for f in $(find . -type f -name '*.zip')
+do 
 
-#	filepattern=$(echo $f | cut -d"_" -f6)"_"$(echo $f | cut -d"_" -f3)"_B"
-#	filename=$OUTDIR$filepattern
-#	echo $filename
-#	if [ ! -e $filename"04.jp2" ]; then
+	filepattern=$(echo $f | cut -d"_" -f6)"_"$(echo $f | cut -d"_" -f3)"_B"
+	filename=$OUTDIR$filepattern
+	echo $filename
+	if [ ! -e $filename"04.jp2" ]; then
 
-#                unzip -n -j $f *B04.jp2 -d $OUTDIR 
-#	fi
+                unzip -n -j $f *B04.jp2 -d $OUTDIR 
+	fi
 
-#	if [ ! -e $filename"08.jp2" ]; then
+	if [ ! -e $filename"08.jp2" ]; then
 
-#		unzip -n -j $f *B08.jp2 -d $OUTDIR 
-#	fi
+		unzip -n -j $f *B08.jp2 -d $OUTDIR 
+	fi
 
-#	
-#	cd $OUTDIR
-#	NIR=$(ls $filepattern"08.jp2" -t |head -1)
-#	VIR=$(ls $filepattern"04.jp2" -t |head -1)
 	
-#	echo "NIR-"$NIR
-#	echo "VIR-"$VIR
+	cd $OUTDIR
+	NIR=$(ls $filepattern"08.jp2" -t |head -1)
+	VIR=$(ls $filepattern"04.jp2" -t |head -1)
 
-#	if [ ! -d "ndvi" ]; then
-#		mkdir ndvi
-#	fi
+	echo "NIR-"$NIR
+	echo "VIR-"$VIR
 
-#	fileout1=$(echo $filepattern | cut -d"_" -f1)
-#        fileout2=$(echo $filepattern | cut -d"_" -f2)
+	if [ ! -d "ndvi" ]; then
+		mkdir ndvi
+	fi
 
-#	if [ ! -e ndvi/$fileout2"_"$fileout1".tiff" ]; then
-#		python3 $HOME/git/GeoSentinel/geosentinel/ndvi.py $NIR $VIR ndvi/$fileout2"_"$fileout1".tiff"
-#	fi
+	fileout1=$(echo $filepattern | cut -d"_" -f1)
+        fileout2=$(echo $filepattern | cut -d"_" -f2)
 
-#        if [ ! -d "video_ndvi" ]; then
-#                mkdir video_ndvi
-#        fi
-#	echo $3
-##	BOX=$(python3 "$HOME/git/GeoSentinel/geosentinel/polygonToBox.py" "$3")
-##	echo $BOX
-##	set GDAL_PAM_ENABLED=NO
-##	gdal_translate -projwin $BOX -projwin_srs WGS84 -ot Byte -scale 0 255 0 255 -of JPEG ndvi/$fileout2"_"$fileout1".tiff" video_ndvi/$fileout2"_"$fileout1".jpg"
-##	set GDAL_PAM_ENABLED=YES
-#	cd ..
+	if [ ! -e ndvi/$fileout2"_"$fileout1".tiff" ]; then
+		python3 $HOME/git/GeoSentinel/geosentinel/ndvi.py $NIR $VIR ndvi/$fileout2"_"$fileout1".tiff"
+	fi
 
-#done
-#cd $OUTDIR"/video_ndvi"
-#ffmpeg  -i %*.jpg -c:v libx264 -vf fps=10 -pix_fmt yuvj422p out.mp4
+        if [ ! -d "video_ndvi" ]; then
+                mkdir video_ndvi
+        fi
+	echo $3
+
+#	set GDAL_PAM_ENABLED=NO
+	gdal_translate -projwin $BOX -projwin_srs WGS84 -ot Byte -scale 0 255 0 255 -of JPEG ndvi/$fileout2"_"$fileout1".tiff" video_ndvi/$fileout2"_"$fileout1".jpg"
+#	set GDAL_PAM_ENABLED=YES
+	cd ..
+
+done
+cd $OUTDIR"/video_ndvi"
+ffmpeg  -i %*.jpg -c:v libx264 -vf fps=10 -pix_fmt yuvj422p out.mp4