CN115048617B - 一种基于信息论的空间因子空间分层异质性强弱的评价方法 - Google Patents

Abstract

本发明属于空间信息技术领域，具体涉及一种基于信息论的空间因子空间分层异质性强弱的评价方法。本发明的目的是为了能够更为有效的定量分析相关因子对目标地理现象的解释能力强弱，解决现有q统计难以处理名义目标变量以及存在误判可能性两方面的不足，提供了一种基于互信息的空间因子空间分层异质性强弱的评价方法。本发明的优点在于无需大量模拟即可根据少量样本点快速有效的得到目标区域各类地物逐像元的概率分布。本发明可有效处理决策变量为类别数据以及连续值数据的情况；只要条件因子同决策变量之间不是相关独立的随机变量，该方法即可度量出条件因子的空间分层异质性。

Description

一种基于信息论的空间因子空间分层异质性强弱的评价方法

技术领域

本发明属于空间信息技术领域，具体涉及一种基于信息论的空间因子空间分层异质性强弱的评价方法。

背景技术

空间异质性是地理相关现象建模的一个重要概念，已被公认为地理科学的第二定律。一般来说，空间异质性是指地理现象在空间上的变异。它是分析地理空间数据(如人口、社区、生态系统和景观分析)的一个重要主题。空间异质性有着不同的形式，其中空间分层异质性是指目标变量在空间因子，如土地利用、气候分带、生态分区和地理区划等，所自然形成的不同分层或区域间的模式变化。

空间分层异质性在分析地理相关现象时是不可或缺的。空间分层异质性对于正确建模和理解目标地理现象非常重要，忽视不同区域间的目标现象形成机制的差异会导致有偏差甚至错误的建模。空间分层异质性可通过一致性程度发现可能存在的因果关系。目前空间分异性已被广泛应用于生产生活实践中。如分析国家级经济技术开发区经济增长率空间分异及影响因素，农地非农化空间非均衡及形成机理的探索，人口老龄化的格局分析，住宅价格的空间差异特征与影响机制的探究，污染物的空间分布特征分析，以及城市人居环境演变及驱动力研究等等。其应用领域十分广泛，并且同国计民生息息相关。

目前q统计是检测和度量空间分层异质性的唯一可用方法。其已广泛应用于许多实际问题的处理中。q统计将空间分层异质性解释为每个分层内的方差与整个研究区域的总方差之间的比率。比率越小，层内方差和总方差之间的差异越大，SSH越显著。然而，在实际应用中q统计量仍存在两个问题。一个是它不能用于处理名义目标变量。二是在某些情况下，q统计可能会忽略空间分层异质性。例如，假设研究区域内有两个分层。每一分层内的目标都服从高斯分布，且均值相同。这两个分层之间的区别是方差不同。在这种情况下，易得q统计值为零，这表示不存在空间分层异质性。然而，这两个分层之间存在着明显的异质性。

本发明从信息论的角度提出了新的空间分层异质性探测方法，以解决q-统计量的上述两个限制。该方法对于名义目标变量和连续值目标变量使用相同的流程但不同的指标进行度量。两个指标都首先使用某种互信息的扩展来度量由空间因子所形成的分层与目标变量之间的一致性程度，并须将其规范化到[0，1]区间，以度量目标变量的空间分层异质性，最后使用置换检验以验证空间分层异质性在统计意义是否显著。面向名义变量，本发明直接将其进行归一化。对于连续值变量，本发明对每个由空间因子所形成的分层和整个研究区域目标变量的概率密度函数的相对熵进行平均归一化，以定量描述其空间分层异质性。

发明内容

本发明的目的是为了能够更为有效的定量分析相关因子对目标地理现象的解释能力强弱，解决现有q统计难以处理名义目标变量以及存在误判可能性两方面的不足，提供一种基于信息论的空间因子空间分层异质性强弱的评价方法。

为解决上述技术问题，本发明采用的技术方案，一种基于信息论的空间因子空间分层异质性强弱的评价方法，包含以下步骤：

步骤1.根据条件因子S将目标地理区域划分为m个互不相交的子区域，分别为s₁,s₂,...,s_m；其中每个子区域内所有对象在条件因子S上的取值完全相同；记为子区域s_i内对象个数所占比例，其中|s_i|为s_i子区域内对象个数，n为目标地理区域内所有对象的总数；

