- •230100 – Информатика и вычислительная техника
- •Часть 2. Алгоритмы цифровой обработки изображений 58
- •Часть 3. Методические указания для практических и лабораторных работ 109
- •Введение
- •Часть 1. Программирование графического адаптера nvidia с использованием технологии cuda
- •Краткая история развития графических адаптеров
- •Архитектура nvidia cuda
- •Архитектура gpu-приложения
- •Основные различия между cpu и cuda-устройством
- •Различие возможностей между поколениями cuda-устройств
- •Инструментарий cuda Toolkit
- •Установка
- •Проверка правильности установки
- •Компиляция примеров
- •Прикладные программные интерфейсы cuda
- •Особенности компиляции
- •Опции компилятора nvcc
- •Профайлер
- •Отладка программ в режиме эмуляции устройства
- •Среда разработки nvidia Parallel Nsight
- •Подключаемые библиотеки
- •Библиотека cublas
- •Библиотека cufft
- •Часть 2. Алгоритмы цифровой обработки изображений
- •Видоизменение гистограммы изображения
- •Основы пространственной фильтрации
- •Сглаживающие пространственные фильтры
- •Сглаживание изображений локальным усреднением
- •Сглаживание по обратному градиенту
- •Сглаживающий сигма-фильтр
- •Сглаживание по наиболее однородным областям
- •Медианный фильтр
- •Описание метода медианной фильтрации
- •Одномерный медианный фильтр
- •Зависимость качества работы медианного фильтра от размера апертуры
- •Двумерный медианный фильтр
- •Пространственные фильтры выделения границ областей
- •Линейные методы выделения границ областей. Лапласиан
- •Нелинейные методы выделения границ областей
- •Оператор Превита
- •Частотные фильтры
- •Сглаживающий частотный фильтр
- •Частотный фильтр повышения резкости
- •Гомоморфная фильтрация
- •Морфологическая фильтрация
- •Операции математической морфологии
- •Дилатация
- •Алгебраические свойства дилатации и эрозии
- •Применение морфологической дилатации
- •Применение морфологической эрозии
- •Размыкание и замыкание
- •Применение морфологических операций размыкания и замыкания
- •Основные морфологические алгоритмы
- •Выделение границ.
- •Заполнение областей.
- •Выпуклая оболочка.
- •Утончение.
- •Утолщение.
- •Построение остова.
- •Часть 3. Методические указания для практических и лабораторных работ
- •Практическая работа №1. Введение в технологию cuda
- •Цель работы
- •Порядок выполнения работы
- •Теоретическая часть
- •Архитектура графических адаптеров nvidia.
- •Архитектура чипа g80.
- •Архитектура чипа Fermi.
- •Вычислительные способности графического процессора.
- •Среда разработки cuda.
- •Язык cuda с.
- •Модель программирования cuda.
- •Теоретическая часть Модель памяти cuda
- •Особенности глобальной памяти.
- •Особенности разделяемой памяти.
- •Теоретическая часть Программный пакет «Sarcis».
- •Команды меню.
- •Панель инструментов.
- •Механизм подключения динамических библиотек.
- •Создание внешней библиотеки
- •Пример листинга программы динамической библиотеки.
- •Медианная фильтрация
- •Варианты индивидуальных заданий
- •Выбор по дискретному критерию.
- •Выбор по среднему контрасту.
- •Выбор по среднему градиенту.
- •Выбор на основе простых статистик.
- •Параметрическая аппроксимация.
- •Моментный подход.
- •Варианты индивидуальных заданий
- •Теоретическая часть
- •Инициализация устройства.
- •Выделение памяти на gpu
- •Пересылка данных между cpu и gpu.
- •Вызов функции ядра.
- •Методические рекомендации
- •Цель работы
- •Варианты индивидуальных заданий
- •Методические рекомендации
- •Защита лабораторной работы
- •Лабораторная работа № 3. Работа с разделяемой памятью в технологии cuda
- •Порядок выполнения работы
- •Методические рекомендации
- •Методические рекомендации
- •Особенности текстурной памяти.
- •Особенности линейной текстурной памяти.
- •Защита лабораторной работы
- •Правила оформления отчета
- •Оформление титульного листа
- •Библиографический список
Выбор по дискретному критерию.
Исходные данные:
гистограмма яркостей
некоторый порог
,
разделяющий точки изображения
на два класса (фон и объект соответственно):
, (3.7)
В дискриминантном методе оптимальным критерием разделимости классов (дискриминантный критерий) является среднеквадратичное расстояние между классами
, (3.8)
где
и
– вероятности принадлежности точки
изображения фону (
)
или объект (
):
,
, (3.9)
а
и
– средние уровни яркости объекта и фона
соответственно:
,
. (3.10)
Оптимальным
является такой порог
,
при котором выполняется условие:
. (3.11)
Очевидно, что
функция имеет равный нулю минимум, когда
или
,
т.е. когда все точки изображения относятся
либо только к объекту, либо только к
фону. В диапазоне
функция принимает положительные
значения, что обеспечивает существование
максимума.
Выбор по среднему контрасту.
В этом алгоритме оптимальным порогом для сегментации изображения считается порог, который выделяет больше высококонтрастных и меньше низко-контрастных перепадов яркости, чем любой другой порог. Количественным выражением критерия является средний контраст, выделяемый данным порогом. Чем больше порог выделяет низко-контрастных перепадов, тем меньше средний контраст. Наоборот, чем больше выделяется высококонтрастных перепадов, тем больше значение среднего контраста. Порог, соответствующий максимальному среднему контрасту, является оптимальным.
Рассмотрим перепад
яркости между двумя смежными точками
и
,
имеющими яркости
,
соответственно. Перепад яркости между
двумя этими точками будет выделяться
порогом
тогда и только тогда, когда
.
Количество перепадов, выделяемых порогом
,
равно
, (3.12)
где
Контраст между
парой точек
,
зададим выражением
. (3.13)
Тогда полный контраст, соответствующий порогу , равен
, (3.14)
где
.
Средний контраст
всех перепадов, обнаруживаемых порогом
,
равен
(равен нулю, если
).
То значение
,
при котором функция
имеет максимальное значение, принимается
за оптимальный порог.
Выбор по среднему градиенту.
Для каждого уровня
яркости
вычисляется среднее значение модуля
градиента:
, (3.15)
где
– значение функции среднего градиента
яркости
;
– число точек
-й
яркости в кадре;
– значение градиента
в
-й
точке яркости
с координатами
.
Расчет градиента производится с помощью оператора Лапласа:
.
В качестве
оптимального порога выбирается уровень
,
для которого
максимально. Поскольку этот уровень
соответствует точкам с наибольшим
градиентом, он должен находиться на
границе между объектом и фоном.
Используя
дополнительную информацию, содержащуюся
в градиентном перепаде изображения и
функции
,
можно производить двухпороговую
сегментацию. В качестве второго порога
выбирается яркость
.
Сегментация производится согласно правилу
.
,
,
.
Выбор на основе простых статистик.
Значение глобального
порога дает отношение двух статик
,
:
, (3.16)
которые вычисляются следующими образом.
Сначала изображение
дифференцируется:
. (3.17)
Сумма абсолютных значений яркости градиента изображения равна
. (3.18)
Другая статистика – сумма произведений уровней яркости на максимальную из производных в точке – вычисляется как
. (3.19)
Эти формулы
справедливы при выполнении ряда условий,
среди которых равномерная яркость
объекта и фона, отсутствие шума на
изображении. Такие ограничения свойств
изображения крайне нереалистичны.
Прежде всего это касается наличия шума.
При его наличии предлагается оценить
параметры распределения шума с целью
выбора порога, удаляющего точки, не
лежащие на границе объект – фон.
Предположив закон распределения шума
нормальным с нулевым средним и дисперсией
,
легко оценить параметры распределения
градиента неграничных точек. Оно также
будет нормальным с нулевым средним и
дисперсией
.
Для удаления 99% неграничных точек
необходимо порог по производной
выбрать
равным
.
Очевидно, что в таком случае превышение яркости объекта над фоном должно быть не меньше этой величины.
Для оценки
среднеквадратического отношения
воспользуемся накапливаемой статистикой
. (3.20)
Таким образом,
процедура определения глобального
порога для сегментации становится
двухэтапной. Сначала по всем точкам
изображения вычисляются
и далее оценка среднеквадратического
отклонения
.
После выбора порога
вновь
просматривается
и в статистики
,
,
включаются только точки со значениями
производной, превышающими
.
Отношение полученных величин дает
несмещенное значение глобального порога
. (3.21)
