• <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ī)療單體的選擇
    女人高潮潮喷娇喘18禁视频| 色综合婷婷激情| 性色av乱码一区二区三区2| 97人妻天天添夜夜摸| 在线看a的网站| 侵犯人妻中文字幕一二三四区| 欧美亚洲 丝袜 人妻 在线| 一区二区日韩欧美中文字幕| ponron亚洲| 色婷婷av一区二区三区视频| 欧美黑人欧美精品刺激| av国产精品久久久久影院| 久久亚洲真实| 丁香欧美五月| 国产激情欧美一区二区| 国产精品av久久久久免费| 亚洲熟女毛片儿| 美女扒开内裤让男人捅视频| 啦啦啦 在线观看视频| 老汉色∧v一级毛片| 精品第一国产精品| 欧美精品人与动牲交sv欧美| 免费黄频网站在线观看国产| 免费在线观看日本一区| tocl精华| 丁香欧美五月| 欧美乱码精品一区二区三区| 在线观看午夜福利视频| 免费黄频网站在线观看国产| 久久婷婷成人综合色麻豆| 女性生殖器流出的白浆| 国产成+人综合+亚洲专区| 美女午夜性视频免费| 国产aⅴ精品一区二区三区波| 久久天躁狠狠躁夜夜2o2o| 久久久久久免费高清国产稀缺| 久久ye,这里只有精品| 亚洲黑人精品在线| 欧美亚洲日本最大视频资源| 无人区码免费观看不卡| 日韩精品免费视频一区二区三区| 亚洲第一欧美日韩一区二区三区| 老司机影院毛片| 999久久久精品免费观看国产| 高清av免费在线| 在线天堂中文资源库| 丝袜美足系列| 免费看十八禁软件| 首页视频小说图片口味搜索| 久久天堂一区二区三区四区| 亚洲自偷自拍图片 自拍| 波多野结衣一区麻豆| 亚洲精品美女久久av网站| 亚洲精品在线美女| www.熟女人妻精品国产| 久久天躁狠狠躁夜夜2o2o| 国产在线一区二区三区精| 日日摸夜夜添夜夜添小说| 国产精品美女特级片免费视频播放器 | 又黄又粗又硬又大视频| 国产亚洲欧美98| 天堂中文最新版在线下载| 天天躁夜夜躁狠狠躁躁| 日韩欧美在线二视频 | 日本a在线网址| 国产精品99久久99久久久不卡| 久久草成人影院| 天天躁狠狠躁夜夜躁狠狠躁| 一级毛片女人18水好多| a级毛片在线看网站| 美女午夜性视频免费| 精品久久久精品久久久| bbb黄色大片| 国产欧美日韩精品亚洲av| 亚洲欧美一区二区三区久久| 欧美色视频一区免费| 国产精品国产av在线观看| 欧美不卡视频在线免费观看 | 啦啦啦在线免费观看视频4| 欧美精品人与动牲交sv欧美| 日韩三级视频一区二区三区| 不卡一级毛片| 啦啦啦在线免费观看视频4| 18禁美女被吸乳视频| 99国产极品粉嫩在线观看| 9191精品国产免费久久| 国产麻豆69| 亚洲视频免费观看视频| 岛国在线观看网站| 在线播放国产精品三级| 国产成人精品久久二区二区91| 日韩 欧美 亚洲 中文字幕| 欧美日韩瑟瑟在线播放| 精品久久蜜臀av无| 国产一区二区三区综合在线观看| 久久婷婷成人综合色麻豆| 国产激情欧美一区二区| 欧美日韩成人在线一区二区| avwww免费| 两个人看的免费小视频| 国产欧美日韩一区二区精品| 久久精品人人爽人人爽视色| 国产又色又爽无遮挡免费看| 黑人巨大精品欧美一区二区mp4| 成熟少妇高潮喷水视频| 久久亚洲真实| av免费在线观看网站| 亚洲成a人片在线一区二区| 欧美激情高清一区二区三区| 纯流量卡能插随身wifi吗| 午夜91福利影院| 男人操女人黄网站| 亚洲av电影在线进入| 大型黄色视频在线免费观看| 欧美性长视频在线观看| 亚洲avbb在线观看| 精品久久久久久久久久免费视频 | 亚洲一码二码三码区别大吗| 欧美黑人欧美精品刺激| a级毛片黄视频| 精品久久蜜臀av无| 交换朋友夫妻互换小说| 国产精品久久久人人做人人爽| 9热在线视频观看99| 精品国产乱子伦一区二区三区| 一进一出抽搐动态| 国产有黄有色有爽视频| 精品国产一区二区三区四区第35| 丝袜美足系列| 婷婷丁香在线五月| 男女下面插进去视频免费观看| 国产xxxxx性猛交| 欧美日韩国产mv在线观看视频| 欧美中文综合在线视频| 91国产中文字幕| 国产精品自产拍在线观看55亚洲 | 亚洲国产欧美网| 天堂动漫精品| 亚洲av成人不卡在线观看播放网| 精品卡一卡二卡四卡免费| 欧美日韩视频精品一区| 国产av精品麻豆| 久久天堂一区二区三区四区| 日韩欧美国产一区二区入口| 香蕉久久夜色| 日日夜夜操网爽| 搡老熟女国产l中国老女人| 两人在一起打扑克的视频| 欧美激情高清一区二区三区| 老汉色∧v一级毛片| 啦啦啦视频在线资源免费观看| 色尼玛亚洲综合影院| 国产极品粉嫩免费观看在线| 天天操日日干夜夜撸| 在线永久观看黄色视频| 免费在线观看日本一区| 熟女少妇亚洲综合色aaa.| 精品久久久久久电影网| 黄频高清免费视频| 女性生殖器流出的白浆| 午夜免费观看网址| 麻豆成人av在线观看| 成人精品一区二区免费| 18禁裸乳无遮挡动漫免费视频| 99热只有精品国产| 99国产精品一区二区三区| 国产精品 欧美亚洲| 成人手机av| 人妻一区二区av| 亚洲中文日韩欧美视频| 女性生殖器流出的白浆| 亚洲精品成人av观看孕妇| 亚洲精品av麻豆狂野| 国产精品秋霞免费鲁丝片| 一进一出好大好爽视频| 一区二区三区精品91| 婷婷成人精品国产| 大码成人一级视频| 18禁黄网站禁片午夜丰满| 亚洲av电影在线进入| 高清欧美精品videossex| 国产日韩一区二区三区精品不卡| 一边摸一边抽搐一进一小说 | 亚洲成人免费av在线播放| 在线av久久热| 国产蜜桃级精品一区二区三区 | 亚洲国产精品sss在线观看 | 性色av乱码一区二区三区2| 黑人巨大精品欧美一区二区蜜桃| 免费在线观看视频国产中文字幕亚洲| 老司机午夜福利在线观看视频| 纯流量卡能插随身wifi吗| 国产亚洲欧美精品永久| 色94色欧美一区二区| 极品少妇高潮喷水抽搐| 人人澡人人妻人| 一进一出抽搐动态| 日韩三级视频一区二区三区| 国产成人欧美| 日日爽夜夜爽网站| 精品无人区乱码1区二区| 国产精品 国内视频| 亚洲精品在线观看二区| 久久精品亚洲精品国产色婷小说| 久久久久久久久免费视频了| 99久久精品国产亚洲精品| 99国产综合亚洲精品| 午夜福利一区二区在线看| 精品亚洲成国产av| 美女国产高潮福利片在线看| 精品国产国语对白av| 欧美日韩亚洲国产一区二区在线观看 | 欧美成人免费av一区二区三区 | 精品一区二区三区视频在线观看免费 | 欧美日韩黄片免| 亚洲 欧美一区二区三区| 大码成人一级视频| 午夜福利免费观看在线| 国产精品国产av在线观看| 日韩 欧美 亚洲 中文字幕| 高潮久久久久久久久久久不卡| 精品国产一区二区三区久久久樱花| 在线观看免费视频网站a站| av线在线观看网站| 国产精品免费视频内射| 高清欧美精品videossex| 日日爽夜夜爽网站| 日日爽夜夜爽网站| 我的亚洲天堂| 女人精品久久久久毛片| 免费看十八禁软件| 午夜福利乱码中文字幕| 国产一区二区三区综合在线观看| 极品教师在线免费播放| 51午夜福利影视在线观看| 精品视频人人做人人爽| 人人妻,人人澡人人爽秒播| 人人妻,人人澡人人爽秒播| 999久久久国产精品视频| 日韩免费av在线播放| 女同久久另类99精品国产91| 老司机亚洲免费影院| 国产成人免费无遮挡视频| 日日摸夜夜添夜夜添小说| 欧美最黄视频在线播放免费 | 老司机午夜十八禁免费视频| 在线观看日韩欧美| videos熟女内射| 涩涩av久久男人的天堂| 伊人久久大香线蕉亚洲五| 国产精品98久久久久久宅男小说| 国产精品亚洲一级av第二区| 国产精品98久久久久久宅男小说| 一本综合久久免费| 国产成+人综合+亚洲专区| 女人爽到高潮嗷嗷叫在线视频| 91九色精品人成在线观看| 欧美午夜高清在线| 在线视频色国产色| 国产精品久久久久久精品古装| 纯流量卡能插随身wifi吗| 国产一区二区三区视频了| 一进一出抽搐gif免费好疼 | 久久精品亚洲av国产电影网| 国产精品av久久久久免费| 五月开心婷婷网| 欧美黄色片欧美黄色片| 国产精品九九99| 夜夜爽天天搞| 国产成人精品在线电影| 99国产精品免费福利视频| 久久国产精品男人的天堂亚洲| 精品人妻在线不人妻| 亚洲片人在线观看| 国产精品国产av在线观看| 午夜福利一区二区在线看| 亚洲精品久久成人aⅴ小说| 99国产极品粉嫩在线观看| 精品亚洲成a人片在线观看| tube8黄色片| 亚洲精品一卡2卡三卡4卡5卡| 高清欧美精品videossex| 欧美乱码精品一区二区三区| 色播在线永久视频| 在线十欧美十亚洲十日本专区| 亚洲精华国产精华精| 中文字幕人妻丝袜一区二区| 老汉色av国产亚洲站长工具| 十八禁网站免费在线| 成人亚洲精品一区在线观看| 久久精品国产清高在天天线| 日本欧美视频一区| 黄色怎么调成土黄色| 淫妇啪啪啪对白视频| 久久热在线av| 看片在线看免费视频| 热99国产精品久久久久久7| 美女国产高潮福利片在线看| 不卡一级毛片| 国产精品久久久久久精品古装| 婷婷精品国产亚洲av在线 | 精品国产乱码久久久久久男人| 91av网站免费观看| 中文字幕人妻丝袜一区二区| 精品国产一区二区三区久久久樱花| 免费观看人在逋| 人成视频在线观看免费观看| 中亚洲国语对白在线视频| 免费女性裸体啪啪无遮挡网站| 国产成人精品久久二区二区91| 亚洲熟妇熟女久久| 国产免费现黄频在线看| 这个男人来自地球电影免费观看| 亚洲一区二区三区欧美精品| 电影成人av| 成人亚洲精品一区在线观看| 久久99一区二区三区| 国产99白浆流出| 99国产极品粉嫩在线观看| 国产xxxxx性猛交| 精品亚洲成国产av| 免费在线观看完整版高清| 久久国产精品大桥未久av| a级毛片在线看网站| 美女 人体艺术 gogo| 国产亚洲欧美98| 黑人欧美特级aaaaaa片| 欧美人与性动交α欧美软件| 亚洲欧美激情在线| 新久久久久国产一级毛片| 亚洲av成人不卡在线观看播放网| 国产野战对白在线观看| 18禁国产床啪视频网站| 人人妻人人澡人人爽人人夜夜| 免费高清在线观看日韩| 国产成人精品在线电影| 亚洲av欧美aⅴ国产| 国产色视频综合| 性色av乱码一区二区三区2| 亚洲五月天丁香| 99热网站在线观看| 99热国产这里只有精品6| 三级毛片av免费| e午夜精品久久久久久久| 性少妇av在线| 色尼玛亚洲综合影院| 欧美大码av| 国产激情欧美一区二区| 麻豆国产av国片精品| 人人妻人人澡人人爽人人夜夜| 国产成人精品在线电影| 久久国产精品男人的天堂亚洲| 国产精品98久久久久久宅男小说| 亚洲精品成人av观看孕妇| 国产精品免费视频内射| 国产亚洲精品久久久久久毛片 | 亚洲 国产 在线| 亚洲色图av天堂| 亚洲欧美一区二区三区久久| 欧美 亚洲 国产 日韩一| www.自偷自拍.com| 精品一品国产午夜福利视频| 精品人妻熟女毛片av久久网站| 欧美日韩成人在线一区二区| 18在线观看网站| 如日韩欧美国产精品一区二区三区| 这个男人来自地球电影免费观看| 性少妇av在线| www.熟女人妻精品国产| 看黄色毛片网站| 日韩免费av在线播放| 国产野战对白在线观看| 丁香欧美五月| 国产精品.久久久| 精品人妻熟女毛片av久久网站| 亚洲熟妇熟女久久| 久久精品国产清高在天天线| 香蕉丝袜av| 自拍欧美九色日韩亚洲蝌蚪91| 黄片播放在线免费| 国产精品电影一区二区三区 | 99国产极品粉嫩在线观看| 免费看a级黄色片| 成年人免费黄色播放视频| 在线观看日韩欧美| 在线观看日韩欧美| 亚洲成a人片在线一区二区| svipshipincom国产片| 制服人妻中文乱码| 侵犯人妻中文字幕一二三四区| 18禁观看日本| 丁香欧美五月| 亚洲久久久国产精品| 免费不卡黄色视频| 男女高潮啪啪啪动态图| 国产日韩一区二区三区精品不卡| 亚洲欧美精品综合一区二区三区| 国产精品综合久久久久久久免费 | 热re99久久国产66热| 欧美日韩国产mv在线观看视频| 国产精品98久久久久久宅男小说| 涩涩av久久男人的天堂| 大香蕉久久网| 亚洲黑人精品在线| 中文字幕人妻熟女乱码| 窝窝影院91人妻| a级毛片黄视频| 丝瓜视频免费看黄片| 黄色视频不卡| 国产精华一区二区三区| 久久久久国产一级毛片高清牌| 亚洲欧美一区二区三区久久| 国产亚洲精品久久久久久毛片 | 久久久久久久精品吃奶| 女警被强在线播放| 黑人欧美特级aaaaaa片| netflix在线观看网站| 成人国产一区最新在线观看| 欧美在线黄色| 久久人妻福利社区极品人妻图片| 视频在线观看一区二区三区| 一级作爱视频免费观看| 免费一级毛片在线播放高清视频 | 一区福利在线观看| 岛国毛片在线播放| 桃红色精品国产亚洲av| 久久精品国产99精品国产亚洲性色 | 99re在线观看精品视频| 捣出白浆h1v1| 在线观看午夜福利视频| 成人18禁在线播放| 热re99久久国产66热| 亚洲国产看品久久| 一本一本久久a久久精品综合妖精| 亚洲七黄色美女视频| 亚洲人成伊人成综合网2020| 国产成人精品久久二区二区免费| 久久久精品区二区三区| 国产主播在线观看一区二区| 日韩欧美三级三区| 日本一区二区免费在线视频| 国产日韩欧美亚洲二区| 精品亚洲成a人片在线观看| 99re6热这里在线精品视频| 高清毛片免费观看视频网站 | 亚洲欧美色中文字幕在线| 国产欧美亚洲国产| 欧美乱妇无乱码| 丰满人妻熟妇乱又伦精品不卡| 免费在线观看影片大全网站| 天堂动漫精品| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品 欧美亚洲| 看免费av毛片| 国产无遮挡羞羞视频在线观看| 手机成人av网站| 亚洲伊人色综图| 精品国产一区二区三区久久久樱花| 多毛熟女@视频| 成年版毛片免费区| 国产精品香港三级国产av潘金莲| 在线看a的网站| 欧美一级毛片孕妇| 精品国产国语对白av| 欧美黑人精品巨大| 精品免费久久久久久久清纯 | 91九色精品人成在线观看| 12—13女人毛片做爰片一| 亚洲一码二码三码区别大吗| 国产精品 欧美亚洲| 欧美+亚洲+日韩+国产| 久久亚洲真实| svipshipincom国产片| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利欧美成人| 黄色片一级片一级黄色片| 一级黄色大片毛片| 日韩熟女老妇一区二区性免费视频| 一区在线观看完整版| 精品国产乱子伦一区二区三区| 精品人妻1区二区| 99精国产麻豆久久婷婷| 在线av久久热| 999精品在线视频| 国产免费现黄频在线看| 久久久国产一区二区| 777久久人妻少妇嫩草av网站| 国产免费av片在线观看野外av| av网站免费在线观看视频| 亚洲综合色网址| a级毛片黄视频| 亚洲成av片中文字幕在线观看| 亚洲av欧美aⅴ国产| 黄色丝袜av网址大全| 日本欧美视频一区| 大香蕉久久成人网| 夜夜夜夜夜久久久久| 精品一区二区三区视频在线观看免费 | 丰满饥渴人妻一区二区三| 他把我摸到了高潮在线观看| 日韩一卡2卡3卡4卡2021年| 女人爽到高潮嗷嗷叫在线视频| 搡老岳熟女国产| 黄色片一级片一级黄色片| 亚洲国产欧美网| 国产亚洲精品第一综合不卡| 又大又爽又粗| 丝袜美足系列| 精品国产美女av久久久久小说| 嫁个100分男人电影在线观看| 免费观看精品视频网站| 久久国产精品男人的天堂亚洲| 天天操日日干夜夜撸| 黄色a级毛片大全视频| 搡老岳熟女国产| 亚洲七黄色美女视频| 丰满人妻熟妇乱又伦精品不卡| 黑人巨大精品欧美一区二区蜜桃| 黄频高清免费视频| 在线播放国产精品三级| 天天躁夜夜躁狠狠躁躁| 精品国产乱码久久久久久男人| 男人操女人黄网站| av不卡在线播放| 午夜福利一区二区在线看| 国产精品成人在线| 亚洲一区中文字幕在线| 国产不卡av网站在线观看| 欧美乱妇无乱码| av免费在线观看网站| 十八禁网站免费在线| 在线播放国产精品三级| 飞空精品影院首页| 欧美激情高清一区二区三区| 我的亚洲天堂| tube8黄色片| 国产有黄有色有爽视频| 黄色视频不卡| 精品人妻在线不人妻| 日韩欧美国产一区二区入口| 成人精品一区二区免费| 亚洲成国产人片在线观看| 久久99一区二区三区| 免费看十八禁软件| 久久久久久亚洲精品国产蜜桃av| 欧美性长视频在线观看| 18禁国产床啪视频网站| 国产精品亚洲av一区麻豆| 国产精品国产高清国产av | 国产又色又爽无遮挡免费看| 亚洲欧美激情综合另类| 欧美日韩瑟瑟在线播放| 精品国产一区二区三区四区第35| 日本撒尿小便嘘嘘汇集6| 国产精品99久久99久久久不卡| 亚洲精品国产精品久久久不卡| 两个人看的免费小视频| 国产成人欧美| 亚洲欧美一区二区三区黑人| netflix在线观看网站| 亚洲国产看品久久| 最近最新免费中文字幕在线| 一边摸一边抽搐一进一出视频| 99精品欧美一区二区三区四区| 免费av中文字幕在线| 久久久久久久国产电影| 人成视频在线观看免费观看| 操美女的视频在线观看| 精品久久久久久久久久免费视频 | av视频免费观看在线观看| 久久久久国内视频| 手机成人av网站| 久久久水蜜桃国产精品网| 精品第一国产精品| 99久久人妻综合| 一边摸一边抽搐一进一小说 | 自线自在国产av| 俄罗斯特黄特色一大片| 亚洲av片天天在线观看| 国产真人三级小视频在线观看| 亚洲精品自拍成人| 欧美日韩瑟瑟在线播放| 99国产极品粉嫩在线观看| 精品免费久久久久久久清纯 | 中文亚洲av片在线观看爽 | 久热这里只有精品99| 国产成人啪精品午夜网站| 中文字幕制服av| 十八禁人妻一区二区| 亚洲一区二区三区不卡视频| 国产精品香港三级国产av潘金莲| 人妻久久中文字幕网| 丝瓜视频免费看黄片| 无人区码免费观看不卡| 老司机在亚洲福利影院| 黄网站色视频无遮挡免费观看| 国产有黄有色有爽视频| 高清av免费在线| 亚洲精品成人av观看孕妇| 人妻一区二区av| av片东京热男人的天堂| 身体一侧抽搐| 亚洲人成电影观看| 男男h啪啪无遮挡| 三上悠亚av全集在线观看| 久久久久久久精品吃奶| 亚洲avbb在线观看| 丝袜在线中文字幕| 大型av网站在线播放| 久久国产亚洲av麻豆专区| 亚洲国产看品久久| 国产精品久久电影中文字幕 | 一级,二级,三级黄色视频| 亚洲第一欧美日韩一区二区三区| 亚洲熟女毛片儿|