Вейвлет-перетворення

Матеріал з testwiki
Перейти до навігації Перейти до пошуку
Приклад 2D Шаблон:Нп, яке використовується в JPEG 2000.

Шаблон:Сирий переклад В математиці, серії вейвлетів є подання Шаблон:Нп (дійсні або комплексні значення) функції за певною ортонормованої серії, породженої вейвлетом. В даний час, вейвлет-перетворення є одним з найпопулярніших тимчасових частот-перетворювань. У даній статті наводиться формальне та математичне визначення ортонормованого вейвлета, та інтегрального вейвлет-перетворення.

Визначення

Функція ψL2() називається ортонормованим вейвлетом, якщо її можна використати для визначення базису Гільберта, тобто повна ортонормована система, для Гільбертого простору L2() Шаблон:Нп. Базис Гільберта будується як сімейство функцій {ψjk:j,k} за допомогою двійкових перенесень і Шаблон:Нп ψ,

ψjk(x)=2j2ψ(2jxk) для цілих j, k. Якщо під стандартним передгільбертовим простором на L2(),
f,g=f(x)g(x)dx то це сімейство ортонормоване і його система теж ортонормована: :ψjk,ψlm=ψjk(x)ψlm(x)dx=δjlδkm

де δjl, є дельта Кронекера. Повнота виконується, якщо кожна функція hL2() може бути розширена в базисі як

h(x)=j,k=cjkψjk(x)

із збіжності рядів розуміється збіжність по нормі. Таке уявлення функції f відомо як серія вейвлета. Це означає, що ортонормований вейвлет Шаблон:Нп. Інтегральне вейвлет-перетворення є інтегральним перетворенням і визначається як : [Wψf](a,b)=1|a|ψ(xba)f(x)dx. Коефіцієнти вейвлетів cjk отримуються з формули: cjk=[Wψf](2j,k2j) Тут, a=2j називається бінарним розширенням або двійковим розширенням, і b=k2j це бінарна або двійкова позиція.

Принцип

Основна ідея вейвлет-перетворень є те, що перетворення має дозволити тільки зміни в продовження часу, але не форму. Це залежить від вибору відповідних базисних функцій, які дозволяють це. Зміни в продовження часу, як очікується, відповідають відповідній частоті аналізу функції базису. Виходячи з принципу невизначеності обробки сигналів:ΔtΔω12, де t представляє час і ω кутову частоту; (ω = 2πf, де f є тимчасова частота). Чим вище необхідний дозвіл в часі, тим нижче має бути дозвіл по частоті. Чим більше вибрано розширення вікна аналізу, тим більше величина .

Коли Δt велике,

  1. Поганий дозвіл за часом.
  2. Гарний дозвіл по частоті.
  3. Низька частота, великий коефіцієнт масштабування.

Коли Δt мала

  1. Гарний дозвіл за часом.
  2. Поганий дозвіл по частоті.
  3. Висока частота, малий коефіцієнт масштабування.

Іншими словами, базову функцію Ψ можна розглядати як імпульсний відгук системи, з якою функція х(t) була відфільтрована. Перетворений сигнал містить інформацію про час і частоту. Таким чином, вейвлет-перетворення містить інформацію, аналогічну до віконного Фур'є-перетворення, але з додатковими спеціальними властивостями вейвлетів, які з'являються в вирішенні під час більш високих частотах аналізу базисної функції. Різниця в дозволі часу на висхідних частотах для перетворення Фур'є і вейвлет-перетворення показано нижче.

Можна побачити, що вейвлет-перетворення кращий в дозволі часу високих частот, в той час як для повільно мінливих функцій, важливий дозвіл по частоті. Інший приклад: Аналіз трьох накладених один на одного синусоїдальних сигналів y(t)=sin(2πf0t)+sin(4πf0t)+sin(8πf0t), з STFT і вейвлет-перетворенням.

Вейвлет стиснення

