[发明专利]基于RGB色彩空间与蒙特卡洛方法的城市公交指标计算方法

权利要求书

1.基于RGB色彩空间与蒙特卡洛方法的城市公交指标计算方法，其特征在于，包括如下步骤：

步骤1，获取城市待研究区域的总体规划，所述待研究区域的总体规划包括待研究区域的土地利用总体规划图和待研究区域的土地面积；

步骤2，获取城市待研究区域中的公共交通站点位置信息和公共交通站点数量，并使用WGS84国际通用地理坐标系的经纬度来表示公共交通站点位置信息；

步骤3，将土地利用总体规划图划分为m×n个像素点，记录各个像素点的RGB向量，得到图二维矩阵

步骤4，根据土地利用总体规划图确定新的坐标系，根据新的坐标系将公共交通站点位置信息转化为公共交通站点位置信息图坐标；具体过程如下：

4.1，选取土地利用总体规划图中的标志性地点作为坐标原点，其WGS84坐标为(i₀,j₀)，土地利用总体规划图的正北方向作为坐标正方向，建立新的坐标系；

4.2，在土地利用总体规划图上选取两个公共交通站点，根据这两个公共交通站点位置信息以及这两个公共交通站点在土地利用总体规划图上的直线距离计算第一转化因子α；

所述第一转化因子α的计算公式为：

第一转化因子α等于两个公共交通站点在土地利用总体规划图上的直线距离除以两个公共交通站点的实际距离再乘以比例尺的分母；

4.3，根据步骤4.2选取的两个公共交通站点在土地利用总体规划图上的直线距离与这两个公共交通站点之间的像素点个数计算第二转化因子β；

所述第二转化因子β的计算公式为：

第二转化因子β等于两个公共交通站点之间的像素点个数除以两个公共交通站点在土地利用总体规划图上的直线距离；

4.4，根据第一转化因子α和第二转化因子β，将公共交通站点位置信息转化为公共交通站点位置信息图坐标

所述公共交通站点位置信息图坐标具体为：

步骤5，确定土地利用总体规划图的空白区域与非城市用地区域的RGB向量，并用一个二维列表记录下来；

所述二维列表具体为：

其中，k＝1,…,t，t为土地利用总体规划图上非城市用地区域及空白区域的RGB个数，b_k为RGB向量，为一个三维向量(R±ε,G±ε,B±ε)，ε为色彩空间误差；

步骤6，基于蒙特卡洛方法，计算待研究区域的用地总面积和公交覆盖总面积；具体过程如下：

6.1，初始化有效研究区域点数M_all＝0和有效公交覆盖点数M_bus＝0；

6.2，确定迭代次数N，N取值为100×m×n，m、n分别代表图二维矩阵的行数、列数；

6.3，确定色彩空间误差ε，ε取值为5；

6.4，每次迭代随机确定位置坐标(x,y)，其中x∈[0,m]，y∈[0,n]，记录a_xy，a_xy代表位置坐标(x,y)的RGB向量；

6.5，若则M_all＝M_all+1；

6.6，若且则M_bus＝M_bus+1，代表公共交通站点位置信息图坐标，K代表转化后的服务半径；

所述转化后的服务半径K＝K₀×α×β，K₀代表公共交通站点服务半径，α为第一转化因子，β为第二转化因子；

6.7，迭代终止后，得到待研究区域的用地总面积公交覆盖总面积A_True代表待研究区域的土地面积；

步骤7，根据公交覆盖总面积和用地总面积计算公交覆盖率，根据公共交通站点数量和用地总面积计算公交站点密度。

2.根据权利要求1所述基于RGB色彩空间与蒙特卡洛方法的城市公交指标计算方法，其特征在于，步骤3所述图二维矩阵具体为：

其中，元素a_ij代表第i行第j列的像素点的RGB向量，为一个三维向量(R,G,B)，i＝1,…,m，j＝1,…,n。

3.根据权利要求1所述基于RGB色彩空间与蒙特卡洛方法的城市公交指标计算方法，其特征在于，步骤7所述公交覆盖率等于公交覆盖总面积除以用地总面积，公交站点密度等于公共交通站点数量除以用地总面积。