夏伊凡 黃曉燕 康晨晨 徐旭

摘要 傳統(tǒng)村落是推進(jìn)鄉(xiāng)村振興和新型城鎮(zhèn)化建設(shè)的重要載體。以陜西省182個(gè)國家級(jí)傳統(tǒng)村落為研究對(duì)象，分析其空間分布特征，并運(yùn)用ArcGIS網(wǎng)絡(luò)分析和地理加權(quán)回歸模型，深入剖析縣級(jí)和鎮(zhèn)級(jí)尺度下傳統(tǒng)村落可達(dá)性特征及影響因素的空間分異。研究結(jié)果表明：① 陜西省國家級(jí)傳統(tǒng)村落呈集聚型空間分布，在空間上形成了4個(gè)主要集聚區(qū)。② 傳統(tǒng)村落可達(dá)性整體水平較低，呈現(xiàn)以關(guān)中為中心的階梯狀遞減特征。③ 不同空間尺度下的傳統(tǒng)村落可達(dá)性與影響因素存在著較強(qiáng)的空間不穩(wěn)定性?？h級(jí)尺度下，高程與傳統(tǒng)村落可達(dá)性基本呈正相關(guān)關(guān)系，但在鎮(zhèn)級(jí)尺度下，兩者正負(fù)向關(guān)系并存，且呈現(xiàn)以榆林西北部和渭南中部為中心向外圍圈層逐漸增強(qiáng)的特點(diǎn)?？h級(jí)尺度下，近水體最短距離與傳統(tǒng)村落可達(dá)性水平呈負(fù)相關(guān)關(guān)系，但在鎮(zhèn)級(jí)尺度下，其對(duì)可達(dá)性水平的影響呈現(xiàn)出從陜南和陜北向中部逐步遞減的趨勢。④高程、坡度和近水體最短距離是影響陜西省國家級(jí)傳統(tǒng)村落可達(dá)性的主要因素，且在不同地區(qū)具有不同的影響效應(yīng)。傳統(tǒng)村落可達(dá)性特征和影響因素的研究可以為不同地區(qū)傳統(tǒng)村落可達(dá)性的改進(jìn)和傳統(tǒng)村落保護(hù)利用策略的制定提供科學(xué)依據(jù)。

關(guān)鍵詞 傳統(tǒng)村落;可達(dá)性;GWR模型;不同尺度;陜西省

Spatial differentiation of accessibility characteristics and influencing

factors of traditional villages in Shaanxi Province

Abstract Traditional villages are an important carrier for promoting rural revitalization and new urbanization construction. The spatial distribution characteristics of 182 national-level traditional villages in Shaanxi Province were analyzed as research objects. Using ArcGIS network analysis and GWR model， the spatial differentiation of traditional village accessibility characteristics and influencing factors at the county and town scales were analyzed in depth. The results of the study show that： ① The national-level traditional villages in Shaanxi Province are characterized by agglomeration spatial distribution， forming four main agglomeration area in space. ② The overall level of accessibility is low， showing a stepwise decline centered on Guanzhong. ③ There is strong spatial instability in the accessibility and influencing factors of traditional village accessibility at different spatial scales. At the county level， the elevation has a positive correlation with the accessibility of traditional villages. However， at the town level， both positive and negative relationships coexist， and are characterized by a gradual increase from the northwestern part of Yulin and the central part of Weinan to the peripheral circles. At the county level， the distance to neighboring rivers is negatively correlated with the level of accessibility of traditional villages. However， at the town level， its impact on accessibility shows a decreasing trend from southern and northern Shaanxi to the central part of the country. ④ Elevation， slope， and the distance to neighboring rivers are the main factors in the accessibility of national-level traditional villages in Shaanxi Province and have different influencing effects in different areas. Research on the accessibility of traditional villages can help reveal the ways to improve the accessibility of traditional villages in different regions， and can provide a scientific basis for the formulation of future conservation and utilization strategies for traditional villages.

Keywords traditional villages; accessibility; geographically weighted regression model; different scales; Shaanxi Province

傳統(tǒng)村落是指形成較早，擁有較豐富的傳統(tǒng)資源，具有一定歷史、文化、科學(xué)、藝術(shù)、社會(huì)和經(jīng)濟(jì)價(jià)值，應(yīng)予以保護(hù)的村落［1］。傳統(tǒng)村落作為傳統(tǒng)文化景觀的重要載體［2］，是推進(jìn)鄉(xiāng)村振興和新型城鎮(zhèn)化建設(shè)極為重要的資源。長期以來，我國高度重視傳統(tǒng)村落的保護(hù)與傳承。自2012年起，全國開展了6次傳統(tǒng)村落摸底調(diào)查，截至2022年底，共認(rèn)定傳統(tǒng)村落8 171個(gè)。然而，在社會(huì)經(jīng)濟(jì)快速發(fā)展的沖擊下，當(dāng)前傳統(tǒng)村落保護(hù)和利用仍面臨著空心化、過度商業(yè)化、民族傳統(tǒng)文化遺失、非物質(zhì)文化景觀衰落和生態(tài)環(huán)境惡化等一系列問題，傳統(tǒng)村落的保護(hù)和利用迫切需要相關(guān)研究的支撐。

我國傳統(tǒng)村落研究起步相對(duì)較晚，但成果豐富。學(xué)者從建筑學(xué)、地理學(xué)、生態(tài)學(xué)、管理學(xué)等不同學(xué)科視角，對(duì)傳統(tǒng)村落開展了大量研究。研究內(nèi)容主要集中在傳統(tǒng)村落空間格局及形態(tài)特征［2-5］、價(jià)值評(píng)價(jià)［6-10］、保護(hù)與利用策略［11-14］、旅游發(fā)展及影響因素［15-18］和景觀基因［19-20］等方面?？傮w而言，現(xiàn)有研究側(cè)重于對(duì)傳統(tǒng)村落的空間分布、價(jià)值解讀和開發(fā)策略的討論，對(duì)傳統(tǒng)村落空間可達(dá)性的深層次分析相對(duì)不足。

