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

    A Study on Hovering Control of Small Aerial Robot by Sensing Existing Floor Features

    2020-08-05 09:40:04ChinthakaPremachandraDangNgocHoangThanhTomotakaKimuraandHiroharuKawanaka
    IEEE/CAA Journal of Automatica Sinica 2020年4期

    Chinthaka Premachandra,,Dang Ngoc Hoang Thanh,Tomotaka Kimura,,and Hiroharu Kawanaka,

    Abstract—Since precise self-position estimation is required for autonomous flight of aerial robots, there has been some studies on self-position estimation of indoor aerial robots. In this study, we tackle the self-position estimation problem by mounting a small downward-facing camera on the chassis of an aerial robot. We obtain robot position by sensing the features on the indoor floor.In this work, we used the vertex points(tile corners) where four tiles on a typical tiled floor connected,as an existing feature of the floor. Furthermore,a small lightweight microcontroller is mounted on the robot to perform image processing for the on board camera. A lightweight image processing algorithm is developed.So, the real-time image processing could be performed by the microcontroller alone which leads to conduct on-board real time tile corner detection. Furthermore,same microcontroller performs control value calculation for flight commanding.The flight commands are implemented based on the detected tile corner information. The above mentioned all devices are mounted on an actual machine,and the effectiveness of the system was investigated.

    I.Introduction

    SMALL aerial robots are capable of motion such as vertical takeoff and landing in narrow spaces, horizontal movements like forward, backward,and hovering.Studies can be performed into autonomous flight considering applications such as room surveillance and transporting cargo by using these motions.Examples of research into autonomous flight include hovering[1],[2], movement to a target location[3],obstacle avoidance[4],[5],and object detection[6]–[8].Furthermore,when performing autonomous flight indoors,since the flight area is limited,only a small robot can be used.Therefore,it is necessary to use the hardware that is mounted on the aerial robot small and lightweight.

    To implement the autonomous flight of an aerial robot,it is essential for the robot to perform self-position estimation.The GPS cannot be used in indoor environments.Indoor selfposition estimation of general robots has been studied using distance sensors such as ultrasonic sensors and infrared sensors[9]–[13],laser range finders[14],and cameras[15]–[21].Although ultrasonic sensors and infrared sensors need only light processing and many lightweight devices are available,but they suffer from low precision.Laser range finders are usually large,and are not suitable for mounting on a small aerial robot.In contrast,cameras can visually capture the surrounding environment,and can implement precision measurement through image processing.

    In this study,we study autonomous flight control of an aerial robot by estimating its self-position using an on-board camera.There are some studies in the literature, regarding self-position estimation and autonomous flight control of aerial robots by using on-board cameras.Those studies include the method of obtaining position information by capturing images of markers embedded in the ground surface,using downward-facing camera mounted on the robot[15]–[18].In those methods,the camera mounted on the robot captures images of the markers.Here,the size and shape of the markers are already known.The robot orientation and position information are obtained from the geometrical relationship between the markers.However, those methods require advance installation of markers inside the flight area.In this case,robot can conduct operation only within the marker existing area.

    Some studies in the literature also perform self-position estimation and robot control by capturing natural features from images[19],[20].Those methods capture natural features those exist on the floor and calculate the optical flow to obtain the robot position information.Since markers are not used,there are no limitations on the flight region,and measurement can also be performed in unknown locations.The common methods of calculating the optical flow include the block matching method and gradient method,and these allow high-precision measurements as long as the distances are short.However,the cumulative error becomes larger as the measurement distance becomes longer.In addition,it is difficult to implement them using light weight hardware like microcontrollers due to complexity of the algorithms.

    Therefore,in order to address these problems,in this study we investigate recognizing tile corners using a camera mounted on an aerial robot to obtain position information and perform autonomous flight.Tiled floors are widely used in indoors.This removes the installing markers in advance because the tile corners can be used as an existing feature.Furthermore, by memorizing the initial position of the tile corner,the error can be reset if the robot returns to the initial position during measurement.Furthermore,since the tile corners are arranged with equal spacing,errors of moving distance can be corrected by using the constant distance between tile corners.

    In this study,the robot position information is obtained based on these tile corners,and flight commands are sent to the robot based on the obtained position information.The final goal of this study is to perform autonomous flight indoors.Since the flight area is limited due to the size of the room and so when performing autonomous flight indoors,it is necessary to use a small robot.Due to the limitations of robot body surface and the weight that robot can carry,a small lightweight camera is mounted on its body surface.Furthermore,small lightweight microcontroller(Raspberry pi3)is mounted on the robot. A lightweight image processing algorithm is developed to perform real-time image processing by this on-board microcontroller alone.In addition,same microcontroller is used to perform flight control values decisions,following detected self-position information by image processing.The entire system is mounted on an actual machine,and the effectiveness of the system was investigated.

    The results showed that it is possible to obtain robot position information by sensing tile corners.Furthermore,we confirmed the possibility of real-time autonomous flight control based on the detected robot positions.

    II.System of Sensing Aerial Robot Position by Using Existing Floor Features

    Fig.1 shows the system configuration of this study.A downward-facing monocular camera is mounted on the robot.In the position information calculation unit,floor images are captured by the camera,and they are sent to the microcontroller mounted on the robot.The microcontroller performs image processing to obtain the robot position information.The control signals are generated by the same microcontroller based on the obtained position information.Ultimately,the robot is controlled by those control signals which is sending to on-board flight controller.All the parts of the system will be detailed in next sub-sections.

    A. Position Information Calculation Unit

    The robot position information is obtained by recognizing the tile corners in images of the floor,captured by the camera.When the position capturing is started,the position of the tile corner first recognized is used as the initial position.Furthermore,the position information is sent to the control value calculation unit.

    Fig.1.System configuration.

    B.Control Value Calculation Unit

    The control values are determined based on the position information from the position information calculation unit.Then determined control values are sent to the flight controller of the aerial robot.In this paper,we target flight control to achieve hovering of the robot.The detailed information on control value calculation is presented in section five.

    III.Estimation of Aerial Robot Position Information

    A. Position Estimation Stream

    Fig.2 shows the procedure of the robot position estimation.Fig.3 shows the floor to be captured,whereas Fig.4 shows the points on the floor to be recognized.In this study,tile corners are extracted by analyzing the directional component of the floor pattern those appear on the images.To achieve this tile corner detection,the general corner detection methods such Harris corner detection and its related methods[22],[23]cannot be applied since they detect all the corner points in the image and the extraction of this specific point(tile corner)from them is very complicated.In this paper,the tile corners are recognized by image processing and then the robot position information is obtained.The methods of each process is present in next sections.In addition, the horizontal translation distance of the robot is measured by comparing the coordinates of the tile corners between consecutive frames.

    Fig.2.Position information estimation procedure.

    Fig.3.Captured image.

    Fig.4.Tile corner.

    Since we select a floor that is more complicated than a regular tiled floor,the system is expected to be able to handle floors with various tile shapes.Furthermore,we develop an algorithm with a small amount of calculations in order to perform the real time image processing on a microcontroller.

    IV.Tile Corner Recognition

    Fig.5 show s the sequence of the image processing during tile corner recognition.In Fig.5,the 80×80 pixel original floor image(a)is converted into a brightness layer image(b),and pattern enhancement and binarization are performed (c).The floor pattern direction components are extracted from the binarized image(d),and image angle correction is performed if needed (e).The direction components are extracted after the rotation correction (f),and noise removal is then performed on the image to obtain a low noise image(g).Finally,processing to recognize the tile corners is performed (h).

    Fig.5.Sequence of image processing.

    A.Conversion to Brightness Layer Image

    In order to perform the calculations more efficiently in the image processing, the color image obtained from the camera is converted into a brightness layer image.A brightness layer image resembles a regular grayscale image,and is converted into an image that represents only the luminosity information from the color RGB image.The brightness layer image colored with 32 levels of luminance as shown in Fig.6.The luminosity valuesYin the brightness layer image are determined following(1).Here,R,G,and B denote red,green,and blue component value of each pixel in the color image,respectively.The color image is converted into a brightness layer image by storing the found luminosity valueYin each of the corresponding pixels.The larger luminosity valueYis denoted by the brighter pixels.

    Fig.6.Brightness layer image.

    B. Pattern Enhancement and Binarization

    In this study,tile corners are recognized based on the regularity of the floor pattern.Because of this,it is important to extract the direction component of the floor pattern.However,the floor pattern is not clear in the brightness layer image,so a pattern enhancement filter is used to make it clearer.In this study,the patterns in the horizontal and vertical directions are enhanced.We design a simple filter to remove unnecessary information and extract only the target information by aligning then×nfilter shown in Fig.7.This can be expressed in the form of an equation as defined by(2),by takingI(x,y)as the image,G(i,j)as then×n(n=3)filter(where 0≤i≤n–1,0 ≤j≤n–1),andH(x,y)as the value after filter processing.TheH(x,y)value for each pixel in the image is calculated.Finally,each pixel is replaced by correspondingH(x,y).We experimentally designed this filter to carry on this work.This filtering can be done with some filters in literature[24]–[28].But, they are comparatively time consuming with compared to above mentioned filter.

    Fig.7.Pattern enhancement filter.

    After the filter processing,an image with the floor pattern enhanced as shown in Fig.5(c).Here binarized image is obtained by setting a threshold.C. Extraction of Direction Components After Pattern Enhancement Filtering and Binarization

    To recognize whether the angle of the binarized image Fig.5(c)needs to be corrected,direction component extraction (DCE)is performed.

    First,a pattern matching process is performed on binarized image Fig.5(c)to achieve directional components.In this process,original pixel pattern is compared with the defined three patterns.Fig.8 illustrates the defined patterns for direction component extraction.They can extractX-direction,Y-direction,and no direction component in the image.This process is performed to entire image to achieve DCE image(Fig.5(d)or 5(f)).We describe about image shown in Fig.5(f)in next sub-section.In Figs.5(d)and 5(f),above three directional components are plotted using three different colors.Fig.9 shows an example of aY-directional component.Based on this processing,whether the image needs to be rotated is judged as detailed in next sub-section.

    Fig.8.Pattern matching.

    Fig.9.Example of direction component pattern.

    D. Angle Correction by Image Rotation

    We take the camera angle shown in Fig.10 as 0°angled image.In reality,the robot will rotate while aerial,and the camera image will be rotated from 0°in the plus or minus directions as shown in Fig.11.Since the recognizedXandYdirection components become smaller due to the rotation,it becomes harder to obtain the floor patterns.As a result, the accuracy of recognizing the tile corners decreases.In order to solve this problem, rotation of the image is calculated,and then angle correction is performed before the tile corner determination.

    Fig.10.At angle0°.

    Fig.11.Rotated image.

    An image,after the pattern enhancement and binarization is considered(e.g.,Fig.5(c)).We determine the total number ofX-direction components in the image asXtotaland the total number ofY-direction components asYtotal.The total of them(Dsum)is calculated following(3).IfDsumexceeds a defined threshold value,we judge that image rotation correction is not needed.IfDsumis less than the threshold value,the following process is conducted to achieve angle correction by rotating the image.

    As shown in Fig.12, the value ofDsumgradually decreases until the camera is rotated by ±45°.Correction is performed to have a 0°angled image following this variation, rotating the image by angler°.The images obtained by rotating in to plus and minus directions are different each other.Furthermore,the direction in which the correction needs to be performed,cannot be judged based on the initially obtainedDsumalone.Therefore,above mentioned pattern matching for directional component extraction is applied to the images,those obtained by rotating the binarized image in to both directions as shown in Fig.13.The value ofDsum+obtained by rotating the image in to plus direction.In the other words,Dsum+is sum of calculatedDsumvalues for each image during the plus direction rotation.The value ofDsum–is obtained by rotating the image in to minus direction in same manner.They are compared based on the rule in(4)to determine the rotation direction.After determining the direction,The image which shows maximumDsumvalue in that direction is selected as angle corrected image.Fig.5(e)is an example image for angle corrected image.In addition,its corresponding DCE image is also generated during this process.following the Fig.5,the image shown in Fig.5(f)is an example for this DCE image.

    After this process,noise removal as described in Section III.Eis performed on DCE image.Furthermore,if it is judged that angle correction is not needed,noise removal is performed on the DCE image without conducting the angle correction.following the Fig.5,step(e)and(f)are not conducted in this case.

    Fig.12.Relationship between camera angle and D sum.

    Fig.13.Deciding the rotation direction.

    E. Noise Removal

    In order to detect the tile corners, noise removal is performed on DCE image since the image that has been processed up to this contains a lot of noise (Figs.5(d) or 5(f)). Fig.14 illustrates the noise removal process.The DCE image is divided into 5×5 pixel segments.The number ofXdirection components(Xi,j)and the number ofYdirection components(Yi,j)in each segment is compared.The direction component that accounts for the highest proportion of each segment is used as the direction componentDk,lfor the entire segment following(5).After this, the segment is filled with defined color for each direction as illustrated in Fig.14.Once it is confirmed that this processing is finished throughout the entire area,the noise removal finishes.The Fig.15 illustrates an image after this noise removal,and this allows the tile corners to be recognized with relatively high precision.

    Fig.14. Noise removal.

    F. Deciding the Tile Corner Coordinates

    The coordinates of the tile corner are decided based on the noise-removed direction component image.Fig.15 shows an example of noise-removed direction component image.Here,the Brown denotesXdirection component while Green denotesYdirection component.Fig.16 shows the tile corner coordinate searching method.The designated row in the image is scanned from the left to search for the coordinates of the boundary between theXdirection component andYdirection component.If a boundary with three consecutiveYdirection components are found they are labeled as①,②,and③.The same process repeat within all designated rows and labeled as illustrated in Fig.16.The line which connects maximum boundaries is selected as border on vertical direction.The border on the horizontal direction is also determined in same manner.Finally,the cross point of these two lines is selected as tile corner.

    Fig.15.Direction component image.

    Fig.16.Searching for the coordinates of the tile corner.

    G. Horizontal Translation Distance Measurement

    The horizontal translation distance of aerial robot is measured by tracking the tile corners.We take the coordinate of the tile corner recognized in then-th frameFnas(xn,yn).The translation distanceLn,n+1betweenFnandFn+1is then given by (6).

    Therefore,the translation distanceL1,nfrom the first frame to then-th frame is given by (7).

    V.Hovering Control

    A. Hovering Control

    There are many studies can be found in the literature regarding the controlling different type of robots[29]–[31].In this paper,we decided to use pulse position modulation (PPM)[32]to transfer the shift between tile corner and aerial robot to achieve hovering control of aerial robot,since it is less time consuming.Furthermore,PPM can easily be implemented on a small microcontroller.The PPM method is a modulation method that converts a signal into time-based pulse position(phase difference)relative to a period.The information is determined by the position of the pulse.The PPM signal used in this study is illustrated in Fig.17.It is composed of 8 channels of pulses and dummy pulse.The dummy pulse is for adjusting the period of the overall signal.This time,hovering control is performed with PID control and PD control following the PPM value, to hover aerial robot over the detected tile corner.Here,PID control is used to achieve the altitude(zdirectional control)and PD control is used to achieve horizontal directional control(xandydirectional control)[33].

    Fig.17.PPM signal.

    We use the PID for thez-axis instead of a PD controller is that the battery voltage/performance decreases with flight time.A full battery needs upwards lift equal to the gravitational force of the aerial robot.A battery at 50%however needs a throttle command in the same motor speeds respectively same upwards force.This loss of “power”,upwards lift over time or acceleration in thez-axis can be compensated by an integrator in the controller.Here the I-part of the PID can constantly increase the throttle command to hold the same height/set point(hover).For thex-andy-axis,we do not need to compensate the loss of battery voltage,since it does not affect thex-yspace positions.The aerial robot would theoretically just stay at the same horizontal position,as long as the roll and pitch angles are 0°.However,we also chose PD-controllers, because the aerial robot only is operated inside a certain range inx-yspace near to a tile corner.

    As mentioned in Section II-B, the control signal is generated by on-board microcontroller based on the position information.Generally,a detected tile corner does not move more than several pixels within a two consecutive frames,during hovering control.Therefore,in this paper,if the detected tile corner in current frame moves irregularly,more than a pre-defined threshold, we ignore the detected position.In this case,hovering control is performed following the detected tile corner in previous frame.

    In this study,we keep both manual control and autonomous control of the aerial robot.Webring the aerial robot to defined altitude by manual control. After that,when the robot detects a tile corner,it starts the autonomous hovering control keeping the current altitude by PID control(zdirectional control).As mentioned above,xandydirectional robot shift from the detected tile corner is modified by PD control.To achieve an effective hovering control, the camera pose has to be horizontal.We have solved this matter applying a pose keeping hardware.

    B. Deciding Flight Commands

    In this study,the center of the 80×80 pixel image at(40,40)is taken as a reference coordinate(Fig.18).When the recognized tile corner is shifted from the reference coordinate,in other words,when the robot is moved from directly above the tile corner,control is performed to remain above the reference coordinate to achieve hovering.Since thex-axis direction in the image corresponds to the forward and backward motion of the robot,if thexcoordinate changes from the reference coordinate,the CH1 pulse width in Fig.17 is changed according toxcoordinate shift value and the aileron is controlled following PD control to return the robot to reference coordinate.Furthermore,since they-axis direction in the image corresponds to the left and right motion of the robot,if theycoordinate changes from the reference coordinate,the CH2 pulse width in Fig.17 is changed according to theycoordinate shift and the elevator is controlled following PD control to return the robot to reference coordinate[33].In this paper,the parameters of PD control is determined by sub-experiments.

    Fig.18.Reference coordinate(The center of the80×80 pixel image).

    VI.Experiments

    A. Investigation of Tile Corner Recognition and Horizontal Translation Distance Measurement

    To investigate the accuracy of the system,we performed accuracy verification experiments on the tile corner recognition and horizontal translation distance measurement by using the on-board downward-facing monocular camera and Raspberry pi3 microcontroller.During the flight,the rotation around thezaxis is assumed to be zero.The camera pixel configuration is set to 80×80 pixels.

    B. Experimental Method of Tile Corner Recognition

    The images of the floor were captured with on-board camera and the accuracy of the tile corner recognition and horizontal translation distance measurement were investigated.Here,the recognition rate was investigated with the approximate camera elevation at 1000 mm.At this camera elevation,there is only one tile corner appear in an image.Furthermore,the machine with camera was approximately moved 200mm horizontally,and the actual horizontal translation distance and measurement result values were compared.The experiment environment used in buildings of the Katsushika Campus,Tokyo University of Science.Fig.19 shows an example image of a floor.In this work,in order to support different kinds of tiled floors,a more complicated tile shape than typical tiled floors were chosen.

    Fig.19.Example of floor image used in experiments.

    C. Results of Tile Corner Detection

    Fig.20 shows the experimental results of tile corner detection.When the tile corner is recognized,a blue point is displayed in the image.Through the experimental results,we confirm that recognition was possible even when the camera was rotated.Furthermore,the recognition success rate as a result of capturing 30 000 frames while flight was 97.6%.This result confirms that the tile corner could be recognized with high accuracy.Furthermore,Fig.21 shows the main stages of the image processing.According to this,it confirms that the subsequent processing was performed after judging whether angle correction was needed from the image.

    Fig.20.Experimental results of tile corner detection.

    Fig.21.Experimental results of tile corner detection with all processing stages.

    Furthermore, the average error of horizontal translation measurement when the machine elevation was approximately at 1000 mm is 2.7%.This shows that high accuracy measurement of horizontal movement was possible.This tile detection could be conducted in real-time since the average computational time was18 ms per frame.

    D. Hovering Control Procedure Following Tile Corner

    To confirm whether the system built in this work can perform robot control,control verification experiments were performed on Raspberry pi3 microcontroller.

    The(x,y)coordinate in the center of the image at (40,40)is treated as a reference coordinate.As shown in Fig.22,it is assumed that the robot will be controlled to return the camera coordinates to the position of the reference coordinate if the robot is moved,and the PPM signal is output to perform this kind of control.In this work, we performed an investigation assuming that the camera coordinates are returned to the position of the reference coordinate when the robot is moved forward,backward,left,and right.For example,a PPM signal is output that controls the elevators to return to(40,40)when(x,y)changes to(20,40)or(60,40)and that controls the ailerons to return to(40,40) when (x,y)changes to(40,20)or(40,60).Fig.23 shows the coordinate relationship between the robot and the camera image.

    Fig.22.Investigation of control command output.

    Fig.23.Coordinate relationship between the robot and camera image.

    E. Results of PPM Generation

    We analyze the PPM variation during the flight movement.Figs.24 and 25 show a comparison of the PPM signal when the camera is at the reference coordinate with the PPM signal when the camera is at(x,y)=(20,40)and(60,40),and Figs.26 and 27 show a comparison of the PPM signal when the camera is at the reference coordinate and the PPM signal when the camera is at (x,y)= (40,20)and (40,60).Fig.24.Comparison of signals between the reference coordinate and (x,y) =(20,40).

    Fig.25.Comparison of signals between the reference coordinate and (x, y) =(60,40).

    Looking at the brown waveform in Fig.24,we can see that when the camera moves from the reference coordinate to(20,40),a PPM signal output becomes shorter following the required minus shift from the original.Similarly,looking at the green waveform in Fig.25,we can see that when the camera moves from the reference coordinate to(60,40),a PPM signal output becomes wider,following the required plus shift.

    Furthermore,looking at the green waveform in Fig.26, we can see that when the camera moves from the reference coordinate to(40,20),a PPM signal output becomes wider following the required plus shift from the original.Similarly,looking at the green waveform in Fig.27,we can see that when the camera moves from the reference coordinate to(40,60),a PPM signal output become smaller according to the required minus shift.

    Fig.26.Comparison of signals between the reference coordinate and (x, y)= (40,20).

    Fig.27.Comparison of signals between the reference coordinate and (x, y)= (40,60).

    This experiment confirms that PPM signals are output based on the obtained position information.

    F. Hovering With a Real Machine

    Fig.28 illustrates the prototype aerial machine prepared for flight experiments.It mainly consists of general parts of a quadcopter.However,to carry on this work,we installed Pixhawk4 flight controller, microcontroller(Raspberry pi3),and small type camera which can be connected to the microcontroller.Tile corner detection is conducted by the microcontroller by processing image from the camera.Pixhawk4 flight controller is connected with microcontroller to communicate each other to achieve data transform.Through the different experiments,we could confirm that hovering control is possible by sensing tile corners on the ground.

    Fig.28.Prepared aerial robot for experiments.

    VII.Conclusion

    In this research,we developed a lightweight image processing algorithm such that image processing can be performed using only a microcontroller with the aim of obtaining the position information of a small aerial robot based on tile corners as an existing floor feature.We then investigated autonomous flight of the robot by sensing the tile corner positions.As a result of the investigation,it was found that the position information for a robot can be obtained by sensing the tile corner using a downward-facing camera mounted on the robot.We also implemented flight commands based on the obtained position information,and confirmed the possibility of autonomous hovering control.

    色综合站精品国产| 美女扒开内裤让男人捅视频| 搞女人的毛片| 韩国精品一区二区三区| 女生性感内裤真人,穿戴方法视频| 一边摸一边抽搐一进一小说| 中文字幕最新亚洲高清| 亚洲三区欧美一区| 亚洲成人国产一区在线观看| 国产精华一区二区三区| 91老司机精品| 精品国产亚洲在线| 一区在线观看完整版| 可以免费在线观看a视频的电影网站| 日日摸夜夜添夜夜添小说| 亚洲国产高清在线一区二区三 | 日韩精品青青久久久久久| 日韩欧美三级三区| 18禁裸乳无遮挡免费网站照片 | 黑人欧美特级aaaaaa片| 亚洲中文字幕日韩| 亚洲男人天堂网一区| 欧美成人一区二区免费高清观看 | 老汉色av国产亚洲站长工具| 亚洲人成电影观看| 亚洲天堂国产精品一区在线| 夜夜看夜夜爽夜夜摸| 色综合亚洲欧美另类图片| 日韩av在线大香蕉| 亚洲国产中文字幕在线视频| 又紧又爽又黄一区二区| 免费在线观看完整版高清| 免费看十八禁软件| 国产在线精品亚洲第一网站| 老司机靠b影院| 国产免费av片在线观看野外av| 国产精品久久久久久人妻精品电影| 久久草成人影院| 女人爽到高潮嗷嗷叫在线视频| 久久久久亚洲av毛片大全| 琪琪午夜伦伦电影理论片6080| 99在线视频只有这里精品首页| 精品一区二区三区av网在线观看| 欧美一区二区精品小视频在线| 日韩三级视频一区二区三区| 极品人妻少妇av视频| 99国产精品99久久久久| 久久人人精品亚洲av| 亚洲狠狠婷婷综合久久图片| 国产成人欧美在线观看| 国产乱人伦免费视频| 国产不卡一卡二| 亚洲aⅴ乱码一区二区在线播放 | 亚洲avbb在线观看| 黑人巨大精品欧美一区二区mp4| 丰满的人妻完整版| 如日韩欧美国产精品一区二区三区| 最新美女视频免费是黄的| 国产精品一区二区精品视频观看| 免费在线观看完整版高清| 国产精华一区二区三区| 看黄色毛片网站| 亚洲精品在线观看二区| 国产成人精品无人区| 午夜福利成人在线免费观看| 国产视频一区二区在线看| 久热这里只有精品99| 啪啪无遮挡十八禁网站| 亚洲在线自拍视频| 国产精品影院久久| 久久人妻熟女aⅴ| 国产精品亚洲av一区麻豆| 欧美乱色亚洲激情| 欧美黑人精品巨大| 后天国语完整版免费观看| 欧美一级a爱片免费观看看 | 黄色 视频免费看| 老汉色av国产亚洲站长工具| 欧美丝袜亚洲另类 | www.www免费av| 欧美国产日韩亚洲一区| 脱女人内裤的视频| 99精品在免费线老司机午夜| 18禁美女被吸乳视频| 欧美成人性av电影在线观看| 午夜精品久久久久久毛片777| 日本在线视频免费播放| 久久精品亚洲熟妇少妇任你| 村上凉子中文字幕在线| 国产一区二区激情短视频| a在线观看视频网站| 久久伊人香网站| 午夜日韩欧美国产| 亚洲精品在线美女| 一区福利在线观看| 欧美日韩亚洲综合一区二区三区_| 大香蕉久久成人网| 成年女人毛片免费观看观看9| 亚洲在线自拍视频| 老司机深夜福利视频在线观看| 18禁观看日本| 午夜免费观看网址| 免费搜索国产男女视频| 成年女人毛片免费观看观看9| 久久人妻av系列| 免费无遮挡裸体视频| 国产在线观看jvid| 欧美日韩亚洲综合一区二区三区_| 国产精华一区二区三区| 久久草成人影院| 精品国产国语对白av| 亚洲中文字幕一区二区三区有码在线看 | 两性午夜刺激爽爽歪歪视频在线观看 | 日日干狠狠操夜夜爽| 一级作爱视频免费观看| 色尼玛亚洲综合影院| 久久久久久久久中文| 中文亚洲av片在线观看爽| 国产午夜福利久久久久久| 啦啦啦免费观看视频1| 国产精品久久久av美女十八| 女生性感内裤真人,穿戴方法视频| 国内精品久久久久久久电影| www日本在线高清视频| 他把我摸到了高潮在线观看| 免费在线观看黄色视频的| 热99re8久久精品国产| 亚洲午夜理论影院| 欧美午夜高清在线| 99久久精品国产亚洲精品| 亚洲最大成人中文| 99国产精品免费福利视频| 非洲黑人性xxxx精品又粗又长| 国产伦一二天堂av在线观看| 亚洲av美国av| 亚洲人成伊人成综合网2020| 一级黄色大片毛片| 91精品三级在线观看| 国产国语露脸激情在线看| 精品久久久精品久久久| 国产成年人精品一区二区| 国产成人av教育| 成人永久免费在线观看视频| 欧美激情 高清一区二区三区| 午夜精品久久久久久毛片777| 久久精品国产综合久久久| 美女免费视频网站| 午夜精品久久久久久毛片777| 香蕉久久夜色| 国产私拍福利视频在线观看| 日韩一卡2卡3卡4卡2021年| 国产精品亚洲美女久久久| 国产精品香港三级国产av潘金莲| 久99久视频精品免费| 香蕉久久夜色| 宅男免费午夜| 日本vs欧美在线观看视频| 久久久久国产精品人妻aⅴ院| 国产一卡二卡三卡精品| 亚洲一区高清亚洲精品| 色哟哟哟哟哟哟| 叶爱在线成人免费视频播放| 国产日韩一区二区三区精品不卡| 亚洲美女黄片视频| 两性夫妻黄色片| 咕卡用的链子| 久久精品国产清高在天天线| 国产单亲对白刺激| 亚洲三区欧美一区| а√天堂www在线а√下载| 91大片在线观看| 亚洲av成人不卡在线观看播放网| 中国美女看黄片| 丁香六月欧美| 日韩精品免费视频一区二区三区| 国产麻豆成人av免费视频| 欧美绝顶高潮抽搐喷水| 欧美激情 高清一区二区三区| 精品无人区乱码1区二区| 亚洲一区中文字幕在线| 午夜福利免费观看在线| 9色porny在线观看| 欧美在线黄色| 好男人在线观看高清免费视频 | 免费看a级黄色片| 久久久久久久久久久久大奶| 亚洲精品av麻豆狂野| 在线视频色国产色| 69精品国产乱码久久久| 久热爱精品视频在线9| 久久久水蜜桃国产精品网| 精品久久蜜臀av无| 午夜精品国产一区二区电影| 18美女黄网站色大片免费观看| 日日夜夜操网爽| 免费少妇av软件| 国产精华一区二区三区| 成年版毛片免费区| 亚洲狠狠婷婷综合久久图片| 啦啦啦观看免费观看视频高清 | 最新美女视频免费是黄的| 亚洲伊人色综图| 欧美成人一区二区免费高清观看 | 久久天躁狠狠躁夜夜2o2o| 丁香六月欧美| 精品国产美女av久久久久小说| 黄色片一级片一级黄色片| 成人国产综合亚洲| 久久中文字幕人妻熟女| 窝窝影院91人妻| 久久国产精品男人的天堂亚洲| 成人手机av| 咕卡用的链子| 亚洲熟妇中文字幕五十中出| 国产激情久久老熟女| 国产精品综合久久久久久久免费 | 好看av亚洲va欧美ⅴa在| 精品国产美女av久久久久小说| 非洲黑人性xxxx精品又粗又长| 又紧又爽又黄一区二区| 日本a在线网址| 精品国产乱子伦一区二区三区| 在线视频色国产色| 国产麻豆69| 在线观看66精品国产| 国产精品精品国产色婷婷| xxx96com| 美女午夜性视频免费| 久久久国产成人免费| 手机成人av网站| 欧美激情极品国产一区二区三区| 色av中文字幕| www.自偷自拍.com| 高清毛片免费观看视频网站| 岛国在线观看网站| 国产精品日韩av在线免费观看 | 亚洲精品国产精品久久久不卡| 美女扒开内裤让男人捅视频| 亚洲欧美日韩无卡精品| 国产麻豆69| 日韩中文字幕欧美一区二区| 日本 欧美在线| 在线观看舔阴道视频| 宅男免费午夜| 亚洲精品国产一区二区精华液| 久久人人97超碰香蕉20202| 久久精品91无色码中文字幕| 国产熟女xx| 亚洲久久久国产精品| 老司机午夜福利在线观看视频| 色哟哟哟哟哟哟| 亚洲国产精品成人综合色| 看片在线看免费视频| 黄色女人牲交| 成年女人毛片免费观看观看9| 午夜福利高清视频| 亚洲色图 男人天堂 中文字幕| 亚洲av五月六月丁香网| 亚洲欧美一区二区三区黑人| 国产亚洲精品av在线| 日韩欧美三级三区| 正在播放国产对白刺激| 精品国产一区二区三区四区第35| av片东京热男人的天堂| 97碰自拍视频| 一级作爱视频免费观看| 亚洲国产高清在线一区二区三 | 麻豆国产av国片精品| 777久久人妻少妇嫩草av网站| 久久中文字幕人妻熟女| 亚洲午夜理论影院| 久久精品人人爽人人爽视色| 欧美最黄视频在线播放免费| 免费在线观看完整版高清| 熟女少妇亚洲综合色aaa.| 女人高潮潮喷娇喘18禁视频| 国产极品粉嫩免费观看在线| 国产高清视频在线播放一区| 黄片大片在线免费观看| 一级a爱视频在线免费观看| 叶爱在线成人免费视频播放| 黑人操中国人逼视频| 久久精品影院6| 天堂影院成人在线观看| 久热这里只有精品99| 亚洲九九香蕉| 色播在线永久视频| 99在线人妻在线中文字幕| svipshipincom国产片| 亚洲精品一区av在线观看| 制服诱惑二区| 亚洲电影在线观看av| 久久久久久久久久久久大奶| 激情在线观看视频在线高清| 搡老妇女老女人老熟妇| 午夜福利影视在线免费观看| 成在线人永久免费视频| 一区二区三区国产精品乱码| 久久香蕉国产精品| 两个人看的免费小视频| 午夜免费观看网址| 午夜久久久久精精品| 国产一区在线观看成人免费| 在线av久久热| 欧美日本亚洲视频在线播放| av片东京热男人的天堂| 男女午夜视频在线观看| 最新在线观看一区二区三区| 午夜免费成人在线视频| 老司机午夜十八禁免费视频| www.自偷自拍.com| 别揉我奶头~嗯~啊~动态视频| 黑人巨大精品欧美一区二区mp4| 久久久久久久精品吃奶| 老熟妇乱子伦视频在线观看| 国产精品久久久久久人妻精品电影| 成人永久免费在线观看视频| 国产成人精品在线电影| 色综合婷婷激情| 亚洲精品国产精品久久久不卡| 精品久久久久久久久久免费视频| 亚洲,欧美精品.| 免费久久久久久久精品成人欧美视频| 午夜福利欧美成人| 午夜福利成人在线免费观看| avwww免费| 亚洲成国产人片在线观看| 国产精品影院久久| 免费不卡黄色视频| 日韩欧美国产一区二区入口| 日韩一卡2卡3卡4卡2021年| 精品电影一区二区在线| 操出白浆在线播放| 天天躁夜夜躁狠狠躁躁| 国产成人一区二区三区免费视频网站| 日韩免费av在线播放| 亚洲成av人片免费观看| 色哟哟哟哟哟哟| 国产精华一区二区三区| 国产午夜福利久久久久久| 亚洲人成电影观看| 精品国产亚洲在线| 好男人在线观看高清免费视频 | 久久精品影院6| 91老司机精品| 亚洲熟妇熟女久久| 国产激情欧美一区二区| 国产精品九九99| 怎么达到女性高潮| 国产一区二区在线av高清观看| 夜夜躁狠狠躁天天躁| 88av欧美| 亚洲 欧美一区二区三区| 99精品欧美一区二区三区四区| 久久热在线av| 亚洲,欧美精品.| 国产日韩一区二区三区精品不卡| 亚洲成av片中文字幕在线观看| 天堂影院成人在线观看| 国产成年人精品一区二区| 两个人看的免费小视频| 老汉色av国产亚洲站长工具| 国产成人免费无遮挡视频| 国产一区二区激情短视频| 国产成人免费无遮挡视频| 国产一区二区激情短视频| 日日爽夜夜爽网站| 成人免费观看视频高清| 亚洲国产日韩欧美精品在线观看 | 中出人妻视频一区二区| 男人舔女人下体高潮全视频| 色尼玛亚洲综合影院| 午夜福利一区二区在线看| 欧美国产日韩亚洲一区| 99热只有精品国产| 一级毛片女人18水好多| 久久久久久大精品| 精品久久久精品久久久| 亚洲av熟女| 岛国在线观看网站| 1024视频免费在线观看| 亚洲人成网站在线播放欧美日韩| 桃红色精品国产亚洲av| 亚洲精品国产区一区二| 国产精品美女特级片免费视频播放器 | 欧美成人一区二区免费高清观看 | 日韩一卡2卡3卡4卡2021年| av中文乱码字幕在线| 操美女的视频在线观看| 美女大奶头视频| 成人三级做爰电影| 他把我摸到了高潮在线观看| 国内精品久久久久精免费| 日韩三级视频一区二区三区| netflix在线观看网站| 国产一区二区三区在线臀色熟女| 又黄又粗又硬又大视频| 满18在线观看网站| 精品少妇一区二区三区视频日本电影| 大香蕉久久成人网| 18禁观看日本| 欧美黄色淫秽网站| 美国免费a级毛片| 亚洲精品在线观看二区| 精品国产一区二区三区四区第35| 亚洲成av片中文字幕在线观看| 国产成人av教育| 99国产精品一区二区三区| 岛国视频午夜一区免费看| 亚洲精品国产色婷婷电影| 黄片小视频在线播放| 黑人欧美特级aaaaaa片| 国产精品日韩av在线免费观看 | 亚洲五月天丁香| 一边摸一边做爽爽视频免费| 亚洲一卡2卡3卡4卡5卡精品中文| 91字幕亚洲| 免费在线观看亚洲国产| 国产一区在线观看成人免费| 中文字幕最新亚洲高清| 19禁男女啪啪无遮挡网站| 亚洲av片天天在线观看| 日本免费a在线| 手机成人av网站| 一边摸一边抽搐一进一小说| 免费在线观看影片大全网站| 精品无人区乱码1区二区| 亚洲欧美日韩无卡精品| 在线播放国产精品三级| 国产欧美日韩精品亚洲av| 一区二区三区国产精品乱码| 亚洲av第一区精品v没综合| 成年版毛片免费区| 两性夫妻黄色片| 亚洲全国av大片| e午夜精品久久久久久久| 美女免费视频网站| 制服人妻中文乱码| 变态另类丝袜制服| 国产免费男女视频| 亚洲国产欧美一区二区综合| 身体一侧抽搐| 亚洲免费av在线视频| 巨乳人妻的诱惑在线观看| 人妻久久中文字幕网| 国产欧美日韩一区二区三| 国内精品久久久久精免费| 国产成人精品在线电影| 国产精品国产高清国产av| 国产成人影院久久av| 女同久久另类99精品国产91| 给我免费播放毛片高清在线观看| 91字幕亚洲| 可以免费在线观看a视频的电影网站| 久热这里只有精品99| 久久久久久久久免费视频了| 91成年电影在线观看| 欧美 亚洲 国产 日韩一| av网站免费在线观看视频| 亚洲精品在线美女| 99久久精品国产亚洲精品| 男人舔女人下体高潮全视频| 欧美一级a爱片免费观看看 | 国产一区二区在线av高清观看| 香蕉久久夜色| a级毛片在线看网站| 999精品在线视频| 免费在线观看亚洲国产| 国产精品香港三级国产av潘金莲| 50天的宝宝边吃奶边哭怎么回事| 国产精品自产拍在线观看55亚洲| 国产精品一区二区精品视频观看| 99精品久久久久人妻精品| 亚洲美女黄片视频| 亚洲精品国产色婷婷电影| 国产午夜精品久久久久久| 久久精品亚洲精品国产色婷小说| www.www免费av| 亚洲 国产 在线| 自拍欧美九色日韩亚洲蝌蚪91| 日韩欧美一区二区三区在线观看| 精品国内亚洲2022精品成人| 中文字幕另类日韩欧美亚洲嫩草| 欧美成人午夜精品| 免费在线观看亚洲国产| 黄色女人牲交| 看黄色毛片网站| 性色av乱码一区二区三区2| 激情在线观看视频在线高清| 国产熟女午夜一区二区三区| 在线视频色国产色| xxx96com| 日韩欧美免费精品| 成人特级黄色片久久久久久久| 三级毛片av免费| 国产不卡一卡二| 欧美一级a爱片免费观看看 | 久久久久久大精品| 亚洲色图av天堂| 别揉我奶头~嗯~啊~动态视频| 欧美av亚洲av综合av国产av| xxx96com| 少妇裸体淫交视频免费看高清 | 久久人人爽av亚洲精品天堂| 国产av精品麻豆| 亚洲欧美日韩另类电影网站| x7x7x7水蜜桃| 丰满的人妻完整版| 精品卡一卡二卡四卡免费| 啪啪无遮挡十八禁网站| 午夜精品国产一区二区电影| 欧美日韩亚洲国产一区二区在线观看| 天堂√8在线中文| 又黄又粗又硬又大视频| 给我免费播放毛片高清在线观看| 亚洲成人精品中文字幕电影| 国产在线精品亚洲第一网站| 制服丝袜大香蕉在线| 亚洲熟女毛片儿| 亚洲欧美激情综合另类| 国产私拍福利视频在线观看| 好看av亚洲va欧美ⅴa在| 99riav亚洲国产免费| 乱人伦中国视频| 丰满人妻熟妇乱又伦精品不卡| 久久久久久久午夜电影| 欧美精品亚洲一区二区| 无人区码免费观看不卡| 少妇被粗大的猛进出69影院| 精品久久久久久久人妻蜜臀av | 搞女人的毛片| tocl精华| 国产又色又爽无遮挡免费看| 国产精品电影一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 精品国产超薄肉色丝袜足j| 夜夜夜夜夜久久久久| 91麻豆精品激情在线观看国产| 亚洲国产日韩欧美精品在线观看 | 欧美一级a爱片免费观看看 | 叶爱在线成人免费视频播放| 侵犯人妻中文字幕一二三四区| 亚洲人成伊人成综合网2020| 91麻豆精品激情在线观看国产| 9191精品国产免费久久| 一本久久中文字幕| 无遮挡黄片免费观看| 久久人妻av系列| 午夜亚洲福利在线播放| 免费在线观看影片大全网站| 精品久久久久久成人av| 色综合站精品国产| 精品不卡国产一区二区三区| 最新在线观看一区二区三区| 国产精品一区二区三区四区久久 | 欧美日韩黄片免| 99久久综合精品五月天人人| 免费高清视频大片| 成年版毛片免费区| 热99re8久久精品国产| 丰满的人妻完整版| 香蕉久久夜色| 亚洲色图 男人天堂 中文字幕| 午夜福利欧美成人| 99国产精品99久久久久| 日韩高清综合在线| 欧美av亚洲av综合av国产av| 亚洲色图综合在线观看| 母亲3免费完整高清在线观看| 亚洲全国av大片| 国产高清videossex| 国产精品自产拍在线观看55亚洲| 精品国产亚洲在线| 成人亚洲精品av一区二区| 久久久久久国产a免费观看| 97人妻天天添夜夜摸| 性少妇av在线| 亚洲成av片中文字幕在线观看| 亚洲一区中文字幕在线| 美女 人体艺术 gogo| 免费av毛片视频| 窝窝影院91人妻| 欧美日韩乱码在线| 97人妻天天添夜夜摸| 男女床上黄色一级片免费看| 69精品国产乱码久久久| 亚洲精品一区av在线观看| 欧美一级毛片孕妇| 免费少妇av软件| 国产高清videossex| 一二三四社区在线视频社区8| 国产精品电影一区二区三区| 亚洲精品中文字幕一二三四区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲五月婷婷丁香| 亚洲电影在线观看av| 午夜精品久久久久久毛片777| 国产亚洲精品久久久久久毛片| 99在线视频只有这里精品首页| 亚洲少妇的诱惑av| 欧美一级毛片孕妇| 一区福利在线观看| 黑人欧美特级aaaaaa片| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品永久免费网站| 长腿黑丝高跟| 亚洲一区二区三区不卡视频| 一级片免费观看大全| 激情视频va一区二区三区| 亚洲人成电影免费在线| 搡老岳熟女国产| 九色国产91popny在线| 亚洲无线在线观看| 美女国产高潮福利片在线看| 欧美乱妇无乱码|