陸丁滒， 吳虹， 郭琪， 陳夢杰

(1.桂林理工大學地球科學學院遙感應用研究所， 桂林　541004； 2.河南省地質(zhì)礦產(chǎn)勘查開發(fā)局第三地質(zhì)勘察院， 洛陽　471023)

基于GoogleEarth影像的漓江水系形態(tài)特征提取與分析

陸丁滒1， 吳虹1， 郭琪2， 陳夢杰1

(1.桂林理工大學地球科學學院遙感應用研究所， 桂林541004； 2.河南省地質(zhì)礦產(chǎn)勘查開發(fā)局第三地質(zhì)勘察院， 洛陽471023)

摘要：為獲取漓江流域水系的空間分布和水系形態(tài)特征，以Google Earth影像為數(shù)據(jù)源，采用人機交互目視解譯法提取了流域內(nèi)的水系形態(tài)特征信息，編制了漓江流域水系分布圖。通過對水系的等級、密度、分支比、各級水道長度以及交匯角等特征信息進行統(tǒng)計，定量化地展示了漓江流域水系的形態(tài)特征； 在此基礎上進行水系形態(tài)特征的控制因素分析，認為地層與斷裂構(gòu)造是控制漓江流域形態(tài)的重要因素，人類活動對漓江流域水系形態(tài)的影響正在日漸增強。該研究成果為漓江流域的綜合治理提供了科學依據(jù)，深化了漓江流域的地貌研究。

關鍵詞：Google Earth； 漓江流域； 水系特征； 信息提?。?控制因素

0引言

漓江是珠江水系在桂北地區(qū)的重要支流，流經(jīng)由加里東期花崗巖、中生代碎屑巖和碳酸鹽巖地層構(gòu)成的斷裂構(gòu)造發(fā)育區(qū)，形成了獨特的巖溶-河流組合地貌。21世紀之前，對漓江的研究主要集中于主干流的水資源與水環(huán)境方面，而對水系地貌特征的研究較少。根據(jù)流域地貌理論，水系地貌特征包括水道級別與流域級別、水道分支比、水道數(shù)量、水道平均長度、水道總長度、水道縱比降、水道交匯角、水道形態(tài)、水道流量和流域面積等28項指標[1]。漓江流域分布范圍廣，形態(tài)與地形復雜，交通不便，依靠地面調(diào)查方式去獲取這些地貌指標信息較為困難。高分辨率衛(wèi)星遙感以其廣視域、大信息量、精定位、同步和重復觀測的巨大優(yōu)勢，使開展這一研究工作成為可能。2000年以來，采用遙感手段研究河流地貌已得到廣泛應用。王磊等[2]采用GIS技術(shù)實現(xiàn)了從DEM獲取黃河源地區(qū)各級水系長度、數(shù)目、密度及流域面積等參數(shù)； 李棟梁[3]結(jié)合形態(tài)學理論，利用TM影像的多光譜信息提取了洛河流域主水系； 劉昱恒等[4]使用基于DEM數(shù)據(jù)的J&D算法及基于輔助數(shù)據(jù)的AGREE方法，完成了對潦河流域水系河網(wǎng)的提取。近年來，專門針對漓江水系地貌的研究也有了一定的進展，吳虹等[5]通過對3期Landsat5/7TM、ETM+遙感影像的信息提取、模式識別和實際調(diào)查，獲得了漓江流域的植被總量、城市化面積和水質(zhì)變化等信息； 秦潤君等[6]采用QuickBird-2和P6衛(wèi)星數(shù)據(jù)對漓江自然地貌破壞現(xiàn)狀進行了調(diào)查； 郝敏[7]基于資源三號衛(wèi)星(ZY-3)數(shù)據(jù)對漓江河床水深反演開展了研究。到目前為止，已開展的漓江水系地貌研究多從植被、土地利用等地理要素著手，水系形態(tài)特征分析較少，且未從整體上對漓江水系形態(tài)的地質(zhì)控制因素開展過專題研究。為此，本文以GoogleEarth影像為數(shù)據(jù)源，在MapGIS軟件的支持下，重點對漓江流域水系形態(tài)的水道級別、水系密度、水道分支比、各級水道長度和水道交匯角5個特征信息進行提取，分析水系空間形態(tài)展布的控制因素，為漓江的綜合治理提供科學依據(jù)。

1水系形態(tài)特征