步骤2.获取目标地理区域内所有n个对象在目标变量d上的测量值，此时目标地理区域内每个对象上有一个其所属子区域的描述以及其目标变量一个测量值；

步骤3.根据目标变量d的测量值以及条件因子S构建的m个互不相交的子区域，计算S与d之间的扩展互信息，具体计算方法为：

(1)若目标变量d为名义变量，则使用公式

计算S与d的扩展互信息；其中，I(d,S)为目标变量d与条件因子S的互信息，H(d)为目标变量d的信息熵；具体的若n个对象在d上的可能取值有w种，这w种取值记为d₁，d₂，...，d_w，且其所占比例分别为p_d1，p_d2，...，p_dw，则同时，记/>为s_i区域内目标变量取值为d_i的对象的个数与该区域所有对象的个数比值，则

(2)若目标变量d为连续值变量，则首先将目标变量d在所有对象上的最大值和最小值分别记为Maxd和Mind，并将区间[Mind,Maxd]等分为k个区间，记为D₁,...,D_k；对于目标地理区域胡每个对象x，若其在d上的取值d(x)满足i<(d(x)-Mind)×k/(Maxd-Mind)≤i+1，则对象x在d上的取值位于区间D_i；当所有对象处理结束后，分别计算每个区间中的对象个数与所有对象个数的比值，并记为并且在每个子区域s_i中计算每个区间中的对象个数与该子区域内所有对象个数的比值，并记为/>随后，对于每个分区s_i计算

最后使用

计算S与d之间的扩展互信息；

步骤4.将原始数据随机置换N次，在每次置换中，首先随机置换所有对象在目标变量d上的取值，然后按照步骤3中的方法计算S与d之间的扩展互信息；所有置换完成后，使用公式(N′+1)/(N+1)计算S与d之间的扩展互信息的统计显著性，其中，N'为所有置换中出现置换后扩展互信息大于等于原始数据扩展互信息的次数；

步骤5.S与d之间的扩展互信息及(N′+1)/(N+1)即可用以度量S因子对于目标变量的空间分层异质性程度。

与现有技术相比，本发明具有以下特点：

(1)可有效处理决策变量为类别数据以及连续值数据的情况；

(2)只要条件因子同决策变量之间不是相关独立的随机变量，该方法即可度量出条件因子的空间分层异质性。

附图说明

图1为本发明的实现流程图。

具体实施方式

如图1所示，本实施例中使用了一个关于房地产价格的相关数据集，其下载地址为：https://geodacenter.github.io/data-and-lab/baltim/。我们选择了该数据集中的三个因子对本专利实施方式进行说明，分别为该房产是否包含一个独立的单元、是否在本地居住、以及房屋出售价格。为方便起见，后面使用DWELL,CITCOU以及PRICE来表示这三个因子。我们按照以下步骤评价DWELL用以对CITCOU进行描述以及DWELL用以对PRICE进行描述时的空间分层异质性的强弱：

步骤1.将DWELL作为条件因子S对研究区域进行划分。由于该实例中所有数据对象在S上的取值为0或1,可以将其划分为两个不同的子区域，子区域内S的取值完全相同，并记这两个子区域记为S＝{s₁,s₂}。

步骤2.在该数据集中，有两个目标变量，分别为名义变量CITCOU和连续值变量PRICE。在该步骤中，我们可以计算得到每个子区域中对象的个数以及其所占比例，具体的s₁中对象的数目为98，s₂中的对象数目为113、s₁中对象个数所占比例为s₂中对象个数所占比例为/>

步骤3.(1)对于目标变量CITCOU，由于其为名义变量，因此使用公式

计算DWELL属性的空间分层异质性。具体的，所有对象在CITCOU上的可能取值有2种，分别为d₁＝0，d₂＝1。两种类别所占比例分别为和/>由此可得

同时，s₁区域内目标变量取值为d₁的对象所占比例为62/211＝0.2939，s₁区域内目标变量取值为d₂的对象所占比例为36/211＝0.1706，s₂区域内目标变量取值为d₁的对象所占比例为21/211＝0.0995，s₂区域内目标变量取值为d₂的对象所占比例为92/211＝0.4360。因此

