Commit dc3c782f7a9071f7fce54b5bad527531d4727bb2

Authored by Jordi Inglada
1 parent 3ae3dbd5
Exists in master

Notes sur la simplification de la generation zonal stats

Showing 1 changed file with 25 additions and 1 deletions   Show diff stats
notebook/zonal_stats.org
... ... @@ -13,13 +13,37 @@ import pandas as pd
13 13 #+RESULTS:
14 14  
15 15  
  16 +* Simplification du fichier TERLAB
  17 +- Eliminer les géométries qui ne contiennent pas d'information de rendement.
  18 +
  19 +* Génération des fichiers de stats zonales
  20 +Le code est ici [[file:/work/OT/theia/oso/jordi/src/maassp/jobs/zonal_stats.pbs]]
  21 +
  22 +Comme l'image pour les stats zonales contient 720 bandes, le fichier vecteur contient > 720 * 4 attributs, ce qui dépasse le nb de colonnes dans un fichier sqlite. Le shapefile n'a pas ce problème, mais est limité à 4GB, ce qui est trop petit.
  23 +
  24 +J'ai ensuite séparé la série en 2 demi-séries pour pouvoir générer 2 shapefiles, mais on a ensuite le pb de faire le lien entre les polygones des 2 shapefiles générés, parce qu'on perd l'identifiant.
  25 +
  26 +On pourrait commencer par ajouter un identifiant unique à chaque polygone du fichier TERLAB de façon à pouvoir faire le lien entre les différents fichiers de stats locales (reflectances, DEM, masques) générés à partir du même fichier TERLAB.
  27 +
  28 +#+begin_src bash :results output :exports code
  29 +export ${iota2dir}=
  30 +export terlabdir='/work/OT/theia/oso/shapes/TERLAB/'
  31 +export terlab31=${terlabdir}'SURFACES-2017-PARCELLES-GRAPHIQUES-CONSTATEES_031_20180210'
  32 +cp ${terlab31}.shp ${terlabdir}'terlab31.shp'
  33 +cp ${terlab31}.dbf ${terlabdir}'terlab31.dbf'
  34 +cp ${terlab31}.shx ${terlabdir}'terlab31.shx'
  35 +cp ${terlab31}.prj ${terlabdir}'terlab31.prj'
  36 +python ${iota2dir}/scripts/VectorTools/AddFieldID.py ${terlabdir}'terlab31.shp'
  37 +#+end_src
  38 +
  39 +
16 40 * Import des fichiers dans des DF
17 41  
18 42 Le fichier TERLAB qui contient le dept 31 est:
19 43  
20 44 #+begin_src python :results output :session :exports code
21 45 terlabdir = '/work/OT/theia/oso/shapes/TERLAB/'
22   -terlab31 = terlabdir+'SURFACES-2017-PARCELLES-GRAPHIQUES-CONSTATEES_031_20180210.shp'
  46 +terlab31 = terlabdir+'terlab31.shp'
23 47 df_terlab = gp.read_file(terlab31)
24 48 print(df_terlab.head())
25 49 #+end_src
... ...