水系作為宏觀的流域地貌組合[8]，一般由主干流和支流河道組成。本文主要對漓江流域的水道級別、水系密度、水道分支比、各級水道總長度和水道交匯角等5種特征指標進行研究和信息提取。

1)水道級別。水道級別是按水系中河道的承繼和匯流關系劃分的水道等級。對水道級別的劃分，目前國際上較受推崇的是Strahler水道級別劃分法，即將位于明顯谷地水流線以內(nèi)的所有間歇性和永久性水道稱為第1級水道，由2個1級水道匯流組成的新水道稱為第2級水道，匯流了2個2級水道的則稱之為第3級水道，依次類推，最終將整個流域的河網(wǎng)水系劃分完畢[9]。

2)水系密度。水系密度表示該地區(qū)水系分布的密集程度。在地學領域，計算河網(wǎng)水系密度有2種方法： 一種是通過計算單位面積內(nèi)河段的總長度來表示，即線密度； 另一種是通過計算單位面積內(nèi)河段的總面積來表示，即面密度[10]。線密度法的計算公式為

D=∑L/∑A，

(1)

式中： D為最終計算出的河網(wǎng)水系線密度； ∑L為一定面積內(nèi)河流的總長度； ∑A為總面積。

3)水道分支比。水道分支比是某一級水道的數(shù)目與比其高一級水道數(shù)目的比值，即

(2)

式中： r為分支比； nx為第x級水道的數(shù)目； nx+1為第x+1級水道的數(shù)目。

4)各級水道總長度。各級水道總長度即某級所有水道長度之和。在劃分了水道級別之后，各級水道總長度便可計算出來，其公式為

(3)

式中： Lu為第u級水道的總長度； n為第u級水道的數(shù)目，i=1,…,n； lui為第u級中第i條水道的長度。

5)水道交匯角。水道交匯角是指2條水道交匯時形成的角度大小。Horton[11]認為2條水道的交匯角與水道縱比降及地表坡度有以下關系

(4)

式中： Z為2條水道的交匯角； Sc為老水道的縱比降，(°)； Sg為新水道發(fā)育于其上地表的平均坡度。

2研究區(qū)概況

漓江發(fā)源于廣西壯族自治區(qū)興安縣華江鄉(xiāng)貓兒山東北面海拔1 732m(黃?；?的老山界南側(cè)，越城嶺主峰貓兒山東麓，由北向南流經(jīng)桂林市轄區(qū)及興安、靈川、臨桂、陽朔、平樂6個縣市，漓江水系由漓江干流與12條次級支流組成，屬復合型水系，干流全長214km。流域范圍為E109°45′18″～111°2′24″，N24°15′28″～25°54′42″，面積約10 733km2，流域內(nèi)多巖溶峰林和中低山丘陵地形(圖1)。

圖1　漓江流域干流及12條次級支流TM假彩色

3研究方案

以GoogleEarth平臺提供的衛(wèi)星遙感影像為數(shù)據(jù)源。具體研究方案流程如圖2所示。

圖2　研究方案流程

首先，采用MapGIS對獲取的影像進行預處理，主要包括影像裁剪、拼接、格式轉(zhuǎn)換、幾何糾正、投影變換和重采樣等； 其次，采用人機交互目視解譯矢量化法，提取漓江水系形態(tài)信息，根據(jù)Strahler水系分級原則，對提取的水道進行分級，在此基礎上開展統(tǒng)計分析，統(tǒng)計內(nèi)容包括各級水道數(shù)目、長度、分支比、交匯角和水系線密度等； 最后，結(jié)合研究區(qū)地質(zhì)圖，從巖性、斷裂構(gòu)造及人類活動的角度分析漓江水系形態(tài)的控制因素。

4漓江水系形態(tài)特征提取

本文使用的GoogleEarth影像分辨率從0.61m到100m不等，其中，桂林市區(qū)有QuickBird影像覆蓋，分辨率達0.61m，可以滿足本次研究的需要?；贕oogleEarth數(shù)據(jù)，采用目視解譯法進行漓江水系特征提取與分析。預處理步驟如下：

1)影像裁剪和拼接。利用GetScreen軟件進行，研究區(qū)長度設置為2.2km，控制單次圖片拼接數(shù)量在800張以內(nèi)。

