[发明专利]基于旋转编码光圈成像系统的深度估计方法

权利要求书

1.一种基于旋转编码光圈成像系统的深度估计方法，其特征是，包括如下步骤：提出旋转编码光圈成像系统的性能评价标准，并依据该评价标准使用遗传算法及坐标下降法优化旋转编码光圈，后借助优化的旋转编码光圈图案实现深度估计及全聚焦图像复原工作。

2.如权利要求1所述的基于旋转编码光圈成像系统的深度估计方法，其特征是，提出旋转编码光圈成像系统的性能评价标准具体步骤是：

(1)失焦模糊成像模型：

f=f0⊗kd+η]]>

其中，f代表成像系统拍摄得到的失焦模糊图像，f₀代表原始清晰图像，k^d代表成像系统在不同深度d处的点扩散函数，描述高斯白噪声，X表示高斯分布，表示高斯噪声的标准差；

(2)三次旋转编码光圈系统的失焦成像模型：

F_i＝F₀·K_i^d+N_i,i＝1,2,3

其中，此模型是(1)中模型在频域中的表达：F_i表示旋转编码光圈成像系统拍摄得到的失焦模糊图像，F₀表示原始图像，K_i^d表示不同深度的点扩散函数，N_i表示第i次旋转编码光圈获得拍摄图像中的高斯白噪声，i表示旋转次数，即三次旋转会产生三次不同的失焦模糊图像；此模型对应于旋转编码光圈成像系统的拍摄所得的三张失焦模糊图像的数学模型；

(3)使用能量函数最小值解决最大后验概率问题(MaximumAPosteriori)以实现对旋转编码光圈的性能优化。

E(d^|F1,F2,F3,σ)=minF0Σi=1,2,3||F^0·Kid^-Fi||2+||C·F^0||2]]>

其中，代表图像先验的噪信比(NoisetoSignalRatios)，可以表示为σ²/A，其中，A使用多张自然图像的平均功率谱来表示，ξ表示频率，使用能量函数最小值意味着使估计模糊图像与原始模糊图像之间的差别最小，即使得估计过程精度最高，深度恢复结果最准确；

(4)使用维纳反卷积法获得估计清晰图像。即对一个既定的估计深度使用能量函数先求解估计清晰图像即在的条件下，可得，

F^0=ΣiFr·K‾id^Σi|Kid^|2+|C|2,i=1,2,3]]>

其中，是K的复共轭，是指在既定估计深度处对应的成像系统点扩散函数的频域表示，本发明中的字母上方的∧均表示该字母所代表变量是一个计算估计值；

(5)将步骤(4)和步骤(2)中的公式带入步骤(3)，重构能量函数，经整理简化得到一个关于估计深度的能量评价函数，

E(d^|K1d*,K2d*,K3d*,σ)=ΣξA|K1dK2d*-K2dK1d*|2+|K1dK3d*-K3dK1d*|2+|K2dK3d*-K3dK2d*|2|K1d|2+|K2d|2+|K3d|2+C2+σ2·Σξ[C2|K1d|2+|K2d|2+|K3d|2+C2+1]]]>

其中，d^*表示位于成像系统所对应的标准深度，d表示可见光图像的各层级的处深度；

(6)在空域中的加性高斯白噪声参数是σ的条件下，对旋转编码光圈位于d^*处的深度估计进行性能评价可得，

R(K1,K2,K3|d*,σ)=mind∈D/d*E(d|K1d*,K2d*,K3d*,σ)-E(d*|K1d*,K2d*,K3d*,σ)]]>

其中，D＝{c₁d^*,c₂d^*,…,c_ld^*}代表系统对应的各深度层级，{c_i}可取为{0.1,0.15,…,1.5}，步长为0.05。

(7)对以上结果进行归一化可得旋转编码光圈的性能评价标准：

M(K1,K2,K3,d,d*)=[1nΣξ|K1dK2d*-K2dK1d*|2+|K1dK3d*-K3dK1d*|2+|K2dK3d*-K3dK2d*|2|K1d|2+|K2d|2+|K3d|2+C2]12R=mind∈D/d*M(K1,K2,K3,d,d*)]]>

R值越大表征评价深度估计的能量函数变化越陡峭，使得深度估计性能对于图像噪声或者彩色纹理的鲁棒性越强。