CN113190794A - 一种新的数据空间离散化算法 - Google Patents

一种新的数据空间离散化算法 Download PDF

Info

Publication number
CN113190794A
CN113190794A CN202110176659.2A CN202110176659A CN113190794A CN 113190794 A CN113190794 A CN 113190794A CN 202110176659 A CN202110176659 A CN 202110176659A CN 113190794 A CN113190794 A CN 113190794A
Authority
CN
China
Prior art keywords
currnode
inval
interval
root
cdcuts
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.)
Pending
Application number
CN202110176659.2A
Other languages
English (en)
Inventor
赵银军
李家旭
丁爱中
邓齐玉
童凯
莫莉萍
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangxi Zhongqingtai Environmental Technology Co ltd
Nanning Normal University
Original Assignee
Guangxi Zhongqingtai Environmental Technology Co ltd
Nanning Normal University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangxi Zhongqingtai Environmental Technology Co ltd, Nanning Normal University filed Critical Guangxi Zhongqingtai Environmental Technology Co ltd
Priority to CN202110176659.2A priority Critical patent/CN113190794A/zh
Publication of CN113190794A publication Critical patent/CN113190794A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0635—Risk analysis of enterprise or organisation activities

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • Computational Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Operations Research (AREA)
  • Strategic Management (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Databases & Information Systems (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Probability & Statistics with Applications (AREA)
  • Evolutionary Biology (AREA)
  • Development Economics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Educational Administration (AREA)
  • Software Systems (AREA)
  • Game Theory and Decision Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Marketing (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Instructional Devices (AREA)

Abstract

本发明涉及空间数据转换技术领域，具体为一种新的数据空间离散化算法，包括初始化候选断点、确定最优候选断点集、确定最优区间数、连续属性离散化等步骤。本算法通过将空间数据的空间特征引入到空间数据的离散化过程中，且兼顾了空间数据的空间特征和属性特征，与现有的离散化算法相比，是一种专门针对地理探测器模型并且充分考虑空间特征的数据离散化算法，弥补传统数据离散化算法的不足且大幅提升了地理探测器q值的精准度，为科研人员准确找出哪些风险因子对实验结果的影响更大提供了有利条件。

Description

一种新的数据空间离散化算法
技术领域
本发明涉及空间数据转换技术领域，具体为一种新的数据空间离散化算法。
背景技术
地理探测器模型是基于空间分层异质性原理来揭示空间分异性背后的驱动因子。现在许多学者将该模型应用于生态环境、区域经济、土地利用、旅游、气象等领域的研究，然而这些领域还有许多如环境因子、社会经济因子等连续型因子数据，例如温度、降水、植被覆盖度、高程、坡度、人口、旅游业总收入等，而地理探测器模型擅长处理的是离散类型因子数据，例如土壤类型、土地利用类型、植被类型、流域类型等，所以连续型因子数据作为空间数据就需要通过数据离散化算法转换为离散型空间数据之后才能导入地理探测器模型进行分析。
现在常用的数据离散化算法有自然断点法、等间距法、几何断点法、分位数法、标准差距离法等，这些算法虽然也能够用于空间数据属性特征的离散化，但均未考虑空间数据的空间特征以及空间实体之间的相关性，空间实体所形成格局的空间分异性，而且各种算法也有各自最适合应用的技术领域，比如分位数法适用于样本容量较大且样本呈线性分布的情况，目前来讲医学领域应用较多，几何断点法是专门为制图而设计的一种数据离散化算法，等间距法对于服从正态分布的数据，该算法效果较好。
目前针对地理探测器模型并且考虑空间特征的数据离散化算法的研究很少，因此从空间数据表现出的空间特征出发，专门为地理探测器模型设计一种空间数据离散化算法，一是能有效弥补传统数据离散化算法的不足，二是能达到提升地理探测器模型精度并降低其应用过程中的不确定性的目的。
发明内容
本发明一种新的数据空间离散化算法，包括初始化候选断点、确定最优候选断点集、确定最优区间数、连续属性离散化等步骤，其将空间数据的空间特征引入到空间数据的离散化过程中，同时兼顾了空间数据的空间特征和属性特征，与传统的离散化算法相比，是一种专门针对地理探测器模型并且充分考虑空间特征的数据离散化算法，弥补传统数据离散化算法的不足且提升了地理探测器q值的精准度，有效解决了上述问题。
本发明采用的技术方案：
1、一种新的数据空间离散化算法，其特征在于，包括以下步骤：
(Ⅰ)初始化候选断点：输入已经处理好的空间数据，指定需要离散化的连续属性，对连续属性的所有取值进行数据预处理，得到候选断点集合；
(Ⅱ)确定最优候选断点集：包括候选区间分割和搜索莫兰指数最优断点两个步骤，根据叠置断点法，用候选断点将连续属性取值集合分割为左右两个区间，分别计算分割后左右两个区间的平均值，再用左右区间各自的平均值代替各自区间的所有属性的取值，则空间数据的连续属性取值集合只包含左右区间的两个平均值，然后计算集合的全局莫兰指数，取莫兰指数最大值对应的分割点作为最优离散化断点，再将分割后的左右区间重复上述过程直至满足叠置断点法的终止条件，最后收集所有最优候选断点，得到最优候选断点集合；
(Ⅲ)确定最优区间数：对于最优候选断点集合，经增量组合最优断点、区间数PD值序列匹配、滑动计算区间比值后可确定最优区间数，即根据每个断点分割的区间长度进行降序排序，将集合中的第一个点作为起点，每次从起点开始增量截取断点，直到从起点截取至最后一个断点，得到与最优断点集合的长度相同的断点组合的集合，计算集合中每个断点组合的因子探测值，按照断点组合长度的降序顺序构造与其对应的因子探测值的序列，计算因子探测值与区间长度的比值，由前向后查找，如果遇到后面比值小于前面比值的，则前面位置对应的组合长度值加一为最优区间数，对应的断点组合为最优断点集合；
(Ⅳ)连续属性离散化：对于步骤(Ⅲ)得到的最优断点集合，将指定的连续属性按照断点进行划分并检查区间是否正常，进而得到带有离散属性的空间数据。
叠置断点法：假设有一个方形的地理空间，里面均匀分布的地理单元既包含自变量的属性值又包含因变量的属性值，对所有地理单元的自变量属性的唯一值进行枚举，如果以某一个唯一值作为断点，自变量属性值小于该断点的地理单元聚集在地理空间的左半部分或右半部分，占空间面积一半的区域，称为1号区域，大于该断点的地理单元聚集在另外一个区域，称为2号区域。如果1号区域内地理单元的因变量同样呈现高度的聚集性，同时2号区域内的因变量也呈现出高度的聚集性，通过计算1号区域内因变量和自变量的莫兰指数的差值为0，2号区域同样的差值也为0，则说明地理单元的因变量和自变量具有高度的空间相关性，则该断点为当前地理空间的最优分割点。通过计算划分出的区域内自变量属性值的莫兰指数，如果小于1 则可以继续划分，否则划分结束。由于通过因变量和自变量进行叠置来获取最优断点，所以这种算法命名为叠置断点法，英文全称为 Overlay Break，简称OB算法。当然这是最理想的状态，实际操作中，地理空间并非是完全符合规则的方形，但OB算法可以通过比较1号区域和2号区域的莫兰指数差值的最小值，接近于理想状态差值的程度来确定最优断点，进而使被划分的地理空间尽可能接近于上述理想状态。
全局莫兰指数作为空间自相关性度量的常用指标，它度量地理空间内某个单元的某种属性值或地理现象，与其空间邻近单元之间的相关程度，其既能处理空间数据的属性值，也能处理空间数据的空间位置，引自《Moran,P.A.P.1950.Notes on continuousstochastic phenomena[J].Biometrika,37(1-2):17-23.》。
进一步的，所述初始化候选断点的具体步骤如下：
(1)首先创建一个存储所有候选断点的数组CdCuts，对于指定需要进行离散化的属性X，读取它的每个属性值Xi，存入一个临时数组 templist；
(2)取templist中所有数值的唯一值即全部的候选断点，每个候选断点值为pi，总个数为u，进行升序排序，将排序结果写入 CdCuts,CdCuts＝{pi|i∈[1,u]}；
(3)初始化一棵空二叉树DivTree，用于离散化过程中收集和组织最优断点集合O；
(4)创建被分割的区间inval，所述inval的初始值为[X1,Xn]，inval左端点L,右端点R，边界点B，区间分割阈值t，其中L＝min(CdCuts)， R＝max(CdCuts)，t＝1，B为空值。
进一步的，所述确定最优CdCuts候选断点集的具体步骤如下：
(1)如果CdCuts中的候选断点个数i大于1则开始遍历CdCuts，取出候选断点pi，判断DivTree中根节点值ROOT是否为空，如果为空，则需要分割的区间inval＝[L,R]，如果不为空，则根据公式1.1确定inval 的端点，
Figure RE-GDA0003047282520000051
如果为空则表明还没有产生最优断点，公式1.1为数学领域常见的分段函数。
(2)计算空间数据图层S中的地理要素属性X和Y，其中X左侧区间invalL的莫兰指数为Glmi(invalL,X)，Y左侧区间invalL的莫兰指数为 Glmi(invalL,X)，根据公式1.2计算左侧区间invalL的地理要素属性X与 Y的差分值dfMiL，同理，X右侧区间invalR的莫兰指数为Glmi(invalR,X)， Y右侧区间invalR的莫兰指数为Glmi(invalR,Y)，根据公式1.3计算右侧区间invalR的地理要素属性X与Y的差分值dfMiR，
dfMiL＝Glmi(invalL,X)-Glmi(invalL,Y) (1.2)，
dfMiR＝Glmi(invalR,X)-Glmi(invalR,Y) (1.3)；
(3)根据公式1.4计算invalL和invalR之间的最小差分值mdfMi，对于从CdCuts中遍历出的任意pi点都对应一个最小差分值mdfMi，
mdfMi＝min(dfMiL,dfMiR) (1.4)；
(4)创建一个数组plist，遍历数组CdCuts的过程中，取出第一个p1点对应的mdfMi作为最小的mdfMimin，将p1点存入数组plist，记p1点分割的区间长度为l'，后面遍历出的pi点对应的mdfMi如果小于mdfMimin，则令mdfMimin＝mdfmi，将pi点存入数组plist，数组plist中的最后一个元素为最优分割点Oi，第一个至倒数第二个元素均为次最优分割点Oi'；
其中数组plist用于存储遍历过程中产生的所有mdfMi的低值点。
(5)如果DivTree的ROOT为空，则直接创建CurrNode得到Oi和对应的l'，然后写入CurrNode并压入DivTree中，此时DivTree的CurrNode根节点值ROOT＝Oi，如果DivTree的ROOT不为空，则先检索以下两个条件：第一、DivTree中是否已经存在Oi，第二、Oi是否为CdCuts中的最后一个候选断点，如果两个条件同时不成立则创建CurrNode，将得到的Oi和对应的l'写入CurrNode中，根据Oi数值与DivTree中已有结点的 root的大小进行比较，然后在相应的位置插入CurrNode，令边界点 B＝Oi，如果两个判断条件只成立其中一个，则表示当前Oi点不符合要求，需要查找次最优分割点O’i，具体过程为：在数组plist中删除掉最后一个元素plist[-1]，原来的倒数第二个元素plist[-2]则为当前的次最优分割点O’i，然后再创建CurrNode，并将得到O’i和对应的l'写入 CurrNode并压入DivTree中，一直到能够压入DivTree为止，则表示找到了符合条件的O’i，此时令边界点B＝Oi，如果O’i仍然无法压入DivTree 中，则返回重新查找次最优分割点O’i，如果一直执行到数组plist中的元素为空，仍然没有找到符合条件的O’i，则说明当前的区间inval没有最优分割点，即区间inval无法再继续分割，此时将表示区间inval无法继续分割的标志NondivFlag赋值为True；
如果两个判断条件只成立其中一个，则表示当前的Oi为重复断点或者Oi到达区间inval的边界。
(6)对于区间inval，首先判断该区间是否同时满足可分割的三个条件：第一个条件Glmi(inval,X)小于t、第二个条件NondivFlag＝False、第三个条件区间inval内的地理要素大于30个，如果不满足则执行步骤(7)，如果同时满足上述三个条件则继续对区间inval进行划分，通过公式1.5计算CdCutsL，令CdCuts＝CdCutsL，再进入区间inval的一次递归过程，一次递归的层数lev初始值为0，递归开始：执行步骤(1) -步骤(6)，一次递归完成后继续判断区间inval是否同时满足上述可分割的三个条件，如果不满足则执行步骤(7)，如果同时满足则 lev＝lev+1进入二次递归，递归开始：执行步骤(1)-步骤(6)，二次递归结束后，确定最优候选断点集的步骤已全部完成；
待步骤(7)执行完成后通过公式1.1确定区间inval，通过公式 1.6计算CdCutsR，令CdCuts＝CdCutsR，进入区间inval的递归过程，当前递归的层数为lev，递归开始：执行步骤(1)-步骤(6)，递归完成后继续判断区间inval是否同时满足上述可分割的三个条件，如果同时满足则lev＝lev+1，进入区间inval的递归过程，递归开始：执行步骤(1)-步骤(6)，递归结束后确定最优候选断点集的步骤已全部完成；
如果区间inval不满足上述可分割的三个条件，则第lev层inval的递归到达出口，然后执行步骤(8)，待步骤(8)完成后如果此时lev＝0 则确定最优候选断点集的步骤已全部完成，否则lev＝lev-1，当前的递归层数又回到lev，然后根据公式1.1确定区间inval，如果 CdCuts＝CdCutsL则执行步骤(9)，待步骤(9)完成后进入区间inval的递归过程，递归开始：执行步骤(1)-步骤(6)，递归结束后，确定最优候选断点集的步骤已全部完成；
如果CdCuts＝CdCutsR则执行步骤(10)，待步骤(10)完成后当前的递归层数还是回到lev，则lev＝lev-1，直到区间inval的递归过程回到lev＝0，确定最优候选断点集的步骤已全部完成，
Figure RE-GDA0003047282520000081
Figure RE-GDA0003047282520000082
公式1.5和公式1.6均为分段函数。
(7)对于inval，如果CurrNode的root等于max(CdCutsL)并且 CurrNode的parentNode不为空，则令L等于CurrNode的root，令CurrNode 指向CurrNode的parentNode，再令B等于CurrNode的root，如果 CurrNode的root等于max(CdCutsR)并且CurrNode的parentNode不为空，令L等于CurrNode的root，令CurrNode指向CurrNode的parentNode，再令B等于CurrNode的root，如果CurrNode的parentNode为空，则说明只产生了一个Oi并且inval的Glmi(inval,X)大于t，inval无法进行继续分割；
(8)对于inval，如果CurrNode的root等于max(CdCutsL)并且 CurrNode的parentNode不为空，令CurrNode指向CurrNode的parentNode，令L等于CurrNode的root，再令B等于CurrNode的parentNode包含的 root，如果CurrNode的root等于max(CdCutsR)并且CurrNode的 parentNode不为空，令L等于CurrNode的root，令CurrNode指向CurrNode 的parentNode，再令B等于CurrNode的root，如果CurrNode的parentNode 为空，则说明只产生了一个Oi并且inval的Glmi(inval,X)大于t，inval无法进行继续分割；
(9)对于inval，同时满足以下两个条件NondivFlag＝True和 CurrNode的root＝max(CdCutsL)，令L等于CurrNode的root，如果 CurrNode的parentNode不为空，则令CurrNode指向CurrNode的 parentNode，再令B等于CurrNode的root，NondivFlag＝False，如果 CurrNode的parentNode为空，令B等于CurrNode的root， NondivFlag＝False，如果NondivFlag＝False并且CurrNode的parentNode 为空，则令L和B都等于CurrNode的root，如果NondivFlag＝False并且 CurrNode的parentNode不为空，判断invalL是否可分，如果可分，令L 等于CurrNode的root，CurrNode指向CurrNode的parentNode，再令B等于CurrNode的root，如果不可分，CurrNode指向CurrNode的parentNode，L等于CurrNode的root，B等于CurrNode的parentNode包含的root；
(10)对于inval，如果满足NondivFlag＝True、CurrNode的parentNode 不为空、CurrNode的root＝max(CdCutsL)这三个条件，则令CurrNode指向 CurrNode的parentNode，L等于CurrNode的root，再令B等于CurrNode的 parentNode包含的root＝False和NondivFlag＝False，如果CurrNode的 parentNode为空，则令NondivFlag＝False。
在确定最优断点集合的过程中，本方案采用递归算法作为算法的主体结构。
进一步的，所述确定最优区间数的具体步骤如下：
(1)首先通过对DivTree进行中根周游，返回其中保存的最优断点集合O以及每个Oi对应的分割区间长度l'；
(2)按照区间长度l'对O降序排序，使排序后每个Oi仍然对应它分割的区间长度l'；
(3)增量截取O中的元素，每次截取的断点数量比前一次多1个，得到{Oi}，{Oi,Oj}，{Oi,Oj,Ok}，{Oi,Oj,Ok,…ON}，并对每个子集合O*进行升序排序，假设
Figure RE-GDA0003047282520000101
为X的子集，元素个数N1,NN+1,NN+1，
Figure RE-GDA0003047282520000102
Figure RE-GDA0003047282520000103
对于图层S中地理要素属性X的取值处于
Figure RE-GDA0003047282520000104
的要素，它们的Y属性的取值集合为
Figure RE-GDA0003047282520000105
方差分别为
Figure RE-GDA0003047282520000106
根据公式1.9计算每个O*的q值，其中n'表示O中的断点对图层S离散化形成的子区域内的要素个数，var 表示子区域内要素的Y取值的方差，
Figure RE-GDA0003047282520000111
Figure RE-GDA0003047282520000112
Figure RE-GDA0003047282520000113
公式1.7和公式1.8为分段函数，公式1.9引自《Moran,P.A. P.1950.Notes oncontinuous stochastic phenomena[J]. Biometrika,37(1-2):17-23》。
(4)Q为q的集合，对Q降序排序，根据公式2.0计算Q中每个元素的qr值，QR为qr的集合，N*表示O*中断点的数目，
qr＝q/N* (2.0)；
(5)通过三点滑动查找QR中所有元素的第一个峰值，如果 qri-1＜qri＜qri+1，则令i＝i+1，如果qri＞qri+1&qri＞qri-1，则确定最优区间数的步骤已完成，qri对应的(BKS'_Len+1)为最优区间数，如果 qri＜qri+1&qri＜qri-1，则qri-1对应的(BKS'_Len+1)为最优区间数。
进一步的，所述连续属性离散化的具体步骤如下：
(1)运用经上述步骤得到的断点集合O对空间数据图层S的属性X 进行划分，首先遍历X，取出Xi；
(2)对于任意的Xi，遍历O，根据公式2.1计算Xj所属的区间编号值VN，
Figure RE-GDA0003047282520000121
公式2.1为分段函数。
(3)把VN写入到空间数据图层S的D属性中；
(4)检查属性D的取值集合中是否存在取值个数小于2的编号值，如果存在将该编号值合并到相邻区间中。
进一步的，所述空间数据为点状矢量数据或面状矢量数据。
进一步的，所述空间数据包含的每个点要素或面要素，同时具有因变量和自变量的属性特征。
综上所述，由于采用了上述技术方案，本发明的有益效果：与传统的离散化算法相比，本申请所提出的一种新的数据空间离散化算法是将空间数据的空间特征引入到空间数据的离散化过程中，同时兼顾了空间数据的空间特征和属性特征，该算法的易实现程度与传统离散化算法相当，但传统的离散化算法需要事先指定离散化的区间数，而本申请提出的算法无需事先指定，通过引入地理探测器q值的优化区间断点组合即可自动计算出最优区间断点，大幅降低了使用地理探测器时空间数据在离散化过程中的不确定性，因此计算结果即地理探测器q值的精准度更高，而且计算更加高效、便捷，是一种专门针对地理探测器模型的数据离散化算法，可在相关领域大规模推广使用。
附图说明
为了更清楚地说明本发明的实例或现有技术中的技术方案，下面将对实施实例或现有技术描述中所需要的附图做简单地介绍，显然，下面描述中的附图仅仅是本发明的一些实例，对于本领域普通技术人员来说，在不付出创造性的前提下，还可以根据这些附图获得其他的附图。
图1一种有效提高地理探测器模型精度的新算法的流程框图；
图2山西和顺县地理位置图；
图3神经管畸形病发病率与本算法得出的高程-坡度图；
图4另外四种常见离散化算法得出的高程图；
图5另外四种常见离散化算法得出的坡度图；
图6坡度与高程根据不同离散化算法得出的q值对比图；
图7坡度与高程在相同区间数根据不同离散化算法得出的q值对比图。
具体实施方式
下面将结合本发明实例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例：
一种新的数据空间离散化算法，如图1所示，包括以下步骤：
(Ⅰ)初始化候选断点：输入已经处理好的空间数据，指定需要离散化的连续属性，对连续属性的所有取值进行数据预处理，得到候选断点集合，具体步骤如下：(1)首先创建一个存储所有候选断点的数组CdCuts，对于指定需要进行离散化的属性X，读取它的每个属性值 Xi，存入一个临时数组templist，(2)取templist中所有数值的唯一值即全部的候选断点，每个候选断点值为pi，总个数为u，进行升序排序，将排序结果写入CdCuts,CdCuts＝{pi|i∈[1,u]}，(3)初始化一棵空二叉树 DivTree，用于离散化过程中收集和组织最优断点集合O，(4)创建被分割的区间inval，所述inval的初始值为[X1,Xn]，inval左端点L,右端点R，边界点B，区间分割阈值t，其中L＝min(CdCuts)，R＝max(CdCuts)， t＝1，B为空值；其中，所述空间数据为点状矢量数据或面状矢量数据，所述空间数据包含的每个点要素或面要素，同时具有因变量和自变量的属性特征；
(Ⅱ)确定最优候选断点集：包括候选区间分割和搜索莫兰指数最优断点两个步骤，根据叠置断点法，用候选断点将连续属性取值集合分割为左右两个区间，分别计算分割后左右两个区间的平均值，再用左右区间各自的平均值代替各自区间的所有属性的取值，则空间数据的连续属性取值集合只包含左右区间的两个平均值，然后计算集合的全局莫兰指数，取莫兰指数最大值对应的分割点作为最优离散化断点，再将分割后的左右区间重复上述过程直至满足叠置断点法的终止条件，最后收集所有最优候选断点，得到最优候选断点集合，具体步骤如下：
(1)如果CdCuts中的候选断点个数i大于1则开始遍历CdCuts，取出候选断点pi，判断DivTree中根节点值ROOT是否为空，如果为空，表明还没有产生最优断点，则需要分割的区间inval＝[L,R]，如果不为空，则根据公式1.1确定inval的端点，
Figure RE-GDA0003047282520000151
(2)计算空间数据图层S中的地理要素属性X和Y，其中X左侧区间invalL的莫兰指数为Glmi(invalL,X)，Y左侧区间invalL的莫兰指数为Glmi(invalL,X)，根据公式1.2计算左侧区间invalL的地理要素属性X与 Y的差分值dfMiL，同理，X右侧区间invalR的莫兰指数为Glmi(invalR,X)， Y右侧区间invalR的莫兰指数为Glmi(invalR,Y)，根据公式1.3计算右侧区间invalR的地理要素属性X与Y的差分值dfMiR，
dfMiL＝Glmi(invalL,X)-Glmi(invalL,Y) (1.2)，
dfMiR＝Glmi(invalR,X)-Glmi(invalR,Y) (1.3)；
(3)根据公式1.4计算invalL和invalR之间的最小差分值mdfMi，对于从CdCuts中遍历出的任意pi点都对应一个最小差分值mdfMi，
mdfMi＝min(dfMiL,dfMiR) (1.4)；
(4)创建一个数组plist用于存储遍历过程中产生的所有mdfMi的低值点遍历数组CdCuts的过程中，取出第一个p1点对应的的mdfMi作为最小的mdfMimin，将p1点存入数组plist，记p1点分割的区间长度为l'，后面遍历出的pi点对应的mdfMi如果小于mdfMimin，则令mdfMimin＝mdfmi，将pi点存入数组plist，数组plist中的最后一个元素为最优分割点Oi，第一个至倒数第二个元素均为次最优分割点Oi'；
(5)如果DivTree的ROOT为空，则直接创建CurrNode得到Oi和对应的l'，然后写入CurrNode并压入DivTree中，此时DivTree的CurrNode根节点值ROOT＝Oi，如果DivTree的ROOT不为空，则先检索以下两个条件：第一、DivTree中是否已经存在Oi，第二、Oi是否为CdCuts中的最后一个候选断点，如果两个条件同时不成立则创建CurrNode，将得到的Oi和对应的l'写入CurrNode中，根据Oi数值与DivTree中已有结点的 root的大小进行比较，然后在相应的位置插入CurrNode，令边界点 B＝Oi，如果两个判断条件只成立其中一个，则表示当前的Oi为重复断点或Oi到达区间inval的边界，即表示当前Oi点不符合要求，需要查找次最优分割点O’i，具体过程为：在数组plist中删除掉最后一个元素 plist[-1]，原来的倒数第二个元素plist[-2]则为当前的次最优分割点O’i，然后再创建CurrNode，并将得到O’i和对应的l'写入CurrNode并压入 DivTree中，一直到能够压入DivTree为止，则表示找到了符合条件的O’i，此时令边界点B＝Oi，如果Oi'仍然无法压入DivTree中，则返回重新查找次最优分割点O’i，如果一直执行到数组plist中的元素为空，仍然没有找到符合条件的Oi'，则说明当前的区间inval没有最优分割点，即区间inval无法再继续分割，此时将表示区间inval无法继续分割的标志 NondivFlag赋值为True；
(6)对于区间inval，首先判断该区间是否同时满足可分割的三个条件：第一个条件Glmi(inval,X)小于t、第二个条件NondivFlag＝False、第三个条件区间inval内的地理要素大于30个，如果不满足则执行步骤(7)，如果同时满足上述三个条件则继续对区间inval进行划分，通过公式1.5计算CdCutsL，令CdCuts＝CdCutsL，再进入区间inval的一次递归过程，一次递归的层数lev初始值为0，递归开始：执行步骤(1) -步骤(6)，一次递归完成后继续判断区间inval是否同时满足上述可分割的三个条件，如果不满足则执行步骤(7)，如果同时满足则 lev＝lev+1进入二次递归，递归开始：执行步骤(1)-步骤(6)，二次递归结束后，确定最优候选断点集的步骤已全部完成；
待步骤(7)执行完成后通过公式1.1确定区间inval，通过公式 1.6计算CdCutsR，令CdCuts＝CdCutsR，进入区间inval的递归过程，当前递归的层数为lev，递归开始：执行步骤(1)-步骤(6)，递归完成后继续判断区间inval是否同时满足上述可分割的三个条件，如果同时满足则lev＝lev+1，进入区间inval的递归过程，递归开始：执行步骤(1)-步骤(6)，递归结束后确定最优候选断点集的步骤已全部完成；
如果区间inval不满足上述可分割的三个条件，则第lev层inval的递归到达出口，然后执行步骤(8)，待步骤(8)完成后如果此时lev＝0 则确定最优候选断点集的步骤已全部完成，否则lev＝lev-1，当前的递归层数又回到lev，然后根据公式1.1确定区间inval，如果 CdCuts＝CdCutsL则执行步骤(9)，待步骤(9)完成后进入区间inval的递归过程，递归开始：执行步骤(1)-步骤(6)，递归结束后，确定最优候选断点集的步骤已全部完成；
如果CdCuts＝CdCutsR则执行步骤(10)，待步骤(10)完成后当前的递归层数还是回到lev，则lev＝lev-1，直到区间inval的递归过程回到lev＝0，确定最优候选断点集的步骤已全部完成，
Figure RE-GDA0003047282520000171
Figure RE-GDA0003047282520000181
(7)对于inval，如果CurrNode的root等于max(CdCutsL)并且CurrNode 的parentNode不为空，则令L等于CurrNode的root，令CurrNode指向CurrNode的parentNode，再令B等于CurrNode的root，如果CurrNode的 root等于max(CdCutsR)并且CurrNode的parentNode不为空，令L等于 CurrNode的root，令CurrNode指向CurrNode的parentNode，再令B等于 CurrNode的root，如果CurrNode的parentNode为空，则说明只产生了一个Oi并且inval的Glmi(inval,X)大于t，inval无法进行继续分割；
(8)对于inval，如果CurrNode的root等于max(CdCutsL)并且CurrNode 的parentNode不为空，令CurrNode指向CurrNode的parentNode，令L等于 CurrNode的root，再令B等于CurrNode的parentNode包含的root，如果 CurrNode的root等于max(CdCutsR)并且CurrNode的parentNode不为空，令 L等于CurrNode的root，令CurrNode指向CurrNode的parentNode，再令B 等于CurrNode的root，如果CurrNode的parentNode为空，则说明只产生了一个Oi并且inval的Glmi(inval,X)大于t，inval无法进行继续分割；
(9)对于inval，同时满足以下两个条件NondivFlag＝True和CurrNode 的root＝max(CdCutsL)，令L等于CurrNode的root，如果CurrNode的 parentNode不为空，则令CurrNode指向CurrNode的parentNode，再令B等于CurrNode的root，NondivFlag＝False，如果CurrNode的parentNode为空，令B等于CurrNode的root，NondivFlag＝False，如果NondivFlag＝False并且CurrNode的parentNode为空，则令L和B都等于CurrNode的root，如果NondivFlag＝False并且CurrNode的parentNode不为空，判断invalL是否可分，如果可分，令L等于CurrNode的root，CurrNode指向CurrNode的 parentNode，再令B等于CurrNode的root，如果不可分，CurrNode指向 CurrNode的parentNode，L等于CurrNode的root，B等于CurrNode的 parentNode包含的root；
(10)对于inval，如果满足NondivFlag＝True、CurrNode的parentNode 不为空、CurrNode的root＝max(CdCutsL)这三个条件，则令CurrNode指向 CurrNode的parentNode，L等于CurrNode的root，再令B等于CurrNode的 parentNode包含的root＝False和NondivFlag＝False，如果CurrNode的 parentNode为空，则令NondivFlag＝False。
(Ⅲ)确定最优区间数：对于最优候选断点集合，经增量组合最优断点、区间数PD值序列匹配、滑动计算区间比值后可确定最优区间数，即根据每个断点分割的区间长度进行降序排序，将集合中的第一个点作为起点，每次从起点开始增量截取断点，直到从起点截取至最后一个断点，得到与最优断点集合的长度相同的断点组合的集合，计算集合中每个断点组合的因子探测值，按照断点组合长度的降序顺序构造与其对应的因子探测值的序列，计算因子探测值与区间长度的比值，由前向后查找，如果遇到后面比值小于前面比值的，则前面位置对应的组合长度值加一为最优区间数，对应的断点组合为最优断点集合，具体步骤如下：
(1)首先通过对DivTree进行中根周游，返回其中保存的最优断点集合O以及每个Oi对应的分割区间长度l'；
(2)按照区间长度l'对O降序排序，使排序后每个Oi仍然对应它分割的区间长度l'；
(3)增量截取O中的元素，每次截取的断点数量比前一次多1个，得到{Oi}，{Oi,Oj}，{Oi,Oj,Ok}，{Oi,Oj,Ok,…ON}，并对每个子集合O*进行升序排序，假设
Figure RE-GDA0003047282520000201
为X的子集，元素个数N1,NN+1,NN+1，
Figure RE-GDA0003047282520000202
Figure RE-GDA0003047282520000203
对于图层S中地理要素属性X的取值处于
Figure RE-GDA0003047282520000204
的要素，它们的Y属性的取值集合为
Figure RE-GDA0003047282520000205
方差分别为
Figure RE-GDA0003047282520000206
根据公式1.9计算每个O*的q值，其中n'表示O中的断点对图层S离散化形成的子区域内的要素个数，var 表示子区域内要素的Y取值的方差，
Figure RE-GDA0003047282520000207
Figure RE-GDA0003047282520000208
Figure RE-GDA0003047282520000209
(4)Q为q的集合，对Q降序排序，根据公式2.0计算Q中每个元素的qr值，QR为qr的集合，N*表示O*中断点的数目，
qr＝q/N* (2.0)；
通过三点滑动查找QR中所有元素的第一个峰值，如果 qri-1＜qri＜qri+1，则令i＝i+1，如果qri＞qri+1&qri＞qri-1，则确定最优区间数的步骤已完成，qri对应的(BKS'_Len+1)为最优区间数，如果 qri＜qri+1&qri＜qri-1，则qri-1对应的(BKS'_Len+1)为最优区间数。
(Ⅳ)连续属性离散化：对于步骤(Ⅲ)得到的最优断点集合，将指定的连续属性按照断点进行划分并检查区间是否正常，进而得到带有离散属性的空间数据，具体步骤如下：首先运用经上述步骤得到的断点集合O对空间数据图层S的属性X进行划分，首先遍历X，取出 Xi，对于任意的Xi，遍历O，然后根据下列公式计算Xj所属的区间编号值VN，
Figure RE-GDA0003047282520000211
把VN写入到空间数据图层S的D属性中，最后检查属性D的取值集合中是否存在取值个数小于2的编号值，如果存在则将该编号值合并到相邻区间中。
神经管出生缺陷(NTD)作为一种出生缺陷疾病，是导致婴儿死亡和残疾的主要原因之一。生理、遗传倾向、自然环境、社会环境都是导致神经管出生缺陷疾病的原因，不同地区由生理、遗传倾向导致的出生缺陷率相差不大，所以自然环境、社会环境可能是影响区域之间神经管出生缺陷疾病发生差异的主要因素，对这些致病的潜在因素进行量化分析，阐明各个因素对疾病发病率的贡献程度，对于神经管出生缺陷的防治具有重要意义。当前地理探测器的分析结果表明，地质岩性、流域类型、土壤类型、高程、坡度等因子均对NTD产生一定的潜在影响，是疾病发生的环境风险因子，而每个环境风险因子都有一个q值，q值为地理探测器的常用函数值。简单点讲，就是通过计算风险因子在地理探测器中的q值，将不同环境风险因子的q值进行比较，q值大的风险因子说明其对导致神经管出生缺陷疾病的影响更大，即找出了哪些环境风险因子是导致神经管出生缺陷疾病的主因，因此q值的精准度就显得至关重要。
本实施例以山西省和顺县为研究对象，该县是世界上NTD高发的地区之一，风险因子选择的是高程和坡度，两者均为连续型风险因子，所述高程和坡度的数据均通过DEM数字高程模型获取。如图2 所示，和顺县位于山西省境内太行山的西侧，地貌主要为山地和丘陵，相对高差在300-500米，平均海拔高度1300米，包含326个行政村，由于难以获得行政村的具体边界，研究人员以泰森多边形作为研究区域各个行政村的边界。
图3中a表示的是山西省和顺县各个区域每1000人当中NTD的发病率，b表示的是连续型风险因子高程经本算法计算后的离散化结果，c表示的是连续型风险因子坡度经本算法计算后的离散化结果。
图4表示的是连续型风险因子高程经另外4种传统离散化算法计算后的结果，图5表示的是连续型风险因子坡度经另外4种传统离散化算法计算后的结果，4种传统离散化算法分别为a几何断点法，b 分位数法，c自然断点法，d等间距法。
将图3的b图与图4的四幅图做对比，明显看到其精度更高，将图3的c图与图5的四幅图做对比，也是c图的精度更高。
图6表示的是风险因子坡度和高程经不同的离散化算法计算后得出的地理探测器q值，由于OB算法无需区分区间数，因此另外4 种算法选择的区间数为2-8。从图中可知坡度的计算中，OB算法计算出的地理探测器q值为0.131，而另外4种传统离散化算法在区间数2-8的计算结果中，仅有两个结果略高于0.131，有两个结果也是 0.131，其余的均小于0.131，高程的计算中，OB算法计算出的地理探测器q值为0.108，而另外4种传统离散化算法在区间数2-8的计算结果全部小于0.108，由此可见OB算法计算出的地理探测器q值相比于另外4种传统离散化算法，精准度明显提高。
图7表示的是风险因子坡度和高程在相同区间数的情况下，经不同的离散化算法计算后得出的地理探测器q值，从图中可知，风险因子坡度在区间数10经OB算法计算出的地理探测器q值为0.131，另外4种传统离散化算法中，只有自然断点法计算出的q值0.132略高出0.001，其余的均小于0.131，而风险因子高程在区间数11经OB 算法计算出的地理探测器q值为0.108，而另外4种传统离散化算法计算出的q值全部小于0.108，由此可见OB算法计算出的地理探测器q值相比于另外4种传统离散化算法，精准度明显提高。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在发明的保护范围之内。

Claims (6)

1.一种新的数据空间离散化算法，其特征在于，包括以下步骤：
(Ⅰ)初始化候选断点：输入已经处理好的空间数据，指定需要离散化的连续属性，对连续属性的所有取值进行数据预处理，得到候选断点集合；
(Ⅱ)确定最优候选断点集：包括候选区间分割和搜索莫兰指数最优断点两个步骤，根据叠置断点法，用候选断点将连续属性取值集合分割为左右两个区间，分别计算分割后左右两个区间的平均值，再用左右区间各自的平均值代替各自区间的所有属性的取值，则空间数据的连续属性取值集合只包含左右区间的两个平均值，然后计算集合的全局莫兰指数，取莫兰指数最大值对应的分割点作为最优离散化断点，再将分割后的左右区间重复上述过程直至满足叠置断点法的终止条件，最后收集所有最优候选断点，得到最优候选断点集合，具体步骤如下：
(1)如果CdCuts中的候选断点个数i大于1则开始遍历CdCuts，取出候选断点pi，判断DivTree中根节点值ROOT是否为空，如果为空，则需要分割的区间inval＝[L,R]，如果不为空，则根据公式1.1确定inval的端点，
Figure RE-FDA0003099018260000011
(2)计算空间数据图层S中的地理要素属性X和Y，其中X左侧区间invalL的莫兰指数为Glmi(invalL,X)，Y左侧区间invalL的莫兰指数为Glmi(invalL,Y)，根据公式1.2计算左侧区间invalL的地理要素属性X与Y的差分值dfMiL，同理，X右侧区间invalR的莫兰指数为Glmi(invalR,X)，Y右侧区间invalR的莫兰指数为Glmi(invalR,Y)，根据公式1.3计算右侧区间invalR的地理要素属性X与Y的差分值dfMiR，
dfMiL＝Glmi(invalL,X)-Glmi(invalL,Y) (1.2)，
dfMiR＝Glmi(invalR,X)-Glmi(invalR,Y) (1.3)；
(3)根据公式1.4计算invalL和invalR之间的最小差分值mdfMi，对于从CdCuts中遍历出的任意pi点都对应一个最小差分值mdfMi，
mdfMi＝min(dfMiL,dfMiR) (1.4)；
(4)创建一个数组plist，遍历数组CdCuts的过程中，取出第一个点p1对应的mdfMi作为最小的mdfMimin，将点p1存入数组plist，记点p1分割的区间长度为l'，后面遍历出的pi点对应的mdfMi如果小于mdfMimin，则令mdfMimin＝mdfMi，将pi点存入数组plist，数组plist中的最后一个元素为最优分割点Oi，第一个至倒数第二个元素均为次最优分割点Oi'；
(5)如果DivTree的ROOT为空，则直接创建CurrNode得到Oi和对应的l'，然后写入CurrNode并压入DivTree中，此时DivTree的CurrNode根节点值ROOT＝Oi，如果DivTree的ROOT不为空，则先检索以下两个条件：第一、DivTree中是否已经存在Oi，第二、Oi是否为CdCuts中的最后一个候选断点，如果两个条件同时不成立则创建CurrNode，将得到的Oi和对应的l'写入CurrNode中，根据Oi数值与DivTree中已有结点的root的大小进行比较，然后在相应的位置插入CurrNode，令边界点B＝Oi，如果两个判断条件只成立其中一个，则表示当前Oi点不符合要求，需要查找次最优分割点O'i，具体过程为：在数组plist中删除掉最后一个元素plist[-1]，原来的倒数第二个元素plist[-2]则为当前的次最优分割点O'i，然后再创建CurrNode，并将得到O'i和对应的l'写入CurrNode并压入DivTree中，一直到能够压入DivTree为止，则表示找到了符合条件的O'i，此时令边界点B＝O'i，如果O'i仍然无法压入DivTree中，则返回重新查找次最优分割点O'i，如果一直执行到数组plist中的元素为空，仍然没有找到符合条件的O'i，则说明当前的区间inval没有最优分割点，即区间inval无法再继续分割，此时将表示区间inval无法继续分割的标志NondivFlag赋值为True；
(6)对于区间inval，首先判断该区间是否同时满足可分割的三个条件：第一个条件Glmi(inval,X)小于t、第二个条件NondivFlag＝False、第三个条件区间inval内的地理要素大于30个，如果不满足则执行步骤(7)，如果同时满足上述三个条件则继续对区间inval进行划分，通过公式1.5计算CdCutsL，令CdCuts＝CdCutsL，再进入区间inval的一次递归过程，一次递归的层数lev初始值为0，递归开始：执行步骤(1)-步骤(6)，一次递归完成后继续判断区间inval是否同时满足上述可分割的三个条件，如果不满足则执行步骤(7)，如果同时满足则lev＝lev+1进入二次递归，递归开始：执行步骤(1)-步骤(6)，二次递归结束后，确定最优候选断点集的步骤已全部完成；
待步骤(7)执行完成后通过公式1.1确定区间inval，通过公式1.6计算CdCutsR，令CdCuts＝CdCutsR，进入区间inval的递归过程，当前递归的层数为lev，递归开始：执行步骤(1)-步骤(6)，递归完成后继续判断区间inval是否同时满足上述可分割的三个条件，如果同时满足则lev＝lev+1，进入区间inval的递归过程，递归开始：执行步骤(1)-步骤(6)，递归结束后确定最优候选断点集的步骤已全部完成；
如果区间inval不满足上述可分割的三个条件，则第lev层inval的递归到达出口，然后执行步骤(8)，待步骤(8)完成后如果此时lev＝0则确定最优候选断点集的步骤已全部完成，否则lev＝lev-1，当前的递归层数又回到lev，然后根据公式1.1确定区间inval，如果CdCuts＝CdCutsL则执行步骤(9)，待步骤(9)完成后进入区间inval的递归过程，递归开始：执行步骤(1)-步骤(6)，递归结束后，确定最优候选断点集的步骤已全部完成；
如果CdCuts＝CdCutsR则执行步骤(10)，待步骤(10)完成后当前的递归层数还是回到lev，则lev＝lev-1，直到区间inval的递归过程回到lev＝0，确定最优候选断点集的步骤已全部完成，
Figure RE-FDA0003099018260000041
Figure RE-FDA0003099018260000042
(7)对于inval，如果CurrNode的root等于max(CdCutsL)并且CurrNode的parentNode不为空，则令L等于CurrNode的root，令CurrNode指向CurrNode的parentNode，再令B等于CurrNode的root，如果CurrNode的root等于max(CdCutsR)并且CurrNode的parentNode不为空，令L等于CurrNode的root，令CurrNode指向CurrNode的parentNode，再令B等于CurrNode的root，如果CurrNode的parentNode为空，则说明只产生了一个Oi并且inval的Glmi(inval,X)大于t，inval无法进行继续分割；
(8)对于inval，如果CurrNode的root等于max(CdCutsL)并且CurrNode的parentNode不为空，令CurrNode指向CurrNode的parentNode，令L等于CurrNode的root，再令B等于CurrNode的parentNode包含的root，如果CurrNode的root等于max(CdCutsR)并且CurrNode的parentNode不为空，令L等于CurrNode的root，令CurrNode指向CurrNode的parentNode，再令B等于CurrNode的root，如果CurrNode的parentNode为空，则说明只产生了一个Oi并且inval的Glmi(inval,X)大于t，区间inval无法再继续分割；
(9)对于inval，同时满足以下两个条件NondivFlag＝True和CurrNode的root＝max(CdCutsL)，令L等于CurrNode的root，如果CurrNode的parentNode不为空，则令CurrNode指向CurrNode的parentNode，再令B等于CurrNode的root，NondivFlag＝False，如果CurrNode的parentNode为空，令B等于CurrNode的root，NondivFlag＝False，如果NondivFlag＝False并且CurrNode的parentNode为空，则令L和B都等于CurrNode的root，如果NondivFlag＝False并且CurrNode的parentNode不为空，判断invalL是否可分，如果可分，令L等于CurrNode的root，CurrNode指向CurrNode的parentNode，再令B等于CurrNode的root，如果不可分，CurrNode指向CurrNode的parentNode，L等于CurrNode的root，B等于CurrNode的parentNode包含的root；
(10)对于inval，如果满足NondivFlag＝True、CurrNode的parentNode不为空、CurrNode的root＝max(CdCutsL)这三个条件，则令CurrNode指向CurrNode的parentNode，L等于CurrNode的root，再令B等于CurrNode的parentNode包含的root＝False和NondivFlag＝False，如果CurrNode的parentNode为空，则令NondivFlag＝False；
(Ⅲ)确定最优区间数：对于最优候选断点集合，经增量组合最优断点、区间数PD值序列匹配、滑动计算区间比值后可确定最优区间数，即根据每个断点分割的区间长度进行降序排序，将集合中的第一个点作为起点，每次从起点开始增量截取断点，直到从起点截取至最后一个断点，得到与最优断点集合的长度相同的断点组合的集合，计算集合中每个断点组合的因子探测值，按照断点组合长度的降序顺序构造与其对应的因子探测值的序列，计算因子探测值与区间长度的比值，由前向后查找，如果遇到后面比值小于前面比值的，则前面位置对应的组合长度值加一为最优区间数，对应的断点组合为最优断点集合；
(Ⅳ)连续属性离散化：对于步骤(Ⅲ)得到的最优断点集合，将指定的连续属性按照断点进行划分并检查区间是否正常，进而得到带有离散属性的空间数据。
2.根据权利要求1所述的一种新的数据空间离散化算法，其特征在于，所述初始化候选断点的具体步骤如下：
(1)创建一个存储所有候选断点的数组CdCuts，对于指定需要进行离散化的属性X，读取它的每个属性值Xi，存入一个临时数组templist；
(2)取templist中所有数值的唯一值即全部的候选断点，每个候选断点值为pi，总个数为u，进行升序排序，将排序结果写入CdCuts,CdCuts＝{pi|i∈[1，u]}；
(3)初始化一棵空二叉树DivTree，用于离散化过程中收集和组织最优断点集合O；
(4)创建被分割的区间inval，所述inval的初始值为[X1，Xn]，inval左端点L,右端点R，边界点B，区间分割阈值t，其中L＝min(CdCuts)，R＝max(CdCuts)，t＝1，B为空值。
3.根据权利要求1所述的一种新的数据空间离散化算法，其特征在于，确定最优区间数的具体步骤如下：
(1)通过对DivTree进行中根周游，返回其中保存的最优断点集合O以及每个Oi对应的分割区间长度l'；
(2)按照区间长度l'对O降序排序，使排序后每个Oi仍然对应它分割的区间长度l'；
(3)增量截取O中的元素，每次截取的断点数量比前一次多1个，得到{Oi}，{Oi，Oj}，{Oi，Oj，Ok}，{Oi，Oj，Ok，···ON}，并对每个子集合O*进行升序排序，假设
Figure RE-FDA0003099018260000071
为X的子集，元素个数N1，Ni+1，NN+1，
Figure RE-FDA0003099018260000072
Figure RE-FDA0003099018260000073
对于图层S中地理要素属性X的取值处于
Figure RE-FDA0003099018260000074
的要素，它们的Y属性的取值集合为
Figure RE-FDA0003099018260000075
方差分别为
Figure RE-FDA0003099018260000076
根据公式1.9计算每个O*的q值，其中n'表示O中的断点对图层S离散化形成的子区域内的要素个数，var表示子区域内要素的Y取值的方差，
Figure RE-FDA0003099018260000081
Figure RE-FDA0003099018260000082
Figure RE-FDA0003099018260000083
(4)Q为q的集合，对Q降序排序，根据公式2.0计算Q中每个元素的qr值，QR为qr的集合，N*表示O*中断点的数目，
qr＝q/N* (2.0)；
(5)通过三点滑动查找QR中所有元素的第一个峰值，如果qri-1＜qri＜qri+1，则令i＝i+1，如果qri＞qri+1&qri＞qri-1，则确定最优区间数的步骤已完成，qri对应的(BKS'_Len+1)为最优区间数，如果qri＜qri+1&qri＜qri-1，则qri-1对应的(BKS'_Len+1)为最优区间数。
4.根据权利要求1所述的一种新的数据空间离散化算法，其特征在于，连续属性离散化的具体步骤如下：
(1)运用经上述步骤得到的断点集合O对空间数据图层S的属性X进行划分，首先遍历X，取出Xi；
(2)对于任意的Xi，遍历O，根据公式2.1计算Xj所属的区间编号值VN，
Figure RE-FDA0003099018260000091
(3)把VN写入到空间数据图层S的D属性中；
(4)检查属性D的取值集合中是否存在取值个数小于2的编号值，如果存在将该编号值合并到相邻区间中。
5.根据权利要求1所述的一种新的数据空间离散化算法，其特征在于：所述空间数据为点状矢量数据或面状矢量数据。
6.根据权利要求5所述的一种新的数据空间离散化算法，其特征在于：所述空间数据包含的每个点要素或面要素，同时具有因变量和自变量的属性特征。
CN202110176659.2A 2021-02-07 2021-02-07 一种新的数据空间离散化算法 Pending CN113190794A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110176659.2A CN113190794A (zh) 2021-02-07 2021-02-07 一种新的数据空间离散化算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110176659.2A CN113190794A (zh) 2021-02-07 2021-02-07 一种新的数据空间离散化算法

Publications (1)

Publication Number Publication Date
CN113190794A true CN113190794A (zh) 2021-07-30

Family

ID=76972921

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110176659.2A Pending CN113190794A (zh) 2021-02-07 2021-02-07 一种新的数据空间离散化算法

Country Status (1)

Country Link
CN (1) CN113190794A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114491416A (zh) * 2022-02-23 2022-05-13 北京百度网讯科技有限公司 特征信息的处理方法、装置、电子设备和存储介质

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114491416A (zh) * 2022-02-23 2022-05-13 北京百度网讯科技有限公司 特征信息的处理方法、装置、电子设备和存储介质

Similar Documents

Publication Publication Date Title
CN109800863B (zh) 一种基于模糊理论和神经网络的测井相识别方法
CN113642849B (zh) 考虑空间分布特征的地质灾害危险性综合评价方法及装置
CN111651545A (zh) 一种基于多源数据融合的城市边缘区提取方法
CN111652198B (zh) 一种城市边缘区识别方法及系统
CN111985389B (zh) 一种基于流域属性距离的流域相似判别方法
CN110929939B (zh) 一种基于聚类-信息量耦合模型下的滑坡灾害易发性空间预测方法
Xiong et al. Drainage basin object-based method for regional-scale landform classification: A case study of loess area in China
CN116205310B (zh) 一种基于可解释集成学习模型的土壤含水量影响因素敏感区间判定方法
ARGIALAS et al. Quantitative description and classification of drainage patterns
CN106326923A (zh) 一种顾及位置重复和密度峰值点的签到位置数据聚类方法
CN116310853A (zh) 一种基于多源数据的中小城市边缘区提取方法
CN110363299A (zh) 面向露头岩层分层的空间案例推理方法
Rocca et al. A continuous scale-space method for the automated placement of spot heights on maps
CN113190794A (zh) 一种新的数据空间离散化算法
CN112634447B (zh) 一种露头岩层分层方法、装置、设备及存储介质
Dawson et al. A method for the automated extraction of environmental variables to help the classification of rivers in Britain
CN110097636B (zh) 一种基于可视域分析的选址规划方法
CN116201535B (zh) 油气藏目标井标志地层自动划分方法、装置及设备
CN112818551A (zh) 一种有效提高地理探测器模型精度的新算法
Hu et al. Using uncertain DM-chameleon clustering algorithm based on machine learning to predict landslide hazards
CN114937215B (zh) 一种城市功能区的识别方法及装置
Hidayat et al. Stream order selection for model generalization of the topographic map of Indonesia
CN118093673B (zh) 一种测绘数据处理方法
Danese et al. Geostatistics in historical macroseismic data analysis
CN111538801A (zh) 多源矢量空间数据中线状要素几何准确性评估方法和系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20210730