• <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
    久久精品国产综合久久久| 亚洲美女黄片视频| 啦啦啦 在线观看视频| 久久99一区二区三区| 久久性视频一级片| 午夜成年电影在线免费观看| 水蜜桃什么品种好| 精品一品国产午夜福利视频| 一级片免费观看大全| 国产极品粉嫩免费观看在线| 一级毛片女人18水好多| 在线观看66精品国产| 亚洲第一av免费看| 国产伦人伦偷精品视频| 一区二区av电影网| 大片免费播放器 马上看| 亚洲 欧美一区二区三区| 亚洲自偷自拍图片 自拍| 国产深夜福利视频在线观看| 精品久久久久久久毛片微露脸| 波多野结衣一区麻豆| 高清在线国产一区| 亚洲午夜精品一区,二区,三区| 国产亚洲精品第一综合不卡| 成在线人永久免费视频| 亚洲成国产人片在线观看| 成人精品一区二区免费| 欧美日韩亚洲高清精品| 成人免费观看视频高清| 欧美日韩福利视频一区二区| 色播在线永久视频| 黑人巨大精品欧美一区二区蜜桃| 亚洲熟妇熟女久久| 欧美日韩福利视频一区二区| 在线av久久热| 欧美另类亚洲清纯唯美| 精品少妇久久久久久888优播| 色综合婷婷激情| 国产精品一区二区精品视频观看| 丝袜在线中文字幕| 少妇粗大呻吟视频| 国产日韩欧美在线精品| 欧美黄色片欧美黄色片| 美女高潮到喷水免费观看| 国产一区二区三区在线臀色熟女 | 一区二区三区乱码不卡18| 亚洲熟女毛片儿| 国产精品麻豆人妻色哟哟久久| 他把我摸到了高潮在线观看 | netflix在线观看网站| 国产精品亚洲一级av第二区| 人人妻,人人澡人人爽秒播| 他把我摸到了高潮在线观看 | 亚洲一区中文字幕在线| 水蜜桃什么品种好| 亚洲美女黄片视频| 国产成人精品在线电影| 在线看a的网站| 黄色片一级片一级黄色片| 久久影院123| 两性午夜刺激爽爽歪歪视频在线观看 | 在线观看免费视频日本深夜| 一区二区三区激情视频| 精品一区二区三区视频在线观看免费 | 国产深夜福利视频在线观看| 国产又色又爽无遮挡免费看| 老汉色av国产亚洲站长工具| 露出奶头的视频| 亚洲一区二区三区欧美精品| 99久久国产精品久久久| 高清av免费在线| 一级毛片精品| 下体分泌物呈黄色| 纯流量卡能插随身wifi吗| 日韩视频在线欧美| 午夜福利视频在线观看免费| 99久久人妻综合| 热99久久久久精品小说推荐| 美女午夜性视频免费| 国产在线观看jvid| 在线看a的网站| 欧美亚洲 丝袜 人妻 在线| 美女扒开内裤让男人捅视频| 一级黄色大片毛片| 老汉色av国产亚洲站长工具| 欧美日韩视频精品一区| 99久久99久久久精品蜜桃| 人人妻,人人澡人人爽秒播| 99国产综合亚洲精品| 99久久国产精品久久久| 五月天丁香电影| 国产又色又爽无遮挡免费看| 欧美成狂野欧美在线观看| 国产精品久久电影中文字幕 | 操出白浆在线播放| 两人在一起打扑克的视频| 精品国产一区二区久久| 亚洲人成电影免费在线| 夫妻午夜视频| 亚洲中文日韩欧美视频| 一级黄色大片毛片| 亚洲国产毛片av蜜桃av| 亚洲人成77777在线视频| 亚洲国产精品一区二区三区在线| 欧美成狂野欧美在线观看| 老司机福利观看| 日日爽夜夜爽网站| 免费日韩欧美在线观看| 在线观看www视频免费| 一个人免费看片子| 捣出白浆h1v1| 国内毛片毛片毛片毛片毛片| 黑人猛操日本美女一级片| 精品少妇一区二区三区视频日本电影| 国产三级黄色录像| 中文字幕色久视频| 色老头精品视频在线观看| 欧美黑人精品巨大| 久久午夜亚洲精品久久| 丰满少妇做爰视频| 一进一出好大好爽视频| 亚洲天堂av无毛| 精品高清国产在线一区| 国产主播在线观看一区二区| 欧美亚洲 丝袜 人妻 在线| 精品国产乱子伦一区二区三区| 国产高清视频在线播放一区| 国产免费福利视频在线观看| 精品一区二区三区av网在线观看 | 男女免费视频国产| 久久人妻av系列| 欧美 日韩 精品 国产| kizo精华| 久9热在线精品视频| 在线观看免费高清a一片| 欧美另类亚洲清纯唯美| 女人精品久久久久毛片| 久久精品亚洲av国产电影网| 亚洲国产精品一区二区三区在线| 国产精品一区二区在线观看99| 亚洲熟妇熟女久久| 亚洲国产看品久久| 首页视频小说图片口味搜索| 亚洲七黄色美女视频| 欧美日韩亚洲国产一区二区在线观看 | 狂野欧美激情性xxxx| 国产精品国产高清国产av | 国产精品 国内视频| 国产免费福利视频在线观看| 国产aⅴ精品一区二区三区波| 国产精品久久久久久精品古装| 国产精品九九99| 极品人妻少妇av视频| 精品久久久精品久久久| a级片在线免费高清观看视频| 中国美女看黄片| 精品少妇久久久久久888优播| 国产精品影院久久| 一边摸一边抽搐一进一小说 | 飞空精品影院首页| 女人久久www免费人成看片| 亚洲成a人片在线一区二区| 99国产精品一区二区三区| h视频一区二区三区| 国产成人影院久久av| 国产99久久九九免费精品| 十八禁高潮呻吟视频| 51午夜福利影视在线观看| 多毛熟女@视频| 国产精品久久久久久精品古装| 亚洲人成77777在线视频| 亚洲精品乱久久久久久| 亚洲精品国产区一区二| 亚洲国产精品一区二区三区在线| 91老司机精品| 嫩草影视91久久| 精品免费久久久久久久清纯 | 一级黄色大片毛片| 99国产综合亚洲精品| 婷婷丁香在线五月| 另类精品久久| 性色av乱码一区二区三区2| 国产高清videossex| 12—13女人毛片做爰片一| 国产av精品麻豆| 美女扒开内裤让男人捅视频| 欧美+亚洲+日韩+国产| 免费观看av网站的网址| 50天的宝宝边吃奶边哭怎么回事| 丁香欧美五月| 中文欧美无线码| 99久久99久久久精品蜜桃| 国产淫语在线视频| 别揉我奶头~嗯~啊~动态视频| 又黄又粗又硬又大视频| 无遮挡黄片免费观看| 男女边摸边吃奶| av国产精品久久久久影院| 国产亚洲精品第一综合不卡| 一边摸一边抽搐一进一小说 | 嫩草影视91久久| 亚洲一码二码三码区别大吗| 人人妻人人添人人爽欧美一区卜| 新久久久久国产一级毛片| 电影成人av| 亚洲欧美日韩高清在线视频 | 熟女少妇亚洲综合色aaa.| 国产av国产精品国产| 精品少妇一区二区三区视频日本电影| 1024香蕉在线观看| 亚洲精品自拍成人| 婷婷丁香在线五月| 97人妻天天添夜夜摸| 精品亚洲乱码少妇综合久久| 国产三级黄色录像| aaaaa片日本免费| 色视频在线一区二区三区| 国产精品99久久99久久久不卡| 国产野战对白在线观看| 国产成人精品在线电影| 国产欧美日韩综合在线一区二区| 日本av手机在线免费观看| 变态另类成人亚洲欧美熟女 | 精品免费久久久久久久清纯 | 高清在线国产一区| 国产无遮挡羞羞视频在线观看| 免费高清在线观看日韩| 久久午夜综合久久蜜桃| 久久久国产欧美日韩av| a在线观看视频网站| 丰满少妇做爰视频| 男女床上黄色一级片免费看| 国产亚洲午夜精品一区二区久久| 亚洲精品在线观看二区| 久久性视频一级片| 一二三四在线观看免费中文在| 久久久久精品人妻al黑| 十分钟在线观看高清视频www| 美女国产高潮福利片在线看| 亚洲中文av在线| 久久久精品94久久精品| 新久久久久国产一级毛片| 国产黄频视频在线观看| 老司机在亚洲福利影院| 肉色欧美久久久久久久蜜桃| 久久青草综合色| 免费看a级黄色片| 欧美日韩成人在线一区二区| 欧美人与性动交α欧美软件| 捣出白浆h1v1| 国产精品 欧美亚洲| 免费在线观看影片大全网站| 欧美日韩成人在线一区二区| 99九九在线精品视频| 欧美变态另类bdsm刘玥| 大码成人一级视频| 99国产极品粉嫩在线观看| 免费人妻精品一区二区三区视频| 亚洲欧美一区二区三区久久| 一二三四社区在线视频社区8| av网站在线播放免费| 国产无遮挡羞羞视频在线观看| 男女免费视频国产| 一本久久精品| 视频区欧美日本亚洲| 久久久久精品国产欧美久久久| 97在线人人人人妻| 亚洲人成电影免费在线| 欧美亚洲日本最大视频资源| 中文字幕人妻丝袜制服| 欧美黄色淫秽网站| 97人妻天天添夜夜摸| 久久精品国产亚洲av香蕉五月 | 美女主播在线视频| 99九九在线精品视频| 国产精品麻豆人妻色哟哟久久| 亚洲av欧美aⅴ国产| 丰满迷人的少妇在线观看| 在线亚洲精品国产二区图片欧美| 欧美激情极品国产一区二区三区| 国产国语露脸激情在线看| 成人精品一区二区免费| 国产精品麻豆人妻色哟哟久久| 男女边摸边吃奶| 怎么达到女性高潮| av天堂在线播放| 肉色欧美久久久久久久蜜桃| 国产成人啪精品午夜网站| 一区二区三区激情视频| 亚洲av日韩在线播放| av国产精品久久久久影院| 久久久精品免费免费高清| 桃花免费在线播放| 另类精品久久| 国产黄色免费在线视频| 在线观看66精品国产| 日本a在线网址| 精品国产国语对白av| 国产在线观看jvid| 少妇的丰满在线观看| 国产精品久久久久久精品古装| 国产精品久久久久久人妻精品电影 | 久久久欧美国产精品| 在线av久久热| 色在线成人网| 成人国语在线视频| 久久人人97超碰香蕉20202| 亚洲av国产av综合av卡| 水蜜桃什么品种好| 亚洲国产欧美一区二区综合| 欧美国产精品va在线观看不卡| 岛国毛片在线播放| 中文字幕av电影在线播放| 日韩欧美一区视频在线观看| 一本—道久久a久久精品蜜桃钙片| 好男人电影高清在线观看| 大香蕉久久成人网| 国产精品98久久久久久宅男小说| 国产成人精品无人区| 精品卡一卡二卡四卡免费| 人妻一区二区av| www.精华液| 国产日韩欧美亚洲二区| 菩萨蛮人人尽说江南好唐韦庄| 国产一区二区三区视频了| a在线观看视频网站| 97人妻天天添夜夜摸| 俄罗斯特黄特色一大片| av电影中文网址| 黄色视频不卡| 中文字幕av电影在线播放| 亚洲成人国产一区在线观看| 成人18禁高潮啪啪吃奶动态图| 男女无遮挡免费网站观看| 狂野欧美激情性xxxx| 女人高潮潮喷娇喘18禁视频| 日韩有码中文字幕| 日韩免费av在线播放| 亚洲熟女精品中文字幕| 欧美 日韩 精品 国产| 成年人黄色毛片网站| 一区二区三区激情视频| 久久久久国内视频| 成人特级黄色片久久久久久久 | 成人永久免费在线观看视频 | 下体分泌物呈黄色| 日韩人妻精品一区2区三区| 日韩免费高清中文字幕av| 91成年电影在线观看| 久久 成人 亚洲| 一二三四社区在线视频社区8| 两个人看的免费小视频| 国产不卡av网站在线观看| 亚洲综合色网址| 精品亚洲乱码少妇综合久久| 夜夜夜夜夜久久久久| 日韩欧美国产一区二区入口| 性少妇av在线| 午夜免费成人在线视频| 免费不卡黄色视频| 亚洲美女黄片视频| 777米奇影视久久| 天天影视国产精品| 久久久久国产一级毛片高清牌| 精品亚洲乱码少妇综合久久| 成人国产av品久久久| 日本一区二区免费在线视频| 91字幕亚洲| 亚洲人成电影免费在线| 亚洲av成人不卡在线观看播放网| 国产精品秋霞免费鲁丝片| 午夜精品国产一区二区电影| 日本wwww免费看| 国产深夜福利视频在线观看| 桃花免费在线播放| 亚洲成人免费av在线播放| 色婷婷av一区二区三区视频| tube8黄色片| 老司机午夜福利在线观看视频 | 国产欧美日韩综合在线一区二区| 亚洲国产成人一精品久久久| 中文字幕人妻熟女乱码| 成人av一区二区三区在线看| 国产精品美女特级片免费视频播放器 | 在线观看免费视频网站a站| 大码成人一级视频| 久热这里只有精品99| 在线播放国产精品三级| 日本精品一区二区三区蜜桃| 国产成人免费无遮挡视频| videosex国产| 午夜日韩欧美国产| 亚洲人成77777在线视频| 99国产精品99久久久久| 成在线人永久免费视频| 欧美精品啪啪一区二区三区| 中亚洲国语对白在线视频| 一级a爱视频在线免费观看| 男女之事视频高清在线观看| 欧美日韩精品网址| 一二三四社区在线视频社区8| 久久99热这里只频精品6学生| 亚洲熟妇熟女久久| 黑人操中国人逼视频| 久久人妻福利社区极品人妻图片| 在线十欧美十亚洲十日本专区| 人成视频在线观看免费观看| 97人妻天天添夜夜摸| 亚洲av欧美aⅴ国产| 侵犯人妻中文字幕一二三四区| 香蕉丝袜av| 天天躁日日躁夜夜躁夜夜| 99国产综合亚洲精品| 搡老熟女国产l中国老女人| 色在线成人网| 国产激情久久老熟女| 天天影视国产精品| 成人三级做爰电影| 80岁老熟妇乱子伦牲交| 9色porny在线观看| 欧美日本中文国产一区发布| 美女福利国产在线| 精品久久久久久久毛片微露脸| 国产亚洲精品一区二区www | 国产精品秋霞免费鲁丝片| 国产精品影院久久| 一级片'在线观看视频| 人妻 亚洲 视频| 丰满迷人的少妇在线观看| 一进一出好大好爽视频| 精品第一国产精品| 一个人免费看片子| 国产成人啪精品午夜网站| 另类亚洲欧美激情| 精品国产乱子伦一区二区三区| 欧美黑人精品巨大| √禁漫天堂资源中文www| 一区福利在线观看| 大码成人一级视频| 91麻豆精品激情在线观看国产 | 十八禁高潮呻吟视频| 这个男人来自地球电影免费观看| 妹子高潮喷水视频| 热99re8久久精品国产| 新久久久久国产一级毛片| 亚洲综合色网址| 精品一品国产午夜福利视频| 国产成人av激情在线播放| 成年人免费黄色播放视频| 丁香六月欧美| 欧美激情高清一区二区三区| 水蜜桃什么品种好| 欧美人与性动交α欧美精品济南到| 91字幕亚洲| 久久久久视频综合| 国产91精品成人一区二区三区 | 精品国产乱码久久久久久男人| 无限看片的www在线观看| 亚洲一码二码三码区别大吗| 不卡av一区二区三区| 亚洲国产精品一区二区三区在线| 一二三四在线观看免费中文在| 欧美亚洲 丝袜 人妻 在线| 午夜视频精品福利| 亚洲精品乱久久久久久| 国产av一区二区精品久久| 久久精品国产99精品国产亚洲性色 | 中文字幕另类日韩欧美亚洲嫩草| 夜夜骑夜夜射夜夜干| 啪啪无遮挡十八禁网站| 精品一区二区三区视频在线观看免费 | 国产一区有黄有色的免费视频| 夫妻午夜视频| 19禁男女啪啪无遮挡网站| 国产成+人综合+亚洲专区| 夜夜爽天天搞| 久久久久久人人人人人| 国产极品粉嫩免费观看在线| 一个人免费看片子| 搡老乐熟女国产| 成年人免费黄色播放视频| 精品国产国语对白av| 欧美日韩亚洲综合一区二区三区_| 97人妻天天添夜夜摸| 老鸭窝网址在线观看| 国产精品免费视频内射| 国产欧美日韩一区二区精品| 色在线成人网| 久久影院123| 露出奶头的视频| 老熟妇乱子伦视频在线观看| 亚洲色图 男人天堂 中文字幕| 九色亚洲精品在线播放| 久久久久国内视频| 黑人巨大精品欧美一区二区蜜桃| av又黄又爽大尺度在线免费看| 亚洲成人国产一区在线观看| 9191精品国产免费久久| 亚洲国产精品一区二区三区在线| 日韩一区二区三区影片| 亚洲国产av新网站| 亚洲精品久久成人aⅴ小说| 欧美另类亚洲清纯唯美| 99re6热这里在线精品视频| 一本久久精品| 一区二区三区精品91| 人人妻,人人澡人人爽秒播| 亚洲欧洲日产国产| 黄色视频不卡| 国产精品国产高清国产av | 色尼玛亚洲综合影院| 久久精品aⅴ一区二区三区四区| 欧美久久黑人一区二区| 俄罗斯特黄特色一大片| 九色亚洲精品在线播放| 最近最新中文字幕大全免费视频| 老熟妇乱子伦视频在线观看| 中文欧美无线码| 精品亚洲成a人片在线观看| 日韩中文字幕视频在线看片| 国产免费av片在线观看野外av| 脱女人内裤的视频| 国产一区二区激情短视频| 日日摸夜夜添夜夜添小说| 肉色欧美久久久久久久蜜桃| 丰满饥渴人妻一区二区三| 1024视频免费在线观看| a级毛片黄视频| 一级,二级,三级黄色视频| 飞空精品影院首页| 两人在一起打扑克的视频| 他把我摸到了高潮在线观看 | 亚洲av欧美aⅴ国产| 亚洲国产成人一精品久久久| 国产男女内射视频| 在线天堂中文资源库| 无限看片的www在线观看| 国产黄频视频在线观看| 一个人免费看片子| 黄色怎么调成土黄色| 国产免费福利视频在线观看| 黄色 视频免费看| av片东京热男人的天堂| 久热这里只有精品99| 美女福利国产在线| 在线看a的网站| 欧美 日韩 精品 国产| 大片电影免费在线观看免费| 久久久水蜜桃国产精品网| 国产精品国产高清国产av | xxxhd国产人妻xxx| 在线天堂中文资源库| 99国产精品一区二区三区| 亚洲色图 男人天堂 中文字幕| 亚洲成人免费电影在线观看| 欧美乱妇无乱码| 美女视频免费永久观看网站| 日韩有码中文字幕| 亚洲少妇的诱惑av| 亚洲成a人片在线一区二区| 免费少妇av软件| 国产成人啪精品午夜网站| 欧美黄色片欧美黄色片| 我的亚洲天堂| 大香蕉久久网| 在线观看舔阴道视频| 大陆偷拍与自拍| 免费日韩欧美在线观看| 男女之事视频高清在线观看| 亚洲一区中文字幕在线| 母亲3免费完整高清在线观看| 又紧又爽又黄一区二区| 欧美日韩成人在线一区二区| a级毛片黄视频| 久久免费观看电影| 国产精品1区2区在线观看. | videosex国产| 免费一级毛片在线播放高清视频 | 热99re8久久精品国产| 国产野战对白在线观看| 亚洲伊人久久精品综合| 国产真人三级小视频在线观看| 色尼玛亚洲综合影院| 天堂动漫精品| 久久精品国产综合久久久| 欧美日韩中文字幕国产精品一区二区三区 | 菩萨蛮人人尽说江南好唐韦庄| 国产野战对白在线观看| 老司机午夜福利在线观看视频 | 别揉我奶头~嗯~啊~动态视频| 亚洲精品久久成人aⅴ小说| 亚洲色图综合在线观看| 久久天躁狠狠躁夜夜2o2o| 亚洲成人免费av在线播放| 久久青草综合色| 一级,二级,三级黄色视频| 三上悠亚av全集在线观看| 美女扒开内裤让男人捅视频| 精品久久久久久久毛片微露脸| 精品国产超薄肉色丝袜足j| 国产单亲对白刺激| 大香蕉久久网| 国产日韩欧美视频二区| e午夜精品久久久久久久| 国产欧美日韩精品亚洲av| av免费在线观看网站| 人人妻人人澡人人看| 国产精品国产av在线观看| 久久精品亚洲精品国产色婷小说| tocl精华| 国产淫语在线视频| 在线观看舔阴道视频| 亚洲一卡2卡3卡4卡5卡精品中文| 日本欧美视频一区| 亚洲av片天天在线观看|