2)幾何精糾正。為了使解譯結(jié)果更準確，以研究區(qū)1∶10萬地形圖為基準，利用MapGIS6.7軟件對所獲取的影像進行幾何糾正，并將其作為研究區(qū)工作底圖。

3)建立水道線文件。根據(jù)遙感解譯原則，建立解譯標志與水道信息分類編碼，矢量化水道。提取的漓江水系河網(wǎng)如圖3所示。

圖3　漓江流域水系形態(tài)信息提取結(jié)果

4.1水道分級信息提取

在解譯出漓江水系分布的基礎上，采用Strahler水道分級法對水道進行級別劃分。經(jīng)統(tǒng)計，漓江流域內(nèi)共有1級水道403條，2級水道75條，3級水道26條，4級水道4條，5級水道1條，其中，漓江水系干流為第5級水道。

4.2水系線密度信息提取

采用公式(1)對漓江水系進行線密度計算。在線密度計算時，用2km×2km的格網(wǎng)對河網(wǎng)矢量數(shù)據(jù)進行網(wǎng)格化，統(tǒng)計各單元網(wǎng)格內(nèi)水系的總長度，利用長度和網(wǎng)格面積計算出線密度值，最后將計算結(jié)果導入Surfer軟件，繪制出的水系密度圖見圖4。

圖4　漓江流域水系密度圖

通過圖4可以發(fā)現(xiàn)，漓江上游水系密度較大，密度值范圍為30～90km/km2，下游水系密度較小，密度值范圍為5～30km/km2。

4.3相鄰水道分支比信息提取

根據(jù)公式(2)，結(jié)合各級水道數(shù)目，計算出漓江水系相鄰2級水道的分支比(表1)。

表1　漓江水系水道分支比

從表1中可以看出，在漓江流域，各級水道之間的分支比有較大差異。這是由于流域內(nèi)不同的自然地理條件造成的。

4.4各級水道總長度信息提取

在MapGIS6.7 軟件屬性庫管理中可以查詢線段的長度，根據(jù)漓江水系分布和分級結(jié)果，經(jīng)統(tǒng)計計算漓江各級水道總長度見圖5。

圖5　漓江水系各級水道長度

從圖5可以看出，除第5級水道——漓江干流外，隨著水道級別的上升，漓江水系支流總長度呈減少的趨勢。

4.5各級水道交匯角信息提取

本研究將相應級別的水道交匯角度劃分為對應級別的交匯角，考慮到低級水道可直接匯入更高級的水道，故在劃分此類交匯角時，以低級水道對應的級別作為交匯角的級別。經(jīng)過統(tǒng)計，漓江水系各級水道平均交匯角見表2。

表2　漓江水系各級水道平均交匯角

從表2中可以看出，漓江水系各級水道呈近垂直狀交匯，這表明漓江各級水道的縱比降與地表平均坡度相差較大。

5控制因素分析

根據(jù)桂林市1∶20萬區(qū)域地質(zhì)圖，結(jié)合漓江水系密度特征，主要從地層、斷裂構(gòu)造以及人類活動的角度對控制漓江水系形態(tài)的因素開展分析。

5.1地層因素

第四紀以來，漓江流域地殼主要以上升為主，出露地層主要有寒武系、奧陶系、泥盆系、石炭系和第四系。筆者基于桂林幅1∶20萬區(qū)域地質(zhì)圖，提取出流域內(nèi)主要巖性信息，同時選取靈川、桂林、陽朔3個地面點作為控制點，并結(jié)合漓江水系密度圖將二者等比例疊放在一起，如圖6所示。

圖6　漓江流域水系密度和地層分布

從圖6可以看出，漓江流域內(nèi)的地層分布以大圩鎮(zhèn)為界，可分為北、南2個部分，大圩鎮(zhèn)以北地層主要走向大致為NE向，以南則以NW向為主，與漓江整體水系形態(tài)展布由下游的NW向到上游逐漸轉(zhuǎn)變?yōu)镹E向一致，說明本區(qū)的地層格局對漓江流域水系形態(tài)的展布與發(fā)展影響明顯。此外，在漓江流域北部，巖性變化由最初的加里東期花崗巖到砂、頁巖，到司門前—溶江段的頁巖、泥質(zhì)灰?guī)r和白云巖，再到溶江—靈川段的白云質(zhì)灰?guī)r、泥質(zhì)灰?guī)r、灰?guī)r和砂巖，具有巖石硬度逐漸降低，抗風化能力逐漸減弱的趨勢； 上游水系形態(tài)表現(xiàn)為徑流發(fā)育，支流繁多，巖石的透水性能差，容易被流水侵蝕，所以水系密度值大。流域南部的巖性為大圩—楊堤段的灰?guī)r、泥質(zhì)灰?guī)r，以及下游的楊堤—陽朔段的泥質(zhì)灰?guī)r、灰?guī)r，其巖性差異較小，故下游水系形態(tài)表現(xiàn)為地表徑流發(fā)育程度不高，地形坡度較小，巖石堅硬且透水性能好，水系長且疏，少有支流發(fā)育，所以水系密度值較小。

