• <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

    一本一本久久a久久精品综合妖精| 日韩成人av中文字幕在线观看| 黄片播放在线免费| 美女福利国产在线| 高清不卡的av网站| 久久97久久精品| 下体分泌物呈黄色| 久久久久久人人人人人| www日本在线高清视频| 亚洲精品国产区一区二| 中文字幕制服av| 欧美成人精品欧美一级黄| 免费看av在线观看网站| 国产成人欧美| 日日啪夜夜爽| 午夜福利一区二区在线看| 女人久久www免费人成看片| 精品国产一区二区三区久久久樱花| 日韩av在线免费看完整版不卡| 婷婷色综合大香蕉| 黑丝袜美女国产一区| 老司机影院毛片| 国产成人精品久久久久久| 久久精品国产亚洲av涩爱| 亚洲成色77777| av.在线天堂| 超色免费av| 亚洲精品,欧美精品| 久久久久精品久久久久真实原创| 大码成人一级视频| 夫妻午夜视频| 精品少妇久久久久久888优播| 王馨瑶露胸无遮挡在线观看| 99久国产av精品国产电影| 欧美成人午夜精品| 亚洲专区中文字幕在线 | 亚洲人成电影观看| 久久久精品国产亚洲av高清涩受| 国产精品免费大片| 亚洲 欧美一区二区三区| 在现免费观看毛片| 欧美亚洲日本最大视频资源| 亚洲国产毛片av蜜桃av| 高清视频免费观看一区二区| 一级毛片黄色毛片免费观看视频| 亚洲人成电影观看| 免费久久久久久久精品成人欧美视频| 男男h啪啪无遮挡| 少妇 在线观看| 国产在线一区二区三区精| 精品酒店卫生间| 波多野结衣一区麻豆| 欧美日韩视频高清一区二区三区二| 亚洲三区欧美一区| 岛国毛片在线播放| 七月丁香在线播放| 免费在线观看视频国产中文字幕亚洲 | bbb黄色大片| 亚洲色图综合在线观看| 精品亚洲成a人片在线观看| 波野结衣二区三区在线| 国产麻豆69| 大片免费播放器 马上看| 嫩草影视91久久| 久久精品人人爽人人爽视色| 少妇人妻 视频| 中文字幕高清在线视频| 丁香六月天网| 天堂中文最新版在线下载| 国产97色在线日韩免费| 99久国产av精品国产电影| 日韩视频在线欧美| 亚洲精品久久久久久婷婷小说| 丰满乱子伦码专区| 在线观看www视频免费| 免费在线观看完整版高清| av在线观看视频网站免费| 黄片播放在线免费| 老汉色∧v一级毛片| 不卡视频在线观看欧美| 精品一区二区三区四区五区乱码 | 2021少妇久久久久久久久久久| 国产一区有黄有色的免费视频| 亚洲欧美精品自产自拍| 夫妻午夜视频| 亚洲国产欧美日韩在线播放| 19禁男女啪啪无遮挡网站| 操出白浆在线播放| 亚洲精品国产av成人精品| 天天影视国产精品| 1024视频免费在线观看| 色视频在线一区二区三区| 999精品在线视频| 亚洲国产日韩一区二区| 亚洲免费av在线视频| 日韩欧美精品免费久久| 久久久国产欧美日韩av| 久久久亚洲精品成人影院| 成人影院久久| 国产爽快片一区二区三区| 国产精品一区二区在线观看99| 亚洲伊人久久精品综合| 久久精品亚洲熟妇少妇任你| 热99久久久久精品小说推荐| 在线观看免费午夜福利视频| 我要看黄色一级片免费的| 男女无遮挡免费网站观看| 男女边摸边吃奶| 国产精品女同一区二区软件| 亚洲精品一区蜜桃| 天天操日日干夜夜撸| www.熟女人妻精品国产| 国产精品蜜桃在线观看| 天天添夜夜摸| 欧美少妇被猛烈插入视频| 亚洲三区欧美一区| 亚洲成色77777| 少妇的丰满在线观看| 亚洲一区中文字幕在线| 91精品三级在线观看| 免费在线观看黄色视频的| 亚洲婷婷狠狠爱综合网| 午夜免费观看性视频| 亚洲精品自拍成人| 高清不卡的av网站| 精品亚洲成国产av| 亚洲精品久久久久久婷婷小说| 久久久久精品国产欧美久久久 | 五月开心婷婷网| 老鸭窝网址在线观看| 美女扒开内裤让男人捅视频| 在现免费观看毛片| 男女无遮挡免费网站观看| 国产精品.久久久| 精品国产一区二区三区四区第35| 成年av动漫网址| 国产精品二区激情视频| 天堂中文最新版在线下载| 在线观看免费日韩欧美大片| 各种免费的搞黄视频| 天美传媒精品一区二区| 欧美日韩av久久| 久久人人97超碰香蕉20202| kizo精华| 男女下面插进去视频免费观看| 99热网站在线观看| 国产成人免费无遮挡视频| 国产片特级美女逼逼视频| 飞空精品影院首页| kizo精华| 久久精品久久久久久噜噜老黄| 人妻人人澡人人爽人人| 午夜福利一区二区在线看| 欧美精品一区二区免费开放| 夜夜骑夜夜射夜夜干| 日韩欧美一区视频在线观看| 国产极品粉嫩免费观看在线| 黄色视频在线播放观看不卡| 精品酒店卫生间| 少妇的丰满在线观看| 亚洲五月色婷婷综合| 国产在线视频一区二区| 高清av免费在线| 久久精品久久久久久久性| 国产男女超爽视频在线观看| 亚洲精品国产色婷婷电影| 午夜精品国产一区二区电影| 人成视频在线观看免费观看| 少妇猛男粗大的猛烈进出视频| 午夜福利视频在线观看免费| 大香蕉久久成人网| 国产免费一区二区三区四区乱码| 午夜免费男女啪啪视频观看| 国产成人免费观看mmmm| 超碰97精品在线观看| 亚洲av电影在线进入| 2018国产大陆天天弄谢| 美女视频免费永久观看网站| 欧美在线一区亚洲| 99久久综合免费| 国产精品香港三级国产av潘金莲 | 制服人妻中文乱码| 波多野结衣av一区二区av| 视频在线观看一区二区三区| 午夜91福利影院| 精品一品国产午夜福利视频| 欧美日韩福利视频一区二区| 黄色视频在线播放观看不卡| 黄色一级大片看看| 搡老乐熟女国产| 久久人人97超碰香蕉20202| 国产女主播在线喷水免费视频网站| 成人漫画全彩无遮挡| 女人高潮潮喷娇喘18禁视频| 国产精品无大码| av电影中文网址| 亚洲视频免费观看视频| 国产成人欧美在线观看 | 国产一卡二卡三卡精品 | 婷婷色av中文字幕| 高清黄色对白视频在线免费看| 午夜免费鲁丝| 男人添女人高潮全过程视频| 视频区图区小说| 国产 精品1| 日韩熟女老妇一区二区性免费视频| 国产精品免费大片| 男女午夜视频在线观看| 国产黄频视频在线观看| 色婷婷久久久亚洲欧美| 亚洲中文av在线| 国产男人的电影天堂91| 丝袜脚勾引网站| 人妻 亚洲 视频| 街头女战士在线观看网站| 观看美女的网站| 成年av动漫网址| 国产高清不卡午夜福利| av片东京热男人的天堂| 好男人视频免费观看在线| 97在线人人人人妻| 性高湖久久久久久久久免费观看| 久久精品久久久久久噜噜老黄| 一区二区三区四区激情视频| 三上悠亚av全集在线观看| 亚洲人成77777在线视频| 亚洲成国产人片在线观看| 97人妻天天添夜夜摸| 久久av网站| 日韩精品有码人妻一区| 国产免费视频播放在线视频| 男人舔女人的私密视频| 欧美激情极品国产一区二区三区| 中文字幕人妻熟女乱码| 91精品国产国语对白视频| 国产爽快片一区二区三区| www.av在线官网国产| 国产精品国产三级专区第一集| 精品少妇久久久久久888优播| 天天躁夜夜躁狠狠躁躁| 国产免费福利视频在线观看| 国产成人欧美在线观看 | 久久久精品国产亚洲av高清涩受| 成人影院久久| 女的被弄到高潮叫床怎么办| 欧美黑人欧美精品刺激| 亚洲激情五月婷婷啪啪| 中文字幕色久视频| 国产精品嫩草影院av在线观看| 亚洲精品久久成人aⅴ小说| 亚洲av欧美aⅴ国产| 国产男人的电影天堂91| 久久精品亚洲熟妇少妇任你| 日韩电影二区| 国产精品一区二区在线观看99| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧洲精品一区二区精品久久久 | 亚洲欧美一区二区三区久久| av在线app专区| 男女之事视频高清在线观看 | 如日韩欧美国产精品一区二区三区| 一区二区三区四区激情视频| 韩国精品一区二区三区| 丝袜美腿诱惑在线| 中文精品一卡2卡3卡4更新| 国产黄色视频一区二区在线观看| 成人亚洲欧美一区二区av| 七月丁香在线播放| 国产精品久久久久久精品电影小说| 日韩 欧美 亚洲 中文字幕| 欧美黄色片欧美黄色片| 国产精品三级大全| 美女中出高潮动态图| 最黄视频免费看| 一级a爱视频在线免费观看| 在线天堂中文资源库| 精品久久蜜臀av无| av在线app专区| 国产黄频视频在线观看| 韩国av在线不卡| 狠狠婷婷综合久久久久久88av| 亚洲,一卡二卡三卡| 国产不卡av网站在线观看| 中文字幕色久视频| 高清在线视频一区二区三区| 国产极品粉嫩免费观看在线| 国产精品蜜桃在线观看| 男人操女人黄网站| 亚洲精品国产区一区二| av网站免费在线观看视频| 9191精品国产免费久久| 午夜福利免费观看在线| 日日摸夜夜添夜夜爱| 黄片小视频在线播放| 久久久久国产一级毛片高清牌| 午夜精品国产一区二区电影| 别揉我奶头~嗯~啊~动态视频 | 亚洲欧美一区二区三区久久| 亚洲av日韩精品久久久久久密 | 99九九在线精品视频| 久久鲁丝午夜福利片| 久久ye,这里只有精品| 搡老岳熟女国产| 亚洲欧美成人综合另类久久久| 国产精品一二三区在线看| 亚洲av欧美aⅴ国产| 高清欧美精品videossex| 19禁男女啪啪无遮挡网站| 制服人妻中文乱码| 1024香蕉在线观看| 在线观看国产h片| 久热这里只有精品99| 777米奇影视久久| 乱人伦中国视频| 欧美日韩福利视频一区二区| 国产精品蜜桃在线观看| 欧美人与性动交α欧美精品济南到| 在线免费观看不下载黄p国产| 亚洲国产精品一区二区三区在线| 在现免费观看毛片| 久久99一区二区三区| 欧美日韩视频高清一区二区三区二| 黑人巨大精品欧美一区二区蜜桃| 男女之事视频高清在线观看 | 国产av码专区亚洲av| 国产人伦9x9x在线观看| 中文字幕制服av| 叶爱在线成人免费视频播放| 亚洲欧洲日产国产| 亚洲国产成人一精品久久久| 国产精品av久久久久免费| 日日撸夜夜添| 一级黄片播放器| 久久久久久久久久久免费av| 国产成人免费观看mmmm| 免费日韩欧美在线观看| 亚洲成人国产一区在线观看 | 少妇人妻 视频| 亚洲综合精品二区| 考比视频在线观看| 欧美人与性动交α欧美精品济南到| 无限看片的www在线观看| 欧美老熟妇乱子伦牲交| 视频区图区小说| 午夜激情av网站| 亚洲精品美女久久av网站| 日韩中文字幕视频在线看片| 超碰97精品在线观看| 国产伦人伦偷精品视频| 丰满迷人的少妇在线观看| tube8黄色片| 国产一区二区三区av在线| av电影中文网址| 人妻一区二区av| 成年女人毛片免费观看观看9 | 亚洲成国产人片在线观看| 一区二区三区激情视频| 国产精品国产三级专区第一集| 侵犯人妻中文字幕一二三四区| 久久免费观看电影| 一边摸一边做爽爽视频免费| 久久99精品国语久久久| 黄网站色视频无遮挡免费观看| 午夜日本视频在线| 秋霞在线观看毛片| 香蕉国产在线看| 啦啦啦啦在线视频资源| 别揉我奶头~嗯~啊~动态视频 | 少妇 在线观看| 麻豆av在线久日| 一级黄片播放器| 婷婷色麻豆天堂久久| 狂野欧美激情性bbbbbb| 国产黄频视频在线观看| 亚洲国产欧美日韩在线播放| 日韩不卡一区二区三区视频在线| 黑人巨大精品欧美一区二区蜜桃| 七月丁香在线播放| 日韩视频在线欧美| www.熟女人妻精品国产| 欧美黄色片欧美黄色片| 久久97久久精品| 国产在线免费精品| 亚洲欧美一区二区三区久久| 亚洲欧美精品自产自拍| 水蜜桃什么品种好| 亚洲精品久久久久久婷婷小说| 黄网站色视频无遮挡免费观看| 在线观看免费日韩欧美大片| 国产精品蜜桃在线观看| 精品久久久久久电影网| 国产亚洲av片在线观看秒播厂| 国产97色在线日韩免费| 美女中出高潮动态图| 免费久久久久久久精品成人欧美视频| 国产成人午夜福利电影在线观看| 综合色丁香网| 亚洲中文av在线| 亚洲国产欧美网| 免费高清在线观看日韩| 国产精品熟女久久久久浪| 欧美在线一区亚洲| 麻豆精品久久久久久蜜桃| 青草久久国产| 三上悠亚av全集在线观看| 午夜老司机福利片| tube8黄色片| 色综合欧美亚洲国产小说| 秋霞伦理黄片| 久久影院123| 亚洲久久久国产精品| 啦啦啦在线免费观看视频4| 看十八女毛片水多多多| av电影中文网址| 亚洲av日韩在线播放| 老司机在亚洲福利影院| 成人手机av| 久久久久久久大尺度免费视频| 日韩精品有码人妻一区| 波多野结衣av一区二区av| 色综合欧美亚洲国产小说| 一区二区日韩欧美中文字幕| 久久久久久久久免费视频了| 五月开心婷婷网| 亚洲国产最新在线播放| 国产成人精品久久久久久| 九色亚洲精品在线播放| 国产黄频视频在线观看| 自线自在国产av| 中国三级夫妇交换| avwww免费| av天堂久久9| 老汉色av国产亚洲站长工具| 欧美日韩国产mv在线观看视频| 香蕉丝袜av| av女优亚洲男人天堂| 又大又爽又粗| 五月天丁香电影| 伊人亚洲综合成人网| 亚洲熟女毛片儿| 69精品国产乱码久久久| av线在线观看网站| 大码成人一级视频| 亚洲第一青青草原| 热re99久久国产66热| 你懂的网址亚洲精品在线观看| 久久久久网色| 亚洲av国产av综合av卡| 亚洲国产欧美日韩在线播放| 日韩一卡2卡3卡4卡2021年| av在线老鸭窝| 久久久久精品国产欧美久久久 | 午夜激情av网站| 黄频高清免费视频| 一边摸一边做爽爽视频免费| 99香蕉大伊视频| 亚洲精品一区蜜桃| 国产亚洲欧美精品永久| kizo精华| 毛片一级片免费看久久久久| 亚洲欧美一区二区三区黑人| 亚洲欧美精品综合一区二区三区| 国产精品一国产av| videos熟女内射| 午夜久久久在线观看| 国产成人一区二区在线| 午夜91福利影院| 亚洲精品aⅴ在线观看| 日韩中文字幕欧美一区二区 | 99久久综合免费| 亚洲欧美精品自产自拍| 国产精品一二三区在线看| 国产国语露脸激情在线看| 国产精品久久久av美女十八| 激情五月婷婷亚洲| 日韩,欧美,国产一区二区三区| 人妻人人澡人人爽人人| 国产黄色免费在线视频| 韩国精品一区二区三区| 免费久久久久久久精品成人欧美视频| 亚洲在久久综合| 一级毛片我不卡| 久久免费观看电影| 18禁观看日本| 成人三级做爰电影| 99re6热这里在线精品视频| 久久精品国产综合久久久| 国产精品久久久久久人妻精品电影 | videosex国产| 曰老女人黄片| 别揉我奶头~嗯~啊~动态视频 | 国产爽快片一区二区三区| 国产一区亚洲一区在线观看| 欧美激情 高清一区二区三区| 黄色毛片三级朝国网站| 精品国产一区二区三区久久久樱花| 欧美日韩国产mv在线观看视频| 精品国产一区二区三区久久久樱花| 精品国产国语对白av| 女性生殖器流出的白浆| 黄色怎么调成土黄色| 亚洲久久久国产精品| 精品一区二区三区av网在线观看 | 久热这里只有精品99| 日本猛色少妇xxxxx猛交久久| 一级爰片在线观看| 日韩中文字幕视频在线看片| 免费少妇av软件| 午夜激情久久久久久久| 精品一区二区三卡| 亚洲国产中文字幕在线视频| 一区二区日韩欧美中文字幕| 午夜久久久在线观看| 日韩电影二区| 色婷婷av一区二区三区视频| 国产乱人偷精品视频| 午夜日本视频在线| 曰老女人黄片| 亚洲国产欧美一区二区综合| 女人爽到高潮嗷嗷叫在线视频| 亚洲成色77777| netflix在线观看网站| 高清视频免费观看一区二区| 成人18禁高潮啪啪吃奶动态图| 国产精品国产av在线观看| 亚洲欧洲精品一区二区精品久久久 | 纯流量卡能插随身wifi吗| 久久精品久久久久久噜噜老黄| 韩国高清视频一区二区三区| 中文欧美无线码| 亚洲综合精品二区| 国产精品一二三区在线看| 男女之事视频高清在线观看 | xxx大片免费视频| 久热这里只有精品99| www.熟女人妻精品国产| 人人妻人人爽人人添夜夜欢视频| 久久综合国产亚洲精品| 久久av网站| 亚洲,一卡二卡三卡| 黑丝袜美女国产一区| 高清欧美精品videossex| 亚洲欧美一区二区三区黑人| 下体分泌物呈黄色| 亚洲七黄色美女视频| 国产精品久久久久久精品电影小说| 国产女主播在线喷水免费视频网站| 99久国产av精品国产电影| 欧美日韩一区二区视频在线观看视频在线| 不卡av一区二区三区| 亚洲精品aⅴ在线观看| 少妇人妻精品综合一区二区| 国产精品国产三级专区第一集| 激情视频va一区二区三区| 国产一区二区三区av在线| 纯流量卡能插随身wifi吗| 亚洲国产av新网站| 国产亚洲精品第一综合不卡| 亚洲国产欧美在线一区| 日韩 欧美 亚洲 中文字幕| 亚洲欧洲日产国产| 亚洲第一青青草原| 国产av国产精品国产| 国产男女超爽视频在线观看| 高清欧美精品videossex| 亚洲国产欧美一区二区综合| 建设人人有责人人尽责人人享有的| 亚洲色图 男人天堂 中文字幕| 天天躁夜夜躁狠狠躁躁| 色婷婷av一区二区三区视频| 久久国产精品男人的天堂亚洲| 黑人猛操日本美女一级片| 深夜精品福利| 91精品伊人久久大香线蕉| 久久人妻熟女aⅴ| 97精品久久久久久久久久精品| 亚洲欧洲国产日韩| a级片在线免费高清观看视频| 精品久久久精品久久久| 成人手机av| 美女中出高潮动态图| 91精品伊人久久大香线蕉| 国产精品国产三级专区第一集| 亚洲国产毛片av蜜桃av| 中文字幕av电影在线播放| 亚洲精品视频女| 国产日韩欧美亚洲二区| 黄色怎么调成土黄色| 欧美精品人与动牲交sv欧美| 日韩制服丝袜自拍偷拍| 国产一区二区 视频在线| 高清黄色对白视频在线免费看| 九九爱精品视频在线观看| 搡老乐熟女国产| 久久精品国产亚洲av高清一级| 成年av动漫网址| 亚洲av电影在线观看一区二区三区| 制服丝袜香蕉在线| 伊人久久大香线蕉亚洲五| 青春草亚洲视频在线观看| 51午夜福利影视在线观看| 晚上一个人看的免费电影| 男女国产视频网站| 国产av码专区亚洲av| 晚上一个人看的免费电影| 精品一区在线观看国产| 肉色欧美久久久久久久蜜桃| 女人被躁到高潮嗷嗷叫费观| 伦理电影免费视频| 大话2 男鬼变身卡| 日韩制服丝袜自拍偷拍| 黄片播放在线免费| 国产精品 国内视频| 波多野结衣av一区二区av| 91精品国产国语对白视频|