便攜的交通是地區(qū)發(fā)展的前提條件。傳統(tǒng)村落大多分布在位置偏遠(yuǎn)、交通基礎(chǔ)設(shè)施相對(duì)較差的地區(qū)，傳統(tǒng)村落的保護(hù)利用高度依賴于交通條件和可進(jìn)入性。傳統(tǒng)村落的可達(dá)性水平反映了其進(jìn)入的交通便利程度，會(huì)對(duì)傳統(tǒng)村落的社會(huì)經(jīng)濟(jì)發(fā)展與文化保護(hù)傳承產(chǎn)生重要影響。準(zhǔn)確把握傳統(tǒng)村落可達(dá)性空間分布特征及其影響因素是對(duì)傳統(tǒng)村落進(jìn)行高效保護(hù)利用的前提。當(dāng)前的可達(dá)性研究主要關(guān)注區(qū)域尺度的可達(dá)性測度、格局變化及其對(duì)經(jīng)濟(jì)、人口或產(chǎn)業(yè)的影響［21-23］，或在城市內(nèi)部分析公共服務(wù)設(shè)施可達(dá)性差異及其效應(yīng)［24-25］。關(guān)于傳統(tǒng)村落交通條件和可達(dá)性的研究相對(duì)有限，且已有研究多側(cè)重于采用緩沖區(qū)或空間句法分析，探討傳統(tǒng)村落可達(dá)性空間分布特征。如張忠訓(xùn)等［26］運(yùn)用空間句法模型分析得出貴州少數(shù)民族聚居區(qū)市域內(nèi)傳統(tǒng)村落全局可達(dá)性與傳統(tǒng)村落分布成正比，局部可達(dá)性高的傳統(tǒng)村落普遍具有極為便利的高等級(jí)公路通達(dá)性。胡金龍等［27］運(yùn)用空間句法及GIS空間分析方法研究了桂林市傳統(tǒng)村落可達(dá)性特征，發(fā)現(xiàn)自然環(huán)境復(fù)雜、社會(huì)發(fā)展滯后區(qū)域的傳統(tǒng)村落可達(dá)性相對(duì)更低。竇銀娣等［28］采用時(shí)間成本分析法研究發(fā)現(xiàn)湖南省傳統(tǒng)村落的整體可達(dá)性較好，市域尺度的傳統(tǒng)村落可達(dá)性差異較大。在傳統(tǒng)村落可達(dá)性的影響因素方面，現(xiàn)有研究主要采用相關(guān)性分析或線性回歸模型，認(rèn)為傳統(tǒng)村落可達(dá)性受海拔、坡度、離水距離、道路等級(jí)、社會(huì)經(jīng)濟(jì)發(fā)展及民族屬性等因素的共同影響［27，29-30］。然而，在傳統(tǒng)村落研究中，缺乏對(duì)可達(dá)性空間分布特征影響因素的深入挖掘，對(duì)影響因素的空間異質(zhì)性更是缺乏討論。

鑒于此，本文以陜西省國家級(jí)傳統(tǒng)村落為研究對(duì)象，運(yùn)用地理加權(quán)回歸（geographically weighted regression， GWR）模型，探究傳統(tǒng)村落可達(dá)性的空間格局、影響因素及各影響因子的空間分異，揭示傳統(tǒng)村落可達(dá)性空間分布規(guī)律和影響因素在不同地理空間背景下的影響差異，以期為傳統(tǒng)村落的保護(hù)開發(fā)提供決策依據(jù)。

1 研究數(shù)據(jù)與方法

1.1 數(shù)據(jù)來源

通過對(duì)中華人民共和國住房和城鄉(xiāng)建設(shè)部網(wǎng)站公布的第一批至第六批中國國家級(jí)傳統(tǒng)村落名錄進(jìn)行梳理，截止2022年，陜西省共有國家級(jí)傳統(tǒng)村落182個(gè)。傳統(tǒng)村落的經(jīng)緯度借助百度地圖API拾取坐標(biāo)系統(tǒng)獲取。研究中的中國海拔高度（DEM）空間分布數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心（http：∥www.resdc.cn／），空間分辨率為1 km。坡度數(shù)據(jù)通過DEM進(jìn)行提取。水系和路網(wǎng)數(shù)據(jù)源于國家基礎(chǔ)地理信息中心（http：∥www.ngcc.cn／ngcc／）提供的最新矢量數(shù)據(jù)。主要數(shù)據(jù)的處理采用了ArcGIS 10.7軟件連接、疊合等空間處理工具。

1.2 研究方法

1.2.1 空間分布特征分析方法

1）最鄰近點(diǎn)指數(shù)。最鄰近距離是表示點(diǎn)狀要素的相互鄰近距離的地理指標(biāo)，反映點(diǎn)狀要素在空間上的分布類型。計(jì)算公式為

2）核密度估計(jì)法。本文采用核密度估計(jì)法判定陜西省國家級(jí)傳統(tǒng)村落在空間上聚集的區(qū)域。核密度估計(jì)法認(rèn)為地理事件可以發(fā)生在任何空間位置上，但地理事件發(fā)生的概率隨空間位置的變化而變化，其中點(diǎn)要素越集聚的區(qū)域，地理事件發(fā)生的概率越高［32］。將陜西省182個(gè)傳統(tǒng)村落抽象成點(diǎn)狀要素，采用核密度分析法對(duì)其進(jìn)行可視化表達(dá)，核密度越高意味著傳統(tǒng)村落在空間上分布越集聚。

式中：p為待估計(jì)傳統(tǒng)村落的位置；pi為落在以p為圓心、r為半徑的圓內(nèi)的第i個(gè)傳統(tǒng)村落的位置;r表示步長，即以p為源點(diǎn)的曲面在空間上的延展的寬度。

1.2.2 可達(dá)性計(jì)算方法

可達(dá)性旨在對(duì)空間上某一要素實(shí)體（點(diǎn)、線或區(qū)域）的位置優(yōu)劣程度進(jìn)行度量。本文采用最短時(shí)間距離模型來測算陜西省國家級(jí)傳統(tǒng)村落到研究范圍內(nèi)不同層級(jí)行政中心（縣和鎮(zhèn)）的便利程度。為便于分析，對(duì)結(jié)果數(shù)值取倒數(shù)，數(shù)值越大，說明該傳統(tǒng)村落可達(dá)性越好。公式為

式中：Ai表示國家級(jí)傳統(tǒng)村落i的可達(dá)時(shí)間;n為不同尺度下的區(qū)域節(jié)點(diǎn)總數(shù)；dij在表示傳統(tǒng)村落i通過交通網(wǎng)絡(luò)到區(qū)域節(jié)點(diǎn)質(zhì)心j的最短時(shí)間。Ai值越小，表征傳統(tǒng)村落i的可達(dá)性越好。反之則可達(dá)性水平相對(duì)較差。

