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

    Real-Time Position and Attitude Estimation for Homing and Docking of an Autonomous Underwater Vehicle Based on Bionic Polarized Optical Guidance

    2020-09-28 00:49:48CHENGHaoyuanCHUJinkuiZHANGRanGUIXinyuanandTIANLianbiao
    Journal of Ocean University of China 2020年5期

    CHENG Haoyuan, CHU Jinkui, ZHANG Ran, GUI Xinyuan, and TIAN Lianbiao

    Real-Time Position and Attitude Estimation for Homing and Docking of an Autonomous Underwater Vehicle Based on Bionic Polarized Optical Guidance

    CHENG Haoyuan, CHU Jinkui*, ZHANG Ran, GUI Xinyuan, and TIAN Lianbiao

    The Key Laboratory for Micro/Nano Technology and System of Liaoning Province, Dalian University of Technology, Dalian 116024, China

    As an important tool for marine exploration, the autonomous underwater vehicle (AUV) must home in and dock at a docking station (DS) to be recharged, repaired, or to exchange information at set intervals. However, the complex and hostile underwater environment makes this process challenging. This study proposes a real-time method based on polarized optical guidance for determining the position and attitude of the AUV relative to its DS. Four polarized artificial underwater landmarks are positioned at the DS, which are recognized by the AUV vision system. Compared with light intensity, the polarization of a light beam is known to be better maintained at greater propagation distances, especially in underwater environments. The proposed method, which is inspired by the ability of marine animals to communicate, calculates the pose parameters in less than 10ms without any other navigational information. The simulation results reveal that the angle errors are small and the position errors are no more than 0.116m within 100m in the coastal ocean. The results of underwater experiments further demonstrate the feasibility of the proposed method, which extends the operating distance of the AUV beyond what is currently possible while maintaining the precision of traditional optical guidance.

    polarization; optical guidance; AUV; underwater docking; position and attitude estimation

    1 Introduction

    As a powerful tool for ocean exploration and the development of sub-sea infrastructure, the autonomous underwater vehicle (AUV) is playing an increasingly important role in military and scientific research. However, due to its limited operational hours, the AUV must be frequently recharged and relaunched and the complex and hostile underwater environment increases the difficulty of this process. Therefore, a reliable homing and docking technology is crucial for the AUV to be able to complete its missions.

    To determine its position and attitude relative to the docking station (DS), the AUV can be equipped with various sensors,, acoustic, electromagnetic, and optical. Compared to other technologies, optical guidance has good directional accuracy, low vulnerability to external detection, and a diverse range of uses (Deltheil, 2002). The optical docking system provides a targeting accuracy on the order of 1cm in real-life conditions, even in turbid bay waters (Cowen, 1997). However, its operating distance is limited because of the underwater conditions. Yu(2001) proposed a navigation method for the AUV based on artificial underwater landmarks (AULs), which can be recognized by a vision system installed in anAUV. Lee(2003) presented a docking system for AUVs that enable them to dock at an underwater station with the aid of a camera installed at the center of the AUV nose. A docking method based on image processing for cruising-type AUVs was proposed by Park(2009), and a similar docking method for hovering-type AUVs was proposed by Kondo(2012). Maki(2013) proposed a docking method for hovering-type AUVs based on both acoustic and visual positioning, and Yahya and Arshad (2015) proposed a tracking method based on the placement of colored markers in the cluttered environment. Vallicrosa(2016) developed a light beacon navigation system for estimating the position of the DS with respect to the AUV.

    All the above mentioned optical guidance systems for AUV homing and docking have established that optical guidance technology is highly feasible and has reliable precision, especially at short distances. However, underwater imaging poses significant challenges at extended ranges compared to similar situations over the ocean surface. Even in the clearest ocean waters, the visibility range, at best, is on the order of 10m owing to the absorption and scattering of light by particles of various origins, including algal cells, detritus, sediments, plankton, and bubbles (Hou, 2009). Therefore, extending the operating range while maintaining the precision of the optical guidance system is the main technological challenge of AUVs. As an information carrier, the characteristics of polarized light, such as the degree of polarization (DOP), have greater robustness than light intensity (Xu, 2015). More than 70 underwater species (Horváth and Varjú, 2004) are polarization-sensitive and some have the ability to navigate, hunt, camouflage themselves, and communicatepolarized vision (Shashar, 2000; Waterman, 2006; Cartron, 2013). Polarization has been found to increase the contrast and detection range in scattering media both above and under water (Schechner, 2003; Schechner and Karpel, 2005). In addition, bionic polarization navigation technology has the potential for application in water (Cheng, 2020).

    In this paper, we propose a real-time method for homing and docking the AUV by determining its position and attitude relative to polarized AULs on the DS, which are recognized by a vision system installed in the AUV. We introduce the theory of polarization imaging, along with the concept of polarized optical guidance. Then, we obtain atmospheric intensity and DOP images and simulate underwater intensity and DOP images at various distances. Pixel-recognition errors are input into the model to calculate the pose errors. We conduct underwater experiments to determine the practicability of the proposed method and the results indicate that the polarized optical guidance system extends the operating distance while maintaining the precision of traditional optical guidance. This biologically inspired polarized optical guidance mechanism provides a novel approach for AUV homing and docking that enables AUVs to perform tasks at much greater distances.

    2 Methods

    2.1 Theory of Polarization Imaging

    Generally, vision systems recognize the intensity characteristic of AULs, whereas the proposed method recognizes the DOP of the polarized AULs. The DOP is the ratio of the component magnitude of the polarized light to the total magnitude of the light, which can vary between 0 (completely unpolarized light) and 100% (completely polarized light).To obtain a DOP image, the DOP of every pixel of the image is calculated and the Stokes vectoris used to describe the polarized light beam. All information about the polarized light is contained in the four components of Eq. (1):

    whereis the total intensity of the light,is the fraction of linear polarization parallel to a reference plane,is the proportion of linear polarization at 45? with respect to the reference plane, andis the fraction of right-handed circular polarization.The linear polarization is used to format the image. The polarization state of the incident light=[,,,]T, which is changed by a polarizing film, can be expressed based on the Muller matrix:

    whereis the angle between the main optical axis and the zero reference line and=[,,,]Tis the polarization state of the emergent light. Here the first row of the Muller matrix is relevant because the light intensity can be obtained directly by the camera:

    Therefore, if the light intensities of the emergent light at three differentvaluesare known, the,,andvalues ofthe incident beam can be calculated. If we setto 0?, 45?, and 90?, the following equation set is obtained:

    The equation set is then transformed as follows:

    Then, the linear DOP is calculated and a DOP image is obtained:

    The precision of the homing and docking process depends on the image quality. To better evaluate the image quality, we introduce the parameter:

    wheretis the intensity or DOP of the target andbis the intensity or DOP of the background in the image.

    2.2 Polarized Optical Guidance

    Polarization imaging, which enables the contrast of target images to be increased and both the illumination effects and backscattering to be decreased, is an ideal technique for underwater target detection, as it is simple to implement, requires little power, and can be used in real-time applications (Dubreuil, 2013).In addition, in the deep sea, there is no daylight, which means point-like bioluminescence is the only light source (Warrant and Locket, 2004). To survive in the dark ocean environment, marine creatures have evolved to acquire polarization sensitivity, which they use for navigation, predation, disguise, and communication in the same way as humans use intensity information. Inspired by these marine creatures, we design a method for determining the position and attitude of the AUV based on polarized AULs, which are recognized by a camera installed in the AUV, as shown in Fig.1. A light source with a wavelength of 532nm is employed because of its greater propagation distance. We use a normal lightsource rather than laser light because AUVs have difficulty recognizing laser light, which is unidirectional, despite its greater propagation distance.

    The proposed system consists of a camera and four AULs, which are polarized light beacons situated in the DS. As the light source is polarized, its DOP is greater than that of the background. The pixel center of each AUL can be identified by finding the maximum DOP point in the image, for which the corresponding pixel coordinates can be obtained. Figs.2 and 3 show the DS, image, and the camera coordinate system at long and short distances, respectively. Fig.4 shows the positions of the four light beacons in the image during the docking process. When the AUV is far away from the DS, the four light beacons appear to be locatedinoneplace on the edge of the AUV vision system, as shown in Figs.2 and 4(a). The angles of rotation,andabout the-axis and-axis, respectively can be calculated as shown in Eq. (8):

    where x0 and y0 are the respective horizontal and vertical coordinates of the camera coordinate system, which can be calculated based on the physical pixel size of the camera and the coordinates of the image coordinate system. The pixel coordinates of the light beacons are obtained by image processing. f is the focal length of the camera.

    Fig.2 Coordinate system of polarized optical guidance at a large distance.

    Fig.3 Coordinate system of polarized optical guidance at a close distance.

    Fig.4 Positions of four light beacons in the image during the docking process. (a), (b), (c), and (d) represent different stages of the docking process.

    At a distance, the AUV approaches the DS based on two angles. As the AUV nears the DS, the four polarized light beacons appear in the AUV vision system, as shown in Figs.3, 4(b), and 4(c). Safe and accurate homing and docking of the AUV depends on the precise determina-tion of the position and attitude of the AUV relative to the DS. It is a perspective--point (PNP) problem, in which the position and attitude of the AUV are determined on the relationship of the positions of the polarized AULs in the image and actual situation whenthe number of AULsisin computervision. Fischler and Bolles (1987) proved that solving a PNP problem and obtaining a unique solution requires at least four coplanar feature points. The PNP problem simplifies the proposed method into a three-dimensional quadratic system of equations:

    where,, andarethe respective anglesamong OcA, OcB, and OcC.,, andare the respective distances between A, B, and C, which are known beforehand. Then,,,and, which are the respective lengths of AOc, BOc, and COc,are known. However, the equation has multiple solutions, and therefore, a fourth point is needed to obtain a unique solution. Coupled with the coordinates for A, B, and C in the camera coordinate system, the coordinates for A, B, and Cin the camera coordinate system are obtained. Finally, we use a matrixforcoordinatetransformation to solve for the six position and attitude parameters of the camera relative to the AULs, as shown below:

    where the positions of the camera relative to the AULs areT,T, andT, and the attitudes are,, and. (,,) and (,,) are the coordinates of the AULs relative to the DS and the camera coordinate system. Fig.5 shows the complete solution process, for which the total computing time is no longer than 10ms. When the AUV is close to the DS, the four light beacons are outside the image, as shown in Fig.4(d). This distance is too short to adjust the attitude of the AUV, so the AUV simply maintains its current direction and homes in.

    The proposed real-time processing algorithm for polarized optical guidance is mainly based on the PNP problem and has a fast calculation speed because of its lower computational complexity. Due to its iterative nature,the algorithm has strong robustness and high precision. It requires only the relative position of each AUL and nothing more. However, the model errors are sensitive to discrepancies between the actual and estimated centers of each AUL.

    2.3 Experiment to Demonstrate the Superiority

    To demonstrate the superior performance of the proposed method, we must prove that the use of polarized optical guidance extends the operating distance while maintaining the precision of traditional optical guidance. In a field experiment, we used a camera to acquire measurements at the roof of a building at night, to obtain intensity and DOP images at 1-m interval at distances less than 30m and 5-m interval at distances greater than 30m. The exposure time was 50ms, which was kept constant for each acquisition. The camera had 2448×2048 pixels and a lens focal length of 10.5mm. The physical size of each pixel was 3.45μm. The AULs comprised four LEDs and linear polarizing films.The power of each AUL was 5mW and the docking target was a square with side lengths of 0.15m. In the experiment, the camera represented the AUV, and the AULs represented the DS. Fig.6 shows a sketch of the optics. The unpolarized light source became polarized by the polarizing film. Intensity data was captured by the camera and the polarization was calculated using a PC. We note that the light pollution in the experimental environment was negligible.

    First, we obtained a series of images of the AULs at different distances up to 45m. Then, we corrected the distortion of the images and simulated the effects of water using an image degradation model (Hufnagel and Stanley, 1964) for oceanic water (Mobley, 1994) with a low concentration of chlorophyl and a moderate concentration of hydrosols. The water type was the coastal ocean for which the absorption coefficient was 0.18m?1and the scattering coefficient was 0.22m?1(Petzold, 1972). Next, the algorithm obtained the pixel-recognition errors of the center of every AUL, which are the errors between the actual and identified pixel centers of the AULs. The averaged pixel-recognition errors were entered into the calculation model, which calculated the errors in the position and attitude of the camera relative to the AULs. To further demonstratethefeasibility of the proposed method underwater, after waterproofing all the equipment, we placed the experimental setup in seawater at a depth of 0.5m.

    Fig.5 Flow chart of solution process for six pose parameters.

    Fig.6 Schematic of the optics of the experiment.

    3 Results and Discussion

    In this section, we describe the simulations and experiments we conducted to illustrate the advantages of polarized versus unpolarized optical guidance. The quality of the DOP images was found to be better than the intensity images both above and under water, so we established that polarization imaging technology extends the operating distance of underwater optical guidance. To simulate the total pose errors of the proposed method, we obtained the pixel-recognition errors of the center of the AULs in the image and proved that underwater polarized optical guidance maintains the precision of traditional optical guidance. We conducted underwater experiments to evaluate the performance of the system and the feasibility of the proposed method.

    3.1 Quality of Intensity and DOP Images

    Fig.7 shows the intensity and DOP images of the atmosphericexperiment. At a distance of 5m, the AULs in the intensity image were barely visible, as depicted in Fig.7(a), whereas in the DOP image, they were clearly visible, as depicted in Fig.7(b). At a distance of 20m, the AULs in the intensity image were invisible, as shown in Fig.7(c), whereas in the DOP image, they were still visible, as shown in Fig.7(d). In the intensity images, the small red circles on the right enclose the light sources, and the big circles on the left are enlarged versions of the same. Table 1 shows the quality of the intensity and DOP images obtained in the atmosphericexperiment, which was calculated using Eq. (7). The above images were then simulated underwater, the quality of which are shown in Table 2. In Tables 1 and 2, we can see improvement with respect to the quality of the DOP images relative to intensity images. The image quality and improvement in Tables1 and 2, as obtained by Eq. (7), are dimensionless. These results establish that the quality of polarization imaging is obviously improved, especially under water and at a large distance. Thus, we have proved that polarization imaging technology hasa longer operating distance than traditional optical guidance technology.

    Fig.7 Atmospheric experiment. (a), intensity image at a distance of 5m; (b), DOP image at a distance of 5m; (c), intensity image at a distance of 20m; and (d), DOP image at a distance of 20m.

    Table 1 Image quality of intensity and DOP above water

    Table 2 Image quality of intensity and DOP under water

    3.2 Precision of Polarized Optical Guidance

    Using the algorithm, we calibrated the atmospheric images and obtained the pixel-recognition errors of the center of the AULs in the images. We acquired 10 images and obtained 40 pixel-recognition errors for every test position. We then calculated the mean errors, and Fig.8 shows the relationship between the distance and the error. At distances less than 3m and greater than 25m, the pixel errors are high due to overexposure and camera performance, respectively. This means the errors obtained at distances between 3m and 25m are valuable, which areshown in Fig.9. Next, we simulated the atmospheric images underwater and obtained the corresponding errors, which are also shown in Fig.9. Both these errors featured some distortion, but the underwater errors were the larger of the two. Owing to the scattering and absorption of particles in water, it is more difficult to recognize the center of the AUL underwater. We found the pixel-recognition error to be inversely proportional to the distance of the beacon from the camera. The closer was the beacon to the camera, the larger was the projected light in the image, which led to greater errors than those of beacons that were far away. These results are consistent with those reported in previous work (Vallicrosa, 2016). At the same time, if the light source had higher power or a bigger luminous area, the results would have been better at greater distances and worse at shorter distances due to overexposure and the large imaging area. In contrast, if the light source had lower power or a smaller luminous area, the camera would experience performance challenges at large distances.

    Fig.8 Mean pixel errors within 45m.

    Fig.9 Pixel errors of atmospheric experiment and underwater simulation between 3 and 25m and their fitting results.

    Fig.9 shows plots in which the atmospheric and underwater errors are fitted, wherein we can see that because of the water degradation, the pixel errors underwater were greater than those in the atmosphere. Lastly, after fitting the pixel errors at different distances, we input them into the AUV pose model to calculate the six parameters of the position and attitude. Due to the iterative nature of the algorithm, the angle and-direction distance errors were small. The total errors mainly depended on the distances along theanddirections. Fig.10 shows the simulation results of the relationship between the distance and the total errors of the AUV, in which we can see that the position errors of both the atmospheric and underwater images increased with the distance between the AUV and the DS. Within 100m, the atmospheric errors were no greater than 0.093m and the underwater errors were no greater than 0.116m. This degree of precision is consistent with the results obtained in previous work (Hong, 2003) using unpolarized optical guidance.

    Fig.10 Total pose errors of the AUV within 100m.

    3.3 Underwater Experiment

    The underwater experiment was conducted in the coastal ocean in the afternoon to be consistent with the conditions used in the simulation. The underwater light field was polarized due to the refracted polarized skylight formedby the atmospheric scattering. However, the DOP of the underwater light field was significantly lower than that of the polarized AULs and had little impact on the experiment. Despite the absence of any interference from fish, seaweed, or other marine creatures, the underwater environment in the experiment was complex, much like the real AUV working environment. The water was semi-turbid and the visibility was poor owing to the absorption and scattering of light by various particles. Long-term underwater work would seriously shorten the life of the equipment because of seawater corrosion and pressure.The ocean waves, shoreline,seabed, and waterproof de-vice also had a negative influence on the experiment due to their impacts on the light. Except for the imaging method, the experimental conditions of the intensity and DOP images were the same. Fig.11 shows the intensity and DOP images of the polarized AULs in the underwater experiment, for which the distance between the beacon and the camera was 2m. In the intensity image, the light source was invisible, but its position is indicated by a red circle.In the DOP image, the light source was clear and successfully recognized, and the polarized AULs were four bright points that had higher DOP values than the background. The underwater DOP image had more noisy points owing to the scattering of various particles, and the calculation accuracy was poorer than the atmospheric image because of the complex underwater optical characteristics. There were also other bright points above the AULs that were darker than the polarized AULs, which were formed by the reflections of the water surface in the DOPimage. The total underwater pose error was 0.011m, which isconsistent with that of the simulation. This result verifies the superior performance of the proposed system and the feasibility of this method underwater.

    Fig.11 Underwater experiment. (a) Intensity image and (b) DOP image.

    4 Conclusions

    In this paper, we proposed a real-time and bioinspired method that uses polarized optical guidance to determine the position and attitude of an AUV relative to a DS for homing and docking. This method, based on polarization imaging technology, extends the operating distance while maintaining the precision of traditional optical guidance. The simulation and experimental results obtained in this research serve as a feasibility study and provide the necessary evidence that such a method is warranted. The obtained angle errors were negligible, and the position errors were no greater than 0.116m within 100m in the coastal ocean. A fast calculation speed enables real-time solutions for the AUV position and attitude. Furthermore, this method is simple and inexpensive, requiring only the relative position of each AUL. Future work will focus on accelerating the calculation speed, developing a dynamic model, and conducting long-distance underwater experiments.

    Acknowledgements

    This work was supported by the National Natural Science Foundation of China (NSFC) (Nos. 51675076, 51505062), the Science Fund for Creative Research Groups of NSFC(No.51621064),and the Pre-Research Foundati on of China (No. 61405180102).

    Cartron, L., Josef, N., Lerner, A., Mccusker, S. D., Darmaillacq, A. S., Dickel, L., and Shashar, N., 2013. Polarization vision can improve object detection in turbid waters by cuttlefish., 447 (3): 80-85.

    Cheng, H. Y., Chu, J. K., Zhang, R., Tian, L. B., and Gui, X. Y., 2020. Underwater polarization patterns considering single Rayleigh scattering of water molecules., 41 (13): 4947-4962.

    Cowen, S., Briest, S., and Dombrowski, J., 1997. Underwater docking of autonomous undersea vehicles using optical terminal guidance.. Halifax, 2: 1143-1147.

    Deltheil, C., Didier, L., Hospital, E., and Brutzman, D. P., 2002. Simulating an optical guidance system for the recovery of an unmanned underwater vehicle., 25 (4): 568-574.

    Dubreuil, M., Delrot, P., Leonard, I., Alfalou, A., Brosseau, C., and Dogariu, A., 2013. Exploring underwater target detection by imaging polarimetry and correlation techniques., 52 (5): 997-1005.

    Fischler, M. A., and Bolles, R. C., 1987. Random sample consensus: A paradigm for model fitting with applications to image analysis and automated cartography., 24 (6), 381-395.

    Hong, Y. H., Kim, J. Y., Lee, P. M., Jeon, B. H., Oh, K. H., and Oh, J. H., 2003. Development of the homing and docking algorithm for AUV.. Honolulu, 205-212.

    Horváth, H. G., and Varjú, D., 2004.. Springer, Berlin, 470pp.

    Hou, W., 2009. A simple underwater imaging model.,34 (17): 2688-2690.

    Hufnagel, R. E., and Stanley, N. R., 1964. Modulation transfer function associated with image transmission through turbulent media., 54 (1): 52.

    Kondo, H., Okayama, K., Choi, J. K., Hotta, T., Kondo, M., Okazaki, T., Singh, H., Chao, Z., Nitadori, K., Igarashi, M., and Fukuchi, T., 2012. Passive acoustic and optical guidance for underwater vehicles.. IEEE, Yeosu, 1: 1-6.

    Lee, P. M., Jeon, B. H., and Kim, S. M., 2003. Visual servoing for underwater docking of an autonomous underwater vehicle with one camera.San Diego,2: 677-682.

    Maki, T., Shiroku, R., Sato, Y., Matsuda, T., Sakamaki, T., and Ura, T., 2013. Docking method for hovering type AUVs by acoustic and visual positioning., Tokyo, 1-6.

    Mobley, C. D., 1994.. Academic, San Diego, 592pp.

    Park, J. Y., Jun, B. H., Lee, P. M., and Oh, J., 2009. Experiments on vision guided docking of an autonomous underwater vehicle using one camera., 36 (1): 48-61.

    Petzold, T. J., 1972.. Scripps Institution of Oceanography La Jolla Ca Visibility Lab, La Jolla, 79pp.

    Schechner, Y. Y., and Karpel, N., 2005. Recovery of underwater visibility and structure by polarization analysis., 30 (3): 570-587.

    Schechner, Y. Y., Narasimhan, S. G., and Nayar, S. K., 2003. Polarization-based vision through haze., 42 (3): 511-525.

    Shashar, N., Hagan, R., Boal, J. G., and Hanlon, R. T., 2000. Cuttlefish use polarization sensitivity in predation on silvery fish., 40 (1): 71-75.

    Vallicrosa, G., Bosch, J., Palomeras, N., Ridao, P., Carreras, M., and Gracias, N., 2016. Autonomous homing and docking for AUVs using range-only localization and light beacons., 49 (23): 54-60.

    Warrant, E. J., and Locket, N. A., 2004. Vision in the deep sea., 79 (3): 671-712.

    Waterman, T. H., 2006. Reviving a neglected celestial underwater polarization compass for aquatic animals., 81 (1): 111-115.

    Xu, Q., Guo, Z. Y., Tao, Q. Q., Jiao, W. Y., Wang, X. S., Qu, S. L., and Gao, J., 2015. Transmitting characteristics of polarization information under seawater., 54 (21): 6584.

    Yahya, M. F., and Arshad, M. R., 2015. Tracking of multiple markers based on color for visual servo control in underwater docking..George Town, 482-487.

    Yu, S. C., Ura, T., Fujii, T., and Kondo, H., 2001. Navigation of autonomous underwater vehicles based on artificial underwater landmarks.Honolulu, 1: 409-416.

    . E-mail: chujk@dlut.edu.cn

    December 8, 2019;

    May 13, 2020;

    May 18, 2020

    (Edited by Xie Jun)

    亚洲国产色片| 深夜精品福利| 精品视频人人做人人爽| 在线观看人妻少妇| 天天影视国产精品| 韩国精品一区二区三区 | 久久久久久久久久久久大奶| 美女大奶头黄色视频| 天天操日日干夜夜撸| 视频在线观看一区二区三区| 成人黄色视频免费在线看| 一区二区三区精品91| 亚洲天堂av无毛| 国产麻豆69| av又黄又爽大尺度在线免费看| 国产日韩一区二区三区精品不卡| 久久99热6这里只有精品| 欧美变态另类bdsm刘玥| 欧美精品一区二区大全| 男女边吃奶边做爰视频| 亚洲精品成人av观看孕妇| 午夜福利乱码中文字幕| 亚洲一区二区三区欧美精品| 久久影院123| 亚洲精品成人av观看孕妇| 午夜福利乱码中文字幕| 最新的欧美精品一区二区| 午夜激情久久久久久久| 水蜜桃什么品种好| 亚洲欧美成人精品一区二区| av在线观看视频网站免费| 欧美另类一区| 久久久久国产精品人妻一区二区| 精品亚洲乱码少妇综合久久| 成人午夜精彩视频在线观看| 人妻少妇偷人精品九色| 五月伊人婷婷丁香| 亚洲精品久久久久久婷婷小说| 波多野结衣一区麻豆| 日日撸夜夜添| 国产精品三级大全| 国产1区2区3区精品| 99久久综合免费| 少妇人妻精品综合一区二区| 99热全是精品| 一级毛片电影观看| 美女国产视频在线观看| 国产一区亚洲一区在线观看| 99久久中文字幕三级久久日本| av一本久久久久| 久久青草综合色| 亚洲av免费高清在线观看| 啦啦啦啦在线视频资源| 精品亚洲成国产av| 国产精品不卡视频一区二区| 看十八女毛片水多多多| 亚洲精品aⅴ在线观看| 成人国产麻豆网| 国产片内射在线| 亚洲成人手机| 97精品久久久久久久久久精品| 黑人巨大精品欧美一区二区蜜桃 | 观看av在线不卡| 国产成人91sexporn| 王馨瑶露胸无遮挡在线观看| 国产一区有黄有色的免费视频| 免费观看在线日韩| 国产精品免费大片| 在线观看www视频免费| 色5月婷婷丁香| 亚洲精品,欧美精品| 在线观看国产h片| 欧美成人午夜精品| 大香蕉久久网| 少妇的丰满在线观看| 五月天丁香电影| 美女中出高潮动态图| 国产女主播在线喷水免费视频网站| 中文精品一卡2卡3卡4更新| av不卡在线播放| 欧美97在线视频| 最近最新中文字幕大全免费视频 | 久久av网站| 亚洲美女黄色视频免费看| 天堂8中文在线网| 丝袜在线中文字幕| 人体艺术视频欧美日本| 激情五月婷婷亚洲| 国产日韩欧美视频二区| videosex国产| 亚洲综合色网址| 人妻少妇偷人精品九色| 日本黄色日本黄色录像| 中国三级夫妇交换| 亚洲,一卡二卡三卡| 国精品久久久久久国模美| 国产精品.久久久| 国产黄色视频一区二区在线观看| 天美传媒精品一区二区| 日日摸夜夜添夜夜爱| 一级毛片我不卡| 亚洲人与动物交配视频| 亚洲av在线观看美女高潮| 日韩人妻精品一区2区三区| av福利片在线| 又粗又硬又长又爽又黄的视频| 国产av国产精品国产| 久久99热6这里只有精品| 热99久久久久精品小说推荐| 久久免费观看电影| 91久久精品国产一区二区三区| 一二三四在线观看免费中文在 | 久久精品国产亚洲av涩爱| 99热国产这里只有精品6| 成人无遮挡网站| 国产高清三级在线| 日产精品乱码卡一卡2卡三| 人人澡人人妻人| 纯流量卡能插随身wifi吗| 久久久久久久久久久免费av| 日韩av免费高清视频| 久久人人爽人人片av| 99re6热这里在线精品视频| 久久精品国产自在天天线| 一级毛片我不卡| 国产淫语在线视频| 波多野结衣一区麻豆| 精品国产国语对白av| 欧美97在线视频| 最新的欧美精品一区二区| 久久99热这里只频精品6学生| 国产在线一区二区三区精| 在线看a的网站| 精品一区在线观看国产| 久久精品国产a三级三级三级| 亚洲第一区二区三区不卡| 一边亲一边摸免费视频| 亚洲精品国产av蜜桃| 捣出白浆h1v1| 天堂中文最新版在线下载| 青春草国产在线视频| 熟女人妻精品中文字幕| 9热在线视频观看99| 亚洲国产欧美日韩在线播放| 青春草视频在线免费观看| 熟女av电影| 欧美xxⅹ黑人| 边亲边吃奶的免费视频| 热99久久久久精品小说推荐| 99热全是精品| 国产一级毛片在线| 免费看不卡的av| 纵有疾风起免费观看全集完整版| 一区二区av电影网| 亚洲成人手机| 高清毛片免费看| 国产一区二区三区综合在线观看 | 男男h啪啪无遮挡| 飞空精品影院首页| 国产不卡av网站在线观看| 边亲边吃奶的免费视频| 五月玫瑰六月丁香| 咕卡用的链子| 精品人妻一区二区三区麻豆| 国产精品一区二区在线观看99| 黄片无遮挡物在线观看| av不卡在线播放| 日韩成人伦理影院| 十分钟在线观看高清视频www| 久久久久精品性色| 国产在线视频一区二区| 精品久久国产蜜桃| 老熟女久久久| 汤姆久久久久久久影院中文字幕| 欧美日韩国产mv在线观看视频| 人妻 亚洲 视频| 欧美日韩精品成人综合77777| 精品久久久精品久久久| 十八禁网站网址无遮挡| 日本猛色少妇xxxxx猛交久久| 亚洲精品456在线播放app| 99国产综合亚洲精品| 欧美国产精品一级二级三级| 久久久精品区二区三区| 国产毛片在线视频| 熟女电影av网| 丝袜喷水一区| 巨乳人妻的诱惑在线观看| 国产男女超爽视频在线观看| 日本与韩国留学比较| 夫妻午夜视频| 老司机亚洲免费影院| 亚洲少妇的诱惑av| 大陆偷拍与自拍| 菩萨蛮人人尽说江南好唐韦庄| 午夜福利视频精品| 九九爱精品视频在线观看| 黄色配什么色好看| 亚洲天堂av无毛| 免费观看性生交大片5| 国产成人欧美| 免费看光身美女| 美女主播在线视频| 黄色怎么调成土黄色| 最近2019中文字幕mv第一页| 久久久精品免费免费高清| 在线看a的网站| 自拍欧美九色日韩亚洲蝌蚪91| 欧美国产精品va在线观看不卡| 又黄又爽又刺激的免费视频.| 日本黄大片高清| 国产欧美另类精品又又久久亚洲欧美| 国产一级毛片在线| 亚洲av中文av极速乱| 久久精品国产综合久久久 | 日韩av免费高清视频| 国产伦理片在线播放av一区| 搡女人真爽免费视频火全软件| 丝袜脚勾引网站| 国产一区二区激情短视频 | 99久国产av精品国产电影| 久久人人97超碰香蕉20202| 有码 亚洲区| 天美传媒精品一区二区| 亚洲精品456在线播放app| 免费看光身美女| 亚洲国产av新网站| 欧美性感艳星| a级毛片黄视频| 新久久久久国产一级毛片| 国产有黄有色有爽视频| 欧美bdsm另类| 丰满少妇做爰视频| 精品亚洲成国产av| 国产免费福利视频在线观看| 国产成人精品婷婷| 亚洲国产精品成人久久小说| 日韩欧美一区视频在线观看| 亚洲激情五月婷婷啪啪| 成人18禁高潮啪啪吃奶动态图| 妹子高潮喷水视频| www.av在线官网国产| 免费在线观看黄色视频的| 亚洲成人手机| 一个人免费看片子| 男的添女的下面高潮视频| 国产精品三级大全| 最近的中文字幕免费完整| 熟女av电影| 日韩三级伦理在线观看| 亚洲伊人色综图| 欧美亚洲 丝袜 人妻 在线| 高清欧美精品videossex| 精品少妇黑人巨大在线播放| 国产高清国产精品国产三级| 日本欧美视频一区| 99热网站在线观看| 欧美成人午夜免费资源| 欧美日韩视频精品一区| a级毛色黄片| 在线看a的网站| av卡一久久| 日韩成人伦理影院| 国产亚洲最大av| 一级爰片在线观看| 日韩大片免费观看网站| 高清在线视频一区二区三区| 久久精品国产自在天天线| 最近的中文字幕免费完整| 狠狠婷婷综合久久久久久88av| 亚洲精品久久午夜乱码| 青春草视频在线免费观看| 美女大奶头黄色视频| 好男人视频免费观看在线| 日韩在线高清观看一区二区三区| 在线观看三级黄色| 久久久久久人妻| 免费av不卡在线播放| 久久人人爽人人爽人人片va| 人妻人人澡人人爽人人| 高清欧美精品videossex| 十分钟在线观看高清视频www| 免费大片18禁| 蜜桃国产av成人99| 国产免费现黄频在线看| 一本—道久久a久久精品蜜桃钙片| 9热在线视频观看99| 深夜精品福利| 精品国产一区二区久久| 大香蕉97超碰在线| 午夜视频国产福利| 亚洲色图 男人天堂 中文字幕 | 欧美xxxx性猛交bbbb| 亚洲在久久综合| www.色视频.com| 美女国产高潮福利片在线看| 国产毛片在线视频| 人人妻人人澡人人看| a级毛片黄视频| 亚洲久久久国产精品| 十分钟在线观看高清视频www| 97精品久久久久久久久久精品| 欧美激情极品国产一区二区三区 | 一级爰片在线观看| 一区二区三区精品91| 久久影院123| 两个人看的免费小视频| 看免费av毛片| 国产淫语在线视频| 国产国语露脸激情在线看| 久久久久精品性色| 免费在线观看完整版高清| 人人妻人人澡人人看| 一边摸一边做爽爽视频免费| 亚洲成人av在线免费| 亚洲一级一片aⅴ在线观看| 曰老女人黄片| 97超碰精品成人国产| 亚洲精品日韩在线中文字幕| 日日爽夜夜爽网站| 国产亚洲av片在线观看秒播厂| 精品熟女少妇av免费看| 精品久久国产蜜桃| 激情视频va一区二区三区| 久久久久国产网址| 色吧在线观看| 亚洲五月色婷婷综合| 国产白丝娇喘喷水9色精品| 亚洲国产欧美在线一区| 国产精品一二三区在线看| 亚洲av电影在线进入| 日本猛色少妇xxxxx猛交久久| 亚洲精品456在线播放app| 精品国产露脸久久av麻豆| 黑人猛操日本美女一级片| 亚洲国产欧美日韩在线播放| 欧美xxⅹ黑人| 青青草视频在线视频观看| 在线观看人妻少妇| 久久99蜜桃精品久久| av电影中文网址| 国产毛片在线视频| av片东京热男人的天堂| 极品少妇高潮喷水抽搐| 久久久久久伊人网av| 韩国精品一区二区三区 | 另类精品久久| 亚洲美女黄色视频免费看| 国产乱来视频区| 1024视频免费在线观看| 日韩一区二区三区影片| 国产精品偷伦视频观看了| 国产日韩欧美视频二区| 久久国产精品大桥未久av| 黑人欧美特级aaaaaa片| 久久精品国产亚洲av天美| 欧美最新免费一区二区三区| 大片电影免费在线观看免费| 中文字幕制服av| av网站免费在线观看视频| 欧美少妇被猛烈插入视频| 欧美 亚洲 国产 日韩一| 国产av一区二区精品久久| 日韩中字成人| a级毛片黄视频| 少妇高潮的动态图| 免费久久久久久久精品成人欧美视频 | 一级毛片电影观看| av线在线观看网站| 男女边摸边吃奶| 99九九在线精品视频| 97超碰精品成人国产| 桃花免费在线播放| 黄网站色视频无遮挡免费观看| 女性生殖器流出的白浆| 国产片内射在线| 日本91视频免费播放| 国产精品一国产av| 下体分泌物呈黄色| 免费大片黄手机在线观看| 黄色怎么调成土黄色| 亚洲精品美女久久久久99蜜臀 | 精品国产露脸久久av麻豆| 精品久久久精品久久久| 永久网站在线| 国产精品久久久久久精品电影小说| 边亲边吃奶的免费视频| 久久人人爽人人片av| 欧美亚洲 丝袜 人妻 在线| 国产精品一二三区在线看| 黄片无遮挡物在线观看| 三上悠亚av全集在线观看| 视频区图区小说| 自线自在国产av| 菩萨蛮人人尽说江南好唐韦庄| 亚洲少妇的诱惑av| 在线观看三级黄色| av视频免费观看在线观看| 国产视频首页在线观看| 亚洲性久久影院| 人人澡人人妻人| 久久久久精品人妻al黑| 美女内射精品一级片tv| 欧美xxⅹ黑人| 久久久久国产精品人妻一区二区| 91aial.com中文字幕在线观看| 国产69精品久久久久777片| 哪个播放器可以免费观看大片| 国产欧美日韩一区二区三区在线| 少妇人妻久久综合中文| 国产又爽黄色视频| 国国产精品蜜臀av免费| 亚洲,欧美,日韩| 一区二区三区四区激情视频| 久久久国产精品麻豆| 日韩欧美精品免费久久| 亚洲欧美一区二区三区黑人 | 一区二区日韩欧美中文字幕 | 男女无遮挡免费网站观看| 免费久久久久久久精品成人欧美视频 | 永久免费av网站大全| 日韩 亚洲 欧美在线| 久久免费观看电影| 亚洲精品中文字幕在线视频| 草草在线视频免费看| 国产69精品久久久久777片| 精品福利永久在线观看| 99热这里只有是精品在线观看| av免费观看日本| 乱人伦中国视频| 日韩大片免费观看网站| 考比视频在线观看| 精品少妇久久久久久888优播| 女性生殖器流出的白浆| 国产亚洲精品久久久com| 母亲3免费完整高清在线观看 | 女的被弄到高潮叫床怎么办| 在线观看免费日韩欧美大片| 午夜老司机福利剧场| 国产精品一区二区在线不卡| 天堂俺去俺来也www色官网| 在线天堂最新版资源| xxxhd国产人妻xxx| 久久狼人影院| 两性夫妻黄色片 | 日日啪夜夜爽| 欧美bdsm另类| 亚洲精品456在线播放app| 蜜臀久久99精品久久宅男| 亚洲国产色片| 亚洲第一区二区三区不卡| 9191精品国产免费久久| 国产一区二区激情短视频 | 成人国产av品久久久| 男人操女人黄网站| 亚洲少妇的诱惑av| 97在线视频观看| 99re6热这里在线精品视频| 热re99久久国产66热| 制服丝袜香蕉在线| 建设人人有责人人尽责人人享有的| 中文字幕人妻熟女乱码| 久久久久久久亚洲中文字幕| 久久久久人妻精品一区果冻| 亚洲经典国产精华液单| 曰老女人黄片| 一区在线观看完整版| av黄色大香蕉| 国产片内射在线| 黄色配什么色好看| 99久国产av精品国产电影| 2022亚洲国产成人精品| 国产欧美日韩一区二区三区在线| 大码成人一级视频| 国产毛片在线视频| 日韩欧美精品免费久久| 亚洲激情五月婷婷啪啪| 中文字幕免费在线视频6| 18禁裸乳无遮挡动漫免费视频| 国产av一区二区精品久久| 91在线精品国自产拍蜜月| 全区人妻精品视频| 一边摸一边做爽爽视频免费| 亚洲国产av新网站| 精品人妻在线不人妻| 久久青草综合色| 国产精品久久久久久久电影| 在线免费观看不下载黄p国产| av福利片在线| 男人添女人高潮全过程视频| xxxhd国产人妻xxx| 午夜福利在线观看免费完整高清在| 成年人午夜在线观看视频| 免费高清在线观看日韩| 欧美xxxx性猛交bbbb| 少妇被粗大猛烈的视频| 午夜影院在线不卡| 久久人人爽av亚洲精品天堂| 最近中文字幕2019免费版| 一二三四在线观看免费中文在 | 麻豆精品久久久久久蜜桃| 日本wwww免费看| 日韩精品免费视频一区二区三区 | 80岁老熟妇乱子伦牲交| 在线观看免费高清a一片| 亚洲第一区二区三区不卡| 国产精品一二三区在线看| 少妇人妻久久综合中文| 国产成人午夜福利电影在线观看| 大片免费播放器 马上看| 亚洲av欧美aⅴ国产| 女的被弄到高潮叫床怎么办| 久久人人爽人人爽人人片va| 日韩,欧美,国产一区二区三区| 2021少妇久久久久久久久久久| 国产黄频视频在线观看| 国产高清三级在线| 在线观看美女被高潮喷水网站| 国产有黄有色有爽视频| 性色av一级| 99热网站在线观看| 18禁国产床啪视频网站| 亚洲婷婷狠狠爱综合网| 男女边摸边吃奶| 韩国精品一区二区三区 | 精品酒店卫生间| 国产精品久久久久久精品电影小说| 999精品在线视频| 一本—道久久a久久精品蜜桃钙片| 欧美老熟妇乱子伦牲交| 亚洲综合精品二区| 黄片无遮挡物在线观看| 一级爰片在线观看| 久久99热这里只频精品6学生| 亚洲欧美一区二区三区黑人 | 久久鲁丝午夜福利片| 在现免费观看毛片| 中文字幕人妻熟女乱码| 久久婷婷青草| 热re99久久国产66热| 91午夜精品亚洲一区二区三区| 黑人巨大精品欧美一区二区蜜桃 | 成人免费观看视频高清| 国产精品一区www在线观看| 久久精品久久久久久噜噜老黄| 亚洲精品久久成人aⅴ小说| 在线观看免费高清a一片| 男女啪啪激烈高潮av片| 国产深夜福利视频在线观看| 狂野欧美激情性xxxx在线观看| 成人亚洲欧美一区二区av| a级毛色黄片| 各种免费的搞黄视频| 欧美+日韩+精品| 老司机亚洲免费影院| 99香蕉大伊视频| 男女无遮挡免费网站观看| 一二三四中文在线观看免费高清| a级毛片在线看网站| 久久精品国产鲁丝片午夜精品| av天堂久久9| 最近中文字幕2019免费版| 激情视频va一区二区三区| 精品久久蜜臀av无| 黄色配什么色好看| 欧美xxxx性猛交bbbb| 亚洲av男天堂| 久久久久久久久久成人| 99久国产av精品国产电影| 中文字幕人妻熟女乱码| 考比视频在线观看| 一区二区三区乱码不卡18| 国产精品无大码| 久久狼人影院| 亚洲高清免费不卡视频| 一区二区三区精品91| 日本欧美视频一区| 91精品国产国语对白视频| 国产黄色免费在线视频| 亚洲国产毛片av蜜桃av| 亚洲,一卡二卡三卡| 男女啪啪激烈高潮av片| 丁香六月天网| 免费观看性生交大片5| 日韩人妻精品一区2区三区| 国产乱人偷精品视频| 黄片播放在线免费| 边亲边吃奶的免费视频| 久久精品国产鲁丝片午夜精品| 久久人人爽av亚洲精品天堂| 午夜福利网站1000一区二区三区| 精品99又大又爽又粗少妇毛片| 伦理电影免费视频| 你懂的网址亚洲精品在线观看| 国产精品蜜桃在线观看| av一本久久久久| 免费观看av网站的网址| 狠狠精品人妻久久久久久综合| 中文字幕最新亚洲高清| 亚洲欧美一区二区三区黑人 | 色网站视频免费| 国产无遮挡羞羞视频在线观看| 少妇的逼好多水| 精品国产一区二区三区久久久樱花| 国产精品久久久久久精品古装| 九色成人免费人妻av| 夫妻性生交免费视频一级片| 亚洲成人一二三区av| 高清av免费在线| 99国产综合亚洲精品| 欧美 日韩 精品 国产| 成人影院久久| 高清黄色对白视频在线免费看| 亚洲精品一区蜜桃| 亚洲精品久久午夜乱码|