wkt to shape

tools/L2AproductListToRGB.sh

@@ -38,7 +38,12 @@ do
 	echo $BOX
         if [  ! -e $JP2DIR"rgbCroped/"$fileout1/$fileout2"_"$fileout1".jpg" ]; then
 
-	        gdal_translate -projwin $BOX -projwin_srs WGS84 -ot Byte -of JPEG $filename"TCI_10m.jp2" $JP2DIR"rgbCroped/"$fileout1/$fileout2"_"$fileout1"_10m.jpg"
+		$OUTIMG=$JP2DIR"rgbCroped/"$fileout1/$fileout2"_"$fileout1"_10m.jpg"
+	        gdal_translate -projwin $BOX -projwin_srs WGS84 -ot Byte -of JPEG $filename"TCI_10m.jp2" $OUTIMG
+		wktToShape.py $POLYGON myShape.shp
+		gdal_rasterize -b 1 -b 2 -b 3 -burn 10 -burn 200 -burn 10 -l wtk myShape.shp $OUTIMG
+
+		rm myShape.*
 ##	        convert rgb/$fileout1/$fileout2"_"$fileout1"_t.jpg" -resize 640x480\!  rgb/$fileout1/$fileout2"_"$fileout1".jpg"
 ##		rm rgb/$fileout1/$fileout2"_"$fileout1"_t.jpg"
 	fi