Мені потрібно зробити мозаїку близько 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 ТБ.
Я шукаю використання трави (ніколи раніше), але i.image.mosaic, здається, обробляє лише 4 файли ... деякі з моїх мають 600 плиток. Будь-які інші варіанти чи програмне забезпечення з відкритим кодом, щоб спробувати?
Вибачте, слід додати, що ми не можемо використовувати мозаїчний набір даних (або еквівалент в іншому програмному забезпеченні), оскільки нам потрібно створити зони з визначеними областями без даних як ecw, щоб їх можна було відкрити в будь-якому програмному забезпеченні ГІС та поєднати з меншою роздільною здатністю / старшою дані, коли нові дані існують непросто.
Приклад того, як деякі мозаїчні файли виглядають у різних програмних засобах. Глобальний Mapper / ERDAS є нормальним, але це не є правильним у аркгізі.
--- СТАРІША ІНФОРМАЦІЯ ---
Вибачте за чорновий малюнок. Тому використання кольорових областей як 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 плиток? -Так ми повинні дивитись на розбиття зон далі?