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

    Golay3型光學(xué)稀疏孔徑系統(tǒng)退化圖像的頻率信息提取及合成研究*1

    2016-08-01 11:44:45姜愛民戴妍峰
    天文研究與技術(shù) 2016年3期

    白 靜,姜愛民,戴妍峰

    (1. 中國科學(xué)院國家天文臺,北京 100012;2. 中國科學(xué)院空間天文與技術(shù)重點(diǎn)實(shí)驗(yàn)室,北京 100012;3. 中國科學(xué)院大學(xué),北京 100049)

    ?

    Golay3型光學(xué)稀疏孔徑系統(tǒng)退化圖像的頻率信息提取及合成研究*1

    白靜1,2,3,姜愛民1,2,戴妍峰1,2

    (1. 中國科學(xué)院國家天文臺,北京100012;2. 中國科學(xué)院空間天文與技術(shù)重點(diǎn)實(shí)驗(yàn)室,北京100012;3. 中國科學(xué)院大學(xué),北京100049)

    摘要:利用組成星座的小衛(wèi)星,分別攜帶分離的子望遠(yuǎn)鏡和合光成像望遠(yuǎn)鏡,構(gòu)成Fizeau型光學(xué)綜合孔徑干涉系統(tǒng),實(shí)現(xiàn)高分辨率的面源目標(biāo)成像是當(dāng)前的研究熱點(diǎn)之一。這種光干涉成像系統(tǒng),由于稀疏度較大,UV覆蓋不全,即空間頻率采樣不連續(xù),表現(xiàn)為系統(tǒng)光學(xué)傳遞函數(shù)有零值存在。要克服UV覆蓋不全的影響,獲得等效的大孔徑望遠(yuǎn)鏡成像效果,需要改變子孔徑的空間排布,獲得不同基線條件下的圖像,進(jìn)行空間頻率信息的提取和合成,最后采用逆濾波的方法達(dá)到提高圖像質(zhì)量的目的。在分析單子孔徑傳遞函數(shù)與系統(tǒng)傳遞函數(shù)關(guān)系的基礎(chǔ)上,優(yōu)化子孔徑的排布方式,采用不同的頻域?yàn)V波器將不同基線獲得的圖像中信噪比高的頻率區(qū)域提取出來進(jìn)行合成,再變換到空域并進(jìn)行逆濾波處理,得到改善的合成圖像。仿真結(jié)果顯示,當(dāng)?shù)玫降母缮鎴D信噪比較低時,該方法可以有效地提高合成圖像的質(zhì)量。

    關(guān)鍵詞:稀疏;孔徑圖像;復(fù)原頻域?yàn)V波;圖像合成

    由于孔徑的衍射效應(yīng),傳統(tǒng)望遠(yuǎn)鏡的成像分辨率受限于入瞳口徑。受加工工藝和成本等因素的制約,光學(xué)系統(tǒng)的口徑不可能無限增大[1]。為了解決這個問題,研究人員提出了光干涉的方法。隨著技術(shù)的發(fā)展,F(xiàn)izeau型光干涉技術(shù)因具有一定的成像視場和對面源短時成像的特點(diǎn),受到越來越多的重視,成為當(dāng)前的研究熱點(diǎn)之一。早在1996年,美國空軍武器實(shí)驗(yàn)室就啟動了用于建立對地觀測的可展開式望遠(yuǎn)鏡陣列UltraLITE[2]項(xiàng)目。

    結(jié)合近年來蓬勃發(fā)展的小衛(wèi)星星座技術(shù),人們提出了基于小衛(wèi)星的Fizeau型光干涉成像系統(tǒng)?;痉椒ㄊ牵涸诿總€小衛(wèi)星上各自攜帶一臺小口徑望遠(yuǎn)鏡,在中心的主衛(wèi)星上攜帶合光成像望遠(yuǎn)鏡,子衛(wèi)星和主衛(wèi)星在空間上按一定的位置排布,來自各子望遠(yuǎn)鏡的光束在合光成像望遠(yuǎn)鏡的焦面上相干疊加,干涉成像。這種形式的光干涉成像系統(tǒng),受小衛(wèi)星避撞設(shè)計(jì)的需要,以及小衛(wèi)星有效載荷質(zhì)量和體積的限制,子望遠(yuǎn)鏡的口徑和要實(shí)現(xiàn)的衍射極限分辨率所對應(yīng)的口徑相比,必然差別很大,屬于稀疏度很大的光學(xué)綜合孔徑系統(tǒng)。此時,由于UV覆蓋不全,系統(tǒng)不可避免地對相應(yīng)的空間頻率在采樣時出現(xiàn)缺失。為了克服UV覆蓋不全的影響,對面源目標(biāo)的空間頻率信息進(jìn)行更好的采樣,可以通過改變子衛(wèi)星和主衛(wèi)星的相對位置來改變系統(tǒng)采樣頻率的覆蓋范圍,獲得之前缺失的空間頻率采樣信息,即獲得不同基線布局下的干涉圖。再經(jīng)過圖像合成和逆濾波處理,最終得到包含目標(biāo)較完整信息的圖像。通過這樣的方法,使合成的圖像分辨率等效為一個更大口徑光學(xué)系統(tǒng)的分辨率,達(dá)到高分辨率成像的目的。

    通常的合成方法是將這些在不同基線下得到的圖像在對齊的基礎(chǔ)上[3]進(jìn)行疊加再解卷積,但是由于采集的圖像受子望遠(yuǎn)鏡口徑較小、積分時間受限等各種因素的影響,信噪比通常都不高。仿真實(shí)驗(yàn)發(fā)現(xiàn),當(dāng)?shù)玫降母缮鎴D在信噪比較低的情況下,采用直接疊加再解卷積的處理方法效果不理想。為了提高圖像質(zhì)量,本文提出將不同基線下獲得的圖像通過相應(yīng)的頻域?yàn)V波器,將高信噪比的頻域信息提取出來再合成,以減少合成圖像中包含的噪聲,從而獲得更好的合成圖像。

    1原理概述

    成像過程可以看作是目標(biāo)物體頻率函數(shù)F(fx,fy)與光學(xué)傳遞函數(shù)OTF(fx,fy)乘積后輸出像頻譜函數(shù)G(fx,fy)。在非相干光照射下,衍射受限系統(tǒng)可認(rèn)為是只改變了輸入的各頻率光強(qiáng)分量而不改變它們的位相,光學(xué)傳遞函數(shù)在數(shù)值上等于調(diào)制傳遞函數(shù)MTF:

    (1)

    其中,?表示自相關(guān)符號;P(λffx,λffy)是孔徑函數(shù)。

    稀疏孔徑成像系統(tǒng)在理想情況下可以看作是由N個直徑為d的圓形孔徑構(gòu)成,在數(shù)學(xué)上可以表示為一個子孔徑光瞳函數(shù)與δ函數(shù)陣列的卷積,即

    (2)

    其中,P(x,y)是光瞳函數(shù);N是子孔徑個數(shù);circ()是圓域孔徑函數(shù);(xn,yn)是第n個子孔徑光瞳的中心坐標(biāo)。

    Golay3是結(jié)構(gòu)最簡單的二維陣列結(jié)構(gòu),由3個圓形子孔徑組成,而且這3個圓心依次位于一個等邊三角形的3個頂點(diǎn)上[4]。

    從(2)式可以推導(dǎo)出Golay3型稀疏孔徑的歸一化調(diào)制傳遞函數(shù)為

    (3)

    其中,(xi-xj,yi-yj)表示子孔徑的相對位置。

    根據(jù)(3)式,可以將稀疏孔徑系統(tǒng)的調(diào)制傳遞函數(shù)看作是由分布在不同位置上,強(qiáng)度為1/3倍(零頻除外)的單個子孔徑的調(diào)制傳遞函數(shù)組合而成。除零頻外,其它調(diào)制傳遞函數(shù)分布的位置都與子孔徑排布有關(guān),也就是說如果已知調(diào)制傳遞函數(shù)的位置分布就可以推算出子孔徑的位置排布。當(dāng)系統(tǒng)調(diào)制傳遞函數(shù)對空間采樣頻率覆蓋的一些區(qū)域出現(xiàn)零值或嚴(yán)重衰減時,就可以找到相應(yīng)的空間位置增加孔徑排布進(jìn)行優(yōu)化,最后將這些在不同排布(即不同基線)下得到的退化圖進(jìn)行合成處理。

    稀疏孔徑系統(tǒng)主要的噪聲有光子噪聲及CCD噪聲等,總噪聲分布近似為高斯分布,一般認(rèn)為是零均值的高斯白噪聲[5],即功率譜密度在整個頻域內(nèi)均勻分布的噪聲。已知單幅退化圖的噪聲功率為ni,如果將退化圖直接疊加就相當(dāng)于直接疊加了∑ini倍的噪聲。

    圖1某子孔徑排布的調(diào)制傳遞函數(shù)示意圖

    Fig.1The MTF of a sub-aperture configuration

    因此,本文針對各種基線排布定義了相應(yīng)的頻域?yàn)V波器Hi,將每幅圖中頻域信噪比較高部分的信號提取出來再合成。由于每幅干涉圖的濾波器Hi加權(quán)為1,而且范圍不完全重復(fù),這樣得到的最終合成圖的噪聲量∑iHi·ni明顯小于通過直接疊加得到的,信噪比有了顯著的提升。

    2圖像頻率提取及合成方法

    設(shè)單個子孔徑直徑為1 m。從系統(tǒng)子孔徑相切開始,以子孔徑直徑為間距向外挪動子孔徑兩次共得到3種不同的排布。此時外接最大包圍圓直徑為6.15 m。根據(jù)這3種排布得到了3種不同分布的調(diào)制傳遞函數(shù)。由于發(fā)現(xiàn)空間頻率信息分別在排布1與排布2、排布2與排布3截止頻率交接處附近明顯衰減,所以增加新的子孔徑排布4和5,對這兩交界附近的調(diào)制傳遞函數(shù)進(jìn)行補(bǔ)償和增強(qiáng)。優(yōu)化后,這5種排布的調(diào)制傳遞函數(shù)在最大截止頻率方向的部分截面示意圖(在同一坐標(biāo)系) 如圖2。

    與稀疏孔徑系統(tǒng)成像質(zhì)量相當(dāng)?shù)膯慰讖街睆椒Q作等效口徑。由于稀疏孔徑系統(tǒng)的調(diào)制傳遞函數(shù)與單孔徑系統(tǒng)相比,分布復(fù)雜,形狀不規(guī)則,目前尚沒有計(jì)算等效口徑的統(tǒng)一方法[6]。從圖2可以看出,當(dāng)前系統(tǒng)最大截止頻率5.4/λf(m)明顯小于系統(tǒng)包圍圓所對應(yīng)的6.15/λf(m)。所以選最大截止頻率對應(yīng)的口徑作為系統(tǒng)等效口徑,此處為5.4 m。

    在這5種不同的子孔徑排布下,獲取5幅不同的圖像(記作g1、g2、g3、g4和g5)。將這5幅圖像通過傅里葉變換到頻域(記作G1、G2、G3、G4和G5)。根據(jù)各自調(diào)制傳遞函數(shù)的分布,分別定義各自的頻域?yàn)V波器,如(4)~(8)式。其中,H1為低通濾波器;H2、H3、H4和H5分別是不同通帶范圍的帶通濾波器。

    (4)

    圖2優(yōu)化后各排布的調(diào)制傳遞函數(shù)部分截面圖

    Fig.2The optimized sectional view of partial MTFfor every sub-aperture configuration

    利用這些頻域?yàn)V波器可以將5幅圖中信噪比較好的頻率信息提出來。提取及合成方法如圖3。圖像頻譜合成方法如(9)式:

    G合成=G1H1+G2H2+G3H3+G4H4+G5H5.

    (9)

    頻譜直接疊加如(10)式:

    (10)

    圖3空間頻率信息提取及合成. (a)原始圖像頻率信息;(b)從孔徑排布1提取的頻率信息;(c)、 (d)、 (e)和(f)分別對應(yīng)從孔徑排布2、3、4和5提取的頻率信息;(g)合成圖像頻率信息

    Fig.3Frequency extraction and synthesis. (a) The frequency of the original objective image; (b) The partial frequency extracted from the degraded image of configuration 1; (c) (d) (e) and (f) The corresponding partial frequency extracted from the degraded image of new configuration 2, 3, 4 and 5; (g) The final synthesized degraded image

    在仿真實(shí)驗(yàn)中引入加性高斯白噪聲,當(dāng)信噪比為35 dB和25 dB時分別得到退化圖如圖4。

    圖4不同方法得到的仿真圖像. (a)等效口徑圖像;(b)直接疊加得到的合成圖像;(c)頻域提取方法得到的合成圖像

    Fig.4Degraded images obtained by different methods. (a) The degraded image obtained by effective Diameter method; (b) The synthesized degraded image obtained by overlapping frequency directly; (c) The synthesized degraded image obtained by frequency filtrations

    從圖4可以看出,當(dāng)圖像信噪比下降后,直接疊加的合成圖像像質(zhì)逐漸變差,此時采用頻域提取的合成圖像像質(zhì)占優(yōu)。

    3圖像復(fù)原(解卷積)及評價

    由不同基線圖像合成的圖像需經(jīng)過解卷積才能更好地恢復(fù)圖像。逆濾波是解卷積最直接的方法,設(shè)G(u,v)為退化圖像的傅里葉譜,H(u,v)為點(diǎn)擴(kuò)散函數(shù)的傅里葉譜,則估計(jì)的原始圖像的傅里葉譜F^(u,v)為

    (11)

    維納濾波(最小均方誤差濾波)是對逆濾波方法的改進(jìn)。維納濾波方法綜合考慮了退化函數(shù)和噪聲的統(tǒng)計(jì)特性,以減少噪聲的影響,如(12)式:

    (12)

    其中,*表示共軛;Sn(u,v)/Sf(u,v)是噪聲功率譜和原始信號功率譜的比值。當(dāng)噪聲的統(tǒng)計(jì)特性未知時,常用一個常數(shù)K近似。若γ=1,W(u,v)是維納濾波器;若γ=0,則系統(tǒng)變成了單純的去卷積濾波器。

    由于圖像經(jīng)過頻率提取和合成后,合成圖的噪信比有較明顯的起伏規(guī)律。因此,在去解卷積濾波時,對維納濾波器也進(jìn)行了相應(yīng)的改進(jìn):在對應(yīng)不同的區(qū)域,根據(jù)合成圖噪信比的不同采用了不同的K值。去卷積后的復(fù)原圖如圖5。

    圖5不同方法得到的復(fù)原效果. (a)等效滿口徑復(fù)原圖;(b)直接疊加得到的復(fù)原圖;(c)頻域提取方法得到的合成復(fù)原圖
    Fig.5Reconstructional images obtained by different methods. (a) The reconstructional image obtained by effective Diameter method; (b) The reconstructional synthetic image obtained by overlapping frequency directly; (c) The reconstructional synthetic image obtained by frequency filtrations

    對于圖像的復(fù)原效果,除了主觀判斷外,采用傳統(tǒng)的均方誤差(Mean Squared Error, MSE)法和當(dāng)前國際上比較流行的SSIM(Structural Similarity)法。均方差法基于誤差理論,將圖像看作是一個個孤立的點(diǎn),通過計(jì)算每個點(diǎn)的灰度值差異而對圖像進(jìn)行評測,見(13)式。SSIM法是對圖像在亮度、對比度、結(jié)構(gòu)3個層面上分別進(jìn)行量化,通過加權(quán)相乘獲得結(jié)構(gòu)相似值,綜合考慮了人類視覺對圖像結(jié)構(gòu)性信息的提取功能,比均方差法更符合人類的視覺特性[7],見(14)式:

    (13)

    (14)

    其中,f(x,y)表示參考圖像;f^(x,y)表示待評價圖像;α、β、γ為加權(quán)系數(shù);l(f,f^)、c(f,f^)、s(f,f^)為相應(yīng)的亮度、對比度、結(jié)構(gòu)比較函數(shù)。

    均方根值越小,表示待評價圖像和參考圖像越相似;結(jié)構(gòu)相似度均值越大,則表示結(jié)構(gòu)相似值越高,說明待評價圖像和參考圖像的相似度越高。

    在不同信噪比下得到的退化圖經(jīng)過解卷積復(fù)原后,分別計(jì)算均方根和結(jié)構(gòu)相似度均值,見表1。

    結(jié)合圖5和表1可以看出,當(dāng)信噪比為35 dB時,對合成圖(c)的客觀評價指標(biāo)比對直接疊加圖像(b)的客觀評價指標(biāo)只是稍微好一點(diǎn),從主觀上觀察,人的眼睛基本分不出哪幅圖更好。

    然而當(dāng)信噪比下降到25 dB時,對合成圖(c)客觀評價指標(biāo)就明顯比對直接疊加圖像(b)的客觀評價值要好一些。在用人眼進(jìn)行主觀評價時,能明顯看出合成圖(c)中的圖像細(xì)節(jié)更清晰。

    表1各退化圖客觀評價計(jì)算結(jié)果

    Table 1Objective evaluation results foreach degraded image

    信噪比35圖(a)圖(b)圖(c)MSSIM0.97680.79190.8003MSE0.00060.00330.0031信噪比25圖(a)圖(b)圖(c)MSSIM0.91660.66300.7441MSE0.00190.00350.0030

    圖(a)是以原圖為評價參考, 圖(b)和圖(c)都是以圖(a)為評價參考
    Fig (a) is with reference to the original image; (b) and (c) are with reference to Fig (a)

    4結(jié)論

    當(dāng)原始圖像的信噪比受限時,利用頻域抽取算法將在不同子孔徑排布下得到的干涉圖合成為一幅包含較完整目標(biāo)信息的退化圖像,再經(jīng)過改進(jìn)的維納濾波算法得到了最終較清晰的復(fù)原圖像。由于保留了各幅干涉圖頻域中信噪比較高的部分,在解卷積復(fù)原后,其圖像質(zhì)量比直接疊加得到的有明顯提升,采用人眼判讀和客觀評價的方法都顯示,該方法能更好地恢復(fù)原圖像的結(jié)構(gòu)信息。

    致謝:本文得到王景宇研究員、薛建偉博士、宋衛(wèi)強(qiáng)工程師的幫助,在此一并致謝。

    參考文獻(xiàn):

    [1]van Belle G T, Meinel A B, Meinel M P. The scaling relationship between telescope cost and aperture size for very large telescopes[C]// Proceedings of SPIE. 2004: 563-570.

    [2]Powers M, Leitner J, Hackney E, et al. Assessment of a large aperture telescope trade space and active opto-mechanical control architecture[C]// Proceedings of IEEE Aerospace Conference.1997: 197-229.

    [3]高莉莉, 劉忠, 金振宇. 綜合孔徑干涉望遠(yuǎn)鏡的高分辨率圖像重建[J]. 天文研究與技術(shù)——國家天文臺臺刊, 2009, 6(4): 327-333.

    Gao Lili, Liu Zhong, Jin Zhenyu. High-resolution reconstruction of images from an interferometry telescope[J]. Astronomical Research & Technology——Publications of National Astronomical Observatories of China, 2009, 6(4): 327-333.

    [4]Golay M J E. Point array having compact, nonredundant autocorrelations[J]. Journal of the Optical Society of America, 1971, 61(2): 272-273.

    [5]劉麗, 江月松. 綜合孔徑成像原理與應(yīng)用[M]. 北京: 國防工業(yè)出版社, 2013.

    [6]Miller N, Duncan B, Dierking M P. Resolution enhanced sparse aperture imaging[C]// Proceedings of IEEE Aerospace Conference. 2006.

    [7]Wang Z, Bovik A C, Sheikh H R, et al. Image quality assessment: from error visibility to structure similarity[J]. IEEE Transactions on Image Processing, 2004, 13(4): 600-612.

    *基金項(xiàng)目:國家高技術(shù)研究發(fā)展計(jì)劃 (863計(jì)劃) (2013AA122603) 資助.

    收稿日期:2015-11-27;

    修訂日期:2016-01-06

    作者簡介:白靜,女,碩士. 研究方向:光學(xué)綜合孔徑圖像處理. Email: jbai@nao.cas.cn

    中圖分類號:TN911.73

    文獻(xiàn)標(biāo)識碼:A

    文章編號:1672-7673(2016)03-0351-07

    Frequency Extraction and Degraded Image Synthesis for Golay3 Type Optical Sparse-Aperture Systems

    Bai Jing1,2,3, Jiang Aimin1,2, Dai Yanfeng1,2

    (1. National Astronomical Observatories, Chinese Academy of Science, Beijing 100012, China; 2. Key Laboratory of Space Astronomy and Technology, National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China; 3. University of Chinese Academy of Sciences, Beijing 100049, China, Email: jbai@nao.cas.cn)

    Abstract:A Fizeau-typesparse-aperture imaging system made of a constellation of telescope-bearing satellites can potentially achieve extremely high resolution. But in such a system, the Modulation Transfer Function (MTF) usually contains a zero value. This means that the system fails to acquire a particular spatial frequency. To overcome this problem, we acquire multiple degraded images by repeatedly changing the sub-aperture configuration and then synthesize their spatial frequencies to create a single improved image. We first analyze the MTFs between each individual sub-aperture and the whole sparse aperture system. Then, we design new sub-aperture configurations that acquire more spatial frequency and use frequency-domain filters to extract the high-SNR part from each degraded image. Finally, we synthesize these frequencies before image reconstruction. Our results demonstrate that this method effectively improves the quality of the final synthesized image by obtaining as much spatial frequency information from all degraded images as possible.

    Key words:Sparse aperture; Image reconstruction; Frequency-domain filtering; Image synthesis

    CN 53-1189/PISSN 1672-7673

    少妇人妻一区二区三区视频| 床上黄色一级片| av在线播放精品| 一本久久精品| 国产综合懂色| 波野结衣二区三区在线| 亚洲欧美成人综合另类久久久 | 国产午夜精品一二区理论片| 日韩人妻高清精品专区| www日本黄色视频网| 黄片wwwwww| 中文精品一卡2卡3卡4更新| 最近中文字幕高清免费大全6| 免费看美女性在线毛片视频| 国产成人aa在线观看| 桃色一区二区三区在线观看| 国产综合懂色| 熟女电影av网| 免费不卡的大黄色大毛片视频在线观看 | 亚洲乱码一区二区免费版| АⅤ资源中文在线天堂| 久久久久九九精品影院| 国产成人a区在线观看| 三级男女做爰猛烈吃奶摸视频| 免费看a级黄色片| 国产亚洲午夜精品一区二区久久 | 欧美zozozo另类| 成年免费大片在线观看| 欧美+日韩+精品| 波多野结衣高清无吗| 亚洲av成人av| 人妻少妇偷人精品九色| 啦啦啦观看免费观看视频高清| av.在线天堂| 亚洲怡红院男人天堂| av在线播放精品| 亚洲精品456在线播放app| 久久亚洲国产成人精品v| 久久热精品热| 成人漫画全彩无遮挡| 中文字幕精品亚洲无线码一区| 啦啦啦韩国在线观看视频| 国产一级毛片七仙女欲春2| 国产成人91sexporn| 久久精品久久精品一区二区三区| 欧美另类亚洲清纯唯美| 久久久久久久久久久免费av| 免费播放大片免费观看视频在线观看 | 日韩 亚洲 欧美在线| 国产成人a区在线观看| 亚洲电影在线观看av| 我的女老师完整版在线观看| 插阴视频在线观看视频| 欧美一级a爱片免费观看看| 久久精品91蜜桃| 高清日韩中文字幕在线| 午夜爱爱视频在线播放| kizo精华| 婷婷色综合大香蕉| videossex国产| 高清毛片免费看| 少妇熟女aⅴ在线视频| 成年版毛片免费区| 中文资源天堂在线| 国产一区二区在线av高清观看| 搡女人真爽免费视频火全软件| 久久久成人免费电影| 99热6这里只有精品| 三级毛片av免费| 一级毛片久久久久久久久女| 久久久久久伊人网av| 久久婷婷人人爽人人干人人爱| 国产亚洲5aaaaa淫片| 亚洲精品乱码久久久久久按摩| 熟妇人妻久久中文字幕3abv| 久久精品熟女亚洲av麻豆精品 | 少妇裸体淫交视频免费看高清| 国产精品国产三级国产av玫瑰| 日本免费在线观看一区| 国内精品美女久久久久久| 国产极品天堂在线| 久久国产乱子免费精品| 久久久成人免费电影| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 春色校园在线视频观看| videos熟女内射| 人人妻人人澡人人爽人人夜夜 | 久久精品熟女亚洲av麻豆精品 | 亚洲国产最新在线播放| 色哟哟·www| 一区二区三区四区激情视频| 毛片一级片免费看久久久久| av卡一久久| 欧美zozozo另类| 一边亲一边摸免费视频| 99在线人妻在线中文字幕| 国产视频首页在线观看| 青青草视频在线视频观看| 国产精品爽爽va在线观看网站| 超碰av人人做人人爽久久| 91在线精品国自产拍蜜月| 亚洲高清免费不卡视频| 18禁裸乳无遮挡免费网站照片| 一级毛片我不卡| 男人和女人高潮做爰伦理| 欧美一级a爱片免费观看看| 99视频精品全部免费 在线| 久久久久久久久大av| 一本久久精品| 国产精品人妻久久久久久| 日产精品乱码卡一卡2卡三| 乱人视频在线观看| 高清毛片免费看| av福利片在线观看| 有码 亚洲区| 在线免费观看不下载黄p国产| 久久久欧美国产精品| 亚洲乱码一区二区免费版| 亚洲av一区综合| 91精品国产九色| 久久久久国产网址| 国产精品人妻久久久影院| 欧美色视频一区免费| 亚洲av中文字字幕乱码综合| 精华霜和精华液先用哪个| 日日摸夜夜添夜夜爱| 国产久久久一区二区三区| 精品国产露脸久久av麻豆 | 精品久久久久久久久av| 麻豆av噜噜一区二区三区| 伊人久久精品亚洲午夜| 国产一区亚洲一区在线观看| 麻豆精品久久久久久蜜桃| 少妇熟女aⅴ在线视频| 免费播放大片免费观看视频在线观看 | 亚洲国产欧美人成| 国产乱来视频区| 国产精品99久久久久久久久| 美女国产视频在线观看| 99热6这里只有精品| 免费看光身美女| 久久这里只有精品中国| 简卡轻食公司| 亚洲国产精品合色在线| 国产黄a三级三级三级人| 麻豆精品久久久久久蜜桃| 久久鲁丝午夜福利片| 久久这里只有精品中国| 久久精品夜色国产| 亚洲最大成人中文| 国产高清视频在线观看网站| 成人鲁丝片一二三区免费| 成人高潮视频无遮挡免费网站| 亚洲av成人精品一二三区| 最近手机中文字幕大全| 在线播放无遮挡| 少妇被粗大猛烈的视频| 精品欧美国产一区二区三| 国产激情偷乱视频一区二区| 久久精品国产自在天天线| av国产久精品久网站免费入址| 国产一区有黄有色的免费视频 | 国产 一区精品| 一个人观看的视频www高清免费观看| 最近最新中文字幕大全电影3| 成人毛片60女人毛片免费| 国产视频首页在线观看| 老女人水多毛片| 色5月婷婷丁香| 久久久久久久久中文| 免费大片18禁| 亚洲欧洲日产国产| 欧美性猛交╳xxx乱大交人| 最近中文字幕高清免费大全6| 一边亲一边摸免费视频| 日本av手机在线免费观看| 亚洲国产精品成人综合色| 在线免费观看不下载黄p国产| 色尼玛亚洲综合影院| 中国美白少妇内射xxxbb| 国产免费福利视频在线观看| 少妇被粗大猛烈的视频| 男女国产视频网站| 18禁裸乳无遮挡免费网站照片| 亚洲不卡免费看| 身体一侧抽搐| 99热网站在线观看| 中文字幕熟女人妻在线| 欧美日韩精品成人综合77777| 精品人妻偷拍中文字幕| 天天躁夜夜躁狠狠久久av| av在线观看视频网站免费| 日韩三级伦理在线观看| 日日摸夜夜添夜夜爱| 亚洲国产精品久久男人天堂| 亚洲丝袜综合中文字幕| 亚洲中文字幕日韩| av免费在线看不卡| 欧美精品国产亚洲| 美女cb高潮喷水在线观看| 日本黄大片高清| 2021天堂中文幕一二区在线观| 亚洲av成人av| 欧美成人一区二区免费高清观看| 国产精品伦人一区二区| 国产成人午夜福利电影在线观看| 成人毛片60女人毛片免费| 国产精品,欧美在线| 91精品伊人久久大香线蕉| 亚洲一区高清亚洲精品| 精品久久久久久久人妻蜜臀av| 国产午夜福利久久久久久| 久久这里只有精品中国| 亚洲欧美日韩卡通动漫| 国产视频内射| 国产在线男女| 日韩欧美 国产精品| 我的女老师完整版在线观看| 亚洲婷婷狠狠爱综合网| 18禁动态无遮挡网站| 91精品一卡2卡3卡4卡| 欧美一区二区精品小视频在线| 亚洲精品乱久久久久久| 亚洲高清免费不卡视频| 欧美另类亚洲清纯唯美| 熟妇人妻久久中文字幕3abv| 亚洲精品,欧美精品| 国产精品爽爽va在线观看网站| 午夜爱爱视频在线播放| 欧美97在线视频| 九九在线视频观看精品| 最新中文字幕久久久久| 久久精品久久精品一区二区三区| 亚洲一区高清亚洲精品| 91狼人影院| av在线观看视频网站免费| 国产高清国产精品国产三级 | 国产精品嫩草影院av在线观看| 亚洲美女搞黄在线观看| 久久鲁丝午夜福利片| 日韩欧美精品v在线| 搡老妇女老女人老熟妇| 国内精品宾馆在线| 日韩欧美国产在线观看| 国产精品乱码一区二三区的特点| 亚洲精品aⅴ在线观看| 亚洲怡红院男人天堂| 女人十人毛片免费观看3o分钟| 免费看av在线观看网站| 日韩欧美在线乱码| 亚洲精品乱码久久久久久按摩| 欧美一区二区国产精品久久精品| 在线免费十八禁| 精品99又大又爽又粗少妇毛片| 女人十人毛片免费观看3o分钟| 国产在线一区二区三区精 | 亚洲国产成人一精品久久久| 欧美日本视频| 国产亚洲5aaaaa淫片| 日本av手机在线免费观看| 国产精品久久久久久久久免| 赤兔流量卡办理| 日韩人妻高清精品专区| 综合色av麻豆| 麻豆乱淫一区二区| 国产精品不卡视频一区二区| 日本与韩国留学比较| 日本爱情动作片www.在线观看| av黄色大香蕉| 天天躁夜夜躁狠狠久久av| 久久这里只有精品中国| 天堂av国产一区二区熟女人妻| 精品久久久久久久久亚洲| 亚洲最大成人手机在线| 久久这里有精品视频免费| 亚洲精品日韩av片在线观看| av卡一久久| 成人午夜精彩视频在线观看| 久久国产乱子免费精品| 欧美成人一区二区免费高清观看| videos熟女内射| 亚洲中文字幕一区二区三区有码在线看| 天堂√8在线中文| 国产黄色视频一区二区在线观看 | 国产成人91sexporn| 亚洲av免费高清在线观看| 少妇的逼水好多| 五月伊人婷婷丁香| 亚洲国产最新在线播放| 国产伦精品一区二区三区四那| 日韩一本色道免费dvd| 91在线精品国自产拍蜜月| 中文字幕免费在线视频6| 国产av码专区亚洲av| АⅤ资源中文在线天堂| 三级国产精品片| 女人久久www免费人成看片 | 国产精品麻豆人妻色哟哟久久 | 日本色播在线视频| 午夜福利在线观看吧| 国产精品麻豆人妻色哟哟久久 | 久久国产乱子免费精品| 18禁在线无遮挡免费观看视频| 国产一区有黄有色的免费视频 | 日本-黄色视频高清免费观看| 男女视频在线观看网站免费| 三级经典国产精品| 99国产精品一区二区蜜桃av| 国产探花极品一区二区| 免费看美女性在线毛片视频| 在线播放无遮挡| 国产一区二区亚洲精品在线观看| 国产精品野战在线观看| 99在线人妻在线中文字幕| 久久人妻av系列| 亚洲美女视频黄频| 亚洲人成网站在线播| 亚洲国产精品sss在线观看| 人妻系列 视频| 中文字幕av成人在线电影| 亚洲欧美日韩高清专用| 我的老师免费观看完整版| 一区二区三区免费毛片| 国产精品一区二区在线观看99 | 国内精品宾馆在线| 亚洲欧美精品自产自拍| 国产亚洲av嫩草精品影院| h日本视频在线播放| 国产精品.久久久| 精品免费久久久久久久清纯| 麻豆成人av视频| 中文字幕熟女人妻在线| 尾随美女入室| 国产精品精品国产色婷婷| 国产激情偷乱视频一区二区| 色综合站精品国产| 18+在线观看网站| 久久久久久国产a免费观看| 国产免费又黄又爽又色| 变态另类丝袜制服| a级一级毛片免费在线观看| 变态另类丝袜制服| 久久这里有精品视频免费| 长腿黑丝高跟| 欧美潮喷喷水| 九九热线精品视视频播放| 一级毛片aaaaaa免费看小| 精品午夜福利在线看| videos熟女内射| 精品久久久久久电影网 | 亚洲国产欧美在线一区| 亚洲在线自拍视频| 精品午夜福利在线看| av在线老鸭窝| 国产国拍精品亚洲av在线观看| 午夜精品一区二区三区免费看| 在线免费观看的www视频| 亚洲精品一区蜜桃| 亚洲综合色惰| 久久久久久久亚洲中文字幕| 国产av不卡久久| 大又大粗又爽又黄少妇毛片口| 在线免费观看的www视频| 老女人水多毛片| 国产精品国产三级国产av玫瑰| 中文在线观看免费www的网站| 日本wwww免费看| 一卡2卡三卡四卡精品乱码亚洲| 日韩中字成人| 亚洲精品一区蜜桃| 亚州av有码| 男女那种视频在线观看| 国产精品一区二区三区四区久久| 亚洲三级黄色毛片| 欧美性猛交黑人性爽| 久久久久久久久中文| 亚洲最大成人av| 日韩在线高清观看一区二区三区| 少妇裸体淫交视频免费看高清| 亚洲,欧美,日韩| 美女xxoo啪啪120秒动态图| 亚洲性久久影院| 日本免费a在线| 麻豆成人av视频| 日韩成人av中文字幕在线观看| 伊人久久精品亚洲午夜| 男女视频在线观看网站免费| 中文亚洲av片在线观看爽| 99久久精品国产国产毛片| 国产精品久久久久久久电影| 欧美色视频一区免费| ponron亚洲| 亚洲人与动物交配视频| 汤姆久久久久久久影院中文字幕 | 日韩欧美国产在线观看| 欧美精品国产亚洲| av黄色大香蕉| eeuss影院久久| 日本wwww免费看| 亚洲av一区综合| 亚洲人成网站高清观看| 精华霜和精华液先用哪个| 最后的刺客免费高清国语| 欧美日韩国产亚洲二区| 久久久成人免费电影| 少妇人妻精品综合一区二区| h日本视频在线播放| 免费观看人在逋| 十八禁国产超污无遮挡网站| www日本黄色视频网| 中文乱码字字幕精品一区二区三区 | 亚洲欧美日韩高清专用| 三级男女做爰猛烈吃奶摸视频| 色综合站精品国产| 99热这里只有精品一区| 久久欧美精品欧美久久欧美| 日韩高清综合在线| 欧美日韩国产亚洲二区| 亚洲国产欧美人成| 久久久a久久爽久久v久久| 人人妻人人澡欧美一区二区| 国语对白做爰xxxⅹ性视频网站| kizo精华| 国产精品国产三级专区第一集| 国产人妻一区二区三区在| 身体一侧抽搐| 国产亚洲av嫩草精品影院| 中国美白少妇内射xxxbb| 有码 亚洲区| 国产av一区在线观看免费| 男女下面进入的视频免费午夜| 国产精品美女特级片免费视频播放器| 成人综合一区亚洲| 在线观看一区二区三区| 韩国高清视频一区二区三区| 日本黄色视频三级网站网址| 我要搜黄色片| 久久久久久久久中文| 久久人人爽人人爽人人片va| 久久久久久久久久久丰满| 婷婷色综合大香蕉| 亚洲一区高清亚洲精品| 亚洲精品亚洲一区二区| 午夜精品在线福利| 美女脱内裤让男人舔精品视频| 丝袜美腿在线中文| 搞女人的毛片| 亚洲av电影在线观看一区二区三区 | 中文字幕av成人在线电影| 日本欧美国产在线视频| 在线观看66精品国产| 天堂中文最新版在线下载 | 午夜福利在线在线| av卡一久久| 亚洲在久久综合| 欧美一区二区国产精品久久精品| 国产 一区 欧美 日韩| 26uuu在线亚洲综合色| 婷婷色av中文字幕| 国产av在哪里看| 精品人妻熟女av久视频| 三级国产精品欧美在线观看| 禁无遮挡网站| 人体艺术视频欧美日本| 免费搜索国产男女视频| 看免费成人av毛片| 嫩草影院新地址| kizo精华| 日韩强制内射视频| 国产淫语在线视频| 日韩三级伦理在线观看| 精品久久久久久久久久久久久| 亚洲欧洲国产日韩| 成年免费大片在线观看| 亚洲精品成人久久久久久| 高清在线视频一区二区三区 | 亚洲av电影不卡..在线观看| 亚洲自偷自拍三级| 欧美日韩国产亚洲二区| 能在线免费看毛片的网站| 我的老师免费观看完整版| 国产亚洲一区二区精品| 亚洲美女视频黄频| 成人欧美大片| 久久精品影院6| 久久久欧美国产精品| 久久6这里有精品| 精品不卡国产一区二区三区| 热99re8久久精品国产| 国产不卡一卡二| 国产 一区 欧美 日韩| 国国产精品蜜臀av免费| 九草在线视频观看| 欧美色视频一区免费| 国产精品无大码| 日韩欧美国产在线观看| 中文字幕熟女人妻在线| 一级黄色大片毛片| 水蜜桃什么品种好| eeuss影院久久| 午夜精品一区二区三区免费看| 偷拍熟女少妇极品色| 一个人看的www免费观看视频| 精品一区二区三区视频在线| 国产亚洲5aaaaa淫片| 成年av动漫网址| 99在线人妻在线中文字幕| 国产午夜精品论理片| 能在线免费看毛片的网站| 男的添女的下面高潮视频| 久久草成人影院| 最近最新中文字幕免费大全7| 免费av观看视频| 亚洲国产精品成人综合色| 永久免费av网站大全| 国产精品久久久久久久久免| 在线免费观看的www视频| 亚洲国产日韩欧美精品在线观看| 嫩草影院入口| 纵有疾风起免费观看全集完整版 | 亚洲一区高清亚洲精品| 亚洲综合色惰| 久久亚洲国产成人精品v| 一边摸一边抽搐一进一小说| 汤姆久久久久久久影院中文字幕 | 赤兔流量卡办理| 毛片一级片免费看久久久久| 日韩人妻高清精品专区| 亚洲精品乱码久久久v下载方式| 一边亲一边摸免费视频| 亚洲国产精品合色在线| 成人性生交大片免费视频hd| 亚洲精品成人久久久久久| 日韩视频在线欧美| av视频在线观看入口| 国产女主播在线喷水免费视频网站 | 特大巨黑吊av在线直播| 国产一区亚洲一区在线观看| 天天一区二区日本电影三级| 久久久国产成人免费| 免费观看a级毛片全部| 亚洲人成网站高清观看| 久久99热这里只频精品6学生 | av卡一久久| 日韩av不卡免费在线播放| 久久精品国产亚洲av天美| 色综合站精品国产| 色网站视频免费| 少妇人妻精品综合一区二区| 久久久久久久久大av| 中国美白少妇内射xxxbb| 最近手机中文字幕大全| 亚洲国产精品sss在线观看| 国产免费又黄又爽又色| 女人十人毛片免费观看3o分钟| 欧美精品一区二区大全| 亚洲美女搞黄在线观看| 噜噜噜噜噜久久久久久91| 亚洲成色77777| 少妇熟女欧美另类| 亚洲av.av天堂| 在线播放无遮挡| 禁无遮挡网站| 黄片wwwwww| 精品久久久久久久人妻蜜臀av| 久久午夜福利片| 欧美日本亚洲视频在线播放| 国产成人精品一,二区| 成人午夜高清在线视频| 中文天堂在线官网| av线在线观看网站| 国产精品爽爽va在线观看网站| 插逼视频在线观看| 嘟嘟电影网在线观看| 少妇的逼好多水| .国产精品久久| 久久人人爽人人片av| av免费在线看不卡| 精品少妇黑人巨大在线播放 | 国产精品无大码| 国产探花在线观看一区二区| 天天躁日日操中文字幕| 久久欧美精品欧美久久欧美| 又黄又爽又刺激的免费视频.| 日日干狠狠操夜夜爽| 男女啪啪激烈高潮av片| 国产av一区在线观看免费| 亚洲av一区综合| 欧美精品一区二区大全| 99久久成人亚洲精品观看| 综合色丁香网| 两个人的视频大全免费| 免费黄网站久久成人精品| 舔av片在线| av.在线天堂| 久热久热在线精品观看| 久久精品综合一区二区三区| 小说图片视频综合网站| 亚洲一区高清亚洲精品| 麻豆国产97在线/欧美| 欧美最新免费一区二区三区| 成人鲁丝片一二三区免费| 有码 亚洲区| 99视频精品全部免费 在线| 国产v大片淫在线免费观看| 成人特级av手机在线观看| 久久精品91蜜桃| 亚洲av电影不卡..在线观看| 我要看日韩黄色一级片| 两个人的视频大全免费| 99久久成人亚洲精品观看| 嫩草影院新地址| 亚洲av中文av极速乱| 亚洲在线观看片| 色播亚洲综合网|