天天看點

新版白話空間統計(25):方向分布(标準差橢圓)

方向分布是蝦神最喜歡的一個空間統計工具,也是最簡單明了,但是用處很廣的一個

點模式的分析中,一般會考察如下五種内容:

新版白話空間統計(25):方向分布(标準差橢圓)

1、點的疏密,包括點資料的分布探索,是否一緻、均勻或者不均勻。

2、點的方位,包括點的分布和方向。

3、點的數量:多少(極值和均值)。

4、點的大小:代表的含義(如點一個點代表多少人口)。

5、其他,如點的一些動态變化等。

我們前面說的好幾種算法,如中心要素、中位數中心和平均中心,都是關于點方位的分析,那麼今天我們要講的這個算法,就是同時對點的方向和分布進行分析的一種經典算法——标準差橢圓。

這算法最早是由美國南加州大學(University of Southern California)社會學教授韋爾蒂.利菲弗(D. Welty Lefever)在1926年提出,是以有的書裡面,也把這個算法稱為Lefever's "Standard Deviational Ellipse"(利菲弗方向性分布)(又到每天的曆史起源科普時間……)。

這個算法最大的特點,就如同他的名詞一樣,是用來度量一組資料的方向和分布的,生成的結果又正如他的别名一樣,會輸出一個橢圓,如下:(在ArcMap中,僅能分析二維資料,在ArcGIS Pro中,支援三維點資料的分析,生成的結果是一個三維體)

新版白話空間統計(25):方向分布(标準差橢圓)

今天我們就看看這個神奇的橢圓代表的啥意思。

橢圓是一個非常神奇的幾何圖形——如果說(哲學意義)圓不僅是幾何圖形,它是未分化的統一性的象征,在神秘主義哲學中,它代表了宇宙創生一分為二時刻之前的那個宇宙。那麼橢圓就是開天辟地之後的萬物運作軌迹——是以開普勒關于行星的運動定律裡面,所有的軌迹都是橢圓形的。

新版白話空間統計(25):方向分布(标準差橢圓)

關于橢圓的方程比圓要麻煩一點,但是隻要有幾個标準參數确定了,也就很容易繪制出來的,繪制一個橢圓需要三個标準的參數:

1、确定圓心。

2、确定旋轉角度。

3、确定XY軸的長度。

算法公式都很清晰,大家自己看看就知道了:

新版白話空間統計(25):方向分布(标準差橢圓)

其中 x 和 y 是 i要素的坐标,{x̄, ȳ} 表示要素的平均中心,n 為要素總數。

确定了中心點和角度之後,計算長短軸的标準差:

新版白話空間統計(25):方向分布(标準差橢圓)

之後就可以直接用橢圓方程繪制出橢圓來了……

好吧,我知道有很多同學有數學恐懼症,是以簡單列一下就好,詳細的請參閱官方文檔,我也就不解釋了

(題外話:新版本的ArcGIS幫助文檔裡面終于改用矩陣方程了……喜大普奔啊,老版本的幫助文檔裡面的公式那個亂的喲)

新版白話空間統計(25):方向分布(标準差橢圓)

如果是學術界或者要寫論文的同學,肯定希望手動演算一下這方程,那麼自己動手豐衣足食……而像蝦神這樣的工程界(數學學渣)人士,隻能通過工具來實作就可以了。

比如,最簡單的用ArcGIS來做,直接使用工具就可以了:

新版白話空間統計(25):方向分布(标準差橢圓)

方向分布工具非常簡單,必選參數就是兩個,一個輸入要素,輸入你要分析的點資料,一個輸出的資料,就是結果的橢圓。

還有三個可選參數,後面再說。

如果我們做預設的分析,僅輸入要分析的要素,不設定權重,得到的結果如下面的圖。下面來簡單解釋一下:

新版白話空間統計(25):方向分布(标準差橢圓)

首先,要分析的資料是是長江流域的傷寒病的案例,在圖上用紅色的點來表示,藍色的河流是長江太湖流域。

當我們使用預設參數進行分析的時候,會得到這樣一個橢圓(上圖粉紅色的橢圓),并且這個橢圓會有如下屬性:

新版白話空間統計(25):方向分布(标準差橢圓)

其中,Shape_Leng和Shape_Area是生成的橢圓的周長和面積,機關與你資料的機關相同,這裡我的資料是經緯度的,是以生成的結果隻能作為相對參考結果。

(空間分析的時候,如果需要比較精确的測量屬性,

一定要使用投影坐标系!

一定要使用投影坐标系!!

一定要使用投影坐标系!!!

重要的事情說三遍)

CenterX和CenterY表示的是橢圓的中心點。

XstdDist和YStdDist表示的X軸的長度和Y軸的長度。

Rotation表示的是橢圓的方向角度。如下:

新版白話空間統計(25):方向分布(标準差橢圓)

結果解讀如下:

1、橢圓的長半軸表示的是資料分布的方向,短半軸表示的是資料分布的範圍,長短半軸的值差距越大(扁率越大),表示資料的方向性越明顯。反之,如果長短半軸越接近,表示方向性越不明顯。如果長短半軸完全相等,就等于是一個圓了,圓的話就表示沒有任何的方向特征。

2、短半軸表示資料分布的範圍,短半軸越短,表示資料呈現的向心力越明顯;反之,短半軸越長,表示資料的離散程度越大。同樣,如果短半軸與長半軸完全相等了,就表示資料沒有任何的分布特征。

3、中心點表示了整個資料的中心位置,一般來說,隻要資料的變異程度不是很大的話,這個中心點的位置大約與算數平均數的位置基本上是一緻的,至于資料變異是什麼情況,請看下面第4點。

4、有的同學會很疑惑,為什麼你畫的這個橢圓,還有很多的點都在外面,沒有把所有的點都包含進去?那麼就是就是“标準差橢圓”這個名詞裡面的“标準差”的含義所在了。

在ArcGIS工具裡面(其他的工具也都差不多),提供了“橢圓大小”(Ellipse_Size)這個參數,這個參數表示你生成的橢圓的級别,一共有三個,如下表:

新版白話空間統計(25):方向分布(标準差橢圓)

三個級别的橢圓,分别表示了你生成的橢圓,能夠包含68%,95%和99%三個級别的資料,我們通過可以指定要表示的标準差數(1、2 或 3)來決定你生成的橢圓包含的資料比例。

當要素具有空間正态分布時(即這些要素在中心處最為密集,而在接近外圍時會逐漸變得稀疏),第一級标準差(預設值)範圍可将約占總數 68% 的輸入要素的質心包含在内。第二級标準差範圍會将約占總數 95% 的要素包含在内,而第三級标準差範圍則會覆寫約占總數 99% 的要素的質心。

是以,當你選擇不同标準差等級的時候,你發現你的中心點的位置也可能不同。

當然,作為空間分析工具,方向分布一樣可以進行權重計算,這個計算主要還是與中心點的位置确定以及橢圓标準差等級生成的橢圓大小有關系。

下面我們來通過一個執行個體來了解方向分布工具的應用:

一共有兩年的傷寒病資料,如下,紅色的是2000年的,藍色是2001年的:

新版白話空間統計(25):方向分布(标準差橢圓)

使用1個标準差的結果,生成的橢圓如上,具體資料如下:

新版白話空間統計(25):方向分布(标準差橢圓)

我們對關鍵性的兩個名額進行對比:

新版白話空間統計(25):方向分布(标準差橢圓)

生成面積表示範圍,可以發現2000的傷寒病發病情況的範圍要大于2001年的,扁率表示他的方向明确性和向心力的程度,2000年生成的橢圓扁率遠高于2001年的,說明了2000的傷寒發病情況,比2001年的方向趨勢更明顯。

2000年生成的橢圓方向與長江的方向基本相相符,是以2000年的傷寒發病源與長江的關系較為密切。而2001年的發病情況雖然方向性并不如2000年的明顯,但是長半軸幾乎與長江垂直,可以基本确認,2001年傷寒發病源與長江關系不大,但是呈現向内陸爆發的趨勢。

方向分布工具在空間統計中是綜合能力最突出的工具之一,有着廣泛的應用,在我們的分析和資料探索的時候,能夠起到非常重大的作用。

最後,列出幾個可能的應用:

1、可用來在地圖上标示一組犯罪行為的分布趨勢,并且能夠确定該行為與特定要素(一系列酒吧或餐館、某條特定街道等)的關系。

2、在地圖上标示地下水井樣本的特定污染,可以訓示毒素的擴散方式,這在部署應急防災政策時非常有用。

3、對各個物種所在區域的橢圓的大小、形狀和重疊部分進行比較可以分析與物種入侵或者隔離相關的深入資訊。

4、繪制一段時間内疾病爆發情況的橢圓可用于建立疾病傳播的模型。

以上案例,以後有機會會進行展示。

最後,說說使用R語言進行方向分布的分析方法,主要使用的是aspace這個包:

新版白話空間統計(25):方向分布(标準差橢圓)

展示結果如下:

新版白話空間統計(25):方向分布(标準差橢圓)

繼續閱讀