Bien que j'utilise python pour créer des scripts / services de géotraitement, j'avais l'impression que l'utilisation d'ArcObjects pour effectuer les opérations équivalentes aurait de meilleures performances.
J'ai publié le service ArcGIS Server GP - RasterIO.dll plante ArcSOC.exe et le script de géotraitement ArcGIS fonctionne correctement sur le bureau, mais se bloque en tant que service de géotraitement? au cours des derniers jours sur l'obtention de scripts de géotraitement qui utilisent les outils Spatial Analyst pour fonctionner en tant que services de géotraitement. Mon échéance approche à grands pas, j'ai donc décidé de suivre la voie SOE pour atteindre la fonctionnalité souhaitée.
Obtenir une analyse du chemin de coût dans ArcObjects était relativement simple à l'aide des méthodes .NET ESRI.ArcGIS.SpatialAnalyst.RasterDistanceOpClass , en particulier les méthodes CostDistanceFull () et CostPath ().
Quelques extraits de code de la façon dont je fais les choses:
Python
# Get Cost Path Origin and Destination Points
inputPointsShp = 'D:/RasterStuff/test_points.shp'
arcpy.MakeFeatureLayer_management(inputPointsShp,"origin",' "TYPE" = \'ORIGIN\' ')
arcpy.MakeFeatureLayer_management(inputPointsShp,"destination",' "TYPE" = \'DESTINATION\' ')
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute CostDistance
outCostDistance = CostDistance("origin",SOURCE_RASTER,"#","backlink")
# Execute CostPath
outCostPath = CostPath("destination", outCostDistance,"backlink")
# Convert Result to Polyline
arcpy.RasterToPolyline_conversion(outCostPath, "leastCostPath")
featSet = arcpy.FeatureSet("leastCostPath")
C #
IDistanceOp distanceOp = new RasterDistanceOpClass();
IRasterBandCollection costDistanceRaster = (IRasterBandCollection)distanceOp.CostDistanceFull((IGeoDataset)sourceFc, (IGeoDataset)raster, true, true, false);
IRasterBand distanceRaster = costDistanceRaster.Item(0);
IRasterBand backLinkRaster = costDistanceRaster.Item(1);
IGeoDataset costPath = distanceOp.CostPath((IGeoDataset)destFc, (IGeoDataset)distanceRaster, (IGeoDataset)backLinkRaster, ESRI.ArcGIS.SpatialAnalyst.esriGeoAnalysisPathEnum.esriGeoAnalysisPathForEachCell);
Une analyse de la trajectoire des coûts dans ArcPy (en utilisant sa.CostDistance et sa.CostPath) prend environ 15-20 secondes. En utilisant exactement les mêmes entrées, la routine basée sur ArcObjects prend 55 à 60 secondes. Même l'utilisation du géoprocesseur .NET est beaucoup plus lente que l'arcpy.
Je suppose que mes questions ici sont:
- Les implémentations ArcPy et ArcObjects pointent-elles vers la même base de code (via leurs wrappers Python et .NET)?
- Des conseils pour optimiser l'analyse des chemins de coûts basée sur ArcObject?