一种基于盲元检测与分析的红外图像质量评价方法与流程

本发明涉及图像质量评价的
技术领域：
，具体地说是一种基于盲元检测与分析的红外图像质量评价方法。
背景技术：
：红外图像是由红外探测器经信号采集和数据处理成像得到，受材料、制造工艺、外部环境等各种因素的限制，红外探测器在制造和使用过程中不可避免地存在盲元：盲元是指探测器响应过高或过低的成像单元，在红外图像中具体表现为无法正常成像，形成过亮、过暗、闪烁或无效的像素单元。盲元的存在严重影响了红外图像质量，甚至会妨碍系统后续目标的探测与识别。因此，迫切需要一种红外图像质量的客观评价方法，准确检测出图像中的盲元并分析其分布情况，判断红外图像是否有效，以满足目标识别系统的实际应用需求。传统的红外图像质量评价方法是在选择理想图像作为参考图像的情况下，比较待评价图像与参考图像之间的差异，峰值信噪比和均方误差是比较常见的两种质量评价方法。但这类方法的局限性在于，它们是从图像像素值的全局统计出发，对于图像局部的质量无从把握。而对于理想图像无法获取的情况，一般基于图像统计特性的方法，采用图像像素均值、标准差、平均梯度等指标对图像质量进行评价，难以准确反映盲元对红外图像质量的影响程度。技术实现要素：针对现有技术的不足，本发明提供一种基于盲元检测与分析的红外图像质量评价方法，解决采用图像像素均值、标准差、平均梯度等指标对图像质量进行评价，难以准确反映盲元对红外图像质量的影响程度的问题。本发明为实现上述目的所采用的技术方案是：一种基于盲元检测与分析的红外图像质量评价方法，包括以下步骤：步骤1：获取红外探测器原始图像iorg，检测出原始图像iorg中的盲元像素，转换成二值化图像ibin；步骤2：对二值化图像ibin进行盲元连通区域标识，并对整个图像进行扫描，得到所有盲元连通区域；步骤3：计算每个连通区域的面积，对其分类统计，得到连通区域的面积分布；步骤4：对二值化图像ibin进行网格划分，统计每个网格内的盲元数量，得到盲元数量的网格分布；步骤5：设定面积分布评价阈值和网格分布评价阈值，根据连通区域的面积分布和盲元数量的网格分布对原始图像iorg进行评价。所述检测出原始图像iorg中的盲元像素，转换成二值化图像ibin，包括：步骤1.1：对原始图像iorg进行5×5中值滤波，得到中值滤波图像imed；将原始图像iorg与中值滤波图像imed进行对减操作，得到过滤图像iflt，即iflt＝|iorg-imed|；步骤1.2：对原始图像iorg进行5×5均值滤波，得到均值滤波图像imean；步骤1.3：将原始图像iorg与均值滤波图像imean分别进行对加操作与对减操作，分别得到对加图像iadd和对减图像isub，即iadd＝iorg+imean，isub＝|iorg-imean|；步骤1.4：将对减图像isub与对加图像iadd进行对除操作，得到对比度图像icon，即icon＝isub./iadd，其中，若当分母为零时，结果为零；步骤1.5：设置对比度阈值t1，若对比度图像icon的值不小于阈值t1，则将对比度图像icon的值置为1，否则置为0；更新对比度图像icon；步骤1.6：将过滤图像iflt与对比度图像icon进行对乘操作，得到二值化图像ibin，即ibin＝iflt.×icon，若ibin的值大于0，则置为1，否则为0；更新二值化图像ibin。所述对二值化图像ibin进行盲元连通区域标识，包括：选取二值化图像ibin中某一盲元点，判断其周围像素点中是否有相邻盲元，如果是，则该盲元点与相邻盲元均为该盲元连通区域中的点；否则，该盲元的连通区域为该盲元点。盲元点周围像素点包括上、下、左、右、左上、右上、左下、右下对应的像素点。所述连通区域的面积为该连通区域的像素数量。所述根据连通区域的面积分布和盲元数量的网格分布对原始图像iorg进行评价，包括：步骤5.1：判断每个连通区域的面积是否不大于面积分布评价阈值，如果是，则执行步骤5.2，否则当前原始图像iorg无效；步骤5.2：判断每个网格中的盲元数量是否不大于网格分布的评价阈值，如果是，则当前原始图像iorg有效；否则当前原始图像iorg无效。本发明具有以下有益效果及优点：1.本发明采用图像处理的方法进行盲元检测，对图像的适应性强，可以有效检测不同型号探测器得到的红外图像；2.本发明通过分析盲元数量和网格化的局部盲元分布，能准确、客观地评价红外图像质量。附图说明图1为本发明的方法流程图；图2为探测器得到的原始红外图像；图3为自动检测出的盲元图像；图4(a)为8连通算法图；图4(b)为一个连通区域的示意图；图5为盲元面积分类的分布情况图；图6(a)为原始图像分成10×10的网格区域图；图6(b)为检测出的盲元图像分成10×10的网格区域图；图7为按照网格统计的局部盲元数量直方图。具体实施方式下面结合附图及实施例对本发明做进一步的详细说明。为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图对本发明的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明。但本发明能够以很多不同于在此描述的其他方式来实施，本领域技术人员可以在不违背发明内涵的情况下做类似改进，因此本发明不受下面公开的具体实施的限制。除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的
技术领域：
的技术人员通常理解的含义相同。本文中在发明的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本发明。如图1所示，本发明方法包括以下步骤：(1)采用图像处理的方法自动检测所获取的红外探测器图像中的盲元像素，并转换生成二值化图像，用1来标识盲元像素，用0来标识非盲元的图像像素；(2)对二值化图像进行盲元连通区域标识，并对整个图像进行扫描，得到所有盲元连通区域；(3)统计每个连通区域的面积，得到图像总体的盲元数量分布情况；(4)将图像网格化，统计每个网格内的盲元数量，得到图像各个局部的盲元数量分布情况；(5)按照总体与局部盲元数量的分类及设定的阈值定量评价盲元对图像质量的影响。自动检测红外图像中盲元像素的图像处理的具体方法为：探测器得到的原始红外图像主要包含的信息有：均匀背景，盲元像素。盲元像素在图像中表现为类似于椒盐噪声的过亮或过暗的像素点，不规则地散落在图像中，如图2所示。对于本发明而言，盲元是待检测目标，相对于背景具有较强的对比度，利用背景抑制进行预处理，同时满足图像中较高对比度的位置来确定盲元。对比度的计算公式如下：其中，f(x,y)代表当前像素，b(x,y)代表当前点所在位置的背景估计值，获取方法为使用5×5大小的中值滤波：对原始图像iorg进行中值滤波，滤波器大小为5×5，得到中值滤波图像imed作为背景估计，将原始图像iorg与中值滤波图像imed进行像素对减操作，得到背景抑制图像iflt，对减操作过程取绝对值以保证结果为非负值，即：iflt＝|iorg-imed|对原始图像iorg进行均值滤波，滤波器大小为5×5，得到均值滤波图像imean；将原始图像iorg与均值滤波图像imean分别进行像素对加操作与像素对减操作，分别得到对加图像iadd和对减图像isub，对减操作过程取绝对值以保证结果为非负值，即：iadd＝iorg+imeanisub＝|iorg-imean|将得到的对减图像isub与对加图像iadd进行像素对除操作，得到对比度图像icon，即：icon＝isub./iadd其中，不同于一般矩阵乘/除法的是，符号的./含义为两个图像矩阵相同行列位置的像素值与像素值进行相除，结果存储在输出图像矩阵的相对应位置，后文的.×符号同理。上式中的特殊情形为，当原始图像iorg和均值滤波图像imean在某同一位置p的像素值均为零时，会出现对加图像相对应位置p的像素值也为零的情况，此时icon＝isub./iadd的分母iadd(p)的值为零，无法进行相除，则不进行除法操作，直接设定icon(p)的值为零。因此，对比度图像icon在点p的像素值的表达式可以写成：设置一定的对比度阈值t1，若对比度图像icon的对比度强度值不小于阈值t1，则认为是盲元点，将对比度图像icon该点的值置为1，否则置为0；用上述策略更新对比度图像icon，表达式可以写成：在本实施例中，阈值t1的取值为0.1。设置一定的过滤阈值t2，若过滤图像iflt的值小于阈值t2，则将过滤图像iflt的值置为0，否则值不变；用上述策略更新对比度图像iflt，表达式可以写成：在本实施例中，阈值t2的取值为3。将上文得到的过滤图像iflt与对比度图像icon进行像素对乘操作，得到二值化图像ibin，即：ibin＝iflt.×icon等价于：ibin(p)＝iflt(p)×icon(p)若ibin的值大于0，则置为1，否则为0；以上述策略更新二值化图像ibin，表达式可以写为：则最终得到的二值化图像ibin如图3所示。接下来对二值化图像进行盲元连通区域标识。对于已检测出盲元的二值化图像，如图4(a)所示，这里选取某一盲元点p周围8个像素(上、下、左、右、左上、右上、左下、右下)，若有相邻盲元，则该盲元点与相邻盲元都被认为是该盲元连通区域中的点。如图4(b)所示，可以看到标记为1的蓝色区域即为该绿色区域中的连通区域。对整个图像进行扫描，得到所有盲元连通区域。以连通区域的像素数量作为该连通区域的面积值；遍历所有盲元连通区域，将连通区域面积值按照数值大小分为8个类别，其面积的范围分别为：1，2，3～4，5～7，8～10，11～20，21～50，51以上。这样即得到了盲元数量在图像整体上的分布情况，从8个面积范围的数量可以反映出盲元在图像中的整体数量和各个盲元连通区域的大小情况。对于该图像而言，其面积分布情况具体列于表1：表1盲元面积分布情况面积范围区域数量1262213～4375～7408～101811～201121～50351及以上0如图5所示，可以看出盲元面积范围的分布情况。最后对图像进行网格化，如图6(a)和(b)所示，将整个图像分成10×10的网格区域；分别统计每个网格内包含的盲元数量，可以得到盲元在整个图像区域的散落情况。对于该图像而言，10×10的网格区域每个网格内的盲元数量列于表2,。将网格化图像的盲元数量分布做直方图，如图7所示，可以直观看到盲元在图像中的散落分布情况。表210×10网格图像中每个网格的盲元数量2466000370062010618066964631127081660107101026011111102162130218006315310961430661641003393871017111300755794063954270934411412008针对于不同应用需求，本发明中8类连通区域面积值和网格大小均可以按需调整。(1)对于连通区域面积值范围而言，每个类别的范围值都可以进行更改，例如对于成像质量较好、盲元较少的图像，可以在1～20的范围内分为8类，而对于成像质量较差，盲元较多的图像，则可以在1～100的范围内分为8类，这样即可对不同成像质量的探测器具有较强的适应能力；(2)对于10×10的网格大小而言，亦可根据对图像关注的细致程度，将网格大小粗化到5×5，或细化到20×20，这样即可重点关注图像感兴趣区域内的盲元数量。综上，通过对图像中盲元的检测与对盲元数量的从整体到局部的分析，从盲元的角度定量、客观地给出图像的质量评价：(1)从面积分布可以看到，对于本例图像，图像中盲元区域面积在5～7范围内的数量最多，面积在3～4范围内的数量次之，面积为1也就是盲元单点的数量亦较为可观；但图像中面积大于20的盲元区域数量较少，且没有面积大于50的盲元区域，也就是说，图像中没有大面积连成一片的盲元区块；(2)从网格分布可以看到，盲元在图像的分布相对比较平均，数量最多的区域主要出现在图像上下两个底边，而数量相对较多的区域主要分布在图像中心区域；虽然整体上没有出现盲元区域面积特别大的情况，但由于盲元分布相对比较分散，即图像各个位置都有盲元存在，对于成像质量要求较高的使用场景会存在干扰，因此该图像质量不合格。当前第1页12