一种单片式宽波段衍射计算成像方法

1.本发明涉及光学计算摄影领域，尤其涉及一种单片式宽波段衍射计算成像方法。


背景技术：

2.降低光学系统复杂性的高质量成像技术一直以来都是工业生产和科研工作的研究目标。为了获得高分辨率图像，传统折/反射式光学系统的做法是增大系统口径或引入多个光学元件来补偿像差，增加了成像系统结构的复杂性、制造成本以及成像系统的质量。衍射光学元件具有超薄和轻巧的物理结构，大而灵活的设计空间和更好的聚焦性能以及离轴成像性能，这些有点使其成为解决高质量成像与轻量化简单光学系统之间矛盾的途径之一。
3.然而，由于衍射本身存在的严重波长依赖性，在应用衍射光学元件成像时会出现强烈的色散，非设计波段的成像会出现大直径模糊，从而大大降低了宽波段成像时的图像质量，甚至在单个色彩通道的波长范围内也会出现色散模糊，因此传统衍射光学元件的应用常局限于窄波段成像系统。
4.随着现代计算机计算能力的飞速发展，计算成像技术已可将图像重建算法作为光学系统的虚拟组件进行引入，将光学系统的负担转移到计算机计算能力上。计算成像技术可以先采集后处理，经光学系统采集的信息需要经过图像端处理才能直接使用，其将光学设计与算法设计相结合，以最终经过处理后的图像质量作为目标来优化光学系统，所以光学部分的设计会根据后端处理方式不同而选择不同的设计目标(不一定以图像清晰度作为目标)，这也是计算成像区别于传统光学设计方法的地方。目前这一技术已使得一些超轻型单片镜成像系统问世，并且为进一步突破衍射光学元件在成像系统上的应用时所面临的带宽、视场、系统质量、复杂性等限制提供了方案。


技术实现要素：

