1.一种基于时空约束的城市车辆可达区域计算方法，其特征在于，包括以下步骤：步骤S1，根据轨迹数据集筛选得到轨迹序列点，以根据所述轨迹序列点计算轨迹可达点；

步骤S2，将存在所述轨迹可达点的轨迹覆盖道路集中，选择道路轨迹抵达数最大的道路，加入边界点计算道路集；以及步骤S3，计算所述边界点计算道路集中所有的边界点和查询时空点连线与x轴正方向的角度，依据所述角度排序并连接所述边界点，获得可达区域。

2.根据权利要求1所述的基于时空约束的城市车辆可达区域计算方法，其特征在于，所述步骤S1包括：步骤S101，在轨迹数据集中筛选出符合查询时空点时空约束的轨迹数据集，以构成所述目标轨迹集；

步骤S102，在所述目标轨迹集中筛选出符合查询时间段约束的轨迹序列点，以构成目标轨迹片段；

步骤S103，记所述目标轨迹片段在道路上的序列点集，以生成轨迹可达点。

3.根据权利要求1所述的基于时空约束的城市车辆可达区域计算方法，其特征在于，所述步骤S2包括：步骤S201，当所述轨迹覆盖道路集中的任意轨迹覆盖道路存在所述轨迹可达点时，记将所述轨迹可达点加入轨迹可达点集中；

步骤S202，遍历所述轨迹覆盖道路集中任意所述轨迹覆盖道路，计算起点和所述查询时空点连线与x轴正方向的第一角度，终点和所述查询时空点S连线与x轴正方向的第二角度，其中，所有轨迹覆盖道路的道路角度构成道路角度集，并规定所述第一角度大于所述第二角度；以及步骤S203，在所述轨迹覆盖道路集中选择道路轨迹抵达数最大的道路，加入边界点计算道路集。

4.根据权利要求3所述的基于时空约束的城市车辆可达区域计算方法，其特征在于，所述步骤S2还包括：步骤S204，遍历所述道路角度集，对除所述轨迹可达点集之外的任意道路角度进行所述轨迹覆盖道路筛选，并更新所述轨迹覆盖道路集属性值；

步骤S205，将所述轨迹覆盖道路集中的所述轨迹覆盖道路删除，并将所述道路角度集中所述道路轨迹抵达数最大的道路的道路角度删除，判断所述轨迹覆盖道路集是否为空。

5.根据权利要求1所述的基于时空约束的城市车辆可达区域计算方法，其特征在于，所述步骤S3包括：步骤S301，遍历所述边界点计算道路集中的所述轨迹可达点，并根据轨迹抵达数对所述轨迹可达点进行降序排列，选择序列中前K个轨迹可达点作为候选轨迹可达点并加入候选可达点集；

步骤S302，依据所述候选可达点集计算边界点，并将所述边界点计算道路集上的所有边界点构成边界点集；以及步骤S303，计算所述边界点集中所有的边界点和所述查询时空点连线与x轴正方向的角度，并将所有角度升序排列，依次连接所述边界点，获得所述可达区域。

6.根据权利要求2所述的基于时空约束的城市车辆可达区域计算方法，其特征在于，所述步骤S103中的所述轨迹可达点Vi(vxi,vyi,numi)计算公式为：其中，vxi为所述轨迹可达点横坐标，STik为所述查询时空点出发的轨迹数据上且在查询时间段T内的轨迹片段， 为轨迹中第m个序列点， 为序列点地理经度转换到平面坐标系的横坐标，vyi为所述轨迹可达点纵坐标， 为序列点地理纬度转换到平面坐标系的纵坐标，numi为抵达该位置点的轨迹数。

7.根据权利要求6所述的基于时空约束的城市车辆可达区域计算方法，其特征在于，当|STia|<|STib|时，则选取起点上的所述目标轨迹片段计算，生成所述轨迹可达点Vi(vxi,i i a bvy ,num )，其中， 当|STi |>|STi |时，则选择终点上的所述轨迹片段计算，生成所述轨迹可达点Vi(vxi,vyi,numi)，其中，

8.根据权利要求4所述的基于时空约束的城市车辆可达区域计算方法，其特征在于，所述步骤S205中判断过程为：将所述轨迹覆盖道路从所述轨迹覆盖道路集中删除，所述道路角度从所述道路角度集中删除，判断所述轨迹覆盖道路集是否为空集，若为空集，则结束，否则执行所述步骤203。

9.根据权利要求5所述的基于时空约束的城市车辆可达区域计算方法，其特征在于，所述依据所述候选可达点集计算边界点的计算公式为：j j

其中，qx为所述边界点的横坐标，qy为所述边界点的纵坐标。

10.根据权利要求5所述的基于时空约束的城市车辆可达区域计算方法，其特征在于，所述计算所述边界点集中所有的边界点和所述查询时空点连线与x轴正方向的角度的计算公式为：其中，Ai为角度值，sy为所述查询时空点的纵坐标，sx为所述查询时空点的横坐标。