Tạo shapefile hiển thị dấu chân của Rasters?

Tôi có khoảng 1.000 hình ảnh vệ tinh ở định dạng tiff và tôi muốn tạo một shapefile sẽ đóng vai trò là một chỉ mục cho các trình quét. Đây là một cái gì đó tương tự như một danh mục raster, nhưng tôi không muốn xây dựng một danh mục raster.

Một số trở ngại tôi có thể thấy trước, đó là hình ảnh được định vị địa lý, vì vậy chúng không có hình chữ nhật (tôi đang nói về vùng dữ liệu).

Để làm rõ, tôi yêu cầu đa giác chỉ bao gồm các pixel khác không (hoặc không gật đầu) của raster chứ không phải toàn bộ raster hình chữ nhật. Hầu hết các câu trả lời cho đến nay, đưa ra một đa giác hình chữ nhật, bao gồm dữ liệu, cũng như các pixel không dữ liệu.

Hình ảnh của tôi

Kết quả được cung cấp bởi các công cụ tôi đã kiểm tra (như catlog raster, Arcscripts khác nhau, tập lệnh Python tùy chỉnh được đưa ra trong một trong các câu trả lời):

Kết quả mà tôi muốn:

Có một plugin trong QGIS được gọi là Image Boundary . Nó là một dụng cụ hữu ích. Trong công cụ này, có một tùy chọn cho "Điểm ảnh hợp lệ", ví dụ như sẽ bỏ qua các cạnh màu đen của hình ảnh vệ tinh.

Đoạn mã sau sẽ lấy một raster đầu vào, lấy phạm vi của nó và chèn phạm vi đó vào một hình đa giác:

import arcpy

r = arcpy.Raster(in_raster)
point = arcpy.Point()
array = arcpy.Array()

corners = ["lowerLeft", "lowerRight", "upperRight", "upperLeft"]

cursor = arcpy.InsertCursor(fc)
feat = cursor.newRow()
for corner in corners:    
    point.X = getattr(r.extent, "%s" % corner).X
    point.Y = getattr(r.extent, "%s" % corner).Y
    array.add(point)
array.add(array.getObject(0))
print len(array)
polygon = arcpy.Polygon(array)
feat.shape = polygon
cursor.insertRow(feat)
array.removeAll()
del feat
del cursor

Bạn có thể chạy nó trong cửa sổ ArcMap Python bằng cách thiết lập in_rastervà fcnhư vậy:

>>> fc = 'r_extent'
>>> in_raster = 'CaliCoast10mNED_HavCurvedPCS'

nơi r_extentcó một đa giác featureclass hiện có. Sau đó, chỉ cần sao chép mã và chạy nó. Tôi nhận được điều này:

Bạn có thể sử dụng gdaltindex cho việc này: http://www.gdal.org/gdaltindex.html

Tuy nhiên, nó vẫn sẽ tạo ra các hình chữ nhật (ví dụ 4 + 1 điểm) trong cùng hệ thống tham chiếu với hình ảnh. Nhưng tôi tự hỏi liệu đó có thực sự là một vấn đề: hình ảnh của bạn lớn đến mức nào?

Tôi đã thử công cụ "Ranh giới hình ảnh" trong QGIS nhưng không thành công với một số trình quét.

Một tùy chọn khác phù hợp với tôi là: - 1. Chuyển đổi raster đầu vào thành một số kỹ thuật số bằng cách sử dụng tiện ích dòng lệnh gdal_translate.

gdal_translate -co nén = lzw -b 1 -ot byte -scale 1 1 inp.tif out.tif

Điều này sẽ tạo raster với các giá trị 1 ở mọi nơi inp.tif có các giá trị khác không.

2. Sau đó, trong ArcMap, bạn có thể sử dụng công cụ "raster to polygon" để lấy đa giác (shapefile).

Shapefile đầu ra sẽ có các ranh giới không phải là tứ giác, với nhiều đỉnh đa giác, theo ranh giới thực của hình ảnh.

Nếu bạn quan tâm đến một giải pháp theo kịch bản, hãy thử lệnh gdal nearblack . Bạn cũng có thể sử dụng cờ - để lấy dữ liệu giới hạn bởi khoảng trắng. Bạn có thể đa giác hóa đầu ra để tạo ra một dấu chân không phải hình chữ nhật của dữ liệu raster.

##run nearblack
shellcmds = "nearblack -white -o " + outpath + " " + inpath
os.system(shellcmds)

##run polygonize
shellcmds = "gdal_polygonize.py " + outpath + " -f" + """ "ESRI Shapefile"  """ + outpath2 + " nearblack"
os.system(shellcmds)

Tôi tin rằng, trong thế giới ArcGIS, bạn có thể đạt được nó bằng công cụ Phân tích 3D \ Chuyển đổi \ Từ Raster \ Raster Domain .

Có một công cụ trong ArcGIS 3D Analyst gọi là Raster Domain , tạo ra dấu chân đa giác hoặc đa giác của các phần dữ liệu của bộ dữ liệu raster - và bỏ qua các ô NoData. Làm việc tuyệt vời cho tôi.

Bạn có thể tạo một danh mục raster với tham chiếu không gian được chỉ định, thêm danh mục raster vào ArcMap và sau đó "Xuất dấu chân". Hãy nhớ rằng có tùy chọn để tạo một danh mục raster "UNMANAGED" để bạn không tạo một bản sao của mỗi bộ dữ liệu raster.