5.本发明要解决的技术问题是：应用传统衍射元件成像时成像系统结构复杂、成本高、波段窄的问题，实现低成本、轻量化、宽波段、保持高色彩保真度的清晰成像。
6.为实现上述目的，本发明采取的技术方案如下：一种单片式宽波段衍射计算成像方法，包括如下步骤：
7.步骤1：设置优化消色差衍射透镜的设计参数，所述的设计参数包括光学孔径、设计波段、优化迭代次数，根据设计波段进行等间隔密集采样得到n个优化波长，记为[λ1,λ2,
…
,λn]，然后生成n个分别以λ1,λ2,
…
,λn作为设计波长的传统位相菲涅尔衍射透镜的台阶结构；
[0008]
步骤2：在光学孔径面按径向划分n个等面积区域，记为s1,s2,
…
,sn，区域si的结构与设计波长为λi的传统位相菲涅尔衍射透镜台阶结构相同，将s1,s2,
…
,sn组合形成分区消色差衍射透镜；
[0009]
步骤3：将所述分区消色差衍射透镜的所有台阶高度作为参数输入，使用粒子群优
化算法进行优化迭代，优化的成本函数为各个波长的点扩散函数与所有波长的平均点扩散函数的差值的加权和，当成本函数小于设定值或达到最大迭代次数时，所述的优化迭代即停止，此时得到的台阶结构即为优化消色差衍射透镜的台阶结构；
[0010]
步骤4：对一幅由所述优化消色差衍射透镜退化的图像，通过双三次样条插值，对其进行0.5倍下采样得到低尺度图像，通过tikhonov正则化对其进行第一次去卷积复原；
[0011]
步骤5：将步骤4中所述的去卷积复原后的图像再经过2倍上采样恢复到原来尺寸，使用l1正则化方法，并将所述上采样后的图像作为先验项加入到l1正则化中，进行第二次去卷积复原，将所述第二次去卷积复原的图像作为最终得到的图像；
[0012]
可选地，步骤1中的传统位相菲涅尔衍射元件是一种台阶型衍射元件，包括但不限于台阶型透射式衍射元件、台阶型反射式衍射元件以及不同台阶数量的衍射元件。
[0013]
可选地，步骤1中所述的设计波段，可以是可见光波段、红外波段等任意光学波段。
[0014]
可选地，步骤3所述的粒子群算法，可以是标准粒子群算法、自适应粒子群算法等经过改进的粒子群算法。
[0015]
可选地，步骤4所述的下采样倍数，其数值不一定为0.5倍，可以是任意小于1的正数。
[0016]
可选地，步骤4所述的去卷积方法，包括但不限于tikhonov正则化方法、totalvariation最小化方法、richard-lucy迭代复原算法等。
[0017]
可选地，步骤5中所述第二次去卷积复原的过程包括：根据半二次惩罚的方法引入一个辅助变量，替换所述l1正则化方法中的正则化项，并生成一个l2范数项，将所述第二次去卷积复原优化函数分解为l1、l2范数部分，根据l2范数部分固定所述辅助变量，直接求解估计清晰图像，根据l1范数部分固定估计图像求解辅助变量，经过交叉迭代优化求解最终估计清晰图像，即所述第二次去卷积复原图像。
[0018]
本发明与现有技术相比的优点在于：
[0019]
本发明的单片式宽波段衍元件计算成像可大幅度减小成像设备上光学系统的质量，简化光学系统结构，具有更小的离轴畸变，并同时拓宽衍射元件的带宽，保持高色彩保真度的清晰成像。在保持相同参数和图像清晰度的情况下，本发明所述的方法能大幅度缩减设备成本，具有容易制备、面型公差宽松、易于装调的优点，适用于集成化成像设备或质量、空间受限的成像系统。相较于传统衍射元件，本发明实施例中的衍射计算成像技术可在可见光波段下将图像的清晰度提高一个数量级以上。
附图说明
[0020]
图1为分区消色差衍射元件的示意图，它是优化消色差衍射元件的初始结构；其中，s1,s2,
…
,sn代表按径向划分的不同等面积区域，λ1,λ2,
…
,λn代表对设计波段密集采样得到的不同波长；
[0021]
图2为本发明实施例中使用的粒子群算法流程图；
[0022]
图3为本实施例中传统菲涅尔衍射透镜、分区消色差衍射透镜、优化消色差衍射透镜分别在450nm、505nm、555nm、605nm、655nm波长下得到的点扩散函数图像；
[0023]
图4为本发明中图像重建算法的框架图；
[0024]
图5为根据本发明设计的实施例的成像效果图，从左至右依次为：原始图像、优化
消色差衍射透镜退化图像、衍射计算成像复原图像；
[0025]
图6为本发明实施例的总体设计流程图。
具体实施方式
[0026]
为了更清楚地展示本发明的技术细节，下面结合附图和具体实施方式来进一步说明本发明。
[0027]
实施例为光学孔径为8毫米的消色差衍射元件计算成像设计，具体方法如下：
[0028]
步骤1：设置优化消色差衍射透镜的设计参数，所述的设计参数包括光学孔径、设计波段、优化迭代次数n，根据设计波段进行等间隔密集采样得到n个优化波长：[λ1,λ2,
…
,λn]，然后生成n个分别以λ1,λ2,
…
,λn作为设计波长的传统位相菲涅尔衍射透镜结构，在光学孔径面按径向划分n个等面积区域：s1,s2,
…
,sn，区域si的结构与设计波长为λi的传统位相菲涅尔衍射透镜结构相同，将s1,s2,
…
,sn组合形成分区消色差衍射透镜；
[0029]
步骤2：构建衍射元件台阶高度δ与点扩散函数psf的映射模型模型：
[0030][0031]
其中，f(
·
)代表傅里叶变换，p为孔径函数，δ为微结构台阶高度，n
λ
为在波长λ下衍射元件基底的折射率，(u,v)为衍射元件面的坐标，再构造成本函数为：
[0032][0033]
式中，pi(δ)为各波长优化过程中的强度psf，t为各波长的平均psf，λ为设计波段，||
·
||1代表l1范数，wi作为成本函数的自适应权重，其取值为：
[0034][0035]
步骤3：获取分区消色差衍射透镜所有微结构台阶高度，将其作为初始参数δ
ini
，使用粒子群优化算法对式(2)进行优化迭代，直到满足或达到迭代次数n终止，此时得到的结构δn即为优化消色差衍射透镜的结构，记为δ
opt
，图3示出了本发明的实施例中粒子群算法的参数设置及优化流程，在迭代的第k 1代时，衍射元件高度参数及高度参数的速度矢量将按照δ
k 1
＝δk vk、v
k 1
＝vk c1·
rand()
·
(p
k-δk) c2·
rand()
·
(gbest-δk)进行更新，其中，c1、c2为学习因子，是大于0的常数，rand()代表随机分布在[0,1]上的随机数，pk来记录每个粒子的历史最优位置，gbest记录所有粒子的全局最优位置；
[0036]
图3示出了本实施例中优化消色差衍射透镜的点扩散函数设计结果，由第一排至第三排分别是：传统菲涅尔衍射透镜、分区消色差衍射透镜、优化消色差衍射透镜在各各波长下的点扩散函数表现；
[0037]
步骤4：图4示出了本发明图像复原算法的框架。对一幅由优化消色差衍射透镜退化的图像，通过双三次样条插值，对其进行0.5倍下采样得到低尺度图像id，通过tikhonov正则化对其进行第一次去卷积复原：
[0038][0039]
其中，f(
·
)
*
为f(
·
)的复共轭，f(
·
)代表傅里叶变换，f-1
(
·
)代表傅里叶逆变换，μ和β是每一项的权重，k是优化消色差衍射透镜的psf，由δ
opt
代入公式(1)得到，d是一阶导数滤波器，bd是经优化消色差衍射透镜退化后的下采样图像；
[0040]
步骤5：将步骤4得到的复原图像再经过去噪后上采样恢复到原来尺寸，记为is，将is作为第二次去卷积的先验项进行复原，复原问题可描述为：
[0041][0042]
其中，i是待估计的原始清晰图像，根据半二次惩罚的方法，引入一个辅助变量q来替换原来的变量di，并将优化函数还写为：
[0043][0044]
其中，γ和ρ是权重因子，对式(5)分成l1、l2范数两个部分，对l2部分，可将q固定，对i直接求取：
[0045][0046]
对l1部分，固定i，根据梯度下降法求q的近端算子对q进行求解：
[0047][0048]
其中，sign(
·
)代表符号函数。
[0049]
对i、q进行交替迭代求解，直到达到最大迭代次数停止，所得到的i即为衍射计算成像的最终图像。
[0050]
图5示出了本发明的成像效果图，图左为原始图像，图中为优化消色差衍射元件直接成像的图像，图右为衍射计算成像实施例最终得到的图像。图6示出了本发明的实施例设计流程，主要分为光学设计和图像重建两个方面。
[0051]
显然，上述实施例仅仅是为清楚地说明所作的举例，而非对实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引申出的显而易见的变化或变动仍处于本发明创造的保护范围之中。