因为需要编写关于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