[发明专利]基于ArcGIS二次开发的线密度统计方法及系统

权利要求书

1.一种基于ArcGIS二次开发的线密度统计方法，其特征在于，所述方法应用在ArcGIS空间分析工具中，所述方法包括：

步骤101，获取简单线数据集，并根据最小搜索半径的数值、半径增加量的数值和循环次数的数值，按照预设半径公式计算第m个搜索半径的数值，所述简单线数据集包括N条具有起点坐标和终点坐标的简单线，所述N为正整数，所述m的取值依次为从1到M，所述M为正整数，所述半径公式为：

r_m＝r₁+(i-1)Δr

其中，r_m表示所述第m个搜索半径，r₁表示所述最小搜索半径，也表示第1个搜索半径，i表示所述循环次数，Δr表示所述半径增加量；

步骤102，统计起点坐标和终点坐标均在以第n简单线的中点坐标为圆心，以第m个搜索半径的数值为半径构成的空间区域内的简单线的数目，得到第m组第n个第一统计结果，所述n的取值依次为从1到所述N；

步骤1031，获取已配置的显著性水平值，并测试第m组中的N个所述第一统计结果形成的数据分布的分布状况；

步骤1032，若第m组中的N个所述第一统计结果形成的数据分布满足正态分布，则基于所述显著性水平值按照正态分布公式计算得到第m个最小聚类结果；

步骤1033，若第m组中的N个所述第一统计结果形成的数据分布满足帕累托分布，则基于所述显著性水平值按照帕累托公式计算得到第m个最小聚类结果；

步骤1041，统计起点坐标和终点坐标均在以未标记的第a条简单线的中点坐标为圆心，以第b个搜索半径的数值为半径构成的空间区域内的未标记的简单线的数目，得到第b组第a个第二统计结果，所述a为正整数且取值依次为从1至A，所述A为所述简单线数据集中未标记的简单线的数值，其中，根据统计学原理和第m组的N个第一统计结果计算得到第m个最大聚类结果，从第1至第M个最大聚类结果中，查找第一个满足第c个最大聚类结果小于第c个最小聚类结果条件的第c个搜索半径，所述b的初始值为c-1，所述c取值1至所述M中的一个；

步骤1042，从第b组A个所述第二统计结果中提取出数值最大的目标统计结果，判断所述目标统计结果是否大于第b个最小聚类结果；

步骤1043，若所述目标统计结果大于所述第b个最小聚类结果，则记i＝i+1，查询与所述目标统计结果对应的目标简单线，标记起点坐标和终点坐标均在以所述目标简单线的中点坐标为圆心，以第b个搜索半径的数值为半径构成的空间区域内的未标记的简单线和所述目标简单线，并记所述目标统计结果为第i个密度结果，所述i的初始值为0；

步骤1044，判断未标记的简单线的数值是否为0，若未标记的简单线的数值不为0，则返回执行所述步骤1041，若未标记的简单线的数值为0，则基于i个密度结果得到密度结果集合；

步骤1045，若所述目标统计结果小于或等于所述第b个最小聚类结果，则当b大于1时，令b＝b-1，返回执行所述步骤1041，当b小于或者等于1时，基于i个密度结果得到密度结果集合。

2.根据权利要求1所述的方法，其特征在于，

所述正态分布公式为：

minlines＝average(Nls)+r*SD(Nls)

其中，minlines表示所述第m个最小聚类结果，average表示平均值函数，Nls表示第m组中的N个所述第一统计结果，r为与所述显著性水平值相关的参数，当所述显著性水平值为99％时，r取值为2.58，当所述显著性水平值为95％时，r取值为1.96，SD表示标准差函数；

所述帕累托公式为：

其中，p为与所述显著性水平值相关的参数，当所述显著性水平值为99％时，p取值为小于0.01，当所述显著性水平值为95％时，p取值为小于0.05，x_m表示第所述m个最小聚类结果，x表示第m组中的N个所述第一统计结果，α表示回归系数，是一个正参数。