一种基于电光晶体偏转的压缩超高速摄影装置的制作方法

本发明属于超快成像技术领域，可用于超快物理、化学、生物等过程的三维x-y-t重建，其中x-y为空间坐标，t为时间坐标，也可用于微纳加工以及一些三维物体的动态观测。

背景技术：

高速成像技术在解决无数科学谜团、推进医学发展、创造艺术效果、揭示细胞甚至机器等事物的内部运作等方面发挥着不可或缺的作用。人们对高速摄影的追求最早的记录可以追溯到1878年记录的奔跑中的马以及1887年对于超声速子弹的摄影。然而，直到二十世纪末，成像速度还始终在每秒10^5帧。后来，基于电荷耦合器件(ccd)和互补金属氧化物半导体(cmos)的电子成像传感器的出现彻底改变了人们对高速成像的认识，使得成像速度可以达到每秒10^7帧。尽管这种传感器得到了非常广泛的应用，但是成像速度受限于ccd和cmos的芯片存储和电子的读出速度，在很多超快领域无法进行直接探测。

电光晶体的普克尔效应可以根据信号电压实现光路的小角度偏转，并且固有响应时间在飞秒量级。主要应用有光学囚禁、q开关激光器、光束偏转、晶体条纹相机。普克尔效应决定了电光晶体在不同的电压下有不同的折射率，类似于光学条纹相机，通过将不同时刻的光信息进行偏转实现将其时间信息转化为空间信息。在成像领域，两幅图像之间的时间间隔决定了是否能很好地区分两帧图像，也即决定了时间分辨能力。在偏转压缩成像领域，时间信息被转换为空间信息，通过对两幅图像之间进行合适的偏转可以实现两帧图像的区分。通常所采集的图像具有一定的空间尺寸，在不偏转的情况下，即使空间二维尺寸非常狭长也会占据很多像素。为了区分相邻帧图像，在不编码的情况下，必须依照每一张图像的本身尺寸作为区分两帧图像的界限；而在编码的情况下，就可以依据面积远小于图像本身尺寸的单像素界限来区分两帧图像，这样可以提高整个系统的时间分辨率。

技术实现要素：

本发明的目的在于提供一种基于电光晶体偏转的压缩超高速摄影装置，该装置基于压缩感知原理和电光晶体的普克尔效应能够克服现有技术中的电光晶体一维成像的限制，对一些超快物理、化学、生物等的动态过程进行成像。

实现本发明目的的具体技术方案是：

一种基于电光晶体偏转的压缩超高速摄影装置，特点是：该装置包括格兰棱镜、第一凸透镜、数字微镜器件、电光晶体、第二凸透镜、第三凸透镜、ccd相机、计算机及数字延时发生器，所述格兰棱镜、第一凸透镜及数字微镜器件依次光路连接，数字微镜器件、电光晶体、第二凸透镜、第三凸透镜、ccd相机依次光路连接，ccd相机与计算机电路连接，数字延时发生器与ccd相机电路连接；其中，所述格兰棱镜用于产生偏振光。

所述电光晶体所匹配的高压电源为100khzq-switchdriver，电压范围为0-2000v,电压上升沿时间为12ns。

所述计算机采用增广拉格朗日(a-l)算法对ccd相机中压缩后的采集结果进行重构，得到高速动态过程的三维数据。

所述增广拉格朗日(a-l)算法，具体包括：

设：被摄对象即动态场景记为x，在ccd相机中得到结果记为y，数据采集过程为y＝lx，l＝ise，其中e为空间编码算符，s为时间剪切算符，i为时空积分算符；解决如下最优解问题：

其中λ为算法乘子，β为正则化参数，φ(x)为全变分函数；

第一步：引入一个新的变量w，其中w＝dx，d为梯度算符，令φ(x)＝||dx||2，对其进行约束再变形后，上述方程就变为：

其中ν是φ(x)的拉格朗日乘子，μ为对应的正则化参数；

第二步：在每次迭代的过程中，将(1)所描述的问题分解为关于变量w和x的两个子问题

w-子问题：

对应的解为：

x-子问题：

对应的解为：

xj＝xj-1-αd(xj-1)(5)

其中d(x)＝(μd^t(dx-w)-d^tν)+βl^t(lx-y)-l^tλ为x的导数，α为迭代优化参数，t为转置矩阵符号；

第三步：将(3)和(5)代入(1)重复第一步、第二步，寻求最优解x。

所述ccd相机(8)的起始工作时间使用数字延时发生器dg645(10)进行同步，其工作时间涵盖动态场景的有效持续时间。

本发明的优点是：

1)首次将电光晶体用于三维(x-y-t)成像，时间分辨率为500ps，空间分辨率在微米量级；

