POI是“POInt of Interest”的縮寫,中文可以翻譯為“興趣點”。在地理信息系統(tǒng)中,一個POI可以是一棟房子、一個商鋪、一個郵筒、一個公交站等。根據(jù)POI點識別功能區(qū)的方法有多種,本文采取最常見的方法,即POI點緩沖區(qū)面積占比大小來識別(具體方法應該根據(jù)研究區(qū)和POI點特征來選?。?/p>
研究區(qū)位置:
數(shù)據(jù)準備:
POI點:由于是長沙市中心,服務業(yè)發(fā)達,興趣點主要以商業(yè)點和居住點為主;
1.各類興趣點點比例差異大;
2.在居住區(qū)內(nèi)商業(yè)服務多,會存在重疊現(xiàn)象;
3.分布不均勻。
經(jīng)去除重復點和按興趣點關鍵字重新調(diào)整后
總POI點:181243個
占比:
商業(yè)點:91725
50.6%
居住點:43139
23.8%
公共服務點:26926
14.8%
工業(yè)點:10649
5.9%
交通設施點:9804
5.4%
該地區(qū)的OSM路網(wǎng)數(shù)據(jù):
分析過程:
圖片
對存在拓撲錯誤的路網(wǎng)進行相關處理后把不符合分區(qū)要求(不能有偽節(jié)點,不能有懸掛點)的線段進行裁剪刪除或者添加(拓撲錯誤修改后面我會專門出一篇文章哦),最后把交通線挖空整理形成936個研究單元(線轉(zhuǎn)面工具)。
存在問題:在對交通線進行裁剪,添加形成新的研究單元時可能會對原先的聚集的興趣點產(chǎn)生一定的分割,影響后續(xù)研究區(qū)的功能分區(qū)分析。
措施:對照原始路網(wǎng)和興趣點分布對一些特征明顯的地區(qū)進行重新劃分研究區(qū)。
路網(wǎng)劃分研究單元如下所示:
圖片
首先對添加進來的5類POI點先進行合并,用研究區(qū)裁剪合并后的POI點(研究區(qū)里面有些道路是被挖空了的,落在空白處的POI點不參與后續(xù)計算,所以要全部落在研究區(qū)范圍內(nèi)的POI點才有效哦)然后用研究區(qū)進行標識:
(注意:標識時連接屬性只要選擇ONLY_FID就行了,因為我標識的目的只要知道我這個POI點屬于哪一個研究區(qū)范圍就行了,其他的研究單元屬性字段我不需要,所以不要連接,如果連接了也沒事,就是看得很多字段很亂哈哈哈哈)
標識的作用就是讓所有落入同一個研究單元內(nèi)的點都知道自己屬于這個研究單元,每個POI點后面都有了所在研究單元的所有屬性字段。
圖片
那么問題來了,怎么知道每個研究單元里面每一類POI點數(shù)量有多少呢?(這一步非常重要)首先我們要在研究區(qū)屬性表添加5個字段用來存放不同的類型的POI數(shù)量,然后打開之前標識點數(shù)據(jù)的屬性表,點擊按屬性選擇,點擊區(qū)分POI點類型的字段“type”,里面只有1,2,3,4,5種選擇
選中其中一種,接下來很關鍵哦?。。?!選中后點擊標識后添加過來的研究區(qū)編號字段(FID_研究區(qū))這是唯一區(qū)分研究單元的編號,右鍵選擇匯總:
這個界面啥也不用管,直接點確定就行了,后面會生成一個表,里面紀錄了每一個研究單元里面含有多少個POI點。這樣分別選中5類POI點,匯總5次就可以得到5個表,每個表對應研究單元里面含有一類POI點的數(shù)量,然后分別連接到研究區(qū)數(shù)據(jù),再分別復制POI數(shù)量字段到之前新添加的字段就行了。
圖片
根據(jù)興趣點公共認知度分別做不同半徑的緩沖區(qū),其中商業(yè)點:200米,居住點:250米,公共服務點:300米,工業(yè)點:350米,交通設施點:400米。
上面文章只是將數(shù)據(jù)給處理好,還沒有開始計算,接下來就是計算為主了,要大量使用字段計算器,如果有小伙伴不懂VB或者Python語法的話可以自行區(qū)百度哦。
功能區(qū)劃分標準:
(其中k:第k個研究區(qū),i:第i類興趣點,Sk:第k個研究區(qū)面積,Si:第i類興趣點緩沖區(qū)在該研究區(qū)內(nèi)的面積,ni:第i類興趣點在該研究區(qū)內(nèi)的數(shù)量)
以單一研究區(qū)內(nèi)的某類興趣點緩沖區(qū)面積占該研究單元內(nèi)所有緩沖區(qū)面積的比重作為該興趣點在該研究區(qū)內(nèi)的權重。由于商業(yè)點比重超過50%,其在單一研究區(qū)內(nèi)H值超過60%則定義該研究區(qū)為商業(yè)區(qū),于此相反,交通設施點占比單一研究區(qū)內(nèi)交通設施點H值超過40%作為交通設施功能區(qū),其他均以超過50%作為劃分單一功能區(qū)的標準,余下的劃分為混合功能區(qū)。
混合功能區(qū)細分:分別計算5大類興趣點在某一單一研究區(qū)內(nèi)任意兩類興趣點的H值之和,選擇其中占比最大的H值,如果其比重高于50%則研究區(qū)劃分為##-##功能區(qū),否則定義該功能區(qū)類型為混合功能區(qū)(這是我針對數(shù)據(jù)的獨特性進行設置的識別閾值,各位小伙伴視不同情況使用不同的劃分標準哦,一定要實事求是!)。
用POI點根據(jù)不同類別建立不同半徑的緩沖區(qū)后,分別用研究區(qū)裁剪這5類緩沖區(qū)(因為有些點在研究區(qū)邊緣,緩沖區(qū)會在研究區(qū)范圍外面,這個沒有意義,還有就是研究區(qū)里面有些道路是被挖空了的,落在空白處的緩沖區(qū)面積也不參與后續(xù)計算),然后再用研究區(qū)分別標識這5個緩沖區(qū)。
(注意:標識時連接屬性只要選擇ONLY_FID就行了,因為我標識的目的只要知道我這個緩沖區(qū)屬于哪一個研究區(qū)范圍就行了,其他的研究單元屬性字段我不需要,所以不要連接,如果連接了也沒事,就是看到很多字段很亂哈哈哈哈)
標識結(jié)果如下所示:
(“FID_研究區(qū)”字段是每個研究單元的標識碼,一般會自動添加面積字段,單位為平方米,如果不確定可以自己添加字段-計算幾何-面積再計算一遍)
圖片
那么問題來了,怎么知道每一個研究單元里面每一類POI點緩沖區(qū)面積有多少呢?(這一步非常重要)首先我們要在研究區(qū)屬性表添加5個字段用來存放不同類型的POI緩沖區(qū)面積,接著打開之前標識緩沖區(qū)數(shù)據(jù)的屬性表,點擊自動添加過來的研究區(qū)編號字段(FID_研究區(qū),這是唯一區(qū)分研究單元的編號),右鍵選擇匯總:
一定要記得在面積字段下拉選項下勾選總和,點擊確定后會生成一個表,里面紀錄了每一個研究單元里面含有該類POI點緩沖區(qū)的面積。這樣用相同的方法分別對其他4類標識緩沖區(qū)進行匯總,最后就可以得到5個表,每個表對應一類POI點緩沖區(qū)在所屬研究單元里面的面積,接著根據(jù)字段(FID_研究區(qū))將5個表分別連接到研究區(qū)屬性表,再分別使用字段計算器復制POI點緩沖區(qū)面積字段到之前新添加的字段就行了。
現(xiàn)在研究區(qū)屬性表里面已經(jīng)有了各類POI點數(shù)量以及各類緩沖區(qū)面積,數(shù)據(jù)已經(jīng)準備齊全了,就差計算了。根據(jù)前面的功能區(qū)劃分標準,在研究區(qū)屬性表添加字段“sum”計算5類POI緩沖區(qū)面積*POI點數(shù)量之和:
圖片
具體代碼如下:
sum([ !sy面積!!商業(yè)! + !jz面積!!居住!+ !gy工業(yè)!* !工業(yè)!+ !gf面積!* !公共服務!+ !jt面積!* !交通設施])
注意:*號前面是緩沖區(qū)面積,后面是對應的POI點數(shù)量!?。。。。。?/p>
下面進行單一功能區(qū)的劃分:添加一個文本類型的字段“功能區(qū)”,右鍵字段計算器,選擇python語言,相關參數(shù)如下:
圖片
具體代碼如下:
def N(a,b,c,d,e,A,B,C,D,E,S):
if(S>0.0):
if(1.0*a*A/S>0.6): return "商業(yè)"if(1.0*b*B/S>0.5): return "居住"if(1.0*c*C/S>0.5): return "工業(yè)"if(1.0*d*D/S>0.5): return "公共服務"if(1.0*e*E/S>0.4): return "交通設備"else: return "混合功能區(qū)"
else:
return "數(shù)據(jù)不足"
N (!商業(yè)!,!居住!,!工業(yè)!,!公共服務!,!交通設施!,!sy面積!,!jz面積!, !gy工業(yè)!,!gf面積!,!jt面積!,!sum!)
接著是混合功能區(qū)的細分:添加一個字段“max”,字段計算器輸入如下代碼:
max([!sy面積!!商業(yè)!+!jz面積!!居住!,!sy面積!!商業(yè)!+ !gy工業(yè)!!工業(yè)!,!sy面積!* !商業(yè)!+!gf面積!!公共服務!,!sy面積!!商業(yè)!+!jt面積!!交通設施!,!jz面積!!居住!+!gy工業(yè)!!工業(yè)!,!jz面積!!居住!+!gf面積!!公共服務!,!jz面積!!居住!+!jt面積!!交通設施!,!gy工業(yè)!!工業(yè)!+!gf面積!!公共服務!,!gy工業(yè)!!工業(yè)!+!jt面積!!交通設施!,!gf面積!!公共服務!+!jt面積!*!交通設施!])
max函數(shù)是計算里面的最大值,這樣可以得到單一研究區(qū)內(nèi)任意兩類興趣點的H值之和的最大值,方便后續(xù)計算。
按屬性選擇“功能區(qū)”字段值為混合功能區(qū)的行,字段計算器設置如下:
圖片
具體代碼如下所示:
def N(a,b,c,d,e,f,g,h,i,j,m,s):
if(1.0*m/s<0.5):
return "混合功能區(qū)"
else:
if(m==j): return "公共服務-交通設備"if(m==i): return "工業(yè)-交通設備"if(m==h): return "工業(yè)-公共服務"if(m==g): return "居住-交通設備"if(m==f): return "居住-公共服務"if(m==e): return "居住-工業(yè)"if(m==d): return "商業(yè)-交通設備"if(m==c): return "商業(yè)-公共服務"if(m==b): return "商業(yè)-工業(yè)"else: return "商業(yè)-居住"
N(!sy面積!!商業(yè)! + !jz面積!!居住!, !sy面積!!商業(yè)! + !gy工業(yè)!!工業(yè)!, !sy面積!* !商業(yè)! + !gf面積!* !公共服務! , !sy面積!* !商業(yè)! + !jt面積!* !交通設施!, !jz面積!!居住! + !gy工業(yè)!!工業(yè)!, !jz面積!* !居住! + !gf面積!!公共服務! , !jz面積!!居住! + !jt面積!!交通設施! , !gy工業(yè)!!工業(yè)! + !gf面積! *!公共服務! , !gy工業(yè)!!工業(yè)! + !jt面積!!交通設施!, !gf面積! *!公共服務! + !jt面積! *!交通設施!,!max!,!sum!)
最后
至此功能區(qū)劃分完畢,最后右鍵研究區(qū)右鍵選擇屬性,設置如下:
符號化后功能區(qū)識別結(jié)果圖:
來源:https://blog.csdn.net/qq_43173805/article/details/119996031
成都途遠GIS是一家專業(yè)致力于無人機航空攝影測繪、航空數(shù)據(jù)處理、GIS地理信息系統(tǒng)研發(fā)、數(shù)字孿生城市制作、數(shù)字沙盤模型等業(yè)務的創(chuàng)新型科技公司,為您提供一站式地理信息服務。
本文鏈接:http://www.lajichun.cn/blog/500.html
本文標簽:ArcGIS