[发明专利]一种基于层次模型与导向先验的单像元相机快速成像方法

说明书

技术领域

本发明属于数字图像处理领域，特别涉及一种基于层次模型与导向先验的单像元相机快速成像方法。

背景技术

随着压缩感知理论的日益成熟，该理论受到越来越多的关注，应用也越来越广泛，典型的应用涉及图像去噪、图像去模糊、单像元相机成像技术，核磁共振成像技术、雷达成像技术等等。但是由于对图像分辨率的要求越来越高，压缩感知成像速度问题日渐突出。

在压缩感知理论框架下，提高图像成像速度有两种方式，一种是优化最终的图像复原算法，另一种则是通过优化测量矩阵的生成，降低计算量从而进行加速。针对第二种思路，现今的大部分方法采用随机方式来产生测量矩阵。随机方式的优点是容易实现，但随着需要的图像的分辨率的提高，计算量大幅增长，导致成像效率显著下降，一定程度上限制了单像元成像技术的应用。

发明内容

针对上述计算量大幅增长，成像效率显著下降技术问题，本发明提供了一种基于层次模型与导向先验的单像元相机快速成像方法，包括如下步骤：

(S1)搭建单像元成像系统，采集成像数据；

(S2)确定最终需要得到的图像的分辨率为M*N，其中，M为所需图像像素的行数，N为所需图像像素的列数；确定起始层次的图像分辨率为m*n，其中m为起始图像像素的行数，n为起始图像像素的列数；M,N,m,n均为整数；

(S3)记起始层次的图像为第0层，测量矩阵采用随机编码的形式，测量矩阵分别率为m*n；通过成像系统采集数据，采用复原算法，复原得到一个m*n分辨率的图像l₀，并提取其导向先验信息G₀；初始化循环次数k＝1；

(S4)进行第k层的图像复原，先计算第k层与第k-1层之间的距离d_k；然后，依据第k-1层的导向先验信息G_k-1，约束生成测量矩阵，利用成像系统采集的成像数据，采用复原算法，复原得到第k层的图像l_k，同时提取导向先验信息G_k；

(S5)判断图像l_k的分辨率是否大于M*N，若是，则输出图像，否则，k增加1，继续步骤S4。

进一步地，所述提取导向先验信息的步骤为：

对当前层次的复原图像进行截取操作，即：设置一个阈值thr，图像像素值低于该阈值的位置数值变0，大于或等于该阈值的位置数值变1，将复原图像变成二值图像；对得到的二值图像进行逐行扫描，检测每个像素点的八个领域，若该点数值为1，且八个领域中的1的个数小于3，则认为该点为噪点，数值重置为0，否则，该点数值不变化；对于二值图像边缘的像素点，不作考虑，均设置为0；经过处理最终得到的二值图像，即为提取的当前层次复原图像的导向先验信息。

进一步地，所述层间距离d_k计算的方法包括：

计算一个图像信息因子F，F＝prior(L_k-1)/L_k-1，其中，L_k-1为第k-1层复原得到的图像的分辨率，即为M_k-1*N_k-1，M_k-1为第k-1层图像像素的行数，N_k-1为第k-1层图像像素的列数，prior(L_k-1)为提取的图像的导向先验信息中数值1的个数。

层间距离d_k为图像信息因子的倒数F，即d_k＝1/F。

进一步地，所述依据第k-1层的导向先验信息G_k-1，约束生成测量矩阵，具体步骤为：

首先随机产生一个与k-1层图像分辨率大小一致的0和1的随机矩阵，0与1个数的比例为1:1；然后将得到的随机矩阵与第k-1层提取的导向先验信息G_k-1做按位置与操作，得到一个新的随机编码矩阵，即为生成的测量矩阵。

进一步地，步骤S3与步骤S4中所述的复原算法，本发明采用l1magic中提供的算法。

采用本发明获得的技术效果：本发明在压缩感知的理论框架下(参考文献：D.Donoho,Compressed sensing,IEEE Transactions on Information Theory 52(4)(2006)1289–1306.)，引入了导向先验，并且结合提出的层次模型指导测量矩阵的生成，从而降低单像元成像复原的算法复杂度，加快单像元相机的成像速度。

附图说明

图1是单像元成像系统示意图；