Россия
Россия
ВАК 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
ВАК 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
ВАК 05.13.10 Управление в социальных и экономических системах
ВАК 05.13.18 Математическое моделирование, численные методы и комплексы программ
ВАК 05.13.19 Методы и системы защиты информации, информационная безопасность
УДК 65 Управление предприятиями. Организация производства, торговли и транспорта
ГРНТИ 20.01 Общие вопросы информатики
ГРНТИ 28.01 Общие вопросы кибернетики
ГРНТИ 49.01 Общие вопросы связи
ГРНТИ 50.01 Общие вопросы автоматики и вычислительной техники
ГРНТИ 82.01 Общие вопросы организации и управления
Рассматривается опыт применения нелинейного сглаживающего фильтра для обработки и анализа изображений. Дано описание модели сглаживающего фильтра, основанного на обобщенном методе наименьших модулей. Фильтр, построенный на основе предложенной модели, эффективно устраняет шум на перепаде. Наряду с подавлением шума в контрастных изображениях данный метод можно использовать для решения задач машинного зрения, медицинской диагностики и в других областях. Установлено, что нелинейная фильтрация на основе обобщенного метода наименьших модулей позволяет решать такие задачи, как уточнение границ контрастных объектов и сегментация изображения. Показана возможность восстановления границ в изображениях, у которых по каким-либо причинам были размыты контрастные границы. В качестве примера использовано изображение рентгеновского снимка кисти животного с расфокусировкой. После фильтрации контрастная граница была восстановлена в том месте, где она и находилась изначально. При обработке снимка флюорографии фильтр позволил удалить из изображения различные артефакты и повысить контраст. Удаление артефактов вместе с уточнением границ контрастных объектов улучшает общую «читаемость» снимка флюорографии, а также позволяет увидеть не различимые ранее детали на изображении. Приведены примеры применения фильтра в задаче кластеризации при помощи алгоритма k-средних. Ввиду недостатка этого алгоритма применение его непосредственно к изображению не дает приемлемого результата. Однако после обработки исходного изображения нелинейным фильтром применение алгоритма k-средних дает требуемый результат.
контрастные изображения, контрастные границы, обработка, подавление шума, сегментация, обобщенный метод наименьших модулей
Введение
В настоящее время цифровые изображения стали повседневным явлением, поэтому к цифровой фото- и видеообработке проявляют все больший интерес. На данный момент разработано множество способов как получения, так и обработки изображений. Вследствие того, что цифровая техника не является идеальной, в некоторых случаях приходится прибегать к обработке изображений для борьбы с шумом, в целях улучшения качества, модификации изображений для последующего анализа и т. д.
Подавляющее большинство изображений сегодня получают при помощи цифровой аппаратуры. Но наряду с преимуществами перед аналоговой аппаратурой цифровая техника также имеет недостатки и, как любое измерительное устройство, имеет погрешности измерений. В результате сформированное изображение состоит как из полезного изображения, так и из цифрового шума в виде наложенной маски из пикселей случайного цвета и яркости [1].
Задача повышения качества цифрового изображения является одной из важнейших в цифровой обработке изображений. Она решается в основном с помощью фильтров изображений [1–4]. Однако, несмотря на многообразие алгоритмов цифровой фильтрации, они, как правило, носят эвристический характер. Особенно это заметно при подавлении шума в контрастных изображениях. Основной причиной является недостаточное исследование свойств цифрового шума в контрастных изображениях.
Многие авторы [5, 6] без проведения исследований и рассмотрения особенностей изображений указывают, что распределение цифрового шума подчиняется гауссовскому закону распределения. В результате предлагаются в основном линейные алгоритмы подавления шума. Однако, как отмечено в [7], из-за математической простоты, которая характерна для работы с моделями гауссовского шума, линейные алгоритмы используются даже в тех ситуациях, когда их применение оправдано, в лучшем случае, лишь частично. Такая неоднозначность высказываний относительно статистических свойств аддитивного шума говорит о недостаточной изученности данного вопроса.
В [8] показано, что аддитивный шум в цифровых контрастных изображениях имеет нелинейный характер. На границе перепадов дисперсия шума имеет разную величину, а на границе минимальной яркости шум имеет негауссовское распределение. Это происходит из-за нелинейного характера формирования шума. Использование линейных цифровых фильтров для контрастных изображений неэффективно, т. к. приводит к размыванию контуров деталей контрастных изображений, что снижает резкость изображения и, следовательно, ухудшает его качество. Применяемые нелинейные алгоритмы подавления шума на основе медианной фильтрации, предложенной Джоном Тьюки [9], в случае контрастных изображений имеют недостатки, т. к. медианные фильтры являются чувствительными к резким перепадам яркости в изображении. Различные модификации медианной фильтрации [10–12] в ряде случаев достигают выигрыша по сравнению с медианной фильтрацией, но они недостаточно формализованы и обычно требуют наличия дополнительной априорной информации.
В [13] была предложена модель нелинейного фильтра на основе обобщенного метода наименьших модулей (ОМНМ) для фильтрации цифровых контрастных изображений. Апробация предложенной модели на простых контрастных изображениях показала ее высокую эффективность – подавление шума происходит с минимальным размытием резких перепадов яркости на границах.
Однако на подавлении шума возможности подложенной модели не заканчиваются и, предположительно, ее можно использовать в обработке изображений для последующего их анализа. Задачи анализа изображений встречаются в разных областях науки и производства и, как часто бывает, требуют предварительной обработки изображений перед применением методов анализа информации, содержащейся в них.
Цель статьи – показать применения ОМНМ в задачах подавления шума, а также обработки и анализа изображений для различных областей знаний.
Описание модели
Функция потерь сглаживающего преобразования была выбрана на основе результатов, описанных в [8], где в качестве функции потерь была использована функция арктангенса
(1)
Однако непосредственное использование (1) не учитывает особенностей изображения и перепада, таких как дисперсия шума и степень контрастности перепада. Так как достижение функцией потерь (1) своего предельного уровня может происходить слишком быстро либо недостаточно быстро, и в том и в другом случае это будет приводить к низкой эффективности оценок и даже к их смещению. Поэтому вместо (1) будем использовать функцию потерь вида
(2)
где α > 0, β > 0 – параметры, подбираемые на основе характеристик изображения и перепада, таких как σ – среднеквадратическое отклонение (уровень) шума в выбранной области изображения; A – нижняя граница уровня яркости перепада; B – верхняя граница уровня яркости перепада; h = A – B – разница между верхней и нижней границами яркости перепада (степень контрастности перепада).
Модель предполагает варьирование двумя параметрами функции потерь (2) в зависимости от степени контрастности и уровня шума. Чтобы учесть негауссовость распределения шума вблизи нижней границы яркости следует использовать свойство вогнутости функции потерь, поэтому задаем параметр α (0; 1], в случае же гауссовости распределения шума использовать выпукло-вогнутую функцию при α
(1; 2]. В случаях, когда требуется обработка изображения без шумовой составляющей, параметр α будет фиксирован. При изменении уровня шума следует менять параметр β, увеличивая его при возрастании уровня шума.
Для правильного выбора параметров функции потерь (2) необходимо подобрать такие функции ,
,
чтобы результат сглаживания ОМНМ-фильтром с функцией потерь (2) был удовлетворительным при любых параметрах изображения [13].
В ходе численного эксперимента, который проводился нами с помощью метода статистических испытаний Монте-Карло, были получены следующие зависимости для a и b:
Соответственно, при отсутствии шума параметр α = 0,7.
С целью упрощения рассматриваем изображение в градациях серого либо яркостную составляющую цветного изображения модели HSV, что даст 256 возможных значений яркости пикселя, где «0» – это черный цвет, а «255» – белый. В этом случае процедура фильтрации описывается формулой
где – фактические и сглаженные с помощью ОМНМ-фильтра значения яркостей пикселей изображений; i – номер пикселя; L – пул точек, входящих в апертуру фильтра, в зависимости от ее вида и размера.
Для оценки качества фильтрации используем среднее линейное отклонение между идеальным изображением и сглаженным, вычисляемое как
,
где fi – идеальные (без шума) значения яркостей пикселей изображений; N – число пикселей в изображении, обработанных при сглаживании.
Применение нелинейной фильтрации в задачах обработки и анализа изображений
При фильтрации реальных контрастных зашумленных изображений (рис. 1) фильтр позволяет сохранить резкую границу без ее размытия.
Рис. 1. Реальное зашумленное изображение здания (а);
увеличенная область, выделенная квадратной рамкой (б)
Чтобы наглядно оценить эффективность сглаживания, отфильтруем реальное зашумленное изображение с помощью различных фильтров: усредняющий фильтр, медианный фильтр, фильтр на основе функции потерь (2). На рис. 2 приведены отфильтрованные различными фильтрами области при увеличении.
Рис. 2. Отфильтрованные изображения выделенной области с помощью линейного усреднения (а);
медианы (б); фильтра на основе функции потерь (2) (в)
Фильтрация зашумленных контрастных изображений на основе ОМНМ позволяет удалять шум из изображения, сохраняя при этом контрастные границы. Также в ходе экспериментов было установлено, что фильтр на основе ОМНМ позволяет восстанавливать контрастные границы в изображениях, в которых по каким-то причинам они были размыты.
Для испытания возьмем исходное незашумленное изображение здания и размоем контрастные границы. Далее проведем фильтрацию медианным и ОМНМ-фильтром. Результаты, в виде увеличенной области изображения, представлены на рис. 3.
Рис. 3. Увеличенные области отфильтрованных изображений: исходное изображение (а);
медианный фильтр (б); фильтр на основе функции потерь (2) (в)
Такое свойство фильтра можно использовать в задачах восстановления дефектных изображений, например рентгеновского снимка с расфокусировкой.
Качество рентгеновского снимка, с технической точки зрения, характеризуется оптической плотностью почернения, резкостью и контрастностью. Кроме того, рентгеновский снимок должен быть правильно выполнен. Правильно выполненным считается тот снимок, который при высоких технических качествах отображает подлинное состояние исследуемого объекта. Но в некоторых случаях по разным причинам получается снимок ненадлежащего качества. Это отсутствие резкости вследствие того, что сбился или неправильно настроен фокус, недоэкспонированный снимок от недостатка мощности, размазывание из-за движения объекта съемки и т. д. Особо важна резкость при съемке скелета, т. к. в этом случае важно четко видеть границу кости.
Некоторые подобные дефектные снимки бывает по разным причинам сложно оперативно переснять. Обработка такого оцифрованного снимка фильтром на основе ОМНМ позволит получить приемлемый результат без проведения повторного сеанса рентгеноскопии.
Возьмем рентгеновский снимок кисти животного в двух вариантах. В первом варианте объект был в фокусе, во втором варианте – не в фокусе. Обработаем изображение, на котором объект не попал в фокус, фильтром. На рис. 4 приведены участки различных изображений в увеличении.
Рис. 4. Рентгеновский снимок кисти: изображение в фокусе (а); изображение не в фокусе (б);
изображение не в фокусе после обработки фильтром (в)
Фильтр восстановил границу контрастного объекта именно на том месте, где она находилась изначально. Это позволит получать необходимые данные даже на некачественных снимках и производить более точные замеры в сравнении с размытым изображением, например, ширины кости.
Оцифрованные снимки флюорографии (ФОГ) в некоторых случаях могут содержать артефакты, мешающие анализировать визуальную информацию. На рис. 5 приведен снимок грудной клетки человека и одна из областей этого снимка при увеличении.
а б
Рис. 5. Исходное изображение снимка ФОГ (а); выделенная область в увеличении (б)
В таких случаях имеет смысл применить фильтр на основе ОМНМ. Используя апертуру типа крест или квадрат с минимальным размером можно избавиться от артефактов. В таком случае мелкие детали размером в половину апертуры могут быть утеряны, но то же случится и при обработке, например, медианой. Однако в отличие от медианного фильтра контраст на границах яркости сохранится. Пример обработки снимка ФОГ фильтром на основе ОМНМ показан на рис. 6.
Рис. 6. Изображение выделенной области снимка ФОГ после обработки
усредняющим фильтром (а); медианным фильтром (б); ОМНМ-фильтром (в)
В машинном зрении [14] часто решается задача сегментации изображения. Смысл этой задачи – выделить определенные объекты на изображении по каким-либо характеристикам, например по цвету или по форме. При сегментации решается задача кластеризации при помощи различных алгоритмов, например алгоритма k-средних [15]. Ввиду недостатков алгоритма k-средних, например высокой чувствительности к различным текстурам, применение его непосредственно к зашумленному изображению не дает приемлемого результата. В случае применения алгоритма k-средних к зашумленному изображению комнаты (рис. 7) обработанное изображение, несмотря на хорошее выделение в кластеры некоторых объектов, содержит большое количество кластеров малого размера.
Рис. 7. Исходное изображение (а); изображение после обработки алгоритмом k-средних (б)
В приведенном примере ярко выраженные текстуры и градиенты мешают алгоритму работать должным образом. На обработанном изображении много отдельных элементов, особенно на полу и стенах. Следовательно, необходимо упростить изображение для дальнейшего применения алгоритма k-средних. Применение фильтра на основе ОМНМ позволяет сильно ослабить контрастность текстур, при этом оставив перепад яркостей контрастным. В данном случае достаточно использовать квадратную апертуру с большим размером, затем к полученному изображению применить алгоритм k-средних. Результат изображен на рис. 8.
Рис. 8. Изображение, обработанное ОМНМ-фильтром (а);
после обработки его алгоритмом k-средних (б)
В результате большое количество кластеров на полу сократилось до одного большого кластера и подавление текстур и шума не привело к растеканию границ контрастных объектов, что могло бы привести к появлению дополнительных ненужных кластеров.
Для применения ОМНМ-фильтра к цветному изображению необходимо перевести изображение из цветового пространства RGB в HSV и обработать отдельно канал V, после чего провести обратное преобразование (HSV в RGB). Далее можно применять алгоритм k-средних.
Другим примером эффективного применения фильтра на основе ОМНМ может служить задача сегментации (выделения) отдельных клеток и клеточных структур [16] в медицине. В приведенном на рис. 9 примере шум, текстура и градиенты мешают алгоритму правильно выделить объекты.
Рис. 9. Исходное изображение (а) и после обработки алгоритмом k-средних (б)
На обработанном участке много отдельных элементов, особенно на темных участках и в месте, где есть межклеточная жидкость. Следовательно, необходимо упростить изображение.
Применение фильтра на основе ОМНМ удалило шум и сильно ослабило контрастность текстур, при этом оставив перепад яркостей контрастным. Далее к получившемуся изображению можно применять алгоритм k-средних. Результат представлен на рис. 10.
Рис. 10. Изображение, обработанное ОМНМ-фильтром (а);
после обработки алгоритмом k-средних (б)
Фильтрация на основе ОМНМ может быть также использована в следующем случае. Традиционный морфометрический подход к определению на изображении площади окрашенной клетки, лежащей на светлом фоне, состоит в выполнении достаточно сложной процедуры выделения границы клетки и подсчета количества пикселей, находящихся в пределах этой границы. Построение полутоновой гистограммы предполагает более простой подход к решению этой задачи [16]. Однако при зашумлении гистограмма «расплывается», что затрудняет последующий анализ. Используя фильтр на основе ОМНМ, можно удалить шум с сохранением контрастной границы, что позволит сократить разброс распределения яркостей на диаграмме (рис. 11).
Рис. 11. Гистограмма исходного изображения (а), зашумленного изображения (б),
изображения после обработки ОМНМ фильтром (в)
В результате гистограммы на рис. 11, а и в, сопоставимы, что позволит получить тот же результат, что и у исходного незашумленного изображения.
Заключение
Фильтрация зашумленных контрастных изображений на основе обобщенного метода наименьших модулей позволяет удалять шум из изображения, сохраняя при этом контрастные границы. Эффект достигается за счет выпукло-вогнутого вида функции потерь и адаптивного управления ее параметрами.
Наряду с подавлением шума нелинейный фильтр на основе ОМНМ может эффективно применяться в других задачах обработки и анализа изображений. Показана возможность восстановления границ в изображениях, у которых по каким-то причинам были размыты контрастные границы. Также установлено, что фильтрация на основе ОМНМ может эффективно применяться в задачах медицинской диагностики и машинного зрения.
1. Ярославский Л. П. Введение в цифровую обработку изображений. М.: Советское радио, 1979. 312 с.
2. Huang T. S., Eklund J.-O., Nussbaumer G. J., et al. Two-Dimensional Digital Signal Processing. Berlin - Heidelberg - New York: Springer-Verlag, 1981. V. 2: Transforms and Median Filters. 224 p.
3. Шеремет И. А., Рукин А. П., Лебедев В. Д. Обработка изображений с помощью целочисленных ортогональных преобразующих матриц // Цифровая обработка сигналов. 2014. № 4. С. 45-53.
4. Клочко В. К., Кузнецов В. П. Методы восстановления изображений и оценивания аппаратной функции по прореженной матрице наблюдений // Автометрия. 2016. Т. 52. № 6. С. 12-21.
5. Лукин А. С., Мамаев Н. В., Юрин Д. В. Локальная оценка уровня шума на цифровых изображениях // Цифровая обработка сигналов и ее применение: тр. 16-й Междунар. конф. 2014. Т. 2. С. 445-449.
6. Бухтояров С. С., Приоров А. Л., Апальков И. В., Хрящев В. В. Переключающийся медианный фильтр с блоком предварительного детектирования // Цифровая обработка сигналов. 2006. № 4. С. 2-8.
7. Gonzalez R. C., Woods R. E. Digital Image Processing. New Jersey: Prentice Hall, 2007. 976 p.
8. Сурин В. А., Тырсин А. Н. Исследование свойств цифрового шума в контрастных изображениях // Изв. высш. учеб. заведений. Физика. 2016. № 8/2 (59). С. 93-96.
9. Tukey J. W. Discussion emphasizing the connection between analysis of variance and spectrum analysis // Technometrics. 1961. V. 3. P. 191-219.
10. Abreu E., Lightstone M., Mitra S., Arakawa K. A new efficient approach for the removal of impulse noise from highly corrupted images // IEEE Trans, on Image Processing. 1996. V. 5. N. 6. P. 1012-1025.
11. Chan R., Ho C., Nikolova M. Convergence of Newton's Method for a Minimization Problem in Impulse Noise Removal // J. Comput. Math. 2004. V. 2. P. 168-177.
12. Schulte S., Nachtegael M., De Witte V., Van der Weken D., Kerre E. A fuzzy impulse noise detection and reduction method // IEEE Trans, on Image Processing. 2006. V. 15. N. 5. P. 1153-1162.
13. Сурин В. А., Тырсин А. Н. Модель нелинейного фильтра для цифровой обработки контрастных изображений // Автометрия. 2018. Т. 54. № 2. С. 54-62.
14. Gorban A. N., Zinovyev A. Y. Principal Graphs and Manifolds // Handbook of Research on Machine Learning Applications and Trends: Algorithms, Methods, and Techniques. IGI Global. Hershey, PA, USA, 2009. Part 2. P. 28-59.
15. MacQueen J. Some methods for classification and analysis of multivariate observations // In Proc. 5th Berkeley Symp. on Math. Statistics and Probability. Univ. of Calif. Press, 1967. V. 1. P. 281-297.
16. Богданов А. К., Проценко В. Д. Практические применения современных методов анализа изображений в медицине: учеб. пособие. М.: Изд-во РУДН, 2008. 77 с.
