А.С.Крылов Лаборатория математических методов обработки изображений imaging.cs.msu.ru

1. Математические методы решения биометрических задачЛекция 2. Идентификация человека по радужной оболочке глаза А.С.Крылов Лаборатория математических методов обработки изображений http://imaging.cs.msu.ru Кафедра математической физики

2. Биометрический погранконтроль: ShenZhen – HongKong • 400,000 пересечений границы каждый день • Два сценария: Пассажиры & Водители машин • Функционируют более 100 автоматов • Три режима: Распознавание лица (инфра-красный; видимый свет) & отпечатки пальцев • 1,800,000 записанных персон. • Скорость проверки: 6 сек / переход • 80,000,000 переходов с июня 2005 г по июнь 2007.

3. Сканирование и распознавание: типы распознавания

4. Уникальные структуры радужки обусловлены радиальной трабекулярной сетью (trabecular meshwork),состоящей из:углублений (крипты, лакуны), гребенчатых стяжек, борозд, колец, морщин, веснушек, короны, иногда пятнышек, сосудов и других черт Радужка: внутренний орган глаза Радужка находится за роговицей и водянистой влагой передней камеры

5. История • XIX век - Альфонс Бертильон -французский антрополог, криминалист, изучал в качестве идентификатора радужку (>цвет, <структуру) • 1936г. - Франк Бурш -известный американский глазной хирург, впервые предложил концепцию использования сложных структур радужной оболочки глаза для распознавания человека как по отпечаткам пальцев • 1978г. - Роберт Хилл -разработал первую глазную систему, запатентовал метод распознавания человека по сетчатке глаза • 1987г. - Леонард Флом и Аран Сафир -два офтальмолога, запатентовали концепцию Франка Бурше, не имея собственных разработок • 1989г. - Доктор Джон Даугман -известный ученый, Леонард Флом и Аран Сафир попросили его разработать математическую теорию и алгоритмы распознавания человека по радужной оболочке, предоставив образцы фотоизображений • 1990, 1994гг. - Доктор Джон Даугман -впервые разработал (1990г.) и запатентовал (1994г.) практический метод кодирования структур радужки.

6. Радужка: формирование, устойчивость 3-й месяц внутриутробного развития- начинает формироваться 8-й месяц - практически сформирована структура Случайность - формируется случайно, даже у однояйцовых близнецов Гены - не влияют на структуры Устойчивость после 1-го года жизни - радужка окончательно сформирована и практически не меняется вплоть до самой смерти, если нет травм или патологий глаза

8. Радужка: наиболее удобный идентификатор Изолирована и защищена от вн. среды Нельзя изменить, не нарушив зрение Реакция на свет и пульсация зрачка используется для защиты от подделок Возможен ненавязчивый, бесконтактный и скрытный метод получения изображений Высокая плотность уникальных структур-3,2 бита/кв. мм. или около 250 независимых характеристик (у других методов - только ~50),30% параметров достаточно, чтобы принять решение о совпадении с вероятностью не более 10E-6

9. Сканирование и распознавание: история с National Geographic Daugman, 2005 12 лет, 1984г. 30 лет, 2002г.

10. Сенсоры 2005

11. Видрадужки Средние размеры по горизонтали: по вертикали: зависит от возраста, состояния здоровья, освещения и др. быстро изменяется форма может достаточно сильно отличаться от круга центр зрачка, как правило, смещен относительно центра радужки по направлению к кончику носа

12. Вариации размеров зрачка Изображения глаза одного человека, полученные на одной установке с интервалом в несколько минут, но при разных условиях освещения

13. Метод оценки качества изображений глаз При получении последовательности изображений, к сожалению, не все изображения могут быть высокого качества. Мы рассмотрим три различного типа изображений. (а) (б) (в) хорошее дефокусированное с ресницами

14. Метод оценки качества изображений глаз преобразование Фурье двух локальных квадратных частей локализация зрачка вычисление значения дескриптора качества выход выведение результата

15. Общая схема задачи идентификации по радужной оболочке глаза Исходное изображение Предобработка изображения Выделение радужной оболочки из изображения Нормализация радужной оболочки Параметризация данных радужной оболочки методом преобразования Эрмита Сравнение изображения с изображениями из базы данных Составление бинарных кодов изображения Вывод о том, что обладатель данной радужной оболочки не может быть найден корректно Идентификация обладателя данной радужной оболочки

16. Локализация зрачка: поиск центра зрачка и радужки

