METHOD OF TEXTURE SEGMENTATION OF AEROSPACE IMAGES

Full Text

Введение В системах аэрокосмического мониторинга земной поверхности для решения задач анализа изображений широко применяется сегментация, которая позволяет выделить на изображениях интересующие области. Все известные методы сегментации изображений можно условно разделить на методы контурной сегментации и методы, основанные на выделении однородных (текстурных) областей. Методы на основе выделения контуров обладают высоким быстродействием, но не всегда оказываются приемлемыми при выделении на изображении областей, не имеющих четко выраженных границ и лишенных существенных деталей. В этом случае необходимо использовать идентификацию областей изображений на основе анализа их текстуры. В [1] показано, что для описания текстуры могут быть использованы марковские случайные поля. Однако при обработке аэрокосмических изображений, представляющих собой полутоновые (или цветные) изображения с числом уровней яркости 2^, возникает проблема хранения в памяти и оперирования с матрицами вероятностей переходов размерностью 28 х28 . Такая обработка изображений требует больших вычислительных ресурсов. Представление ^-разрядных полутоновых изображений g-разрядными двоичными изображениями (РДИ) в [2-3] позволило «Инфокоммуникационные технологии» Том 10, № 4, 2012 42 Карлушин К.А., Курбатова Е.Е., Медведева Е.В. снизить вычислительные ресурсы за счет оперирования с матрицами вероятностей переходов размером 2x2. В статье предложен метод текстурной сегментации аэрокосмических изображений на основе марковских случайных полей, позволяющий при меньших вычислительных ресурсах эффективно выделять протяженные области с однородными статистическими характеристиками. Метод текстурной сегментации Пусть g-разрядное цифровое полутоновое изображение (ЦПИ) представлено набором g-РДИ. Учитывая характер статистической связи между элементами изображения, будем считать, что элементы ЦПИ образуют двумерную цепь Маркова с несколькими состояниями, а РДИ - двумерную цепь Маркова с двумя равновероятными состояниями м[1), мр и матрицами вероятностей переходов по горизонтали 'П = 2x2 и вертикали П = 2я« пу 2x2 На рис. 1 представлен фрагмент двумерного /-го РДИ, соответствующего марковскому случайному полю. Элементами ближайшей окрес- (0 тности анализируемого элемента vy являются (0 (0 элементы V} v^'. Так как однородная двумерная цепь Маркова полностью определяется матрицами вероятностей переходов вместе с вектором начальных вероятностей, то в качестве текстурного признака предлагается использовать оценку вероятностей перехода в двумерной цепи Маркова. Предполагается, что текстурный признак внутри одной области изменяется незначительно, а на разных областях принимает существенно различные значения. Учитывая, что основные детальные области выражены на старших разрядах ЦПИ, для их выделения необходимо обработать РДИ старших разрядов ЦПИ с наиболее выраженными текстурными признаками. тт (0 Для определения вероятности перехода пу в двумерной цепи Маркова воспользуемся многомерной математической моделью изображения, синтезируемой в [2]. Вероятность состояния элемента (см. рис. 1) определится энтропией (0 элемента ' относительно состояний элементов ближайшей окрестности v^j: #(V^ v(2,),v|')j = -log (1) где j, j - одномерные плот ности вероятностей перехода между соседними ( (0 (0 (0^ значениями; wlv^' v}',vv2' I - плотность вероятности перехода в двумерной цепи Маркова. Плотность вероятности перехода в двумерной двоичной цепи Маркова можно выразить как hj,q=\ xSCvj0 -Mf )x5(v^ -M<°), (2) где 6(0 - дельта функция. На основании аргумента (1) матрица для сочетаний состояний эле-А«: ментов окрестности Ау имеет вид: П (0= af <4° af ”52 л2 а« af л2 л2 а(4г) af i,j = 1,2; і* Г, / = l,g; (3) элементы которой симметричны и удовлетворя- (0 , '(0 і ют условию нормировки ' + a^w = 1. Элементы матрицы (3) связаны с элементами матриц *11^ и 2П^ как где - элемент матрицы вероятностей переходов 3П^ = х 2П^. Вычисление оценки вероятности перехода по горизонтали сводится к вычислению средней длины последовательности элементов одного знака [3]: Рис. 1. Фрагмент двоичного изображения (4) «Инфокоммуникационные технологии» Том 10, № 4, 2012 Карлушин К.А., Курбатова Е.Е., Медведева Е.В. 43 лС'.'О где X - оценка средней длины последовательности одинаковых элементов 7-го разряда на г -ом шаге уточнения оценки. Начиная со второй строки РДИ, используя множество элементов i|/=и ра- 1 «(Л нее вычисленную оценку л;/ вероятности переходов по горизонтали, вычисляются оценки 2п\0 и л \lj}. Учитывая локальные изменения статистических характеристик на изображении, для их вычисления используется метод «сканирующего окна» [4]. Размер окна задается исходя из заданной точности оценки и минимальных вычислительных затрат. При обработке первых строк и столбцов изображения размер окна будет увеличиваться по мере увеличения номера обрабатываемого элемента, пока не достигнет заданных размеров. Размер окна на каждом этапе обработки определяется из пары соотношений: т = 1,(2/-1), г = 2,(М-1)/2; п = 1,(27-!), j = 2,{N-\)/2, (5) где i, j - координаты обрабатываемого изображения, M,N - заданные ширина и высота сканирующего окна. Усредняя оценки вероятностей переходов в пределах окна, вычисляется оценка „ _(/) средней вероятности перехода для элемента, соответствующего центральному элементу окна: яш 1 тхп т п Г=1 к=1 (6) На каждом последующем шаге окно перемещается на один элемент (слева-направо, сверху-вниз) и вычисляется оценка средней вероятности перехода для каждого элемента РДИ. Для выделения областей с разной тексту- „ ~(Л рой вычисленная оценка я),./ сравнивается с порогом. Пороговым значением h между двумя разными текстурными областями является оценка йда , которая для текстурных областей с разными вероятностными характеристиками может быть выбрана на основе анализа гистограмм изображения. В случае когда априорно известно, что изображение содержит области двух текстур с разными статистическими характеристиками, для их выделения будет достаточно двух меток (0 и 1). Всем элементам, для которых оценка nf) превышает порог, присваивается метка 1, остальным элементам - метка 0. В случае когда ЦПИ содержит несколько текстур, необходимо каждой из них присвоить свою метку. При этом используется несколько пороговых значений, соответствующих нескольким текстурам на изображении. Если ЦПИ содержит несколько текстур, необходимо каждой из них присвоить свою метку. При этом используется несколько пороговых значений, соответствующих нескольким текстурам на изображении. На рис. 2 приведен пример гистограммы текстурного признака изображения, содержащего области четырех текстур с матрицами вероятно сти переходов равными: *n = 2n = 0,05 0,95 *П = 2П = 1П = 2П = 0,75 0,25 0,25 0,75 0,5 0,5 0,5 0,5 , ‘n = 2n = 0,6 0,4 0,4 0,6 и Пикселей, тыс. 20 15 10 1 f "Г 1 1 1 т 1 Г 1 і і у — і : і _L_ 1 _j_і і і 7” 1 і 1 1 ! 1/ 1 \---h- \ 1 і і і і 1 7 Н - 0,546 [’ * = °Д \ 1 7 т \ і і і Т 1 і і L [ 4 11 ИІУ! -1—і--1— L_l Т і \| -L 1 н = о|?зб j 0,5 0,6 0,7 0,8 0,9 Рис. 2. Гистограмма текстурного признака Число пиков на гистограмме текстурного признака изображения соответствует количеству текстур на текстурном изображении. Пороговое значение выбирается как минимальное значение между двумя соседними пиками гистограммы. Выбрав порог по тестовому набору изображений с разными статистическими характеристиками, данное пороговое значение можно использовать для сегментации других изображений, содержащих области с такими же статистическими характеристиками. Результаты моделирования Моделирование проводилось на реальных и искусственных изображениях, имитирующих реальные аэрокосмические изображения и содержащее естественные текстурные области в среде Matlab. «Инфокоммуникационные технологии» Том 10, № 4, 2012 44 Карлушин К.А., Курбатова Е.Е., Медведева Е.В. В ч а) б) шш % ч. В) е) Щ Щ ж) з) и) к) Рис. 3. Результаты сегментации искусственного изображения: а) исходное изображение; б) разметка; в) РДИ; г) результат сегментации разработанным методом; д,ж,и) контуры, выделенные методами Канни, Собела, Лапласиана гауссиана, соответственно; е, з, к) сегментированные изображения, полученные методами Канни, Собела, Лапласиана гауссиана, соответственно Для получения оценки качества сегментированного изображения искусственное изображение формировалось по заданной разметке. На рис.3 представлены результаты сегментации искусственного аэрокосмического изображения (размер 1024*1024 пикселей) разработанным и известными контурными методами. Сегментация разработанным и известными методами осуществлялась по 6-му РДИ. Для закрашивания контурных областей (см. рис. 3 д, ж, з), полученных известными методами, применялся построчный алгоритм заполнения с затравкой [5]. Из приведенных на рис. 3 данных следует, что разработанный метод эффективно выделяет на изображении два сегмента, один из которых относится к области леса, а другой к области озера (см. рис. 3г). Область леса однородна с точки зрения текстуры, но при этом внутри этой области присутствует большое число перепадов яркости, которые воспринимаются методами, основанными на вычислении градиента как границы областей. Известные методы выделяют значительно большее количество сегментов. Таблица 1 Тип ЦПИ Разрабо танный метод Метод Канни Метод Собела Метод Лапла сиана гаус сиана Искус ственное 0,0079 0,7665 0,4811 0,5956 Реальное 0,0388 0,7268 0,4707 0,5969 Для оценки качества работы метода выделения текстурных областей и сравнения его с известными контурными методами вычислялось количество ошибочно сегментированных элементов (ESE), для определения которых сегментированное изображение сравнивалось с разметкой. Результаты вычисленного критерия ESE приведены в таблице 1. Из приведенных в таблице 1 данных следует, что ошибка сегментации разработанным методом более чем в 10 раз ниже по сравнению с известными методами, основанными на вычислении градиента. Таким образом, предложенный метод текстурной сегментации аэрокосмических изображений на основе марковских случайных полей при меньших вычислительных ресурсах эффективно выделяет протяженные области с однородными статистическими характеристиками. «Инфокоммуникационные технологии» Том 10, № 4, 2012

About the authors

K. A Karlushin

Email: karlushin-work@mail.ru

E. E Kurbatova

Email: kurbatovae@gmail.com

E. V Medvedeva

Email: EMedv@mail.ru

References

Supplementary files