CN108960254B - 一种基于mscn信息识别模糊图像的方法和装置、计算设备及存储介质 - Google Patents
一种基于mscn信息识别模糊图像的方法和装置、计算设备及存储介质 Download PDFInfo
- Publication number
- CN108960254B CN108960254B CN201810678755.5A CN201810678755A CN108960254B CN 108960254 B CN108960254 B CN 108960254B CN 201810678755 A CN201810678755 A CN 201810678755A CN 108960254 B CN108960254 B CN 108960254B
- Authority
- CN
- China
- Prior art keywords
- image
- mean
- value
- identified
- impairment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/50—Extraction of image or video features by performing operations within image blocks; by using histograms, e.g. histogram of oriented gradients [HoG]; by summing image-intensity values; Projection analysis
- G06V10/507—Summing image-intensity values; Histogram projection analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/74—Image or video pattern matching; Proximity measures in feature spaces
- G06V10/75—Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
- G06V10/755—Deformable models or variational models, e.g. snakes or active contours
- G06V10/7553—Deformable models or variational models, e.g. snakes or active contours based on shape, e.g. active shape models [ASM]
Abstract
本发明提供了一种能够基于MSCN信息快速地识别模糊图像的方法和装置、计算设备及存储介质。所述方法包括以下步骤:对待识别图像进行预处理以得到MSCN图像;对所述MSCN图像的统计直方图中的‑0.5至0.5范围内的MSCN系数进行统计,以得到MSCN系数和;将所述MSCN系数和与预定值进行比较;如果所述MSCN系数和大于等于所述预定值,则判别所述待识别图像为模糊图像,否则判别所述待识别图像为清晰图像。本发明的方法和装置不仅识别速度快,而且准确率高,从而可以有效应用到手机终端。
Description
技术领域
本发明涉及图像识别技术领域,特别涉及一种基于MSCN信息识别模糊图像的方法和装置、计算设备及存储介质。
背景技术
图像质量客观评价是指通过设计合理的算法使计算机自动精确地预测图像的感知质量。图像质量的评价结果可以为图像处理算法的参数优化、图像处理系统的性能评估和图像处理设备的质量检测提供重要的指标和依据,已成为图像处理领域的研究热点之一。其中,盲图像质量评价(Blind Image Quality Assessment,BIQA)是指在没有参考图像的情况下对任意输入图像的视觉质量进行精确预测。由于在绝大多数的实际应用场景中,测试图像对应的参考图像是无法或难以得到的,因此对于盲图像质量评价的研究至关重要。
通过对现有文献的检索,目前无参考图像质量评价的代表性方法有两种。一种是AnishMittal等人在2012年IEEE Transactions on Image Processing,vol.21(12),pp.4695-4708(2012年IEEE图像处理会刊第21卷12期,4695至4708页)上发表的“No-reference imagequality assessment in the spatial domain(空间域无参考图像质量评价)”一文中提出的用于无参考图像质量评价的自然图像统计特征构建方法(简称为BRISQUE)。该方法直接对多个不同尺度或方位的预处理图像拟合广义高斯分布模型,将模型参数共计36个系数值作为图像自然统计特征,采用支持向量机(Support VectorMachine,简称SVM)进行模型训练和测试。然而拟合广义高斯分布模型的步骤对输入图像做了过强的假设,不可避免地降低了图像的原始信息量,进而影响了模型的精度。而且,该方法速度慢、需要用到svm分类器导致模型大。
另一种代表性方法是Wufeng Xue等人在2014年IEEE Transactions on ImageProcessing,vol.23(11),pp.4850-4862(2014年IEEE图像处理会刊第23卷11期,4850至4862页)上发表的“Blind Image Quality Assessment Using Joint Statistics ofGradient Magnitude and Laplacian Features(使用梯度幅度和拉普拉斯特征的联合统计量进行盲图像质量评估)”一文中提出的联合梯度和拉普拉斯-高斯信号的盲图像质量评价方法。该方法从人类视觉的处理机制出发,首先提取图像的梯度模值和LOG信号,通过对这两种信号进行联合自适应归一化消除高阶冗余性,并且,在此基础上分别计算二者的边缘分布概率直方图和条件分布概率直方图作为图像感知特征。该方法符合人类视觉早期处理过程,并且具有更好的区分自然图像与非自然图像的能力。然而,该方法同样地速度慢,且需要用到svm分类器导致模型大。
发明内容
为了克服现有技术中的一个或多个缺陷,本发明提供一种能够基于均值减损对比归一化(mean subtracted contrast normalized,MSCN)信息快速地识别模糊图像的方法和装置、计算设备及存储介质,以便于在人脸检测前剔除模糊图像,提高人脸识别率。
根据本发明的第一方面,提供一种基于均值减损对比归一化信息来识别模糊图像的方法。所述方法包括以下步骤:对待识别图像进行预处理以得到均值减损对比归一化图像;对所述均值减损对比归一化图像的统计直方图中的-0.5至0.5范围内的均值减损对比归一化系数进行统计,以得到均值减损对比归一化系数和;将所述均值减损对比归一化系数和与预定值进行比较;如果所述均值减损对比归一化系数和大于等于所述预定值,则判别所述待识别图像为模糊图像,否则判别所述待识别图像为清晰图像。
优选地,所述预处理采用如下公式:
其中,I(i,j)是所述待识别图像;i,j是像素坐标,i∈1,2,...M,j∈1,2,...N,M和N分别为所述待识别图像的长宽;C=1,为防止分母为零的常数;I'(i,j)为所述均值减损对比归一化图像;μ(i,j)和σ(i,j)通过下式来得到:
其中,w={wk,l|k=-K,...K,l=-L,...L}是二维高斯窗口,K=L=3,μ(i,j)是在窗口内的局部均值,σ(i,j)是窗口内的局部方差。
优选地,所述预定值为0.9。
根据本发明的第二方面,提供一种基于梯度幅值和方向值稀疏特征信息以及均值减损对比归一化信息来识别模糊图像的方法。所述方法包括以下步骤:计算待识别图像的梯度幅值;将梯度幅值归一化为0~1;以第一间隔gap0对归一化后的梯度幅值进行量化,即,量化后的梯度幅度范围是0~1/gap0;统计量化后的梯度幅度的直方图中的0~rate*1/gap0范围内的非空点数N0,其中,rate为比率常数;计算所述待识别图像在0~360度范围内的梯度方向值;将梯度方向值归一化为0~1;以第二间隔gap1对归一化后的梯度方向值进行量化,即,量化后的梯度方向值范围是0~1/gap1;统计量化后的梯度方向值的直方图中的所有的非空点数N1;对所述待识别图像进行预处理以得到均值减损对比归一化图像;对所述均值减损对比归一化图像的统计直方图中的-0.5至0.5范围内的均值减损对比归一化系数进行统计,以得到均值减损对比归一化系数和N2;基于非空点数N0和N1以及均值减损对比归一化系数和N2来判别所述待识别图像为模糊图像还是为清晰图像。
优选地,所述基于非空点数N0和N1以及均值减损对比归一化系数和N2来判别所述待识别图像为模糊图像还是为清晰图像的步骤包括:
计算下式:rate0*|N0-center0|/(|N0-center0|+|N0-center0’|)+rate1*|N1-center1|/(|N1-center1|+|N1-center1’|)+(1-rate0-rate1)*|N2-center2|/(|N2-center2|+|N2-center2’|),其中,rate0和rate1是交叉验证参数,center0和center1分别为对训练集中的所有模糊样本计算而得到的非空点数N0和N1的均值,center0’和center1’分别为对训练集中的所有清晰样本计算而得到的非空点数N0和N1的均值,center2为对训练集中的所有模糊样本计算而得到的均值减损对比归一化系数和N2的均值,center2’为对训练集中的所有清晰样本计算而得到的均值减损对比归一化系数和N2的均值;
将计算出的上式与预定值进行比较;
如果计算出的上式小于所述预定值,则判别所述待识别图像为模糊图像,否则判别所述待识别图像为清晰图像。
优选地,所述预处理采用如下公式:
其中,I(i,j)是所述待识别图像;i,j是像素坐标,i∈1,2,...M,j∈1,2,...N,M和N分别为所述待识别图像的长宽;C=1,为防止分母为零的常数;I'(i,j)为所述均值减损对比归一化图像;μ(i,j)和σ(i,j)通过下式来得到:
其中,w={wk,l|k=-K,...K,l=-L,...L}是二维高斯窗口,K=L=3,μ(i,j)是在窗口内的局部均值,σ(i,j)是窗口内的局部方差。
优选地,所述预定值为0.5。
优选地,所述归一化为最大最小归一化。
优选地,所述第一间隔gap0为0.001,所述比率常数rate为0.01,并且,所述第二间隔gap1为0.00025。
根据本发明的第三方面,提供一种用于基于均值减损对比归一化信息来识别模糊图像的装置。所述装置包括:预处理模块,用于对待识别图像进行预处理以得到均值减损对比归一化图像;均值减损对比归一化系数统计模块,用于对所述均值减损对比归一化图像的统计直方图中的-0.5至0.5范围内的均值减损对比归一化系数进行统计,以得到均值减损对比归一化系数和;比较模块,用于将所述均值减损对比归一化系数和与预定值进行比较;判别模块,用于如果所述均值减损对比归一化系数和大于等于所述预定值,则判别所述待识别图像为模糊图像,否则判别所述待识别图像为清晰图像。
优选地,所述预处理采用如下公式:
其中,I(i,j)是所述待识别图像;i,j是像素坐标,i∈1,2,...M,j∈1,2,...N,M和N分别为所述待识别图像的长宽;C=1,为防止分母为零的常数;I'(i,j)为所述均值减损对比归一化图像;μ(i,j)和σ(i,j)通过下式来得到:
其中,w={wk,l|k=-K,...K,l=-L,...L}是二维高斯窗口,K=L=3,μ(i,j)是在窗口内的局部均值,σ(i,j)是窗口内的局部方差。
优选地,所述预定值为0.9。
根据本发明的第四方面,提供一种用于基于梯度幅值和方向值稀疏特征信息以及均值减损对比归一化信息来识别模糊图像的装置。所述装置包括:梯度幅值计算模块,用于计算待识别图像的梯度幅值;梯度幅值归一化模块,用于将梯度幅值归一化为0~1;梯度幅值量化模块,用于以第一间隔gap0对归一化后的梯度幅值进行量化,即,量化后的梯度幅度范围是0~1/gap0;梯度幅值统计模块,用于统计量化后的梯度幅度的直方图中的0~rate*1/gap0范围内的非空点数N0,其中,rate为比率常数;梯度方向值计算模块,用于计算所述待识别图像在0~360度范围内的梯度方向值;梯度方向值归一化模块,用于将梯度方向值归一化为0~1;梯度方向值量化模块,用于以第二间隔gap1对归一化后的梯度方向值进行量化,即,量化后的梯度方向值范围是0~1/gap1;梯度方向值统计模块,用于统计量化后的梯度方向值的直方图中的所有的非空点数N1;预处理模块,用于对所述待识别图像进行预处理以得到均值减损对比归一化图像;均值减损对比归一化系数统计模块,用于对所述均值减损对比归一化图像的统计直方图中的-0.5至0.5范围内的均值减损对比归一化系数进行统计,以得到均值减损对比归一化系数和N2;判别模块,用于基于非空点数N0和N1以及均值减损对比归一化系数和N2来判别所述待识别图像为模糊图像还是为清晰图像。
优选地,所述判别模块包括:判别式计算单元,用于计算下式:
rate0*|N0-center0|/(|N0-center0|+|N0-center0’|)+rate1*|N1-center1|/(|N1-center1|+|N1-center1’|)+(1-rate0-rate1)*|N2-center2|/(|N2-center2|+|N2-center2’|),其中,rate0和rate1是交叉验证参数,center0和center1分别为对训练集中的所有模糊样本计算而得到的非空点数N0和N1的均值,center0’和center1’分别为对训练集中的所有清晰样本计算而得到的非空点数N0和N1的均值,center2为对训练集中的所有模糊样本计算而得到的均值减损对比归一化系数和N2的均值,center2’为对训练集中的所有清晰样本计算而得到的均值减损对比归一化系数和N2的均值;比较单元,用于将判别式计算单元计算出的上式与预定值进行比较;综合判别单元,用于如果判别式计算单元计算出的上式小于所述预定值,则判别所述待识别图像为模糊图像,否则判别所述待识别图像为清晰图像。
优选地,所述预处理采用如下公式:
其中,I(i,j)是所述待识别图像;i,j是像素坐标,i∈1,2,...M,j∈1,2,...N,M和N分别为所述待识别图像的长宽;C=1,为防止分母为零的常数;I'(i,j)为所述均值减损对比归一化图像;μ(i,j)和σ(i,j)通过下式来得到:
其中,w={wk,l|k=-K,...K,l=-L,...L}是二维高斯窗口,K=L=3,μ(i,j)是在窗口内的局部均值,σ(i,j)是窗口内的局部方差。
优选地,所述预定值为0.5。
优选地,所述归一化为最大最小归一化。
优选地,所述第一间隔gap0为0.001,所述比率常数rate为0.01,并且,所述第二间隔gap1为0.00025。
根据本发明的第五方面,提供一种计算设备。所述计算设备包括:处理器;以及存储器,其上存储有可执行代码,当所述可执行代码被所述处理器执行时,使所述处理器执行上述方法。
根据本发明的第六方面,提供一种非暂时性机器可读存储介质,其上存储有可执行代码,当所述可执行代码被电子设备的处理器执行时,使所述处理器执行上述方法。
与现有技术相比,本发明具有以下有益效果:本发明提供一种新的基于MSCN信息来识别模糊图像的方法和装置。所述方法和装置只需要统计MSCN图像的直方图中的-0.5至0.5范围内的MSCN系数,并将得到的MSCN系数之和与预定值进行比较,即可判别待识别图像是模糊图像还是清晰图像。由此可知,所述方法和装置不仅识别速度快,而且准确率高,从而可以有效应用到手机终端。
本领域技术人员应当理解的是,能够用本发明实现的目的和优点不限于以上具体所述,并且根据以下详细说明将更清楚地理解本发明能够实现的上述和其他目的。
并且,应当理解,前述大体的描述和后续详尽的描述均为示例性说明和解释,并不应当用作对本发明所要求保护内容的限制。
附图说明
参考随附的附图,本发明更多的目的、功能和优点将通过本发明实施方式的如下描述得以阐明,其中:
图1为根据本发明的第一示例性实施例的模糊图像识别方法的流程图;
图2为根据本发明的第一示例性实施例的模糊图像识别装置的框图;
图3为根据本发明的第二示例性实施例的模糊图像识别方法的流程图;
图4为根据本发明的第二示例性实施例的模糊图像识别方法中的判别步骤的流程图;
图5为根据本发明的第二示例性实施例的模糊图像识别方法的处理的示意图;
图6为根据本发明的第二示例性实施例的模糊图像识别装置的框图;
图7为根据本发明的第二示例性实施例的模糊图像识别装置中的判别模块的框图;
图8为根据本发明示例性实施例的可用于实现模糊图像识别方法的计算设备的结构示意图。
具体实施方式
通过参考示范性实施例,本发明的目的和功能以及用于实现这些目的和功能的方法将得以阐明。然而,本发明并不受限于以下所公开的示范性实施例;可以通过不同形式来对其加以实现。说明书的实质仅仅是帮助相关领域技术人员综合理解本发明的具体细节。
为了使本发明的技术方案更加清楚、明白,下面将参照附图并结合具体实施例对本发明的模糊图像识别方法和装置进行详细描述。
图1示出了根据本发明的第一示例性实施例的模糊图像识别方法的流程图。如图1所示,所述模糊图像识别方法包括预处理步骤S1、MSCN系数统计步骤S2、比较步骤S3和判别步骤S4。
首先,在步骤S1中,采用如下公式对待识别图像进行预处理,以得到MSCN图像:
其中,I(i,j)是待识别图像;i,j是像素坐标,i∈1,2,...M,j∈1,2,...N,M和N分别为待识别图像的长宽;C=1,为防止分母为零的常数;I'(i,j)为MSCN图像;μ(i,j)和σ(i,j)通过下式来得到:
其中,w={wk,l|k=-K,...K,l=-L,...L}是二维高斯窗口,K=L=3,μ(i,j)是在窗口内的局部均值,σ(i,j)是窗口内的局部方差。
接着,在步骤S2中,利用下式对MSCN图像的统计直方图中的-0.5至0.5范围内的MSCN系数进行统计,以得到MSCN系数和N2:
其中,mscn表示MSCN图像的统计直方图中的系数值。
接下来,在步骤S3中,将步骤S2中获得的MSCN系数和N2与预定值T进行比较。经验上,通常可以将预定值T设为0.9。然后,在步骤S4中,基于步骤S3的比较结果来判别待识别图像为模糊图像还是清晰图像。具体地,如果MSCN系数和大于等于预定值T,则判别待识别图像为模糊图像。相反地,如果MSCN系数和小于预定值T,则判别待识别图像为清晰图像。
发明人采用上述模糊图像识别方法进行了试验。试验结果是:可达90.1%的识别率。由此可知,本示例性实施例的模糊图像识别方法不仅运算速度快,而且准确率高。
另外,本示例性实施例还提供一种用于实现上述模糊图像识别方法的模糊图像识别装置。图2示出了该模糊图像识别装置的框图。如图2所示,模糊图像识别装置100包括预处理模块101,MSCN系统统计模块102,比较模块103和判别模块104。
预处理模块101用于通过如下公式对待识别图像进行预处理以得到MSCN图像:
其中,I(i,j)是待识别图像;i,j是像素坐标,i∈1,2,...M,j∈1,2,...N,M和N分别为待识别图像的长宽;C=1,为防止分母为零的常数;I'(i,j)为MSCN图像;μ(i,j)和σ(i,j)通过下式来得到:
其中,w={wk,l|k=-K,...K,l=-L,...L}是二维高斯窗口,K=L=3,μ(i,j)是在窗口内的局部均值,σ(i,j)是窗口内的局部方差。
MSCN系统统计模块102用于通过下式对MSCN图像的统计直方图中的-0.5至0.5范围内的MSCN系数进行统计,以得到MSCN系数和:
其中,mscn表示MSCN图像的统计直方图中的系数值。
比较模块103用于将MSCN系统统计模块102所得到的MSCN系数和N2与预定值T进行比较。经验上,通常可以将预定值T设为0.9。判别模块104用于基于比较模块103所得到的比较结果来判别待识别图像为模糊图像还是清晰图像。具体地,如果MSCN系数和大于等于预定值T,则判别待识别图像为模糊图像。相反地,如果MSCN系数和小于预定值T,则判别待识别图像为清晰图像。
图3示出了根据本发明的第二示例性实施例的模糊图像识别方法的流程图,并且,图5给出了该模糊图像识别方法的处理的示意图。下面参照图3和图5来描述根据本发明的第二示例性实施例的模糊图像识别方法。
如图3所示,所述模糊图像识别方法包括梯度幅值计算步骤S101、梯度幅值归一化步骤S102、梯度幅值量化步骤S103、梯度幅值统计步骤S104、梯度方向值计算步骤S105、梯度方向值归一化步骤S106、梯度方向值量化步骤S107、梯度方向值统计步骤S108、预处理步骤S109、MSCN系数统计步骤S110和判别步骤S111。
首先,在步骤S101中,计算待识别图像的梯度幅值。具体地,利用下式获得待识别图像中每个像素点的梯度幅值:
其中,A(x,y)是待识别图像上的坐标为(x,y)的像素点的梯度幅值,I(x,y)是待识别图像上的坐标为(x,y)的像素点的像素值。
接下来,在步骤S102中,将步骤S101中计算得到的梯度幅值归一化为0~1。具体地,利用下式对梯度幅值进行最大最小归一化:
其中,A*(x,y)为归一化后的[0,1]范围内的梯度幅值,Amax和Amin分别是步骤S101中计算得到的梯度幅值中的最大值和最小值。
接着,在步骤S103中,以间隔gap0对归一化后的梯度幅值进行量化,即,量化后的梯度幅度范围是0~1/gap0。例如,假设取间隔gap0为0.001,则量化后的梯度幅度范围是0~1000。然后,在步骤S104中,统计量化后的梯度幅度的直方图中的0~rate*1/gap0范围内的非空点数N0,其中,rate为比率常数。例如,假设rate为0.01,即,提取量化后的梯度幅度的直方图中的前1%梯度幅值并统计其中的非空点数。
在步骤S105中,计算待识别图像在0~360度范围内的梯度方向值。具体地,利用下式获得待识别图像中每个像素点在各个方向上的梯度方向值:
其中,T(x,y)是待识别图像上的坐标为(x,y)的像素点的梯度方向值,I(x,y)是待识别图像上的坐标为(x,y)的像素点的像素值。
接下来,在步骤S106中,将步骤S105中计算得到的梯度方向值归一化为0~1。具体地,利用下式对梯度方向值进行最大最小归一化:
其中,T*(x,y)为归一化后的[0,1]范围内的梯度方向值,Tmax和Tmin分别是步骤S105中计算得到的梯度方向值中的最大值和最小值。
接着,在步骤S107中,以间隔gap1对归一化后的梯度方向值进行量化,即,量化后的梯度方向值范围是0~1/gap1。例如,假设取间隔gap1为0.00025,则量化后的梯度幅度范围是0~4000。然后,在步骤S108中,统计量化后的梯度方向值的直方图中的所有的非空点数N1。
此后,在步骤S109中,采用如下公式对待识别图像进行预处理,以得到MSCN图像:
其中,I(i,j)是待识别图像;i,j是像素坐标,i∈1,2,...M,j∈1,2,...N,M和N分别为待识别图像的长宽;C=1,为防止分母为零的常数;I'(i,j)为MSCN图像;μ(i,j)和σ(i,j)通过下式来得到:
其中,w={wk,l|k=-K,...K,l=-L,...L}是二维高斯窗口,K=L=3,μ(i,j)是在窗口内的局部均值,σ(i,j)是窗口内的局部方差。
接着,在步骤S110中,利用下式对MSCN图像的统计直方图中的-0.5至0.5范围内的MSCN系数进行统计,以得到MSCN系数和N2:
其中,mscn表示MSCN图像的统计直方图中的系数值。
最后,在步骤S111中,基于步骤S104中得到的非空点数N0、步骤S108中得到的非空点数N1以及步骤S110中得到的MSCN系数和N2来判别待识别图像为模糊图像还是为清晰图像。具体地,首先,在判别式计算步骤S111A中,计算下式:
rate0*|N0-center0|/(|N0-center0|+|N0-center0’|)+rate1*|N1-center1|/(|N1-center1|+|N1-center1’|)+(1-rate0-rate1)*|N2-center2|/(|N2-center2|+|N2-center2’|),
其中,rate0和rate1是交叉验证参数,center0和center1分别为对训练集中的所有模糊样本计算而得到的非空点数N0和N1的均值,center0’和center1’分别为对训练集中的所有清晰样本计算而得到的非空点数N0和N1的均值,center2为对训练集中的所有模糊样本计算而得到的均值减损对比归一化系数和N2的均值,center2’为对训练集中的所有清晰样本计算而得到的均值减损对比归一化系数和N2的均值。接着,在比较步骤S111B中,将步骤S111A中计算出的上式与预定值进行比较。经验上,通常可以将预定值设为0.5。然后,在综合判别步骤S111C中,基于步骤S111B的比较结果来判别待识别图像为模糊图像还是清晰图像。具体地,如果步骤S111A中计算出的上式小于预定值,则判别待识别图像为模糊图像。相反地,如果步骤S111A中计算出的上式大于等于预定值,则判别待识别图像为清晰图像。
发明人采用本示例性实施例的模糊图像识别方法进行了试验。试验结果是:可高达93.0%的识别率。相比于第一示例性实施例的模糊图像识别方法,本示例性实施例的模糊图像识别方法还考虑到了梯度幅值和方向值稀疏特征信息,因而保留了原始图像的更多信息,从而能够获得更好的评价效果。
另外,本示例性实施例还提供一种用于实现上述模糊图像识别方法的模糊图像识别装置。图6示出了该模糊图像识别装置的框图。如图6所示,所述模糊图像识别装置200包括梯度幅值计算模块201、梯度幅值归一化模块202、梯度幅值量化模块203、梯度幅值统计模块204、梯度方向值计算模块205、梯度方向值归一化模块206、梯度方向值量化模块207、梯度方向值统计模块208、预处理模块209、MSCN系数统计模块210和判别模块211。
梯度幅值计算模块201用于计算待识别图像的梯度幅值。具体地,利用下式获得待识别图像中每个像素点的梯度幅值:
其中,A(x,y)是待识别图像上的坐标为(x,y)的像素点的梯度幅值,I(x,y)是待识别图像上的坐标为(x,y)的像素点的像素值。
梯度幅值归一化模块202用于将梯度幅值计算模块201所得到的梯度幅值归一化为0~1。具体地,利用下式对梯度幅值进行最大最小归一化:
其中,A*(x,y)为归一化后的[0,1]范围内的梯度幅值,Amax和Amin分别是梯度幅值计算模块201计算出的梯度幅值中的最大值和最小值。
梯度幅值量化模块203用于以间隔gap0对由梯度幅值归一化模块202归一化后的梯度幅值进行量化,即,量化后的梯度幅度范围是0~1/gap0。例如,假设取间隔gap0为0.001,则量化后的梯度幅度范围是0~1000。
梯度幅值统计模块204用于统计梯度幅值量化模块203量化后的梯度幅度的直方图中的0~rate*1/gap0范围内的非空点数N0,其中,rate为比率常数。例如,假设rate为0.01,即,提取量化后的梯度幅度的直方图中的前1%梯度幅值并统计其中的非空点数。
梯度方向值计算模块205用于计算待识别图像在0~360度范围内的梯度方向值。具体地,利用下式获得待识别图像中每个像素点在各个方向上的梯度方向值:
其中,T(x,y)是待识别图像上的坐标为(x,y)的像素点的梯度方向值,I(x,y)是待识别图像上的坐标为(x,y)的像素点的像素值。
梯度方向值归一化模块206用于将梯度方向值计算模块205所得到的梯度方向值归一化为0~1。具体地,利用下式对梯度方向值进行最大最小归一化:
其中,T*(x,y)为归一化后的[0,1]范围内的梯度方向值,Tmax和Tmin分别是梯度方向值计算模块205计算出的梯度方向值中的最大值和最小值。
梯度方向值量化模块207用于以间隔gap1对梯度方向值归一化模块206归一化后的梯度方向值进行量化,即,量化后的梯度方向值范围是0~1/gap1。例如,假设取间隔gap1为0.00025,则量化后的梯度幅度范围是0~4000。
梯度方向值统计模块208用于统计梯度方向值量化模块207量化后的梯度方向值的直方图中的所有的非空点数N1。
预处理模块209用于采用如下公式对待识别图像进行预处理,以得到MSCN图像:
其中,I(i,j)是待识别图像;i,j是像素坐标,i∈1,2,...M,j∈1,2,...N,M和N分别为待识别图像的长宽;C=1,为防止分母为零的常数;I'(i,j)为MSCN图像;μ(i,j)和σ(i,j)通过下式来得到:
其中,w={wk,l|k=-K,...K,l=-L,...L}是二维高斯窗口,K=L=3,μ(i,j)是在窗口内的局部均值,σ(i,j)是窗口内的局部方差。
MSCN系数统计模块210用于通过下式对MSCN图像的统计直方图中的-0.5至0.5范围内的MSCN系数进行统计,以得到MSCN系数和N2:
其中,mscn表示MSCN图像的统计直方图中的系数值。
判别模块211用于基于梯度幅值统计模块204统计得到的非空点数N0、梯度方向值统计模块208统计得到的非空点数N1以及MSCN系数统计模块210统计得到的MSCN系数和N2来判别待识别图像为模糊图像还是为清晰图像。如图7所示,判别模块211包括判别式计算单元211A、比较单元211B和综合判别单元211C。
判别式计算单元211A用于计算下式:
rate0*|N0-center0|/(|N0-center0|+|N0-center0’|)+rate1*|N1-center1|/(|N1-center1|+|N1-center1’|)+(1-rate0-rate1)*|N2-center2|/(|N2-center2|+|N2-center2’|),其中,rate0和rate1是交叉验证参数,center0和center1分别为对训练集中的所有模糊样本计算而得到的非空点数N0和N1的均值,center0’和center1’分别为对训练集中的所有清晰样本计算而得到的非空点数N0和N1的均值,center2为对训练集中的所有模糊样本计算而得到的均值减损对比归一化系数和N2的均值,center2’为对训练集中的所有清晰样本计算而得到的均值减损对比归一化系数和N2的均值。
比较单元211B用于将判别式计算单元211A计算出的上式与预定值进行比较。经验上,通常可以将预定值设为0.5。
综合判别单元211C用于基于比较单元211B的比较结果来判别待识别图像为模糊图像还是清晰图像。具体地,如果判别式计算单元211A计算出的上式小于预定值,则判别待识别图像为模糊图像。相反地,如果判别式计算单元211A计算出的上式大于等于预定值,则判别待识别图像为清晰图像。
图8示出了根据本发明示例性实施例的可用于实现上述模糊图像识别方法的计算设备的结构示意图。
参见图8,计算设备1000包括存储器1010和处理器1020。
处理器1020可以是一个多核的处理器,也可以包含多个处理器。在一些实施例中,处理器1020可以包含一个通用的主处理器以及一个或多个特殊的协处理器,例如图形处理器(GPU)、数字信号处理器(DSP)等等。在一些实施例中,处理器1020可以使用定制的电路实现,例如特定用途集成电路(ASIC,Application Specific Integrated Circuit)或者现场可编程逻辑门阵列(FPGA,Field Programmable GateArrays)。
存储器1010可以包括各种类型的存储单元,例如系统内存、只读存储器(ROM),和永久存储装置。其中,ROM可以存储处理器1020或者计算机的其他模块需要的静态数据或者指令。永久存储装置可以是可读写的存储装置。永久存储装置可以是即使计算机断电后也不会失去存储的指令和数据的非易失性存储设备。在一些实施方式中,永久性存储装置采用大容量存储装置(例如磁或光盘、闪存)作为永久存储装置。另外一些实施方式中,永久性存储装置可以是可移除的存储设备(例如软盘、光驱)。系统内存可以是可读写存储设备或者易失性可读写存储设备,例如动态随机访问内存。系统内存可以存储一些或者所有处理器在运行时需要的指令和数据。此外,存储器1010可以包括任意计算机可读存储媒介的组合,包括各种类型的半导体存储芯片(DRAM,SRAM,SDRAM,闪存,可编程只读存储器),磁盘和/或光盘也可以采用。在一些实施方式中,存储器1010可以包括可读和/或写的可移除的存储设备,例如激光唱片(CD)、只读数字多功能光盘(例如DVD-ROM,双层DVD-ROM)、只读蓝光光盘、超密度光盘、闪存卡(例如SD卡、min SD卡、Micro-SD卡等等)、磁性软盘等等。计算机可读存储媒介不包含载波和通过无线或有线传输的瞬间电子信号。
存储器1010上存储有可处理代码,当可处理代码被处理器1020处理时,可以使处理器1020执行上文述及的基于机器学习的色情应用识别方法。
上文中已经参考附图详细描述了根据本发明的基于机器学习的色情应用识别方法和装置。
此外,根据本发明的方法还可以实现为一种计算机程序或计算机程序产品,该计算机程序或计算机程序产品包括用于执行本发明的上述方法中限定的上述各步骤的计算机程序代码指令。
或者,本发明还可以实施为一种非暂时性机器可读存储介质(或计算机可读存储介质、或机器可读存储介质),其上存储有可执行代码(或计算机程序、或计算机指令代码),当所述可执行代码(或计算机程序、或计算机指令代码)被电子设备(或计算设备、服务器等)的处理器执行时,使所述处理器执行根据本发明的上述方法的各个步骤。
本领域技术人员还将明白的是,结合这里的公开所描述的各种示例性逻辑块、模块、电路和算法步骤可以被实现为电子硬件、计算机软件或两者的组合。
附图中的流程图和框图显示了根据本发明的多个实施例的系统和方法的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现中,方框中所标记的功能也可以以不同于附图中所标记的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或操作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中的技术的改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (18)
1.一种基于均值减损对比归一化信息来识别模糊图像的方法,其特征在于,所述方法包括以下步骤:
对待识别图像进行预处理以得到均值减损对比归一化图像;
对所述均值减损对比归一化图像的统计直方图中的-0.5至0.5范围内的均值减损对比归一化系数进行统计,以得到均值减损对比归一化系数和;
将所述均值减损对比归一化系数和与预定值进行比较;
如果所述均值减损对比归一化系数和大于等于所述预定值,则判别所述待识别图像为模糊图像,否则判别所述待识别图像为清晰图像。
3.根据权利要求1所述的方法,其特征在于,所述预定值为0.9。
4.一种基于梯度幅值和方向值稀疏特征信息以及均值减损对比归一化信息来识别模糊图像的方法,其特征在于,所述方法包括以下步骤:
计算待识别图像的梯度幅值;
将梯度幅值归一化为0~1;
以第一间隔gap0对归一化后的梯度幅值进行量化,即,量化后的梯度幅度范围是0~1/gap0;
统计量化后的梯度幅度的直方图中的0~rate*1/gap0范围内的非空点数N0,其中,rate为比率常数;
计算所述待识别图像在0~360度范围内的梯度方向值;
将梯度方向值归一化为0~1;
以第二间隔gap1对归一化后的梯度方向值进行量化,即,量化后的梯度方向值范围是0~1/gap1;
统计量化后的梯度方向值的直方图中的所有的非空点数N1;
对所述待识别图像进行预处理以得到均值减损对比归一化图像;
对所述均值减损对比归一化图像的统计直方图中的-0.5至0.5范围内的均值减损对比归一化系数进行统计,以得到均值减损对比归一化系数和N2;
基于非空点数N0和N1以及均值减损对比归一化系数和N2来判别所述待识别图像为模糊图像还是为清晰图像;其中,所述基于非空点数N0和N1以及均值减损对比归一化系数和N2来判别所述待识别图像为模糊图像还是为清晰图像的步骤包括:
计算下式:rate0*|N0-center0|/(|N0-center0|+|N0-center0’|)+rate1*|N1-center1|/(|N1-center1|+|N1-center1’|)+(1-rate0-rate1)*|N2-center2|/(|N2-center2|+|N2-center2’|),其中,rate0和rate1是交叉验证参数,center0和center1分别为对训练集中的所有模糊样本计算而得到的非空点数N0和N1的均值,center0’和center1’分别为对训练集中的所有清晰样本计算而得到的非空点数N0和N1的均值,center2为对训练集中的所有模糊样本计算而得到的均值减损对比归一化系数和N2的均值,center2’为对训练集中的所有清晰样本计算而得到的均值减损对比归一化系数和N2的均值;
将计算出的上式与预定值进行比较;
如果计算出的上式小于所述预定值,则判别所述待识别图像为模糊图像,否则判别所述待识别图像为清晰图像。
6.根据权利要求4所述的方法,其特征在于,所述预定值为0.5。
7.根据权利要求4所述的方法,其特征在于,所述归一化为最大最小归一化。
8.根据权利要求4所述的方法,其特征在于,所述第一间隔gap0为0.001,所述比率常数rate为0.01,并且,所述第二间隔gap1为0.00025。
9.一种用于基于均值减损对比归一化信息来识别模糊图像的装置,其特征在于,所述装置包括:
预处理模块,用于对待识别图像进行预处理以得到均值减损对比归一化图像;
均值减损对比归一化系数统计模块,用于对所述均值减损对比归一化图像的统计直方图中的-0.5至0.5范围内的均值减损对比归一化系数进行统计,以得到均值减损对比归一化系数和;
比较模块,用于将所述均值减损对比归一化系数和与预定值进行比较;
判别模块,用于如果所述均值减损对比归一化系数和大于等于所述预定值,则判别所述待识别图像为模糊图像,否则判别所述待识别图像为清晰图像。
11.根据权利要求9所述的装置,其特征在于,所述预定值为0.9。
12.一种用于基于梯度幅值和方向值稀疏特征信息以及均值减损对比归一化信息来识别模糊图像的装置,其特征在于,所述装置包括:
梯度幅值计算模块,用于计算待识别图像的梯度幅值;
梯度幅值归一化模块,用于将梯度幅值归一化为0~1;
梯度幅值量化模块,用于以第一间隔gap0对归一化后的梯度幅值进行量化,即,量化后的梯度幅度范围是0~1/gap0;
梯度幅值统计模块,用于统计量化后的梯度幅度的直方图中的0~rate*1/gap0范围内的非空点数N0,其中,rate为比率常数;
梯度方向值计算模块,用于计算所述待识别图像在0~360度范围内的梯度方向值;
梯度方向值归一化模块,用于将梯度方向值归一化为0~1;
梯度方向值量化模块,用于以第二间隔gap1对归一化后的梯度方向值进行量化,即,量化后的梯度方向值范围是0~1/gap1;
梯度方向值统计模块,用于统计量化后的梯度方向值的直方图中的所有的非空点数N1;
预处理模块,用于对所述待识别图像进行预处理以得到均值减损对比归一化图像;
均值减损对比归一化系数统计模块,用于对所述均值减损对比归一化图像的统计直方图中的-0.5至0.5范围内的均值减损对比归一化系数进行统计,以得到均值减损对比归一化系数和N2;
判别模块,用于基于非空点数N0和N1以及均值减损对比归一化系数和N2来判别所述待识别图像为模糊图像还是为清晰图像;
其中,所述判别模块包括:
判别式计算单元,用于计算下式:
rate0*|N0-center0|/(|N0-center0|+|N0-center0’|)+rate1*|N1-center1|/(|N1-center1|+|N1-center1’|)+(1-rate0-rate1)*|N2-center2|/(|N2-center2|+|N2-center2’|),其中,rate0和rate1是交叉验证参数,center0和center1分别为对训练集中的所有模糊样本计算而得到的非空点数N0和N1的均值,center0’和center1’分别为对训练集中的所有清晰样本计算而得到的非空点数N0和N1的均值,center2为对训练集中的所有模糊样本计算而得到的均值减损对比归一化系数和N2的均值,center2’为对训练集中的所有清晰样本计算而得到的均值减损对比归一化系数和N2的均值;
比较单元,用于将判别式计算单元计算出的上式与预定值进行比较;
综合判别单元,用于如果判别式计算单元计算出的上式小于所述预定值,则判别所述待识别图像为模糊图像,否则判别所述待识别图像为清晰图像。
14.根据权利要求12所述的装置,其特征在于,所述预定值为0.5。
15.根据权利要求12所述的装置,其特征在于,所述归一化为最大最小归一化。
16.根据权利要求12所述的装置,其特征在于,所述第一间隔gap0为0.001,所述比率常数rate为0.01,并且,所述第二间隔gap1为0.00025。
17.一种计算设备,包括:
处理器;以及
存储器,其上存储有可执行代码,当所述可执行代码被所述处理器执行时,使所述处理器执行如权利要求1至8中任何一项所述的方法。
18.一种非暂时性机器可读存储介质,其上存储有可执行代码,当所述可执行代码被电子设备的处理器执行时,使所述处理器执行如权利要求1至8中任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810678755.5A CN108960254B (zh) | 2018-06-27 | 2018-06-27 | 一种基于mscn信息识别模糊图像的方法和装置、计算设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810678755.5A CN108960254B (zh) | 2018-06-27 | 2018-06-27 | 一种基于mscn信息识别模糊图像的方法和装置、计算设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108960254A CN108960254A (zh) | 2018-12-07 |
CN108960254B true CN108960254B (zh) | 2021-12-10 |
Family
ID=64487180
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810678755.5A Active CN108960254B (zh) | 2018-06-27 | 2018-06-27 | 一种基于mscn信息识别模糊图像的方法和装置、计算设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108960254B (zh) |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9277148B2 (en) * | 2012-06-06 | 2016-03-01 | Board Of Regents, The University Of Texas System | Maximizing perceptual quality and naturalness of captured images |
CN103778636B (zh) * | 2014-01-22 | 2017-10-31 | 上海交通大学 | 一种用于无参考图像质量评价的特征构建方法 |
CN105007488A (zh) * | 2015-07-06 | 2015-10-28 | 浙江理工大学 | 基于变换域和空间域的通用无参考图像质量评价方法 |
CN105049851B (zh) * | 2015-07-06 | 2018-01-19 | 浙江理工大学 | 基于色彩感知的通用无参考图像质量评价方法 |
CN105957036B (zh) * | 2016-05-06 | 2018-07-10 | 电子科技大学 | 一种加强字符先验的视频去运动模糊方法 |
CN106447646A (zh) * | 2016-06-28 | 2017-02-22 | 中国人民解放军陆军军官学院 | 一种用于无人机图像的质量盲评价方法 |
CN106485702B (zh) * | 2016-09-30 | 2019-11-05 | 杭州电子科技大学 | 基于自然图像特征统计的图像模糊检测方法 |
-
2018
- 2018-06-27 CN CN201810678755.5A patent/CN108960254B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108960254A (zh) | 2018-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110020592B (zh) | 物体检测模型训练方法、装置、计算机设备及存储介质 | |
US10319096B2 (en) | Automated tattoo recognition techniques | |
WO2019011165A1 (zh) | 人脸识别方法、装置、电子设备及存储介质 | |
US20130039588A1 (en) | Image processing method and apparatus for tamper proofing | |
US20210232846A1 (en) | Image processing method and apparatus for object detection | |
CN109815797B (zh) | 活体检测方法和装置 | |
CN111611863B (zh) | 车牌图像质量评价方法、装置以及计算机设备 | |
CN111144425B (zh) | 检测拍屏图片的方法、装置、电子设备及存储介质 | |
CN108985350B (zh) | 一种基于梯度幅值稀疏特征信息识别模糊图像的方法和装置、计算设备及存储介质 | |
CN109712134B (zh) | 虹膜图像质量评价方法、装置及电子设备 | |
CN116612306B (zh) | 基于计算机视觉的法兰盘智能对位方法及系统 | |
CN108960254B (zh) | 一种基于mscn信息识别模糊图像的方法和装置、计算设备及存储介质 | |
CN111163332A (zh) | 视频色情度检测方法、终端及介质 | |
WO2019177130A1 (ja) | 情報処理装置および情報処理方法 | |
CN113822361B (zh) | 一种基于汉明距离的sar图像相似程度度量方法和系统 | |
CN115830351A (zh) | 图像处理方法、设备以及存储介质 | |
CN108985351B (zh) | 一种基于梯度方向稀疏特征信息识别模糊图像的方法和装置、计算设备及存储介质 | |
Saparudin et al. | Segmentation of fingerprint image based on gradient magnitude and coherence | |
CN113850252A (zh) | 防振锤螺栓检测方法、装置和电子设备 | |
CN111985423A (zh) | 活体检测方法、装置、设备及可读存储介质 | |
Kalaiselvi et al. | A Study on Validation Metrics of Digital Image Processing | |
CN116071658B (zh) | 一种基于深度学习的sar图像小目标检测识别方法及设备 | |
CN113205082B (zh) | 基于采集不确定性解耦的鲁棒虹膜识别方法 | |
CN112329572B (zh) | 一种基于边框和闪光点的快速静态活体检测方法及装置 | |
Paithane et al. | Performance Analysis of Face Detection system using HOG and QualHOG Features |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A method and device, computing device, and storage medium for recognizing fuzzy images based on MSCN information Effective date of registration: 20230324 Granted publication date: 20211210 Pledgee: Haidian Beijing science and technology enterprise financing Company limited by guarantee Pledgor: BEIJING SINOSECU TECHNOLOGY CO.,LTD. Registration number: Y2023110000122 |