CN106951828A - 一种基于卫星影像和网络的城市区域功能属性的识别方法 - Google Patents

一种基于卫星影像和网络的城市区域功能属性的识别方法 Download PDF

Info

Publication number
CN106951828A
CN106951828A CN201710097718.0A CN201710097718A CN106951828A CN 106951828 A CN106951828 A CN 106951828A CN 201710097718 A CN201710097718 A CN 201710097718A CN 106951828 A CN106951828 A CN 106951828A
Authority
CN
China
Prior art keywords
urban area
user
network
access
group
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
Application number
CN201710097718.0A
Other languages
English (en)
Other versions
CN106951828B (zh
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.)
Tsinghua University
Original Assignee
Tsinghua 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 Tsinghua University filed Critical Tsinghua University
Priority to CN201710097718.0A priority Critical patent/CN106951828B/zh
Publication of CN106951828A publication Critical patent/CN106951828A/zh
Application granted granted Critical
Publication of CN106951828B publication Critical patent/CN106951828B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • G06V20/176Urban or other man-made structures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/25Fusion techniques
    • G06F18/253Fusion techniques of extracted features
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Biology (AREA)
  • Biophysics (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Computational Linguistics (AREA)
  • Biomedical Technology (AREA)
  • Health & Medical Sciences (AREA)
  • Multimedia (AREA)
  • Navigation (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本发明涉及一种基于卫星影像和网络的城市区域功能属性的识别方法,属于多模态传感器数据处理技术领域。首先,收集在城市区域内分布的路由器上活动的用户的连接行为信息以及整个城市区域对应的卫星影像。其次,基于收集到的数据进行城市区域的特征定义与提取,即对网络数据进行用户行为特征的提取,对卫星影像特征通过卷积神经网络进行视觉特征提取。最后,采用AdaBoost算法整合一系列的决策树模型进行城市区域功能属性识别,从而增加在训练的迭代过程中发现关键特征的稳定性。本发明的优点是融合了两种模态的传感器数据,提高了城市区域功能属性识别的灵活性和准确率。本方法采用的数据易获取,算法的计算复杂度低,可扩展性强,适合于实际工业应用。

Description

一种基于卫星影像和网络的城市区域功能属性的识别方法
技术领域
本发明涉及一种基于卫星影像和网络的城市区域功能属性的识别方法,属于多模态传感器数据处理技术领域。
背景技术
具备不同功能属性的城市区域是随着城市的自然演化自然形成的,并为满足居民日常生活所需提供不同的功能服务。对城市区域功能的获取可以方便很多有价值的服务,例如城市规划、商业选址、出行建议等。目前一类自动进行城市区域功能判别的方法大多依赖于对用户在区域间移动的数据分析。虽然用户在区域间转移的行为模式在一定程度上反映了区域的功能特征,但是用户在区域内部的行为模式在当前的研究中并没有充分利用,以更好地理解用户的行为如驻留时长特征是如何受到区域功能影响的。
通过在城市范围内部署的大量Wi-Fi无线接入点,可以获得大量用户的Wi-Fi连接记录。相比于GPS轨迹具有在建筑内部信号缺失的问题,以及手机基站定位精度不高的问题,这种数据形式具有体量大、精度高、富含用户行为特征的优势。除了Wi-Fi访问提供的用户行为特征,城市基础设施的地理位置分布,如公园、建筑、公路等,可以提供城市功能的静态视觉特征反映。采用高清的卫星影像,可以基于公共基础设施的俯视图特征对城市区域进行分析。通过对Wi-Fi用户行为和卫星影像这两种城市传感器数据进行融合,可以提供更好的功能属性识别方法。
综上所述,融合卫星影像和Wi-Fi连接进行城市区域功能属性的识别方法是利用多种传感器数据的识别能力对城市区域功能属性进行有效识别的关键技术。这种城市区域功能属性的识别方法的作用是给出城市范围内的大量用户Wi-Fi的连接记录以及城市卫星影像数据,通过对行为特征和视觉特征进行提取,并设计有效分类器进行实现。随着城市化进程的快速发展,以及城市日常运营过程中产生的大量数据,这种城市区域功能属性的识别方法不可避免地需要满足如下要求:1)对多种数据类型应用的灵活性;2)对不同数量和质量的城市数据的适应性;3)对多种区域功能类型识别的扩展性。在城市区域功能属性识别的过程中,应该尽量保证上述三点,这样才能保证区域功能属性的识别在视城市大数据上的灵活应用,否则在区域功能属性的识别的过程中,很难利用好多种数据来源的识别能力,从而影响最终功能识别的效果。因此,城市区域功能属性的识别算法需要尽可能的保证以上三点要求。
发明内容
本发明提出了一种基于卫星影像和网络的城市区域功能属性的识别方法,其核心在于通过分析易收集的大规模的城市Wi-Fi上用户连接信息以及城市区域的地理卫星影像,来自动地感知城市中不同区域的功能特征,从而让城市居民以及城市的管理者更好地理解城市区域分布的现状以及发展的进程,为包括旅游路线选择、商业选址以及城市建设规划等方面的实际应用做出建议和贡献。
本发明提出的基于卫星影像和网络的城市区域功能属性的识别方法,包括以下步骤:
(1)在城市区域内布置供网络用户免费使用的路由器,利用路由器收集用户访问网络的时间TA和离开网络的时间TL,或者城市区域内网络用户利用手机自动数据收集程序发送访问网络的时间TA和离开网络的时间TL,每次用户访问的持续时间记为D,则:
D=TL-TA
记用户第i次访问网络的行为Bi,Bi用一个四元组(TA,i,TL,i,Di,Uidi)表示,其中Uidi代表该用户的代码ID,每个用户对应一个单独的代码ID,记所有用户访问的总次数为Q,记所有用户访问行为构成的集合为Gall
(2)利用电子地图软件下载城市区域的相应卫星影像,卫星影像的大小为256*256像素,每个像素用RGB颜色值表示;
(3)根据上述步骤(1)的用户访问网络时间TA,离开网络时间TL以及持续时间D,计算城市区域的网络特征向量VN,该特征向量VN包括36维,计算过程如下:
(3-1)根据用户访问网络时间TA处于一天24小时中的某一小时,将所有用户的访问行为分为24组,记为G[1,2,…,24],其中G[a]组中的用户访问网络时间TA满足a-1≤TA<a,各组的用户访问行为的数量作为所述的特征向量VN的1-24维VN[1,2,…,24];
(3-2)根据用户访问网络时间TA处于周中或周末,将所有用户的访问行为分为2组,记为G[25,26],其中G[25]组中的用户访问发生在周中,G[26]组中的用户访问发生在周末,各组的用户访问行为数量作为特征向量VN的25-26维VN[25,26];
(3-3)分别计算上述步骤(3-2)中G[25]组和G[26]组中用户访问网络的平均持续时间AvgD[25,26]作为特征向量VN的27-28维VN[27,28]:
(3-4)设置一个用户访问网络的持续时间阈值Dth,计算用户访问网络持续时间大于该阈值的用户访问行为的数量占所有用户访问行为总数量的比率Rlong作为特征向量VN的29维VN[29]:
其中,I(·)为指示函数,当条件·为真时,I(·)取值为1,当条件·为假时,I(·)取值为0;
(3-5)将G[25]组,G[26]组以及Gall中的用户数量Nu作为特征向量VN的30-32维VN[30,31,32];
(3-6)分别计算G[25]组用户访问行为、G[26]组用户访问行为和Gall用户访问行为的平均重连时间Tre作为特征向量VN的33-35维VN[33,34,35],平均重连时间Tre的计算步骤如下(以G[25]为例):
(3-6-1)根据用户的代码IDUid,将G[25]组的用户访问行为分成Nu个部分,Nu为与该组用户访问行为相对应的用户数量,即每个部分中的所有用户访问行为对应同一个Uid,对于Nu个部分中的第c部分用户访问行为,将上述步骤(1)的访问网络时间TA从小到大排序,得到序列B[1,2,…,p],其中p为第c部分用户访问行为的数量,则第c部分用户访问行为的平均重连时间Tre,c为:
其中Maximum为设定的平均重连时间Tre极大值,取值为168小时;
(3-6-2)利用下式计算G[25]组用户访问行为的平均重连时间Tre为:
(3-6-3)重复步骤(3-6-1)和步骤(3-6-2),分别计算G[26]组用户访问行为和Gall用户访问行为的平均重连时间;
(3-7)根据上述步骤(1)的访问网络时间TA,将所有用户访问行为按照日期顺序记为G′[1,2,…,ND],其中每组用户访问行为发生在同一天,计算G′[1,2,…,ND]中第e天的用户群变化率CRCe
其中N′u,e为与G′[e]组访问行为相对应的用户数量,G′[e]为第e天的的用户访问行为,com(e-1,e)为与G′[e-1]组访问行为和G′[e]组访问行为相对应的用户中共同重复的用户数量;
利用下式计算平均用户群变化率CRC作为特征向量VN的36维VN[36]:
(4)将上述步骤(2)的卫星影像中输入到卷积神经网络中,将卷积神经网络的第二个全连接层的输出作为影像特征向量VI[1,2,…,4096];
(5)将上述步骤(3)得到的城市区域的网络特征向量VN[1,2,…,36]和上述步骤(4)得到的城市区域的卫星影像特征向量Vi[1,2,…,4096]连接,合并成向量V[1,2,…,4132],作为城市区域的特征向量;
(6)建立功能属性分类模型,过程如下:
(6-1)将具有已知功能属性的N个城市区域作为训练样本,记N个城市区域中任意一个城市区域为Rf,1≤f≤N,具有已知的功能属性tf,利用上述步骤(5)得到特征向量Vf,遍历N个城市区域中的所有城市区域,得到具有已知功能属性的N个城市区域的特征向量;
(6-2)利用AdaBoost算法,训练基于决策树的功能属性分类模型步骤如下:
(6-2-1)定义一个N维样本权重向量WN,权重向量wN中的每个元素对应一个具有已知功能属性的城市区域训练样本,将每个元素初始化为
(6-2-2)将误差函数εm最小化,得到决策树函数ym(V):
其中,I(·)为指示函数,当条件·为真时,I(·)取值为1,当条件·为假时,I(·)取值为0,m为迭代次数,V为上述步骤(5)中的城市区域的特征向量,Vf为上述步骤(6-1)中的N个城市区域中任意一个城市区域Rf的特征向量;
利用下式,计算决策树函数ym(V)的权重αm
利用下式调整权重向量wN
其中,Zm是归一化因子,保证权重向量wN的和为1;
(6-2-3)重复上述步骤(6-2-2)M次,得到M个决策树模型;
(6-2-4)将步骤(6-2-3)得到的M个决策树模型加权,得到如下功能属性分类模型:
(7)对未知功能属性的城市区域R,重复步骤(1)-(5),得到特征向量V,将V输入步骤(6-2-4)的功能属性分类模型YM(V),输出即为城市区域R的功能属性。
本发明提出的基于卫星影像和网络的城市区域功能属性的识别方法,其优点是:
1、本发明基于卫星影像和网络的城市区域功能属性的识别方法中,充分使用城市区域的卫星影像和群体Wi-Fi连接数据,对两种传感器的隐含特征进行有机融合,从而增强了城市区域功能属性识别的准确率。由于本方法所需要的传感器数据易于获取,方法直接应用于对城市功能属性的识别,具有很强的灵活性和适应性。
2、本发明的用于城市区域功能属性的识别方法中,从卫星影像和Wi-Fi连接记录中所提取的特征能够更好反映城市区域的功能属性特征。从Wi-Fi连接记录中提取的群体行为特征受区域本身提供的功能所影响,群体用户在区域间转移及在单个区域的驻留时空特征等能更好地反映城市的功能属性。对区域地理卫星影像特征的提取采用最新的卷积神经网络的方法,所采用的视觉特征对不同的城市区域具有很强的区分度。本发明算法将两者结合起来,从而极大提高了对城市区域功能属性的识别效果。
3、本发明的用于城市区域功能属性的识别方法中,提出了采用AdaBoost算法来整合一系列的决策树模型组成一个优化的决策树森林,从而增加在训练的迭代过程中发现关键特征的稳定性。该方法将多个弱分类器有机结合,从而提高了最终城市区域功能属性的识别方法的效果。
附图说明
图1是本发明方法的流程框图。
图2是本发明方法中使用的卷积神经网络的结构图。
图3是本发明方法中模型训练算法流程图。
具体实施方式
本发明提出的基于卫星影像和网络的城市区域功能属性的识别方法,其流程框图如图1所示,该方法包括以下步骤:
(1)在城市区域内布置供网络用户免费使用的路由器,利用路由器收集用户访问网络的时间TA和离开网络的时间TL,或者城市区域内网络用户利用手机自动数据收集程序发送访问网络的时间TA和离开网络的时间TL,每次用户访问的持续时间记为D,则:
D=TL-TA
记用户第i次访问网络的行为Bi,Bi用一个四元组(TA,i,TL,i,Di,Uidi)表示,其中Uidi代表该用户的代码ID,每个用户对应一个单独的代码ID,记所有用户访问的总次数为Q,记所有用户访问行为构成的集合为Gall
(2)利用电子地图软件(例如腾讯地图和谷歌地图)下载城市区域的相应卫星影像,卫星影像的大小为256*256像素,每个像素用RGB颜色值表示;
(3)根据上述步骤(1)的用户访问网络时间TA,离开网络时间TL以及持续时间D,计算城市区域的网络特征向量VN,该特征向量VN包括36维,计算过程如下:
(3-1)根据用户访问网络时间TA处于一天24小时中的某一小时,将所有用户的访问行为分为24组,记为G[1,2,…,24],其中G[a]组中的用户访问网络时间TA满足a-1≤TA<a,各组的用户访问行为的数量作为所述的特征向量VN的1-24维VN[1,2,…,24];
(3-2)根据用户访问网络时间TA处于周中或周末,将所有用户的访问行为分为2组,记为G[25,26],其中G[25]组中的用户访问发生在周中,G[26]组中的用户访问发生在周末,各组的用户访问行为数量作为特征向量VN的25-26维VN[25,26];
(3-3)分别计算上述步骤(3-2)中G[25]组和G[26]组中用户访问网络的平均持续时间AvgD[25,26]作为特征向量VN的27-28维VN[27,28]:
(3-4)设置一个用户访问网络的持续时间阈值Dth,持续时间阈值Dth的取值可以为30分钟,计算用户访问网络持续时间大于该阈值的用户访问行为的数量占所有用户访问行为总数量的比率Rlong作为特征向量VN的29维VN[29]:
其中,I(·)为指示函数,当条件·为真时,I(·)取值为1,当条件·为假时,I(·)取值为0;
(3-5)将G[25]组,G[26]组以及Gall中的用户数量Nu作为特征向量VN的30-32维VN[30,31,32];
(3-6)分别计算G[25]组用户访问行为、G[26]组用户访问行为和Gall用户访问行为的平均重连时间Tre作为特征向量VN的33-35维VN[33,34,35],平均重连时间Tre的计算步骤如下(以G[25]为例):
(3-6-1)根据用户的代码IDUid,将G[25]组的用户访问行为分成Nu个部分,Nu为与该组用户访问行为相对应的用户数量,即每个部分中的所有用户访问行为对应同一个Uid,对于Nu个部分中的第c部分用户访问行为,将上述步骤(1)的访问网络时间TA从小到大排序,得到序列B[1,2,…,p],其中p为第c部分用户访问行为的数量,则第c部分用户访问行为的平均重连时间Tre,c为:
其中Maximum为设定的平均重连时间Tre极大值,取值为168小时;
(3-6-2)利用下式计算G[25]组用户访问行为的平均重连时间Tre为:
(3-6-3)重复步骤(3-6-1)和步骤(3-6-2),分别计算G[26]组用户访问行为和Gall用户访问行为的平均重连时间;
(3-7)根据上述步骤(1)的访问网络时间TA,将所有用户访问行为按照日期顺序记为G′[1,2,…,ND],其中每组用户访问行为发生在同一天,计算G′[1,2,…,ND]中第e天的用户群变化率CRCe
其中N′u,e为与G′[e]组访问行为相对应的用户数量,G′[e]为第e天的的用户访问行为,com(e-1,e)为与G′[e-1]组访问行为和G′[e]组访问行为相对应的用户中共同重复的用户数量;
利用下式计算平均用户群变化率CRC作为特征向量VN的36维VN[36]:
(4)将上述步骤(2)的卫星影像中输入到卷积神经网络中,网络结构如图2所示,将卷积神经网络的第二个全连接层(fc7)的输出(4096维)作为影像特征向量VI[1,2,…,4096];
选择这一层的原因是该层距离最后输出层较近,隐含特征已经在之前的网络层被充分挖掘融合,同时特征维数相比最后一个全连接层(fc8)的1000层要高,可以保留更多特征;
(5)将上述步骤(3)得到的城市区域的网络特征向量VN[1,2,…,36]和上述步骤(4)得到的城市区域的卫星影像特征向量VI[1,2,…,4096]连接,合并成向量V[1,2,…,4132],作为城市区域的特征向量;
(6)建立功能属性分类模型,过程如下:
(6-1)将具有已知功能属性的N个城市区域作为训练样本,记N个城市区域中任意一个城市区域为Rf,1≤f≤N,具有已知的功能属性tf,利用上述步骤(5)得到特征向量Vf,遍历N个城市区域中的所有城市区域,得到具有已知功能属性的N个城市区域的特征向量;
(6-2)利用AdaBoost算法,训练基于决策树的功能属性分类模型,具体流程如图3所示,步骤如下:
(6-2-1)定义一个N维样本权重向量wN,权重向量wN中的每个元素对应一个具有已知功能属性的城市区域训练样本,将每个元素初始化为
(6-2-2)将误差函数εm最小化,得到决策树函数ym(V):
其中,I(·)为指示函数,当条件·为真时,I(·)取值为1,当条件·为假时,I(·)取值为0,m为迭代次数,V为上述步骤(5)中的城市区域的特征向量,Vf为上述步骤(6-1)中的N个城市区域中任意一个城市区域Rf的特征向量;
利用下式,计算决策树函数ym(V)的权重αm
利用下式调整权重向量wN
其中,Zm是归一化因子,保证权重向量wN的和为1;
在此步骤中,4132维特征向量V中对城市区域的功能属性分类贡献较大的维度将被自动赋予更高的权重,从而实现网络特征部分和影像特征部分的有效结合;
(6-2-3)重复上述步骤(6-2-2)M次,得到M个决策树模型;
(6-2-4)将步骤(6-2-3)得到的M个决策树模型加权,得到如下功能属性分类模型:
(7)对未知功能属性的城市区域R,重复步骤(1)-(5),得到特征向量V,将V输入步骤(6-2-4)的功能属性分类模型YM(V),输出即为城市区域R的功能属性。

