Commit 35985b4b98ee371d424c712ff41d1d1a9c7c9130

Authored by Jordi Inglada
1 parent badd3f72
Exists in master

zonal stats for dem

Showing 1 changed file with 12 additions and 3 deletions   Show diff stats
jobs/zonal_stats.pbs
... ... @@ -8,13 +8,22 @@
8 8  
9 9 TERLAB_DIR=/work/OT/theia/oso/shapes/TERLAB/
10 10 DPT=31
  11 +TILE=T31TCJ
  12 +DEM_DIR=/work/OT/theia/oso/sensorsDatas/MNT_S2_N2/S2__TEST_AUX_REFDE2_${TILE}_0001.DBL.DIR
11 13 TERLAB_FILE=RPG_TERLAB_DEP${DPT}_2017
12 14 TERLAB_ZIP=$TERLAB_DIR/$TERLAB_FILE.7z
13 15 TERLAB_SHP=$TERLAB_DIR/SURFACES-2017-PARCELLES-GRAPHIQUES-CONSTATEES_0${DPT}_20180210.shp
14 16 MAASSP_DIR=/work/OT/theia/oso/jordi/maassp_data
15 17  
16   -7z x $TERLAB_ZIP -aoa -o$TERLAB_DIR
17   -#otbcli_ZonalStatistics -in ${MAASSP_DIR}/tiles/T31TCJ/Sentinel2_ST_MASK.tif -inbv -1000 -inzone.vector.in $TERLAB_SHP -inzone.vector.reproject 1 -out.vector.filename ${MAASSP_DIR}/zonal_stats_mask.shp -ram 30000
  18 +#7z x $TERLAB_ZIP -aoa -o$TERLAB_DIR
  19 +#otbcli_ZonalStatistics -in ${MAASSP_DIR}/tiles/T31TCJ/Sentinel2_ST_MASK.tif -inbv -10000 -inzone.vector.in $TERLAB_SHP -inzone.vector.reproject 1 -out.vector.filename ${MAASSP_DIR}/zonal_stats_mask.shp -ram 30000
18 20  
19   -otbcli_ZonalStatistics -in ${MAASSP_DIR}/tiles/T31TCJ/Sentinel2_ST_REFL.tif -inbv -1000 -inzone.vector.in $TERLAB_SHP -inzone.vector.reproject 1 -out.vector.filename ${MAASSP_DIR}/zonal_stats_refl.shp -ram 30000
  21 +#otbcli_ZonalStatistics -in ${MAASSP_DIR}/tiles/T31TCJ/Sentinel2_ST_REFL.tif -inbv -10000 -inzone.vector.in $TERLAB_SHP -inzone.vector.reproject 1 -out.vector.filename ${MAASSP_DIR}/zonal_stats_refl.shp -ram 30000
20 22  
  23 +#### DEM
  24 +ALTIMG=${DEM_DIR}/S2__TEST_AUX_REFDE2_${TILE}_0001_ALT_R1.TIF
  25 +ASPIMG=${DEM_DIR}/S2__TEST_AUX_REFDE2_${TILE}_0001_ASP_R1.TIF
  26 +SLPIMG=${DEM_DIR}/S2__TEST_AUX_REFDE2_${TILE}_0001_SLP_R1.TIF
  27 +DEMIMG=${MAASSP_DIR}/tiles/T31TCJ/DEM_${TILE}.tif
  28 +otbcli_ConcatenateImages -il ${ALTIMG} ${ASPIMG} ${SLPIMG} -out ${DEMIMG} -ram 30000
  29 +otbcli_ZonalStatistics -in ${DEMIMG} -inbv -10000 -inzone.vector.in $TERLAB_SHP -inzone.vector.reproject 1 -out.vector.filename ${MAASSP_DIR}/zonal_stats_dem.shp -ram 30000
... ...