压缩感知(Compressed sensing),也被称为压缩采样(Compressive sampling)或稀疏采样(Sparse sampling),是一种寻找欠定线性系统的稀疏解的技术。压缩感知被应用于电子工程尤其是信号处理中,用于获取和重构稀疏或可压缩的信号。這個方法用到訊號稀疏的特性,得以從相對較少的測量值還原出原來整個欲得知的訊號。MRI就是一個可能使用此方法的應用。这一方法至少已经存在了四十年,由于David Donoho、Emmanuel Candès和陶哲轩的工作,最近这个领域有了长足的发展。

目录

  • 1 概览
  • 2 取样方法
    • 2.1 结构简化采样矩阵
    • 2.2 数值简化采样矩阵
  • 3 欠定线性系统
  • 4 求解/重构方法
    • 4.1 基追踪
    • 4.2 匹配追踪
  • 5 应用
    • 5.1 摄影学
    • 5.2 全息成像
    • 5.3 面部识别
    • 5.4 核磁共振成像
    • 5.5 网络诊断
    • 5.6 短红外摄影
  • 6 參考資料
  • 7 外部链接

概览

信号处理领域中的一个常见问题就是从一系列的采样中重建原本的信号。大致而言,由于不可能重建出未被采样的部分信号,这一任务是无法完全实现的。然而通过借助对于信号(性质)的预先了解或合理假设,完美地通过一系列采样重建原信号就成为了可能。随着科学的发展,数学家们逐步增进了如何作出合理假设的能力,并慢慢了解到在何种情况下可将这些假设一般化、推广化。

信号处理领域中的一次较早的突破是采样定理的提出。这一定理证明了若信号的最高频率小于采样频率的一半,便可完美地从采样结果中恢复原本信号。

在2004年左右,Emmanuel Candès、陶哲轩和David Donoho证明了在已知信号稀疏性的情况下,可能凭借较采样定理所规定更少的采样数重建原信号。这一理论也是压缩感知的基石。

取样方法

在理论上,为了确保信号重建的准确度,需要令所采用的取样矩阵各行列之间相干性([[:en:Coherence_(signal_processing)}]])尽量地低,且须矩阵元素(element)取值随机性尽量地高。

考虑到以上原因,最为标准的做法是采用相同独立分布随机高斯矩阵(identical independent distributed random Gaussian matrix)对待处理信号进行取样,即可确保在信号具有足够稀疏性的前提下得到较佳的压缩感知重建效果。

但是在采用相同独立分布随机高斯矩阵时,所需要的数据存储量过于庞大——每个矩阵元素都要单独记录,且数据类型一般为浮点数——这就促进了简化取样矩阵的研究;目前被提出的简化取样矩阵主要包括两种:结构简化采样矩阵与数值简化采样矩阵。

结构简化采样矩阵

一种基本的结构简化采样矩阵是Toeplitz矩阵:

其中

在采用Toeplitz矩阵进行压缩感知取样时,便可以只在系统内存中储存第一行与第一列的矩阵元素,大大降低了储存成本,为算法的硬件实现降低了门槛。

数值简化采样矩阵

数值简化采样矩阵普遍将原有的、按照高斯分布随机取值的采样矩阵元素以数值上更为简单的元素取代,但在分布上维持矩阵元素的分布随机性——即缩减了储存浮点数这一方面的成本。

一种较为基本的数值简化采样矩阵是0-1伯努利矩阵,其中的元素仅有0和1两种,分布模式为相同独立的伯努利分布(identical independent Bernoulli distribution)。

对于每一个矩阵元素,该分布的概率质量函数为:

欠定线性系统

如果一个线性方程组未知数的数目超过方程的数目,这个方程组被称为欠定,并且一般而言有无数个解。但是,如果这个欠定系统只有唯一一个稀疏解,那么我们可以利用压缩感知理论和方法来寻找这个解。值得注意的是,不是所有欠定线性方程组都有稀疏解。

求解/重构方法

压缩感知利用了很多信号中所存在的冗余(换言之,这些信号并非完全是噪声)。具体而言,很多信号都是“稀疏”的;在适当的表示域中,它们的很多系数都是等于或约等于零。

在信号获取阶段,压缩感知在信号并不稀疏的域对信号进行线性测量。

为了从线性测量中重构出原来的信号,压缩感知求解一个称为L1-范数正则化的最小二乘问题。从理论上可以证明,在某些条件下,这个正则化最小二乘问题的解正是原欠定线性系统的稀疏解。

基追踪

基追踪(basis pursuit[1])是一种信号重建演算法,被广泛地用于压缩感知领域,属于数学最优化问题的范畴,形式为 。,

其中xN×1尺寸向量,表示输出(信号),yM×1的测量向量,AM × N的“转换矩阵”或曰“测量矩阵”;其中M < N

