[发明专利]一种基于卫星影像和网络的城市区域功能属性的识别方法

权利要求书

1.一种基于卫星影像和网络的城市区域功能属性的识别方法，其特征在于该方法包括以下步骤：

(1)在城市区域内布置供网络用户免费使用的路由器，利用路由器收集用户访问网络的时间T_A和离开网络的时间T_L，或者城市区域内网络用户利用手机自动数据收集程序发送访问网络的时间T_A和离开网络的时间T_L，每次用户访问的持续时间记为D，则：

D＝T_L-T_A

记用户第i次访问网络的行为B_i，B_i用一个四元组(T_A，i，T_L，i，D_i，Uid_i)表示，其中Uid_i代表该用户的代码ID，每个用户对应一个单独的代码ID，记所有用户访问的总次数为Q，记所有用户访问行为构成的集合为G_all；

(2)利用电子地图软件下载城市区域的相应卫星影像，卫星影像的大小为256*256像素，每个像素用RGB颜色值表示；

(3)根据上述步骤(1)的用户访问网络时间T_A，离开网络时间T_L以及持续时间D，计算城市区域的网络特征向量V_N，该特征向量V_N包括36维，计算过程如下：

(3-1)根据用户访问网络时间T_A处于一天24小时中的某一小时，将所有用户的访问行为分为24组，记为G[1，2，...，24]，其中G[a]组中的用户访问网络时间T_A满足a-1≤T_A＜a，各组的用户访问行为的数量作为所述的特征向量V_N的1-24维V_N[1，2，...，24]；

(3-2)根据用户访问网络时间T_A处于周中或周末，将所有用户的访问行为分为2组，记为G[25，26]，其中G[25]组中的用户访问发生在周中，G[26]组中的用户访问发生在周末，各组的用户访问行为数量作为特征向量V_N的25-26维V_N[25，26]；

(3-3)分别计算上述步骤(3-2)中G[25]组和G[26]组中用户访问网络的平均持续时间AvgD[25，26]作为特征向量V_N的27-28维V_N[27，28]：

b＝25，26

(3-4)设置一个用户访问网络的持续时间阈值D_th，计算用户访问网络持续时间大于该阈值的用户访问行为的数量占所有用户访问行为总数量的比率R_long作为特征向量V_N的29维V_N[29]：

其中，I(·)为指示函数，当条件·为真时，I(·)取值为1，当条件·为假时，I(·)取值为0；

(3-5)将G[25]组，G[26]组以及G_all中的用户数量N_u作为特征向量V_N的30-32维V_N[30，31，32]；

(3-6)分别计算G[25]组用户访问行为、G[26]组用户访问行为和G_all用户访问行为的平均重连时间T_re作为特征向量V_N的33-35维V_N[33，34，35]，平均重连时间T_re的计算步骤，以G[25]为例，介绍如下，：

(3-6-1)根据用户的代码IDUid，将G[25]组的用户访问行为分成N_u个部分，N_u为与该组用户访问行为相对应的用户数量，即每个部分中的所有用户访问行为对应同一个Uid，对于N_u个部分中的第c部分用户访问行为，将上述步骤(1)的访问网络时间T_A从小到大排序，得到序列B[1，2，...，p]，其中p为第c部分用户访问行为的数量，则第c部分用户访问行为的平均重连时间T_re，c为：

其中Maximum为设定的平均重连时间T_re极大值，取值为168小时；

(3-6-2)利用下式计算G[25]组用户访问行为的平均重连时间T_re为：

(3-6-3)重复步骤(3-6-1)和步骤(3-6-2)，分别计算G[26]组用户访问行为和G_all用户访问行为的平均重连时间；

(3-7)根据上述步骤(1)的访问网络时间T_A，将所有用户访问行为按照日期顺序记为G′[1，2，...，N_D]，其中每组用户访问行为发生在同一天，计算G′[1，2，...，N_D]中第e天的用户群变化率CRC_e：

其中N′_u，e为与G′[e]组访问行为相对应的用户数量，G′[e]为第e天的的用户访问行为，com(e-1，e)为与G′[e-1]组访问行为和G′[e]组访问行为相对应的用户中共同重复的用户数量；

利用下式计算平均用户群变化率CRC作为特征向量V_N的36维V_N[36]：

(4)将上述步骤(2)的卫星影像中输入到卷积神经网络中，将卷积神经网络的第二个全连接层的输出作为影像特征向量V_I[1，2，...，4096]；

(5)将上述步骤(3)得到的城市区域的网络特征向量V_N[1，2，...，36]和上述步骤(4)得到的城市区域的卫星影像特征向量V_I[1，2，...，4096]连接，合并成向量V[1，2，...，4132]，作为城市区域的特征向量；

(6)建立功能属性分类模型，过程如下：

(6-1)将具有已知功能属性的N个城市区域作为训练样本，记N个城市区域中任意一个城市区域为R_f，1≤f≤N，具有已知的功能属性t_f，利用上述步骤(5)得到特征向量V_f，遍历N个城市区域中的所有城市区域，得到具有已知功能属性的N个城市区域的特征向量；

(6-2)利用AdaBoost算法，训练基于决策树的功能属性分类模型步骤如下：

(6-2-1)定义一个N维样本权重向量w_N，权重向量w_N中的每个元素对应一个具有已知功能属性的城市区域训练样本，将每个元素初始化为

(6-2-2)将误差函数ε_m最小化，得到决策树函数y_m(V)：

其中，I(·)为指示函数，当条件·为真时，I(·)取值为1，当条件·为假时，I(·)取值为0，m为迭代次数，V为上述步骤(5)中的城市区域的特征向量，V_f为上述步骤(6-1)中的N个城市区域中任意一个城市区域R_f的特征向量；

利用下式，计算决策树函数y_m(V)的权重α_m：

利用下式调整权重向量w_N：

其中，Z_m是归一化因子，保证权重向量w_N的和为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)的功能属性分类模型Y_M(V)，输出即为城市区域R的功能属性。