离散可分离剪切波变换(DSST)及其数值计算.docVIP

离散可分离剪切波变换及其数值计算 1、离散可分离剪切波变换DSST Wang-Q Lim在2010年提出了离散可分离剪切波变换（Discrete Separable Shearlet Transform，DSST），其中一个主要特征就是可以选择可分离的尺度函数及剪切波生成函数（ ），即函数可表示为 ，， 下面将在水平锥上构造可分离剪切波生成函数以及与之相关的尺度函数，垂直锥同理。令是一维紧支撑尺度函数，并选择某个相适的滤波器（所需条件将在后面讨论）使其满足： （1） 如果与之相关的一维紧支撑小波函数，可用相适的滤波器表示为： （2） 那么，此剪切波生成函数可表示为 （3） 尺度函数可表示为 （4） 对于固定的J 0，假设函数可表示为 （5） 这是一个数字实现的常规假设，尺度系数可被看做的抽样值，事实上，通过选择合适的可使。由上面的讨论，可知剪切波系数可通过下式计算 （6） 若不为整数，则需选取。式（6）显示了剪切波系数的离散化方法：首先应用与各向异性抽样矩阵相关联的离散可分离小波变换计算。与式（）相比，要求包含在以下尺度空间中 由上式很容易看出，如果剪切系数为非整数，则剪切矩阵不能将规则网格保留在VJ中，即 为了解决这个问题，定义一个新的尺度空间 可以看出由仿射规则网格沿着水平方向以为因子，可得到尺度空间。基于这种改变，新网格在剪切算子下是不变的，由于 其中，Q=diag（2,1）。因此，对于式（6）中的表示我们有如下引理。 引理1 保留本部分的符号和定义。令表示以为因子的一维上采样算子，表示沿水平方向的一维卷积算子，是三角多项式的傅里叶系数 （7） 可得 其中， 此引理的证明需要以下结论，它来自小波理论的级联算法 命题1 假定和分别满足等式(3)和(4)，对正整数，有 （8） 和 （9） 其中，和分别是三角多项式和的傅里叶系数，的定义如式（7）所示。对固定的j0，的定义为 引理1的证明：令，将其带入式（8），有 （10） 由于是一个二维可分离函数，即，因此有 由式（10）可得 其中，Q=diag（2,1）。使用，最终可表示为 证明完毕。 式（6）的第二步数字化是剪切波本身的离散化，由命题1可得到如下的结论。 引理2 保留此部分的符号和定义。 如前所述，充分应用与各向异性抽样矩阵相关联的离散可分离小波变换，对与，，定义线性变换 （11） 引理1和2共同完成了剪切波系数的离散化算术实现。因此，有如下定理： 定理1 保留此部分的符号和定义。令表示以2j/2为因子的一维上采样算子，表示沿水平方向的一维卷积算子，剪切波系数可化为 其中 图1（a） （12） 在详细描述实现步骤之前，先进一步研究尺度系数。可以看做是通过数字剪切操作在整数格上对的采样。图1（b）时的基本过程。 图1 （a）计算剪切波系数的步骤：沿水平轴细化（上行），与剪切矩阵相关的重采样（中间）和可分离小波变换（下行）；（b），时的水平方向细化过程 事实上，对于任意剪切系数k，式（12）中的滤波器系数都可以很容易预先计算得到。在实际计算中经常假设，有时甚至可以跳过这一卷积步骤。 离散可分离剪切波变换（DSST）的计算方法可通过如下步骤表示： Step 1：在精细尺度以2j/2为因子对其给定的数据一维上采样。 Step 2：在精细尺度，计算上采样后的数据和一维低通滤波器的一维卷积，得到。 Step 3：在精细尺度，根据剪切抽样矩阵重采样得到。由于整数格在剪切矩阵下是不变的，因此重采样的步骤也很简单。 Step 4：在精细尺度，以2j/2为因子对一维下采样后与进行一维卷积。 Step 5：应用可分离小波变换遍历尺度。 2、冗余度分析 实用性要求的主要问题之一就是可控的冗余。为了能够定量的分析冗余离散剪切波变换的冗余度，假设输入数据f是一个由二维尺度函数转换的有限的线性组合，在尺度J表示如下： 上式满足式（5）的假设。为了使结果更具普遍性，在变换中使用任意的抽样矩阵，剪切元素的形式如下所示 那么可以得到如下结论： 命题2 离散可分离剪切波变换（DSST）的冗余度为 证明：首先考虑在固定尺度水平域的剪切元素，观查可知，剪切指数k有个，与尺度矩阵和抽样矩阵相关的变换指数有个。因此，在水平域表示f需要个剪切元素。同理，在垂直域也需要相同个数的剪切元素。最后，在最粗糙尺度尺度函数需要个变换。 遍历所有尺度所需必要剪切元素的个数为 冗余度为此数据与原始