一种基于微透镜阵列与光子集成芯片干涉成像方法及系统与流程

1.本发明涉及光学领域，具体而言，涉及一种基于微透镜阵列与光子集成芯片干涉成像方法及系统。


背景技术：

2.现有技术中，普通基于精确折光原理的成像系统由成像物镜组和光电转换器件阵列构成，物体光进入成像物镜组，形成的像由光电转换器件阵列接收，从而得到物体图像，通过成像系统可实时记录和观测物体形貌。传统光学望远镜的基本设计原理仍基于精确折光的设计理念，受衍射极限约束，为提高分辨率，必须增大系统口径，为实现更高的分辨率，传统空间望远镜的体积、重量将极其庞大。例如哈勃空间望远镜主镜口径为2.4m， jwst望远镜采用拼接主镜为6m，正在论证的atlast望远镜将达到惊人的 8m口径。其发射成本和在轨维护难度具有极大的挑战性。
3.光波干涉是一种基本的光学物理现象，光学干涉测量是通过产生相干的两路光，一路作为参考光，一路作为测量光，参考光照射到位置固定的参考反射镜，测量光照射到被测对象，由参考反射镜反射回的参考光与被测对象反射回的物光光波叠加发生干涉，干涉条纹反映两路光光程差信息，从而根据干涉条纹的变化获得被测对象的特征信息。在干涉成像中，两路干涉光都来自观测目标，远场成像目标的光到达孔径平面，孔径平面光分布即为成像目标的频谱，两个干涉孔径接收到频域信息并耦合到光子集成芯片，通过设计好的光子集成芯片结构来解调出频域信息的幅值和相位部分。
4.在richard.l.kendrick等人提出了一种用于光电侦察的分段平面成像探测器(segmented planar imaging detector for electro
‑
opticalreconnaissance，spider)的概念中，即是干涉成像的一种形式，richard. l.kendrick等人采用在光子集成芯片前段放置小透镜，其光路的共相调整难度大，每个小透镜都有6个调整的自由度，实现干涉成像往往需要大量用来进行两两干涉的小透镜，因此其调整难度很大，急需一种容易进行调节的且可以提供聚光功能的器件代替分立的小透镜。


技术实现要素：

