Czy Clip (zarządzanie danymi) lub Extract By Mask (Spatial Analyst) jest bardziej wydajny?

Próbuję wyciąć ortofoto do granicy hrabstwa. Korzystam z narzędzia Wyodrębnij według maski, ale proces ten trwa około 12 godzin! Przeszły już 2 cykle, więc zakładam, że to każdy zespół. Czy ktoś może mnie poprawić, jeśli się mylę.

Czy klip w zarządzaniu danymi działałby lepiej (szybciej)? Jaka byłaby dokładność?

Zawsze próbowałem użyć Extract By Mask, ale to trwa zbyt długo.

Przeprowadziłem test, aby ustalić, w jaki sposób prędkość i jakość różnią się między dwiema metodami, oto wyniki:

Dane wejściowe

1. 4-pasmowy obraz NAIP DOQQ w formacie .img (349,34 MB)
2. Klasa funkcji używana jako maska ​​/ maszynka do strzyżenia

Występ

Przeprowadzono trzy próby i przeprowadzono testy porównawcze. Metoda Clip (zarządzanie danymi) jest znacznie szybsza niż metoda Extract by Mask (Spatial Analyst).

Jakość

Oba zakresy były identyczne, podobnie jak wartości NoData. Ocena wizualna wykazała jednak, że metoda ekstrakcji metodą maski nieznacznie zmieniła układ pikseli - prawdopodobnie wynik pewnego rodzaju ponownego próbkowania. Układ pikseli w operacji klipu był identyczny z oryginalnym obrazem wejściowym.

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