5.2斷裂構(gòu)造因素

在現(xiàn)代構(gòu)造運動明顯和活斷層附近地區(qū)，水系形態(tài)的發(fā)展主要受構(gòu)造活動的控制[12]。筆者將漓江流域地質(zhì)圖中的線性斷裂構(gòu)造信息單獨提取出來，以靈川、桂林、陽朔3個地面點作為控制點，等比例疊放在漓江流域水系線密度圖上，如圖7。

圖7　漓江流域水系密度和斷層構(gòu)造分布

從圖7可以看出，在老人山以北地區(qū)分布的斷裂自北向南，主要為NNE向的資源斷裂、SN向的龍勝—永福斷裂、NE向的桂林—來賓斷裂，整體上斷層方向有從NE向到NNE向再到SN向的逆時針轉(zhuǎn)動趨勢。漓江流域水系高密度帶呈NE向展布，并且自北向南由NE向向NNE向逐漸轉(zhuǎn)變，這與該地區(qū)的斷裂構(gòu)造方向一致。老人山以南地區(qū)，分布的斷裂帶有觀音閣斷裂和白石斷裂，整體由南及北呈NW向、NNW向再到SN向的順時針轉(zhuǎn)動趨勢。雖然該區(qū)水系密度較低，但其展布趨勢依然明顯，整體為由陽朔—遇龍河一帶的NW向陽朔—楊堤—奇峰鎮(zhèn)一帶的NNW方向轉(zhuǎn)動，這依然與該地區(qū)的斷裂構(gòu)造方向一致。

5.3人類活動因素

水系形態(tài)的展布除了受巖性、構(gòu)造等自然地質(zhì)因素的影響外，還在人類開發(fā)利用流域內(nèi)自然資源以及抗擊自然災害的過程中，不斷地被人為活動改造。早在2 200多a前，秦王朝出于統(tǒng)一中國政治和軍事行動的需要，鑿建了現(xiàn)今位于興安縣境內(nèi)的靈渠。改革開放近40a以來，隨著桂林城市化進程的加速及自然資源的過度開發(fā)利用，人類活動對漓江流域的主觀改造日益明顯，這在一定程度上都影響著漓江水系的自然形態(tài)。

根據(jù)上文解譯得到的水系形態(tài)信息，以司門前鎮(zhèn)為起點，每隔2km取1個干流水道寬度值，可得到漓江干流寬度變化趨勢，如圖8所示。

圖8　漓江干流寬度變化趨勢(作圖縱橫比例尺為1∶200)

從圖8可以看出，漓江上、中、下游的河流寬度值范圍差異較大。在漓江上游，由于匯入的支流不多，河流水量少，表現(xiàn)為主流寬度較窄，最窄僅40m左右。到了河流中游地區(qū)，由于有多條支流的匯入，河流的寬度明顯增大，并且由于存在規(guī)模比較大的江心洲，使河流的寬度范圍存在幾十m到200多m不等的變化。在距源頭60km附近，由于桂林市及靈川縣的城市用水，使河流水量急劇減少，河流寬度開始出現(xiàn)減小趨勢； 隨著距離河流源頭越來越遠，河流的動能減少，河流寬度在下游地區(qū)漸趨穩(wěn)定； 距源頭140km以后，由于流域水量的不斷匯聚，主干道又逐步變寬。

