苑惠麗 楊亞楠 陳文君

摘要 基于多時相中分辨率遙感數(shù)據(jù)，結(jié)合野外實地采樣數(shù)據(jù)和面向?qū)ο蟮倪b感影像分析方法，開展黃河故道地區(qū)果園信息提取，結(jié)合地形因子、海岸線分布等基礎(chǔ)地理數(shù)據(jù)，探討果園空間分布格局，并分析不同空間尺度下果園的聚集度特征。結(jié)果表明：黃河故道果園面積約1 268.30 km，其中碭山縣、豐縣和蕭縣為果園種植的主要地區(qū)；果園主要分布于高程＞20～60 m、坡度≤4°的范圍內(nèi)的地勢平坦地帶，分別占果園總面積的86.84%和88.43%；果園主要分布于距離海岸線＞280 km的區(qū)域內(nèi)，占果園總面積的91.20%；縣（區(qū)）域、10 km×10 km、5 km×5 km和1 km×1 km格網(wǎng)尺度下的果園種植密度呈顯著空間正相關(guān)，屬于空間聚集分布模式，隨空間尺度降低，聚集度越高?？梢婞S河故道地區(qū)果園逐漸呈現(xiàn)出種植規(guī)?；?、連片化的特征。

關(guān)鍵詞 果園；空間分布特征；黃河故道地區(qū)

中圖分類號 S181.3 文獻標識碼 A

文章編號 0517-6611（2024）02-0063-06

doi：10.3969/j.issn.0517-6611.2024.02.013

開放科學(xué)（資源服務(wù)）標識碼（OSID）：

Spatial Distribution Characteristics of Orchards in the Ancient Canal of Yellow River

YUAN Hui-li，YANG Ya-nan，CHEN Wen-jun

（1.Architectural Engineering Institute，Jinling Institute of Technology，Nanjing，Jiangsu 211169; 2.State Key Laboratory of Lake Science and Environment，Nanjing Institute of Geography and Limnology，Chinese Academy of Science，Nanjing，Jiangsu 210008;3.Software Engineering Institute，Jinling Institute of Technology，Nanjing，Jiangsu 211169）

Abstract Based on multi temporal and medium resolution remote sensing data，combined with field sampling data and object-oriented remote sensing image analysis methods，information extraction of orchards in the ancient canal of Yellow River was carried out.Combining basic geographic data such as terrain factors and coastline distribution，the spatial distribution pattern of orchards was explored，and the clustering characteristics of orchards at different spatial scales were analyzed.The results showed that the area of orchards along the ancient canal of Yellow River was about 1 268.30 km.Dangshan County，F(xiàn)eng County and Xiao County were the major planting areas.Orchards were mainly distributed in flat terrain areas with an elevation of >20-60 m and a slope of ≤4°，accounting for 86.84% and 88.43% of the total orchard area，respectively;orchards were mainly distributed in areas more than 280 km away from the coastline，accounting for 91.20% of the total area of orchards.The planting density of orchards at the county （district） scale，10 km×10 km，5 km×5 km and 1 km×1 km grid scale showed a significant spatial positive correlation，belonging to a spatial aggregation distribution pattern.As the spatial scale decreased，the aggregation degree increased.It could be seen that orchards in the ancient canal Yellow River were gradually showing the characteristics of large-scale cultivation and contiguity.

Key words Orchard;Spatial distribution characteristic;Ancient canal of Yellow River

基金項目 國家自然科學(xué)基金項目（42101476，42201295）；江蘇省高等學(xué)校自然科學(xué)研究面上項目（21KJB170025）。

作者簡介 苑惠麗（1981—），女，山東濱州人，講師，博士，從事流域地表要素遙感反演研究。

收稿日期 2022-11-21；修回日期 2022-12-13

黃河故道又稱為故黃河、淤黃河，是歷史上黃河流域的重要組成部分，位于京津冀、黃河流域和長三角3個區(qū)域的連接地帶，地理位置重要，發(fā)展?jié)摿Υ蟆?019年，習(xí)近平總書記在黃河流域生態(tài)保護和高質(zhì)量發(fā)展座談會上的講話指出，黃河流域是打贏脫貧攻堅戰(zhàn)的重要區(qū)域；黃河流域生態(tài)環(huán)境持續(xù)明顯向好，下游河口濕地面積逐年回升，生物多樣性明顯增加，但流域生態(tài)環(huán)境脆弱。然而，由黃河水淤積形成的土壤有機質(zhì)含量低，農(nóng)業(yè)生產(chǎn)能力弱，沿線貧困化程度較高，但其獨特的氣候特征和砂質(zhì)土壤，適宜果樹生長。20世紀50年代開始建起了近百個國營園藝場，有力推動了黃河故道地區(qū)果園產(chǎn)業(yè)的快速發(fā)展。然而，經(jīng)過近20年的快速發(fā)展，在2010年左右果園的面積及產(chǎn)量增長速度變緩甚至有的區(qū)域開始減少種植面積，勢必會影響故道周邊的生態(tài)環(huán)境和經(jīng)濟發(fā)展。因此，有必要對黃河故道地區(qū)果園的種植現(xiàn)狀進行調(diào)查，獲取準確的果園種植情況，從而為改善生態(tài)環(huán)境和保證周邊經(jīng)濟的可持續(xù)發(fā)展提供數(shù)據(jù)支持。

黃河故道果園的研究多集中于果樹的栽培、病蟲害防治、土壤性狀及不同耕作和覆蓋條件下對果園產(chǎn)量的影響等方面。常規(guī)的統(tǒng)計方法很難及時準確地得到精確的農(nóng)作物種植面積，更無法實時獲取農(nóng)作物的空間分布狀況。隨著遙感技術(shù)的出現(xiàn)，快速準確地獲取果園的空間分布情況成為可能，尤其是大范圍果園信息提取。

因此，該研究擬基于多時相中分辨率遙感數(shù)據(jù)，結(jié)合野外實地采樣數(shù)據(jù)和面向?qū)ο蟮倪b感影像分析方法，開展黃河故道地區(qū)果園信息提取，結(jié)合地形因子、海岸線分布等基礎(chǔ)地理數(shù)據(jù)，探討果園空間分布格局，并分析不同空間尺度下果園的聚集度特征，為黃河故道地區(qū)果園產(chǎn)業(yè)化規(guī)模健康發(fā)展提供參考。

1 資料與方法

1.1 研究區(qū)概況

黃河故道自河南省蘭考縣向東南，途徑3個省的21個縣（區(qū)），自西向東依次經(jīng)過民權(quán)縣、商丘市區(qū)、安徽省碭山縣、江蘇省豐縣、徐州市區(qū)、睢寧縣、宿遷市區(qū)、泗陽縣、淮安市區(qū)，過漣水縣后，由濱?？h的大淤尖村套子口入黃海，全長約728.3 km（圖1）。除安徽蕭縣境內(nèi)、徐州市市轄區(qū)和徐州銅山區(qū)境內(nèi)有小面積山地海拔較高，區(qū)域內(nèi)總體地勢西高東低，平均海拔31.52 m，該區(qū)域?qū)贇夂蚪缬谂瘻貛Ш捅眮啛釒еg，屬于季風(fēng)半濕潤氣候區(qū)，冬秋無霜期長，光熱充足，東西跨度較大，由內(nèi)陸到沿海，氣候特征稍有差異，年均氣溫14.0～14.7 ℃，年無霜期213～220 d，年降雨量678～949 mm。其獨特的氣候特征和砂質(zhì)土壤，適宜果樹生長。

1.2 基礎(chǔ)地理數(shù)據(jù)

該研究使用的數(shù)據(jù)包括：①黃河故道地區(qū)30 m空間分辨率ASTER GDEM地形數(shù)據(jù)，下載于美國航天局（NASA）數(shù)據(jù)統(tǒng)一發(fā)布平臺Reverb （http：//reverb.echo.nasa.gov/reverb）；②Landsat 8 OLI（operational land image）L1級標準產(chǎn)品，該級別數(shù)據(jù)產(chǎn)品使用地面控制點和數(shù)字高程模型數(shù)據(jù)進行了幾何校正，坐標精度能滿足該研究要求。

經(jīng)過統(tǒng)計發(fā)現(xiàn)，覆蓋整個研究區(qū)需要5景Landsat 8 OLI影像，條帶號分別為120/36、120/37、121/36、122/36、123/36。據(jù)旱地作物的生長特征可知，10月中下旬農(nóng)作物已全部收獲，此時是提取果園的最佳時相。其中，位于條帶號為122/36的安徽省碭山縣梨園具有大面積連片分布的特征，且緊鄰農(nóng)村居民點，上覆有落葉林的農(nóng)村居民點與果園有相似的光譜特征和時相特征，利用IOPT指數(shù)可以將其有效區(qū)分，選取4月初、10月初及12月份影像進行該區(qū)域果園提取。因此，該研究中共選取出2017—2020年7景無云影像作為數(shù)據(jù)源（表1）。利用2020年4月獲取的野外調(diào)查樣點進行精度的驗證與評價。

1.3 基于NDVI的果園信息提取

黃河故道地區(qū)東西跨度較大，果園在內(nèi)陸和沿海地區(qū)有明顯不同的分布特點和影像特征。內(nèi)陸地區(qū)果園一般集中分布、面積大，與耕地鄰接，兩者在遙感影像上植被特征極其相似。但因其在生長期內(nèi)與耕地有明顯的差異，10月中下旬是耕地內(nèi)莊稼的收獲期，NDVI最低，而果園仍然茂盛，因此選取此期間獲取的影像可以有效區(qū)分耕地和果園。沿海地區(qū)多是小型采摘園，面積較小，中分辨率遙感影像上識別困難，僅利用光譜特征難以從遙感影像上將其提取出來，需結(jié)合Google Earth平臺的高分辨率衛(wèi)星遙感影像進行驗證并糾錯。

1.4 地形與區(qū)位特征

在同一空間參考系統(tǒng)下，將同一地區(qū)具有不同地理特征和屬性的數(shù)據(jù)進行疊加，能夠得到該地區(qū)的多重屬性特征，或建立地理對象之間的空間對應(yīng)關(guān)系。該研究利用空間疊加分析方法，對黃河故道地區(qū)果園的分布特征進行研究，包括：①將提取的果園數(shù)據(jù)與黃河故道地區(qū)縣、區(qū)分布圖進行疊加，統(tǒng)計得到各縣、區(qū)果園的種植面積；②將果園數(shù)據(jù)分別與高程和坡度2個地形因子進行疊加，統(tǒng)計分析果園在不同高程和坡度下的空間分布情況；③利用海岸線數(shù)據(jù)構(gòu)建多環(huán)緩沖區(qū)與果園數(shù)據(jù)進行疊加，統(tǒng)計分析果園在不同緩沖區(qū)范圍內(nèi)的分布特征。

1.5 基于Moran’s I指數(shù)的果園空間聚集度

相關(guān)位置上的地理數(shù)據(jù)由于空間相互作用和擴散的影響，彼此之間相互依賴，具有空間自相關(guān)性。定量化探究空間自相關(guān)性的方法按功能一般分為全局空間自相關(guān)和局部空間自相關(guān)。前者是從整個研究區(qū)域內(nèi)探測變量在空間分布上的聚集性，后者是從特定局部區(qū)域內(nèi)探測變量在空間分布上的聚集性。

