idwt3

Одноуровневое обратное дискретное 3-D вейвлет

Синтаксис

X = idwt3(WT)
C = idwt3(WT,TYPE)

Описание

idwt3 команда выполняет одноуровневую трехмерную реконструкцию вейвлета, начиная с одноуровневого трехмерного разложения вейвлет.

X = idwt3(WT) вычисляет одноуровневый восстановленный трехмерный массив X на основе трехмерного разложения вейвлет, сохраненного в структуре WT. Эта структура содержит следующие поля.

sizeINI

Размер трехмерного массива X.

mode

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

filters

Структура с 4 полями, LoD, HiD, LoR, HiR, которые содержат фильтры, используемые для DWT.

dec

Массив ячеек 2 x 2 x 2, содержащий коэффициенты разложения.

dec{i,j,k}, i, j, k = 1 или 2 содержит коэффициенты, полученные путем lowpass фильтрации (для i или j или k = 1) или высокочастотную фильтрацию (для i или j или k = 2).

C = idwt3(WT,TYPE) вычисляет одноуровневый восстановленный компонент на основе трехмерного разложения вейвлет. Допустимые значения для TYPE являются:

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

  • The char 'd' (или 'h' или 'D' или 'H'), которая задает сумму всех компонентов, отличных от компонента lowpass.

Примеры

свернуть все

Определите исходные данные 3-D.

X  = reshape(1:64,4,4,4)
X = 
X(:,:,1) =

     1     5     9    13
     2     6    10    14
     3     7    11    15
     4     8    12    16


X(:,:,2) =

    17    21    25    29
    18    22    26    30
    19    23    27    31
    20    24    28    32


X(:,:,3) =

    33    37    41    45
    34    38    42    46
    35    39    43    47
    36    40    44    48


X(:,:,4) =

    49    53    57    61
    50    54    58    62
    51    55    59    63
    52    56    60    64

Разложение X использование 'db1'.

wt = dwt3(X,'db1');

Восстановление X от коэффициентов. Проверьте, что восстановленные данные согласуются с исходными данными для точности машины.

XR = idwt3(wt);

dff = max(abs(X-XR))
dff = 
dff(:,:,1) =

   1.0e-13 *

    0.0266    0.0355    0.0888    0.1066


dff(:,:,2) =

   1.0e-13 *

    0.1066    0.1066    0.2132    0.2132


dff(:,:,3) =

   1.0e-13 *

    0.1421    0.1421    0.2132    0.2132


dff(:,:,4) =

   1.0e-13 *

    0.3553    0.3553    0.2842    0.2842

Вычислите восстановленное приближение, которое состоит из lowpass компонента.

A  = idwt3(wt,'aaa');

Вычислите сумму всех компонентов, отличных от lowpass компонента.

D  = idwt3(wt,'d');

Восстановите компонент, сопоставленный с lowpass в направлениях x и z, и highpass в направлении y.

ADA  = idwt3(wt,'ada');

См. также

| |

Введенный в R2010a