压缩感知算法在细胞图像压缩重建中研究.docVIP

压缩感知算法在细胞图像压缩重建中研究 　　摘 要：由于传统的压缩技术局限于奈奎斯特采样定理的限制，其对于存储空间及数据的处理能力都要求比较高，这就给硬件的实现带来了很大困难。压缩感知技术与传统的压缩技术相比，其优势是：数据采样与数据压缩是同时进行的。压缩感知理论的提出很好地克服了资源浪费的现象。该文利用了Contourlet变换工具，对细胞图像进行了边缘轮廓的提取。提出了一种新的采样方法：对细胞图像的不同频率部分，采取不同的采样方法，为了使细胞图像的压缩比能够提高，该文结合细胞图像的特点，对Contourlet算法进行了改进，然后采用阈值迭代的方法，准确、快速地来对细胞图像进行重构。 　　关键词：压缩感知 细胞图像压缩 稀疏表示 分类字典 　　中图分类号：TP391 文献标识码：A 文章编号：1672-3791（2016）02（b）-0157-02 　　图像压缩技术诞生于1948年，但是在现实的应用中，为了节约图像传输、存储及处理方面的成本，通常情况下是对图像进行压缩处理。在2006年，美国著名的数学教授Candes与Donho提出了压缩感知理论[1]，此理论的提出为上述问题的解决带来新的思路。 　　1 稀疏性定义 　　信号能够进行稀疏表示是压缩感知能够进行的前提条件。信号在时域中是不存在稀疏性，但是在某些域中存在，人多眼睛识别东西是有限度的，有些情况视觉是无法察觉到的，这样，就可以在视觉察觉不到信号损失时，丢掉小值系数，用留下的大值系数进行编码及传输的处理。 　　2 KSVD字典 　　设置Ψ作为初始化字典，设置J作为训练次数，设置Y作为训练集，图像中的各块表示为，对于任意的i，j，KSVD字典的训练方法就可以表述为：首先，执行步骤（1），（2）J次。 　　（1）利用下面的式子（1）并利用OMP算法进行求解各个图像块在字典Ψ上的稀疏表示系数； 　　对于任意的i，j ； 　　0 st.||≤（Cσ）2 （1） 　　式子中：C为噪声增益；σ为误差。 　　（2）原子字典需要进行更新，需要设定字典中原子的个数，假如其中有K个原子。 　　第一步：寻找原子图像块，并找到第一个使用的集合，l=1，2，……，k； 　　第二步：进行计算并对误差进行表示： 　　（2） 　　其中， 原子进行图像表示的误差。 　　第三步：对E1进行设置，它的列为{}（i，j）∈W1 　　第四步：将E1进行SVD处理： 　　2.1 细胞图像边缘的检测方法 　　从细胞图像中能够很明显的看出：图像的边缘含有较多的信息，毫无疑问，图像的边缘必定成为图像重构的重点关注对象，对于图像边缘所进行的各种操作都成为影响最后重构效果的重要因素。假如一个图像的像素点坐标为（x，y），其对应的灰度值为g（x，y），那么其相应的梯度计算为： 　　g= （3） 　　|g|= （4） 　　β= （5） 　　其中，点（x，y）处灰度变化最快的方向用β来表示，最大变化率的数值是用|g|来表示的。对图像进行灰度跃变检测就是常说的边缘检测，并将图像进行转化，将其变成二值图像，经常使用的算法有：Canny，Sobel以及Robert。它主要是在Canny准则的基础进行了一定的发展而形成的。著名的Canny准则为： 　　（1）较小的边缘检测误差； 　　（2）边缘定位很准确； 　　（3）较低的边缘响应次数； 　　Canny算子处理过程为：先进行滤波处理，然后求导，最后求导，获得在最大梯度的边缘检测： 　　（1）为了减少噪声的干扰，对图像数据进行了二维高斯滤波处理； 　　（2）计算灰度值为g（x，y）像素点的梯度； 　　（3）利用上的（4）式子进行梯度方向的求解； 　　（4）将求解的方向梯度进行分类，分成4个种类，并找出其与灰度值最为接近的邻接像素点； 　　（5）用最大值抑制的方法进行边缘像素点的判断，假如是边缘点就置1，反之，则为0； 　　受噪声影响，可能出现斑纹现象。但是，使用Canny算子就不会出现上述情况。因为其使用的是双阈值（N1和N2，并且N1 N2），如大于N2肯定是边缘，小于N1，肯定不是边缘；介于两者的中间，就进一步考察像素点是否有比N1高且相邻的像素，有的话，那就是边缘点，反之，则不是。 　　为了获得更精确的边缘点需要对图像边缘进行细化。 　　2.2 图像中的非边缘点 　　图像中的其他点主要是比较平滑的，并且在整图中所占的比例是很大的，假如对其采用高采样频率，就会浪费大量的数据，所以该章采用低频率。要尽可能多的获取有用的信息，需要保证选取像素点的位置是随机的，随机是对于全图像来说的，并不是BCS[2]中的K个图像块。 　　通过实验发现：由于对非边缘像素点的选择是随机性的，没有办法保证其位置，及它的数量，这就加大了