該研究使用全局空間自相關(guān)指標Moran’s I指數(shù)描述黃河故道地區(qū)果園種植密度的空間相關(guān)性及其分布模式。Moran’s I指數(shù)的取值為[-1，1]，Moran’s I<0表示負相關(guān)，值越小，空間差異越大；Moran’s I=0表示不相關(guān)，空間呈隨機性；Moran’s I>0表示正相關(guān)，值越大，空間相關(guān)性越大。針對Moran’s I指數(shù)，使用標準化統(tǒng)計量Z來檢驗多個區(qū)域是否存在空間自相關(guān)關(guān)系。Moran’s I指數(shù)計算公式如下：

式中：I為Moran’s I指數(shù)；n為研究區(qū)的空間單元總數(shù)；x、x分別為i、j空間單元的屬性值；x為所有空間單元屬性值的均值；w為空間權(quán)重；S為空間權(quán)重矩陣之和。

利用Z檢驗進行驗證：

式中：E（I）為 Moran’s I指數(shù)的期望值；VAR（I）為 Moran’s I指數(shù)的方差。為研究不同空間尺度下黃河故道地區(qū)果園的空間聚集特征，構(gòu)建10 km×10 km、5 km×5 km、1 km×1 km 3種尺度的網(wǎng)格，計算包括縣（區(qū)）域的4種空間尺度下果園種植密度的全局自相關(guān)指數(shù)Moran’s I。該方法將網(wǎng)格圖與果園分布數(shù)據(jù)相疊加，統(tǒng)計每個網(wǎng)格的果園面積，計算得到每個網(wǎng)格果園面積百分比，將計算結(jié)果賦給對應(yīng)網(wǎng)格的中心點。

2 結(jié)果與分析

2.1 黃河故道地區(qū)果園種植規(guī)模

基于遙感影像提取果園結(jié)果（表2和圖2）表明：黃河故道地區(qū)果園種植總面積為1 268.30 km，其中，種植面積超過40 km的縣（區(qū)）有碭山縣、豐縣、蕭縣、虞城縣、寧陵縣和民權(quán)縣，占總種植面積的83.20%，其中碭山縣、豐縣和蕭縣種植面積較大，分別占總種植面積的44.26%、15.46%和8.31%。主要的果園類型包括蘋果園、梨園、桃園和葡萄園，其中梨園面積最大，為1 068.68 km，占果園總面積的84.26%；蘋果園面積為87.06 km，占果園總面積的6.86%；桃園面積為68.06 km，占果園總面積的5.37%；葡萄園面積為44.50 km，占果園總面積的3.51%。

2.2 果園分布的地形特征

地形為土地利用格局的形成提供基礎(chǔ)并影響著土地利用的演變過程。高程和坡度2個地形因子直接影響著農(nóng)林牧業(yè)的用地分布、水土流失等情況，是土地資源2個重要的環(huán)境因子。果園數(shù)據(jù)與地形因子疊加結(jié)果（表3）表明，黃河故道地區(qū)果園主要種植的高程為>20～60 m，占總種植面積的86.84%。其中，在>20～30 m高程帶內(nèi)，果園種植面積為193.53 km，占總種植面積的15.26%；>30～40 m高程帶內(nèi)，果園種植面積為428.98 km，占總種植面積的33.82%；>40～50 m高程帶內(nèi)，種植面積為329.22 km，占總種植面積的25.96%；>50～60 m高程帶內(nèi)，果園種植面積為149.68 km，占總種植面積的11.80%。如圖3所示，≤20 m高程范圍主要位于研究區(qū)東部，耕地類型主要為水田，幾乎沒有連片種植果園，果園種植多是采摘園形式，零星分布，面積較小，所占比重不大。>70 m高程范圍主要包含了民權(quán)縣部分園地、安徽蕭縣和徐州市市轄區(qū)范圍內(nèi)的部分山地果園。

黃河故道地區(qū)總體地勢平坦，統(tǒng)計表明（圖4），≤4°坡度范圍內(nèi)果園面積有1 121.56 km，占果園總面積的88.43%。

2.3 果園分布的區(qū)位特征

