基于核密度估計(jì)的點(diǎn)群密度制圖應(yīng)用研究

盧 敏,楊 柳,王金茵,黃 煌,王結(jié)臣

(1.江蘇省地理信息技術(shù)重點(diǎn)實(shí)驗(yàn)室，江蘇 南京 210023;2.南京大學(xué) 地理信息科學(xué)系，江蘇 南京 210023;3.江蘇省基礎(chǔ)地理信息中心，江蘇 南京 210013)

基于核密度估計(jì)的點(diǎn)群密度制圖應(yīng)用研究

盧 敏1,2,楊 柳3,王金茵1,2,黃 煌1,2,王結(jié)臣1,2

(1.江蘇省地理信息技術(shù)重點(diǎn)實(shí)驗(yàn)室，江蘇 南京 210023;2.南京大學(xué) 地理信息科學(xué)系，江蘇 南京 210023;3.江蘇省基礎(chǔ)地理信息中心，江蘇 南京 210013)

隨著信息技術(shù)的發(fā)展，點(diǎn)狀專題數(shù)據(jù)激增，需要突破傳統(tǒng)制圖方法的表達(dá)局限，利用海量點(diǎn)群數(shù)據(jù)對(duì)專題信息進(jìn)行快速、直觀、綜合的制圖表達(dá)。文中以南京市銀行的ATM網(wǎng)點(diǎn)數(shù)據(jù)以及全國KTV營業(yè)點(diǎn)分布為制圖數(shù)據(jù)源，探討了點(diǎn)群數(shù)據(jù)核密度估計(jì)圖以及三維核密度曲面圖的制作方法，并對(duì)其中的關(guān)鍵技術(shù)作了深入的研究。研究表明，基于核密度估計(jì)的網(wǎng)點(diǎn)密度制圖方法，消除了邊界分隔對(duì)密度計(jì)算精度的影響，能更好反映點(diǎn)群數(shù)據(jù)的空間分布密度、區(qū)域差異、動(dòng)態(tài)變化等特征，具備更豐富的表現(xiàn)力。

核密度估計(jì)；點(diǎn)群密度制圖；最佳帶寬計(jì)算；三維核密度曲面圖

點(diǎn)狀分布是專題要素空間分布的最常見狀態(tài)之一，比如銀行、學(xué)校、超市、停車場、餐廳等[1]，這些點(diǎn)狀分布的專題要素種類繁多、分布廣泛，點(diǎn)群目標(biāo)的制圖表達(dá)適用范圍廣泛，通過合理運(yùn)用豐富的專題地圖直觀、有效的可視化表達(dá)，可以綜合多樣地展現(xiàn)研究成果、直觀地傳達(dá)專題信息，并反映其空間分布的分布密度、區(qū)域差異、動(dòng)態(tài)變化等特征，使讀圖者一目了然。從整體上看，點(diǎn)群要素最重要的特征是空間分布特征，如何直觀、快速、有效地表達(dá)點(diǎn)群目標(biāo)在空間中的數(shù)量、質(zhì)量、密度和相互關(guān)系，是制圖表達(dá)的重點(diǎn)[2]。隨著信息技術(shù)的發(fā)展，數(shù)據(jù)獲取能力的增強(qiáng)，許多領(lǐng)域、行業(yè)的點(diǎn)狀專題數(shù)據(jù)激增，如何突破傳統(tǒng)制圖方法的表達(dá)局限，利用海量點(diǎn)群數(shù)據(jù)對(duì)專題信息進(jìn)行快速、直觀、綜合的制圖表達(dá)，成為了制圖領(lǐng)域亟待解決的問題。

對(duì)點(diǎn)群要素空間分布特征的表達(dá)，應(yīng)力求展現(xiàn)其分布密度的區(qū)域差異，進(jìn)而揭示某些潛在規(guī)律，如居民行為習(xí)慣、事件發(fā)生熱點(diǎn)、商圈吸引力等[3-4]。利用離散點(diǎn)要素進(jìn)行制作密度圖時(shí)，通常有兩種思路。一種是將整個(gè)區(qū)域按指定的邊界(如行政區(qū))進(jìn)行分割，統(tǒng)計(jì)各個(gè)區(qū)域內(nèi)落入的點(diǎn)的數(shù)量，計(jì)算各區(qū)域的要素密度，從而將點(diǎn)要素轉(zhuǎn)換為具有指定邊界和不同值的面要素。在該思路中可采用量底法、點(diǎn)值法等方法制圖，但是在這類依邊界分隔的密度圖中，制圖表達(dá)的精度與區(qū)域劃分的精細(xì)程度密切相關(guān)，區(qū)域劃分越小，精度越高，區(qū)域劃分越大，精度也就越低。因此，另一種思路是不進(jìn)行區(qū)域分割，根據(jù)離散的點(diǎn)要素?cái)?shù)據(jù)內(nèi)插計(jì)算整個(gè)區(qū)域的值，從而生成一個(gè)連續(xù)的密度表面，通過顏色、符號(hào)等在底圖上的疊加來展現(xiàn)專題要素的聚集特點(diǎn)[5]。根據(jù)內(nèi)插原理的不同可分為簡單點(diǎn)密度制圖和核函數(shù)密度制圖，可以采用等值線法、量底法、三角形圖表法等方法制圖。點(diǎn)密度制圖與核密度制圖均是通過將點(diǎn)狀要素轉(zhuǎn)換為連續(xù)曲面，并以面上某點(diǎn)的值代表該點(diǎn)的要素密度，來消除邊界分隔對(duì)密度計(jì)算精度的影響，較好地展現(xiàn)了點(diǎn)群數(shù)據(jù)的聚集特點(diǎn)。

本次研究將以南京市所有銀行的ATM(Automatic Teller Machine，自動(dòng)柜員機(jī))營業(yè)網(wǎng)點(diǎn)數(shù)據(jù)，以及中國大陸地區(qū)KTV(Karaoke Television，通常指提供卡拉OK影音設(shè)備與視唱空間的場所)營業(yè)網(wǎng)點(diǎn)數(shù)據(jù)為例，探討核密度估計(jì)方法在密度制圖中的應(yīng)用。

1 核密度估計(jì)

1.1 核密度估計(jì)的基本原理

事件在空間上的出現(xiàn)具有一定的隨意性，但這種隨意性在一定的空間過程作用下會(huì)受到影響，從而導(dǎo)致其在不同位置上出現(xiàn)的概率或高或低。如果某一區(qū)域出現(xiàn)更多的事件，則可以認(rèn)為事件在受到某種空間作用下在此處出現(xiàn)的概率會(huì)更高，反之則更低。結(jié)合地理學(xué)第一定律空間自相似的規(guī)律特征，這種密度(或者概率)上的變化可以用來區(qū)分事件在空間上的分布模式[6-8]。核密度估計(jì)(Kernel Density Estimation, KDE)方法就反映了這一思想，它認(rèn)為區(qū)域內(nèi)任意位置都有一個(gè)可測(cè)度的事件密度(也稱強(qiáng)度)，該位置的事件密度可以通過其周圍單位面積區(qū)域內(nèi)的事件點(diǎn)數(shù)量來估計(jì)[9-10]。核密度估計(jì)可以想象成一個(gè)移動(dòng)的三維函數(shù)(內(nèi)核)，其影響范圍內(nèi)的事件點(diǎn)到s的距離視為事件點(diǎn)對(duì)s影響的權(quán)重值，進(jìn)而計(jì)算出點(diǎn)s處的事件密度[11]。

(1)

式中：k( )表示的是核的權(quán)重函數(shù)，τ>0，稱作帶寬，即以s為原點(diǎn)的曲面在空間上延展的寬度，τ值的選擇會(huì)影響到分布密度估計(jì)的光滑程度；s-si是密度估值點(diǎn)s和si之間的距離。對(duì)于核密度估計(jì)的結(jié)果，可以采用合適的等值線算法或者以柵格的形式，將研究區(qū)域R以密度變化的連續(xù)表面來表示。

圖1 對(duì)點(diǎn)模式的核密度估計(jì)

(2)

式中：τ為核函數(shù)的帶寬，di為事件點(diǎn)到s的距離。研究表明，不同數(shù)學(xué)形式和核函數(shù)對(duì)密度估計(jì)的影響很小，對(duì)KDE的曲面特征影響較大是帶寬τ的選取。因此,帶寬的選擇比較重要，而核函數(shù)可以根據(jù)實(shí)際需求選擇即可。

1.2 最適帶寬的計(jì)算方法

在KDE方法中，帶寬的確定至關(guān)重要，在很大程度上會(huì)影響到分析的結(jié)果。實(shí)際的經(jīng)驗(yàn)表明，如果想要更光滑的密度表面，τ的值可以設(shè)的大一些，然而此時(shí)就有可能掩蓋想要探測(cè)的熱點(diǎn)區(qū)域；而如果將τ的值設(shè)的太小，所得到的密度表面變化就會(huì)變得突兀不平，即局部的特征會(huì)得到凸顯，而忽略了稍大尺度上的整體關(guān)聯(lián)性，使得熱點(diǎn)區(qū)域過于分散。因此，在KDE方法的實(shí)際應(yīng)用中，需要采用合適的方法計(jì)算得到τ的取值，探索相應(yīng)帶寬值τ下密度表面的光滑程度，以檢驗(yàn)該帶寬下的KDE結(jié)果是否合理。

一些學(xué)者經(jīng)過研究總結(jié)出一些最適值的計(jì)算方法，其中最簡單的一種是根據(jù)事件點(diǎn)數(shù)據(jù)的最小外接矩形確定，即用最小外接矩形的較短邊除以30來確定τ值，即

(3)

這種方法缺點(diǎn)是僅考慮了研究區(qū)的形狀特征，而沒要考慮事件點(diǎn)數(shù)量的影響。

一種改進(jìn)的τ值確定方式是，首先確定整個(gè)事件點(diǎn)的平均中心，然后計(jì)算平均中心到各個(gè)事件點(diǎn)的距離，取這些距離的中位數(shù)Dm，并計(jì)算事件點(diǎn)的標(biāo)準(zhǔn)距離SD(所謂標(biāo)準(zhǔn)距離與標(biāo)準(zhǔn)差類似，是對(duì)事件點(diǎn)的一種概括性描述)，則τ值滿足公式[12]

(4)

式中,n為事件點(diǎn)數(shù)量。

Bailey和Gatrell給出τ的建議值[13]，其采用事件點(diǎn)數(shù)量的0.68倍的-0.2次方與研究區(qū)面積的平方根的乘積確定，即

(5)

式中:n為事件點(diǎn)數(shù)量，A為研究區(qū)面積。

上述幾種方法中均沒有考慮事件點(diǎn)之間的空間關(guān)系，因此，有學(xué)者提出一種基于事件點(diǎn)之間的k階最鄰近距離方法來確定τ值[14-15]，即

(6)

其中,dij表示k階最鄰近距離，即從一個(gè)事件點(diǎn)一直到第k個(gè)最鄰近點(diǎn)的距離均值。k值決定了密度表面的光滑程度，k值越大帶寬τ也就越大，生成的密度表面也就越光滑。

2 核密度估計(jì)的密度制圖

核密度制圖的思路是，給每個(gè)點(diǎn)要素定義一個(gè)圓形鄰域，在其上方覆蓋一個(gè)平滑曲面，依據(jù)二次核函數(shù)運(yùn)算使曲面值呈正態(tài)分布，在該點(diǎn)所在位置處曲面的值最高，隨著與點(diǎn)的距離的增大曲面值逐漸減小，到鄰域邊界處減小為零。其中，曲面最高值所在點(diǎn)到圓形鄰域值為零處的距離，就是帶寬。每個(gè)輸出像元的密度均為疊加在該柵格像元中心的所有密度曲面值之和。

平面視角下的核密度估計(jì)，涉及的主要參數(shù)為帶寬。根據(jù)已有文獻(xiàn)的研究成果[16-17]，總結(jié)出幾種帶寬的計(jì)算方法。但是由于這些帶寬的計(jì)算方法對(duì)于數(shù)據(jù)集的側(cè)重點(diǎn)各不相同，因此分別計(jì)算所得到的帶寬值差異較大。在具體研究中，通常的做法是：首先利用這幾種方法分別計(jì)算來得到可選的帶寬，然后將其作為參數(shù)分別進(jìn)行核函數(shù)估計(jì)，最后根據(jù)曲面的平滑程度，以及是否能兼顧全局趨勢(shì)信息與局部熱點(diǎn)信息，從而最終確立合適的核密度估計(jì)帶寬。

本次研究以南京市主城區(qū)的ATM網(wǎng)點(diǎn)數(shù)據(jù)為例，研究區(qū)包括了南京市的主要城區(qū)范圍(面積約為254.86 km2)，銀行ATM數(shù)據(jù)來源于網(wǎng)絡(luò)電子地圖，研究中主要選擇了南京市主城區(qū)范圍內(nèi)網(wǎng)點(diǎn)數(shù)量大于50的商業(yè)銀行ATM網(wǎng)點(diǎn)，包括工商銀行、建設(shè)銀行、農(nóng)業(yè)銀行等合計(jì)749家ATM網(wǎng)點(diǎn)。其他主要參數(shù)和帶寬的計(jì)算結(jié)果如表1所示。

表1 研究區(qū)帶寬計(jì)算結(jié)果

根據(jù)表1的計(jì)算結(jié)果分別就4種不同的帶寬進(jìn)行密度估計(jì)計(jì)算，其結(jié)果如圖2所示。由于不同帶寬下計(jì)算所得到密度估計(jì)值的方法各不相同，為了更好地比較4種帶寬密度估計(jì)的影響，在分級(jí)渲染的過程中并不能采用統(tǒng)一的分級(jí)斷點(diǎn)，這里統(tǒng)一采用“自然斷點(diǎn)”進(jìn)行分級(jí)。

圖2 不同帶寬下的核密度估計(jì)結(jié)果

由圖2可以看出，隨著帶寬的增大，密度估計(jì)表面越來越趨向于光滑，特別是當(dāng)τ= 4 590 m時(shí)，整個(gè)密度曲面基本上呈現(xiàn)一種由中心向外圍擴(kuò)散的圈層結(jié)構(gòu)，雖然這在一定程度上反應(yīng)了ATM網(wǎng)點(diǎn)在南京主城區(qū)范圍內(nèi)分布的大體趨勢(shì)(由中心城區(qū)向外圍遞減)，但對(duì)網(wǎng)點(diǎn)分布的熱點(diǎn)區(qū)域探測(cè)并沒有太大的用處。對(duì)于其他3種方法而言，采用τ=643 m和τ=1 022 m這兩種帶寬的效果均較為理想，在反映網(wǎng)點(diǎn)分布的整體趨勢(shì)的同時(shí)，也能夠很好地呈現(xiàn)局部的熱點(diǎn)區(qū)域信息；k階最鄰近距離方法這里取k=3時(shí)計(jì)算所得的帶寬，即τ=345 m，可以看出其探測(cè)出的熱點(diǎn)區(qū)域最為分散，對(duì)網(wǎng)點(diǎn)分布熱點(diǎn)區(qū)域的探測(cè)而言尺度過小。但是事實(shí)上，k階最鄰近距離方法具有很強(qiáng)的靈活性，上述其他幾種方法計(jì)算得到的帶寬都能在k階最鄰近距離方法相應(yīng)的階數(shù)下找到對(duì)應(yīng)的值，因此，一種選擇合適帶寬的方式可以是：首先通過默認(rèn)方法計(jì)算其帶寬，然后在k階最鄰近距離方法中找到對(duì)應(yīng)的階數(shù)，并以此為基準(zhǔn)，增加或減少階數(shù)直到找到適合的帶寬，得到較好的視覺效果。

KDE的方法消除了樣方分析中由于樣方的尺寸和形狀等對(duì)局部密度的影響，并且具有更好的可視化效果，在探索事件分布熱點(diǎn)、測(cè)度局部密度概括指標(biāo)上具有一定的優(yōu)勢(shì)。然而KDE雖然能夠?qū)哂幸浑A密度或平穩(wěn)性的分布模式做出有效且實(shí)用的檢驗(yàn)，但對(duì)于二階距離特性的探測(cè)還顯得不夠。

3 三維核密度曲面圖

將普通二維專題地圖進(jìn)行三維顯示，通過對(duì)專題信息的立體化表達(dá)，能給讀圖者帶來更直接更具沖擊力的視覺感受。同時(shí)，三維制圖還有助于減弱專題地圖中對(duì)色彩的依賴，使地圖適應(yīng)傳統(tǒng)印刷工藝對(duì)色彩的限制。在立體的現(xiàn)實(shí)世界中，具有高低起伏特點(diǎn)的景觀有很多，縱橫交錯(cuò)的山脈、鱗次櫛比的樓宇，都可以作為三維制圖的效果參照。受人們閱讀習(xí)慣和設(shè)備制作工藝的限制，目前常用的圖形顯示設(shè)備依然采用平面顯示方式，三維效果通常需要依靠物體形狀、顏色、陰影、輔助線等的輔助渲染來體現(xiàn)[18]。

相對(duì)于平面地圖而言，三維模擬更形象、直觀，也更符合人們?nèi)粘Ｓ^察事物的習(xí)慣。三維表面圖的應(yīng)用較為廣泛，不僅可以直觀地展現(xiàn)高程、氣溫、人口密度等連續(xù)變量，還可以通過疊加其它要素參數(shù)來表達(dá)綜合專題信息。在點(diǎn)群要素三維制圖表達(dá)中，可根據(jù)數(shù)據(jù)特點(diǎn)和地圖用途的不同，選取適當(dāng)?shù)闹茍D表達(dá)方法，結(jié)合構(gòu)圖和色彩搭配，可以形成較好的效果。較高的區(qū)域會(huì)遮擋后面較矮的區(qū)域，若要較好地掌握數(shù)據(jù)的整體狀況，需從多個(gè)角度進(jìn)行觀察。對(duì)地形的三維模擬多從這兩種視角出發(fā)，或是如身臨其境般的虛擬現(xiàn)實(shí)場景，或是用沙盤等立體模型展現(xiàn)地勢(shì)起伏。為了更好地展現(xiàn)制圖區(qū)域內(nèi)的總體情況，地形模擬更多地采用了俯瞰的視角。然而，從地圖制圖的角度來看，在現(xiàn)有的專題地圖制圖表達(dá)方法中，三維地圖制圖方法一般用于自然地理要素的表達(dá)，如地形的高低起伏、溫度的變化、降雨量地區(qū)差異的展示等[19]，但對(duì)于人文地理要素中社會(huì)經(jīng)濟(jì)現(xiàn)象的空間分布特征而言，三維的制圖表達(dá)手段更加直觀、豐富多樣，但卻是少有人運(yùn)用該方法來呈現(xiàn)此類人文地理要素的空間差異特征。劉銳、胡偉平等人[20]利用核密度估計(jì)模型與基尼系數(shù)為切入點(diǎn)，以我國沿海11個(gè)省份2000—2012年的數(shù)據(jù)為例，繪制了海洋經(jīng)濟(jì)發(fā)展的動(dòng)態(tài)演變趨勢(shì)與規(guī)律，但也依舊沒有使用三維制圖進(jìn)行專題制圖展示。

本次制圖所用的KTV營業(yè)網(wǎng)點(diǎn)數(shù)據(jù)是通過網(wǎng)絡(luò)電子地圖整理得到的。作為一種典型的娛樂休閑行業(yè)，KTV行業(yè)在中國的發(fā)展并不均衡。經(jīng)統(tǒng)計(jì)，在占大陸地區(qū)全部國土面積13.5%的東部地區(qū)12省(或自治區(qū)、直轄市)中，共有KTV營業(yè)網(wǎng)點(diǎn)44 637個(gè)，占大陸地區(qū)KTV網(wǎng)點(diǎn)總量的50.9%；在占大陸地區(qū)國土面積29.3%的中部地區(qū)9省(或自治區(qū))中，擁有KTV營業(yè)網(wǎng)點(diǎn)24 883個(gè)，占大陸地區(qū)KTV網(wǎng)點(diǎn)總量的28.4%；西部地區(qū)的10個(gè)省(或自治區(qū)、直轄市)占據(jù)大陸地區(qū)國土面積的57.2%，共有KTV營業(yè)網(wǎng)點(diǎn)18 210個(gè)，僅占大陸地區(qū)KTV網(wǎng)點(diǎn)總量的20.7%。利用前文所描述的核密度最佳帶寬確定方法，首先分別計(jì)算出4種方法對(duì)應(yīng)的帶寬，然后分別用這4個(gè)帶寬作出核密度估計(jì)圖，通過視覺的對(duì)比，選出其中能兼顧全局趨勢(shì)信息與局部熱點(diǎn)信息的密度圖，從而確定制圖的最佳帶寬。通過該最佳帶寬做出平面圖后添加俯瞰視角的三維顯示效果，可以得到如圖3所示的效果圖。

可以看到，將平面的核密度圖拉伸成立體曲面后，密度較高的區(qū)域得到有效拉伸，輔以色彩上的區(qū)分，能夠有效突出了區(qū)域間的差異。經(jīng)過三維立體化并拉伸顯示后的核密度曲面，高低值之間的差異被適當(dāng)?shù)乜浯螅@樣可以有效地凸顯出熱點(diǎn)，如圖3所示的上海、北京、廣州等，均是KTV密度最高的幾個(gè)城市。此外，在三維核密度曲面上疊加了注記、圖例、行政邊界等輔助讀圖的地圖要素之外，由于三維透視中存在顯示比例近大遠(yuǎn)小的變化，圖中無法添加比例尺，為了彌補(bǔ)這一缺憾，在制圖中還添加了規(guī)則的格網(wǎng)，它不僅使讀圖者對(duì)空間位置的把握更準(zhǔn)確，并且有指示方向的作用。因此在制圖的過程中，應(yīng)當(dāng)根據(jù)制圖數(shù)據(jù)的特點(diǎn)和制圖需要選取合適的角度，全面、清晰、直觀地展現(xiàn)空間和屬性信息，提升地圖的可讀性和可用性。

圖3 中國大陸KTV營業(yè)網(wǎng)點(diǎn)三維密度圖

4 結(jié)束語

大多數(shù)專題制圖要素在空間上的分布是不均衡的，使用傳統(tǒng)的專題地圖制圖表達(dá)手法往往受區(qū)域幾何面積、分區(qū)單元大小等的影響，不能直觀地展現(xiàn)局部統(tǒng)計(jì)信息的密度差異。在諸如量底法這類依邊界分隔的密度圖中，制圖表達(dá)的精度與區(qū)域劃分的精細(xì)程度相關(guān)，區(qū)域劃分越小，精度越高，就越能消除分區(qū)間的密度差異，分區(qū)間的密度分布越平滑；反之區(qū)域劃分越大，精度也就越低，分區(qū)間的密度差異越大，密度分布越粗糙。

本文將核密度估計(jì)的相關(guān)理論運(yùn)用到空間點(diǎn)群數(shù)據(jù)的制圖上，分析了不同數(shù)據(jù)源的情況下KDE函數(shù)帶寬的變化對(duì)核密度制圖的影響。研究表明，利用核密度估計(jì)的方法能夠計(jì)算得到較為平滑的密度分布，能有效消除邊界分隔對(duì)密度計(jì)算造成的影響。而三維曲面圖用立體模型展現(xiàn)其凸起和凹陷，從而使專題信息的展示更加直觀多樣。不僅有效地反映了專題要素的空間形態(tài)、數(shù)量特征，而且具有方向性和對(duì)比性，表現(xiàn)力出色。這兩種方法是專題制圖方法的擴(kuò)充與發(fā)展，為點(diǎn)群數(shù)據(jù)的制圖表達(dá)提供了一種新的思路。

[1] 周愛華,張寶秀,董恒年,等. 基于GIS的北京城區(qū)老字號(hào)小吃店空間分布特征研究[J].北京聯(lián)合大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,29(4):29-35.

[2] 曹一冰,付洪貝,潘志超. 地圖設(shè)計(jì)中的美學(xué)研究[J]. 測(cè)繪與空間地理信息,2011,34(1): 231-233.

[3] 王勁峰. 地圖的定性和定量分析[J]. 地球信息科學(xué)學(xué)報(bào),2009,11(2): 169-175.

[4] 禹文豪,艾廷華. 核密度估計(jì)法支持下的網(wǎng)絡(luò)空間POI點(diǎn)可視化與分析[J]. 測(cè)繪學(xué)報(bào),2015,44(1):82-90.

[5] 祝國瑞,郭禮珍,尹貢白,等. 地圖設(shè)計(jì)與編繪[M]. 武漢：武漢大學(xué)出版社,2001.

[6] ANDERSON T K. Kernel density estimation and K-means clustering to profile road accident hotspots[J]. Accident Analysis & Prevention,2009,41(3): 359-364.

[7] 閆慶武,卞正富,趙華. 人口密度空間化的一種方法[J]. 地理與地理信息科學(xué),2005，21(5): 45-48.

[8] 張桂銘,朱阿興,楊勝天,等. 基于核密度估計(jì)的動(dòng)物生境適宜度制圖方法[J].生態(tài)學(xué)報(bào),2013,33(23):7590-7600.

[9] XIE Z,YAN J. Kernel Density Estimation of traffic accidents in a network space[J]. Computers,Environment and Urban Systems. 2008,32(5): 396-406.

[10] XIE Z,YAN J. Detecting traffic accident clusters with network kernel density estimation and local spatial statistics: an integrated approach[J]. Journal of Transport Geography. 2013,31(0): 64-71.

[11] OKABE A,SATOH T,SUGIHARA K. A kernel density estimation method for networks,its computational method and a GIS-based tool[J]. International Journal of Geographical Information Science. 2009,23(1): 7-32.

[12] ArcGIS10.2在線幫助文檔中核密度估計(jì)函數(shù)的默認(rèn)帶寬設(shè)定[EB/OL]. http://resources.arcgis.com/en/help/main/10.2/index.html#/How_Kernel_Density _works/009z00000011000000/ . (2014,10).

[13] BAILEY T C,GATRELL A C. Interactive spatial data analysis[M]. Longman Scientific & Technical Essex,1995.

[14] 王遠(yuǎn)飛,何洪林. 空間數(shù)據(jù)分析方法[M]. 北京：科學(xué)出版社,2007.

[15] WILLIAMSON D,MCLAFFERTY S,GOLDSMITH V,et al. A better method to smooth crime incident data[J]. ESRI ArcUser Magazine. 1999: 1-5.

[16] SCOTT D W. Multivariate density estimation: theory,practice,and visualization[M]. John Wiley & Sons,2015.

[17] SIMONOFF J S. Smoothing methods in statistics[M]. Springer Science & Business Media,2012.

[18] 徐智勇,艾廷華,危擁軍，等. 三維地圖符號(hào)視覺參量研究[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2006，31(6): 557-560.

[19] 葉思源,吳樹仁,歐陽永龍. 地質(zhì)三維可視化建模與其剖面自動(dòng)制圖應(yīng)用研究[J]. 地質(zhì)與勘探,2011,47(3):498-504.

[20] 劉銳,胡偉平,王紅亮,等. 基于核密度估計(jì)的廣佛都市區(qū)路網(wǎng)演變分析[J]. 地理科學(xué),2011,31(1):81-86.

[責(zé)任編輯：劉文霞]

Applications of pointgroup density cartography based on kernel density estimation

LU Min1,2, YANG Liu3, WANG Jinyin1,2, HUANG Huang1,2, WANG Jiechen1,2

(1. Jiangsu Province Key Laboratory of Geographic Information Science and Technology, Nanjing University, Nanjing 210023, China;2. Dept.of Geographic Information Science, Nanjing University, Nanjing 210023,China;3. Jiangsu Province Surveying & Mapping Engineering Institute, Nanjing 210013,China)

With the development of information technology and the increase of thematic point data, it is necessary to break through the limitation of the traditional cartography and apply mass point group data to express thematic information more quickly, intuitively and comprehensively. This paper, taking ATM network data of banks in Nanjing and KTV network data in the whole country as given data, studies the methods to make kernel density map and 3D kernel density surface map and presents the key technology further. The study shows that the point density cartography based on the kernel density eliminates the influence of edge on the accuracy of calculating density and reflects the density of point distribution, the regional difference, dynamic change and other characteristics, so which has more abundant expression.

kernel density estimation; point group density cartography; optimal bandwidth computation; 3-D kernel density surface map

引用著錄：盧 敏,楊 柳,王金茵，等.基于核密度估計(jì)的點(diǎn)群密度制圖應(yīng)用研究[J].測(cè)繪工程，2017，26(4)：70-74，80.

10.19349/j.cnki.issn1006-7949.2017.04.013

2016-02-25；

2016-03-19

國家自然科學(xué)基金資助項(xiàng)目(41571377)；國家自然科學(xué)基金資助項(xiàng)目(41401450)

盧 敏(1991-)，男，碩士研究生.

P283

A

1006-7949(2017)04-0070-05