基于张量动态纹理模型的极光视频分类方法
【技术领域】
[0001] 本发明属于图像处理技术领域，涉及四种形态极光视频的计算机分类方法，可用 于极光视频的特征提取和计算机图像识别。
【背景技术】
[0002] 极光是太阳风通过日侧极隙区注入到地球磁层时，沉降粒子沿磁力线与地球高层 大气相互作用而产生的绚丽光辉。极光是极区空间天气物理过程的观测窗口，直观地反映 了太阳风与地磁层的耦合过程，蕴含着大量日地空间的电磁活动信息，有着深刻的研宄意 义。
[0003] 中国北极黄河站的全天空成像系统（All-skyCamera)同时对极光的三个典 型谱段427. 8nm，557. 7nm和630.Onm进行连续观测，产生数以万计的极光图像，数据量 巨大。WangQ等人在文章"Spatialtexturebasedautomaticclassificationof daysideaurorainall-skyimages.JournalofAtmosphericandSolar-Terrestrial Physics，2010, 72(5) : 498-508. "中将极光按形态分为弧状、辐射状、帷幔状和热点状四类， 并得出了四种极光类型的统计分布规律。Pedersen等人发表的文章"PedersenTR，Gerken EA.Creationofvisibleartificialopticalemissionintheaurorabyhigh-power radiowaves.Nature，2005, 433 (7025) : 498-500"，胡泽俊等人发表文章"HuZJ，Yang H，HuangD，etal.Synopticdistributionofdaysideaurora:Multiple-wavelength all-skyobservationatYellowRiverStationinNy-Alesund,Svalbard.J.Atmos. Sol. -Terr.Phys.，2009, 71 (89) : 794-804" 以及文章"LorentzenDA，MoenJ，Oksavik K，etal.Insitumeasurementofanewlycreatedpolarcappatch.J.Geophys. Res.，2010, 115(A12). "，提供了大量的研宄材料，证明不同形态的极光对应着不同的磁层 边界层动力学过程。如何准确高效地对极光视频进行分类，既是揭示其磁层源区动力学过 程的一个关键，也是其发生机制研宄的重要环节，然而极光形态和动态变化复杂，无疑为极 地科研工作者带来了巨大的困难.
[0004] 计算机图像识别和分析技术的发展为海量极光数据分类研宄提供了可 能。2004 年，Syrjiisuo等人在文章"SyrjasuoM，PartamiesN.Numericimage featuresfordetectionofaurora[J] ?GeoscienceandRemoteSensing Letters,IEEE, 2012, 9(2) : 176-179. "中将计算机视觉的方法引入极光图像自动分类系 统中，该方法从分割后的极光区域中提取Fourier算子作为特征，通过最近邻方法实现 极光图像的自动分类，由于受到分割算法的影响，该方法仅对形状特征明显的弧状极光 图像分类效果良好；Wang等人于2007年在文章"WangQian,LiangJimin,HuZejun,Hu HaiHong，ZhaoHeng，HuHongQiao,GaoXinbo,YangHuigen.Spatialtexturebased automaticclassificationofdaysideaurorainall-skyimages.Journalof AtmosphericandSolar-TerrestrialPhysics，2010, 72 (5): 498 - 508. " 中使用主成分 分析法PCA对极光图像的灰度特征进行提取，提出了一种基于表象的极光分类方法，在冕 状极光分类研究方向取得了一定进展；2008年，Gao等人发表文章"L.Gao,X.B.Gao,and J.M.Liang.Daysidecoronaautoradetectionbasedonsampleselectionand adaBoostalgorithm.J.ImageGraph，2010, 15 (1) : 116-121. "，提出基于Gabor变换的 极光图像分类方法，采用了局部Gabor滤波器提取图像特征，在确保计算精度的情况下降 低了特征冗余信息，取得了较好的分类效果；2009年，Fu等人在文章"FuRu，JieLiand X.B.Gao. .Automaticauroraimagesclassificationalgorithmbasedonseparated texture.Proc.Int.Conf.RoboticsandBiomimetics, 2009:1331-1335. " 中将形态学成 分分析（MCA)与极光图像处理相结合，从经过MCA分离后所得到的极光纹理子图中提取 特征，用于弧冕两类极光图像的分类，提高了弧冕极光分类的正确率。后续的相关研宄 还有：Han等人在文章"BingHan,XiaojingZhao,DachengTao,etal.Daysideaurora classificationviaBIFs-basedsparserepresentationusingmanifoldlearning. InternationalJournalofComputerMathematics.Publishedonline:12Nov2013. " 中 又提出基于BIFs特征和C均值分类的极光分类法；Yang等人在文章"YangXi，LiJie，Han Bing,GaoXinbo.Wavelethierarchicalmodelforauroraimagesclassification. JournalofXidianUniversity, 2013, 40 (2): 18-24."中提出多层小波变换来表示极光图 像特征，取得了较高的分类准确率；2013年，Han等人在文章"HanB，YangC，GaoXB.Aurora imageclassificationbasedonLDAcombiningwithsaliencyinformation.Ruanjian XueBao/JournalofSoftware, 2013, 24(11) :2758-2766?"中引入隐含狄利克雷分布模型 LDA，并结合图像显著性信息，又进一步提高了极光图像的分类准确率.
[0005] 然而，上述这些已有的极光图像分析多是基于单幅图像和静态特征，关于极光 序列的自动分类的相关工作仍然还比较少。相关的进展主要有：Yang等人于2013年在 文章"YangQiuju.AuroralEventsDetectionandAnalysisBasedonASIandUVI Images[D].Xi'an:Xidianuniversity, 2013. "中提出基于隐马尔科夫模型的极光序列 分类理论，但是该方法本质仍是基于单幅图像特征；另外，Han在文章"HanB,LiaoQ,Gao XB.Spatial-Temporalpolewardvolumelocalbinarypatternsforaurorasequences eventdetection.RuanJianXueBao/JournalofSoftware, 2014, 25(9) :2172_2179?"中 构造了空时极向LBP表征算子提取极光序列的STP-LBP特征，用于检测极光视频中的极向 运动现象。但是，该算法只针对极向弧状极光视频的特征，并不具备普适性。由于目前没有 对于四类形态极光视频的普适性建模和提取视频动态特征的方法，计算机无法直接对四类 极光视频进行自动分类，只能进行单帧操作，或者单一性地对某一类极光视频进行分析，分 类准确率和分类效率都比较低。
【发明内容】
[0006] 本发明的目的在于针对上述已有技术的不足，提出一种基于张量动态纹理模型的 极光视频分类方法，以利用极光视频帧之间的重复相关性，实现对四类极光视频的计算机 自动分类，并提高分类准确率。
[0007] 为实现上述目的本发明的技术方案包括如下步骤：
[0008] 1)从已标记类别的极光视频数据库中任取N个视频组成训练集 {yi，y2,...，yk，...，yiJ，yk是第k个训练集样本，k= 1，2,...，N，将剩余极光视频组成测试 集，从测试集中取一个样本作为测试极光视频ytest;
[0009] 2)提取极光视频的混合观测矩阵C' ：
[0010] 2a)将测试极光视频ytest表示为y'⑴，/(〇eZ?/x/:，L为当前视频帧像素矩阵的 行数，12为当前视频帧像素矩阵的列数，t= 1，. . .，T，T是视频总帧数；
[0011] 2b)用观测到的极光视频帧y'（t)组成视频张量Y' ：
[0012] r= [..v'(l)，.v-'(2),....v-V)，... 0，(0]， K'eA5/|X/：xr；
[0013] 2c)令i= 1，2, 3,将极光视频张量Y'按第i维展开，得到展开矩阵Y'⑴，即 F'⑴e妒'r(2) ，r(3>eiT'其中e=I2XT，f=IiXT，g=IiXI2;
[0014] 2d)对展开矩阵Y'⑴进行SVD分解，将其分为三个矩阵，即Y'⑴=U⑴S(i)V(i)H，其 中U(i)GR为左基矩阵，S(i)为奇异值矩阵，V(i)为右基矩阵；
[0015] 当i= 1时，左基矩阵t/(1) 为行观测矩阵，对应的奇异值矩阵为，ef' 右基矩阵为F(1)eiTlXe;
[0016] 当i= 2时，左基矩阵t/(2> 为列观测矩阵，对应的奇异值矩阵为eF' 右基矩阵为F(21eiTx/;
[0017] 当i= 3时，左基矩阵t/(3>ei??5为张量动态纹理特征矩阵，对应的奇异值矩阵为 W