以70 km為間隔，對黃河故道地區(qū)海岸線構(gòu)建多環(huán)緩沖區(qū)，與果園數(shù)據(jù)疊加后統(tǒng)計結(jié)果如表4所示。分析表明，隨離海岸線距離的增加，果園種植面積逐漸增大，距離海岸線>280 km范圍內(nèi)果園種植面積為1 156.67 km，占果園總面積的91.20%。

2.4 果園分布的空間聚集特征

構(gòu)建10 km×10 km、5 km×5 km、1 km×1 km 3種不同的空間尺度格網(wǎng)，通過與黃河故道地區(qū)果園數(shù)據(jù)得到的不同空間尺度下的種植密度，分別計算包括縣（區(qū)）域的4種空間尺度下果園種植密度的全局自相關(guān)指數(shù)Moran’s I，得到研究區(qū)果園種植密度在不同空間尺度下的空間聚集特征（圖5）。統(tǒng)計表明，黃河故道地區(qū)縣（區(qū)）域、10 km×10 km、5 km×5 km、1 km×1 km 4種尺度下的Moran’s I指數(shù)分別為0.14、0.66、0.79、0.83，均大于0，Z檢驗值分別為2.39、24.61、56.84和300.41，均顯著。可見，在這4種空間尺度下，研究區(qū)內(nèi)的果園種植密度在整體空間分布上呈現(xiàn)出明顯的正相關(guān)關(guān)系，表現(xiàn)出顯著的空間聚集模式。通過比較發(fā)現(xiàn)，Moran’s I 指數(shù)隨著空間尺度的變小而逐漸增大，空間分布也更為聚集，相鄰空間單元的果園種植密度差異降低，相似性增大，表明果園種植密度在空間分布的聚集程度與空間尺度的大小有強烈的相關(guān)性。

3 結(jié)論

該研究利用中分辨率遙感影像，基于面向?qū)ο蟮倪b感影像分析方法，進行黃河故道地區(qū)果園信息提?。徊⒗肎IS空間分析方法，結(jié)合研究區(qū)地形因子與海岸線分布情況，對果園的空間分布格局進行了分析，結(jié)論如下：

（1）黃河故道地區(qū)果園種植的總面積約為1 268.30 km ，碭山縣、豐縣和蕭縣為該地區(qū)果園種植的主要區(qū)域。

（2）高程>20～60 m、坡度≤4°的區(qū)域是黃河故道地區(qū)果園主要分布地帶，具有典型的平原區(qū)果園種植特征。

（3）黃河故道地區(qū)果園主要分布于距離海岸線>280 m的內(nèi)陸區(qū)域，體現(xiàn)出內(nèi)陸自然條件依賴性。

（4）不同空間尺度下，黃河故道地區(qū)果園種植密度在整體空間分布上呈現(xiàn)出明顯的正相關(guān)關(guān)系，屬于空間聚集分布模式。隨著空間尺度的變小，種植密度的全局自相關(guān)性 Moran’s I 指數(shù)也逐漸增大，空間聚集程度逐漸增高，果園種植密度在空間內(nèi)分布的聚集程度與空間尺度的大小有強烈的相關(guān)性。

參考文獻

[1]趙瑜.讓黃河成為造福人民的幸福河：全國政協(xié)提案委員會“將黃河故道納入黃河流域生態(tài)保護和高質(zhì)量發(fā)展國家戰(zhàn)略”重點提案辦理協(xié)商會側(cè)記[J].中國政協(xié)，2020（22）：34-35.

[2]朱永增，丁繼輝，李倩倩，等.黃河故道土地生態(tài)整治與高標準農(nóng)田建設(shè)模式研究[J].江蘇水利，2020（8）：38-41.

[3]馬凌珂.黃河故道地區(qū)梨產(chǎn)業(yè)發(fā)展淺析[J].果農(nóng)之友，2020（11）：42-43.

[4]王寶秋.應(yīng)當開展果園面積變動統(tǒng)計[J].浙江統(tǒng)計，1999（1）：34.

[5]盛炳成.黃河故道地區(qū)果樹綜論[M].北京：中國農(nóng)業(yè)出版社，2008.