17. Локализация зрачка: поиск радиуса зрачка Применяя метод выделения контуров (Edge Detection), мы четко выделяем границу зрачка.

18. Выравнивание освещения

20. Внешняя граница При поиске внешней границы интегрирование ведется не по всему круговому контору, а только по его части, соответствующей боковым частям радужной оболочки , .

21. Преобразование Эрмита Нормализация радужной оболочки Вычисляются знаки сверток

22. Функции Эрмита где – полиномы Эрмита: Функции Эрмита: • Образуют полную и ортонормированную в систему функций • Локализованы как в пространственном, так и в частотном пространствах • Имеют среднее значение в нуле

23. Преобразование Эрмита - функция интенсивности изображения - m,n-я двумерная функция Эрмита Преобразование Эрмита определяется значением сверток Далее используются знаки этих сверток

24. Поиск оптимальных номеров функций Эрмита База данных CASIA-IrisV3 (2655 изображений) Метрика сравнения изображений – расстояние Хэмминга Ham(L) Обозначим Наилучший результат - используя значение (1, 0) + (2, 0) .

25. Сравнение изображений • Сортировка изображений из базы данных по значениям и • Циклический сдвиг нормализованного изображения на 3, 6, 9, 12, 15, 18, 21пикселей как вправо, так и влево - учитываются углы поворота глаза от ‑14º до 14º. • FAR = 0, FRR= 082%

26. Определение наложения нижнего века на радужную оболочку

27. Результаты

28. Введение шума в изображения • Исходное изображение

29. Введение шума в изображения • Равномерно распределенный шум. Порог = 15

30. Введение шума в изображения • Равномерно распределенный шум. Порог = 30

31. Размытие изображений • Исходное изображение

32. Размытие изображений • Размытие Гауссом с sigma = 1

33. Проекционный метод Эрмита Проекционный метод Эрмита основан на свойствах функций Эрмита, являющихся собственными функциями преобразования Фурье.Аппроксимация сигнала рядом этих функций позволяет осуществлять локализацию информации одновременно в пространственной и частотной областях.

34. Функции Эрмита: A) Б) Образуют полную ортонормированную систему функций в . В)

35. Общая форма разложения: Быстрый метод: где и - нули полинома Эрмита

36. Иерархический проекционный метод Эрмита Первый уровень иерархии: для каждой точки оси OX среднее значение интенсивности вдоль оси OY раскладывается в ряд по функциям Эрмита Второй уровень иерархии: Четвертый уровень иерархии:

37. Сравнение • Метрика - сумма квадратов отклонения коэффициентов. • Сортировка изображений из базы данных по расстоянию до взятого изображения происходит на всех уровнях иерархии • Циклический сдвиг нормализованного изображения на 3, 6, 9, 12, 15пикселей как вправо, так и влево - учитываются углы поворота глаза от ‑10º до 10º.

38. Полярное преобразование Эрмита • Полярные коэффициенты Эрмита - декартовые коэффициенты Эрмита

39. Полярное преобразование Эрмита • При повороте исходного изображения на угол полярные коэффициенты пересчитываются по формуле: • Для вычисления полярных коэффициентов Эрмита для повернутого изображения на заданный угол , достаточно умножить полярные коэффициенты исходного неповернутого изображения на матрицу поворота.

40. Полярный метод Эрмита • Нормализация изображения глаза • Считаются полярные коэффициенты Эрмита для входного нормализованного изображения и для этого же изображения, повернутого на углы • Сравнение с полярными коэффициентами всех других изображений глаз. • Сортировка изображений из базы данных по расстоянию до взятого изображения.

41. Результаты работы алгоритма определения угла поворота глаза (в) ближайшее к нему изображение из базы данных. Угол поворота равен 8º (а) исходное изображение (б) нормализованное изображение (г) изображение (б), повернутое на 8º.

42. Изображение 1

43. Изображение 2

44. Вспомогательный метод определения обладателя радужной оболочки глаза • В изображении радужной оболочки с нормализованными границами берется прямоугольная область: • В этой области используется проекционный метод Эрмита.

45. Алгоритм сравнения результатов методов параметризации

46. Результаты работы проекционного метода БЕЗ ШУМА • FAR = 0.1% • FRR = 9% • С ШУМОМ • Уровень шума 25 - FRR = 20% • Уровень шума 30 - FRR = 25%

47. Результаты работы алгоритма на зашумленных изображениях • Равномерно распределенный шум различных уровней

48. СПАСИБО ЗА ВНИМАНИЕ