Добрый день, в данный момент я работаю над плагином для Metashape 2.2.1.
1) Чтобы избежать экспортирования ортофотопланов, плагин принимает на вход тайлы из проекта с помощью Python кода приведенного ниже и создает из них виртуальный датасет(GDAL .vrt). При попытке прочитать эти тайлы с помощью GDAL или rasterio считывается только overview самого низкого качества. Поизучав причины выяснил, что это связанно с тем, что тайлы в проекте хранят самое низкое разрешение на 1-й странице мультистраничного .tif, а GDAL ожидает, что 1-я страница хранит высокое качество, а остальные overview. Предусмотрена ли возможность чтения тайлов высокого разрешения из проекта без экспорта совместимая с GDAL?
app = Metashape.app
doc = app.document
tiles = [[tile['tile_path'] for tile in getattr(orthomosaic, '__tiles') if 'l5' not in Path(tile['tile_path']).stem] for orthomosaic in self.chunk.orthomosaics]
2) В тайлах проекта в GeoTransform каждого тайла внутри ортофотоплана записано одинаковое географическое преобразование (похоже преобразование родительского ортофотоплана). Есть ли возможность считывать правильное географическое преобразование для тайлов из проекта?
Мне удалось получить решение, при котором виртуальный датасет создается по 9-й странице каждого тайла, а географические преобразования для тайлов пересчитываются исходя из их x, y, позоции указанной в имени и orthomosaic/doc.xml, но это требует дополнительных вычислений и не очень интуитивно.
С уважением,
Гонтарь Денис
Topic: Несовместимость внутреннего формата хранения тайлов ортофотоплана с GDAL (Read 270 times)