為了解決城市用水導致漓江水量急劇減少的問題，桂林市政府提出了枯季補水方案，先后分3期完成。第一期修建青獅潭水庫，第二期修建斧子口水庫、川江水庫和小溶江水庫，第三期則主要保持水面景觀。目前，第一期青獅潭水庫已建成并投入使用，其年內(nèi)調(diào)節(jié)能力對緩解漓江的枯水徑流問題具有積極的作用[13]。但是，由于水系形態(tài)的變化和人類活動都具有動態(tài)的不穩(wěn)定性，目前人類活動的影響主要表現(xiàn)在對河道的改造方面，它的影響既有積極的作用，又有消極的作用。因此，漓江流域水系形態(tài)的展布、發(fā)展與演變是一項需要長期觀察研究和對比分析的課題。

6結(jié)論

本文基于GoogleEarth影像，綜合運用遙感和GIS技術(shù)解譯提取了漓江流域的水系形態(tài)特征信息。

1)利用MapGIS軟件完成了對漓江流域水系形態(tài)信息的提取，獲得了漓江流域水系分布圖，可作為下一步研究的基礎數(shù)據(jù)。

2)通過對水系等級、支流數(shù)目、水系密度、分支比、交匯角和水道長度等信息的統(tǒng)計，認為漓江流域為5級水系，流域內(nèi)共有大小河道509條，流域北部水系密度較大，密度值范圍為30～90km/km2，南部水系密度較小，密度值范圍為5～30km/km2； 流域內(nèi)各級水道呈近90°交匯，整個水系總長度為1 946.63km。

3)通過分析漓江流域水系形態(tài)的控制因素，認為地層與斷裂構(gòu)造對其形態(tài)的控制較為明顯，人類活動對漓江流域自然地貌的影響也日益顯著，漓江流域的綜合治理已刻不容緩。

需要指出的是，自然形成的水系是由多種因素共同控制和影響的，本文僅從地層、斷裂構(gòu)造及人類活動的角度分析漓江水系形態(tài)的控制因素，分析結(jié)果有一定的局限性； 此外，GoogleEarth影像自身存在分辨率的差異，在一定程度上也影響了解譯的精度，這些問題有待于在今后的研究中進一步完善和解決。

參考文獻(References)：

[1]承繼成,江美球.流域地貌數(shù)學模型[M].北京:科學出版社,1986.

ChengJC,JiangMQ.TheMathematicalModelofWatershedGeomorphology[M].Beijing:SciencePress,1986.

[2]王磊,胡伍生,吳波.基于GIS的水系提取方法與分析[J].現(xiàn)代測繪,2004,27(6):45-46.

WangL,HuWS,WuB.BaseonGISwatersystemdistillmethodandanalyse[J].ModernSurveyingandMapping,2004,27(6):45-46.

[3]李棟梁.基于TM影像的水系信息提取及變化制圖研究[D].南京:河海大學,2006.

LiDL.DrainageSystemBasedonTMImageInformationExtractionandChangeCartographyResearch[D].Nanjing:HohaiUniversity,2006.

[4]劉昱恒,徐宏根.兩種不同DEM處理方法提取的水系比較[J].資源調(diào)查與環(huán)境,2012,33(3):206-210.

LiuYH,XuHG.ComparisonsofwatersystemsextractedfromDEMwithtwodifferentprocessingmethods[J].ResourcesSurveyandEnvironment,2012,33(3):206-210.

[5]吳虹,郭遠飛,郭建東,等.基于TM/ETM+和MODIS的漓江流域生態(tài)環(huán)境遙感動態(tài)監(jiān)測[J].國土資源遙感,2007,19(3):56-61.doi:10.6046/gtzyyg.2007.03.13.

WuH,GuoYF,GuoJD,etal.Long-termandshort-termdynamicmonitoringofecologicalenvironmentvariationinthedrainageareaoftheLijiangRiverbasedonTM/ETM+andMODIS[J].RemoteSensingforLandandResources,2007,19(3):56-61.doi:10.6046/gtzyyg.2007.03.13.

[6]秦潤君,吳虹,郭琪,等.基于遙感和GIS技術(shù)的漓江自然地貌破壞現(xiàn)狀調(diào)查[J].國土資源遙感,2013,25(1):160-164.doi:10.6046/gtzyyg.2013.01.28.

QinRJ,WuH,GuoQ,etal.InvestigationofdamagesituationofthenaturallandformalongLijiangRiverbasedonGISandRS[J].RemoteSensingforLandandResources,2013,25(1):160-164.doi:10.6046/gtzyyg.2013.01.28.

[7]郝敏.基于TM&QB-2數(shù)據(jù)的漓江河床水深遙感反演研究[D].桂林:桂林理工大學,2014.

