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

    基于Google Earth Engine的洱海流域建成區(qū)范圍及生態(tài)質(zhì)量變化分析

    2024-04-26 23:47:56鄭舒元陳星宇海燕王建雄
    湖北農(nóng)業(yè)科學(xué) 2024年2期

    鄭舒元 陳星宇 海燕 王建雄

    鄭舒元,陳星宇,海 燕,等. 基于Google Earth Engine的洱海流域建成區(qū)范圍及生態(tài)質(zhì)量變化分析[J]. 湖北農(nóng)業(yè)科學(xué),2024, 63(2):199-204.

    摘要:以洱海流域?yàn)檠芯繀^(qū),借助Google Earth Engine云平臺(tái),利用Landsat 8 OLI系列影像,合成研究區(qū)域的增強(qiáng)的指數(shù)型建筑指數(shù)(EIBI)并對(duì)洱海流域建成區(qū)進(jìn)行提取,再對(duì)研究區(qū)域遙感生態(tài)指數(shù)進(jìn)行合成,通過對(duì)建成區(qū)范圍與流域范圍內(nèi)遙感生態(tài)指數(shù)進(jìn)行空間相關(guān)分析,從而對(duì)2014—2021年洱海流域及其建成區(qū)進(jìn)行生態(tài)評(píng)估。結(jié)果表明,EIBI可有效增強(qiáng)建筑信息,使用EIBI可有效、準(zhǔn)確地提取研究區(qū)域建成區(qū),經(jīng)過EIBI提取,洱海流域2014—2021年建成區(qū)擴(kuò)張面積共14.712 km2;2014—2021年洱海流域生態(tài)明顯改善,遙感生態(tài)指數(shù)增加0.008,其上升速率為0.001/年,流域內(nèi)超過65%的區(qū)域生態(tài)得到改善;洱海流域建成區(qū)與生態(tài)未得到改善區(qū)域大致匹配,建成區(qū)經(jīng)過生態(tài)修復(fù)后生態(tài)質(zhì)量未見顯著提升,且建成區(qū)生態(tài)質(zhì)量呈逐年下降趨勢(shì),說明建成區(qū)對(duì)環(huán)境改善具有阻力。

    關(guān)鍵詞:洱海流域;增強(qiáng)的指數(shù)型建筑指數(shù);遙感生態(tài)指數(shù);Google Earth Engine;流域環(huán)境

    中圖分類號(hào):P237;X321? ? ? ? ?文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):0439-8114(2024)02-0199-06

    DOI:10.14088/j.cnki.issn0439-8114.2024.02.030 開放科學(xué)(資源服務(wù))標(biāo)識(shí)碼(OSID):

    Analysis of built-up area range and ecological quality changes in the Erhai Basin?based on Google Earth Engine

    ZHENG Shu-yuan1,2, CHEN Xing-yu3, HAI Yan1,2, WANG Jian-xiong1,2

    (1. Yunnan Agricultural University, Kunming? 650201,China; 2. Agricultural Remote Sensing and Precision Agriculture Engineering Research Center of Yunnan Universities, Kunming? 650201,China;3. College of Economics and Management,Southwest Forestry University, Kunming? 650224,China)

    Abstract: Taking Erhai Basin as the research area, with the help of Google Earth Engine cloud platform and Landsat 8 OLI series images, the enhanced index-based build-up index (EIBI) of the research area was synthesized and the built-up areas in Erhai Basin was extracted. Then, the remote sensing-based ecological index of the study area was synthesized. Through the spatial correlation analysis of the remote sensing-based ecological index of the built-up area and the basin, the ecological assessment of Erhai Basin and its built-up area from 2014 to 2021 was carried out. The results showed that the EIBI could effectively enhance the building information. Using EIBI could effectively and accurately extract the built-up area of the study area. After EIBI extraction, the expansion area of the built-up area of Erhai Basin from 2014 to 2021 was 14.712 km2. From 2014 to 2021, the ecology of Erhai Basin was significantly improved, the remote sensing ecological index increased by 0.008, with an increase rate of 0.001/year, and more than 65% of the regional ecology in the basin was improved. The built-up area in Erhai Basin roughly matched the area where the ecology had not been improved. The ecological quality of the built-up area had not been significantly improved after ecological restoration, and the ecological quality of the built-up area showed a downward trend year by year, indicating that the built-up area had resistance to environmental improvement.

    Key words: Erhai Basin; enhanced index-based built-up index; remote sensing-based ecological index; Google Earth Engine; basin environment

    洱海作為云南省九大高原湖泊之一,其流域?yàn)樯鷳B(tài)系統(tǒng)的重要組成部分。21世紀(jì)以來,洱海流域周邊土地進(jìn)行了不同程度的開發(fā)利用,導(dǎo)致建成區(qū)不斷擴(kuò)張,人口密度的上升、建成區(qū)不斷擴(kuò)張以及對(duì)植被的破壞使得洱海流域生態(tài)質(zhì)量下降問題日漸突出,對(duì)洱海流域生態(tài)質(zhì)量進(jìn)行評(píng)估對(duì)其生態(tài)建設(shè)具有重要意義。

    當(dāng)前對(duì)生態(tài)環(huán)境遙感監(jiān)測(cè)的研究有了較大發(fā)展,大多研究從土地利用變化、城市熱島效應(yīng)以及不透水面蓋度等角度對(duì)生態(tài)效應(yīng)影響進(jìn)行評(píng)價(jià),其多為單角度及單一生態(tài)因子對(duì)生態(tài)效應(yīng)進(jìn)行評(píng)價(jià),徐涵秋[1,2]于2013年提出遙感生態(tài)指數(shù)(Remote sensing based ecological index,RSEI),其利用主成分分析法對(duì)綠度、濕度、干度及熱度生態(tài)指標(biāo)進(jìn)行合成,規(guī)避了人為確定權(quán)重而導(dǎo)致評(píng)價(jià)結(jié)果的不確定性,其計(jì)算結(jié)果較為穩(wěn)定客觀,能夠用于對(duì)研究區(qū)的生態(tài)環(huán)境遙感綜合評(píng)價(jià)[3-10],近年來,吳志杰等[11]以福建省龍巖市為例,利用遙感生態(tài)指數(shù)對(duì)山地城市進(jìn)行生態(tài)遙感變化分析,對(duì)其環(huán)境質(zhì)量做出評(píng)估;張娟等[12]以河南省焦作市為例,利用遙感生態(tài)指數(shù)對(duì)資源枯竭型城市生態(tài)環(huán)境變化影響進(jìn)行評(píng)估;王志超等[13]利用遙感生態(tài)指數(shù)與植被覆蓋度相結(jié)合,對(duì)成都市錦江區(qū)的生態(tài)質(zhì)量做出評(píng)估并證明RSEI能夠反映2002—2018年錦江區(qū)生態(tài)質(zhì)量及變化情況;劉栩位等[14]利用遙感生態(tài)指數(shù)對(duì)三峽庫(kù)區(qū)重慶段水土保持生態(tài)功能區(qū)生態(tài)環(huán)境質(zhì)量狀況進(jìn)行了評(píng)估;范德芹等[15]利用改進(jìn)的遙感生態(tài)指數(shù)對(duì)神府礦區(qū)的生態(tài)環(huán)境變化進(jìn)行綜合評(píng)價(jià),得到目標(biāo)區(qū)域生態(tài)逐漸得到修復(fù),生態(tài)環(huán)境質(zhì)量逐漸提升的結(jié)論。

    傳統(tǒng)的遙感數(shù)據(jù)處理方法對(duì)長(zhǎng)時(shí)序研究而言工作量較大,近年來云平臺(tái)的發(fā)展使得遙感數(shù)據(jù)的存儲(chǔ)計(jì)算及可視化得到了新的處理途徑[16-19],Google Earth Engine(下文簡(jiǎn)稱GEE)平臺(tái)由Google公司提供技術(shù)支撐,其提供大量全球尺度不同分辨率衛(wèi)星影像數(shù)據(jù),能夠?yàn)殚L(zhǎng)期監(jiān)測(cè)和長(zhǎng)時(shí)序數(shù)據(jù)處理提供有力支撐[20-25],Hamud等[26]利用GEE對(duì)索馬里巴納迪爾30年的城市擴(kuò)展及土地利用進(jìn)行了提取分析。本研究借助GEE獲取2014—2021年洱海流域的Landsat 8 OLI系列影像,通過GEE計(jì)算研究區(qū)域增強(qiáng)的指數(shù)型建筑指數(shù)(Enhanced index-based built-up index,EIBI)提取該流域建成區(qū),再分別計(jì)算流域綠度指標(biāo)、濕度指標(biāo)、干度指標(biāo)及熱度指標(biāo),通過以上計(jì)算結(jié)果對(duì)RSEI指數(shù)進(jìn)行合成,并與提取的建成區(qū)進(jìn)行疊加統(tǒng)計(jì)分析,從而對(duì)洱海流域生態(tài)環(huán)境質(zhì)量進(jìn)行評(píng)估。

    1 研究區(qū)概況

    研究區(qū)域位于云南省中部偏西洱海流域,地處北緯25°25′—26°10′,東經(jīng)99°32′—100°27′,流域面積2 607.75 km2,湖面高程1 966 m時(shí)湖面面積為252 km2,一般湖水面積約246 km2,蓄水量約29.5×108 m3,呈狹長(zhǎng)形,流域內(nèi)河流呈現(xiàn)聚集狀匯入中心湖泊,境內(nèi)有大小河溪共117條,主要入湖溝渠125條,轄大理市、洱源縣2市縣16個(gè)鄉(xiāng)鎮(zhèn)。具有供水、灌溉、發(fā)電、調(diào)節(jié)氣候、漁業(yè)、航運(yùn)及旅游七大功能,是云南省九大高原湖泊之一,干濕季分明,5—10月為雨季,11月至次年4月為旱季,多年年均溫度及降雨量分別為15.1 ℃和1 048 mm。

    2 數(shù)據(jù)來源與研究方法

    2.1 數(shù)據(jù)來源及預(yù)處理

    研究采用Landsat 8地表反射率(Surface reflectance,SR)數(shù)據(jù)作為數(shù)據(jù)源,其從Google Earth Engine(GEE,https://earthengine.google.com/)平臺(tái)直接獲Landsat 8 OLI數(shù)據(jù),時(shí)間跨度為2014—2021年,選取2014年及2021年1月相同時(shí)段影像數(shù)據(jù)進(jìn)行對(duì)比分析。數(shù)據(jù)受到的云、雨雪等不利因素利用CFMASK方法去除,再利用頭文件信息對(duì)影像進(jìn)行輻射定標(biāo)及大氣校正,最終得到可用于處理的影像數(shù)據(jù)。

    2.2 建筑擴(kuò)張信息提取

    不同地物具有不同光譜特征,且不同地物對(duì)太陽(yáng)光的反射和吸收的程度不同,其反射光譜曲線具有顯著差異,利用這一差異可區(qū)別土地利用類型,從而進(jìn)一步提取感興趣區(qū),通過分析并找出城市用地與其他地物之間的差距,可將建筑用地信息與其他地物相區(qū)別開來,從而提取建成區(qū)。

    Zha等[27]改進(jìn)了歸一化建筑指數(shù)(Normalized differences built-up index,NDBI)用于提取建筑用地,該指數(shù)思想與歸一化植被指數(shù)NDVI相似,經(jīng)過對(duì)光譜的進(jìn)一步分析而得來,其表達(dá)式為:

    該公式由近紅外與紅光波段計(jì)算得來,ρnir與ρred分別表示近紅外波段與紅波段的反射率,在這2個(gè)波段之間,建筑用地、裸土與沙地的灰度值逐漸增加,其他值則逐漸下降,因此圖像上大于0的像元值可歸為建筑用地及沙地裸土。

    徐涵秋[28]提出的基于壓縮數(shù)據(jù)維方法提取建筑用地信息的指數(shù),又稱指數(shù)型建筑用地指數(shù)(Index-based Built-up Index,IBI),采用原始波段衍生出的3個(gè)指數(shù)波段來構(gòu)建新的影像對(duì)建筑用地進(jìn)行提取,其采用了土壤植被調(diào)節(jié)指數(shù)(Soil adjusted vegetation index,SAVI)、修正歸一化水體指數(shù)(Modified normalized differences water index,MNDWI)和歸一化建筑指數(shù)(NDBI),模型如下所示。

    式中,ρnir、ρred、ρgreen、ρmir1分別表示近紅外波段、紅波段、綠波段、中紅外波段的反射率;L代表常量,默認(rèn)值為0.5。

    由于研究區(qū)域包含大量山地、裸土區(qū)域,在分析時(shí)需要剔除大部分干擾因素的影響從而準(zhǔn)確提取目標(biāo)區(qū)域,為準(zhǔn)確提取研究區(qū)域建成區(qū),選用吳志杰等[29]提出的增強(qiáng)的指數(shù)型建筑指數(shù)(Enhanced index-based built-up index,EIBI)來對(duì)研究區(qū)建成區(qū)進(jìn)行提取,其對(duì)富含黏土型礦物的地表裸土有很好的抑制作用,經(jīng)過對(duì)漳州市及福州市研究區(qū)提取結(jié)果對(duì)比驗(yàn)證,證實(shí)其能夠準(zhǔn)確區(qū)分研究區(qū)域內(nèi)裸土和建筑物,精度可達(dá)90%以上;代穎懿等[30]基于Landsat-7 影像,利用EIBI提取福州市區(qū)建成區(qū)的精度可達(dá)92%。因此,選擇EIBI來提取研究區(qū)的建成區(qū)信息,其公式模型如下:

    其中,EIBI由歸一化差值裸地與建筑用地指數(shù)(Normalized difference bareness and built-up Index,NDBBI)、裸土指數(shù)(Bare soil index,BSI)、增強(qiáng)型裸土指數(shù)(Enhanced bare soil index,EBSI)計(jì)算得到,ρblue為藍(lán)光波段反射率,ρgreen為綠光波段反射率,ρred為紅光波段反射率,ρnir為近紅外波段反射率,ρswir1和 ρswir2為短波紅外波段反射率。且在使用 NDBBI、SAVI、MNDWI、BSI、EBSI等指數(shù)模型生成指數(shù)波段時(shí)不能直接使用歸一化數(shù)據(jù)進(jìn)行計(jì)算,其中NDBBI要經(jīng)過0~255的灰度值拉伸。運(yùn)用EIBI建筑指數(shù)對(duì)2014—2021年洱海流域建成區(qū)進(jìn)行提取結(jié)果如圖1所示,通過將EIBI提取結(jié)果與原始影像進(jìn)行疊加,并進(jìn)行目視解譯,該指數(shù)能夠準(zhǔn)確剔除研究區(qū)內(nèi)山體及裸土,提取建成區(qū)。

    2.3 RSEI指數(shù)計(jì)算

    遙感生態(tài)指數(shù)(Remote sensing based ecological index,RSEI)是一種基于濕度、綠度、干度和熱度4種自然因子為主的評(píng)價(jià)生態(tài)狀況的體系,可以快速評(píng)價(jià)生態(tài)環(huán)境狀況。需利用綠度指標(biāo)、濕度指標(biāo)、干度指標(biāo)及熱度指標(biāo)4個(gè)波段來進(jìn)行合成,再對(duì)合成的多波段影像進(jìn)行主成分分析得到,其中綠度指標(biāo)為歸一化植被指數(shù)(Normalized difference vegetation index,NDVI);濕度指標(biāo)為地表濕度(Wet),具體表現(xiàn)為纓帽變換的第三分量;干度指標(biāo)為建筑指數(shù)與裸土指數(shù)合成干度指數(shù)(Normalized difference soil index,NDSI),由建筑指數(shù)(Index-based built-up index,IBI)與土壤指數(shù)(Soil index,SI)合成;熱度指標(biāo)為地表溫度(Land surface temperaturae,LST)。以上指標(biāo)計(jì)算公式如下所示。

    式中,LST為地表溫度;ε為地表比輻射率;λ為熱紅外波段的中心波長(zhǎng);ρ=14 380 umk,T=k2-ln (k1/L+1),L=gainDN+bias,其中DN表示衛(wèi)星傳感器記錄下的數(shù)字信號(hào)值轉(zhuǎn)化來的反射率,gain為傳感器增益,bias為傳感器偏移。k1、k2、gain、bias 均可在影像頭文件中獲取。計(jì)算完4個(gè)波段的指數(shù)后,將其進(jìn)行多波段合成,其合成順序?yàn)镹DVI、Wet、NDSI、LST,對(duì)合成后的新影像進(jìn)行主成分分析(PCA),提取主成分分析后所得的多波段影像的第一波段,再將第一主成分進(jìn)行歸一化處理,歸一化后的第一主成分即為遙感生態(tài)指數(shù)(RSEI),其歸一化的表達(dá)式為:

    其中,PC1min為第一主成分最小值,PC1max為第一主成分最大值。通過對(duì)數(shù)據(jù)進(jìn)行輻射定標(biāo)、大氣校正、幾何校正、地形校正等預(yù)處理后,對(duì)綠度指標(biāo)、濕度指標(biāo)、干度指標(biāo)、熱度指標(biāo)進(jìn)行提取并進(jìn)行歸一化處理與主成分分析,可得到4個(gè)指標(biāo)分量的主成分分析結(jié)果(表1)。

    由表1可知,2期影像的第一主成分貢獻(xiàn)率分別為86.72%、84.73%,說明第一主成分PC1包含了影像絕大部分特征信息,且濕度指標(biāo)與綠度指標(biāo)的荷載值為正值,其對(duì)遙感生態(tài)指數(shù)為正面貢獻(xiàn);干度指標(biāo)與熱度指標(biāo)的荷載值為負(fù)值,說明其對(duì)遙感生態(tài)指數(shù)為負(fù)面貢獻(xiàn)。綜上,PC1可充分集成4個(gè)指標(biāo)分量的信息,該研究區(qū)域可采用第一主成分PC1來提取遙感生態(tài)指數(shù)(RSEI)。

    RSEI越大表明所處地區(qū)環(huán)境越好,其通過不同閾值(間隔為0.2)分級(jí)可分為1—5級(jí),其生態(tài)狀況根據(jù)分級(jí)可分為差、較差、中等、良、優(yōu)5個(gè)等級(jí),通過對(duì)2014年與2020年洱海流域建成區(qū)提取結(jié)果與重分類后所得的洱海流域RSEI指數(shù)進(jìn)行疊加分析,可掌握洱海流域建成區(qū)的生態(tài)情況變化,經(jīng)計(jì)算后的洱海流域2014年、2020年RSEI如圖2所示。

    3 結(jié)果與分析

    3.1 2014—2021年洱海流域建筑擴(kuò)張變化及影響

    利用2021年建成區(qū)與2014年建成區(qū)相比較,可以計(jì)算得到8年間增長(zhǎng)的建成區(qū)面積與建成區(qū)區(qū)域的遙感生態(tài)指數(shù),通過對(duì)不同年份之間的建成區(qū)指數(shù)之間作差即可得出2014—2021年建成區(qū)范圍變化,8年間建成區(qū)面積與遙感生態(tài)指數(shù)變化結(jié)果如表2所示,建成區(qū)范圍變化如圖1所示。從表2可看出,2014—2021年,建成區(qū)面積增加了14.712 km2,而建成區(qū)遙感生態(tài)指數(shù)降低了0.008,洱海流域內(nèi)建筑面積年均擴(kuò)展為1.839 km2,其建成區(qū)的遙感生態(tài)指數(shù)年均降低0.001 3。

    由圖1可知,洱海流域建成區(qū)擴(kuò)展主要集中于洱海流域下游大理市及杭瑞高速途經(jīng)流域路段,流域上游洱源縣及流域西側(cè)喜洲鎮(zhèn)、古城區(qū)以及西景線、大麗線沿途。由于政策原因,洱海東岸海東鎮(zhèn)開發(fā)速度有所減緩,2014—2021年該區(qū)域未見顯著擴(kuò)張。對(duì)2014年建成區(qū)與2021年建成區(qū)遙感生態(tài)指數(shù)進(jìn)行提取,2014年建成區(qū)遙感生態(tài)指數(shù)為0.811,2021年建成區(qū)遙感生態(tài)指數(shù)為0.800,可見建成區(qū)遙感生態(tài)指數(shù)明顯低于洱海流域整體遙感生態(tài)指數(shù),且2014—2021年建成區(qū)遙感生態(tài)指數(shù)呈逐年下降趨勢(shì),通過對(duì)2014年與2021年建成區(qū)RSEI做差可分析建成區(qū)內(nèi)RSEI整體變化趨勢(shì),其結(jié)果如圖3所示。

    由圖3可知,RSEI呈上升趨勢(shì)的區(qū)域占整個(gè)建成區(qū)的47.54%,RSEI呈下降趨勢(shì)的區(qū)域占整個(gè)建成區(qū)的52.06%,其中呈下降趨勢(shì)的區(qū)域,遙感生態(tài)指數(shù)下降的速率為0.043/年,呈上升趨勢(shì)的區(qū)域上升速率為0.029/年,綜合說明建成區(qū)整體生態(tài)呈下降趨勢(shì),生態(tài)質(zhì)量下降來源于建成區(qū)的擴(kuò)展導(dǎo)致的環(huán)境破壞及原有建成區(qū)生態(tài)質(zhì)量的下滑,證明建成區(qū)擴(kuò)展對(duì)洱海流域生態(tài)有負(fù)面影響。

    3.2 洱海流域生態(tài)質(zhì)量時(shí)空變化分析

    為能夠更好反映洱海流域的生態(tài)質(zhì)量變化,根據(jù)徐涵秋[1]對(duì)RSEI的評(píng)估等級(jí)將研究區(qū)內(nèi)RSEI進(jìn)行分級(jí),生態(tài)質(zhì)量差為1級(jí)(0~0.2),生態(tài)質(zhì)量較差為2級(jí)(0.2~0.4),生態(tài)質(zhì)量中等為3級(jí)(0.4~0.6),生態(tài)質(zhì)量良好為4級(jí)(0.6~0.8),生態(tài)質(zhì)量?jī)?yōu)為5級(jí)(0.8~1.0),對(duì)研究區(qū)RSEI進(jìn)行分級(jí)統(tǒng)計(jì),其結(jié)果如表3所示。

    由表3可知,2014—2021年洱海流域生態(tài)指數(shù)等級(jí)以優(yōu)(5級(jí))為主,其占比均超過85%,該時(shí)間段內(nèi)1級(jí)、2級(jí)、3級(jí)占比均不足1%;研究時(shí)間段內(nèi)洱海流域生態(tài)指數(shù)大幅上升,超過65%的4級(jí)地升為5級(jí)地。圖4從空間上直觀體現(xiàn)洱海流域的RSEI指數(shù)等級(jí)變化,可見研究區(qū)整體生態(tài)環(huán)境有巨大改善,大部分區(qū)域RSEI評(píng)級(jí)得以上升;未見顯著改善區(qū)域經(jīng)與建成區(qū)提取結(jié)果疊加分析可得,研究區(qū)生態(tài)未見顯著改善區(qū)域基本與建成區(qū)提取結(jié)果相吻合,2014年建成區(qū)所在的區(qū)域生態(tài)評(píng)價(jià)等級(jí)處于較低水平,說明建成區(qū)擴(kuò)張對(duì)研究區(qū)域的環(huán)境有負(fù)面影響。

    綜上,2014—2021年洱海流域的遙感生態(tài)指數(shù)從0.869上升到0.877,共上升了0.008,表明研究區(qū)域生態(tài)質(zhì)量呈逐年改善趨勢(shì),年均上升速度為0.001/年,與建成區(qū)生態(tài)呈下降趨勢(shì)不同,由于國(guó)家于研究時(shí)間段內(nèi)對(duì)洱海流域生態(tài)進(jìn)行整治與保護(hù)以及對(duì)洱海流域提出轉(zhuǎn)型發(fā)展要求,對(duì)洱海周邊的建成區(qū)進(jìn)行了轉(zhuǎn)型及調(diào)整,于2018—2019年時(shí)間段內(nèi)關(guān)停了洱海流域內(nèi)所有礦山,以及對(duì)原有礦山設(shè)施進(jìn)行拆除,對(duì)洱海流域?qū)嵤┝饔蛎嫔骄G化建設(shè)與退耕還林和陡坡地生態(tài)治理工程,實(shí)施西湖居民區(qū)生態(tài)搬遷工程及分類清理建設(shè)用地,對(duì)整個(gè)洱海流域生態(tài)質(zhì)量有了較大提高,使得整個(gè)研究區(qū)生態(tài)呈改善趨勢(shì)。

    4 結(jié)論

    利用Google Earth Engine云平臺(tái),并利用Landsat8 OLI系列遙感數(shù)據(jù),對(duì)洱海流域建成區(qū)擴(kuò)張分析后可得以下結(jié)論。

    1)使用增強(qiáng)的指數(shù)型建筑指數(shù)(EIBI)可有效抑制富含黏土礦物的裸土信息,而增強(qiáng)建成區(qū)信息,使用EIBI可有效、準(zhǔn)確地提取研究區(qū)域建成區(qū),經(jīng)過EIBI提取,洱海流域2014—2021年建筑擴(kuò)張面積共14.712 km2

    2)使用遙感生態(tài)指數(shù)(RSEI)對(duì)洱海流域建筑擴(kuò)張地域進(jìn)行生態(tài)評(píng)估,評(píng)估結(jié)果為2014—2021年洱海流域生態(tài)環(huán)境得到顯著改善,遙感生態(tài)指數(shù)逐年上升,其上升速率為0.001/年,2014—2021年流域內(nèi)超過65%的區(qū)域生態(tài)得到改善。

    3)洱海流域建成區(qū)與生態(tài)未得到改善區(qū)域大致匹配,建成區(qū)經(jīng)過生態(tài)修復(fù)后生態(tài)質(zhì)量未見顯著提升,且建成區(qū)生態(tài)質(zhì)量呈逐年下降趨勢(shì),說明建成區(qū)對(duì)環(huán)境改善具有阻力。

    參考文獻(xiàn):

    [1] 徐涵秋.城市遙感生態(tài)指數(shù)的創(chuàng)建及其應(yīng)用[J].生態(tài)學(xué)報(bào),2013,33(24):7853-7862.

    [2] 徐涵秋. 區(qū)域生態(tài)環(huán)境變化的遙感評(píng)價(jià)指數(shù)[J]. 中國(guó)環(huán)境科學(xué),2013,33(5):889-897.

    [3] HU X S, XU H Q. A new remote sensing index for assessing the spatial heterogeneity in urban ecological quality: A case from Fuzhou City, China[J]. Ecological indicators, 2018, 89:11-21.

    [4] 王東升, 王小磊, 雷澤勇. 基于遙感生態(tài)指數(shù)的阜新市生態(tài)質(zhì)量評(píng)估[J]. 生態(tài)科學(xué), 2020, 39(3):88-94.

    [5] 李 妍,張國(guó)欽,吝 濤,等. 鄉(xiāng)鎮(zhèn)遙感生態(tài)指數(shù)時(shí)空變化及影響因子研究——以天津市薊州區(qū)為例[J]. 生態(tài)學(xué)報(bào),2022,42(2):474-486.

    [6] ZHU D Y,CHEN T,ZHEN N, et al. Monitoring the effects of open-pit mining on the eco-environment using a moving window-based remote sensing ecological index[J]. Environmental science and pollution research, 2020, 27(13):15716-15728.

    [7] 范小安, 王金輝, 徐美娜,等. 典型入??诤0稁鷳B(tài)環(huán)境質(zhì)量遙感評(píng)價(jià)研究——以寶山區(qū)為例[J]. 海洋學(xué)研究,2020,38(2):66-73.

    [8] 喬 敏,張永彬,宇林軍.基于遙感生態(tài)指數(shù)的北京市生態(tài)環(huán)境質(zhì)量動(dòng)態(tài)監(jiān)測(cè)[J].華北理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2020,? ? ? 42(4):17-23.

    [9] 張 寒,潘 安,羅 芳.基于遙感生態(tài)指數(shù)的南充市生態(tài)環(huán)境質(zhì)量評(píng)價(jià)[J].西華師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2020,41(4):395-403.

    [10] 項(xiàng) 頌,吳 越,呂興菊,等. 洱海流域農(nóng)業(yè)面源污染空間分布特征及分類控制策略[J].環(huán)境科學(xué)研究,2020,33(11):2474-2483.

    [11] 吳志杰,鄒 丹,黃紹霖.山地城市擴(kuò)展及生態(tài)變化遙感分析[J].亞熱帶資源與環(huán)境學(xué)報(bào),2016,11(3):81-87,94.

    [12] 張 娟,文廣超,王恩營(yíng),等.基于遙感生態(tài)指數(shù)的焦作市生態(tài)環(huán)境動(dòng)態(tài)監(jiān)測(cè)與評(píng)價(jià)[J].水土保持通報(bào),2020,40(6):107-114.

    [13] 王志超,何新華.基于植被覆蓋度和遙感生態(tài)指數(shù)的成都市錦江區(qū)生態(tài)質(zhì)量評(píng)估[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2021,37(4):492-500.

    [14] 劉栩位,周啟剛,周 浪,等.基于RSEI的三峽庫(kù)區(qū)重慶段水土保持生態(tài)功能區(qū)生態(tài)環(huán)境質(zhì)量動(dòng)態(tài)監(jiān)測(cè)[J].水土保持研究,2021,28(5):278-286.

    [15] 范德芹,邱 玥,孫文彬,等. 基于遙感生態(tài)指數(shù)的神府礦區(qū)生態(tài)環(huán)境評(píng)價(jià)[J].測(cè)繪通報(bào),2021(7):23-28.

    [16] SHETTY S,GUPTA P K,BELGIU M,et al. Assessing the effect of training sampling design on the performance of machine learning classifiers for land cover mapping using multi-temporal remote sensing data and Google Earth Engine[J]. Remote sensing,2021,13(8):1-22.

    [17] LI K,WANG L,YIN D M. Deriving corn and soybeans fractions with Land Remote-Sensing Satellite(System, Landsat) imagery by accounting for endmember variability on Google Earth Engine[J]. International journal of remote sensing,2021,42(12):4493-4513.

    [18] ZHANG Z M,WEI M Y,PU D C,et al. Assessment of annual composite images obtained by Google Earth Engine for urban areas mapping using random forest[J]. Remote sensing,2021,13(4):748.

    [19] 鄭子豪,吳志峰,陳穎彪,等.基于Google Earth Engine的長(zhǎng)三角城市群生態(tài)環(huán)境變化與城市化特征分析[J].生態(tài)學(xué)報(bào),2021,41(2):717-729.

    [20] 廖文秀,陳奕云,趙 曦,等. 基于GEE的湖北省近30年湖泊及其岸線演變分析[J].湖北農(nóng)業(yè)科學(xué),2021,60(10):46-54,59.

    [21] 張禎祺,蔡惠文,張平平,等.基于GEE遙感云平臺(tái)的三江源植被碳源/匯時(shí)空變化研究[J].自然資源遙感,2023,35(1):231-242.

    [22] 谷 雷,岳彩榮,張國(guó)飛,等.基于Google Earth Engine平臺(tái)的大湄公河次區(qū)域2001—2019年植被NPP時(shí)空變化分析[J].西部林業(yè)科學(xué),2021,50(2):132-139.

    [23] 楊坤士,盧 遠(yuǎn),翁月梅,等.Google Earth Engine平臺(tái)支持下的南流江流域生態(tài)環(huán)境質(zhì)量動(dòng)態(tài)監(jiān)測(cè)[J].農(nóng)業(yè)資源與環(huán)境學(xué)報(bào),2021,38(6):1112-1121.

    [24] 皮青山,席 茜,周玉科.Google Earth Engine云平臺(tái)及植被遙感案例研究[J].測(cè)繪與空間地理信息,2021,44(2):38-41.

    [25] 姜伊蘭,陳保旺,黃玉芳,等.基于Google Earth Engine和NDVI時(shí)序差異指數(shù)的作物種植區(qū)提?。跩].地球信息科學(xué)學(xué)報(bào),2021,23(5):938-947.

    [26] HAMUD A M,SHAFRI H Z M,SHAHARUM N S N. Monitoring urban expansion and land use/land cover changes in Banadir, Somalia using Google Earth Engine (GEE)[J]. IOP Conference Series: Earth and environmental science,2021,767(1):012041.

    [27] ZHA Y, GAO J, NI S. Use of normalized difference built-up index in automatically mapping urban areas from TM imagery[J]. International journal of remote sensing,2003,24(3):583-594.

    [28] 徐涵秋.基于壓縮數(shù)據(jù)維的城市建筑用地遙感信息提?。跩].中國(guó)圖象圖形學(xué)報(bào),2005,10(2):223-229.

    [29] 吳志杰,趙書河.基于TM圖像的“增強(qiáng)的指數(shù)型建筑用地指數(shù)”研究[J].國(guó)土資源遙感,2012(2):50-55.

    [30] 代穎懿,劉 輝,黎啟燃,等. 基于中等空間分辨率遙感影像的建筑用地信息提?。跩].亞熱帶資源與環(huán)境學(xué)報(bào),2013,8(4):68-74.

    收稿日期:2022-09-13

    基金項(xiàng)目:國(guó)家自然科學(xué)基金項(xiàng)目(41867040)

    作者簡(jiǎn)介:鄭舒元(1997-),男,云南玉溪人,碩士,主要從事資源與環(huán)境遙感研究,(電話)15261805767(電子信箱)952012141@qq.com;通信作者,王建雄(1975-),男,云南玉溪人,教授,博士,主要從事資源與環(huán)境遙感研究,(電話)13608873306(電子信箱)jianxiongw@126.com。

    女警被强在线播放| 亚洲专区中文字幕在线| 国产老妇伦熟女老妇高清| 国产精品久久久久久精品电影小说| 国产精品久久久久成人av| 中文乱码字字幕精品一区二区三区| 99九九在线精品视频| 侵犯人妻中文字幕一二三四区| 精品国产超薄肉色丝袜足j| 秋霞在线观看毛片| a 毛片基地| 丝袜脚勾引网站| 亚洲欧美中文字幕日韩二区| 久久人人爽人人片av| 精品国产超薄肉色丝袜足j| 999久久久国产精品视频| 久久人人爽人人片av| 午夜视频精品福利| 国产片特级美女逼逼视频| 91精品国产国语对白视频| 久久99热这里只频精品6学生| 青青草视频在线视频观看| 久久人人爽av亚洲精品天堂| 人人澡人人妻人| 高潮久久久久久久久久久不卡| 免费一级毛片在线播放高清视频 | 悠悠久久av| 黑人巨大精品欧美一区二区蜜桃| 免费少妇av软件| 国产av精品麻豆| 日本wwww免费看| 国产99久久九九免费精品| 亚洲精品第二区| 精品一区二区三卡| 亚洲国产欧美一区二区综合| 久久午夜综合久久蜜桃| 婷婷成人精品国产| 女性生殖器流出的白浆| 成人国语在线视频| 在线观看免费高清a一片| 天天躁夜夜躁狠狠躁躁| 成在线人永久免费视频| 大码成人一级视频| 97人妻天天添夜夜摸| 男女无遮挡免费网站观看| 亚洲欧美一区二区三区黑人| 午夜久久久在线观看| 五月天丁香电影| av在线app专区| 国产三级黄色录像| 成人黄色视频免费在线看| 婷婷丁香在线五月| 曰老女人黄片| 免费在线观看完整版高清| 精品福利永久在线观看| 欧美精品av麻豆av| 亚洲一区中文字幕在线| 18禁观看日本| 午夜福利免费观看在线| 国产极品粉嫩免费观看在线| av国产精品久久久久影院| 精品少妇一区二区三区视频日本电影| 黄片小视频在线播放| 亚洲黑人精品在线| 激情五月婷婷亚洲| 日韩一区二区三区影片| 久久久欧美国产精品| 高清黄色对白视频在线免费看| 日韩大片免费观看网站| 国产精品久久久久成人av| 97在线人人人人妻| 中文字幕另类日韩欧美亚洲嫩草| 精品熟女少妇八av免费久了| 亚洲,一卡二卡三卡| 超碰97精品在线观看| 18禁国产床啪视频网站| 亚洲美女黄色视频免费看| 制服人妻中文乱码| 亚洲国产欧美日韩在线播放| 99re6热这里在线精品视频| videosex国产| 国产亚洲精品久久久久5区| 制服诱惑二区| 欧美人与性动交α欧美精品济南到| 嫁个100分男人电影在线观看 | www日本在线高清视频| 国产精品99久久99久久久不卡| 美女视频免费永久观看网站| 天天影视国产精品| 久久亚洲精品不卡| 一区福利在线观看| 亚洲成人国产一区在线观看 | 性少妇av在线| 又粗又硬又长又爽又黄的视频| 两个人看的免费小视频| 国产成人精品久久二区二区免费| www日本在线高清视频| 国产精品久久久久久精品古装| 天天躁夜夜躁狠狠躁躁| 久久久精品免费免费高清| 色94色欧美一区二区| 在线观看免费日韩欧美大片| 婷婷色综合大香蕉| 亚洲色图 男人天堂 中文字幕| 这个男人来自地球电影免费观看| 亚洲精品美女久久av网站| 亚洲av美国av| 免费女性裸体啪啪无遮挡网站| 中文欧美无线码| 美女扒开内裤让男人捅视频| 男男h啪啪无遮挡| 国产无遮挡羞羞视频在线观看| 免费少妇av软件| 国产在线免费精品| 亚洲一区二区三区欧美精品| 十分钟在线观看高清视频www| 91老司机精品| 日本欧美国产在线视频| 欧美日韩黄片免| 国产精品香港三级国产av潘金莲 | 久久ye,这里只有精品| 国产在视频线精品| 久久久精品区二区三区| 丁香六月天网| 又大又爽又粗| 男男h啪啪无遮挡| 国产一区有黄有色的免费视频| 欧美精品一区二区大全| 青春草亚洲视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 狠狠婷婷综合久久久久久88av| 午夜激情久久久久久久| 国产亚洲一区二区精品| 亚洲成色77777| 日韩一卡2卡3卡4卡2021年| av不卡在线播放| 午夜福利一区二区在线看| 久久九九热精品免费| 久久 成人 亚洲| 男女午夜视频在线观看| 久久久久久久久久久久大奶| 91精品国产国语对白视频| 午夜福利视频在线观看免费| 久久99精品国语久久久| 少妇猛男粗大的猛烈进出视频| 午夜精品国产一区二区电影| 人成视频在线观看免费观看| 日日夜夜操网爽| 十八禁高潮呻吟视频| 黄频高清免费视频| 久久久久久久精品精品| 成人国产一区最新在线观看 | 日韩大码丰满熟妇| 国产男女超爽视频在线观看| 老鸭窝网址在线观看| 欧美精品人与动牲交sv欧美| 亚洲欧洲国产日韩| av天堂久久9| 一边摸一边抽搐一进一出视频| 丝瓜视频免费看黄片| 亚洲精品美女久久久久99蜜臀 | 欧美日韩一级在线毛片| 99香蕉大伊视频| 国产日韩一区二区三区精品不卡| 国产片特级美女逼逼视频| 精品亚洲成a人片在线观看| www.av在线官网国产| 亚洲精品中文字幕在线视频| 超碰97精品在线观看| 男女床上黄色一级片免费看| 日韩大片免费观看网站| 老司机靠b影院| 美女脱内裤让男人舔精品视频| 亚洲国产中文字幕在线视频| 好男人电影高清在线观看| 国产福利在线免费观看视频| 国产精品一区二区免费欧美 | 一二三四社区在线视频社区8| av片东京热男人的天堂| 国产亚洲一区二区精品| 人体艺术视频欧美日本| 搡老乐熟女国产| 成年av动漫网址| 亚洲七黄色美女视频| 久久午夜综合久久蜜桃| 国产日韩欧美在线精品| 精品熟女少妇八av免费久了| 国产老妇伦熟女老妇高清| 国产精品二区激情视频| 国产精品三级大全| 人人妻,人人澡人人爽秒播 | 亚洲欧美中文字幕日韩二区| 波野结衣二区三区在线| 国产成人av激情在线播放| 97在线人人人人妻| 久9热在线精品视频| 老司机深夜福利视频在线观看 | 久久精品亚洲av国产电影网| 男人爽女人下面视频在线观看| 人人妻人人爽人人添夜夜欢视频| 色综合欧美亚洲国产小说| 欧美亚洲日本最大视频资源| 亚洲欧美精品自产自拍| 大片免费播放器 马上看| 考比视频在线观看| 一区二区av电影网| 亚洲欧美中文字幕日韩二区| 黄色一级大片看看| 一级片'在线观看视频| 妹子高潮喷水视频| 永久免费av网站大全| 校园人妻丝袜中文字幕| 午夜免费男女啪啪视频观看| 十八禁人妻一区二区| 99久久99久久久精品蜜桃| 日韩精品免费视频一区二区三区| 亚洲av成人精品一二三区| 久久久久国产精品人妻一区二区| 美女高潮到喷水免费观看| 亚洲精品国产一区二区精华液| 亚洲国产成人一精品久久久| 只有这里有精品99| 中文字幕精品免费在线观看视频| 亚洲精品美女久久久久99蜜臀 | 亚洲精品一二三| av片东京热男人的天堂| 亚洲国产成人一精品久久久| 亚洲人成电影观看| 丝袜脚勾引网站| 视频区欧美日本亚洲| 一级毛片电影观看| 在线观看一区二区三区激情| 女警被强在线播放| 亚洲av综合色区一区| 亚洲欧美一区二区三区黑人| 99久久99久久久精品蜜桃| 欧美老熟妇乱子伦牲交| 精品欧美一区二区三区在线| 亚洲av成人不卡在线观看播放网 | 夫妻午夜视频| 成人国产一区最新在线观看 | 国产欧美日韩一区二区三区在线| 满18在线观看网站| 亚洲情色 制服丝袜| 亚洲国产精品国产精品| 99香蕉大伊视频| 国产精品免费视频内射| 亚洲人成77777在线视频| 伦理电影免费视频| 国产真人三级小视频在线观看| 制服诱惑二区| 午夜福利影视在线免费观看| www.999成人在线观看| 美女扒开内裤让男人捅视频| 亚洲少妇的诱惑av| 欧美日韩福利视频一区二区| 韩国高清视频一区二区三区| 国产成人精品久久二区二区91| 国产av国产精品国产| 午夜福利乱码中文字幕| 色综合欧美亚洲国产小说| 亚洲国产成人一精品久久久| 亚洲美女黄色视频免费看| av不卡在线播放| 国产精品九九99| 欧美少妇被猛烈插入视频| 亚洲精品一区蜜桃| 最黄视频免费看| 亚洲专区国产一区二区| 国产福利在线免费观看视频| 免费在线观看视频国产中文字幕亚洲 | 黑人巨大精品欧美一区二区蜜桃| 免费观看a级毛片全部| www日本在线高清视频| 90打野战视频偷拍视频| 免费在线观看日本一区| 考比视频在线观看| 欧美中文综合在线视频| 国产亚洲精品第一综合不卡| 亚洲国产av影院在线观看| 热re99久久精品国产66热6| 高清黄色对白视频在线免费看| 亚洲av片天天在线观看| videosex国产| 真人做人爱边吃奶动态| 欧美人与善性xxx| svipshipincom国产片| 精品亚洲成a人片在线观看| 欧美日韩一级在线毛片| 悠悠久久av| 日韩大码丰满熟妇| 精品少妇一区二区三区视频日本电影| 母亲3免费完整高清在线观看| 丰满少妇做爰视频| 黑丝袜美女国产一区| 一级毛片女人18水好多 | 亚洲免费av在线视频| 校园人妻丝袜中文字幕| 99国产精品一区二区三区| 亚洲图色成人| 国产av一区二区精品久久| av网站免费在线观看视频| 亚洲国产中文字幕在线视频| 一级a爱视频在线免费观看| 99国产精品免费福利视频| 日本wwww免费看| 国产老妇伦熟女老妇高清| 国产91精品成人一区二区三区 | 久久狼人影院| 欧美日韩成人在线一区二区| 久久国产亚洲av麻豆专区| 无遮挡黄片免费观看| 亚洲精品日韩在线中文字幕| 亚洲黑人精品在线| av在线播放精品| 免费黄频网站在线观看国产| 亚洲av国产av综合av卡| 中文乱码字字幕精品一区二区三区| 国产高清不卡午夜福利| 中文字幕高清在线视频| videosex国产| 99久久综合免费| 久久久久精品人妻al黑| 久久精品久久精品一区二区三区| 欧美人与善性xxx| 搡老乐熟女国产| 人妻 亚洲 视频| 自拍欧美九色日韩亚洲蝌蚪91| 新久久久久国产一级毛片| bbb黄色大片| 超碰97精品在线观看| 十八禁人妻一区二区| 天天躁日日躁夜夜躁夜夜| 久久亚洲精品不卡| 麻豆国产av国片精品| 一本久久精品| 亚洲成人国产一区在线观看 | 手机成人av网站| 天天躁狠狠躁夜夜躁狠狠躁| 狠狠婷婷综合久久久久久88av| www.熟女人妻精品国产| 亚洲熟女毛片儿| 欧美黄色片欧美黄色片| 国产精品 国内视频| 久久久久久久久免费视频了| 9热在线视频观看99| 国产精品久久久av美女十八| 国产福利在线免费观看视频| 天堂8中文在线网| 无遮挡黄片免费观看| 两性夫妻黄色片| 亚洲自偷自拍图片 自拍| 中文字幕人妻丝袜一区二区| 日韩一卡2卡3卡4卡2021年| 亚洲专区国产一区二区| 亚洲av电影在线进入| 男女边吃奶边做爰视频| 天天躁夜夜躁狠狠躁躁| 天天躁夜夜躁狠狠久久av| 国产成人a∨麻豆精品| 两个人免费观看高清视频| 日本一区二区免费在线视频| 亚洲精品日本国产第一区| 亚洲成av片中文字幕在线观看| 中文精品一卡2卡3卡4更新| 国产高清视频在线播放一区 | 成人黄色视频免费在线看| 亚洲熟女毛片儿| 嫁个100分男人电影在线观看 | 国产麻豆69| 久久人人97超碰香蕉20202| 69精品国产乱码久久久| 纯流量卡能插随身wifi吗| 美女主播在线视频| 天天躁夜夜躁狠狠躁躁| 色婷婷av一区二区三区视频| 在线观看人妻少妇| 婷婷丁香在线五月| 两个人免费观看高清视频| svipshipincom国产片| av视频免费观看在线观看| 欧美另类一区| kizo精华| 免费在线观看日本一区| 国产成人精品在线电影| 一级毛片我不卡| 国产熟女欧美一区二区| 中文乱码字字幕精品一区二区三区| 国产成人精品无人区| 色婷婷久久久亚洲欧美| 电影成人av| 亚洲欧美精品综合一区二区三区| 丝袜人妻中文字幕| 老汉色∧v一级毛片| 如日韩欧美国产精品一区二区三区| 啦啦啦啦在线视频资源| 亚洲激情五月婷婷啪啪| 免费日韩欧美在线观看| 亚洲欧美精品自产自拍| 九色亚洲精品在线播放| 久久精品国产亚洲av高清一级| av在线播放精品| 九草在线视频观看| 咕卡用的链子| 99国产精品一区二区蜜桃av | 如日韩欧美国产精品一区二区三区| 亚洲精品中文字幕在线视频| 亚洲av成人精品一二三区| 亚洲少妇的诱惑av| 一级片'在线观看视频| 久久精品aⅴ一区二区三区四区| 国产精品三级大全| 人妻一区二区av| av国产精品久久久久影院| 女性生殖器流出的白浆| 中文字幕精品免费在线观看视频| 满18在线观看网站| 欧美成人午夜精品| 99香蕉大伊视频| 亚洲av成人不卡在线观看播放网 | 我要看黄色一级片免费的| 在线看a的网站| 国产精品99久久99久久久不卡| 最新的欧美精品一区二区| 韩国高清视频一区二区三区| bbb黄色大片| 精品亚洲成a人片在线观看| 大码成人一级视频| 亚洲视频免费观看视频| 亚洲,一卡二卡三卡| 波野结衣二区三区在线| 色婷婷久久久亚洲欧美| 男女边摸边吃奶| 精品国产超薄肉色丝袜足j| 精品熟女少妇八av免费久了| 亚洲成人免费av在线播放| 久久久国产一区二区| 久久久欧美国产精品| 我的亚洲天堂| 午夜福利免费观看在线| 老司机影院成人| 电影成人av| 在线观看www视频免费| 最黄视频免费看| 午夜免费成人在线视频| www.精华液| kizo精华| 在线观看免费高清a一片| 色综合欧美亚洲国产小说| 亚洲精品久久成人aⅴ小说| 亚洲欧洲精品一区二区精品久久久| 国产伦理片在线播放av一区| 美女福利国产在线| 男女免费视频国产| 黄网站色视频无遮挡免费观看| 欧美人与性动交α欧美精品济南到| 99久久人妻综合| 色综合欧美亚洲国产小说| 国产又色又爽无遮挡免| 黄色毛片三级朝国网站| 亚洲欧美精品自产自拍| 黄网站色视频无遮挡免费观看| 99精品久久久久人妻精品| 精品亚洲成a人片在线观看| 美女扒开内裤让男人捅视频| 在线精品无人区一区二区三| 精品卡一卡二卡四卡免费| 免费观看av网站的网址| 亚洲精品第二区| 中文字幕人妻熟女乱码| 久久久国产欧美日韩av| 亚洲精品国产色婷婷电影| www.999成人在线观看| 99九九在线精品视频| 交换朋友夫妻互换小说| 老司机靠b影院| 桃花免费在线播放| 午夜福利一区二区在线看| 看免费av毛片| 免费观看人在逋| 水蜜桃什么品种好| 亚洲天堂av无毛| 人人妻,人人澡人人爽秒播 | 波多野结衣av一区二区av| 欧美另类一区| 五月开心婷婷网| 在现免费观看毛片| 国产色视频综合| 亚洲专区中文字幕在线| 50天的宝宝边吃奶边哭怎么回事| avwww免费| 免费久久久久久久精品成人欧美视频| 美女中出高潮动态图| xxx大片免费视频| tube8黄色片| 一级毛片 在线播放| 伊人久久大香线蕉亚洲五| 欧美人与性动交α欧美精品济南到| 日韩中文字幕欧美一区二区 | 国产精品免费视频内射| 老鸭窝网址在线观看| 欧美日韩av久久| 永久免费av网站大全| 国产老妇伦熟女老妇高清| 久久99一区二区三区| 我要看黄色一级片免费的| 免费观看av网站的网址| 又黄又粗又硬又大视频| 你懂的网址亚洲精品在线观看| 欧美乱码精品一区二区三区| 成在线人永久免费视频| 亚洲国产最新在线播放| 久久女婷五月综合色啪小说| 午夜福利视频精品| 国产精品一区二区在线观看99| 亚洲av片天天在线观看| 99国产综合亚洲精品| 国产在视频线精品| 大话2 男鬼变身卡| 色婷婷久久久亚洲欧美| 亚洲一码二码三码区别大吗| 一边摸一边做爽爽视频免费| 精品久久久久久电影网| 精品人妻在线不人妻| 亚洲 欧美一区二区三区| 日本a在线网址| 婷婷色av中文字幕| 爱豆传媒免费全集在线观看| 亚洲,欧美,日韩| 一本大道久久a久久精品| 国产黄频视频在线观看| 精品人妻一区二区三区麻豆| 国产精品一区二区在线不卡| 最近中文字幕2019免费版| 欧美成人午夜精品| 一级毛片 在线播放| a级片在线免费高清观看视频| 久久精品国产综合久久久| 日韩av在线免费看完整版不卡| 纯流量卡能插随身wifi吗| 日韩av免费高清视频| 看免费成人av毛片| 9热在线视频观看99| 操出白浆在线播放| 久热这里只有精品99| 99国产精品99久久久久| 麻豆av在线久日| 激情视频va一区二区三区| 欧美日韩视频精品一区| 久久久精品区二区三区| 国产精品成人在线| 久久久久久人人人人人| 亚洲欧美成人综合另类久久久| 国产深夜福利视频在线观看| 精品久久久久久电影网| 伊人久久大香线蕉亚洲五| 又大又黄又爽视频免费| 久久久久久久大尺度免费视频| 麻豆国产av国片精品| 成年av动漫网址| 女人精品久久久久毛片| 国产日韩欧美在线精品| 国产亚洲欧美精品永久| 午夜福利一区二区在线看| 亚洲av电影在线观看一区二区三区| 人人妻人人爽人人添夜夜欢视频| 日日夜夜操网爽| 精品国产一区二区三区久久久樱花| 超碰成人久久| 色精品久久人妻99蜜桃| 热re99久久精品国产66热6| 国产熟女午夜一区二区三区| 亚洲图色成人| 日韩中文字幕视频在线看片| 久久亚洲精品不卡| 久久热在线av| 日韩熟女老妇一区二区性免费视频| 咕卡用的链子| 纵有疾风起免费观看全集完整版| 亚洲欧美中文字幕日韩二区| 91国产中文字幕| 亚洲精品国产av蜜桃| 国产精品人妻久久久影院| 久久精品aⅴ一区二区三区四区| 亚洲国产精品一区三区| 亚洲欧美一区二区三区黑人| 午夜精品国产一区二区电影| av一本久久久久| 精品国产乱码久久久久久男人| 国产一级毛片在线| 国产精品欧美亚洲77777| 大陆偷拍与自拍| 999精品在线视频| 国产精品三级大全| 国产av精品麻豆| 欧美亚洲日本最大视频资源| 国产精品三级大全| 欧美日韩综合久久久久久| 亚洲精品国产区一区二| 一边亲一边摸免费视频| 色94色欧美一区二区| 一区二区三区激情视频| av又黄又爽大尺度在线免费看| 女人久久www免费人成看片| 亚洲欧美激情在线| 国产一区亚洲一区在线观看| 在线天堂中文资源库| 精品少妇一区二区三区视频日本电影| 一本—道久久a久久精品蜜桃钙片| 日韩,欧美,国产一区二区三区| 久久久久国产精品人妻一区二区| 成人国语在线视频| 大香蕉久久网| 一区二区三区激情视频| 美女扒开内裤让男人捅视频|