Ist Clip (Datenverwaltung) oder Nach Maske extrahieren (Spatial Analyst) effizienter?

Ich versuche, ein Orthofoto an einer Kreisgrenze zu befestigen. Ich habe das Extract By Mask Tool verwendet, aber dieser Vorgang wurde ungefähr 12 Stunden lang ausgeführt! Es hat bereits 2 Zyklen durchlaufen, also gehe ich davon aus, dass dies jede Band ist. Kann mich jemand korrigieren, wenn ich mich irre?

Wäre der Clip in der Datenverwaltung besser (schneller)? Wie wäre die Genauigkeit?

Ich habe immer versucht, Extract By Mask zu verwenden, aber das dauert zu lange.

Ich habe einen Test durchgeführt, um festzustellen, wie sich die Geschwindigkeit und die Qualität der beiden Methoden unterscheiden. Hier sind die Ergebnisse:

Eingabedaten

1. 4-Band-NAIP-DOQQ-Bild im .img-Format (349,34 MB)
2. Eine Feature-Class, die als Maske / Clipper verwendet wird

Performance

Drei Versuche wurden durchgeführt und bewertet. Die Methode "Clip" (Datenverwaltung) ist erheblich schneller als die Methode "Nach Maske extrahieren" (Spatial Analyst).

Qualität

Beide Bereiche waren identisch, ebenso wie die NoData-Werte. Eine visuelle Beurteilung ergab jedoch, dass die Pixelanordnung durch das Maskenextraktverfahren leicht verändert wurde - wahrscheinlich das Ergebnis einer Art von Resampling. Die Pixelanordnung bei der Clip-Operation war identisch mit dem ursprünglichen Eingabebild.

# 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)