一种空间数据访问规律动态统计及信息压缩的方法及系统的制作方法
【专利摘要】本发明针对空间数据访问规律统计信息具有的离散性、无结构性，设计了一种空间数据访问规律动态统计及信息压缩的方法及系统，通过空间数据访问规律统计信息的点云表达、聚类梯度计算与聚类区域提取，实现偶发性访问的无效空间数据访问规律统计信息剔除，减少信息基量；通过对每个聚类点云集区域，按照点云差数和访问次数的点云序列用1个字节编码形成最终压缩统计数据输出，减少信息总量。本发明相对于现有技术，解决了由于全局动态统计带来的不可承载的网络流量问题，实现了实时动态跟踪空间数据访问全局分布规律，大大改善了空间数据系统公众服务能力。
【专利说明】一种空间数据访问规律动态统计及信息压缩的方法及系统
【技术领域】
[0001]本发明属于大规模分布式环境下空间地理信息系统【技术领域】，尤其涉及海量空间信息存储、访问、缓存和控制应用【技术领域】，具体涉及一种云计算环境下空间数据访问规律动态统计及信息压缩的方法及系统。
【背景技术】
[0002]随着高精度遥感影像技术的提高，对空间数据(D0M/DEM/DLG/DRG)的开发和利用，已经被广泛重视并应用到众多领域。数字地球中的空间信息服务质量问题，即如何支持大规模、高强度的并发用户访问以提供高可靠、高可用的空间信息服务以及如何最快速地将空间数据传输并显示在数字地球的客户端，成为工业界和学术界的关注热点。随着互联网络技术的发展和宽带网络的普及，空间数据的服务形式从桌面转变到分布式网络环境中，新的服务形式迫切需要将空间数据服务和网络通信技术融合到一起，从而在分布式网络环境下为用户各种空间数据浏览、分析、查询、漫游等服务提供高效的平台支持。空间信息服务的质量，对于加速推广和普及空间信息服务和提升数字地球用户的服务体验，具有重要意义，相关的研究也就成为热点。
[0003]空间数据访问规律的动态统计，可有效的为空间数据的存储组织、副本数量及其存储分布提供依据，而动态统计所产生的海量空间数据访问规律统计信息将带来不可承载的网络流量，既而有算法研究如何通过P2P的方式，减少统计信息融合的副本数量。但由于统计信息基本量巨大，简单的代理分发和融合依然难以满足要求。
[0004]研究表明，空间用户访问统计信息存在由于用户兴趣、地区热点等引发的用户访问特征行为，也存在用户随意访问等偶发性行为。针对用户特征访问行为，其对应瓦片被再次访问的可能性大，而对偶发性访问行为，其行为准则难以判断，对应瓦片被再次访问的可能性不可预测。针对偶发性访问而进行空间数据瓦片的存储组织调整，将变得非常频繁和低效，不但很难有效提高空间信息服务系统的性能，相反，这种频繁的无序调整还将额外增加系统的计算和传输开销因此，如何减少偶发性访问形成的空间数据访问规律信息，成为问题的关键之一。
[0005]针对空间数据的压缩有很多研究，包括删格数据的压缩、矢量数据的压缩等。将空间数据映射成空间点云进行数据压缩处理的研究相对较少。纯粹的基于点云处理的研究包括边界提取、表面重建以及模型简化处理等，而点云去噪是基于点云测量和处理的重要过程之一，应用领域较广，包括地形测量，如利用地面三维激光扫描进行地形测量数据去噪，三维重构和曲面拟合等。
[0006]不同应用目的，点去噪算法也存在差异，对于有序或者部分有序的点云来说，可以采用的去噪算法有很多，如基于最小二乘估计的维纳滤波、卡尔曼滤波等。针对无序点云，一般是先利用相关规则建立其相互之间的拓扑关系，如八叉树法、k-d tree法以及空间单元格法等。海量空间统计数据点云规模巨大,数据映射到点云后,其点云基本是无序的，相互之间没有拓扑关系或很难确定其拓扑关系，针对这种无序海量点云，目前的研究算法较少。
[0007]所以，解决由于海量空间数据访问规律统计信息传输带来的不可承载的网络流量问题，从而实现能实时动态跟踪空间数据访问全局分布规律、大大改善空间数据系统公众服务能力的方法及系统，是目前该领域的重要研究方向。
【发明内容】
[0008]本发明目的在于解决现有技术不足，提供一种能解决由于海量空间数据访问规律统计信息传输带来的不可承载的网络流量问题、能实时动态跟踪空间数据访问全局分布规律、大大改善空间数据系统公众服务能力的空间数据访问规律动态统计及信息压缩的方法及系统。
[0009]本发明的方法所采用的技术方案是:一种空间数据访问规律动态统计及信息压缩的方法，其特征在于，包括以下步骤:
[0010]步骤1:空间数据访问规律统计信息点云输出:将空间数据访问规律信息按照每块瓦片数据的空间位置及其被访问下载次数，形成空间数据访问规律统计信息点云输出；
[0011]步骤2:空间数据访问规律统计信息聚类计算:根据获取的聚类计算半径参数计算空间数据访问规律统计信息点云聚类梯度，根据计算得到的聚类梯度剔除偶发性的离散访问点云，保留有效特征访问点云，减少空间数据访问规律统计信息基量；
[0012]步骤3:空间数据访问规律统计信息区域提取:根据获取的聚类区域提取扩散半径、聚类区域内点云差数量参数计算聚类内点云欧式距离，将特征访问点云进行聚类分组，形成聚类点云集，对空间数据访问规律统计信息进行聚类区域划分和分类提取；
[0013]步骤4:空间数据访问规律统计信息聚类压缩:根据每个聚类点云集区域，计算聚类区域内点云差，按照点云差数和访问次数的点云序列用I个字节编码形成最终压缩统计数据输出。
[0014]作为优选，所述的空间数据按照金子塔瓦片的形式存储，每块瓦片的被访问次数为空间数据访问规律统计信息；空间数据访问规律统计信息包括空间数据的位置坐标和访问次数。
[0015]作为优选，所述的空间数据访问规律统计信息空间坐标范围、空间数据块数量参数通过系统计算自动获取；聚类梯度计算参数、聚类区域提取扩散半径、聚类区域内点云差数量通过配置参数设置。
[0016]本发明的系统所采用的技术方案是:一种空间数据访问规律动态统计及信息压缩的系统，其特征在于，包括:空间数据访问规律统计信息表达单元、空间数据访问规律统计信息聚类计算单元、空间数据访问规律统计信息区域提取单元和空间数据访问规律统计信息聚类压缩单元；
[0017]所述的空间数据访问规律统计信息表达单元用于将空间数据访问规律信息按照每块瓦片数据的空间位置及其被访问下载次数，形成空间数据访问规律统计信息点云表达；
[0018]所述的空间数据访问规律统计信息聚类计算单元用于根据获取的聚类计算半径参数计算空间数据访问规律统计信息点云聚类梯度，剔除偶发性的离散访问点云，保留有效特征访问点云，减少空间数据访问规律统计信息基量；[0019]所述的空间数据访问规律统计信息区域提取单元用于根据获取的聚类区域提取扩散半径、聚类区域内点云差数量参数计算聚类内点云欧式距离，将特征访问点云进行聚类分组，形成聚类点云集，达到对空间数据访问规律统计信息进行聚类区域划分和提取；
[0020]所述的空间数据访问规律统计信息聚类压缩单元用于对提取的每个聚类点云集区域，按照点云差数和访问次数的点云序列用I个字节编码形成最终压缩统计数据输出。
[0021]作为优选，所述的空间数据访问规律统计信息表达单元包括空间参数获取模块、点云表达调节因子计算模块和空间点云输出模块；
[0022]所述的空间参数获取模块用于获取空间数据访问规律统计信息空间坐标范围、空间数据块数量、空间数据被访问或下载最大次数参数；
[0023]所述的点云表达调节因子计算模块用于计算空间数据访问规律统计信息中空间数据点云输出调节因子；
[0024]所述的空间点云输出模块用于根据空间数据位置和调节因子输出将空间数据访问规律输出成空间点云。
[0025]作为优选，所述的空间数据访问规律统计信息聚类计算单元包括:聚类梯度计算参数输入模块和聚类梯度计算模块；
[0026]所述的聚类梯度计算参数输入模块用于控制聚类计算半径；
[0027]所述的聚类梯度计算模块用于区域内聚类梯度计算。
[0028]作为优选，所述的空间数据访问规律统计信息区域提取单元包括聚类区域提取参数输入模块和聚类区域提取模块；
[0029]所述的聚类区域提取参数输入模块用于限制聚类区域提取扩散半径、限制聚类区域内点云差数量；
[0030]所述的聚类区域提取模块用于提取聚类区域。
[0031]作为优选，所述的空间数据访问规律统计信息聚类压缩单元包括聚类区域点云差计算模块、点云序列组织模块和点云序列游程编码压缩模块；
[0032]所述的聚类区域点云差计算模块用于计算聚类区域内点云差；
[0033]所述的点云序列组织模块用于将聚类区域按照坐标顺序和点云差数量依序组织；
[0034]所述的点云序列游程编码压缩模块用于按照游程编码规则，将按序组织的点云序列进行编码压缩。
[0035]本发明具有的有益效果是:空间数据访问规律统计信息存在由于用户兴趣、地区热点等引发的用户访问特征行为，也存在用户随意访问等偶发性行为；针对用户特征访问行为，其对应空间数据瓦片被再次访问的可能性大，而对偶发性访问行为，其行为准则难以判断，对应瓦片被再次访问的可能性不可预测；本发明能在空间数据访问规律统计信息中，剔除偶发性访问形成的统计信息，并对点云数据根据欧式距离形成聚类集后压缩编码，形成最终统计信息输出，大大减少实际需要传输分发的空间数据访问规律统计信息总量；解决了由于海量空间数据访问规律统计信息传输带来的不可承载的网络流量问题，能实时动态跟踪空间数据访问全局分布规律，大大改善空间数据系统公众服务能力。
【专利附图】
【附图说明】[0036]附图1:是本发明的方法流程图。
[0037]附图2:是本发明实施例中系统的结构示意图。
[0038]附图3:是本发明实施例中空间数据访问规律统计信息表达单元的功能模块图。
[0039]附图4:是本发明实施例中空间数据访问规律统计信息聚类计算单元的功能模块图。
[0040]附图5:是本发明实施例中空间数据访问规律统计信息区域提取单元的功能模块图。
[0041]附图6:是本发明实施例中空间数据访问规律统计信息聚类压缩单元的功能模块图。
【具体实施方式】
[0042]以下结合附图和【具体实施方式】对本发明做进一步的阐述，本实施例主要以地理信息系统作为载体进行具体说明。
[0043]针对空间数据访问分布规律，微软研究院通过Hotmap模型得到数字地球系统中的影像瓦片请求符合幂律(Power-1aw)访问分布,但这种局部、静态的访问分布规律,不能真正实时、动态的反映访问分布规律的全局信息，从而制约了系统服务能力的提高。一方面，由于存在地形区域的差异(海洋、沙漠、城市)、时事热点的差异(海湾战争、朝韩冲突、台海问题)等，针对空间数据的访问是高度分布不均的，而且国际时事风云变幻，热点区域经常转移，导致空间数据的访问规律也动态变化；另一方面，分布式空间信息服务环境下的异构服务器也存在更新替换、性能老化，导致这种分布不均和设备性能差异(处理能力的差异(内存、CPU)、存储能力的差异(存储聚合速度、存储容量、容错能力)、有效传输带宽的差异)是经常变化的；因此需要对这些要素信息进行动态统计和分发，以便跟踪这些变化，及时调整地形数据的存储组织策略和预取、缓存策略。
[0044]通过简单的同步、复制机制进行异构服务器的全局、动态统计和获取空间数据访问分布规律，将带来的不可承载的网络流量。以SRTM90地形数据为例，共3538890块瓦片，若每个字节表示一个瓦片的访问次数，则每个字节可以代表最大255次访问，可以计算，当服务器数量达到1000个时，一次统计分发，由于没有信息压缩，每个服务器产生的数据都将分发999次，总共将分发3292.55GB的总数据量，通过这种方式实时、动态统计将持续产生大量的统计流量，从而挤占用户服务带宽，影响服务体验。因此，如何有效的减少空间数据访问规律动态统计信息总量成为关键。
[0045]请见图1、图2、图4、图5和图6，本发明的方法所采用的技术方案是:一种空间数据访问规律动态统计及信息压缩的方法，包括以下步骤:
[0046]步骤1:空间数据访问规律统计信息点云输出:将空间数据访问规律信息按照每块瓦片数据的空间位置及其被访问下载次数，形成空间数据访问规律统计信息点云输出；
[0047]步骤2:空间数据访问规律统计信息聚类计算:根据获取的聚类计算半径参数计算空间数据访问规律统计信息点云聚类梯度，根据计算得到的聚类梯度剔除偶发性的离散访问点云，保留有效特征访问点云，减少空间数据访问规律统计信息基量；
[0048]步骤3:空间数据访问规律统计信息区域提取:根据获取的聚类区域提取扩散半径、聚类区域内点云差数量参数计算聚类内点云欧式距离，将特征访问点云进行聚类分组，形成聚类点云集，对空间数据访问规律统计信息进行聚类区域划分和分类提取；
[0049]步骤4:空间数据访问规律统计信息聚类压缩:根据每个聚类点云集区域，计算聚类区域内点云差，按照点云差数和访问次数的点云序列用I个字节编码形成最终压缩统计数据输出。
[0050]其中，空间数据按照金子塔瓦片的形式存储，每块瓦片的被访问次数为空间数据访问规律统计信息；空间数据访问规律统计信息包括空间数据的位置坐标和访问次数；空间数据访问规律统计信息空间坐标范围、空间数据块数量参数通过系统计算自动获取；聚类梯度计算参数、聚类区域提取扩散半径、聚类区域内点云差数量通过配置参数设置。
[0051]请见图2、本发明的系统所采用的技术方案是:一种空间数据访问规律动态统计及信息压缩的系统，包括:空间数据访问规律统计信息表达单元100、空间数据访问规律统计信息聚类计算单元200、空间数据访问规律统计信息区域提取单元300和空间数据访问规律统计信息聚类压缩单元400 ;空间数据访问规律统计信息表达单元100用于将空间数据访问规律信息按照每块瓦片数据的空间位置及其被访问下载次数，形成空间数据访问规律统计信息点云表达；空间数据访问规律统计信息聚类计算单元200用于根据获取的聚类计算半径参数计算空间数据访问规律统计信息点云聚类梯度，剔除偶发性的离散访问点云，保留有效特征访问点云，减少空间数据访问规律统计信息基量；空间数据访问规律统计信息区域提取单元300用于根据获取的聚类区域提取扩散半径、聚类区域内点云差数量参数计算聚类内点云欧式距离，将特征访问点云进行聚类分组，形成聚类点云集，达到对空间数据访问规律统计信息进行聚类区域划分和提取；空间数据访问规律统计信息聚类压缩单元400用于对提取的每个聚类点云集区域，按照点云差数和访问次数的点云序列用I个字节编码形成最终压缩统计数据输出。
[0052]空间数据访问规律统计信息表达单元100包括空间参数获取模块101、点云表达调节因子计算模块102和空间点云输出模块103 ;空间参数获取模块101用于获取空间数据访问规律统计信息空间坐标范围、空间数据块数量、空间数据被访问或下载最大次数参数；点云表达调节因子计算模块102用于计算空间数据访问规律统计信息中空间数据点云输出调节因子；
[0053]空间点云输出模块103用于根据空间数据位置和调节因子输出将空间数据访问规律输出成空间点云。
[0054]其中，空间数据访问规律统计信息聚类计算单元200包括:聚类梯度计算参数输入模块201和聚类梯度计算模块202 ;聚类梯度计算参数输入模块201用于控制聚类计算半径；聚类梯度计算模块202用于区域内聚类梯度计算。空间数据访问规律统计信息区域提取单元300包括聚类区域提取参数输入模块301和聚类区域提取模块302 ;聚类区域提取参数输入模块301用于限制聚类区域提取扩散半径、限制聚类区域内点云差数量；聚类区域提取模块302用于提取聚类区域。空间数据访问规律统计信息聚类压缩单元400包括聚类区域点云差计算模块401、点云序列组织模块402和点云序列游程编码压缩模块403 ；聚类区域点云差计算模块401用于计算聚类区域内点云差；点云序列组织模块402用于将聚类区域按照坐标顺序和点云差数量依序组织；点云序列游程编码压缩模块403用于按照游程编码规则，将按序组织的点云序列进行编码压缩。
[0055]本系统实际工作时的具体流程为:[0056](I)空间数据访问规律统计信息点云输出:通过所述空间数据访问规律统计信息表达单元100将空间数据访问规律信息按照每块瓦片数据的空间位置及其被访问下载次数，形成空间数据访问规律统计信息点云输出；在这一过程，包括以下方面:
[0057]①利用所述的空间参数获取模块101获取空间数据访问规律统计信息空间坐标范围、空间数据块数量、空间数据被访问或下载最大次数参数；
[0058]空间数据访问规律统计信息的空间坐标范围根据所有统计的空间数据坐标，分别选择最大坐标和最小坐标作为空间坐标范围；
[0059]空间数据块数量N依据空间坐标范围计算面积可得，即连续空间坐标范围内，每个点代表一块空间数据；
[0060]空间数据被访问或下载最大次数参数只需要依次获取每块空间数据块i的被访问或下载次数ki，取最大值即可。
[0061]②利用所述的点云表达调节因子计算模块102计算空间数据访问规律统计信息中空间数据点云输出调节因子；
[0062]若空间瓦片总个数为N，每个瓦片的访问次数分别为L则调节因子α要求满足式(I)限制；
[0063]a ≤ 255-max(kj), i=l, 2, — ,N (I)
[0064]其中255表示最大访问次数必须不大于255，对于大于255的空间数据，通过拆分成多份数据进行压缩处理，最大访问次数必须小于255的限制，是为了满足后续处理中游程编码运算效率要求。
[0065]③利用所述的空间点云输出模块103将空间数据访问规律信息根据空间数据位置和调节因子输出成空间点云；
[0066]首先将地理信息系统的每个瓦片看成一个点，所有瓦片形成一个空间连续点集，没有被访问的空间瓦片称为“虚点”，被访问的瓦片称“实点”；
[0067]同时，“虚点”以颜色RGB (255，255，255)表示，实点i根据其访次数，分别按照颜色RGB (ki+α，0，0)表示，其中ki为空间数据实点i被访问或下载次数，最大不超过255次；
[0068]最后，按照“虚点”、“实点”的相对坐标，将空间数据访问映射成空间点云，其中“实点”代表被访问的瓦片，点云颜色深浅表示被访问的次数。
[0069](2)空间数据访问规律统计信息聚类计算:通过所述空间数据访问规律统计信息聚类计算单元200，根据获取的聚类计算半径参数计算空间数据访问规律统计信息点云聚类梯度,根据计算得到的聚类梯度剔除偶发性的离散访问点云,保留有效特征访问点云,减
少空间数据访问规律统计信息基量；在这一过程，包括以下方面:
[0070]①利用所述的聚类梯度计算参数输入模块201获取用于控制聚类计算半径参数；
[0071]控制聚类计算半径参数r和空间访问聚类梯度阀值T通过配置管理输入。
[0072]②利用所述的聚类梯度计算模块202进行区域内聚类梯度计算；
[0073]考察一段时间t，设某瓦片“实点”Pi的访问次数为ki;以Pi为中心，边长r的区域内，其ni个邻居“实点”瓦片p' i；J(j=l, 2,..., Iii)的访问次数分别为k'
访问聚类梯度表示某“实点”瓦片被访问的聚集程度，聚集程度越低，是偶发性访问的可能性越大，定义瓦片Pi的访问聚类梯度为T(Pi)，则有:
【权利要求】
1.一种空间数据访问规律动态统计及信息压缩的方法，其特征在于，包括以下步骤: 步骤1:空间数据访问规律统计信息点云输出:将空间数据访问规律信息按照每块瓦片数据的空间位置及其被访问下载次数，形成空间数据访问规律统计信息点云输出； 步骤2:空间数据访问规律统计信息聚类计算:根据获取的聚类计算半径参数计算空间数据访问规律统计信息点云聚类梯度，根据计算得到的聚类梯度剔除偶发性的离散访问点云,保留有效特征访问点云,减少空间数据访问规律统计信息基量； 步骤3:空间数据访问规律统计信息区域提取:根据获取的聚类区域提取扩散半径、聚类区域内点云差数量参数计算聚类内点云欧式距离，将特征访问点云进行聚类分组，形成聚类点云集，对空间数据访问规律统计信息进行聚类区域划分和分类提取； 步骤4:空间数据访问规律统计信息聚类压缩:根据每个聚类点云集区域，计算聚类区域内点云差，按照点云差数和访问次数的点云序列用I个字节编码形成最终压缩统计数据输出。
2.根据权利要求1所述的空间数据访问规律动态统计及信息压缩的方法，其特征在于:所述的空间数据按照金子塔瓦片的形式存储，每块瓦片的被访问次数为空间数据访问规律统计信息；空间数据访问规律统计信息包括空间数据的位置坐标和访问次数。
3.根据权利要求1所述的空间数据访问规律动态统计及信息压缩的方法，其特征在于:所述的空间数据访问规律统计信息空间坐标范围、空间数据块数量参数通过系统计算自动获取；聚类梯度计算参数、聚类区域提取扩散半径、聚类区域内点云差数量通过配置参数设置。
4.一种空间数据访问规律动态统计及信息压缩的系统，其特征在于，包括:空间数据访问规律统计信息表达单元(100)、空间数据访问规律统计信息聚类计算单元(200)、空间数据访问规律统计信息区域提取单元(300)和空间数据访问规律统计信息聚类压缩单元(400)； 所述的空间数据访问规律统计信息表达单元(100)用于将空间数据访问规律信息按照每块瓦片数据的空间位置及其被访问下载次数，形成空间数据访问规律统计信息点云表达； 所述的空间数据访问规律统计信息聚类计算单元(200)用于根据获取的聚类计算半径参数计算空间数据访问规律统计信息点云聚类梯度，剔除偶发性的离散访问点云，保留有效特征访问点云，减少空间数据访问规律统计信息基量； 所述的空间数据访问规律统计信息区域提取单元(300)用于根据获取的聚类区域提取扩散半径、聚类区域内点云差数量参数计算聚类内点云欧式距离，将特征访问点云进行聚类分组，形成聚类点云集，达到对空间数据访问规律统计信息进行聚类区域划分和提取； 所述的空间数据访问规律统计信息聚类压缩单元(400)用于对提取的每个聚类点云集区域，按照点云差数和访问次数的点云序列用I个字节编码形成最终压缩统计数据输出。
5.根据权利要求4所述的空间数据访问规律动态统计及信息压缩的系统，其特征在于:所述的空间数据访问规律统计信息表达单元(100)包括空间参数获取模块(101)、点云表达调节因子计算模块(102)和空间点云输出模块(103)； 所述的空间参数获取模块(101)用于获取空间数据访问规律统计信息空间坐标范围、空间数据块数量、空间数据被访问或下载最大次数参数；所述的点云表达调节因子计算模块(102)用于计算空间数据访问规律统计信息中空间数据点云输出调节因子； 所述的空间点云输出模块(103)用于根据空间数据位置和调节因子输出将空间数据访问规律输出成空间点云。
6.根据权利要求4所述的空间数据访问规律动态统计及信息压缩的系统，其特征在于:所述的空间数据访问规律统计信息聚类计算单元(200)包括:聚类梯度计算参数输入模块(201)和聚类梯度计算模块(202)； 所述的聚类梯度计算参数输入模块(201)用于控制聚类计算半径； 所述的聚类梯度计算模块(202 )用于区域内聚类梯度计算。
7.根据权利要求4所述的空间数据访问规律动态统计及信息压缩的系统，其特征在于:所述的空间数据访问规律统计信息区域提取单元(300)包括聚类区域提取参数输入模块(301)和聚类区域提取模块(302); 所述的聚类区域提取参数输入模块(301)用于限制聚类区域提取扩散半径、限制聚类区域内点云差数量； 所述的聚类区域提取模块(302)用于提取聚类区域。
8.根据权利要求4所述的空间数据访问规律动态统计及信息压缩的系统，其特征在于:所述的空间数据访问规律统计信息聚类压缩单元(400)包括聚类区域点云差计算模块(401)、点云序列组织模块(402)和点云序列游程编码压缩模块(403)； 所述的聚类区域点云差计算 模块(401)用于计算聚类区域内点云差； 所述的点云序列组织模块(402)用于将聚类区域按照坐标顺序和点云差数量依序组织； 所述的点云序列游程编码压缩模块(403)用于按照游程编码规则，将按序组织的点云序列进行编码压缩。
【文档编号】H04L12/24GK103458032SQ201310396525
【公开日】2013年12月18日 申请日期:2013年9月3日 优先权日:2013年9月3日
【发明者】潘少明, 李红, 汤戈 申请人:武汉大学