Est-ce que Clip (Gestion des données) ou Extraire par masque (Spatial Analyst) est plus efficace?

J'essaie de couper une orthophoto à la limite d'un comté. J'utilise l'outil Extraire par masque, mais ce processus fonctionne depuis environ 12 heures! Il a déjà traversé 2 cycles, donc je suppose que c'est chaque groupe. Quelqu'un peut-il me corriger si je me trompe à ce sujet.

Le clip dans la gestion des données fonctionnerait-il mieux (plus rapidement)? Quelle serait la précision?

J'ai toujours essayé d'utiliser Extract By Mask, mais cela prend trop de temps.

J'ai effectué un test pour déterminer en quoi la vitesse et la qualité diffèrent entre les deux méthodes, voici les résultats:

Des données d'entrée

1. Image NAIP DOQQ 4 bandes au format .img (349,34 Mo)
2. Une classe d'entités utilisée comme masque / tondeuse

Performance

Trois essais ont été réalisés et comparés. La méthode Clip (Data Management) est nettement plus rapide que la méthode Extract by Mask (Spatial Analyst).

Qualité

Les deux étendues étaient identiques, tout comme les valeurs NoData. Cependant, une évaluation visuelle a montré que la méthode d'extrait par masque modifiait légèrement la disposition des pixels - probablement le résultat d'un certain type de rééchantillonnage. La disposition des pixels dans l'opération de clip était identique à l'image d'entrée d'origine.

# Import system modules
import arcpy, time, os
from arcpy import env
from arcpy.sa import *

env.overwriteOutput = 1

# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")

# Set local variables
inRaster = r'C:\temp\naip2011.img'
inMaskData = r'C:\temp\fgdb.gdb\clipper'
outws = r'C:\temp'
out1 = os.path.join(outws, 'extractbymask.img')
out2 = os.path.join(outws, 'clip.img')

#############TRIAL 1: EXTRACT BY MASK#########################################
start = time.clock()

# ExtractByMask
outExtractByMask = ExtractByMask(inRaster, inMaskData)

# Save the output
outExtractByMask.save(out1)

end = time.clock()
total = end - start

print "The extract by mask method took:  %s seconds" % round(total, 3)

###############TRIAL 2: CLIP#################################################
start = time.clock()

# Clip
arcpy.Clip_management(inRaster, "#", out2, inMaskData, "", "ClippingGeometry")

end = time.clock()
total = end - start

print "The clip method took:             %s seconds" % round(total, 3)