本方法涉及图像压缩技术领域，特别涉及一种自适应的多尺度分块压缩感知采样方法。
背景技术：
随着现代社会的信息化发展进程，人们对接受信息的需求越来越迫切，多媒体作为信息技术的一种，在社会生活中的作用越来越突出。迄今为止，包括汉字、动画、声音、图像、影像等都是多媒体领域相关的元素，其中，图像的数据量特别巨大，给信号处理过程中的存储空间、信号带宽和计算速度各方面带来很大的压力。为了解决上述问题，对图像的压缩处理势在必行，图像压缩技术的研究势必对人们的现代化生活产生重要意义。
传统的图像压缩标准，如JPEG(Joint Photographic Experts Group)、JPEG2000、MPEG4(Moving Pictures Experts Group)等都是较为成熟的图像压缩标准，然而这些标准都基于香农采样定理，要求采样频率高于信号带宽的两倍，对数据的存储、处理和传输造成巨大挑战。压缩感知理论突破了奈奎斯特采样率的限制，采样和压缩同步进行，大大降低了传感器的采样和计算成本，因而在众多图像压缩技术中脱颖而出，在科学和工程的许多领域有着重要的理论与实际应用价值。
目前，多尺度分块压缩感知作为压缩感知的拓展研究，其根据信号的多分辨分析特性，基于分块压缩感知对采样率的分配应用了小波域内分尺度加权调整策略，取得了不错的图像处理效果。然而，多尺度分块压缩感知忽略了小波域图像中包含大量先验信息的低频系数，且同一层小波系数采用一致的采样率。事实上，即使是同一层小波系数内部，各系数块之间含有的纹理信息和方向性也是不相同的，这势必会因为采样率分配不当而导致重构效果不尽理想。
技术实现要素：
本发明目的在于克服上述现有技术的不足，提供一种结合纹理信息和方向估计的自适应多尺度分块压缩感知采样方法，实现自适应采样，提高图像的重构质量。
为实现上述目的，本发明的技术方案包括如下步骤：
步骤1，对原图像进行三层小波分解，得到小波域图像；
所述的小波分解采用9/7正交小波分解，分解方式是：
首先，将原图像分解为第一级的低频系数子带和三个高频系数子带；
接着，将第一级分解的低频子带进行进一步分解，得到下一层的低频及高频子带，而由第一级分解获得的三个高频子带的数据将不再参与分解；
最后，按照此方式再进行一步分解，即可得到三层小波分解后的小波域图像；
步骤2，提取步骤1得到的小波域图像中的高频系数进行分层、分子带、分块；
所述的分层是指将小波域图像按照三层分解方式分为三层；
所述的分子带是指除低频部分外，得到九个水平、垂直及对角高频子带；
所述的分块是指将各子带分成64个大小相同且不重叠的图像块；
步骤3，提取步骤1得到的小波域图像中的低频系数，对低频系数进行小波逆变换得到预估计图像，将预估计图像分块分成64个大小相同且不重叠的图像块；
步骤4，计算步骤3得到的预估计图像各块的灰度熵和所有图像块的灰度熵之和；
各图像块的灰度熵用Hj表示，j表示第j(j＝1,2，…64)个图像块；灰度熵的计算公式为：
其中，pi表示图像中灰度值为i的像素在整个图像中所占的比例，由灰度直方图获得；
步骤5，利用如下经验公式将各块灰度熵转换为初步自适应采样率：
其中，rj表示第j个图像块的初步自适应采样率，S表示本文中预设的目标采样率，h表示分块数目，为了防止初步自适应采样率过低，最低采样率的阈值设置为Smin，取值为：
步骤6，对步骤3得到的预估计图像各块利用图像块的梯度信息进行方向估计；
所述的方向估计的详细步骤为：
(6.1)对每个图像块的梯度矩阵G进行奇异值分解：
其中U是N×2维的正交矩阵，S是2×2维矩阵，表示主方向上的能量，v1与v2相互正交，v1表示梯度场主方向，v2表示图像块的主方向；
(6.2)利用奇异值s1和s2之间的差值对图像块的方向显著性进行判断：
R∈[0,1]，R的值大表示图像块的方向性较显著；反之，R的值小则说明图像的方向性并不明显；
当R→0时，表示图像块不具有明显的方向性；当R≥τ(τ为判断方向显著性的阈值)时，图像块有明确的主方向；
(6.3)对那些主方向非常明确的图像块，利用下式计算各块的主方向角θ：
θ＝180/π·arctan(v22/v21)
(6.4)按照(6.3)得到的主方向角对各图像块的主方向进行判断，方向角和方向的关系规定如下：
步骤7，将步骤5得到的初步自适应采样率按照步骤6得到的方向信息对各层各子带各图像块进行采样率的分配：主方向明确的图像块按主方向分配采样率，没有明确主方向的图像块平均分配采样率；低频系数全采样；
所述的采样率分配的详细步骤为：
(7.1)求各层子采样率Sl,j，即第l层各子带的第j个图像块的采样率：
利用已求得的第j个图像块的初步自适应采样率rj作为该图像块各层子采样率的目标采样率，应用如下公式求解Sl,j：
Sl,j＝WlS' (1)
Wl＝16L-l+1 (3)
其中，rj和Wl已知，通过式(2)可以很容易地解出S'，再利用式(1)就能求出Sl,j，这个求解过程会产生一个或多个Sl,j＞1；因此，需要修改求解方式和步骤来保证各层的Sl,j≤1；特别的，由式(1)和(2)求得S'和S1,j后，检查是否有S1,j＞1，若是，设置S1,j＝1，若否，保留S1,j的值；按此法继续依次求解S2,j、S3,j，使得所有解均小于等于1；
(7.2)求各层各子带中各图像块的采样率Sl,s,j，即第l(l＝1,2,3)层s(s∈H,V,D)子带第j个图像块的采样率：
利用(7.1)中获得的各层子采样率Sl,j和步骤6得到的方向性，对各层各子带中各图像块分配采样率，假设第j个图像块的方向性是水平方向，那么水平方向的子带应多采样，垂直和对角方向的子带应少采样，若第j个图像块没有明确的主方向，那么各个方向的图像块平均分配该层采样率；
步骤8，计算各层各子带各图像块的观测矩阵Φl,s,j；
首先，利用(7.2)获得的各层各子带中各图像块的采样率Sl,s,j计算出每个对应块的自适应观测数目，计算方法是Ml,s,j＝Sl,s,j×Bl2；
接着，由高斯分布N(0,1)构造的随机高斯矩阵正交化得到一个Bl×Bl维正交变换的随机矩阵
最后，Φl,s,j由从矩阵中随机抽取的Ml,s,j个行向量构成；
步骤9，将步骤8得到的观测矩阵Φl,s,j代入式yl,s,j＝Φl,s,jxl,s,j，计算各层各子带各图像块的观测集合；
其中，xl,s,j表示小波域内l层子带s的第j个图像块，yl,s,j表示相应的观测集合；
步骤10，采用平滑投影Landweber重构方法恢复图像；
重构过程由投影和滤波交替迭代进行，直到达到停止迭代的条件。
有益效果
本发明充分利用代表先验信息的小波域低频系数，首先将灰度熵应用于衡量图像的纹理信息，并以此量化值作为确定初步自适应采样率的依据；然后利用图像自身具有的方向性，采用估计主方向的方法对图像进行了更深入和细致的衡量。通过结合图像的纹理信息和方向性克服了采样率分配不合理的现象，其明显的优点是：重构图像的峰值信噪比(Peak Signal to Noise Ratio,PSNR)、结构相似度(Structural Similarity Index Measurement,SSIM)以及视觉效果比原算法均有提升，在相同的采样率下能够提高重构图像质量，广泛应用于医学领域、军事领域和公共安全领域，降低图像处理设备成本。
下面结合实施例附图对本发明作进一步说明。
附图说明
图1是本发明方法的整体流程框图；
图2是以Lena为例预估计图像的分块情况；
图3是以Lena为例原图像经过三层小波分解后的分块情况；
图4是以Lena为例各图像块初步自适应采样率与固定采样率的对比图；
图5是以Barbara为例各图像块的方向估计情况；
图6是各测试图像重构性能和视觉效果对比图。
具体实施方式
参照图1，本发明的实现步骤如下：
步骤1，对原图像进行三层小波分解，得到小波域图像；
本实施例中将三幅512×512的标准测试图像Lena、Barbara、Airplane作为原图像，所述的小波分解采用9/7正交小波分解，分解方式是：
首先，将原图像分解为第一级的低频系数子带和三个高频系数子带；
接着，将第一级分解的低频子带进行进一步分解，得到下一层的低频及高频子带，而由第一级分解获得的三个高频子带的数据将不再参与分解；
最后，按照此方式再进行一步分解，即可得到三层小波分解后的小波域图像；
步骤2，提取步骤1得到的小波域图像中的高频系数进行分层、分子带、分块；
所述的分层是指将小波域图像按照三层分解方式分为三层；
所述的分子带是指除低频部分外，各层得到三个水平、垂直及对角高频子带；
所述的分块是指将各子带分成64个大小相同且不重叠的图像块，则各层分块大小分别为32,16,8；
分层、分子带和分块具体情况参照图2；
步骤3，提取步骤1得到的小波域图像中的低频系数，对低频系数进行小波逆变换得到预估计图像，将预估计图像分块分成64个大小相同且不重叠的图像块，分块情况参照图3；
步骤4，计算步骤3得到的预估计图像各块的灰度熵和所有图像块的灰度熵之和；
各图像块的灰度熵用Hj表示，j表示第j(j＝1,2，…64)个图像块；灰度熵的计算公式为：
其中，pi表示图像中灰度值为i的像素在整个图像中所占的比例，由灰度直方图获得；
步骤5，利用如下经验公式将各块灰度熵转换为初步自适应采样率：
其中，rj表示第j个图像块的初步自适应采样率，S表示本文中预设的目标采样率，h(本实施例中h＝64)表示分块数目，为了防止初步自适应采样率过低，最低采样率的阈值设置为Smin，取值为：
步骤6，对步骤3得到的预估计图像各块利用图像块的梯度信息进行方向估计；
所述的方向估计的详细步骤为：
(6.1)对每个图像块的梯度矩阵G进行奇异值分解：
其中U是N×2维的正交矩阵，S是2×2维矩阵，表示主方向上的能量，v1与v2相互正交，v1表示梯度场主方向，v2表示图像块的主方向；
(6.2)利用奇异值s1和s2之间的差值对图像块的方向显著性进行判断：
R∈[0,1]，R的值大表示图像块的方向性较显著；反之，R的值小则说明图像的方向性并不明显；
当R→0时，表示图像块不具有明显的方向性；当R≥τ(本实施例中τ＝0.1)时，图像块有明确的主方向；
(6.3)对那些主方向非常明确的图像块，利用下式计算各块的主方向角θ：
θ＝180/π ·arctan(v22/v21)
(6.4)按照(6.3)得到的主方向角对各图像块的主方向进行判断，方向角和方向的关系规定如下：
步骤7，将步骤5得到的初步自适应采样率按照步骤6得到的方向信息对各层各子带各图像块进行采样率的分配：主方向明确的图像块按主方向分配采样率，没有明确主方向的图像块平均分配采样率；低频系数全采样；
所述的采样率分配的详细步骤为：
(7.1)求各层子采样率Sl,j，即第l层各子带的第j个图像块的采样率：
利用已求得的第j个图像块的初步自适应采样率rj作为该图像块各层子采样率的目标采样率，应用如下公式求解Sl,j：
Sl,j＝WlS' (1)
Wl＝16L-l+1 (3)
其中，rj和Wl已知，通过式(2)可以很容易地解出S'，再利用式(1)就能求出Sl,j，这个求解过程会产生一个或多个Sl,j＞1；因此，需要修改求解方式和步骤来保证各层的Sl,j≤1；特别的，由式(1)和(2)求得S'和S1,j后，检查是否有S1,j＞1，若是，设置S1,j＝1，若否，保留S1,j的值；按此法继续依次求解S2,j、S3,j，使得所有解均小于等于1；
(7.2)求各层各子带中各图像块的采样率Sl,s,j，即第l(l＝1,2,3)层s(s∈H,V,D)子带第j个图像块的采样率：
利用(7.1)中获得的各层子采样率Sl,j和步骤6得到的方向性，对各层各子带中各图像块分配采样率；
假设第j个图像块的方向性是水平方向，那么水平方向的子带应多采样，垂直和对角方向的子带应少采样，本实施例中设置采样率Sl,s,j为：
若第j个图像块没有明确的主方向，那么各个方向的图像块平均分配该层采样率；
步骤8，计算各层各子带各图像块的观测矩阵Φl,s,j；
首先，利用(7.2)获得的各层各子带中各图像块的采样率Sl,s,j计算出每个对应块的自适应观测数目，计算方法是Ml,s,j＝Sl,s,j×Bl2；
接着，由高斯分布N(0,1)构造的随机高斯矩阵正交化得到一个Bl×Bl维正交变换的随机矩阵
最后，Φl,s,j由从矩阵中随机抽取的Ml,s,j个行向量构成；
步骤9，将步骤8得到的观测矩阵Φl,s,j代入式yl,s,j＝Φl,s,jxl,s,j，计算各层各子带各图像块的观测集合；
其中，xl,s,j表示小波域内l层子带s的第j个图像块，yl,s,j表示相应的观测集合；
步骤10，采用平滑投影Landweber重构方法恢复图像；
首先由Φl,s,j和yl,s,j求出小波域图像，经小波逆变换和维纳滤波去块效应得到图像的初始解后，投影和维纳滤波交替迭代，直到达到停止迭代的条件，得到最优解，即恢复的图像。
本发明的效果可以通过以下仿真内容和结果进一步说明：
仿真1，图像的纹理信息转换成的初步自适应采样率，图4是以Lena为例，各图像块初步自适应采样率与固定采样率的对比图。各图像块的初步自适应采样率是由代表图像纹理信息的灰度熵决定的，纹理信息越复杂，灰度熵越大，该图像块的初步自适应采样率的值越高；反之，纹理平滑图像块的初步自适应采样率的值低；
仿真2，图像各块的方向估计结果，图5是以Barbara为例，各图像块的方向估计情况。其中，N表示该块不满足方向显著性条件R≥τ，没有明确的主方向，H，D，V分别表示该块主方向为水平、对角和垂直方向，即表征了该图像块的方向细节；
仿真3，重构图像的峰值信噪比PSNR和结构相似度SSIM，下述表1和表2分别描述针对三幅标准测试图像，对比分块压缩感知采样方法(Block Compressed Sensing with Smoothed Projected Landweber Reconstruction,BCS-SPL)、基于小波域的多尺度分块压缩感知方法(Multiscale Block Compressed Sensing with Smoothed Projected Landweber Reconstruction,MS-BCS-SPL)和本发明的自适应多尺度分块压缩感知采样方法(Adaptive Multiscale Block Compressed Sensing AMS-BCS-SPL)，在采样率0.3～0.6下的峰值信噪比PSNR和结构相似度SSIM的数值对比结果(较高的数值用粗体标记)；
表1，各重构图像的峰值信噪比PSNR(dB)
表2，各重构图像的结构相似度SSIM
由表1和表2可以看出，在大部分采样率下，本发明的自适应多尺度分块压缩感知采样方法的重构性能优于其他压缩感知采样方法，所以本发明在相同采样率下提高了图像的重构质量；
仿真4，重构图像的视觉效果，参照图6各测试图像重构性能和视觉效果对比图。其中图6(a)是采样率为0.3时，Lena整体图像的重构性能和视觉效果；图6(b)是采样率为0.6时，Barbara图像的重构性能和局部放大图的视觉效果对比图；图6(c)是采样率为0.5时，Airplane图像的重构性能和局部放大图的视觉效果对比图。均给出了相应的峰值信噪比PSNR和结构相似度SSIM值。由图6可见，本发明的重构图像性能和视觉效果明显提高。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的思想和原理之内,所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。