[6]李曉敏.河南省蘋果品種發(fā)展及栽培利用研究[D].楊凌：西北農(nóng)林科技大學(xué)，2018.

[7]邵礫群.中國蘋果矮化密植集約栽培模式技術(shù)經(jīng)濟評價研究[D].楊凌：西北農(nóng)林科技大學(xué)，2015.

[8]趙增鋒.蘋果病蟲害種類、地域分布及主要病蟲害發(fā)生趨勢研究[D].保定：河北農(nóng)業(yè)大學(xué)，2012.

[9]符丹丹.中國蘋果炭疽病病原菌的遺傳多樣性[D].楊凌：西北農(nóng)林科技大學(xué)，2014.

[10]簡興，陳鴻，肖昕迪.黃河故道土壤有機碳及易氧化有機碳含量分析[J].安徽農(nóng)學(xué)通報，2017，25（14）：62-64.

[11]高凝.不同地面覆蓋對黃河故道地區(qū)蘋果樹體生長及果實品質(zhì)的影響[D].泰安：山東農(nóng)業(yè)大學(xué)，2020.

[12]李芳，李姝蓉，章敏.顧及地形特征的山區(qū)農(nóng)作物種植面積測量與空間分布特征分析：以云南省文山州丘北縣為例[C]∥2016年全國統(tǒng)計建模大賽論文集.北京：國家統(tǒng)計局，2016：362-401.

[13]YUAN H L，MA R H，LUO J H.Mapping orchards on plain terrains using multi-temporal medium-resolution satellite imagery[J].Applied engineering in agriculture，2015，31（3）：351-362.

[14]徐晗澤宇，劉沖，齊述華，等.基于隨機森林算法的贛南柑橘果園遙感信息提取[J].江西師范大學(xué)學(xué)報（自然科學(xué)版），2018，42（4）：434-440.

[15]陳優(yōu)良，胡錦景，王兆茹，等.贛南柑橘果園面積動態(tài)變化分析：以信豐縣為例[J].江西農(nóng)業(yè)大學(xué)學(xué)報，2017，39（5）：884-894.

[16]陳優(yōu)良，胡錦景，王兆茹，等.贛南柑橘果園面積動態(tài)變化分析：以尋烏縣為例[J].江西理工大學(xué)學(xué)報，2017，38（5）：19-25.

[17]董芳.基于多源數(shù)據(jù)的丘陵區(qū)蘋果園地信息遙感提取技術(shù)研究[D].泰安：山東農(nóng)業(yè)大學(xué)，2012.

[18]RANCHIN T，NAERT B，ALBUISSON M，et al.An automatic method for vine detection in airborne imagery using wavelet transform and multiresolution analysis[J].Photogrammetric engineering and remote sensing，2001，67（1）：91-98.

[19]TRIAS-SANZ R.Texture orientation and period estimator for discriminating between forests，orchards，vineyards，and tilled fields[J].IEEE transactions on geoscience and remote sensing，2006，44（10）：2755-2760.

[20]喬斌，祝存兄，曹曉云，等.格網(wǎng)尺度下青海瑪多縣土地利用及生態(tài)系統(tǒng)服務(wù)價值空間自相關(guān)分析[J].應(yīng)用生態(tài)學(xué)報，2020，31（5）：1660-1672.

[21]劉洋，余建新，向冬蕾，等.基于局部空間自相關(guān)的思茅區(qū)耕地利用保護綜合分區(qū)[J].水土保持研究，2020，27（1）：183-188，196.

[22]張起明，徐晗澤宇，江豐，等.贛南地區(qū)柑橘果園空間分布特征分析[J].江西科學(xué)，2018，36（4）：591-598.

[23]徐建華.計量地理學(xué)[M].2版.北京：高等教育出版社，2014：84-91.

[24]陳運強，王榮遠.土地利用空間格局與坡度坡向的關(guān)系[J].測繪與空間地理信息，2015，38（10）：171-174，177.