Как организовать построение пирамид изображений?

Question

[luzhskij]

Здравия.
Возникла необходимость в построении пирамид изображений.
И...кто может посоветовать с какой стороны подступиться? По идее, opencv умеет делать изменение размера и интерполяцию, но мои файлы могут быть размером и 10+гб.
В своем вопросе ниже как раз и выяснял, каким объемом можно оперировать. И если говорить об открытии по кускам... То непонятно, как это вообще может выглядеть%)
Работаю под windows.

Answer 1

[Koss1024]

Это уже очень большое изображение. Стоит его разбить на тайлы со стороной степени двойки
И их использовать как основание пирамиды.

Элементарной интерполяций получить тайны вчетверо меньше и т/д

Собсно все

Но только мне кажется что я неправильно понял вопрос поскольку это весьма очевидно :)

Answer 2

[Михаил Григорьев]

Попробуйте воспользоваться библиотекой GDAL - в ней есть средства для работы с большими изображениями.

Например, для организации просмотра больших изображений можно через GDAL запрашивать тайлы нужного масштаба (LOD-а) изображения и отображать их (GDAL самостоятельно будет масштабировать изображение и строить тайлы). Для повышения производительности этого процесса есть возможность построения overviews (версий изображения меньших масштабов) с помощью GDAL – рядом с файлом GDAL помещает файл с пирамидой и потом сам использует его, когда у него запрашивают данные в меньшем разрешении.

Кроме того, можно сохранять overviews в отдельные файлы, потом резать их на тайлы (тоже GDAL-ом), получится требуемая пирамида масштабов.

У GDAL-а есть с++ и cи api. Еще c# и питон, cmd, ну и что там до кучи.

Минусы: мало примеров использования в интернете, придется повозиться с документацией и поэкспериментировать.

Comments

[luzhskij]

Да, соглашусь с Вами - gdal хорош. И именно с его помощью и зачитываю файлы. Он умеет очень много. Но весь потенциал его пока не выходит раскрыть. Например, как нарезать тайлы по уровням и расположить их в своем порядке? В порядке похожем на openstreetmap. Т.е. уровень пирамиды(папка)->номер столбика(папка)->номер строчки(файл). Или как-то можно обращаться к overviews в процессе тайлирования?

[Михаил Григорьев]

Указанную нарезку тайлов несложно сделать через API вручную. То есть рассчитывать пирамиду вручную, запрашивать данные нужного размера с нужным прореживанием и сохранять их в требуемую структуру папок. Возможно ли это сделать через cmd - не знаю. С overviews, если не стоит вопрос о производительности, можно вообще не заморачиваться (а можно построить их один раз и GDAL будет сам их автоматически использовать).

[luzhskij]

"То есть рассчитывать пирамиду вручную, запрашивать данные нужного размера с нужным прореживанием" - вот очень любопытно. Я на C++ веду разработку. Сейчас через opencv делаю пирамиды, т.е. отрезается тайл - интерполируется, уменьшается, сохраняется и так пока всё изображение не обработается на уровне пирамиды, потом переходим на второй уровень и опять.... И это очень медленно.
Наверное, то как Вы написали - будет быстрее. Был бы очень благодарен за указание где можно об этом прочесть.

[Михаил Григорьев]

Вам нужна функция RasterIO www.gdal.org/classGDALDataset.html#ae077c53268d227...
Она весьма замороченная. Вот пример использования:

Rectangle rect ( 0, 0, resultBitmap->Width, resultBitmap->Height );
BitmapData^ bmpData = resultBitmap->LockBits(rect, ImageLockMode::ReadWrite,
PixelFormat::Format24bppRgb);

// Get the address of the first line.
int bands[3] = { 3, 2, 1 };
IntPtr ptr = bmpData->Scan0;
pDataset->RasterIO( GF_Read, 0, 0, srcSize.Width, srcSize.Height, (unsigned char*)ptr.ToInt32(),
resultBitmap->Width, resultBitmap->Height, GDT_Byte, 3, bands, 3, bmpData->Stride, 1 );

Здесь все большое изображение считывается в bitmap какого-то (меньшего размера).

[Михаил Григорьев]

int bands[3] = { 3, 2, 1 };
означает, что нужно цветовые каналы брать в порядке BGR.

[luzhskij]

Longy: Спасибо, буду пробовать. Попутно нашел GDAL Warp. Но тоже непонятно, как применить.