发明名称 基于分层高斯混合模型的统计压缩感知图像重构方法
摘要 本发明公开了一种基于分层高斯混合模型的统计压缩感知图像重构方法,具体涉及一种基于图像块数据的统计模型学和信号重构算法,主要解决现有统计压缩感知图像重构方法中使用的单高斯先验模型不能灵活准确地刻画子图像块的非高斯统计性质,从而导致重构图像质量不高的缺点。本发明在模型第一层学全局高斯混合模型对子图像块进行硬聚类,从全局上充分利用图像块之间的相似性结构,在模型第二层分类学子图像块的局部高斯混合模型,从局部上对子图像块之间的差异性进行区分和建模。与现有传统压缩感知和统计压缩感知重构技术相比,能够以更少的测量获得更高的重构精度,适合对自然图像的重构。
申请公布号 CN105787895A 申请公布日期 2016.07.20
申请号 CN201610110804.6 申请日期 2016.02.29
申请人 中国计量学院 发明人 武娇;曹飞龙
分类号 G06T5/00(2006.01)I;G06T5/50(2006.01)I;G06T3/40(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 杭州浙科专利事务所(普通合伙) 33213 代理人 杜立
主权项 基于分层高斯混合模型的统计压缩感知图像重构方法,包括如下步骤:(1) 将图像分成<img file="dest_path_image001.GIF" wi="36" he="16" />的无重叠的子图像块,对每个子图像块以<img file="dest_path_image002.GIF" wi="35" he="20" />的采样率进行压缩采样,得到测量<img file="dest_path_image003.GIF" wi="87" he="20" />:<img file="dest_path_image004.GIF" wi="69" he="17" />,其中<img file="dest_path_image005.GIF" wi="18" he="14" />是第<img file="dest_path_image006.GIF" wi="11" he="14" />个子图像块像素值,<img file="dest_path_image007.GIF" wi="13" he="17" />是<img file="dest_path_image008.GIF" wi="46" he="18" />维高斯随机矩阵,<img file="dest_path_image009.GIF" wi="51" he="18" />;(2) 由黑白边缘图像生成<img file="dest_path_image010.GIF" wi="42" he="16" />个初始的方向高斯分布<img file="dest_path_image011.GIF" wi="112" he="22" />,利用DCT变换生成第<img file="dest_path_image012.GIF" wi="16" he="17" />个方向高斯分布<img file="dest_path_image013.GIF" wi="80" he="20" />,<img file="dest_path_image014.GIF" wi="17" he="14" />的初始值设定为<img file="246975dest_path_image015.GIF" wi="9" he="17" />,<img file="632957dest_path_image012.GIF" wi="16" he="17" />设定为19;(3) 在混合模型第一层,由测量<img file="dest_path_image016.GIF" wi="18" he="14" />计算子图像块<img file="740590dest_path_image005.GIF" wi="18" he="14" />在<img file="604641dest_path_image012.GIF" wi="16" he="17" />个高斯模型下的MAP估计和对数后验概率:<img file="dest_path_image017.GIF" wi="239" he="22" />,<img file="dest_path_image018.GIF" wi="138" he="19" /><img file="dest_path_image019.GIF" wi="166" he="21" />,<img file="dest_path_image020.GIF" wi="186" he="17" />比较<img file="dest_path_image021.GIF" wi="16" he="21" />将子图像块聚成<img file="638587dest_path_image012.GIF" wi="16" he="17" />类:<img file="dest_path_image022.GIF" wi="134" he="30" />,第<img file="523367dest_path_image006.GIF" wi="11" he="14" />个子图像块被分到第<img file="dest_path_image023.GIF" wi="28" he="24" />类,<img file="dest_path_image024.GIF" wi="121" he="24" />;(4) 在混合模型第二层,对每个聚类中的子图像块,选择具有最大对数后验概率的<img file="dest_path_image025.GIF" wi="11" he="16" />个高斯分布的混合模型作为先验分布,利用EM算法从压缩测量迭代地学习高斯混合模型,并重构属于每个聚类的子图像块<img file="dest_path_image026.GIF" wi="59" he="20" />,其中<img file="dest_path_image027.GIF" wi="16" he="18" />是第<img file="dest_path_image028.GIF" wi="9" he="18" />个聚类中子图像块的标识集,<img file="478816dest_path_image029.GIF" wi="88" he="17" />;(5) 使用重构的子图像块<img file="dest_path_image030.GIF" wi="59" he="20" />更新<img file="5612dest_path_image012.GIF" wi="16" he="17" />个方向高斯分布的均值向量和协方差矩阵:<img file="dest_path_image031.GIF" wi="108" he="33" />,<img file="dest_path_image032.GIF" wi="272" he="34" />,<img file="205649dest_path_image029.GIF" wi="88" he="17" />其中<img file="dest_path_image033.GIF" wi="16" he="25" />是正则化参数,设定为<img file="dest_path_image034.GIF" wi="31" he="19" />;(6) 如果未达到指定的迭代重构次数<img file="dest_path_image035.GIF" wi="46" he="18" />,返回步骤(3);否则,将重构的子图像块拼在一起,得到原始图像的重构图像。
地址 310018 浙江省杭州市下沙高教园区学源街258号