HaoM.RemoteSensingInversionStudyofLijiangRiverWaterDepthBasedonTMandQB-2[D].Guilin:GuilinUniversityofTechnology,2014.

[8]曹伯勛.地貌學及第四紀地質(zhì)學[M].武漢:中國地質(zhì)大學出版社,1995.

CaoBX.GeomorphologyandQuaternaryGeology[M].Wuhan:ChinaUniversityofGeosciencesPress,1995.

[9]StrahlerAN.Quantitativeanalysisofwatershedgeomorphology[J].TransactionsAmericanGeophysicalUnion,1957,38(6):913-920.

[10]王金艷,黃永林,譚慧明,等.利用ArcGIS生成水系密度的方法[J].地理空間信息,2010,8(6):101-102.

WangJY,HuangYL,TanHM,etal.NewmethodfordrawingriverdensitymapwithArcGIS[J].GeospatialInformation,2010,8(6):101-102.

[11]HortonRE.Erosionaldevelopmentofstreamsandtheirdrainagebasins:Hydrophysicalapproachtoquantitativemorphology[J].GeologicalSocietyofAmericaBulletin,1945,56(3):275-370.

[12]丁國瑜.活動走滑斷裂帶的斷錯水系與地震[J].地震,1982(1):3-8.

DingGY.Activestrike-slipfaultzoneofthedislocationdrainageandearthquake[J].Earthquake,1982(1):3-8.

[13]蔣亞萍.漓江枯水、功能性斷流與水利工程的關系[J].桂林工學院學報,2005,25(4):426-431.

JiangYP.LowWaterandfunctionalcutoffandreservoirsofLijiangRiver[J].JournalofGuilinUniversityofTechnology,2005,25(4):426-431.

(責任編輯： 陳理)

FeatureextractionandanalysisoftheLijiangRiverwatersystemformbasedontheGoogleEarthimage

LUDingge1,WUHong1,GUOQi2,CHENMengjie1

(1. Guilin University of Technology Remote Sensing Institute， Guilin 541004， China; 2. The Third Geological Prospecting Institute of Geology and Mineral Resources Exploration Development Authorities， Luoyang 471023， China)

Abstract:In order to obtain the river system space distribution and river system feature information of the Lijiang River Basin, the authors, with Google Earth images as information source and by using the method of man-machine interactive visual interpretation, extracted basin river system configuration information and compiled water distribution map of Lijiang River basin. On the basis of information extraction and statistics of such factors as the water level, drainage density, stream tributaries branching ratio, length, and intersection angle, the morphological characteristics of Lijiang River Basin system was quantitatively demonstrated. An analysis of control factors based on river system morphology characteristics revealed that the strata and faults are the important controlling factors of Lijiang River morphology, and that the influence of human activity is growing. The research results provide objective scientific basis for the Lijiang River comprehensive control and treatment and also fill the blank in the study of the Lijiang River basin landform.

Keywords:Google Earth； Lijiang River basin; drainage characteristics; information extraction; control factor

doi:10.6046/gtzyyg.2016.02.25

收稿日期：2014-10-23；

修訂日期：2014-12-11

基金項目：國家科技支撐計劃項目“漓江流域遙感動態(tài)監(jiān)測與應用示范關鍵技術(shù)”(編號： 2012BAC16B01-2)和廣西科技廳科技攻關項目“漓江流域生態(tài)環(huán)境保護與可持續(xù)發(fā)展研究”(編號： 桂科攻1298006-1)共同資助。

中圖法分類號：TP 79

文獻標志碼：A

文章編號：1001-070X(2016)02-0161-07

第一作者簡介：陸丁滒(1986-)，男，碩士研究生。主要從事遙感技術(shù)與應用方面的研究。Email： ldg8677@163.com。

通信作者：吳虹(1947-)，教授，博士研究生導師。Email： wuhong@glut.edu.cn。

引用格式： 陸丁滒,吳虹,郭琪,等.基于GoogleEarth影像的漓江水系形態(tài)特征提取與分析[J].國土資源遙感,2016,28(2):161-167.(LuDG,WuH,GuoQ,etal.FeatureextractionandanalysisoftheLijiangRiverwatersystemformbasedontheGoogleEarthimage[J].RemoteSensingforLandandResources,2016,28(2):161-167.)