Створення формфайлу, що показує сліди растрових?

У мене є близько 1000 супутникових знімків у форматі tiff, і я хочу створити файл форми, який буде слугувати в якості покажчика растрам. Це щось подібне до растрового каталогу, але я не хочу будувати растровий каталог.

Деякі перешкоди, які я можу передбачити, полягають у тому, що зображення є георефрендами, тому вони не мають прямокутної форми (я кажу про область даних).

Для уточнення я вимагаю, щоб багатокутник охоплював лише ненульові (або не нодатні) пікселі растра, а не весь прямокутний растр. Більшість відповідей поки що дає прямокутний багатокутник, який охоплює дані, а також непікселі.

Мій образ

Результат, отриманий інструментами, які я перевірив (на кшталт растрових каталогів, різних архівних скриптів, користувацького сценарію Python, наведеного в одній з відповідей):

Результат, який я хочу:

У QGIS є плагін, який називається Image Boundary . Це чудовий інструмент. У цьому інструменті є опція "Дійсні пікселі", яка опускає, наприклад, чорні краї супутникового зображення.

Наступний код візьме растровий вхід, отримає його масштаб і вставить цю міру в клас класів багатокутника:

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

Ви можете запустити його у вікні ArcMap Python, встановивши in_rasterі fcтак:

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

де r_extentє існуючий багатокутний клас класу. Потім просто скопіюйте код і запустіть його. Я отримую це:

Ви можете використовувати gdaltindex для цього: http://www.gdal.org/gdaltindex.html

Однак він все одно створюватиме прямокутники (наприклад, 4 + 1 бали) у тій же системі відліку, що і зображення. Але мені цікаво, чи це насправді проблема: наскільки великі ваші образи?

Я спробував інструмент "Межі зображення" в QGIS, але він не вдається з деякими растрами.

Інший варіант, який працював для мене, це: - 1. Перетворити растровий вхід в єдине цифрове число за допомогою утиліти командного рядка gdal_translate.

gdal_translate -co компрес = lzw -b 1 -ot byte -scale 1 1 inp.tif out.tif

Це створить растр зі значеннями 1 скрізь, коли .p.tif мав ненульові значення.

2. Потім в ArcMap ви можете скористатися інструментом «растр до багатокутника», щоб отримати полігон (shapefile).

Форма виводу файлу матиме не чотирикутні межі з декількома вершинами багатокутника за фактичною межею зображення.

Якщо вас цікавить сценарій рішення, спробуйте команду gdal nearblack . Ви можете використовувати прапор -белі для отримання даних, обмежених пробілом. Ви можете полігонізувати вихід, щоб створити не прямокутний слід растрових даних.

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

Я вважаю, що в світі ArcGIS ви можете досягти цього за допомогою інструмента 3D Analyst \ Conversion \ From Raster \ Raster Domain Domain .

У ArcGIS 3D Analyst є інструмент, який називається Raster Domain, який створює полігон або полілінійний відбиток частин даних растрового набору даних і ігнорує комірки NoData. Для мене чудово працювали.

Ви можете створити растровий каталог із присвоєною просторовою посиланням, додати каталог растра до ArcMap, а потім "Експорт сліду". Пам'ятайте, що є можливість створити растровий каталог "НЕМАГАЛЬНИЙ", щоб ви не робили дубліката кожного растрового набору даних.