Вейвлет стиснення є формою стиснення даних, яка добре підходить для стиснення зображень (іноді також стиснення відео і аудіо стиснення). Помітні реалізацій JPEG 2000, DjVu і нерухомих зображень, Шаблон:Нп і BBC's Шаблон:Нп. Мета полягає в тому, щоб зберігати дані зображення, зайнявши, при цьому як менше місця, наскільки це можливо в файлі. Вейвлет стиснення може бути або без втрат, або з втратами.[1]

За допомогою вейвлет-перетворення, методи вейвлет стиснення є достатніми для подання перехідних процесів, таких як ударні звуки в аудіо, або високочастотні компоненти в двовимірних зображень, наприклад, зображення зірок на нічному небі. Це означає, що перехідні елементи даних сигналу можуть бути представлені меншою кількістю інформації, ніж було б у разі використання якої-небудь іншої трансформації, наприклад широкого дискретного косинусного перетворення.

Дискретне вейвлет-перетворення було успішно застосовано для стиснення електрокардіографічних (ЕКГ) сигналів.[2] У даній роботі, висока кореляція між відповідними вейвлет-коефіцієнтами сигналів послідовних серцевих циклів і використовуються лінійне передбачення з використанням.

Вейвлет стиснення не підходить для всіх видів даних: перехідні характеристики сигналу означають гарне стиснення вейвлета, в той час як гладкі, періодичні сигнали краще стискаються за допомогою інших методів, зокрема, традиційна гармоніка стиснення (частотна область, за допомогою перетворень Фур'є і пов'язаними з ними).

Дивись Щоденник x264 Розробник: Проблеми з вейвлетами (2010) (Diary Of An x264 Developer: The problems with wavelets) для обговорення практичних питань існуючих методів з використанням вейвлетів для стиснення відео.

Метод

Перше вейвлет-перетворення було прикладним. Він виробляє стільки коефіцієнтів, скільки пікселів у зображенні (тобто, немає ніякого стиснення, оскільки це тільки перетворення). Ці коефіцієнти можна легше стиснути, тому що загальна матриця містить надлишкову інформацію, а значення їх — статистично залежні. Цей принцип називається Шаблон:Нп. Після цього коефіцієнти квантуються і квантовані значення ентропійно кодуються і / або кодується довжина послідовності. Декілька 1D і 2D застосувань вейвлет-стиснення використовують технологію, яка називається «вейвлетний слід» (Шаблон:Lang-en).[3][4]

Порівняння з перетворенням Фур'є та частотно-часового аналізу

Перетворення Представлення Вхідні дані
Перетворення Фур'є f^(ξ)=f(x)e2πixξdx ξ, частота
Частотно-часовий аналіз X(t,f) t, time; f, частота
Вейвлет-перетворення X(a,b)=1aΨ(tba)x(t)dt a, масштаб; b, час

Вейвлети мають деякі незначні переваги в порівнянні з Фур'є, наприклад, в скороченні обчислень при розгляді конкретних частот. Проте, вони рідко бувають більш чутливими, і, дійсно, спільний Шаблон:Нп — це математично ідентичний до віконного перетворення Фур'є з використанням функції вікна Гаусса.[5] Виняток є при пошуку сигналів відомої синусоїдальної форми (наприклад, серцебиття); в цьому випадку, використовуючи відповідний вейвлет можна перевершити стандартний аналіз STFT / Морлет.[6]

Інші практичні використання

Вейвлет-перетворення може дати нам з частотою сигналів час, пов'язаним з цими частотами, робить його дуже зручним для його застосування в різних областях. Так, наприклад, обробка сигналів прискорень для аналізу ходи,[7] для виявлення несправностей,[8] для дизайну електрокардіостимуляторів низької потужності, а також в надширокосмугових (Сніп) бездротового зв'язку.[9]

1.Дискретизація c-τ осі

Застосовують наступну дискретизацію частоти і часу:

cn=c0nτm=mTc0n

Призводить до форми вейвлетів, яка є дискретна формулою для базису вейвлета: Ψ(k,n,m)=1c0nΨ[kmc0nc0nT]=1c0nΨ[(kc0nm)T]

