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

    基于高分六號影像的杭州市黑臭水體遙感識別研究

    2022-10-09 07:12:52黃祺宇于之鋒彭曉雪袁小紅
    關(guān)鍵詞:黑臭反射率波段

    黃祺宇,肖 晗,于之鋒,3,周 斌,3,彭曉雪,袁小紅,3

    (1. 杭州師范大學(xué)遙感與地球科學(xué)研究院,浙江 杭州311121; 2. 武漢大學(xué)遙感信息工程學(xué)院,湖北 武漢 430079;3. 浙江省城市濕地與區(qū)域變化研究重點實驗室,浙江 杭州 311121)

    城市黑臭水體是指城市建成區(qū)內(nèi)呈現(xiàn)令人不悅的顏色和散發(fā)令人不適氣味的水體[1].城市黑臭水體不僅對城市水環(huán)境造成惡劣影響,而且損害居民身體健康,降低生活品質(zhì)[2].我國政府針對城市黑臭水體采取了一系列措施,其中快速準(zhǔn)確地獲取城市黑臭水體的時空分布是及時治理黑臭水體或阻止正常水體黑臭化的關(guān)鍵.黑臭水體時空分布信息來源于城市水體監(jiān)測,相比人工巡航采樣結(jié)合實驗室化學(xué)分析,遙感監(jiān)測城市黑臭水體速度更快、范圍更大且時效性更強,已被世界各地用作長期監(jiān)測城市黑臭水體的主要技術(shù)手段[3].隨著我國遙感與衛(wèi)星技術(shù)的不斷革新突破,高分(GF)系列、資源(ZY)系列等衛(wèi)星的陸續(xù)發(fā)射及其數(shù)據(jù)應(yīng)用的推廣,許多國產(chǎn)高分辨率衛(wèi)星影像已經(jīng)能夠支持并實現(xiàn)對城市細(xì)小、零碎河流的觀測,提高了城市黑臭水體的監(jiān)測精度[4].

    在對既定研究區(qū)黑臭水體識別應(yīng)用方面,國內(nèi)外已經(jīng)建立并驗證諸多適用性強且精度高的定量化監(jiān)測評價模型.Bai等[5]基于SeaWiFS衛(wèi)星的原位數(shù)據(jù)和海洋顏色遙感數(shù)據(jù),監(jiān)測到長江河口存在黑水區(qū)域并解釋了其異常光學(xué);Leif等[6]利用Landsat影像反演了美國明尼蘇達州各個區(qū)位湖泊的透明度,證實了從遙感影像中可以獲得時間和空間全覆蓋的水質(zhì)特征.Zou等[7]通過對典型黑臭水聚集區(qū)分析,開發(fā)了基于DN值和研究區(qū)域方差的黑臭水聚集區(qū)經(jīng)驗檢索算法并證明其精度.此外,有部分學(xué)者基于某一城市的黑臭水體構(gòu)建影像監(jiān)測算法,利用實測水體光譜結(jié)合遙感影像,分析黑臭水體與一般水體的光譜差異,再通過波段組合構(gòu)建黑臭水體判別模型.例如,溫爽等[8]分析了南京地區(qū)的黑臭水體光譜與一般水體光譜的特征區(qū)別,構(gòu)建了基于高分二號影像的單波段閾值、差值、比值、色度值的黑臭識別模型.胡國慶等[9]基于高分二號,利用單波段閾值法、波段差值法、歸一化指數(shù)法、斜率指數(shù)法對安徽省蕪湖市城市黑臭水現(xiàn)狀進行識別.紀(jì)剛[10]基于高分一號和高分二號衛(wèi)星影像,構(gòu)建了北京市黑臭水體識別系列模型(HeiChou Identification, HCI).曹紅業(yè)[11]綜合多個城市各類水體的固有光學(xué)量、表觀光學(xué)量和水質(zhì)參數(shù)特性,提出了基于實測遙感反射率的飽和度法和光譜指數(shù)法.張雪等[12]基于高分一號衛(wèi)星,利用單波段法、差值法、比值法、HCI1和HCI2算法對深圳市黑臭水體進行識別.韓文聰?shù)萚13]基于高分二號,使用營養(yǎng)狀態(tài)指數(shù)、歸一化黑臭水體指數(shù)(normalized difference black-odorousous water index, NDBWI)、反射率光譜指數(shù)(black-odorous water index, BOI)3種方法對浙江省寧??h進行疑似黑臭水體的判讀和識別,取得較好成果;姚月等[14]使用高分二號衛(wèi)星影像,基于遙感反射率、瑞利散射校正反射率遙感識別沈陽市黑臭水體,取得較好的識別精度.

    目前,遙感水色識別方法及城市黑臭水體識別模型雖在實際應(yīng)用中具有較高的精度,但仍存在不足,如何選定最優(yōu)閾值或自動化提取水體是各算法模型需要克服的重大難題之一[15].不僅如此,算法應(yīng)用于不同衛(wèi)星影像的反演效果不盡相同,將算法移植至其他研究區(qū)或應(yīng)用于其他衛(wèi)星影像的反演效果仍有待檢驗.本研究基于高分六號衛(wèi)星影像,利用杭州市城區(qū)實測黑臭水體光譜,應(yīng)用5種黑臭水體識別模型并進行精度檢驗,探究每種黑臭水體識別算法在杭州地區(qū)的適用性.

    1 研究區(qū)與數(shù)據(jù)

    1.1 研究區(qū)概況

    杭州市區(qū)河道縱橫交錯、河網(wǎng)密布,共有河道1 400余條,主城區(qū)內(nèi)河道461條[16].依照浙江省制定的黑臭河標(biāo)準(zhǔn)和地表水環(huán)境質(zhì)量標(biāo)準(zhǔn),杭州市在開展“五水共治”之前有348條垃圾河和黑臭河,在納入監(jiān)測范圍的1 334條河道中,劣Ⅴ類水河道約占40%,主城區(qū)劣Ⅴ類河道占河道總數(shù)的2/3以上[17].目前,主要河道已基本消除了黑臭,但由于歷史因素,加上杭州河道數(shù)量眾多,市區(qū)范圍內(nèi)河道仍存在漏測、“復(fù)黑”、富營養(yǎng)化等潛在風(fēng)險,水質(zhì)亟待改善[18].

    1.2 數(shù)據(jù)采集

    2019年3月23日—24日,對杭州市拱墅區(qū)后橫港河、陶家圩,西湖區(qū)族濱漾、西行河、余杭塘河、十字港橋,余杭區(qū)蓬橋河港、紅衛(wèi)港,濱江區(qū)陸家水潭、永久河、傅家峙河、善慶莊橫河,共12條河道展開黑臭水體的野外實地測量.測量時,每條河道選擇上、中、下游3部分至少一個樣本點,并且要求該樣本點能代表較大面積河段的水質(zhì),在遇到有明顯水質(zhì)變化的河段時增加采樣,最終共選擇了62個樣本點:余杭區(qū)點簇7個樣點,拱墅區(qū)點簇7個樣點,西湖區(qū)從左至右,從上至下點簇分別為4個、7個、5個、7個、5個樣點,濱江區(qū)從左至右,從上至下點簇分別為4個、7個、5個、4個樣點(圖1).

    圖1 實測水體點位分布Fig.1 Point distribution of measured water body

    水體光譜數(shù)據(jù)采集時,使用手持式光譜儀(analytical spectral devices,ASD)和30%反射率的標(biāo)準(zhǔn)板,依據(jù)水面以上測量法[19]獲取各站點的水體光譜數(shù)據(jù),剔除各觀測點中偏差較大的異常光譜,計算剩余光譜數(shù)據(jù)的水體遙感反射率.溶解氧與氧化電位采用多參數(shù)水質(zhì)儀獲取,依照水質(zhì)儀操作步驟[20],設(shè)定待測水質(zhì)參數(shù)后,開啟自動測量模式獲取相應(yīng)水質(zhì)參數(shù)信息.使用透明度儀測量水體的透明度信息,取2次測定平均值作為水體透明度值.

    (1)

    式中,Rrs為遙感反射率,Lu、Lsky、Lp分別為光譜儀面向水體、天空和灰板時的輻亮度測量信號,ρp為灰板反射率,ρf為菲涅爾反射率.依據(jù)測站點的風(fēng)速設(shè)置ρf的數(shù)值,平靜水面可取ρf=0.022;在5 m/s左右風(fēng)速下,ρf可取0.025.

    1.3 高分六號影像處理

    本文使用最貼近實測日期的可用影像數(shù)據(jù)(2019年3月12日的兩幅高分六號影像數(shù)據(jù)),數(shù)據(jù)來源于中國資源衛(wèi)星應(yīng)用中心,影像無云覆蓋,可用于杭州市黑臭水體的識別.利用ENVI軟件對高分六號影像進行數(shù)據(jù)預(yù)處理,包括正射校正、輻射定標(biāo)、大氣校正等.大氣校正采用FLAASH大氣校正模塊,由于高分六號缺少大氣校正氣溶膠反演相應(yīng)的波段[21-22],因此不采用氣溶膠反演.為減少建筑物高亮部分對水體識別的影響,對影像進行水體掩膜,使用NCWI[23]提取影像水體部分,并對感興趣區(qū)(region of interest,ROI)進行細(xì)節(jié)修正排除建筑物高亮部分,保留細(xì)小水體,剔除混合像元.

    2 城市黑臭水體遙感識別算法

    2.1 黑臭水體光譜特征分析

    依據(jù)曹紅業(yè)[11]提出的新黑臭水體判別指標(biāo)(以下簡稱新指標(biāo)),對一般水體與黑臭水體進行判別.新指標(biāo)基于《城市黑臭水體整治工作指南》[24](以下簡稱舊指標(biāo))且更加成熟:溶解氧(dissolved oxygen, DO)水平由2 mg/L提升至3.5 mg/L,氧化還原電位(oxidation reduction potential, ORP)由50 mV提升至60 mV,氨氮水平由8 mg/L降至6 mg/L.

    依照城市黑臭水體判別標(biāo)準(zhǔn)中新指標(biāo)(表1),根據(jù)各采樣點的氨氮水平、透明度、溶解氧水平、氧化還原電位,判別各采樣點處水體是否為黑臭水體.最終62個樣本點的數(shù)據(jù)中,黑臭水體樣36個、正常水體樣26個.依據(jù)分類結(jié)果對正常河段和黑臭河段的樣本點光譜數(shù)據(jù)進行處理,計算得到遙感反射率,繪制光譜曲線(圖2),對比分析黑臭水體與正常水體的光譜特征.

    表1 城市黑臭水體判別標(biāo)準(zhǔn)[12]Tab.1 Discrimination standard of urban black-odorous water body

    a:黑臭水體; b:正常水體.

    在波長400~550 nm范圍內(nèi),黑臭水體與正常水體的遙感反射率數(shù)值均上升,但黑臭水體的整體上升趨勢更加平緩.大部分黑臭水體含氧量低,藻類植物數(shù)量少,因此相較于正常水體,其670 nm附近吸收谷不明顯,700 nm附近反射峰不明顯[25].此外,760 nm處正常水體具有明顯的液態(tài)水特征(吸收形成的吸收谷),而黑臭水體無該特征.總體而言,杭州黑臭水體的遙感反射率(<0.038 sr-1)低于正常水體的遙感反射率.永久河與十字港橋水體清澈,受水底河床干擾,其遙感反射率較低.杭州黑臭水體的遙感反射率走勢較正常水體平緩,波峰波谷不明顯,數(shù)值與正常水體有明顯的分界閾值.利用高分六號光譜響應(yīng)函數(shù)和野外實際測量計算得到的水體遙感反射率均值可生成高分六號模擬結(jié)果,對比圖3顯示,高分六號影像大大縮減了水體的光譜特征細(xì)節(jié)差異[26].在斜率與數(shù)值上,正常水體與黑臭水體存在較大差異,這與余佳龍等[17]研究結(jié)果一致.因此,利用高分六號監(jiān)測黑臭水體時可通過構(gòu)建藍光、綠光、紅光和近紅外波段的相關(guān)算法識別黑臭水體.

    圖3 基于兩類水體反射率均值的高分六號模擬結(jié)果Fig.3 Simulated average reflectivity of two kinds of water body multi-spectrum of GF-6

    利用采樣點的光譜數(shù)據(jù)建立單波段法、差值法、比值法、HCI1、HCI2[10,12,27]5種黑臭水體識別算法,對基于高分六號影像的黑臭水體識別算法的精度進行驗證和對比分析.由于不同地區(qū)的地理位置和水質(zhì)環(huán)境等具有差異,不同城市的正常水體和黑臭水體的光譜有差異[12].因此,上述黑臭水體識別算法雖在某一地區(qū)能有效區(qū)分黑臭水體和正常水體,但在不同地區(qū)應(yīng)用時識別精度易受影響,故需修正提取算法的閾值,并進行精度評價[28-29].從高分六號模擬結(jié)果中,隨機篩選出15個黑臭水體樣本和10個正常水體樣本,對算法的閾值進行修正.如圖4所示,紅線處為表2所列的水體區(qū)分閾值,RR為紅光波段、RG為綠光波段、RB為藍光波段、RN為近紅外波段,N1、N2、N3、N4、N5、N6分別為各方法中使用到的閾值.

    圖4 各算法建模樣點Fig.4 Modeling sample points of each algorithm

    表2 各算法閾值修正結(jié)果Tab.2 Threshold correction results of each algorithm

    2.2 精度分析

    將表2中的5種黑臭水體識別算法應(yīng)用于經(jīng)過處理的高分六號影像,提取影像中的黑臭水體.根據(jù)與衛(wèi)星同步的實測水體進行精度評價,評價指標(biāo)采用識別正確率,結(jié)果如表3所示.

    表3 黑臭水體的影像提取算法識別正確率Tab.3 Correct recognition rate of image extraction algorithm for black-odorous water body

    (2)

    式中,N為該水體類別中總樣本個數(shù),M為應(yīng)用判別模型后識別正確的樣本個數(shù).

    結(jié)果表明,5種識別算法中,比值法和HCI2方法的精度較低,且各類算法將正常水體錯分為黑臭水體的現(xiàn)象較為普遍,可能原因如下:

    1)杭州地區(qū)部分水體水質(zhì)較好、水體透明度高、水體光譜吸收強、水體反射率值較低,在水體識別時易被錯分為黑臭水體.如永久河水質(zhì)較好,水體反射率反而較低,且在綠光和紅光波段上反射率值較為接近,因而在比值法和HCI2中的區(qū)分度較差.水體很淺而導(dǎo)致底泥影響水體光譜特性較大的情況,在水體提取階段已被剔除,因此不需考慮底泥光譜信息導(dǎo)致的誤判.

    2)在測量的黑臭水體中,善慶莊橫河的水體表面存在嚴(yán)重油污,油污遮蓋水體表面導(dǎo)致水體遙感反射率偏高,因此該水體易被錯分為正常水體.

    3)進行大氣校正時,由于高分六號影像不存在大氣校正氣溶膠反演相應(yīng)的波段,大氣校正結(jié)果可能與實際情況存在一定偏差.

    4)杭州市長期以來積極開展水體的清潔工作,在實際測量時也發(fā)現(xiàn)有清潔人員在對水體中的異物進行清除,影像成像時間與實際測量時間具有一定時間間隔,可能存在短時間內(nèi)水體遙感反射率發(fā)生變化的情況,導(dǎo)致數(shù)據(jù)誤差.

    在5種識別算法中,波段差值方法識別正常水體的正確率相對較高,正常水體識別正確點位為23個,識別正確率達88.5%.HCI1方法識別黑臭水體的正確率較高,黑臭水體識別正確點位為34個,識別正確率達94.4%.

    采用不同算法提取的黑臭水體如圖5所示.結(jié)果顯示,單波段法提取出黑臭河段13條,總長度97.30 km;波段差值法提取出黑臭河段17條,總長度124.433 km;波段比值法提取出黑臭河段13條,總長度89.194 km;3種算法在提取拱墅區(qū)黑臭水體時基本一致,提取上城區(qū)、濱江區(qū)和蕭山區(qū)黑臭水體時有較小差異,提取余杭區(qū)黑臭水體差異較大.

    (a)單波段法

    2.3 5種算法適用性分析與評價

    1)單波段法:單波段算法最大的優(yōu)勢在于運算快、所需數(shù)據(jù)量少,并且具有較高的正常水體識別精度,但對黑臭水體識別精度明顯不如另外4種算法,在以精度為目的的反演情況下不推薦使用.在隨機抽樣調(diào)查城市河道水質(zhì)或判別中大型城市湖泊等水體時,單波段法可以對黑臭水體進行快速定性.

    2)差值法:差值法對于黑臭水體識別精度處于中等水平,差值法可將正常水體與黑臭水體光譜特征差異放大.另外,差值法在識別正常水體時正確率最高,可以通過同一水體在時間序列上的差異來監(jiān)測整治進度.

    3)比值法:比值法對黑臭水體的正確識別率較高,達88.9%,僅次于HCI1法.但其錯分現(xiàn)象最為嚴(yán)重,推測原因可能與清澈水體在紅光、綠光波段反射率差距不大有較大關(guān)聯(lián),因此不建議在識別少數(shù)幾處黑臭水體時使用.

    4)HCI1法:HCI1法對黑臭水體識別正確率最高,達94.4%,基本可以滿足對黑臭水體的精確監(jiān)測要求.因此在需要準(zhǔn)確識別城市中所有黑臭水體時,HCI1法為最佳算法.

    5)HCI2法:HCI2法對于黑臭水體和正常水體識別準(zhǔn)確度均較差,且錯分率較高,對于實際監(jiān)測城市黑臭水體的價值不如其他4種.總體而言不推薦使用該方法進行黑臭水體遙感識別.

    3 結(jié)論與展望

    基于高分六號影像以及現(xiàn)有的黑臭水體識別算法,對杭州市黑臭水體識別算法的適用性進行探究.依據(jù)實測的杭州市黑臭水體與正常水體的光譜數(shù)據(jù),對各算法進行閾值修正,并進行精度驗證,得到以下結(jié)論:

    1)與正常水體相比,杭州市黑臭水體反射率相對平緩且數(shù)值較低,但西湖、西溪濕地等較為清潔的湖泊水體因光譜吸收強、反射率低,易被識別為黑臭水體.

    2)現(xiàn)有黑臭水體算法對基于高分六號的黑臭水體識別均有一定的應(yīng)用價值,大部分算法模型對黑臭水體與正常水體均有較高的識別準(zhǔn)確率.

    3)對差值法或HCI1法進一步優(yōu)化,有可能按黑臭程度對黑臭水體進行分級識別,為城市黑臭水整治提供更多科學(xué)依據(jù).

    4)基于高分六號影像反演識別城市黑臭水體的5種方法中,單波段法效率高、速度快;HCI1法對黑臭水體識別精度最高;差值法對正常水體識別精度最高.

    另外,后續(xù)還需改進高分六號影像的大氣校正,提高對黑臭水體的識別精度,對黑臭水體識別精度較高的HCI1方法和差值法進行進一步研究.

    猜你喜歡
    黑臭反射率波段
    春日暖陽
    影響Mini LED板油墨層反射率的因素
    近岸水體異源遙感反射率產(chǎn)品的融合方法研究
    海洋通報(2022年4期)2022-10-10 07:40:26
    具有顏色恒常性的光譜反射率重建
    基于“高分二號”衛(wèi)星影像的黑臭水體識別
    應(yīng)用生態(tài)修復(fù)技術(shù)治理黑臭河道研究
    綠色科技(2019年16期)2019-11-22 14:39:17
    利用微生物技術(shù)修復(fù)黑臭和富營養(yǎng)水體的研究
    科技資訊(2018年6期)2018-07-26 10:40:56
    溫州市甌海區(qū)排污口整治與河道生態(tài)復(fù)合處理方法及技術(shù)
    綠色科技(2017年10期)2017-07-05 16:17:22
    化學(xué)腐蝕硅表面結(jié)構(gòu)反射率影響因素的研究*
    電子器件(2017年2期)2017-04-25 08:58:37
    M87的多波段輻射過程及其能譜擬合
    两个人的视频大全免费| 免费人成视频x8x8入口观看| 免费高清视频大片| 欧美黑人欧美精品刺激| 国产精品98久久久久久宅男小说| 男女那种视频在线观看| 亚洲国产中文字幕在线视频| 久久午夜亚洲精品久久| 亚洲黑人精品在线| 90打野战视频偷拍视频| 88av欧美| 一二三四社区在线视频社区8| 亚洲美女视频黄频| 黄色女人牲交| 精华霜和精华液先用哪个| 色在线成人网| 男女午夜视频在线观看| 人人妻,人人澡人人爽秒播| 999精品在线视频| 国产精品久久久久久亚洲av鲁大| 1000部很黄的大片| 999精品在线视频| 91字幕亚洲| 好男人电影高清在线观看| 国产高清视频在线播放一区| 伊人久久大香线蕉亚洲五| 午夜福利高清视频| 一进一出抽搐gif免费好疼| 亚洲中文字幕日韩| 国产精品亚洲一级av第二区| 国产视频一区二区在线看| 国产黄片美女视频| 女同久久另类99精品国产91| 欧美三级亚洲精品| 黄色 视频免费看| 老司机在亚洲福利影院| 白带黄色成豆腐渣| 亚洲人成网站在线播放欧美日韩| 一级黄色大片毛片| 变态另类成人亚洲欧美熟女| 国产亚洲av高清不卡| 亚洲电影在线观看av| 一区二区三区高清视频在线| 男女午夜视频在线观看| 欧美+亚洲+日韩+国产| 一级黄色大片毛片| 小说图片视频综合网站| 午夜两性在线视频| 日本熟妇午夜| 日本熟妇午夜| 国产精品99久久99久久久不卡| 手机成人av网站| 亚洲va日本ⅴa欧美va伊人久久| 日本熟妇午夜| 欧美日韩综合久久久久久 | 床上黄色一级片| 亚洲熟妇中文字幕五十中出| 国产91精品成人一区二区三区| 国产精品一区二区三区四区免费观看 | 在线免费观看不下载黄p国产 | 久久这里只有精品19| 亚洲精品久久国产高清桃花| 国产亚洲av嫩草精品影院| 国产淫片久久久久久久久 | 亚洲av电影不卡..在线观看| 国产97色在线日韩免费| 国产真实乱freesex| 中国美女看黄片| 嫁个100分男人电影在线观看| 欧美激情在线99| 亚洲精品中文字幕一二三四区| 国产精品精品国产色婷婷| 国内毛片毛片毛片毛片毛片| 免费大片18禁| 少妇裸体淫交视频免费看高清| 免费看美女性在线毛片视频| 免费看美女性在线毛片视频| 99国产精品一区二区蜜桃av| 又大又爽又粗| 欧美黑人巨大hd| 欧美激情在线99| 桃红色精品国产亚洲av| 亚洲成人久久爱视频| 亚洲欧美激情综合另类| 成年女人毛片免费观看观看9| 97超级碰碰碰精品色视频在线观看| 可以在线观看的亚洲视频| 久久久久性生活片| 性欧美人与动物交配| 欧美日韩亚洲国产一区二区在线观看| 国产视频内射| 少妇的丰满在线观看| 国产精品亚洲av一区麻豆| 国产97色在线日韩免费| 久久久国产成人精品二区| 国产精品永久免费网站| 脱女人内裤的视频| 亚洲av片天天在线观看| 亚洲一区二区三区不卡视频| 欧美激情久久久久久爽电影| 日本成人三级电影网站| 亚洲国产高清在线一区二区三| 精品一区二区三区视频在线观看免费| 我的老师免费观看完整版| 999久久久国产精品视频| 国产在线精品亚洲第一网站| 久久伊人香网站| 午夜精品一区二区三区免费看| 中文字幕最新亚洲高清| 免费在线观看日本一区| 丁香六月欧美| 国产成人精品久久二区二区免费| 国产伦一二天堂av在线观看| 亚洲最大成人中文| 亚洲国产高清在线一区二区三| 日本黄大片高清| 亚洲 国产 在线| 97碰自拍视频| 亚洲成人精品中文字幕电影| 在线永久观看黄色视频| 亚洲中文av在线| 偷拍熟女少妇极品色| 啪啪无遮挡十八禁网站| 国产亚洲精品久久久com| 国内精品久久久久久久电影| 亚洲av电影不卡..在线观看| 中文字幕人妻丝袜一区二区| 婷婷六月久久综合丁香| 国产伦精品一区二区三区视频9 | 日韩精品中文字幕看吧| 两个人的视频大全免费| 亚洲国产精品sss在线观看| 男人舔女人下体高潮全视频| 哪里可以看免费的av片| 亚洲七黄色美女视频| 99riav亚洲国产免费| aaaaa片日本免费| 美女扒开内裤让男人捅视频| 日韩欧美在线二视频| 中文字幕高清在线视频| 免费在线观看视频国产中文字幕亚洲| 国产毛片a区久久久久| 我的老师免费观看完整版| 在线视频色国产色| 长腿黑丝高跟| 69av精品久久久久久| 亚洲色图 男人天堂 中文字幕| 国产综合懂色| 9191精品国产免费久久| www日本在线高清视频| 老熟妇乱子伦视频在线观看| 成年女人看的毛片在线观看| 午夜福利在线在线| 草草在线视频免费看| 国产一区二区激情短视频| 成人无遮挡网站| 国产成人aa在线观看| svipshipincom国产片| 啦啦啦观看免费观看视频高清| 亚洲av成人精品一区久久| 国产精品亚洲美女久久久| 午夜福利在线观看免费完整高清在 | 国产av麻豆久久久久久久| 白带黄色成豆腐渣| 亚洲专区中文字幕在线| 亚洲成人精品中文字幕电影| 久久精品国产清高在天天线| or卡值多少钱| 丝袜人妻中文字幕| 国产成人av教育| 校园春色视频在线观看| 国产探花在线观看一区二区| 欧美不卡视频在线免费观看| 97碰自拍视频| 欧美午夜高清在线| 偷拍熟女少妇极品色| 久久久久九九精品影院| e午夜精品久久久久久久| 精品久久久久久成人av| 国产亚洲精品综合一区在线观看| 黄色成人免费大全| 2021天堂中文幕一二区在线观| 久久这里只有精品中国| 国内毛片毛片毛片毛片毛片| www.www免费av| 亚洲一区高清亚洲精品| 亚洲中文日韩欧美视频| av福利片在线观看| 日本五十路高清| 久久中文看片网| 嫩草影院入口| 狂野欧美激情性xxxx| 免费观看精品视频网站| 夜夜看夜夜爽夜夜摸| 少妇丰满av| 一个人免费在线观看电影 | 国产精品乱码一区二三区的特点| 国产亚洲av高清不卡| 亚洲精品在线观看二区| 国产精品免费一区二区三区在线| av在线天堂中文字幕| 国产成人一区二区三区免费视频网站| 亚洲avbb在线观看| 国产亚洲精品一区二区www| 两个人的视频大全免费| 亚洲精品一区av在线观看| 99热只有精品国产| 51午夜福利影视在线观看| 欧美色视频一区免费| 欧美国产日韩亚洲一区| 黑人操中国人逼视频| 老熟妇仑乱视频hdxx| 亚洲av美国av| 亚洲av第一区精品v没综合| 两人在一起打扑克的视频| 午夜久久久久精精品| 欧美精品啪啪一区二区三区| 亚洲精品乱码久久久v下载方式 | 美女高潮的动态| 看免费av毛片| 国产精品99久久久久久久久| av在线蜜桃| 宅男免费午夜| 亚洲精品中文字幕一二三四区| 给我免费播放毛片高清在线观看| av片东京热男人的天堂| 欧美国产日韩亚洲一区| 欧美乱码精品一区二区三区| 精品不卡国产一区二区三区| 波多野结衣高清无吗| 一区二区三区高清视频在线| 亚洲国产看品久久| 久久欧美精品欧美久久欧美| 亚洲18禁久久av| 一边摸一边抽搐一进一小说| x7x7x7水蜜桃| 欧美最黄视频在线播放免费| 在线观看66精品国产| 国产精品久久久av美女十八| 美女黄网站色视频| 国产精品久久久久久人妻精品电影| 成人性生交大片免费视频hd| 亚洲人成网站高清观看| 757午夜福利合集在线观看| 精品久久久久久,| 国产一区二区三区在线臀色熟女| 国产精品 国内视频| 亚洲人成伊人成综合网2020| 啦啦啦观看免费观看视频高清| 韩国av一区二区三区四区| 日韩欧美国产一区二区入口| 一二三四社区在线视频社区8| 欧美激情在线99| 免费看日本二区| 一区二区三区国产精品乱码| 午夜福利在线观看吧| 亚洲欧美激情综合另类| 亚洲午夜理论影院| 最好的美女福利视频网| 嫩草影院精品99| 麻豆国产av国片精品| 熟女少妇亚洲综合色aaa.| 久久久久久久久中文| 亚洲专区中文字幕在线| 中文字幕人妻丝袜一区二区| 禁无遮挡网站| 国产精品永久免费网站| 欧美日韩精品网址| 久久久久亚洲av毛片大全| 久久久久久久午夜电影| 亚洲色图 男人天堂 中文字幕| 国产日本99.免费观看| 91在线观看av| 成人国产一区最新在线观看| 久久久水蜜桃国产精品网| tocl精华| 国产成人啪精品午夜网站| 国产又黄又爽又无遮挡在线| 亚洲精品一卡2卡三卡4卡5卡| 嫩草影视91久久| 美女黄网站色视频| 18美女黄网站色大片免费观看| 亚洲精品456在线播放app | 99精品在免费线老司机午夜| 12—13女人毛片做爰片一| 亚洲美女视频黄频| 99riav亚洲国产免费| 精品国产乱子伦一区二区三区| 亚洲专区国产一区二区| 成人三级黄色视频| 嫩草影院入口| 国产在线精品亚洲第一网站| 精品一区二区三区视频在线观看免费| 午夜免费观看网址| 两性夫妻黄色片| 欧美国产日韩亚洲一区| av天堂中文字幕网| 夜夜看夜夜爽夜夜摸| 天堂动漫精品| 别揉我奶头~嗯~啊~动态视频| 精品久久久久久久久久久久久| 亚洲国产精品合色在线| 久久精品夜夜夜夜夜久久蜜豆| а√天堂www在线а√下载| 亚洲国产中文字幕在线视频| 好看av亚洲va欧美ⅴa在| 中文字幕人成人乱码亚洲影| 欧美黑人巨大hd| 国产 一区 欧美 日韩| 亚洲成人久久性| 久久久精品欧美日韩精品| 91麻豆av在线| 欧美高清成人免费视频www| 国产男靠女视频免费网站| 叶爱在线成人免费视频播放| 后天国语完整版免费观看| 黑人欧美特级aaaaaa片| 色哟哟哟哟哟哟| 亚洲精品粉嫩美女一区| 国产成人一区二区三区免费视频网站| 欧美一级a爱片免费观看看| 麻豆久久精品国产亚洲av| 欧美一区二区国产精品久久精品| 18禁裸乳无遮挡免费网站照片| 国产激情久久老熟女| 国产熟女xx| 欧美在线黄色| av视频在线观看入口| 99久久成人亚洲精品观看| 三级国产精品欧美在线观看 | 国产精品久久电影中文字幕| 久久香蕉精品热| 国产精品九九99| 欧美日韩福利视频一区二区| 一级毛片精品| 免费看a级黄色片| 成人午夜高清在线视频| 欧美国产日韩亚洲一区| 国产伦在线观看视频一区| 中文字幕人成人乱码亚洲影| 禁无遮挡网站| 国产精品野战在线观看| 午夜福利18| 国产高清有码在线观看视频| 观看免费一级毛片| 男人舔女人的私密视频| 久久久水蜜桃国产精品网| 三级毛片av免费| 熟妇人妻久久中文字幕3abv| 草草在线视频免费看| 亚洲成av人片免费观看| 丁香六月欧美| 久久久久久人人人人人| 19禁男女啪啪无遮挡网站| 国产淫片久久久久久久久 | 成年版毛片免费区| 久久精品国产亚洲av香蕉五月| 免费搜索国产男女视频| 日韩欧美国产在线观看| 欧美三级亚洲精品| 国产午夜精品久久久久久| 国产日本99.免费观看| 手机成人av网站| 国产主播在线观看一区二区| 后天国语完整版免费观看| 啦啦啦观看免费观看视频高清| 此物有八面人人有两片| 琪琪午夜伦伦电影理论片6080| bbb黄色大片| 亚洲午夜理论影院| av国产免费在线观看| 中亚洲国语对白在线视频| 制服人妻中文乱码| 老熟妇仑乱视频hdxx| av视频在线观看入口| 国产 一区 欧美 日韩| 午夜视频精品福利| 色综合站精品国产| 国产精品亚洲av一区麻豆| 一个人看视频在线观看www免费 | 在线国产一区二区在线| 免费看a级黄色片| 又大又爽又粗| 男女那种视频在线观看| 国产av在哪里看| 国产精品 国内视频| 国产1区2区3区精品| 成人av一区二区三区在线看| 久久中文字幕人妻熟女| 国产精品 欧美亚洲| h日本视频在线播放| 中文字幕人妻丝袜一区二区| 午夜精品久久久久久毛片777| 波多野结衣巨乳人妻| 日韩 欧美 亚洲 中文字幕| 欧美成人一区二区免费高清观看 | 色噜噜av男人的天堂激情| 美女cb高潮喷水在线观看 | 亚洲国产看品久久| 国产成年人精品一区二区| www.www免费av| 一区二区三区国产精品乱码| 一a级毛片在线观看| 国产淫片久久久久久久久 | 日韩中文字幕欧美一区二区| 国产av不卡久久| 狂野欧美激情性xxxx| 白带黄色成豆腐渣| АⅤ资源中文在线天堂| 欧美色欧美亚洲另类二区| 在线观看免费视频日本深夜| a级毛片a级免费在线| 伦理电影免费视频| 神马国产精品三级电影在线观看| 中文字幕熟女人妻在线| 最近最新中文字幕大全电影3| 悠悠久久av| 国产精品女同一区二区软件 | 亚洲avbb在线观看| 国产高清视频在线播放一区| 麻豆国产97在线/欧美| 麻豆成人午夜福利视频| 国产69精品久久久久777片 | 亚洲国产日韩欧美精品在线观看 | 首页视频小说图片口味搜索| 神马国产精品三级电影在线观看| 成人三级黄色视频| 色老头精品视频在线观看| 曰老女人黄片| 国产精品av视频在线免费观看| 国产成人精品无人区| www.精华液| av在线蜜桃| 日韩高清综合在线| 真人做人爱边吃奶动态| 99re在线观看精品视频| e午夜精品久久久久久久| 亚洲av五月六月丁香网| 欧美黄色片欧美黄色片| 欧美日韩亚洲国产一区二区在线观看| 九九在线视频观看精品| 亚洲av日韩精品久久久久久密| 美女高潮的动态| 999久久久国产精品视频| 欧美日韩乱码在线| 欧美大码av| 国产高清三级在线| 久久99热这里只有精品18| av女优亚洲男人天堂 | 一级a爱片免费观看的视频| 给我免费播放毛片高清在线观看| 在线观看66精品国产| 极品教师在线免费播放| 亚洲自偷自拍图片 自拍| 宅男免费午夜| 午夜成年电影在线免费观看| 女人高潮潮喷娇喘18禁视频| 美女 人体艺术 gogo| 日韩中文字幕欧美一区二区| 在线观看免费午夜福利视频| 国产精品98久久久久久宅男小说| 岛国在线免费视频观看| 国产熟女xx| 日本a在线网址| svipshipincom国产片| 久久伊人香网站| 欧美激情久久久久久爽电影| 搡老岳熟女国产| 免费人成视频x8x8入口观看| 成熟少妇高潮喷水视频| 一二三四在线观看免费中文在| 国产精品永久免费网站| 三级国产精品欧美在线观看 | 性色av乱码一区二区三区2| 久久九九热精品免费| 亚洲精品乱码久久久v下载方式 | 国内久久婷婷六月综合欲色啪| 亚洲专区中文字幕在线| 免费观看人在逋| 国产精品 欧美亚洲| 亚洲中文字幕日韩| 亚洲av第一区精品v没综合| 成人国产综合亚洲| 黄频高清免费视频| 岛国在线观看网站| 极品教师在线免费播放| 三级男女做爰猛烈吃奶摸视频| 色综合站精品国产| 国产伦精品一区二区三区四那| 久久国产乱子伦精品免费另类| 美女cb高潮喷水在线观看 | 欧美性猛交黑人性爽| 人妻久久中文字幕网| 在线观看一区二区三区| 国产精品久久视频播放| 久久草成人影院| 久久久久国内视频| 久久久国产成人精品二区| 午夜激情欧美在线| 成人欧美大片| 嫩草影视91久久| 久久精品人妻少妇| a级毛片a级免费在线| 亚洲色图 男人天堂 中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 亚洲色图av天堂| 他把我摸到了高潮在线观看| 十八禁人妻一区二区| 狂野欧美白嫩少妇大欣赏| 国内精品一区二区在线观看| 午夜福利18| 十八禁网站免费在线| 久久久国产成人免费| 亚洲国产欧美网| 999久久久国产精品视频| 久久精品综合一区二区三区| 两人在一起打扑克的视频| 午夜福利成人在线免费观看| 国产男靠女视频免费网站| 午夜福利成人在线免费观看| 男女下面进入的视频免费午夜| 国产久久久一区二区三区| 欧美在线黄色| 国产aⅴ精品一区二区三区波| 亚洲国产看品久久| 国产日本99.免费观看| 欧美在线黄色| 日本撒尿小便嘘嘘汇集6| 午夜福利在线观看免费完整高清在 | 制服人妻中文乱码| 国产91精品成人一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 在线免费观看的www视频| 法律面前人人平等表现在哪些方面| 男人舔奶头视频| 色尼玛亚洲综合影院| 久久亚洲精品不卡| 九色国产91popny在线| 97超视频在线观看视频| 丰满人妻熟妇乱又伦精品不卡| 1024香蕉在线观看| 亚洲国产精品久久男人天堂| 在线视频色国产色| 欧美成人免费av一区二区三区| 日韩欧美 国产精品| 热99re8久久精品国产| www.自偷自拍.com| 亚洲一区高清亚洲精品| 美女免费视频网站| 操出白浆在线播放| 村上凉子中文字幕在线| xxxwww97欧美| 日韩欧美一区二区三区在线观看| 欧美一级毛片孕妇| 亚洲欧美日韩无卡精品| 国产高清三级在线| 免费电影在线观看免费观看| 中文资源天堂在线| 叶爱在线成人免费视频播放| 久久久久国产精品人妻aⅴ院| 欧美性猛交黑人性爽| 精品福利观看| 黄频高清免费视频| 老司机深夜福利视频在线观看| 最好的美女福利视频网| 中文字幕熟女人妻在线| 欧美丝袜亚洲另类 | 日本成人三级电影网站| 久久中文字幕一级| 在线免费观看不下载黄p国产 | 成年人黄色毛片网站| 天天一区二区日本电影三级| 毛片女人毛片| 人妻夜夜爽99麻豆av| 久久久久国产一级毛片高清牌| 999久久久国产精品视频| 色av中文字幕| 日韩欧美精品v在线| 久久午夜综合久久蜜桃| 久久中文字幕一级| 国产亚洲精品一区二区www| 欧美成人一区二区免费高清观看 | 久久久久久久精品吃奶| 国产精品一区二区三区四区免费观看 | 久久欧美精品欧美久久欧美| 狂野欧美白嫩少妇大欣赏| 亚洲精品在线观看二区| 少妇的丰满在线观看| 校园春色视频在线观看| 久久久久久久久中文| 成年免费大片在线观看| 中文资源天堂在线| av视频在线观看入口| 亚洲精品一卡2卡三卡4卡5卡| 亚洲七黄色美女视频| 成年版毛片免费区| 性色avwww在线观看| 亚洲美女黄片视频| 久久久久国产一级毛片高清牌| 嫩草影视91久久| 亚洲一区高清亚洲精品| 在线永久观看黄色视频| 特级一级黄色大片| 亚洲国产欧美一区二区综合| 午夜福利在线观看免费完整高清在 | 欧美午夜高清在线| 国产成人av激情在线播放| 欧美日韩乱码在线| 久久香蕉精品热| 欧美一级毛片孕妇| АⅤ资源中文在线天堂| 日韩欧美在线乱码| 真人一进一出gif抽搐免费| 日本与韩国留学比较| 国产成人av激情在线播放|