本发明涉及一种基于服务时间的区域共享自行车有效车辆数获取方法，属于城市公共自行车和共享单车运营
技术领域：
。
背景技术：
：：随着我国城市化不断发展，城市人口增加，机动化水平提高，带来了拥堵、环境，能源等一系列问题，同时居民出行需求也存在多样化和个性化的发展趋势。而以公共自行车和共享单车为代表的共享自行车，作为城市公共交通网络的补充，较好解决“最后一公里”问题，满足中短距离的差异化出行需求，有助于形成“多模式、一体化、人性化”的城市综合交通系统。公共自行车和共享单车的区别主要在于信息共享设备分别是车桩和智能锁，在服务对象、出行距离和使用功能等方面都具有明显的一致性，两者可以统称为共享自行车。现阶段共享自行车服务已经在世界各地得到广泛推广，但是规划管理和运营调度主要是依靠实地经验和传统理论，未能利用共享自行车系统的运营服务数据来进行更加科学高效的经营管理。在利用运营数据精细化指导企业经营的过程中，一个关键问题就是一个区域内的共享自行车有效车辆数量问题。因为共享自行车服务的流动性本质，所以一个区域内车辆数量是处在实时动态变化的过程中，这一点不利于计算不同区域的车辆数量、车辆密度、车辆每日使用次数等运营指标，阻碍了企业基于不同区域运营指标分析来实现精细化经营管理。目前尚无确定区域内共享自行车实际服务车辆数的有效方法。技术实现要素：本发明所要解决的技术问题是克服共享自行车流动性对服务车辆数量统计带来的困难，提供一种基于服务时间的区域共享自行车有效车辆数获取方法。本发明划分共享自行车分析区域，利用企业的运营数据提取出车辆的时空定位数据，根据车辆在区域内的服务时间，确定各区域的有效车辆数量，为共享自行车精细化运营调度奠定指标分析基础。本发明为解决上述技术问题采用以下技术方案：本发明提供一种基于服务时间的区域共享自行车有效车辆数获取方法，包括以下步骤：步骤1，根据实际情况划分共享自行车分析的各个区域，区域数量为m，研究周期为t；步骤2，将订单数据处理转换为车辆的定位数据，并判断每辆车每次定位所在的区域，车辆总数为n，第i辆车的定位次数为di；步骤3，通过判断每辆车的每次定位与下次定位是否在同一区域，计算每辆车在各区域的累计服务时间，来确定各区域的有效车辆数量。作为本发明的进一步技术方案，步骤1中根据土地利用特征、社会经济属性、人口分布特征、自然边界、行政区划、城市道路网络和轨道交通网络等影响因素划分共享自行车分析的各个区域。作为本发明的进一步技术方案，步骤2中订单数据处理转换为车辆的定位数据。定位数据包括车辆id、定位时间、经度、纬度等4个字段，不需要处理转换，可以直接使用。订单数据包括订单id、用户id、车辆id、借车定位时间、借车经度、借车纬度、还车定位时间、还车经度、还车纬度等9个字段，一条订单数据可以拆分为借车、还车的两次定位数据，两次定位的车辆id不变，定位时间、经度、纬度分别是借车、还车的定位时间、经度、纬度。作为本发明的进一步技术方案，步骤3具体为：3.1，计算每辆车在各区域的累计服务时间。其中，tij为第i辆车在第j个区域的累计服务时间，是计算第i辆车多大概率上在第j个区域提供服务的依据。满足约束条件：0≤tij≤t，tid为第i辆车第d次定位到第d+1次定位的时段长度。当d＝0时，tid为第i辆车研究周期开始(认为是第0次定位)到第1次定位的时段长度；当d＝di时，tid为第i辆车第di次定位到研究周期结束(认为是第di+1次定位)的时段长度；当0<d<di时，tid为第i辆车第d次定位到第d+1次定位的时段长度。kdj为在第i辆车第d次定位到第d+1次定位的时段中是否位于第j个区域的系数，取值为0、0.5或1。当d＝0时，如果第i辆车第1次定位是位于第j个区域，kdj＝1，否则kdj＝0；当d＝di时，如果第i辆车第di次定位是位于第j个区域，kdj＝1，否则kdj＝0；当0<d<di时，第i辆车第d次定位和第d+1次定位都是位于第j个区域则kdj＝1，第d次定位和第d+1次定位都不是位于第j个区域则kdj＝0，否则kdj＝0.5。3.2，计算各区域的有效车辆数量。对于第j个区域，基于所有车辆服务时间的有效车辆数量(为第j个区域提供服务的共享自行车车辆数量)为：其中，bj是第j个区域的有效车辆数量。满足约束条件：即所有区域的有效车辆数量之和等于车辆总数n。本发明采用以上技术方案具有以下技术效果：本发明划分共享自行车分析区域，利用企业的运营数据提取出车辆的时空定位数据，根据车辆在区域内的服务时间，确定各区域的有效车辆数量，为共享自行车精细化运营调度奠定指标分析基础。本发明能科学地计算共享自行车在相应区域内的有效车辆数量，克服车辆流动性对各区域运营指标计算的不利影响，为通过运营指标空间研究引导企业高效经营提供了分析基础。附图说明图1是本发明基于服务时间分析获取区域内共享自行车有效车辆数量的流程图。图2是本发明实例的共享自行车分析的区域划分示意图。具体实施方式下面结合附图对本发明的技术方案做进一步的详细说明：图2所示为本发明确定区域内共享自行车有效车辆数量流程图，下面结合图2对本发明作进一步的说明。本发明基于服务时间的区域共享自行车有效车辆数获取方法，主要包含以下步骤：(1)，根据实际情况划分共享自行车分析的各个区域，区域数量为m，研究周期为t。根据土地利用特征、社会经济属性、人口分布特征、自然边界、行政区划、城市道路网络和轨道交通网络等影响因素划分共享自行车分析的各个区域。(2)，将运营数据处理转换为车辆的定位数据，并判断每辆车每次定位所在的区域，车辆总数为n，第i辆车的定位次数为di。定位数据包括车辆id、定位时间、经度、纬度等4个字段。定位数据或订单数据等运营数据处理转换为车辆的定位数据。定位数据包括车辆id、定位时间、经度、纬度等4个字段，不需要处理转换，可以直接使用。订单数据包括订单id、用户id、车辆id、借车定位时间、借车经度、借车纬度、还车定位时间、还车经度、还车纬度等9个字段，一条订单数据可以拆分为借、还车的两次定位数据，两次定位的车辆id不变，定位时间、经度、纬度分别是借、还车的定位时间、经度、纬度。(3)，通过判断每辆车的每次定位与下次定位是否在同一区域，计算每辆车在各区域的累计服务时间，来确定各区域的有效车辆数量。所述步骤(3)中根据服务时间确定各区域的有效车辆数量的具体步骤包括：(31)计算每辆车在各区域的累计服务时间。其中，tij为第i辆车在第j个区域的累计服务时间，是计算第i辆车多大概率上在第j个区域提供服务的依据。满足约束条件：0≤tij≤t。其中，tid为第i辆车第d次定位到第d+1次定位的时段长度。当d＝0时，tid为第i辆车研究周期开始(认为是第0次定位)到第1次定位的时段长度；当d＝di时，tid为第i辆车第di次定位到研究周期结束(认为是第di+1次定位)的时段长度；当0<d<di时，tid为第i辆车第d次定位到第d+1次定位的时段长度。kdj为在第i辆车第d次定位到第d+1次定位的时段中是否位于第j个区域的系数，取值为0、0.5或1。当d＝0时，如果第i辆车第1次定位是位于第j个区域，kdj＝1，否则kdj＝0；当d＝di时，如果第i辆车第di次定位是位于第j个区域，kdj＝1，否则kdj＝0；当0<d<di时，第i辆车第d次定位和第d+1次定位都是位于第j个区域则kdj＝1，第d次定位和第d+1次定位都不是位于第j个区域则kdj＝0，否则kdj＝0.5。(32)计算各区域的有效车辆数量。对于第j个区域，基于所有车辆服务时间的有效车辆数量(真正在第j个区域内服务的共享自行车车辆数量)为：其中，bj是第j个区域的有效车辆数量。满足约束条件：即所有区域的有效车辆数量之和等于车辆总数n。下面用具体实施例来进一步说明本发明的技术方案。案例研究对象为共享单车，共有两辆车，获取的订单数据如表1所示。表1案例订单数据1、根据南京新街口-大行宫附近地区实际情况划分共享自行车分析的区域，分为两个区域，区域id分别为1(代表新街口区域)和2(代表大行宫区域)，研究周期为一天(2017年8月30日)。根据土地利用特征、社会经济属性、人口分布特征、自然边界、行政区划、城市道路网络和轨道交通网络等影响因素划分共享自行车分析的两个个区域,即新街口区域和大行宫区域。2、将运营数据处理转换为车辆的定位数据，并判断每辆车每次定位所在的区域，车辆总数为n，第i辆车的定位次数为di。定位数据包括车辆id、定位时间、经度、纬度等4个字段。定位数据或订单数据等运营数据处理转换为车辆的定位数据。定位数据包括车辆id、定位时间、经度、纬度等4个字段，不需要处理转换，可以直接使用。订单数据包括订单id、用户id、车辆id、借车定位时间、借车经度、借车纬度、还车定位时间、还车经度、还车纬度等9个字段，一条订单数据可以拆分为借、还车的两次定位数据，两次定位的车辆id不变，定位时间、经度、纬度分别是借、还车的定位时间、经度、纬度。表2处理转换后的定位数据车辆id定位时间经度纬度12017/8/306:00:00118.780564°32.041535°12017/8/3018:00:00118.775388°32.046271°22017/8/306:00:00118.780445°32.044703°22017/8/3018:00:00118.799943°32.041301°(3)，通过判断每辆车的每次定位与下次定位是否在同一区域，计算每辆车在各区域的累计服务时间，来确定各区域的有效车辆数量。表2处理转换后的定位数据车辆id定位时间所在区域id定位次序与下次定位是否在同一区域12017/8/306:00:0011是12017/8/3018:00:0012—(最后一次定位)22017/8/306:00:0011否22017/8/3018:00:0022—(最后一次定位)所述步骤(3)中根据服务时间确定各区域的有效车辆数量的具体步骤包括：(31)计算每辆车在各区域的累计服务时间。其中，tij为第i辆车在第j个区域的累计服务时间，是计算第i辆车多大概率上在第j个区域提供服务的依据。满足约束条件：0≤tij≤t。其中，tid为第i辆车第d次定位到第d+1次定位的时段长度。当d＝0时，tid为第i辆车研究周期开始(认为是第0次定位)到第1次定位的时段长度；当d＝di时，tid为第i辆车第di次定位到研究周期结束(认为是第di+1次定位)的时段长度；当0<d<di时，tid为第i辆车第d次定位到第d+1次定位的时段长度。kdj为在第i辆车第d次定位到第d+1次定位的时段中是否位于第j个区域的系数，取值为0、0.5或1。当d＝0时，如果第i辆车第1次定位是位于第j个区域，kdj＝1，否则kdj＝0；当d＝di时，如果第i辆车第di次定位是位于第j个区域，kdj＝1，否则kdj＝0；当0<d<di时，第i辆车第d次定位和第d+1次定位都是位于第j个区域则kdj＝1，第d次定位和第d+1次定位都不是位于第j个区域则kdj＝0，否则kdj＝0.5。表3每辆车在各区域的服务时间计算过程及结果(32)计算各区域的有效车辆数量。对于第j个区域，基于所有车辆服务时间的有效车辆数量(为第j个区域服务的共享自行车车辆数量)为：其中，bj是第j个区域的有效车辆数量。满足约束条件：即所有区域的有效车辆数量之和等于车辆总数n。表4各区域的有效车辆数量计算过程及结果以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属
技术领域：
的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替代，都应当视为属于本发明的保护范围。当前第1页12