基追踪常在需要完美满足欠定线性方程组系统中y = Ax时被应用,且要求解在L1基准下为最稀疏的。

若应用情景允许降低对完美恢复的要求,以换取更加稀疏的解x,降噪基追踪(basis pursuit denoising[2])更为适用。

匹配追踪

匹配追踪(Matching pursuit[3])是一种稀疏近似运算,旨在找到多维数据在某个超完备字典(dictionary)上映射的“最佳匹配”。其基本的概念就是要用来自的函数组(称为基元,或atom)的加权和来表示希尔伯特空间上的信号

其中表示被选取基元的序数,是每个基元的加权常数。对于固定的字典,匹配追踪会先找到与信号间内积最大的一个基元,再减去该基元带来的影响,然后重复找寻影响力次大的基元直到信号被较好地分解。

相较而言,以傅里叶级数表示的信号来说,字典是构建于正弦基函数之上的。信号处理领域中傅里叶分析的主要缺点就在于它只能提取出信号中常存的特征,而不能适应当前的分析目标信号

若采用一组极端保险、带有大量冗余的字典,我们就能够找到可以完美匹配信号的函数组。在对信号进行编码与压缩时,最好是能够找到一组映射,使加权和中的大部分系数都接近零(具有“稀疏性”)。

应用

压缩感知与包括欠定系统、群验、稀疏编码、复用、稀疏采样等。在成像科技领域,亦有许多技术如(编码孔与计算摄影学)与压缩感知相关。亦有许多在不同技术完成度下将压缩感知实现的案例。

摄影学

压缩感知技术被用于手机摄像头传感器设计中。这一技术使得在获取图像时所耗费的能量大大降低,可达原先耗费的15分之一——当然,需要引入复杂的解压算法;这一运算可能会需要脱机状态下的预先安装、设置。

压缩感知亦被用于实现单像素摄影。贝尔实验室运用了这一技术,用无镜片单像素相机在栅格中随机选取的孔隙处拍照,以达到成像效果。随着拍照次数的逐渐增多,照片质量也会逐步提高;一般来说,这种技术较之传统的摄影成像技术仅仅需要一小部分的数据占用量,且能完全避免与镜片或聚焦相关的故障或失常;参见[4] 。

全息成像

压缩感知可被用于增加通过单幅全息图中所能得到的立体像素的数量改进全息摄影技术中的图像重建问题。在光学全息图或毫米波全息图采样率不足的情况下,这一技术也能够被用于图像检索以作出改善。

面部识别

压缩感知目前被用于面部识别应用之中,参加Engineers Test Highly Accurate Face Recognition。

核磁共振成像

压缩感知被用于缩短核磁共振成像中的扫描时间,参加[5] 、[6] 与[7];其中涉及的重建方法包括ISTA、FISTA、SISTA等。

网络诊断

压缩感知在被用于旨在利于网络管理的网络诊断应用中时带来了极佳成效。对网络延时的估计和网络拥塞的探知均可被归纳、建模为非定性的线性方程组系统,其中的参数矩阵正是所分析网络的路由选择矩阵。此外,在互联网情景中,网络路由矩阵常常能够满足压缩感知技术所要求的几个基本要素:低相关性、稀疏性及(可能的)R.I.P条件,请参见[8] 。

短红外摄影

目前,基于压缩感知技术的商用短红外相机已被推出[9] 。这些相机的光敏度大约从0.9µm到1.7µm,在上述波段上,人类的肉眼是无效的。

參考資料

  • Candès, E.J., & Wakin, M.B., An Introduction To Compressive Sampling, IEEE Signal Processing Magazine, V.21, March 2008

外部链接

  • "The Fundamentals of Compressive Sensing" Part 1, Part 2 and Part 3: video tutorial by Mark Davenport, Georgia Tech. at SigView, the IEEE Signal Processing Society Tutorial Library.
  • Using Math to Turn Lo-Res Datasets Into Hi-Res Samples Wired Magazine article
  • Compressive Sensing Resources at Rice University.
  • Compressed Sensing Makes Every Pixel Count – article in the AMS What's Happening in the Mathematical Sciences series
  • Wiki on sparse reconstruction

