• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于小型無人機(jī)與MODIS數(shù)據(jù)的草地植被覆蓋度研究
    ——以甘南州為例

    2017-03-02 08:29:05宋清潔張瑤瑤孟寶平高金龍向宇軒
    草業(yè)科學(xué) 2017年1期
    關(guān)鍵詞:甘南州蓋度植被指數(shù)

    宋清潔,崔 霞,張瑤瑤,孟寶平,高金龍,向宇軒

    (1.蘭州大學(xué)資源環(huán)境學(xué)院西部環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,甘肅 蘭州 730000; 2.草地農(nóng)業(yè)生態(tài)系統(tǒng)國(guó)家重點(diǎn)實(shí)驗(yàn)室,蘭州大學(xué)草地農(nóng)業(yè)科技學(xué)院,甘肅 蘭州 730020; 3.蘭州大學(xué)資源環(huán)境學(xué)院,甘肅 蘭州 730000)

    基于小型無人機(jī)與MODIS數(shù)據(jù)的草地植被覆蓋度研究
    ——以甘南州為例

    宋清潔1,崔 霞1,張瑤瑤1,孟寶平2,高金龍2,向宇軒3

    (1.蘭州大學(xué)資源環(huán)境學(xué)院西部環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,甘肅 蘭州 730000; 2.草地農(nóng)業(yè)生態(tài)系統(tǒng)國(guó)家重點(diǎn)實(shí)驗(yàn)室,蘭州大學(xué)草地農(nóng)業(yè)科技學(xué)院,甘肅 蘭州 730020; 3.蘭州大學(xué)資源環(huán)境學(xué)院,甘肅 蘭州 730000)

    以甘南州為研究區(qū),采用Canon數(shù)碼相機(jī)和小型無人機(jī)搭載相機(jī)在不同大小樣方上獲取草地植被覆蓋度數(shù)碼照片,結(jié)合2015年5月-10月的Terra/MODIS植被指數(shù)產(chǎn)品MOD13Q1,分析了增強(qiáng)型植被指數(shù)(EVI)和歸一化植被指數(shù)(NDVI)與草地植被覆蓋度之間的相關(guān)性,建立了研究區(qū)草地植被覆蓋度的回歸模型,并對(duì)模型進(jìn)行了精度評(píng)價(jià),篩選出甘南州草地植被覆蓋度最優(yōu)遙感反演模型,并對(duì)草地生長(zhǎng)季時(shí)期覆蓋度時(shí)空上的動(dòng)態(tài)特征進(jìn)行分析。結(jié)果表明,1)利用小型無人機(jī)搭載相機(jī)獲取草地大樣方植被數(shù)碼照片的方法能用于地面草地覆蓋度數(shù)據(jù)的獲??;2)與NDVI相比,用EVI估算草地覆蓋度更優(yōu),因此確定基于EVI構(gòu)建的對(duì)數(shù)模型為甘南州草地植被覆蓋度最優(yōu)反演模型,模型精度可達(dá)88.00%;3)研究區(qū)2015年生長(zhǎng)季草地植被覆蓋度除了低平地草甸在8月達(dá)到最大值外,其它草地類型均在7月達(dá)到最大值;4)甘南州以中高植被覆蓋度為主,主要分布在瑪曲、碌曲、夏河以及合作四縣市。整體而言,中西部和西南部區(qū)域草地覆蓋度高于東部。通過精確草地植被覆蓋度模型的建立,不僅有利于及時(shí)準(zhǔn)確的了解草地植被覆蓋度的時(shí)空分布狀況和季節(jié)性動(dòng)態(tài)變化,也有利于維護(hù)甘南州草地生態(tài)系統(tǒng)的持續(xù)穩(wěn)定發(fā)展。

    Terra/MODIS;NDVI;EVI;相關(guān)性分析;遙感反演模型;精度評(píng)價(jià);草地覆蓋度動(dòng)態(tài)分析

    草地是陸地生態(tài)系統(tǒng)的重要組成部分[1],具有水源涵養(yǎng)、水土保持和氣候調(diào)節(jié)等重要功能,也是畜牧業(yè)的經(jīng)濟(jì)支柱[2]。描述草地的生物物理參數(shù)主要有植被類型、蓋度、生物量、反照率、粗糙度、葉面積指數(shù)等[3-4],其中植被覆蓋度是指樣地中地上部分的全部植物個(gè)體(包括葉、莖、枝等)的垂直投影面積占樣地總面積的百分比[5-6],是表征生態(tài)系統(tǒng)植被群落生長(zhǎng)狀況及生態(tài)環(huán)境質(zhì)量的重要參數(shù),其值大小能直接反映出一個(gè)地區(qū)的生態(tài)狀況[7]。所以,精準(zhǔn)草地蓋度遙感模型的建立,不但能及時(shí)準(zhǔn)確地掌握草地植被覆蓋度的時(shí)空分布狀況和季節(jié)性動(dòng)態(tài)變化規(guī)律,而且有利于維護(hù)草地生態(tài)系統(tǒng)的持續(xù)穩(wěn)定發(fā)展。

    目前,對(duì)于草地植被覆蓋度的估測(cè)主要有兩種方法[8]:地面測(cè)量法和遙感反演法。地面測(cè)量法主要包括拍照法、采樣法、儀器法和目估法等[9]。地面觀測(cè)方法往往受到不同研究區(qū)空間異質(zhì)性的影響,若擴(kuò)展到整個(gè)區(qū)域有一定的難度,并且受人為因素影響大,對(duì)于草地植被覆蓋度的大范圍精確估測(cè)較難??墒谴朔ㄓ休^高精度的樣方尺度,因此經(jīng)常將此法用來作為遙感反演的地面驗(yàn)證數(shù)據(jù)。遙感技術(shù)具有監(jiān)測(cè)范圍廣、空間連續(xù)、周期短、受地面條件限制少、信息量大、資源消耗少等特點(diǎn),有利于快速獲取大區(qū)域草地植被覆蓋度信息[10]。在草地植被覆蓋度的遙感研究方面,早期多用NOAA/AVHRR(advanced very high resolution radiometer)遙感資料進(jìn)行研究[11-13],隨著1999年MODIS(moderate-resolution imaging spectroradiometer)衛(wèi)星的發(fā)射成功,MODIS數(shù)據(jù)在土地覆蓋、冰雪變化、監(jiān)測(cè)荒漠化程度以及植被蓋度等領(lǐng)域應(yīng)用廣泛[14-20],利用MODIS植被指數(shù)數(shù)據(jù)估算草地植被覆蓋度也成為研究熱點(diǎn)[21-23]。植被指數(shù)(vegetation index,VI)是由遙感衛(wèi)星探測(cè)數(shù)據(jù)的線性和非線性組合而構(gòu)成的可以反映綠色植被生長(zhǎng)狀況和分布特征的指數(shù)[24],可以很好地反映蓋度和生物量、綠色植被生長(zhǎng)狀況及空間分布格局等特征。其中增強(qiáng)型植被指數(shù)(enhanced vegetation index,EVI)和歸一化植被指數(shù)(normalized difference vegetation index,NDVI)[25]已廣泛應(yīng)用于草地植被覆蓋度及草地地上生物量、草地退化狀況分析、草地利用強(qiáng)度分析等領(lǐng)域[26-28]。

    在草地植被覆蓋度野外調(diào)查中,地表樣方大小設(shè)置與對(duì)應(yīng)的衛(wèi)星影像分辨率大小之間的空間匹配問題一直是遙感反演中的難題。傳統(tǒng)的利用照相法獲取地表樣方植被覆蓋度基本都在較小樣方上進(jìn)行[7,9,29-30],但對(duì)于中分辨率遙感數(shù)據(jù)來說,這種地面測(cè)量樣方是否能很好地與遙感產(chǎn)品象元匹配,測(cè)量樣方是否能體現(xiàn)對(duì)應(yīng)遙感象元的地物狀況還有待商榷。

    針對(duì)以上問題,本研究利用小型無人機(jī)(大疆精靈Phantom 3 Professional:1/2.3英寸CMOS影像傳感器,有效像素1 240 萬;照片最大分辨率4 000×3 000)及Canon數(shù)碼相機(jī)(EOS 5D Mark Ⅱ:2110萬有效像素;最高分辨率5 616×3 744)獲取草地植被覆蓋度大樣方的小樣方實(shí)測(cè)數(shù)據(jù),結(jié)合高時(shí)間分辨率的Terra/MODIS的植被指數(shù)產(chǎn)品MOD13Q1,旨在建立基于大樣方實(shí)測(cè)數(shù)據(jù)和MODIS植被指數(shù)數(shù)據(jù)的甘南草地植被覆蓋度反演模型,并分析甘南州2015年草地生長(zhǎng)季(5月-10月)期間覆蓋度時(shí)間與空間上的動(dòng)態(tài)變化特征,為甘南州草地資源的可持續(xù)利用和動(dòng)態(tài)監(jiān)測(cè)提供科學(xué)合理的依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    甘南藏族自治州位于甘肅省西南,處于黃土高原的西部和青藏高原東北邊緣的過渡地段(33°06′-35°44′ N,100°46′-104°44′ E)[29,31-33],境內(nèi)地形復(fù)雜,山地、高原占土地總面積的98%以上,從西北到東南地勢(shì)逐漸降低,平均海拔3 000 m。甘南州屬于高原大陸性氣候,年降水量在400~800 mm;年均溫1~3 ℃;全年平均日照時(shí)數(shù)2 200~2 400 h,自西北向東南逐漸降低[31]。甘南州草原牧草豐茂,是甘肅省重要的畜牧業(yè)生產(chǎn)基地,被譽(yù)為“亞洲最優(yōu)良的牧場(chǎng)”之一,不僅具有很高的經(jīng)濟(jì)價(jià)值,而且有利于維護(hù)長(zhǎng)江和黃河源頭地區(qū)的生態(tài)安全。甘南州草地面積為2.602 5×106hm2,其占甘南州土地總面積的67.64%,其中可以利用的草地面積大概為2.495×106hm2,占草地總面積的95.86%;耕地面積為7.66×104hm2,占甘南州土地總面積的2.01%[34]。高寒草甸、高寒灌叢草甸、溫性草甸草原、溫性草原、沼澤、暖性草叢、低平地草甸地7種是甘南州的主要草地類型[31,35-36],瑪曲、碌曲、夏河、合作4個(gè)縣市為甘南州的主要牧區(qū)。近幾十年來,在自然和人為兩種因素共同作用下,我國(guó)草原生態(tài)環(huán)境問題仍較突出[37],甘南州草地植被出現(xiàn)退化狀況[29,38],在很大程度上阻礙了甘南州的進(jìn)一步發(fā)展,與此同時(shí)也給當(dāng)?shù)氐淖匀簧鷳B(tài)和社會(huì)經(jīng)濟(jì)系統(tǒng)的可持續(xù)性發(fā)展帶來了不利影響[29,38-39]。

    1.2 數(shù)據(jù)獲取

    1.2.1 地面實(shí)測(cè)數(shù)據(jù)的獲取 選擇草地植被空間分布上較為一致、能夠代表較大區(qū)域植被覆蓋狀況的地區(qū)設(shè)置30 m×30 m的樣地,每塊樣地設(shè)置兩種不同尺度的樣方:小樣方5個(gè)(1 m×1 m)[9]和大樣方1個(gè)(30 m×30 m)。對(duì)于小樣方,用Canon數(shù)碼相機(jī)垂直對(duì)樣方進(jìn)行拍攝,獲取覆蓋小樣方的草地植被照片;對(duì)于大樣方,用小型無人機(jī)在樣地中心上空一定高度垂直對(duì)地拍攝,獲取覆蓋大樣地的草地植被照片。每塊樣地具體調(diào)查內(nèi)容包括草地植被覆蓋度、生物量、牧草高度、土壤質(zhì)地等指標(biāo),并使用GPS記錄采樣點(diǎn)的經(jīng)緯度及海拔信息。Canon數(shù)碼相機(jī)和小型無人機(jī)所拍攝的植被照片根據(jù)任世龍等[40]對(duì)植被照片的處理方法進(jìn)行處理,采用三人處理結(jié)果的均值作為該樣方的草地植被覆蓋度實(shí)測(cè)數(shù)據(jù)。甘南州牧區(qū)主要集中在瑪曲、碌曲、夏河和合作4個(gè)縣市,迭部、卓尼、臨潭和舟曲主要為林地和農(nóng)業(yè)用地,故樣地主要分布在瑪曲、碌曲、夏河和合作4個(gè)縣市。2015年7月-9月草地生長(zhǎng)季期間,在研究區(qū)共計(jì)調(diào)查樣地86塊,包含430個(gè)小樣方和86個(gè)大樣方(圖1)。

    1.2.2 遙感數(shù)據(jù)的獲取 本研究所用遙感數(shù)據(jù)為MOD13Q1產(chǎn)品,屬于NASA MODIS陸地三級(jí)標(biāo)準(zhǔn)數(shù)據(jù)產(chǎn)品,包括NDVI和EVI,空間分辨率250 m。下載的MOD13Q1數(shù)據(jù)時(shí)間序列為2015年生長(zhǎng)季(5月-10月)。此數(shù)據(jù)可在NASA對(duì)地觀測(cè)系統(tǒng)數(shù)據(jù)共享平臺(tái)上(EOS data gateway)下載(http://delenn.gsfc.nasa.gov/-imswww/pub/imswelcome/)[29,32]。

    圖1 研究區(qū)及采樣點(diǎn)分布圖Fig.1 Study area and distribution of sampling points

    1.2.3 其它數(shù)據(jù) 甘南州行政邊界圖、甘南州草地類型圖等相關(guān)空間數(shù)據(jù)庫(kù),數(shù)據(jù)來源于草地農(nóng)業(yè)科技學(xué)院草地農(nóng)業(yè)生態(tài)系統(tǒng)國(guó)家重點(diǎn)實(shí)驗(yàn)室。

    1.3 研究方法

    1.3.1 草地植被覆蓋度信息獲取 參考任世龍等[40]對(duì)植被照片的處理方法對(duì)本研究中Canon數(shù)碼相機(jī)和小型無人機(jī)所拍攝的植被照片進(jìn)行處理。為了減少草地植被覆蓋度計(jì)算中人為因素產(chǎn)生的誤差,用三人處理結(jié)果的均值作為該樣方的草地植被覆蓋度數(shù)據(jù)。對(duì)于設(shè)置的5個(gè)小樣方,用其均值作為該樣地草地植被覆蓋度驗(yàn)證數(shù)據(jù)。由于小型無人機(jī)拍攝照片覆蓋面積較大,有部分無人機(jī)照片的草地植被覆蓋度無法識(shí)別,因而剔除相應(yīng)樣地小樣方獲取的草地植被覆蓋度數(shù)據(jù),研究區(qū)共計(jì)調(diào)查樣地73塊,包含365個(gè)小樣方和73個(gè)大樣方。

    1.3.2 衛(wèi)星遙感數(shù)據(jù)的處理 運(yùn)用MRT(modis reprojection tool)進(jìn)行數(shù)據(jù)讀取[25],同時(shí)轉(zhuǎn)換遙感影像的數(shù)據(jù)格式和地圖投影;然后運(yùn)用ArcGIS軟件掩膜提取研究區(qū)數(shù)據(jù)。

    在使用MOD13Q1數(shù)據(jù)之前將每個(gè)象元值縮小10 000倍,以得到真實(shí)的EVI、NDVI值,并將16天植被指數(shù)數(shù)據(jù)采用最大值合成法(maximum value composite,MVC)[41]生成5月-10月生長(zhǎng)季的月植被指數(shù)及2015年年最大植被指數(shù)。

    1.3.3 不同樣方大小獲取草地植被覆蓋度結(jié)果相關(guān)性分析 為了檢驗(yàn)利用小型無人機(jī)對(duì)大樣方(30 m×30 m)獲取地面草地植被覆蓋度方法的準(zhǔn)確性,對(duì)大樣方的草地植被覆蓋度與同一樣地所設(shè)置的5個(gè)小樣方(1 m×1 m)的覆蓋度均值之間進(jìn)行相關(guān)性分析。

    1.3.4 草地植被覆蓋度反演模型的構(gòu)建及精度評(píng)價(jià) 草地遙感植被覆蓋度研究中用到的模型一般有線性和非線性兩種。但是由于不同研究區(qū)水熱條件、植被特征、土壤等條件不同,最終確定采用的模型也不同[42]。

    采用統(tǒng)計(jì)分析軟件的回歸分析方法,分別以每塊樣地的EVI、NDVI兩種植被指數(shù)作為自變量,對(duì)應(yīng)的小型無人機(jī)獲取的草地植被覆蓋度數(shù)據(jù)作為因變量,建立草地植被覆蓋度與EVI、NDVI兩種植被指數(shù)之間的回歸模型。

    均方根誤差(root mean square error,RMSE)[25,29,32]常用來衡量模型的估測(cè)能力、描述模型的不確定性以及量化模型精度[43],RMSE數(shù)值越低,表明回歸模型精確度越高[44]。以Canon數(shù)碼相機(jī)獲取的樣地草地植被覆蓋度為真實(shí)值,對(duì)所建立的回歸模型進(jìn)行精度評(píng)價(jià),分別計(jì)算各反演模型估測(cè)的草地植被覆蓋度與真實(shí)值之間的平均相對(duì)誤差(P)(式1)、均方根誤差(RMSE)(式2),分析各回歸模型的草地植被覆蓋度估測(cè)精度(Accuracy)(式3),篩選出研究區(qū)草地植被覆蓋度監(jiān)測(cè)的最佳遙感反演模型。

    (1)

    (2)

    Accuracy=1-P/100

    (3)

    1.3.5 草地植被覆蓋度空間反演 利用最終確定的研究區(qū)遙感監(jiān)測(cè)模型及2015年5月-10月生長(zhǎng)季月植被指數(shù)及2015年年最大植被指數(shù)分析甘南州不同草地類型植被覆蓋度月季變化特征及全州2015年草地植被覆蓋度空間分布特征。

    1.3.6 草地植被覆蓋度分級(jí) 利用2015年MODIS年最大植被指數(shù)數(shù)據(jù)及甘南州草地植被覆蓋度最優(yōu)反演模型計(jì)算得出甘南州2015年草地植被覆蓋度年最大數(shù)據(jù)。根據(jù)甘南州2015年生長(zhǎng)季草地植被覆蓋度[32]、國(guó)家《土地利用現(xiàn)狀調(diào)查技術(shù)規(guī)程》、全國(guó)《草場(chǎng)資源調(diào)查技術(shù)規(guī)程》以及前人關(guān)于植被劃分的標(biāo)準(zhǔn)[45-47],同時(shí)結(jié)合甘南州植被特有的生態(tài)特征,將研究區(qū)的草地植被覆蓋度分為五級(jí)。一級(jí)植被:植被覆蓋度大于80%,屬高覆蓋度植被,相當(dāng)于高產(chǎn)草地、密灌木地;二級(jí)植被:植被覆蓋度為60%~80%,屬中高覆蓋度植被,相當(dāng)于中高產(chǎn)草地;三級(jí)植被:植被覆蓋度為40%~60%,屬中覆蓋度植被,相當(dāng)于中產(chǎn)草地;四級(jí)植被:植被覆蓋度為20%~40%,屬低覆蓋度植被,相當(dāng)于中低產(chǎn)草地;五級(jí)植被:植被覆蓋度小于20%,屬低覆蓋度植被,相當(dāng)于強(qiáng)度、中度沙漠化土地,低產(chǎn)草地。

    2 結(jié)果與分析

    2.1 不同樣方大小獲取的草地植被覆蓋度相關(guān)性分析

    利用小型無人機(jī)對(duì)大樣方(30m×30m)獲取地面草地植被覆蓋度與同一樣地所設(shè)置的5個(gè)小樣方(1m×1m)草地植被覆蓋度(Canon數(shù)碼相機(jī)獲取的數(shù)據(jù))的平均值進(jìn)行相關(guān)性分析(圖2)。結(jié)果表明,兩種不同樣方大小獲取的草地植被覆蓋度之間具有很高的相關(guān)性(R2=0.963,P<0.01)。說明利用小型無人機(jī)對(duì)地面大樣方拍攝的方法可以用于地面草地植被覆蓋度數(shù)據(jù)的獲取。

    2.2 草地植被覆蓋度反演模型構(gòu)建及精度評(píng)價(jià)

    通過統(tǒng)計(jì)分析軟件分析出小型無人機(jī)獲取的地面草地植被覆蓋度與MODIS兩種植被指數(shù)(EVI,NDVI)之間的線性、對(duì)數(shù)、乘冪、增長(zhǎng)及指數(shù)回歸模型(表1)可以看出,草地植被覆蓋度與EVI、NDVI之間均具有顯著的相關(guān)性。對(duì)于EVI而言,對(duì)數(shù)模型(R2=0.707)優(yōu)于其它4種模型;對(duì)于NDVI而言,線性模型(R2=0.694)優(yōu)于其它4種模型。從精度評(píng)價(jià)結(jié)果(表2)可以看出,EVI對(duì)數(shù)模型的精度最高,達(dá)88.00%,平均相對(duì)誤差(P)和均方根誤差(RMSE)最低。由此可見,基于MODIS-EVI對(duì)數(shù)模型反演的甘南州草地植被覆蓋度與真實(shí)值最接近,能較好地模擬生長(zhǎng)季期間草地生長(zhǎng)狀態(tài)和覆蓋度動(dòng)態(tài)變化。因此,基于MODIS-EVI的對(duì)數(shù)模型為甘南州草地植被覆蓋度的最優(yōu)模型:

    圖2 Canon數(shù)碼相機(jī)和小型無人機(jī)獲取草地 植被覆蓋度相關(guān)性分析Fig. 2 Correlation analysis between the grassland fractional vegetation cover (FVC) as measured by Canon digital camera and small unmanned aerial vehicles (UAV)

    y=65.317lnx+111.964R2=0.707

    (4)

    式中:y為草地植被覆蓋度預(yù)測(cè)值(%);x為MODIS-EVI。

    2.3 草地植被覆蓋度月度變化動(dòng)態(tài)

    從甘南州2015年草地植被生長(zhǎng)季植被覆蓋度變化動(dòng)態(tài)分析可以看出,草地植被覆蓋度在7月達(dá)到最大值(66.29%),最低值(14.21%)出現(xiàn)在5月。在生長(zhǎng)季期間,草地植被覆蓋度呈先增長(zhǎng)再降低的單峰拋物線型(圖3)。草地于5月進(jìn)入返青期且快速生長(zhǎng);6月份的草地蓋度較5月份增長(zhǎng)了大約34.07%;6月到7月草地蓋度繼續(xù)上升;7月-8月為草地生長(zhǎng)的旺盛期,在此時(shí)草地蓋度達(dá)到了最大值,8月以后氣候漸冷,草地生長(zhǎng)逐漸停止,因此9月的覆蓋度較8月下降了22.98%,到10月時(shí)草地植被覆蓋度下降了20.26%。

    從甘南州不同草地類型植被覆蓋度的月變化動(dòng)態(tài)可以看出,低平地草甸8月達(dá)到最大值,其它類型草地均在7月達(dá)到最大值(圖4)。雖然沼澤類的覆蓋度在7種草地類型中是最大的,但是沼澤在甘南州占的比例很少,所以對(duì)整個(gè)甘南州草地植被覆蓋度影響不大。高寒草甸和高寒灌叢草甸是甘南州最主要的草地類型,植被覆蓋度也相對(duì)較高,在很大程度上影響著甘南州草地植被覆蓋度的動(dòng)態(tài)變化。

    2.4 草地植被覆蓋度空間反演及分級(jí)

    對(duì)甘南州2015年草地植被覆蓋度進(jìn)行分級(jí),得到草地植被覆蓋度空間分布圖(圖5)。從中可以看出,甘南州中西部和西南部地區(qū)的草地植被覆蓋度較東部高,以中高植被覆蓋度為主,主要分布在瑪曲、碌曲、夏河以及合作4個(gè)縣(市)(表3),草地植被覆蓋度的等級(jí)都為二級(jí),其它等級(jí)的草地植被覆蓋度所占比例相對(duì)較少。

    表1 甘南州不同草地植被指數(shù)與草地植被覆蓋度回歸模型比較Table 1 Comparison of regression models between fractional vegetation cover (%) and vegetation index in Gannan Prefecture (N=73)

    注:** 表示極顯著(P<0.01);EVI為增強(qiáng)型植被指數(shù),NDVI為歸一型植被指數(shù);表2同。

    Note:** means significant at the 0.01 level;EVI, enhanced vegetation index, NDVI, normalized difference vegetation index; Similarly for Table 2.

    表2 甘南州草地植被指數(shù)與植被覆蓋度回歸模型精度評(píng)價(jià)(N=73)Table 2 Accuracy assessment of regression models between fractional vegetation cover and vegetation index in Gannan Prefecture

    圖3 甘南州2015年生長(zhǎng)季草地植被覆蓋度月度變化動(dòng)態(tài)Fig.3 Monthly variation of grassland fractional vegetation cover in 2015 growing season in Gannan Prefecture

    3 討論與結(jié)論

    3.1 討論

    本研究中草地蓋度的調(diào)查方法采用小型無人機(jī)所拍攝的照片計(jì)算而得到,相較于傳統(tǒng)地面樣方尺度蓋度數(shù)據(jù),其主要特點(diǎn)包括:1)覆蓋范圍大(大小約900 m2),具有更好的空間代表性,樣方覆蓋范圍是傳統(tǒng)調(diào)查樣方的900倍,能夠更好地反映草地地表蓋度情況,傳統(tǒng)的點(diǎn)測(cè)量數(shù)據(jù)缺乏中尺度甚至更大尺度的代表性,如何取得具有代表性的數(shù)據(jù),也就是獲得中尺度數(shù)據(jù)的真實(shí)值,是一個(gè)非常棘手的問題[48],已有的對(duì)于大中尺度植被蓋度的研究,通常借助較高分辨率遙感數(shù)據(jù)(如TM、ASTER)改進(jìn)較低分辨率遙感數(shù)據(jù)(如MODIS)對(duì)植被覆蓋度監(jiān)測(cè)的精度[49-50],但較高分辨率遙感數(shù)據(jù)的可獲得性往往受到限制,采用無人機(jī)獲取地面高分辨率蓋度圖像受天氣干擾較小,可獲得性高;2)操作快捷方便,效率高。在傳統(tǒng)地表樣方布置方法中,受樣方大小的限制,往往通過增加樣地內(nèi)樣方數(shù)目和改進(jìn)采樣布置方案來增加樣地的空間代表性,在提高草地蓋度反演模型精度的同時(shí)增加了工作量。而利用小型無人機(jī)攝相技術(shù)估測(cè)草地蓋度大量減少地表采樣工作(由原來的3~5個(gè)樣方減少為1個(gè)樣方)[9,28],無需進(jìn)入草地便可獲取草地蓋度信息,節(jié)省了大量時(shí)間和人力。3)估測(cè)精度較傳統(tǒng)方法高,本研究中使用小型無人機(jī)照片所計(jì)算的蓋度與傳統(tǒng)方法所計(jì)算的蓋度之間有較高的相關(guān)性,其相關(guān)系數(shù)高達(dá)0.963。與NDVI、EVI所構(gòu)建的草地蓋度反演模型精度較高,精度高達(dá)81.99%和88.00%。

    圖4 甘南州2015年生長(zhǎng)季各類草地月最大植被覆蓋度動(dòng)態(tài)變化Fig. 4 Maximum values of monthly fractional vegetation cover for different grassland types in 2015 growing season in Gannan

    圖5 甘南州2015年生長(zhǎng)季最大草地植被 覆蓋度分級(jí)空間分布圖Fig.5 The largest classification of grassland fractional vegetation cover spatial distribution in 2015 growing season表3 甘南州2015年生長(zhǎng)季主要牧業(yè)縣(市) 草地植被覆蓋度及等級(jí)Table 3 Value and level of grassland fractional vegetation cover in 2015 growing season in the main pastoral county of Gannan Prefecture

    縣(市)County(City)面積Area/km2草地植被覆蓋度均值TheaverageofgrasslandFVC/%植被覆蓋度等級(jí)ThelevelofgrasslandFVC瑪曲Maqu1019074.85二級(jí)2nd碌曲Luqu529874.84二級(jí)2nd夏河Xiahe627470.37二級(jí)2nd合作Hezuo267070.70二級(jí)2nd

    近年來發(fā)展無人機(jī)為遙感平臺(tái)的近地(50~100 m)遙感信息獲取技術(shù)已成為構(gòu)建星空一體化遙感系統(tǒng)的理想銜接點(diǎn)[51],同時(shí)隨著數(shù)碼相機(jī)的普及,攝影技術(shù)已成為一種炙手可熱的數(shù)字圖像信息獲取手段,基于數(shù)字圖像技術(shù)的植物生物量、葉面積指數(shù)(LAI)估測(cè)研究也得到了長(zhǎng)足發(fā)展[52-54],本研究使用的小型旋翼無人機(jī)自帶的相機(jī)不是專門為量測(cè)目的設(shè)計(jì)的,其邊緣畸變不能被忽視,在今后利用無人機(jī)進(jìn)行植被監(jiān)測(cè)的研究中需要充分認(rèn)識(shí)圖像邊緣畸變對(duì)地面實(shí)測(cè)值的影響。

    3.2 結(jié)論

    植被覆蓋度是衡量植被狀況的重要定量指標(biāo),本研究利用小型無人機(jī)獲取的地面草地植被覆蓋度數(shù)據(jù)結(jié)合高分辨率的MODIS植被指數(shù)數(shù)據(jù)(EVI,NDVI)確定甘南州草地覆蓋度最優(yōu)反演模型,與此同時(shí)分析了2015年生長(zhǎng)季草地植被覆蓋度的時(shí)空變化特征,得出如下主要結(jié)論:

    1)利用小型無人機(jī)獲取的地面大樣方草地植被覆蓋度數(shù)據(jù)與傳統(tǒng)方法獲取的小樣方數(shù)據(jù)之間具有很高的相關(guān)性,說明小型無人機(jī)可以用于獲取地面草地植被覆蓋度數(shù)據(jù),小型無人機(jī)的使用增大了監(jiān)測(cè)樣方的面積,不僅準(zhǔn)確快捷而且使樣方設(shè)置更具代表性,為快速獲取地面草地植被狀況提供了新的途徑。

    2)MODIS產(chǎn)品的植被指數(shù)數(shù)據(jù)(EVI/NDVI)與甘南州草地植被覆蓋度具有較強(qiáng)的相關(guān)性(表1),說明MODIS產(chǎn)品在研究區(qū)具有很好的適用性,可以及時(shí)有效地監(jiān)測(cè)甘南州草地動(dòng)態(tài)變化,對(duì)促進(jìn)當(dāng)?shù)啬翗I(yè)發(fā)展具有重大意義。與NDVI相比,EVI可以更好地反映研究區(qū)草地植被覆蓋度的時(shí)空差異,因此確定基于EVI構(gòu)建的草地植被覆蓋度最優(yōu)反演模型為y=65.317lnx+111.964,該模型可以解釋研究區(qū)草地植被覆蓋度變化的70.7%,精度可達(dá)88.00%,能很好地模擬草地生長(zhǎng)季期間(5月-10月)的草地植被覆蓋度的動(dòng)態(tài)變化。

    3)研究區(qū)2015年生長(zhǎng)季草地植被覆蓋度在7月達(dá)到最大值(66.29%),在5月出現(xiàn)最低值(14.21%)。在生長(zhǎng)季期間,草地植被覆蓋度具有先上升后下降的特點(diǎn),其變化趨勢(shì)呈單峰拋物線型,草地植被覆蓋度年最大值主要集中在7月-8月,除了低平地草甸在8月達(dá)到最大值外,其它草地類型均在7月達(dá)到最大值。

    4)甘南州2015年生長(zhǎng)季草地植被覆蓋度狀況較好,中西部和西南部地區(qū)的草地植被覆蓋度較東部高,以中高植被覆蓋度為主(表3),主要分布在瑪曲、碌曲、夏河以及合作4個(gè)牧業(yè)縣(市),其它等級(jí)的草地植被覆蓋度所占比例相對(duì)較少,說明甘南州牧業(yè)縣(市)草地狀況良好。使用遙感技術(shù)進(jìn)行大范圍草地蓋度分析可以為草地資源管理與合理規(guī)劃放牧活動(dòng)提供指導(dǎo)意見,甘南州主要牧業(yè)縣(市)可實(shí)行劃區(qū)輪牧,根據(jù)草地植被生長(zhǎng)狀況,合理配置草地。

    References:

    [1] 趙同謙,歐陽志云,賈良青,鄭華.中國(guó)草地生態(tài)系統(tǒng)服務(wù)功能間接價(jià)值評(píng)價(jià).生態(tài)學(xué)報(bào),2004,24(6):1101-1110. Zhao T Q,Ouyang Z Y,Jia L Q,Zheng H.Ecosystem services and their valuation of China grassland.Acta Ecologica Sinica,2004,24(6):1101-1110.(in Chinese)

    [2] 查勇,Gao J,倪紹祥.國(guó)際草地資源遙感研究新進(jìn)展.地理科學(xué)進(jìn)展,2003,22(6):607-617. Zha Y,Gao J,Ni S X.Most recent progress of international research on remote sensing of grassland resources.Progress in Geography,2003,22(6):607-617.(in Chinese)

    [3] 王寶元.保護(hù)草地生態(tài)環(huán)境走可持續(xù)發(fā)展之路.四川草原,2003(6):1-3. Wang B Y.Protection of the eco-environment on grasslands and the paths to sustainable development.Sichuan Caoyuan,2003(6):1-3.(in Chinese)

    [4] 高飛,邢文淵,李大平,趙睿,師慶東.和靜草地覆蓋度變化遙感監(jiān)測(cè)及分析.草業(yè)科學(xué),2007,24(4):27-30. Gao F,Xing W Y,Li D P,Zhao R,Shi Q D.Remote sensing monitoring and analysis of grassland coverage in Hejing.Pratacultural Science,2007,24(4):27-30.(in Chinese)

    [5] 溫慶可,張?jiān)鱿?劉斌,喬竹萍.草地覆蓋度測(cè)算方法研究進(jìn)展.草業(yè)科學(xué),2009,26(12):30-36. Wen Q K,Zhang Z X,Liu B,Qiao Z P.Research progress in grassland fractional coverage estimation methods.Pratacultural Science,2009,26(12):30-36.(in Chinese)

    [6] Gitelson A A,Kaufman Y J,Stark R,Rundquist D.Novel algorithms for remote estimation of vegetation fraction.Remote Sensing of Environment,2002,80:76-87.

    [7] 陳建軍,宜樹華,任世龍,秦彧,王曉云.疏勒河上游高寒草地植被蓋度反演及精度評(píng)價(jià).草業(yè)科學(xué),2014,31(1):56-65. Chen J J,Yi S H,Ren S L,Qin Y,Wang X Y.Retrieval of fractional vegetation cover of alpline grassland and the efficiency of remote sensing retrieval in the upper of Shule River Basin.Pratacultural Science,2014,31(1):56-65.(in Chinese)

    [8] 顧祝軍,曾志遠(yuǎn).遙感植被蓋度研究.水土保持研究,2005,12(2):18-21. Gu Z J,Zeng Z Y.Overview of researches on vegetation coverage in remote sensing.Research of Soil and Water Conservation,2005,12(2):18-21.(in Chinese)

    [9] 宜樹華,王增如,謝霞,楊時(shí)海,黃磊,葉柏生.高寒草地植被蓋度估算及其與凍土的關(guān)系.草業(yè)科學(xué),2011,28(3):353-358. Yi S H,Wang Z R,Xie X,Yang S H,Huang L,Ye B S.Estimation of fractional vegetation cover and its relation with permafrost in the upstream regions of Shule River Basin.Pratacultural Science,2011,28(3):353-358.(in Chinese)

    [10] 郭述茂.基于3S技術(shù)的高寒草地植被覆蓋度分布特征及動(dòng)態(tài)變化研究.蘭州:蘭州大學(xué)碩士學(xué)位論文,2009. Guo S M.Research of vegetation coverage distribution changes of alpine grassland based on 3S technology.Master Thesis.Lanzhou:Lanzhou University,2009.(in Chinese)

    [11] Purevdorj T S,Tateishi R,Ishiyama T,Honda Y.Relationships between percent vegetation cover and vegetation indices.International Journal of Remote Sensing,1998,19(18):3519-3535.

    [12] Zeng X B,Dickinson R E,Walker A,Shaikh M,Defries R S,Qi J G.Derivation and evaluation of global 1-km fractional vegetation cover data for land modeling.Journal of Applied Meteorology,2000,39(6):826-839.

    [13] Wittich K P,Hansing O.Area-averaged vegetative cover fraction estimated from satellite data.International Journal of Biometeorology,1995,38:209-215.

    [14] 趙彩霞,施昆,寧平.EOS-MODIS在環(huán)境科學(xué)中的應(yīng)用與研究進(jìn)展.環(huán)境科學(xué)導(dǎo)刊,2008,27(2):15-20. Zhao C X,Shi K,Ning P.Research progress of application of EOS-MODIS on environmental science.Environmental Science Survey,2008,27(2):15-20.(in Chinese)

    [15] 劉闖,文洪濤,趙立成,張瑋.我國(guó)EOS-MODIS地面站建設(shè)的現(xiàn)狀、問題與對(duì)策.遙感信息,2003(4):42-47. Liu C,Wen H T,Zhao L C,Zhang W.Retrieving foliar biochemical components with statistical and physical models.Remote Sensing Information,2003(4):42-47.(in Chinese)

    [16] 王瑋,馮琦勝,郭鈮,沙莎,胡蝶,王麗娟,李耀輝.基于長(zhǎng)時(shí)間序列NDVI資料的我國(guó)西北干旱區(qū)植被覆蓋動(dòng)態(tài)監(jiān)測(cè).草業(yè)科學(xué),2015,32(12):1969-1979. Wang W,Feng Q S,Guo N,Sha S,Hu D,Wang L J,Li Y H.Dynamic monitoring of vegetation coverage based on long time-series NDVI data sets in northwest arid region of China.Pratacultural Science,2015,32(12):1969-1979.(in Chinese)

    [17] 嚴(yán)建武,李春娥,袁雷,陳全功.EOS-MODIS數(shù)據(jù)在草地資源監(jiān)測(cè)中的應(yīng)用進(jìn)展綜述.草業(yè)科學(xué),2008,25(4):1-9. Yan J W,Li C E,Yuan L,Chen Q G.Application summary of EOS-MODIS data in the monitoring of grassland resources.Pratacultural Science,2008,25(4):1-9.(in Chinese)

    [18] Maffei C,Leone A P,Vella M,Guiseppe M,Menenti M.The use of MODIS-simulated spectral bands for monitoring plant water stress as a help for dynamic fire risk assessment.Remote Sensing for Agriculture,Ecosystems,and Hydrology Ⅷ,2006 doi:10.1117/12.689905.

    [19] 李聰,肖繼東,宮恒瑞,石玉,馮志敏.MODIS數(shù)據(jù)在烏魯木齊地區(qū)植被景觀動(dòng)態(tài)監(jiān)測(cè)中的應(yīng)用.新疆氣象,2005,28(3):27-29. Li C,Xiao J D,Gong H R,Shi Y,Feng Z M.The application of MODIS data to plant view dynamic monitoring in Urumqi.Bimonthly of Xinjiang Meteorology,2005,28(3):27-29.(in Chinese)

    [20] 宮攀,陳仲新.基于MODIS數(shù)據(jù)的東北地區(qū)植被物候參數(shù)提取.土壤通報(bào),2009,40(2):213-217. Gong P,Chen Z X.Regional vegetation phenology monitoring based on MODIS.Chinese Journal of Soil Science,2009,40(2):213-217.(in Chinese)

    [21] 王春香,張滌非,任萬輝.MODIS數(shù)據(jù)植被覆蓋度提取算法比較.大氣與環(huán)境光學(xué)學(xué)報(bào),2010,5(6):457-462. Wang C X,Zhang D F,Ren W H.Comparison of vegetation coverage extracting based on MODIS data.Journal of Atmospheric and Environmental Optics,2010,5(6):457-462.(in Chinese)

    [22] Lehnert L W,Meyer H,Wang Y,Miehe G,Thies B.Retrieval of grassland plant coverage on the Tibetan Plateau based on a multi-scale,multi-sensor and multi-method approach.Remote Sensing of Environment,2015,164:197-207.

    [23] 劉春靜,張麗,周宇,張炳華,侯小麗.中國(guó)新疆及中亞五國(guó)干旱區(qū)草地覆蓋度反演與分析.草業(yè)科學(xué),2016,33(5):861-870. Liu C J,Zhang L,Zhou Y,Zhang B H,Hou X L.Retrieval and analysis of grassland coverage in arid Xinjiang,China and five countries of Centrial Asia.Pratacultural Science,2016,33(5):861-870.(in Chinese)

    [24] 程乾,黃敬峰,王人潮.MODIS和NOAA/AVHRR植被指數(shù)差異初步分析.科技通報(bào),2005,21(2):205-209. Cheng Q,Huang J F,Wang R C.Preliminary analysis of difference between MODIS and NOAA/AVHRR vegetation indices.Bulletin of Science and Technology,2005,21(2):205-209.(in Chinese)

    [25] 方金.基于多源遙感數(shù)據(jù)的甘南牧區(qū)草地生物量遙感監(jiān)測(cè)研究.蘭州:蘭州大學(xué)碩士學(xué)位論文,2013. Fang J.Monitoring grassland biomass based on multi-source remote sensing data in Gannan pastoral area.Master Thesis.Lanzhou:Lanzhou University,2013.(in Chinese)

    [26] 郭芬芬,范建容,邊金虎,劉飛,張懷珍.基于MODIS NDVI時(shí)間序列數(shù)據(jù)的藏北草地類型識(shí)別.遙感技術(shù)與應(yīng)用,2011,26(6):821-826. Guo F F,Fan J R,Bian J H,Liu F,Zhang H Z.Grassland types identification based on time-series MODIS NDVI data in northern Tibet.Remote Sensing Technology and Application,2011,26(6):821-826.(in Chinese)

    [27] 范瑛,李小雁,李廣泳.基于MODIS/EVI的內(nèi)蒙古高原西部植被變化.中國(guó)沙漠,2014,34(6):1671-1677. Fan Y,Li X Y,Li G Y.Spatio-temporal vegetation variation in the western Inner Mongolia Plateau based on MODIS/EVI.Journal of Desert Research,2014,34(6):1671-1677.(in Chinese)

    [28] 除多,次仁多吉,王彩云,姬秋梅,德央.利用MODIS數(shù)據(jù)估算西藏高原地表植被覆蓋度.遙感技術(shù)與應(yīng)用,2010,25(5):707-713. Chu D,Cirenduoji,Wang C Y,Ji Q M,De Y.Vegetation coverage estimation over Tibetan Plateau using MODIS data.Remote Sensing Technology and Application,2010,25(5):707-713.(in Chinese)

    [29] 馬琳雅,崔霞,馮琦勝,梁天剛.2001-2011年甘南草地植被覆蓋度動(dòng)態(tài)變化分析.草業(yè)學(xué)報(bào),2014,23(4):1-9. Ma L Y,Cui X,Feng Q S,Liang T G.Dynamic changes of grassland vegetation coverage from 2001 to 2011 in Gannan Prefecture.Acta Prataculturae Sinica,2014,23(4):1-9.(in Chinese)

    [30] Curran P J,Williamson H D.Sample size for ground and remotely sensed data.Remote Sensing of Environment,1986,20(1):31-41.

    [31] 梁天剛,崔霞,馮琦勝,王鶯,夏文韜.2001-2008年甘南牧區(qū)草地地上生物量與載畜量遙感動(dòng)態(tài)監(jiān)測(cè).草業(yè)學(xué)報(bào),2009,18(6):12-22. Liang T G,Cui X,Feng Q S,Wang Y,Xia W T.Remotely sensed dynamics monitoring of grassland aboveground biomass and carrying capacity during 2001-2008 in Gannan pastoral area.Acta Prataculturae Sinica,2009,18(6):12-22.(in Chinese)

    [32] 馬琳雅.甘南州草地植被覆蓋度與物候期時(shí)空變化動(dòng)態(tài)特征.蘭州:蘭州大學(xué)碩士學(xué)位論文,2013. Ma L Y.Spatio-temporal dynamic changes of grassland vegetation cover and phenology in Gannan Prefecture.Master Thesis.Lanzhou:Lanzhou University,2013.(in Chinese)

    [33] 郭正剛,高新華,劉興元,梁天剛.甘南草地類型的生態(tài)經(jīng)濟(jì)價(jià)值與功能及其分類經(jīng)營(yíng).山地學(xué)報(bào),2004,22(6):655-660. Guo Z G,Gao X H,Liu X Y,Liang T G.Ecological economic value and functions and classification management for grassland in Gannan Prefecture,Gansu Province.Journal of Mountain Science,2004,22(6):655-660.(in Chinese)

    [34] 甘肅省草原總站.甘肅草地資源.蘭州:甘肅科技出版社,1999:258-261.

    [35] Turner D P,Ritts W D,Cohen W B,Gower S T,Running S W,Zhao M S,Costa M H,Kirschbaum A A,Ham J M,Saleska S R,Ahl D E.Evaluation of MODIS NPP and GPP products across multiple biomes.Remote Sensing of Environment,2006,102(3):282-292.

    [36] 梁天剛,馮琦勝,夏文韜,崔霞.甘南牧區(qū)草畜平衡優(yōu)化方案與管理決策.生態(tài)學(xué)報(bào),2011,31(4):1111-1123. Liang T G,Feng Q S,Xia W T,Cui X.Optimization strategy and management decision-making in balancing forage and livestock in Gannan pastoral area.Acta Ecologica Sinica,2011,31(4):1111-1123.(in Chinese)

    [37] 楊旭東,楊春,孟志興.我國(guó)草原生態(tài)保護(hù)現(xiàn)狀、存在問題及建議.草業(yè)科學(xué),2016,33(9):1901-1909. Yang X D,Yang C,Meng Z X.The current situation,problems and suggestions of grassland ecological protection in China.Pratacultural Science,2016,33(9):1901-1909.(in Chinese)

    [38] 張維云,公保才讓,仁青.海南州發(fā)展草產(chǎn)業(yè)的前景與對(duì)策.四川草原,2006(4):57-60. Zhang W Y,Gongbaocairang,Ren Q.Prospects and countermeasures for the development of grass industry in Hainan Prefecture.Sichuan Caoyuan,2006(4):57-60.(in Chinese)

    [39] 朱建國(guó),袁翀.甘南州發(fā)展草產(chǎn)業(yè)的前景與對(duì)策.草業(yè)科學(xué),2002,19(2):26-28. Zhu J G,Yuan C.Prospects and countermeasures for the development of grass industry in Gannan Prefecture.Pratacultural Science,2002,19(2):26-28.(in Chinese)

    [40] 任世龍,宜樹華,陳建軍,秦彧,王曉云.基于不同數(shù)碼相機(jī)和圖像處理方法的高寒草地植被蓋度估算的比較.草業(yè)科學(xué),2014,31(6):1007-1013. Ren S L,Yi S H,Chen J J,Qin Y,Wang X Y.Comparisons of alpine grassland fractional vegetation cover estimation using different cameras and different image analysis methods.Pratacultural Science,2014,31(6):1007-1013.(in Chinese)

    [41] 陳學(xué)兄,張小軍,陳永貴,常慶瑞.陜西省1998-2008年植被覆蓋度的時(shí)空變化研究.武漢大學(xué)學(xué)報(bào)·信息科學(xué)版,2013,38(6):674-678. Chen X X,Zhang X J,Chen Y G,Chang Q R.Spatiotemporal change of vegetation coverage in Shaanxi Province from 1998 to 2008.Geomatics and Information Science of Wuhan University,2013,38(6):674-678.(in Chinese)

    [42] 崔霞,梁天剛,劉勇.基于MOD09GA產(chǎn)品的草地生物量遙感估算模型.蘭州大學(xué)學(xué)報(bào):自然科學(xué)版,2009,45(5):79-87. Cui X,Liang T G,Liu Y.Modeling of aboveground biomass of grassland using remotely sensed MOD09GA data.Journal of Lanzhou University:Natural Sciences,2009,45(5):79-87.(in Chinese)

    [43] 馮琦勝,高新華,黃曉東,于惠,梁天剛.2001-2010年青藏高原草地生長(zhǎng)狀況遙感動(dòng)態(tài)監(jiān)測(cè).蘭州大學(xué)學(xué)報(bào):自然科學(xué)版,2011,47(4):75-81. Feng Q S,Gao X H,Huang X D,Yu H,Liang T G.Remote sensing dynamic monitoring of grass growth in Qinghai-Tibet plateau from 2001 to 2010.Journal of Lanzhou University:Natural Sciences,2011,47(4):75-81.(in Chinese)

    [44] 孟寶平,陳思宇,崔霞,馮琦勝,梁天剛.基于多源遙感數(shù)據(jù)的高寒草地生物量反演模型精度——以夏河縣??撇菰囼?yàn)區(qū)為例.草業(yè)科學(xué),2015,32(11):1730-1739. Meng B P,Chen S Y,Cui X,Feng Q S,Liang T G.The accuracy of grassland vegetation biomass estimated model based on multi-source remote sensing data——As a case of experimental area in Sangke grassland in Xiahe County.Pratacultural Science,2015,32(11):1730-1739.(in Chinese)

    [45] 趙軍,王小敏,李東成.基于MODIS的民勤綠洲植被覆蓋變化定量分析.干旱區(qū)資源與環(huán)境,2012,26(10):91-96. Zhao J,Wang X M,Li D C.Quantitative analysis on the vegetation coverage changes of Minqin Oasis based on MODIS.Journal of Arid Land Resources and Environment,2012,26(10):91-96.(in Chinese)

    [46] 王浩,李文龍,杜國(guó)禎,朱曉麗.基于3S技術(shù)的甘南草地覆蓋度動(dòng)態(tài)變化研究.草業(yè)學(xué)報(bào),2012,21(3):26-37. Wang H,Li W L,Du G Z,Zhu X L.Research on the change of grassland vegetation coverage using 3S technology in Gannan.Acta Prataculturae Sinica,2012,21(3):26-37.(in Chinese)

    [47] 葉貴祥,李維青,田源.基于NDVI的干旱區(qū)典型綠洲植被覆蓋動(dòng)態(tài)變化分析——以策勒綠洲為例.干旱區(qū)資源與環(huán)境,2009,23(9):128-133. Ye G X,Li W Q,Tian Y.Dynamic changes of vegetation cover in typical oasis of arid areas based on NDVI——A case study in Cele Oasis.Journal of Arid Land Resources and Environment,2009,23(9):128-133.(in Chinese)

    [48] 張仁華,田靜,李召良,蘇紅波,陳少輝.定量遙感產(chǎn)品真實(shí)性檢驗(yàn)的基礎(chǔ)與方法.中國(guó)科學(xué):地球科學(xué),2010,40(2):211-222. Zhang R H,Tian J,Li Z L,Su H B,Chen S H.Principles and methods for the validation of quantitative remote sensing products.Scientia Sinica Terrae,2010,40(2):211-222.(in Chinese)

    [49] 張?jiān)葡?張?jiān)骑w,李曉兵.地面測(cè)量與ASTER影像綜合計(jì)算植被蓋度.生態(tài)學(xué)報(bào),2007,27(3):964-976. Zhang Y X,Zhang Y F,Li X B.The synthetically estimating vegetation fractional coverage of grassland using field data and ASTER remote sensing imagine.Acta Ecologica Sinica,2007,27(3):964-976.(in Chinese)

    [50] 孟寶平,崔霞,楊淑霞,高金龍,胡遠(yuǎn)寧,陳思宇,梁天剛.基于Landsat 8 OLI和MODIS數(shù)據(jù)的高寒草地蓋度升尺度效應(yīng)研究——以夏河縣??撇菰囼?yàn)區(qū)為例.草業(yè)學(xué)報(bào),2016,25(7):1-12. Meng B P,Cui X,Yang S X,Gao J L,Hu Y N,Chen S Y,Liang T G.Scaling-up methodology for alpine grassland coverage monitoring based on Landsat 8 OLI and MODIS remote sensing data——A case study in Xiahe Sangke grassland.Acta Prataculturae Sinica,2016,25(7):1-12.(in Chinese)

    [51] 高林,楊貴軍,李紅軍,李振海,馮海寬,王磊,董錦繪,賀鵬.基于無人機(jī)數(shù)碼影像的冬小麥葉面積指數(shù)探測(cè)研究.中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2016,24(9):1254-1264. Gao L,Yang G J,Li H J,Li Z H,Feng H K,Wang L,Dong J H,He P.Winter wheat LAI estimation using unmanned aerial vehicle RGB-imaging.Chinese Journal of Eco-Agriculture,2016,24(9):1254-1264.(in Chinese)

    [52] 單成鋼.數(shù)字圖像技術(shù)在冬小麥群體監(jiān)測(cè)中的應(yīng)用研究.北京:中國(guó)農(nóng)業(yè)大學(xué)博士學(xué)位論文,2006. Shan C G.Application of digital image technology in winter wheat population monitoring.PhD Thesis.Beijing:China Agricultural University,2006.(in Chinese)

    [53] 李明,張長(zhǎng)利,房俊龍.基于圖像處理技術(shù)的小麥葉面積指數(shù)的提取.農(nóng)業(yè)工程學(xué)報(bào),2010,26(1):205-209. Li M,Zhang C L,Fang J L.Extraction of leaf area index of wheat based on image processing technique.Transactions of the CSAE,2010,26(1):205-209.(in Chinese)

    [54] 孫濤,劉振波,葛云健,顧祝軍.基于數(shù)碼相片Gamma校正的水稻葉面積指數(shù)估算.生態(tài)學(xué)報(bào),2014,34(13):3548-3557. Sun T,Liu Z B,Ge Y J,Gu Z J.Estimation of paddy rice leaf area index based on photo gamma correction.Acta Ecologica Sinica,2014,34(13):3548-3557.(in Chinese)

    (責(zé)任編輯 茍燕妮)

    Grassland fractional vegetation cover analysis using small UVAs and MODIS——A case study in Gannan Prefecture

    Song Qing-jie1, Cui Xia1, Zhang Yao-yao1, Meng Bao-ping2, Gao Jin-long2, Xiang Yu-xuan3

    (1.Key Laboratory of Western China’s Environmental Systems (Ministry of Education), College of Earth and Environmental Sciences, Lanzhou University, Lanzhou 730000, China;2.State Key Laboratory of Grassland Agro-ecosystems, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou 730020, China; 3.College of Earth and Environmental Sciences, Lanzhou University, Lanzhou 730000, China)

    Gannan Prefecture was selected as the study area. Photos of different sized subplots of grasslands were taken with a camera (Canon) and small UAVs (unmanned aerial vehicles). Terra/MODIS MOD13Q1 vegetation indices were obtained from May to October 2015. Correlation analyses were conducted between vegetation indices [enhanced vegetation index (EVI), normalized difference vegetation index (NDVI)] and grassland fractional vegetation cover (FVC), and regression models were determined. Based on the accuracy of different regression models, the best estimation model for grassland FVC in Gannan Prefecture was selected, and grassland spatio-temporal characteristics of FVC during the growing season were analysed. The results showed that 1) small UAVs carrying cameras could be used to obtain digital photos containing data on grassland FVC; 2) EVI was more correlated with grassland FVC than NDVI, and the logarithmic model based on EVI was the most suitable model for monitoring grassland FVC in Gannan Prefecture, with model accuracy reaching 88.00%; 3) grassland FVC of flat meadows in 2015 reached the maximum in August, and the remaining grassland types reached the maximum in July; 4) grasslands in Gannan Prefecture were dominated by high and mid-grassland FVC levels, and these grasslands were mainly distributed in Maqu County, Luqu County, and Xiahe County and Hezuo City. Overall, the grassland FVC in the mid-western and south-western areas was better than that in the eastern area. Establishment of the best model for grassland FVC was beneficial not only in elucidating the spatio-temporal distribution and seasonal dynamic variation of grassland FVC, but also to the maintenance of stable development of grassland ecosystems in Gannan Prefecture.

    Terra/MODIS; NDVI; EVI; correlation analysis; remote sensing inversion model; accuracy assessment; dynamic analysis of grassland fractional vegetation cover

    Cui Xia E-mail:xiacui@lzu.edu.cn

    10.11829/j.issn.1001-0629.2016-0046

    2016-01-22 接受日期:2016-04-05

    蘭州大學(xué)中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助(lzujbky-2015-140);國(guó)家自然科學(xué)基金項(xiàng)目(41401472、31372367);蘭州大學(xué)地理學(xué)基地科研訓(xùn)練及科研能力提高項(xiàng)目(J1210065)

    宋清潔(1992-),女,河南睢縣人,在讀碩士生,主要從事草地遙感研究。E-mail:songqj15@lzu.edu.cn

    崔霞(1984-),女,甘肅民勤人,講師,博士,主要從事環(huán)境遙感應(yīng)用研究。E-mail:xiacui@lzu.edu.cn

    S812-05;S127

    A

    1001-0629(2017)1-0040-11*

    宋清潔,崔霞,張瑤瑤,孟寶平,高金龍,向宇軒.基于小型無人機(jī)與MODIS數(shù)據(jù)的草地植被覆蓋度研究——以甘南州為例.草業(yè)科學(xué),2017,34(1):40-50.

    Song Q J,Cui X,Zhang Y Y,Meng B P,Gao J L,Xiang Y X.Grassland fractional vegetation cover analysis using small UVAs and MODIS——A case study in Gannan Prefecture.Pratacultural Science,2017,34(1):40-50.

    猜你喜歡
    甘南州蓋度植被指數(shù)
    黃山市近30 a植被蓋度時(shí)空變化遙感分析
    基于GIS的甘肅省甘南州垃圾填埋場(chǎng)適宜性分析
    甘肅科技(2020年21期)2020-04-13 00:33:40
    黃土高原地區(qū)植被蓋度對(duì)產(chǎn)流產(chǎn)沙的影響
    AMSR_2微波植被指數(shù)在黃河流域的適用性對(duì)比與分析
    河南省冬小麥產(chǎn)量遙感監(jiān)測(cè)精度比較研究
    Reflections on the Teaching Mode of English Reading in Senior High Schools
    基于DEM的1983-2012年甘南州降水空間變化趨勢(shì)及特征
    主要植被指數(shù)在生態(tài)環(huán)評(píng)中的作用
    西藏科技(2015年1期)2015-09-26 12:09:29
    坡位與植被蓋度對(duì)楓香造林成效的影響
    基于MODIS數(shù)據(jù)的植被指數(shù)與植被覆蓋度關(guān)系研究
    免费搜索国产男女视频| 国产精品美女特级片免费视频播放器| 日韩高清综合在线| 99久久人妻综合| 久久久久免费精品人妻一区二区| 久久热精品热| АⅤ资源中文在线天堂| 国产精品一区www在线观看| av在线老鸭窝| 我的老师免费观看完整版| 亚洲av成人精品一二三区| 乱人视频在线观看| 国产69精品久久久久777片| 成人一区二区视频在线观看| 22中文网久久字幕| 国产一区二区亚洲精品在线观看| 免费黄色在线免费观看| 色5月婷婷丁香| 亚洲精品自拍成人| 国产精品精品国产色婷婷| 国产精品美女特级片免费视频播放器| 97超碰精品成人国产| 亚洲av日韩在线播放| 色综合亚洲欧美另类图片| 日韩av在线大香蕉| 尤物成人国产欧美一区二区三区| 一个人看的www免费观看视频| av免费观看日本| 偷拍熟女少妇极品色| 国产中年淑女户外野战色| 免费人成在线观看视频色| 一个人看的www免费观看视频| 日韩精品有码人妻一区| 日本免费a在线| av女优亚洲男人天堂| 老女人水多毛片| 黄色日韩在线| 美女黄网站色视频| 欧美色视频一区免费| 久久精品国产鲁丝片午夜精品| 免费无遮挡裸体视频| 亚洲真实伦在线观看| 少妇的逼好多水| 欧美成人一区二区免费高清观看| 国产精品麻豆人妻色哟哟久久 | 黄色一级大片看看| 久久久久精品久久久久真实原创| 国模一区二区三区四区视频| 天堂网av新在线| 精品熟女少妇av免费看| 国产男人的电影天堂91| 午夜激情欧美在线| 人人妻人人看人人澡| 精品人妻偷拍中文字幕| 一级二级三级毛片免费看| 三级国产精品片| 久久久欧美国产精品| 午夜福利高清视频| 超碰av人人做人人爽久久| 国产黄色视频一区二区在线观看 | 一级二级三级毛片免费看| 建设人人有责人人尽责人人享有的 | 欧美日本视频| 黄色日韩在线| 欧美又色又爽又黄视频| 婷婷六月久久综合丁香| 有码 亚洲区| 国产精品精品国产色婷婷| 一个人观看的视频www高清免费观看| 2021少妇久久久久久久久久久| 国产黄色视频一区二区在线观看 | 久久国产乱子免费精品| 国产精品av视频在线免费观看| 长腿黑丝高跟| 最近中文字幕高清免费大全6| av又黄又爽大尺度在线免费看 | 一级av片app| 欧美激情国产日韩精品一区| 99久久精品热视频| 99在线视频只有这里精品首页| 99热精品在线国产| 村上凉子中文字幕在线| 毛片一级片免费看久久久久| videos熟女内射| 全区人妻精品视频| 色综合站精品国产| 精品无人区乱码1区二区| 中国美白少妇内射xxxbb| 欧美日本亚洲视频在线播放| 国产精品久久久久久av不卡| 啦啦啦韩国在线观看视频| 91久久精品电影网| 校园人妻丝袜中文字幕| 中文字幕人妻熟人妻熟丝袜美| 人体艺术视频欧美日本| 亚洲精品国产成人久久av| 久久久久性生活片| h日本视频在线播放| 久久久久久久久大av| 性色avwww在线观看| 99热这里只有是精品50| 久久99精品国语久久久| 99国产精品一区二区蜜桃av| 国产成人a区在线观看| 亚洲av日韩在线播放| 亚洲国产欧美人成| 蜜桃亚洲精品一区二区三区| 熟女电影av网| 久久精品影院6| 久久精品夜色国产| 国产免费男女视频| av免费在线看不卡| 韩国av在线不卡| 国产成人freesex在线| 热99在线观看视频| 禁无遮挡网站| 中文资源天堂在线| 美女国产视频在线观看| 好男人在线观看高清免费视频| 可以在线观看毛片的网站| ponron亚洲| 超碰av人人做人人爽久久| 亚洲人成网站在线播| 全区人妻精品视频| 色综合站精品国产| 别揉我奶头 嗯啊视频| 高清午夜精品一区二区三区| 嫩草影院精品99| 国产一级毛片在线| 日本-黄色视频高清免费观看| 99在线人妻在线中文字幕| 能在线免费观看的黄片| 伦理电影大哥的女人| 国产成人免费观看mmmm| 在线免费十八禁| 中文资源天堂在线| 国产精品,欧美在线| 简卡轻食公司| 黄色一级大片看看| 亚洲国产欧洲综合997久久,| 国产精品伦人一区二区| av.在线天堂| 久久久久久久午夜电影| av线在线观看网站| 日韩一本色道免费dvd| 免费看美女性在线毛片视频| 国产久久久一区二区三区| 日本-黄色视频高清免费观看| 欧美一区二区国产精品久久精品| 午夜福利成人在线免费观看| 国产爱豆传媒在线观看| 免费大片18禁| 国产色爽女视频免费观看| 亚洲精品乱码久久久久久按摩| 精品免费久久久久久久清纯| 一级爰片在线观看| 国产极品天堂在线| 性色avwww在线观看| 亚洲av中文字字幕乱码综合| 日本wwww免费看| 国产老妇伦熟女老妇高清| 97超碰精品成人国产| 日韩精品有码人妻一区| 成人av在线播放网站| 久久精品国产亚洲av涩爱| 免费av毛片视频| 最近2019中文字幕mv第一页| 国内精品一区二区在线观看| 成年女人永久免费观看视频| 欧美极品一区二区三区四区| 在线观看av片永久免费下载| 一级毛片久久久久久久久女| 午夜久久久久精精品| 国内少妇人妻偷人精品xxx网站| 久久99热这里只频精品6学生 | 精华霜和精华液先用哪个| 最近手机中文字幕大全| a级一级毛片免费在线观看| av播播在线观看一区| 国产精品久久视频播放| 91精品一卡2卡3卡4卡| 日本wwww免费看| 熟女人妻精品中文字幕| 99久久九九国产精品国产免费| 国产亚洲91精品色在线| 国产精品久久电影中文字幕| 男的添女的下面高潮视频| 日本av手机在线免费观看| av黄色大香蕉| 欧美日本亚洲视频在线播放| 日本一二三区视频观看| 欧美人与善性xxx| 亚洲一级一片aⅴ在线观看| 少妇人妻精品综合一区二区| 网址你懂的国产日韩在线| 变态另类丝袜制服| 亚洲人成网站在线播| ponron亚洲| 亚洲国产精品成人综合色| 大又大粗又爽又黄少妇毛片口| 美女大奶头视频| 一本一本综合久久| 国产精品乱码一区二三区的特点| 97在线视频观看| 高清毛片免费看| 国产一级毛片七仙女欲春2| 亚洲内射少妇av| 国产成人aa在线观看| 精品国内亚洲2022精品成人| 欧美三级亚洲精品| 亚洲国产精品专区欧美| 一区二区三区免费毛片| 国产乱人偷精品视频| av福利片在线观看| 亚洲国产精品专区欧美| 亚洲综合精品二区| 男女那种视频在线观看| 国内精品美女久久久久久| 国产精品一区二区三区四区久久| 一本一本综合久久| 国产淫语在线视频| 欧美精品一区二区大全| 亚洲真实伦在线观看| 久久久精品欧美日韩精品| 麻豆av噜噜一区二区三区| av在线亚洲专区| 97超视频在线观看视频| 如何舔出高潮| 一区二区三区高清视频在线| 亚洲综合精品二区| 网址你懂的国产日韩在线| 91久久精品国产一区二区三区| 亚洲美女搞黄在线观看| 久久99蜜桃精品久久| 亚洲欧美清纯卡通| 久久久久国产网址| 丰满乱子伦码专区| 午夜精品一区二区三区免费看| 亚洲乱码一区二区免费版| 毛片女人毛片| 亚洲av二区三区四区| 中文字幕熟女人妻在线| 啦啦啦观看免费观看视频高清| 日韩精品有码人妻一区| 免费看日本二区| 身体一侧抽搐| 91久久精品国产一区二区成人| 免费观看精品视频网站| 免费av毛片视频| 男人狂女人下面高潮的视频| 国产av在哪里看| 午夜精品在线福利| 成人亚洲欧美一区二区av| 一级二级三级毛片免费看| 久久久精品大字幕| 在线观看一区二区三区| 午夜激情福利司机影院| 亚洲欧美精品自产自拍| 一边摸一边抽搐一进一小说| 国产精品国产高清国产av| 欧美性猛交黑人性爽| 国产亚洲av片在线观看秒播厂 | 久久久久久伊人网av| 色视频www国产| 国产精品一区二区性色av| 天天躁日日操中文字幕| 伦精品一区二区三区| 国产成人福利小说| 国产 一区 欧美 日韩| 美女大奶头视频| 51国产日韩欧美| 色网站视频免费| 国产黄片视频在线免费观看| 国产精品不卡视频一区二区| 欧美一区二区国产精品久久精品| 国产成人精品一,二区| 欧美bdsm另类| 国产v大片淫在线免费观看| 天天躁夜夜躁狠狠久久av| 一边摸一边抽搐一进一小说| 国产精品一区二区在线观看99 | 午夜福利在线观看吧| 内射极品少妇av片p| 日本猛色少妇xxxxx猛交久久| 99久久成人亚洲精品观看| 日韩欧美精品v在线| 亚洲人与动物交配视频| 永久免费av网站大全| 国产精品久久视频播放| 亚洲经典国产精华液单| 少妇猛男粗大的猛烈进出视频 | 久久这里有精品视频免费| 国产精品一及| 国产av一区在线观看免费| 亚洲欧美精品专区久久| 国产精品一二三区在线看| 麻豆久久精品国产亚洲av| 欧美成人a在线观看| 不卡视频在线观看欧美| 视频中文字幕在线观看| 国产成人一区二区在线| 深爱激情五月婷婷| 国产伦在线观看视频一区| 麻豆国产97在线/欧美| 国产大屁股一区二区在线视频| 国产不卡一卡二| 国产精品一二三区在线看| 精品人妻偷拍中文字幕| 国产美女午夜福利| 成人毛片a级毛片在线播放| 禁无遮挡网站| 高清午夜精品一区二区三区| 国产在视频线精品| 日本黄色视频三级网站网址| 国语对白做爰xxxⅹ性视频网站| 伦理电影大哥的女人| 视频中文字幕在线观看| 国产探花极品一区二区| 免费看美女性在线毛片视频| 成人一区二区视频在线观看| 久久久久精品久久久久真实原创| 国产中年淑女户外野战色| 精品熟女少妇av免费看| 一级毛片久久久久久久久女| 免费观看精品视频网站| 国产成人精品婷婷| av在线老鸭窝| 亚洲熟妇中文字幕五十中出| 一级爰片在线观看| 亚洲av福利一区| 中文在线观看免费www的网站| 少妇丰满av| 99久国产av精品| 日本五十路高清| 51国产日韩欧美| 国内少妇人妻偷人精品xxx网站| 精品欧美国产一区二区三| 在线观看66精品国产| 日日摸夜夜添夜夜爱| 不卡视频在线观看欧美| 亚洲av熟女| 欧美+日韩+精品| 舔av片在线| 人人妻人人看人人澡| 久久精品综合一区二区三区| 成人午夜精彩视频在线观看| 一级毛片电影观看 | 色综合亚洲欧美另类图片| 亚洲不卡免费看| 久久精品夜夜夜夜夜久久蜜豆| 成年免费大片在线观看| 精品人妻视频免费看| 久99久视频精品免费| 国产极品精品免费视频能看的| 日韩高清综合在线| 你懂的网址亚洲精品在线观看 | 欧美精品一区二区大全| 亚洲成人久久爱视频| 国产爱豆传媒在线观看| 99热全是精品| 欧美日韩在线观看h| 久久精品久久久久久久性| 91午夜精品亚洲一区二区三区| 免费搜索国产男女视频| 日韩在线高清观看一区二区三区| 国产精品久久久久久精品电影| 国产69精品久久久久777片| 日本免费一区二区三区高清不卡| 精品免费久久久久久久清纯| 99热精品在线国产| 热99在线观看视频| 国产片特级美女逼逼视频| 天天一区二区日本电影三级| 天堂√8在线中文| 亚洲av不卡在线观看| 人妻夜夜爽99麻豆av| 久久6这里有精品| 亚洲av中文字字幕乱码综合| www.色视频.com| 午夜激情福利司机影院| 国产熟女欧美一区二区| 青春草亚洲视频在线观看| av在线老鸭窝| 日本午夜av视频| 午夜爱爱视频在线播放| 国产白丝娇喘喷水9色精品| 男女边吃奶边做爰视频| 午夜久久久久精精品| av国产免费在线观看| 免费黄色在线免费观看| 成人毛片60女人毛片免费| 久久久久九九精品影院| 一个人看的www免费观看视频| 舔av片在线| 99热这里只有是精品50| 色尼玛亚洲综合影院| 在线免费观看的www视频| a级毛色黄片| 22中文网久久字幕| 久久午夜福利片| 国产色婷婷99| 日本五十路高清| 午夜免费激情av| 国产黄片美女视频| 中国国产av一级| 欧美潮喷喷水| 高清日韩中文字幕在线| 99热这里只有是精品在线观看| 亚洲av日韩在线播放| 亚洲精品影视一区二区三区av| 日日摸夜夜添夜夜添av毛片| 热99在线观看视频| 99久国产av精品| 免费看a级黄色片| 精品一区二区三区视频在线| 人人妻人人澡人人爽人人夜夜 | 中文亚洲av片在线观看爽| 综合色丁香网| 国产在线男女| 亚洲怡红院男人天堂| 一级av片app| 一区二区三区高清视频在线| 亚洲精品影视一区二区三区av| 国产伦在线观看视频一区| 一边亲一边摸免费视频| 国产精品久久电影中文字幕| 免费播放大片免费观看视频在线观看 | 国产精品久久视频播放| 天天一区二区日本电影三级| 又粗又硬又长又爽又黄的视频| 成人二区视频| 中文天堂在线官网| 一级毛片久久久久久久久女| 又黄又爽又刺激的免费视频.| 国产精品一及| 高清毛片免费看| 亚洲真实伦在线观看| 51国产日韩欧美| 网址你懂的国产日韩在线| 噜噜噜噜噜久久久久久91| av专区在线播放| 久久精品国产自在天天线| 久久久久网色| 欧美日本视频| 少妇丰满av| 欧美3d第一页| 成人性生交大片免费视频hd| 美女高潮的动态| 亚洲国产日韩欧美精品在线观看| 99久久九九国产精品国产免费| 九九久久精品国产亚洲av麻豆| 久久国内精品自在自线图片| 丰满少妇做爰视频| 精品国内亚洲2022精品成人| 色噜噜av男人的天堂激情| 男人和女人高潮做爰伦理| 麻豆成人av视频| 亚洲一级一片aⅴ在线观看| 成年免费大片在线观看| av.在线天堂| 观看美女的网站| 国产免费又黄又爽又色| 国产精品国产三级专区第一集| 欧美性感艳星| 春色校园在线视频观看| 性插视频无遮挡在线免费观看| 乱码一卡2卡4卡精品| 欧美日韩一区二区视频在线观看视频在线 | 国产探花极品一区二区| 国产私拍福利视频在线观看| 免费看日本二区| 22中文网久久字幕| 亚洲精品影视一区二区三区av| 一级毛片aaaaaa免费看小| 蜜臀久久99精品久久宅男| 亚洲最大成人手机在线| 国产精品综合久久久久久久免费| 啦啦啦啦在线视频资源| 久久99热这里只有精品18| 国产免费又黄又爽又色| 青春草亚洲视频在线观看| 99久久九九国产精品国产免费| 日本一二三区视频观看| 天天躁夜夜躁狠狠久久av| 韩国av在线不卡| 最新中文字幕久久久久| 亚洲在线观看片| 久久精品综合一区二区三区| 一个人看视频在线观看www免费| 嫩草影院入口| 亚洲av日韩在线播放| 亚洲欧美一区二区三区国产| 国产真实伦视频高清在线观看| 亚洲不卡免费看| 国产精品一区二区三区四区久久| 欧美高清成人免费视频www| 欧美性感艳星| av天堂中文字幕网| 久99久视频精品免费| 亚洲欧美日韩无卡精品| 午夜久久久久精精品| 国产亚洲午夜精品一区二区久久 | 国产又黄又爽又无遮挡在线| 麻豆国产97在线/欧美| 国产伦理片在线播放av一区| 一区二区三区免费毛片| 纵有疾风起免费观看全集完整版 | 国产在线一区二区三区精 | 99久久成人亚洲精品观看| 国产精品一区二区三区四区久久| 欧美+日韩+精品| 嫩草影院入口| 亚洲精品自拍成人| 一卡2卡三卡四卡精品乱码亚洲| 中文精品一卡2卡3卡4更新| 欧美高清成人免费视频www| 91久久精品国产一区二区成人| 亚洲欧美精品自产自拍| a级毛色黄片| 男女下面进入的视频免费午夜| av在线老鸭窝| 青青草视频在线视频观看| 国产黄片美女视频| 成人性生交大片免费视频hd| 婷婷色综合大香蕉| 久久久久久久亚洲中文字幕| 欧美一区二区亚洲| 插逼视频在线观看| 韩国av在线不卡| 亚洲国产精品专区欧美| 波多野结衣巨乳人妻| 天天躁日日操中文字幕| 22中文网久久字幕| 高清毛片免费看| 99在线人妻在线中文字幕| 国产成人freesex在线| 亚洲欧美清纯卡通| 国产av在哪里看| 高清毛片免费看| 成人特级av手机在线观看| 精品少妇黑人巨大在线播放 | 又粗又爽又猛毛片免费看| 如何舔出高潮| 亚洲一级一片aⅴ在线观看| 日日摸夜夜添夜夜添av毛片| 免费搜索国产男女视频| 亚洲激情五月婷婷啪啪| 91精品国产九色| 亚洲av福利一区| 国产亚洲精品av在线| 18禁在线无遮挡免费观看视频| 22中文网久久字幕| 亚洲av不卡在线观看| 一边亲一边摸免费视频| 人妻制服诱惑在线中文字幕| 久久久国产成人免费| 国产一区二区三区av在线| 精品久久久噜噜| 深爱激情五月婷婷| 国产极品精品免费视频能看的| 国产视频内射| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品国内亚洲2022精品成人| 级片在线观看| 少妇丰满av| 午夜老司机福利剧场| 老司机福利观看| 天天躁夜夜躁狠狠久久av| 国产老妇伦熟女老妇高清| 日本免费在线观看一区| 99久久无色码亚洲精品果冻| 亚洲精华国产精华液的使用体验| 亚洲自偷自拍三级| 中文在线观看免费www的网站| 午夜a级毛片| 亚洲国产成人一精品久久久| 国产一区二区三区av在线| 男人狂女人下面高潮的视频| 国产一区二区三区av在线| 国产69精品久久久久777片| 亚洲国产精品sss在线观看| 国产亚洲一区二区精品| 国产av在哪里看| 欧美3d第一页| 97人妻精品一区二区三区麻豆| 欧美3d第一页| 五月伊人婷婷丁香| 午夜福利成人在线免费观看| 五月伊人婷婷丁香| 亚洲av电影不卡..在线观看| 有码 亚洲区| 国产精品一区二区三区四区免费观看| 亚洲欧美中文字幕日韩二区| 国产一区二区在线观看日韩| 国产精品电影一区二区三区| 69av精品久久久久久| 午夜福利高清视频| av黄色大香蕉| 偷拍熟女少妇极品色| 国产免费福利视频在线观看| 99久国产av精品| 一边亲一边摸免费视频| 看十八女毛片水多多多| 国产麻豆成人av免费视频| 搡女人真爽免费视频火全软件| 欧美高清成人免费视频www| 男人的好看免费观看在线视频| 91av网一区二区| 秋霞伦理黄片| 91精品伊人久久大香线蕉| 久久99精品国语久久久| 又黄又爽又刺激的免费视频.| 欧美激情在线99| 亚洲精品久久久久久婷婷小说 | 人体艺术视频欧美日本|