processing

tools/L2ASCL_AreaProcessing.sh

@@ -15,8 +15,8 @@ NC='\033[0m' # No Color
 echo $USERDIR
 cd $USERDIR
 
-#Link L2A products 
-echo "Linking Products"
+#1.- Link L2A products 
+echo "Linking Products..."
 if [ ! -d "L2A" ]; then
 	mkdir L2A
 fi
@@ -25,13 +25,34 @@ cd L2A
 cd ..
 echo yes | findProducts.py $USERDIR"findProducts.json"
 
-#Extract Images
-echo "Extracting JP2 Images"
+#2.- Extract Images
+echo "Extracting JP2 Images..."
 if [ ! -d "jp2" ]; then
 	mkdir jp2
 fi
 L2AProductListExtractData.sh $USERDIR"L2A/" $USERDIR"jp2/" 4
 
-#Merge Images
-echo "Merging Images"
+#3.- Merge Images
+echo "Merging Images..."
 mergeImagesByDirectory.sh $USERDIR"jp2/" "$WKT" 4
+
+#4.- Delete images with few data
+cd jp2/merged_out/
+if [ ! -d "SCL" ]; then
+	mkdir SCL
+fi
+cp *SCL_20m.tif SCL
+cd SCL
+ls *.tif | parallel --jobs 4 imageMissingData.py {} 50
+
+##for tci in $($ls *.tif)
+##do
+##	fileprefix=$(echo tci )
+##done
+cd $USERDIR
+##Extract SCL information
+#if [ ! -d "json" ]; then
+#	mkdir json
+#fi
+#mv *SCL_20m.tif json
+#cd json