压缩感知(Compressed sensing)from wiki相关推荐

  1. matlab压缩感知,压缩感知(Compressed Sensing, CS)matlab代码

    [实例简介] 压缩感知(Compressed Sensing, CS)matlab代码.实现多个正弦信号的随机欠采样,通过压缩感知恢复.两个m文件分别是两个算法,正交匹配追踪(OMP)算法和SPGL1 ...

  2. 压缩感知 (Compressive sensing)之公式及关键点说明

    前言 压缩感知(Compressive sensing)和压缩采样(Compressive sampling)是一样的. 本想细更一下这个内容,但细看完"参考资料1"后,我决定放弃 ...

  3. 压缩感知 compressive sensing(转)

    信号处理新理论 采样定理(又称取样定理.抽样定理)是采样带限信号过 程所遵循的规律, 1928年由美国电信工程师H.奈奎斯特首先提出来的,因此称为奈奎斯特采样定理.1948年信息论的创始人C.E.香农 ...

  4. 初识压缩感知 compressive sensing

    压缩感知是近年来极为热门的研究前沿,在若干应用领域中都引起瞩目.最近粗浅地看了这方面一些研究,对于Compressive Sensing有了初步理解,在此分享一些资料与精华.本文针对陶哲轩和Emman ...

  5. 初识压缩感知Compressive Sensing

    压缩感知是近年来极为热门的研究前沿,在若干应用领域中都引起瞩目.最近粗浅地看了这方面一些研究,对于Compressive Sensing有了初步理解,在此分享一些资料与精华.本文针对陶哲轩和Emman ...

  6. 科学松鼠会压缩感知科普文章两篇:“压缩感知与单像素相机(陶哲轩)”“填补空白:用数学方法将低分辨率图像变成高分辨率图像(Jordan Ellenberg)

    题目:科学松鼠会压缩感知科普文章两篇:"压缩感知与单像素相机(陶哲轩)""填补空白:用数学方法将低分辨率图像变成高分辨率图像(Jordan Ellenberg)" ...

  7. 压缩感知算法matlab,压缩感知算法matlab

    MATLAB代码如下: 压缩感知仿真实例 压缩感知仿真实例在MATLAB R2001b中的计算结果如下: 原图像 采样率0.7 采样率0.5 采样率0.3 压缩感知仿真实例采用均方误差...... M ...

  8. matlab21世纪论坛,compressive sensing 压缩感知(转) 21世纪最火的研究方向

    转载于http://to-cs.blog.sohu.com/110268260.html 整个人类匍匐在上帝的脚下,艰难并努力地生存着.有多少次,被上帝无情的戏弄:也有多少次,尝试着和上帝较量.每当想 ...

  9. 压缩感知进阶——有关稀疏矩阵

    上一篇<初识压缩感知Compressive Sensing>中我们已经讲过了压缩感知的作用和基本想法,涉及的领域,本文通过学习陶哲轩对compressive sensing(CS)的课程, ...

  10. 压缩感知先进——关于稀疏矩阵

    前<初识压缩感知Compressive Sensing>中我们已经讲过了压缩感知的作用和基本想法,涉及的领域,本文通过学习陶哲轩对compressive sensing(CS)的课程,对压 ...

最新文章

  1. 知识图谱学习笔记(1)
  2. jQuery笔记---选择器
  3. Change FZU - 2277(线段树+dfs序)
  4. 真想学习?扔掉你详细的计划吧!
  5. linux自带浏览器如何更新,在deepin上安装edge浏览器后更新系统报错的解决方案
  6. php 文章浏览量 缓存,WordPress缓存文章浏览量访问不自动增加怎么办?WordPress缓存导致文章阅读数点赞数不更新...
  7. 每天进步一点点《ML - 异常点检测》
  8. Video-Swin-Transformer 跑通遇到的问题
  9. iOS 使用NJKWebViewProgress做webview进度条(加载网页时获取加载进度)
  10. 下载mysql的驱动包
  11. Step7 V5.6 win10版本安装包下载
  12. 热敏打印机ESCPOS指令封装类,票据打印从此轻松,佳博系列实测通过
  13. 更强大、更灵活、更全面丨一文搞懂DolphinDB窗口计算
  14. 一个小白的自渡-Git 仓库基础操作
  15. 小程序推送代码到远程库
  16. 微信投放卡券-支付即会员
  17. 建模方法(十)-灰色预测模型GM(1,1)
  18. MySQL第一讲 一遍让你彻底掌握MVCC多版本并发控制机制原理
  19. 5800计算机交点法程序大全,FX5800计算器交点法正反算程序加入高程及隧道超欠挖...
  20. AM335X BeagleBone 之格式化SD卡

热门文章

  1. 马化腾:互联网企业竞争要注重什么
  2. web安全day37:Linux脚本判断和循环,编写脚本实现内网主机存活批量检测
  3. OSPFv3中LSA详解(一)——概述
  4. 计算机网络误区——为什么我的PING包发不出去?
  5. MySQL 第五次练习(MySQL备份恢复)
  6. Linux之LVM(逻辑卷管理)、分层存储Stratis、VDO、SWAP分区及相应案例
  7. numpy 中的nan和常用的统计方法
  8. ko.js循环绑定值问题(工作遇见)
  9. 写表单验证等页面的总结
  10. 亚马逊EC2服务器使用Rsync+Inotify实时同步