REFAC: set vector data extensions, prefix, suffix, and SRS in environment variables

scripts/LoadEnvironmentVariables.sh

@@ -38,7 +38,7 @@ export ENVELOPES_DIR=/home/coupesrases/ENVELOPES
 export ENVELOPES_PREFIX=TILE_
 export ENVELOPES_SUFFIX=.shp
 
-# Tiles to download
+# Tiles to download and process
 export TILES=(T30UUU T30UVV T30UVU T30TVT T30UWA T30UWV T30UWU T30TWT T30TWS T30TWP T30TWN T30UXA T30UXV T30UXU T30TXT T30TXS T30TXR T30TXQ T30TXP T30TXN T30UYA T30UYV T30UYU T30TYT T30TYS T30TYR T30TYQ T30TYP T30TYN T31UCS T31UCR T31UCQ T31UCP T31TCN T31TCM T31TCL T31TCK T31TCJ T31TCH T31UDS T31UDR T31UDQ T31UDP T31TDN T31TDM T31TDL T31TDK T31TDJ T31TDH T31TDG T31UES T31UER T31UEQ T31UEP T31TEN T31TEM T31TEL T31TEK T31TEJ T31TEH T31UFR T31UFQ T31UFP T31TFN T31TFM T31TFL T31TFK T31TFJ T31TFH T31UGR T31UGQ T31UGP T31TGN T31TGM T31TGL T31TGK T31TGJ T31TGH T32ULV T32ULU T32TLT T32TLS T32TLR T32TLQ T32TLP T32UMV T32UMU T32TMN T32TMM T32TNN T32TNM T32TNL)
 
 export MAX_CLOUD_COVERAGE=100
@@ -48,6 +48,16 @@ export VEGETATION_PERIOD_START_DAY=01
 export VEGETATION_PERIOD_END_MONTH=08
 export VEGETATION_PERIOD_END_DAY=31
 
+# Suffixes for produced data (NDVI tiles, Output Vector Data)
+export NDVI_DIRNAME=NDVI
+export NDVI_SUFFIX=_NDVI
+export NDVI_COMPLETE_FILE=ndvi_complete.txt
+export OUTPUT_DIRNAME=CLEARCUTS
+export OUTPUT_VECTORDATA_PREFIX=vector_data_
+export OUTPUT_VECTORDATA_SUFFIX=_clip_epsg2154
+export OUTPUT_VECTORDATA_EXTENSION=.shp
+export OUTPUT_VECTORDATA_SRS=epsg:2154
+
 # OTB environment
 export OTB_INSTALL_DIR=/home/coupesrases/OTB/superbuild_install
 export LD_LIBRARY_PATH=${OTB_INSTALL_DIR}/lib/:${OTB_INSTALL_DIR}/lib/otb/applications
@@ -67,9 +77,3 @@ CheckDirectory $ENVELOPES_DIR
 export LOG_FILE=${DATA_DIR}/log.txt
 touch $LOG_FILE
 CheckFile $LOG_FILE
-
-# Suffixes
-export NDVI_DIRNAME=NDVI
-export NDVI_SUFFIX=_NDVI
-export NDVI_COMPLETE_FILE=ndvi_complete.txt
-export OUTPUT_DIRNAME=CLEARCUTS

scripts/ProcessTimeSeries.sh

@@ -81,8 +81,8 @@ do
   logging "Start otbapp. See $OUTPUT_LOG_FILE"
 
   # Names of produced files
-  OUTPUT_VECTORDATA=${CURRENT_TILE_OUTPUT_DIR}/vector_data_${TILE}.shp
-  OUTPUT_VECTORDATA_FINAL=${CURRENT_TILE_OUTPUT_DIR}/vector_data_cliped_${TILE}.shp
+  OUTPUT_VECTORDATA=${CURRENT_TILE_OUTPUT_DIR}/${OUTPUT_VECTORDATA_PREFIX}${TILE}${OUTPUT_VECTORDATA_EXTENSION}
+  OUTPUT_VECTORDATA_FINAL=${CURRENT_TILE_OUTPUT_DIR}/${OUTPUT_VECTORDATA_PREFIX}${TILE}${OUTPUT_VECTORDATA_SUFFIX}${OUTPUT_VECTORDATA_EXTENSION}
 
   # Compute time series
   otbcli_ClearCutsMultitemporalDetection -il $IMAGES_FILENAMES_LIST -dates $OUTPUT_DATES_FILE -masksdir $MASKS_DIR -outvec $OUTPUT_VECTORDATA -inputscale 0.0001 -forward.nmin 3 -sfilter connect >> $OUTPUT_LOG_FILE 2>&1
@@ -90,6 +90,6 @@ do
   # Clip vector data on the tile envelope
   CURRENT_TILE_ENVELOPE=${ENVELOPES_DIR}/${ENVELOPES_PREFIX}${TILE}${ENVELOPES_SUFFIX}
   CheckFile $CURRENT_TILE_ENVELOPE
-  ogr2ogr -clipsrc $CURRENT_TILE_ENVELOPE $OUTPUT_VECTORDATA_FINAL $OUTPUT_VECTORDATA
+  ogr2ogr -clipdst $CURRENT_TILE_ENVELOPE $OUTPUT_VECTORDATA_FINAL $OUTPUT_VECTORDATA -t_srs $OUTPUT_VECTORDATA_SRS
 
 done # next tile