Le portail francophone de la géomatique

#Import modules
import arcgisscripting, sys, os, string

# Create the geoprocessor object
gp = arcgisscripting.create()

# Check out any necessary licenses
gp.CheckOutExtension("Spatial")
gp.cellSize = 1

gp.Workspace = "D:\\Test"

try:

    #----------BUFFER------------------------------------
    # Zone des Buffers 
    InRasterBuffer = "D:\\Test\\buffer\\EssaiNDVI_2.tif"
    print InRasterBuffer

    zone_field = "Value"
    
    #---------RASTER-----------------------------------
    # Get a list of tif in the workspace.
    InputRasters = gp.ListRasters("", "TIF")
    Raster_uniq = InputRasters.next()
    i=0
    # boucle de lecture des raster
    while Raster_uniq:
        print Raster_uniq
        i=i+1
        Nom_Image = Raster_uniq[0:8] # sera ultérieurement rajouté dans la table pour identifier les résultats
        print Nom_Image

        output = "ZonalStats" + str(i) + ".dbf"
        print output
        
        # Process: Zonal Statistics as Table
        gp.ZonalStatisticsAsTable_sa(InRasterBuffer, zone_field, Raster_uniq, output, "DATA")
        print "image " + Nom_Image + " traitée"

        # retour boucle : Raster suivant
        Raster_uniq = InputRasters.next()
        print "C'est fini"
        
except:
    gp.AddMessage(gp.GetMessages())
    print gp.GetMessages()