Такі дискретні вейвлети можуть бути використані для трансформації:

YDW(n,m)=1c0nk=0K1y(k)Ψ[(kc0nm)T]

2.Реалізація через БПФ (швидке перетворення Фур'є)

Видно з подання вейвлет-перетворення (як показано нижче) :YW(c,τ)=1cy(t)Ψ(tτc)dt

де с — коефіцієнт масштабування, τ являє собою фактор зсуву часу і, як вже згадувалися в цьому контексті, вейвлет-перетворення відповідає згортку функції у(t) і вейвлет-функцію]. Згортка може бути реалізована як множення в частотній області. При цьому наступний підхід результатів реалізації в:

  1. Фур'є-перетворення сигналу у(к) з FFT.
  2. Вибір дискретного коефіцієнта масштабування cn.
  3. Масштабування базису вейвлет-функції за допомогою цього чинника cn і результат швидкого перетворення Фур'є цієї функції.
  4. Множення з перетвореним сигналом YFFT з першого кроку.

№Зворотне перетворення продукту в результати тимчасової області в YW(c,τ) для різних дискретних значень τ і дискретне значення cn. №Повертаємося до другого кроку, до тих пір, поки всі значення дискретного масштабування для cn обробляються.

Є багато різних типів вейвлет-перетворення для конкретних цілей. Див. також: повний список вейвлета пов'язаних перетворень, але загальні з них перераховані нижче: Шаблон:Нп, вейвлет Хаара, Вейвлет, вейвлет Добеші, трикутний імпульс.

Див. також

  • Continuous wavelet transform
  • Шаблон:Нп
  • Complex wavelet transform
  • Шаблон:Нп
  • Multiresolution analysis
  • MrSID, the image format developed from original wavelet compression research at Los Alamos National Laboratory (LANL).
  • ECW (file format), a wavelet-based geospatial image format designed for speed and processing efficiency
  • JPEG 2000, вейвлет оснований на стандарті стиснення зображень
  • DjVu format uses wavelet-based IW44 algorithm for image compression
  • scaleograms, a type of spectrogrm generated using wavelets instead of a short-time Fourier transformтип спектрограми, яка генерується з використанням вейвлетів, замість віконного перетворення Фур'є.
  • Шаблон:Нп
  • Chirplet transform
  • Time-frequency representation
  • S transform
  • Віконне перетворення Фур'є
  • Ів Мейєр
  • Інгрід Добеші
  • Шаблон:Нп
  • Шаблон:Нп

Література

  • Вейвлет-перетворення у компресії та попередній обробці зображень / О. В. Капшій, О. І. Коваль, Б. П. Русин ; НАН України, Фіз.-мех. ін-т ім. Г. В. Карпенка. − Л. : Сполом, 2008. − 208 с. : іл. − Бібліогр. : с. 187−203 (238 назв). − ISBN 978-966-665-554-0.

Примітки

Шаблон:Reflist

  1. JPEG 2000, for example, may use a 5/3 wavelet for lossless (reversible) transform and a 9/7 wavelet for lossy (irreversible) transform.
  2. A. G. Ramakrishnan and S. Saha, "ECG coding by wavelet-based linear prediction, " IEEE Trans. Biomed. Eng., Vol. 44, No. 12, pp. 1253—1261, 1977.
  3. N. Malmurugan, A. Shanmugam, S. Jayaraman and V. V. Dinesh Chander. "A New and Novel Image Compression Algorithm Using Wavelet Footprints" Шаблон:Webarchive
  4. Ho Tatt Wei and Jeoti, V. "A wavelet footprints-based compression scheme for ECG signals". Шаблон:Cite book
  5. Шаблон:Cite journal
  6. Шаблон:Cite book
  7. «Novel method for stride length estimation with body area network accelerometers» Шаблон:Webarchive, IEEE BioWireless 2011, pp. 79-82
  8. Шаблон:Cite journal
  9. Шаблон:Cite journal