本文在網(wǎng)絡(luò)數(shù)據(jù)集的基礎(chǔ)上對(duì)高速公路、國道、省道、縣道和鄉(xiāng)道等不同等級(jí)的道路設(shè)定不同的行駛速度。道路行駛速度結(jié)合陜西省道路通行情況，根據(jù)《中華人民共和國公路工程技術(shù)標(biāo)準(zhǔn)（JTGB—2006）》，按照高速公路120 km／h、國道100 km／h、省道80 km／h、縣道50 km／h、鄉(xiāng)道40 km／h進(jìn)行設(shè)定。

1.2.3 GWR模型

GWR模型是在最小二乘法（ordinary least squares，OLS）基礎(chǔ)上的進(jìn)一步拓展，其考慮了空間數(shù)據(jù)的異質(zhì)性特征，能夠?qū)⒛Ｐ偷膮?shù)估計(jì)值反映在地圖上，有助于探索空間變異特征和規(guī)律［31］。GWR模型考慮了數(shù)據(jù)的空間地理變化，在回歸參數(shù)中納入對(duì)數(shù)據(jù)地理位置的考慮，使回歸參數(shù)變成空間可觀測，從而解釋因變量隨自變量空間位置的變化［32］。計(jì)算公式如下

式中：yi為傳統(tǒng)村落i的因變量；xik為第i個(gè)傳統(tǒng)村落點(diǎn)上第k個(gè)變量的觀測值；（ui，vi）為第i個(gè)傳統(tǒng)村落的地理坐標(biāo)；β0（ui，vi）是回歸方程的截距；βk（ui，vi）表示為第i個(gè)傳統(tǒng)村落的回歸參數(shù)；k為傳統(tǒng)村落數(shù)量；εi則為誤差項(xiàng)。

根據(jù)地理學(xué)第一定律，兩個(gè)地物之間存在空間相關(guān)性，一般來說，空間相關(guān)性與兩地物間距離正相關(guān)。因此，對(duì)樣本點(diǎn)進(jìn)行加權(quán)處理，“加權(quán)”通過權(quán)重矩陣來完成，利用加權(quán)最小二乘法估計(jì)參數(shù)的數(shù)學(xué)表達(dá)式為

式中：dij為點(diǎn)（ui，yi）到點(diǎn)（uj，yj）的歐氏距離;b為帶寬（bandwidth），是用來控制權(quán)重隨距離dij之間衰減的程度的參數(shù)。選擇AIC信息準(zhǔn)則（akaike information criterion，AIC）最小化來確定最佳帶寬，公式為

式中：AAIC值表示模型的赤池信息值，根據(jù)AIC準(zhǔn)則，AIC最小的值對(duì)應(yīng)的帶寬即模型中的最優(yōu)帶寬;n表示傳統(tǒng)村落點(diǎn)的數(shù)量;tr（S）表示GWR的S矩陣的跡，是帶寬的函數(shù);σ表示誤差項(xiàng)估計(jì)的標(biāo)準(zhǔn)離差。

2 傳統(tǒng)村落空間分布及可達(dá)性特征

2.1 空間分布特征

通過計(jì)算陜西省國家級(jí)傳統(tǒng)村落的最鄰近點(diǎn)指數(shù)，發(fā)現(xiàn)呈現(xiàn)集聚型空間分布特征。計(jì)算得出的理論鄰近距離rE=25.52 km，實(shí)際鄰近距離r1=16.50 km。最鄰近指數(shù)R約為0.65<1，表明陜西省國家級(jí)傳統(tǒng)村落在空間分布上成集聚狀態(tài)。

選取帶寬150 km的核密度工具生成陜西省國家級(jí)傳統(tǒng)村落空間分布的核密度熱點(diǎn)圖。如圖1所示，由于地區(qū)地理環(huán)境不同，陜西國家級(jí)傳統(tǒng)村落在空間分布上呈現(xiàn)出較強(qiáng)的不均衡性，形成了4個(gè)主要的空間集聚區(qū)，分別是陜北中東部傳統(tǒng)村落集聚區(qū)、關(guān)中東北部傳統(tǒng)村落集聚區(qū)、陜南中南部傳統(tǒng)村落集聚區(qū)以及陜南西北部傳統(tǒng)村落集聚區(qū)。

呈現(xiàn)該空間格局分布特征主要有兩方面原因：對(duì)于關(guān)中平原地區(qū)，由于其社會(huì)經(jīng)濟(jì)發(fā)展水平較高，具有良好交通和區(qū)位條件，傳統(tǒng)村落的保護(hù)利用水平相對(duì)較高;對(duì)于陜北黃土高原及陜南秦巴山地，地形條件導(dǎo)致對(duì)外交通聯(lián)系較少，傳統(tǒng)村落受外來影響小，眾多的傳統(tǒng)村落得以保留。

2.2 可達(dá)性空間特征

傳統(tǒng)村落的可達(dá)性可衡量外界進(jìn)入村落的難易程度。 本文采用最短時(shí)間距離模型計(jì)算可達(dá)性， 分別測算182個(gè)陜西省國家級(jí)傳統(tǒng)村落到研究范圍內(nèi)107個(gè)縣城、 1 398個(gè)鄉(xiāng)鎮(zhèn)的可達(dá)性水平。

為分析傳統(tǒng)村落可達(dá)性的空間分布特征及其差異，按照最短時(shí)間距離劃分可達(dá)時(shí)間等級(jí)，利用ArcGIS軟件對(duì)陜西省國家級(jí)傳統(tǒng)村落的可達(dá)性進(jìn)行可視化處理。通過測算得出縣級(jí)和鎮(zhèn)級(jí)尺度下各傳統(tǒng)村落的可達(dá)性均值為363.97 min和377.86 min，可達(dá)性空間格局基本相似，陜西省國家級(jí)傳統(tǒng)村落的可達(dá)性等級(jí)總體上呈現(xiàn)以關(guān)中地區(qū)為中心的階梯狀向四周遞減的圈層式結(jié)構(gòu)（見圖2）。具備較高可達(dá)性的區(qū)域主要分布在關(guān)中中部的西安市、咸陽市和渭南市，而陜北地區(qū)和陜南地區(qū)的可達(dá)性相對(duì)較差。這與自然地理?xiàng)l件、道路基礎(chǔ)設(shè)施和經(jīng)濟(jì)發(fā)展水平密切相關(guān)。

