Le portail francophone de la géomatique

Salut

Je ne suis pas sur de comprendre...ce que tu appelles des tables sont des tables attributaires dbf ou des couches de donnees geometriques et attributaires (shapefile ou autre). Si tu as 70 shapefiles de bassin versant, pourquoi ne pas les combiner en utilisant l'outil combiner/merge, ce n'est pas un probleme s'ils se superposent, puis utilises l'outil "zonal statistic as table".

A+

Rodolphe

Simon,

Bien sur que tu peux utiliser l'outil combiner et sans rien perdre en plus. Tu peux combiner tous tes polygones de bassins versant dans un seul shapefile, meme s'ils se superposent entierement ou partiellement. Bien sur le resultat sera une couche avec une topologie incorrecte, mais ca n'est pas ton probleme et ca repondra parfaitement a ta demande ce qui est l'essentiel.

A+

Rodolphe

Simon,

OK, desole ! C'est vrai qu'en y regardant de pres, en fait l'outil statistique zonale rasterize le shapefile en un grid qu'il utilise comme raster de zones d'ou le fait que ca ne marche pas pour les polygones qui se superposent...
Une solution est d'utiliser un script python.
Voici un script que j'ai fait il y a quelques temps pour faire ce genre de choses.

#Script for Miguel to extract zonal
# statistics as table in batch from a folder
# path even if this one contain subfolders
# and a raster list
#created on: Tue 16 Sept 2008
#By Rodolphe Bernard
#From Imperial College
#     Silwood park Campus
 
#container for boolbox path list
tbx_path = []
 
#--variables to change before running this script to suits your needs and your computer settings--
#---------------------Make sure that all pathes contains slash instead of backslash------------------------
 
#COMPUTER SETTINGS
#Data management toolbox path
tbx_path.append("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx")
#Spatial Analyst toolbox path
tbx_path.append("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Spatial Analyst Tools.tbx")  
 
#YOUR SETTINGS
#path to the root folder containing all folders with shapefiles
folderpath = "C:/Documents and Settings/rbernard/My Documents/DATA/GIS_data/temp/Miguel_files/test"
#common shapefiles fieldname to use to identify the entity to group within your shapefile (like if you do a dissolve) if necessary  before performing zonal statistics
fieldname = "SHAPENAME"
#raster path list like ["rasterpath1","rasterpath2",...]
rasterpathlist = ["C:/Documents and Settings/rbernard/My Documents/DATA/GIS_data/temp/Miguel_files/ataet00"]
#scratch folder
scratch_folder = "C:/Documents and Settings/rbernard/My Documents/DATA/GIS_data/temp/Miguel_files/scratch"
#final output table
final_output = "C:/Documents and Settings/rbernard/My Documents/DATA/GIS_data/temp/Miguel_files/final_ouput_table.dbf"
 
#-------------------------------------------------------------------------------------------------
 
# Import system modules
import sys, string, os, arcgisscripting
 
# Create the Geoprocessor object
gp = arcgisscripting.create()
 
# Check out ArcGIS Spatial Analyst extension license
gp.CheckOutExtension("Spatial")
 
# Load required toolboxes...
for tbx in tbx_path:
    gp.AddToolbox(tbx)
 
#extension file to look for
searched_extension = "shp"
  
#container for filepath
fileslist = []
 
#go through all folders and subfolders and return a list of shapefiles path
for root, dirs, files in os.walk(folderpath):
    for name in files:
        if name.endswith(searched_extension) or name.endswith(searched_extension.upper()):
                #root is the dir path, root + os.sep + name is the absolute path
                fileslist.append(root + "/" + name)
#container for string list of tables
tablelist=""
#for each shapefile
for shapefile in fileslist:
    #for each raster in rasterpathlist
    for raster in rasterpathlist:
        #perform zonal statistics as table between shapefile and raster
        output_table = scratch_folder + "/" + shapefile[shapefile.rfind("/")+1:len(shapefile)-4] + "_x_" + raster[raster.rfind("/")+1:len(raster)] + ".dbf"
        try:
            #gp.ZonalStatisticsAsTable_sa(outputshapefile, fieldname, raster, output_table, "DATA")
            gp.ZonalStatisticsAsTable_sa(shapefile, fieldname, raster, output_table, "DATA")
            # Add a field to this table to be able to know with wich raster the statistics have been calculated, of type TEXT
            gp.addfield (output_table, "rastername", "TEXT", "20")
            gp.CalculateField_management(output_table, "rastername", "\""+raster[raster.rfind("/")+1:len(raster)]+"\"", "VB")
            #create a list of output table for the next operation in the following format: "table1;table2;table3;..."
            if tablelist == "":
                tablelist = output_table
            else:
                tablelist = tablelist + ";" + output_table
 
        except:
            print "Error:",
            print "could not perform the zonal statistic calculation between layers:" + shapefile[shapefile.rfind("/")+1:len(shapefile)] + " and " + raster[raster.rfind("/")+1:len(raster)]
            print "n"
        
