Exponenta Banner
exponenta event banner

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 равно 7*WDEC.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

Последовательные размеры компонентов разложения

Примеры

свернуть все

Скрипт 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
Скрипт 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

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

| |

Введен в R2010a

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

Поддержка