Collecter des événements avec rendu (Statistiques spatiales)
Récapitulatif
Convertit les données d'événement en données ponctuelles pondérées, puis applique un rendu par cercles gradués au champ de total résultant.
Héritage :En raison des nouvelles fonctionnalités d'ArcGIS qui permettent d'associer la sortie d'outils de script et de modèle à un rendu par défaut, cet outil sera déconseillé après la version 10.1 d'ArcGIS.
Illustration
 |
Utilisation
-
L'outil Collecter des événements avec rendu est un modèle qui combine l'outil Collect Events avec l'outil Rendu d'effectif.
Depuis la version ArcGIS 9.3, la classe d'entités en sortie de l'outil Collect Events est automatiquement ajoutée à la table des matières avec un rendu par cercles gradués par défaut appliqué au champ ICOUNT.
Syntaxe
| Paramètre | Explication | Type de données |
Input_Incident_Features |
Entités représentant les données d'événement ou d'incident. | Feature Layer |
Output_Layer |
Fichier de couche (.lyr) permettant de stocker les informations de rendu par cercles gradués. | Layer File |
Output_Weighted_Point_Feature_Class |
Classe d'entités en sortie contenant les données ponctuelles pondérées. | Feature Class |
Exemple de code
Le script de fenêtre Python ci-dessous illustre l'utilisation de l'outil Collecter des événements avec rendu.
import arcpy
arcpy.env.workspace = "C:/Data"
arcpy.CollectEventsRendered_stats("911.shp", "911Count_rendered.lyr", "911Count.shp")
Le script Python autonome ci-dessous illustre l'utilisation de l'outil Collecter des événements avec rendu.
# Analyze the spatial distribution of 911 calls in a metropolitan area
# using the Hot-Spot Analysis Tool (Local Gi*)
# Import system modules
import arcpy
# Set geoprocessor object property to overwrite existing output, by default
arcpy.gp.overwriteOutput = True
# Local variables...
workspace = "C:/Data"
try:
# Set the current workspace (to avoid having to specify the full path to the feature classes each time)
arcpy.env.workspace = workspace
# Copy the input feature class and integrate the points to snap
# together at 500 feet
# Process: Copy Features and Integrate
cf = arcpy.CopyFeatures_management("911Calls.shp", "911Copied.shp",
"#", 0, 0, 0)
integrate = arcpy.Integrate_management("911Copied.shp #", "500 Feet")
# Use Collect Events to count the number of calls at each location
# Process: Collect Events
ce = arcpy.CollectEventsRendered_stats("911Copied.shp", "911Count_rendered.lyr", "911Count.shp")
# Add a unique ID field to the count feature class
# Process: Add Field and Calculate Field
af = arcpy.AddField_management("911Count.shp", "MyID", "LONG", "#", "#", "#", "#",
"NON_NULLABLE", "NON_REQUIRED", "#",
"911Count.shp")
cf = arcpy.CalculateField_management("911Count.shp", "MyID", "[FID]", "VB")
# Create Spatial Weights Matrix for Calculations
# Process: Generate Spatial Weights Matrix...
swm = arcpy.GenerateSpatialWeightsMatrix_stats("911Count.shp", "MYID",
"euclidean6Neighs.swm",
"K_NEAREST_NEIGHBORS",
"#", "#", "#", 6,
"NO_STANDARDIZATION")
# Hot Spot Analysis of 911 Calls
# Process: Hot Spot Analysis (Getis-Ord Gi*)
hs = arcpy.HotSpots_stats("911Count.shp", "ICOUNT", "911HotSpots.shp",
"GET_SPATIAL_WEIGHTS_FROM_FILE",
"EUCLIDEAN_DISTANCE", "NONE",
"#", "#", "euclidean6Neighs.swm")
except:
# If an error occurred when running the tool, print out the error message.
print arcpy.GetMessages()