CN105447100B - 一种基于形状特征的云图检索方法 - Google Patents

Abstract

本发明公开了一种基于形状特征的云图检索方法，有效结合云图云顶亮温信息，采用迭代阈值分割方法分割云系，同时采用计算复杂度较低且鲁棒性强的几何不变矩提取云系形状特征克服了传统的基于形状特征的云图检索方法在分析云图时的困难，特别是方法普适性、计算复杂度、鲁棒性等问题，此外，本发明将底层云图特征从欧氏空间转换到了稀疏空间，利用稀疏空间的分布规律对云图数据库进行检索，有效地缓解了传统检索方法的语义鸿沟问题。

Description

一种基于形状特征的云图检索方法

技术领域

本发明涉及一种云图的检索方法，尤其是涉及一种基于形状特征的云图检索方法。

背景技术

卫星云图能从多角度展示各类云系特征及其演变过程，在天气实况监测、气候研究或灾难天气救灾决策等方面具有重要的意义。卫星云图中各云系往往对应着不同的天气信息，一般若两幅云图在灰度、形状和纹理等特征上相似，则两幅云图在该阶段的天气状况存在一定相似度并具有借鉴意义。同时，随着气象卫星技术的发展，各资料站每天能够接收几乎覆盖全球的GB数量级的海量云图数据，传统的云图人工标注方法已捉襟见肘。如何更为有效借助历史云图信息并投入实际领域，吸引了大批气象卫星工作研究者进行科学研究。因此，如何设计一个高效的基于内容的云图检索系统将成为破解当前难题的关键技术，已引起了许多学者的关注。

在基于内容的云图检索研究中要实现云图较为精确的检索，必须找到能有效逼近云图高层语义内容的底层特征。就图像特征而言，灰度、纹理和形状特征是图像最重要的三类特征，图像的形状特征尽管有别于灰度和纹理特征，却是对灰度或纹理特征的几何重现，认知心理学研究发现，形状在检索过程中显现出来的作用比灰度和纹理更为重要。同时云形状多种多样，变化也较复杂，如高积云形状差异大，常以椭圆形或水波状密集云条存在，卷云一般呈现丝缕状、团簇状，台风等热带气旋常呈现巨大云团涡旋等。因此，云图在形状方面具有较强的区分度，研究一种基于形状特征的云图检索方法具有重要的意义。

当前基于内容的卫星云图检索研究处于起步阶段，国内外尚未形成可靠的利用形状特征的云图检索方案。日本研究了可变形的椭圆来提取台风云系的云团形状特征，并根据椭圆位置的变化来表示运动轨迹，设计了一种面向台风云图的检索系统：Kitamoto，但该系统只能针对台风云系，不具有一般性；意大利的Acqua等人利用点扩散技术，通过位置、旋转度和尺度来刻画云图的形状特征，研究了针对飓风与非飓风云系的检索系统，该系统同样存在应用面较窄的问题且计算复杂；印度的Deepak利用云团的面积和周长作为形状特征，实现了一个云图检索系统，但实验表明该类特征也仅对台风云系适应性较好。在国内，李艳兵等人利用圆形的“可变形模型”描述云团形状，提取多云块参数特征，但圆形提取和匹配的复杂性都很大，计算比较费时。哈尔滨工程大学的上官伟采用变分法提取云图形状特征，通过模糊相似度计算，实现了一种云图检索系统，但系统的检索精度不高；南京航空航天大学融合粒子群优化方法和FCM方法，首先对云系进行前期聚类得到积云云系，再由几何不变矩提取形状特征，进行相似度匹配，实现云图检索，但当该方法面对复杂云系时，难于提取全面的形状特征，并且相似度匹配也不易实现。

