waverec3

3-D реконструкция вейвлета

Синтаксис

X = waverec3(WDEC) C = waverec3(WDEC,TYPE,N) X = waverec3(WDEC,'a',0) X = waverec3(WDEC,'ca',0) C = waverec3(WDEC,TYPE) C = waverec3(WDEC,TYPE,N)

Описание

waverec3 выполняет многоуровневую 3-D реконструкцию вейвлета, начинающую с многоуровневого 3-D разложения вейвлета.

X = waverec3(WDEC) восстанавливает трехмерный массив X на основе многоуровневой структуры разложения вейвлета WDEC. Можно также использовать waverec3 извлекать коэффициенты из 3-D разложения вейвлета.

WDEC структура с полями, показанными в таблице.

C = waverec3(WDEC,TYPE,N) восстанавливает многоуровневые компоненты на уровне N из 3-D разложения вейвлета. N должно быть положительное целое число, меньше чем или равное уровню разложения.

Допустимые значения для TYPE :

Группа из трех символов 'xyz', один на направление, с 'x'Y и 'z' выбранный в наборе {'a'DLH} или в соответствующем прописном наборе {'A'DLH}), где 'A' (или 'L') фильтр lowpass и 'D' (или 'H') фильтр высоких частот.
char 'd' (или 'h' или 'D' или 'H') дает сумму всех компонентов, отличающихся от lowpass.
Символьный 'a' (или 'l' или 'A' или 'L') дает компонент lowpass (приближение на уровне N).

Для экстракции допустимые значения для ВВОДЯТ, то же самое, но снабженный префиксом 'c' или 'C'.

X = waverec3(WDEC,'a',0) или X = waverec3(WDEC,'ca',0) эквивалентно X = waverec3(WDEC)X реконструкция коэффициентов в WDEC на уровне 0.

C = waverec3(WDEC,TYPE) эквивалентно C = waverec3(WDEC,TYPE,N) с N равняйтесь уровню разложения.

`sizeINI`	Размер 3D массива `X`
`level`	Уровень разложения
`mode`	Имя вейвлета преобразовывает дополнительный режим
`filters`	Структура с 4 полями, `LoD`, `HiD`, `LoR`, и `HiR`, которые содержат фильтры, используемые в DWT
`dec`	`N` x 1 массив ячеек, содержащий коэффициенты разложения. `N` равно `7WDEC.level+1`. `dec{1}` содержит компонент lowpass (приближение) на уровне разложения. Приближение эквивалентно операциям `'LLL'` фильтрации. `dec{k+2},...,dec{k+8}` с `k = 0,7,14,...,7(WDEC.level-1)` содержите 3-D коэффициенты вейвлета для мультиразрешения начиная с самого грубого уровня когда `k=0`. Например, если `WDEC.level=3`, `dec{2},...,dec{8}` содержите коэффициенты вейвлета для уровня 3 (`k=0`), `dec{9},...,dec{15}` содержите коэффициенты вейвлета для уровня 2 (`k=7`), и `dec{16},...,dec{22}` содержите коэффициенты вейвлета для уровня 1 (`k=7*(WDEC.level-1)`). На каждом уровне, коэффициентах вейвлета в `dec{k+2},...,dec{k+8}` в следующем порядке: `'HLL'`, `'LHL'`, `'HHL'`, `'LLH'`, `'HLH'`, `'LHH'`, `'HHH'`. Последовательность букв дает распоряжение, в котором отделимые операции фильтрации применяются слева направо. Например, `'LHH'` средние значения, что lowpass (масштабирование) фильтр с субдискретизацией применяется к строкам `X`, сопровождаемый highpass (вейвлет) фильтр с субдискретизацией применился к столбцам `X`. Наконец, фильтр highpass с субдискретизацией применяется к 3-й размерности `X`.
`sizes`	Последовательные размеры компонентов разложения

Примеры

свернуть все

Совершенная реконструкция с 3-D дискретным вейвлетом преобразовывает

Скрипт Open Live Script

Создайте 3-D матрицу, получите вейвлет, преобразовывают вниз к уровню 2 с помощью 'db2' вейвлет, и восстанавливает матрицу, чтобы проверить совершенную реконструкцию.

Создайте 3-D матрицу.

M = magic(8);
X = repmat(M,[1 1 8]);

Получите 3-D дискретное преобразование вейвлета матрицы и восстановите вход на основе 3-D приближения и детализируйте коэффициенты.

wd = wavedec3(X,2,'db2');
XR = waverec3(wd);

Проверьте совершенную реконструкцию с помощью разложения вейвлета вниз для уровня 2.

err1 = max(abs(X(:)-XR(:)))

err1 = 8.6057e-11

Проверьте, что матрица данных является суммой приближения и деталей от уровней 2 и 1. Восстановите сумму компонентов, отличающихся от компонента lowpass, и проверяйте тот X = A + D.

A = waverec3(wd,'LLL');
D = waverec3(wd,'d');
err2 = max(abs(X(:)-A(:)-D(:)))

err2 = 8.6054e-11

Сравните waverec3 и idwt3

Скрипт Open Live Script

Сравните реконструкции уровня 1 на основе операций 'LLH' фильтрации использование idwt3 и waverec3.

M = magic(8);
X = repmat(M,[1 1 8]);
wd = wavedec3(X,2,'db2','mode','per');
dwtOut = dwt3(X,'db2');
Xr = idwt3(dwtOut,'LLH');
Xrec = waverec3(wd,'LLH',1);
norm(Xr(:)-Xrec(:))

ans = 2.7511e-14

Смотрите также

idwt3 | wavedec3 | waveinfo

Документация