同時(shí)，高可達(dá)性的傳統(tǒng)村落與陜西省經(jīng)濟(jì)地理空間格局具有明顯耦合性，這些傳統(tǒng)村落所在地區(qū)的地理區(qū)位、道路基礎(chǔ)設(shè)施和經(jīng)濟(jì)發(fā)展水平具有明顯優(yōu)勢，人口密度較高、經(jīng)濟(jì)活動(dòng)相對(duì)較為活躍。

3 傳統(tǒng)村落可達(dá)性的影響因素及空間分異

3.1 OLS回歸結(jié)果

陜西省各國家級(jí)傳統(tǒng)村落交通可達(dá)性為被解釋變量。參考已有研究［29-30］，選取傳統(tǒng)村落所在地高程、坡度、近水體最短距離、路網(wǎng)密度和河網(wǎng)密度等5項(xiàng)指標(biāo)作為解釋變量。為避免指標(biāo)出現(xiàn)多重共線性，引起模型擬合錯(cuò)誤，首先，計(jì)算回歸模型參數(shù)的方差膨脹因子（variance inflation factor， VIF）。結(jié)果表明，VIF最大值是1.28，遠(yuǎn)小于7.5，變量未出現(xiàn)多重共線性問題，5項(xiàng)指標(biāo)均可作為解釋變量。其次，建立OLS經(jīng)典線性回歸模型初步測度解釋變量對(duì)被解釋變量的影響。由表1可知，傳統(tǒng)村落縣級(jí)和鎮(zhèn)級(jí)可達(dá)性均受到高程和近水體最短距離的顯著影響，其中高程均通過了1％的顯著性檢驗(yàn)。此外，傳統(tǒng)村落縣級(jí)可達(dá)性還受到坡度的顯著影響。

本文使用Jarque-Bera檢驗(yàn)方法考察樣本數(shù)據(jù)是否具有符合正態(tài)分布的偏度和峰度的擬合優(yōu)度。縣級(jí)和鎮(zhèn)級(jí)尺度Jarque-Bera檢驗(yàn)結(jié)果不顯著，可知縣級(jí)和鎮(zhèn)級(jí)尺度下的回歸方程均服從正態(tài)分布假設(shè)，可以進(jìn)行下一步的GWR模型參數(shù)回歸。

3.2 GWR回歸結(jié)果

借助GWR模型進(jìn)一步分析地理位置對(duì)傳統(tǒng)村落可達(dá)性水平的影響。運(yùn)算得出縣級(jí)和鎮(zhèn)級(jí)尺度下傳統(tǒng)村落可達(dá)性水平的Morans I值分別為0.72和0.65，Z得分為16.32和14.68，p值均為0.00，可得知不同尺度下傳統(tǒng)村落的可達(dá)性水平存在顯著的空間集聚特征。

運(yùn)用GWR 4軟件構(gòu)建GWR模型，選擇Adptive bi-square為權(quán)重函數(shù)，帶寬確定準(zhǔn)則選擇AICc，輸出得到回歸結(jié)果（見表2）?？h級(jí)尺度和鎮(zhèn)級(jí)尺度GWR模型的AICc值分別為575.854和575.332，遠(yuǎn)小于一般OLS回歸產(chǎn)生的AICc值695.869和676.804。縣級(jí)和鎮(zhèn)級(jí)尺度的矯正模型優(yōu)度（Adjusted R2）分別為0.603和0.549，均高于OLS回歸的總體擬合優(yōu)度0.111和0.105。因此，本研究使用的GWR模型結(jié)果更優(yōu)。

為探究各變量對(duì)傳統(tǒng)村落可達(dá)性的影響，對(duì)GWR模型分析結(jié)果進(jìn)行整理，分別選取最小值、上四分位數(shù)、中位數(shù)、下四分位數(shù)、最大值和平均值6個(gè)統(tǒng)計(jì)項(xiàng)目繪制表3。GWR各變量的回歸系數(shù)的平均值反映了各因素對(duì)可達(dá)性的平均貢獻(xiàn)程度。變量的回歸系數(shù)有正有負(fù)，回歸系數(shù)為正說明該變量對(duì)傳統(tǒng)村落可達(dá)性產(chǎn)生正向影響，反之亦然。

GWR回歸結(jié)果表明， 縣級(jí)和鎮(zhèn)級(jí)尺度下， 高程、 坡度和近水體最短距離的回歸系數(shù)較大， 且回歸系數(shù)的波動(dòng)較大， 表明3個(gè)變量對(duì)傳統(tǒng)村落可達(dá)性空間分布有較大影響。 同時(shí)高程、 坡度和近水體最短距離的回歸系數(shù)有正有負(fù)， 表明這3個(gè)變量對(duì)部分地區(qū)的傳統(tǒng)村落可達(dá)性有正向促進(jìn)作用， 對(duì)部分地區(qū)的傳統(tǒng)村落可達(dá)性有負(fù)向牽制作用。

3.3 影響因素的空間異質(zhì)性

采用ArcGIS軟件對(duì)GWR模型結(jié)果中各變量的回歸系數(shù)進(jìn)行可視化分析， 得到各變量系數(shù)的空間分布圖， 進(jìn)而研究各變量影響效應(yīng)的空間差異。

3.3.1 高程

繪制山體高程與傳統(tǒng)村落的疊加圖〔見圖3（a）〕。運(yùn)用柵格提取技術(shù)得到各個(gè)傳統(tǒng)村落的高程，最高地勢屬西安市周至縣厚畛子鄉(xiāng)老縣城村，高程達(dá)1 758 m。而渭南市華陰市岳廟街道雙泉村的高程只有344 m。得出陜西省國家級(jí)傳統(tǒng)村落的高程極值差異較大。

縣級(jí)尺度下，高程的回歸系數(shù)有正有負(fù)，但高程與傳統(tǒng)村落可達(dá)性基本呈現(xiàn)正相關(guān)關(guān)系〔見圖3（b）〕。正值區(qū)集中分布在關(guān)中西部及關(guān)中渭南東北部、陜北延安市及榆林東部、陜南。負(fù)值區(qū)集中分布在關(guān)中渭南中西部、陜南安康西部、陜北榆林西部。正值區(qū)地形類型多樣，包括平原、高原及山地。

