【摘要】 本发明一种复杂背景下红外背景预测算子的自动选择方法，可概括为：首先对单帧图像中的每个点进行邻域操作得到以该点为中心的子图像邻域，并计算该子图像的平滑度，然后根据子图像的平滑度选择对应的背景预测算子，以获得背景图像。对具有不同图像平滑度的局部邻域自动采用不同的预测算子进行背景预测，比单独采用某一种预测算子计算得到的背景图像更接近原图，能够更大程度上的抑制背景噪声。 【专利类型】发明授权 【申请人】北京航空航天大学 【申请人类型】学校 【申请人地址】100191 北京市海淀区学院路37号北京航空航天大学电子信息工程学院 【申请人地区】中国 【申请人城市】北京市 【申请人区县】海淀区 【申请号】CN200810222206.3 【申请日】2008-09-11 【申请年份】2008 【公开公告号】CN101359402B 【公开公告日】2010-12-15 【公开公告年份】2010 【授权公告号】CN101359402B 【授权公告日】2010-12-15 【授权公告年份】2010.0 【IPC分类号】G06T7/20 【发明人】毛峡; 黄康; 胡海勇; 梁晓庚 【主权项内容】1.一种复杂背景下红外背景预测算子的自动选择方法，其特征在于，包括下述几个步骤： 步骤1：获取子图像邻域： 采用红外成像系统获得红外弱小目标图像并对其直接进行处理，即对红外成像系统产生的单帧图像进行处理，来获取子图像邻域，该红外成像系统的成像参数如下：光学系统有效面积A0＝12.0cm2，探测器星响应度 探测最小信噪比SNR＝4，探测器面积Ad＝0.4mm2，等效噪声带宽Δf＝2000.0Hz，光学系统红外辐射透过率τ0＝0.8，信号衰减因子δ≈1；“目标成像系统”和“辅助成像系统”的成像波段为红外中波波段3μm-5μm； 步骤2：计算局部图像的平滑度 复杂背景下红外弱小运动目标中的图像目标信息很弱，对获得的子图像邻域首先计算局部图像的平滑度，然后根据计算结果判断是否需要进一步计算局部图像灰度均值，以选择最适合的背景预测算子，具体步骤如下： 计算局部图像的平滑度： 平滑度用于描述图像灰度的偏移，设方差为v，则平滑度R表示为： 显然R∈[0，1]，且R越小表示图像越平滑，方差v的计算公式如下： 式中，L表示图像的灰度级，zi表示图像的灰度随机变量，P(zi)表示该灰度随机变量的概率，m表示局部图像灰度均值； 步骤3：根据局部图像的平滑度计算结果选择背景预测算子 3.1计算局部图像灰度均值m 当平滑度小于平滑度阈值TR时，需要计算局部图像灰度均值m，局部图像灰度均值m表示图像的整体灰度，其计算公式如下： 式中，L表示图像的灰度级，zi表示图像的灰度随机变量，P(zi)表示该灰度随机变量的概率； 3.2基本加权背景预测 当平滑度R大于平滑度阈值TR时，选择基本加权背景预测算子，基本加权背景预测算子计算较为简单，具体计算方法如式4-6所示； 所有的背景预测算子均由同一个背景预测的算法模型衍生出来，此背景预测的算法模型表示为： 式中，(i，j)表示图像中某点的坐标，Y(i，j)表示该点的背景预测值，S表示以该点为中心的邻域，X(i-l，j-k)表示原图像中该邻域中某点的像素值，Wn(l，k)表示该点的权重； 基本加权背景预测算子在背景预测算法模型的基础上定义了邻域各点的权重，如下所 示： 式中，Wn(l，k)表示邻域S中的某一点的权重，r(i，j)表示该点到中心点的距离，本发明采用7×7的模板进行背景预测，权重Wn(l，k)用矩阵形式Wn表示为： 将式6代入式4，即求出基本加权背景预测算子的计算结果； 步骤4：选择背景预测算子 根据局部图像平滑度和局部图像灰度均值的计算结果，在下列4种背景预测算子中选择合适的算子作为当前邻域的背景预测算子：即基本加权背景预测算子、最大化背景预测算子、最小化背景预测算子以及最相似背景预测算子； 以图像中某点为中心邻域S，以该点为中心划分为4个象限，用数字i表示为i＝1，2，3，4，根据基本加权背景预测的计算方法，将式6代入式4，计算4个象限对应的预测值，分别为：Y1，Y2，Y3，Y4，最小背景预测算子取其中最小值作为背景预测值，用Ymin表示最小背景预测算子计算得到的背景预测值，用于计算背景灰度均值较小区域的背景预测值，如下式所示： 最大背景预测算子取其中最大值作为背景预测值，用Ymax表示最大背景预测算子计算得到的背景预测值，如下式所示： 最相似背景预测算子取Y1～Y4中与原图像的点灰度值最接近的值作为背景预测值，如下式所示： 式9中abs(Yi-X)分别表示4个象限的预测值和原图像素值X差的绝对值； 如果当前领域S的局部图像灰度均值m小于Tm1，则中心点的灰度值低的可能性大，因此选择最小化背景预测算子；同理，如果当前领域S的局部图像灰度值m大于Tm2，则采用最大化背景预测算子；如果当前领域S的局部图像灰度值m属于[Tm1，Tm2]，则选择最相似背景预测算子； 根据蒙特卡洛实验结果，灰度阈值取值如下：Tm1＝0.2，Tm2＝0.5。 【当前权利人】北京航空航天大学 【当前专利权人地址】北京市海淀区学院路37号北京航空航天大学电子信息工程学院 【统一社会信用代码】12100000400011227Y 【引证次数】3.0 【被引证次数】2 【自引次数】1.0 【他引次数】2.0 【被他引次数】2.0 【家族引证次数】3.0 【家族被引证次数】9