5.本发明实施例提供了一种基于微透镜阵列与光子集成芯片干涉成像方法及系统，以至少解决现有干涉成像系统调整难度大的技术问题。
6.根据本发明的一实施例，提供了一种基于微透镜阵列与光子集成芯片干涉成像方法，包括以下步骤：
7.利用微透镜阵列接收光信号，同时采用vhc的光耦合方式接收从微透镜阵列出射的光信号，对光耦合后的光信号进行处理获取观测目标频域信息的解算信号：
8.解算出解算信号中观测目标的频域信息，在频域信息中恢复观测目标。
9.进一步地，利用微透镜阵列接收光信号，同时采用vhc的光耦合方式接收从微透镜阵列出射的光信号，对光耦合后的光信号进行处理获取观测目标频域信息的解算信号包
括：
10.排布微透镜阵列；
11.对频域信息进行采样；
12.位于物体平面的光信号通过自由空间传播至微透镜阵列所在的孔径平面，微透镜阵列将光信号进行汇聚；
13.将光信号从孔径平面传输到光子集成芯片的光电探测器上。
14.进一步地，位于物体平面的光信号通过自由空间传播至微透镜阵列所在的孔径平面，微透镜阵列将光信号进行汇聚包括：
15.物体平面的光信号传播到孔径平面后得到孔径平面上的互谱密度函数观测目标的傅里叶变化；孔径平面上的互谱密度函数表示为：
[0016][0017]
其中空间频率为α'表示观测平面上某一点的坐标；表示孔径平面上任意一点的坐标，观测光的频率为ν、波长为λ，非相干源的归一化频谱可以表示为g(ν)，并有∫g(ν)dν＝1，观测目标光强分布为i
object
(α')；孔径平面上，任意两点之间的互谱密度函数只与孔径平面上的两个干涉位置差有关，与物体的绝对位置无关；互谱密度函数是物体信号的傅里叶变换，f1‑
f2越大观测目标越高频的信号。
[0018]
进一步地，将光信号从孔径平面传输到光子集成芯片的光电探测器上包括：
[0019]
光信号从孔径平面到光子集成芯片的光电探测器传播，微透镜阵列将光汇聚到光子集成芯片上的vhc结构上，vhc结构为垂直转水平耦合，是一种片上光栅结构，vhc结构将汇聚光耦合入光子集成芯片内部，经过光子集成芯片上的输入波导阵列、波分解复用结构、频谱信号解算结构，最后被光子集成芯片的光电探测器接收，并转换为电信号；此部分过程描述为：
[0020]
i
k
＝|m
k
|2＝s2 l2 2slcos[(φ
s
‑
φ
l
) φ
k
]
[0021]
其中i
k
表示光子集成芯片输出的光强信号，m
k
表示光子集成芯片输出的振幅；s、φ
s
表示其中一个输入波导的幅值和相位；l、φ
l
表示另一个输入波导的幅值和相位；φ
k
表示光子集成芯片引入的相位差，φ
k
的值通过光子集成芯片可设为
[0022]
进一步地，对频域信息进行采样包括：
[0023]
频域得到的观测目标的空间频率f与基线b的关系表示为：
[0024]
f＝λz/(d
i
‑
d
j
)
[0025]
＝λz/b
[0026]
λ表示成像波长，z表示物体平面到孔径平面的距离；b＝b
i,j
＝d
i
‑
d
j
表示径向上两个不同的微透镜中心位置之差，称为基线baseline；空间频率f表示频域坐标，α'表示空域坐标；观测目标i
object
(α')通过fouriertransform即可得到频域分布信息，频谱信息表示为：
[0027]
v(f)＝ft{i
object
(α')}
[0028]
其频域的采样方式表示为：
[0029][0030]
其中h(f)表示在空间频率f处的对频谱信息v(f)的采样权重。
[0031]
进一步地，解算出解算信号中观测目标的频域信息包括：
[0032]
经过微透镜汇聚的光耦合到输入波导阵列，在经过波分解复用结构可以得到多个准单色光λ
k
信号，这些信号通过片上频谱信号解算结构对应基线长度b
ij
可以得到四个光强信号输出i
k
(k＝1,2,3,4)，这些光强信号通过abcd方法即可解算对应空间频率为f的频域信号：
[0033][0034][0035]
v(f)表示空间频率f对应的频谱信号，abs{v(f)}表示其幅值部分，arg{v(f)}表示其相位部分。
[0036]
进一步地，在频域信息中恢复观测目标包括：
[0037]
对观测目标频域信息降质影响表示为：
[0038]
v
sample
＝v
·
h
[0039]
离散情况下为：
[0040]
v
sample
(f
kij
)＝v
object
(f
kij
)
·
h(f
kij
)
[0041]
其中v表示观测目标的频域信息分布，v(f)＝ft{i
object
(α')}是观测目标i
object
的fouriertransform，v
sample
表示采样得到的频域信息，
·
表示点乘操作；对于v、v
sample
离散的分量为v
ij
、降质过程是和采样的离散分量h
ij
的对应值相乘得到:
[0042][0043]
进一步地，频域信息恢复观测目标为对于已采样的频谱信息部分做最大程度的保留，添加datafidelity项；对于未采样部分选择使用惩罚函数来添加对观测目标的先验；
[0044]
最终需要优化的组合函数包括两个部分datafidelity项α
·
χ2及penaltyfunction项β
·
p，恢复观测目标采用如下公式：
[0045]
j＝α
·
χ2 β
·
p
[0046][0047][0048]
v
sample
＝v(f)＝abs{v(f)}exp(i
·
arg{v(f)})
[0049]
其中，α，β表示各项的系数，表示对矩阵m中的所有元素求平方和，即表示对矩阵m中的所有元素求平方和，即表示对观测目标的估计值，a表示系统对观测目标的降质作用，v
sample
表示采样得到的频域信息；当j的值收敛到最小值，即可得；p为惩罚函数，采用最
大熵maximum entropy惩罚，此外添加tp限制即：
[0050][0051]
其中x
ij
表示x中的分量；通过采用不同的优化算法计算迭代，使得组合函数j的值收敛到最小值，得到惩罚函数为：
[0052][0053]
其中α
ij
、a
ij
、x
ij
、v
ijsample
、β'
ij
、β”ij
表示α、a、x、v
sample
、β'、β”的分量；x表示观测目标的估计值选用梯度下降进行x的更新可以得到：
[0054][0055]
通过多次的迭代得到频域信息恢复观测目标。
[0056]
根据本发明的另一实施例，提供了一种基于微透镜阵列与光子集成芯片干涉成像系统，包括：
[0057]
解算信号获取单元，用于利用微透镜阵列接收光信号，同时采用vhc 的光耦合方式接收从微透镜阵列出射的光信号，对光耦合后的光信号进行处理获取观测目标频域信息的解算信号：
[0058]
观测目标恢复单元，用于解算出解算信号中观测目标的频域信息，在频域信息中恢复观测目标。
[0059]
进一步地，解算信号获取单元包括：
[0060]
排布单元，用于排布微透镜阵列；
[0061]
采样单元，用于对频域信息进行采样；
[0062]
光信号汇聚单元，用于位于物体平面的光信号通过自由空间传播至微透镜阵列所在的孔径平面，微透镜阵列将光信号进行汇聚；
[0063]
光信号接收单元，用于将光信号从孔径平面传输到光子集成芯片的光电探测器上。
[0064]
本发明实施例中的基于微透镜阵列与光子集成芯片干涉成像方法及系统，利用微透镜阵列接收光信号，同时采用vhc的光耦合方式接收从微透镜阵列出射的光信号，对光耦合后的光信号进行处理获取观测目标频域信息的解算信号，解算出解算信号中观测目标的频域信息，在频域信息中恢复观测目标。本发明利用微透镜阵列代替其分立的小透镜，共相调整难度可以得到很大程度的降低，从而进一步优化系统更有利于集成，提升系统的稳定性。采用基于微透镜阵列与干涉成像系统还可以进一步降低系统的体积。采用vhc的耦合方式可以降低空间光耦合入芯片的对光难度，还可以大大压缩系统的长度，依靠超薄的微透镜阵列和光子集成芯片即可完成成像，从而实现平板成像。
附图说明
[0065]
此处所说明的附图用来提供对本发明的进一步理解，构成本技术的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
[0066]
图1是本发明实施例提供的基于微透镜阵列与光子集成芯片干涉平板成像系统的
流程图；
[0067]
图2是本发明实施例提供的光子集成芯片的结构示意图；
[0068]
图3是本发明实施例提供的单个微透镜后端光路示意图；
[0069]
图4是本发明实施例提供的单个微透镜后焦面结构；
[0070]
图5是本发明实施例提供的光栅阵列中单个光栅的结构示意图；
[0071]
图6是本发明实施例提供的光路模型图；
[0072]
图7是本发明实施例提供的微透镜与光子集成芯片的结构示意图；
[0073]
图8是本发明实施例提供的频谱信号解算结构示意图；
[0074]
图9是本发明实施例提供的频谱信号解算结构和光电探测器的组合示意图；
[0075]
图10是本发明实施例提供的仿真的观测目标重建结果。
[0076]
附图标号说明：
[0077]
10、微透镜阵列；11、第一基板；12、微透镜列；
[0078]
20、光子集成芯片；21、光栅阵列；22、波分解复用结构；23、频谱信号解算结构；231、移相结构；232、波导耦合结构；24、光电探测器； 25、第二基板；26、背板。
具体实施方式
[0079]
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
[0080]
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
[0081]
实施例1
[0082]
根据本发明一实施例，提供了一种基于微透镜阵列与光子集成芯片干涉成像方法，参见图1
‑
10，包括以下步骤：
[0083]
利用微透镜阵列10接收光信号，同时采用vhc的光耦合方式接收从微透镜阵列10出射的光信号，对光耦合后的光信号进行处理获取观测目标频域信息的解算信号：
[0084]
解算出解算信号中观测目标的频域信息，在频域信息中恢复观测目标。
[0085]
本发明实施例中的基于微透镜阵列与光子集成芯片干涉成像方法，利用微透镜阵列10接收光信号，同时采用vhc的光耦合方式接收从微透镜阵列10出射的光信号，对光耦合后的光信号进行处理获取观测目标频域信息的解算信号，解算出解算信号中观测目标的频域信息，在频域信息中恢复观测目标。本发明利用微透镜阵列10代替其分立的小透镜，共相调整难度可以得到很大程度的降低，从而进一步优化系统更有利于集成，提升系统的稳定
性。采用基于微透镜阵列10与干涉成像系统还可以进一步降低系统的体积。采用vhc的耦合方式可以降低空间光耦合入芯片的对光难度，还可以大大压缩系统的长度，依靠超薄的微透镜阵列10和光子集成芯片20即可完成成像，从而实现平板成像。
[0086]
其中，利用微透镜阵列10接收光信号，同时采用vhc的光耦合方式接收从微透镜阵列10出射的光信号，对光耦合后的光信号进行处理获取观测目标频域信息的解算信号包括：
[0087]
排布微透镜阵列10；
[0088]
对频域信息进行采样；
[0089]
位于物体平面的光信号通过自由空间传播至微透镜阵列10所在的孔径平面，微透镜阵列10将光信号进行汇聚；；
[0090]
将光信号从孔径平面传输到光子集成芯片20的光电探测器24上。
[0091]
其中，位于物体平面的光信号通过自由空间传播至微透镜阵列10所在的孔径平面，微透镜阵列10将光信号进行汇聚包括：
[0092]
物体平面的光信号传播到孔径平面后得到孔径平面上的互谱密度函数观测目标的傅里叶变化；孔径平面上的互谱密度函数表示为：
[0093][0094]
其中空间频率为α'表示观测平面上某一点的坐标；表示孔径平面上任意一点的坐标，观测光的频率为ν、波长为λ，非相干源的归一化频谱可以表示为g(ν)，并有∫g(ν)dν＝1，观测目标光强分布为i
object
(α')；孔径平面上，任意两点之间的互谱密度函数只与孔径平面上的两个干涉位置差有关，与物体的绝对位置无关；互谱密度函数是物体信号的傅里叶变换，f1‑
f2越大观测目标越高频的信号。
[0095]
其中，将光信号从孔径平面传输到光子集成芯片20的光电探测器24 上包括：
[0096]
光信号从孔径平面到光子集成芯片20的光电探测器24传播，微透镜阵列10将光汇聚到光子集成芯片20上的光栅阵列21即vhc结构上(垂直转水平耦合，vertical
‑
to
‑
horizontal converter，是一种片上光栅结构： surface grating)，vhc结构可以将汇聚光耦合入光子集成芯片20内部，经过光子集成芯片20上的输入波导(surface grating)阵列、波分解复用结构22、频谱信号解算结构23，最后被光子集成芯片20的光电探测器 24接收，并转换为电信号；此部分过程描述为：
[0097]
i
k
＝|m
k
|2＝s2 l2 2slcos[(φ
s
‑
φ
l
) φ
k
]
[0098]
其中i
k
表示光子集成芯片20输出的光强信号，m
k
表示光子集成芯片 20输出的振幅；s、φ
s
表示其中一个输入波导的幅值和相位；l、φ
l
表示另一个输入波导的幅值和相位；φ
k
表示光子集成芯片20引入的相位差，φ
k
的值通过光子集成芯片20可设为
[0099]
其中，对频域信息进行采样包括：
[0100]
频域得到的观测目标的空间频率f与基线b的关系表示为：
[0101]
f＝λz/(d
i
‑
d
j
)
[0102]
＝λz/b
[0103]
λ表示成像波长，z表示物体平面到孔径平面的距离；b＝b
i,j
＝d
i
‑
d
j
表示径向上两个不同的微透镜中心位置之差，称为基线baseline；空间频率f表示频域坐标，α'表示空域坐标；观测目标i
object
(α')通过fouriertransform即可得到频域分布信息，频谱信息表示为：
[0104]
v(f)＝ft{i
object
(α')}
[0105]
其频域的采样方式表示为：
[0106][0107]
其中h(f)表示在空间频率f处的对频谱信息v(f)的采样权重。
[0108]
其中，解算出解算信号中观测目标的频域信息包括：
[0109]
经过微透镜汇聚的光耦合到输入波导(surfacegrating)阵列，在经过波分解复用结构22可以得到多个准单色光λ
k
信号，这些信号通过片上频谱信号解算结构23对应基线长度b
ij
可以得到四个光强信号输出i
k
(k＝1,2,3,4)，这些光强信号通过abcd方法即可解算对应空间频率为f
nij
的频域信号：
[0110][0111][0112]
v(f)表示空间频率f对应的频谱信号，abs{v(f)}表示其幅值部分，arg{v(f)}表示其相位部分。
[0113]
其中，在频域信息中恢复观测目标包括：
[0114]
对观测目标频域信息降质影响表示为：
[0115]
v
sample
＝v
·
h
[0116]
离散情况下为：
[0117]
v
sample
(f
kij
)＝v
object
(f
kij
)
·
h(f
kij
)
[0118]
其中v表示观测目标的频域信息分布，v(f)＝ft{i
object
(α')}是观测目标i
object
的fouriertransform，v
sample
表示采样得到的频域信息，
·
表示点乘操作；对于v、v
sample
离散的分量为v
ij
、降质过程是和采样的离散分量h
ij
的对应值相乘得到:
[0119][0120]
其中，频域信息恢复观测目标为对于已采样的频谱信息部分做最大程度的保留，添加datafidelity项；对于未采样部分选择使用惩罚函数来添加对观测目标的先验；
[0121]
最终需要优化的组合函数包括两个部分datafidelity项α
·
χ2及penaltyfunction项β
·
p，恢复观测目标采用如下公式：
[0122]
j＝α
·
χ2 β
·
p
[0123]
[0124][0125]
v
sample
＝v(f)＝abs{v(f)}exp(i
·
arg{v(f)})
[0126]
其中，α，β表示各项的系数，表示对矩阵m中的所有元素求平方和，即表示对矩阵m中的所有元素求平方和，即表示对观测目标的估计值，a表示系统对观测目标的降质作用，v
sample
表示采样得到的频域信息；当j的值收敛到最小值，即可得；p为惩罚函数，采用最大熵maximum entropy惩罚，此外添加tp限制即：
[0127][0128]
其中x
ij
表示x中的分量；通过采用不同的优化算法计算迭代，使得组合函数j的值收敛到最小值，得到惩罚函数为：
[0129][0130]
其中α
ij
、a
ij
、x
ij
、v
ijsample
、β'
ij
、β”ij
表示α、a、x、v
sample
、β'、β”的分量；x表示观测目标的估计值选用梯度下降进行x的更新可以得到：
[0131][0132]
通过多次的迭代得到频域信息恢复观测目标。
[0133]
下面以具体实施例，对本发明的基于微透镜阵列与光子集成芯片干涉成像方法进行详细说明：
[0134]
具体参见图1
‑
10，本发明使用阵列的微透镜形式代替每个分立的小透镜，其光路的调整难度大大降低，可以更有利于集成，提升系统的稳定性。微透镜是一种直径为微米级的小透镜，微透镜按一定的周期以一维或二维的方式排列在基底上构成微透镜阵列10。微透镜阵列10利用其表面波长量级的浮雕结构对光波进行调制，将一个完整的激光波前在空间上分成许多微小的部分，每一部分都被相应的小透镜聚焦在焦平面上，一系列微透镜就可以得到由很多均匀、规则排布的焦点组成的平面。将该结构引入干涉成像代替原来的分立透镜可以大大降低系统的调整难度。
[0135]
如图3所示为空间光耦合波导的两种形式：端面耦合以及垂直转水平耦合(vhc，vertical
‑
to
‑
horizontal converter，是一种片上光栅结构： surface grating)。可以将垂直或者有一定偏角照射的光转换到水平输出，spider采用端面耦合的形式其对空间光对准波导的角度精度以及位置精度要求很高，本发明提出采用vhc的耦合方式可以降低空间光耦合入芯片的对光难度，还可以大大压缩系统的长度，依靠超薄的微透镜阵列10和光子集成芯片20即可完成成像，从而实现平板成像。
[0136]
因此基于已有的精确折光传统大口径成像系统存在着体积、功耗和重量大的缺点，本发明提供了一种基于微透镜阵列与光子集成芯片干涉成像方法，使用多个小口径透镜代替大口径成像系统，利用波导作为导波和光干涉介质，可以大大减小系统的体积、功耗和重量；其次使用了光子集成芯片20的移相结构231组件(例如热控或者电控移相器)可以降低成像系统的共相调整难度。
[0137]
相比于spider系统，本发明利用微透镜阵列10代替其分立的小透镜，共相调整难
度可以得到很大程度的降低，从而进一步优化系统更有利于集成，提升系统的稳定性。采用基于微透镜阵列与光子集成芯片干涉成像方方法还可以进一步降低系统的体积。采用vhc的耦合方式可以降低空间光耦合入芯片的对光难度，还可以大大压缩系统的长度，依靠超薄的微透镜阵列10和光子集成芯片20即可完成成像，从而实现平板成像。
[0138]
本发明公开的一种基于微透镜阵列与光子集成芯片干涉成像方法中，远场目标发出的光经过自由传播到达微透镜阵列10，经过每个微透镜汇聚后耦合入每个微透镜后端对应的光子集成芯片20的输入波导(surfacegrating)阵列中(如图4
‑
5所示，每个输入波导都有一定的视场，将输出波导做成阵列结构可以增加其成像视场，光子集成芯片20利用波分解复用结构22(例如arrayedwaveguidegrating结构)、频谱信号解算结构23(如图9，包含移相结构231、波导耦合结构232)对耦合入芯片光分别分频、移相和干涉操作，光子集成芯光子集成芯片20的光电探测器24(photoelectricdetector，pd)将光信号转换为电信号输出，可得到光强值，从这些光强值可以计算出观测目标的频域信息的幅值和相位值。
[0139]
本发明公开的一种基于微透镜阵列与光子集成芯片干涉成像方法，其整体流程如图1所示，具体包括如下步骤：
[0140]
步骤1：获取观测目标频域信息的解算信号
[0141]
如图2所示该系统的单个结构示意图，图1为成像系统孔径排布方式示例，本发明的干涉成像方法通过利用微透镜阵列10接收光信号，通过光子集成芯片20进行信号的分频干涉及恢复算法等操作，最后通过光子集成芯片20上的pd获得光强信号，该信号即可解算观测目标的频域信息。该系统结构能够获取观测目标频域信息的解算信号。
[0142]
观测目标频域信息的解算信号的获取包括如下步骤：
[0143]
步骤1.1：光信号从物体平面传播到孔径平面
[0144]
如图6所示，物体平面的光信号传播到孔径平面，这个过程可以应用部分相干理论描述，可以得到孔径平面上的互谱密度函数(mutualspectraldensityfunction，msdf)观测目标的傅里叶变化。该自由传播过程不受系统配置的影响。
[0145]
具体地，远场非相干源的光传播到孔径平面，这个过程可以应用部分相干理论描述，孔径平面上的互谱密度函数(mutualspectraldensityfunction，msdf)可以表示为：
[0146][0147]
其中空间频率为α'表示观测平面上某一点的坐标；表示孔径平面上任意一点的坐标，观测光的频率为ν、波长为λ，非相干源的归一化频谱可以表示为g(ν)，并有∫g(ν)dν＝1，观测目标光强分布为i
object
(α')。孔径平面上，任意两点之间的互谱密度函数只与孔径平面上的两个干涉位置差有关，与物体的绝对位置无关。互谱密度函数是物体信号的傅里叶变换，f1‑
f2越大可以观测目标越高频的信号，因此孔径平面上的两个干涉位置差越大就可以接收到更高频的物体信号，从而得到更好的细节信息，低频包含了物体的轮廓信息。
[0148]
步骤1.2：光信号从孔径平面到光子集成芯片20的光电探测器24
[0149]
如图2所示，光信号从孔径平面到光子集成芯片20的pd(光电探测器24)传播，微透镜阵列10将光汇聚到光子集成芯片20上的光栅阵列21即vhc结构上(垂直转水平耦合，vertical
‑
to
‑
horizontal converter，是一种片上光栅结构：surface grating)，vhc结构
可以将汇聚光耦合入光子集成芯片20内部，需要经过光子集成芯片20上的输入波导(surfacegrating)阵列、波分解复用结构22(例如arrayedwaveguidegrating结构)、频谱信号解算结构23(如图9，包含移相结构231、波导耦合结构232)，最后被光子集成芯片20的光电探测器24接收，并转换为电信号。此部分过程可以描述为：
[0150]
i
k
＝|m
k
|2＝s2 l2 2slcos[(φ
s
‑
φ
l
) φ
k
]
[0151]
其中i
k
表示光子集成芯片20输出的光强信号，m
k
表示光子集成芯片20输出的振幅。s、φ
s
表示其中一个输入波导的幅值和相位，即l、φ
l
表示另一个输入波导的幅值和相位，即φ
k
表示光子集成芯片20引入的相位差，为了解算观测目标的频率信息，φ
k
的值通过光子集成芯片20可设为
[0152]
注:文中加粗的量表示矢量，具有数值大小以及方向的量。
[0153]
具体地，如图2所示，光从孔径平面到光子集成芯片20的pd传播，需要经过微透镜，每个微透镜将光汇聚到其对应的后端光子集成芯片20上的输入波导(surfacegrating、vhc结构)阵列，经过片上波分解复用结构22，将宽谱段光转换为多个准单色光λ
k
输出。径向上两个不同中心位置(对应基线b
ij
＝d
i
‑
d
j
)的微透镜得到准单色光λ
k
的光信号输入到同一个频谱信号解算结构23(如图9，包含移相结构231、波导耦合结构232)，既可以得到该基线下b
ij
的四个有不同相位差的信号，四个信号的相位差每个信号最后被对应光子集成芯片20上的光电探测器24接收，并转换为电信号。单个空间频率可以得到的四个光强信号如下，此部分过程可以描述为：
[0154]
i
k
＝|m
k
|2＝s2 l2 2slcos[(φ
s
‑
φ
l
) φ
k
]
[0155]
其中i
k
(k＝1,2,3,4)表示单个频谱信号解算结构23的光强信号，这些信号对应于确定的基线b
ij
，波长λ
k
及空间频率f。根据傅里叶光学，频域信息是空域信息的傅里叶变换，高空间频率信息代表了目标的细节信息部分，低空间频率信息代表了目标的轮廓信息，低频和高频对成像都很重要，因此需要增加基线b
ij
数量及波分解复用得到的准单色光λ
k
的个数，来增加对不同空间频率f对应频域信号v(f)的采样。m
k
表示光子集成芯片20输出的振幅。s，φ
s
表示其中一个输入波导的幅值和相位，即l，φ
l
表示另一个输入波导的幅值和相位，即φ
k
表示光子集成芯片20中频谱信号解算结构23引入的相位差，为了解算观测目标的频率信息，φ
k
的值通过光子集成芯片20可设为
[0156]
步骤1.3：微透镜阵列10排布与频域信息采样
[0157]
不同的孔径阵列排布会得到不同观测目标的频域信号分布，通过合适的设计可以得到较好的频域信号分布，本发明为了展示该系统的可行性，选择如图7的孔径排布方式(采用和spider相同的孔径排布方式)，光子集成芯片20设置在背板26上。频域得到的观测目标的空间频率f与基线b的关系可以表示为：
[0158]
f＝λz/(d
i
‑
d
j
)
[0159]
＝λz/b
[0160]
λ表示成像波长，z表示物体平面到孔径平面的距离。b＝b
i,j
＝d
i
‑
d
j
表示径向上(如图7)两个不同的微透镜中心位置之差，称为基线 (baseline)。空间频率f表示频域坐标，α'表示空域坐标。观测目标 i
object
(α')通过fourier transform(ft)即可得到频域分布信息，频谱信息可表示为：
[0161]
v(f)＝ft{i
object
(α')}
[0162]
其频域的采样方式可以表示为：
[0163][0164]
其中h(f)表示在空间频率f处的对频谱信息v(f)的采样权重。
[0165]
具体地，不同的孔径阵列排布会得到不同的观测目标的频域信号分布，通过合适的设计可以得到较好的频域信号分布，本发明为了展示该系统的可行性，选择如图7的孔径排布方式(采用和spider相同的空间排布方式)。其有多个(spider设计的是37个)如图2所示的结构组合而成，将图2结构按照径向进行排布，即可得到如图7的孔径排布方式。由于采用离散的孔径排布，因此频域采样也是离散采样，下面对系统的各个参量进行离散处理，频域得到的观测目标的空间频率f与基线b的关系可以表示为：
[0166]
f＝λz/(d
‑
d)
[0167]
＝λz/b
[0168]
离散表示主要是两个部分，由于波分解复用结构22的存在可以将宽谱段的光转换为多个准单色光输出λ
n
(n＝1,2,3,...,t
n
)；每个微透镜阵列10和光子集成芯片20都是按照径向排布，因此干涉基线也是沿着径向的分布，用i(i＝1,2,3,4,...,t
i
)表示径向坐标，j(j＝1,2,3,...,t
j
)表示横向坐标，则d
i'j
、d
i”j
表示径向上两个不同微透镜的中心位置矢量，b
ij
＝d
i'j
‑
d
i”j
表示径向上两个不同的微透镜中心位置之差，称为基线(baseline)，由于采用的是两两干涉因此第j个径向上的基线b
ij
个数可以达到因此频域得到的空间频率采样中f
kij
与基线b
ij
的关系可以表示为：
[0169]
f
kij
＝λ
k
z/(d
i'j
‑
d
i”j
)
[0170]
＝λ
k
z/b
ij
[0171]
空间频率f
kij
是通过fourier transform(ft)之后空域坐标与频域坐标的转换关系，也表征观测目标细节信息分布情况；本发明用α'表征空域坐标。观测目标i
object
(α')通过fouriertransform(ft)即可得到频域分布信息，频谱信息可表示为：
[0172]
v(f)＝ft{i
object
(α')}
[0173]
其离散的表达式可表示为：
[0174]
v(f
kij
)＝ft{i
object
(α')}
[0175]
其频域的采样方式可以表示为：
[0176][0177]
其中h(f
kij
)表示在空间频率f
kij
处的采样权重。
[0178]
步骤2：目标图像重建
[0179]
通过光子集成芯片20上移相结构231的共相调整后，使得两路干涉的微透镜之间达到较好的共相关系，此时通过光子集成芯片20的pd可以得到观测目标频域信息的解算信号，pd只能接受光强信号，无法接受到振幅信号，因此需要在芯片上设计振幅信号的解算结构，如图9所示为片上频谱信号解算结构23(包含移相结构231、波导耦合结构232)的简单示意图，移相结构231可以采用片上热调移相或者电调移相技术，波导耦合接口包括yjunction和xcoupler结构来实现对光信号的分光干涉。通过片上频谱信号解算结构23可以得到单个波长特定基线下的四个光强信号i
k
(k＝1,2,3,4)，这四个光强信号可以解算出该基线b
ij
下对应的频谱信号由于微透镜阵列10排布稀疏，观测目标频域信息不能被充分采集到，因此如果直接进行inversefouriertransform(ift)恢复的效果会很差，需要借助特殊设计的算法，消除这些频域信息未充分采集对最终成像的影响，最终可以获得高分辨的图像。
[0180]
步骤2.1：解算观测目标的频域信息
[0181]
如图2，经过微透镜汇聚的光耦合到输入波导(surfacegrating)阵列，在经过波分解复用结构22可以得到多个准单色光λ
k
信号，这些信号通过片上频谱信号解算结构23对应基线长度b
ij
可以得到四个光强信号输出i
k
(k＝1,2,3,4)，这些光强信号通过abcd方法即可解算对应空间频率为f
nij
的频域信号：
[0182][0183][0184]
v(f)表示空间频率f对应的频谱信号，abs{v(f)}表示其幅值部分，arg{v(f)}表示其相位部分。
[0185]
步骤2.2：频域信息恢复观测目标
[0186]
本发明中的成像系统对观测目标频域信息降质影响可以简单表示为：
[0187]
v
sample
＝v
·
h
[0188]
离散情况下可以写为：
[0189]
v
sample
(f
kij
)＝v
object
(f
kij
)
·
h(f
kij
)
[0190]
其中v表示观测目标的频域信息分布，v(f)＝ft{i
object
(α')}是观测目标i
object
的fouriertransform(ft)，v
sample
表示采样得到的频域信息，
·
表示点乘操作。对于v、v
sample
离散的分量为v
ij
、降质过程实际上是和采样的离散分量h
ij
的对应值相乘得到:
[0191][0192]
频域信息恢复观测目标基本思想是对逆向问题的求解，其基本思路是对于已采样的频谱信息部分做最大程度的保留，即添加datafidelity项；对于未采样部分选择使用惩罚函数来添加对观测目标的先验，而且较好的惩罚函数(penaltyfunction)可以使病态问题存在唯一解，也可以让解得到快速的收敛。
[0193]
最终需要优化的组合函数包括两个部分datafidelity项α
·
χ2及penalty
function项β
·
p，即，恢复观测目标采用如下公式：
[0194]
j＝α
·
χ2 β
·
p
[0195][0196][0197]
v
sample
＝v(f)＝abs{v(f)}exp(i
·
arg{v(f)})
[0198]
其中，α，β表示各项的系数(此系数可以通过人为设定或者自更新得到，相关操作可以参考凸优化算法)；p为惩罚函数；表示对矩阵 m中的所有元素求平方和，即m中的所有元素求平方和，即表示对观测目标的估计值，a表示系统对观测目标的降质作用，v
sample
表示采样得到的频域信息；当j的值收敛到最小值，即可得
[0199]
惩罚函数p有很多形式，例如总变分惩罚，熵惩罚等。本发明中为了展示成像的可行性，采用最大熵(maximum entropy)惩罚，此外还添加了 tp(total power)限制即
[0200][0201]
其中x
ij
表示x中的分量。通过采用不同的优化算法计算迭代，使得组合函数j的值收敛到最小值，即可得到恢复的观测目标。因此可以得到其惩罚函数：
[0202][0203]
其中α
ij
、a
ij
、x
ij
、v
ijsample
、β'
ij
、β”ij
表示α、a、x、v
sample
、β'、β”的分量。x表示观测目标的估计值选用梯度下降进行x的更新可以得到：
[0204][0205]
通过多次的迭代即可得到频域信息恢复观测目标。
[0206]
实施例2
[0207]
根据本发明的另一实施例，提供了一种基于微透镜阵列与光子集成芯片干涉成像系统，参见图1
‑
10，包括：
[0208]
解算信号获取单元，用于利用微透镜阵列10接收光信号，同时采用 vhc的光耦合方式接收从微透镜阵列10出射的光信号，对光耦合后的光信号进行处理获取观测目标频域信息的解算信号：
[0209]
观测目标恢复单元，用于解算出解算信号中观测目标的频域信息，在频域信息中恢复观测目标。
[0210]
本发明实施例中的基于微透镜阵列与光子集成芯片干涉成像系统，利用微透镜阵列10接收光信号，同时采用vhc的光耦合方式接收从微透镜阵列10出射的光信号，对光耦合后的光信号进行处理获取观测目标频域信息的解算信号，解算出解算信号中观测目标的频域信息，在频域信息中恢复观测目标。本发明利用微透镜阵列10代替其分立的小透镜，共相调整难度可以得到很大程度的降低，从而进一步优化系统更有利于集成，提升系统的稳定性。采用基于微透镜阵列10与干涉成像系统还可以进一步降低系统的体积。采用vhc的耦合方式可以降低空间光耦合入芯片的对光难度，还可以大大压缩系统的长度，依靠超薄的微
透镜阵列10和光子集成芯片20即可完成成像，从而实现平板成像。
[0211]
其中，解算信号获取单元包括：
[0212]
排布单元，用于排布微透镜阵列10；
[0213]
采样单元，用于对频域信息进行采样；
[0214]
光信号汇聚单元，用于位于物体平面的光信号通过自由空间传播至微透镜阵列10所在的孔径平面，微透镜阵列10将光信号进行汇聚；
[0215]
光信号接收单元，用于将光信号从孔径平面传输到光子集成芯片20 的光电探测器24上；
[0216]
采样单元，用于对频域信息进行采样。
[0217]
下面以具体实施例，对本发明的基于微透镜阵列与光子集成芯片干涉成像系统进行详细说明：
[0218]
具体参见图1
‑
10，本发明使用阵列的微透镜形式代替每个分立的小透镜，其光路的调整难度大大降低，可以更有利于集成，提升系统的稳定性。微透镜是一种直径为微米级的小透镜，微透镜按一定的周期以一维或二维的方式排列在基底上构成微透镜阵列10。微透镜阵列10利用其表面波长量级的浮雕结构对光波进行调制，将一个完整的激光波前在空间上分成许多微小的部分，每一部分都被相应的小透镜聚焦在焦平面上，一系列微透镜就可以得到由很多均匀、规则排布的焦点组成的平面。将该结构引入干涉成像代替原来的分立透镜可以大大降低系统的调整难度。
[0219]
如图3所示为空间光耦合波导的两种形式：端面耦合以及垂直转水平耦合(vhc，vertical
‑
to
‑
horizontal converter，是一种片上光栅结构： surface grating)。可以将垂直或者有一定偏角照射的光转换到水平输出， spider采用端面耦合的形式其对空间光对准波导的角度精度以及位置精度要求很高，本发明提出采用vhc的耦合方式可以降低空间光耦合入芯片的对光难度，还可以大大压缩系统的长度，依靠超薄的微透镜阵列10和光子集成芯片20即可完成成像，从而实现平板成像。
[0220]
因此基于已有的精确折光传统大口径成像系统存在着体积、功耗和重量大的缺点，本发明提供了一种基于微透镜阵列与光子集成芯片干涉成像系统，使用多个小口径透镜代替大口径成像系统，利用波导作为导波和光干涉介质，可以大大减小系统的体积、功耗和重量；其次使用了光子集成芯片20的移相结构231组件(例如热控或者电控移相器)可以降低成像系统的共相调整难度。
[0221]
相比于spider系统，本发明利用微透镜阵列10代替其分立的小透镜，共相调整难度可以得到很大程度的降低，从而进一步优化系统更有利于集成，提升系统的稳定性。采用基于微透镜阵列与光子集成芯片干涉成像系统还可以进一步降低系统的体积。采用vhc的耦合方式可以降低空间光耦合入芯片的对光难度，还可以大大压缩系统的长度，依靠超薄的微透镜阵列10和光子集成芯片20即可完成成像，从而实现平板成像。
[0222]
本发明公开的一种基于微透镜阵列与光子集成芯片干涉成像系统，包括微透镜阵列10、光子集成芯片20，微透镜阵列10包括第一基板11以及若干个刻蚀在第一基板11上的微透镜列12，每个微透镜列12包括若干个微透镜，光子集成芯片20包括第二基板25以及集成在第二基板25上的光栅阵列21(输入波导阵列)、波分解复用结构22、频谱信号解算结构 23和光电探测器24，微透镜阵列10、光栅阵列21、波分解复用结构22、频谱信号解算结构23
和光电探测器24沿光路依次设置，微透镜阵列10用于将入射光波汇聚，光栅阵列21用于将垂直或者有一定偏角照射的入射光波转换到水平输出，经过微透镜汇聚的光波耦合入光栅阵列21后通过波分解复用结构22和频谱信号解算结构23实现分频、相位调整以满足干涉条件，光电探测器24用于获取频谱信号解算结构23输出的光强信号(即观测目标的频域信息解算信号)并将光强信号转换为电信号，并对电信号进行计算处理输出重建图像。
[0223]
本发明公开的一种基于微透镜阵列与光子集成芯片干涉成像系统中，远场目标发出的光经过自由传播到达微透镜阵列10，经过每个微透镜汇聚后耦合入每个微透镜后端对应的光子集成芯片20的输入波导(surfacegrating)阵列中(如图2所示，每个输入波导都有一定的视场，将输出波导做成阵列结构可以增加其成像视场，光子集成芯片20利用波分解复用结构22(例如arrayedwaveguidegrating结构)、频谱信号解算结构23(如图9，包含移相结构231、波导耦合结构232)对耦合入芯片光分别分频、移相和干涉操作，光子集成芯光子集成芯片20的光电探测器24(photoelectricdetector，pd)将光信号转换为电信号输出，可得到光强值，从这些光强值可以计算出观测目标的频域信息的幅值和相位值。
[0224]
本发明公开的一种基于微透镜阵列与光子集成芯片干涉成像系统，其整体流程如图1所示，具体包括如下步骤：
[0225]
步骤1：解算信号获取单元：获取观测目标频域信息的解算信号
[0226]
如图2所示该系统的单个结构示意图，图1为成像系统孔径排布方式示例，本发明的干涉成像方法通过利用微透镜阵列10接收光信号，通过光子集成芯片20进行信号的分频干涉及恢复算法等操作，最后通过光子集成芯片20上的pd获得光强信号，该信号即可解算观测目标的频域信息。该系统结构能够获取观测目标频域信息的解算信号。
[0227]
观测目标频域信息的解算信号的获取包括如下步骤：
[0228]
步骤1.1：光信号汇聚单元：光信号从物体平面传播到孔径平面
[0229]
如图6所示，物体平面的光信号传播到孔径平面，这个过程可以应用部分相干理论描述，可以得到孔径平面上的互谱密度函数(mutualspectraldensityfunction，msdf)观测目标的傅里叶变化。该自由传播过程不受系统配置的影响。
[0230]
具体地，远场非相干源的光传播到孔径平面，这个过程可以应用部分相干理论描述，孔径平面上的互谱密度函数(mutualspectraldensityfunction，msdf)可以表示为：
[0231][0232]
其中空间频率为α'表示观测平面上某一点的坐标；表示孔径平面上任意一点的坐标，观测光的频率为ν、波长为λ，非相干源的归一化频谱可以表示为g(ν)，并有∫g(ν)dν＝1，观测目标光强分布为i
object
(α')。孔径平面上，任意两点之间的互谱密度函数只与孔径平面上的两个干涉位置差有关，与物体的绝对位置无关。互谱密度函数是物体信号的傅里叶变换，f1‑
f2越大可以观测目标越高频的信号，因此孔径平面上的两个干涉位置差越大就可以接收到更高频的物体信号，从而得到更好的细节信息，低频包含了物体的轮廓信息。
[0233]
步骤1.2：光信号接收单元：光信号从孔径平面到光子集成芯片20的光电探测器24
[0234]
如图2所示，光信号从孔径平面到光子集成芯片20的pd(光电探测器 24)传播，需
要经过光子集成芯片20上的输入波导(surface grating) 阵列、波分解复用结构22(例如arrayed wave guide grating结构)、频谱信号解算结构23(如图9，包含移相结构231、波导耦合结构232)，最后被光子集成芯片20的光电探测器24接收，并转换为电信号。此部分过程可以描述为：
[0235]
i
k
＝|m
k
|2＝s2 l2 2slcos[(φ
s
‑
φ
l
) φ
k
]
[0236]
其中i
k
表示光子集成芯片20输出的光强信号，m
k
表示光子集成芯片 20输出的振幅。s、φ
s
表示其中一个输入波导的幅值和相位，即l、φ
l
表示另一个输入波导的幅值和相位，即φ
k
表示光子集成芯片20 引入的相位差，为了解算观测目标的频率信息，φ
k
的值通过光子集成芯片20可设为
[0237]
注:文中加粗的量表示矢量，具有数值大小以及方向的量。
[0238]
具体地，如图2所示，光从孔径平面到光子集成芯片20的pd传播，需要经过微透镜，每个微透镜将光汇聚到其对应的后端光子集成芯片20 上的输入波导(surface grating、vhc结构)阵列，经过片上波分解复用结构22，将宽谱段光转换为多个准单色光λ
k
输出。径向上两个不同中心位置(对应基线b
ij
＝d
i
‑
d
j
)的微透镜得到准单色光λ
k
的光信号输入到同一个频谱信号解算结构23(如图9，包含移相结构231、波导耦合结构232)，既可以得到该基线下b
ij
的四个有不同相位差的信号，四个信号的相位差每个信号最后被对应光子集成芯片20上的光电探测器 24接收，并转换为电信号。单个空间频率可以得到的四个光强信号如下，此部分过程可以描述为：
[0239]
i
k
＝|m
k
|2＝s2 l2 2slcos[(φ
s
‑
φ
l
) φ
k
]
[0240]
其中i
k
(k＝1,2,3,4)表示单个频谱信号解算结构23的光强信号，这些信号对应于确定的基线b
ij
，波长λ
k
及空间频率f。根据傅里叶光学，频域信息是空域信息的傅里叶变换，高空间频率信息代表了目标的细节信息部分，低空间频率信息代表了目标的轮廓信息，低频和高频对成像都很重要，因此需要增加基线b
ij
数量及波分解复用得到的准单色光λ
k
的个数，来增加对不同空间频率f对应频域信号v(f)的采样。m
k
表示光子集成芯片20输出的振幅。s，φ
s
表示其中一个输入波导的幅值和相位，即l，φ
l
表示另一个输入波导的幅值和相位，即φ
k
表示光子集成芯片20中频谱信号解算结构23引入的相位差，为了解算观测目标的频率信息，φ
k
的值通过光子集成芯片20可设为
[0241]
步骤1.3：排布单元及采样单元：微透镜阵列10排布与频域信息采样
[0242]
不同的孔径阵列排布会得到不同观测目标的频域信号分布，通过合适的设计可以得到较好的频域信号分布，本发明为了展示该系统的可行性，选择如图7的孔径排布方式(采用和spider相同的孔径排布方式)。频域得到的观测目标的空间频率f与基线b的关系可以表示为：
[0243]
f＝λz/(d
i
‑
d
j
)
[0244]
＝λz/b
[0245]
λ表示成像波长，z表示物体平面到孔径平面的距离。b＝b
i,j
＝d
i
‑
d
j
表示径向上(如图7)两个不同的微透镜中心位置之差，称为基线 (baseline)。空间频率f表示频域坐标，α'表示空域坐标。观测目标 i
object
(α')通过fourier transform(ft)即可得到频域分布信息，频谱信息可表示为：
[0246]
v(f)＝ft{i
object
(α')}
[0247]
其频域的采样方式可以表示为：
[0248][0249]
其中h(f)表示在空间频率f处的对频谱信息v(f)的采样权重。
[0250]
具体地，不同的孔径阵列排布会得到不同的观测目标的频域信号分布，通过合适的设计可以得到较好的频域信号分布，本发明为了展示该系统的可行性，选择如图7的孔径排布方式(采用和spider相同的空间排布方式)。其有多个(spider设计的是37个)如图2所示的结构组合而成，将图2结构按照径向进行排布，即可得到如图7的孔径排布方式。由于采用离散的孔径排布，因此频域采样也是离散采样，下面对系统的各个参量进行离散处理，频域得到的观测目标的空间频率f与基线b的关系可以表示为：
[0251]
f＝λz/(d
‑
d)
[0252]
＝λz/b
[0253]
离散表示主要是两个部分，由于波分解复用结构22的存在可以将宽谱段的光转换为多个准单色光输出λ
n
(n＝1,2,3,...,t
n
)；每个微透镜阵列10和光子集成芯片20都是按照径向排布，因此干涉基线也是沿着径向的分布，用i(i＝1,2,3,4,...,t
i
)表示径向坐标，j(j＝1,2,3,...,t
j
)表示横向坐标，则d
i'j
、d
i”j
表示径向上两个不同微透镜的中心位置矢量，b
ij
＝d
i'j
‑
d
i”j
表示径向上两个不同的微透镜中心位置之差，称为基线(baseline)，由于采用的是两两干涉因此第j个径向上的基线b
ij
个数可以达到因此频域得到的空间频率采样中f
kij
与基线b
ij
的关系可以表示为：
[0254]
f
kij
＝λ
k
z/(d
i'j
‑
d
i”j
)
[0255]
＝λ
k
z/b
ij
[0256]
空间频率f
kij
是通过fourier transform(ft)之后空域坐标与频域坐标的转换关系，也表征观测目标细节信息分布情况；本发明用α'表征空域坐标。观测目标i
object
(α')通过fouriertransform(ft)即可得到频域分布信息，频谱信息可表示为：
[0257]
v(f)＝ft{i
object
(α')}
[0258]
其离散的表达式可表示为：
[0259]
v(f
kij
)＝ft{i
object
(α')}
[0260]
其频域的采样方式可以表示为：
[0261][0262]
其中h(f
kij
)表示在空间频率f
kij
处的采样权重。
[0263]
步骤2：目标图像重建
[0264]
通过光子集成芯片20上移相结构231的共相调整后，使得两路干涉的微透镜之间
达到较好的共相关系，此时通过光子集成芯片20的pd可以得到观测目标频域信息的解算信号，pd只能接受光强信号，无法接受到振幅信号，因此需要在芯片上设计振幅信号的解算结构，如图9所示为片上频谱信号解算结构23(包含移相结构231、波导耦合结构232)的简单示意图，移相结构231可以采用片上热调移相或者电调移相技术，波导耦合接口包括yjunction和xcoupler结构来实现对光信号的分光干涉。通过片上频谱信号解算结构23可以得到单个波长特定基线下的四个光强信号i
k
(k＝1,2,3,4)，这四个光强信号可以解算出该基线b
ij
下对应的频谱信号由于微透镜阵列10排布稀疏，观测目标频域信息不能被充分采集到，因此如果直接进行inversefouriertransform(ift)恢复的效果会很差，需要借助特殊设计的算法，消除这些频域信息未充分采集对最终成像的影响，最终可以获得高分辨的图像。
[0265]
步骤2.1：解算观测目标的频域信息
[0266]
如图2，经过微透镜汇聚的光耦合到输入波导(surfacegrating)阵列，在经过波分解复用结构22可以得到多个准单色光λ
k
信号，这些信号通过片上频谱信号解算结构23对应基线长度b
ij
可以得到四个光强信号输出i
k
(k＝1,2,3,4)，这些光强信号通过abcd方法即可解算对应空间频率为f
kij
的频域信号：
[0267][0268][0269]
v(f)表示空间频率f对应的频谱信号，abs{v(f)}表示其幅值部分，arg{v(f)}表示其相位部分。
[0270]
步骤2.2：频域信息恢复观测目标
[0271]
本发明中的成像系统对观测目标频域信息降质影响可以简单表示为：
[0272]
v
sample
＝v
·
h
[0273]
离散情况下可以写为：
[0274]
v
sample
(f
kij
)＝v
object
(f
kij
)
·
h(f
kij
)
[0275]
其中v表示观测目标的频域信息分布，v(f)＝ft{i
object
(α')}是观测目标i
object
的fouriertransform(ft)，v
sample
表示采样得到的频域信息，
·
表示点乘操作。对于v、v
sample
离散的分量为v
ij
、降质过程实际上是和采样的离散分量h
ij
的对应值相乘得到:
[0276][0277]
频域信息恢复观测目标基本思想是对逆向问题的求解，其基本思路是对于已采样的频谱信息部分做最大程度的保留，即添加datafidelity项；对于未采样部分选择使用惩罚函数来添加对观测目标的先验，而且较好的惩罚函数(penaltyfunction)可以使病态问题存在唯一解，也可以让解得到快速的收敛。
[0278]
最终需要优化的组合函数包括两个部分datafidelity项α
·
χ2及penaltyfunction项β
·
p，恢复观测目标采用如下公式：
[0279]
j＝α
·
χ2 β
·
p
[0280][0281][0282]
v
sample
＝v(f)＝abs{v(f)}exp(i
·
arg{v(f)})
[0283]
其中，α，β表示各项的系数，表示对矩阵m中的所有元素求平方和，即表示对矩阵m中的所有元素求平方和，即表示对观测目标的估计值，a表示系统对观测目标的降质作用，v
sample
表示采样得到的频域信息；当j的值收敛到最小值，即可得；惩罚函数p有很多形式，例如总变分惩罚，熵惩罚等。本发明中为了展示成像的可行性，采用最大熵(maximum entropy)惩罚，此外还添加了tp(total power)限制即
[0284][0285]
其中x
ij
表示x中的分量。通过采用不同的优化算法计算迭代，使得组合函数j的值收敛到最小值，即可得到恢复的观测目标。因此可以得到其惩罚函数：
[0286][0287]
其中α
ij
、a
ij
、x
ij
、v
ijsample
、β'
ij
、β”ij
表示α、a、x、v
sample
、β'、β”的分量。x表示观测目标的估计值选用梯度下降进行x的更新可以得到：
[0288][0289]
通过多次的迭代即可得到频域信息恢复观测目标。
[0290]
本发明的仿真条件为：
[0291]
本发明在inter(r)cpu3.00ghz，8g的pc机上，使用matlabr2018a 平台，对不同特征的目标图像进行了成像仿真。表一是仿真的系统配置。目标图像重建选用最大熵方法(maximum entropy method)。
[0292][0293][0294]
表1、仿真系统相关参数
[0295]
其仿真结果如下：
[0296]
图1为本发明对不同特征在表1的系统配置下的成像结果。利用峰值信噪比(peak signal
‑
to
‑
noise ratio，psnr)和结构相似度指数 (structural similarity index，ssim)进一步评估恢复的目标分布与观测目标之间的相似性。结果如表2所示，其中添加的噪声是均值为0、方差为0.004的高斯噪声。可以看出，恢复后的图像更接近观测目标值。
[0297][0298]
表二、干涉成像系统成像的像质评价
[0299]
从图1可以看出，本发明可以使用稀疏的子孔径干涉成像即可实现二维目标成像，且对不同特征的目标图像重建质量都很高。本发明的技术方案已经过仿真实验验证，效果好。
[0300]
本发明具有如下有益效果：
[0301]
(1)成像质量高
[0302]
由仿真结果可以看出，重建图像的ssim(结构相似性系数)均达到 0.890以上，psnr(峰值信噪)可以达到12.8db以上。相比于spider(用于光电侦察的分段平面成像探测器)系统，本发明将高精度的超表面加工技术及光子芯片技术组合，将孔径平面的原本的分立小透镜转换为集成在同一个基板上的微透镜可以很大程度上提升系统的集成度。
[0303]
(2)系统的调整难度低
[0304]
本发明对比传统的基于精确折光的成像系统而言，其共相调整是借助于波导的移相器件(例如热控或者电控移相器)，可以很高的调整精度而且调整方便。目前波导的热控移相器以及光纤移相器可以达到很高的移相精度，借助干涉的信号以及这些移相器件可以得到高精度的共相调整。相比于spider(用于光电侦察的分段平面成像探测器)系统，本发明将高精度的超表面加工技术及光子芯片技术组合，大大减低了芯片与聚光系统 (即微透镜阵列10)之间光路调整调整量，从而减小光路调整难度。
[0305]
(3)体积、功耗和重量上的优势
[0306]
与具有类似有效孔径和空间分辨率的传统望远镜相比，本发明多孔径干涉成像技术由于采用多个小孔径以及光子波导(光纤和光子集成芯片20 等)在整个系统的尺寸、重量和功率(size，weight and power，swap)可以得到10倍及以上的减小，对于应用于车载或者星载，可以大大减少载荷的负担。采用vhc的耦合方式可以降低空间光耦合入芯片的对光难度，还可以大大压缩系统的长度，依靠超薄的微透镜阵列10和光子集成芯片 20即可完成成像，从而实现平板成像。
[0307]
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
[0308]
在本发明的上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见其他实施例的相关描述。
[0309]
在本技术所提供的几个实施例中，应该理解到，所揭露的技术内容，可通过其它的
方式实现。其中，以上所描述的系统实施例仅仅是示意性的，例如单元的划分，可以为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，单元或模块的间接耦合或通信连接，可以是电性或其它的形式。
[0310]
作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
[0311]
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
[0312]
集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例方法的全部或部分步骤。而前述的存储介质包括：u盘、只读存储器(rom，read
‑
onlymemory)、随机存取存储器(ram，random access memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
[0313]
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。