鎮(zhèn)級(jí)尺度下，高程對(duì)傳統(tǒng)村落可達(dá)性是正負(fù)向影響并存，且市際差異較大，形成以陜北榆林西北部和關(guān)中渭南中部為核心向外逐漸遞增的圈層結(jié)構(gòu)〔見圖3（c）〕。高程回歸系數(shù)負(fù)值集中在關(guān)中渭南中西部、陜北榆林南部，表明上述兩個(gè)地區(qū)的絕大部分傳統(tǒng)村落可達(dá)性受到高程的阻礙。其原因可能是榆林南部地處黃土高原丘陵溝壑區(qū)，較高的周邊地勢影響了傳統(tǒng)村落可達(dá)性，而渭南受到華山山脈影響，致使部分傳統(tǒng)村落可達(dá)性的水平較低。

將傳統(tǒng)村落可達(dá)性數(shù)值進(jìn)行對(duì)比可發(fā)現(xiàn)，村落的高程與傳統(tǒng)村落的可進(jìn)入程度具有明顯耦合關(guān)系，即高程較高的傳統(tǒng)村落的可進(jìn)入程度較差。高程對(duì)傳統(tǒng)村落可達(dá)性的影響程度較大，且影響程度與高程大小相關(guān)。

3.3.2 坡度

坡度對(duì)傳統(tǒng)村落可達(dá)性具有重要影響。將坡度與傳統(tǒng)村落空間數(shù)據(jù)在ArcGIS軟件中進(jìn)行疊加，進(jìn)一步分析坡度對(duì)傳統(tǒng)村落可達(dá)性的影響，如圖4（a）。

縣級(jí)尺度下，坡度對(duì)傳統(tǒng)村落可達(dá)性水平的影響呈現(xiàn)較為顯著的正向與負(fù)向空間分區(qū)特征〔見圖4（b）〕。回歸系數(shù)呈現(xiàn)以陜南商洛為核心向外逐漸遞減的圈層結(jié)構(gòu)，即陜南商洛受坡度的影響最強(qiáng)。陜南商洛的坡度最高為16.19°，較高的坡度降低了傳統(tǒng)村落可達(dá)性水平。由于高值區(qū)的坡度相對(duì)較高，坡度的降低對(duì)提高傳統(tǒng)村落的可進(jìn)入性水平效果顯著?；貧w系數(shù)最高值每提升1個(gè)坡度水平，傳統(tǒng)村落交通可達(dá)時(shí)間數(shù)值將減少0.088 min，即降低坡度，傳統(tǒng)村落的可進(jìn)入程度將會(huì)得到大幅度提升。

對(duì)比圖4（a）和4（b）可知，回歸系數(shù)高值區(qū)部分集中在坡度相對(duì)較高的地區(qū)?？傮w來說，分布在關(guān)中地區(qū)的傳統(tǒng)村落坡度與傳統(tǒng)村落可達(dá)性水平相耦合，即交通出行便利的關(guān)中地區(qū)傳統(tǒng)村落具有較高的可達(dá)性水平。

3.3.3 近水體最短距離

陜西省地跨黃河和長江兩大流域，全省境內(nèi)河流眾多。利用ArcGIS 10.7軟件繪制各水體與傳統(tǒng)村落的疊加圖〔見圖5（a）〕。從圖中可觀察到陜西省傳統(tǒng)村落的空間分布存在明顯的依水系發(fā)展的特征，大多傳統(tǒng)村落分布在各水系兩側(cè)，已有研究也有類似發(fā)現(xiàn)［33］。陜北、關(guān)中的傳統(tǒng)村落分布在黃河支流附近，陜南傳統(tǒng)村落聚集在長江支流——漢江、嘉陵江周圍。陜西省傳統(tǒng)村落大多集中在河邊并沿河散布。沿河流分布的傳統(tǒng)村落，可獲得河流提供的水資源進(jìn)行農(nóng)業(yè)生產(chǎn)，因而獲得了更好的生產(chǎn)和發(fā)展機(jī)會(huì)。

縣級(jí)尺度下，近水體最短距離對(duì)傳統(tǒng)村落可達(dá)性水平的影響基本呈現(xiàn)負(fù)相關(guān)效應(yīng)〔見圖5（b）〕。即河流的存在對(duì)傳統(tǒng)村落的可進(jìn)入性具有限制作用。從回歸系數(shù)的空間分布來看，正向影響主要分布在關(guān)中中部，負(fù)向影響分布在關(guān)中東部、陜南漢中及安康、陜北，主要原因是河流能為村民生計(jì)發(fā)展提供充足水源，成為傳統(tǒng)村落選址的必要條件，但河流往往是高原和山地地區(qū)道路建設(shè)的阻礙。河流的阻隔加上所處的地形導(dǎo)致道路系統(tǒng)的修筑成本增加，進(jìn)而影響可達(dá)性。因此，陜北和陜南大部分距河流較近的傳統(tǒng)村落可達(dá)性較差。

鎮(zhèn)級(jí)尺度下，近水體最短距離對(duì)傳統(tǒng)村落可達(dá)性影響的正負(fù)向效應(yīng)并存〔見圖5（c）〕?；貧w系數(shù)基本呈現(xiàn)從陜南和陜北向中部逐步遞減的趨勢。關(guān)中中部及東部、陜北延安東南部、陜南漢中中部具有正相關(guān)效應(yīng)。而關(guān)中寶雞、陜北榆林以及延安北部、陜南安康及商洛具有負(fù)相關(guān)效應(yīng)。正向效應(yīng)集中分布區(qū)的回歸系數(shù)絕對(duì)值較低，即隨著與水體最短距離的增加，該區(qū)域的可進(jìn)入性水平對(duì)近水體最短距離的敏感性減弱。

4 結(jié)論與討論

4.1 結(jié)論

傳統(tǒng)村落作為重要?dú)v史文化遺產(chǎn)和文化旅游資源，分析其可達(dá)性及影響因素對(duì)于促進(jìn)傳統(tǒng)村落所在地區(qū)的社會(huì)經(jīng)濟(jì)發(fā)展和文化傳承具有重要意義。本文以182個(gè)陜西省國家級(jí)傳統(tǒng)村落為研究對(duì)象，分析了傳統(tǒng)村落的空間分布和可達(dá)性特征，運(yùn)用回歸模型分析了傳統(tǒng)村落可達(dá)性的影響因素，并基于GWR模型深入揭示了各影響因子的空間異質(zhì)性。主要結(jié)論如下。