#copy row of the first table of the list to a new table and create this table (final output table)
gp.maketableview(tablelist[0:tablelist.find(";")], "temp_table")
gp.copyrows(tablelist[0:tablelist.find(";")],final_output,"temp_table")        
#Append all other tables rows into the final output table
gp.Append_management(tablelist[tablelist.find(";")+1:len(tablelist)], final_output, "NO_TEST")

Copies et colles ce script dans un editeur de texte, sauve le avec l'extension .py. Tu as juste a changer les parametres au debut du script, mettre les bon chemins vers tes donnees, puis l'ouvrir dans Python IDLE et l'executer.
Il te faut le chemin vers ton dossier contenant tes shapefiles, le champ commun a tous tes shapefile qui identifie tes bassins versants, le chemin complet d'acces a ton raster, le chemin vers un dossier temporaire pour stocker les dbf temp (ce dossier doit exister) et pour finir le chemin vers la table dbf finale.
Le script execute ensuite en boucle le zonal stat pour chacun de tes shapefiles puis rassemble le tout dans un dbf final. Le script  est livre tel quel sans erreur handling sauf si un zonal stat ne peux etre realiser entre un shapefile et un raster (different datum, differentes emprises,...) la un message s'affiche.

A+

Rodolphe

Désolé de relancer cette dscussion...j ai lancé le script sur plusieurs bassins versants mais j'obtiens toujours ce type d'erreur :

"Error: could not perform the zonal statistic calculation between layers:bretagne5.shp and mass_mvt
n

Traceback (most recent call last):
  File "C:\Documents and Settings\barreau\Bureau\breta.py", line 87, in -toplevel-
    gp.maketableview(tablelist[0:tablelist.find(";")], "temp_table")
RuntimeError:
Echec de l'exécution. Les paramètres ne sont pas valides.
Table en entrée: valeur requise
Echec de l'exécution de (MakeTableView)."

dans le script voici les éléments que j'ai modifiés :

#YOUR SETTINGS
#path to the root folder containing all folders with shapefiles
folderpath = "D:/Donnees-sig/bassins_cotiers/bretagne"
#common shapefiles fieldname to use to identify the entity to group within your shapefile (like if you do a dissolve) if necessary  before performing zonal statistics
fieldname = "CD_STATION"
#raster path list like ["rasterpath1","rasterpath2",...]
rasterpathlist = ["D:/Donnees-sig/indicateurs/mass_mvt"]
#scratch folder
scratch_folder = "D:/Donnees-sig/bassins_cotiers/bretagne/script"
#final output table
final_output = "D:/Donnees-sig/bassins_cotiers/bretagne/script/essai.dbf"

je ne vois pas du tout pourquoi ça ne marche pas...CD station à un format texte (erreur?) j'ai créé le dossier script pour le srcatch folder mais il est vide, faut il créer la table essai.dbf même si elle est vide et ensuite lancer le script???
Donc pour résumer merci d'avance si vous avez quelques idées qui pourrait expliquer pourquoi ça ne fonctionne pas...

Simon

encore une fois merci pour l'aide,

pour le shape files j'ai :
Projection: Lambert_Conformal_Conic
Geographic Coordinate System: GCS_User_Defined
Datum: D_User_Defined

pour le raster j'ai :
Référence_spatiale>User_Defined_Lambert_Conformal_Conic
<Datum>D_User_Defined

donc je sais pas trop si c'est correct...
je précise que les stats zonales fonctionnait correctement avec le raster et les shapefiles que j'ai choisis pour le test...

@+
Simon

Pierre, ESRI a prévu les cas non directement prévu. C'est pour ça que dans ArcGIS il y a des "Model Builder" et des "langages de scripts". Cela permet de réaliser un outil spécialisé, ce n'est pas exactement une astuce, c'est juste une utilisation avancée.

Utiliser ces outils de programmation haut niveau permet de bénéficier du moteur SIG qu'est ArcGIS.

Il est très rare que pour une tache de géomatique il faut se décider à abandonner une solution SIG bureautique du style d'ArcGIS et créer un outil ex nihilo.

Grace à l'aide d un collègue j'ai enfin trouvé d'où venait le problème...
la ligne output_table = scratch_folder + "/" + shapefile[shapefile.rfind("/")+1:len(shapefile)-4] + "_x_" + raster[raster.rfind("/")+1:len(raster)] + ".dbf"
bloquait les calculs car la taille du fichier dbf crée était trop longue.
Donc encore merci à Rodolphe pour son script et ses conseils !

@+

Simon