[發明專利]基于ArcGIS二次開發的線密度統計方法及系統

有效

• 申請號: 201810803200.9
• 申請日: 2018-07-20
• 公開（公告）號: CN110019633B
• 公開（公告）日: 2021-08-10
• 發明（設計）人: 辜智慧;張艷
• 申請（專利權）人: 深圳大學
• 主分類號: G06F16/29
• 分類號: G06F16/29;G06F17/18;G06K9/62
• 代理公司: 深圳市恒申知識產權事務所(普通合伙) 44312
• 代理人: 彭海民

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個所述第一統計結果，α表示回歸系數，是一個正參數。