1）陜西省國家級(jí)傳統(tǒng)村落呈集聚型空間分布特征，在空間上形成了4個(gè)主要的空間集聚區(qū)，分別是陜北中東部傳統(tǒng)村落集聚區(qū)、關(guān)中東北部傳統(tǒng)村落集聚區(qū)、陜南中南部傳統(tǒng)村落集聚區(qū)以及陜南西北部傳統(tǒng)村落集聚區(qū)。

2）陜西省國家級(jí)傳統(tǒng)村落的可達(dá)性呈現(xiàn)出以關(guān)中為中心階梯狀遞減的趨勢。各國家級(jí)傳統(tǒng)村落的可達(dá)時(shí)間等級(jí)在空間上呈連續(xù)的圈層式分布，高可達(dá)時(shí)間等級(jí)的傳統(tǒng)村落主要分布在關(guān)中中部地區(qū)，而其他可達(dá)時(shí)間等級(jí)的傳統(tǒng)村落則依次分布于外圍地區(qū)。

3）不同空間尺度下，傳統(tǒng)村落可達(dá)性與影響因素間存在著較強(qiáng)的空間不穩(wěn)定性?？h級(jí)尺度下，高程對(duì)傳統(tǒng)村落可達(dá)性基本呈現(xiàn)正相關(guān)關(guān)系;坡度對(duì)傳統(tǒng)村落可達(dá)性水平的影響呈現(xiàn)以陜南商洛為核心向外逐漸遞減的圈層結(jié)構(gòu);近水體最短距離對(duì)傳統(tǒng)村落可達(dá)性水平的影響基本呈現(xiàn)負(fù)相關(guān)效應(yīng)。鎮(zhèn)級(jí)尺度下，高程對(duì)傳統(tǒng)村落可達(dá)性水平的影響呈現(xiàn)以榆林西北部和渭南中部為中心向外圍圈層逐漸增強(qiáng)的特點(diǎn);近水體最短距離對(duì)傳統(tǒng)村落可達(dá)性水平的影響基本呈現(xiàn)從陜南和陜北向中部逐步遞減的趨勢。

4）高程、坡度和近水體最短距離是陜西省國家級(jí)傳統(tǒng)村落可達(dá)性的主要影響因素，且在不同地區(qū)具有不同的影響效應(yīng)。高程是影響關(guān)中傳統(tǒng)村落可達(dá)性的主要因素;影響陜北地區(qū)傳統(tǒng)村落可達(dá)性的因素較多，包括高程以及近水體最短距離;坡度是影響陜南傳統(tǒng)村落可達(dá)性的主要因素。

4.2 討論

生態(tài)文明建設(shè)和鄉(xiāng)村振興背景下，傳統(tǒng)村落的保護(hù)和利用面臨著前所未有的機(jī)遇和挑戰(zhàn)。交通作為地區(qū)發(fā)展的前提條件，是支撐傳統(tǒng)村落保護(hù)和利用的基礎(chǔ)。完善傳統(tǒng)村落交通基礎(chǔ)設(shè)施，提升可達(dá)性水平對(duì)于傳統(tǒng)村落的保護(hù)利用至關(guān)重要。結(jié)合本文發(fā)現(xiàn)，從可達(dá)性的影響因素入手，可探究不同傳統(tǒng)村落保護(hù)與振興的差異化路徑。

高程是影響關(guān)中傳統(tǒng)村落可達(dá)性的主要因素。對(duì)于高程較低、容易進(jìn)入的傳統(tǒng)村落，可以采取優(yōu)先發(fā)展策略。相關(guān)管理部門應(yīng)積極整合各方資源，充分挖掘傳統(tǒng)村落的文化資源，形成有機(jī)整體，全力打造生態(tài)旅游、全域旅游。對(duì)高程較高的傳統(tǒng)村落，需要加強(qiáng)對(duì)其獨(dú)特優(yōu)勢的深入挖掘與宣傳，同時(shí)通過旅游線路規(guī)劃設(shè)計(jì)，充分依托周邊景點(diǎn)的帶動(dòng)輻射作用，與周邊區(qū)域景點(diǎn)連接成旅游產(chǎn)品共同推向市場。

影響陜北地區(qū)傳統(tǒng)村落可達(dá)性的因素較多，需制定有層次的發(fā)展策略來改善該地區(qū)傳統(tǒng)村落可達(dá)性。坡度是影響陜南傳統(tǒng)村落可達(dá)性的主要因素，建議依據(jù)坡度的高低，差異化構(gòu)建傳統(tǒng)村落的重點(diǎn)保護(hù)系統(tǒng)，從而為有目標(biāo)、有層級(jí)地改善傳統(tǒng)村落的可達(dá)性夯實(shí)基礎(chǔ)。

本文還存在以下不足。在探究陜西省國家級(jí)傳統(tǒng)村落可達(dá)性時(shí)，將傳統(tǒng)村落抽象為點(diǎn)要素，忽視了傳統(tǒng)村落的內(nèi)在價(jià)值和外在吸引力。未來研究可以對(duì)陜西省國家級(jí)傳統(tǒng)村落的價(jià)值進(jìn)行梳理，從不同價(jià)值層面對(duì)傳統(tǒng)村落可達(dá)性及其影響因素進(jìn)行研究。

參考文獻(xiàn)

［1］ 中華人民共和國住房和城鄉(xiāng)建設(shè)部. 住房城鄉(xiāng)建設(shè)部文化部財(cái)政部關(guān)于公布第一批列入中國傳統(tǒng)村落名錄村落名單的通知［EB／OL］.（2012-12-20）［2015-06-19］.http：∥www.mohurd.gov.cn／wjfb／201212／t20121219-212340. html.

［2］ 陶偉， 陳紅葉， 林杰勇. 句法視角下廣州傳統(tǒng)村落空間形態(tài)及認(rèn)知研究［J］. 地理學(xué)報(bào)， 2013， 68（2）： 209-218.

TAO W， CHEN H Y， LIN J Y. Spatial form and spatial cognition of traditional village in syntactical view： A case study of Xiaozhou Village， Guangzhou［J］. Acta Geographica Sinica， 2013， 68（2）： 209-218.

［3］ 康璟瑤， 章錦河， 胡歡， 等. 中國傳統(tǒng)村落空間分布特征分析［J］. 地理科學(xué)進(jìn)展， 2016， 35（7）： 839-850.

KANG J Y， ZHANG J H， HU H， et al. Analysis on the spatial distribution characteristics of Chinese traditional villages［J］. Progress in Geography， 2016， 35（7）： 839-850.