2)单次拍照再现过程，可以拍摄不重复或者不可逆事件的发生；

3)可接收式成像，在自发光场景的探测中不需要进行主动探测光照明；

附图说明

图1为本发明结构示意图；

图2为本发明使用状态示意图；

图3为本发明的实验原理图；

图4为本发明拍摄的罗丹明b溶液的荧光衰减过程示意图。

具体实施方式

以下结合附图和实施例进一步详细阐述本发明。

参阅图1，本发明包括格兰棱镜2、第一凸透镜3、数字微镜器件4、电光晶体5、第二凸透镜6、第三凸透镜7、ccd相机8、计算机9、数字延时发生器10，所述格兰棱镜2、第一凸透镜3及数字微镜器件4依次光路连接，数字微镜器件4、电光晶体5、第二凸透镜6、第三凸透镜7、ccd相机8依次光路连接，ccd相机8与计算机9电路连接，数字延时发生器10与ccd相机8电路连接。

参阅图2，本发明使用时，动态场景1为包含时间信息的x-y-t三维矩阵的数据立方，光信息进入格兰棱镜2产生偏振的光信息，然后进入第一凸透镜3将不同时刻的图像信息进行缩小，使得像的尺寸小于数字微镜器件4的有效面积，然后数字微镜器件4对到达其上的光信息进行0-1伪随机二元化编码，将反射回的光信息进一步传送至电光晶体5中。电光晶体5在高压电源的控制下，不同时刻具有不同的电压值，由于电光晶体5的普克尔效应，不同时刻就会具有不同的折射率，折射角不同，不同时刻进入电光晶体5的光信息的传播方向就会进行发散，不同时刻的折射光在电光晶体5之后进入第二凸透镜6、第三凸透镜7中进行像的放大，然后进入ccd相机8中。ccd相机8的起始工作时间使用数字延时发生器dg64510进行同步，其工作时间涵盖动态场景的有效持续时间。ccd相机8的曝光时间内就会产生一个时间偏转的光信息的积分效应。

本发明包括正演过程和反演过程。参阅图3，正演过程是高速动态场景x由数字微镜器件4进行编码，然后反射进入电光晶体5进行偏转，偏转之后的图像叠加显示在ccd相机8中得到拍摄结果y。整个过程可概括为y＝lx，l＝ise，其中e为空间编码算符，s为时间剪切算符，i为时空积分算符。反演过程基于压缩感知算法，采用增广拉格朗日方法解决如下欠定最优解问题：

其中λ为算法乘子，β为正则化参数，φ(x)为全变分函数。计算机对ccd相机中采集到的数据进行处理后可得到高速动态过程的三维数据。

在具体实现中，基于可分离变量的原理，利用tv项(φ(x))引入一个新的变量w，其中w＝dx，d为梯度算符，令φ(x)＝||dx||2，对其进行约束再变形后，上述方程就变为：

其中ν是tv项的拉格朗日乘子，μ为对应的正则化参数。

在每次迭代的过程中，式(1)所描述的问题可分解为关于变量w和x的两个子问题。

w-子问题：

对应的解为：

x-子问题：

对应的解为：

xj＝xj-1-αd(xj-1)(5)

其中d(x)＝(μd^t(dx-w)-d^tν)+βl^t(lx-y)-l^tλ为x的导数，α为迭代优化参数，t为转置矩阵符号；

每次迭代的过程就是将(3)和(5)代入(1)反复寻求最优解x的过程。

实施例

本实施例各部件参阅图2设置，包括动态场景1、格兰棱镜2、第一凸透镜3、数字微镜器件4、电光晶体5、第二凸透镜6、第三凸透镜7、ccd相机8及计算机9。动态场景1发出的散射光依次进入格兰棱镜2、第一凸透镜3、数字微镜器件4，然后经过数字微镜器件4编码后进入电光晶体5、第二凸透镜6、第三凸透镜7和ccd相机8。ccd相机8接收到数字延时发生器10的信号后开始工作。计算机9对ccd相机中采集到的数据利用增广拉格朗日算法进行重构。

在本实施例中，如图4(a)所示，使用50fs的激光脉冲打在罗丹明b溶液中，产生其荧光衰减的动态过程，然后用本发明重现整个动态过程。图4(b)为使用增广拉格朗日算法重构出来的荧光光斑演化图像，图4(c)为(b)中对应的不同时刻的光强总和拟合出的荧光强度变化曲线，发现其可以和e指数衰减函数相吻合(图中小方框为实验数据，线条为拟合结果)。本发明对荧光捕捉的是二维信号，是一个二维信号演化的过程，并且拟合出其寿命为2.84ns，和已报道的荧光寿命相近。

于此，本实施例完成了三维动态场景的拍摄。