Exponenta Banner
exponenta event banner

waverec3

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

свернуть все на странице

    Синтаксис

    x = waverec3 (wdec)
    c = waverec3 (wdec, тип, n)

    Описание

    x = waverec3(wdec) восстанавливает массив 3-D x на основе многоуровневой структуры вейвлет-декомпозиции wdec.

    пример

    c = waverec3(wdec,type,n) реконструирует или извлекает на уровне n многоуровневые компоненты, указанные type. Если type начинается с 'c' или 'C', waverec3 извлекает указанные компоненты. В противном случае waverec3 восстанавливает компоненты.

    x = waverec3(wdec,'a',0) и x = waverec3(wdec,'ca',0) эквивалентны x = waverec3(wdec), где 'a' задает компонент нижних частот. x - реконструкция коэффициентов в wdec на уровне 0.

    c = waverec3(wdec,type) эквивалентно c = waverec3(wdec,type,wdec.level).

    Примеры

    свернуть все

    Открыть сценарий в реальном времени

    Создайте матрицу 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. Реконструируйте сумму компонентов, отличных от компонента нижних частот, и проверьте, что X = A + D.

    A = waverec3(wd,'LLL');
D = waverec3(wd,'d');
err2 = max(abs(X(:)-A(:)-D(:)))
    err2 = 8.6054e-11
    Открыть сценарий в реальном времени

    Сравнение реконструкций уровня 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

    Входные аргументы

    свернуть все

    Вейвлет-декомпозиция, заданная как структура. Структура является выходом wavedec3 и имеет следующие поля:

    Размер массива 3-D, заданный как вектор 1 на 3.

    Уровень разложения, заданный как целое число.

    Имя режима расширения вейвлет-преобразования, заданного как символьный вектор.

    Вейвлет-фильтры, используемые для разложения, задаются как структура со следующими полями:

    • LoD - фильтр разложения нижних частот

    • HiD - фильтр разложения верхних частот

    • LoR - фильтр разложения нижних частот

    • HiR - фильтр разложения верхних частот

    Коэффициенты разложения, заданные как массив N-by-1 ячеек, где N равно 7 ×wdec.level+1.

    dec{1} содержит компонент нижних частот (аппроксимацию) на уровне разложения. Аппроксимация эквивалентна операциям фильтрации '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' означает, что фильтр нижних частот (масштабирование) с понижающей дискретизацией применяется к строкам x, за которым следует фильтр верхних частот (вейвлет) с понижающей дискретизацией, применяемой к столбцам x. Наконец, фильтр верхних частот с понижающей дискретизацией применяется к 3-му размеру x.

    Последовательные размеры компонентов разложения, указанные как wdec.levelМатрица + 1 на 2.

    Тип реконструкции или извлечения, заданный как вектор символа или скаляр строки. Для реконструкции действительные значения type являются:

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

    • Случайная работа 'd' (или 'h' или 'D' или 'H') дает сумму всех компонентов, отличных от компонента нижних частот.

    • Случайная работа 'a' (или 'l' или 'A' или 'L') дает компонент нижних частот (аппроксимация на уровне n).

    Для извлечения коэффициентов действительные значения для type те же, но с префиксом 'c' или 'C'.

    Уровень разложения, заданный как целое число.

    Выходные аргументы

    свернуть все

    Реконструкция, возвращенная как 3-D массив размера sz(1)около-sz(2)около-sz(3), где sz = wpdec.sizeINI.

    Извлеченные коэффициенты, возвращенные в виде массива 3-D.

    См. также

    | |

    Представлен в R2010a

    Документация инструментария вейвлет

    Поддержка