［4］ 劉大均， 胡靜， 陳君子， 等. 中國傳統(tǒng)村落的空間分布格局研究［J］. 中國人口·資源與環(huán)境， 2014， 24（4）： 157-162.

LIU D J， HU J， CHEN J Z， et al. The study of spatial distribution of traditional villages in China［J］. China Population， Resources and Environment， 2014， 24（4）： 157-162.

［5］ 陳志文， 胡希軍， 葉向陽，等. 中國傳統(tǒng)村落有機(jī)體生長內(nèi)在邏輯研究［J］. 經(jīng)濟(jì)地理， 2020， 40（11）： 225-232.

CHEN Z W， HU X J， YE X Y， et al. Internal logic of organic growth in Chinese traditional villages［J］. Economic Geography， 2020， 40（11）： 225-232.

［6］ 楊麗婷， 曾禎. 古村落保護(hù)與開發(fā)綜合價(jià)值評(píng)價(jià)研究：以浙江省磐安縣為例［J］. 地域研究與開發(fā)， 2013， 32（4）： 112-116，122.

YANG L T， ZENG Z. Research on evaluation of ancient villages integrated value of the protection and development： A case study of Panan County of Zhejiang Province ［J］. Areal Research and Development， 2013， 32（4）： 112-116，122.

［7］ 汪清蓉， 李凡. 古村落綜合價(jià)值的定量評(píng)價(jià)方法及實(shí)證研究：以大旗頭古村為例［J］.旅游學(xué)刊， 2006， 21（1）： 19-24.

WANG Q R， LI F. On the quantitative evaluation method of ancient villages comprehensive value and demonstration study［J］. Tourism Tribune， 2006， 21（1）： 19-24.

［8］ 何艷冰， 張彤， 熊冬梅. 傳統(tǒng)村落文化價(jià)值評(píng)價(jià)及差異化振興路徑：以河南省焦作市為例［J］. 經(jīng)濟(jì)地理， 2020， 40（10）： 230-239.

HE Y B， ZHANG T， XIONG D M. Evaluation on cultural value of traditional villages and differential revitalization： A case study of Jiaozuo City， Henan Province［J］. Economic Geography， 2020， 40（10）： 230-239.

［9］ 康晨晨，黃曉燕，夏伊凡.傳統(tǒng)村落文化遺產(chǎn)價(jià)值分級(jí)分類評(píng)價(jià)體系構(gòu)建及實(shí)證：以陜西省國家級(jí)傳統(tǒng)村落為例［J］.陜西師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2023，51（2）：84-96.

KANG C C， HUANG X Y， XIA Y F. Construction and empirical evidence of graded and classified evaluati on system of tradition alvillage cultural herita gevalue： A case of national traditional villages in Shaanxi Province［J］.Journal of Shaanxi Normal University（Natural Science Edition），2023，51（2）：84-96.

［10］王云才，郭煥成，楊麗.北京市郊區(qū)傳統(tǒng)村落價(jià)值評(píng)價(jià)及可持續(xù)利用模式探討：以北京市門頭溝區(qū)傳統(tǒng)村落的調(diào)查研究為例［J］.地理科學(xué)， 2006， 26（6）： 735-742.

WANG Y C， GUO H C， YANG L. The discussion on evaluation and sustainable development models of traditional villages in the western mountainous area of Beijing： A case study of traditional villages in Mentougou District of Beijing［J］. Scientia Geographica Sinica， 2006， 26（6）： 735-742.

［11］劉春臘，劉沛林.北京山區(qū)溝域經(jīng)濟(jì)建設(shè)背景下的古村落保護(hù)與開發(fā)研究［J］.經(jīng)濟(jì)地理， 2011， 31（11）： 1923-1929.

LIU C L， LIU P L. Probe into protecting and developing ancient villages while constructing valley economy in the mountain area of Beijing［J］. Economic Geography， 2011， 31（11）： 1923-1929.

［12］祁嘉華，鄭曄梅.新農(nóng)村建設(shè)語境中的古村落保護(hù)與發(fā)展：以陜西為例［J］.西安建筑科技大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版）， 2011， 30（6）： 36-41.

QI J H， ZHENG Y M. The preservation and development of ancient villages in the new rural construction： Taking Shaanxi as the example［J］.Journal of Xian University of Architecture ＆ Technology （Social Science Edition）， 2011， 30（6）： 36-41.

［13］吳琳，余建忠.浙江傳統(tǒng)村落空間特征與保護(hù)利用探索［J］.城市發(fā)展研究， 2021， 28（3）： 32-39.

WU L， YU J Z. Exploration on the spatial characteristics and protection utilization of traditional villages in Zhejiang Province［J］. Urban Development Studies， 2021， 28（3）： 32-39.

［14］許建和，樂詠梅，毛洲，等.湖南省傳統(tǒng)村落空間格局影響因素與保護(hù)模式［J］.經(jīng)濟(jì)地理， 2020， 40（10）： 147-153.

XU J H， LE Y M， MAO Z， et al. Influencing factors and protection modes of the spatial pattern of traditional villages in Hunan Province［J］. Economic Geography， 2020， 40（10）： 147-153.

［15］邵秀英，劉亞玲，王向東，等.黃河流域傳統(tǒng)村落旅游響應(yīng)度及影響因素研究［J］.干旱區(qū)資源與環(huán)境， 2021， 35（6）： 200-208.

SHAO X Y， LIU Y L， WANG X D， et al. Tourism responsiveness of traditional villages in the Yellow River Basin［J］.Arid Land Resources and Environment， 2021， 35（6）： 200-208.

［16］陳剛.發(fā)展人類學(xué)視角下西部民族地區(qū)傳統(tǒng)村落旅游開發(fā)［J］.旅游學(xué)刊， 2017， 32（2）： 11-12.

CHEN G. Spatial form and spatial cognition of traditional village in syntactical view： A case study of Xiaozhou Village， Guangzhou［J］.Tourism Tribune， 2017， 32（2）： 11-12.

［17］李連璞.基于多維屬性整合的古村落旅游發(fā)展模式研究：以歷史文化名村為例［J］.人文地理， 2013， 28（4）： 155-160.

LI L P. Tourism sustainable development of villages based on multiattributes： A case study of national historical cultural village［J］. Human Geography， 2013， 28（4）： 155-160.

［18］黃卉潔，苗紅，王云.西北地區(qū)傳統(tǒng)村落旅游開發(fā)潛力評(píng)價(jià)研究［J］.地域研究與開發(fā)， 2021， 40（6）： 106-110，117.

