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

    基于MERIS遙感圖像的太湖葉綠素濃度反演

    2017-05-30 10:10:08王根深王得玉
    安徽農(nóng)業(yè)科學(xué) 2017年30期
    關(guān)鍵詞:太湖反射率波段

    王根深 王得玉

    摘要以太湖為試驗(yàn)區(qū),基于MERIS遙感圖像數(shù)據(jù)以及同步實(shí)測的太湖水質(zhì)參數(shù)數(shù)據(jù),應(yīng)用歸一化葉綠素指數(shù)算法(NDCI),對太湖水體葉綠素a濃度進(jìn)行反演,得到了太湖區(qū)域的水體葉綠素a反演結(jié)果,并對反演結(jié)果進(jìn)行了驗(yàn)證和分析。結(jié)果表明:歸一化葉綠素指數(shù)反演算法能夠精確地反演太湖區(qū)域的葉綠素a濃度值,模型的決定系數(shù)(R2)為0.881 2,反演精度優(yōu)于傳統(tǒng)經(jīng)驗(yàn)統(tǒng)計(jì)模型,可為今后更精確地反演內(nèi)陸水體的葉綠素a濃度提供參考依據(jù)。

    關(guān)鍵詞MERIS遙感圖像;葉綠素a;歸一化葉綠素指數(shù)反演算法;太湖

    中圖分類號X87文獻(xiàn)標(biāo)識(shí)碼

    A文章編號0517-6611(2017)30-0071-04

    AbstractA novel index,normalized difference chlorophyll index(NDCI) was introduced in this study to invert the chlorophyll a(Chl.a) concentration form MERIS remote sensing image and the data of water quality parameters of Taihu Lake,and the inversion results were verified and analyzed.The inversion results showed that normalized Difference chlorophyll Index could accurately reflect the chlorophyll a concentration in the Taihu Lake, the coefficient of decision was 0.881 2, which was better than traditional empirical model,and provided a reference for the more accurate inversion of chlorophyll a concentration in inland water.

    Key wordsMERIS remote sensing image;Chlorophylla;NDCI;Taihu Lake

    內(nèi)陸水體,特別是位于經(jīng)濟(jì)發(fā)達(dá)、人口密集地區(qū)的湖泊和河口,一般受到人類活動(dòng)影響較大,湖泊水污染和水體富營養(yǎng)程度日益加重。葉綠素a是重要的水色參數(shù)之一,也是評價(jià)水體富營養(yǎng)化程度的重要指標(biāo)[1]。通過測定葉綠素a濃度,能夠直接反映水體的營養(yǎng)化程度,這是目前水環(huán)境監(jiān)測中常用并行之有效的方法之一。

    遙感技術(shù)的出現(xiàn)和應(yīng)用為水體污染檢測和研究開辟了新的途徑。利用遙感技術(shù)監(jiān)測水體,可快速獲得較為準(zhǔn)確的大面積水體信息及空間分布特征。

    目前,利用遙感技術(shù)估算葉綠素a濃度主要有3種方法:經(jīng)驗(yàn)法、半分析法和分析法。經(jīng)驗(yàn)方法主要對大量的地面實(shí)測數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析,應(yīng)用簡便,運(yùn)算快捷。楊偉等[2]驗(yàn)證了TM二三四波段反演葉綠素a濃度效果較好,但經(jīng)驗(yàn)方法由于缺乏一定的物理基礎(chǔ),過多依賴于實(shí)測數(shù)據(jù),有較大的時(shí)空限制。分析法以光在水體中的傳輸機(jī)理為理論基礎(chǔ)。李云梅等[3]利用Gordon模型對太湖的水體反射波譜進(jìn)行模擬,并反演水體的葉綠素濃度。但是分析法需要大量準(zhǔn)同步的氣象和水體固有光學(xué)特性等數(shù)據(jù),而這些數(shù)據(jù)的獲取又很困難,因此,分析法很難在實(shí)際中得到應(yīng)用。目前應(yīng)用最廣泛的是半分析法。半分析法通過分析水體組分與固有光學(xué)量、固有光學(xué)量與表觀光學(xué)量之間的關(guān)系來反演水質(zhì)參數(shù)濃度,具有很強(qiáng)的實(shí)用性。Gitelson等[4]利用MERIS數(shù)據(jù)和MODIS數(shù)據(jù)驗(yàn)證了兩波段法、三波段法在渾濁水體中的估算精度;Le等[5]利用半分析法估算了太湖的葉綠素a濃度。

    筆者以太湖為研究區(qū)域,利用MERIS遙感數(shù)據(jù)及水面實(shí)測數(shù)據(jù),對太湖葉綠素a濃度的估算方法進(jìn)一步深入研究,引入一種估算精度高、計(jì)算簡單、適用性強(qiáng)的葉綠素a濃度估算方法——?dú)w一化葉綠素指數(shù)法(NDCI),并對該方法進(jìn)行精度評價(jià)。

    1研究區(qū)概況

    研究區(qū)選擇位于長江三角洲的太湖。太湖是我國第三大淡水湖泊,地處長江三角洲,流域范圍包括江蘇、上海、浙江、安徽等地區(qū),水資源總量豐富,地勢平坦,土壤肥沃,氣候溫?zé)?,水陸交通方便,是著名的魚米之鄉(xiāng)、旅游勝地,太湖為當(dāng)?shù)亟?jīng)濟(jì)發(fā)展提供了良好的條件。太湖水面面積約2 338.1 km2,位于119°52′32″~120°36′10″ E,30°55′40″~31°32′58″ N,區(qū)域位置見圖1。近年來,隨著太湖地區(qū)經(jīng)濟(jì)發(fā)展,水體污染嚴(yán)重程度增加,水體富營養(yǎng)化程度日益加重。

    2試驗(yàn)數(shù)據(jù)及預(yù)處理

    2.1水面實(shí)測數(shù)據(jù)

    2004年10月20—29日,在太湖67個(gè)測點(diǎn)(圖2)進(jìn)行了野外測量,測點(diǎn)均勻分布于整個(gè)太湖水域,涵蓋了從極度到中度富營養(yǎng)化水體,從渾濁到清澈水體。首先根據(jù)NASA SeaWiFS海洋光學(xué)規(guī)范,測量了水體反射率光譜和水體的后向散射系數(shù)。接著在該采樣點(diǎn)取水樣,帶回實(shí)驗(yàn)室進(jìn)行相應(yīng)水質(zhì)參數(shù)和水色組分吸收系數(shù)的測量。

    2.2衛(wèi)星數(shù)據(jù)獲取

    內(nèi)陸水體光學(xué)特性復(fù)雜,對傳感器提出了較高要求,需要較高光譜分辨率、空間分辨率和時(shí)間分辨率等[6]。MERIS是搭載在ENVISAT-1上的中等分辨率成像頻譜儀,是目前最具優(yōu)勢的水色傳感器之一,主要用于海洋海岸帶的水色監(jiān)測,具有高時(shí)間、高空間、高光譜分辨率及高靈敏度的優(yōu)點(diǎn)。MERIS傳感器在光譜分辨率、水色波段設(shè)置以及輻射靈敏度等方面優(yōu)于SeaWIFS、MODIS。其空間分辨率為300 m,波譜范圍是412~1 050 nm,波譜分辨率范圍為7.5~20.0 nm,有15個(gè)波段(390~1 040 nm),覆蓋可見光到近紅外,可以應(yīng)用于復(fù)雜二類水體的監(jiān)測。

    MERIS數(shù)據(jù)產(chǎn)品分為3級:Level 1b數(shù)據(jù)(經(jīng)過輻射校正后的大氣層頂輻射亮度);Level 2b數(shù)據(jù)(經(jīng)過大氣校正后的產(chǎn)品數(shù)據(jù),主要包括葉綠素濃度、各波段歸一化離水輻射率、氣溶膠輻射等產(chǎn)品信息);Level 3b數(shù)據(jù)(經(jīng)地理校正及平均后的海洋產(chǎn)品數(shù)據(jù),分辨率降低)[7]。該研究選用的MERIS遙感數(shù)據(jù)是2004年10月20日衛(wèi)星經(jīng)過太湖境內(nèi)實(shí)時(shí)數(shù)據(jù),可用于太湖葉綠素濃度的反演分析。

    2.3數(shù)據(jù)預(yù)處理

    該研究所選用的MERIS 2b遙感數(shù)據(jù)是在1b數(shù)據(jù)的基礎(chǔ)上經(jīng)過粗略幾何校正和大氣校正后的數(shù)據(jù),帶有衛(wèi)星方位角、時(shí)間信息、經(jīng)緯度信息、太陽方位角和太陽天頂角等附屬信息,圖像較為清晰,但其精度還需要進(jìn)一步校正。

    BEAM軟件是一款開源集成處理遙感柵格數(shù)據(jù)的軟件,由BEAM軟件處理器和經(jīng)過可擴(kuò)展預(yù)先定義的參數(shù)文件構(gòu)成,還由經(jīng)過幾何校正的二級產(chǎn)品子集處理器組成。該軟件對于MERIS、AATSR、ASAR、MODIS等數(shù)據(jù)格式的處理比其他遙感處理軟件高效。

    首先利用BEAM 5.0軟件,對影像進(jìn)行預(yù)處理,由于MERIS 2b產(chǎn)品已經(jīng)經(jīng)過了輻射能量轉(zhuǎn)換,即將大氣頂層輻射亮度轉(zhuǎn)換為大氣層頂?shù)姆瓷渎?,故無需進(jìn)行輻射能量轉(zhuǎn)換,可以直接進(jìn)行其他預(yù)處理流程,流程如下:

    (1)輻射校正。輻射校正就是將遙感圖像上的灰度值轉(zhuǎn)化為入瞳輻亮度。該研究采用輻射傳輸模型對數(shù)據(jù)進(jìn)行輻射校正。

    (2)幾何校正。幾何校正是指消除或改正遙感影像幾何誤差的過程。MERIS遙感數(shù)據(jù)的經(jīng)緯度信息存儲(chǔ)在Tie Point Grid中。宋瑜等[8]通過對比使用BEAM軟件中的幾何校正模塊和ENVI軟件中的Georeference MERIS對MERIS數(shù)據(jù)進(jìn)行幾何校正的結(jié)果,發(fā)現(xiàn)前者的校正精度更好。因此,該研究采用BEAM軟件進(jìn)行幾何校正,選擇一幅影像作為基準(zhǔn)圖像,將其他影像與基準(zhǔn)圖像進(jìn)行幾何配準(zhǔn),幾何校正的同時(shí)進(jìn)行投影轉(zhuǎn)換,所選擇的控制點(diǎn)應(yīng)該盡可能均勻地分布在校正區(qū)域內(nèi),且數(shù)量要合理,使相同地標(biāo)出現(xiàn)在校正后的圖像相同位置[9]。幾何校正的同時(shí)進(jìn)行重投影轉(zhuǎn)換,坐標(biāo)系為WGS84,經(jīng)緯度投影,采用最近鄰方法重采樣,選擇3次多項(xiàng)式進(jìn)行校正,目的是提高幾何校正精度。

    (3)圖像裁剪。裁剪的目的是將研究區(qū)域之外的圖像去除,按照太湖所在的經(jīng)緯度范圍進(jìn)行裁剪,使裁剪成的太湖圖像的邊界范圍變成一個(gè)矩形。按照太湖所在的經(jīng)緯度范圍進(jìn)行裁剪,設(shè)置數(shù)據(jù)格式:North latitude bound:31.615;West longitude bound:119.997;South latitude bound:30.824;East longitude bound:120.541。

    (4)掩膜處理。裁剪完成后還需進(jìn)行一定的掩膜處理,掩膜(mask)掉陸地和云信息,主要剔除被標(biāo)識(shí)為云層(cloud)和陸地(land)的像元、含有耀斑的像元(bright)、雜散光污染的像元(glint)、大氣校正算法失效的像元(invalid)、產(chǎn)品值溢出或明顯異常的像元等(duplicated),但是云層邊緣和薄云部分難以去除[10]。

    3歸一化葉綠素指數(shù)法

    3.1葉綠素的光譜特征分析

    在內(nèi)陸水體葉綠素的遙感反演中,葉綠素a的光譜特征相對較為復(fù)雜,在藍(lán)、綠、紅、近紅外波段都存在能反映其濃度高低的特征區(qū)間。對于野外實(shí)測水體光譜反射率值,利用Matlab軟件繪制出水體的反射率光譜,以便找出指示葉綠素含量的重要敏感性波段。圖3給出了部分采樣點(diǎn)的反射率光譜曲線,實(shí)測的水體反射率光譜中,剔除數(shù)據(jù)異常的17#、48#和57#樣點(diǎn),有效數(shù)據(jù)樣點(diǎn)為64個(gè)。

    在可見光范圍內(nèi),一般水體的反射率總體上比較低,不超過10%,一般為4%~5%[11],且隨著波長的增大逐漸降低,到600 nm處為2%~3%,超過800 nm,水體幾乎為全吸收體,因此在近紅外的遙感影像上,清澈的水體呈現(xiàn)黑色[12]。但是在富營養(yǎng)化水體中,水體葉綠素濃度增加,藍(lán)光波段的反射率下降,綠光波段的反射率增高,水體反射率光譜特征其結(jié)論并不完全相同,近紅外波段仍存在一定的反射率,該波段水體影像不呈現(xiàn)黑色,而是呈現(xiàn)灰色,甚至是淺灰色[13]。圖3葉綠素光譜曲線呈現(xiàn)明顯的“峰和谷特征”,特有的峰谷特征為其在遙感影像上提供了有效的信息,具體葉綠素光譜特征:①波長540~557、697~703、666~693 nm處光譜反射率的大小受葉綠素濃度變化的影響較大,對探測葉綠素具有重要作用。②波長400~500、605~620、690~700 nm處,由于葉綠素吸收作用明顯,反射率低,對葉綠素含量的變化最敏感。③波長675 nm附近的吸收谷,葉綠素在紅波波段吸收最強(qiáng),葉綠素吸收和細(xì)胞壁散射均衡,對藻類密度和葉綠素濃度的反射敏感度最低[14]。④波長700 nm附近的反射峰是由于浮游植物分子吸收光能后,受到太陽激發(fā)進(jìn)行光合作用,激發(fā)出的能量產(chǎn)生熒光化的結(jié)果,它是判斷湖泊水體是否含有葉綠素的重要光譜特征,反射峰的位置和數(shù)值是葉綠素濃度的指示,熒光反射峰近似滿足高斯正態(tài)分布[15]。⑤在波長706 nm附近,反射曲線的斜率開始發(fā)生變化,近紅外波段處水體強(qiáng)烈吸收,反射率逐步下降,該處的反射依賴于有機(jī)和無機(jī)懸浮物的濃度,對藻類色素的反應(yīng)最不敏感[16]。

    3.2歸一化葉綠素指數(shù)法

    水中葉綠素a濃度影響到水的光譜響應(yīng)。葉綠素濃度上升時(shí),藍(lán)光波段的響應(yīng)下降,綠光、紅光波段的響應(yīng)上升。并且當(dāng)葉綠素a的濃度增加到一定值時(shí),葉綠素a的診斷波段向長波段方向移動(dòng),因此在進(jìn)行水質(zhì)葉綠素a濃度檢測時(shí),核心問題在于建立水體反射率和葉綠素a濃度之間的定量關(guān)系。因此,在進(jìn)行葉綠素a濃度反演時(shí),常用方法是對葉綠素a的有效響應(yīng)波段建立最佳比值或各種組合模型。

    Mishra[17]在對比了兩波段、三波段等模型的優(yōu)缺點(diǎn)之后,提出一種估算葉綠素濃度的新指數(shù)——?dú)w一化葉綠素指數(shù),并且驗(yàn)證了該方法在渾濁 Ⅱ 類水體中的適用性。其模型公式為

    CChl.a∞[Rrs(λ2)-Rrs(λ1)]

    [Rrs(λ2)+Rrs(λ1)] (1)

    式中,CChl.a是葉綠素a濃度,Rrs(λ1)、Rrs(λ2)分別為波段λ1、λ2處的反射率。NDCI指數(shù)的構(gòu)建采用歸一化的形式,部分消除由于太陽高度角的變化、大氣輻射等因素所帶來的影響。

    λ1和λ2的選擇有所限定:λ1選在665~675 nm附近的吸收峰處,該處的吸收峰主要是葉綠素吸收所致;

    λ2位于700 nm附近的反射峰,該反射峰對水體中葉綠素a濃度含量變化敏感。假定CDOM和TSS在這兩波段處吸收近似相等,因此在選擇λ1和λ2的波段時(shí)距離不宜過遠(yuǎn)。λ1、λ2均選擇光學(xué)特性由葉綠素a主導(dǎo)的特征光譜處,也在一定程度上減少了其他水體組分的影響。

    4建模與分析

    4.1模型構(gòu)建實(shí)測數(shù)據(jù)隨機(jī)選取15個(gè)采樣點(diǎn)(表1),用于模型構(gòu)建。

    首先在BEAM軟件中選擇第7、9波段,對應(yīng)的中心波長分別為665和708 nm。然后在MERIS影像中得到對應(yīng)點(diǎn)的反射率因子,再將[R(708)-R(665)]/[R(708)+R(665)]作為自變量,實(shí)測的樣本點(diǎn)的葉綠素a濃度作為因變量,在Matlab中作圖進(jìn)行回歸分析,得到每個(gè)遙感反射率因子與葉綠素a濃度之間的回歸關(guān)系以及反演模型(圖4)。

    根據(jù)圖4,建立的反演模型如下:

    y=27 084 x2+1 440.3x+20.636 (2)

    從圖4可以看出,葉綠素a濃度與水體反射率因子之間的關(guān)系是非線性的,數(shù)據(jù)樣點(diǎn)基本均勻分布在曲線兩側(cè),NDCI與葉綠素a濃度實(shí)測值的相關(guān)性較好,R2=0.881 2。

    4.2模型驗(yàn)證

    為了驗(yàn)證歸一化葉綠素指數(shù)法的反演精度,隨機(jī)挑選9個(gè)樣點(diǎn)對該模型進(jìn)行驗(yàn)證。由表2可知,反演結(jié)果最小相對誤差為0.63%,最大相對誤差為52.33%,平均相對誤差為14.40%。

    結(jié)果表明:歸一化葉綠素指數(shù)法反演得到的太湖水域葉綠素a濃度值與對應(yīng)的葉綠素a濃度實(shí)測值大致相同,估算結(jié)果較為準(zhǔn)確,相對誤差較小,基本可滿足目前水質(zhì)參數(shù)估算的精度要求。

    4.3反演結(jié)果誤差分析部分樣點(diǎn)葉綠素a濃度與實(shí)測值相差較大,可能由以下原因引起。

    (1)湖泊水體存在較長的時(shí)間、空間差異性,不同環(huán)境、氣候、測量條件的數(shù)據(jù)之間很難保持在統(tǒng)一的影響因子水平,使得利用簡單的相關(guān)分析方法得到的結(jié)果往往具有一定的不確定性。

    (2)由于內(nèi)陸湖泊水體的葉綠素光譜特征受懸浮顆粒物、黃色物質(zhì)的影響,其光譜特征明顯高于一類水體葉綠素

    光譜特征。葉綠素通常會(huì)被懸浮顆粒物、黃色物質(zhì)所掩蓋。

    (3)實(shí)測葉綠素a濃度之間相差不大,影響光譜特征的光學(xué)物質(zhì)含量相似,導(dǎo)致水體光譜具有相似特征。

    5結(jié)語

    該研究基于太湖實(shí)測葉綠素濃度數(shù)據(jù)和MERIS遙感數(shù)據(jù),應(yīng)用歸一化葉綠素指數(shù)法,反演出太湖葉綠素濃度,并利用太湖實(shí)測數(shù)據(jù)進(jìn)行驗(yàn)證,結(jié)果顯示,歸一化葉綠素指數(shù)法是一種精度高、計(jì)算簡單、適用性強(qiáng)的葉綠素a濃度估算方法,為內(nèi)陸渾濁水體葉綠素a濃度遙感監(jiān)測的業(yè)務(wù)運(yùn)行提供了有力例證。

    參考文獻(xiàn)

    [1]

    尹球,疏小舟,徐兆安,等.湖泊水環(huán)境指標(biāo)的超光譜響應(yīng)特征分析[J].紅外與毫米波學(xué)報(bào), 2004, 23(6):427-430, 435.

    [2] 楊偉, 陳晉, 松下文經(jīng).基于生物光學(xué)模型的水體葉綠素濃度反演算法[J].光譜學(xué)與光譜分析, 2009, 29(1): 38-42.

    [3] 李云梅, 黃家柱, 韋玉春,等.用分析模型方法反演水體葉綠素的濃度[J].遙感學(xué)報(bào), 2006, 10(2):169-175.

    [4] GITELSON A A,DALLOLMO G,MOSES W, et al.A simple semianalytical model for remote estimation of chlorophylla in turbid waters: Validation[J].Remote sensing of environment, 2008, 112(9): 3582-3593.

    [5] LE C F, LI Y M, ZHA Y, et al.A fourband semianalytical model for estimating chlorophyll a in highly turbid lakes:The case of Taihu Lake, China[J].Remote sensing of environment, 2009, 113(6):1175-1182.

    [6] 姜廣甲, 周琳, 馬榮華,等.渾濁 Ⅱ 類水體葉綠素a濃度遙感反演(Ⅱ):MERIS遙感數(shù)據(jù)的應(yīng)用[J].紅外與毫米波學(xué)報(bào), 2013, 32(4): 372-378.

    [7] 王云飛.東海赤潮監(jiān)測衛(wèi)星遙感方法研究[D].青島:中國海洋大學(xué), 2009.

    [8] 宋瑜, 宋曉東, 郭照冰,等.利用MERIS產(chǎn)品數(shù)據(jù)反演太湖葉綠素a濃度研究[J].遙感信息, 2009(4): 19-24.

    [9] 馮登超, 陳剛, 肖楷樂,等.遙感圖像的幾何精校正研究[J].國外電子測量技術(shù), 2012, 31(5): 41-43.

    [10] 李莉, 尹球, 鞏彩蘭,等.太湖不同葉綠素a濃度水體熒光特征分析[J].光譜學(xué)與光譜分析,2011,31(1):36-40.

    [11] 王志輝, 易善楨.不同指數(shù)模型法在水體遙感提取中的比較研究[J].科學(xué)技術(shù)與工程, 2007, 7(4): 534-537.

    [12] 韓麗君.遙感在資源環(huán)境調(diào)查中的應(yīng)用[J].太原理工大學(xué)學(xué)報(bào)(社會(huì)科學(xué)版), 2005, 23(S1): 69-71.

    [13] 熊旭平.基于遙感的公路域生態(tài)環(huán)境質(zhì)量評價(jià)[D].長沙:長沙理工大學(xué), 2008.

    [14] 吳傳慶, 楊志峰, 王橋,等.葉綠素a濃度的動(dòng)態(tài)峰反演方法[J].湖泊科學(xué), 2009, 21(2): 223-227.

    [15] 劉英.千島湖水體水質(zhì)參數(shù)遙感及其估測模型研究[D].杭州:浙江大學(xué), 2003.

    [16] 展建敏.實(shí)測高光譜數(shù)據(jù)歸一化區(qū)間變化對太湖葉綠素a含量估算的影響[D].南京:南京師范大學(xué), 2009.

    [17] MISHRA S,MISHRA D R.Normalized difference chlorophyll index: A novel model for remote estimation of chlorophylla, concentration in turbid productive waters[J].Remote sensing of environment, 2012, 117(2): 394-406.

    猜你喜歡
    太湖反射率波段
    春日暖陽
    影響Mini LED板油墨層反射率的因素
    近岸水體異源遙感反射率產(chǎn)品的融合方法研究
    具有顏色恒常性的光譜反射率重建
    太湖思變2017
    玩具世界(2017年4期)2017-07-21 13:27:24
    化學(xué)腐蝕硅表面結(jié)構(gòu)反射率影響因素的研究*
    電子器件(2017年2期)2017-04-25 08:58:37
    太湖攬春
    寶藏(2017年2期)2017-03-20 13:16:42
    太湖
    中亞信息(2016年3期)2016-12-01 06:08:24
    M87的多波段輻射過程及其能譜擬合
    日常維護(hù)對L 波段雷達(dá)的重要性
    西藏科技(2015年4期)2015-09-26 12:12:58
    一二三四在线观看免费中文在| 国产精品久久久久成人av| av在线播放精品| 日日爽夜夜爽网站| netflix在线观看网站| 欧美日韩国产mv在线观看视频| 国语对白做爰xxxⅹ性视频网站| 国产精品久久久久久精品电影小说| 亚洲五月色婷婷综合| 久久精品久久久久久久性| 日韩 欧美 亚洲 中文字幕| 老司机影院毛片| 国产精品国产三级专区第一集| 午夜久久久在线观看| 亚洲国产欧美网| 999久久久国产精品视频| 国产精品国产av在线观看| 亚洲一码二码三码区别大吗| 爱豆传媒免费全集在线观看| 欧美日韩国产mv在线观看视频| 18禁动态无遮挡网站| 国产成人一区二区在线| 免费观看a级毛片全部| 国产一区二区 视频在线| 老司机亚洲免费影院| 亚洲av在线观看美女高潮| 欧美日本中文国产一区发布| 又大又黄又爽视频免费| 亚洲婷婷狠狠爱综合网| 久久久久久久国产电影| 久久女婷五月综合色啪小说| 国产97色在线日韩免费| 老司机亚洲免费影院| 一级爰片在线观看| 天天躁夜夜躁狠狠久久av| 国产亚洲欧美精品永久| 2021少妇久久久久久久久久久| 十八禁高潮呻吟视频| 国产精品一区二区在线观看99| 国产午夜精品一二区理论片| 久久精品国产亚洲av高清一级| 一本—道久久a久久精品蜜桃钙片| 99九九在线精品视频| 欧美成人精品欧美一级黄| 国产成人一区二区在线| 男人操女人黄网站| 最近中文字幕高清免费大全6| 免费看av在线观看网站| 久久久久久久久久久免费av| 少妇的丰满在线观看| 夜夜骑夜夜射夜夜干| 国产深夜福利视频在线观看| 最近中文字幕2019免费版| 日本vs欧美在线观看视频| 亚洲欧美一区二区三区久久| 亚洲少妇的诱惑av| 精品国产露脸久久av麻豆| 少妇人妻久久综合中文| 一本大道久久a久久精品| 91成人精品电影| av.在线天堂| 国产97色在线日韩免费| 亚洲欧美一区二区三区国产| 国产成人午夜福利电影在线观看| av免费观看日本| 啦啦啦中文免费视频观看日本| 国产一级毛片在线| 国产高清不卡午夜福利| 成人三级做爰电影| 亚洲欧美激情在线| 黄色一级大片看看| 亚洲成国产人片在线观看| 高清欧美精品videossex| 丝袜喷水一区| 最近最新中文字幕大全免费视频 | 亚洲精品乱久久久久久| 777久久人妻少妇嫩草av网站| 日韩电影二区| 日本av免费视频播放| 最近中文字幕高清免费大全6| 在线亚洲精品国产二区图片欧美| 国产伦人伦偷精品视频| 美女扒开内裤让男人捅视频| 妹子高潮喷水视频| 亚洲av电影在线进入| 亚洲欧洲国产日韩| av电影中文网址| 国产精品麻豆人妻色哟哟久久| 免费女性裸体啪啪无遮挡网站| 中文字幕色久视频| 侵犯人妻中文字幕一二三四区| 热99久久久久精品小说推荐| 在线观看免费日韩欧美大片| 美女福利国产在线| 男女无遮挡免费网站观看| 母亲3免费完整高清在线观看| √禁漫天堂资源中文www| 一二三四在线观看免费中文在| 熟女av电影| 免费观看a级毛片全部| 卡戴珊不雅视频在线播放| 自线自在国产av| 欧美日韩综合久久久久久| 色吧在线观看| 日韩 亚洲 欧美在线| 欧美黑人精品巨大| 午夜福利在线免费观看网站| 超色免费av| 色婷婷av一区二区三区视频| 国产高清国产精品国产三级| 少妇猛男粗大的猛烈进出视频| 亚洲精品国产av成人精品| 韩国高清视频一区二区三区| 午夜影院在线不卡| 精品酒店卫生间| 1024视频免费在线观看| 亚洲成国产人片在线观看| 丁香六月欧美| 亚洲精品第二区| 少妇精品久久久久久久| 看免费成人av毛片| 电影成人av| 亚洲av福利一区| 丝袜美足系列| 蜜桃国产av成人99| 免费高清在线观看日韩| 黑人猛操日本美女一级片| 国产乱来视频区| 99re6热这里在线精品视频| 欧美日韩福利视频一区二区| 美国免费a级毛片| 亚洲国产看品久久| 哪个播放器可以免费观看大片| 在线观看免费高清a一片| 国产野战对白在线观看| 又大又爽又粗| 好男人视频免费观看在线| 国产又爽黄色视频| 国产精品无大码| 国产成人a∨麻豆精品| 91成人精品电影| 中文天堂在线官网| 91精品国产国语对白视频| videos熟女内射| 亚洲精品aⅴ在线观看| 伦理电影免费视频| 亚洲美女黄片视频| 黄片大片在线免费观看| 麻豆一二三区av精品| 亚洲欧美日韩无卡精品| av超薄肉色丝袜交足视频| 高潮久久久久久久久久久不卡| 久久婷婷成人综合色麻豆| 亚洲久久久国产精品| av免费在线观看网站| 黄色a级毛片大全视频| 亚洲午夜理论影院| 日韩精品中文字幕看吧| 在线观看66精品国产| 麻豆久久精品国产亚洲av| 亚洲欧美日韩无卡精品| 亚洲,欧美精品.| 男女之事视频高清在线观看| 麻豆国产av国片精品| 国产成人av教育| 精品久久久久久久久久免费视频| 亚洲av五月六月丁香网| 国产精品影院久久| 亚洲精品国产精品久久久不卡| 成人手机av| 亚洲精品国产色婷婷电影| 女性被躁到高潮视频| 国产伦人伦偷精品视频| 后天国语完整版免费观看| 午夜福利成人在线免费观看| 国产精品一区二区精品视频观看| 国产三级黄色录像| 日韩欧美在线二视频| 每晚都被弄得嗷嗷叫到高潮| av网站免费在线观看视频| 人人澡人人妻人| 人人妻人人澡欧美一区二区 | 在线国产一区二区在线| 久久婷婷人人爽人人干人人爱 | 日韩成人在线观看一区二区三区| 一区在线观看完整版| 国产免费男女视频| 长腿黑丝高跟| 法律面前人人平等表现在哪些方面| 亚洲成人久久性| 制服丝袜大香蕉在线| 欧美大码av| 日本五十路高清| 久久天躁狠狠躁夜夜2o2o| 黄色毛片三级朝国网站| 一进一出好大好爽视频| 精品欧美一区二区三区在线| 国产欧美日韩一区二区三区在线| 亚洲aⅴ乱码一区二区在线播放 | 可以在线观看毛片的网站| 亚洲少妇的诱惑av| 久久精品91无色码中文字幕| 中文字幕精品免费在线观看视频| 欧美激情极品国产一区二区三区| 亚洲av电影不卡..在线观看| 亚洲国产欧美日韩在线播放| 国产精品久久视频播放| 亚洲自偷自拍图片 自拍| 欧美激情 高清一区二区三区| 麻豆国产av国片精品| 国产成人免费无遮挡视频| 亚洲国产高清在线一区二区三 | 黄网站色视频无遮挡免费观看| 老汉色av国产亚洲站长工具| 国产av一区在线观看免费| 国产精品影院久久| 久久香蕉激情| 亚洲国产高清在线一区二区三 | 男人舔女人下体高潮全视频| 女生性感内裤真人,穿戴方法视频| 日韩大码丰满熟妇| 午夜久久久久精精品| 激情视频va一区二区三区| 亚洲男人天堂网一区| 香蕉国产在线看| 午夜免费鲁丝| 999精品在线视频| 国产黄a三级三级三级人| 国产极品粉嫩免费观看在线| 亚洲 欧美一区二区三区| 精品电影一区二区在线| av有码第一页| 国产麻豆成人av免费视频| av电影中文网址| 午夜a级毛片| 免费人成视频x8x8入口观看| 我的亚洲天堂| 日韩欧美一区二区三区在线观看| 两人在一起打扑克的视频| 久久香蕉精品热| 欧美日韩亚洲综合一区二区三区_| xxx96com| 一进一出好大好爽视频| 亚洲国产中文字幕在线视频| 欧美av亚洲av综合av国产av| 国产成人av教育| 女同久久另类99精品国产91| 精品免费久久久久久久清纯| 黄色毛片三级朝国网站| 国产精品二区激情视频| 色尼玛亚洲综合影院| 神马国产精品三级电影在线观看 | 亚洲av第一区精品v没综合| 精品国内亚洲2022精品成人| 亚洲第一青青草原| 自拍欧美九色日韩亚洲蝌蚪91| 成人av一区二区三区在线看| 日本 欧美在线| 亚洲一区二区三区不卡视频| 19禁男女啪啪无遮挡网站| 丰满人妻熟妇乱又伦精品不卡| 人人妻,人人澡人人爽秒播| 国产精华一区二区三区| 国产精品免费视频内射| 日韩欧美三级三区| 国产麻豆69| 午夜免费鲁丝| 亚洲电影在线观看av| 欧美黑人欧美精品刺激| 夜夜爽天天搞| 伦理电影免费视频| 悠悠久久av| 亚洲av日韩精品久久久久久密| 欧美丝袜亚洲另类 | netflix在线观看网站| 精品久久久久久,| 一个人观看的视频www高清免费观看 | 国产精品久久久久久人妻精品电影| 最新美女视频免费是黄的| 精品久久久久久,| 国产麻豆成人av免费视频| 亚洲 欧美 日韩 在线 免费| 一级,二级,三级黄色视频| 欧美绝顶高潮抽搐喷水| 香蕉丝袜av| 99国产精品99久久久久| 日韩三级视频一区二区三区| 怎么达到女性高潮| 亚洲av美国av| 精品第一国产精品| 国产91精品成人一区二区三区| 精品午夜福利视频在线观看一区| av天堂在线播放| 午夜精品久久久久久毛片777| 亚洲色图av天堂| 俄罗斯特黄特色一大片| 一级毛片女人18水好多| 亚洲av美国av| www.熟女人妻精品国产| 久久精品aⅴ一区二区三区四区| e午夜精品久久久久久久| 在线观看66精品国产| 亚洲av熟女| 高潮久久久久久久久久久不卡| 国产精品久久久久久精品电影 | 长腿黑丝高跟| 国产一卡二卡三卡精品| 亚洲七黄色美女视频| 操美女的视频在线观看| 看黄色毛片网站| 免费一级毛片在线播放高清视频 | 国产成人精品久久二区二区91| 欧美中文综合在线视频| 久久婷婷成人综合色麻豆| 99国产综合亚洲精品| 成人18禁在线播放| 午夜精品国产一区二区电影| 久久精品国产亚洲av高清一级| 精品第一国产精品| 精品人妻1区二区| 久久天躁狠狠躁夜夜2o2o| 亚洲中文日韩欧美视频| 国产蜜桃级精品一区二区三区| 无人区码免费观看不卡| 亚洲欧美精品综合一区二区三区| 国产精品二区激情视频| e午夜精品久久久久久久| 天天躁夜夜躁狠狠躁躁| 麻豆国产av国片精品| 日韩中文字幕欧美一区二区| 国产精华一区二区三区| 身体一侧抽搐| 精品熟女少妇八av免费久了| 色哟哟哟哟哟哟| 一区二区三区精品91| 欧美日韩精品网址| av视频在线观看入口| 久久人人爽av亚洲精品天堂| 国产一区在线观看成人免费| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲中文字幕一区二区三区有码在线看 | 久久久久久亚洲精品国产蜜桃av| 色在线成人网| 精品久久蜜臀av无| 亚洲国产精品sss在线观看| e午夜精品久久久久久久| 婷婷精品国产亚洲av在线| bbb黄色大片| 纯流量卡能插随身wifi吗| 99精品久久久久人妻精品| 亚洲国产欧美日韩在线播放| 女警被强在线播放| 婷婷精品国产亚洲av在线| 三级毛片av免费| 91字幕亚洲| bbb黄色大片| 91字幕亚洲| 国产成人av教育| 一区二区三区高清视频在线| 操出白浆在线播放| 好男人电影高清在线观看| 女警被强在线播放| 亚洲精品一区av在线观看| 三级毛片av免费| 黑人巨大精品欧美一区二区蜜桃| 国内精品久久久久久久电影| 天天躁狠狠躁夜夜躁狠狠躁| 在线av久久热| 久久国产乱子伦精品免费另类| 欧美黄色淫秽网站| 搡老岳熟女国产| 自拍欧美九色日韩亚洲蝌蚪91| 国产xxxxx性猛交| 大码成人一级视频| 男女做爰动态图高潮gif福利片 | 日本五十路高清| 精品无人区乱码1区二区| 亚洲全国av大片| tocl精华| 精品电影一区二区在线| 亚洲成a人片在线一区二区| 91成人精品电影| 极品人妻少妇av视频| 18禁观看日本| 在线av久久热| 精品一品国产午夜福利视频| 久久影院123| 在线免费观看的www视频| 黄片播放在线免费| 久久久久久久午夜电影| 91九色精品人成在线观看| 欧美激情高清一区二区三区| 国产精品国产高清国产av| 电影成人av| 在线观看免费视频日本深夜| 久久九九热精品免费| 免费在线观看影片大全网站| 亚洲人成伊人成综合网2020| 一级a爱片免费观看的视频| 校园春色视频在线观看| 国产私拍福利视频在线观看| 亚洲色图综合在线观看| 久久久精品欧美日韩精品| 波多野结衣巨乳人妻| 国产av精品麻豆| 两性午夜刺激爽爽歪歪视频在线观看 | 精品人妻1区二区| 国产精品亚洲av一区麻豆| 搡老岳熟女国产| www.精华液| 宅男免费午夜| 99精品在免费线老司机午夜| 俄罗斯特黄特色一大片| 欧美不卡视频在线免费观看 | 日日干狠狠操夜夜爽| 国产乱人伦免费视频| 又大又爽又粗| 午夜福利高清视频| 桃色一区二区三区在线观看| 在线观看免费视频网站a站| 中文字幕最新亚洲高清| 国产欧美日韩精品亚洲av| 91精品国产国语对白视频| 99国产综合亚洲精品| www日本在线高清视频| 母亲3免费完整高清在线观看| 国产亚洲精品久久久久久毛片| 给我免费播放毛片高清在线观看| ponron亚洲| 97超级碰碰碰精品色视频在线观看| 亚洲国产精品sss在线观看| 国产成人啪精品午夜网站| 久久精品国产清高在天天线| 亚洲色图综合在线观看| 91av网站免费观看| 黄色成人免费大全| 亚洲一卡2卡3卡4卡5卡精品中文| av天堂在线播放| 婷婷丁香在线五月| 1024香蕉在线观看| 99国产精品99久久久久| 午夜亚洲福利在线播放| 黑人巨大精品欧美一区二区mp4| 欧美成狂野欧美在线观看| 97碰自拍视频| 欧美一级毛片孕妇| 国产三级在线视频| 亚洲成国产人片在线观看| 国产主播在线观看一区二区| 高潮久久久久久久久久久不卡| 国产精品爽爽va在线观看网站 | 亚洲精品国产区一区二| 97超级碰碰碰精品色视频在线观看| 国产xxxxx性猛交| 亚洲天堂国产精品一区在线| 女性生殖器流出的白浆| 超碰成人久久| 大香蕉久久成人网| 韩国精品一区二区三区| 大码成人一级视频| 亚洲男人天堂网一区| 国产激情欧美一区二区| 免费观看人在逋| 午夜激情av网站| 青草久久国产| 亚洲 国产 在线| 久久草成人影院| 18禁国产床啪视频网站| 色综合婷婷激情| 日韩av在线大香蕉| 欧美亚洲日本最大视频资源| 亚洲一区高清亚洲精品| 日韩欧美一区二区三区在线观看| 免费一级毛片在线播放高清视频 | 亚洲色图av天堂| 真人一进一出gif抽搐免费| av中文乱码字幕在线| www.精华液| 天堂√8在线中文| 成人永久免费在线观看视频| 国产成人影院久久av| 伊人久久大香线蕉亚洲五| 欧美成人免费av一区二区三区| 在线永久观看黄色视频| 日日摸夜夜添夜夜添小说| 亚洲av第一区精品v没综合| 亚洲成人免费电影在线观看| 日韩一卡2卡3卡4卡2021年| 中文字幕人成人乱码亚洲影| 制服人妻中文乱码| 黑丝袜美女国产一区| 国产精品98久久久久久宅男小说| 免费在线观看黄色视频的| 欧美av亚洲av综合av国产av| av天堂久久9| 国产精品日韩av在线免费观看 | 两个人视频免费观看高清| 动漫黄色视频在线观看| 操美女的视频在线观看| 男女下面进入的视频免费午夜 | 久久 成人 亚洲| 男女之事视频高清在线观看| 动漫黄色视频在线观看| 精品一区二区三区四区五区乱码| 日韩精品免费视频一区二区三区| 一级毛片女人18水好多| 久久精品国产亚洲av高清一级| 在线观看www视频免费| 精品无人区乱码1区二区| 日本 av在线| 50天的宝宝边吃奶边哭怎么回事| av在线天堂中文字幕| 国产视频一区二区在线看| 三级毛片av免费| 日韩精品免费视频一区二区三区| 成人欧美大片| 色播亚洲综合网| 日韩免费av在线播放| 欧美另类亚洲清纯唯美| 黄色视频不卡| 亚洲最大成人中文| 一级作爱视频免费观看| 午夜福利欧美成人| 成人欧美大片| 亚洲aⅴ乱码一区二区在线播放 | 一边摸一边抽搐一进一出视频| 欧美色视频一区免费| 搡老岳熟女国产| 一级毛片女人18水好多| 大型av网站在线播放| 99久久久亚洲精品蜜臀av| 91麻豆av在线| 国产午夜精品久久久久久| 久久久久国产精品人妻aⅴ院| 后天国语完整版免费观看| 亚洲第一av免费看| 国产精品电影一区二区三区| av网站免费在线观看视频| 天天躁夜夜躁狠狠躁躁| 丰满人妻熟妇乱又伦精品不卡| 91成年电影在线观看| 亚洲欧美日韩另类电影网站| 久99久视频精品免费| 久热这里只有精品99| 色播在线永久视频| 高清毛片免费观看视频网站| 日韩成人在线观看一区二区三区| 欧美一级a爱片免费观看看 | 叶爱在线成人免费视频播放| 激情视频va一区二区三区| www国产在线视频色| 国产欧美日韩精品亚洲av| 日韩精品中文字幕看吧| 国产麻豆成人av免费视频| 国产精品久久久久久精品电影 | 国产成人系列免费观看| 人妻丰满熟妇av一区二区三区| 高清黄色对白视频在线免费看| 欧美一区二区精品小视频在线| 欧美中文综合在线视频| 精品久久久久久久毛片微露脸| 国产精品亚洲一级av第二区| 亚洲成国产人片在线观看| 国产欧美日韩精品亚洲av| 老熟妇乱子伦视频在线观看| 精品一区二区三区四区五区乱码| 国产精品免费一区二区三区在线| 成人永久免费在线观看视频| 国产欧美日韩一区二区精品| 国产主播在线观看一区二区| 9191精品国产免费久久| 青草久久国产| av在线播放免费不卡| 色综合婷婷激情| 亚洲五月婷婷丁香| 99re在线观看精品视频| 国产精品野战在线观看| 一级,二级,三级黄色视频| 最好的美女福利视频网| 日韩大尺度精品在线看网址 | 亚洲成人久久性| 多毛熟女@视频| 叶爱在线成人免费视频播放| 国产成人欧美| 丰满人妻熟妇乱又伦精品不卡| 亚洲九九香蕉| 女人被躁到高潮嗷嗷叫费观| 午夜成年电影在线免费观看| 99国产综合亚洲精品| 老司机靠b影院| 亚洲国产中文字幕在线视频| 国产精品精品国产色婷婷| 精品一区二区三区视频在线观看免费| 色av中文字幕| 日韩视频一区二区在线观看| 亚洲中文字幕一区二区三区有码在线看 | 人人妻,人人澡人人爽秒播| 国产在线精品亚洲第一网站| 最近最新中文字幕大全免费视频| 日韩精品青青久久久久久| 黄片播放在线免费| 母亲3免费完整高清在线观看| 人人妻人人澡欧美一区二区 | а√天堂www在线а√下载| 日本免费一区二区三区高清不卡 | 黑丝袜美女国产一区| 欧美亚洲日本最大视频资源| 国内久久婷婷六月综合欲色啪| 免费看美女性在线毛片视频| 老熟妇仑乱视频hdxx| 制服人妻中文乱码| 69av精品久久久久久| 人人妻人人澡欧美一区二区 | 高清黄色对白视频在线免费看| 国产精品久久久久久亚洲av鲁大|