[发明专利]星空图像坏点鲁棒检测方法

摘要

本发明公开了一种星空图像坏点鲁棒检测方法，用于解决现有图像坏点检测方法局限性大的技术问题。技术方案是充分考虑星空图像的特性和坏点在图像像素统计中的特性，利用星空图像局部区域统计稳定的特性、中值统计量和平均绝对离差统计量等鲁棒统计量排除样本中的外点对统计量估计的干扰，并将外点位置的像素当做局部像素的统计外点，设计基于统计外点检测思想的统计方法，并利用连续序列图像中充分的样本对每个位置进行判断从而检测出图像中外点。该方法参数少且统计意义明确，易于根据具体应用调整，并且适用于不同程度坏点的检测。

主权项

一种星空图像坏点鲁棒检测方法，其特征在于包括以下步骤：步骤一、首先在多帧图像序列S_I＝{I₁,I₂,...,I_t,...I_T}中随机选取一帧图像I_R作为参考图像,其中t＝1,...,T，T表示序列中图像帧数,I_t表示任一帧图像；对于该单帧参考星空图像I_R，利用阈值对星空图像进行二值化分割，即根据像素值是否大于阈值ε_bw而进行分割，图像I_t中(i,j)处像素点灰度值为且有1≤i≤h,1≤j≤w，图像中像素灰度值满足的像素被分割为星点光斑区域，而满足的像素点分割为背景区域；通过对星图数据进行鲁棒统计得到分割阈值ε_bw＝μ_bw+δ_bwσ_bw           (1)其中，μ_bw为图像灰度中值，σ_bw为图像平均绝对偏差δ_bw为阈值参数，δ_bw＝6；利用二值化图像对星体区域进行剔除，得到星空图像背景像素点集合其中i_B＝1,...,N_B,表示背景像素集合中的单个像素值，N_B表示背景像素点集合中的像素个数，该数目与二值化分割结果相关；利用剩余的背景像素对星空图像背景噪声程度进行鲁棒估计，对B计算统计中值得到背景灰度统计中值μ_B,并计算B的平均绝对离差这两个值将作为后续坏点检测的参数；步骤二、对星空图像中的坏点位置进行检测；该过程是一个滑动窗口的过程，以每个位置的多帧中像素点及其在多帧中的邻域内像素点进行鲁棒统计外点检测，通过判断该像素点是否为外点进而判断该点是否为图像中的坏点；首先人工设定一个局部选取尺寸参数值为d,以此参数对原图像序列S_I＝{I₁,I₂,...,I_t,...I_T}进行边缘扩展得到扩展图像集合其中t＝1,...,T，T表示序列中图像帧数，为原图像I_t边缘扩展后得到扩展后图像，即若原图像I_t为$I_t=\left[\begin{array}{ccc}{x}_{\mathrm{1,1}}& ...& x_{1,w}\\ ...& ...& ...\\ x_{h,1}& ...& x_{h,w}\end{array}\right]---\left(2\right)$其中，h为图像高度，w为图像宽度，x_i,j为图像中的像素，有1≤i≤h,1≤j≤w；扩展后的对应图像为$I_t^{\mathrm{ep}}=\left[\begin{array}{ccccccc}{x}_{d+1,d+1}& ...& x_{d+\mathrm{1,1}}& ...& x_{d+1,w}& ...& x_{d+1,w-d}\\ ...& ...& ...& ...& ...& ...& ...\\ x_{\mathrm{1,1}+d}& ...& x_{\mathrm{1,1}}& ...& x_{1,w}& ...& x_{1,w-d}\\ ...& ...& ...& ...& ...& ...& ...\\ x_{h,1+d}& ...& x_{h,1}& ...& x_{h,w}& ...& x_{h,w-d}\\ ...& ...& ...& ...& ...& ...& ...\\ x_{h-d,1+d}& ...& x_{h-d,1}& ...& x_{h-d,w}& ...& x_{h-d,w-d}\end{array}\right]---\left(3\right)$其中，扩展部分看作是对原图像以边缘进行镜面对称扩展的结果；扩展后图像的高度为h+2d,宽度为w+2d；用(i′,j′)表示扩展后图像的序号，有1≤i′≤(h+2d)与1≤j′≤(w+2d)，中位于(i′,j′)处的像素为根据镜面对称扩展关系，中像素与I_t中像素对应关系表示为：$I_t^{\mathrm{ep}}(i^{\′},j^{\′})=\left\{\begin{array}{c}{I}_t^{(d-i^{\′}+1,d-j^{\′}+1)},1\≤i^{\′}\≤d,1\≤j^{\′}\≤d\\ I_t^{(d-i^{\′}+1,j^{\′})},1\≤i^{\′}\≤d,d+1\≤j^{\′}\≤w+d\\ I_t^{(d-i^{\′}+\mathrm{1,2}h+d-j^{\′}-1)},1\≤i^{\′}\≤d,h+d+1\≤j^{\′}\≤w+2d\\ I_t^{(i^{\′}-d,d-j^{\′}+1)},d+1\≤i^{\′}\≤h+d,1\≤j^{\′}\≤d\\ I_t^{(i^{\′}-d,j^{\′}-d)},d+1\≤i^{\′}\≤h+d,d+1\≤j^{\′}\≤w+d\\ I_t^{(i^{\′}-d,2h+d-j^{\′}-1)},d+1\≤i^{\′}\≤h+d,w+d+1\≤j^{\′}\≤w+2d\\ I_t^{(2w+d-i^{\′}-1,d-j^{\′}+1)},h+d+1\≤i^{\′}\≤h+2d,1\≤j^{\′}\≤d\\ I_t^{(2w+d-i^{\′}-1,j^{\′}-d)},h+d+1\≤i^{\′}\≤h+2d,d+1\≤j^{\′}\≤w+d\\ I_t^{(2w+d-i^{\′}-\mathrm{1,2}h+d-j^{\′}-1)},h+d+1\≤i^{\′}\≤h+2d,w+d+1\≤j^{\′}\≤w+2d\end{array}\right.---\left(4\right)$中对应于原图像的I_t中的像素为满足d+1≤i′≤h+d,d+1≤j′≤w+d范围内的像素；通过对图像进行扫描检测坏点，具体过程为：(a)对于对d+1≤i′≤h+d且1+d≤j′≤w+d范围内的像素进行检测，即对对应于原图像中像素的位置的像素进行扫描；对于任意位置(i′,j′)处的像素，在每帧中以相同位置为中心各选取一个长宽大小均为d′＝2×d+1的正方形局部图像块，得到一个局部图像块序列$S_L^{(i^{\′},j^{\′})}=\{L_1^{(i^{\′},j^{\′})},L_2^{(i^{\′},j^{\′})},...,L_t^{(i^{\′},j^{\′})},...,L_T^{(i^{\′},j^{\′})}\},$对序列中图像(i′,j′)处的像素是否为坏点的判断将依据中的数据；(b)对中心位置(i′,j′)处像素统计特性进行鲁棒估计；$S_L^{(i^{\′},j^{\′})}=\{L_1^{(i^{\′},j^{\′})},L_2^{(i^{\′},j^{\′})},...,L_t^{(i^{\′},j^{\′})},...,L_T^{(i^{\′},j^{\′})}\}$中的每个局部图像块中的中心像素值，即位于每个扩展图像中(i′,j′)处或每个原图像I_t中(i′‑d,j′‑d)位置处的像素值，组成的序列$S_c^{(i^{\′},j^{\′})}=\{c_1^{(i^{\′},j^{\′})},c_2^{(i^{\′},j^{\′})},...,c_t^{(i^{\′},j^{\′})},...,c_T^{(i^{\′},j^{\′})}\},$其中为图像块中心像素值，即扩展后图像中(i′,j′)处的像素值；对像素集合计算统计中值得到中心点像素的统计中值μ_(i′,j′),并计算的平均绝对离差这两个值作为通过多帧数据鲁棒统计得到的(i′,j′)处的像素统计特性；(c)基于)对中心位置(i′,j′)邻域像素统计特性进行鲁棒估计；$S_L^{(i^{\′},j^{\′})}=\{L_1^{(i^{\′},j^{\′})},L_2^{(i^{\′},j^{\′})},...,L_t^{(i^{\′},j^{\′})},...,L_T^{(i^{\′},j^{\′})}\}$中每个局部图像块包含d′×d′个像素，除去中心像素之后中所有像素个数为T′＝d′×d′×T‑T，这些像素组成像素集合$S_a^{(i^{\′},j^{\′})}=\{a_1^{(i^{\′},j^{\′})},a_2^{(i^{\′},j^{\′})},...,a_{t^{\′}}^{(i^{\′},j^{\′})},...,a_{T^{\′}}^{(i^{\′},j^{\′})}\},$其中表示中序号为t′的像素值，序号t′满足1≤t′≤T′；对中样本计算统计中值得到邻域像素的统计中值并计算的平均绝对离差这两个值作为(i′,j′)处领域的像素统计特性；(d)根据已经估计得到的(i′,j′)处与邻域的像素值统计特性，基于统计外点检测方法判定(i′,j′)像素是否为坏点；μ_(i′,j′)与σ_(i′,j′)表示(i′,j′)处像素的统计量，而与表示(i′,j′)处领域像素的统计量，表示局部与内容相关的像素值的鲁棒统计量；为了判断(i′,j′)处像素是否为坏点，首先通过经验阈值对该位置像素统计量进行判断，若满足：μ_(i′,j′)∈[α,β]                (5)与表示该位置处的多帧中的像素点亮度不满足坏点具有的——亮度过亮或过暗且相邻多帧中同一位置处像素值一致性好——的特性，该点不被判断为坏点,其中α表示过暗坏点的亮度上限，β分别过亮坏点的亮度下限，表示多帧稳定的坏点像素的平均绝对离差的上限；利用步骤(a)中计算得到背景像素的统计值对参数进行赋值，α＝μ_B+θ₁σ_B,β＝μ_B‑θ₂σ_B,θ₁、θ₂和θ₃为根据经验给与的控制参数；若以上条件不满足，基于公式(7)判断(i′,j′)处像素是否为不符合局部观测像素值分布的统计外点：$|{\mathrm{\μ}}_{(i^{\′},j^{\′})}-{\mathrm{\μ}}_{(i^{\′},j^{\′})}^a|\<\mathrm{\λ}{\mathrm{\σ}}_{(i^{\′},j^{\′})}^a---\left(7\right)$其中，λ为置信度控制参数；若该条件满足，则(i′,j′)处像素被判定为不符合局部像素分布特性的统计外点，被判定为坏点；对中每个位置的像素重复以上的步骤(a)、步骤(b)和步骤(c),对每个位置处的像素进行判断；由于为S_I的镜面对称边缘扩展版本，在中检测得到坏点的位置(i′,j′)转换到原始图像中，位置为(i′‑d,j′‑d),d为镜面扩展参数；对图像序列S_I进行坏点检测得到其中i_D＝1,...,N_D,表示坏点集合中的单个坏点位置，N_D表示坏点集合中的坏点个数；N_D用来衡量图像传感器在当前情况下质量状态和工作情况。