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

    Measurement of particle size based on digital imaging technique*

    2013-06-01 12:29:57CHENHong陳紅TANGHongwu唐洪武LIUYun劉云WANGHao王浩LIUGuiping劉貴平
    水動力學研究與進展 B輯 2013年2期
    關(guān)鍵詞:劉云陳紅王浩

    CHEN Hong (陳紅), TANG Hong-wu (唐洪武), LIU Yun (劉云), WANG Hao (王浩), LIU Gui-ping (劉貴平)

    State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing 210098, China, E-mail: chh_hhu@hhu.edu.cn

    Measurement of particle size based on digital imaging technique*

    CHEN Hong (陳紅), TANG Hong-wu (唐洪武), LIU Yun (劉云), WANG Hao (王浩), LIU Gui-ping (劉貴平)

    State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing 210098, China, E-mail: chh_hhu@hhu.edu.cn

    (Received May 23, 2012, Revised August 22, 2012)

    To improve the analysis methods for the measurement of the sediment particle sizes with a wide distribution and of irregular shapes, a sediment particle image measurement, an analysis system, and an extraction algorithm of the optimal threshold based on the gray histogram peak values are proposed. Recording the pixels of the sediment particles by labeling them, the algorithm can effectively separate the sediment particle images from the background images using the equivalent pixel circles with the same diameters to represent the sediment particles. Compared with the laser analyzer for the case of blue plastic sands, the measurement results of the system are shown to be reasonably similar. The errors are mainly due to the small size of the particles and the limitation of the apparatus. The measurement accuracy can be improved by increasing the Charge-Coupled Devices (CCD) camera resolution. The analysis method of the sediment particle images can provide a technical support for the rapid measurement of the sediment particle size and its distribution.

    particle size, extraction algorithm of optimal threshold, equivalent circle transformation

    Introduction

    Recently, many methods were developed to determine the sediment particle size, including the original manual method and the automated instrumental method. These methods can be classified as the sieving method, the sedimentation method, the resistance method, the microscope method, and the laser particle sizer method. The traditional analysis methods are usually characterized by the complex operations, the low automation levels, the time-consuming operations, and the low speed. On the other hand, the laser particle sizer method has a high degree of automation, but it is not suitable for a large-scale use because of its high price and the maintenance cost. The particle size should be confined within 2 mm, otherwise, the instrument would fail to work properly. The measuring ranges of these methods are all limited. Thus, a variety of methods are usually adopted jointly to analyze a wide-gradation sediment. Due to the differences of the measuring principles and the representations of different methods, there are certain link errors, and the grading curve is not smoothly continuous.

    With the development of computer technologies, the chips, Charge-Coupled Devices (CCD), and especially, the image processing technologies are widely applied in the particle analysis fields, such as in the powder, mineral, chemical and medical industries. Misirli et al.[4]and Apostolopoulos et al.[5]used the image processing technology to identify cells, and put forward the “best rounding fitting” method to measure the cell size. Souza and Menegalli[6]used the digital image processing technology to measure the particle size, and analyzed the particle morphology charcteristics. Carter et al.[7,8]developed a set of particle size measurement devices, and realized the particle sizemeasurement in pipelines. The rapid development of image processing technologies improves the measurement of the particle size, the particle morphology, and other particle parameters. The commercial particle analysis systems were successfully applied to the particle analysis fields[9,10]..However, the image technology of the sediment particle is far behind other industries. Using the digital image processing technology, Li et al.[11]obtained the representative diameter, the flat index, the triaxial shape parameters, the roundness coefficient, and other geometric parameters of pebbles, and provided a new method for the study of the mechanism of the pebble movement. As a whole, the applications in this area are quite limited, and the main reasons are as follows: (1) The sediment involves a wide range of particle sizes, the poor uniformity of the particle distribution, and various morphologies of particles, which increases the difficulty in the image processing. (2) For the sediment, not only the morphological parameters of a single particle, but also the group characteristics are important. For the analysis of sediment particles, more samples are required, and the measured results should be representative. In general, it is difficult for particle analysis systems with microscope imaging to collect enough sediment samples for analysis of sediment gradation characteristics. (3) To meet the need of a wide gradation, neither the processing accuracy nor the system automation can be satisfactory. As a result, the integration level of the system is low and the operation is complex. These facts affect the development of image processing technologies in the sediment particle field. Liu[12]developed a simple sediment particle size analysis system, but the automation level and the measurement accuracy for a small-scale sediment remain to be improved.

    This study focuses on how to obtain accurate image information of wide-gradation sediment particles and to accurately measure the particle size through integrating the image acquisition device, the image processing, and the fast statistics technologies of sediment particle pixels. Based on the existing particle analysis system, a suitable system for the sediment particle image analysis is developed. Through the image acquisition device with an optimal design, a high-quality and high-stability Light Emitting Diode (LED) light source array is built, to improve the accuracy of the sediment particle projection and to avoid the submerging small-scale sediment. Based on the peak value of the gray histogram, an optimal threshold extraction algorithm is designed to improve the extraction accuracy of the wide-gradation sediment particle. Through the labeling records and the statistics of the sediment particle pixels, the processing speed and the measurement accuracy of the system are improved, and the operation process is simplified.

    1. Experimental devices

    The system includes the measurement point locating device, a black and white CCD camera, LED light sources, computer, and the image processing device, as shown in Fig.1. The uniform sediment particles are scattered at the measurement positions under stable LED light sources. The sediment particle images are obtained by the digital CCD camera and sent to the computer. The image information is extracted through the image processing. Finally, the pixels of the sediment particle images are transformed to equivalent pixel circles with the same diameters.

    With the black and white CCD camera adopted here, the amount of image data can effectively reduced and the processing speed is improved. The resolution of the CCD chip is 1 628 × 1 236. The pixel size is 4.4 μm × 4.4 μm. The data interface is set to be 1 394, and the output is an 8-bit grayscale image. A lens with a focal length of 16 mm is used. The measurement point locating device includes four LED light source linear arrays, each of which consists of five 2 W LED lights (the power could be changed by current). The LED lights are arranged at equal distances from each other to ensure the sediment particles under a uniform and relatively stable light intensity to obtain an accurate projection of the sediment particles.

    The full dispersion of the sediment particles is under an accurate measurement condition. The large particles such as the pebbles can be dispersed easily, but it is difficult to disperse small particles, which are always dispersed by artificial means. This means that a particle lap would appear, and the lap image will be determined as a single particle in the image processing.

    For dispersing the sediment, a special oscillator is designed as shown in Fig.1. This oscillator includes several oscillating plates, and they are divided into four groups, numbered as O1, O2, O3and O4. Every group has two oscillating plates. The oscillation frequency and amplitude can be adjusted according to the sediment particle sizes. For a particle size of lessthan 0.1 mm, the oscillation frequency can be set at 3 Hz to 5 Hz, and the amplitude can be set at 2 mm to 5 mm. Under the action of the oscillator, the sediment particles can be fully dispersed.

    Fig.1 Structure of experimental devices

    2. Information extraction of sediment particles

    2.1 Sediment particle image acquisition

    The sediment particle images are directly related to the measurement precision of the sediment particle diameter and the analysis of the group features. The number of the sediment particles in the conventional particle image analysis methods of microscope imaging is very small, and it is difficult to obtain the group characteristics of the sediment. The ordinary lighting is a solid light source, which emits light in all directions. When the number of the sediment particles is small, the light in different directions passes through the sediment, which greatly reduces the sediment projection resolution. Thence a high light intensity is required to cover up the tiny sediment particles, and clear sediment particle images can not be obtained by a light source of low intensity. Therefore, it is necessary to adjust the light intensity according to the light resistance characteristics of the sediment for the sediment particle image acquisition.

    Fig.2 Image acquisition process

    This experimental system chooses the LED as the light sources, with which, the light intensity can be adjusted by the current. To obtain accurate particle images, the image acquisition process as shown in Fig.2 is adopted. First, with the image acquisition, 20 frames of the sediment particle images are captured with the LED current of 0.8 A. Then the LED current is adjusted, and 20 frames are captured with the LED current of 0.75 A. These two sets of 20 frame images are processed image by image, respectively. Through the threshold transformation, the sediment particle image pixels are counted, and the results for the two sets are compared.n0is the number of sediment image pixels obtained with the LED current of 0.8 A, n1is the number of sediment image pixels obtained with the LED current of 0.75 A. When the difference between n0and n1is less than 0.5%, all sediment particles are considered to be counted. If it is more than 0.5%, some small particles might be submerged in the light. Thus, we should adjust the LED current to reduce the light intensity until the difference between ni-1and ni(iis the adjusting steps for the LED current) is less than 0.5%.

    2.2 Optimal threshold algorithm

    When the control parameters are kept constant, the image quality is relatively stable because of the low noise-signal ratio of the uniform light sources, and the camera and other devices in the closed measurement positions. The collected images are grayscale images with gray-scale values ranging from 0 to 255. The sediment particle images are separated from the background images by the threshold transform, which could be expressed as

    whereT is the threshold value, and f( x, y)is the gray-scale value at the point(x, y)in the pixel lattice of images.Aand Bare the gray-scale values obtained by the threshold transform, which could be termed as the binarization transform when A =1and B=0.

    Recently, the threshold transform is widely applied to the gray-scale distribution of images. The grayscale distribution diagrams, also named the histograms, have three types: the unimodal, the bimodal, and the multi-peak types. The unimodal histogram describes the non-uniform distribution of the background and the objects, with little to do with the gray-scale distribution, and partially shows the information of the objects. The bimodal histogram shows the existence of the gray level area of the objects and the background in the images, focusing on the difference between the objects and the background. The troughs of the gray values between two peaks are considered as the threshold. The threshold transform is used to separate the background from the objects. The multipeak histogram is a complicated image containing three or more peaks, with each peak corresponding to an object separated in each zone before the image is divided into zones.

    The sediment particles are represented by different types of histograms according to different sizes and materials. The sediments with the uniform material and well-distributed sizes of particles are to be represented by unimodal histograms, while the sediments with the non-uniform material and a particle distri-bution are to be represented by multi-peak histograms. Thus, it is difficult to meet the actual needs just by adopting a single sediment segmentation.

    Fig.3 Original image

    Fig.4 Result image with threshold of 190

    Fig.5 Segmented image with threshold of 182

    With the threshold segmentation as shown in Fig.3, the segmented images are obtained with the threshold values of 190 and 182, respectively. Compared with the image in Fig.4, the image in Fig.5 is shown to be better segmented, therefore, the separation result between the sediment particles and the background is directly affected by the threshold value. The experiments show that there exists a certain relationship between the optimal threshold value of the sediment and the peak value of the gray-scale, which may be expressed based on the statistics as follows

    where Tbestis the optimal threshold value represented by an integer,G1is the peak value of the grayscale. The peak value in the gray histogram of the original image shown in Fig.3 is 223, and the optimal threshold value calculated by Eq.(2) is 182, as well consistent with the treatment effect. Using the extraction algorithm of the optimal threshold based on the peak value of the gray histogram, the dynamic threshold segmentation based on images of different sediment particles could be achieved, with results better and more flexible compared with the fixed threshold value.

    Fig.6 Single threshold value segmentation

    Fig.7 Original image

    Fig.8 Partitioned threshold value segmentation

    Although a homogeneous LED linear light source is used, it is still difficult to guarantee a consistent optical path for each sediment particle, especially, as it is affected by the reflection from the wall at the measurement positions in the area near the wall, resulting in differences in the light intensity received by different sediment particles. Under such a condition, the treatment is not very effective (Fig.6) in segmenting the particles from the background (Fig.7) with a single threshold value. Thus, the image should be divided into various regions according to the distribution of the gray histogram in each region, and thenthe threshold segmentation is conducted in each region of the image based on Eq.(2). The results of the treatment of the partitioned threshold segmentation are shown in Fig.8, which implies that the partitioned treatment is more efficient for the segmentation between the particles and the background. After the threshold segmentation, the gray value is 0 for the sediment particles and 255 for the background.

    2.3 Fast pixel statistics of particle image

    To calculate the area of the particle image, the pixels in the image should be counted. The labeling method is applied in this system to obtain the pixel statistics. The principle of the method is as follows: we mark the same label for the pixel points of the same target in the image while different labels are for different images, and the size of the target is obtained by using the statistics of the pixels with the same label.

    The processing steps are as follows: (1) Determining the pixel characteristics of the sediment particles, setting the gray value 0 for the sediment particles, and 255 for the background, and making the pixel points of the same particle interconnected. (2) Determining the starting mark point and the initial tagged value, which are, in this system, respectively, the zero point of the image coordinate (0,0)and –1, indicating that the labeling has not started. (3) Determining the labeling principle that only the pixels of the sediment particles but not the background are labeled, with the interconnected pixels taking the same label. (4) Establishing three one-dimensional tag arrays, PX[ a],PY[ a]and PV[ a], for the x coordinate, the ycoordinate, and the corresponding tagged value of the pixels, respectively, and a is the total number of the pixels of the sediment particles. (5) Interpolating the gray values for each pixel along the row direction in turn, and recording the tagged value according to the sediment particle pixel points. The tagged value increases by 1 when a particle pixel point is shifted to a background one, while it keeps invariant when the background pixel point is shifted to a particle one or a particle pixel point is shifted to another. (6) Calculating the number of the pixels of the same tagged value in the tag array to obtain the sediment particle size.

    Fig.9 Sketch of pixel statistics

    Figure 9, as a processing example, shows the processing program: as f(0,0)=255, the pixel point (0, 0) will not be marked, and the tagged value keeps invariant,f(2,0)=0, so the pixel point (2, 0) is a sediment particle pixel point, and the tagged value should increase by 1 and be saved to the tag array: PX [1]=2,PY [1]=0and PV[1]=0, hence the number of pixels will be 1,f (3,0)=225, so the pixel point (3, 0) is not a particle pixel point with no mark, and the tagged value increases by 1; the points are analyzed in turns,f (8,0)=0,(8,0)is a sediment particle pixel point, hencePX [2]=8,PY[2]= 0,PV [2]=1.f(14,0)=0, (14, 0) is a sediment particle pixel point,PX [3]=14,PY [3]=0,PY [3]=2. f (15,0)=0, (15, 0) is a sediment particle pixel point, PX [4]=15,PY [4]=0,PV[4]=2. After analyzing the pixels in the first row, the second row will be considered, in whichf (2,1)=0, so (2,1)is a particle pixel point, now the distance between thePX[ a], PY[ a]of the first four pixel points and (2,1)could be calculated. If the distance is less than or equal to 1, it means that (2,1)is interconnected with this particle pixel and would take the same tagged value,PX[5]= 2,PY [5]=1and PV[4]=0. The other pixel points are analyzed similarly, and a value of 2 is marked, which means that there are three sediment particles in Fig.9. Statistics of the tagged values show that the three particles occupy 11, 12 and 15 pixel points, respectively.

    Fig.10 Equivalent area transformation

    2.4 Equivalent transformation

    The sediment particle images obtained through the image analysis are the largest cross-section projections with various shapes. Thus, it is necessary to establish a unified method for the sediment particle size calculation through an equivalent transformation between an image of the sediment particles and a circle with the same value of pixels and diameter as the sediment particles. The equivalent principle is shown in Fig.10, and the equivalent circle diameter is

    where n is the number of pixels of the sediment par-ticle image,αis the pixel size, which is 4.4 μm for MVC1000MF black and white camera adopted in this system,Lis the vertical distance from the camera lens to the plane where the sediment particles are measured, andfis the focal length of the camera lens. The minimum size measured by this system is determined byLandα.Lis adjustable within the range from 100 mm to 400 mm, and the corresponding minimum diameter measurable ranges from 12.42 μm to 24.83 μm. The measurement of a smaller particle size is also possible by using a higher-resolution camera and a lens with a greater focal length. With the smaller sediment particle size, the camera range drops as well, resulting in the decrease in the number of sediment particles measurable simultaneously.Lshould be reasonably set according to the sediment sample to be measured in order to meet the measurement needs.

    IfLandfare set at 280 mm and 16 mm, respectively, the equivalent circle diameters for the three sediment particles in Fig.4 will be as follows

    3. Application and error analysis

    The gradation curve is usually used to represent the sediment particle size information in the analysis of sediment particles, for the number of particles is huge, and the sediment usually moves as a group. The abscissa represents the particle size in a gradation curve, while the vertical axis represents the progressively accumulated percentage value of the quality of sediment particles of multi-size levels.

    To test the measurement effect of this system, a comparison is made between the results of the Malvern laser particle size analyzer and the image analysis measurement by selecting spherical and irregularly shaped plastic sands as measurement objects, as shown in Fig.11. Where Cpis the cumulative percent andDis the partice size. It is shown that the shapes of gradation curves obtained by the two methods are similar. The median diameter of the laser particle method is 0.33 mm, while that of the image measurement method is 0.34 mm. The deviation of the results between the two methods mainly comes from small particles. The quality percentage of the sediment particles measured by the image measurement is 0 when the particle size is less than 0.04 mm, and the cumulative quality percentage by the image measurement method is smaller than that of the laser particle size analyzer when the particle size is less than 0.7 mm. The image measurement method is limitedly applicable to the minimum measurable size, and part of the small particles could not be correctly measured, so the quality percentage of large-size and mid-size particles are larger than the real ones.

    Fig.11 Comparison of results of spherical plastic sands obtained by Malvern laser particle size analyzer and image measurement method

    Fig.12 Comparison of results of irregular shaped plastic sands obtained by Malvern laser particle size analyzer and image measurement method

    The comparison of the results of irregular shaped plastic sands obtained by the Malvern laser particle size analyzer and the image analysis measurement method is shown in Fig.12. In order to reduce the impact of small particles on the measurement results, the particles with sizes less than 0.06 mm are left out by using a screen. The gradation curves of the two methods are reasonably similar in shape, and the median diameters are, respectively, 0.33 mm and 0.34 mm for the laser particle size method and the image measurement method. The measured quality percentage of particles obtained by the image measurement method is larger than that of the laser particle size method when the particle size is larger than 0.2 mm, mainly due to the irregular shape. The imagemeasurement mainly concerns with the cross-section projection information of the sediment particles, and the measured size is larger than the actual one for the largest cross-section projection of the particle.

    For the laser particle size analyzer and the image analysis method, different definitions of equivalent diameter are used. For the former, the particle diameters are based on the scattering and diffraction theory[13-15], and the diameter is equivalent to the diameter of a ball with the same surface area as the sediment particles. For the latter, the diameter is equivalent to that of a circle with the same area as the maximal projection in the light propagation direction. Therefore, for irregular particles, the particle size measurement results show greater difference, while for spherical particles, the measurement results tend to be more consistent.

    4. Conclusions

    The grain size information is extracted by the threshold segmentation of the sediment particle images obtained by a camera. The pixels of each sediment particle are counted using the labeling method, and the diameter of a circle with an equivalent area of the sediment particle is calculated based on the equivalent area method. This system improves the sediment particle image segmentation by an extracting algorithm of the optimal threshold based on the gray histogram peak value, and makes the threshold transform more accurate. Theoretical calculations show that the minimum measurable size reaches 24.83 μm, enough to meet the requirement of the general sediment sample analysis, when the pixel size, the height of the camera, and the focal length are set at 4.4 μm, 400 mm, and 16 mm, respectively. This indicates that the sediment particle image analysis method is practical.

    The grading form of the image analysis method is basically similar to but steeper than that of the laser particle sizer. This system is limited by the minimum measurable size, providing a larger median particle size, which means that the deviation would be a little bit larger for small sediment particles. To improve the measurement accuracy for small sediment particles, a high-resolution CCD camera should be used.

    In the image measurements, only the particle information in the cross-section projection can be obtained, therefore, a multi-camera technology should be developed to reconstruct the three-dimensional shape of the sediment particles for further study of the physical properties of sediment particles.

    [1] SHAO Xue-jun, WANG Xing-kui. Introduction to river mechanics[M]. Beijing: Tsinghua University Press, 2005(in Chinese).

    [2] QIAN Ning, WAN Zhao-hui. Sediment dynamics[M]. Beijing, China: Science Press, 2003(in Chinese).

    [3] QIAN Jin, ZHENG Sha-sha and WANG Pei-fang et al. Experimental study on sediment resuspension in Taihu Lake under different hydrodynamic disturbances[J]. Journal of Hydrodynamics, 2011, 23(6): 826-833.

    [4] MISIRLI Z., ONER E. T. and KIRDAR B. Real imaging and size values of Saccharomyces cerevisiae cells with comparable contrast tuning to two environmental scanning electron microscopy modes[J]. Scanning, 2007, 29(1): 11-19.

    [5] APOSTOLOPOULOS G., TSINOPOULOS S. and DERMATAS E. Estimation of human red blood cells size using light scattering images[J]. Journal of Computational Methods in Sciences and Engineering, 2009, 9(1-2): 19-30.

    [6] SOUZA D. O. C., MENEGALLI F. C. Statistical study of particle size distribution and shape characterization[J]. Powder Technology, 2011, 214(1): 57-63

    [7] CARTER R. M., YAN Y. Measurement of particle shape using digital imaging techniques[J]. Journal of Physics: Conference Series, 2005, 15: 177-182.

    [8] CARTER R. M., YAN Y. and LEE P. On-line nonintrusive measurement of particle size distribution through digital imaging[J]. IEEE Transactions on Instrumentation and Measurement, 2006, 55(6): 2034-2038.

    [9] TYSMANS D., CLAEYS P. and DERIEMAEKER L. et al. Size and shape analysis of sedimentary grains by automated dynamic image analysis[J]. Particle and Particle Systems Characterization, 2007, 23(5): 381-387.

    [10] BOSCHETTO A., GIORDANO V. Powder sampling and characterization by digital image analysis[J]. Measurement: Journal of the International Measurement Confederation, 2012, 45(5): 1023-1038.

    [11] LI Dan-xun, XU Si-tao and YU Ming-zhong et al. Application of image processing to the analysis of gravel shape in Jialing River[J]. Journal of Sediment Research, 2000, (6): 6-9(in Chinese).

    [12] LIU Yun. Study on the sediment particle analysis system based on digital image process[D]. Doctoral Thesis, Nanjing, China: Hohai University, 2011(in Chinese).

    [13] ZOBECK T. M. Rapid soil particle size analyses using laser diffraction[J]. Applied Engineering in Agriculture, 2004, 20(5): 633-639.

    [14] VEIHELMANN B., KONERT M. and Van Der ZANDE W. J. Size distribution of mineral aerosol: Using lightscattering models in laser particle sizing[J]. Applied Optics, 2006, 45(23): 6022-6029.

    [15] PIERI L., BITTELLI M. and PISA P. R. Laser diffraction, transmission electron microscopy and image analysis to evaluate a bimodal Gaussian model for particle size distribution in soils[J]. Geoderma, 2006, 135: 118-132.

    10.1016/S1001-6058(13)60359-5

    * Project supported by the National Key Basic Research and Development of China (973 Program, Grant No. 2011CB403303), the China National Funds for Distinguished Young Scientists (Grant No. 51125034) and the National Natural Science Foundation of China (Grant Nos. 50909036, 50879019).

    Biography: CHEN Hong (1981-), Male, Ph. D. Candidate

    TANG Hong-wu,

    E-mail: hwtang@hhu.edu.cn

    猜你喜歡
    劉云陳紅王浩
    更正
    民警查逆行,他突然說自己殺了人
    方圓(2021年15期)2021-09-07 09:21:01
    詩與遠方
    陳紅作品
    喜龍?zhí)栵w船
    變臉的媽媽
    “想象魔法堡” 第41期來稿選登
    神秘的房東
    故事會(2016年8期)2016-05-06 09:20:04
    買走你的余生
    故事林(2012年4期)2012-05-14 17:29:57
    av在线蜜桃| 亚洲国产精品合色在线| 亚洲天堂国产精品一区在线| 国产精品久久久久久久久免 | 好男人在线观看高清免费视频| 久久久色成人| 亚洲五月婷婷丁香| 精品人妻熟女av久视频| 亚洲欧美日韩高清专用| 欧美xxxx性猛交bbbb| 极品教师在线视频| 夜夜躁狠狠躁天天躁| 熟妇人妻久久中文字幕3abv| 在线观看免费视频日本深夜| 亚洲av二区三区四区| 麻豆成人av在线观看| 欧美最新免费一区二区三区 | 国产蜜桃级精品一区二区三区| 久久国产乱子免费精品| 亚洲乱码一区二区免费版| 亚洲专区中文字幕在线| 2021天堂中文幕一二区在线观| 久久国产乱子伦精品免费另类| 久久人人精品亚洲av| 12—13女人毛片做爰片一| 日韩欧美国产一区二区入口| 色在线成人网| 国产黄a三级三级三级人| 免费搜索国产男女视频| 桃色一区二区三区在线观看| 真实男女啪啪啪动态图| 精品一区二区三区av网在线观看| 一本综合久久免费| 成人特级黄色片久久久久久久| 每晚都被弄得嗷嗷叫到高潮| 中文字幕av成人在线电影| 亚洲人成网站在线播放欧美日韩| 少妇熟女aⅴ在线视频| 网址你懂的国产日韩在线| 99久久无色码亚洲精品果冻| 毛片一级片免费看久久久久 | 国产欧美日韩一区二区三| 欧美一区二区国产精品久久精品| 亚洲性夜色夜夜综合| 我要看日韩黄色一级片| 少妇裸体淫交视频免费看高清| 国产三级黄色录像| 国产精品三级大全| 色播亚洲综合网| 日韩有码中文字幕| 久久精品国产亚洲av天美| 中文字幕高清在线视频| 日本熟妇午夜| 亚洲最大成人中文| 久久久久免费精品人妻一区二区| 成人精品一区二区免费| 国产激情偷乱视频一区二区| av在线老鸭窝| 国产综合懂色| 久久久久久久久久成人| 一本综合久久免费| 欧美+日韩+精品| 看黄色毛片网站| 欧洲精品卡2卡3卡4卡5卡区| 级片在线观看| 国产综合懂色| 日日夜夜操网爽| 91九色精品人成在线观看| 成人国产一区最新在线观看| 国产精品久久久久久亚洲av鲁大| 欧美bdsm另类| 亚洲熟妇中文字幕五十中出| 免费大片18禁| 久久九九热精品免费| 精品午夜福利视频在线观看一区| 级片在线观看| 搡老熟女国产l中国老女人| 国产精品一区二区性色av| 丁香欧美五月| 婷婷精品国产亚洲av在线| 亚洲自拍偷在线| 欧美黑人欧美精品刺激| 午夜福利免费观看在线| 久久久久国内视频| 蜜桃亚洲精品一区二区三区| 乱码一卡2卡4卡精品| 高潮久久久久久久久久久不卡| 精品国内亚洲2022精品成人| 国产精品一区二区三区四区免费观看 | 黄片小视频在线播放| 久久久久亚洲av毛片大全| 又黄又爽又免费观看的视频| 欧美三级亚洲精品| 亚洲精品一区av在线观看| 亚洲片人在线观看| 麻豆国产97在线/欧美| 亚洲aⅴ乱码一区二区在线播放| 麻豆成人午夜福利视频| 91麻豆av在线| 精品人妻一区二区三区麻豆 | 99久久久亚洲精品蜜臀av| 日韩中字成人| 琪琪午夜伦伦电影理论片6080| 91久久精品国产一区二区成人| 精品熟女少妇八av免费久了| 亚洲人成网站在线播| 亚洲av中文字字幕乱码综合| 亚洲av二区三区四区| 亚洲精品一区av在线观看| 男女床上黄色一级片免费看| 日本 欧美在线| 国产成人av教育| 国产综合懂色| 网址你懂的国产日韩在线| 精品人妻1区二区| 精华霜和精华液先用哪个| 真人做人爱边吃奶动态| 国产精品国产高清国产av| 久久久久亚洲av毛片大全| 免费av毛片视频| 中文字幕熟女人妻在线| 性插视频无遮挡在线免费观看| 97超级碰碰碰精品色视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 级片在线观看| 欧美日本视频| 激情在线观看视频在线高清| 亚洲 欧美 日韩 在线 免费| 色在线成人网| 有码 亚洲区| 在线观看午夜福利视频| 久久久国产成人精品二区| 成人美女网站在线观看视频| 日本熟妇午夜| 日韩精品青青久久久久久| 99热只有精品国产| 少妇人妻一区二区三区视频| 两性午夜刺激爽爽歪歪视频在线观看| 美女免费视频网站| 九九久久精品国产亚洲av麻豆| 18禁黄网站禁片免费观看直播| 亚洲精品在线美女| 亚洲天堂国产精品一区在线| 精品欧美国产一区二区三| 日本免费a在线| 欧美精品国产亚洲| 精品一区二区免费观看| 人人妻,人人澡人人爽秒播| 国产黄片美女视频| 老司机深夜福利视频在线观看| 人妻丰满熟妇av一区二区三区| 午夜精品久久久久久毛片777| 国产免费av片在线观看野外av| bbb黄色大片| 中文资源天堂在线| 亚洲精华国产精华精| 一进一出好大好爽视频| 人人妻,人人澡人人爽秒播| 999久久久精品免费观看国产| 99国产精品一区二区三区| 国产精品久久电影中文字幕| 两人在一起打扑克的视频| 丝袜美腿在线中文| 中文字幕人成人乱码亚洲影| 国产精品乱码一区二三区的特点| 国产成人a区在线观看| 欧美三级亚洲精品| 亚洲av.av天堂| 久久午夜亚洲精品久久| 伊人久久精品亚洲午夜| 欧美黄色片欧美黄色片| 色视频www国产| 国产精品不卡视频一区二区 | 在线观看舔阴道视频| 怎么达到女性高潮| 在线播放国产精品三级| 国产精品久久电影中文字幕| 国产人妻一区二区三区在| 18禁裸乳无遮挡免费网站照片| 在线看三级毛片| 日本撒尿小便嘘嘘汇集6| 天堂影院成人在线观看| 日韩欧美在线乱码| 在线a可以看的网站| 一本综合久久免费| 欧美日韩瑟瑟在线播放| 国产成+人综合+亚洲专区| 免费一级毛片在线播放高清视频| xxxwww97欧美| 国产亚洲欧美在线一区二区| 99精品在免费线老司机午夜| 国产av一区在线观看免费| 久久国产精品人妻蜜桃| 国产色爽女视频免费观看| 欧美高清成人免费视频www| 757午夜福利合集在线观看| 国产熟女xx| 老熟妇乱子伦视频在线观看| ponron亚洲| 国产主播在线观看一区二区| 一进一出好大好爽视频| 国内少妇人妻偷人精品xxx网站| 麻豆一二三区av精品| 亚洲人成网站在线播| 亚洲欧美精品综合久久99| 最近在线观看免费完整版| 99热6这里只有精品| 非洲黑人性xxxx精品又粗又长| 女人被狂操c到高潮| 人妻久久中文字幕网| 国产精品久久电影中文字幕| 高清在线国产一区| av专区在线播放| 国产综合懂色| 一级av片app| 搡老妇女老女人老熟妇| 国产一区二区三区在线臀色熟女| 久久久精品欧美日韩精品| 国产三级黄色录像| 成人特级av手机在线观看| 国产黄a三级三级三级人| 欧美精品国产亚洲| 亚洲国产欧洲综合997久久,| 亚洲精品色激情综合| 国产欧美日韩精品亚洲av| 国产精品三级大全| 亚洲欧美清纯卡通| 特级一级黄色大片| 变态另类丝袜制服| 久久精品国产亚洲av天美| 精品人妻偷拍中文字幕| 一区二区三区免费毛片| 成人永久免费在线观看视频| 91麻豆精品激情在线观看国产| 午夜激情欧美在线| 永久网站在线| 日韩欧美免费精品| 窝窝影院91人妻| 欧美在线黄色| 亚洲人与动物交配视频| 欧美区成人在线视频| 丰满乱子伦码专区| 日本一本二区三区精品| 久久亚洲精品不卡| 97超级碰碰碰精品色视频在线观看| 免费av不卡在线播放| 国产黄a三级三级三级人| 一卡2卡三卡四卡精品乱码亚洲| 永久网站在线| 中亚洲国语对白在线视频| 麻豆一二三区av精品| 国产伦精品一区二区三区视频9| 简卡轻食公司| 国产老妇女一区| 久久久精品欧美日韩精品| 中文在线观看免费www的网站| 国产一区二区三区视频了| 1000部很黄的大片| 婷婷丁香在线五月| 草草在线视频免费看| 日韩欧美 国产精品| 欧美在线黄色| 夜夜爽天天搞| 亚洲精品乱码久久久v下载方式| 一个人看的www免费观看视频| 日本 av在线| 久久天躁狠狠躁夜夜2o2o| 69av精品久久久久久| 亚洲男人的天堂狠狠| 欧美日韩综合久久久久久 | 欧美日韩瑟瑟在线播放| 他把我摸到了高潮在线观看| 最近在线观看免费完整版| 赤兔流量卡办理| 婷婷精品国产亚洲av在线| 嫩草影院精品99| 91麻豆av在线| 一个人观看的视频www高清免费观看| 欧洲精品卡2卡3卡4卡5卡区| 我的老师免费观看完整版| 露出奶头的视频| 日本黄色片子视频| 午夜两性在线视频| 久久久精品欧美日韩精品| 国产av一区在线观看免费| 亚洲真实伦在线观看| 欧美zozozo另类| 国产色爽女视频免费观看| 男人舔奶头视频| 一区二区三区四区激情视频 | 老熟妇乱子伦视频在线观看| 欧美一区二区国产精品久久精品| 精品99又大又爽又粗少妇毛片 | 深夜a级毛片| 欧美成人一区二区免费高清观看| 精品一区二区三区av网在线观看| 人人妻,人人澡人人爽秒播| 日日干狠狠操夜夜爽| 成人亚洲精品av一区二区| 少妇熟女aⅴ在线视频| 久久精品久久久久久噜噜老黄 | 五月玫瑰六月丁香| 亚洲成人久久爱视频| 欧美+日韩+精品| 最新中文字幕久久久久| 欧美乱妇无乱码| 亚洲成av人片在线播放无| www日本黄色视频网| 桃红色精品国产亚洲av| 成人精品一区二区免费| 国产不卡一卡二| 午夜激情欧美在线| 此物有八面人人有两片| 51国产日韩欧美| 日本黄大片高清| 男女视频在线观看网站免费| 啦啦啦观看免费观看视频高清| 久久精品91蜜桃| 人妻制服诱惑在线中文字幕| 91av网一区二区| 亚洲成av人片免费观看| 午夜日韩欧美国产| 嫩草影视91久久| 亚洲av一区综合| 99久久精品国产亚洲精品| 欧美激情国产日韩精品一区| 嫩草影视91久久| 91在线精品国自产拍蜜月| 十八禁国产超污无遮挡网站| 久久热精品热| 少妇裸体淫交视频免费看高清| 精品久久久久久久久久久久久| 人妻夜夜爽99麻豆av| 亚洲欧美清纯卡通| 亚洲第一区二区三区不卡| 中国美女看黄片| 国产精品乱码一区二三区的特点| 欧美成人性av电影在线观看| 成人国产综合亚洲| 天天躁日日操中文字幕| 网址你懂的国产日韩在线| 日日干狠狠操夜夜爽| 老司机深夜福利视频在线观看| 极品教师在线免费播放| 内射极品少妇av片p| 超碰av人人做人人爽久久| 国产色婷婷99| 日韩中字成人| 一夜夜www| 99久久成人亚洲精品观看| 日韩精品中文字幕看吧| 国产精品久久久久久精品电影| 亚洲人成电影免费在线| 久久久久久久久久成人| 宅男免费午夜| 99精品久久久久人妻精品| 男女那种视频在线观看| 深夜a级毛片| 欧美丝袜亚洲另类 | 精品人妻偷拍中文字幕| x7x7x7水蜜桃| 最好的美女福利视频网| 色视频www国产| 日本与韩国留学比较| 日本黄色片子视频| 亚洲精品成人久久久久久| 成人鲁丝片一二三区免费| 精品人妻1区二区| 日韩 亚洲 欧美在线| 日本黄大片高清| 精品欧美国产一区二区三| 精品一区二区三区人妻视频| 国产精品美女特级片免费视频播放器| 国产精品影院久久| 欧美丝袜亚洲另类 | а√天堂www在线а√下载| 亚洲,欧美,日韩| 日本免费a在线| 免费看光身美女| a级毛片a级免费在线| 色av中文字幕| 国语自产精品视频在线第100页| 久久精品国产自在天天线| 给我免费播放毛片高清在线观看| 特级一级黄色大片| 亚洲av免费在线观看| 亚洲国产色片| 免费av观看视频| 日韩欧美一区二区三区在线观看| 最近视频中文字幕2019在线8| 国产精品久久久久久久久免 | 在线观看免费视频日本深夜| 一个人免费在线观看的高清视频| 又爽又黄a免费视频| 欧美日本亚洲视频在线播放| 最后的刺客免费高清国语| 91狼人影院| 日日摸夜夜添夜夜添小说| 小说图片视频综合网站| 欧美精品国产亚洲| 国产精品亚洲av一区麻豆| 国产午夜福利久久久久久| 两个人视频免费观看高清| av专区在线播放| 久久婷婷人人爽人人干人人爱| 99久久精品国产亚洲精品| 亚洲第一区二区三区不卡| 女人十人毛片免费观看3o分钟| 久久久久久久久中文| 乱人视频在线观看| 欧美绝顶高潮抽搐喷水| 亚洲aⅴ乱码一区二区在线播放| 亚洲av二区三区四区| 久久天躁狠狠躁夜夜2o2o| 麻豆成人av在线观看| 午夜亚洲福利在线播放| 久久人妻av系列| 精品人妻视频免费看| 亚洲国产色片| 亚洲人成网站在线播放欧美日韩| 日本a在线网址| 可以在线观看毛片的网站| 亚洲av美国av| 久久久国产成人精品二区| 两个人的视频大全免费| 亚洲五月天丁香| 色5月婷婷丁香| 麻豆国产97在线/欧美| 欧美国产日韩亚洲一区| 此物有八面人人有两片| 国产亚洲欧美98| 久久精品国产自在天天线| 99久久成人亚洲精品观看| 亚洲精品色激情综合| 国产精品野战在线观看| 国产黄色小视频在线观看| 黄色配什么色好看| 97热精品久久久久久| 久久精品91蜜桃| 熟女人妻精品中文字幕| 亚洲欧美日韩无卡精品| 内地一区二区视频在线| 亚洲国产欧洲综合997久久,| 国产精品不卡视频一区二区 | 欧美日本视频| 久久久久久久精品吃奶| 亚洲电影在线观看av| 1024手机看黄色片| 亚洲人成网站在线播放欧美日韩| 老鸭窝网址在线观看| 日本熟妇午夜| www.999成人在线观看| 九色国产91popny在线| 久久99热这里只有精品18| 天天躁日日操中文字幕| 午夜久久久久精精品| 男女做爰动态图高潮gif福利片| 又爽又黄无遮挡网站| 欧美日本亚洲视频在线播放| 亚洲成人久久爱视频| 国产免费一级a男人的天堂| 性色avwww在线观看| 成年版毛片免费区| 99国产精品一区二区蜜桃av| 亚洲欧美日韩东京热| 欧美性猛交╳xxx乱大交人| 久久久久久久久中文| 精品国产三级普通话版| 他把我摸到了高潮在线观看| 久久精品国产99精品国产亚洲性色| 此物有八面人人有两片| 美女 人体艺术 gogo| 免费av不卡在线播放| 亚洲最大成人av| 亚洲欧美日韩高清在线视频| avwww免费| 国产视频内射| 日本黄色视频三级网站网址| 夜夜爽天天搞| 国产精品一区二区三区四区免费观看 | 亚洲真实伦在线观看| 1024手机看黄色片| 亚洲中文字幕日韩| 天天躁日日操中文字幕| 国产精品1区2区在线观看.| 久久九九热精品免费| 精品人妻1区二区| 欧洲精品卡2卡3卡4卡5卡区| 有码 亚洲区| 特级一级黄色大片| 国产精华一区二区三区| 婷婷精品国产亚洲av在线| 亚洲美女搞黄在线观看 | 九九久久精品国产亚洲av麻豆| 啦啦啦观看免费观看视频高清| 18禁在线播放成人免费| 国产视频一区二区在线看| 亚洲精品久久国产高清桃花| 国语自产精品视频在线第100页| 亚洲精华国产精华精| 在线十欧美十亚洲十日本专区| 欧美性猛交╳xxx乱大交人| a级一级毛片免费在线观看| 久久人人爽人人爽人人片va | 9191精品国产免费久久| 国产欧美日韩精品亚洲av| 久久人妻av系列| 国产男靠女视频免费网站| 久99久视频精品免费| 人妻夜夜爽99麻豆av| 欧美成狂野欧美在线观看| 国产精品国产高清国产av| 欧美成人性av电影在线观看| 女同久久另类99精品国产91| 老司机午夜十八禁免费视频| 亚洲美女视频黄频| 99久久精品一区二区三区| 久久久久久久精品吃奶| 日韩欧美 国产精品| 欧美黄色淫秽网站| 美女免费视频网站| 国产精品久久电影中文字幕| 国产蜜桃级精品一区二区三区| 日韩欧美国产在线观看| 他把我摸到了高潮在线观看| 国产成人啪精品午夜网站| 免费电影在线观看免费观看| 精品人妻一区二区三区麻豆 | 免费av不卡在线播放| 久久6这里有精品| 国产精品98久久久久久宅男小说| 亚洲五月天丁香| av福利片在线观看| 中国美女看黄片| 我的老师免费观看完整版| 久久久国产成人免费| 欧美+亚洲+日韩+国产| 久久久久性生活片| 神马国产精品三级电影在线观看| 好男人在线观看高清免费视频| 国产白丝娇喘喷水9色精品| 欧美精品国产亚洲| 成年女人毛片免费观看观看9| 婷婷精品国产亚洲av| 亚洲内射少妇av| 丰满人妻熟妇乱又伦精品不卡| 麻豆久久精品国产亚洲av| 久久天躁狠狠躁夜夜2o2o| 老司机深夜福利视频在线观看| 国产精品爽爽va在线观看网站| 亚洲在线观看片| 亚洲片人在线观看| 首页视频小说图片口味搜索| 国产高潮美女av| 高清毛片免费观看视频网站| 欧美绝顶高潮抽搐喷水| 日韩欧美 国产精品| 精品一区二区三区人妻视频| 99热精品在线国产| 成人av一区二区三区在线看| 观看美女的网站| 啦啦啦观看免费观看视频高清| 国产一区二区在线观看日韩| 很黄的视频免费| 人妻夜夜爽99麻豆av| 亚洲专区中文字幕在线| 毛片女人毛片| 久久热精品热| 美女高潮的动态| 嫩草影视91久久| 美女高潮的动态| 成人国产综合亚洲| 精品国产三级普通话版| 欧美中文日本在线观看视频| 亚洲精品亚洲一区二区| 在线免费观看不下载黄p国产 | 亚洲成av人片在线播放无| 午夜亚洲福利在线播放| 精品福利观看| 成年人黄色毛片网站| 又粗又爽又猛毛片免费看| 成年免费大片在线观看| eeuss影院久久| 欧美国产日韩亚洲一区| 美女黄网站色视频| 国产伦精品一区二区三区四那| 深爱激情五月婷婷| 国产白丝娇喘喷水9色精品| 男人的好看免费观看在线视频| 成人毛片a级毛片在线播放| 亚洲成a人片在线一区二区| 国产不卡一卡二| 日本一本二区三区精品| 综合色av麻豆| 757午夜福利合集在线观看| 十八禁人妻一区二区| 国产精品综合久久久久久久免费| 一级黄片播放器| 亚洲成人中文字幕在线播放| www.熟女人妻精品国产| 99在线人妻在线中文字幕| 亚洲国产日韩欧美精品在线观看| 天美传媒精品一区二区| 免费av不卡在线播放| 午夜精品久久久久久毛片777| 国产男靠女视频免费网站| 老司机福利观看| 久久久久久国产a免费观看| 在线天堂最新版资源| 99久国产av精品| 一个人免费在线观看的高清视频| 亚洲第一区二区三区不卡| 国产精品久久电影中文字幕| 美女xxoo啪啪120秒动态图 | 熟女人妻精品中文字幕|