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

    基于單體擴張的端元提取算法

    2016-03-27 08:04:32董安國龔文娟
    紅外技術(shù) 2016年11期
    關(guān)鍵詞:端元情形單體

    董安國,韓 雪,龔文娟

    ?

    基于單體擴張的端元提取算法

    董安國,韓 雪,龔文娟

    (長安大學(xué) 理學(xué)院,陜西 西安 710064)

    端元提取是高光譜混合像元分解的重要環(huán)節(jié)。為了提取高光譜圖像的端元,本文基于線性表示理論與凸錐模型理論,論證了:與單體共面的單體外向量被單體的頂點向量線性表示時,表示系數(shù)必有負值,從而給出了理想情形下判別端元的充要條件,并在此基礎(chǔ)上,針對非理想情形提出了一種提取端元的迭代算法。實驗結(jié)果表明,算法提取端元的精度優(yōu)于VCA算法、效率高于搜索算法,算法穩(wěn)定性好,對噪聲的敏感性低。

    高光譜;凸錐模型;端元提取;單體

    0 引言

    高光譜遙感圖像因具有較高光譜分辨率而成為當(dāng)今遙感領(lǐng)域研究的熱點,但其較低的空間分辨率和復(fù)雜的地物分布,導(dǎo)致影像中普遍存在混合像元,所以,基于高光譜圖像進行精確的地物識別和分類,首先要對混合像元進行分解,而端元提取是混合像元分解的關(guān)鍵環(huán)節(jié)[1]。目前,端元提取算法較多,總體可以分為3類:窮舉算法,投影算法與迭代算法。窮舉是一種重復(fù)型算法,它的基本思想是,對問題的所有可能狀態(tài)一一進行測試,直到找到解或全部可能的狀態(tài)都測試過為止,窮舉算法提取端元的精度高,但是算法運算量大,包括N-finder[2-3]算法,單體增長(simplex growing algorithm,SGA)[4-5]算法,其中N-FINDR算法,盡管無參數(shù)選擇效果好,但因其搜索凸面體最大體積需遍歷所有可能像元,計算量巨大。投影算法在進行端元提取時,計算速度相對較快,但是算法受噪聲影響較大,精度不高。例如:像元純凈指數(shù)(pixel purity index,PPI)[6-7]算法,頂點成分分析(vertex component analysis,VCA)[8-9]算法,它們都是基于投影的算法,只是投影的方式不同。迭代是不斷利用新值取代舊值,或者由舊值遞推出變量新值的過程,此法運算速度快,工作量少,精度取決于算法的本身。例如,序列最大角凸錐(sequential maximum angle convex cone,SMACC)[10-11]算法,它雖然計算速度快,但不能保證獲取端元的精度。

    為此,論文在兼顧效率和精度的同時,基于線性表示理論,提出了一種迭代算法,即基于單體擴張的端元提取算法。首先,隨機選取已知數(shù)目的像元作為初始端元集,用初始端元線性表示其他像元,得到一系列表示系數(shù),確定出負系數(shù)的絕對值最大的像元,表明此像元離單體的“距離”最遠,將其替換初始端元中的某一個,使得替換后的單體體積最大;重復(fù)上述過程直到無負系數(shù)。從而找到使凸面單體體積達到最大的像元集,即端元。

    1 端元提取的相關(guān)理論

    相關(guān)理論是基于理想情形下給出,即數(shù)據(jù)無噪聲,且完全滿足線性混合模型。

    1.1 端元提取的理論模型

    本文將基于線性混合模型展開討論。

    1.2 基于線性表示的端元判定

    在中任取個無關(guān)點1,2,…,,并記={1,2,…,},={1,2,…,}是端元集,則有以下結(jié)論。

    定理1的凸面單體()中的任意向量均可被1,2,…,線性表示,且系數(shù)和為1。

    定理2的凸面單體包含于的凸面單體中,即()ì()。

    定理1表明,對任意?,均可被唯一地線性表示,且系數(shù)之和為1,若系數(shù)全非負,則?(),否則?()。定理2表明,將?均由來線性表示,若所有的表示系數(shù)全非負,則就是端元。

    定理3={1,2,…,}是端元集的充分必要條件是高光譜向量集被線性表示的系數(shù)均非負。

    必要性,由定理1和定理2顯然可得。

    2 理想情形下端元提取算法

    基于上章給出的理論,端元提取的實現(xiàn)涉及3個內(nèi)容,分別是單體體積計算、線性表示的系數(shù)確定以及算法流程。

    2.1 單體體積計算

    根據(jù)文獻[13],以1,2,…,為頂點的-1維凸面單體的體積(())的計算公式為:

    式中:={2-1,3-1,…,-1}。

    2.2 線性表示的系數(shù)確定

    若端元個數(shù)為,則在理想情形下,ì(),高光譜向量集的秩為,由于1,2,…,線性無關(guān),所以?(=1,2,…,),均可由1,2,…,唯一地線性表示,={1,2,…,}其系數(shù)={1,2,…,},由線性方程組確定:

    (=1,2,…,) (3)

    2.3 端元提取算法流程

    算法采用迭代思想,首先從中任意選取個無關(guān)點作為初始端元={1,2,…,},再從中找出一個向量0來代替中的一個向量0,形成迭代過程。

    確定0的規(guī)則是使其離單體()盡量“遠”,考慮到當(dāng)用1,2,…,來表示x時,其負系數(shù)的絕對值越大,表明離單體()的“距離”越遠,所以可用線性表示的系數(shù)來確定0。確定被替換對象0的規(guī)則是使替換后的單體體積盡量大。具體流程如下:

    輸入:高光譜向量={1,2,…,},端元個數(shù);

    輸出:端元向量1,2,…,。

    第一步:在中任取個無關(guān)點1,2,…,

    第二步:對?,求線性方程組c得到系數(shù),并記=min(1,2,…,)(=1, 2,…,);

    第三步:若≥0(=1, 2,…,),則1,2,…,是端元向量并輸出,結(jié)束;

    盡管未能證明算法第三步中((0))>(())一定成立,但大量數(shù)值實驗均未出現(xiàn)算法失效的情形。

    3 非理想情形的端元提取算法

    上述理論和算法均基于理想情形,而實際數(shù)據(jù)一定不是理想狀況,疊加了各種噪聲,理論條件均不滿足,幾個關(guān)鍵的問題及算法改進如下:

    1)數(shù)據(jù)集不在同一-1維超平面內(nèi),而分布在超平面“兩側(cè)”附近,ì()只在“距離很近”的意義下成立。

    2)的秩遠遠大于,線性方程組(3)無解,但由于“近似”落在維超平面內(nèi),所以,可以求出線性方程組(3)的最小二乘解,作為“近似線性表示”系數(shù)。

    3)由于噪聲的存在,算法步驟第三步中不能保證((0))>(()),每次迭代得到的單體體積不單調(diào),不能保證算法收斂,出現(xiàn)算法失效的情況。對此,將算法第三步找出一個離()“最遠”點0的步驟,嵌套入一個內(nèi)層迭代,這樣,算法第三步改進如下:

    第三步:若≥0(=1,2,…,),則1,2,…,是端元向量并輸出,結(jié)束;

    理論上,可能存在算法失效的情形,但大量數(shù)值實驗均未出現(xiàn)失效情況。

    4 實驗與分析

    實驗首先通過人造數(shù)據(jù),在端元已知的情況下,運用算法提取端元,再對真實高光譜數(shù)據(jù)進行實驗。

    4.1 仿真實驗

    首先針對不添加噪聲的理想情形,再次針對添加噪聲的非理想情形。

    4.1.1 理想情形下

    為驗證本算法的有效性,實驗首先用人工合成的三維數(shù)據(jù),點(1,1.2,3),(0.4,2,1),(3,2,1)為理論端元,被隨機產(chǎn)生的歸一化的非負系數(shù)合成2000個點,即這2000個點在3個端點向量組成的三角形內(nèi),此時不添加任何的噪聲。實驗結(jié)果如圖1所示,圖中黑色部分圈出的點為提取出的“端元”,圖1(a)為算法初始時隨機選取的3個“端元”,圖1(b)為算法經(jīng)過第一次迭代,可觀察到此時提取出其中的一個端元,圖1(d)為算法最終迭代,提取出的端元,恰為理論端元。

    圖1 理想情形下的端元提取過程圖

    4.1.2 非理想情形下

    實驗依舊用人工合成的三維數(shù)據(jù),點(1,1.2,3),(0.4,2,1),(3,2,1)為理論端元,被隨機產(chǎn)生的歸一化的非負系數(shù)合成2000個點,再附加0均值的高斯噪聲,實驗結(jié)果如圖2所示,圖中黑色部分,圈出的點為提取出的“端元”,圖2(a)為算法初始時隨機選取的3個“端元”,圖2(b)為算法第一次迭代,可觀察到此時提取出其中的一個端元,圖2(d)為算法最終迭代提取出的端元,恰為理論端元。計算出本算法提取的端元單體體積最大,如表1。

    圖2 非理想情形下的端元提取過程圖

    表1 三種算法所得端元單體的體積

    當(dāng)合成數(shù)據(jù)為20維,3個端元時,被隨機產(chǎn)生歸一化的非負系數(shù)合成10000個像元,分別加入信噪比為15dB,25dB,35dB,45dB,55dB的高斯白噪聲時,用3種方法提取出端元,計算平均光譜角距離(mSAD)的值,得到的結(jié)果如圖3(a)所示,本文算法對噪聲的敏感性低。

    當(dāng)端元數(shù)目分別為3,4,5,6,7個時,維數(shù)為20,分別被隨機產(chǎn)生歸一化的非負系數(shù)合成10000個像元,加入信噪比為35dB的高斯白噪聲,用3種算法提取端元計算提取的端元與真實端元的mSAD值,如圖3(b)所示,本文算法得到mSAD值為0,說明端元提取精度最好的是本文算法,其次是VCA算法,PPI算法最差。

    4.2 高光譜數(shù)據(jù)的算法實驗分析

    本文用美國內(nèi)華達州Cuprite地區(qū)的AVIRIS數(shù)據(jù)進行分解實驗,Cuprite地區(qū)AVIRIS原始數(shù)據(jù),共224個波段,光譜分辨率為10nm,尺寸為350×350像元,刪除噪聲較大和光譜吸收較大的波段,本文選擇172~221光譜區(qū)間的50個波段用于算法測試,已知端元數(shù)目為8[14],算法迭代16次,便可提取出最終端元,用提取出的端元對高光譜影像進行豐度反演,反演時運用基于幾何距離的全約束混合像元分解算法(DGAE)[15],得到的豐度圖如圖4(a)~圖4(h)所示,將各豐度圖與圖4(j)的調(diào)查分類結(jié)果圖比較,從而確定端元對應(yīng)的礦物。圖4(i)為端元的位置分布。圖4(k)為用于處理的區(qū)域的AVIRIS假彩色合成圖像。最終將本文算法與VCA,PPI算法得到的端元體積進行比較,如表2所示,本文算法得到的體積最大。

    圖3 3種算法性能比較

    表2 3種算法所得端元單體的體積

    為了進一步比較,我們將USGS庫中對應(yīng)的礦物光譜作為參考光譜,將參考光譜在相應(yīng)的波段范圍內(nèi)能量求積分,轉(zhuǎn)換為光譜向量,并求取解混所得的端元與它們之間的光譜角,如表3所示,可以看出,用本文方法提取的端元與USGS光譜庫中的標準光譜相似度較高,總體上優(yōu)于VCA與PPI算法提取的端元光譜。

    5 結(jié)論

    本論文是在線性混合模型的基礎(chǔ)上,假設(shè)純凈像元存在的前提下提出的,是一種基于單體擴張的端元提取算法,它的主要思想就是基于線性表示的端元判定,給出了端元判別的定理,并將定理進行證明,是另外一種迭代算法,并且,通過迭代,優(yōu)化其時間效率。文中用人工仿真的數(shù)據(jù)與真實的高光譜數(shù)據(jù)進行實驗,表明本算法精度高,穩(wěn)定性好。算法是在端元數(shù)目已知的情況下給出的,后續(xù)工作將針對端元數(shù)目未知時進行研究。

    圖4 美國Cuprite地區(qū)AVIRIS數(shù)據(jù)端元分解結(jié)果(350像元×350像元)

    表3 提取的端元與標準光譜庫中的標準光譜的光譜角距離

    [1] 李二森, 朱述龍, 周曉明, 等. 高光譜圖像端元提取算法研究進展與比較[J]. 遙感學(xué)報, 2011,15(4): 659-679.

    LI Ersen, ZHU Shulong, ZHOU Xiaoming, et al. The development and comparison of endmember extraction algorithms using hyperspectral imagery[J]., 2011, 15(4): 659-679.

    [2] JI Luyan, GENG Xiurui, Sun Kang, et al. Modified N-FINDR endmember extraction algorithm for remote-sensing imagery[J]., 2015, 36(8): 2148-2162.

    [3] 楊可明, 魏華鋒, 劉飛, 等. 以光譜信息熵改進的N-FINDR端元提取算法[J]. 地球信息科學(xué)學(xué)報, 2015, 17(8): 980-985.

    YANG Keming, WEI Huafeng, LIU Fei, et al. Improved N-FINDR algorithm on hyperspectral endmember extraction based on spectral Shannon entropy[J]., 2015, 17(8): 980-985.

    [4] ZHAO L, FAN M, WANG L. Fast implementation of linear and nonlinear simplex growing algorithm for hyperspectral endmember extraction[J].,2015, 126(23): 4072-4077.

    [5] LIU Junmin, ZHANG Jiangshe. A new maximum simplex volume method based on householder transformation for endmember extraction[J]., 2012, 50(1): 104-111.

    [6] Rob Heylen, Paul Scheunders. Multidimensional pixel purity index for convex Hull estimation and endmember extraction[J]., 2013, 51(7): 4059-4069.

    [7] 徐軍, 徐富紅, 蔡體健, 等. 一種基于最大距離的純像元指數(shù)端元提取算法[J].地球信息科學(xué)學(xué)報, 2015, 17(1): 86-90.

    XU Jun, XU Fuhong, CAI Tijian, et al. A novel pure pixel endmember extraction algorithm based on the maximum distance[J]., 2015, 17(1): 86-90.

    [8] Jose M, Rodriguez Alves, Jose M, et al. Vertex component analysis Gpu-based implementation for hyperspectral unimixing[C]//, 2012: 10.1109/WHISPERS.2012.6874337.

    [9] 何曉寧, 曹建農(nóng), 高坡, 等. 基于頂點成分分析法的端元提取改進算法[J]. 測繪通報, 2013 (7): 30-34.

    HE Xiaoning, CAO Jiannong, GAO Po, et al. An improved algorithm of endmember extraction based on vertex component analysis[J]., 2013(7): 30-34.

    [10] John Gruninger, Anthony J Ratkowski, Michael L Hoke. The sequential maximum angle Convex Cone(SMACC) endmember model[C]//,[C]//2004, 5425: doi:10.1117/12.543794.

    [11] LIU C, LI J, WANG G. Endmember extraction algorithm for hyperspectral image based on PCA-SMACC[C]//,, 2014, 9142:91421A.

    [12] Ambikapathi A M, Chan T H, Ma W K. A robust alternating volume maximization algorithm for endmember extraction in hyperspectral images[C]// 2on, 2010: 10.1109/WHISPERS.2010.5594862.

    [13] SUN Kang, GENG Xiurui, WANG Panshi, et al. A fast endmember extraction algorithm based on gram determinant[J].,2014, 11(6): 1124-1127.

    [14] 田玉剛, 楊貴. 端元快速提取的光譜梯度特征搜索法[J]. 測繪學(xué)報, 2015, 44(2): 214-219.

    TIAN Yugang, YANG Gui. A fast endmember extraction algorithmusing spectrum gradient features[J]., 2015,44(2): 214-219.

    [15] PU H, XIA W, WANG B, et al. A fully constrained linear spectral unmixing algorithm based on distance geometry[J].,2014, 52(2): 1157-1176.

    Endmember Extraction Algorithm Based on the Simplex Expansion

    DONG Anguo,HAN Xue,GONG Wenjuan

    (School of Science, Chang’an University, Xi’an 710064, China)

    Endmember extraction is a key step in unmixing hyperspectral mixed pixels. In order to extract endmembers of hyperspectral image, this paper proves that if a vector in vitro and vivo is represented by a vertex vector, there will be a negative coefficient based on the theory of linear representation and the theory of convex cone model. The theory gives a necessary and sufficient condition for the endmembers identification. A new iterative algorithm is proposed under non-ideal situation. The experimental results show that the precision of the endmember extraction using the algorithm proposed in this paper is better than VCA algorithm. And this algorithm has better efficiency than the search algorithm. It has good stability and a low sensitivity to noise.

    hyperspectral,convex cone,endmember extraction,simplex

    TP7

    A

    1001-8891(2016)11-0947-06

    2016-05-24;

    2016-08-05.

    董安國(1964-),男,浙江象山人,教授,碩士生導(dǎo)師,主要研究領(lǐng)域為數(shù)值代數(shù),數(shù)字圖像處理,E-mail:donganguo@chd.edu.cn。

    國家自然科學(xué)基金項目(41571346),國家自然科學(xué)基金項目(40971217),國家自然科學(xué)基金項目(11201038)。

    猜你喜歡
    端元情形單體
    現(xiàn)代黃河三角洲沉積物粒度特征及其來源
    避免房地產(chǎn)繼承糾紛的十二種情形
    四種情形拖欠勞動報酬構(gòu)成“拒不支付”犯罪
    公民與法治(2020年4期)2020-05-30 12:31:34
    南昌地區(qū)不透水面遙感估算研究
    單體光電產(chǎn)品檢驗驗收方案問題探討
    兩種基于異常權(quán)重的N-FINDR端元提取算法
    出借車輛,五種情形下須擔(dān)責(zé)
    公民與法治(2016年9期)2016-05-17 04:12:18
    基于Gram行列式的快速端元提取方法
    相變大單體MPEGMA的制備與性能
    巨無霸式醫(yī)療單體的選擇
    久久人人爽av亚洲精品天堂| 国内毛片毛片毛片毛片毛片| 黄片大片在线免费观看| 欧美日韩亚洲国产一区二区在线观看 | 欧美精品av麻豆av| 老司机午夜福利在线观看视频 | videos熟女内射| 国产精品久久久av美女十八| 伊人久久大香线蕉亚洲五| 国产高清videossex| 欧美日韩一级在线毛片| 久久性视频一级片| 69精品国产乱码久久久| 免费女性裸体啪啪无遮挡网站| 国产精品偷伦视频观看了| 各种免费的搞黄视频| 国产一区二区 视频在线| 9191精品国产免费久久| 亚洲精品中文字幕在线视频| 日韩精品免费视频一区二区三区| av有码第一页| a级毛片在线看网站| 两性夫妻黄色片| 亚洲国产av影院在线观看| 欧美少妇被猛烈插入视频| 美女福利国产在线| 欧美成人午夜精品| 国产xxxxx性猛交| 国产一区二区三区综合在线观看| 人成视频在线观看免费观看| 18禁黄网站禁片午夜丰满| 国产高清国产精品国产三级| 99久久综合免费| 亚洲欧美激情在线| tube8黄色片| 国产欧美日韩一区二区三 | kizo精华| 亚洲精品美女久久久久99蜜臀| 精品人妻一区二区三区麻豆| 搡老岳熟女国产| 成人国语在线视频| svipshipincom国产片| 国产成人a∨麻豆精品| 亚洲精品成人av观看孕妇| 精品亚洲成a人片在线观看| 欧美老熟妇乱子伦牲交| 亚洲人成77777在线视频| 久久久久国内视频| 中文字幕色久视频| 一区二区三区四区激情视频| 国产精品99久久99久久久不卡| 老熟妇乱子伦视频在线观看 | 国产精品99久久99久久久不卡| 国产欧美日韩一区二区三区在线| 91成人精品电影| 9热在线视频观看99| 侵犯人妻中文字幕一二三四区| 亚洲三区欧美一区| 天天躁夜夜躁狠狠躁躁| 国产又爽黄色视频| 久久中文看片网| 90打野战视频偷拍视频| 男女免费视频国产| 国产成人精品无人区| 在线看a的网站| 免费高清在线观看视频在线观看| 一区二区av电影网| 午夜成年电影在线免费观看| 高潮久久久久久久久久久不卡| 美国免费a级毛片| 99热全是精品| 国产又爽黄色视频| 亚洲久久久国产精品| 最新的欧美精品一区二区| 午夜久久久在线观看| 一边摸一边抽搐一进一出视频| 精品国产一区二区三区四区第35| 国产福利在线免费观看视频| 成人三级做爰电影| 亚洲熟女毛片儿| 欧美黄色片欧美黄色片| 一区二区av电影网| 9热在线视频观看99| 80岁老熟妇乱子伦牲交| 亚洲人成电影观看| 高潮久久久久久久久久久不卡| 欧美日韩亚洲高清精品| 日韩熟女老妇一区二区性免费视频| 一区二区av电影网| 欧美人与性动交α欧美软件| 少妇被粗大的猛进出69影院| 91国产中文字幕| 亚洲色图综合在线观看| 亚洲少妇的诱惑av| 肉色欧美久久久久久久蜜桃| 久久ye,这里只有精品| 欧美成人午夜精品| 一个人免费在线观看的高清视频 | 另类精品久久| 日韩视频一区二区在线观看| 美女扒开内裤让男人捅视频| 亚洲精品中文字幕在线视频| 中文字幕人妻丝袜制服| 欧美大码av| 动漫黄色视频在线观看| 亚洲男人天堂网一区| 成年人黄色毛片网站| 亚洲 欧美一区二区三区| 午夜精品国产一区二区电影| 狠狠狠狠99中文字幕| 国产激情久久老熟女| 精品少妇一区二区三区视频日本电影| 啦啦啦 在线观看视频| 精品国产一区二区久久| 丰满少妇做爰视频| 国产精品国产av在线观看| 国产黄频视频在线观看| 国产成人精品久久二区二区免费| 精品一品国产午夜福利视频| 国产色视频综合| 99国产精品99久久久久| 久久精品熟女亚洲av麻豆精品| 久久天躁狠狠躁夜夜2o2o| 黄片播放在线免费| 91字幕亚洲| 欧美日韩亚洲综合一区二区三区_| 可以免费在线观看a视频的电影网站| 久久久久网色| 老汉色av国产亚洲站长工具| 桃花免费在线播放| 午夜两性在线视频| 亚洲九九香蕉| 免费久久久久久久精品成人欧美视频| 久久久久久久大尺度免费视频| 国产精品久久久人人做人人爽| 青草久久国产| 日本av手机在线免费观看| 男人舔女人的私密视频| 一级片免费观看大全| 91大片在线观看| 精品国产超薄肉色丝袜足j| 亚洲精品粉嫩美女一区| 可以免费在线观看a视频的电影网站| 国产精品免费视频内射| 国产成人精品在线电影| 欧美大码av| 亚洲专区中文字幕在线| 国产在线视频一区二区| 狠狠婷婷综合久久久久久88av| 久久国产精品影院| 亚洲国产欧美日韩在线播放| 精品亚洲成a人片在线观看| 在线永久观看黄色视频| 香蕉丝袜av| 日韩有码中文字幕| 久久久久精品国产欧美久久久 | 成人国产一区最新在线观看| 深夜精品福利| 视频在线观看一区二区三区| 亚洲精品国产一区二区精华液| 亚洲一区中文字幕在线| 久久人人97超碰香蕉20202| 亚洲伊人色综图| 美女高潮喷水抽搐中文字幕| 国产精品.久久久| 操出白浆在线播放| 国产欧美亚洲国产| 青草久久国产| 免费日韩欧美在线观看| 精品少妇黑人巨大在线播放| 狂野欧美激情性xxxx| 亚洲熟女毛片儿| 国产精品一区二区在线观看99| av网站在线播放免费| 国产精品国产三级国产专区5o| 女性生殖器流出的白浆| 制服人妻中文乱码| 免费黄频网站在线观看国产| 国产成人影院久久av| 免费一级毛片在线播放高清视频 | 亚洲av欧美aⅴ国产| 高潮久久久久久久久久久不卡| 在线观看人妻少妇| 欧美日韩一级在线毛片| 国产精品99久久99久久久不卡| 久久久久精品人妻al黑| 老司机午夜福利在线观看视频 | 丝袜喷水一区| 国产极品粉嫩免费观看在线| 亚洲熟女毛片儿| 久久国产精品大桥未久av| 久久久精品94久久精品| 国产亚洲av片在线观看秒播厂| 久久免费观看电影| 日韩免费高清中文字幕av| 国产一卡二卡三卡精品| 日本撒尿小便嘘嘘汇集6| 侵犯人妻中文字幕一二三四区| 久久精品人人爽人人爽视色| 日本猛色少妇xxxxx猛交久久| 午夜视频精品福利| 久久99一区二区三区| 国产成人a∨麻豆精品| 欧美亚洲日本最大视频资源| 亚洲av欧美aⅴ国产| 啦啦啦在线免费观看视频4| 在线观看一区二区三区激情| 人人妻人人澡人人看| 亚洲国产精品成人久久小说| 大香蕉久久网| 久久国产精品人妻蜜桃| 欧美97在线视频| 99香蕉大伊视频| www日本在线高清视频| 99精品欧美一区二区三区四区| 桃花免费在线播放| 午夜精品国产一区二区电影| 亚洲五月婷婷丁香| 色老头精品视频在线观看| 日韩,欧美,国产一区二区三区| 日韩中文字幕视频在线看片| 999久久久国产精品视频| 纵有疾风起免费观看全集完整版| 黑人操中国人逼视频| 他把我摸到了高潮在线观看 | 亚洲专区字幕在线| 亚洲精品国产精品久久久不卡| 亚洲avbb在线观看| 涩涩av久久男人的天堂| 国产国语露脸激情在线看| 丝袜脚勾引网站| 五月开心婷婷网| 欧美中文综合在线视频| 少妇精品久久久久久久| 中文字幕人妻熟女乱码| 亚洲成人免费电影在线观看| av超薄肉色丝袜交足视频| 精品国产国语对白av| 国产欧美日韩一区二区三 | 天天躁夜夜躁狠狠躁躁| 中文欧美无线码| 亚洲欧美清纯卡通| 国产一区二区三区在线臀色熟女 | 黑丝袜美女国产一区| 两个人看的免费小视频| 少妇被粗大的猛进出69影院| av网站在线播放免费| 男人舔女人的私密视频| 免费观看a级毛片全部| 下体分泌物呈黄色| 51午夜福利影视在线观看| 一级毛片女人18水好多| 亚洲精品一二三| 一本色道久久久久久精品综合| 婷婷色av中文字幕| 男女边摸边吃奶| 一本综合久久免费| 老司机影院毛片| 国产精品一区二区免费欧美 | 久久青草综合色| av欧美777| 99热网站在线观看| 老熟女久久久| bbb黄色大片| 国产亚洲av高清不卡| 狠狠婷婷综合久久久久久88av| 啦啦啦在线免费观看视频4| 考比视频在线观看| 国产成人啪精品午夜网站| 国产亚洲精品第一综合不卡| 一区在线观看完整版| 一级a爱视频在线免费观看| 日韩视频一区二区在线观看| 久久久久久久精品精品| 亚洲第一青青草原| 黑人猛操日本美女一级片| 国产亚洲精品第一综合不卡| 少妇人妻久久综合中文| 999久久久国产精品视频| 黄片大片在线免费观看| 99热全是精品| 久久精品成人免费网站| 午夜激情av网站| 99精国产麻豆久久婷婷| 国产高清videossex| 一级a爱视频在线免费观看| 久久久久久久久免费视频了| 热re99久久精品国产66热6| 人人妻人人爽人人添夜夜欢视频| 欧美亚洲 丝袜 人妻 在线| 久久午夜综合久久蜜桃| 欧美 亚洲 国产 日韩一| 精品一品国产午夜福利视频| 黑丝袜美女国产一区| 国产精品亚洲av一区麻豆| 69精品国产乱码久久久| 最黄视频免费看| 青青草视频在线视频观看| 另类精品久久| 99re6热这里在线精品视频| bbb黄色大片| 五月天丁香电影| 99精品欧美一区二区三区四区| 老司机午夜十八禁免费视频| www日本在线高清视频| 美女高潮喷水抽搐中文字幕| 日韩制服骚丝袜av| 国产精品99久久99久久久不卡| 一二三四社区在线视频社区8| 热99re8久久精品国产| 在线十欧美十亚洲十日本专区| 午夜福利视频在线观看免费| 久久久久久久久久久久大奶| 淫妇啪啪啪对白视频 | 亚洲 国产 在线| 狠狠狠狠99中文字幕| 99久久99久久久精品蜜桃| 欧美日本中文国产一区发布| 高潮久久久久久久久久久不卡| 18禁国产床啪视频网站| 成人影院久久| 国产精品香港三级国产av潘金莲| 亚洲综合色网址| 久久久国产欧美日韩av| 夜夜夜夜夜久久久久| 91字幕亚洲| 亚洲人成电影免费在线| 免费在线观看影片大全网站| 后天国语完整版免费观看| 精品少妇一区二区三区视频日本电影| 国产黄色免费在线视频| 亚洲成人免费电影在线观看| 国产精品成人在线| 老汉色∧v一级毛片| 亚洲五月婷婷丁香| 国产精品久久久久久精品古装| 亚洲五月婷婷丁香| 欧美日韩黄片免| 久久性视频一级片| 老汉色∧v一级毛片| 国产精品成人在线| 免费在线观看影片大全网站| 亚洲人成电影免费在线| 老熟女久久久| 最新在线观看一区二区三区| 国产精品久久久av美女十八| 成人三级做爰电影| 久久精品aⅴ一区二区三区四区| 老司机午夜福利在线观看视频 | 久久天堂一区二区三区四区| 国产一区二区 视频在线| 国产91精品成人一区二区三区 | 99热全是精品| 国产色视频综合| 激情视频va一区二区三区| 午夜福利在线观看吧| 午夜激情av网站| 美女大奶头黄色视频| 欧美日韩精品网址| 色精品久久人妻99蜜桃| 欧美成人午夜精品| 久久性视频一级片| 久久久久久人人人人人| 久久这里只有精品19| 亚洲 国产 在线| 99久久国产精品久久久| 国产日韩欧美亚洲二区| 成年女人毛片免费观看观看9 | 欧美 亚洲 国产 日韩一| 91av网站免费观看| 亚洲午夜精品一区,二区,三区| 狠狠狠狠99中文字幕| 久久久久网色| 亚洲一区二区三区欧美精品| 午夜两性在线视频| 欧美人与性动交α欧美精品济南到| 欧美 日韩 精品 国产| 国产一区二区三区av在线| 日韩一区二区三区影片| 久久九九热精品免费| 国产视频一区二区在线看| 最新的欧美精品一区二区| 日韩中文字幕欧美一区二区| 69精品国产乱码久久久| 乱人伦中国视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲少妇的诱惑av| av片东京热男人的天堂| 天天影视国产精品| 欧美另类一区| 国产xxxxx性猛交| 亚洲欧美成人综合另类久久久| 一边摸一边抽搐一进一出视频| cao死你这个sao货| 少妇 在线观看| 麻豆av在线久日| 日韩一区二区三区影片| 国产熟女午夜一区二区三区| 韩国精品一区二区三区| 国产成人欧美在线观看 | 精品人妻1区二区| 桃花免费在线播放| 日韩精品免费视频一区二区三区| 搡老乐熟女国产| 精品乱码久久久久久99久播| xxxhd国产人妻xxx| 国产男女超爽视频在线观看| www日本在线高清视频| 一级片免费观看大全| 久久人人爽人人片av| 国产精品成人在线| 成人国产一区最新在线观看| 18禁裸乳无遮挡动漫免费视频| 国产一区有黄有色的免费视频| 国产亚洲欧美精品永久| 欧美激情 高清一区二区三区| 人人妻人人添人人爽欧美一区卜| 国产三级黄色录像| 高清欧美精品videossex| 国产成人精品久久二区二区91| a级毛片黄视频| 视频在线观看一区二区三区| 飞空精品影院首页| 精品久久久久久电影网| 欧美在线黄色| 久久久精品94久久精品| 老汉色av国产亚洲站长工具| 咕卡用的链子| 老鸭窝网址在线观看| 又紧又爽又黄一区二区| 国产一区二区三区av在线| 91字幕亚洲| 国产精品亚洲av一区麻豆| 亚洲国产看品久久| 久久午夜综合久久蜜桃| 国产精品国产av在线观看| 国产男人的电影天堂91| 91字幕亚洲| 午夜福利在线免费观看网站| 黄色怎么调成土黄色| 日韩欧美一区视频在线观看| 国产极品粉嫩免费观看在线| 久久久水蜜桃国产精品网| 亚洲欧美一区二区三区久久| 免费在线观看完整版高清| 777米奇影视久久| 老熟女久久久| av福利片在线| 12—13女人毛片做爰片一| 1024视频免费在线观看| 国产老妇伦熟女老妇高清| 男女之事视频高清在线观看| 国产1区2区3区精品| 99re6热这里在线精品视频| 日本av手机在线免费观看| a级片在线免费高清观看视频| 国产精品欧美亚洲77777| 国产日韩欧美在线精品| 十分钟在线观看高清视频www| 国产在线一区二区三区精| 久久性视频一级片| 一级黄色大片毛片| 亚洲黑人精品在线| 搡老熟女国产l中国老女人| 国产亚洲欧美在线一区二区| av在线播放精品| 国产欧美亚洲国产| 国产一区二区 视频在线| 日韩欧美一区二区三区在线观看 | 亚洲精品成人av观看孕妇| 两性午夜刺激爽爽歪歪视频在线观看 | 国产人伦9x9x在线观看| 18禁观看日本| 91国产中文字幕| 亚洲熟女毛片儿| 久久久久国内视频| 99热全是精品| 欧美精品人与动牲交sv欧美| 国产精品熟女久久久久浪| 性色av乱码一区二区三区2| cao死你这个sao货| 午夜91福利影院| 国产精品二区激情视频| 精品一区在线观看国产| 青春草视频在线免费观看| 日韩一卡2卡3卡4卡2021年| 欧美成人午夜精品| 大型av网站在线播放| 国产亚洲精品一区二区www | 国产一区二区三区在线臀色熟女 | 丝袜在线中文字幕| 麻豆av在线久日| 欧美性长视频在线观看| 美女脱内裤让男人舔精品视频| 啦啦啦在线免费观看视频4| 19禁男女啪啪无遮挡网站| 人人妻人人添人人爽欧美一区卜| 国产精品久久久久久精品古装| 建设人人有责人人尽责人人享有的| 高清黄色对白视频在线免费看| 欧美精品啪啪一区二区三区 | 高清黄色对白视频在线免费看| 久久精品亚洲熟妇少妇任你| 中文字幕另类日韩欧美亚洲嫩草| 久久精品国产综合久久久| 亚洲九九香蕉| 亚洲情色 制服丝袜| 精品福利观看| 国产精品 国内视频| 另类精品久久| 国产极品粉嫩免费观看在线| av线在线观看网站| 婷婷丁香在线五月| 久久久欧美国产精品| av又黄又爽大尺度在线免费看| 精品视频人人做人人爽| 欧美国产精品一级二级三级| 久久天堂一区二区三区四区| 日本一区二区免费在线视频| 天天躁夜夜躁狠狠躁躁| 久久久久久亚洲精品国产蜜桃av| 亚洲国产毛片av蜜桃av| 一本—道久久a久久精品蜜桃钙片| 精品人妻1区二区| 日韩制服骚丝袜av| 精品亚洲乱码少妇综合久久| 视频区图区小说| 伊人久久大香线蕉亚洲五| 18禁黄网站禁片午夜丰满| 老熟女久久久| 亚洲伊人久久精品综合| 午夜免费成人在线视频| 色婷婷av一区二区三区视频| 亚洲成人免费av在线播放| 亚洲av男天堂| 91成人精品电影| 日韩中文字幕视频在线看片| 国产又爽黄色视频| 婷婷丁香在线五月| 天天影视国产精品| 午夜久久久在线观看| 精品人妻在线不人妻| 国产人伦9x9x在线观看| 国产在线视频一区二区| 午夜日韩欧美国产| 一个人免费在线观看的高清视频 | 啦啦啦啦在线视频资源| 国产黄频视频在线观看| 久久久久久人人人人人| 国产麻豆69| 嫩草影视91久久| 国产97色在线日韩免费| av在线老鸭窝| 2018国产大陆天天弄谢| 日韩欧美一区视频在线观看| 黑人猛操日本美女一级片| 国产91精品成人一区二区三区 | 天天躁日日躁夜夜躁夜夜| 欧美xxⅹ黑人| 亚洲欧美日韩另类电影网站| 久久久久国内视频| 丝袜在线中文字幕| 性少妇av在线| 丁香六月天网| 午夜成年电影在线免费观看| 人妻一区二区av| 丰满迷人的少妇在线观看| 国产亚洲欧美在线一区二区| 一本综合久久免费| 伊人亚洲综合成人网| 麻豆乱淫一区二区| 别揉我奶头~嗯~啊~动态视频 | 国产伦人伦偷精品视频| 亚洲一区中文字幕在线| 国产黄频视频在线观看| 在线亚洲精品国产二区图片欧美| 在线观看免费午夜福利视频| 国产亚洲精品一区二区www | 亚洲精华国产精华精| 韩国高清视频一区二区三区| 久久这里只有精品19| 在线观看一区二区三区激情| 国产成人啪精品午夜网站| 亚洲第一青青草原| 精品国内亚洲2022精品成人 | 久久久久国产一级毛片高清牌| 国产在线观看jvid| 成人三级做爰电影| 人人妻人人添人人爽欧美一区卜| 精品卡一卡二卡四卡免费| 啦啦啦 在线观看视频| 中国国产av一级| 18禁国产床啪视频网站| 日韩有码中文字幕| 国产免费一区二区三区四区乱码| 一个人免费在线观看的高清视频 | 男人操女人黄网站| 性色av一级| 日韩大码丰满熟妇| 飞空精品影院首页| 最黄视频免费看| av在线老鸭窝| 精品一区二区三区四区五区乱码| 十八禁高潮呻吟视频| 亚洲精品乱久久久久久| 国产av一区二区精品久久| 脱女人内裤的视频| 久久女婷五月综合色啪小说| 老司机影院成人| 欧美日韩视频精品一区| 亚洲精品国产av蜜桃| 日本vs欧美在线观看视频| 亚洲国产看品久久| 狂野欧美激情性xxxx| 另类亚洲欧美激情|