为克服传统基于形状特征的云图检索方法在分析云图时的困难，特别是方法普适性、计算复杂度、鲁棒性等问题，有必要引入信号与图像处理领域的新理念、新技术。近年来，随着压缩感知技术在信号处理及成像等领域的成功应用，稀疏表示理论受到了各国学者的广泛关注。研究表明，稀疏表示符合灵长类动物大脑视觉皮层对复杂刺激的感知过程，可作为神经信息群体分布式表达的有效策略，更加符合人眼的视觉特性，正是由于稀疏表示所具有的特征保持性和稀疏性，以冗余字典为代表的稀疏表示理论被应用于信号处理的各个领域，特别是其在人脸识别中的成功应用，为在图像分类和检索领域的应用提供了重要参考。

发明内容

本发明所要解决的技术问题是提供一种计算复杂度较低、鲁棒性强且检索精度高的基于形状特征的云图检索方法。

本发明解决上述技术问题所采用的技术方案为：一种基于形状特征的云图检索方法，包括以下步骤：

(1)将静止卫星云图接收系统所接收到的灰度区间为[0,1024]的云图，依据云图接收系统所对应的云图灰度与云顶亮温关系对照表，将云图灰度转化为云顶亮温；

(2)完成云图亮温空间转换后，再采用迭代阈值分割方法进行云图分割；

(3)通过开、闭运算处理，去除云图的云系空洞和边界毛刺；

(4)对开、闭运算处理后的整幅云图采用100*100大小进行不重叠分块处理，得到云图各子块的区域信息；

(5)对整幅云图分块处理后，采用几何不变矩对各子块进行区域特征描述，将提取的所有子块的特征向量依次排为一长列向量作为该幅云图的形状特征向量，表示为K＝[f₁…f_i…f_m]^T，其中K∈R^m，表示K是实数空间的m维向量，m表示形状特征列向量的维数，T为向量的转置，f_i为K中的分量，i为区间[1，m]中的正整数；

(6)对于该幅云图，首先统计得到形状特征向量K中的最大值和最小值，分别记为F_max和F_min，再对形状特征向量K中的任意分量f_i进行归一化处理，表示如下：

