conturlet变换-nsct_toolbox工具箱源码简要分析

因为需要编写关于contourlet的相关代码，苦于网上没有太多的相关资料，无奈，硬着头皮强行啃食nsct_toolbox的源码，现总结出一些经验和自己的理解，供苦×的广大学子参考，里面还有很多疑问点，希望看到该篇文章的朋友不吝赐教，有错误之处还请指正。

• y = nsctdec(x, levels, [dfilt, pfilt] )
  这就是非下采样（也可选择采样滤波器）的contourlet分解的主函数

y = nsctdec(x, levels, [dfilt, pfilt] )

x:图像的输入矩阵，双精度数据类型
levels:方向滤波器组分解层数向量,注意是一组向量值
dfilt:方向滤波器组
pfilt:塔式分解的滤波器组
out:输出的结果是一组向量集合，集合第一个向量是低频分量的二维数组，接下来是从小到大排列高频子带。
首先是一些滤波器的准备：

filters = cell(4) ;//初始化一个数组存放滤波器
[h1, h2] = dfilters(dfilt, 'd');//得到方向二维滤波器组，一低一高
//得到比例尺度
h1 = h1./sqrt(2) ;
h2 = h2./sq