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