Мне нужно объединить около 550 Гб изображений TIF вместе, и программное обеспечение, которое я пробовал, продолжает работать. Область была разделена на зоны, так что наименьшая имеет около 200 плиток.
Я использовал последние версии ERDAS (Imagine and Mapper), ArcINFO и Global Mapper на 3,30 ГГц Intel Xeon E31245, DELL, 16 ГБ оперативной памяти, 64-разрядной версии Win 7 Professional. Многоядерный (всего 4), Гиперпоточный (всего 8). Мой C имеет 700 ГБ свободного и D имеет 1,5 ТБ.
Я изучаю использование Grass (никогда раньше), но i.image.mosaic, кажется, обрабатывает только 4 файла ... некоторые из моих имеют 600 плиток. Любые другие варианты или программное обеспечение с открытым исходным кодом, чтобы попробовать?
Извините, следует добавить, что мы не можем использовать набор данных мозаики (или эквивалент в другом программном обеспечении), так как нам нужно создавать зоны с определенными областями без данных в виде ecw, чтобы их можно было открывать в любом программном обеспечении ГИС и объединять с более низким разрешением / более старым данные, когда новые данные не существуют без проблем.
Пример того, как некоторые мозаичные файлы выглядят в разных программах. Global Mapper / ERDAS - это хорошо, но это не правильно в arcgis.
--- более старая информация ---
Извините за грубый рисунок. Таким образом, наличие цветных областей в виде 5 зон минимизирует области без данных в более крупном AOI.
В arcgis код выглядит следующим образом (он запускается как модель, а не в python, так как я не могу получить его для ввода tifList).
arcpy.MosaicToNewRaster_management(tifList+";" +mask,RootOutput,"Tile1.tif","PROJCS['GDA_1994_MGA_Zone_55',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',147.0],PARAMETER['Scale_Factor',0.9996],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]","16_BIT_UNSIGNED","0.5","3","MAXIMUM","#")
# Replace a layer/table view name with a path to a dataset (which can be a layer file) or create the layer/table view within the script
# The following inputs are layers or table views: "test2"
arcpy.CopyRaster_management(OutputFile,RootOutput+"Tile1b.tif","#","256","256","NONE","NONE","16_BIT_UNSIGNED")
где tifList должен быть прочитан из CSV-файла, но это не работает в Python, поэтому я запускаю выше в модели вместо ...
У меня есть 1,5 ТБ + свободного места на диске, но процесс завершается с ошибкой 9999.
Будет ли обрабатываться даже 100 плиток? -должны ли мы смотреть на разделение зон дальше?