CN117761896A - 一种单次曝光压缩超快相干调制成像装置 - Google Patents
一种单次曝光压缩超快相干调制成像装置 Download PDFInfo
- Publication number
- CN117761896A CN117761896A CN202311844219.5A CN202311844219A CN117761896A CN 117761896 A CN117761896 A CN 117761896A CN 202311844219 A CN202311844219 A CN 202311844219A CN 117761896 A CN117761896 A CN 117761896A
- Authority
- CN
- China
- Prior art keywords
- phase
- light
- camera
- imaging
- ultrafast
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 38
- 230000001427 coherent effect Effects 0.000 title claims abstract description 34
- 230000006835 compression Effects 0.000 title claims abstract description 17
- 238000007906 compression Methods 0.000 title claims abstract description 17
- 238000000034 method Methods 0.000 claims abstract description 39
- 230000008569 process Effects 0.000 claims abstract description 30
- 238000001514 detection method Methods 0.000 claims abstract description 20
- 238000009826 distribution Methods 0.000 claims abstract description 17
- 238000005259 measurement Methods 0.000 claims abstract description 6
- 238000012545 processing Methods 0.000 claims description 21
- 238000005286 illumination Methods 0.000 claims description 20
- 238000013461 design Methods 0.000 claims description 8
- 238000005457 optimization Methods 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 6
- HPTJABJPZMULFH-UHFFFAOYSA-N 12-[(Cyclohexylcarbamoyl)amino]dodecanoic acid Chemical compound OC(=O)CCCCCCCCCCCNC(=O)NC1CCCCC1 HPTJABJPZMULFH-UHFFFAOYSA-N 0.000 claims description 5
- 238000012546 transfer Methods 0.000 claims description 4
- 230000003190 augmentative effect Effects 0.000 claims description 2
- 238000013527 convolutional neural network Methods 0.000 claims description 2
- 238000011084 recovery Methods 0.000 claims description 2
- 230000009467 reduction Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 13
- 238000012634 optical imaging Methods 0.000 abstract description 9
- 230000001052 transient effect Effects 0.000 abstract description 5
- 230000003287 optical effect Effects 0.000 description 14
- 239000000523 sample Substances 0.000 description 11
- 238000002679 ablation Methods 0.000 description 7
- 239000010408 film Substances 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 239000000758 substrate Substances 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 239000010409 thin film Substances 0.000 description 5
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000000608 laser ablation Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000005684 electric field Effects 0.000 description 3
- 238000005530 etching Methods 0.000 description 3
- 239000005350 fused silica glass Substances 0.000 description 3
- 230000001788 irregular Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000000644 propagated effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000001747 exhibiting effect Effects 0.000 description 2
- AMGQUBHHOARCQH-UHFFFAOYSA-N indium;oxotin Chemical compound [In].[Sn]=O AMGQUBHHOARCQH-UHFFFAOYSA-N 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- VYZAMTAEIAYCRO-UHFFFAOYSA-N Chromium Chemical compound [Cr] VYZAMTAEIAYCRO-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000010587 phase diagram Methods 0.000 description 1
- 238000007747 plating Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 239000005361 soda-lime glass Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种单次曝光压缩超快相干调制成像装置。通过结合相干调制成像和压缩超快成像技术，使用相干光照明待测物体，在样品对光波衍射后经过随机相位调制板的相位调制，利用狭缝完全打开的条纹相机在探测平面上偏转记录调制波前的积分强度分布，得到单张压缩衍射图案，再开发基于即插即用框架和深度图像先验的空间二维至时空三维图像重构算法，使用该算法对采集到的图像进行数据恢复。相比于传统的超快光学成像与相位成像技术，本发明可以实现在单次测量过程中定量的获取瞬态场景完整的强度和相位信息，具有单次曝光相位成像下皮秒量级时间分辨率和实时定量强度‑相位成像能力的特点，可广泛应用于各种需要观测相位的超快现象。
Description
技术领域
本发明涉及超快光学成像技术领域，特别是一种单次曝光压缩超快相干调制成像装置(CS-CMUI)。
背景技术
单次曝光超快成像技术是超快光学成像的一个重要分支，其中压缩超快成像(CUP)适用于记录不可重复或不可逆的瞬态事件，并通过图像重建算法从二维图像中恢复原始动态场景。相位作为一个重要的参量，在光学研究中描述了光波的波前形状和光程差，对于光学成像具有重要的影响。相位可以提供光强无法提供的关于物体的折射率、厚度等信息，这对于待测物理现象的理解具有更重要的作用。相干调制成像(CMI)是相位成像中的一种无透镜衍射成像技术，通过在被测样品和图像传感器之间引入一个随机相位或振幅调制器，通过对物波的多次调制来增加约束条件，提升迭代算法收敛速度，可实现从单幅衍射强度图中重建复杂的光场。
然而，超快光学成像为强度记录方式，无法获取物体的相位信息；相位成像技术仅能实现静态记录，无法获取时间分辨信息。目前还没有一种技术可以实现在单次测量过程中定量的获取瞬态场景完整的强度和相位信息。超快相位成像的难点在于，目前的传感器或其他记录介质只能感知样品的强度信息，无法直接获取相位信息，因此光学相位成像需要通过特殊的光学设计和信号处理方法进行，同时相位信息容易受到光路中的折射、散射和吸收等失真的影响，成像质量和应用效果不佳。
发明内容
为了克服传统光学器件无法探测相位信息的难点与突破传统相位成像技术无法获取时间分辨信息的瓶颈，本发明提出了一种单次曝光下同时具有皮秒量级时间分辨和定量强度-相位成像能力的时空三维成像装置，其能够克服当前超快光学成像技术无法实时定量探测相位变化的缺陷，可以同时获取超快瞬态场景的强度和相位随时间的变化，提高超快光学成像的性能，拓展光学成像技术的应用领域。本发明是一种结合相干调制成像(CMI)和压缩超快成像(CUP)的单次曝光压缩超快相干调制成像(CS-CMUI)，使用无透镜设计降低光路结构失真的影响，并基于CMI实现相位和强度信息的定量探测，基于压缩感知原理和条纹相机实现单次曝光下时空全息信息的压缩记录，最后再利用基于即插即用框架(plug and play,PnP)和深度图像先验(deep image prior,DIP)的迭代重构算法PnP-DIP对测量数据进行还原，最终实现具有皮秒量级时间分辨率的相位和强度定量超快成像。
实现本发明目的的具体技术方案是：
一种单次曝光压缩超快相干调制成像装置，特点是所述装置包括探测光相位调制光路、条纹相机与数据处理系统；
所述的探测光相位调制光路包括依次光路连接的相干照明光源、含相位信息的待测物体(可为动态场景)和随机相位调制板；
所述的探测光相位调制光路采用无透镜设计，降低了透镜造成的光路结构相位失真的影响。
所述的条纹相机内部包括依次光路连接的入射狭缝、光阴极、微通道板、扫描电极、荧光屏与内置相机；
所述的数据处理系统由计算机构成，所述计算机搭载计算显卡和CUDA通用并行计算架构；其中：
所述相干照明光源的镜头中心、待测物体的几何中心、随机相位调制板的几何中心与条纹相机的镜头中心在同一条直线上，此直线即为光的传播方向；
所述相干照明光源、含相位信息的待测物体、随机相位调制板和条纹相机分别对应着光路中的四个平面，依次分别为出射平面、入射平面、调制平面和探测平面；四个平面均垂直于光的传播方向，且光路与平面的交点在平面的几何中心。
所述CUDA通用并行计算架构中含一种图像重建算法，其基于即插即用框架和深度图像先验的一种空间二维至时空三维图像重构算法即PnP-DIP算法；数据处理系统在信号采集完成后，PnP-DIP算法对记录到的压缩衍射全息图案进行数据处理，进而重建出被测物体的相位和振幅信息。所述PnP-DIP具体数据处理过程如下：
设待测物体的复振幅为：
U0(x,y,t；0)＝I(x,y,t；0)exp[iφ(x,y,t；0)] (1)
其中i代表虚数单位，x、y是物体的空间二维坐标，t代表时间，分号“；”后的0代表当前的传播距离；I(x,y,t；0)表示物体的振幅信息，φ(x,y,t；0)表示物体的相位信息；经过d1距离的传播后，衍射光的复振幅表示为：
其中代表d1距离的传递函数，k代表波数，λ代表相干照明光源的波长；/>是U0的二维傅里叶变换，fx和fy分别是图像在x和y方向上的空间频率；经过随机相位调制板后的复振幅表示为：
其中m(x,y)代表相位调制板的空间相位分布；UM继续经过d2距离的传播后，衍射光的复振幅表示为：
此时，相机接收到的光强信息为：
在条纹相机内部，时变光强信息经过偏转后被积分叠加，最终得到的压缩结果为：
其中v代表条纹相机的偏转速度。上述所有过程可以用矩阵的形式表示为：
I(x,y)＝TS|Hd2MHd1U(x,y,t)|2 (7)
其中，U(x,y,t)代表待测物体的复振幅，Hd1和Hd2为空间传递算符，M代表相位调制算符，S代表条纹相机的时空偏转算符，T表示条纹相机CMOS在单次曝光过程中的积分算符；
以上正向采集过程进一步的简略表示为：
y＝Ax (8)
其中表示待测物体的复振幅分布，/>表示测量矩阵，表示最终获得的压缩数据结果；这里，变量W、H和B分别表示动态场景的x、y和t方向上的像素数，变量N表示压缩快照坐标系中y方向上的像素数。通常，N的值大于H的值，但比H×B的值小得多；数据恢复的过程即是从已知信号y计算得到待测数据x的过程；从数学角度出发，已知y和A求解原始信号x是一个明显的欠定问题，因为已知数据的数量远小于未知数据的数量，所以要准确恢复信号必须增加其他的正则约束条件，缩小解空间的大小；同时，由于A的矩阵形式过于庞大，求解A-1的过程相当耗费计算资源和时间，在此，通过引入DIP来代替求逆过程，将逆问题改写为如下形式的最优化问题：
其中R(x)是一个先验项，λ为正则化参数，用来平衡保真项和先验项；G(θ,z)代表一个可更新参数的深度卷积神经网络，其输入是一个与原始数据x同等大小的标准高斯分布白噪声，并且在优化过程中保持不变，θ为网络中可更新的权值参数，将G(θ,z)简写为G(θ)；进一步，通过引入一个可更新的辅助参数u和一个手动设置的惩罚参数ρ，并构造增广拉格朗日函数Lρ，将上述带约束的优化问题改写为无约束优化问题：
利用ADMM框架，的最小化问题分解为以下三个子问题:
1)
2)
3)uk+1＝uk+(xk+1-G(θk+1)). (3)
对于x的更新，观察到它是输入G(θk)-uk关于的近端算子的形式，在此使用即插即用设计利用降噪算法代替近端算子，以防止深度图像先验过拟合；对于θ的更新，在每次迭代中采用一个1阶梯度下降步骤；
从式(10)、式(11)、式(12)与式(13)中不难看出，随着迭代次数k的增加，求解对象xk+1会不断趋近于真实待测物体的复振幅分布x，最后在k无限大的时候收敛于x，根据式(8)，x即为待测物体的复振幅信息U0(x,y,t；0)＝I(x,y,t；0)exp[iφ(x,y,t；0)]；通过手动选择合适的迭代次数k，从获得的压缩数据结果重建出被测物体的相位和振幅信息。
本发明的相干照明光源为532纳米调Q激光器，脉冲宽度为10纳秒(半峰全宽)，最大脉冲能量为5毫焦。
所述随机相位调制板在熔融石英衬底上577纳米深度处通过不规则蚀刻形成，单位尺寸为50微米×50微米，在532纳米波长处呈现0和π相调制的二元随机分布。
所述条纹相机是一种时间偏转成像设备，光阴极将入射光转换为与光强度成正比的电子，电子根据光强度的不同在扫描电极的作用下发生位移不同的垂直偏转，从而将光强度转换为电强度。
本发明的光从相干照明光源出发，经过待测物体、随机相位调制板，最终被条纹相机接收，形成一条完整的光路。
所述条纹相机将接收到的电信号输入到数据处理系统中的计算机中。
本发明的PnP-DIP算法利用高性能计算显卡和CUDA通用并行计算架构对训练与推理过程进行加速，显卡的大规模并行计算能力可以有效提升算法运行速度，减少计算时间。
本发明除数据处理系统外其他元件须在光学平台上搭建，同时环境光要求尽可能少，减少对系统探测结果的影响。
所述探测光相位调制光路采用无透镜设计，降低了透镜造成的光路结构相位失真的影响。
本发明的有益效果是，兼具超灵敏相位探测能力和超高时间分辨能力的光学成像技术可以实现对样品在超短时间尺度上的动态过程进行成像和诊断。本发明结合相干调制成像技术和压缩超快成像技术，提出一种新型的单次曝光压缩超快相干调制成像技术，实现对超短时间尺度上的动态过程进行定量强度和相位成像，为探索研究相位型超快瞬态事件的动力学过程提供强有力的技术支撑，有望促进诸多相关研究的进展。
附图说明
图1是本发明结构示意图；
图2为照明激光脉冲通过E形调制图案后得到E形激光脉冲的示意图；
图3为本发明捕获E形纳秒激光脉冲得到的强度图与相位图；
图4为通过本发明实时检测ITO薄膜(铟锡氧化物半导体透明导电膜)激光烧蚀过程的示意图；
图5为本发明检测ITO薄膜激光烧蚀得到的强度图与相位图。
具体实施方式
以下结合附图及实施例对本发明做详细描述。
参阅图1，本发明包括探测光相位调制光路100、条纹相机200与数据处理系统300三个部分；探测光相位调制光路100包括相干照明光源1、含相位信息的待测物体2(可为动态场景)和随机相位调制板3三个部分；条纹相机200内部由入射狭缝4、光阴极5、微通道板6、扫描电极7、荧光屏8与内置相机9六个部分组成；数据处理系统300由计算机10组成。
在系统搭建好之后，首先，使用相干照明光源1的超快相干脉冲光照明待测物体2，基于CMI原理，由于待测物体2对相干光的衍射作用，光束会在传播方向上扩散，传播70mm距离后，被随机相位调制板3调制，形成一个大发散角、颗粒状的衍射光斑，实现对探测器动态响应范围要求的降低；继续传播30mm距离之后，基于压缩超快成像的原理，携带物体复振幅信息的光束进入一台狭缝4完全打开的条纹相机200，条纹相机内部的光阴极5将光信号转换为电子信号通过微通道板6，电子脉冲在传输过程中受到扫描电极7偏转板内线性电场的作用产生偏转，使得不同时间到达的电子在空间上分开，最后被偏转的电子轰击荧光屏8后重新产生光信号并被内置相机图像采集9芯片记录下来，得到一张压缩衍射全息图案，完成正向探测过程。整个光路采用无透镜设计，避免了光学元件本身产生的或制造缺陷引起的波前形变。
参阅图1，待测物体2与随机相位调制板3之间的距离为70mm，随机相位调制板3与条纹相机200入射平面之间的距离为30mm。
参阅图1，相干照明光源1为532纳米调Q激光器，脉冲宽度为10纳秒(半峰全宽)，最大脉冲能量为5毫焦；随机相位调制板3在熔融石英衬底上577纳米深度处通过不规则蚀刻形成，单位尺寸为50微米×50微米，在532纳米波长处呈现0和π相调制的二元随机分布。
参阅图1，为了捕捉完整的超快过程，条纹相机200以一维模式运作，其选择与待测事件相近的时间挡位进行记录。根据需求，可考虑使用线阵超高速相机替代。典型的线阵超高速相机具有百kHz级别的连续曝光速度和微秒级别的时间分辨率。
参阅图1，在条纹相机200记录到压缩衍射全息图案后，使用基于即插即用框架和深度图像先验的PnP-DIP算法对记录到的压缩衍射全息图案进行数据处理，通过对图像配准、变换、裁剪和归一化，进而重建出被测物体的相位和振幅信息。PnP-DIP算法本质上是一种从获得的空间二维数据还原至时空三维图像的重构算法。CS-CMUI装置具有较高的灵活性和可调谐性，可以通过扫描或者多角度成像的方式，实现对物体的三维成像，获得更加全面和准确的信息，从而提高光学系统的相位精度和稳定性。
实施例1
本实施例通过测量调制纳秒激光脉冲验证本发明的强度和相位检测能力。参阅图1及图2，在相干照明光源1和待测物体2之间增加一个“E”形调制图案，该图案由空心E形石英玻璃基板上的镀铬薄膜形成，并通过该图案进行空间调制照明激光脉冲，从而产生波前为E形的激光脉冲。
参阅图1，E形脉冲光照明待测物体2，基于CMI原理，由于待测物体2对相干光的衍射作用，光束会在传播方向上扩散，传播70mm距离后，被随机相位板3调制，形成一个大发散角、颗粒状的衍射光斑，继续传播30mm距离之后，基于压缩超快成像的原理，携带物体复振幅信息的光束进入一台狭缝4完全打开的条纹相机200，条纹相机内部的光阴极5将光信号转换为电子信号通过微通道板6，电子脉冲在传输过程中受到扫描电极7偏转板内线性电场的作用产生偏转，使得不同时间到达的电子在空间上分开，最后被偏转的电子轰击荧光屏8后重新产生光信号并被内置相机图像采集9芯片记录下来，得到一张压缩衍射全息图案，完成正向探测过程。
参阅图1，待测物体2与随机相位调制板3之间的距离为70mm，随机相位调制板3与条纹相机200入射平面之间的距离为30mm。
参阅图1及图2，相干照明光源1为532纳米调Q激光器，脉冲宽度为10纳秒(半峰全宽)，最大脉冲能量为5毫焦；“E”形调制图案的中心点与光束相交，且图案所在平面与光传播方向垂直；随机相位调制板3在熔融石英衬底上577纳米深度处通过不规则蚀刻形成，单位尺寸为50微米×50微米，在532纳米波长处呈现0和π相调制的二元随机分布。
参阅图1，为了捕捉完整的超快过程，条纹相机200以一维模式运作，其选择与待测事件相近的时间挡位进行记录。
参阅图1及图3，在条纹相机200记录到压缩衍射全息图案后，数据处理系统300对采集的数据进行配准和变换，再经过裁剪和归一化处理后，使用英伟达RTX3090高性能显卡与PnP-DIP算法对记录到的压缩衍射全息图案并行计算，进行数据处理，进而重建出被测物体2的相位和振幅信息。图3即为数据处理系统300中的计算机10采集到的E形激光脉冲的相位和振幅数据。
参阅图3，CS-CMUI装置恢复了帧间隔为2ns的-5至5ns的6个强度和6相位图像。在这里，时间零点被设置为与最大激光脉冲强度相对应的时刻。
参阅图3，为了获得相位梯度，对所有相位图像进行一阶差分，可以很好地反映激光腔镜内部曲率的变化。空间相位分布呈现出一致的凹陷趋势，这与调Q激光器中后腔镜的物理特性相吻合。
实施例2
本实施例旨在观测纳秒激光烧蚀ITO薄膜(铟锡氧化物半导体透明导电膜)超快动态过程中的强度与相位变化。无论是基础研究还是应用方面，科学界广泛关注纳秒激光与材料相互作用的领域。超快激光能够改变材料的状态和特性，适用于几乎所有材料的高质量、高精度和复杂三维结构加工。本发明以观测纳秒激光烧蚀ITO薄膜的超快动态过程为例，旨在深入研究这一领域。
参阅图1及图4，探测光相位调制光路100中的相干照明光源1被图4中的烧蚀后透过的探测激光脉冲取代。
参阅图4，将纳秒泵浦激光脉冲斜射到ITO薄膜表面，并使用垂直于表面的另一个纳秒探测激光脉冲检测激光诱导烧蚀动力学。通过测量ITO薄膜透射探测激光脉冲的时空强度和相位演化，可以观察到泵浦激光脉冲引起的烧蚀过程。同时，还记录了未进行泵浦激光脉冲照射的透射探头激光脉冲的时空强度和相位分布，作为参考。
参阅图5，通过将被测激光脉冲的时空强度演变归一化为参考激光脉冲的时空强度演变，得到被测激光脉冲的时空强度变化。通过从被测激光脉冲的时空相位分布中减去参考激光脉冲的时空相位分布，可以确定被测激光脉冲的时空相位变化。
参阅图4，本实施例选用的样品为ITO薄膜，其涂覆在钠钙玻璃基板上；泵浦光用于在ITO样品上诱导烧蚀动态场景，其从玻璃基板侧入射，以增强烧蚀效果。
参阅图4，ITO薄膜的厚度为200±50nm；泵浦光与探测光均来自于532nm相干照明光源1的同两部分，测得入射泵浦激光脉冲的单脉冲能量为1.6mJ。
参阅图1及图4，烧蚀后透过的探测激光脉冲取代相干照明光源1对待测物体2进行照明，经过待测物体2衍射后传播70mm，被随机相位板3调制，形成一个大发散角、颗粒状的衍射光斑，继续传播30mm距离之后，携带物体复振幅信息的光束进入一台狭缝4完全打开的条纹相机200，条纹相机内部的光阴极5将光信号转换为电子信号通过微通道板6，电子脉冲在传输过程中受到扫描电极7偏转板内线性电场的作用产生偏转，使得不同时间到达的电子在空间上分开，最后被偏转的电子轰击荧光屏8后重新产生光信号并被内置相机图像采集9芯片记录下来，得到一张压缩衍射全息图案，完成正向探测过程。
参阅图1、图4及图5，此测量物体2是纳秒量级超快变化的动态事件。在条纹相机200记录到压缩衍射全息图案后，数据处理系统300对采集的数据进行配准和变换，再经过裁剪和归一化处理后，使用英伟达RTX3090高性能显卡与PnP-DIP算法对记录到的压缩衍射全息图案并行计算，进行数据处理，进而重建出超快动态事件的相位和振幅信息。图5即为数据处理系统300中的计算机10采集到的此纳秒激光烧蚀ITO薄膜过程的相位和振幅数据。
参阅图5，所选的-9—5ns强度和相位图像中，时间零点对应于最大泵浦激光脉冲强度。
参阅图5，在泵浦激光脉冲引起的烧蚀发生之前，相位分布几乎保持不变。随着泵浦激光脉冲内瞬时强度的增加，烧蚀区域的相位值开始减小。在大约-2ns时观察到突然的相位变化，然后该区域随着相位值的下降而扩大。
参阅图5，烧蚀坑的相形受激光通量的影响。当激光通量低时，相形为圆形或近圆形。然而，当激光通量高时，相形变得非常复杂，激光烧蚀导致的最大相位差可达0.516弧度。本发明对烧蚀动力学的详细观察为分析激光-物质相互作用机理和优化激光参数提供了关键解决方案，这反过来又促进了激光加工在制造、医药、航空航天等各个领域的进一步发展和应用。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (2)
1.一种单次曝光压缩超快相干调制成像装置，其特征在于，所述装置包括探测光相位调制光路(100)、条纹相机(200)与数据处理系统(300)；
所述的探测光相位调制光路(100)包括依次光路连接的相干照明光源(1)、含相位信息的待测物体(2)和随机相位调制板(3)；
所述的探测光相位调制光路(100)采用无透镜设计，降低了透镜造成的光路结构相位失真的影响。
所述的条纹相机(200)内部包括依次光路连接的入射狭缝(4)、光阴极(5)、微通道板(6)、扫描电极(7)、荧光屏(8)与内置相机(9)；
所述的数据处理系统(300)由计算机(10)构成，所述计算机(10)搭载计算显卡和CUDA通用并行计算架构；其中：
所述相干照明光源(1)的镜头中心、待测物体(2)的几何中心、随机相位调制板(3)的几何中心与条纹相机(200)的镜头中心在同一条直线上，此直线即为光的传播方向；
所述相干照明光源(1)、含相位信息的待测物体(2)、随机相位调制板(3)和条纹相机(200)分别对应着光路中的四个平面，依次分别为出射平面、入射平面、调制平面和探测平面；四个平面均垂直于光的传播方向，且光路与平面的交点在平面的几何中心。
2.根据权利要求1所述的单次曝光压缩超快相干调制成像装置，其特征在于，所述CUDA通用并行计算架构中含一种图像重建算法，所述算法基于即插即用框架和深度图像先验的一种空间二维至时空三维图像重构算法即PnP-DIP算法，其具体数据处理过程如下：
设待测物体的复振幅为：
U0(x,y,t；0)＝I(x,y,t；0)exp[iφ(x,y,t；0)] (1)
其中i代表虚数单位，x、y是物体的空间二维坐标，t代表时间，分号后的0代表当前的传播距离；I(x,y,t；0)表示物体的振幅信息，φ(x,y,t；0)表示物体的相位信息；经过d1距离的传播后，衍射光的复振幅表示为：
其中代表d1距离的传递函数，k代表波数，λ代表相干照明光源的波长；/>是U0的二维傅里叶变换，fx和fy分别是图像在x和y方向上的空间频率；经过随机相位调制板后的复振幅表示为：
其中m(x,y)代表相位调制板的空间相位分布；UM继续经过d2距离的传播后，衍射光的复振幅表示为：
此时，相机接收到的光强信息为：
在条纹相机内部，时变光强信息经过偏转后被积分叠加，最终得到的压缩结果为：
其中v代表条纹相机的偏转速度；上述所有过程用矩阵的形式表示为：
其中，U(x,y,t)代表待测物体的复振幅，Hd1和Hd2为空间传递算符，M代表相位调制算符，S代表条纹相机的时空偏转算符，T表示条纹相机CMOS在单次曝光过程中的积分算符；
正向采集过程进一步的简略表示为：
y＝Ax (8)
其中表示待测物体的复振幅分布，/>表示测量矩阵，表示最终获得的压缩数据结果；这里，变量W、H和B分别表示动态场景的x、y和t方向上的像素数，变量N表示压缩快照坐标系中y方向上的像素数；N的值大于H的值，但比H×B的值小；数据恢复的过程即是从已知信号y计算得到待测数据x的过程；从数学角度出发，已知y和A求解原始信号x是一个明显的欠定问题，因为已知数据的数量远小于未知数据的数量，所以要准确恢复信号必须增加其他的正则约束条件，缩小解空间的大小；同时，由于A的矩阵形式过于庞大，求解A-1的过程相当耗费计算资源和时间，在此，通过引入DIP来代替求逆过程，将逆问题改写为如下形式的最优化问题：
s.t.x＝G(θ,z)
其中R(x)是一个先验项，λ为正则化参数，用来平衡保真项和先验项；G(θ,z)代表一个可更新参数的深度卷积神经网络，其输入是一个与原始数据x同等大小的标准高斯分布白噪声，并且在优化过程中保持不变，θ为网络中可更新的权值参数，将G(θ,z)简写为G(θ)；进一步，通过引入一个可更新的辅助参数u和一个手动设置的惩罚参数ρ，并构造增广拉格朗日函数Lρ，将带约束的优化问题改写为无约束优化问题：
利用ADMM框架，的最小化问题分解为以下三个子问题:
1.
2.
3.uk+1＝uk+(xk+1-G(θk+1)). (13)
对于x的更新，观察到它是输入G(θk)-uk关于的近端算子的形式，在此使用即插即用设计利用降噪算法代替近端算子，以防止深度图像先验过拟合；对于θ的更新，在每次迭代中采用一个1阶梯度下降步骤；
从式(10)、式(11)、式(12)与式(13)中不难看出，随着迭代次数k的增加，求解对象xk+1会不断趋近于真实待测物体的复振幅分布x，最后在k无限大的时候收敛于x，根据式(8)，x即为待测物体的复振幅信息U0(x,y,t；0)＝I(x,y,t；0)exp[iφ(x,y,t；0)]；通过手动选择迭代次数k，从获得的压缩数据结果重建出被测物体的相位和振幅信息。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311844219.5A CN117761896A (zh) | 2023-12-29 | 2023-12-29 | 一种单次曝光压缩超快相干调制成像装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311844219.5A CN117761896A (zh) | 2023-12-29 | 2023-12-29 | 一种单次曝光压缩超快相干调制成像装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117761896A true CN117761896A (zh) | 2024-03-26 |
Family
ID=90316220
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311844219.5A Pending CN117761896A (zh) | 2023-12-29 | 2023-12-29 | 一种单次曝光压缩超快相干调制成像装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117761896A (zh) |
-
2023
- 2023-12-29 CN CN202311844219.5A patent/CN117761896A/zh active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Luo et al. | Computational imaging without a computer: seeing through random diffusers at the speed of light | |
| US10194100B2 (en) | Glare suppression through fog by optical phase conjugation assisted active cancellation | |
| CN108332866B (zh) | 基于振幅调制的波前在线快速重建装置和重建方法 | |
| CN106990694B (zh) | 一种部分相干光照明下的非迭代相位恢复装置及方法 | |
| CN106094487B (zh) | 基于多个记录距离的太赫兹同轴全息成像方法 | |
| CN111366557A (zh) | 一种基于薄散射介质的相位成像方法 | |
| CN110455834B (zh) | 基于光强传输方程的x射线单次曝光成像装置及方法 | |
| CN113340570B (zh) | 高功率脉冲激光动态光场及动态焦斑的测量装置和测量方法 | |
| Berdeu et al. | Reconstruction of in-line holograms: combining model-based and regularized inversion | |
| CN115372339A (zh) | 一种等离子体三维特性的重构方法 | |
| CN110048293A (zh) | 对畸变飞秒激光脉冲同步聚焦整形的装置及控制方法 | |
| Su et al. | Physics-Guided Terahertz Computational Imaging: A tutorial on state-of-the-art techniques | |
| Zhai et al. | End-to-end infrared radiation sensing technique based on holography-guided visual attention network | |
| Wen et al. | Phase retrieval with dynamic linear combination in multiple intensity measurements | |
| CN117761896A (zh) | 一种单次曝光压缩超快相干调制成像装置 | |
| CN108918529A (zh) | 一种透明介质的表面划痕成像检测装置及方法 | |
| CN110160663A (zh) | 一种高分辨率的近场波前测量装置和测量方法 | |
| CN114965364A (zh) | 单脉冲瞬态非线性折射率测量装置和测量方法 | |
| Chen et al. | 3PointTM: faster measurement of high-dimensional transmission matrices | |
| CN208255081U (zh) | 一种透明介质的表面划痕成像检测装置 | |
| Jin et al. | Single-Shot Intensity-and Phase-Sensitive Compressive Sensing-Based Coherent Modulation Ultrafast Imaging | |
| Hu et al. | Reflective coherent diffraction imaging with binary random sampling and updated support constraints | |
| US20240310283A1 (en) | Super-resolution lens-free microscopy | |
| CN117571140A (zh) | 基于单次曝光的激光脉冲时空复振幅测量装置及方法 | |
| Su et al. | Physics-guided terahertz computational imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |