余卓芮 孫怡恒 劉鈺 沈國強 劉岱寧

關鍵詞：開源；R語言；空間分析；教學改革；廣東省

1 空間分析教學的現(xiàn)狀及突出問題

空間分析是地理信息系統(tǒng)課程教學的核心環(huán)節(jié)。通過學習空間數(shù)據(jù)的分析方法，學生可以揭示地理對象在空間分布、空間形態(tài)、空間演變等方面的特征，并基于分析結果和專題地圖為現(xiàn)實問題的決策提供技術支持。本環(huán)節(jié)擬達到的教學目標為，通過對一系列空間分析模型的講授，使學生在理解模型算法的基礎上，具備根據(jù)現(xiàn)實問題選取恰當分析工具的能力，并且能夠對分析結果進行科學解讀，從而為現(xiàn)實問題的解決提供技術支持。但是，課堂測驗和期末考試結果顯示空間分析環(huán)節(jié)的教學效果往往較差。其中的突出問題包括：學生不能深入理解空間分析工具的計算原理，導致在分析工具的選擇上具有盲目性，進而在分析結果的解讀上存在較大偏差，最終無法將教學內容應用于對現(xiàn)實問題的分析。

與相關院校教學團隊的深入交流以及與學生的廣泛溝通中發(fā)現(xiàn)，空間分析教學中出現(xiàn)的上述問題與教學使用的分析軟件有密切關系。ESRI公司出品的ArcGIS系列軟件是全球地理信息系統(tǒng)領域的領軍產(chǎn)品，在市場占有率和行業(yè)認可度上都遠遠高于同類產(chǎn)品，因此也被廣泛用于空間分析教學中。但ArcGIS應用于教學存在明顯的局限性。除了價格高昂，還突出表現(xiàn)在作為商業(yè)軟件，其源代碼是不公開的。這使得學生只能被動地使用分析模型，模型中的很多參數(shù)不能根據(jù)現(xiàn)實需求進行調整、也很難根據(jù)現(xiàn)實需求引入或構建新的分析模型，更不能通過查看源代碼來深入理解分析模型是如何運行的。這就導致學生在分析模型的選擇上具有盲目性，進而在分析結果的解讀上也往往存在偏差，甚至影響學生深入探索的積極性。由此可見，使用不開放源代碼的商業(yè)軟件，是造成空間分析教學效果較差的重要原因。

為解決這一問題，開放源代碼的R語言被嘗試引入空間分析教學。R語言是一種面向數(shù)據(jù)分析的計算機編程語言，具有免費開源、擴展性強和簡便易學等優(yōu)點，已被用于統(tǒng)計學等專業(yè)的教學中[1]。相對于不開放源代碼的商業(yè)軟件，R語言具有以下突出優(yōu)勢：1） R語言社區(qū)提供了sf和terra等擴展程序包，支持空間數(shù)據(jù)的讀取、分析和可視化，在功能上完全滿足空間分析的教學需要；2） R語言及其擴展程序包的源代碼是開放的，學生可以自由地閱讀和修改空間分析模型代碼，從而可以深入地理解分析模型的原理和計算過程，也可以針對現(xiàn)實需求自行調用、修改和編寫分析模型。這使得學生由被動地使用軟件轉為主動地對分析模型進行探索。由于R語言具有明顯優(yōu)勢，已有教師嘗試將R語言引入地理信息專業(yè)的教學中[2]。教學實踐表明，通過將開源R語言引入空間分析教學，可以提高學生學習的積極性，教學效果也得到了明顯提升。

2 開源R 語言建模引入空間分析教學的主要途徑

如前所述，開源R語言作為空間分析教學工具有著突出的優(yōu)勢。將R語言建模引入空間分析教學，主要包括以下三種典型途徑：

1） 將其他領域的分析模型引入空間分析教學。在理論和現(xiàn)實需求的推動下，很多其他學科的模型被不斷引入空間分析領域。由于商業(yè)軟件需要控制軟件的規(guī)模和維持軟件的穩(wěn)定性，不會將所有的相關模型都納入軟件產(chǎn)品。由此導致，很多空間分析實踐中經(jīng)常使用的模型，并沒有納入以ArcGIS為代表的商業(yè)軟件中。以Theil指數(shù)為例，這是經(jīng)濟學領域用以測度收入不平等程度的基礎模型。近年來它被引入空間分析領域，成為空間差異分析中經(jīng)常使用的模型[3]，但是卻沒有被納入ArcGIS軟件中。基于R語言編寫Theil指數(shù)的計算函數(shù)并引入教學過程，可以有效地支持空間差異分析的案例教學。

2） 開放空間分析模型的參數(shù)?？臻g分析模型的運行往往需要輸入多個參數(shù)，如果將所有參數(shù)都開放，分析軟件的界面將變得臃腫不堪，因此相關軟件往往會將部分參數(shù)進行封閉。以測度空間中心聚集性的標準距離模型為例，ArcGIS只能輸入標準距離倍數(shù)、權重字段和分組字段三個參數(shù)，包括區(qū)域中心坐標在內的其他參數(shù)是封閉的，用戶無法進行調整。對于特定的現(xiàn)實問題，區(qū)域中心的坐標往往需要由用戶來進行指定，否則分析會產(chǎn)生偏差。但是，這在Arc?GIS中是無法實現(xiàn)的。通過將開源R語言引入教學，就可以很好的解決這一難題。R語言的源代碼是公開的，通過調整代碼，可以將封閉的參數(shù)開放給學生，從而支持空間分析的教學需求。

3） 針對具體場景開發(fā)新的分析模型。不同的應用場景對空間分析模型的需求是不同的，因此經(jīng)常需要針對現(xiàn)實需求定制適用的分析模型。雖然ArcGIS 等成熟的商業(yè)軟件提供了腳本語言，可以通過編程定制新的分析模型。但是，腳本語言屬于“膠水”性質的語言，只能用以“黏合”已有的模型，并不適于創(chuàng)建新的分析模型。而在R編程語言的支持下，學生可以針對具體的應用場景，構建適用的新型模型，并對分析結果進行解讀和驗證。這可以極大地提高學生的參與度和積極性，從而顯著地提高教學效果。

3 基于R 語言建模的空間分析教學案例

如前節(jié)所述，將R語言引入空間分析教學有三種途徑?；诮虒W實踐中使用的廣東省土地城市化案例，現(xiàn)對這三種途徑進行具體介紹。本案例基于Zhang等提取的不透水地表數(shù)據(jù)[4]，采用Li等提出的技術框架獲得了廣東省的城市用地范圍[5]，并通過楊喜等提出的土地城市化計算方法[6]，計算了廣東省124 個縣級行政單元（包括縣、縣級市及市轄區(qū)）的土地城市化率。基于R語言建模，可對土地城市化的空間格局進行深入分析。

3.1 引入經(jīng)濟學領域模型的空間差異分析

空間差異是指一定地域內的實體普遍存在差異性。例如，廣東省一百多個縣級行政單元在土地城市化水平上存在巨大差異，很少出現(xiàn)兩個或多個行政單元的土地城市化水平完全相同的情況。

空間差異的程度可以用Theil指數(shù)來測度。Theil 指數(shù)早期被經(jīng)濟學家用來測度收入不平等的程度，近年來被引入空間分析領域[3]，其公式為：

式中xi 和-x分別為各行政單元的土地城市化率及其均值，n 為行政單元的數(shù)量。Theil指數(shù)的取值范圍為0～ln（n），數(shù)值越大說明各行政單元之間的差異越大。

Theil指數(shù)計算結果為0.33，說明廣東省縣級行政單元的土地城市化水平存在一定的差異。為分析空間差異的來源，進一步按照行政等級將縣級單元分為兩組，其中市轄區(qū)為一組，縣和縣級市為一組。兩組的Theil指數(shù)分別為0.13和0.15，遠小于廣東省整體的Theil指數(shù)。這說明兩個行政等級分組內部的差異較小，廣東省土地城市化在空間上的差異主要來自行政等級分組之間的差異。這可歸因于我國城市的用地指標是按照行政等級分配的，由此導致行政等級較高的城市的土地開發(fā)強度也相應較高[7]。

