Практика цифровой обработки оптических сигналов

1. Практика цифровой обработки оптических сигналов Учебное пособие (краткий курс лекций) Лычагов В.В., Рябухо В.П. ГОУ ВПО «Саратовский государственный университет им. Н.Г. Чернышевского» Кафедра оптики и биофотоники Саратов, 2010

2. Постановка задачи

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

4. Интерферометрия (лазерная, низкокогерентная, спекл-, полного поля, микро-, голографическая и т.д.) • Измерения: Толщин; Профиля поверхности; Деформаций; Смещений; Вибраций; Параметров движения; Оптических напряжений; Скрытых дефектов; … Цифровая фильтрация, морфологический анализ, корреляционный анализ, частотно-временные и частотно-пространственные преобразования…

5. Оптическая когерентная томография (ОКТ) • Восстановление внутренней оптической 2-D и 3-D структуры объекта по одномерным сечениям (формирование искусственных изображений) Цифровая фильтрация, визуализация, частотно-временные и пространственно-временные преобразования…

6. Фурье-спектроскопия • Расчет спектра излучения источника по его функции автокорреляции посредством Фурье-преобразования Фурье- преобразование Фурье-преобразования, цифровая фильтрация…

7. Спектральный метод ОКТ Обратное преобразование Фурье Фурье-преобразования, цифровая фильтрация, корреляционный анализ, морфологический анализ…

8. Аналого-цифровыеи цифро-аналоговые преобразования

9. Преобразование сигналов: аналоговый – цифровойДискретизация по времени: теорема отсчетов • T – период дискретизации; • Fd=1/T – частота дискретизации; • Критерий (Теорема) Найквиста (Шеннона, Котельникова): Частота дискретизации сигнала должна быть как минимум в два раза больше максимальной частоты, содержащейся в сигнале Fd ≥ 2*Fmax Fn = Fd / 2 Fmax ≤ Fn

10. Преобразование сигналов: аналоговый – цифровойДискретизация по времени: периодичность спектра

11. Преобразование сигналов: аналоговый – цифровойДискретизация по времени: наложение спектров • Частота дискретизации достаточная, чтобы адекватно разрешить все частоты сигнала • Частота дискретизации мала – происходит наложение спектров

12. Преобразование сигналов: аналоговый – цифровойКвантование по уровню: ошибка квантования • x(n) = x(t=nT) – дискретный сигнал • xq(n) – цифровой сигнал • e(n) = x(n) - xq(n)– шум квантования • xq(n) = x(n) + e(n) • T –период дискретизации • Q – шаг квантования

13. Преобразование сигналов: аналоговый – цифровойКвантование по уровню: метод округления • Ошибка метода округления -Q/2 ≤ e(n) ≤ Q/2 Среднее значение ошибки округления равно 0

14. Преобразование сигналов: аналоговый – цифровойКвантование по уровню: метод усечения • Ошибка метода усечения 0 ≤ e(n) ≤ Q Среднее значение ошибки усечения равно Q/2 Диапазон входных напряжений АЦП +/-5В Разрядность АЦП 12 бит Шаг квантования 2.4мВ

15. Преобразование сигналов: цифровой – аналоговый Интерполятор нулевого порядка

16. Ряд ФурьеПреобразование ФурьеДискретное преобразование ФурьеБыстрое преобразование ФурьеСпектральный анализОценка спектра

17. Преобразование Фурье

18. Разложение в ряд Фурье Коэффициенты ряда Фурье t – здесь - время, но может быть чем угодно w – циклическая частота первой гармоники T – период изменения сигнала T = 2p/w k = 1, 2, 3, 4, …

19. Разложение в ряд Фурьечетных и нечетных функций cos(- kwt) = cos(kwt) Ряд Фурье четной функции содержит только косинусоидальные члены: sin(- kwt) = - sin(kwt) Ряд Фурье нечетной функции содержит только синусоидальные члены:

20. Экспоненциальное представление ряда Фурье связь с коэффициентами действительного ряда Фурье: амплитудный и фазовый спектр:

21. Спектр последовательности прямоугольных импульсов • действительная часть коэффициентов ряда Фурье – коэффициенты при косинусе • мнимая часть коэффициентов ряда Фурье – коэффициенты при синусе • амплитудный спектр сигнала • спектр мощности сигнала

22. Разложение непериодических сигналов Интегральное преобразование Фурье

23. Свойства преобразования Фурье • Линейность: • Сдвиг сигнала во времени:

24. Свойства преобразования Фурье • Подобие: • Теорема Парсеваля:

25. Дискретное преобразование Фурье k, i = 0, 1, 2,…, N-1 x1, x2, x3,…xN-1– отсчеты дискретного сигнала

26. Быстрое Преобразование Фурье k = 0, 1, 2,… N-1 x(i), i = 0, 1, 2,... N-1 x1(i) = x(2i), i = 0, 1, 2,... N/2 – 1 x2(i) = x(2i+1), i = 0, 1, 2,... N/2 – 1

27. Быстрое преобразование ФурьеВычисление 8-точечного БПФ X(0) = X1(0) + W0X2(0); X(4) = X1(0) – W0X2(0); … и так далее

28. Быстрое преобразование ФурьеВычисление 8-точечного БПФ операция инверсии битов:

29. Проблемы спектрального анализа Эффекты конечного размера реализации Составляющая сигнала на частоте f?не может быть представлена в спектре. Возможный вариант решения проблемы – использование дополняющих нулей. Исходная реализация содержит Nточек отсчитанных через T секунд. Дополненная реализация – N+N`точек через те же Tсекунд. • «гребешковое искажение»

30. Эффекты конечного размера реализации • Просачивание спектральных составляющих и размывание спектра.Выборка данных в течение некоторого времени эквивалентна умножению сигнала на функцию (временное окно)спектр которой (спектральное окно) равен

31. Проблемы спектрального анализа Случайный характер измеряемых величин • Каждой выборке или реализации случайного процесса соответствует выборочный спектр, так же имеющий случайный характер. • Для детерминированных сигналов выборочный спектр при увеличении времени измерения сходится кистинному спектру сигнала. • Для случайных сигналов выборочный спектр не сходится к какому-либо предельному значению. • Для случайных сигналов можно говорить только об оценке спектра • Оценка спектра, являясь статистической величиной, может быть охарактеризована смещением и дисперсией.

32. Немного статистики • Математическое ожидание величины x(i): • Дисперсия ряда x(i): • Если оценка спектра величина статистическая, то она должна иметь функцию распределения плотности вероятности • Смещение оценки спектра: – истинный спектр случайного процесса – статистическая оценка этого спектра

33. Метод периодограмм x(i), i = 0 ... N-1 –дискретный временной ряд • Величина называется периодограммой и представляет собой оценку спектра мощности случайного сигнала в дискретном представлении • Для больших N смещение периодограммы является незначительным и сходится к • Дисперсия периодограммы не равна 0 и не стремится к нему при увеличении N – коэффициенты дискретного преобразования Фурье

34. Модифицированные периодограммы метод Бартлетта метод Уэлча • дисперсия оценки по методу Уэлча меньше, но перекрывающиеся сегменты коррелируют между собой.

35. Взвешивание • узкое временное окно (широкое спектральное окно, малое М) – дисперсия оценки малая, но большое смещение; • широкое временное окно (узкое спектральное окно, большое М) – дисперсия увеличивается, но смещение оценки становится малым

36. Небольшой пример реализация шумового сигнала... ...он же, растянутый во времени

37. Небольшой пример

38. Корреляция и свертка:функция корреляции,коэффициент корреляции,циклическая корреляция,линейная корреляция,быстрая корреляция,импульсная характеристика системы,теорема о свертке,обращение свертки...

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

40. Нормированная функция корреляцииКоэффициенты корреляции коэффициент взаимной корреляции функция автокорреляции коэффициент автокорреляции

41. Циклическая и линейная корреляции x={x1,x2,x3,x4} y={y1,y2,y3} Циклическая корреляция периодична с периодом, равным длине более короткой реализации Для получения линейной корреляции последовательности необходимо дополнить нулями: x={x1,x2,x3,x4,0,0} y={y1,y2,y3,0,0,0} Длина конечного вектора: N1+N2-1 N1 – длина первого вектора N2 – длина второго вектора

42. Быстрая корреляция Расчет быстрой корреляции методом БПФ Теорема о корреляции: Это выражение для циклической корреляции. Для получения линейной корреляции нужно использовать дополняющие нули. Расчет быстрой корреляции рекурсивным методом (новое значение) = (предыдущее значение) + +1/N(произведение двух новых членов)- -1/N(произведение первых двух членов)

43. Применение корреляции в обработке одномерных данных Расчет спектра мощности сигнала. Теорема Винера-Хинчина Выделение периодических составляющих из зашумленных сигналов

44. Свертка. Связь входного сигнала, выходного сигнала и импульсной характеристики системы y(0) = h(0)x(0); y(1) = h(1)x(0) + h(0)x(1); y(2) = h(2)x(0) + h(1)x(1) + h(0)x(2); … y(n) = h(n)x(0) + h(n-1)x(1) + … + h(0)x(n)

45. Графическая интерпретация операции свертки x(m) h(m) x(m) h(-m) x(m) h(t-m)

46. Свертка во временной и частотной областях Теорема о сверткеСвойства свертки • Коммутативность • Дистрибутивность • Ассоциативность

47. Обращение сверткиКруговая и линейная свертка Обращение свертки Идентификация системы N M-1 M N-1

48. Вычисление свертки и корреляции методом сегментации

49. Литература, которая может понадобиться: • Бендат Дж., Пирсол А. Применение корреляционного и спектрального анализа: Пер. с англ. – М.:Мир, 1983. – 312 с. • Макс Ж. Методы и техника обработки сигналов при физических измерениях: В 2-х томах. Пер. с франц. – М.:Мир, 1983. – 568 с. • Айфичер Э.С., Джервис Б.У. Цифровая обработка сигналов: практический подход, 2-е издание: Пер. с англ. – М.: Издательский дом «Вильямс», 2004. – 992 с. • Оппенгейм А.В., Шафер Р.В. Цифровая обработка сигналов: Пер. с англ. – М.: Связь, 1979. – 416 с. • Рабинер Л., Гоулд Б. Теория и применение цифровой обработки сигналов: Пер. с англ. – М.:Мир, 1978. – 849 с. • Гонсалес Р., Вудс Р. Цифровая обработка изображений. – М.: Техносфера, 2005. – 1072 с. • Прэтт Э. Цифровая обработка изображений: Пер. с англ. – М.:Мир, 1982. – 312 с.