Clip (Quản lý dữ liệu) hay Extract By Mask (Spatial Analyst) hiệu quả hơn?

Tôi đang cố gắng cắt một Orthophoto đến một ranh giới quận. Tôi đã sử dụng Extract By Mask Tool, nhưng quá trình này đã chạy được khoảng 12 giờ! Nó đã trải qua 2 chu kỳ rồi, vì vậy tôi cho rằng đây là mỗi nhóm. Ai đó có thể sửa cho tôi nếu tôi sai về điều đó.

Clip trong Quản lý dữ liệu có hoạt động tốt hơn (nhanh hơn) không? Độ chính xác sẽ như thế nào?

Tôi đã luôn cố gắng sử dụng Extract By Mask, nhưng việc này mất quá nhiều thời gian.

Tôi đã chạy thử nghiệm để xác định tốc độ và chất lượng khác nhau như thế nào giữa hai phương pháp, đây là kết quả:

Dữ liệu đầu vào

1. Hình ảnh NAIP DOQQ 4 băng tần ở định dạng .img (349,34 MB)
2. Một lớp tính năng được sử dụng làm mặt nạ / clipper

Hiệu suất

Ba thử nghiệm đã được thực hiện và điểm chuẩn. Phương pháp Clip (Quản lý dữ liệu) nhanh hơn đáng kể so với phương pháp Trích xuất bằng mặt nạ (Phân tích không gian).

Chất lượng

Cả hai phạm vi đều giống hệt như giá trị NoData. Tuy nhiên, một đánh giá trực quan cho thấy rằng trích xuất bằng phương pháp mặt nạ đã thay đổi một chút cách sắp xếp pixel - có thể là kết quả của một số kiểu lấy mẫu lại. Bố trí pixel trong thao tác clip giống hệt với hình ảnh đầu vào ban đầu.

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