一种基于非负矩阵分解的手机图像特征提取方法
【专利摘要】本发明公开了一种基于非负矩阵分解的手机图像特征提取方法，特点是包括以下步骤：①获取用于表示手机任意一个主功能的界面图片，在该界面图片上截取多幅用于表示主功能的不同子功能的特征子图并进行预处理，得到对应的预处理后的尺寸大小一致的图像，然后将每幅预处理后的图像按列扫描形成一个M维的列向量，接着将所有预处理后的图像对应的列向量按序排列构成一个待分解矩阵V；②对V进行稀疏非负矩阵分解，得到一个基矩阵W和一个系数矩阵H；优点是通过对预处理后的图像进行稀疏非负矩阵分解，降低了数据计算量，分解得到的矩阵是具有合适的稀疏性的稀疏矩阵，因此在保留数据的主要特征的基础上能够减少储存空间并提高运算效率。
【专利说明】一种基于非负矩阵分解的手机图像特征提取方法
【技术领域】
[0001]本发明涉及一种图像特征提取方法，尤其是涉及一种基于非负矩阵分解的手机图像特征提取方法。
【背景技术】
[0002]随着科学技术的快速发展，非负矩阵分解已逐渐成为信号处理、生物医学工程、模式识别、计算机视觉和图像工程等研究领域中最受欢迎的多维数据处理工具之一。非负矩阵分解使得分解后的所有分量均为非负值，并且同时实现了非线性的维数约减。传统的非负矩阵分解算法包括LNMF法、CNMF法、NSC法和SNMF法。然而，这些传统的非负矩阵分解算法并不适用于手机图像的处理，这是因为手机图像数据量大，各种应用图标特征比较特另O，需要能保留图像的大致特征、得出来的数据具有可观性且效率高的算法。因此，需要研究一种基于非负矩阵分解算法的手机图像处理方法。
【发明内容】
[0003]本发明所要解决的技术问题是提供一种运算效率高、存储空间小且对大规模数据处理速度快的基于非负矩阵分解的手机图像特征提取方法。
[0004]本发明解决上述技术问题所采用的技术方案为:一种基于非负矩阵分解的手机图像特征提取方法，包括以下步骤:
[0005]①获取用于表示手机任意一个主功能的界面图片，在该界面图片上截取多幅用于表示主功能的不同子功能的特征子图，对这些特征子图进行预处理，得到对应的预处理后的图像，且这些特征子图各自对应的预处理后的图像的尺寸大小一致，然后将每幅预处理后的图像按列扫描形成一个M·维的列向量，接着将所有预处理后的图像对应的列向量按序排列构成一个待分解矩阵，记为V，V=[Vl, - ,Vi,…，vN]，其中，M=W’ XH’，W’表示预处理后的图像的宽度，H’表示预处理后的图像的高度，V的维数为MXN，Vl表示第I幅预处理后的图像对应的列向量，Vi表示第i幅预处理后的图像对应的列向量，vN表示第N幅预处理后的图像对应的列向量，N，N表示截取的特征子图的总幅数；
[0006]②对V进行稀疏非负矩阵分解，得到一个MX R维的包含图像特征的基矩阵W和一
M X N
个RXN维的系数矩阵H，V=WXH,其中，O < Λ* < ——。
M + JN
[0007]所述的步骤②的具体过程为:
[0008]②-1、令W表示V稀疏非负矩阵分解后得到的包含图像特征的基矩阵，令H表示V稀疏非负矩阵分解后得到的系数矩阵，并令Sw表示W的稀疏度，令Sh表示H的稀疏度；
[0009]②-2、给定Sw和Sh的初始值，然后判断Sw的初始值是否为0，如果是，则将W初始化为一个稀疏度任意的随机矩阵，否则，将W初始化为一个稀疏度为给定的初始值的随机矩阵；同样判断Sh的初始值是否为0，如果是，则将H初始化为一个稀疏度任意的随机矩阵，否则，将H初始化为一个稀疏度为给定的初始值的随机矩阵；[0010]②-3、根据步骤②-2确定的W和H，计算初始目标函数值，记为Etjld,Eold= V-WHl I2，其中，符号“I I 112”为2-范数符号；
[0011]②-4、判断W的稀疏度是否为0，如果是，则根据W’ =WX (VHt) / (WHHt)更新W，得到更新后的基矩阵W’，然后执行步骤②-5 ;否则，根据W’=W-Uw(WH-V)Ht更新W，得到更新后的基矩阵W’，然后对W’进行映射，得到映射后的基矩阵，记为W’ ’，接着令W’ =W’ ’，再执行步骤②-5 ;其中，W’=W’’中的“=”为赋值符号，HtSH的转置矩阵，μ w表示步长，μ?>0且μ w是一个较小的常数；
[0012]②-5、判断H的稀疏度是否为0，如果是，则根据H’ =HX (WtV) / (WtWH)更新H，得到更新后的系数矩阵H’，然后执行步骤②_6 ;否则，根据H’=H-μ hWt(WH-V)更新H，得到更新后的系数矩阵H’，然后对H’进行映射，得到映射后的系数矩阵，记为H’’，接着令H’ =Η’’，再执行步骤②_6 ;其中，H’ =Η’ ’中的“=”为赋值符号，Wt为W的转置矩阵，μ η表示步长，μ Η>0且μ η是一个较小的常数；
[0013]②-6、令W=W’，并令Η=Η’，然后计算新目标函数值，记为Enew，Enew= I V-WH| |2，其中，W=W'和H=H’中的“=”为赋值符号；
[0014]②-7、计算Enew与Etjld的差值，记为ε，ε =Enew-Eold,然后判断ε是否达到预先设置的表示无穷小的值，如果达到，则结束，得到最终的MXR维的包含图像特征的基矩阵W和最终的RXN维的系数矩阵H ，否则，返回步骤②-4继续执行。
[0015]所述的步骤②_4中μ w的取值范围为0.01~10 ;所述的步骤②_5中μ H的取值范围为0.01~10。
[0016]所述的步骤②-4中对W’进行映射的具体过程为:
[0017]②-4-1、分别获取与W’中的每列列向量欧式距离最近的列向量:设置一个初始列向量 SpSi=Wi' +(λ「Σ Wi')/dim (Wi'),其中，I ≤ i ≤ R, λ 工=| w/ I1, Wi' | I1 表示 Wi'的1-范数，即 λ J= Σ I Wi' I, dim (w/)表示 Wi'的维数，即 dim (w/) =M ;
[0018]②-4-2、设置一个初始状态为空集的集合Z，并设置一个R维的中间列向量m，将该中间列向量m中的第j个元素记为Hij,其中，I≤j≤R ;
[0019]②-4-3、令所,=j。1、’如果/eZ，其中，dim(Z)表示集合Z的维数；
[0020]②-4-4、根据s/=Iiij+a (Wi’-mj)更新Si，得到更新后的列向fiSi’，其中，α为迭
代系数，且 <2 = (-/rhrcal(V/)3 -Aac))12a , real 为取实部函数，a= Σ (stones (M, DXA1/
(M-dim(Z)))2，ones (Μ, I)为 MX I 阶元素全部为 I 的矩阵，b=2X [stones (M，I) X λ /(M-dim(Z)) ]TX Si, c= Σ sit2- λ 2，sit 为 Si 中的第 t 个元素，I ^ t ^ Μ, λ 2 表示 w/ 的 2-范数；然后令Si=Si'，其中，Si=Si'中的“=”为赋值符号；
[0021]②-4-5、判断更新后的Si中的所有元素是否均为非负，如果是，则将更新得到的Si作为需要求得的W’ ’的第i列列向量，并结束；否则，将更新后的Si中所有小于O的元素中的第h个元素记为sih,并令Z=Z U {h:sih〈0},且令sih=0,然后执行步骤②-4-6,其中，Z=Z U {j:SiJ<0}中的“=”为赋值符号，“U”为并集运算符号，I≤h≤H，H表示Si中所有小于O的元素的总个数，{h:sih<0}表示Si中所有小于O的元素的集合；[0022]②-4-6、将更新后的Si中所有大于O的元素中的第k个元素记为sik，并根据sik’ =Sik-S'更新sik，其中，1≤k≤K，K表示Si中所有大于O的元素的总个数，s' = ( Σ sik- A1)/ (dim(w/ ) -dim(Z));然后令 sik=sik’，其中，sik=sik’ 中的“=”为赋值符号;再返回步骤②-4-3继续迭代。
[0023]所述的步骤②-5中对H’进行映射的具体过程与所述的步骤②_4中对W’进行映射的具体过程相同。该映射保证了能够得到响应稀疏度的最近向量，使得求解目标的稀疏度得到精确地控制，且计算过程简单易实现，提升了算法结果的灵活性，可根据需要进行合
理设置。
[0024]所述的步骤②_7中预先设置的表示无穷小的值为1 X 10_7。
[0025]与现有技术相比，本发明的优点在于:对于手机任意一个主功能的界面图片，在每张界面图片内截取多幅用于表示主功能的不同子功能的特征子图，对这些特征子图进行预处理，得到对应的预处理后的图像，再对预处理后的图像进行稀疏非负矩阵分解，降低了数据计算量，分解得到的结果比较直观，而且采用稀疏非负矩阵分解的方法使得分解得到的矩阵是具有合适的稀疏性的稀疏矩阵，因此在保留数据的主要特征的基础上能够减少储存空间、提高运算效率，并且有利于对大规模数据进行快速处理。
【专利附图】
【附图说明】
[0026]图1为本发明的处理过程流程图。
【具体实施方式】
[0027]以下结合附图实施例对本发明作进一步详细描述。
[0028]一种基于非负矩阵分解的手机图像特征提取方法，包括以下步骤:
[0029]①获取用于表示手机任意一个主功能的界面图片，在该界面图片上截取多幅用于表示主功能的不同子功能的特征子图，对这些特征子图进行预处理，得到对应的预处理后的图像，且这些特征子图各自对应的预处理后的图像的尺寸大小一致，然后将每幅预处理后的图像按列扫描形成一个M维的列向量，接着将所有预处理后的图像对应的列向量按序排列构成一个待分解矩阵，记为V，V=[Vl, - ,Vi,…，vN]，其中，M=W’ XH’，W’表示预处理后的图像的宽度，H’表示预处理后的图像的高度，V的维数为MXN，Vl表示第I幅预处理后的图像对应的列向量，Vi表示第i幅预处理后的图像对应的列向量，vN表示第N幅预处理后的图像对应的列向量，1≤i≤N，N表示截取的特征子图的总幅数。
[0030]②对V进行稀疏非负矩阵分解，得到一个MX R维的包含图像特征的基矩阵W和一
个RXN维的系数矩阵H，V=WXH,其中，
【权利要求】
1.一种基于非负矩阵分解的手机图像特征提取方法，其特征在于包括以下步骤: ①获取用于表示手机任意一个主功能的界面图片，在该界面图片上截取多幅用于表示主功能的不同子功能的特征子图，对这些特征子图进行预处理，得到对应的预处理后的图像，且这些特征子图各自对应的预处理后的图像的尺寸大小一致，然后将每幅预处理后的图像按列扫描形成一个M维的列向量，接着将所有预处理后的图像对应的列向量按序排列构成一个待分解矩阵，记为V，V=[Vl, - ,Vi,…，vN]，其中，M=W' XH’，W’表示预处理后的图像的宽度，H’表示预处理后的图像的高度，V的维数为MXN，Vl表示第I幅预处理后的图像对应的列向量，Vi表示第i幅预处理后的图像对应的列向量，vN表示第N幅预处理后的图像对应的列向量，I ^ i ^ N, N表示截取的特征子图的总幅数； ②对V进行稀疏非负矩阵分解，得到一个MXR维的包含图像特征的基矩阵W和一个
MxNRXN维的系数矩阵H，V=WXH，其中，0<R< ^ 。
M + N
2.根据权利要求1所述的一种基于非负矩阵分解的手机图像特征提取方法，其特征在于所述的步骤②的具体过程为: ②-1、令W表示V稀疏非负矩阵分解后得到的包含图像特征的基矩阵，令H表示V稀疏非负矩阵分解后得到的系数矩阵，并令Sw表示W的稀疏度，令Sh表示H的稀疏度； ②_2、给定Sw和Sh的初始值，然后判断Sw的初始值是否为O，如果是，则将W初始化为一个稀疏度任意的随机矩阵，否则，将W初始化为一个稀疏度为给定的初始值的随机矩阵；同样判断Sh的初始值是否为O，`如果是，则将H初始化为一个稀疏度任意的随机矩阵，否则，将H初始化为一个稀疏度为给定的初始值的随机矩阵； ②-3、根据步骤②-2确定的W和H，计算初始目标函数值，记为Etjld, Etjld= I V-WHl |2，其中，符号“II M2 ”为2-范数符号； ②_4、判断W的稀疏度是否为O，如果是，则根据W’ =WX (VHt) / (WHHt)更新W，得到更新后的基矩阵W’，然后执行步骤②-5 ;否则，根据W’=W-Uw(WH-V) Ht更新W，得到更新后的基矩阵W’，然后对W’进行映射，得到映射后的基矩阵，记为W’ ’，接着令W’ =W’ ’，再执行步骤②-5 ;其中，W’=W’’中的“=”为赋值符号，HtSH的转置矩阵，μ?表示步长，是一个较小的常数； ②-5、判断H的稀疏度是否为O，如果是，则根据H’ =HX (WtV) / (WtWH)更新H，得到更新后的系数矩阵H’，然后执行步骤②-6 ;否则，根据H’ =H- μ HffT (WH-V)更新H，得到更新后的系数矩阵H’，然后对H’进行映射，得到映射后的系数矩阵，记为H’’，接着令H’ =H’’，再执行步骤②-6;其中，H’=H’’中的“=”为赋值符号，Wt为W的转置矩阵，μ H表示步长，μΗ>0且μ η是一个较小的常数； ②-6、令W=W’，并令Η=Η’，然后计算新目标函数值，记为En?，Enew= I V-WHl |2，其中，W=W’和H=H’中的“=”为赋值符号； ②-7、计算Enew与Etjld的差值，记为ε，ε =Enew-Etjld，然后判断ε是否达到预先设置的表示无穷小的值，如果达到，则结束，得到最终的MXR维的包含图像特征的基矩阵W和最终的RXN维的系数矩阵H，否则，返回步骤②-4继续执行。
3.根据权利要求2所述的一种基于非负矩阵分解的手机图像特征提取方法，其特征在于所述的步骤②-4中μ w的取值范围为0.01~10 ;所述的步骤②-5中μ H的取值范围为.0.01 ~10。
4.根据权利要求2或3所述的一种基于非负矩阵分解的手机图像特征提取方法，其特征在于所述的步骤②-4中对W’进行映射的具体过程为: ②-4-1、分别获取与W’中的每列列向量欧式距离最近的列向量:设置一个初始列向量
5.根据权利要求4所述的一种基于非负矩阵分解的手机图像特征提取方法，其特征在于所述的步骤②-7中预先设置的表示无穷小的值为1X10'
【文档编号】G06T7/00GK103679715SQ201310651941
【公开日】2014年3月26日 申请日期:2013年12月5日 优先权日:2013年12月5日
【发明者】吴月, 叶庆卫, 周宇, 王晓东 申请人:宁波大学