原标题:Advanced Photonics|深度神经网络实现透过厚散射介质成像

现有的解决光散射后成像的方法主要是采用相干光成像,该方法由于噪声高导致对比度很低,不能得到清晰的图像。近日,中国科学院上海光学精密机械研究所司徒国海研究员课题组提出了一种混合神经网络,成功实现了通过厚散射介质成像。

背景

当光在厚的散射介质(如皮肤等生物组织、浓雾等自然环境)中传播时会发生许多次散射,波前会受到许多次随机调制,导致无法形成清晰的像。针对这一问题,传统思路是通过门控技术选通只经过很少次散射的弹道光,如基于光学相干成像(OCT)即是利用相干门选通弹道光的方法。但弹道光的功率随着光在散射介质中的传播而呈指数衰减(衰减到1/e为一个光学厚度)。虽然在一般条件下弹道光能成像,但由于在“门”打开的短暂时间内仍然有大量的散射光到达像面,导致所成的像对比度极低,甚至会导致及其微弱的弹道光淹没在散射光的“噪声”之中。最近,随着计算光学成像技术的发展,研究人员提出了多种利用散射光而不是弹道光的成像方法,如传输矩阵法、相位共轭法、波前整形法和记忆效应法等等。但这些方法都各有优缺点和适用条件。例如,基于散射光关联特性的记忆效应法在面对厚的散射介质时会失效;同时,由于它是一种介观效应,当被观察物体尺寸较大时也会失效。

从记录到的散斑恢复物体图像是一个欠定的逆问题,尤其在被厚散射体散射的情况下,弹道光的信号极弱而散射“噪声”很强,导致经典的解逆问题的最优化算法都无能为力。中国科学院上海光学精密机械研究所司徒国海研究员课题组发现以深度学习为代表的人工智能技术可以用来解决这一难题。

深度学习方法通过对一系列原始数据的学习,自动提取数据中的抽象信息,并将抽象信息进行组合得到复杂抽象的特征,再使用这些特征解决最终的问题。这个方法已经广泛地应用在计算机视觉、图像处理、语音识别、自动驾驶、目标识别等领域,并取得了十分瞩目的成果。

该课题组于2017年提出了将深度神经网络用于解决透过厚散射体成像的想法。此后,MIT的George Barbastathis教授课题组提出了利用卷积神经网络通过薄的毛玻璃恢复物体的相位,Boston大学的Lei Tian教授课题组利用通过不同毛玻璃之后采集到的散斑图对卷积神经网络进行训练,提高了神经网络的稳健性,从而能够在一定程度上实现通过不同的毛玻璃成像。

最近,该课题组在2017年工作的基础上,针对卷积神经网络对通过强散射介质之后从散斑图中恢复物体时拟合能力不强等问题,提出了一种混合神经网络,成功实现了通过厚散射介质成像。

实验原理及结果

课题组成功实现了用波长为532 nm的连续激光透过一块光学厚度达13.4的聚苯乙烯板获得成像。该方法将卷积神经网络模型与全连接网络模型相结合:首先用2个卷积层提取散斑图中的结构信息,随后用3个全连接层增加网络的拟合能力,最后用3个卷积层对图像进行恢复重建。这样的方法既提高了神经网络的拟合能力,又尽量避免了全连接网络中物体结构信息丢失等问题,成功地实现了通过厚散射介质成像。在实验中,3990张MNIST手写字符及其对应的散斑图被用于训练网络,训练网络的时间大约为194 s。网络训练完成后,20张不包括在训练图像集合中的MNIST手写字符所对应的散斑图被用来测试网络,结果表明,从每张散斑图恢复出物体的时间大约为0.78 s。

实验光路如图1所示,被空间光调制器(SLM)反射的激光通过厚的散射介质之后被相机接收,图1(b)与图1(c)分别表示通过散射介质之前与之后采集到的物体图片。将采集到的散斑图与加载在SLM上的图片输入神经网络进行训练,网络训练完成之后就可以利用深度神经网络实现通过散射介质成像。实验恢复结果如图2所示。

