Klip (Veri yönetimi) veya Maske İle Ayıkla (Uzamsal Analist) daha verimli mi?

Bir Ortofotoyu bir ilçe sınırına indirmeye çalışıyorum. Maske İle Ayıkla aracını kullanıyorum, ancak bu işlem yaklaşık 12 saattir devam ediyor! Zaten 2 döngüden geçti, bu yüzden her grubun bu olduğunu varsayıyorum. Bu konuda yanılıyorsam birisi beni düzeltebilir mi?

Veri Yönetimindeki Klip daha iyi (daha hızlı) çalışır mı? Doğruluk nasıl olurdu?

Hep Mask By Extract'ı kullanmaya çalıştım, ama bu çok uzun sürüyor.

Hız ve kalitenin iki yöntem arasında nasıl farklı olduğunu belirlemek için bir test yaptım, işte sonuçlar:

Veri girişi

1. .İmg formatında 4 bantlı NAIP DOQQ görüntüsü (349.34MB)
2. Maske / kesme makinesi olarak kullanılan bir özellik sınıfı

Verim

Üç deneme yapıldı ve kıyaslandı. Klip (Veri Yönetimi) yöntemi Maske ile Ayıkla (Uzamsal Analist) yönteminden önemli ölçüde daha hızlıdır.

Kalite

Her iki uzantı da NoData değerleri ile aynıdır. Bununla birlikte, görsel bir değerlendirme, maske yöntemiyle ekstraktın piksel düzenlemesini biraz değiştirdiğini - muhtemelen bir tür yeniden örneklemenin sonucunu gösterdi. Klip işlemindeki piksel düzenlemesi orijinal giriş görüntüsüyle aynıydı.

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