一种基于图像特征的PM2.5污染指数的识别方法与流程

本发明属于图像处理及气象学交叉领域，具体涉及一种基于图像特征的pm2.5污染指数的识别方法。
背景技术：
：空气中漂浮着的颗粒物大小不同，pm2.5是其中较细小的那部分。测定pm2.5的浓度一般分为两步：(1)把pm2.5与较大的颗粒物分离，几乎所有的测定方法都需要这一步；(2)测定分离出来的pm2.5的重量。目前，各国环保部门广泛采用的pm2.5测定方法有三种：重量法、β射线吸收法和微量振荡天平法。pm2.5手动监测方法又称标准称重法(滤膜称重法)属于手工监测方法，主要用于研究，或对自动分析方法的校准。这种方法较简便，只需要一个pm2.5切割头、一台泵和膜架及其滤膜，采集24小时样品后，取下滤膜称重即可，必要时可以平行采集3个样品，经恒温恒湿后再称重。滤膜称重法的优点是经济成本低，容易实施。缺点：(1)气流长时间不断通过采样滤膜，滤膜上采集到的物质随着气流和温度的变化也会造成挥发性和半挥发性物质的损失；(2)一些极细小的颗粒穿过滤膜，造成结果偏低；(3)气态物质可能被滤膜吸附，造成结果偏高。自动监测方法的测站点无人值守，数据通过网络直接获得，包括β射线吸收法和微量振荡天平法。β射线吸收法的基本原理是利用堆积在石英滤膜上的颗粒物对c14释放的β射线衰减量的变化检测大气颗粒物质量的变化。环境空气由采样泵经切割器吸入采样管，经过滤膜后排出。颗粒物沉淀在条状石英滤膜上，当β射线通过沉积着颗粒物的滤膜时，β射线强度发生衰减，通过对衰减量的测定计算出颗粒物的浓度。微量振荡天平(teom)法又称微量石英振荡天平法，该方法是在质量传感器内使用一个石英空心锥形管，在空心锥形管振荡端上安放可更换的滤膜，振荡频率取决于石英锥形管特性和它的质量。当采样气流通过滤膜，其中的颗粒物沉积在滤膜上，滤膜质量变化导致振荡频率变化，通过测量振荡频率的变化计算出沉积在滤膜上颗粒物的质量，再根据采样流量、采样现场环境温度和气压计算出该时段的颗粒物标态质量浓度。振荡天平法的优点是定量关系明确。但其存在下面两个缺点：(1)样品加热后挥发性和半挥发性物质的损失，导致测定结果偏低；(2)需要加装膜动态测量系统(filterdynamicmeasurementsystem，)对偏低的结果进行校准。β射线法是基于两个假设：仪器的石英采样滤膜条带均一；采集下来的pm2.5粒子物理特性均一，对β射线强度衰减率相同。但是现实条件下，以上两个假设通常很难成立，因此测定数据一般被认为也存在偏差。这种方法的缺点是在相对干净和干燥的地区故障率低，而在潮湿高温区域故障率很高。技术实现要素：本发明的目的是解决现有pm2.5的测定方法存在操作复杂和测定结果不精确的技术问题，提供一种基于图像特征的pm2.5污染指数的识别方法。为解决上述技术问题，本发明采用的技术方案为：一种基于图像特征的pm2.5污染指数的识别方法，包括以下步骤：步骤1.搜集标注了pm2.5实测值的历史图像，将搜集的历史图像全部统一缩放为xs×ys像素大小；如果图像是彩色的，那么就将其转换为灰度图；步骤2.计算每个缩放并转换为灰度图后图像的图像特征；步骤3.将每个图像的特征作为输入，使用回归方法得到pm2.5的预测模型；步骤4.采集目标图像，按照步骤1中同样的方法将目标图像缩放为xs×ys像素大小，如果图像是彩色的，将其转换为灰度图；步骤5.针对缩放后的目标图像，按照步骤2中的同样的方法，计算缩放并转换为灰度图后图像的图像特征；步骤6.将目标图像的图像特征作为步骤3中已经训练好的预测模型的输入，计算得到目标图像的pm2.5值。所述步骤2中在计算每个缩放并转换为灰度图后图像的图像特征之后,还需要对每个图像的图像特征以及pm2.5的值进行归一化；所述步骤5中针对缩放后的目标图像，在计算缩放并转换为灰度图后图像的图像特征之后，还需要对目标图像的图像特征按照步骤2中同样的归一化方法，进行归一化。在机器学习中，使用同训练时一样的方法对数据进行预处理，即使用同训练时一样的方法将图像转换为灰度值，计算图像特征并进行归一化，在使用先前训练好的预测模型进行预测时才能得到最佳的预测结果。所述步骤1或步骤4中采用公式gray＝0.299×r+0.587×g+0.114×b将彩色图像转换为灰度图像，其中r、g、b分别表示彩色图像的红色、绿色和蓝色波段；所述步骤2或步骤5中的图像特征为图像的平均像素值、局部极值点的个数nmax、局部方差的均值μlvar、图像水平方向和垂直方向二阶差分加权平均值、图像水平方向和垂直方向的二阶梯度和的加权平均值、qx、qy和qg折线差分的加权平均值以及图像中的斑点数；之所以选择这些特征，是因为在实验过程中发现这些特征能够很好的反映实际拍摄图像中pm2.5高低。所述图像的平均像素值为所有像元的像素值加起来之后除以xs×ys；所述图像局部极值点的个数nmax和局部方差的均值μlvar的计算方法为：(1)设定nmax和μlvar的初值为0；(2)对于图像中的每一个非边界像元找出邻近的8个像元，记8个像元的构成集合为pixa；如果当前像元坐标为(i,j),i和j分别表示该像元和图像最左上角像元分别在垂直和水平方向同当前像元间隔的像元个数，那么其8个近邻的坐标分别为(i+1,j)、(i,j+1)、(i-1,j)、(i,j-1)、(i-1,j-1)、(i-1,j+1)、(i+1,j-1)和(i+1,j+1)；这里的非边界像元是指图像中所有能够找到上述8个近邻，并且任一近邻像元像素值不为空的像元；(3)如果当前像元值大于周边的8个像元，那么nmax值加1；并且这里，p表示当前像元，i(p)表示像元p的灰度值，表示pixa中的8个像元和当前像元像素值的均值；(4)当所有像元都处理完之后，μlvar＝μlvar/((xs-2)×(ys-2))；所述图像水平方向和垂直方向二阶差分加权平均值、图像水平方向和垂直方向二阶梯度和的加权平均值的计算方法为：(1)分别计算图像水平方向x和垂直方向y每个像元p的二阶差分,分别记作和并分别将两个方向的二阶差分按照从小到大的序列排序，排序过程中去掉重复值，并过滤掉小于水平方向给定阈值αx或垂直方向给定阈值αy的所有差分值，由此得到两个大小分别为nx和ny的水平方向队列qx和垂直方向队列qy，两个队列中第i的元素分别为和即为队列中排序为i的像元所对应的二阶差分值；(2)计算图像每个像元p的二阶梯度和并将二阶梯度和值按照从小到大的序列排序，排序过程中去掉重复值，并过滤掉小于给定阈值αg的所有二阶梯度和值，由此得到大小为ng的二阶梯度和队列qg，队列中第i的元素分别为qg(i)＝grad(pi)，即为队列中排序为i的像元所对应的二阶梯度和值；(3)分别对三个队列qx、qy和qg计算对应折线总长度和以及总体加权变化度和(4)分别计算三个队列qx、qy和qg的加权平均值wvlx＝vx/lx、wvly＝vy/ly和wvlg＝vg/lg，这三个值就分别为水平方向二阶差分加权平均值、垂直方向二阶差分加权平均值以及二阶梯度和的加权平均值；所述图像中的斑点数的计算方法为：(1)给定的一系列有序尺度值σ1,σ2,…,σs,s≥3，并选定某个半径r，根据规范化拉普拉斯高斯算子生成对应各个尺度的卷积模板，其中每个卷积模板的中心坐标为(0，0)，所有和卷积模板中心水平或者垂直方向相差r个像元的所有网格都属于模板；模板中所有网格的坐标集合为{(tx,ty)|tx＝-r,-r+1,…,r-1,r∧ty＝-r,-r+1,…,r-1,r}，对应尺度σ∈{σ1,σ2,…,σs}的卷积模板中每个网格(tx,ty)上的值使用下面公式计算：对每个尺度σ∈{σ1,σ2,…,σs}，使用对应的卷积模板对图像做卷积运算，得到一系列卷积图像i1,i2,…,is；(2)对于每个卷积图像ii,i＝2,…,s-1中的每个像元p，假设其在图像中的坐标为(i,j)，其8个最近邻的像元坐标分别为(i+1,j)、(i,j+1)、(i-1,j)、(i,j-1)、(i-1,j-1)、(i-1,j+1)、(i+1,j-1)和(i+1,j+1)；如果ii(p)是所有近邻中最小的，并且ii(p)也比ii-1和ii+1中所有坐标为(i,j)、(i+1,j)、(i,j+1)、(i-1,j)、(i,j-1)、(i-1,j-1)、(i-1,j+1)、(i+1,j-1)和(i+1,j+1)的像元的值都小，那么该像元对应的就是一个斑点；如果ii(p)是所有近邻中最大的，并且ii(p)也比ii-1和ii+1中所有坐标为(i,j)、(i+1,j)、(i,j+1)、(i-1,j)、(i,j-1)、(i-1,j-1)、(i-1,j+1)、(i+1,j-1)和(i+1,j+1)的像元的值都大，那么该像元对应的也是一个斑点；在所有卷积图像ii,i＝2,…,s-1上找到的斑点总数即为图像中的斑点数。所述qx、qy和qg折线差分的加权平均值的计算方法为：(1)首先针对qx、qy和qg的折线分别计算相邻两个元素之间的差值，并且将这一差值按照从小到大排序，分别形成一个队列，分别记为qqx、qqy和qqg，其第i个元素分别为qqx(i)、qqy(i)和qqg(i)。(2)分别对三个队列qqx、qqy和qqg计算对应折线总长度和以及总体加权变化度和(3)分别计算三个队列qqx、qqy和qqg的加权平均值wvlqx＝vqx/lqx、wvlqy＝vqy/lqy和wvlqg＝vqg/lqg，这三个值就分别为qx、qy和qg折线差分的加权平均值。所述步骤2和步骤5中按照如下公式对图像的特征和pm2.5的值进行归一化：其中f表示某个特征，f(i)表示缩放并转换为灰度图后的历史图像或者目标图像i在特征f上的取值，minf表示特征f在所有缩放并转换为灰度图后的历史图像上的最小值，maxf表示特征f在所有缩放并转换为灰度图后的历史图像上的最大值,fnorm(i)表示i在特征f上归一化后的值，并且在对转换为灰度图后的目标图像进行归一化时若fnorm(i)＞1时将fnorm(i)置为1，若fnorm(i)＜0时将fnorm(i)置为0,inew为待归一化的新值。这一方法实现起来比较简单高效，因此选择这一方法进行归一化。所述步骤3中使用回归方法得到pm2.5的预测模型的方法为前馈神经网络、高斯过程或常用的最小二乘等方法。所述采用前馈神经网络方法得到pm2.5的预测模型的具体操作步骤为：(1)在前馈神经网络中设定好激活函数，这一步可以使用逻辑斯特函数或者双曲正切函数；(2)设定神经网络的层次和每层的神经元数量；(3)使用反向传播算法学习神经网络；(4)将学习到的网络保存下来，需要保存的内容有使用了何种激活函数，神经网络的层数、每层的神经元数量以及每个神经元的权重；所述步骤6中采用前馈神经网络方法得到pm2.5的预测模型时，需要读入保存的基于前馈神经网络的预测模型参数，包括激活函数，神经网络的层数、每层的神经元数量以及每个神经元的权重，然后按照前馈神经网络的前馈传播方法将归一化后的图像特征作为输入计算出inew所对应的归一化后的pm2.5值pm2.5norm(inew)，最后将其按照下面公式反归一化后作为新拍摄图像的pm2.5值pm2.5(inew)：pm2.5(inew)＝pm2.5norm(inew)×(maxf-minf)+minf。本发明充分利用历史pm2.5的测量值和对应的实景照片，综合使用图像处理技术和数据挖掘技术，只需要使用目标区域景物照片就能够快速估计其pm2.5指数，克服了传统方法通常只能定点测量的缺点，提供了一种能够快速估计大范围区域pm2.5的技术方法。本发明采用以上技术方案，与现有技术相比的优点在于：(1)不需要制作或者使用特别的仪器以及化学试剂，仅使用目标区域的照片就能够识别pm2.5指数；(2)当只有目标区域照片时，无法重建或者重复观察目标区域的环境因素，本方法仍然可以识别其pm2.5指数；(3)能够直接识别大面积区域的平均pm2.5值，而不是定点观测。附图说明图1为本发明的实现流程图；图2-5为本发明实例中的四个输入图像。具体实施方式本实施例中的一种基于图像特征的pm2.5污染指数的识别方法，包括以下步骤：步骤1.选择一个区域，使用摄像头拍摄不同时间的彩色照片并记录这一时间的pm2.5值以搜集标注了pm2.5实测值的历史图像，根据季节和时间将所有的照片区分为白天和晚上两大类，下面计算都以白天为例；将搜集的历史图像全部统一缩放为xs×ys像素大小，这里使用imagemagick软件的convert命令：“convert–resize640x320”将搜集的历史图像全部统一缩放为640×320像素大小；这里使用只使用彩色照片的红色波段的方式将图像转换为灰度图，也可以采用公式gray＝0.299×r+0.587×g+0.114×b将图像转换为灰度图，其中r、g、b分别表示彩色图像的红色、绿色和蓝色波段；步骤2.计算每个缩放并转换为灰度图后图像的平均像素值、局部极值点的个数nmax、局部方差的均值μlvar、图像水平方向和垂直方向二阶差分加权平均值、图像水平方向和垂直方向的二阶梯度和的加权平均值、qx、qy和qg折线差分的加权平均值以及图像中的斑点数；所述图像的平均像素值为所有像元的像素值加起来之后除以xs×ys；所述图像局部极值点的个数nmax和局部方差的均值μlvar的计算方法为：(1)设定nmax和μlvar的初值为0；(2)对于图像中的每一个非边界像元找出邻近的8个像元，记8个像元的构成集合为pixa；如果当前像元坐标为(i,j),i和j分别表示该像元和图像最左上角像元分别在垂直和水平方向同当前像元间隔的像元个数，那么其8个近邻的坐标分别为(i+1,j)、(i,j+1)、(i-1,j)、(i,j-1)、(i-1,j-1)、(i-1,j+1)、(i+1,j-1)和(i+1,j+1)；这里的非边界像元是指图像中所有能够找到上述8个近邻，并且任一近邻像元像素值不为空的像元；(3)如果当前像元值大于周边的8个像元，那么nmax值加1；并且这里，p表示当前像元，i(p)表示像元p的灰度值，表示pixa中的8个像元和当前像元像素值的均值；(4)当所有像元都处理完之后，μlvar＝μlvar/((xs-2)×(ys-2))；所述图像水平方向和垂直方向二阶差分加权平均值、图像水平方向和垂直方向二阶梯度和的加权平均值的计算方法为：(1)分别计算图像水平方向x和垂直方向y每个像元p的二阶差分,分别记作和并分别将两个方向的二阶差分按照从小到大的序列排序，排序过程中去掉重复值，并过滤掉小于水平方向给定阈值αx或垂直方向给定阈值αy的所有差分值，由此得到两个大小分别为nx和ny的水平方向队列qx和垂直方向队列qy，两个队列中第i的元素分别为和即为队列中排序为i的像元所对应的二阶差分值；(2)计算图像每个像元p的二阶梯度和并将二阶梯度和值按照从小到大的序列排序，排序过程中去掉重复值，并过滤掉小于给定阈值αg的所有二阶梯度和值，由此得到大小为ng的二阶梯度和队列qg，队列中第i的元素分别为qg(i)＝grad(pi)，即为队列中排序为i的像元所对应的二阶梯度和值；(3)分别对三个队列qx、qy和qg计算对应折线总长度和以及总体加权变化度和(4)分别计算三个队列qx、qy和qg的加权平均值wvlx＝vx/lx、wvly＝vy/ly和wvlg＝vg/lg，这三个值就分别为水平方向二阶差分加权平均值、垂直方向二阶差分加权平均值以及二阶梯度和的加权平均值；所述qx、qy和qg折线差分的加权平均值的计算方法为：(1)首先针对qx、qy和qg的折线分别计算相邻两个元素之间的差值，并且将这一差值按照从小到大排序，分别形成一个队列，分别记为qqx、qqy和qqg，其第i个元素分别为qqx(i)、qqy(i)和qqg(i)。(2)分别对三个队列qqx、qqy和qqg计算对应折线总长度和以及总体加权变化度和(3)分别计算三个队列qqx、qqy和qqg的加权平均值wvlqx＝vqx/lqx、wvlqy＝vqy/lqy和wvlqg＝vqg/lqg，这三个值就分别为qx、qy和qg折线差分的加权平均值。所述图像中的斑点数的计算方法为：(1)给定的一系列有序尺度值σ1,σ2,…,σs,s≥3，并选定某个半径r，根据规范化拉普拉斯高斯算子生成对应各个尺度的卷积模板，其中每个卷积模板的中心坐标为(0，0)，所有和卷积模板中心水平或者垂直方向相差r个像元的所有网格都属于模板；模板中所有网格的坐标集合为{(tx，ty)|tx＝-r，-r+1，…，r-1，r∧ty＝-r，-r+1，…，r-1，r}，对应尺度σ∈{σ1,σ2,…,σs}的卷积模板中每个网格(tx，ty)上的值使用下面公式计算：对每个尺度σ∈{σ1,σ2,…,σs}，使用对应的卷积模板对图像做卷积运算，得到一系列卷积图像i1,i2,…,is；(2)对于每个卷积图像ii,i＝2,…,s-1中的每个像元p，假设其在图像中的坐标为(i,j)，其8个最近邻的像元坐标分别为(i+1,j)、(i,j+1)、(i-1,j)、(i,j-1)、(i-1,j-1)、(i-1,j+1)、(i+1,j-1)和(i+1,j+1)；如果ii(p)是所有近邻中最小的，并且ii(p)也比ii-1和ii+1中所有坐标为(i,j)、(i+1,j)、(i,j+1)、(i-1,j)、(i,j-1)、(i-1,j-1)、(i-1,j+1)、(i+1,j-1)和(i+1,j+1)的像元的值都小，那么该像元对应的就是一个斑点；如果ii(p)是所有近邻中最大的，并且ii(p)也比ii-1和ii+1中所有坐标为(i,j)、(i+1,j)、(i,j+1)、(i-1,j)、(i,j-1)、(i-1,j-1)、(i-1,j+1)、(i+1,j-1)和(i+1,j+1)的像元的值都大，那么该像元对应的也是一个斑点；在所有卷积图像ii,i＝2,…,s-1上找到的斑点总数即为图像中的斑点数；本实例选择了计算局部极值点的个数nmax、局部方差的均值μlvar、平均像素值，基于拉普拉斯算子的二阶梯度和加权平均值和qg折线差分的加权平均值、基于sobel算子的qy折线差分的加权平均值以及基于高斯拉普拉斯算子的图像中的斑点数这些特征，计算出的图像部分特征如下表所示：计算高斯拉普拉斯图像斑点数时使用了opencv库，其部分代码如下：matimage＝img；matdescriptors；vector<keypoint>keypoints；simpleblobdetector::paramsparams；params.minthreshold＝10；params.maxthreshold＝100；params.thresholdstep＝10；params.minarea＝10；params.minconvexity＝0.3；params.mininertiaratio＝0.01；params.maxarea＝8000；params.maxconvexity＝10；params.filterbycolor＝false；params.filterbycircularity＝false；simpleblobdetectorblobdetector(params)；blobdetector.create("simpleblob")；blobdetector.detect(image,keypoints)；drawkeypoints(image,keypoints,image,scalar(255,0,0))；cout<<keypoints.size()<<""；计算基于拉普拉斯算子的二阶梯度和加权平均值和qg折线差分的加权平均值也使用了opencv库，其部分代码如下：这里αg阈值设定为50。并对每个图像的各图像特征以及pm2.5的值按照如下公式进行归一化：其中f表示某个特征，f(i)表示缩放并转换为灰度图后的历史图像或者目标图像i在特征f上的取值，minf表示特征f在所有缩放并转换为灰度图后的历史图像上的最小值，maxf表示特征f在所有缩放并转换为灰度图后的历史图像上的最大值,fnorm(i)表示i在特征f上归一化后的值，并且在对转换为灰度图后的目标图像进行归一化时若fnorm(i)＞1时将fnorm(i)置为1，若fnorm(i)＜0时将fnorm(i)置为0,inew为待归一化的新值。以平均像素值为例，训练数据中这一特征的最大值为147.323，最小值为49.6862，因此，当这一值为113.21时，其归一化后的值为(113.21-49.6862)/(147.323-49.6862)＝0.181324。下表为各个特征和pm2.5归一化之后的结果：步骤3.将每个图像的特征作为输入，使用前馈神经网络方法得到pm2.5的预测模型，其具体操作步骤为：(1)在前馈神经网络中设定好激活函数，这一步可以使用逻辑斯特函数或者双曲正切函数，这里设置激活函数为逻辑斯特函数，也叫做sigmoid函数；(2)设定神经网络的层次和每层的神经元数量，这里设置网络层数为3，其中输入层神经元个数为7,输出层神经元个数为1，中间层神经元个数为1000；(3)使用反向传播算法学习神经网络，这里使用的参数最低训练速率为0.001，最高训练速率为0.01，误差区间为[0.001,0.02],迭代次数为10000；神经网络采用了lwneuralnetplus0.998版本的库，部分训练代码如下：iomanager＝newiomanagelwnnfann()；iomanager->info_from_file(argv[1],&npatterns,&ninput,&noutput)；net＝newnetwork(network::logistic,3,ninput,1000,noutput)；net->set_momentum(0)；net->set_learning_rate(0.001)；net->jolt(0.03,0.22)；train＝newtrainer(net,"","")；train->set_iomanager(iomanager)；cout<<"loadingtrainingdata..."<<endl；train->load_training(argv[1])；train_go(argv,train,net)；其中train_go部分代码如下所示：(4)将学习到的网络保存下来，需要保存的内容有使用了何种激活函数，神经网络的层数、每层的神经元数量以及每个神经元的权重；步骤4.采集目标图像，按照步骤1中同样的方法将目标图像缩放为xs×ys像素大小，如果图像是彩色的，将其转换为灰度图；步骤5.针对缩放后的目标图像，按照步骤2中的同样的方法，计算缩放并转换为灰度图后图像的平均像素值、局部极值点的个数nmax、局部方差的均值μlvar、图像水平方向和垂直方向二阶差分加权平均值、图像水平方向和垂直方向的二阶梯度和的加权平均值、qx、qy和qg折线差分的加权平均值以及图像中的斑点数，并对目标图像的各图像特征按照步骤2中同样的归一化方法，进行归一化；步骤6.将目标图像的图像特征作为步骤3中已经训练好的预测模型的输入，读入保存的基于前馈神经网络的预测模型参数，包括激活函数，神经网络的层数、每层的神经元数量以及每个神经元的权重，然后按照前馈神经网络的前馈传播方法将归一化后的图像特征作为输入计算出inew所对应的归一化后的pm2.5值pm2.5norm(inew)，最后将其按照下面公式反归一化后作为新拍摄图像的pm2.5值pm2.5(inew)：pm2.5(inew)＝pm2.5norm(inew)×(maxf-minf)+minf。我们选用了下图中的四个图像作为新拍摄图像进行识别。这里需要先对下图中的各个图像的特征进行归一化，然后带入先前训练好的神经网络计算归一化后的pm2.5值，最后进行反归一化即可。识别效果如下表所示：实测值预测值图像13749图像2223247图像3323306图像4117103所述步骤2和步骤5中不论使用何种归一化方法，只要能够将特征值或者pm2.5值变为[0,1]区间或者(0,1)区间的数，就属于本专利权利要求范围；所述步骤3中的预测模型的可选建立方法包括神经网络、高斯过程以及常用的最小二乘等方法等，不论使用何种方法，只要是在预测过程中使用了神经网络、高斯过程以及常用的最小二乘等方法中的一种或者多种特征都在本专利的权利要求范围；所述步骤2和步骤5中，计算图像水平方向和垂直方向二阶差分加权平均值以及二阶梯度和的加权平均值可以使用多种方法，重复计算并获得多种特征，例如分别使用sobel算子和拉普拉斯算子计算这三个特征；所述步骤2和步骤5的归一化为可选步骤，如果不使用步骤2和步骤5，则需要再步骤6中预测出的pm2.5值无需进行反归一化操作。本专利中的预测目标pm2.5的值也可以为能见度等可以使用实数进行度量的其他实际预测目标，不论使用何种预测目标，只要使用了图像水平方向和垂直方向二阶差分加权平均值、二阶梯度和的加权平均值以及qx、qy和qg折线差分的加权平均值这几个特征就都属于本专利的权利要求范围。以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。当前第1页12