HUANG H J， MIAO H， WANG Y. Research on the tourism development potential of traditional villages in Northwest China［J］. Areal Research and Development， 2021， 40（6）： 106-110，117.

［19］楊立國，劉沛林，林琳.傳統(tǒng)村落景觀基因在地方認(rèn)同建構(gòu)中的作用效應(yīng)：以侗族村寨為例［J］.地理科學(xué)， 2015， 35（5）： 593-598.

YANG L G， LIU P L， LIN L. The effect of traditional village landscape genes in construction place identity： Taking the Dong minority village as an example［J］. Scientia Geographica Sinica， 2015， 35（5）： 593-598.

［20］李伯華，石凱霞.傳統(tǒng)村落景觀基因識(shí)別及其特征分析：以張谷英村為例［J］.衡陽師范學(xué)院學(xué)報(bào)， 2018， 39（3）： 1-7.

LI B H， SHI K X. Landscape gene identification and its characteristic analysis of traditional villages： A case study of Zhangguiying Villages［J］. Journal of Hengyang Normal University， 2018， 39（3）： 1-7.

［21］陳娛，金鳳君，陸玉麒，等.京津冀地區(qū)陸路交通網(wǎng)絡(luò)發(fā)展過程及可達(dá)性演變特征［J］.地理學(xué)報(bào)， 2017， 72（12）： 2252-2264.

CHEN Y， JIN F J， LU Y Q， et al. Development history and accessibility evolution of land transportation network in Beijing-Tianjin-Hebei region［J］. Acta Geographica Sinica， 2017， 72（12）： 2252-2264.

［22］路海艷，趙鵬軍，董亞寧， 等.高鐵站可達(dá)性與城鄉(xiāng)居民收入差異空間分異研究［J］.地理科學(xué)進(jìn)展， 2022， 41（1）： 131-142.

LU H Y， ZHAO P J， DONG Y N， et al. Accessibility of high-speed rail stations and spatial disparity of urban-rural income gaps［J］. Progress in Geography， 2022， 41（1）： 131-142.

［23］陳潔， 陸鋒， 程昌秀. 可達(dá)性度量方法及應(yīng)用研究進(jìn)展評(píng)述［J］. 地理科學(xué)進(jìn)展， 2007， 26（5）： 100-110.

CHEN J， LU F， CHENG C X. Advance in accessibility evaluation approaches and applications［J］. Progress in Geography， 2007， 26（5）： 100-110.

［24］翟石艷，何新新，孔云峰， 等.基于最優(yōu)供需分配的公共設(shè)施空間可達(dá)性分析［J］.地理學(xué)報(bào)， 2022， 77（4）： 1028-1038.

ZHAI S Y， HE X X， KONG Y F， et al. Measuring the spatial accessibility of public services by optimal supply-demand allocation［J］. Acta Geographica Sinica， 2022， 77（4）： 1028-1038.

［25］仝德，孫裔煜，謝苗苗.基于改進(jìn)高斯兩步移動(dòng)搜索法的深圳市公園綠地可達(dá)性評(píng)價(jià)［J］.地理科學(xué)進(jìn)展， 2021， 40（7）： 1113-1126.

TONG D， SUN Y Y， XIE M M. Evaluation of green space accessibility based on improved Gaussian two-step floating catchment area method： A case study of Shenzhen City， China［J］. Progress in Geography， 2021， 40（7）： 1113-1126.

［26］張忠訓(xùn)，楊慶媛，王立，等.少數(shù)民族地區(qū)傳統(tǒng)村落交通可達(dá)性分析——以貴州省銅仁市為例［J］.資源科學(xué)， 2018， 40（11）： 2296-2306.

ZHANG Z X， YANG Q Y， WANG L， et al. Traffic accessibility analysis of traditional villages in minority areas： A case study of Tongren City， Guizhou Province［J］. Resources Science， 2018， 40（11）： 2296-2306.

［27］胡金龍，艾燁，鄭文俊，等.基于空間句法的桂林市傳統(tǒng)村落可達(dá)性研究［J］.西安建筑科技大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版）， 2022， 41（1）： 38-46.

HU J L， AI Y， ZHENG W J， et al. A Research on the accessibility of traditional villages in Guilin based on space syntax［J］. Journal of Xian University of Architecture ＆ Technology（Social Science Edition）， 2022， 41（1）： 38-46.

［28］竇銀娣，彭姍姍，李伯華，等.湖南省傳統(tǒng)村落空間可達(dá)性研究［J］.資源開發(fā)與市場， 2015， 31（5）： 554-558，583，641.

DOU Y D， PENG S S， LI B H， et al. Study on spatial accessibility of traditional village in Human Province［J］. Resource Development ＆ Market， 2015， 31（5）： 554-558，583，641.

［29］付強(qiáng)， 楊壯， 董鎖成， 等. 河南省國家級(jí)傳統(tǒng)村落空間可達(dá)性及影響因素研究［J］. 河南師范大學(xué)學(xué)報(bào)（自然科學(xué)版）. 2021， 49（6）： 82-90.

FU Q， YANG Z， DONG S C， et al. Research on spatial accessibility and influencing factors of national traditional villages in Henan Province［J］. Journal of Henan Normal University（Natural Science Edition）， 2021， 49（6）： 82-90.

［30］馬勇，黃智洵.基于GWR模型的長江中游城市群傳統(tǒng)村落空間格局及可達(dá)性探究［J］.人文地理， 2017， 32（4）：78-85.

MA Y， HUANG Z X. Research on spatial distribution and accessibility of the traditional village in the urban agglomeration on the middle reaches of the Yangtze River based on GWR model［J］. Human Geography， 2017， 32（4）： 78-85.

［31］FOTHERINGHAM A S，CHARLTON M E，BRUNSDON C. Geographically weighted regression： A natural evolution of the expansion method for spatial data analysis［J］. Environment and Planning A： Economy and Space， 1998， 30（11）： 1905-1927.

［32］BRUNSDON C， FOTHERINGHAM A S， CHARLTON M E. Geographically weighted regression： A method for exploring spatial nonstationarity［J］. Geographical Analysis， 1996， 28（4）： 281-298.

［33］MCMILLEN D P. Geographically weighted regression： The analysis of spatially varying relationships［J］. American Journal of Agricultural Economics， 2004， 86（2）： 554-556.