NSCT变换定义及部分物理意义（附Matlab代码）

NSCT（非下采样轮廓波变换）

一、定义

NSCT变换由非下采样金字塔(non-subsampled pyramid，NSP)和非下采样方向滤波器组(non-subsampled directional filter bank，NSDFB)构成。它与轮廓波变换类似，包括子带分解和方向变换两个步骤。先由拉普拉斯金字塔变换对图像进行多分尺度分解，捕捉高频信息，再通过方向滤波器组进一步分解得到不同方向子带。NSCT变换分解过程见图1。

         图2 非下采样金字塔基本结构

（1）NSCT分解在轮廓波变换的基础上取消了其在金字塔分解和方向分解和重构中的上下采样，图2为NSP的分解和重建示意图，其中Hk(Z)(k=0,1)表示分解滤波器，Gk(Z)(k=0,1)表示重构滤波器。一次NSP滤波可产生一个低通子带图像（y0）和一个带通方向子带图像(y1)，尺寸均与源图像大小相同。
（2）NSDFB第一级分解由水平和垂直两个方向的扇形滤波器对图像做分解，得到两个方向上的高频子带。在第二级分解中级联象限滤波器，分解为4个方向，在第三级分解中再级联锲型滤波器进行多方向分解。图3为NSDFB三级分解结构图和频率分割图。将频域平面分成多方向的楔形块状结构，代表源图像在响应方向上的细节信息，最终形成一个由多个双通道NSDFB构成的滤波结构。

               图3 非下采样方向滤波器组基本结构

以经典的zoneplate图像为例，对其进行3层NSCT分解。如图4©所示为NSCT一级分解的垂直和水平两个方