剪切波的构造Shearlet是一类新的多尺度几何分析方法,该方法通过对基本函数的膨胀、剪切和平移变换来构造,体现了函数的几何和数学特性,如近几年来许多领域的研究学者所强调的函数的方向性、尺度和振荡等。Shearlet可以在广义多分辨率分析的框架中研究,这样就可以获得像小波一样的迭代算法,并推广到经典的级联算法。因此Shearlet变换作为一种新型的多尺度几何分析工具为图像处理领域的研究人员所广泛接受。1、Shearlet及其变换的定义函数f(x)的连续Shearlet变换为:,,(,,),fastSHastf(1.1)其中,3/411,,()(())astxaABxt(1.2)为剪切波函数,aR为尺度参数,sR为剪切参数,2tR为平移参数,12,0;0,aa是各向异性膨胀矩阵,1,;0,1s是剪切矩阵。对任何2121ˆ(,),0R,令满足21121ˆˆˆ()()()(1.3)其中,ˆ为,,ast的傅里叶变换,1为连续小波函数,1ˆRC(),1ˆ54,1414,54supp,2为bump函数,2ˆRC(),2ˆ1,1supp,在区间1,1上20且21。由以上定义可得,,()astx的傅里叶变换为342122,,1121ˆˆˆitastaeaas显然,剪切波的几何性质在频域上更为直观。由1ˆ和2ˆ的支撑条件很容易看到,,ˆast有如下的频域支撑:,,121212112supp,:,,,/22astsaaaaa由上式可知,在不同尺度,,,ˆast支撑在以原点对称、以s为斜率的梯形对上;改变剪切参数s,支撑区域可获得保持面积不变的旋转;旋转区域由尺度参数a控制,随着a→0,支撑区间逐渐变窄。图1.1给出了)3,1(),0,4/1(),0,1(sasasa时的,,ˆast频域支撑。a=1,s=-3a=1/4,s=0a=1,s=012图1.1不同a和s值时,,ast的频域支撑连续剪切波变换的平移参数可检测到所有奇异点的位置,而剪切参数则可显示出奇异曲线的方向。2、Shearlet的离散化为了实现剪切波的离散化,令尺度参数2jjajZ,剪切参数122,2jjkskakkZ,以及平移参数2,,,,jkmasjjmtDmZ。假定22101ˆ21,8jjww(2.1)和222112ˆ21,1jjjkwkw(2.2)由式(2.1)和式(2.2)可知:对任何120,C,有1222122(0)221200100=-221ˆ()(2)(2)1jjjjjkjjjjkkkAB其中,2012121ˆ,:18,1RC,如图2.1(a)所示,即函数00()jkAB形成0C的一个剖分,如图2.1(b)所示。由以上的讨论,可知集合3(0)(0)2224,,00()2():0,22,jkjjjjkmxxmjkmBAZ是22200ˆ(){():supp}VLfLfCC得一个紧框架。其中,1022002Α,01101Β由图2.1(b)可以看出,剪切波的每个元素,,jkm支撑在梯形对上,每一个梯形包含在一个大小近似为222jj的盒子里,方向沿着斜率为2jl的直线。图2.1(a)水平锥0C和垂直锥1C(b)剪切波的频域剖分图同样可以构造一个21()VLC的紧框架,其中,1C是垂直锥2112212ˆ,:18,1RC。由下式给定(1)11222ˆˆˆ()()()(2.3)则集合3(1)(1)2224,,11()2():0,22,jkjjjjkmxxmjkmBAZ是21()VLC的一个紧框架。其中,1212002Α,11011Β。最后,令22LR满足:对任何2R,有2121222(0)(1)001100k=-2k=-2ˆˆˆ()()1jjjjjkjkjjABAB上式暗含2ˆsupp18,18,且2ˆ=1。因此剪切波集合为2()222,,()():():0,22,,0,1mdjjjkmxxmmxjkmdZZ3、剪切波的主要性质从上述可以看出,剪切波具有以下良好特性:(1)剪切波具有非常好的局部化特性。在上述剪切波的构造中,剪切波在频域内是紧支撑的,并且在空域内具有快速的衰减特性。(2)剪切波满足抛物线尺度化特性。每一个元素mkj,,ˆ支撑在一个梯形对上,且每个梯形对包含在一个大小近似为jj222的盒子内,如图2.1(b)所示。这是因为剪切波具有非常好的局部化特性,在空域内每个mkj,,本质上支撑在一个大小为jj222的盒子里面。当j时,元素的支撑区间会逐渐变窄。(3)剪切波体现了非常高的方向敏感性。元素mkj,,ˆ的方向是沿着斜率为jl2的直线。相应的元素mkj,,的方向是沿着斜率为jl2的直线,并且方向的数目随着尺度的不断细化而逐层加倍。(4)剪切波是空域局部化的。对任意固定的尺度和方向,剪切波可以通过在格2上平移来获得。(5)剪切波具有接近最优的稀疏表示性能。