在本案例中，基于R語言編寫了Theil指數(shù)的計算函數(shù)，將經(jīng)濟學領域的分析工具引入空間分析教學，彌補了商業(yè)軟件ArcGIS中空間差異分析工具缺失的不足。此外，通過實例分析，學生不僅測度了空間差異的程度，還探索了空間差異的來源，這對學生科學解讀空間差異的分析結果具有明顯的助益。

3.2 開放標準距離模型參數(shù)的空間中心聚集分析

空間中心聚集是指區(qū)域要素有向區(qū)域中心聚集的趨勢。例如，距離廣東省的中心城市廣州越近，各行政單元的土地開發(fā)強度也會越高，其土地城市化水平也會相應上升。

區(qū)域要素的空間中心聚集程度常采用加權標準距離模型進行分析[8]，其公式為：

式中xi 和yi 為行政單元i 在東西和南北方向上的坐標；xc 和yc 為區(qū)域中心的坐標；wi 為計算權重，在本案例中為行政單元i 的土地城市化水平。SDw 越小，說明指向中心城市的聚集程度越高。

由公式（2） 可見，區(qū)域中心的坐標對SDw 的計算結果有重要影響，所以應針對應用場景選擇適宜的區(qū)域中心。以廣東省為例，廣州市作為省會城市，是廣東省政治和文化中心，也是城市用地規(guī)模最大的城市；而且城市間的交通組織也是以廣州市為中心節(jié)點?？梢?，廣州市是廣東省城市發(fā)展的中心，在做中心聚集分析時應以廣州作為區(qū)域中心來進行計算。但是，在商業(yè)軟件ArcGIS中，區(qū)域中心的坐標被指定為輸入數(shù)據(jù)點位的加權重心坐標，是不能由用戶指定和更改的。以廣東省縣級行政單元的土地城市化分析為例，ArcGIS分析工具使用的區(qū)域中心坐標位于東莞市境內（圖1中B點），相對廣州市中心（圖1中A 點）偏移了45公里。顯然，這種簡單依賴算法推斷出來的“區(qū)域中心”與廣東省城市發(fā)展的空間結構是嚴重不符的。而且SDw 的計算結果也顯示，ArcGIS的計算結果對比以廣州市為中心的計算結果，偏差了5190 米。本案例中，區(qū)域中心廣州市臨近廣東省的重心位置，因此SDw 計算結果的偏差相對較小。不難推斷，當區(qū)域中心遠離區(qū)域重心位置時（例如，長三角地區(qū)的區(qū)域中心上海市），SDw 計算結果會出現(xiàn)更大的偏差。

從本案例可見，在商業(yè)軟件ArcGIS的空間分析工具中，部分參數(shù)是不對用戶開放的，由此導致其分析工具并不適用于所有的應用場景，在某些場景中其分析結果會存在偏差。與之相對，基于開源R語言，可以將空間分析模型中的任意參數(shù)開放給學生。學生可以根據(jù)應用需要對分析模型中的任意參數(shù)進行合理設定，從而獲得科學的分析結果?？梢姡赗語言的支持下，可以引導學生分析應用場景的實際需求并合理設定模型的參數(shù)，從而提升教學效果。

3.3 開發(fā)新模型的空間異常分區(qū)分析

空間異常分區(qū)是指明顯偏離全域整體水平的特定局域。例如，相對于廣東全省，中心城市廣州在城市規(guī)模和土地開發(fā)強度等方面都處于較高水平，屬于具有明顯優(yōu)勢的異常分區(qū)。

對特定區(qū)域的異常水平進行科學測度，是識別異常分區(qū)、分析其形成機制和制訂引導政策的基礎。但是，商業(yè)軟件ArcGIS提供的熱點分析工具只能對異常點位進行識別，缺乏分區(qū)尺度上的異常水平測度工具?；陂_源R語言，在教學中可以開發(fā)針對特定區(qū)域的異常水平測度模型，公式為：

式中l(wèi)i為特定分區(qū)內行政單元i 的屬性特征，在本案例中為土地城市化水平；-l 為特定分區(qū)內各個行政單元屬性的均值；-g為全域各個行政單元屬性的均值。式（3） 借鑒通行的歐幾里得距離測算方法，其分子和分母分別測度了特定分區(qū)內各個行政單元相對于分區(qū)均值和全域均值的偏離程度。異常水平指數(shù)Di 的取值范圍為0～1；特定分區(qū)內部各單元與分區(qū)均值越接近、與全域均值的偏離越大，Di 的取值越低，表明特定分區(qū)的異常水平越高。

