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

    Compressive imaging based on multi-scale modulation and reconstruction in spatial frequency domain*

    2021-01-21 02:07:36FanLiu劉璠XueFengLiu劉雪峰RuoMingLan藍(lán)若明XuRiYao姚旭日ShenChengDou竇申成XiaoQingWang王小慶andGuangJieZhai翟光杰
    Chinese Physics B 2021年1期
    關(guān)鍵詞:旭日雪峰

    Fan Liu(劉璠), Xue-Feng Liu(劉雪峰),?, Ruo-Ming Lan(藍(lán)若明), Xu-Ri Yao(姚旭日),Shen-Cheng Dou(竇申成), Xiao-Qing Wang(王小慶), and Guang-Jie Zhai(翟光杰)

    1Key Laboratory of Electronics and Information Technology for Space Systems,National Space Science Center,Chinese Academy of Sciences,Beijing 100190,China

    2School of Physics and Electronics,Shandong Normal University,Ji’nan 250014,China

    3University of Chinese Academy of Sciences,Beijing 100049,China

    Keywords: compressed sensing,imaging quality,spatial frequency domain,multi-scale modulation

    1. Introduction

    Compressed sensing (CS)-based optical imaging has attracted significant interests in recent years. The CS theory is a new sampling theorem proposed by Donoho, Tao, and Cand`es,[1,2]which breaks the Nyquist sampling theorem,with sub-sampling in the signal measurement.A great benefit of CS to imaging is single-pixel imaging feasibility, with the image spatially modulated and then measured by a point detector.[3]After several modulations and measurements,the image is reconstructed based on a CS algorithm using the modulation matrix and point-detection results. Since CS imaging dramatically lowers the array detector requirement, it is already widely applied in many areas such as imaging spectrum,[4,5]3D Lidar,[6,7]Terahertz,[8]and biological imaging.[9]

    A practical problem for CS imaging, however, is the imaging quality. Since the image is obtained by postprocessing of the CS algorithm rather than direct measurement, a reconstruction error is introduced. Although the CS theorem mathematically proves that the signal is accurately recovered under restricted conditions,[10]achieving CS imaging with identical quality as traditional imaging is difficult in practice. Many approaches to resolve this problem exist,with optimizing the measurement matrix design among the most important. The CS theorem requires the measurement matrix to satisfy the restricted isometry property (RIP),[10,11]which implies uncorrelation of the measurement matrix and sparse representation basis. At first,the Gaussian random matrix and the binary random matrix were often employed as measurement matrix to satisfy the RIP criterion. However, the general random matrix was later realized to involve limitations.Through proper design of statistical properties of the measurement matrix, the CS imaging quality can be improved.Therefore, the Optimized Hadamard matrix, Fourier matrix,and other structured matrices were advanced to replace the Gaussian random measurement matrix.[12–15]Besides, Yu et al. derived a complementary modulation matrix that improved the quality of CS telescope imaging.[16]Czajkowski et al. designed the sampling functions by convolving Morlet wavelets with white Gaussian noise.[17]

    The measurement matrix resolution is also related to the CS imaging quality. The CS imaging involves only a point detector, with the imaging resolution actually provided by spatial modulation. Therefore, in a general situation, high spatial resolution modulation is preferred to increase the imaging resolution. However, some studies reveal a trade-off between imaging resolution and reconstruction quality.[18,19]For a fixed scene, finer spatial modulation means more pixels are needed for reconstruction. In CS imaging using light intensities measured through analog-to-digital conversion, a high imaging pixel number amplifies quantization noise and decreases the imaging quality, and in photon-counting region the increase of pixel size has the similar influence. Therefore, some previous studies considered measurement matrix with different spatial frequencies.[20–22]For instance, Fourier single-pixel imaging method[23,24]is equivalent to use a lowpass spatial filter to measure the low frequency part of the image, and in Hadamard single-pixel imaging,[24]the blocks of different sizes in matrix correspond to measurement with different spatial resolutions. However, these methods tend to collect low frequency information if sub-sampling is needed,which means the highest resolution cannot be achieved. The matrix in adaptive CS imaging can be optimized to have varying resolutions in different positions according to the target features,[25]but it requires different spatial frequency components to be separated in space, which may not be applicable to many real scenes. Overall, there is no multi-scale design which can achieve global highest resolution with subsampling.

    In this study, we examine the relationship between the measurement matrix resolution and reconstruction qualities of different image spatial frequencies. We also improve the design of the complementary matrix to achieve multi-scale modulation in the spatial frequency domain. Further, an improved reconstruction method is proposed for image reconstruction for different spatial frequency components. Through the multi-scale reconstruction, the imaging quality is effectively improved compared to traditional single-scale CS imaging,with global highest resolution guaranteed.

    The remainder of this paper is organized as follows: Section 2 contains information on how the resolution of the measurement matrix affects the quality of the CS imaging. In Section 3,the multi-scale measurement matrix design and reconstruction method are proposed and results from numerical simulations are presented. In Section 4, an experimental system set to demonstrate the effect of multi-scale CS imaging is reported,while the main conclusions of the study are summarized in Section 5.

    2. Relationship between measurement matrix resolution and CS imaging quality

    Basic CS imaging can be described through the simple relationship expressed as

    where x is the target,Φ is the measurement matrix with each row representing a spatial modulation of the target x, and y is a column matrix denoting the measurement results. After sub-sampling,the original signal x is recovered by solving the equation according to the measurement result y and measurement matrix Φ.

    In the traditional CS imaging,a measurement matrix of 1-pixel resolution is employed to guarantee the imaging resolution.However,in the present research,we design the measurement matrices with different resolutions(e.g.,n=1–4 pixels)and investigate their effects on the imaging quality through the numerical simulations. A matrix resolution of n pixels means that in each spatial modulation,the matrix is divided into several n×n blocks,and the values in each block are identical.To enhance the modulation randomness,the block division has a random starting point in the matrix. Several examples of measurement matrices with 2-pixel resolution are shown in the following Fig.1.

    Fig. 1. Examples of measurement matrices with 2-pixel resolution. The initial positions of the 2×2 blocks in panels(a),(b),and(c)are different.

    The imaging results for measurement matrices of different resolutions are shown in Fig.2. Rather than simply calculating the quality of the entire reconstructed image as in traditional CS imaging, the imaging quality is considered in the spatial frequency domain. Figure 3 shows the image divided into four spatial frequencies by extracting different regions of the Fourier spectrum. The terms Fre 1 to Fre 4 represent the highest to lowest spatial frequency parts, the boundaries of which are determined as follows. If the width of a 1-pixel point is broadened to n pixels after low-pass filters with the cutoff frequency of fn, then fnis the boundary of Fre n and Fre (n-1). By this means, the terms Fre 1 to Fre 4 corresponding to images with resolutions of 1–2 pixels,2–3 pixels,3–4 pixels, and >4 pixels. The reconstructed errors for the spatial frequencies are presented in Table 1. The maximum of the original image has been normalized to 1, so the RMSEs can be expressed as percentages.

    Fig.2. Imaging results involving measurement matrices of different resolutions. The measurement matrix resolutions for panels(a)–(d)are 1,2,3,and 4 pixels,respectively.

    Fig.3.The reconstructed image is divided into 4 spatial frequencies through Fourier transform.The terms Fre 1 to Fre 4 are in the order of high frequency to low frequency.

    Table 1. Reconstructed errors of different spatial frequencies.

    From the simulation results, the 1-pixel resolution measurement matrix allows adequate reconstruction of the high frequency component but exhibits relatively poor quality in the low frequency component. In contrast, decreasing the measurement matrix resolution is beneficial for proper reconstruction of the low frequency image. Precisely, the spatial frequency Fre n is best reconstructed when the measurement matrix resolution is n pixel.This observation differs from conventional understanding of measurement matrix design,as the single-scale matrix cannot ensure the images have the best reconstruction qualities in both low and high frequencies.

    The regularization revealed by the simulation is further explained below. First, the 1-pixel resolution modulation is obviously necessary for the reconstruction of a high spatial frequency image. Second, when the low spatial frequency component of the image is reconstructed, the effective pixels are lower than in the original image. For example, for a 2-pixel resolution 256×256 image,the effective pixels recovered are 128×128 pixels. Therefore, a low-resolution measurement matrix is applicable for reconstructing a low spatial frequency image at a higher sampling rate compared with the high-resolution matrix. Specifically,the n-pixel measurement matrix achieves the highest sampling rate,with enough imaging resolution in reconstructing an n-pixel resolution image.We notice that the advantage of 1-pixel matrix in high frequency reconstruction is not particular obvious.That is mainly because the Lenna image is heavily represented by lower frequency coefficients and the high frequency information is relative less.

    3. Multi-scale modulation and reconstructionbased CS imaging

    3.1. Multi-scale measurement matrix design

    In this study,we improve the measurement matrix design based on complementary compressed imaging,to increase the CS imaging quality. The complementary compressed imaging model can be expressed as

    Let y=y1-y2,Φ=Φ1-Φ2,and e=e1-e2,equation(4)can be written using the typical compressed sampling formula y =Φx+e. The measurement matrix is changed from the 0–1 binary matrix to the±1 binary matrix with mean of 0 to satisfy the RIP criterion.

    A problem with the traditional complementary CS imaging is that the final measurement matrix is the single-scale 1-pixel resolution matrix. This is only suitable for reconstructing the high frequency image but not the low frequency component according to the regularization mentioned in Section 2. Here, we improve this by designing a multi-scale measurement matrix suitable for simultaneous high-and lowresolution modulations. We initially constructed a 2-pixel resolution p–q random binary matrix A2with identical values in each 2×2 area.Without loss of generality,we selected p=0.5 and q=1 in our study as displayed in Fig. 4, with the initial positions of the 2×2 blocks remaining random as in Fig. 1.Then, A2was split into complementary matrices B1and B2that obeyed the 0–p or 0–q Bernoulli distribution in every 2×2 area. Repeating this process with n random matrices A2produced n pairs of matrices B1and B2.

    Fig.4. Design of multi-scale matrices A1 and A2 representing 1-pixel and 2-pixel resolutions,respectively,which are realized by summation and subtraction of complementary matrix pairs B1 and B2.

    This means that the 1-pixel and 2-pixel resolution measurements are achieved concurrently. Finally,we used the 1-pixel and 2-pixel resolution matrices to reconstruct the target and obtained two values ?x1and ?x2, respectively. According to Section 2,the high frequency reconstruction quality is better in ?x1,while ?x2shows higher quality in the low frequency reconstruction.By extracting and superimposing the high frequency component from the 1-pixel matrix reconstruction result ?x1and low frequency component from the 2-pixel matrix reconstruction result ?x2, a new image reconstruction is obtained,with overall quality better than ?x1and ?x2.

    Different from traditional complementary single-pixel CS imaging using matrices which are complementary in the entire region,the matrices in the multi-scale modulation are complementary only in each 2×2 region. Therefore,a grayscale matrix can be obtained by adding the complementary matrix pair.On the contrary,the addition of traditional complementary matrices can only generate an all 1s matrix which provides little information for the reconstruction. In fact, in the traditional complementary CS imaging,2k times sampling is needed for the measurement matrices B1and B2,with the final effective sampling number only k times. Obviously, some information on the 2k times measurements is wasted. Through the multiscale matrix design, both the addition and subtraction of the k pairs measurements contribute to the image reconstruction,and the information of complementary modulations is fully utilized.

    3.2. Numerical simulations of multi-scale reconstructionbased CS imaging

    In the simulations, we used the images sized 256×256,measuring the images using the multi-scale matrices designed in Section 3.1.Then,the entire image was reconstructed based on Eqs. (6a) and (6b). The sampling rate for each image is 25%.

    The original 256×256 sized images are displayed in Figs. 5(a)–5(c), while the reconstruction results with the 1-pixel resolution measurement matrix according to Eq. (6a)are exhibited in Figs. 5(d)–5(f). Also, the reconstruction results for the 2-pixel resolution measurement matrix according to Eq. (6b) are displayed in Figs. 5(g)–5(i). Finally, the multi-scale reconstruction images are shown in Figs. 5(j)–5(l). The red regions in the images contain abundant lowfrequency components,and these are reconstructed better with the 2-pixel resolution matrix. The blue regions represent more high spatial frequency components,which exhibit better quality with the 1-pixel resolution matrix.

    The RMSEs of the reconstructed images for different spatial frequency components are presented in Table 2. The effects of the multi-scale reconstruction on different targets look similar, coinciding with our theoretical analysis. Clearly, the images reconstructed by the 1-pixel matrix produce lower RMSEs than those reconstructed by the 2-pixel matrix at high frequency, whereas images reconstructed by the 2-pixel matrix show lower RMSEs in the low frequency components. By extracting and superimposing the high and low spatial frequency components,the RMSE of multi-scale reconstructed image is lower than those from single-scale reconstruction for the 1-pixel and 2-pixel resolutions.

    Fig. 5. Simulation results using the proposed method for (a)–(c) original images;(d)–(f)images reconstructed by the 1-pixel resolution measurement matrix;(g)–(i)images reconstructed by the 2-pixel resolution measurement matrix;and(j)–(l)multi-scale reconstruction images.

    Table 2. Reconstructed errors of different spatial frequencies with singlescale and multi-scale reconstruction.

    4. Experiment results

    The CS imaging setup we built based on a digital micromirror device(DMD)to confirm the practical impact of our proposed multi-scale reconstruction is depicted in Fig. 6. To quantitatively evaluate the imaging quality,a digital target instead of a real target is employed to accurately retain the original image. The light source is a halogen lamp and the DMD contains 1024×768 micromirrors, with each mirror of size 13.68 μm×13.68 μm. Each micromirror can be modulated to alter its flip angle separately,reflecting light in different directions to achieve spatial modulation of the light intensity.In the experiment,we used 2×2 micromirrors as a unit to modulate a pixel of the image,so the actual size of the modulation pixel is 27.36 μm×27.36 μm. The DMD flips according to the dot multiplication of the measurement matrix and digital target,equivalent to that the target is imaged on the DMD and then modulated by the measurement matrix. The reflected light is collected by a charge-coupled device (CCD) via a collecting lens L1. However, the CCD is only used as a bucket, that is,as a point detector, as the gray values of every pixels in each sampling are added to obtain the total light intensity of once measurement. The original image sized 64×64 utilized in the experiment is shown in Fig. 7(a). The lines at the bottom of the image involve width and separation of 1 pixel, which can emerge as high spatial frequency targets. In the top of the image, lines with width and separation of 2 pixels are present,and display more low spatial frequency components.

    Fig. 6. Experimental setup for verification of multi-scale CS imaging. The dot multiplication of the measurement matrix and digital target is loaded on the DMD to realize imaging of a virtual digital target.

    In the experiment,we used 1000 pairs of complementary measurement matrices with the sampling rate of 25%. The image reconstructed by the 1-pixel resolution matrix is shown in Fig.7(b),while that reconstructed by the 2-pixel resolution matrix is displayed in Fig. 7(c), and the superimposed multiscale reconstructed image is exhibited in Fig. 7(d). Here, we use“1-pixel reconstruction”for the image reconstructed by 1-pixel measurement matrix, and 2-pixel reconstruction means the image reconstructed by 2-pixel matrix.

    From the experimental results of the 2-pixel measurement matrix, the 1-pixel line is indistinguishable, so obviously the reconstruction quality is poor. However, the reconstructed 2-pixel lines and Tai Chi are smoother, with less noise compared with the 1-pixel reconstruction result. This is because the 1-pixel lines represent almost pure high frequency information, while the 2-pixel lines and Tai Chi contain abundant low frequency components. Although the 1-pixel matrix provides enough resolution for distinguishing the 1-pixel and 2-pixel lines, the overall imaging quality is unsatisfactory. By combining the high frequency component of the 1-pixel reconstruction with the low frequency component of the 2-pixel reconstruction,the image with the best reconstruction quality for high and low spatial frequencies is obtained.

    Fig. 7. Experimental results involving 1000 pairs of complementary measurements for(a)original image;(b)image reconstructed by the 1-pixel matrix; (c) image reconstructed by the 2-pixel matrix; and (d) superimposed multi-scale image.

    For objective evaluation,data for the RMSEs for images in Fig. 7 are presented in Table 3. The RMSE of the high frequency component in the 1-pixel reconstruction is 11.09%,while that for the 2-pixel reconstruction is 14.59%. In the low frequency component, the RMSE of the 2-pixel component is 6.45%,which is much lower than the RMSE of the 1-pixel reconstruction. These results are consistent with expectations. The total RMSE of the multi-scale reconstruction is only 12.82%,which is dramatically improved compared with that for the traditional complementary CS imaging based on the 1-pixel reconstruction(15.75%).

    Table 3. Reconstruction RMSEs of the experimental results.

    Furthermore,we compared the traditional single-scale CS imaging with multi-scale CS imaging under different sampling rates,with the reconstruction results presented in Table 4. The sampling rate is calculated according to the pair number of the complementary matrix,same as in traditional complementary CS imaging. The data reveals that our method is more effective at a low sampling rate. When the sampling rate η is 15%to 25%,the RMSEs of the low frequency component for the 2-pixel reconstruction changes from 8.88%to 6.45%,and is onsistently lower than those for the 1-pixel reconstruction.When η increases to 50%,the RMSEs of the 2-pixel reconstruction in the low frequency component decreases slowly from 6.45%to 5.66%, but the RMSEs of the 1-pixel reconstruction decrease significantly from 11.20% to 2.41%, becoming lower than those for the 2-pixel reconstruction. This is mainly because the 1-pixel measurement matrix already contains enough information for reconstructing all spatial frequency components of the image when the sampling rate increases to a certain degree. Therefore,the advantage displayed by the 2-pixel matrix for restoring the low frequency component is inexistent. When the imaging is performed at a low sampling rate,the information measured by the 1-pixel matrix is not enough,and thus,the 2-pixel matrix exhibits more advantages in lowfrequency reconstruction.

    Table 4. RMSEs of the experimental results under different sampling rates.

    The relationship between the sampling rate and reconstruction quality in different frequency components using traditional 1-pixel reconstruction,2-pixel reconstruction,and the multi-scale method are displayed in Fig. 8. The reconstruction RMSEs of high- and low- frequency components under different sampling rates are shown in Figs. 8(a) and 8(b), respectively. The reconstruction errors of total images containing both high- and low- frequencies are show in Fig. 8(c).The blue, red, and green curves represent the images reconstructed by 1-pixel resolution measurement matrix, 2-pixel resolution matrix,and multi-scale matrix,respectively.For the reconstruction error of the high frequency component shown in Fig. 8(a), the 2-pixel reconstruction is always higher than that for the 1-pixel,especially when the sampling rate is high.This is easily understood as the 2-pixel matrix cannot provide enough resolution to recover details of 1 pixel. In Fig. 8(b),the advantage of the 2-pixel matrix in restoring the low frequency component is more obvious for the lower sampling rate,whereas as η increases gradually,the 1-pixel reconstruction exhibits a lower RMSE limit than that of the 2-pixel reconstruction. Therefore, the advantages of our multi-scale method differ from those of the traditional single-scale 1-pixel reconstruction, according to the sampling rate. As shown in Fig.8(c),the traditional method is more effective at high sampling rates, while our method is better if the sampling rate is low. For different targets, although the critical sampling rate may vary,the qualitative trend will remain unchanged.

    Fig. 8. Relationship between the sampling rate and reconstruction RMSE for (a) the high-frequency component reconstruction; (b) the low-frequency component reconstruction; and (c) the total reconstruction containing both high- and low-frequencies. The blue, red, and green curves represent the measurement matrix of 1-pixel resolution,2-pixel resolution,and multi-scale matrix,respectively.

    5. Conclusion

    In summary, we developed a new method for CS imaging that involved modulating and reconstructing an image in multi-scale spatial frequencies. We showed that the resolution of the measurement matrix was related to the imaging quality of different spatial frequency components. High-resolution modulation was useful for high frequency reconstruction and vice versa. Based on these results, we designed an improved multi-scale complementary measurement matrix capable of achieving high- and low-resolution modulations simultaneously. Compared with the traditional single-pixel imaging which usually reconstructs an image only in the single scale of one pixel, the multi-scale matrix ensured the best reconstruction qualities for the low- and high-frequency components.The simulations and experiment results proved that the proposed multi-scale reconstruction method effectively decreased the RMSE of the reconstructed image.

    Imaging time is a key performance of compressive imaging. To realize multi-scale modulation,the measurement matrix needs to be a grayscale matrix,for example,a three-value matrix in our design. This can be realized by adjusting the residence times of micromirrors in two directions. As in most applications the limited target intensity restricts the modulation rate, the high speed of DMD (up to 20 kHz) is redundant and the grayscale modulation does not actually affect the sampling time. The algorithm for multi-scale reconstruction can be any ones used in traditional CS imaging,so the stability and reconstruction time under different sampling rates are similar to traditional CS imaging. Because the reconstructions need to be performed with 1-pixel and 2-pixel resolution matrices separately in our method,the total reconstruction time is about double of that in traditional CS imaging,but this can be resolved by increasing computing resource or parallel processing. Moreover, the multi-scale reconstruction exhibits more benefits at low sampling rate,so it is significant to increase the sampling speed in compressive imaging. In practice,this can be widely exploited in biological and dynamic compressive imaging involving limited sampling times.

    In this study, we achieved two-scale modulation and reconstruction to verify the idea of improving measurement matrix from the perspective of spatial frequency. Achieving additional scales of modulation at the same sampling number will further enhance the imaging quality. This is the focus of our next study.

    Acknowledgment

    Xue-Feng Liu thanks Amit Ashok for helpful discussions.

    猜你喜歡
    旭日雪峰
    Single exposure passive three-dimensional information reconstruction based on an ordinary imaging system
    雨過(guò)冬晴見(jiàn)黃金一葉獨(dú)掛枝頭
    FSAE賽車(chē)轉(zhuǎn)向系統(tǒng)優(yōu)化設(shè)計(jì)
    要退休了
    雜文月刊(2019年19期)2019-12-04 07:48:34
    日本將刊登“旭日旗”韓文說(shuō)明
    晨釣
    看山是山?看山非山?
    雪峰下的草場(chǎng)
    王雪峰國(guó)畫(huà)
    歌海(2016年1期)2016-03-28 10:08:55
    韓雪峰的“臺(tái)賬”
    av免费观看日本| 久久久久久久久中文| 建设人人有责人人尽责人人享有的 | 亚洲图色成人| 欧美一区二区国产精品久久精品| 色5月婷婷丁香| 美女被艹到高潮喷水动态| 日本黄色片子视频| 久99久视频精品免费| 国产成人aa在线观看| 精品无人区乱码1区二区| 亚洲欧美精品综合久久99| 中文字幕亚洲精品专区| 一级毛片aaaaaa免费看小| 成人高潮视频无遮挡免费网站| av在线蜜桃| 久久精品91蜜桃| 赤兔流量卡办理| 欧美另类亚洲清纯唯美| 精品国内亚洲2022精品成人| 欧美日韩国产亚洲二区| 一级黄色大片毛片| 夜夜看夜夜爽夜夜摸| av在线观看视频网站免费| 婷婷色综合大香蕉| 少妇的逼好多水| 欧美成人a在线观看| 国产美女午夜福利| 嘟嘟电影网在线观看| 夜夜爽夜夜爽视频| 久久久久久久久久久免费av| 久久热精品热| 精品人妻视频免费看| 97人妻精品一区二区三区麻豆| 自拍偷自拍亚洲精品老妇| 男人舔女人下体高潮全视频| 亚洲国产精品专区欧美| 亚州av有码| 97在线视频观看| 青青草视频在线视频观看| 亚洲高清免费不卡视频| 久久久久网色| 国产精品一二三区在线看| 欧美极品一区二区三区四区| 午夜视频国产福利| 91午夜精品亚洲一区二区三区| 亚洲国产欧美人成| 成年免费大片在线观看| 白带黄色成豆腐渣| 国产 一区精品| 日本免费a在线| 建设人人有责人人尽责人人享有的 | 麻豆一二三区av精品| 久久久久久九九精品二区国产| 少妇的逼好多水| 在线免费十八禁| 天堂网av新在线| 亚洲中文字幕一区二区三区有码在线看| 亚洲精品乱久久久久久| 51国产日韩欧美| 国产伦精品一区二区三区四那| 午夜福利网站1000一区二区三区| 久久精品综合一区二区三区| 成人鲁丝片一二三区免费| 熟妇人妻久久中文字幕3abv| 日韩精品有码人妻一区| 91久久精品国产一区二区成人| 久热久热在线精品观看| 国产精品1区2区在线观看.| 亚洲美女视频黄频| 国产成人福利小说| 亚洲国产精品成人综合色| 久久久色成人| 久久久a久久爽久久v久久| 亚洲,欧美,日韩| 人人妻人人澡人人爽人人夜夜 | 成年av动漫网址| 国产又黄又爽又无遮挡在线| 欧美日本视频| 久久久久久久久久久免费av| 久久久午夜欧美精品| 亚洲综合色惰| 淫秽高清视频在线观看| 好男人视频免费观看在线| 国产黄色视频一区二区在线观看 | 日本午夜av视频| 在线免费观看的www视频| 舔av片在线| 午夜福利视频1000在线观看| 欧美成人精品欧美一级黄| 国产免费又黄又爽又色| 99久久九九国产精品国产免费| 少妇人妻一区二区三区视频| 有码 亚洲区| 亚洲熟妇中文字幕五十中出| 美女大奶头视频| 亚洲精品乱久久久久久| 国产亚洲午夜精品一区二区久久 | 精品久久久久久久久亚洲| 中文字幕av成人在线电影| 国产av在哪里看| 亚洲精品乱久久久久久| 夜夜爽夜夜爽视频| 蜜桃久久精品国产亚洲av| 午夜视频国产福利| 免费av观看视频| a级一级毛片免费在线观看| 不卡视频在线观看欧美| 精品久久久噜噜| 亚洲国产成人一精品久久久| 欧美成人午夜免费资源| 最近最新中文字幕大全电影3| 国产探花在线观看一区二区| 久久久久久久久久成人| 十八禁国产超污无遮挡网站| 久久99热这里只频精品6学生 | 秋霞在线观看毛片| 亚洲美女搞黄在线观看| 国产午夜精品一二区理论片| 精品一区二区三区人妻视频| 国内揄拍国产精品人妻在线| 在线a可以看的网站| 97在线视频观看| 色综合亚洲欧美另类图片| 亚洲综合色惰| 亚洲激情五月婷婷啪啪| 精品久久久久久久久亚洲| 天堂√8在线中文| 女人十人毛片免费观看3o分钟| 午夜老司机福利剧场| 国产69精品久久久久777片| 国产色爽女视频免费观看| www.av在线官网国产| 两个人的视频大全免费| 看非洲黑人一级黄片| 我的女老师完整版在线观看| 听说在线观看完整版免费高清| 99久国产av精品国产电影| 精品久久久久久久久久久久久| 国产成人一区二区在线| 欧美三级亚洲精品| 国产精品一及| 午夜久久久久精精品| 国产黄片美女视频| 成人国产麻豆网| av国产久精品久网站免费入址| 国产精品久久电影中文字幕| 久久鲁丝午夜福利片| 2022亚洲国产成人精品| 亚洲美女视频黄频| 九色成人免费人妻av| 亚洲婷婷狠狠爱综合网| 蜜桃亚洲精品一区二区三区| 久久亚洲国产成人精品v| 欧美区成人在线视频| 日韩视频在线欧美| 大又大粗又爽又黄少妇毛片口| 亚洲色图av天堂| av免费观看日本| 综合色av麻豆| 狂野欧美白嫩少妇大欣赏| 联通29元200g的流量卡| 夫妻性生交免费视频一级片| 精品久久久久久电影网 | 亚洲av成人精品一区久久| 不卡视频在线观看欧美| 91aial.com中文字幕在线观看| 变态另类丝袜制服| 男人和女人高潮做爰伦理| 最近最新中文字幕大全电影3| 中文亚洲av片在线观看爽| 亚洲av.av天堂| 激情 狠狠 欧美| 一区二区三区乱码不卡18| 99热这里只有精品一区| 大香蕉久久网| 26uuu在线亚洲综合色| 天美传媒精品一区二区| 男插女下体视频免费在线播放| 欧美丝袜亚洲另类| 午夜免费男女啪啪视频观看| 91久久精品电影网| 又粗又硬又长又爽又黄的视频| 春色校园在线视频观看| 99久久人妻综合| 九草在线视频观看| 少妇被粗大猛烈的视频| 激情 狠狠 欧美| 久久人人爽人人片av| 秋霞伦理黄片| 国内精品美女久久久久久| 亚洲av男天堂| 嫩草影院入口| 亚洲精品乱码久久久v下载方式| 国产毛片a区久久久久| 亚洲欧洲国产日韩| 一个人免费在线观看电影| 精品久久久久久成人av| 国产精品国产三级专区第一集| 男女视频在线观看网站免费| 国产亚洲精品av在线| 婷婷六月久久综合丁香| 非洲黑人性xxxx精品又粗又长| 一个人看的www免费观看视频| 床上黄色一级片| 亚洲欧美精品综合久久99| 三级国产精品欧美在线观看| 一区二区三区免费毛片| 内射极品少妇av片p| 国产精品,欧美在线| 桃色一区二区三区在线观看| 欧美日韩综合久久久久久| 国产视频首页在线观看| 国产 一区 欧美 日韩| 久久99精品国语久久久| 春色校园在线视频观看| 网址你懂的国产日韩在线| 成年女人永久免费观看视频| videos熟女内射| 国产亚洲5aaaaa淫片| 97热精品久久久久久| 久久久欧美国产精品| 精品一区二区免费观看| 久热久热在线精品观看| 亚洲精品久久久久久婷婷小说 | 国产又色又爽无遮挡免| 国产三级中文精品| 成人美女网站在线观看视频| 晚上一个人看的免费电影| 色视频www国产| 99久久成人亚洲精品观看| 亚洲国产精品合色在线| 国产伦在线观看视频一区| 亚洲欧美日韩卡通动漫| 日本欧美国产在线视频| 我要搜黄色片| 三级国产精品片| 青春草国产在线视频| 亚洲最大成人av| 久久久国产成人精品二区| 内射极品少妇av片p| 亚洲成色77777| 内地一区二区视频在线| 午夜视频国产福利| 大香蕉久久网| 欧美日韩在线观看h| 超碰av人人做人人爽久久| 国产v大片淫在线免费观看| 亚洲欧美成人精品一区二区| 熟女电影av网| 我的老师免费观看完整版| 久久久久久久国产电影| 免费一级毛片在线播放高清视频| 精品少妇黑人巨大在线播放 | 欧美精品国产亚洲| 午夜免费激情av| 国产成人a区在线观看| 欧美一区二区精品小视频在线| 精品久久国产蜜桃| 99热6这里只有精品| 国产免费男女视频| 久久久久精品久久久久真实原创| 欧美成人一区二区免费高清观看| 久久人妻av系列| 中文天堂在线官网| 99九九线精品视频在线观看视频| 久久久久久久国产电影| 中文亚洲av片在线观看爽| 日本猛色少妇xxxxx猛交久久| 久久精品国产亚洲网站| 看十八女毛片水多多多| 村上凉子中文字幕在线| 国产黄色视频一区二区在线观看 | 中文字幕精品亚洲无线码一区| 你懂的网址亚洲精品在线观看 | 寂寞人妻少妇视频99o| 丝袜美腿在线中文| 亚洲国产精品合色在线| 久久99热这里只频精品6学生 | 亚洲国产欧美在线一区| 亚洲天堂国产精品一区在线| 菩萨蛮人人尽说江南好唐韦庄 | 男人舔奶头视频| 亚洲精品国产av成人精品| 少妇熟女aⅴ在线视频| 国产精品爽爽va在线观看网站| 99久久人妻综合| 哪个播放器可以免费观看大片| 51国产日韩欧美| 寂寞人妻少妇视频99o| 欧美xxxx黑人xx丫x性爽| 久久鲁丝午夜福利片| 乱人视频在线观看| 国产色婷婷99| 国产精品人妻久久久影院| 在线播放国产精品三级| 搡女人真爽免费视频火全软件| 国产亚洲一区二区精品| 亚洲真实伦在线观看| 在线免费十八禁| av.在线天堂| 好男人视频免费观看在线| 久久鲁丝午夜福利片| 欧美性感艳星| av在线亚洲专区| 成年女人永久免费观看视频| 免费一级毛片在线播放高清视频| 麻豆国产97在线/欧美| 久久久a久久爽久久v久久| 久久婷婷人人爽人人干人人爱| 久热久热在线精品观看| 哪个播放器可以免费观看大片| a级毛色黄片| 男女啪啪激烈高潮av片| 国产黄色视频一区二区在线观看 | 又爽又黄无遮挡网站| 神马国产精品三级电影在线观看| 国产精品美女特级片免费视频播放器| 国产女主播在线喷水免费视频网站 | 国产亚洲精品久久久com| 国语自产精品视频在线第100页| 国产成人精品婷婷| 久久久久性生活片| 真实男女啪啪啪动态图| 中文字幕亚洲精品专区| 成人美女网站在线观看视频| 欧美日本亚洲视频在线播放| 最近手机中文字幕大全| 99热全是精品| 日韩欧美在线乱码| 国产精品一及| 日本欧美国产在线视频| 国产伦精品一区二区三区四那| 日本黄色片子视频| 亚洲精华国产精华液的使用体验| 国产不卡一卡二| 成人性生交大片免费视频hd| 成人高潮视频无遮挡免费网站| 久久久久久久午夜电影| 国产在视频线在精品| 三级国产精品欧美在线观看| 综合色av麻豆| 国产精品蜜桃在线观看| 午夜福利网站1000一区二区三区| 女人十人毛片免费观看3o分钟| 久久精品熟女亚洲av麻豆精品 | 国产精品.久久久| 中文精品一卡2卡3卡4更新| 亚洲人与动物交配视频| 男的添女的下面高潮视频| 亚洲成人中文字幕在线播放| 国产一区二区在线av高清观看| 国产亚洲一区二区精品| 最近中文字幕2019免费版| 亚洲av一区综合| 中文字幕久久专区| 看非洲黑人一级黄片| 亚洲精品影视一区二区三区av| 日韩一区二区视频免费看| 日本午夜av视频| 1024手机看黄色片| 亚洲欧美精品专区久久| 日本av手机在线免费观看| 久久人妻av系列| 欧美日韩在线观看h| 久久久色成人| 日韩,欧美,国产一区二区三区 | 久久久久性生活片| 美女脱内裤让男人舔精品视频| 激情 狠狠 欧美| 国产在视频线在精品| 国产av在哪里看| av女优亚洲男人天堂| 直男gayav资源| 亚洲精品日韩在线中文字幕| 国产精品久久久久久av不卡| 国产亚洲精品av在线| 一边亲一边摸免费视频| 国产大屁股一区二区在线视频| 亚洲激情五月婷婷啪啪| 国产黄色小视频在线观看| 亚洲av男天堂| 成年免费大片在线观看| 国产成人精品久久久久久| 毛片一级片免费看久久久久| 国产精品人妻久久久久久| 99久久中文字幕三级久久日本| 午夜激情福利司机影院| 身体一侧抽搐| 能在线免费观看的黄片| av在线观看视频网站免费| 久久精品久久久久久久性| 久久99蜜桃精品久久| 啦啦啦啦在线视频资源| 3wmmmm亚洲av在线观看| av线在线观看网站| 亚洲18禁久久av| 欧美97在线视频| 麻豆国产97在线/欧美| 波多野结衣巨乳人妻| 乱人视频在线观看| 欧美激情久久久久久爽电影| 性色avwww在线观看| 国产极品精品免费视频能看的| 99在线视频只有这里精品首页| 天天躁夜夜躁狠狠久久av| 少妇熟女欧美另类| 亚洲精品成人久久久久久| 午夜福利成人在线免费观看| 高清日韩中文字幕在线| 99热精品在线国产| 青春草亚洲视频在线观看| 搞女人的毛片| 性插视频无遮挡在线免费观看| 最近中文字幕2019免费版| 一本一本综合久久| 波多野结衣高清无吗| 建设人人有责人人尽责人人享有的 | 亚洲精品色激情综合| 亚洲在线自拍视频| 狂野欧美激情性xxxx在线观看| 高清视频免费观看一区二区 | 亚洲av不卡在线观看| 亚洲欧美成人综合另类久久久 | 欧美性猛交黑人性爽| 色5月婷婷丁香| 亚洲av熟女| 人人妻人人澡人人爽人人夜夜 | 99热这里只有是精品在线观看| 国产 一区精品| 九九爱精品视频在线观看| 我的女老师完整版在线观看| a级一级毛片免费在线观看| 欧美性感艳星| 黄色日韩在线| 一本一本综合久久| 国产av在哪里看| 中文字幕亚洲精品专区| 亚洲一级一片aⅴ在线观看| 日韩视频在线欧美| 欧美一区二区亚洲| 国产精品一区二区在线观看99 | 成人鲁丝片一二三区免费| 看非洲黑人一级黄片| 中文在线观看免费www的网站| 国产精品国产三级国产av玫瑰| 成人欧美大片| 日本爱情动作片www.在线观看| 波野结衣二区三区在线| kizo精华| 两个人视频免费观看高清| 亚洲成人av在线免费| 亚洲欧美日韩高清专用| 日产精品乱码卡一卡2卡三| 亚洲国产色片| 91午夜精品亚洲一区二区三区| 国产69精品久久久久777片| 久久精品夜夜夜夜夜久久蜜豆| 十八禁国产超污无遮挡网站| 欧美成人午夜免费资源| 内射极品少妇av片p| 天天躁日日操中文字幕| 纵有疾风起免费观看全集完整版 | 亚洲国产精品成人久久小说| 日韩欧美精品v在线| av线在线观看网站| 男女边吃奶边做爰视频| 搞女人的毛片| 国产精华一区二区三区| 亚洲真实伦在线观看| 亚洲久久久久久中文字幕| 亚洲精品乱久久久久久| 日韩成人伦理影院| av国产久精品久网站免费入址| 男的添女的下面高潮视频| 国产单亲对白刺激| 七月丁香在线播放| 久久久欧美国产精品| 69av精品久久久久久| 你懂的网址亚洲精品在线观看 | 精品国产三级普通话版| 青春草视频在线免费观看| 免费看美女性在线毛片视频| 不卡视频在线观看欧美| 一个人免费在线观看电影| 欧美日本视频| 男女国产视频网站| 国产色爽女视频免费观看| 最新中文字幕久久久久| 晚上一个人看的免费电影| 神马国产精品三级电影在线观看| 日本欧美国产在线视频| 成人毛片60女人毛片免费| 欧美色视频一区免费| 大话2 男鬼变身卡| 午夜精品在线福利| 午夜福利在线观看免费完整高清在| 美女cb高潮喷水在线观看| 人体艺术视频欧美日本| 亚洲国产最新在线播放| 一级黄色大片毛片| av国产久精品久网站免费入址| av.在线天堂| 国产视频内射| 日韩大片免费观看网站 | 中文精品一卡2卡3卡4更新| 少妇丰满av| 午夜福利成人在线免费观看| 在线观看一区二区三区| 国产成人a∨麻豆精品| 日韩一区二区三区影片| 亚洲精品乱码久久久久久按摩| 国产在视频线精品| 特级一级黄色大片| 国产亚洲最大av| 少妇的逼水好多| 久久久精品94久久精品| 91狼人影院| 嫩草影院新地址| av天堂中文字幕网| 国内精品宾馆在线| 在线观看av片永久免费下载| 三级国产精品片| 色噜噜av男人的天堂激情| 精品人妻一区二区三区麻豆| 国产一区亚洲一区在线观看| 精品久久久久久久人妻蜜臀av| 免费看美女性在线毛片视频| 国内精品一区二区在线观看| 亚洲国产精品成人综合色| 国产老妇女一区| 人妻少妇偷人精品九色| 久久久午夜欧美精品| 亚洲一区高清亚洲精品| 边亲边吃奶的免费视频| 最近2019中文字幕mv第一页| 1000部很黄的大片| 午夜老司机福利剧场| 白带黄色成豆腐渣| 久久综合国产亚洲精品| 国产伦在线观看视频一区| 国产黄色小视频在线观看| 看免费成人av毛片| 国产成人a∨麻豆精品| 精品不卡国产一区二区三区| 亚洲精品国产av成人精品| 精品国产一区二区三区久久久樱花 | 国产美女午夜福利| 久久久久久久久久久免费av| 毛片一级片免费看久久久久| 国语自产精品视频在线第100页| 丝袜美腿在线中文| 国产成人freesex在线| 一区二区三区四区激情视频| 国产精品乱码一区二三区的特点| 有码 亚洲区| a级一级毛片免费在线观看| 亚洲内射少妇av| 秋霞伦理黄片| 18禁在线播放成人免费| 亚洲最大成人中文| 六月丁香七月| av.在线天堂| 春色校园在线视频观看| 国产精品无大码| 亚洲人成网站在线播| 亚洲欧美成人综合另类久久久 | 免费播放大片免费观看视频在线观看 | 插阴视频在线观看视频| 夜夜看夜夜爽夜夜摸| 亚洲国产精品成人久久小说| 亚洲av中文av极速乱| 国产精品一二三区在线看| 97人妻精品一区二区三区麻豆| 亚洲图色成人| 99久久中文字幕三级久久日本| 国产又黄又爽又无遮挡在线| 亚洲内射少妇av| 国产在线男女| 亚洲av一区综合| 亚洲最大成人中文| 欧美成人a在线观看| 夫妻性生交免费视频一级片| 91精品伊人久久大香线蕉| 亚洲色图av天堂| 亚洲va在线va天堂va国产| 国产精品一区二区性色av| 国产精品久久久久久精品电影小说 | 国产探花在线观看一区二区| 国产精品精品国产色婷婷| 亚洲成色77777| 乱系列少妇在线播放| 亚洲国产精品sss在线观看| 天堂影院成人在线观看| 又粗又硬又长又爽又黄的视频| 国产亚洲精品av在线| 毛片一级片免费看久久久久| 国产精品一区二区性色av| 非洲黑人性xxxx精品又粗又长| 国产精品无大码| 国产精品永久免费网站| 久久精品影院6| 亚洲性久久影院| 91精品伊人久久大香线蕉| 亚洲精品一区蜜桃| 两个人的视频大全免费| 51国产日韩欧美| 亚洲自拍偷在线| 国产在视频线在精品| 爱豆传媒免费全集在线观看| 久久久精品大字幕| 最近最新中文字幕免费大全7| 国产色爽女视频免费观看| 国产成人福利小说| 亚洲,欧美,日韩|