用于交通需求预测的时空地理加权回归模型及应用 周禹佳,沈辛夷,金盛 (漸江大学建筑工程学院,浙江杭州310058) 摘要:交通需求预测对交通规划具有重要作用,但当对交通需求进行预测,需同 时考虑时空因素对其造成的影响时则具有较大挑战性,本文选取时空地理加权回归 模型(GT瞅R模型〕进行预测,该模型能较好的捕捉交通需求的时空特性以及交通需 求与建成环境之间的关系,运用2016年杭州市车牌数据验证该模型的精确性,本 文运用一般最小二乘回归模型(0LS模型)、地理加权回归模型(GR模型)以及 GTwR模型分别对各个交通小区的交通需求进行预测,结果可以发现:GTwR模型在 需同时考虑时空非平稳性方面有较大优势,且就拟合度方面来说:0LS模型、GWR 模型以及θTR模型的精度分别为12.90%,51.04%以及91.85%,即GTR模型在描 述交通需求预测方面优于传统模型 关键词 交通需求预测;LS模型;GR模型;GTwR模型;时空非平稳性 中图分类号:U238 Geographically and Temporally weighted regression Model for Traffic demand Forecasting and Its application Zhouru-jia, Shen Xin-yi, Jin Sheng College of Civil Enginee ing and Architecture, Zhejiang University, Hangzhou 310058) Abstract: Traffic demand estimation is of great importance to transportation planning. Both the spatial and temporal dependences need to be considered simultaneously, which makes traffic demand estimation challenging. This paper introduces a Geographically and Temporally Weighted Regression(GTWR) model to capture the spatiotemporal characteristics of traffic demand and correlations between traffic demand and built environment. Experimental data, the license plate data in Hangzhou, 2016, is utilized to evaluate the accuracy of the GTwr model. This paper uses the global Ordinary Least Squares(OLS)model, Geographically ighted Regression(GWR)model and GTwR model to predict the traffic demand of each raffic Analysis Zone (TAZ) respectively. The results show that the GTwr model has substantial benefits in modeling both spatial and temporal non-stationarity simultaneously. In the test sample, in terms of goodness-of-fit, 12.90% of the variation in the traffic demand can be explained by the OLs model, while 51.04% by the GwR model, 91.85% by the Gtwr 通讯作者:金盛(1982—),男,浙江温州人,副教授,博士生 (c)1994-2019ChinaAcademicJournalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
用于交通需求预测的时空地理加权回归模型及应用1 周禹佳,沈辛夷,金盛 (浙江大学建筑工程学院，浙江 杭州 310058） 摘 要： 交通需求预测对交通规划具有重要作用，但当对交通需求进行预测，需同 时考虑时空因素对其造成的影响时则具有较大挑战性，本文选取时空地理加权回归 模型（GTWR 模型）进行预测，该模型能较好的捕捉交通需求的时空特性以及交通需 求与建成环境之间的关系，运用 2016 年杭州市车牌数据验证该模型的精确性，本 文运用一般最小二乘回归模型（OLS 模型）、地理加权回归模型（GWR 模型）以及 GTWR 模型分别对各个交通小区的交通需求进行预测，结果可以发现：GTWR 模型在 需同时考虑时空非平稳性方面有较大优势，且就拟合度方面来说：OLS 模型、GWR 模型以及 GTWR 模型的精度分别为 12.90%，51.04%以及 91.85%，即 GTWR 模型在描 述交通需求预测方面优于传统模型。 关键词： 交通需求预测；OLS 模型；GWR 模型；GTWR 模型；时空非平稳性 中图分类号：U238 Geographically and Temporally Weighted Regression Model for Traffic Demand Forecasting and Its Application ZhouYu-jia, ShenXin-yi, JinSheng ( College of Civil Engineering and Architecture,Zhejiang University, Hangzhou 310058） Abstract: Traffic demand estimation is of great importance to transportation planning. Both the spatial and temporal dependences need to be considered simultaneously, which makes traffic demand estimation challenging. This paper introduces a Geographically and Temporally Weighted Regression (GTWR) model to capture the spatiotemporal characteristics of traffic demand and correlations between traffic demand and built environment. Experimental data, the license plate data in Hangzhou, 2016, is utilized to evaluate the accuracy of the GTWR model. This paper uses the global Ordinary Least Squares (OLS) model, Geographically Weighted Regression (GWR) model and GTWR model to predict the traffic demand of each Traffic Analysis Zone (TAZ) respectively. The results show that the GTWR model has substantial benefits in modeling both spatial and temporal non-stationarity simultaneously. In the test sample, in terms of goodness-of-fit, 12.90% of the variation in the traffic demand can be explained by the OLS model, while 51.04% by the GWR model, 91.85% by the GTWR 通讯作者：金盛（1982—），男，浙江温州人，副教授，博士生
model That is, the GTWR model is better than the traditional model in describing traffic demand forecasting key words: Traffic Demand Estimation, OLSModel, GWR Model, GTWR Model, Spatial and Temporal Non-Stationarity 、引言 近年来,随着机动车保有量不断增加,城市交通需求日益增长,仅挖掘现有道路的通行 能力远远不够,应结合交通规划的方法由供给端向需求端改进,对未来道路的交通需求进行 合理的预测从而为道路网规划建设提供较好的支撑,因此作为交通规划理论核心的交通需求 预测,其预测模型的准确度直接影响交通需求预测结果,选择能够较好的指导城市交通规划 实践的交通需求预测模型至关重要 在过去的几年里已经吸引较多的学者研究交通需求预测模型,例如:四阶段法、时间序 列法等。以往的方法较少考虑建成环境对于交通需求的影响,但众所周知,交通小区土地利 用密度及混合使用对其交通需求量具有较大影响,以往有关的研究均说明该种现象:首先在 土地利用密度方面:在考虑了家庭小汽车拥有量及收支情况等因素后,Liu和Shen对美国巴 尔的摩市居住地建成环境与机动车行驶里程之间的关系进行了研究,结果表明居住密度对机 动车出行有显著的总影响。Hong和Shen2,Hong和 Goodchild分别研究了居住密度对 机动车出行的影响,研究结果均表明在控制了家庭收入及小汽车拥有等因素后,居住密度与 机动车出行成负相关。Ding等人的研究证实,工作地的职住关系对于通勤人员机动车出行 距离的影响比居住地的影响更为显著;在土地性质混和使用方面: Acker和wtox分析 了居住地和工作地的建成环境对居民出行方式选择的影响,研究表明随着居住地土地使用混 合度的提高,选择机动车出行的人数显著降低。 Antipova等人分析了美国巴吞鲁日的居住 地的职住平衡对居民机动车出行的影响,研究结果表明较高的职住平衡能够显著降低居民机 动车行驶里程,他们认为较高的职住平衡意味着短距离就业的可能性较大。由此,建成环 境与机动车出行密切相关 近几年随着空间计量模型的发展,在研究交通需求预测时,也将土地利用密度及混合使 用因素进行考量,其中OLS模型是一种较为广泛应用于分析多种因素混杂作用关系的方法 在OLS模型中,假定了各交通小区的数据是相互独立没有空间联系、且在各个小区之间的数 值是固定的。然而,机动化出行需求并不满足独立假定,由于小区与小区之间的差异性,各 交通小区间都存在较大的相互作用。因此,将OLS模型用于此研究的做法由于忽视了空间差 异而有很大弊端。 学者对传统OLS模型进行了改进,纳入了对空间差异性的考虑。其中GW模型应用最为 广泛,GWR通过构建基于距离的权重模型,以分析影响因素与出行需求之间的贡献程度关 系。GWR可以很好地刻画影响因素在空间维度上的作用,因此被广泛应用于交通研究中 但在使用GwR方法进行分析时,时间维度是GWR方法无法处理的范畴。然而时间维度 (c)1994-2019ChinaAcademicJournalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
model.That is, the GTWR model is better than the traditional model in describing traffic demand forecasting. key words: Traffic Demand Estimation, OLSModel, GWR Model, GTWR Model, Spatial and Temporal Non-Stationarity 一、引言 近年来，随着机动车保有量不断增加，城市交通需求日益增长，仅挖掘现有道路的通行 能力远远不够，应结合交通规划的方法由供给端向需求端改进，对未来道路的交通需求进行 合理的预测从而为道路网规划建设提供较好的支撑，因此作为交通规划理论核心的交通需求 预测，其预测模型的准确度直接影响交通需求预测结果，选择能够较好的指导城市交通规划 实践的交通需求预测模型至关重要。 在过去的几年里已经吸引较多的学者研究交通需求预测模型，例如：四阶段法、时间序 列法等。以往的方法较少考虑建成环境对于交通需求的影响，但众所周知，交通小区土地利 用密度及混合使用对其交通需求量具有较大影响，以往有关的研究均说明该种现象：首先在 土地利用密度方面：在考虑了家庭小汽车拥有量及收支情况等因素后，Liu和Shen对美国巴 尔的摩市居住地建成环境与机动车行驶里程之间的关系进行了研究，结果表明居住密度对机 动车出行有显著的总影响[1] 。Hong 和Shen[2] ，Hong 和 Goodchild[3]分别研究了居住密度对 机动车出行的影响，研究结果均表明在控制了家庭收入及小汽车拥有等因素后，居住密度与 机动车出行成负相关。Ding 等人的研究证实，工作地的职住关系对于通勤人员机动车出行 距离的影响比居住地的影响更为显著[4]；在土地性质混和使用方面：Acker 和 Witlox 分析 了居住地和工作地的建成环境对居民出行方式选择的影响，研究表明随着居住地土地使用混 合度的提高，选择机动车出行的人数显著降低[5]。Antipova 等人分析了美国巴吞鲁日的居住 地的职住平衡对居民机动车出行的影响，研究结果表明较高的职住平衡能够显著降低居民机 动车行驶里程，他们认为较高的职住平衡意味着短距离就业的可能性较大[6]。由此，建成环 境与机动车出行密切相关。 近几年随着空间计量模型的发展，在研究交通需求预测时，也将土地利用密度及混合使 用因素进行考量，其中OLS模型是一种较为广泛应用于分析多种因素混杂作用关系的方法。 在OLS模型中，假定了各交通小区的数据是相互独立没有空间联系、且在各个小区之间的数 值是固定的。然而，机动化出行需求并不满足独立假定，由于小区与小区之间的差异性，各 交通小区间都存在较大的相互作用。因此，将OLS模型用于此研究的做法由于忽视了空间差 异而有很大弊端。 学者对传统OLS模型进行了改进，纳入了对空间差异性的考虑。其中GWR模型应用最为 广泛，GWR通过构建基于距离的权重模型，以分析影响因素与出行需求之间的贡献程度关 系。GWR可以很好地刻画影响因素在空间维度上的作用，因此被广泛应用于交通研究中。 但在使用GWR方法进行分析时，时间维度是GWR方法无法处理的范畴。然而时间维度
在交通分析中尤为重要,例如,早晚高峰的出行需求往往比平峰多出50%,这说明出行需求 是具有时间差异性的。出行需求的时间差异特性代表了它对历史数据具有敏感性。因此,出 行需求数据的时空特性分析在交通领域中一直是炙手可热的研究方向。 尽管已有不少研究分析过出行需求的时空特性,这些研究大多着力于对个体出行者的出 行行为研究,少有研究基于交通小区对出行需求的时空特性进行分析。同时,多数在交通小 区基础上分析各种因素对出行需求影响的研究,仅基于GWR模型研究了空间差异性对出行 需求的影响,而少有将时间差异性与空间差异性结合的研究。为了弥补这一缺憾,本文采用 GTwR模型对机动化出行需求进行研究和分析。 本文余下部分研究结构如下所示:方法部分对于GTWR模型进行具体介绍:案例部分描 述了运用杭州市车牌识别数据进行交通需求预测的过程:结果对比与分析部分介绍了三种预 测模型的具体应用过程,并将三种模型预测精度进行对比 二、方法 黄波在GWR方法的基础上提出了考虑时空因素的GTWR模型,其用于对与时空特性相 关的房价问题的建模;随后,GTWR模型被大量运用于有关时空问题的研究中。韩兆洲等 人基于2006-2014年省级时空数据,分别运用时间加权回归模型(TWR)、GWR模型及GTWR 模型分析我国各个地区的最低工资增长率的非平稳性;丁焕峰等人运用GTWR模型探索了 城乡区域间收入水平的时空差异性;马晓磊等人首次将GTWR模型运用于交通领域,分析 建成环境对于公交客流的影响田。 本文应用GIWR模型,根据杭州市各个交通小区的空间位置,各种类型的PO密度以及 居民出行信息进行交通需求预测,模型的具体步骤如下所示: ①获取数据 获取杭州市各个交通小区的地理坐标(u),Y表示i小区的交通需求,t表示采样时间 (本文取1个小时),表示各个交通小区各种类型的PO1密度矩阵。 ②计算每个交通小区与其他交通小区的时空距离 由于时间通常以小时或天为测量单位,空间通常为m或km为测量单位,二者测量标 准的差异会对结果造成较大偏差,因此需要设置参数平衡时空差异,交通小区i和j之间的 时空距离公式如下所示: (d)2=1(n1-1)2+("-")]+(1-1)2 μ和λ为平衡时空差异的参数,μ和λ值的设定对于研究结果具有较大影响 ③构建交通小区之间的时空相关性矩阵W H= diag(W1…,Wn)为对角矩阵,对角元素为两个交通小区之间相互影响的权重值, (c)1994-2019ChinaAcademicJournalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
在交通分析中尤为重要，例如，早晚高峰的出行需求往往比平峰多出50%，这说明出行需求 是具有时间差异性的。出行需求的时间差异特性代表了它对历史数据具有敏感性。因此，出 行需求数据的时空特性分析在交通领域中一直是炙手可热的研究方向。 尽管已有不少研究分析过出行需求的时空特性，这些研究大多着力于对个体出行者的出 行行为研究，少有研究基于交通小区对出行需求的时空特性进行分析。同时，多数在交通小 区基础上分析各种因素对出行需求影响的研究，仅基于GWR模型研究了空间差异性对出行 需求的影响，而少有将时间差异性与空间差异性结合的研究。为了弥补这一缺憾，本文采用 GTWR模型对机动化出行需求进行研究和分析。 本文余下部分研究结构如下所示：方法部分对于GTWR模型进行具体介绍；案例部分描 述了运用杭州市车牌识别数据进行交通需求预测的过程；结果对比与分析部分介绍了三种预 测模型的具体应用过程，并将三种模型预测精度进行对比。 二、方法 黄波在GWR方法的基础上提出了考虑时空因素的GTWR模型，其用于对与时空特性相 关的房价问题的建模[7]；随后，GTWR模型被大量运用于有关时空问题的研究中。韩兆洲等 人基于2006—2014年省级时空数据,分别运用时间加权回归模型(TWR)、GWR模型及GTWR 模型分析我国各个地区的最低工资增长率的非平稳性[8]；丁焕峰等人运用GTWR模型探索了 城乡区域间收入水平的时空差异性[9]；马晓磊等人首次将GTWR模型运用于交通领域，分析 建成环境对于公交客流的影响[11]。 本文应用GTWR模型，根据杭州市各个交通小区的空间位置，各种类型的POI密度以及 居民出行信息进行交通需求预测，模型的具体步骤如下所示： ①获取数据 获取杭州市各个交通小区的地理坐标(u,v)，Yi 表示 i 小区的交通需求，t 表示采样时间 （本文取 1 个小时），Xik表示各个交通小区各种类型的 POI 密度矩阵。 ②计算每个交通小区与其他交通小区的时空距离 由于时间通常以小时或天为测量单位，空间通常为 m 或 km 为测量单位，二者测量标 准的差异会对结果造成较大偏差，因此需要设置参数平衡时空差异，交通小区 i 和 j 之间的 时空距离公式如下所示： 2 2 2 2 ( ) [( ) ( ) ] ( ) i j i j i j ST d u u v v t t （1） μ 和 λ 为平衡时空差异的参数，μ 和 λ 值的设定对于研究结果具有较大影响。 ③构建交通小区之间的时空相关性矩阵 W ( ,...., ) W diag W1 Wn 为对角矩阵，对角元素为两个交通小区之间相互影响的权重值
般来说,与交通小区i距离越近的交通小区,对该小区交通需求的影响越大,时空距离矩 阵一般运用高斯函数构建: (d)2 W=exp(-=12) ④对带宽h,时空参数μ和进行标定 为了简化上述模型,令τ=uA,则距离公式可以表示为 (d)2 之=4-4)+(-1)]+(-1)2 (3) 利用交叉验证计算残差平方和,绘制以h为横坐标,以残差平方和为纵坐标来确定参数 h的最优值。 CRSS(h)=∑(y-y2(h) ⑤回归参数估计 GTWR模型基本形式如下所示 Y=B(u1,",)+∑B(u1,,)X从+E 其中,G为误差项,(xy2b为不同类型PO对于交通小区交通需求的时空相关程度, 使用最小二乘估计方法,通过使残差平方和最小来估计相关系数 因此,各种类型POI密度与交通需求之间的参数估计可以表示为: B(u v, =[Xw(u v,oXrY'w(u, v, L, )Y 三、案例分析 在本部分,杭州市是我们选择验证GTWR模型是否适用于交通需求预测研究的城市, 我们应用杭州市车牌数据及杭州市路网进行探究,将杭州市划分为126个交通小区,并且建 立杭州市交通需求影响因素的数据集 )车牌数据 车牌数据相对于其他数据来说具有样本量大和精确度高的特点,车牌数据包括车辆的 车牌号、运行时间、地点以及车辆类别等,因此我们能从数据中获取所有车辆的出行链。杭 州市共有1399个检测器,且检测器覆盖率在90%以上,因此据此获取的车牌数据能够较完 整的描述杭州市机动化出行,表1显示了检测设备的记录表,检测设备的详细经纬度坐标便 (c)1994-2019ChinaAcademicJournalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
一般来说，与交通小区 i 距离越近的交通小区，对该小区交通需求的影响越大，时空距离矩 阵一般运用高斯函数构建： exp( ) 2 2 h d W ST ij （ ） （2） ④对带宽 h，时空参数 μ 和 λ 进行标定 为了简化上述模型，令 τ=μ/λ，则距离公式可以表示为 2 2 2 2 [( ) ( ) ] ( ) ( ) i j i j i j ST u u v v t t d （3） 利用交叉验证计算残差平方和，绘制以 h 为横坐标，以残差平方和为纵坐标来确定参数 h 的最优值。 2 1 CVRSS(h) ( y y (h)) i i （4） ⑤回归参数估计 GTWR 模型基本形式如下所示： i i i ik i k Yi (ui ,vi ,ti) k (u ,v ,t )X 0 （5） 其中，εi 为误差项，βk(xi,yi,ti)为不同类型 POI 对于交通小区交通需求的时空相关程度， 使用最小二乘估计方法，通过使残差平方和最小来估计相关系数。 因此，各种类型 POI 密度与交通需求之间的参数估计可以表示为： u v t X W u v t X X W ui vi ti Y T i i i T i i i ( , , ) [ ( , , ) ] ( , , ) 1 （6） 三、案例分析 在本部分，杭州市是我们选择验证 GTWR 模型是否适用于交通需求预测研究的城市， 我们应用杭州市车牌数据及杭州市路网进行探究，将杭州市划分为 126 个交通小区，并且建 立杭州市交通需求影响因素的数据集。 （一）车牌数据 车牌数据相对于其他数据来说具有样本量大和精确度高的特点[12]，车牌数据包括车辆的 车牌号、运行时间、地点以及车辆类别等，因此我们能从数据中获取所有车辆的出行链。杭 州市共有 1399 个检测器，且检测器覆盖率在 90%以上，因此据此获取的车牌数据能够较完 整的描述杭州市机动化出行，表 1 显示了检测设备的记录表，检测设备的详细经纬度坐标便
于区分车辆所在的地点。表2展示了车牌数据的相关信息表 表1设备点位记录表 设备表主要信息 数据形式 数据说明 DEV ID 50091 备点位编码 DEV NAME 玉古路求是路北口 设备点位名称 120.127709,30.26547 备点位经纬度坐标 表2过车记录表 车牌识别主要信息 数据形式 数据说明 过车记录编号 DEV_ID 2147726 采集设备点位编号 DEV NAME 中河高架西湖大道东下匝 采集设备点位名称 WAY ID 车辆匹配时所处的车道 车牌识别数据,存在未识别和无牌 CAR TYPE 车辆类型,不同数字代表不同类型 CAP DATE 2016/6/130:09:1 采集时间,精确到秒 (二)建成环境 随着计算机技术、数据传输等技术的不断发展,各种交通数据传感器应运而生,使得研 究者能够大范围、自动化地获取机动车出行路径,而且越来越多的网络地图数据也能为研究 所应用,比如:百度地图AP12,谷歌地图AP1,高德地图开放平台等,本文利用高德 地图开放平台并根据高德地图自定义编码表爬取所需类型的POI数据,见表3:每种类型的 POI记录都包含该种类型PO1的名称以及经纬度坐标,运用 AEGIS软件得到每个交通小区 各种类型POI的密度,第k种类型POI密度的计算公式如下所示 Vk第k种类型的PO1数量/交通小区i的面积 (7) 据此,交通小区i的POl向量被表示为x=(v,2…,vF,此处F表示为POl类别的数量。 表3P0|类别 Po|种类 Po|种类 院 银行和金融机构 旅店 14 生活服务 停车场 789 对外交通 政府机构 公交地铁站 便利店 咖啡厅 (三)交通需求预测 我们利用杭州市车牌数据计算交通需求,根据车牌数据中的车牌号以及检测时间,可以 得到每辆机动车的一系列出行点,因此每辆车的出行链被获取,通过识别出行链中每个相邻 点之间时空关系,我们能够识别出行链中的间断点,根据检测器编号的地理信息,将其聚类 (c)1994-2019ChinaAcademicJournalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
于区分车辆所在的地点。表 2 展示了车牌数据的相关信息表。 表 1 设备点位记录表 设备表主要信息 数据形式 数据说明 DEV_ID 2150091 设备点位编码 DEV_NAME 玉古路求是路北口 设备点位名称 XY 120.127709,30.265471 设备点位经纬度坐标 表 2 过车记录表 车牌识别主要信息 数据形式 数据说明 ID 8071219675 过车记录编号 DEV_ID 2147726 采集设备点位编号 DEV_NAME 中河高架西湖大道东下匝道 采集设备点位名称 WAY_ID 1 车辆匹配时所处的车道 CAR_NUM 浙 AT5738 车牌识别数据，存在未识别和无牌 CAR_TYPE 1 车辆类型，不同数字代表不同类型 CAP_DATE 2016/6/13 0:09:15 采集时间，精确到秒 （二）建成环境 随着计算机技术、数据传输等技术的不断发展，各种交通数据传感器应运而生，使得研 究者能够大范围、自动化地获取机动车出行路径，而且越来越多的网络地图数据也能为研究 所应用，比如：百度地图 API[12]，谷歌地图 API[13]，高德地图开放平台[14]等，本文利用高德 地图开放平台并根据高德地图自定义编码表爬取所需类型的 POI 数据，见表 3：每种类型的 POI 记录都包含该种类型 POI 的名称以及经纬度坐标，运用 ACGIS 软件得到每个交通小区 各种类型 POI 的密度，第 k 种类型 POI 密度的计算公式如下所示： Vk=第 k 种类型的 POI 数量/交通小区 i 的面积 (7) 据此，交通小区 i 的 POI 向量被表示为 xi=(v1,v2,...,vF),此处 F 表示为 POI 类别的数量。 表 3 POI 类别 编号 POI 种类 编号 POI 种类 1 医院 11 购物中心 2 学校 12 饭店 3 银行和金融机构 13 景区 4 旅店 14 运动 5 住宅 15 影剧院 6 生活服务 16 停车场 7 对外交通 17 政府机构 8 公交地铁站 18 便利店 9 公司 19 咖啡厅 10 娱乐设施 （三）交通需求预测 我们利用杭州市车牌数据计算交通需求，根据车牌数据中的车牌号以及检测时间，可以 得到每辆机动车的一系列出行点，因此每辆车的出行链被获取，通过识别出行链中每个相邻 点之间时空关系，我们能够识别出行链中的间断点，根据检测器编号的地理信息，将其聚类
在已经划分好的126个交通小区中,即可获取每个交通小区在每个时间段内的交通需求。各 个交通小区在早晚高峰期间的交通需求计算结果如图1所示: mm on g amJ (17pm 图1交通需求的时空分布 (四)时空非平稳性判断 在应用GTwR模型之前,首先应分析GIWR模型相对于OLS与GWR模型是否能较 好地对交通需求进行预测,意味着我们应先判断研究的问题是否具有明显的时空非平稳性, 在早期的研究中, Brunsdon等人提出运用假设检验的方法来解决,在本文中,基于方差 计算的F检验被运用,结果如表4所示:GWR模型相对于OLS模型的残差平方和有所改善, 表明交通需求预测具有空间非平稳性,进一步将运用GIWR模型与GWR模型的结果进行 比较可以发现该问题也具有时间非平稳性。因此,运用基于GWR模型对交通需求问题进行 预测是合适的,而且,从结果可以看出运用传统的GWR模型描述时空非平稳性的问题是不 合适的,因此,我们提出了在GWR模型的基础上考虑了时间因素的更精确的GTWR模型 表4在GwR/叫LS,GTwR/GN模型之间的方差分析对比 P值 OLS残差 2351400671 117570033.56 GWR残差 1321361282513.4252571246 TWR残差 219834211.23 2537.42 86636.90 (c)1994-2019ChinaAcademicJournalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
在已经划分好的 126 个交通小区中，即可获取每个交通小区在每个时间段内的交通需求。各 个交通小区在早晚高峰期间的交通需求计算结果如图 1 所示： 图 1 交通需求的时空分布 （四）时空非平稳性判断 在应用 GTWR 模型之前，首先应分析 GTWR 模型相对于 OLS 与 GWR 模型是否能较 好地对交通需求进行预测，意味着我们应先判断研究的问题是否具有明显的时空非平稳性， 在早期的研究中，Brunsdon 等人提出运用假设检验的方法来解决[15]，在本文中，基于方差 计算的 F 检验被运用，结果如表 4 所示：GWR 模型相对于 OLS 模型的残差平方和有所改善， 表明交通需求预测具有空间非平稳性，进一步将运用 GTWR 模型与 GWR 模型的结果进行 比较可以发现该问题也具有时间非平稳性。因此，运用基于 GWR 模型对交通需求问题进行 预测是合适的，而且，从结果可以看出运用传统的 GWR 模型描述时空非平稳性的问题是不 合适的，因此，我们提出了在 GWR 模型的基础上考虑了时间因素的更精确的 GTWR 模型。 表 4 在 GWR/OLS, GTWR/GWR 模型之间的方差分析对比 变量 RSS DF MS F 值 P 值 OLS 残差 2351400671.28 20 117570033.56 GWR 残差 1321336213.28 2513.42 525712.46 GTWR 残差 219834211.23 2537.42 86636.90
GWR/OLS改善 10300945460052.5819591541.1337 GTWR/OLS改善 2131566050 76.58 27835596.33321.29 GTWR/GWR改善 1101502002.05 45895916.75529.750.00 四、结果对比和分析 )OLS模型结果 运用OLS模型对杭州市2016年6月1日各个交通小区交通需求进行预测,结果如表5所示, 结果显示进行交通需求预测可以利用建成环境建模,根据t检验结果,有13个决策变量在以 95%为置信区间中具有统计学意义,而且,风景名胜、对外交通等类型PO与交通需求成负 相关,而公司和饭店等类型的PO与交通需求成正相关,即公司与饭店等类型的PO的密度 越大,风景区以及对外交通等类型的PO密度越小,交通需求越大 表50Ls模型参数估计总结 参数 系数 744.044 28.735 教育 312 -2.190 0.029 银行与金融机构 0.000 916 055* 住宅 -0.010 0.560 0.57 生活服务 0.122 0.000 公交地铁站 0.93 0.352* 娱乐 -5.782 物中心 0.635 0.525* 4.83 风景名胜 0.122* 0.000 0.857* 判断信息 0.1290 标准残差 949,532 残差平方和 2351400671.278 (二)GWR与GTWR模型结果 本部分运用GWR和GTwR模型进行交通需求的估计,结果如表6所示。估计系数表 示每种类型的POI密度对于交通需求预测结果的影响程度,表6给出了估计系数有特征性 的值,LQ代表下四分位数,UQ代表上四分位数 表6也给出了详细的统计对比,可以看出:OLS模型、GWR模型与 GTWR模型的R2 分别为12.90%、51.0%以及91.85%。我们推断因为GTwR模型能够更好地处理时空非平稳 (c)1994-2019ChinaAcademicJournalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
GWR/OLS 改善 1030064458.00 52.58 19591541.13 37.27 0.00 GTWR/OLS 改善 2131566460.05 76.58 27835596.33 321.29 0.00 GTWR/GWR 改善 1101502002.05 24 45895916.75 529.75 0.00 四、结果对比和分析 （一）OLS 模型结果 运用OLS模型对杭州市2016年6月1日各个交通小区交通需求进行预测，结果如表5所示， 结果显示进行交通需求预测可以利用建成环境建模，根据t检验结果，有13个决策变量在以 95%为置信区间中具有统计学意义，而且，风景名胜、对外交通等类型POI与交通需求成负 相关，而公司和饭店等类型的POI与交通需求成正相关，即公司与饭店等类型的POI的密度 越大，风景区以及对外交通等类型的POI密度越小，交通需求越大。 表 5 OLS 模型参数估计总结 参数 系数 T 值 P 值 截距 744.044 28.735 0.000 医院 -0.443 -3.171 0.002 教育 -0.312 -2.190 0.029 银行与金融机构 - 0.381 -5.319 0.000 旅店 -0.091 -1.916 0.055* 住宅 -0.010 -0.560 0.575* 生活服务 0.122 4.909 0.000 对外交通 -1.839 -5.223 0.000 公交地铁站 0.362 0.931 0.352* 公司 0.235 7.732 0.000 娱乐 -1.677 -5.782 0.000 购物中心 0.025 0.635 0.525* 饭店 0.203 4.831 0.000 风景名胜 -0.378 -4.841 0.000 运动 0.389 1.547 0.122* 影剧院 -0.560 -0.574 0.566* 停车场 -0.190 -3.581 0.000 政府机构 -0.200 -5.174 0.000 便利店 0.992 4.599 0.000 咖啡厅 -0.021 -0.180 0.857* 判断信息 R 2 0.1290 标准残差 949.532 残差平方和 2351400671.278 AIC 16.564 （二）GWR 与 GTWR 模型结果 本部分运用 GWR 和 GTWR 模型进行交通需求的估计，结果如表 6 所示。估计系数表 示每种类型的 POI 密度对于交通需求预测结果的影响程度，表 6 给出了估计系数有特征性 的值，LQ 代表下四分位数，UQ 代表上四分位数。 表 6 也给出了详细的统计对比，可以看出：OLS 模型、GWR 模型与 GTWR 模型的 R2 分别为 12.90%、51.04%以及 91.85%。我们推断因为 GTWR 模型能够更好地处理时空非平稳
性事件,而交通需求预测问题是明显的时空问题 表6R与R模型参数估计总结 参数 GWR(带宽=0.3756) GTWR(带宽=0.237 LQ 截距 0.0 0.00.00.0 0.00.00.00.4 医院 143.3-6.40.14.9405.1 4.7|-3.80.03.8721.7 教育 64.6 0.04.5160.2-113.4-4.60.03.4270.2 银行与金融-47.8-6.60.01.9 -65.7-4.00.02.3102.2 机构 旅店-125.2-2.30.23.4164.7-230.7-2.20.13.11120.5 住宅 30.1-1.00.25.5587-68 0.70.22.3 79,6 生活服务66.0-2.60.42.996.0-67.3-2.60.32.8179.7 对外交通-161.70.00.00.01526-53.8-1.20.00.0845.9 公交地铁站-2630.13.116.64594-10760.12.814.2764 公司 娱乐-50.3-940.26.11163-6:5-3.510.144405.8 购物中 0.70.9 饭店 4.1 35 风景名胜-63.5-2.40.03.715.9-2450-2.50.027158.3 124.5-8.00.04.4164.4-81.5-5.50.02.3 影剧院 0.60.0 -470. 0.0 786 停车场-487-3.20.57.5743-70.9-2.50.35.2 政府机构-75.5-3.50.33 939-2.60.2 便利店-131.2-1.01.614.212.-156.3-0.11.912.4456.7 咖啡厅|-275.3-3.20.110.04.6-313.3-2.90.14.1315.5 判断信息 0.5104 0.9185 残差 711,792 290.331 平方和 1321336213.278 219834211.231 AIC 15.988 14.194 说明:*参数x10 为了避免过拟合现象的发生,我们随机的删除总体样本中一部分数据验证模型的精度, 结果如表7所示,可以看出:数据的随机缺失没有显著降低模型的精度,表明GIwR模型 能较好地适用于交通需求的预测 表7不同数据比例的模型精度对比 数据比例 70% 60% 0.918 0.918 0.9125 0.8859 0.9140 五、结论 以前的研究显示,交通小区的交通需求不随时间和空间变化的假设是不成立的,本文采 取了一种分析杭州市各个交通小区交通需求的非静态的方法,在对各个交通小区的交通需求 进行预测时同时包含了时空因素,并将其与没有包含时空因素的OLS方法及仅包含空间因 素的GWR方法进行对比,结果显示OLS,GWR,GTWR模型的拟合度分别为12.90%,51.04% 以及91.85%。因此我们得出GTWR模型更适合对交通小区的发生量进行建模。 然而我们的研究还有一些缺陷,在模型建立过程中,仅使用一个参数平衡时空关系,且 (c)1994-2019ChinaAcademicJournalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
性事件，而交通需求预测问题是明显的时空问题。 表 6 GWR 与 GTWR 模型参数估计总结 参数 GWR (带宽=0.3756) GTWR (带宽=0.2337) Min LQ Med UQ Max Min LQ Med UQ Max 截距 0.0 0.0 0.0 0.0 0.3 0.0 0.0 0.0 0.0 0.4 医院 -143.3 -6.4 0.1 4.9 405.1 -74.7 -3.8 0.0 3.8 721.7 教育 -64.6 -9.5 0.0 4.5 160.2 -113.4 -4.6 0.0 3.4 270.2 银行与金融 机构 -47.8 -6.6 0.0 1.9 83.5 -65.7 -4.0 0.0 2.3 102.2 旅店 -125.2 -2.3 0.2 3.4 164.7 -230.7 -2.2 0.1 3.1 120.5 住宅 -30.1 -1.0 0.2 5.5 58.7 -68.9 -0.7 0.2 2.3 79.6 生活服务 -66.0 -2.6 0.4 2.9 96.0 -67.3 -2.6 0.3 2.8 179.7 对外交通 -161.7 0.0 0.0 0.0 1526 -533.8 -1.2 0.0 0.0 845.9 公交地铁站 -266.3 0.1 3.1 16.6 459.4 -1076 0.1 2.8 14.2 768.4 公司 -39.5 0.1 1.5 4.1 126.1 -55.2 0.0 1.2 4.4 236.0 娱乐 -50.3 -9.4 0.2 6.1 116.3 -61.5 -3.5 0.1 4.4 405.8 购物中心 -34.9 -0.7 0.9 5.3 51.7 -66.9 -1.1 0.7 4.7 77.9 饭店 -42.6 -4.1 0.2 8.1 104.9 -115.9 -3.7 0.1 3.5 135.6 风景名胜 -63.5 -2.4 0.0 3.7 115.9 -245.0 -2.5 0.0 2.7 158.3 运动 -124.5 -8.0 0.0 4.4 164.4 -81.5 -5.5 0.0 2.3 105.4 影剧院 -544.5 -0.6 0.0 5.1 320.3 -470.7 -0.4 0.0 4.0 1786 停车场 -48.7 -3.2 0.5 7.5 74.3 -70.9 -2.5 0.3 5.2 139.2 政府机构 -75.5 -3.5 0.3 3.1 52.6 -93.9 -2.6 0.2 3.4 97.2 便利店 -131.2 -1.0 1.6 14.2 112.9 -156.3 -0.1 1.9 12.4 456.7 咖啡厅 -275.3 -3.2 0.1 10.0 44.6 -313.3 -2.9 0.1 4.1 315.5 判断信息 R 2 0.5104 0.9185 标准残差 711.792 290.331 残差平方和 1321336213.278 219834211.231 AIC 15.988 14.194 说明：*参数 x10-4 . 为了避免过拟合现象的发生，我们随机的删除总体样本中一部分数据验证模型的精度， 结果如表 7 所示，可以看出：数据的随机缺失没有显著降低模型的精度，表明 GTWR 模型 能较好地适用于交通需求的预测。 表 7 不同数据比例的模型精度对比 数据比例 100% 90% 80% 70% 60% R 2 0.9185 0.918 0.9125 0.8859 0.9140 五、结论 以前的研究显示，交通小区的交通需求不随时间和空间变化的假设是不成立的，本文采 取了一种分析杭州市各个交通小区交通需求的非静态的方法，在对各个交通小区的交通需求 进行预测时同时包含了时空因素，并将其与没有包含时空因素的 OLS 方法及仅包含空间因 素的 GWR 方法进行对比，结果显示 OLS,GWR,GTWR 模型的拟合度分别为 12.90%，51.04% 以及 91.85%。因此我们得出 GTWR 模型更适合对交通小区的发生量进行建模。 然而我们的研究还有一些缺陷，在模型建立过程中，仅使用一个参数平衡时空关系，且
仅使用一般的线性回归的方法得到交通小区的交通需求的预测值,因此更多精确有效的研究 方法需要进一步研究。 参考文献 [1] Liu, C,& Shen, Q.(2011). An empirical analysis of the influence of urban form or household travel and energy consumption. Computers, Environment and Urban Systems 35(5),347-357 2 Hong, J, Shen, Q(2013). Residential density and transportation emissions: examining the connection by addressing spatial autocorrelation and self-selection. Transportation Research Part D: Transport and Environment, 22, 75-79 3 Hong, J.,& Goodchild, A(2014). Land use policies and transport emissions: Modeling the impact of trip speed, vehicle characteristics and residential location. Transportation Research Part D: Transport and Environment, 26, 47-51 []Ding, C, Wang, Y,, Xie, B,& Liu, C.(2014). Understanding the role of built environment in reducing vehicle miles traveled accounting for spatial heterogeneity Sustainability, 6(2), 589-601 5 Van Acker, V,& witlox, F(2011). Commuting trips within tours: how is commuting related to land use?. Transportation, 38(3), 465-486 6 Antipova, A, Wang, F,& wilmot, C.(2011). Urban land uses, socio-demographic attributes and commuting: A multilevel modeling approach. Applied Geography, 31(3), 1010-1018 [7 Huang, B, Wu, B,& Barry, M.(2010). Geographically and temporally weighted regression for modeling spatio-temporal variation in house prices. International Journal of Geographical Information Science, 24(3), 383-401 [8]韩兆洲,林仲源.我国最低工资增长机制时空非平稳性测度研究[].统计研究, 2017(06):40-53 [9]丁焕峰,刘心怡.城镇化背景下城乡收入差距的时空演化[J.经济地理,2017(4). [10] Ma, X, Zhang, J, Ding, C,& Wang, Y,(2018). A geographically and temporally weighted regression model to explore the spatiotemporal influence of built environment on transit ridership. Computers, Environment and Urban Systems, 70, 113-124 [11]Sun, J, Feng, Y(2011). A new method of OD estimation based on automatic vehicle dentification data. Journal of Tongji University. Natural Science, 39(12), 1800-1804 [12]baiduMapOpenPlatformhttp://ibSyun.baidu.com.AccessedJun23,2018 [13]googleMapsApis.https://developers.googlecom/maps.acCessedJun28,2018 [14amap.http://ibS.amap.com.AccessedJuly3,2018 [15 Brunsdon, C, Fotheringham, A. S,& Charlton, M.(1999). Some notes on parametric ignificance tests for geographically weighted regression. Journal of Regional science 39(3),497-524 作者简介 周禹佳,工作单位:浙江大学建筑工程学院智能交通研究所,电话:08657188208968, (c)1994-2019ChinaAcademicJournalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
仅使用一般的线性回归的方法得到交通小区的交通需求的预测值，因此更多精确有效的研究 方法需要进一步研究。 参考文献 [1] Liu, C., & Shen, Q. (2011). An empirical analysis of the influence of urban form on household travel and energy consumption. Computers, Environment and Urban Systems, 35(5), 347-357. [2] Hong, J., & Shen, Q. (2013). Residential density and transportation emissions: examining the connection by addressing spatial autocorrelation and self-selection. Transportation Research Part D: Transport and Environment, 22, 75-79. [3] Hong, J., & Goodchild, A. (2014). Land use policies and transport emissions: Modeling the impact of trip speed, vehicle characteristics and residential location. Transportation Research Part D: Transport and Environment, 26, 47-51. [4] Ding, C., Wang, Y., Xie, B., & Liu, C. (2014). Understanding the role of built environment in reducing vehicle miles traveled accounting for spatial heterogeneity. Sustainability, 6(2), 589-601. [5] Van Acker, V., & Witlox, F. (2011). Commuting trips within tours: how is commuting related to land use?. Transportation, 38(3), 465-486. [6] Antipova, A., Wang, F., & Wilmot, C. (2011). Urban land uses, socio-demographic attributes and commuting: A multilevel modeling approach. Applied Geography, 31(3), 1010-1018. [7] Huang, B., Wu, B., & Barry, M. (2010). Geographically and temporally weighted regression for modeling spatio-temporal variation in house prices. International Journal of Geographical Information Science, 24(3), 383-401. [8] 韩兆洲, 林仲源. 我国最低工资增长机制时空非平稳性测度研究[J].统计研究, 2017(06):40-53. [9] 丁焕峰, 刘心怡. 城镇化背景下城乡收入差距的时空演化[J]. 经济地理, 2017(4). [10] Ma, X., Zhang, J., Ding, C., & Wang, Y. (2018). A geographically and temporally weighted regression model to explore the spatiotemporal influence of built environment on transit ridership. Computers, Environment and Urban Systems, 70, 113-124. [11] Sun, J., & Feng, Y. (2011). A new method of OD estimation based on automatic vehicle identification data. Journal of Tongji University. Natural Science, 39(12), 1800-1804. [12] Baidu Map Open Platform. http://Ibsyun.baidu.com. Accessed Jun 23, 2018. [13] Google Maps APIs. https://developers.google.com/maps. Accessed Jun 28, 2018. [14] AMAP. http://Ibs.amap.com. Accessed July 3, 2018. [15] Brunsdon, C., Fotheringham, A. S., & Charlton, M. (1999). Some notes on parametric significance tests for geographically weighted regression. Journal of Regional science, 39(3), 497-524. 作者简介 周禹佳，工作单位：浙江大学建筑工程学院智能交通研究所，电话：086 571 88208968
传真:08657188208968,电子信箱:2396256725 gg. com 沈辛夷,工作单位:浙江大学建筑工程学院智能交通研究所,电话:08657188208968, 传真:08657188208968,电子信箱:shenxi@ziu.edu.cn 金盛,工作单位:浙江大学建筑工程学院智能交通研究所,电话:08657188208968, 传真:08657188208968,电子信箱:jinsheng@ziu.edu.cn (c)1994-2019ChinaAcademicJournalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
传真：086 571 88208968，电子信箱：2396256725@qq.com 沈辛夷，工作单位：浙江大学建筑工程学院智能交通研究所，电话：086 571 88208968， 传真：086 571 88208968，电子信箱：shenxy@zju.edu.cn 金盛，工作单位：浙江大学建筑工程学院智能交通研究所，电话：086 571 88208968， 传真：086 571 88208968，电子信箱：jinsheng@zju.edu.cn