分析結果顯示，包含11個市轄區(qū)單元的廣州市，其Di 為0.34。由此推斷，廣州市的土地城市化水平與廣東省的整體水平之間存在較大偏差，廣州市為土地城市化的異常區(qū)域。由圖2可見，廣州各市轄區(qū)土地城市化率的均值為69%，是廣東全省均值33%的兩倍有余，廣州市的土地城市化水平明顯高于廣東省的平均水平。具體來看，廣州市11個市轄區(qū)中，僅有2個低于廣東省均值，其余9個市轄區(qū)均高于廣東省均值。按照公式（3） 中分母和分子的分別測算，廣州各市轄區(qū)相對于廣東省均值的偏離水平（221%） 遠高于其相對于廣州市均值的偏離水平（76%） ?？梢?，廣州市各市轄區(qū)的土地城市化水平存在向自身均值集中的態(tài)勢，而與廣東省的平均水平產(chǎn)生了較大偏離；而且其土地城市化水平明顯高于廣東省的平均水平。因此，廣州市是廣東省城市土地開發(fā)的顯著熱點區(qū)域。

以上結果，從城市土地發(fā)展機制上看，反映了廣州市作為廣東省政治文化和人口聚集的中心，對城市土地有巨大的現(xiàn)實需求，并且在城市土地開發(fā)上得到了廣東省政府的政策支持。從政策制定上看，廣州市應適當控制城市土地開發(fā)強度，避免過度開發(fā)導致的生態(tài)損失和省內發(fā)展不平衡問題。

在本教學案例中，采用R語言開發(fā)了針對特定分區(qū)的異常水平測度模型，豐富了空間分析的模型庫。該模型本質上是將分區(qū)內各行政單元對本分區(qū)均值的偏離程度，與這些行政單元相對于全域均值的偏離程度進行比較。如果分區(qū)內部各單元趨向于向本分區(qū)的均值集中，而與全域均值形成較大的偏離，這說明該分區(qū)具有較高的異常水平。在指導學生進行新型模型開發(fā)的過程中，學生掌握了如何將傳統(tǒng)統(tǒng)計學中對均值的偏離程度測度方法引入空間分析。相對于對新型模型開發(fā)有諸多限制的封閉的商業(yè)軟件，開源R語言為引導學生深入理解分析模型的算法原理，為他們提高學習的參與度和積極性提供了有力支持。

4 結論

從以上教學案例可見，商業(yè)軟件的封閉性導致學生在深入理解和合理選擇空間分析工具時存在困難，甚至難以科學解讀分析結果，從而導致課堂教學效果較差。與之相對，通過開放源代碼的R語言，學生能夠自由地查看、修改和編寫分析工具的代碼，從而得以深入理解空間分析模型的算法，并能夠合理選擇和科學解讀其分析結果，從而極大地提高學生學習的參與度和積極性，使得教學效果得到了顯著提升。

教學實踐表明，基于R語言可以將其他領域的分析工具引入空間分析教學，以及改進現(xiàn)有的空間分析模型，還可以根據(jù)應用場景構建新的分析模型。通過這三種途徑，R語言可為空間分析教學提供全面和有力的支持。除了以空間分析為核心的地理信息系統(tǒng)課程教學，R語言強大的建模能力還可以進一步應用于城市地理、區(qū)域規(guī)劃等課程，為這些課程的定量化教學改革提供技術支持。

隨著大數(shù)據(jù)時代的到來，從海量數(shù)據(jù)中揭示研究對象的特征和演化規(guī)律已成為解決現(xiàn)實問題需具備的核心技能。在經(jīng)濟生產(chǎn)和社會發(fā)展中，對具備相關技能的人才具有廣泛而迫切的需求[1]。因此，廣大教師應積極將R語言引入一線教學工作，結合現(xiàn)實需求推進教學改革，引導學生成為大數(shù)據(jù)時代的棟梁之才。