将形状特征向量K归一化到[0,1]范围，得到新的云图形状特征向量，表示为；H＝[f₁'…f_i'…f_m']^T，H∈R^m；

(7)将通过上述步骤提取的所有云图的形状特征，建立一个云图形状特征库A＝[H₁,H₂,...H_j...,H_w]，A∈R^m×w，其中w表示云图库中总的云图数量，H_j表示任意一幅云图的特征向量，j为区间[1，w]中的正整数；

(8)对于一幅待检索云图，采用步骤(1)～(6)的方法提取形状特征向量，表示为：y∈R^m；然后采用稀疏分类方法将待检索云图的形状特征与云图形状特征库进行相似度计算，得到最接近的9幅云图，完成云图检索。

采用迭代阈值分割方法进行云图分割的具体方法为：

①先用如下公式计算云图初始亮温阈值iThreshold：

其中，max_value_gray表示为整幅云图中的最大亮温值，min_value_gray表示整幅云图中的最小亮温值；

②根据初始亮温阈值iThreshold，采用下式统计云图的重要云系成分Z₀和云系背景Z₁：

式中Z₀表示重要云系成分，Z₁表示云系背景，(x,y)表示云图中的像素坐标点，f(x,y)表示云图在坐标点(x,y)处亮温大小，N(x,y)表示权重大小，设定为1，并设定初始亮温阈值iThreshold大小为175；

③根据Z₀和Z₁重新求解新亮温阈值iNewThreshold：

④如果新亮温阈值iNewThreshold与初始亮温阈值iThreshold差值在50之内，则继续下一步骤；否则用iNewThreshold作为初始亮温阈值iThreshold，转到步骤②；

⑤根据新亮温阈值iNewThreshold对原始云图进行分割，即当云图中像素的亮温大于iNewThreshold时设定为1，否则 设定为0，得到黑白二值图像。

去除云图的云系空洞和边界毛刺的具体步骤为：

①对阈值分割后的云图，先采用半径为5的平坦型圆盘结构元素进行膨胀处理，再通过半径为7的平坦型圆盘结构元素进行腐蚀处理，通过闭运算去除云系内部大部分细小空洞；

②对闭运算处理后的云图，先采用半径为6的平坦型圆盘结构元素进行腐蚀处理，再通过半径为7的平坦型圆盘结构元素进行膨胀处理，通过开运算去除边界毛刺。

采用几何不变矩对各子块进行区域特征描述的具体方法为：

①首先计算各云图子块的区域几何中心坐标

②计算云图各子块的中心矩μ_pq：式中p，q满足p +q＝2或者p+q＝3，且p，q为正整数；

③对中心矩μ_pq进行规格化处理，式中满足γ＝(p+q)/2；

④根据规格化后的中心矩η_pq，提取各子块的7个形状不变矩I_k，其中k＝1，2，3，4，5，6，7：I₁＝η₂₀+η₀₂，I₃＝(η₃₀-3η₁₂)²+(3η₂₁-η₀₃)²，I₄＝(η₃₀+η₁₂)²+(η₂₁+η₀₃)²，

I₅＝(η₃₀-3η₁₂)(η₃₀+η₁₂)[(η₁₂+η₃₀)²-3(η₂₁+η₃₀)²]

+(3η₂₁-η₀₃)(η₂₁+η₀₃)[3(η₁₂+η₃₀)²-(η₂₁+η₃₀)²]，

I₆＝(η₂₀-η₀₂)[(η₁₂+η₃₀)²-(η₂₁+η₀₃)²]+4η₁₁(η₁₂+η₃₀)(η₂₁+η₀₃)，

I₇＝(3η₂₁-η₀₃)(η₃₀+η₁₂)[(η₁₂+η₃₀)²-3(η₂₁+η₀₃)²]

+(η₃₀-2η₁₂)(η₂₁+η₀₃)[3(η₁₂+η₃₀)²-(η₂₁+η₀₃)²]，

将所得的7个形状不变矩排为一列作为该云图子块的形状不变矩特征；

⑤当完成整幅云图所有子块形状不变矩特征的提取后，将所有子块的特征向量依次排为一长列向量作为该幅云图的形状特征向量，表示为向量K，K∈R^m。

采用稀疏分类方法进行云图检索相似度计算的具体步骤为：

①首先对待检索云图和云图特征库进行空间投影，采用正态分布随机矩阵R∈R^n×m(n<m)对云图特征库进行空间映射如下：同样对待检索云图特征也进行映射处理如下：

②根据稀疏分类算法，计算残差，进行类别判断：其中a_s表示稀疏系数中第s幅云图数据的对应系数；

③重复第①和②步骤，通过多次迭代，求取重构残差均值：其中l表示迭代次数，表示第t次迭代时属于s幅云图样本的残差值；

④对残差均值E[r_s]进行大小排序，获取与待检索云图的相似度信息，其中残差越小，相似性越大，最后返回残差最小的9幅云图。

与现有技术相比，本发明的优点在于克服了传统的基于形状特征的云图检索方法在分析云图时的困难，特别是方法普适性、计算复杂度、鲁棒性等问题，有效结合云图云顶亮温信息，采用迭代阈值分割方法分割云系，同时采用计算复杂度较低且鲁棒性强的几何不变矩提取云系形状特征；此外，本发明将底层云图特征从欧氏空间转换到了稀疏空间，利用稀疏空间的分布规律对云图数据库进行检索，有效地缓解了传统检索方法的语义鸿沟问题。

附图说明

图1为针对普通云图本发明与现有技术的检索结果示意图；其中(a)为本发明方法，(b)为现有技术的HOG特征检索算法，(c)为现有技术的Zernike矩形状描述的检索算法；

图2为针对特殊天气本发明与现有技术的检索结果示意图；其中(a)为本发明方法，(b)为现有技术的HOG特征检索算法，(c)为现有技术的Zernike矩形状描述的检索算法；

图3为本发明与现有技术云图旋转后的检索结果示意图；其中(a)为本发明方法，(b)为现有技术的HOG特征检索算法，(c)为现有技术的Zernike矩形状描述的检索算法；

图4为本发明与现有技术云图平移后的检索结果示意图；其中(a)为本发明方法，(b)为现有技术的HOG特征检索算法，(c)为现有技术的Zernike矩形状描述的检索算法。

具体实施方式

以下结合附图实施例对本发明作进一步详细描述。

一种基于形状特征的云图检索方法，包括以下步骤：

(1)将静止卫星云图接收系统所接收到的灰度区间为[0,1024]的云图，依据云图接收系统所对应的云图灰度与云顶亮温关系对照表，将云图灰度转化为云顶亮温；

(2)完成云图亮温空间转换后，再采用迭代阈值分割方法进行云图分割，具体方法为：

(2)－①先用如下公式计算云图初始亮温阈值iThreshold：

其中，max_value_gray表示为整幅云图中的最大亮温值，min_value_gray表示整幅云图中的最小亮温值；

(2)－②根据初始亮温阈值iThreshold，采用下式统计云图的重要云系成分Z₀和云系背景Z₁：

式中Z₀表示重要云系成分，Z₁表示云系背景，(x,y)表示云图中的像素坐标点，f(x,y)表示云图在坐标点(x,y)处亮温大小，N(x,y)表示权重大小，设定为1，并设定初始亮温阈值iThreshold大小为175；

(2)－③根据Z₀和Z₁重新求解新亮温阈值iNewThreshold：

(2)－④如果新亮温阈值iNewThreshold与初始亮温阈值iThreshold差值在50之内，则继续下一步骤；否则用iNewThreshold作为初始亮温阈值iThreshold，转到步骤②；

(2)－⑤根据新亮温阈值iNewThreshold对原始云图进行分割，即当云图中像素的亮温大于iNewThreshold时设定为1，否则 设定为0，得到黑白二值图像；

(3)通过开、闭运算处理，去除云图的云系空洞和边界毛刺，具体步骤为：

(3)－①对阈值分割后的云图，先采用半径为5的平坦型圆盘结构元素进行膨胀处理，再通过半径为7的平坦型圆盘结构元素进行腐蚀处理，通过闭运算去除云系内部大部分细小空洞；

(3)－②对闭运算处理后的云图，先采用半径为6的平坦型圆盘结构元素进行腐蚀处理，再通过半径为7的平坦型圆盘结构元素进行膨胀处理，通过开运算去除边界毛刺；

(4)对开、闭运算处理后的整幅云图采用100*100大小进行不重叠分块处理，得到云图各子块的区域信息；

(5)对整幅云图分块处理后，采用几何不变矩对各子块进行区域特征描述，具体方法为：

(5)－①首先计算各云图子块的区域几何中心坐标

(5)－②计算云图各子块的中心矩μ_pq：式中p，q 满足p+q＝2或者p+q＝3，且p，q为正整数；

(5)－③对中心矩μ_pq进行规格化处理，式中满足γ＝(p+q)/2；

(5)－④根据规格化后的中心矩η_pq，提取各子块的7个形状不变矩I_k，其中k＝1，2，3，4，5，6，7：I₁＝η₂₀+η₀₂，I₃＝(η₃₀-3η₁₂)²+(3η₂₁-η₀₃)²，

I₄＝(η₃₀+η₁₂)²+(η₂₁+η₀₃)²，

I₅＝(η₃₀-3η₁₂)(η₃₀+η₁₂)[(η₁₂+η₃₀)²-3(η₂₁+η₃₀)²]

+(3η₂₁-η₀₃)(η₂₁+η₀₃)[3(η₁₂+η₃₀)²-(η₂₁+η₃₀)²]，

I₆＝(η₂₀-η₀₂)[(η₁₂+η₃₀)²-(η₂₁+η₀₃)²]+4η₁₁(η₁₂+η₃₀)(η₂₁+η₀₃)，

I₇＝(3η₂₁-η₀₃)(η₃₀+η₁₂)[(η₁₂+η₃₀)²-3(η₂₁+η₀₃)²]

+(η₃₀-2η₁₂)(η₂₁+η₀₃)[3(η₁₂+η₃₀)²-(η₂₁+η₀₃)²]，

将所得的7个形状不变矩排为一列作为该云图子块的形状不变矩特征；

(5)－⑤当完成整幅云图所有子块形状不变矩特征的提取后，将所有子块的特征向量依次排为一长列向量作为该幅云图的形状特征向量，表示为向量K，K∈R^m，将提取的所有子块的特征向量依次排为一长列向量作为该幅云图的形状特征向量，表示为K＝[f₁…f_i…f_m]^T，其中K∈R^m，表示K是实数空间的m维向量，m表示形状特征列向量的维数，T为向量的转置，f_i为K中的分量，i为区间[1，m]中的正整数；

(6)对于该幅云图，首先统计得到形状特征向量K中的最大值和最小值，分别记为F_max和F_min，再对形状特征向量K中的任意分量f_i进行归一化处理，表示如下：

将形状特征向量K归一化到[0,1]范围，得到新的云图形状特征向量，表示为；H＝[f₁'…f_i'…f_m']^T，H∈R^m；

(7)将通过上述步骤提取的所有云图的形状特征，建立一个云图形状特征库A＝[H₁,H₂,...H_j...,H_w]，A∈R^m×w，其中w表示云图库中总的云图数量，H_j表示任意一幅云图的特征向量，j为区间[1，w]中的正整数；

(8)对于一幅待检索云图，采用步骤(1)～(6)的方法提取形状特征向量，表示为：y∈R^m；然后采用稀疏分类方法将待检索云图的形状特征与云图形状特征库进行相似度计算，具体步骤为：

①首先对待检索云图和云图特征库进行空间投影，采用正态分布随机矩阵R∈R^n×m(n<m)对云图特征库进行空间映射如下：同样对待检索云图特征也进行映射处理如下：

②根据稀疏分类算法，计算残差，进行类别判断：其中a_s表示稀疏系数中第s幅云图数据的对应系数；

③重复第①和②步骤，通过多次迭代，求取重构残差均值：其中l表示迭代次数，表示第t次迭代时属于s幅云图样本的残差值；

④对残差均值E[r_s]进行大小排序，获取与待检索云图的相似度信息，其中残差越小，相似性越大，最后返回残差最小的9幅云图。

为检验本发明基于形状特征的云图检索方法的性能，以下设置了多组对比试验，分别采用经典描述形状特征的HOG特征检索方法^[1]，记为HOG_CBIR；采用Zernike矩形状描述的检索方法^[2]，记为Zernike_CBIR；本发明基于形状特征的检索方法，记为HU_SRC。

在图像检索中，为了评价各类检索算法的有效性，需要算法性能评价准则。图像检索算法中公认度较高的性能评价准则有查全率、查准率和算法时间复杂度。查全率(Recall)是指检索所得的图像中与检索图像相关的图像数占整个数据库中相关图像数的比例，查准率(Precision)是指检索所得的图像中与检索图像相关的图像数占返回图像数目的比例，一般而言，Recall和Precision的值越高，代表检索算法更优异，但两值为矛盾值，当返回图像数越多，查全率越高，查准率反而越低，一般存在一定条件使两者达到最大。假设m₁代表检索数据库中与检索图像相关的图像个数，m₂代表相关但并未检索到的图像个数，m₃代表查询返回的图像中不相关的图像个数，则查全率、查准率可表示如下：

由上可知，当返回图像数越多，查全率越高，而此时的查准率会呈下降趋势，因此一般在检索系统中存在查全率和查准率都较高的平衡点，此时也是该检索算法的最高性能。为此，为平衡查全率、查准率，同时使本发明的图像检索算法性能更为精确，本发明又增加了信息检索中的综合评价指标F-Measure，又称F-Score，可表示为：

此外，本发明采用了算法的时间复杂度进行检索性能的评价，即所需的查询时间。

例一：采用2014年04月10日00点32分普通天气数据作为待检索图片，其中与该天云图比较相似的云图共计15幅。检索效果如图1所示。

由检索的实验结果可以看出，相对于其他形状检索方法，本发明检索方法的效果更为明显，检索查准率和综合指标较高，可见本发明基于形状的云图检索方法具有可行性。

例二：由于台风天气云图可以检验本发明方法的普适性，同时减少主观误差，利于观察，因此选用了2013年11月01日08点30分台风“罗莎”天气数据作为待检索图片，其中云图数据库中与该天云图比较相似的云图共计12幅，实验结果如下。

表1不同方法的检索性能比较

TBB:表示云顶亮温。

如下显示为各类方法的检索结果，其中返回的9幅云图相似度从大到小排列，如图2所示。

例三：为评价本发明方法旋转、平移鲁棒性能，对云图分别采取旋转10°、右下平移70个像素后进行检索比较，各方法的检索结果如图3和图4所示，由检索结果可知，所有检索方法基本能检索到旋转、平移处理后的云图，并作为最相似图像返回，这与事实一致；此外由显示的云图检索结果可知，相对于其他检索方法，本发明方法检索效果更好，检索图像基本为与待检索云图时间间隔相邻且相似的台风天气，其他方法在云图旋转平移下也能大致检索到相似云图，但检索失误率较大。总体来说，本发明方法具有较好的抗平移、旋转特性。

参考文献：

[1]董俊杰.基于HOG和SVM的服装图像检索系统的设计与实现[D].中山大学.2014.

[2]曹闻,万明英,李润生,等.基于Zernike形状矩的地图匹配算法[J].计算机应用研究,2011,28(7):2786-2792.

Claims (5)

1.一种基于形状特征的云图检索方法，其特征在于包括以下步骤：

(1)将静止卫星云图接收系统所接收到的灰度区间为[0,1024]的云图，依据云图接收系统所对应的云图灰度与云顶亮温关系对照表，将云图灰度转化为云顶亮温；

(2)完成云图亮温空间转换后，再采用迭代阈值分割方法进行云图分割；

(3)通过开、闭运算处理，去除云图的云系空洞和边界毛刺；

(4)对开、闭运算处理后的整幅云图采用100*100大小进行不重叠分块处理，得到云图各子块的区域信息；

(5)对整幅云图分块处理后，采用几何不变矩对各子块进行区域特征描述，将提取的所有子块的特征向量依次排为一长列向量作为该幅云图的形状特征向量，表示为K＝[f₁…f_i…f_m]^T，其中K∈R^m，表示K是实数空间的m维向量，m表示形状特征列向量的维数，T为向量的转置，f_i为形状特征向量K中的第i个分量，i为区间[1，m]中的正整数；

(6)对于该幅云图，首先统计得到形状特征向量K中的最大值和最小值，分别记为F_max和F_min，再对形状特征向量K中的任意第i个分量f_i进行归一化处理，表示如下：

将形状特征向量K归一化到[0,1]范围，得到新的云图形状特征向量，表示为；H＝[f₁'Kf_i'K f′_m]^T，H∈R^m；

(7)将通过上述步骤提取的所有云图的形状特征，建立一个云图形状特征库A＝[H₁,H₂,...H_j...,H_w]，A∈R^m×w，表示A是m行w列的矩阵，其中w表示云图库中总的云图数量，云图形状特征库H_j表示任意一幅云图的特征向量，j为区间[1，w]中的正整数；

(8)对于一幅待检索云图，采用步骤(1)～(6)的方法提取形状特征向量，表示为：y∈R^m，表示y是维度为m的列向量；然后采用稀疏分类方法将待检索云图的形状特征与云图形状特征库进行相似度计算，得到最接近的9幅云图，完成云图检索。

2.如权利要求1所述的一种基于形状特征的云图检索方法，其特征在于采用迭代阈值分割方法进行云图分割的具体方法为：

①先用如下公式计算云图初始亮温阈值iThreshold：

其中，max_value_gray表示为整幅云图中的最大亮温值，min_value_gray表示整幅云图中的最小亮温值；

②根据初始亮温阈值iThreshold，采用下式统计云图的重要云系成分Z₀和云系背景Z₁：

式中Z₀表示重要云系成分，Z₁表示云系背景，(x,y)表示云图中的像素坐标点，f(x,y)表示云图在坐标点(x,y)处亮温大小，N(x,y)表示权重大小，设定为1，并设定初始亮温阈值iThreshold大小为175；

③根据Z₀和Z₁重新求解新亮温阈值iNewThreshold：

④如果新亮温阈值iNewThreshold与初始亮温阈值iThreshold差值在50之内，则继续下一步骤；否则用iNewThreshold作为初始亮温阈值iThreshold，转到步骤②；

⑤根据新亮温阈值iNewThreshold对原始云图进行分割，即当云图中像素的亮温大于iNewThreshold时设定为1，否则设定为0，得到黑白二值图像。

3.如权利要求1所述的一种基于形状特征的云图检索方法，其特征在于去除云图的云系空洞和边界毛刺的具体步骤为：

①对阈值分割后的云图，先采用半径为5的平坦型圆盘结构元素进行膨胀处理，再通过半径为7的平坦型圆盘结构元素进行腐蚀处理，通过闭运算去除云系内部大部分细小空洞；

②对闭运算处理后的云图，先采用半径为6的平坦型圆盘结构元素进行腐蚀处理，再通过半径为7的平坦型圆盘结构元素进行膨胀处理，通过开运算去除边界毛刺。

4.如权利要求1所述的一种基于形状特征的云图检索方法，其特征在于采用几何不变矩对各子块进行区域特征描述的具体方法为：

①首先计算各云图子块的区域几何中心坐标

②计算云图各子块的中心矩μ_pq：式中p，q为正整数1或2，且满足p+q＝2或者p+q＝3，f(x,y)表示云图子块中坐标(x,y)处的像素值；

③对中心矩μ_pq进行规格化处理，式中为零阶中心矩的γ次方，γ为与p，q相关的指数变量，满足γ＝(p+q)/2；

④根据规格化后的中心矩η_pq，提取各子块的7个形状不变矩I_k，其中k＝1，2，3，4，5，6，7：I₁＝η₂₀+η₀₂，I₃＝(η₃₀-3η₁₂)²+(3η₂₁-η₀₃)²，

I₄＝(η₃₀+η₁₂)²+(η₂₁+η₀₃)²，

I₅＝(η₃₀-3η₁₂)(η₃₀+η₁₂)[(η₁₂+η₃₀)²-3(η₂₁+η₃₀)²]

+(3η₂₁-η₀₃)(η₂₁+η₀₃)[3(η₁₂+η₃₀)²-(η₂₁+η₃₀)²]，

I₆＝(η₂₀-η₀₂)[(η₁₂+η₃₀)²-(η₂₁+η₀₃)²]+4η₁₁(η₁₂+η₃₀)(η₂₁+η₀₃)，

I₇＝(3η₂₁-η₀₃)(η₃₀+η₁₂)[(η₁₂+η₃₀)²-3(η₂₁+η₀₃)²]

+(η₃₀-2η₁₂)(η₂₁+η₀₃)[3(η₁₂+η₃₀)²-(η₂₁+η₀₃)²]，

将所得的7个形状不变矩排为一列作为该云图子块的形状不变矩特征；

⑤当完成整幅云图所有子块形状不变矩特征的提取后，将所有子块的特征向量依次排为一长列向量作为该幅云图的形状特征向量，表示为向量K，K∈R^m。

5.如权利要求1所述的一种基于形状特征的云图检索方法，其特征在于采用稀疏分类方法进行云图检索相似度计算的具体步骤为：

①首先对待检索云图和云图特征库进行空间投影，采用正态分布随机矩阵R∈R^n×m(n<m)对云图特征库进行空间映射如下：同样对待检索云图特征也进行映射处理如下： 其中R是n行m列的矩阵；

②根据稀疏分类算法，计算残差，进行类别判断： 其中a_s表示稀疏系数中第s幅云图数据的对应系数；

③重复第①和②步骤，通过多次迭代，求取重构残差均值：其中l表示迭代次数，表示第t次迭代时属于s幅云图样本的残差值；

④对残差均值E[r_s]进行大小排序，获取与待检索云图的相似度信息，其中残差越小，相似性越大，最后返回残差最小的9幅云图。