Exponenta Banner
exponenta event banner

wavedec3

3-D вейвлет

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

    Синтаксис

    wdec = wavedec3(x,n,wname)
    wdec = wavedec3(x,n,wname,'mode',extmode)
    wdec = wavedec3(x,n,{LoD,HiD,LoR,HiR})

    Описание

    пример

    wdec = wavedec3(x,n,wname) возвращает вейвлет трехмерного массива x на уровне n, с использованием вейвлета, заданного как wname. wavedec3 использует режим расширения по умолчанию 'sym'.

    пример

    wdec = wavedec3(x,n,wname,'mode',extmode) использует заданный режим расширения extmode.

    пример

    wdec = wavedec3(x,n,{LoD,HiD,LoR,HiR}) использует указанные фильтры разложения и реконструкции LoD,HiD и LoR,HiR, соответственно.

    Примеры

    свернуть все

    Открыть Live Script

    Найдите 3-D DWT тома. Постройте матрицу целых чисел 8 на 8 на 8 1 - 64 и сделайте данные 3D. 

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

    Получите 3-D дискретное вейвлет-преобразование на уровне 1 с помощью вейвлета Haar и симметричного режима расширения всей точки по умолчанию. 

    wd1 = wavedec3(X,1,'db1');
    Открыть Live Script

    Сравните выходные данные wavedec3 и dwt3 для иллюстрации упорядоченного расположения 3-D вейвлет, описанных в dec описание поля. 

    X = reshape(1:512,8,8,8);
dwtOut = dwt3(X,'db1','mode','per');
wdec = wavedec3(X,1,'db1','mode','per');
max(abs((wdec.dec{4}(:)-dwtOut.dec{2,2,1}(:))))
    ans = 0

    max(abs((wdec.dec{5}(:)-dwtOut.dec{1,1,2}(:))))
    ans = 0
    Открыть Live Script

    Задайте фильтры разложения и реконструкции как массив ячеек. Постройте матрицу целых чисел 8 на 8 на 8 1 - 64 и сделайте данные 3D.

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

    Получите 3-D дискретное вейвлет-преобразование до уровня 2 с помощью вейвлета конечной фазы Daubechies с двумя моментами исчезновения. Введите фильтры разложения и реконструкции как массив ячеек. Используйте периодический режим расширения. 

    [LoD,HiD,LoR,HiR] = wfilters('db2');
wd2 = wavedec3(X,2,{LoD,HiD,LoR,HiR},'mode','per');

    Входные параметры

    свернуть все

    Входные данные, заданные как трехмерный массив.

    Типы данных: double

    Уровень разложения, заданный как положительное целое число. wavedec3 не устанавливает ограничение максимального уровня. Посмотрите wmaxlev.

    Типы данных: double

    Анализ вейвлета, заданный как вектор символов или строковый скаляр.

    Примечание

    wavedec3 Поддержки только Тип 1 (ортогональный) или Тип 2 (биортогональные) вейвлеты. Посмотрите wfilters для получения списка ортогональных и биортогональных вейвлеты.

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

    mode

    Режим расширения DWT

    'zpd'

    Нулевой внутренний номер

    'sp0'

    Сглаживайте расширение порядка 0

    'spd' (или 'sp1')

    Сглаживайте расширение порядка 1

    'sym' или 'symh'

    Симметричное расширение (половина точки): граничное значение симметричной репликации

    'symw'

    Симметричное расширение (вся точка): граничное значение симметричной репликации

    'asym' или 'asymh'

    Антисимметричное расширение (половина точки): краевое значение антисимметричной репликации

    'asymw'

    Антисимметричное расширение (вся точка): краевое значение антисимметричной репликации

    'ppd', 'per'

    Периодизированное расширение

    Если длина сигнала нечетная и mode является 'per'дополнительная выборка, равная последнему значению, добавляется вправо, и расширение выполняется в 'ppd' режим. Если длина сигнала четная, 'per' эквивалентно 'ppd'. Это правило также применяется к изображениям.

    Глобальная переменная, управляемая dwtmode задает режим расширения по умолчанию. Посмотрите dwtmode для описаний режима расширения.

    Фильтры вейвлет-разложения, сопоставленные с ортогональным или биортогональным вейвлетом, заданы как реальные векторы четной длины. LoD - lowpass разложения, и HiD - фильтр высокочастотного разложения. Посмотрите wfilters для получения дополнительной информации.

    Фильтры вейвлет-реконструкции, сопоставленные с ортогональным или биортогональным вейвлетом, заданы как реальные векторы четной длины. LoR - lowpass фильтр реконструкции, и HiR - фильтр реконструкции верхних частот. Посмотрите wfilters для получения дополнительной информации.

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

    свернуть все

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

    Размер входных данных, возвращенный как вектор 1 на 3.

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

    Имя режима расширения преобразования вейвлета, возвращаемое как вектор символов.

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

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

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

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

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

    Коэффициенты разложения, возвращенные как N-на-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. Наконец, фильтр высокой частоты с понижающей дискретизацией применяется к 3-ьей размерности x.

    Последующие размеры компонентов разложения, возвращенные как n+ 1 на 2 матрица.

    См. также

    | | | | |

    Введенный в R2010a

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

    Поддержка