CN105528423A - 顾及空间实例距离权重的自适应同位模式获取方法和装置 - Google Patents
顾及空间实例距离权重的自适应同位模式获取方法和装置 Download PDFInfo
- Publication number
- CN105528423A CN105528423A CN201510891697.0A CN201510891697A CN105528423A CN 105528423 A CN105528423 A CN 105528423A CN 201510891697 A CN201510891697 A CN 201510891697A CN 105528423 A CN105528423 A CN 105528423A
- Authority
- CN
- China
- Prior art keywords
- distance
- pattern
- antithesis
- normalization
- instance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
Landscapes
- Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Image Generation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明提供一种顾及空间实例距离权重的自适应同位模式获取方法和装置;所述方法,包括:将同一空间特征的实例作为目标图层,构建泰森多边形,基于当前目标图层的实例与其它特征实例在泰森多边形影响范围的连接关系,构建归一化实例对偶表;对归一化实例对偶表的距离集合进行分布参数计算,获得归一化距离截断参数;利用归一化距离截断参数构造距离回馈函数,计算任意两个不同特征的连接实例对的距离回馈值;在所述距离回馈阈值下,获取最终的基于距离回馈函数计算的同位模式。解决了常规同位模式挖掘方法在区域适应性、高效性和有效性方面的不足。
Description
技术领域
本发明涉及空间数据模式挖掘领域,尤其涉及一种针对空间点要素数据的自适应的,并且考虑点要素邻近距离权重的同位模式挖掘方法和装置。
背景技术
空间同位模式是指一系列空间特征的组合,这些组合的点实例要素在同一地点频繁的出现,往往表征着它们在空间上有某种普遍的依赖关系。同位模式挖掘目前已经在种群分布、移动通讯、公共安全、环境管理、城市规划等领域有广泛的应用,其挖掘的结果可以为学者进一步研究空间特征分布的内在机理提供线索,也可以辅助决策者进行科学的预判。以城市规划为例,有人发现某一新城大部分的学校周边1000米没有邮政点,而利用北京核心功能区的设施数据进行同位模式挖掘,其结果表明学校和邮政点属于显著度为89%的同位模式,基于该事实,决策者可适量的将邮政引入学校周边,以减小人口流动半径,缓建交通压力,降低各方面的资源浪费。
传统的同位模式挖掘方法普遍需要用户自行设定统一的距离阈值,即使其中有一些由于所处理的数据类型和挖掘目标有所拓展而在原来距离度量的指标上做了改进而换成了其它指标,如把距离阈值改为密度阈值。然而,这种仅凭用户经验来设定距离阈值的方法很不现实:区域具有特异性,区域上的数据分布也有疏密之分,偏小的距离阈值难以挖出有用的模式,偏大的距离阈值会导致大量无用的模式产生。虽然目前已有一些改进方法不需要用户设置距离阈值,但他们普遍利用最近邻思路或迭代优化的方法对算法进行改造,很难对实例间的距离特征有全局的认识,而考虑全局连接关系的效率问题会限制算法在大数据集上的计算。其次,传统方法基于参与度的流行模式判定方法,在计算中并没有考虑实例之间的连接紧密程度,即没有讨论实例之间的距离权重对模式流行程度的影响。而实际上,根据地理学第一定律,相近的事物联系更紧密,所以近距离实例关系对模式流行程度大小的影响,要大于远距离的实例关系。因此,需要构建一种新型有效的计算体系来评价模式的显著性。
总之,目前没有相关的文献或公开的方法,能够在同时考虑区域数据疏密和实例之间距离权重影响的情况下,高效快速的发现空间上的同位模式。
发明内容
本发明旨在解决上面描述的问题。本发明的一个目的是提供一种解决以上问题中的任何一个的顾及空间实例距离权重的自适应同位模式获取方法和装置。具体为:
本发明的一方面提供了一种顾及空间实例距离权重的自适应同位模式获取方法,包括:
将同一空间特征的实例作为目标图层,构建泰森多边形,基于当前目标图层的实例与其它特征实例在泰森多边形影响范围的连接关系,构建归一化实例对偶表;
对归一化实例对偶表的距离集合进行分布参数计算,获得归一化距离截断参数;
利用归一化距离截断参数构造距离回馈函数,计算任意两个不同特征的连接实例对的距离回馈值;
在所述距离回馈阈值下,获取最终的基于距离回馈函数计算的同位模式。
本发明提供的方法还具有如下特点:所述构建归一化实例对偶表的步骤包括:
针对每一类空间特征ei,将其实例作为目标图层,创建泰森多边形图;
基于泰森多边形图所划定的目标实例的影响区间,建立每一个实例与其它特征ej的实例的关联关系,并将这些关联的实例对和它们的距离存储在二维哈希表table_instancevor中;其中,table_instancevor(ei,ej)为二阶有序模式(ei,ej)的Cell信息,表示以ei为目标图层,特征为ej的实例与目标图层实例进行连接所获取的行实例信息,表示特征为ei的第k个实例,Cell信息包括二阶有序模式(ei,ej)的可连接实例对insPair和其对应的距离值dis;
通过公式(1),计算不同特征的连接实例对之间的归一化距离,替换所述二维哈希表的距离值:
其中,dis(a,b)表示两个实例a,b之间的欧氏距离,
DisVP表示table_instancevor中所有的实例对的距离集合,
min(T)和max(T)分别表示取集合T的最小和最大值;
对所述二维哈希表进行无序处理,得到归一化实例对偶表步骤包括:
把满足j>i的table_instancevor(ej,ei)中的任意有序二阶模式Cell信息的实例对顺序互换后,得到新的二维哈希表
将所述新的二维哈希表中满足j>i的的Cell信息与原table_instancevor(ei,ej)对应的Cell信息进行求并处理,得到最终的归一化实例对偶表
本发明提供的方法还具有如下特点:对于实例数大于12000的数据集,利用所述归一化实例对偶表中的归一化距离集合,在显著水平α=0.05的前提下,进行广义极大值分布函数的参数拟合,对于χ2检验中拟合优度大于0.8的情况,设定所述的归一化距离截断为μ±stderr;其中μ表示广义极值极大值拟合函数的位置参数,stderr表示拟合的残差;或者,
对于其它不符合实例数和拟合优度条件的情况,将归一化距离的均值作为截断值。
本发明提供的方法还具有如下特点:利用所述归一化距离截断参数,构建距离回馈函数,计算任意两个不同特征的连接实例对的距离回馈值:
其中,lclose=μ-stderr,且lfar=μ+stderr。
本发明提供的方法还具有如下特点:在所述距离回馈阈值下,获取最终的基于距离回馈函数计算的同位模式的步骤包括:
以所述归一化实例对偶表的行列位置所对应的特征对{ei,ej}为初始的候选同位模式c,迭代以下子步骤1-5,直到|c|阶同位模式没有大于距离回馈阈值的成员,其中,j>i,
当|c|=2时,省略子步骤1-2,直接跳转到子步骤3;
子步骤1:对于候选模式c={ei...ek},根据所述归一化实例对偶表中的实例对,求出候选模式的团实例insCliquec,其中,候选模式c中包含的特征按照字典序进行排列,即满足k≥j>i,insCliquec是指基于泰森多边形,在空间上有全连接关系的不同特征的实例的组合,即团实例中抽取的任意两个实例对,在归一化实例对偶表中都能找到对应项;
子步骤2:根据公式(3),求算候选模式c中包含的任意两个二阶模式{ei,ej}的实例对在insCliquec约束下的距离回馈值:
其中,表示求取二阶模式{ei,ej}在insCliquec上的投影;
子步骤3:通过公式(4),求算二阶模式{ei,ej}在c={ei...ek}约束下的距离回馈值:
其中,表示二阶模式{ei,ej}对应的归一化实例对偶表中的所有实例对集合,avg表示求算平均值;
子步骤4:通过公式(5),求算模式c={ei...ek}的基于距离回馈函数计算的流行程度:
Rewardc=mini<j≤k{Rewardc(ei,ej)}(5),
将大于给定距离回馈阈值θ的候选模式判定为最终的|c|阶同位模式,即最终的同位模式需满足Rewardc>θ;
子步骤5:以|c|阶同位模式为依据,将前|c|-1位具有相同特征的同位模式求并,获取|c|+1阶候选同位模式c’,令c=c’。
本发明的另一个方面还提供了一种顾及空间实例距离权重的自适应同位模式获取装置,包括:
归一化实例对偶表构建模块,用于将同一空间特征的实例作为目标图层,构建泰森多边形,基于当前目标图层的实例与其它特征实例在泰森多边形影响范围的连接关系,构建归一化实例对偶表;
归一化距离截断参数获取模块,用于对归一化实例对偶表的距离集合进行分布参数计算,获得归一化距离截断参数;
实例对距离回馈值计算模块,用于利用归一化距离截断参数构造距离回馈函数,计算任意两个不同特征的连接实例对的距离回馈值;
同位模式获取模块,用于在所述距离回馈阈值下,获取最终的基于距离回馈函数计算的同位模式。
本发明提供的装置还具有如下特点:所述归一化实例对偶表构建模块具体用于:
泰森多边形创建单元,用于针对每一类空间特征ei,将其实例作为目标图层,创建泰森多边形图;
二维哈希表创建单元,用于基于泰森多边形图所划定的目标实例的影响区间,建立每一个实例与其它特征ej的实例的关联关系,并将这些关联的实例对和它们的距离存储在二维哈希表table_instancevor中;其中,table_instancevor(ei,ej)为二阶有序模式(ei,ej)的Cell信息,表示以ei为目标图层,特征为ej的实例与目标图层实例进行连接所获取的行实例信息,表示特征为ei的第k个实例,Cell信息包括二阶有序模式(ei,ej)的可连接实例对insPair和其对应的距离值dis;
归一化距离计算单元,用于计算不同特征的连接实例对之间的归一化距离,替换所述二维哈希表的距离值,见公式(1):
其中,dis(a,b)表示两个实例a,b之间的欧氏距离,
DisVP表示table_instancevor中所有的实例对的距离集合,
min(T)和max(T)分别表示取集合T的最小和最大值;
无序处理单元,用于对所述二维哈希表进行无序处理,得到归一化实例对偶表
所述无序处理单元,具体用于:
把满足j>i的table_instancevor(ej,ei)中的任意有序二阶模式Cell信息的实例对顺序互换后,得到新的二维哈希表
将所述新的二维哈希表中满足j>i的的Cell信息与原table_instancevor(ei,ej)对应的Cell信息进行求并处理,得到最终的归一化实例对偶表
本发明提供的装置还具有如下特点:所述归一化距离截断参数获取模块具体用于:
对于实例数大于12000的数据集,利用所述归一化实例对偶表中的归一化距离集合,在显著水平α=0.05的前提下,进行广义极大值分布函数的参数拟合,对于χ2检验中拟合优度大于0.8的情况,设定所述的归一化距离截断为μ±stderr;其中μ表示广义极值极大值拟合函数的位置参数,stderr表示拟合的残差;或者,
对于其它不符合实例数和拟合优度条件的情况,将归一化距离的均值作为截断值。
本发明提供的装置还具有如下特点:所述实例对距离回馈值计算模块具体用于,计算任意两个不同特征的连接实例对的距离回馈值,见公式(2):
其中,lclose=μ-stderr,且lfar=μ+stderr。
本发明提供的装置还具有如下特点:所述同位模式获取模块,具体用于:
以所述归一化实例对偶表的行列位置所对应的特征对{ei,ej}为初始的候选同位模式c,迭代以下子步骤1-5,直到|c|阶同位模式没有大于距离回馈阈值的成员,其中,j>i,
当|c|=2时,省略子步骤1-2,直接跳转到子步骤3;
子步骤1:对于候选模式c={ei...ek},根据所述归一化实例对偶表中的实例对,求出候选模式的团实例insCliquec,其中,候选模式c中包含的特征按照字典序进行排列,即k≥j>i,insCliquec是指基于泰森多边形,在空间上有全连接关系的不同特征的实例的组合,即团实例中抽取的任意两个实例对,在归一化实例对偶表中都能找到对应项;
子步骤2:根据公式(3),求算候选模式c中包含的任意两个二阶模式{ei,ej}的实例对在insCliquec约束下的距离回馈值:
其中,表示求取二阶模式{ei,ej}在insCliquec上的投影;
子步骤3:通过公式(4),求算二阶模式{ei,ej}在c={ei...ek}约束下的距离回馈值:
其中,表示二阶模式{ei,ej}对应的归一化实例对偶表中的所有实例对集合,avg表示求算平均值;
子步骤4:通过公式(5),求算模式c={ei...ek}的基于距离回馈函数计算的流行程度:
Rewardc=mini<j≤k{Rewardc(ei,ej)}(5),
将大于给定距离回馈阈值θ的候选模式判定为最终的|c|阶同位模式,即最终的同位模式需满足Rewardc>θ;
子步骤5:以|c|阶同位模式为依据,将前|c|-1位具有相同特征的同位模式求并,获取|c|+1阶候选同位模式c’,令c=c’。
同位模式是一组空间特征的子集,表示他们的实例组合在空间上频繁的出现。现有相关技术中,很难在同时考虑区域数据疏密和实例之间距离权重影响的情况下,高效快速的发现空间上的同位模式。
本发明针对该问题,提出了一种基于泰森多边形和距离回馈函数的同位模式挖掘方式,利用泰森多边形的特性,在实例的影响区域内确立不同特征实例间的领域关系,并以此为依据,获得实例之间的极大距离统计特征,而后将特征参数用于构建候选模式的距离回馈函数,来度量模式的流行程度。
该方法的优点在于:1)省去了用户预先设定距离阈值的操作,规避了在未知区域上进行挖掘的不确定性;
2)考虑了实例之间距离权重对模式流行程度的影响,更符合实际应用场景;
3)泰森多边形的使用,极大降低了实例对连接的计算量,使得运算效率大幅提高。
参照附图来阅读对于示例性实施例的以下描述,本发明的其它特性特征和优点将变得清晰。
附图说明
并入到说明书中并且构成说明书的一部分的附图示出了本发明的实施例,并且与描述一起用于解释本发明的原理。在这些附图中,类似的附图标记用于表示类似的要素。下面描述中的附图是本发明的一些实施例,而不是全部实施例。对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,可以根据这些附图获得其它的附图。
图1为本发明的实施例一提供的一种顾及空间实例距离权重的自适应同位模式获取方法的流程图;
图2为构建归一化实例对偶表的方法示意图;
图3为从二阶候选同位模式,获取最终的基于距离回馈函数计算的同位模式的迭代流程图;
图4为本发明的实施例二提供的一种顾及空间实例距离权重的自适应同位模式获取装置的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互任意组合。
为了更好的对本发明的实施例所提供的技术方案进行阐述,首先对如下概念进行说明:
泰森多边形:是由一组连接两邻接点的直线的垂直平分线组成的连续多边形几何。它按照最邻近原则划分平面,每个点与它的最近邻区域相关联。
SHP数据类型:shapefile的简称,是ESRI公司开发的一种空间数据开放格式。
有序模式:是指一系列有排序的空间特征的向量。
广义极值分布:符合高峰胖尾的曲线特征,是极值理论的一个重要分支,统一了Gumbel,Fréchet和Weibull三种极值分布模型,用于模拟变量为极值序列的应用场景,在气象、经济学上已有广泛的应用。该分布函数涉及三个关键参数,即位置参数μ,尺度参数σ,和形状参数k。
下面结合附图,对本发明的实施例一进行说明。
本发明提供一种顾及空间实例距离权重的自适应同位模式获取方法,输入为欧式空间的点状SHP数据P={p1,p2,...pn},以及这些数据所涉及到的空间特征类型E={e1,e2,...em},其中,每个点实例包含编号、类型和X、Y坐标(为了强调实例的空间特征,以下将特征为ei的第j个实例表示为)。该方法的流程如图1所示,包括:
步骤101、将同一空间特征的实例作为目标图层,构建泰森多边形,基于当前目标图层的实例与其它特征实例在泰森多边形影响范围的连接关系,构建归一化实例对偶表;
归一化实例对偶表的构建理念源于自然界的一个性质:大多数物体在空间上会处于它们形成的场力的平衡点。两种类型(A、B)的物体若在空间上存在某种关系,那属于类型A的物体更倾向于受到它周边场力最大的B类型的物体的影响,反之亦然。由此可推理出,若A,B类型相互吸引,那么它们的实例大多会分布在彼此影响区域的核地带,反之则会分布在彼此场力最薄弱的边缘地带。由于其连接操作被限制在单个实例的影响范围内,因此其连接数相比传统的方法,有很大程度的减少。该步骤涉及三个子步骤:
子步骤1:利用泰森多边形来确定所有特征对应实例的领域对象。对于每一类空间特征ei,将该特征下的实例作为目标图层,创建泰森多边形图。这样,每一个实例相关的其它特征的实例可表示为:
上式中,可以在地理信息系统中通过叠加分析进行计算;表示特征为非ei的实例集合;dis(a,b)表示两个空间实例a,b的欧氏距离。
子步骤2:将上一步中获取的实例对和它们之间的距离值存储在m*m的二维哈希表中,用table_instancevor表示。table_instancevor(ei,ej)即表示以ei为目标图层,特征为ej的实例与目标图层实例进行连接所获取的行实例信息,以下简称为二阶有序模式(ei,ej)的Cell信息。Cell信息包括二阶有序模式(ei,ej)的可连接实例对insPair和其对应的距离值dis。
通过公式(2),计算不同特征的连接实例对之间的归一化距离,替换所述二维哈希表的距离值。
上式中,dis(a,b)表示两个实例a,b之间的欧氏距离;DisVP表示table_instancevor中所有的实例对的距离集合;min(T)和max(T)分别表示取集合T的最小和最大值。
子步骤3:将table_instancevor中有序模式的实例对变为无序,即把满足j>i的table_instancevor(ej,ei)中的任意有序二阶模式Cell信息的实例对顺序互换后,得到新的二维哈希表然后,将新的二维哈希表中满足j>i的的Cell信息与原table_instancevor(ei,ej)对应的Cell信息进行求并处理,得到最终的归一化实例对偶表
为了更清楚的说明步骤101的过程,图2给出了图形化的示例。假设空间X中有一系列特征为E={A,B,...,et}的实例P={A1,A2,...A5,B1,...B9,C1,...C5,...},图2(a)展示了以特征A的实例为目标图层构建的泰森多边形(为了简化说明,只以阴影区域为提取示例),其中,实例A5的对应连接实例为B4、B8、B7和C3。对所有特征做完实例对连接和距离归一化处理后,形成了图2(b)所示的二维哈希表。最后,对table_instancevor进行无序化处理,形成图3(c)的
步骤102、对归一化实例对偶表的距离集合进行分布参数计算,获得归一化距离截断参数;
该步骤通过计算空间实例与其影响范围内的其它特征实例的距离的统计分布规律,来获得指代这个区域大概数据密度特征的两个“截断参数”,用以取代传统同位模式挖掘的距离阈值设定过程。
在归一化实例对偶表的构建过程中,基于泰森多边形求取当前实例最靠近目标特征的实例,是一种区域极大值的求值策略,因此,可以推测出的“归一化距离集合(记作NordisVP)”的分布更倾向于广义极值分布规律。实际上,本发明利用北京的城市设施点状数据进行了分区同位模式挖掘的实验,8个数据量都在1.2万以上的连续平面区域的NordisVP,在显著水平α=0.05的χ2拟合检测中,都符合广义极大值的统计分布规律,拟合优度均在0.89以上,最高达到0.9839,满足了人类活动现象的统计拟合要求。如此一来,在中,那些具有潜在邻近关系的特征对,通常包含更多相对“变量的拟合峰值μ±stderr(α=0.05)”而偏小的归一化距离值,反之,则包含更多相对于拟合峰值而偏大的归一化距离值。其中,μ表示广义极值极大值拟合函数的位置参数,stderr表示拟合的残差。
根据上述推理和所述的实验结果,在实例数大于12000,且显著性水平α=0.05的χ2拟合检测中,若拟合优度大于0.8,本发明将μ±stderr(α=0.05)的两个边界值作为截断参数,这两个值把[0,1]距离轴剖分为三部分:1)>μ+stderr(α=0.05),指代远距离关系;2)<μ-stderr(α=0.05),指代近距离关系;3)[μ-stderr,μ+stderr],表征归一化距离在该区间出现的概率最大,反映的是分布的平均特征。对于其它不符合实例数和拟合优度条件的情况,将归一化距离的均值作为截断值。
步骤103、利用归一化距离截断参数构造距离回馈函数,计算任意两个不同特征的连接实例对的距离回馈值;
根据地理学第一定律,近距离关系的实例对,对模式流行程度的贡献,要大于远距离关系的实例对。基于这一设定,本发明构建了一个单调递减函数,见公式(3),用以计算两个不同类型的实例对的距离回馈值,该函数的取值范围为[-1,1]。
上式中,lclose=μ-stderr,且lfar=μ+stderr。
步骤104、在给定的距离回馈阈值下,获取最终的基于距离回馈函数计算的同位模式。
本步骤以归一化实例对偶表的行列位置所对应的特征对{ei,ej}为初始的候选同位模式c,迭代子步骤1-5,直到|c|阶同位模式没有大于距离回馈阈值的成员,该步骤的迭代流程如图3所示。其中,j>i。当|c|=2时,省略子步骤1-2,直接跳转到子步骤3。
子步骤1:对于候选模式c={ei...ek},根据所述归一化实例对偶表中的实例对,求出候选模式的团实例insCliquec。其中,候选模式c中包含的特征按照字典序进行排列,即k≥j>i,insCliquec是指基于泰森多边形,在空间上有全连接关系的不同特征的实例的组合,即团实例中抽取的任意两个实例对,在归一化实例对偶表中都能找到对应项。
子步骤2:根据公式(3),求算候选模式c中包含的任意两个二阶模式{ei,ej}的实例对在insCliquec约束下的距离回馈值。参与团实例构建的实例对,保留其原来的距离回馈值;不能参与构成团实例的实例对,对模式的流行程度会有一个逆向的作用,而这种逆向作用,通过求取其距离回馈值的绝对值负数来定义。
上式中,表示求取二阶模式{ei,ej}在insCliquec上的投影。
子步骤3:通过公式(4),求算二阶模式{ei,ej}在c={ei...ek}约束下的距离回馈值。该值的计算同时顾忌了远距离和近距离实例对,对模式流行程度的影响,因此在模式的显著度判定上,比传统的仅依赖“团实例所属模式的所有单个特征上的实例投影个数和其相同特征的实例个数的比值的最小值”的流行指数计算方法,更有说服力。
上式中,表示二阶模式{ei,ej}对应的归一化实例对偶表中的所有实例对集合,avg表示求算平均值。
子步骤4:通过公式(5),求算模式c={ei...ek}的基于距离回馈函数计算的流行程度。
Rewardc=mini<j≤k{Rewardc(ei,ej)}(5)
从公式(5)可以看出,该值的取值范围为(-1,1),因此用户可以模拟流行阈值的设定方式,为模式的距离回馈值设定下限,以获得层次性的同位模式挖掘结果,即将大于给定距离回馈阈值θ的候选模式判定为最终的|c|阶同位模式。
子步骤5:以|c|阶同位模式为依据,将前|c|-1位具有相同特征的同位模式求并,获取|c|+1阶候选同位模式c’,令c=c’。
图4为本发明提供的顾及空间实例距离权重的自适应同位模式获取装置。图4所示装置,包括:
归一化实例表构建模块401,用于将同一空间特征的实例为目标图层,构建泰森多边形,基于当前目标图层的实例与其它特征实例在泰森多边形影响范围的连接关系,构建归一化实例对偶表;
归一化距离截断参数获取模块402,用于对归一化实例对偶表的距离集合进行分布参数计算,获得归一化距离截断参数;
实例对距离回馈值计算模块403,用于利用归一化距离截断参数构造的距离回馈函数,计算任意两个不同特征的连接实例对的距离回馈值;
同位模式获取模块404,用于在所述距离回馈阈值下,获取最终的基于距离回馈函数计算的同位模式。
其中,所述构建模块401具体用于:
泰森多边形创建单元,用于针对每一类空间特征ei,将其实例作为目标图层,创建泰森多边形图;
二维哈希表创建单元,用于基于泰森多边形图所划定的目标实例的影响区间,建立每一个实例与其它特征ej的实例的关联关系,并将这些关联的实例对和它们的距离存储在二维哈希表table_instancevor中;其中,table_instancevor(ei,ej)为二阶有序模式(ei,ej)的Cell信息,表示以ei为目标图层,特征为ej的实例与目标图层实例进行连接所获取的行实例信息,Cell信息包括二阶有序模式(ei,ej)的可连接实例对insPair和其对应的距离值dis;
归一化距离计算单元,用于计算不同特征的连接实例对之间的归一化距离,替换所述二维哈希表的距离值,见公式(1):
上式中,dis(a,b)表示两个实例a,b之间的欧氏距离,
DisVP表示table_instancevor中所有的实例对的距离集合,
min(T)和max(T)分别表示取集合T的最小和最大值;
无序处理单元,用于对所述二维哈希表进行无序处理,得到归一化实例对偶表
其中,所述无序处理单元,具体用于:
把满足j>i的table_instancevor(ej,ei)中的任意有序二阶模式Cell信息的实例对顺序互换后,得到新的二维哈希表
将所述新的二维哈希表中满足j>i的中的Cell信息与原table_instancevor(ei,ej)对应的Cell信息进行求并处理,得到最终的归一化实例对偶表
其中,所述归一化截断参数获取模块402,具体用于:
对于实例数大于12000的数据集,利用所述归一化实例对偶表中的归一化距离集合,在显著水平α=0.05的前提下,进行广义极大值分布函数的参数拟合,对于χ2检验中拟合优度大于0.8的情况,设定所述的归一化距离截断为μ±stderr;其中μ表示广义极值极大值拟合函数的位置参数,stderr表示拟合的残差;或者,
对于其它不符合实例数和拟合优度条件的情况,将归一化距离的均值作为截断值。
其中,所述实例对距离回馈值计算模块403,具体用于,计算任意两个不同特征的连接实例对的距离回馈值,见公式(2):
上式中,lclose=μ-stderr,且lfar=μ+stderr。
其中,所述同位模式获取模块404,具体用于:
以所述归一化实例对偶表的行列位置所对应的特征对{ei,ej}为初始的候选同位模式c,迭代以下子步骤1-5,直到|c|阶同位模式没有大于距离回馈阈值的成员,其中,j>i,
当|c|=2时,省略子步骤1-2,直接跳转到子步骤3;
子步骤1:对于候选模式c={ei...ek},根据所述归一化实例对偶表中的实例对,求出候选模式的团实例insCliquec,其中,候选模式c中包含的特征按照字典序进行排列,即k≥j>i,insCliquec是指基于泰森多边形,在空间上有全连接关系的不同特征的实例的组合,即团实例中抽取的任意两个实例对,在归一化实例对偶表中都能找到对应项;
子步骤2:根据公式(3),求算候选模式c中包含的任意两个二阶模式{ei,ej}的实例对在insCliquec约束下的距离回馈值:
上式中,表示求取二阶模式{ei,ej}在insCliquec上的投影;
子步骤3:通过公式(4),求算二阶模式{ei,ej}在c={ei...ek}约束下的距离回馈值:
上式中,表示二阶模式{ei,ej}对应的归一化实例对偶表中的所有实例对集合,avg表示求算平均值;
子步骤4:通过公式(5),求算模式c={ei...ek}的基于距离回馈函数计算的流行程度:
Rewardc=mini<j≤k{Rewardc(ei,ej)}(5),
将大于给定距离回馈阈值θ的候选模式判定为最终的|c|阶同位模式,即最终的同位模式需满足Rewardc>θ;
子步骤5:以|c|阶同位模式为依据,将前|c|-1位具有相同特征的同位模式求并,获取|c|+1阶候选同位模式c’,令c=c’。
本发明提供的方法和装置实施例,与现有技术中的同位模式挖掘方法相比,在区域适应性、高效性和有效性方面享有优势:第一,不需要用户预先设定距离阈值,而是从区域实例与其影响范围内的实例的距离分布规律,得出相应的距离截断参数来替代距离阈值,因此挖掘不受限于用户对数据的了解程度,具有更强的区域适用性;第二,由于使用了泰森多边形,使得连接计算限制在实例的影响范围内,取代了传统的单一距离阈值限定的连接思路,因此大幅度减小了连接计算量,改善了挖掘的效率;第三,考虑了实例连接的距离权重,尊重了地理学第一定律,因此是一种更切合实际的同位模式挖掘方法,说服力也更强。上面描述的内容可以单独地或者以各种方式组合起来实施,而这些变型方式都在本发明的保护范围之内。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制。尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种顾及空间实例距离权重的自适应同位模式获取方法,其特征在于,包括:
将同一空间特征的实例作为目标图层,构建泰森多边形,基于当前目标图层的实例与其它特征实例在泰森多边形影响范围的连接关系,构建归一化实例对偶表;
对归一化实例对偶表的距离集合进行分布参数计算,获得归一化距离截断参数;
利用归一化距离截断参数构造距离回馈函数,计算任意两个不同特征的连接实例对的距离回馈值;
在所述距离回馈阈值下,获取最终的基于距离回馈函数计算的同位模式。
2.根据权利要求1所述的方法,其特征在于,所述基于泰森多边形,构建归一化实例对偶表的步骤包括:
针对每一类空间特征ei,将其实例作为目标图层,创建泰森多边形图;
基于泰森多边形图所划定的目标实例的影响区间,建立每一个实例与其它特征ej的实例的关联关系,并将这些关联的实例对和它们的距离存储在二维哈希表table_instancevor中;其中,table_instancevor(ei,ej)为二阶有序模式(ei,ej)的Cell信息,表示以ei为目标图层,特征为ej的实例与目标图层实例进行连接所获取的行实例信息,表示特征为ei的第k个实例,Cell信息包括二阶有序模式(ei,ej)的可连接实例对insPair和其对应的距离值dis;
通过公式(1),计算不同特征的连接实例对之间的归一化距离,替换所述二维哈希表的距离值:
其中,dis(a,b)表示两个实例a,b之间的欧氏距离,
DisVP表示table_instancevor中所有的实例对的距离集合,
min(T)和max(T)分别表示取集合T的最小和最大值;
对所述二维哈希表进行无序处理,得到归一化实例对偶表步骤包括:
把满足j>i的table_instancevor(ej,ei)中的任意有序二阶模式Cell信息的实例对顺序互换后,得到新的二维哈希表
将所述新的二维哈希表中满足j>i的的Cell信息与原table_instancevor(ei,ej)对应的Cell信息进行求并处理,得到最终的归一化实例对偶表
3.根据权利要求1所述的方法,其特征在于:
对于实例数大于12000的数据集,利用所述归一化实例对偶表中的归一化距离集合,在显著水平α=0.05的前提下,进行广义极大值分布函数的参数拟合,对于χ2检验中拟合优度大于0.8的情况,设定所述的归一化距离截断为μ±stderr;其中μ表示广义极值极大值拟合函数的位置参数,stderr表示拟合的残差;或者,
对于其它不符合实例数和拟合优度条件的情况,将归一化距离的均值作为截断值。
4.根据权利要求1所述的方法,其特征在于,利用所述归一化距离截断参数,构建距离回馈函数,计算任意两个不同特征的连接实例对的距离回馈值:
其中,lclose=μ-stderr,且lfar=μ+stderr。
5.根据权利要求1所述的方法,其特征在于,在所述距离回馈阈值下,获取最终的基于距离回馈函数计算的同位模式的步骤包括:
以所述归一化实例对偶表的行列位置所对应的特征对{ei,ej}为初始的候选同位模式c,迭代以下子步骤1-5,直到|c|阶同位模式没有大于距离回馈阈值的成员,其中,j>i,
当|c|=2时,省略子步骤1-2,直接跳转到子步骤3;
子步骤1:对于候选模式c={ei...ek},根据所述归一化实例对偶表中的实例对,求出候选模式的团实例insCliquec,其中,候选模式c中包含的特征按照字典序进行排列,即满足k≥j>i,insCliquec是指基于泰森多边形,在空间上有全连接关系的不同特征的实例的组合,即团实例中抽取的任意两个实例对,在归一化实例对偶表中都能找到对应项;
子步骤2:根据公式(3),求算候选模式c中包含的任意两个二阶模式{ei,ej}的实例对在insCliquec约束下的距离回馈值:
其中,表示求取二阶模式{ei,ej}在insCliquec上的投影;
子步骤3:通过公式(4),求算二阶模式{ei,ej}在c={ei...ek}约束下的距离回馈值:
其中,表示二阶模式{ei,ej}对应的归一化实例对偶表中的所有实例对集合,avg表示求算平均值;
子步骤4:通过公式(5),求算模式c={ei...ek}的基于距离回馈函数计算的流行程度:
Rewardc=mini<j≤k{Rewardc(ei,ej)}(5),
将大于给定距离回馈阈值θ的候选模式判定为最终的|c|阶同位模式,即最终的同位模式需满足Rewardc>θ;
子步骤5:以|c|阶同位模式为依据,将前|c|-1位具有相同特征的同位模式求并,获取|c|+1阶候选同位模式c’,令c=c’。
6.一种顾及空间实例距离权重的自适应同位模式获取装置,其特征在于,包括:
归一化实例对偶表构建模块,用于将同一空间特征的实例作为目标图层,构建泰森多边形,基于当前目标图层的实例与其它特征实例在泰森多边形影响范围的连接关系,构建归一化实例对偶表;
归一化距离截断参数获取模块,用于对归一化实例对偶表的距离集合进行分布参数计算,获得归一化距离截断参数;
实例对距离回馈值计算模块,用于利用归一化距离截断参数构造距离回馈函数,计算任意两个不同特征的连接实例对的距离回馈值;
同位模式获取模块,用于在所述距离回馈阈值下,获取最终的基于距离回馈函数计算的同位模式。
7.根据权利要求6所述的装置,其特征在于,所述归一化实例对偶表构建模块具体用于:
泰森多边形创建单元,用于针对每一类空间特征ei,将其实例作为目标图层,创建泰森多边形图;
二维哈希表创建单元,用于基于泰森多边形图所划定的目标实例的影响区间,建立每一个实例与其它特征ej的实例的关联关系,并将这些关联的实例对和它们的距离存储在二维哈希表table_instancevor中;其中,table_instancevor(ei,ej)为二阶有序模式(ei,ej)的Cell信息,表示以ei为目标图层,特征为ej的实例与目标图层实例进行连接所获取的行实例信息,表示特征为ei的第k个实例,Cell信息包括二阶有序模式(ei,ej)的可连接实例对insPair和其对应的距离值dis;
归一化距离计算单元,用于计算不同特征的连接实例对之间的归一化距离,替换所述二维哈希表的距离值,见公式(1):
其中,dis(a,b)表示两个实例a,b之间的欧氏距离,
DisVP表示table_instancevor中所有的实例对的距离集合,
min(T)和max(T)分别表示取集合T的最小和最大值;
无序处理单元,用于对所述二维哈希表进行无序处理,得到归一化实例对偶表
所述无序处理单元,具体用于:
把满足j>i的table_instancevor(ej,ei)中的任意有序二阶模式Cell信息的实例对顺序互换后,得到新的二维哈希表
将所述新的二维哈希表中满足j>i的的Cell信息与原table_instancevor(ei,ej)对应的Cell信息进行求并处理,得到最终的归一化实例对偶表
8.根据权利要求6所述的装置,其特征在于,所述归一化距离截断参数获取模块具体用于:
对于实例数大于12000的数据集,利用所述归一化实例对偶表中的归一化距离集合,在显著水平α=0.05的前提下,进行广义极大值分布函数的参数拟合,对于χ2检验中拟合优度大于0.8的情况,设定所述的归一化距离截断为μ±stderr;其中μ表示广义极值极大值拟合函数的位置参数,stderr表示拟合的残差;或者,
对于其它不符合实例数和拟合优度条件的情况,将归一化距离的均值作为截断值。
9.根据权利要求6所述的装置,其特征在于,所述实例对距离回馈值计算模块具体用于,计算任意两个不同特征的连接实例对的距离回馈值,见公式(2):
其中,lclose=μ-stderr,且lfar=μ+stderr。
10.根据权利要求6所述的装置,其特征在于,所述同位模式获取模块,具体用于:
以所述归一化实例对偶表的行列位置所对应的特征对{ei,ej}为初始的候选同位模式c,迭代以下子步骤1-5,直到|c|阶同位模式没有大于距离回馈阈值的成员,其中,j>i,
当|c|=2时,省略子步骤1-2,直接跳转到子步骤3;
子步骤1:对于候选模式c={ei...ek},根据所述归一化实例对偶表中的实例对,求出候选模式的团实例insCliquec,其中,候选模式c中包含的特征按照字典序进行排列,即k≥j>i,insCliquec是指基于泰森多边形,在空间上有全连接关系的不同特征的实例的组合,即团实例中抽取的任意两个实例对,在归一化实例对偶表中都能找到对应项;
子步骤2:根据公式(3),求算候选模式c中包含的任意两个二阶模式{ei,ej}的实例对在insCliquec约束下的距离回馈值:
其中,表示求取二阶模式{ei,ej}在insCliquec上的投影;
子步骤3:通过公式(4),求算二阶模式{ei,ej}在c={ei...ek}约束下的距离回馈值:
其中,表示二阶模式{ei,ej}对应的归一化实例对偶表中的所有实例对集合,avg表示求算平均值;
子步骤4:通过公式(5),求算模式c={ei...ek}的基于距离回馈函数计算的流行程度:
Rewardc=mini<j≤k{Rewardc(ei,ej)}(5),
将大于给定距离回馈阈值θ的候选模式判定为最终的|c|阶同位模式,即最终的同位模式需满足Rewardc>θ;
子步骤5:以|c|阶同位模式为依据,将前|c|-1位具有相同特征的同位模式求并,获取|c|+1阶候选同位模式c’,令c=c’。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510891697.0A CN105528423B (zh) | 2015-12-07 | 2015-12-07 | 顾及空间实例距离权重的自适应同位模式获取方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510891697.0A CN105528423B (zh) | 2015-12-07 | 2015-12-07 | 顾及空间实例距离权重的自适应同位模式获取方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105528423A true CN105528423A (zh) | 2016-04-27 |
CN105528423B CN105528423B (zh) | 2018-08-17 |
Family
ID=55770646
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510891697.0A Active CN105528423B (zh) | 2015-12-07 | 2015-12-07 | 顾及空间实例距离权重的自适应同位模式获取方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105528423B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106780262A (zh) * | 2017-01-13 | 2017-05-31 | 中国科学院遥感与数字地球研究所 | 一种考虑城市道路网络约束的同位模式发现方法及装置 |
CN107194541A (zh) * | 2017-04-24 | 2017-09-22 | 国网江西省电力公司赣州供电分公司 | 一种基于自适应权重Voronoi图的配电网供电分区方法 |
CN113157717A (zh) * | 2021-05-26 | 2021-07-23 | 深圳平安智汇企业信息管理有限公司 | 长数据链路的缓存刷新方法、装置、设备及存储介质 |
CN116229442A (zh) * | 2023-01-03 | 2023-06-06 | 武汉工程大学 | 一种文本图像合成和实例化权重的迁移学习方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103631882A (zh) * | 2013-11-14 | 2014-03-12 | 北京邮电大学 | 基于图挖掘技术的语义化业务生成系统和方法 |
CN104794201A (zh) * | 2015-04-22 | 2015-07-22 | 哈尔滨工业大学深圳研究生院 | 一种基于Voronoi图与地标的层次化空间位置描述方法及装置 |
-
2015
- 2015-12-07 CN CN201510891697.0A patent/CN105528423B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103631882A (zh) * | 2013-11-14 | 2014-03-12 | 北京邮电大学 | 基于图挖掘技术的语义化业务生成系统和方法 |
CN104794201A (zh) * | 2015-04-22 | 2015-07-22 | 哈尔滨工业大学深圳研究生院 | 一种基于Voronoi图与地标的层次化空间位置描述方法及装置 |
Non-Patent Citations (4)
Title |
---|
SEUNG KWAN KIM ET AL: "Maximal Cliques Generating Algorithm for Spatial Co-location Pattern Mining", 《INTERNATIONAL CONFERENCE ON SECURE AND TRUST COMPUTING, DATA MANAGEMENT, AND APPLICATION》 * |
VENKATESAN MEENAKSHI SUNDARAM ET AL: "Discovering co-location patterns from spatial domain using a Delaunay approach", 《PROCEDIA ENGINEERING》 * |
周剑云 等: "特定空间对象同位模式挖掘算法研究", 《电脑知识与技术》 * |
李中元 等: "空间同位模式挖掘研究进展", 《地理空间信息》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106780262A (zh) * | 2017-01-13 | 2017-05-31 | 中国科学院遥感与数字地球研究所 | 一种考虑城市道路网络约束的同位模式发现方法及装置 |
CN107194541A (zh) * | 2017-04-24 | 2017-09-22 | 国网江西省电力公司赣州供电分公司 | 一种基于自适应权重Voronoi图的配电网供电分区方法 |
CN113157717A (zh) * | 2021-05-26 | 2021-07-23 | 深圳平安智汇企业信息管理有限公司 | 长数据链路的缓存刷新方法、装置、设备及存储介质 |
CN116229442A (zh) * | 2023-01-03 | 2023-06-06 | 武汉工程大学 | 一种文本图像合成和实例化权重的迁移学习方法 |
CN116229442B (zh) * | 2023-01-03 | 2024-05-28 | 武汉工程大学 | 一种文本图像合成和实例化权重的迁移学习方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105528423B (zh) | 2018-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107862738B (zh) | 一种基于移动激光测量点云进行室内结构化三维重建方法 | |
Hansen | GIS-based multi-criteria analysis of wind farm development | |
CN107113764B (zh) | 提高人工神经网络定位性能的方法和装置 | |
CN103247041B (zh) | 一种基于局部采样的多几何特征点云数据的分割方法 | |
Liu et al. | On coverage of wireless sensor networks for rolling terrains | |
CN103401922B (zh) | 无线传感器网络中基于博弈方法的分布式定位装置与方法 | |
CN105528423A (zh) | 顾及空间实例距离权重的自适应同位模式获取方法和装置 | |
CN105512759A (zh) | 一种基于生物地理学优化算法的城市ca模型参数优化方法 | |
CN103648164B (zh) | 一种基于到达时间差和Gossip算法的无线传感器网络分布式定位方法 | |
CN105761312A (zh) | 一种微地形表面重建方法 | |
CN110346654B (zh) | 基于普通克里金插值的电磁频谱地图构建方法 | |
CN106600617A (zh) | 基于曲率从Lidar点云数据提取建筑物轮廓线的方法 | |
CN102708327A (zh) | 一种基于谱优化的网络社区发现方法 | |
CN103729846B (zh) | 基于不规则三角网的LiDAR点云数据边缘检测方法 | |
CN109947879A (zh) | 一种基于复杂网络的海洋观测大数据可视化分析方法 | |
Thong | Soft computing methods for WiMax network planning on 3D geographical information systems | |
CN116518960B (zh) | 路网更新方法、装置、电子设备和存储介质 | |
Amirabadizadeh et al. | Downscaling daily precipitation and temperatures over the Langat River Basin in Malaysia: a comparison of two statistical downscaling approaches | |
Dimitrios | Urban growth prediction modelling using fractals and theory of chaos | |
CN105825430A (zh) | 一种基于异构社会网络的检测方法 | |
CN107506513A (zh) | 一种基于a*算法的兵棋六角格地图路径规划方法 | |
CN105678844A (zh) | 一种基于地物散乱点逐点增加的轮廓构建算法 | |
CN103916953A (zh) | 目标定位的方法、系统及检测节点 | |
Courtat et al. | Stochastic simulation of urban environments. application to path-loss in wireless systems | |
Vuckovik et al. | Durkin’s propagation model based on triangular irregular network terrain |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |