关于nsct

最近在学习nsct，一点点的过来总结，先贴从前辈那看来有兴趣的
        众所周知，所谓NSCT就是取消了contourlet中对于图像的上下采样，而改用对滤波器进行上采样，接下来问题就来了。用于塔式分解的滤波器的上采样比较好理解。但是对于方向滤波器的上采样就不太理解了。
        有的文献上说每一级通道的方向滤波器都是对上一级的方向滤波器进行五株上采样（梅花上采样Q0=[1 -1;1 1]或Q1=[1 1;-1 1]）；有的说第一级用扇型滤波器，第二级用棋盘型滤波器，第三级及第三级之后的都用平行四边形的滤波器。
        按说这两种描述应该是等价的，而且用toolbox中的代码实现，第二级的棋盘型滤波器确实是对第一级扇形进行五株上采样得到的，但之后的就有点不对了，toolbox代码中不再是是对棋盘滤波器进行五株采样，而是用推挤变换（平行四边形变换R0=[1 1;0 1]等），虽说得到了平行四边形滤波器组，但出来的频谱并不是预期的那样，但toolbox代码另外编写的C程序的代码进行卷积运算，很难理解，不过它确实是实现了相应的滤波。

        我现在就只想用matlab来编写，对于上采样我用的是仿射imtransform变换，但是出来的结果和toolbox的不一样，原因可能是imtransform不是像toolbox用的零插值，而是用的最近邻等方法。所以现在就卡在这个地方了呢，求大神帮忙！！！

--------------------------http://www.ilovematlab.cn/thread-246372-1-1.html

下采样：

　　定义：对于一个样值序列间隔几个样值取样一次，这样得到新序列就是原序列　　的下采样。

　　采样率变化主要是由于信号处理的不同模块可能有不同的采样率要求。下采样相对于最初的连续时间信号而言，还是要满足采样定理才行，否则这样的下采样会引起信号成分混叠。

　　下采样就是抽取，是多速率信号处理中的基本内容之一。在不同应用场合，下采样可以带来许多相应的好处。就以在最常见的数字接收机中为例，最后要得到的基带信号的采样率等于符号速率，这个速率是比较低的，但通常的做法并不是直接以这个采样率对模拟信号进行采样，而是采用高的多（几十甚至上百倍）的采样率，这样可以提高采样得到的信号的信噪比，然后再用数字的方法对信号进行多级的滤波和抽取，直到最后信号的采样率与符号速率相等。这样处理可以获得的信噪比增益为最初采样率与最后输出信号采样率之比。

　　不同的采样率之前，是有一个带宽与该采样率对应的滤波器的，

　　采样率越高，滤波器带宽就越大，对于宽带噪声而言（噪声带宽高于最高的采样率），

　　通过的噪声功率就越高（噪声功率即功率谱密度乘上带宽，也即是每采样值中噪声分量的平方取均值。）

　　信号功率在采样前后始终是没有变化的（信号功率即是每采样值中信号分量的平方取均值）。

　　对于窄带噪声或者窄带干扰（噪声或者干扰带宽低于最高采样率），下采样获得“信噪比增益

　　为最初采样率与最后输出信号采样率之比”