Claims (1)

1.一种基于卫星影像和网络的城市区域功能属性的识别方法,其特征在于该方法包括以下步骤:
(1)在城市区域内布置供网络用户免费使用的路由器,利用路由器收集用户访问网络的时间TA和离开网络的时间TL,或者城市区域内网络用户利用手机自动数据收集程序发送访问网络的时间TA和离开网络的时间TL,每次用户访问的持续时间记为D,则:
D=TL-TA
记用户第i次访问网络的行为Bi,Bi用一个四元组(TA,i,TL,i,Di,Uidi)表示,其中Uidi代表该用户的代码ID,每个用户对应一个单独的代码ID,记所有用户访问的总次数为Q,记所有用户访问行为构成的集合为Gall
(2)利用电子地图软件下载城市区域的相应卫星影像,卫星影像的大小为256*256像素,每个像素用RGB颜色值表示;
(3)根据上述步骤(1)的用户访问网络时间TA,离开网络时间TL以及持续时间D,计算城市区域的网络特征向量VN,该特征向量VN包括36维,计算过程如下:
(3-1)根据用户访问网络时间TA处于一天24小时中的某一小时,将所有用户的访问行为分为24组,记为G[1,2,…,24],其中G[a]组中的用户访问网络时间TA满足a-1≤TA<a,各组的用户访问行为的数量作为所述的特征向量VN的1-24维VN[1,2,…,24];
(3-2)根据用户访问网络时间TA处于周中或周末,将所有用户的访问行为分为2组,记为G[25,26],其中G[25]组中的用户访问发生在周中,G[26]组中的用户访问发生在周末,各组的用户访问行为数量作为特征向量VN的25-26维VN[25,26];
(3-3)分别计算上述步骤(3-2)中G[25]组和G[26]组中用户访问网络的平均持续时间AvgD[25,26]作为特征向量VN的27-28维VN[27,28]:
A v g D &lsqb; b &rsqb; = &Sigma; u &Element; G &lsqb; b &rsqb; D u V N &lsqb; b &rsqb; , b = 25 , 26
(3-4)设置一个用户访问网络的持续时间阈值Dth,计算用户访问网络持续时间大于该阈值的用户访问行为的数量占所有用户访问行为总数量的比率Rlong作为特征向量VN的29维VN[29]:
R l o n g = &Sigma; j = 1 Q I ( D j > D t h ) Q
其中,I(·)为指示函数,当条件·为真时,I(·)取值为1,当条件·为假时,I(·)取值为0;
(3-5)将G[25]组,G[26]组以及Gall中的用户数量Nu作为特征向量VN的30-32维VN[30,31,32];
(3-6)分别计算G[25]组用户访问行为、G[26]组用户访问行为和Gall用户访问行为的平均重连时间Tre作为特征向量VN的33-35维VN[33,34,35],平均重连时间Tre的计算步骤如下(以G[25]为例):
(3-6-1)根据用户的代码IDUid,将G[25]组的用户访问行为分成Nu个部分,Nu为与该组用户访问行为相对应的用户数量,即每个部分中的所有用户访问行为对应同一个Uid,对于Nu个部分中的第c部分用户访问行为,将上述步骤(1)的访问网络时间TA从小到大排序,得到序列B[1,2,…,p],其中p为第c部分用户访问行为的数量,则第c部分用户访问行为的平均重连时间Tre,c为:
T r e , c = &Sigma; k = 1 p - 1 ( T A , k + 1 - T L , k ) p - 1 , p &GreaterEqual; 2 M a x i m u m , p = 1
其中Maximum为设定的平均重连时间Tre极大值,取值为168小时;
(3-6-2)利用下式计算G[25]组用户访问行为的平均重连时间Tre为:
T r e = &Sigma; c = 1 N u T r e , c N u
(3-6-3)重复步骤(3-6-1)和步骤(3-6-2),分别计算G[26]组用户访问行为和Gall用户访问行为的平均重连时间;
(3-7)根据上述步骤(1)的访问网络时间TA,将所有用户访问行为按照日期顺序记为G′[1,2,…,ND],其中每组用户访问行为发生在同一天,计算G′[1,2,…,ND]中第e天的用户群变化率CRCe
CRC e = 1 - c o m ( e - 1 , e ) N u , e - 1 &prime; + N u , e &prime; - c o m ( e - 1 , e ) , 2 &le; e &le; N D
其中N′u,e为与G′[e]组访问行为相对应的用户数量,G′[e]为第e天的的用户访问行为,com(e-1,e)为与G′[e-1]组访问行为和G′[e]组访问行为相对应的用户中共同重复的用户数量;
利用下式计算平均用户群变化率CRC作为特征向量VN的36维VN[36]:
C R C = &Sigma; e = 2 N D CRC e N D - 1
(4)将上述步骤(2)的卫星影像中输入到卷积神经网络中,将卷积神经网络的第二个全连接层的输出作为影像特征向量VI[1,2,…,4096];
(5)将上述步骤(3)得到的城市区域的网络特征向量VN[1,2,…,36]和上述步骤(4)得到的城市区域的卫星影像特征向量VI[1,2,…,4096]连接,合并成向量V[1,2,…,4132],作为城市区域的特征向量;
(6)建立功能属性分类模型,过程如下:
(6-1)将具有已知功能属性的N个城市区域作为训练样本,记N个城市区域中任意一个城市区域为Rf,1≤f≤N,具有已知的功能属性tf,利用上述步骤(5)得到特征向量Vf,遍历N个城市区域中的所有城市区域,得到具有已知功能属性的N个城市区域的特征向量;
(6-2)利用AdaBoost算法,训练基于决策树的功能属性分类模型步骤如下:
(6-2-1)定义一个N维样本权重向量wN,权重向量wN中的每个元素对应一个具有已知功能属性的城市区域训练样本,将每个元素初始化为
(6-2-2)将误差函数εm最小化,得到决策树函数ym(V):
&epsiv; m = &Sigma; f = 1 N w N , f ( m ) I ( y m ( V f ) &NotEqual; t f )
其中,I(·)为指示函数,当条件·为真时,I(·)取值为1,当条件·为假时,I(·)取值为0,m为迭代次数,V为上述步骤(5)中的城市区域的特征向量,Vf为上述步骤(6-1)中的N个城市区域中任意一个城市区域Rf的特征向量;
利用下式,计算决策树函数ym(V)的权重αm
&alpha; m = l n { 1 - &epsiv; m &epsiv; m }
利用下式调整权重向量wN
w N , f ( m + 1 ) = w N , f ( m ) Z m * exp ( - &alpha; m ) , y m ( V f ) = t f exp ( &alpha; m ) , y m ( V f ) &NotEqual; t f , 1 &le; f &le; N ,
其中,Zm是归一化因子,保证权重向量wN的和为1;
(6-2-3)重复上述步骤(6-2-2)M次,得到M个决策树模型;
(6-2-4)将步骤(6-2-3)得到的M个决策树模型加权,得到如下功能属性分类模型:
Y M ( V ) = &Sigma; m = 1 M &alpha; m y m ( V )
(7)对未知功能属性的城市区域R,重复步骤(1)-(5),得到特征向量V,将V输入步骤(6-2-4)的功能属性分类模型YM(V),输出即为城市区域R的功能属性。
CN201710097718.0A 2017-02-22 2017-02-22 一种基于卫星影像和网络的城市区域功能属性的识别方法 Expired - Fee Related CN106951828B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710097718.0A CN106951828B (zh) 2017-02-22 2017-02-22 一种基于卫星影像和网络的城市区域功能属性的识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710097718.0A CN106951828B (zh) 2017-02-22 2017-02-22 一种基于卫星影像和网络的城市区域功能属性的识别方法

