• <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
    国产精品 国内视频| 亚洲成人免费电影在线观看| www.熟女人妻精品国产| 成人手机av| 午夜精品在线福利| 高清黄色对白视频在线免费看| 变态另类丝袜制服| 国产精品二区激情视频| 19禁男女啪啪无遮挡网站| videosex国产| 人成视频在线观看免费观看| 69av精品久久久久久| 亚洲色图 男人天堂 中文字幕| 国内久久婷婷六月综合欲色啪| 国产免费男女视频| 欧美日本亚洲视频在线播放| 免费观看人在逋| 一级作爱视频免费观看| 日日爽夜夜爽网站| 此物有八面人人有两片| 久久精品人人爽人人爽视色| 黄色片一级片一级黄色片| 91老司机精品| 亚洲成人久久性| 久久中文看片网| 香蕉丝袜av| 日韩一卡2卡3卡4卡2021年| 99精品在免费线老司机午夜| 亚洲最大成人中文| av欧美777| 亚洲精品中文字幕一二三四区| 欧美色欧美亚洲另类二区 | 日韩精品免费视频一区二区三区| 老鸭窝网址在线观看| 老司机午夜福利在线观看视频| 9191精品国产免费久久| 日韩成人在线观看一区二区三区| 黄色视频,在线免费观看| 一级作爱视频免费观看| 丁香六月欧美| 少妇粗大呻吟视频| 国产精品久久久久久亚洲av鲁大| 黄色片一级片一级黄色片| 69精品国产乱码久久久| 一进一出抽搐动态| 国产成年人精品一区二区| 免费在线观看影片大全网站| 在线观看66精品国产| 欧美最黄视频在线播放免费| 最近最新中文字幕大全电影3 | 身体一侧抽搐| 少妇 在线观看| 久久香蕉精品热| 国产精品影院久久| 国产一级毛片七仙女欲春2 | 一边摸一边抽搐一进一出视频| 天堂影院成人在线观看| 亚洲精品国产色婷婷电影| av天堂久久9| 人妻久久中文字幕网| 国产激情久久老熟女| 亚洲男人的天堂狠狠| 亚洲色图综合在线观看| 9热在线视频观看99| 好男人电影高清在线观看| 男人舔女人下体高潮全视频| 国产一区二区激情短视频| 日日干狠狠操夜夜爽| 女人被狂操c到高潮| 搞女人的毛片| 婷婷六月久久综合丁香| 淫秽高清视频在线观看| 夜夜夜夜夜久久久久| 香蕉久久夜色| 久久亚洲精品不卡| 国产精华一区二区三区| 怎么达到女性高潮| 亚洲欧美激情在线| 少妇粗大呻吟视频| 一进一出抽搐动态| 国产午夜精品久久久久久| 国产午夜精品久久久久久| 99国产精品免费福利视频| 亚洲男人天堂网一区| 一二三四社区在线视频社区8| 99精品欧美一区二区三区四区| 97人妻精品一区二区三区麻豆 | 搞女人的毛片| 不卡一级毛片| 欧美黑人精品巨大| 国产av精品麻豆| 大型av网站在线播放| 国产成人一区二区三区免费视频网站| 国产人伦9x9x在线观看| 国产精品久久电影中文字幕| 欧美乱妇无乱码| 欧美中文日本在线观看视频| 亚洲午夜理论影院| 免费无遮挡裸体视频| 99国产精品一区二区三区| 纯流量卡能插随身wifi吗| 午夜免费观看网址| 色播亚洲综合网| 很黄的视频免费| 黄色毛片三级朝国网站| 老司机深夜福利视频在线观看| 91麻豆av在线| 亚洲av片天天在线观看| 精品国产超薄肉色丝袜足j| 欧美av亚洲av综合av国产av| 精品人妻在线不人妻| 精品国产乱码久久久久久男人| 国产午夜精品久久久久久| 一边摸一边做爽爽视频免费| 黄网站色视频无遮挡免费观看| 久久国产乱子伦精品免费另类| 99在线人妻在线中文字幕| 男人的好看免费观看在线视频 | av天堂在线播放| 欧美午夜高清在线| 性色av乱码一区二区三区2| 如日韩欧美国产精品一区二区三区| 色av中文字幕| 搡老熟女国产l中国老女人| 国产三级在线视频| 亚洲性夜色夜夜综合| 欧美色视频一区免费| 国产高清videossex| 欧美日韩瑟瑟在线播放| 神马国产精品三级电影在线观看 | 美女扒开内裤让男人捅视频| 欧美乱码精品一区二区三区| 嫩草影视91久久| 国产精品二区激情视频| 18美女黄网站色大片免费观看| 非洲黑人性xxxx精品又粗又长| 免费在线观看亚洲国产| 成人18禁高潮啪啪吃奶动态图| 成人永久免费在线观看视频| 国产精品1区2区在线观看.| 国产一级毛片七仙女欲春2 | 精品国产一区二区三区四区第35| 人人妻人人澡欧美一区二区 | 久久天躁狠狠躁夜夜2o2o| 一区二区三区精品91| 国产成人欧美在线观看| 黑人欧美特级aaaaaa片| 国产免费男女视频| 精品一区二区三区视频在线观看免费| 香蕉久久夜色| 国产区一区二久久| 国产精品 欧美亚洲| 欧美中文综合在线视频| 亚洲欧美精品综合久久99| 人妻丰满熟妇av一区二区三区| 亚洲色图av天堂| 999久久久国产精品视频| 亚洲成人国产一区在线观看| av超薄肉色丝袜交足视频| 国产单亲对白刺激| 一区二区三区激情视频| 淫秽高清视频在线观看| 久久国产精品人妻蜜桃| 97人妻精品一区二区三区麻豆 | 国内毛片毛片毛片毛片毛片| 国产激情久久老熟女| 身体一侧抽搐| 欧美成人性av电影在线观看| 一级毛片精品| 久久精品国产亚洲av香蕉五月| 亚洲一卡2卡3卡4卡5卡精品中文| 9热在线视频观看99| 国产xxxxx性猛交| 久久久水蜜桃国产精品网| 国产成人影院久久av| 免费在线观看亚洲国产| 黄色毛片三级朝国网站| 一区二区三区激情视频| 婷婷六月久久综合丁香| 免费高清在线观看日韩| 真人一进一出gif抽搐免费| 成人亚洲精品av一区二区| 欧美乱码精品一区二区三区| 日本免费一区二区三区高清不卡 | 操出白浆在线播放| 国产精品日韩av在线免费观看 | 50天的宝宝边吃奶边哭怎么回事| 香蕉久久夜色| 老汉色∧v一级毛片| 人人妻人人澡欧美一区二区 | 亚洲国产精品久久男人天堂| 搡老熟女国产l中国老女人| 日本免费一区二区三区高清不卡 | 啦啦啦韩国在线观看视频| 国产午夜福利久久久久久| 亚洲av美国av| 一本久久中文字幕| 欧美日本视频| 一二三四社区在线视频社区8| 国产成人欧美在线观看| 不卡一级毛片| 国产成人精品无人区| 精品国产乱子伦一区二区三区| 色在线成人网| 国产精华一区二区三区| 久久久久国产精品人妻aⅴ院| 亚洲成av片中文字幕在线观看| 欧美色欧美亚洲另类二区 | 久久久精品国产亚洲av高清涩受| 精品国产亚洲在线| 亚洲国产精品成人综合色| 国产激情欧美一区二区| 久久中文看片网| 国产97色在线日韩免费| 日韩免费av在线播放| 色播亚洲综合网| 日日爽夜夜爽网站| 亚洲,欧美精品.| 操美女的视频在线观看| 少妇粗大呻吟视频| www.www免费av| 国产亚洲精品一区二区www| 国产精品,欧美在线| 桃色一区二区三区在线观看| 12—13女人毛片做爰片一| cao死你这个sao货| 久久国产亚洲av麻豆专区| 久久婷婷人人爽人人干人人爱 | 男人舔女人下体高潮全视频| e午夜精品久久久久久久| 亚洲电影在线观看av| 搡老妇女老女人老熟妇| 精品国产乱子伦一区二区三区| 亚洲精品中文字幕在线视频| 亚洲视频免费观看视频| 亚洲av第一区精品v没综合| x7x7x7水蜜桃| av电影中文网址| 日韩欧美一区视频在线观看| 免费高清在线观看日韩| 免费不卡黄色视频| 欧美不卡视频在线免费观看 | 一本综合久久免费| 国产成人系列免费观看| 色精品久久人妻99蜜桃| 国产精品一区二区在线不卡| 久久婷婷成人综合色麻豆| 久久精品国产亚洲av香蕉五月| 18美女黄网站色大片免费观看| 国产激情久久老熟女| 国产精品99久久99久久久不卡| 久久精品aⅴ一区二区三区四区| 亚洲av第一区精品v没综合| 久久久国产成人精品二区| 久久精品国产综合久久久| 91字幕亚洲| 国产精品av久久久久免费| 99久久综合精品五月天人人| 色婷婷久久久亚洲欧美| 自线自在国产av| 国产三级黄色录像| 精品久久久久久成人av| 日本一区二区免费在线视频| 人妻丰满熟妇av一区二区三区| 妹子高潮喷水视频| 俄罗斯特黄特色一大片| 一区二区三区国产精品乱码| 久久精品影院6| 校园春色视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品秋霞免费鲁丝片| 欧美性长视频在线观看| 99国产综合亚洲精品| 淫妇啪啪啪对白视频| 男女下面进入的视频免费午夜 | 亚洲精品久久国产高清桃花| 一边摸一边抽搐一进一出视频| 啦啦啦 在线观看视频| 成人18禁在线播放| 国产高清激情床上av| 在线av久久热| 国产三级在线视频| 老司机午夜十八禁免费视频| 老司机靠b影院| 久久精品91无色码中文字幕| 手机成人av网站| 国产一区二区在线av高清观看| 在线av久久热| 在线十欧美十亚洲十日本专区| 91麻豆精品激情在线观看国产| 人人妻人人澡人人看| 欧美国产精品va在线观看不卡| 18禁国产床啪视频网站| 操出白浆在线播放| 久久精品亚洲精品国产色婷小说| 亚洲视频免费观看视频| 国产成人一区二区三区免费视频网站| 国产成人精品久久二区二区免费| 波多野结衣av一区二区av| 亚洲国产精品成人综合色| 黑人欧美特级aaaaaa片| 亚洲精品在线观看二区| 黄色视频,在线免费观看| 国产激情久久老熟女| 欧美日韩黄片免| 国产精品av久久久久免费| 久久久久精品国产欧美久久久| 精品国产亚洲在线| 制服丝袜大香蕉在线| 久久午夜综合久久蜜桃| 亚洲国产高清在线一区二区三 | 丝袜人妻中文字幕| 多毛熟女@视频| 伊人久久大香线蕉亚洲五| 久久久久久久精品吃奶| 国内久久婷婷六月综合欲色啪| 国产精品综合久久久久久久免费 | 午夜久久久在线观看| 久久婷婷成人综合色麻豆| 亚洲av第一区精品v没综合| 国产精品一区二区免费欧美| 伊人久久大香线蕉亚洲五| 丝袜人妻中文字幕| 亚洲中文日韩欧美视频| 久久精品aⅴ一区二区三区四区| 免费在线观看黄色视频的| 成人国产综合亚洲| 国产欧美日韩精品亚洲av| 无人区码免费观看不卡| 日韩国内少妇激情av| 亚洲美女黄片视频| 国产欧美日韩综合在线一区二区| 久久久久国产精品人妻aⅴ院| АⅤ资源中文在线天堂| 日韩成人在线观看一区二区三区| 亚洲国产欧美一区二区综合| 真人做人爱边吃奶动态| 亚洲精品久久成人aⅴ小说| 免费在线观看黄色视频的| 国产精品国产高清国产av| 日韩三级视频一区二区三区| 成熟少妇高潮喷水视频| ponron亚洲| 一个人观看的视频www高清免费观看 | 国产成人系列免费观看| 久久久久国内视频| 99久久99久久久精品蜜桃| 亚洲激情在线av| 日本撒尿小便嘘嘘汇集6| 亚洲欧美日韩另类电影网站| 国产精品日韩av在线免费观看 | 黄色女人牲交| 操出白浆在线播放| 女性被躁到高潮视频| 亚洲av成人av| 成熟少妇高潮喷水视频| 国产成人欧美| 日韩av在线大香蕉| 91老司机精品| 97人妻精品一区二区三区麻豆 | 搡老妇女老女人老熟妇| 国产精品一区二区免费欧美| 精品久久久久久,| 精品欧美国产一区二区三| 国产精品乱码一区二三区的特点 | 99精品在免费线老司机午夜| 巨乳人妻的诱惑在线观看| 69精品国产乱码久久久| 激情在线观看视频在线高清| 成人三级黄色视频| 大香蕉久久成人网| 91在线观看av| 日韩欧美三级三区| 免费在线观看亚洲国产| 国产精品一区二区精品视频观看| 宅男免费午夜| 狂野欧美激情性xxxx| 侵犯人妻中文字幕一二三四区| 中文字幕人妻丝袜一区二区| 亚洲欧美精品综合久久99| 午夜精品国产一区二区电影| 亚洲精品中文字幕一二三四区| 免费观看精品视频网站| 精品国产一区二区三区四区第35| 欧美日韩福利视频一区二区| 一级a爱视频在线免费观看| or卡值多少钱| 两个人看的免费小视频| 国产成人系列免费观看| 中文字幕人妻熟女乱码| 色老头精品视频在线观看| 99在线人妻在线中文字幕| 国产av一区在线观看免费| 日韩精品青青久久久久久| 一边摸一边抽搐一进一小说| 成人永久免费在线观看视频| 久久精品影院6| 成熟少妇高潮喷水视频| 精品国产亚洲在线| 十八禁网站免费在线| 亚洲精品国产一区二区精华液| 亚洲 国产 在线| 精品欧美国产一区二区三| 国语自产精品视频在线第100页| 亚洲狠狠婷婷综合久久图片| 九色国产91popny在线| 国产xxxxx性猛交| 久久午夜亚洲精品久久| 女性被躁到高潮视频| 免费高清在线观看日韩| 日韩大尺度精品在线看网址 | 欧美日韩亚洲国产一区二区在线观看| 久久欧美精品欧美久久欧美| 久9热在线精品视频| 亚洲午夜精品一区,二区,三区| 国产黄a三级三级三级人| av天堂久久9| 国产精品香港三级国产av潘金莲| 三级毛片av免费| 变态另类成人亚洲欧美熟女 | 91在线观看av| 免费高清视频大片| 大香蕉久久成人网| 亚洲色图av天堂| 亚洲av熟女| 国产亚洲精品一区二区www| 久久久久久久午夜电影| 国产精品一区二区在线不卡| av在线播放免费不卡| 中文字幕色久视频| 黑人巨大精品欧美一区二区蜜桃| 大香蕉久久成人网| 欧美乱色亚洲激情| 久久香蕉国产精品| 国产精品一区二区在线不卡| 日本三级黄在线观看| 久久中文看片网| 美女扒开内裤让男人捅视频| 一区二区三区高清视频在线| 午夜精品在线福利| 中文亚洲av片在线观看爽| 日日干狠狠操夜夜爽| 国产伦一二天堂av在线观看| 美女高潮到喷水免费观看| 精品久久久久久久久久免费视频| 婷婷六月久久综合丁香| 欧美中文日本在线观看视频| 亚洲精品一区av在线观看| 好男人在线观看高清免费视频 | 午夜两性在线视频| 男女之事视频高清在线观看| 亚洲av片天天在线观看| 久久精品影院6| 国产精品一区二区在线不卡| 日本欧美视频一区| 在线观看日韩欧美| 久久香蕉精品热| 一级,二级,三级黄色视频| 亚洲va日本ⅴa欧美va伊人久久| 日韩欧美免费精品| 成人欧美大片| 国产不卡一卡二| e午夜精品久久久久久久| 久久欧美精品欧美久久欧美| 中文字幕人妻熟女乱码| 老司机在亚洲福利影院| 亚洲成人国产一区在线观看| 美女扒开内裤让男人捅视频| 国产xxxxx性猛交| 十分钟在线观看高清视频www| 亚洲精品中文字幕在线视频| 51午夜福利影视在线观看| 黑人欧美特级aaaaaa片| 成人国产一区最新在线观看| 午夜福利在线观看吧| 香蕉久久夜色| 亚洲久久久国产精品| 亚洲人成伊人成综合网2020| 一区二区三区激情视频| 狠狠狠狠99中文字幕| 国产精品野战在线观看| 男女下面进入的视频免费午夜 | 波多野结衣av一区二区av| 好男人在线观看高清免费视频 | 欧美 亚洲 国产 日韩一| 欧美黄色片欧美黄色片| 色综合欧美亚洲国产小说| 日韩欧美国产一区二区入口| 在线免费观看的www视频| 一卡2卡三卡四卡精品乱码亚洲| 村上凉子中文字幕在线| 亚洲一区高清亚洲精品| 亚洲国产精品sss在线观看| 精品欧美一区二区三区在线| 日韩 欧美 亚洲 中文字幕| 老司机午夜十八禁免费视频| 国产野战对白在线观看| 国产麻豆69| 亚洲美女黄片视频| 成人国语在线视频| 在线视频色国产色| 在线av久久热| 色在线成人网| 首页视频小说图片口味搜索| 亚洲欧美日韩无卡精品| 十分钟在线观看高清视频www| 国产野战对白在线观看| 亚洲一区高清亚洲精品| 美女大奶头视频| 久久久国产成人精品二区| 国产精品久久久av美女十八| 狂野欧美激情性xxxx| 中文字幕人成人乱码亚洲影| 校园春色视频在线观看| 一区福利在线观看| 午夜久久久在线观看| 两个人免费观看高清视频| 欧美成狂野欧美在线观看| 亚洲电影在线观看av| 免费女性裸体啪啪无遮挡网站| 99久久99久久久精品蜜桃| 99国产综合亚洲精品| 无限看片的www在线观看| 亚洲一码二码三码区别大吗| 日韩国内少妇激情av| 亚洲av熟女| 欧美日韩精品网址| 国产精品 国内视频| 精品福利观看| 这个男人来自地球电影免费观看| 一级黄色大片毛片| 夜夜爽天天搞| 国产乱人伦免费视频| 电影成人av| 国产精品亚洲一级av第二区| 一本久久中文字幕| 亚洲一区二区三区色噜噜| 别揉我奶头~嗯~啊~动态视频| 天天躁夜夜躁狠狠躁躁| 欧美色欧美亚洲另类二区 | 无人区码免费观看不卡| 亚洲av成人不卡在线观看播放网| aaaaa片日本免费| 欧美av亚洲av综合av国产av| 欧美中文日本在线观看视频| 日本欧美视频一区| 熟女少妇亚洲综合色aaa.| 两个人免费观看高清视频| 黑人欧美特级aaaaaa片| 在线天堂中文资源库| 久久精品国产综合久久久| 亚洲九九香蕉| 亚洲自偷自拍图片 自拍| 亚洲国产精品成人综合色| 黑人巨大精品欧美一区二区mp4| 桃红色精品国产亚洲av| 老汉色∧v一级毛片| 国产欧美日韩综合在线一区二区| 伊人久久大香线蕉亚洲五| x7x7x7水蜜桃| 日本欧美视频一区| 精品久久久久久成人av| 少妇被粗大的猛进出69影院| 9色porny在线观看| 成人免费观看视频高清| 国产成+人综合+亚洲专区| 国产精品 国内视频| 国产午夜福利久久久久久| 国产精品美女特级片免费视频播放器 | 在线十欧美十亚洲十日本专区| 亚洲国产欧美一区二区综合| 亚洲成人免费电影在线观看| 国产色视频综合| 亚洲专区中文字幕在线| 狂野欧美激情性xxxx| 国产免费av片在线观看野外av| 淫妇啪啪啪对白视频| 亚洲av日韩精品久久久久久密| 亚洲一区二区三区不卡视频| 欧美成人性av电影在线观看| 夜夜看夜夜爽夜夜摸| 国产成人av教育| 久久人人精品亚洲av| 91麻豆av在线| 久久中文看片网| 久久国产亚洲av麻豆专区| 国产亚洲精品久久久久5区| 欧美日韩精品网址| 亚洲成av人片免费观看| 波多野结衣高清无吗| 成人亚洲精品av一区二区| 欧美黄色片欧美黄色片| 热re99久久国产66热| 免费女性裸体啪啪无遮挡网站| 精品福利观看| 久久午夜亚洲精品久久| 婷婷精品国产亚洲av在线| 校园春色视频在线观看| 美女高潮喷水抽搐中文字幕| 一级毛片精品| 中文字幕另类日韩欧美亚洲嫩草| 9热在线视频观看99| 香蕉国产在线看| 亚洲欧美精品综合一区二区三区| 在线观看免费日韩欧美大片| www.自偷自拍.com| 久久久久国产一级毛片高清牌| 午夜免费观看网址| 国产激情欧美一区二区| av在线天堂中文字幕| av免费在线观看网站| 69精品国产乱码久久久| 在线视频色国产色| 亚洲欧洲精品一区二区精品久久久| 亚洲中文av在线|