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

    Light pen whole space coordinate measuring system based on a tracking turntable

    2021-12-21 13:34:26WANGSenLIUShuguiMAOQing

    WANG Sen, LIU Shugui, MAO Qing

    (State Key Laboratory of Precision Measuring Technology and Instruments, Tianjin University, Tianjin 300072, China)

    Abstract: Light pen coordinate measuring system (LPCMS) is a kind of portable coordinate measuring technique based on vision metrology. In classical LPCMS, the measuring range is limited to the camera’s field of view. To overcome this defect, a new LPCMS is designed in this paper to fulfil whole space coordinate measurement. The camera is installed on a turntable instead of a tripod, so that the camera can rotate to track the movement of the light pen. The new system can be applied to large scale onsite measurement, and therefore it notably extends the application of LPCMS. To guarantee the accuracy of the new system, a method to calibrate the parameters of the tracking turntable is also proposed. Fixing the light pen at a stationary position, and changing the azimuth angles of the turntable’s two shafts, so that the camera can capture the images of the light pen from different view angles. According to the invariant spatial relationship between the camera and the pedestal of the tracking turntable, a system of nonlinear equations can be established to solve the parameters of the turntable. Experimental results show that the whole space coordinate measuring accuracy of the new system can reach 0.25 mm within 10 m. It can be concluded that the newly designed system can significantly expand the measuring range of LPCMS without losing too much accuracy.

    Key words: light pen; whole space coordinate measuring system; tracking turntable; vision metrology; parameter calibration; homogeneous coordinate transformation matrix

    0 Introduction

    Coordinate measuring is a very important technique in manufacturing industry and reverse engineering[1-3]. With the continuous development of the manufacturing industry, coordinate measuring technique is increasingly applied to on-site measurement of workpieces. To meet this demand, several kinds of portable coordinate measuring systems have been designed in the last three decades[4-7]. Among these systems, light pen coordinate measuring system is indispensable for its ability to measure the hidden features of workpieces[8-9].

    Light pen coordinate measuring system (LPCMS) is developed on the basis of vision metrology principle. It has the advantages of good portability, high measuring accuracy and low environmental requirements, which makes it competent for on-site measurement.

    When three dimensional object space is mapped to two dimensional image space through the camera’s perspective projection model, the depth information is lost[10]. Therefore, it is impossible to recover the three dimensional information of an object just from one single image. LPCMS utilizes control point technique to provide additional constraints to the perspective projection model, so that it is able to measure three dimensional coordinates from the image of the light pen.

    Despite classical LPCMS has many good features, it has its own limitation. The most salient one is that its measuring range is limited to the camera’s field of view, whose shape is a rectangular pyramid. In practical application, the position of the camera needs to change frequently according to the positon of the measured object, so that the object is guaranteed in the camera’s field of view all the time. Once the camera is moved, the reference coordinate system of LPCMS changes. Therefore these measuring results cannot be unified in the same reference coordinate system. In other words, it is impossible for the classical LPCMS to fulfil whole space coordinate measurement.

    To overcome the defect of the classical LPCMS, a new light pen whole space coordinate measuring system based on a tracking turntable is designed in this paper. By mounting it to a tracking turntable, the camera can rotate driven by the turntable to track the movement the light pen[11-13]. Thus the measuring range of the newly designed system is no longer confined to the camera’s field of view, and therefore can realize whole space coordinate measurement. This new system can be applied to on-site measurement of large workpieces such as ship hull, train body and plane shell, which exceeds the measuring range of the classical LPCMS.

    In the newly designed light pen whole space coordinate measuring system, the positon of the camera is constantly changing with the movement of the light pen during the measuring process. Therefore, the reference coordinate system of the new system cannot be fixed to the camera any more. Instead, the coordinate system which is fasten to the pedestal of the turntable can serve as the new reference coordinate system. To transform the final measuring results to the coordinates in the new reference coordinate system, the coordinate transformation matrix between the camera and the turntable’s pedestal must be determined. For this purpose, the parameters of the tracking turntable must be carefully calibrated to guarantee the measuring accuracy of the new system in the whole space.

    To calibrate the parameters of the tracking turntable, a method that utilizes a stationary light pen is also proposed in this paper. By changing the azimuth angles of the turntable’s two rotation shafts, the camera can capture the image of the stationary light pen from different view directions. A series of nonlinear equations can be established according to the invariance of the coordinate transformation matrix between the light pen and the turntable’s pedestal. The parameters of the tracking turntable can be solved from this equation system by nonlinear least squares algorithm.

    1 Detailed system description

    1.1 System configuration

    Fig.1 shows the configuration of the light pen whole space coordinate measuring system based on a tracking turntable designed in this paper. From Fig.1, it can be seen that the system mainly consists of 4 components: a light pen, a laptop, an industrial camera and a tracking turntable. In addition to these, there are some wires and accessories as well.

    Fig.1 Configuration of light pen whole space coordinate measuring system based on a tracking turntable

    Fig.2 shows the structure of the light pen adopted in this paper. It uses luminous infrared LEDs as control points, and there are 13 control points distributing on the light pen. The relative spatial positions between these control points need to be accurately calibrated in advance to provide accurate additional constraints to the camera’s perspective projection model. Besides the control points, there is a contact probe at the end of the light pen as well. The positon of the center of the probe tip also needs to be calibrated in advance to obtain its position relative to all the control points.

    The laptop is used to run the measuring software. Besides the mathematical calculation of the whole measuring system, the measuring software is also in charge of the communication with the light pen and the turntable. According to actual requirements, the laptop can control the rotation angles of the turntable’s two shafts through RS232 serial communication, and control the luminescence time and the electric current of the infrared LEDs through wireless communication.

    The industrial camera is used to capture the image of the light pen. Both the intrinsic parameter matrix of the camera and the distortion factors of the lens need to be calibrated in advance to guarantee the accuracy of the camera’s imaging model. To track the movement of the light pen, the camera is installed on a tracking turntable instead of a tripod.

    Fig.3 gives the structure of the tracking turntable. The turntable consists of a stationary pedestal, two precision rotation shafts and two high precision angular encoders. The camera is installed on a rectangular metal frame which can rotate driven by these two shafts. The vertical shaft can drive the camera conduct left and right deflection in the horizontal plane, and it is calledAaxis in this paper. The horizontal shaft can drive the camera conduct up and down pitching in the vertical plane, and it is calledBaxis in this paper. Under ideal condition, theAaxis and theBaxis of the turntable are strictly perpendicularly intersecting. However, because of the existence of manufacture error and assembly error, they are neither perpendicular nor intersecting. In actual application, theAaxis and theBaxis are treated as two skew lines in different planes, which has no special geometric relationship.

    Fig.3 Structure of tracking turntable

    1.2 Measurement principle

    The new whole space coordinate measuring system designed in this paper is developed on the basis of traditional LPCMS, so it is necessary to give a detailed introduction of the measurement principle of classical LPCMS first.

    For convenience of description, two coordinate systems are established in LPCMS. They are fixed to the light pen and the camera, respectively, and they can be called the light pen coordinate system and the camera coordinate system severally. The light pen coordinate system can be denoted asOl-xlylzl, and the camera coordinate system can be denoted asOc-xcyczc.

    (1)

    whereAl,cis the 4×4 homogeneous coordinate transformation matrix betweenOl-xlylzlandOc-xcyczc;Rl,cis the 3×3 orthogonal rotation matrix, andtl,cis the three dimensional translation vector.

    (2)

    whereKis the intrinsic parameter matrix of the camera;λis an arbitrary scale factor that is nonzero;fuandfvare the effective focal lengths along the two directions of the imaging plane;cuandcvare the pixel coordinates of the principal point, andsis a parameter that reflects the non-perpendicularity of the two axes of the imaging plane.

    What needs to be emphasized here is that the lens distortion is not taken into consideration, because it is not the main issue of this paper. In practice, the lens distortions in the homogeneous pixel coordinates need to be compensated first to get the undistorted pixel coordinates, so that the linear perspective projection model in Eq.(2) is satisfied.

    When LPCMS is used to measure the coordinates of an object, the camera needs to be located at a proper stationary position, so that the whole object is in the camera’s field of view. The operator holds the light pen and makes the probe tip contact the surface of the object. According to the image of the light pen captured by the camera, the coordinates of the point touched by the light pen can be measured.

    According to Eqs.(1) and (2), the equation satisfying for each control point can be expressed as

    (3)

    After eliminating the scale factorλ, two nonlinear equations can be derived from Eq.(3) for each control point. The unknown variables in these equations are the rotation matrixRl,cand the translation vectortl,cbetweenOl-xlylzlandOc-xcyczc. Merging these equations together, if the number of the control points is big enough, an overdetermined nonlinear equation system can be established.Rl,candtl,ccan be solved from this equation system by various nonlinear least squares algorithms. This is a well-known pose estimation issue in vision metrology, which is named as perspective-n-point (PnP) problem.

    For the reason that the position of the camera is fixed during the whole measuring process, the camera coordinate systemOc-xcyczccan serve as the reference coordinate system of LPCMS. The final measuring result of LPCMS is the coordinates of the measured point in the reference coordinate systemOc-xcyczc.

    (4)

    From Eq.(4), it can be seen that how to determine the homogeneous transformation matrixAc,tbetweenOc-xcyczcandOt-xtytztis very important for the new system. For the reason that theAaxis and theBaxis of the turntable are treated as skew lines in different planes, this transformation process can be decomposed into 4 steps, and two intermediate coordinate systems need to be established to make each step simple and clear. The two intermediate coordinate systems are fixed to and rotate with theAaxis and theBaxis of the tracking turntable, respectively, and they can be denoted as theA-axis coordinate systemOA-xAyAzAand theB-axis coordinate systemOB-xByBzBseverally.

    The two intermediate coordinate systems are established when theAaxis and theBaxis are both at zero position. The way to establishOA-xAyAzAis: take theAaxis as thexaxisxA; take the common perpendicular of theAaxis and theBaxis as thezaxiszA; take the intersection of theAaxis and the common perpendicular as the originOA; theyaxisyAcan be determined according to the right hand rule. The way to establishOB-xByBzBis: take theBaxis as theyaxisyB; take the common perpendicular of theAaxis and theBaxis as thezaxiszB; take the intersection of theBaxis and the common perpendicular as the originOB; thexaxisxBcan be determined according to the right hand rule.

    From the structure of the turntable shown in Fig.3, it can be inferred thatOA-xAyAzAonly rotates with theAaxis andOB-xByBzBrotates with theAaxis and theBaxis at the same time. Suppose the azimuth angles of theAaxis and theBaxis can be denoted asθandφ, respectively, and the two intermediate coordinate systems can be further denoted asOA-xAyAzA(θ) andOB-xByBzB(θ,φ) when the two shafts of the turntable are at any azimuth angles. When theAaxis is at zero position, theA-axis coordinate system can be denoted asOA-xAyAzA(θ=0), and it can serve as the reference coordinate system of the newly designed measuring system. In other words, the coordinate systemsOA-xAyAzA(θ=0) is exactly the turntable coordinate systemOt-xtytzt.

    After finishing the establishment of the relevant coordinate systems, the 4 steps to deduce the expression form ofAc,tare shown as follows.

    (5)

    Because there are no special restrictions on the way that the camera is installed on the turntable, the matrixA1represents a general three dimensional rigid body transformation. The orthogonal rotation matrixR1has only 3 degrees of freedom, and it can be parameterized with 3 rotation angles around the coordinate axes, which can be denoted as (α,β,γ). The translation vectort1also has 3 degrees of freedom, and it is composed of 3 translational components along the coordinate axes, which can be denoted as (r,s,t). These 6 parameters describes how the camera is installed on the tracking turntable, and therefore they are called the assembly parameters of the camera in this paper[15-16].

    Fig.4 Rotation of OB-xByBzB(θ,φ) around B axis

    (6)

    where

    (7)

    Fig.5 Spatial relationship of two shafts of turntable

    (8)

    where

    (9)

    The two parametersdandεare determined by the structure of the tracking turntable, and therefore they are called the structure parameters of the tracking turntable in this paper.

    Fig.6 Rotation of OA-xAyAzA(θ) around A axis

    (10)

    where

    (11)

    By establishing two intermedia coordinate systemsOA-xAyAzAandOB-xByBzB, the 4 steps of coordinate transformation all become straight for ward. According to the coordinate transformation chain listed in Eqs.(5), (6), (8) and (10), the homogeneous coordinate transformation matrix between the camera coordinate systemOc-xcyczcand the turntable coordinate systemOt-xtytztcan be calculated from the product of the transformation matrix in each step:

    Ac,t=A4A3A2A1.

    (12)

    From the expression forms of the 4 matrixes in Eqs.(5), (7), (9) and (11), it can be seen that the matrixAc,tcontains 10 variables in total. Among them, the two azimuth anglesθandφof theAaxis and theBaxis are directly given by the readings of the two high precision angular encoders installed on the tracking turntable. The rest 8 variables (α,β,γ,r,s,t,d,ε) are the parameters of the tracking turntable. As is mentioned in Steps 1) and 3), they can be divided into two parts: the assembly parameters of the camera (α,β,γ,r,s,t), and the structure parameters of the tracking turntabledandε. The values of these parameters need to be determined by the calibration of the tracking turntable. This will be introduced in the Section 2.

    To sum up, the measurement principle of the newly designed light pen whole space coordinate measuring system based on a tracking turntable can be summarized as

    (13)

    The matrixAl,cis calculated from the image of the light pen by solving PnP problem, and the matrixAc,tis calculated according to Eq.(12).

    1.3 Tracking equation of turntable

    In the newly designed light pen whole space coordinate measuring system, the camera can track the movement of the light pen automatically. The general principle of the tracking process is to make the light pen always located at the center of the camera’s field of view approximately. As long as there is at least one control point can be seen by the camera, the rotation angles of the two shafts of the turntable can be calculated according to the current positon of the light pen in the image.

    (14)

    The objective of the tracking process is to make the centroid of the control points coincide with the principal point of the camera’s imaging plane. According to the pinhole imaging model, the rotation angles of the turntable’s two shafts, which can be denoted as Δθand Δφ, can be calculated by

    (15)

    where (cu,cv,fu,fv) are the intrinsic parameters of the camera as shown in Eq.(3).

    (16)

    Once the movement of the light pen makes the condition in Eq.(16) satisfied, the two shafts need to rotate by the angles calculated according to Eq.(15), so that the moving light pen is always in the center of the camera’s field of view approximately. Therefore, the equation in Eq.(15) is called the tracking equation of the turntable.

    2 System calibration

    In order to calculate the matrixAc,tin Eq.(13) accurately, the parameters (α,β,γ,r,s,t,d,ε) of the tracking turntable must be calibrated as well. In this paper, a method is proposed to calibrate the parameters of the tracking turntable.

    2.1 Calibration procedure

    The overall idea of the proposed calibration method is utilizing the camera to capture the image of the stationary light pen from different view angles. The calibration procedure is as follows:

    1) Fix the pedestal of the tracking turntable on a stable platform and make sure that its position remains unchanged during the whole calibration process.

    2) Set the azimuth angles of theAaxis and theBaxis of the turntable both to zero, i.e.θ=0°,φ=0°.

    3) As shown in Fig.7(a), locate the light pen at a proper distance from the camera and make it at the center of the camera’s field of view. Fix the current position of the light pen and make sure its position remains unchanged during the whole calibration process as well.

    Fig.7 Schematic diagram of calibration procedure of turntable’s parameters

    4) As shown in Fig.7(b), let theAaxis deflect to the right and theBaxis pitch down respectively to change the view angle of the camera, so that the light pen becomes at the top left corner of the camera’s field of view. Record the current azimuth angles of theAaxis and theBaxis.

    5) As shown in Figs.7(c)-(e), similar to the operation in Step 4), change the azimuth angles of theAaxis and theBaxis sequentially, so that the light pen becomes at the other three corners of the camera’s field of view, respectively. Record the azimuth angles of theAaxis and theBaxis that correspond to these 3 situations.

    6) Through the operations in Steps 4) and 5), the variation ranges of the two azimuth angles during the whole calibration process are determined, which can be denoted asθ∈[θb,θe] andφ∈[φb,φe]. According to these two ranges, a two dimensional subdivision calibration grid is generated by segmenting both the two ranges into several equal intervals.

    7) Let the two shafts of the tracking turntable rotate automatically to scan line by line according to the generated calibration grid, so that the azimuth angles of theAaxis and theBaxis can traverse every node of the calibration grid. At each node of the grid, use the camera to capture the image of the stationary light pen from different view angles.

    8) With the captured light pen images, a series of nonlinear equations can be established according to the invariable coordinate transformation between the light pen and the pedestal of the turntable. The parameters of the tracking turntable can be solved from the established equation system with nonlinear least squares method.

    2.2 Equation establishment

    For the reason that the positions of the light pen and the pedestal of the tracking turntable both remain unchanged during the whole calibration process, the homogeneous coordinate transformation matrix between the light pen coordinate systemOl-xlylzland the turntable coordinate systemOt-xtytztis constant, which can be denoted asAl,t. According to this invariant relation, a series of nonlinear equations can be established.

    Suppose that the numbers of intervals in the azimuth angles of theAaxis and theBaxis arepandq, respectively, then the azimuth angles of the two shafts at each node of the calibration grid can be denoted as (θi,φk), (i=1,2,…,p;k=1,2,…,q). They can be calculated by

    (17)

    (18)

    wherei1,i2=1,2,…,p;k1,k2=1,2,…,q;i1≠i2;k1≠k2.

    From Eq.(18), 12 nonlinear equations which take the 8 parameters of the tracking turntable as unknown variables can be deduced. The concrete expression forms of these equations are quite complicated. For convenience, they are expressed in a simplified form as

    (19)

    whereh=1,2,…,12. When all the nodes of the calibration grid are taken into consideration at the same time, an overdetermined nonlinear equation system can be established. The parameters of the tracking turntable can be solved from this equation system by nonlinear least squares algorithm.

    2.3 Solving algorithm

    1) Give the initial solution vectorX(0). Set the threshold value of the convergence conditionσ1andσ2. Set the iteration index tol=0.

    2) In thel-th iteration, the current solution vector isX(l). Calculate the value vectorG(X(l)) and the Jacobian matrixJ(X(l)) according to this solution vector.

    3) Solve the Gauss-Newton equation to obtain the iteration direction at the current solution vectorX(l):Z(l)=-(J(X(l))TJ(X(l)))-1J(X(l))G(X(l)).

    4) Calculate the step length along iteration directionZ(k)by using one dimensional line search, which can be formulized as the optimization problem:a(l)=argmin‖G(X(l)+λZ(l))‖2. This problem can be solved with the rational extremum method.

    5) With the calculated iteration direction and step length, update the solution vector according to the iterative criterion:X(l+1)=X(l)+a(l)Z(l).

    6) Check the threshold condition and judge whether the iteration process converges: if ‖d(l)Z(l)‖≤σ1‖X(l)‖ or ‖G(X(l+1))-G(X(l))‖≤σ2‖G(X(l))‖, the iteration terminates and the solution vectorX(l+1)is the final solution; otherwise go back to Step 2) and continue the iteration with the updated solution vectorX(l+1). At the same time, update the iteration indexl=l+1.

    3 Results and discussion

    To verify the precision of the proposed calibration method, the calibration experiments are repeated 10 times. The calibration results of the parameters of the turntable are listed in Table 1. In the camera’s assembly parameters, the repeatability of the three rotation angles (α,β,γ) is no more than 0.02°, and the repeatability of the three translational components (r,s,t) is no more than 0.05 mm. In the structure parameters of the turntable, the repeatability of the distancedbetween the two shafts can reach 0.023 8 mm, and the repeatability of the non-perpendicularityεbetween the two shafts can reach 0.001 19°. From Table 1, it can be concluded that the proposed method to calibrate the parameters of the tracking turntable has good repeatability.

    To verify the measuring accuracy of the light pen whole space coordinate measuring system based on a tracking turntable, the system is used to measure a standard gauge block whose length is 1 000 mm. In order to highlight the superiority of the newly designed whole space coordinate measuring system, the gauge block is placed in many different directions around the turntable, and located at several different distances from the turntable. Because the camera can rotate driven by the turntable to track the movement of the light pen, the new coordinate measuring system can complete all these measurements through setting up it only once. This is impossible for the classical LPCMS. Table 2 shows the results of the measuring experiments at 7 different distances, and the measurements are repeated 10 times at each distance. It can be seen that the accuracy of the newly designed system can reach 0.25 mm within 10 m, which can satisfy the accuracy requirement of most in situ measurement of large workpieces such as ship hull, plane shell, etc.

    Table 1 Results of 10 calibration experiments (PAR, AVE, STD, REP are abbreviations of parameter,

    From the measuring results in Table 2, it can be seen that the worst accuracy appears when the measuring distance is 1.5 m. The reason is that the major error source at this distance is the rotation of the tracking turntable. The smaller the measuring distance is, the more widely the two shafts of the turntable rotate. Although calibrated, the parameters of the turntable inevitably involve some errors. The measuring errors will magnify rapidly with the enlargement of the rotation angles. As the measuring distance increases, the major error source gradually becomes the decrease of the spatial resolution of the camera. According to the perspective projection model, in the premise that the image resolution is unchanged, the spatial resolution decreases with the enlargement of distance. As a result, for the 1 000 mm standard gauge block, the measuring accuracy first decreases rapidly and then gradually increases with the enlargement of the measuring distance within a certain distance range.

    Table 2 Measuring results of 1 000 mm gauge block (DIS, AVE, RAN, STD are abbreviations of distance,

    4 Conclusions

    In this paper, a new kind of light pen whole space coordinate measuring system based on a tracking turntable is designed. By introducing a turntable to the classical LPCMS, the camera can track the movement of the light pen during the whole measuring process. Benefiting from this, the measuring range is no longer confined to the view field of the camera. As a result, the newly designed system can accomplish whole space coordinate measuring, and therefore can be applied to large scale on-site measurement of workpieces such as ship hull, train body, plane shell, etc. To guarantee the accuracy of the new system, a method to calibrate the parameters of the tracking turntable is also proposed. After calibration, the real-time homogeneous coordinate transformation between the camera and the pedestal of the turntable can be accurately determined according to the turntable’s parameters and the azimuth angles of the two shafts. After experimental verification, the whole space coordinate measuring accuracy can reach 0.25 mm within 10 m. It can be concluded that the newly designed light pen whole space coordinate measuring system significantly expands the measuring range of classical LPCMS without losing too much accuracy. The research significance of the new system is extending the application scope of classical LPCMS and providing a brand new solution to large scale on-site coordinate measurement.

    免费不卡黄色视频| 一级片免费观看大全| 欧美激情 高清一区二区三区| 欧美黑人欧美精品刺激| 美女大奶头黄色视频| 国产极品粉嫩免费观看在线| 国产精品一区二区在线不卡| 777久久人妻少妇嫩草av网站| 悠悠久久av| 日韩电影二区| 狂野欧美激情性bbbbbb| 亚洲中文字幕日韩| 男女之事视频高清在线观看| 精品卡一卡二卡四卡免费| 男人操女人黄网站| 久久久精品免费免费高清| 男人添女人高潮全过程视频| 亚洲精品一卡2卡三卡4卡5卡 | 午夜激情久久久久久久| 69精品国产乱码久久久| 欧美精品一区二区大全| e午夜精品久久久久久久| 国产精品九九99| 国产精品麻豆人妻色哟哟久久| 建设人人有责人人尽责人人享有的| 欧美乱码精品一区二区三区| 国产99久久九九免费精品| 欧美日韩精品网址| 国产不卡av网站在线观看| cao死你这个sao货| 夫妻午夜视频| 国产视频一区二区在线看| 午夜福利乱码中文字幕| 欧美激情高清一区二区三区| 十分钟在线观看高清视频www| 国产99久久九九免费精品| 黄色 视频免费看| 亚洲第一青青草原| 免费观看a级毛片全部| 色视频在线一区二区三区| 黑人巨大精品欧美一区二区mp4| 精品视频人人做人人爽| 一级,二级,三级黄色视频| 多毛熟女@视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品久久久久久精品古装| 久久久精品免费免费高清| 男女床上黄色一级片免费看| 国产成人影院久久av| 欧美另类一区| 麻豆乱淫一区二区| 十八禁网站免费在线| 狠狠狠狠99中文字幕| 亚洲免费av在线视频| netflix在线观看网站| 久久人人97超碰香蕉20202| www.999成人在线观看| 中文欧美无线码| 久久国产精品人妻蜜桃| 黄色 视频免费看| 欧美日本中文国产一区发布| 久久性视频一级片| 日日爽夜夜爽网站| 精品亚洲成a人片在线观看| 天天添夜夜摸| 成人国语在线视频| 丝袜脚勾引网站| 日韩精品免费视频一区二区三区| 欧美精品人与动牲交sv欧美| 正在播放国产对白刺激| 侵犯人妻中文字幕一二三四区| av欧美777| 男女免费视频国产| 国产精品 欧美亚洲| 老司机在亚洲福利影院| 久久亚洲国产成人精品v| 久久久久国内视频| 亚洲精品久久午夜乱码| 极品少妇高潮喷水抽搐| 亚洲七黄色美女视频| 成人三级做爰电影| 久久久久网色| 亚洲精品一二三| 日韩大码丰满熟妇| 精品一区二区三区av网在线观看 | 午夜视频精品福利| 黄色视频,在线免费观看| 99热网站在线观看| 成年人午夜在线观看视频| 国产主播在线观看一区二区| 久久精品亚洲av国产电影网| 国产成人免费无遮挡视频| 天天躁夜夜躁狠狠躁躁| 90打野战视频偷拍视频| 久久久久视频综合| 少妇 在线观看| 他把我摸到了高潮在线观看 | 中国国产av一级| 国产亚洲av高清不卡| 18禁国产床啪视频网站| 免费观看a级毛片全部| 在线观看免费高清a一片| 韩国精品一区二区三区| 国产男女超爽视频在线观看| 午夜成年电影在线免费观看| 欧美日韩黄片免| 男女床上黄色一级片免费看| 桃花免费在线播放| 狂野欧美激情性bbbbbb| 国产亚洲欧美在线一区二区| 大片电影免费在线观看免费| 日韩一区二区三区影片| 亚洲国产日韩一区二区| 国产精品 国内视频| av片东京热男人的天堂| 99久久综合免费| 9191精品国产免费久久| 国产精品久久久久成人av| 免费在线观看完整版高清| 新久久久久国产一级毛片| 女人爽到高潮嗷嗷叫在线视频| 99九九在线精品视频| 亚洲欧洲精品一区二区精品久久久| 热99re8久久精品国产| 欧美黄色淫秽网站| 免费在线观看黄色视频的| 久久久国产一区二区| 9色porny在线观看| 欧美日韩亚洲综合一区二区三区_| 天堂俺去俺来也www色官网| 夜夜骑夜夜射夜夜干| 欧美 亚洲 国产 日韩一| 精品第一国产精品| 99国产精品99久久久久| 精品少妇内射三级| 手机成人av网站| 18禁观看日本| 天天添夜夜摸| 国产激情久久老熟女| 麻豆乱淫一区二区| 久久人妻福利社区极品人妻图片| 欧美精品av麻豆av| 午夜免费成人在线视频| 色精品久久人妻99蜜桃| 91精品国产国语对白视频| 精品亚洲成国产av| 视频在线观看一区二区三区| 在线观看免费高清a一片| 午夜影院在线不卡| 桃红色精品国产亚洲av| 黄色a级毛片大全视频| 成年av动漫网址| 精品免费久久久久久久清纯 | 国产黄频视频在线观看| 99精国产麻豆久久婷婷| 天堂中文最新版在线下载| 麻豆av在线久日| 国产成人精品无人区| 日韩人妻精品一区2区三区| 一进一出抽搐动态| av天堂久久9| 777久久人妻少妇嫩草av网站| 啦啦啦免费观看视频1| 国产在线视频一区二区| 午夜福利免费观看在线| 精品福利永久在线观看| av福利片在线| 精品乱码久久久久久99久播| 久久精品国产亚洲av香蕉五月 | 午夜成年电影在线免费观看| 亚洲欧美清纯卡通| 999久久久国产精品视频| kizo精华| 久久久久久久久久久久大奶| 亚洲第一欧美日韩一区二区三区 | 日韩一卡2卡3卡4卡2021年| 成人18禁高潮啪啪吃奶动态图| 国产亚洲精品一区二区www | 国产在线视频一区二区| 欧美日韩亚洲综合一区二区三区_| 国产精品久久久久久人妻精品电影 | 亚洲精品美女久久av网站| 色精品久久人妻99蜜桃| 免费在线观看日本一区| 青春草视频在线免费观看| 日本av免费视频播放| 热re99久久国产66热| 午夜福利乱码中文字幕| e午夜精品久久久久久久| 亚洲精品成人av观看孕妇| www.av在线官网国产| 狂野欧美激情性xxxx| 国产老妇伦熟女老妇高清| 久久av网站| bbb黄色大片| 丝袜脚勾引网站| 在线天堂中文资源库| 免费在线观看黄色视频的| 久久久久网色| 美女高潮到喷水免费观看| 亚洲精品国产一区二区精华液| 12—13女人毛片做爰片一| 亚洲九九香蕉| 国产成人免费无遮挡视频| 美女中出高潮动态图| 亚洲色图 男人天堂 中文字幕| 国产免费一区二区三区四区乱码| 人人妻人人爽人人添夜夜欢视频| 午夜激情久久久久久久| 男女午夜视频在线观看| 亚洲色图 男人天堂 中文字幕| 免费在线观看视频国产中文字幕亚洲 | 日本wwww免费看| 王馨瑶露胸无遮挡在线观看| 人人妻人人爽人人添夜夜欢视频| 少妇的丰满在线观看| 一区二区三区激情视频| 男女下面插进去视频免费观看| 男女国产视频网站| 午夜激情av网站| 老汉色∧v一级毛片| 天天躁狠狠躁夜夜躁狠狠躁| 香蕉丝袜av| 亚洲欧美一区二区三区黑人| 亚洲熟女精品中文字幕| 国产极品粉嫩免费观看在线| 真人做人爱边吃奶动态| 999久久久国产精品视频| 他把我摸到了高潮在线观看 | 99久久人妻综合| 精品视频人人做人人爽| 中国美女看黄片| 日韩制服骚丝袜av| 首页视频小说图片口味搜索| 国产成+人综合+亚洲专区| 国产精品久久久av美女十八| 国产精品自产拍在线观看55亚洲 | 岛国在线观看网站| 国产欧美亚洲国产| 中文字幕最新亚洲高清| 国产淫语在线视频| 国产麻豆69| 国产精品一二三区在线看| 香蕉丝袜av| 一级毛片女人18水好多| 精品人妻1区二区| 超色免费av| 欧美乱码精品一区二区三区| 女人久久www免费人成看片| 国产精品免费大片| 免费少妇av软件| 人人妻人人澡人人看| 亚洲自偷自拍图片 自拍| 亚洲男人天堂网一区| 亚洲综合色网址| 纯流量卡能插随身wifi吗| 黑丝袜美女国产一区| 淫妇啪啪啪对白视频 | 女人精品久久久久毛片| 每晚都被弄得嗷嗷叫到高潮| 中文欧美无线码| 亚洲专区字幕在线| 老司机靠b影院| 亚洲国产欧美网| 午夜福利视频在线观看免费| 久久精品成人免费网站| 青青草视频在线视频观看| 国产成人啪精品午夜网站| 国产成人av激情在线播放| 亚洲国产精品成人久久小说| 狠狠精品人妻久久久久久综合| 欧美日本中文国产一区发布| 91国产中文字幕| 十分钟在线观看高清视频www| 又黄又粗又硬又大视频| 一个人免费在线观看的高清视频 | 国产欧美日韩一区二区三 | 日韩熟女老妇一区二区性免费视频| 亚洲视频免费观看视频| 欧美变态另类bdsm刘玥| 欧美xxⅹ黑人| 大香蕉久久成人网| 久久精品国产a三级三级三级| 搡老岳熟女国产| 国产成人精品无人区| 欧美另类亚洲清纯唯美| 女性被躁到高潮视频| 捣出白浆h1v1| 精品一区在线观看国产| 精品一品国产午夜福利视频| 亚洲免费av在线视频| 午夜视频精品福利| 男女之事视频高清在线观看| 我要看黄色一级片免费的| 日韩三级视频一区二区三区| 国产一区二区 视频在线| 亚洲第一av免费看| 三级毛片av免费| 黄网站色视频无遮挡免费观看| 亚洲黑人精品在线| 99久久99久久久精品蜜桃| 国产又色又爽无遮挡免| 日本一区二区免费在线视频| 老司机亚洲免费影院| 国产在线观看jvid| 久久久精品免费免费高清| 别揉我奶头~嗯~啊~动态视频 | 十八禁高潮呻吟视频| 99国产精品免费福利视频| 免费人妻精品一区二区三区视频| 黄片播放在线免费| 亚洲精品一区蜜桃| 欧美在线一区亚洲| 国产高清videossex| 美女福利国产在线| 在线av久久热| 中文字幕精品免费在线观看视频| 国产亚洲av片在线观看秒播厂| 亚洲国产欧美一区二区综合| 免费在线观看日本一区| 无遮挡黄片免费观看| 母亲3免费完整高清在线观看| 黑人猛操日本美女一级片| 91字幕亚洲| 香蕉丝袜av| 久久久久精品国产欧美久久久 | 亚洲成人免费电影在线观看| 美女高潮喷水抽搐中文字幕| 18禁黄网站禁片午夜丰满| av电影中文网址| 亚洲专区中文字幕在线| 人人妻人人澡人人爽人人夜夜| 菩萨蛮人人尽说江南好唐韦庄| 国产日韩一区二区三区精品不卡| 成年av动漫网址| 天天操日日干夜夜撸| av国产精品久久久久影院| 国产日韩欧美在线精品| 超碰成人久久| 久久国产精品男人的天堂亚洲| 亚洲美女黄色视频免费看| www.熟女人妻精品国产| 午夜日韩欧美国产| 国产精品麻豆人妻色哟哟久久| 国产精品.久久久| 日本av免费视频播放| 欧美日韩亚洲综合一区二区三区_| 国产麻豆69| 黄色视频在线播放观看不卡| 亚洲精品久久成人aⅴ小说| 热99re8久久精品国产| 国产在视频线精品| 亚洲精品一二三| 操出白浆在线播放| 国产精品久久久久成人av| 成年人免费黄色播放视频| 女人精品久久久久毛片| 丝袜脚勾引网站| 性高湖久久久久久久久免费观看| 可以免费在线观看a视频的电影网站| 亚洲中文字幕日韩| 19禁男女啪啪无遮挡网站| 人妻久久中文字幕网| 少妇粗大呻吟视频| 交换朋友夫妻互换小说| 一区二区三区四区激情视频| 亚洲第一av免费看| 桃红色精品国产亚洲av| xxxhd国产人妻xxx| 免费黄频网站在线观看国产| 欧美精品一区二区免费开放| 两性午夜刺激爽爽歪歪视频在线观看 | 国产区一区二久久| 国产视频一区二区在线看| 黄色视频,在线免费观看| 亚洲av成人一区二区三| 狠狠婷婷综合久久久久久88av| 国精品久久久久久国模美| 在线天堂中文资源库| 日日夜夜操网爽| 亚洲中文av在线| 中文字幕高清在线视频| 黄色视频不卡| 午夜91福利影院| 老司机影院毛片| 亚洲av成人一区二区三| 精品免费久久久久久久清纯 | netflix在线观看网站| 久久国产精品大桥未久av| 丰满少妇做爰视频| 满18在线观看网站| 侵犯人妻中文字幕一二三四区| 丝袜人妻中文字幕| 少妇 在线观看| 欧美亚洲 丝袜 人妻 在线| 久热爱精品视频在线9| 久久久久久久久久久久大奶| 天堂俺去俺来也www色官网| 香蕉国产在线看| 丰满饥渴人妻一区二区三| 一区在线观看完整版| 在线亚洲精品国产二区图片欧美| 午夜老司机福利片| 欧美av亚洲av综合av国产av| 国产亚洲午夜精品一区二区久久| 亚洲国产精品999| 国产精品亚洲av一区麻豆| 国产在线观看jvid| 一级毛片电影观看| 亚洲精品一二三| 人人妻人人添人人爽欧美一区卜| 国产欧美亚洲国产| 搡老熟女国产l中国老女人| 大陆偷拍与自拍| 高清欧美精品videossex| 国产成人精品无人区| 少妇 在线观看| 久久九九热精品免费| 在线亚洲精品国产二区图片欧美| 亚洲免费av在线视频| 欧美在线黄色| 亚洲av电影在线观看一区二区三区| 91精品三级在线观看| 国产免费视频播放在线视频| 国产亚洲精品久久久久5区| 美女国产高潮福利片在线看| 在线观看免费日韩欧美大片| 啦啦啦啦在线视频资源| 一区二区av电影网| 日韩视频在线欧美| 午夜激情久久久久久久| 亚洲,欧美精品.| 中文字幕av电影在线播放| av有码第一页| 99国产精品一区二区三区| 一级片免费观看大全| 国产亚洲一区二区精品| 黄色片一级片一级黄色片| 狂野欧美激情性bbbbbb| 麻豆国产av国片精品| 动漫黄色视频在线观看| av福利片在线| av在线老鸭窝| 丰满人妻熟妇乱又伦精品不卡| 精品乱码久久久久久99久播| 中文字幕高清在线视频| 久久99一区二区三区| 国产视频一区二区在线看| 亚洲一卡2卡3卡4卡5卡精品中文| 丁香六月天网| 丝袜喷水一区| 精品久久蜜臀av无| 国产免费av片在线观看野外av| 91老司机精品| 精品国产一区二区三区久久久樱花| 无限看片的www在线观看| 亚洲一码二码三码区别大吗| 精品国产超薄肉色丝袜足j| 精品熟女少妇八av免费久了| 一个人免费看片子| 中文字幕最新亚洲高清| 久久久国产精品麻豆| 99国产精品一区二区三区| 免费看十八禁软件| 男人舔女人的私密视频| 可以免费在线观看a视频的电影网站| 久久精品aⅴ一区二区三区四区| 青青草视频在线视频观看| 91成年电影在线观看| 丰满少妇做爰视频| 中亚洲国语对白在线视频| 一级a爱视频在线免费观看| 国产精品成人在线| 97在线人人人人妻| 久久久精品94久久精品| 亚洲精品美女久久av网站| 国产亚洲精品一区二区www | 国产成人精品久久二区二区91| 亚洲欧美激情在线| 女警被强在线播放| 免费少妇av软件| 久久国产精品人妻蜜桃| 日韩视频一区二区在线观看| 亚洲精品一卡2卡三卡4卡5卡 | 久久人人爽人人片av| av福利片在线| 日韩 欧美 亚洲 中文字幕| 亚洲情色 制服丝袜| 捣出白浆h1v1| 啦啦啦中文免费视频观看日本| 亚洲第一av免费看| 大香蕉久久网| 久久国产精品人妻蜜桃| av视频免费观看在线观看| 欧美黄色淫秽网站| 午夜精品国产一区二区电影| 69精品国产乱码久久久| 啦啦啦啦在线视频资源| 国产黄频视频在线观看| 亚洲精品第二区| 久久人妻熟女aⅴ| 午夜福利影视在线免费观看| 黄色视频,在线免费观看| 99九九在线精品视频| 亚洲av电影在线进入| 午夜日韩欧美国产| 日韩欧美免费精品| 两性午夜刺激爽爽歪歪视频在线观看 | 别揉我奶头~嗯~啊~动态视频 | 9色porny在线观看| 国产有黄有色有爽视频| 丝袜脚勾引网站| av福利片在线| 成年动漫av网址| 日日爽夜夜爽网站| 啪啪无遮挡十八禁网站| 97精品久久久久久久久久精品| 国产极品粉嫩免费观看在线| 亚洲精品久久成人aⅴ小说| 亚洲色图综合在线观看| 狠狠婷婷综合久久久久久88av| 亚洲三区欧美一区| 国产一区二区在线观看av| 青春草视频在线免费观看| 免费在线观看日本一区| 国产成人精品无人区| 飞空精品影院首页| 高清av免费在线| 啦啦啦 在线观看视频| 成年人午夜在线观看视频| 69精品国产乱码久久久| 日韩精品免费视频一区二区三区| 国产成人免费无遮挡视频| 国产精品久久久久久人妻精品电影 | 亚洲精品一区蜜桃| av天堂久久9| 国产精品av久久久久免费| 国产精品1区2区在线观看. | 欧美老熟妇乱子伦牲交| av电影中文网址| 国产黄频视频在线观看| 亚洲国产毛片av蜜桃av| 久久久久久久精品精品| 欧美在线一区亚洲| 好男人电影高清在线观看| 欧美另类一区| 国产精品自产拍在线观看55亚洲 | a 毛片基地| 99久久国产精品久久久| 他把我摸到了高潮在线观看 | 免费黄频网站在线观看国产| 50天的宝宝边吃奶边哭怎么回事| 高清黄色对白视频在线免费看| 久久久久久久大尺度免费视频| 可以免费在线观看a视频的电影网站| 99re6热这里在线精品视频| 老司机亚洲免费影院| 最黄视频免费看| 久久久精品免费免费高清| 亚洲av片天天在线观看| 大香蕉久久成人网| 国产在线观看jvid| 亚洲精品av麻豆狂野| a级毛片在线看网站| 国产在线一区二区三区精| 久久影院123| 黄片小视频在线播放| 91老司机精品| 日韩中文字幕欧美一区二区| 在线天堂中文资源库| 精品第一国产精品| 精品人妻熟女毛片av久久网站| 一级毛片精品| av有码第一页| 久久99热这里只频精品6学生| 亚洲视频免费观看视频| 黄片播放在线免费| 两个人免费观看高清视频| 日韩三级视频一区二区三区| 一区二区三区四区激情视频| 久久女婷五月综合色啪小说| 日日爽夜夜爽网站| 精品免费久久久久久久清纯 | 啪啪无遮挡十八禁网站| 国产精品1区2区在线观看. | 欧美人与性动交α欧美软件| 亚洲av日韩在线播放| 亚洲精品一卡2卡三卡4卡5卡 | 脱女人内裤的视频| 久久精品人人爽人人爽视色| 国产老妇伦熟女老妇高清| 国产成人一区二区三区免费视频网站| 宅男免费午夜| 亚洲自偷自拍图片 自拍| 欧美日韩福利视频一区二区| 国产淫语在线视频| 欧美日韩国产mv在线观看视频| 亚洲色图 男人天堂 中文字幕| 午夜福利乱码中文字幕| 咕卡用的链子| 久9热在线精品视频| 国产主播在线观看一区二区| 女人爽到高潮嗷嗷叫在线视频| 久9热在线精品视频| 国产精品国产av在线观看| 久久午夜综合久久蜜桃| 亚洲一码二码三码区别大吗| 亚洲伊人色综图| 人人妻人人添人人爽欧美一区卜| 天天操日日干夜夜撸| 男女午夜视频在线观看| 免费在线观看日本一区| 这个男人来自地球电影免费观看| 国产又爽黄色视频| 亚洲av成人一区二区三|