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) возвращает вейвлет-декомпозицию массива 3-D 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соответственно.

    Примеры

    свернуть все

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

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

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

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

    wd1 = wavedec3(X,1,'db1');
    Открыть сценарий в реальном времени

    Сравнение выходных данных 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
    Открыть сценарий в реальном времени

    Укажите фильтры разложения и реконструкции как массив ячеек. Постройте матрицу целых чисел 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');

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

    свернуть все

    Входные данные, указанные как массив 3-D.

    Типы данных: 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 является фильтром разложения нижних частот, и HiD - фильтр разложения верхних частот. Посмотрите wfilters для получения подробной информации.

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

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

    свернуть все

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

    Размер входных данных, возвращаемых в виде вектора 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.

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

    См. также

    | | | | |

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

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

    Поддержка