Publications (2)

Publication Number Publication Date
CN106951828A true CN106951828A (zh) 2017-07-14
CN106951828B CN106951828B (zh) 2020-02-21

Family

ID=59467526

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710097718.0A Expired - Fee Related CN106951828B (zh) 2017-02-22 2017-02-22 一种基于卫星影像和网络的城市区域功能属性的识别方法

Country Status (1)

Country Link
CN (1) CN106951828B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107392463A (zh) * 2017-07-18 2017-11-24 中国地质大学(武汉) 一种城市功能区识别方法、模块、装置及其存储装置
CN108287958A (zh) * 2018-01-15 2018-07-17 北京航空航天大学 人工智能程序员书写数字飞行器源代码有限选择决策方法
CN110442662A (zh) * 2019-07-08 2019-11-12 清华大学 一种确定用户属性信息的方法以及信息推送方法
CN110659427A (zh) * 2019-09-06 2020-01-07 北京百度网讯科技有限公司 基于多源数据的城市功能划分方法、装置和电子设备
CN111126120A (zh) * 2018-11-01 2020-05-08 百度在线网络技术(北京)有限公司 城市区域分类方法、装置、设备和介质
CN112364776A (zh) * 2020-11-12 2021-02-12 中国地质大学(武汉) 一种城市功能区识别模型建立、分类方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102572799A (zh) * 2011-12-29 2012-07-11 华为终端有限公司 获取接入Wi-Fi热点的客户终端的网络连接信息的方法和终端
CN103035036A (zh) * 2012-12-05 2013-04-10 福建省计量科学研究院 基于无线网络和卫星定位实现远程监控出租车计价器方法
US20150025914A1 (en) * 2013-07-16 2015-01-22 Esurance Insurance Services, Inc. Property inspection using aerial imagery
CN104968027A (zh) * 2015-05-29 2015-10-07 努比亚技术有限公司 Wi-Fi控制方法及装置
CN106028283A (zh) * 2016-05-26 2016-10-12 上海网数信息科技有限公司 通过无线网络数据分析用户线下行为的方法及系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102572799A (zh) * 2011-12-29 2012-07-11 华为终端有限公司 获取接入Wi-Fi热点的客户终端的网络连接信息的方法和终端
CN103035036A (zh) * 2012-12-05 2013-04-10 福建省计量科学研究院 基于无线网络和卫星定位实现远程监控出租车计价器方法
US20150025914A1 (en) * 2013-07-16 2015-01-22 Esurance Insurance Services, Inc. Property inspection using aerial imagery
CN104968027A (zh) * 2015-05-29 2015-10-07 努比亚技术有限公司 Wi-Fi控制方法及装置
CN106028283A (zh) * 2016-05-26 2016-10-12 上海网数信息科技有限公司 通过无线网络数据分析用户线下行为的方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ZHANG MINGGUO ETAL.: ""UAV Image Identification in Urban Region Satellite Image Using Global Feature and Local Feature"", 《IEEE》 *
王鹏 等: ""面向城市规划编制的大数据类型及应用方式研究"", 《规划师论坛》 *
钱家航 等: ""基于QuickBird影像城市道路特征语义信息提取"", 《遥感技术与应用》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107392463A (zh) * 2017-07-18 2017-11-24 中国地质大学(武汉) 一种城市功能区识别方法、模块、装置及其存储装置
CN107392463B (zh) * 2017-07-18 2021-01-01 中国地质大学(武汉) 一种城市功能区识别方法、模块、装置及其存储装置
CN108287958A (zh) * 2018-01-15 2018-07-17 北京航空航天大学 人工智能程序员书写数字飞行器源代码有限选择决策方法
CN108287958B (zh) * 2018-01-15 2021-05-11 北京航空航天大学 人工智能程序员书写数字飞行器源代码有限选择决策方法
CN111126120A (zh) * 2018-11-01 2020-05-08 百度在线网络技术(北京)有限公司 城市区域分类方法、装置、设备和介质
CN111126120B (zh) * 2018-11-01 2024-02-23 百度在线网络技术(北京)有限公司 城市区域分类方法、装置、设备和介质
CN110442662A (zh) * 2019-07-08 2019-11-12 清华大学 一种确定用户属性信息的方法以及信息推送方法
CN110442662B (zh) * 2019-07-08 2022-05-20 清华大学 一种确定用户属性信息的方法以及信息推送方法
CN110659427A (zh) * 2019-09-06 2020-01-07 北京百度网讯科技有限公司 基于多源数据的城市功能划分方法、装置和电子设备
CN112364776A (zh) * 2020-11-12 2021-02-12 中国地质大学(武汉) 一种城市功能区识别模型建立、分类方法及装置

Also Published As

Publication number Publication date
CN106951828B (zh) 2020-02-21

Similar Documents

Publication Publication Date Title
CN106951828A (zh) 一种基于卫星影像和网络的城市区域功能属性的识别方法
CN104239556B (zh) 基于密度聚类的自适应轨迹预测方法
Hu et al. Federated region-learning: An edge computing based framework for urban environment sensing
CN105045858B (zh) 基于投票的出租车载客点推荐方法
CN109409261B (zh) 一种农作物分类方法及系统
CN107247938A (zh) 一种高分辨率遥感影像城市建筑物功能分类的方法
CN111723257B (zh) 基于用水规律的用户画像方法与系统
CN108717676A (zh) 基于多数据融合的不同尺度下的职住空间评价方法及系统
CN113887517B (zh) 基于并行注意力机制的农作物遥感图像语义分割方法
CN112598438A (zh) 一种基于大规模用户画像的户外广告推荐系统及方法
CN112579922B (zh) 一种基于用户轨迹的兴趣点推荐方法
CN109143408A (zh) 基于mlp的动态区域联合短时降水预报方法
CN110245773A (zh) 一种多源实况时空预报因子提取及纳入模式解释应用的方法
Shi et al. Mapping vegetation and land cover in a large urban area using a multiple classifier system
CN116796904A (zh) 一种轨道交通新线客流预测方法、系统、电子设备及介质
CN110457706A (zh) 兴趣点名称选择模型训练方法、使用方法、装置及存储介质
CN110035393A (zh) 移动终端关系的识别方法
CN114462572A (zh) 一种基于时空融合图的空气质量预测方法及装置
CN114329240A (zh) 选址特征筛选方法、装置、电子设备及存储介质
CN113807278A (zh) 一种基于深度学习的土地利用分类及变化预测方法
CN116528282B (zh) 覆盖场景识别方法、装置、电子设备和可读存储介质
CN109063026A (zh) 一种地理信息数据采集方法及系统
CN114707785A (zh) 基于深度学习的农村居民点多尺度空间特征分析方法
CN108399413A (zh) 一种图片拍摄区域识别及地理定位方法及装置
Nice et al. The nature of human settlement: building an understanding of high performance city design

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
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200221

Termination date: 20210222