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

    Detection of broken manhole coverusing improved Hough and image contrast

    2015-03-01 08:07:31ZhangFengyanChenRongbaoLiYangGuoXiucheng

    Zhang Fengyan  Chen Rongbao  Li Yang  Guo Xiucheng

    (1School of Transportation, Southeast University, Nanjing 210096, China)(2School of Electrical Engineering and Automation, Hefei University of Technology, Hefei 230009, China)(3Jiangsu Transportation Institute Co., Ltd., Nanjing 210017, China)(4Department of Civil Engineering, Anhui Communications Vocational and Technical College, Hefei 230051, China)

    ?

    Detection of broken manhole coverusing improved Hough and image contrast

    Zhang Fengyan1,3,4Chen Rongbao2Li Yang2Guo Xiucheng1

    (1School of Transportation, Southeast University, Nanjing 210096, China)(2School of Electrical Engineering and Automation, Hefei University of Technology, Hefei 230009, China)(3Jiangsu Transportation Institute Co., Ltd., Nanjing 210017, China)(4Department of Civil Engineering, Anhui Communications Vocational and Technical College, Hefei 230051, China)

    Abstract:The damage or loss of urban road manhole covers may cause great risk to residents’ lives and property if they cannot be discovered in time. Most existing research recommendations for solving this problem are difficult to implement. This paper proposes an algorithm that combines the improved Hough transform and image comparison to identify the damage or loss of the manhole covers in complicated surface conditions by using existing urban road video images. Focusing on the pre-processed images, the edge contour tracking algorithm is applied to find all of the edges. Then with the improved Hough transformation, color recognition and image matching algorithm, the manhole cover area is found and the change rates of the manhole cover area are calculated. Based on the threshold of the change rates, it can be determined whether there is potential damage or loss in the manhole cover. Compared with the traditional Hough transform, the proposed method can effectively improve the processing speed and reduce invalid sampling and accumulation. Experimental results indicate that the proposed algorithm has the functions of effective positioning and early warning in the conditions of complex background, different perspectives, and different videoing time and conditions, such as when the target is partially covered.

    Key words:manhole cover; edge tracking; improved Hough transform; shape detection; image contrast

    Received 2015-01-20.

    Biographies:Zhang Fengyan(1972—),doctor, male, zhangleike@126.com; Guo Xiucheng (corresponding author), male, doctor, professor, seuguo@163.com.

    Foundation item:The Natural Science Fundation of Education Department of Anhui Province (No.KJ2012B051).

    Citation:Zhang Fengyan, Chen Rongbao, Li Yang, et al. Detection of broken manhole cover using improved Hough and image contrast[J].Journal of Southeast University (English Edition),2015,31(4):553-558.[doi:10.3969/j.issn.1003-7985.2015.04.021]

    The manhole cover is essential to any city. There are typically a large number of manhole covers in a city which are difficult to manage and monitor. Once a manhole cover is damaged but not dealt with timely, it will bring great hidden dangers to the lives, property, or safety of people. Prasad et al.[1]used the edge curvature to determine the edge based on the traditional Hough transform method, which searches each edge profile in the area and makes a grouping calculation based on the conditions of the edge contour to determine the elliptical shape. However, this method is limited by the standard elliptical contour; therefore, it cannot be used on the shape detection of broken and damaged elliptical manhole covers. Kasemir et al.[2]proposed a detection algorithm for the elliptical shape in a high noise level image. The algorithm combines the improved Hough transform and the genetic algorithm. However, we need to design a novel algorithm that synthesizes the improved Hough transform and image comparison for solving different technical problems; so it is difficult to generalize this algorithm for the detection of the broken manhole covers. Wei et al.[3]developed a detection method on the incomplete elliptical image in strong noise environments based on the iterative randomized Hough transform (IRHT). The IRHT method iteratively applies the randomized Hough transforms, constantly adjusting the iterative parameter and repeatedly amplifying the target curve by using image space and parameter space. Due to the fact that the core of this method is the iterative computation, it is time and computer memory demanding, and thus its adaptability for moving vehicles is poor in dynamic environments. Much research has been done on the manhole cover, but the research mainly focuses on manufacturing intelligence, by which a timely alarm can be set if the manhole cover is damaged. Although this method can effectively monitor the manhole cover, the implementation process of these measures requires extensive human and material resources, which will be very expense. Since a large number of manhole covers need to be maintained, there will be a great resistance for the application of manufacturing intelligence[4]. In recent years, many cities have established municipal video monitoring networks. Therefore, using existing video surveillance networks, and the image acquisition and processing technology of computer vision, we develop a method to monitor the manhole covers and raise alarm in time when there are missing and broken manhole covers. This can achieve the aim of monitoring manhole covers without any hardware investment and repetitive construction.

    Usually, very complex images can be taken by cameras with a complex background, but people are often interested in a particular region of the image with special characteristics, which is the target area. It is necessary to extract the target area from the original image to further identify and analyze the target. This paper focuses on identifying the manhole cover region with more clear geometric features and computational geometry characteristics. Hence, compared with other spatial invariance features, calculating geometric features is more simple and convenient[5], and with the help of the colour recognition algorithm, it can better meet the requirements of accuracy and real-time. Therefore, we select the geometry features and colors of the manhole cover as the positioning method to contrast the images and detect manhole cover breakage.

    1Overall Scheme

    The manhole cover in the image can be viewed approximately as the distortion of ellipses. The recognition of the manhole cover is to determine the region[4]and position of the manhole cover with edge curves calculated by the image processing algorithm. The Hough transform is not sensitive to image noise, and it easily conducts parallel processing and convenient for handling characteristics when images are covered locally. Therefore, it is more suitable for the elliptical extraction of the manhole cover image with complex backgrounds[5-6]. However, the traditional Hough transform has several limitations:

    1) It requires a large amount of computation. For linear, the algorithm complexity isO(n2); for circle and ellipse, the algorithm complexities areO(n3) andO(n5), respectively.

    2) It needs a large CPU memory.

    3) The parameters extracted by quantization have interval constraints of parameter space, so it is difficult to meet high precision and real-time requirements.

    Therefore, conducting an elliptical traditional Hough transform directly in the image with complicated background is not feasible[7]. In order to overcome the limitations of the traditional Hough transform, this paper first conducts the edge search to find the edges above all the image and stores them in the image list. Then we divide them into each sub image which is transformed by Hough. This can reduce the noise influence to a great extent and reduce invalid image accumulation. Also, the space of the Hough transform and time consumption are reduced to meet the system’s high precision and real-time requirements. Due to the interference from other elliptic objects[8-10], the road condition is very complex. So this paper combines the color recognition and image alignment algorithm to recognize and localize the manhole cover, which can greatly improve the identification accuracy of the manhole cover.

    As shown in Fig.1, the method mainly consists of four parts: pre-processing, contour extraction, target localization and identification, and detecting suspected damage. In order to ensure an accurate extraction of the target region, before positioning, it is necessary to do some pre-processing operations such as brightness and contrast adjustment, smoothing, and denoising and edge detection. Then we search the edge, delete the useless contour, extract the useful contour, and recognize and localize the manhole cover by the improved Hough transform and the color recognition algorithm. Finally, the image matching algorithm is used to calculate the suspected discrimination by the change rate.

    Fig.1 The process of the proposed method

    2Identification and Localization of the Manhole Cover

    Since images are affected by acquisition conditions such as illumination uniformity and the time period, we need to carry out image preprocessing, which mainly includes the adjustment of brightness, equalization and smoothing of the contrast histogram, the denoising method of mathematical morphology and the detection of Canny edge[11]. These preprocessing algorithms are classic approaches, which are not focused on in this paper.

    2.1 Edge tracking

    In order to analyze the shape of the collected image, we need to find edge points of the target and then connect the edge points into one continuous boundary for boundary tracking. As shown in Fig.2, the specific steps are as follows:

    1) Set a discontinuous reference distance “gap”, establish a two-dimensional image arrayI(store two images after binaryzation) and a dynamic two-dimensional array “cur” (coordinate storage of edge point).

    2) Scan a two-dimensional arrayIof an image from top to bottom and from left to right, save the point on the left side ofA0after locating the edge pointA0(x0,y0) in the array “cur”, and setA0(x0,y0) to correspond to the coordinate pointI(x,y)=0 in theIarray, showing that the point has been recorded.

    3) Scan the local window [x-gap,x+gap]×[y-gap,y+gap] centered by coordinates (x,y) in the imageIfrom top to bottom and from left to right, and seek the nearest edge pointA1to point (x,y). IfA1exists, then recordA1for coordinates (x,y) and repeat Step 3), until no edge points can be found. Denote edge pointsAe(xe,ye). It should be noted that in Step 3) the tracking edge contour is carried out only in one direction.

    4)A0(x0,y0) is regarded as the starting point coordinates. Set (x,y)=(x0,y0), repeat Step 3), and the end point coordinate is denoted asAs(xs,ys).

    5) If there arenlines of continuous edge in the images, repeat Steps 2) to 4) and produce a closed edge point set. The corresponding edge point set is {cur(i),i=1,2,…,n}.

    Fig.2 The schematic of edge tracking

    2.2 Elliptic automatic recognition

    Among a plurality of edge point set boundaries calculated by the edge tracking algorithm, there is the elliptic boundary of the manhole cover and other interfering boundaries. To distinguish the target boundary from the boundaries, we need to recognize the elliptical objects. The camera angle is different, so the manhole cover in the image will stretch and distort to a different elliptic according to a certain angle, that is, an oblique ellipse with a certain included angle existing in the direction of the main shaft and theXshaft. The standard equation of the ellipse is

    (1)

    In the actual situation, the axis of the elliptical image of the manhole cover is not parallel to the coordinate axes and can be approximately seen as translation and rotation.

    Also, the major and minor axis of ellipse may not be parallel to the coordinate axes, which can be viewed as the translation and rotation of the coordinates. As shown in Fig.3, thexoyaxis is translated tox′o′y′, and thenx′o′y′is rotated tox″o″y″. The ellipse in the coordinatesx″o″y″ is finished. Supposing that the rotation angle isθ, the elliptic equations in thexoycoordinates can be written as

    (2)

    The expression of the short half axisbis

    (3)

    Fig.3 The schematic of elliptic rotation and translation

    Based on contour tracking, each profile in the list is processed in turn. The steps are as follows:

    1) For each profile in the list, initialize the horizontal and vertical coordinates of the initial and end points from the neighborhood control structure.

    2) Do circulation processing of each pixel point in the current contour neighborhood supposing that the point is the center pointCof the ellipse. Calculate the Euclidean distance of every point in the current profile to pointC, find the maximum distance value, set it as the semi major axis of the elliptic, and preserve the corresponding pointP.

    3) If the maximum distance value found in Step 2) is less than the specified threshold, then switch to Step 2) and treat the next pixel of neighborhood structure, otherwise, calculate the angle between pointPandCconnection andxaxis.

    4) Do circular processing of every pixel in the current profile. According to Eq.(3), calculate the ellipse semi minor axisb, and plus 1 to the value of each pixel in the accumulator of the semi minor axis.

    5) Carry out circular treatment for each element of the accumulator of the semi minor axis. If its value is between the upper and lower bounds of the prescribed threshold, set the currenta,b,C,θ, respectively, as the semi major axis, ellipse semi minor axis, center, rotation angle of coordinate value, to obtain all the parameters of elliptic equations, and add them to the list of ellipses.

    6) Empty the semi minor axis accumulator array. If there are any pixels of the current neighborhood profile not disposed of, turn to Step 2), otherwise, turn to Step 7).

    7) If there are some outlines left in the list, turn to Step 1), otherwise, the algorithm terminates.

    Performing the Hough transform to each scattered contour field and reducing the processing range can improve the processing speed effectively, so the system can satisfy the real-time requirement[12].

    2.3 The image matching algorithm

    Although the improved Hough transform algorithm can recognize ellipses on the surface of the road, the average case is very complex due to many other interfering elliptic objects; therefore, the manhole cover cannot be distinguished from the elliptical interference. To solve this problem, this paper adopts the color recognition and image alignment algorithm after staining the manhole cover at the same time. Besides the improved Hough method, the comprehensive discrimination method combining these two algorithms can improve the accuracy of the recognition. At the same time, the image matching algorithm can distinguish between suspect broken or missing manhole covers by different rates[13].

    For the image matching algorithm,it requires a stored image of the standard pavement of the manhole cover part as a standard reference, and then the image is pretreated. Since the detected pavement image and standard image are different, after pretreatment, a registration of the standard pavement image and the detected pavement image is required. In the process of detecting, compared with the input image and the standard image, if the difference is greater than a certain threshold, the manhole cover is suspected to be broken or missing. As shown in Fig.4, the image in Fig.4 is the manhole cover pavement image.

    Fig.4 The process of inspecting manhole cover’s breakage and loss

    Image registration is a fundamental task in image processing and is used to match two or more images acquired at different time, different sensors, different angles and different shooting conditions[13]. Considering the actual situation of the defective detection system, the difference between the collected detected image and standard image is rigid deformation, so we can use the registration method based on gray level information.

    Assume that the standard reference image isR; the detecting image isS, the size ofRism×n, the size ofSisM×N, as shown in Fig.5. The basic process of the image registration method based on gray level information is as follows: superimpose the reference imageRon the imageSand translate; let the area covered by the reference image on the searching imageSto be the subgraph ofSij; (i,j) is the coordinates of the upper left corner of subgraph in the image on theS. The search range is

    1≤i≤M-m, 1≤j≤N-n

    (4)

    Fig.5 The basic process of image registration method

    The registration process is completed by analyzing the similarity betweenRandSij.

    According to the similarity measure function, the registration methods based on gray level information can be classified into the cross-correlation algorithm, the maximum mutual information registration method etc. The cross-correlation algorithm is used to realize the registration of two images, which is the most basic image registration method based on gray level statistics. It requires the reference image and the matching image with a similar scale and gray information. The cross-correlation algorithm is used to match the selected region of the image on the reference template and to calculate the cross correlation coefficient for each position on the reference image and the matching image. The choice of the method is due to that before calculating the cross correlation coefficient of the reference image and the test image, a rotating operation is needed to ensure that the maximum cross correlation coefficient is achieved when the axes of two plates are parallel. To maintain an unchanged standard plate direction, we only move the test plate, which can reduce error and improve registration accuracy. The registration process is shown in Fig.6.

    1) According to the preset step value, rotate the measured pavement image within a certain angle range (including the forward rotation and the backward rotation). After each rotation, conduct a cross correlation calculation of the selected area of two pictures.

    2) Choose the corresponding rotation angle when the

    Fig.6 The process of image registration by cross-correlation algorithm

    cross-correlation coefficient reaches the maximum; modify the measured pavement image according to the angle and the axial of two images will become parallel.

    3) Select some area on the revised image of the pavement surface, and then calculate the cross-correlation coefficient. This is mainly in order to obtain the corresponding position when the coefficient reaches its maximum.

    4) Using the position of the maximum value, we can deduce the relationship between pixel values of the corresponding points on the two images. Through translation and a cutting operation, complete the registration.

    Conduct the XOR Boolean calculation after matching the detecting manhole cover image with the standard image of the two values, and then create the outline of image change. After some morphological treatment, we can obtain the suspected damaged area of the manhole cover. We find the manhole cover area through the method of Hough and color recognition and the alarm will sound when the manhole cover area is beyond a certain threshold by calculating the suspected damage zone size in its region.

    3Program Verification

    A simulation platform was developed based on the Matlab. The experimental image is obtained by the color video camera. Figs.7 to 9 show the key steps of image processing. In each of Figs.7 to 9, the subfigure (a) is the original image of a manhole cover, subfigure (b) is a manhole cover image after dyeing, subfigure (c) is an image of the manhole cover after damage or loss, subfigure (d) is the preprocessed image of a manhole cover after damage or loss, subfigure (e) is an alarming manhole cover image of the suspected defect. Fig.7 shows a manhole cover with a small area of the damage; Fig.8 presents a manhole cover with a large area of damage; and Fig.9 shows the pictures of a completely lost manhole cover. It is seen that this system can accurately sound the alarm for manhole cover damage or loss in the above three situations.

    4Conclusion

    In this paper, the edge contour tracking algorithm is first used to find all the edges of broken manhole cover based on the pre-processing image. Then, the improved Hough transform is adopted to find all the elliptic regions

    (a) (b) (c) (d) (e)

    Fig.7The manhole cover with a small area of damage. (a) Original image; (b) Dyed image; (c) Image of damaged cover; (d) Preprocessed image; (e) Alarm of damage

    (a) (b) (c) (d) (e)

    Fig.8The manhole cover with a large area of damage. (a) Original image; (b) Dyed image; (c) Image of damaged cover; (d) Preprocessed image; (e) Alarm of damage

    (a) (b) (c) (d) (e)

    Fig.9The completely lost manhole cover. (a) Original image; (b) Dyed image; (c) Image of damaged cover; (d) Preprocessed image; (e) Alarm of damage

    and the color recognition algorithm is adopted to distinguish the manhole cover area from the elliptic regions. The proposed synthesis algorithm can not only eliminate the interference from other elliptic objects, but also exclude the interference of other color approximate regions. Therefore, it can accurately determine the manhole cover area. After the manhole cover area is identified, the damage rate can be calculated by the image matching algorithm. Whether the change rate reaches the threshold value can be used to determine whether it is a suspected damaged manhole cover or not. This synthesis algorithm can effectively improve the processing speed of the traditional Hough transform and reduce inefficient sampling accumulation. The experimental results show that the synthesis algorithm has a good early warning effect in complex backgrounds from different perspectives, and in different videoing time and conditions, such as when the target is partially covered.

    References

    [1]Prasad D K, Leung M K H, Cho S Y. Edge curvature and convexity based ellipse detection method[J].PatternRecognition, 2012, 45(9): 3204-3221.

    [2]Kasemir K U, Betzler K. Detecting ellipses of limited eccentricity in images with high noise levels[J].ImageandVisionComputing, 2003, 21(2): 221-227.

    [3]Wei H, Tan J L. Detection of incomplete ellipse in images with strong noise by iterative randomized Hough transform (IRHT) [J].PatternRecognition, 2007, 41(4): 1268-1279.

    [4]Chen Ying, Xiong Zhang, Wang Yunhong. Location and target detection of manhole cover based on improved Hough[J].OpticalTechnique, 2006, 8(4): 463-467. (in Chinese)

    [5]Liu Jianhua. Automatic location and recognition of manhole cover based on driving on the road with sensor[J].ApplicationResearchofComputers, 2011, 28(8): 3138-3140. (in Chinese)

    [6]Wang Zheshen, Li Cuihua. Building object search and recognition based on improved Hough[J].JournalofImageandGraphics, 2005, 10(4): 463-468. (in Chinese)

    [7]Zhou Xiang, Kong Xiaodong, Zeng Guihua. A new ellipse contour detection method based on Hough[J].ComputerEngineering, 2007, 33(16): 166-171. (in Chinese)

    [8]Li Ziqiang, Teng Hongfei. Generalized Hough fast randomized oval detection[J].JournalofComputer-AidedDesign&ComputerGraphics, 2006, 18(1): 27-33. (in Chinese)

    [9]Chen L, Chen H J, Pan Y, et al. A fast efficient parallel Hough algorithm on LARPBS[J].TheJournalofSuperComputing, 2004, 32(5): 759-762.

    [10]Ou Xiangyin. Image identification and feature point calculation based on Hough to the known figure[D]. Changchun: School of Computer Science and Technology, Jilin University, 2005. (in Chinese)

    [11]Chen R B, Wei H, Li Y, et al. Recognition of lanes on complex roads based on video image[J].AdvancedMaterialsResearch, 2014, 950: 298-305.

    [12]Li Q, Xie Y. Randomised Hough transform with error propagation for line and circle detection[J].PatternAnalysis&Applications, 2003, 6(1): 55-64.

    [13]Yang Fan.Digitalimageprocessingandapplication(Matlab) [M]. Beijing: Chemical Industry Press, 2003. (in Chinese)

    doi:10.3969/j.issn.1003-7985.2015.04.021

    成人美女网站在线观看视频| 精品久久久久久成人av| 国国产精品蜜臀av免费| 亚洲国产精品sss在线观看| 日本免费a在线| 校园人妻丝袜中文字幕| av卡一久久| 午夜老司机福利剧场| 国产在视频线精品| av国产久精品久网站免费入址| 69人妻影院| 日日摸夜夜添夜夜添av毛片| 国产黄色小视频在线观看| 在线观看一区二区三区| 美女主播在线视频| 亚洲伊人久久精品综合| 六月丁香七月| 好男人视频免费观看在线| 成年人午夜在线观看视频 | 亚洲人与动物交配视频| 亚洲精品成人av观看孕妇| 久久精品久久久久久久性| 国产久久久一区二区三区| 免费观看精品视频网站| 亚洲性久久影院| 久久久久免费精品人妻一区二区| 精品国产露脸久久av麻豆 | 深爱激情五月婷婷| 嫩草影院新地址| 国语对白做爰xxxⅹ性视频网站| 亚洲美女搞黄在线观看| 欧美区成人在线视频| 久久国产乱子免费精品| 亚洲精品色激情综合| 日本一本二区三区精品| 精品欧美国产一区二区三| 麻豆成人午夜福利视频| 亚洲四区av| 简卡轻食公司| 一级二级三级毛片免费看| 免费观看无遮挡的男女| 午夜福利视频1000在线观看| 婷婷色综合大香蕉| 午夜精品一区二区三区免费看| 熟妇人妻久久中文字幕3abv| 3wmmmm亚洲av在线观看| 大片免费播放器 马上看| 中文字幕免费在线视频6| 国产亚洲精品av在线| 亚洲成人av在线免费| 18+在线观看网站| 国产伦在线观看视频一区| 国产黄片美女视频| 午夜福利网站1000一区二区三区| av免费观看日本| 免费不卡的大黄色大毛片视频在线观看 | 国产伦在线观看视频一区| 欧美最新免费一区二区三区| 国产精品综合久久久久久久免费| 日产精品乱码卡一卡2卡三| 国产伦精品一区二区三区视频9| 亚洲熟妇中文字幕五十中出| 卡戴珊不雅视频在线播放| 男人和女人高潮做爰伦理| 亚洲人成网站在线播| 女人被狂操c到高潮| 七月丁香在线播放| 2018国产大陆天天弄谢| 欧美激情在线99| 干丝袜人妻中文字幕| 国产成人精品一,二区| 日本熟妇午夜| 精品一区二区三卡| 三级经典国产精品| 免费高清在线观看视频在线观看| 七月丁香在线播放| av国产久精品久网站免费入址| 久久草成人影院| 欧美日本视频| 人体艺术视频欧美日本| 日韩一本色道免费dvd| 久久久久久久久大av| 午夜亚洲福利在线播放| 国产av国产精品国产| 99久久精品国产国产毛片| 狂野欧美白嫩少妇大欣赏| 成人亚洲精品av一区二区| 免费黄频网站在线观看国产| 午夜免费观看性视频| 女人久久www免费人成看片| 国产精品一区二区性色av| 汤姆久久久久久久影院中文字幕 | 国产男人的电影天堂91| 亚洲av一区综合| av在线播放精品| 亚洲va在线va天堂va国产| 秋霞在线观看毛片| 乱码一卡2卡4卡精品| 国产精品99久久久久久久久| 肉色欧美久久久久久久蜜桃 | 国产伦理片在线播放av一区| 国产精品99久久久久久久久| 99久久九九国产精品国产免费| 国产成人精品久久久久久| 国产爱豆传媒在线观看| 久久99精品国语久久久| 美女国产视频在线观看| 狂野欧美激情性xxxx在线观看| 69人妻影院| 在线免费观看不下载黄p国产| 国产大屁股一区二区在线视频| 特级一级黄色大片| 欧美xxⅹ黑人| 超碰97精品在线观看| 高清毛片免费看| av在线亚洲专区| 成人毛片60女人毛片免费| 一级毛片 在线播放| 欧美bdsm另类| av国产免费在线观看| 久久久久久国产a免费观看| 精品久久久久久久末码| 超碰97精品在线观看| 精品久久久久久久人妻蜜臀av| 免费av不卡在线播放| 国产一区二区在线观看日韩| 亚洲av成人精品一二三区| 男女那种视频在线观看| 高清视频免费观看一区二区 | 综合色av麻豆| 国产一区二区在线观看日韩| 色网站视频免费| 你懂的网址亚洲精品在线观看| 国产午夜精品久久久久久一区二区三区| 爱豆传媒免费全集在线观看| 丝袜美腿在线中文| 国产亚洲一区二区精品| 亚洲丝袜综合中文字幕| 九九在线视频观看精品| 十八禁国产超污无遮挡网站| 一个人看视频在线观看www免费| 亚洲精品乱久久久久久| 九九爱精品视频在线观看| 国产黄色免费在线视频| 卡戴珊不雅视频在线播放| 国产综合懂色| 国产精品国产三级专区第一集| 日本爱情动作片www.在线观看| 欧美成人午夜免费资源| 男女啪啪激烈高潮av片| 97热精品久久久久久| 黄色欧美视频在线观看| 国产一区二区在线观看日韩| 日本欧美国产在线视频| 波多野结衣巨乳人妻| 国产一区二区三区av在线| 91精品伊人久久大香线蕉| 亚洲在线观看片| 亚洲精品,欧美精品| 丝袜美腿在线中文| 午夜免费男女啪啪视频观看| www.色视频.com| 中文乱码字字幕精品一区二区三区 | 蜜臀久久99精品久久宅男| 色播亚洲综合网| 又黄又爽又刺激的免费视频.| 国产永久视频网站| 亚洲国产最新在线播放| 国产精品不卡视频一区二区| 亚洲av成人精品一二三区| 好男人在线观看高清免费视频| 男女那种视频在线观看| 人体艺术视频欧美日本| 亚洲成人中文字幕在线播放| 亚洲精品自拍成人| 国产精品国产三级国产av玫瑰| 亚洲欧美日韩卡通动漫| 亚洲一区高清亚洲精品| 久久精品国产鲁丝片午夜精品| 黄色一级大片看看| 国产免费一级a男人的天堂| 一级毛片 在线播放| 插阴视频在线观看视频| 丝瓜视频免费看黄片| 免费看光身美女| 午夜日本视频在线| 2018国产大陆天天弄谢| 亚洲四区av| 久久精品国产亚洲av天美| 精品少妇黑人巨大在线播放| 日韩一本色道免费dvd| 国产av不卡久久| 亚洲va在线va天堂va国产| 国产精品人妻久久久久久| 91久久精品电影网| xxx大片免费视频| 国模一区二区三区四区视频| 免费看av在线观看网站| 国产伦理片在线播放av一区| 日日啪夜夜撸| 搡老乐熟女国产| 午夜激情欧美在线| 看黄色毛片网站| 国内少妇人妻偷人精品xxx网站| 又粗又硬又长又爽又黄的视频| 国产精品美女特级片免费视频播放器| 午夜精品一区二区三区免费看| 国产一区二区三区av在线| 高清欧美精品videossex| 91狼人影院| 在线免费观看的www视频| 亚洲一区高清亚洲精品| 国产爱豆传媒在线观看| 男插女下体视频免费在线播放| 尾随美女入室| 3wmmmm亚洲av在线观看| 久久韩国三级中文字幕| 国产 一区 欧美 日韩| 三级国产精品欧美在线观看| 国产精品日韩av在线免费观看| 国产麻豆成人av免费视频| 波多野结衣巨乳人妻| 中文欧美无线码| 69人妻影院| 高清在线视频一区二区三区| 麻豆久久精品国产亚洲av| 国产亚洲最大av| 亚洲成人一二三区av| 大香蕉久久网| 又大又黄又爽视频免费| 午夜福利在线观看吧| 亚洲在线自拍视频| 中文字幕av成人在线电影| 少妇裸体淫交视频免费看高清| 久久精品久久精品一区二区三区| 欧美xxxx黑人xx丫x性爽| 日日啪夜夜爽| 欧美日韩综合久久久久久| 最近最新中文字幕免费大全7| 一二三四中文在线观看免费高清| 如何舔出高潮| 成年版毛片免费区| 成人国产麻豆网| 嫩草影院入口| 色尼玛亚洲综合影院| 一个人看的www免费观看视频| 国产伦精品一区二区三区视频9| 午夜福利在线观看吧| 亚洲精品一区蜜桃| 亚洲成人av在线免费| av免费在线看不卡| 高清欧美精品videossex| 丝袜喷水一区| 肉色欧美久久久久久久蜜桃 | 日本猛色少妇xxxxx猛交久久| 高清毛片免费看| 日产精品乱码卡一卡2卡三| 禁无遮挡网站| 亚洲欧美日韩东京热| 亚洲经典国产精华液单| 亚洲成人精品中文字幕电影| 99久久人妻综合| 国产一区有黄有色的免费视频 | 99视频精品全部免费 在线| 91午夜精品亚洲一区二区三区| 我的女老师完整版在线观看| 天美传媒精品一区二区| 性色avwww在线观看| av黄色大香蕉| 日本黄大片高清| 一级a做视频免费观看| 国产视频内射| 精品少妇黑人巨大在线播放| 日韩精品有码人妻一区| 能在线免费看毛片的网站| 老师上课跳d突然被开到最大视频| 午夜福利网站1000一区二区三区| 嫩草影院入口| 日韩精品青青久久久久久| 久久午夜福利片| 特级一级黄色大片| 偷拍熟女少妇极品色| 成人美女网站在线观看视频| 日韩大片免费观看网站| 又爽又黄a免费视频| av国产免费在线观看| 非洲黑人性xxxx精品又粗又长| 国产高清不卡午夜福利| 欧美激情在线99| 91精品一卡2卡3卡4卡| 久久人人爽人人爽人人片va| 免费少妇av软件| 成人国产麻豆网| 日本黄色片子视频| 免费看a级黄色片| 国产精品熟女久久久久浪| av福利片在线观看| 亚洲婷婷狠狠爱综合网| h日本视频在线播放| 免费大片18禁| 联通29元200g的流量卡| 国产成年人精品一区二区| 99热这里只有精品一区| 2021天堂中文幕一二区在线观| 高清欧美精品videossex| 少妇高潮的动态图| 啦啦啦中文免费视频观看日本| kizo精华| 床上黄色一级片| 99热这里只有是精品50| 爱豆传媒免费全集在线观看| 十八禁国产超污无遮挡网站| 亚洲国产av新网站| 99久久人妻综合| 日韩中字成人| 中文字幕制服av| 2021天堂中文幕一二区在线观| 国产美女午夜福利| 日韩人妻高清精品专区| 中文乱码字字幕精品一区二区三区 | 国产精品国产三级国产av玫瑰| 又爽又黄无遮挡网站| 亚洲精品乱码久久久久久按摩| 亚洲精品国产av成人精品| 欧美不卡视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 久久久久免费精品人妻一区二区| 一级二级三级毛片免费看| 亚洲人成网站在线观看播放| 有码 亚洲区| 又黄又爽又刺激的免费视频.| 精品国产三级普通话版| 国产麻豆成人av免费视频| 免费大片黄手机在线观看| 中文欧美无线码| 午夜亚洲福利在线播放| 免费黄频网站在线观看国产| 亚洲精品日韩av片在线观看| 亚洲精品成人av观看孕妇| 亚洲av二区三区四区| av在线天堂中文字幕| 日本三级黄在线观看| 久久99精品国语久久久| 2021少妇久久久久久久久久久| 国产一区二区三区综合在线观看 | 国产综合懂色| 美女主播在线视频| 在线免费观看不下载黄p国产| 99视频精品全部免费 在线| 免费观看性生交大片5| av专区在线播放| 99久久中文字幕三级久久日本| 秋霞在线观看毛片| 免费av不卡在线播放| 亚州av有码| 成人二区视频| 天天躁夜夜躁狠狠久久av| 婷婷色综合www| 内射极品少妇av片p| 九色成人免费人妻av| 如何舔出高潮| 国产精品国产三级专区第一集| 国产亚洲av嫩草精品影院| 亚洲最大成人中文| or卡值多少钱| 亚洲人成网站在线播| 久久久久久久久久久免费av| 色吧在线观看| 成人二区视频| ponron亚洲| 两个人视频免费观看高清| 啦啦啦韩国在线观看视频| 草草在线视频免费看| 国产乱来视频区| 国产免费一级a男人的天堂| 毛片一级片免费看久久久久| 天堂av国产一区二区熟女人妻| 99久久中文字幕三级久久日本| 女人被狂操c到高潮| 国产女主播在线喷水免费视频网站 | 波野结衣二区三区在线| 色吧在线观看| 美女高潮的动态| 麻豆久久精品国产亚洲av| 美女cb高潮喷水在线观看| 国产久久久一区二区三区| 国产亚洲5aaaaa淫片| 亚洲av免费高清在线观看| 黄色配什么色好看| 激情 狠狠 欧美| 99热这里只有是精品50| 亚洲精品,欧美精品| 国产亚洲av嫩草精品影院| 亚洲色图av天堂| 欧美zozozo另类| 女人久久www免费人成看片| 性插视频无遮挡在线免费观看| 看十八女毛片水多多多| 精品久久久精品久久久| 国产午夜福利久久久久久| 久久国内精品自在自线图片| 男的添女的下面高潮视频| 亚洲人与动物交配视频| 国产精品国产三级专区第一集| 国产麻豆成人av免费视频| 肉色欧美久久久久久久蜜桃 | 亚洲婷婷狠狠爱综合网| 麻豆成人av视频| 欧美精品国产亚洲| 精品99又大又爽又粗少妇毛片| 亚洲精品自拍成人| 亚洲精品国产av成人精品| 国内精品一区二区在线观看| 一级毛片电影观看| 国产亚洲av片在线观看秒播厂 | 少妇裸体淫交视频免费看高清| 夜夜爽夜夜爽视频| 99热网站在线观看| 国产精品蜜桃在线观看| 男女边摸边吃奶| 亚洲精品国产av成人精品| 看免费成人av毛片| 淫秽高清视频在线观看| 黄色欧美视频在线观看| www.色视频.com| 天堂俺去俺来也www色官网 | 成人特级av手机在线观看| 亚洲在线自拍视频| 久久久成人免费电影| 女人被狂操c到高潮| 亚洲精品亚洲一区二区| 精品欧美国产一区二区三| 精品国产露脸久久av麻豆 | 亚洲精品国产av蜜桃| 美女主播在线视频| 国产大屁股一区二区在线视频| 国产男人的电影天堂91| 91精品一卡2卡3卡4卡| 国产黄片视频在线免费观看| 老司机影院毛片| 一本一本综合久久| 日产精品乱码卡一卡2卡三| 蜜桃亚洲精品一区二区三区| 精品一区二区免费观看| 一级爰片在线观看| 搡老乐熟女国产| 国产一区二区三区综合在线观看 | 亚洲精品影视一区二区三区av| 99视频精品全部免费 在线| 成人毛片a级毛片在线播放| 精品久久久久久久末码| 成人毛片60女人毛片免费| 精品国产露脸久久av麻豆 | 色网站视频免费| 五月天丁香电影| 婷婷色麻豆天堂久久| 久久草成人影院| 亚洲伊人久久精品综合| 亚洲三级黄色毛片| 青春草亚洲视频在线观看| 免费看光身美女| 成人无遮挡网站| 国内精品美女久久久久久| 非洲黑人性xxxx精品又粗又长| 日本色播在线视频| 真实男女啪啪啪动态图| 综合色av麻豆| 亚洲欧美日韩东京热| 最后的刺客免费高清国语| av网站免费在线观看视频 | 男的添女的下面高潮视频| 成年免费大片在线观看| 亚洲美女视频黄频| 大又大粗又爽又黄少妇毛片口| 激情 狠狠 欧美| 色综合色国产| 精品欧美国产一区二区三| 国产亚洲5aaaaa淫片| 亚洲人成网站在线播| 国产一区亚洲一区在线观看| 看黄色毛片网站| 亚洲欧美日韩卡通动漫| 日本黄大片高清| 国产精品一二三区在线看| 黄片wwwwww| 乱码一卡2卡4卡精品| 国产亚洲午夜精品一区二区久久 | 亚洲成人中文字幕在线播放| 又黄又爽又刺激的免费视频.| 国产精品日韩av在线免费观看| 老师上课跳d突然被开到最大视频| 国产亚洲午夜精品一区二区久久 | 美女主播在线视频| 一级爰片在线观看| 五月天丁香电影| 全区人妻精品视频| 97热精品久久久久久| 国产高清不卡午夜福利| 午夜老司机福利剧场| 网址你懂的国产日韩在线| 日日干狠狠操夜夜爽| 久久久久国产网址| 国产高清有码在线观看视频| 精品酒店卫生间| 天天躁日日操中文字幕| 大香蕉久久网| 久久6这里有精品| 免费无遮挡裸体视频| 日日摸夜夜添夜夜爱| 久久精品国产亚洲av涩爱| 91久久精品国产一区二区三区| 日韩电影二区| 亚洲熟妇中文字幕五十中出| 欧美日韩精品成人综合77777| 国产精品久久久久久精品电影小说 | 美女主播在线视频| 国产精品久久久久久久久免| 91aial.com中文字幕在线观看| 色视频www国产| 亚洲欧美精品专区久久| 成人鲁丝片一二三区免费| 只有这里有精品99| 麻豆av噜噜一区二区三区| av播播在线观看一区| 久热久热在线精品观看| 免费观看a级毛片全部| 特大巨黑吊av在线直播| 3wmmmm亚洲av在线观看| 特大巨黑吊av在线直播| 久久精品夜夜夜夜夜久久蜜豆| av线在线观看网站| 搡女人真爽免费视频火全软件| 久久精品国产亚洲av天美| 免费av毛片视频| av黄色大香蕉| av又黄又爽大尺度在线免费看| 国产免费又黄又爽又色| 一边亲一边摸免费视频| 啦啦啦韩国在线观看视频| 国国产精品蜜臀av免费| 国产在线男女| 男女国产视频网站| 亚洲国产日韩欧美精品在线观看| 18禁在线播放成人免费| 天天一区二区日本电影三级| 欧美激情在线99| 精品不卡国产一区二区三区| 国产精品一二三区在线看| 免费电影在线观看免费观看| 国产午夜精品论理片| 亚洲最大成人手机在线| 少妇人妻一区二区三区视频| 成人美女网站在线观看视频| 久久97久久精品| 一区二区三区高清视频在线| 高清毛片免费看| 国产视频首页在线观看| 九九在线视频观看精品| 女人被狂操c到高潮| 国产男人的电影天堂91| 少妇高潮的动态图| 少妇人妻精品综合一区二区| 国产精品一及| 我的老师免费观看完整版| 欧美日韩亚洲高清精品| 高清av免费在线| 看十八女毛片水多多多| 久久人人爽人人片av| 日日撸夜夜添| 熟妇人妻久久中文字幕3abv| 成年人午夜在线观看视频 | 国产视频内射| 亚洲精品色激情综合| 黄色配什么色好看| 51国产日韩欧美| 99热网站在线观看| 欧美激情国产日韩精品一区| 亚洲国产色片| 精品久久久久久电影网| 日日摸夜夜添夜夜添av毛片| 国产精品熟女久久久久浪| 日韩中字成人| 国产伦精品一区二区三区四那| 国产精品久久视频播放| 91aial.com中文字幕在线观看| 日韩视频在线欧美| 欧美成人a在线观看| 亚洲精品,欧美精品| 夜夜爽夜夜爽视频| 毛片一级片免费看久久久久| 免费看a级黄色片| 天堂网av新在线| 亚洲欧美日韩卡通动漫| 午夜日本视频在线| 国产一区有黄有色的免费视频 | 中文资源天堂在线| 2021少妇久久久久久久久久久| 3wmmmm亚洲av在线观看| 人人妻人人看人人澡| 国产亚洲精品久久久com| 男女下面进入的视频免费午夜| 午夜免费男女啪啪视频观看| 午夜免费激情av| av国产免费在线观看| 人人妻人人看人人澡| 亚洲最大成人手机在线| 亚洲欧美成人综合另类久久久| 日韩在线高清观看一区二区三区| 国产精品熟女久久久久浪| 成人国产麻豆网| 国产伦在线观看视频一区| 麻豆成人午夜福利视频| 久久久久久伊人网av| 人妻制服诱惑在线中文字幕| 韩国av在线不卡|