图1 (a)利用深度神经网络实现通过厚散射介质成像实验光路示意图,(b)通过散射介质之前图像,(c)通过散射介质之后采集到的散斑图,(d)散射介质聚苯乙烯板,(e)散射介质侧视图。

图2 利用深度神经网络恢复的结果。(a)散斑图,(b)利用深度神经网络实现的恢复结果,(c)物体图片,(d)利用记忆效应恢复的结果。

从图2(b)的结果可以看出,利用深度神经网络可以实现记忆效应(图2(d))无法实现的效果,这证明了深度神经网络通过厚散射介质成像的有效性。

此外,实验发现,利用基于深度神经网络的通过散射介质成像技术,可以实现仅仅利用原始散斑图0.1%的信息对物体进行恢复。实验中,原始采集到的散斑图为512 pixel×512 pixel、16位深,而利用深度神经网络,可以利用从采集到的散斑图中选取的任意64 pixel×64 pixel二值化图像实现物体成像恢复,相关的实验结果如图3所示。

图3 利用深度神经网络实现图像压缩,第一行为采集到的散斑图,第二行为利用深度学习恢复结果,第三行为物体图片。(G1)为灰度区间为[25000,30000]的散斑图,(G2)为灰度区间为[35000,40000]的散斑图,(G3)为灰度区间为[40000,45000]的散斑图,(G4)为以35000为阈值的二值化的散斑图。

结论

由此可见,利用深度学习方法可以实现通过厚的光学散射介质成像,并可以减少恢复物体需要的散斑图的信息量:仅仅利用散斑图0.1%的信息就能够实现物体图像的恢复。这种基于深度学习的散射成像技术今后有望在通过生物组织成像等领域得到广泛应用,其对于散斑图像冗余性的讨论也对更深入理解散斑图像所包含的信息具有一定的启发意义。

End

如需转载,请直接留言。

免责声明

本文注明来源为其他媒体或网站的文/图等稿件均为转载,如涉及版权等问题,请作者在20个工作日之内来电或来函联系,我们将协调给予处理(按照法规支付稿费或删除)。返回搜狐,查看更多

责任编辑:

透过散射介质的成像matlab,Advanced Photonics|深度神经网络实现透过厚散射介质成像...相关推荐

  1. 透过散射介质的成像matlab,利用空间光调制器实现基于相位多样性的散射介质成像的制作方法...

    本发明属于散射介质成像领域,具体说是一种利用空间光调制器实现基于相位多样性的散射介质成像方法. 背景技术: 传统的光学成像方法通常无法直接获得隐藏在散射介质后方的物体像,因此如何利用光学技术实现穿透散 ...

  2. 基于学习的透过厚散射介质的无透镜光学成像

    基于学习的透过厚散射介质的无透镜光学成像 Learning-based lensless imaging through optically thick scattering media 文献引文信息 ...

  3. SAR成像系列:【8】合成孔径雷达(SAR)成像算法-压缩感知(Compressed Sensing,CS)成像算法(附Matlab代码)

    压缩感知(Compressed Sensing,CS)该理论指出:对于满足约束等距条件(Restricted Isometry Property,RIP)的稀疏或可压缩信号,通过低于(甚至远低于)Ny ...

  4. x射线 相衬成像 matlab模拟,X射线相衬成像系统与成像方法

    X射线相衬成像系统与成像方法 [技术领域] [0001] 本发明涉及X射线光栅成像技术,特别涉及基于分布式X射线源的X射线相衬成 像系统与成像方法. [背景技术] [0002] 在现有技术例如CT扫描 ...

  5. SAR成像系列:【7】合成孔径雷达(SAR)成像算法-后向投影(Back Projecting)算法(附Matlab代码)

    前面介绍了SAR成像的RD算法和CS算法,接下来介绍两种时域成像算法,其一就是后向投影(BP)算法. BP成像的优点:成像算法简单,鲁棒性好,分辨率高,适用于任何轨道或飞行轨迹模型,不存在斜距近似假设 ...

  6. SAR成像系列:【9】合成孔径雷达(SAR)成像算法-波数域(omega-K)成像算法[也叫距离徙动(RM)算法](附Matlab代码)

    波数域()成像算法作为本系列的最后一种成像算法介绍.关于SAR成像的其他的各种改进算法就不一一列举了.在实际成像中,万变不离其踪,最主要的是关注成像的几何模型,再根据指标选择不同的基础成像算法,然后进 ...

  7. matlab 合成孔径,逆合成孔径成像MATLAB算法)配书源代码

    [实例简介] Inverse Synthetic Aperture Radar Imaging With MATLAB Algorithms(逆合成孔径成像MATLAB算法)配书源代码,可用于学习和理 ...

  8. 合成孔径雷达成像matlab,合成孔径雷达成像Matlab仿真研究.doc

    合成孔径雷达成像Matlab仿真研究.doc 合成孔径雷达成像Matlab仿真研究 摘 要: 计算机仿真是现代雷达研究中的重要技术之一,针对合成孔径雷达(SAR)成像中影响仿真结果的因素,从分析SAR ...

  9. mimo 雷达成像 matlab,MIMO雷达成像算法研究

    摘要: 作为一种新兴雷达技术,多输入多输出(MIMO)雷达受到国内外科研人员的广泛关注.该雷达综合利用阵列与分集技术,能够引入远多于实际物理阵元数目的观测通道和自由度,并在信号层联合进行多通道回波数据 ...

最新文章

  1. TiDB 在大型互联网的深度实践及应用--大数据平台技术栈08
  2. 量子信息技术研究现状与未来——郭光灿
  3. 对软件工程与计算机科学之间区别的看法
  4. 【Python】青少年蓝桥杯_每日一题_4.03_求偶数
  5. js模拟form表单提交数据, js模拟a标签点击跳转,避开使用window.open引起来的浏览器阻止问题...
  6. Python面试题总结(4)--数据类型(列表)
  7. 中国搪瓷板(VE+Panel)市场趋势报告、技术动态创新及市场预测
  8. 速修复!NSA 报告四个严重和高危 Exchange Server RCE 漏洞
  9. iphone 开发第四天 - 字符串
  10. 7时过2小时是几时_小姐姐花1万7租了一个2小时男友,全程上演偶像剧......
  11. 取消计算机硬盘自检,教你电脑开机自检如何关闭
  12. 一般时序电路设计(次态卡诺图化简法)
  13. 黑马程序员rocketmq第一章
  14. java nio rewind_java.nio.ByteBuffer中的flip()、rewind()、compact()等方法的使用和区别
  15. 基于K均值算法的鸢尾花聚类实验(Sklearn实现)
  16. Java 进制转换工具类
  17. FLV Binder下载
  18. C++ 的静态成员变量为什么一定要在类外定义
  19. html制作百度页面
  20. 基于C的双声道PCM文件数据处理

热门文章

  1. spring中自带的缓存springcache
  2. 对于vue项目中报Disallow self-closing on HTML void elements格式错误时的解决方案
  3. java中GUI中显示当前时间_【java编程】返回当地时间Clock到GUI面板Panel上
  4. SYBASE 12.5 维护手册 3. 技术专刊
  5. 【R语言】ggplot2作图补充(1)
  6. 2021年,我们还需要入门深度学习吗?
  7. uboot环境下mmc操作_uboot mmc命令详解
  8. 2-- 关机/重启 、 创建 / 切换目录
  9. Java8新特性(一)—————Lambda表达式
  10. 计算机二级心得体会论文,学习计算机二级心得体会.docx