买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：北京交通大学

摘要：本发明涉及一种基于轨迹数据的不确定性医疗可达性计算方法，具体包括如下步骤：步骤1，数据获取和预处理；步骤2，空间处理与网格生成；步骤3，提取历史出行时间；步骤4，构建OD矩阵；步骤5，计算出行完成概率；步骤6，计算可靠性约束下的可达性；本发明提供的方法，具有以下优点：1得到更加符合现实情况的医疗可达性；2模型简单，易于理解和计算，应用型较强；3相比于正方形，六边形具有更少的邻近图形、更小的边面积比；同时，具有各向同性的几何性质，能够减少可达性计算中的误差。

主权项：1.一种基于轨迹数据的不确定性医疗可达性计算方法，其特征在于，所述方法步骤如下：步骤1，数据获取和预处理；首先，通过使用互联网地图公司提供的API接口爬取医疗设施POI数据；然后，剔除因建筑遮蔽、GPS信号丢失而产生的异常信息，并对时间字段进行标准化处理；同时，以轨迹为单位，计算相邻轨迹点之间的时间差、轨迹点与起点之间的时间差、轨迹的行驶距离和平均速度，完成轨迹数据的预处理工作，所述轨迹点与起点之间的时间差为行程时间；步骤2，空间处理与网格生成；在ArcGIS中，设置新建数据框的地理坐标系及空间投影，添加研究区域的行政区划图和城市主要道路，利用ThiessenPolygon工具箱将研究区域划分为若干边长确定的六边形网格；添加步骤1处理完成的轨迹数据，按照位置将六边形网格与轨迹数据进行空间连接；添加医疗设施POI并与六边形网格进行空间连接，参考医疗设施的床位数，设置不同类型医疗设施的权重，计算得到每个网格的医疗设施得分，得到公式1：设Gj为格子j的医疗设施得分，为格子j内设施h的权重，则有： 式中，s为格子j内的医疗设施总数；步骤3，提取历史行程时间；在步骤2处理得到的轨迹数据的基础上，提取每条轨迹上所有的OD出行，表示为i，j，假设轨迹τ由k个轨迹点组成，Qτ表示轨迹τ经过的格子编号的集合，m表示轨迹τ经过的格子数，得到公式2： 式中，表示轨迹τ上的轨迹点k所在的格子编号，则该轨迹上包含的OD对个数为t1，t2，t3，…，tk-1，tk表示轨迹点1，2，3，…，k-1，k的打点时刻；假设轨迹τ在格子i中的轨迹点数量为ni，得到公式3： 则Tdi＝maxtni表示格子i的出发时刻，其值为轨迹τ在格子i中的最后一个点的打点时刻；因此，格子i，j间的行程时间表示为Tdj-Tdi，且Tdj-Tdi＞0；步骤4，构建OD矩阵；在步骤3的基础上，对相同OD的行程时间进行合并，得到该OD出行的历史行程时间矩阵：ODij＝{Tdj-Tdi|i＝0，1，2…，1101，j＝0，1，2，…，1101，i≠j}4其中，ODij为ij之间的历史行程时间集合，0，1，2…，1101为六边形网格编号；步骤5，计算出行完成概率；假设行程时间服从正态分布，则出行完成概率如式5所示： 式中，CDFijT为ij之间行程时间的累积分布函数，σ为ij之间行程时间的样本标准差，μ为ij之间行程时间的样本均值；步骤6，计算可靠性约束下的可达性；累计机会模型表达如式6所示： 式中，CUMi为格子i未考虑可靠性的空间可达性，为格子的机会数，是一个0-1变量，当格子i，j间的行程时间Tdj-Tdi小于等于设定的时间阈值T时，否则在6中引入可靠性约束，则得到7： 式中，PCUMiT，r0为可靠性要求为r0，时间阈值为T时的可达性，为格子的机会数，取值为0或1，当i，j间的行程时间小于等于时间阈值T时，i处居民获得j格子的机会，否则，在时间阈值和可靠性约束下，居民能够获得的机会总数为格子i的可达性。

全文数据：

权利要求：

百度查询： 北京交通大学 一种基于轨迹数据的不确定性医疗可达性计算方法