<\/HEAD>Les ensembles de données raster contiennent une grande variété d'informations au-delà de l'affichage basique des pixels. Ces informations sont stockées dans les propriétés et sont utiles pour mieux comprendre les données. Localiser les propriétés d'un ensemble de données raster peut être un processus fastidieux si vous essayez de compiler des informations pour de nombreux ensembles de données. En effet, le flux de travail pour trouver les propriétés d'un ensemble de données raster dans le Catalog est de faire un clic droit sur chaque ensemble de données et choisir Propriétés > onglet Général.
<\/A>
Par exemple, vous travaillez sur un projet qui nécessite la création d'un service d'image d'élévation. La base du service d'image est un
mosaic dataset qui référence les rasters source pour chaque comté de l'état. Vous créez le mosaic dataset, ajoutez les rasters, construisez les aperçus, et l'ajoutez au document cartographique. Le mosaic dataset semble parfait jusqu'à ce que vous commenciez à zoomer sur les données sources. À un endroit, vous constatez que l'élévation apparaît très sombre au lieu d'une échelle de gris cohérente.
<\/A>
Après avoir examiné le raster source, vous découvrez que l'ensemble de données est un entier 4 bits au lieu du flottant 32 bits normal. Cela vous fait vous demander s'il existe d'autres ensembles de données créés incorrectement. Pour trouver ces ensembles, vous pouvez soit zoomer sur des zones à petite échelle à travers tout l'état, soit ouvrir les propriétés pour chaque raster individuel. Dans tous les cas, le processus serait très chronophage car chaque raster devrait être examiné individuellement.
Une autre option est d'accéder aux propriétés raster via une boucle Python. Ce processus offre une solution plus efficace et simplifiée. Les informations peuvent être interrogées via l'objet python Raster object et l'outil de géotraitement
Get Raster Properties.
[python]
import arcpy
arcpy.env.workspace = r"C:\\temp"
rasterList = arcpy.ListRasters()
for raster in rasterList:
rasterObj = arcpy.Raster(raster)
print raster
bands = arcpy.GetRasterProperties_management(raster, "BandCount")
print "Nombre de bandes : %s" %bands
noData = rasterObj.noDataValue
print "Valeur NoData : %s \n" %noData
[\/python]
<\/A>
Les valeurs retournées peuvent être écrites dans un tableau pour identifier facilement quels ensembles de données sont incorrects.
<\/A>Ressources supplémentaires<\H2>
Un exemple d'outil script Python pour obtenir les propriétés raster peut être téléchargé depuis ArcGIS.com :
Write Raster Properties to Table.
Timothy H. – Analyste Support Senior<\EM>
<\BODY><\HTML>