最终计算I_N(d,S)＝0.1077/0.6702＝0.1606

(2)目标变量PRICE为连续值变量，其最大值和最小值分别为165和3.5。本实施实例中，首先将区间[3.5,165]划分为k＝6个区间，分别为D₁＝[3.5,30.42)、D₂＝[30.42,57.34)、D₃＝[57.34,84.26)、D₄＝[84.26,111.18)、D₅＝[111.18,138.1)、D₆＝[138.1,165)。然后，将每个对象PRICE上的值映射到这6个区间上。例如，若对象x的PRICE的取d(x)为92，由于(92-3.5)×6/(165-3.5)在4和5之间，其取值改为D₄。之后，对于每个子区域s_i使用下述公式

计算可得R(s₁,d)＝0.2076,R(s₂,d)＝0.1714。最后使用

计算S与d之间的扩展互信息I_c(d,S)＝0.1428。

步骤4.进行置换检验。将所有对象的目标变量CITCOU和PRICE分别进行随机置换，基于置换后的数据，按照步骤3中的方法重新计算扩展互信息，我们将原始数据随机置换N＝1000次，所有置换完成后，记录打乱之后的扩展互信息大于打乱之前的扩展互信息的次数N'＝0，因此(N′+1)/(N+1)＝1/1001<0.01。

通过上述步骤可得，DWELL用以对CITCOU进行描述以及DWELL用以对PRICE进行描述时的空间分层异质性程度分别为0.1606和0.1428，且二者在置信度为0.01时，均统计显著。

Claims (1)

1.一种基于信息论的空间因子空间分层异质性强弱的评价方法，其特征在于:包括以下步骤：

步骤1.根据条件因子S将目标地理区域划分为m个互不相交的子区域，分别为s₁,s₂,...,s_m；其中每个子区域内所有对象在条件因子S上的取值完全相同；记为子区域s_i内对象个数所占比例，其中|s_i|为s_i子区域内对象个数，n为目标地理区域内所有对象的总数；

步骤2.获取目标地理区域内所有n个对象在目标变量d上的测量值，此时目标地理区域内每个对象上有一个其所属子区域的描述以及其目标变量一个测量值；

步骤3.根据目标变量d的测量值以及条件因子S构建的m个互不相交的子区域，计算S与d之间的扩展互信息，具体计算方法为：

(1)若目标变量d为名义变量，则使用公式

计算S与d的扩展互信息；其中，I(d,S)为目标变量d与条件因子S的互信息，H(d)为目标变量d的信息熵；具体的若n个对象在d上的可能取值有w种，这w种取值记为d₁，d₂，...，d_w，且其所占比例分别为则/>同时，记/>为s_i区域内目标变量取值为d_i的对象的个数与该区域所有对象的个数比值，则

(2)若目标变量d为连续值变量，则首先将目标变量d在所有对象上的最大值和最小值分别记为Maxd和Mind，并将区间[Mind,Maxd]等分为k个区间，记为D₁,...,D_k；对于目标地理区域的每个对象x，若其在d上的取值d(x)满足i<(d(x)-Mind)×k/(Maxd-Mind)≤i+1，则对象x在d上的取值位于区间D_i；当所有对象处理结束后，分别计算每个区间中的对象个数与所有对象个数的比值，并记为并且在每个子区域s_i中计算每个区间中的对象个数与该子区域内所有对象个数的比值，并记为/>随后，对于每个分区s_i计算

最后使用

计算S与d之间的扩展互信息；

步骤4.将原始数据随机置换N次，在每次置换中，首先随机置换所有对象在目标变量d上的取值，然后按照步骤3中的方法计算S与d之间的扩展互信息；所有置换完成后，使用公式(N′+1)/(N+1)计算S与d之间的扩展互信息的统计显著性，其中，N'为所有置换中出现置换后扩展互信息大于等于原始数据扩展互信息的次数；

步骤5.S与d之间的扩展互信息及(N′+1)/(N+1)即可用以度量S因子对于目标变量的空间分层异质性程度。