• <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
    亚洲三级黄色毛片| 国产毛片a区久久久久| 麻豆精品久久久久久蜜桃| 狂野欧美激情性xxxx在线观看| 美女大奶头视频| 亚洲精品乱码久久久久久按摩| 精品一区二区三区人妻视频| 黄色一级大片看看| 久久精品国产亚洲av涩爱| 免费观看精品视频网站| 免费黄网站久久成人精品| 亚洲国产成人一精品久久久| 欧美丝袜亚洲另类| 精品欧美国产一区二区三| 乱人视频在线观看| 色吧在线观看| 久久久久精品久久久久真实原创| 久久99精品国语久久久| 国产精品综合久久久久久久免费| 精品国产一区二区三区久久久樱花 | www.色视频.com| 国产乱人视频| 神马国产精品三级电影在线观看| 神马国产精品三级电影在线观看| 国产爱豆传媒在线观看| 免费观看在线日韩| 国产男女超爽视频在线观看| av在线亚洲专区| 高清视频免费观看一区二区 | 在线免费观看的www视频| 色吧在线观看| 七月丁香在线播放| 免费电影在线观看免费观看| 国产伦精品一区二区三区视频9| 久久久久精品久久久久真实原创| 街头女战士在线观看网站| 国产女主播在线喷水免费视频网站 | 国内精品美女久久久久久| 全区人妻精品视频| 国产熟女欧美一区二区| 麻豆久久精品国产亚洲av| 自拍偷自拍亚洲精品老妇| 免费观看a级毛片全部| 日日撸夜夜添| 美女xxoo啪啪120秒动态图| 美女xxoo啪啪120秒动态图| 一级毛片久久久久久久久女| 一级二级三级毛片免费看| 男女边吃奶边做爰视频| 又粗又硬又长又爽又黄的视频| 免费不卡的大黄色大毛片视频在线观看 | 综合色av麻豆| 一级毛片久久久久久久久女| 观看美女的网站| 久久韩国三级中文字幕| 插阴视频在线观看视频| 国产老妇伦熟女老妇高清| 一级黄片播放器| 特大巨黑吊av在线直播| 又爽又黄a免费视频| 亚洲在线观看片| 最近最新中文字幕大全电影3| 99热这里只有是精品50| 国产在线男女| 麻豆成人av视频| 在线a可以看的网站| 人人妻人人澡人人爽人人夜夜 | 亚洲欧美成人综合另类久久久| 国产精品久久久久久精品电影| 国产一区亚洲一区在线观看| 男女边吃奶边做爰视频| 日韩成人av中文字幕在线观看| 亚洲高清免费不卡视频| 在线免费十八禁| 成年免费大片在线观看| 精品人妻一区二区三区麻豆| 夜夜爽夜夜爽视频| 亚洲人成网站在线观看播放| 中国美白少妇内射xxxbb| 伊人久久精品亚洲午夜| 久久久久久久久久久免费av| 天堂中文最新版在线下载 | 亚洲性久久影院| 成人高潮视频无遮挡免费网站| 国产成人精品福利久久| 人体艺术视频欧美日本| 亚洲精品自拍成人| 日韩精品有码人妻一区| 国产成人精品婷婷| 国产精品爽爽va在线观看网站| 两个人视频免费观看高清| 午夜免费激情av| 99热6这里只有精品| 亚洲不卡免费看| 中文字幕久久专区| 亚洲无线观看免费| 看十八女毛片水多多多| 91精品一卡2卡3卡4卡| 天堂av国产一区二区熟女人妻| 国产成人一区二区在线| 99视频精品全部免费 在线| 午夜福利网站1000一区二区三区| 只有这里有精品99| 日本wwww免费看| 日韩三级伦理在线观看| 丰满乱子伦码专区| 深爱激情五月婷婷| 毛片一级片免费看久久久久| 一本久久精品| 麻豆av噜噜一区二区三区| 精品久久久久久久久久久久久| 免费观看精品视频网站| 亚洲精华国产精华液的使用体验| av天堂中文字幕网| 1000部很黄的大片| 只有这里有精品99| 欧美日韩综合久久久久久| 内地一区二区视频在线| 亚洲国产精品成人综合色| 777米奇影视久久| 国产一区二区三区av在线| 身体一侧抽搐| 成年版毛片免费区| av在线亚洲专区| 精品久久久噜噜| 国产亚洲av片在线观看秒播厂 | 久久精品久久久久久久性| 人体艺术视频欧美日本| 国产亚洲91精品色在线| 国产成人福利小说| 18禁在线无遮挡免费观看视频| 亚洲国产精品成人综合色| 最近视频中文字幕2019在线8| 日韩av在线大香蕉| 精品欧美国产一区二区三| 国产v大片淫在线免费观看| 天堂网av新在线| 日韩,欧美,国产一区二区三区| 两个人的视频大全免费| 精华霜和精华液先用哪个| 男女国产视频网站| 久久97久久精品| 久久久久九九精品影院| 国产69精品久久久久777片| 丝瓜视频免费看黄片| 国产淫语在线视频| 黄色一级大片看看| 成人毛片a级毛片在线播放| 精华霜和精华液先用哪个| 2018国产大陆天天弄谢| av卡一久久| 人妻少妇偷人精品九色| 麻豆av噜噜一区二区三区| 亚洲av在线观看美女高潮| 国产激情偷乱视频一区二区| 免费无遮挡裸体视频| 美女脱内裤让男人舔精品视频| 男女边摸边吃奶| 日韩制服骚丝袜av| 偷拍熟女少妇极品色| 久久久久九九精品影院| 久久久久性生活片| 国产免费又黄又爽又色| 国产精品一区二区三区四区免费观看| 成人毛片a级毛片在线播放| 26uuu在线亚洲综合色| 国产视频首页在线观看| 久久午夜福利片| 免费在线观看成人毛片| 国产乱人偷精品视频| 一个人看视频在线观看www免费| 国产女主播在线喷水免费视频网站 | 久久久久久国产a免费观看| 99久国产av精品国产电影| 伊人久久国产一区二区| 看免费成人av毛片| 亚洲精品乱码久久久久久按摩| 又粗又硬又长又爽又黄的视频| 天天躁日日操中文字幕| 最近2019中文字幕mv第一页| 高清在线视频一区二区三区| 国产精品.久久久| 国产黄色免费在线视频| 亚洲av免费在线观看| 国产黄a三级三级三级人| 国产精品一及| 男女下面进入的视频免费午夜| 欧美日韩在线观看h| ponron亚洲| 寂寞人妻少妇视频99o| av福利片在线观看| 国内精品一区二区在线观看| 在线播放无遮挡| 99热这里只有是精品50| 亚洲国产精品国产精品| 韩国高清视频一区二区三区| 国产精品久久视频播放| 欧美区成人在线视频| 久久久久久久亚洲中文字幕| 51国产日韩欧美| 久99久视频精品免费| 18禁动态无遮挡网站| 日韩欧美一区视频在线观看 | 丝袜美腿在线中文| 春色校园在线视频观看| 免费在线观看成人毛片| av线在线观看网站| 啦啦啦韩国在线观看视频| 最近中文字幕2019免费版| 亚洲精品自拍成人| 日韩三级伦理在线观看| 大片免费播放器 马上看| 99热全是精品| 国产在线一区二区三区精| 小蜜桃在线观看免费完整版高清| 国产免费福利视频在线观看| 日产精品乱码卡一卡2卡三| 日韩欧美三级三区| 搡老乐熟女国产| 狂野欧美白嫩少妇大欣赏| 国产在视频线在精品| 久久久久久久久久人人人人人人| 日韩一本色道免费dvd| 亚洲自拍偷在线| 久久精品综合一区二区三区| 欧美xxxx黑人xx丫x性爽| 亚洲精品第二区| 国产 一区 欧美 日韩| 中国美白少妇内射xxxbb| 大陆偷拍与自拍| 免费高清在线观看视频在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人美女网站在线观看视频| 免费av不卡在线播放| 夜夜看夜夜爽夜夜摸| 男的添女的下面高潮视频| 少妇丰满av| 噜噜噜噜噜久久久久久91| 在线免费十八禁| 晚上一个人看的免费电影| 久久久久久久亚洲中文字幕| 亚洲在久久综合| 国产精品久久久久久久电影| 男人舔女人下体高潮全视频| 国精品久久久久久国模美| 性色avwww在线观看| 亚洲精品乱码久久久v下载方式| 国产午夜精品一二区理论片| 99视频精品全部免费 在线| 日韩中字成人| 美女xxoo啪啪120秒动态图| 午夜福利视频1000在线观看| 亚洲一区高清亚洲精品| 欧美 日韩 精品 国产| 高清视频免费观看一区二区 | 有码 亚洲区| 国产老妇伦熟女老妇高清| 亚洲av中文字字幕乱码综合| 久久久久久久久久久免费av| 在线观看一区二区三区| 中国国产av一级| 人妻少妇偷人精品九色| 午夜免费观看性视频| 女人久久www免费人成看片| 汤姆久久久久久久影院中文字幕 | 免费大片18禁| 性色avwww在线观看| 亚洲国产成人一精品久久久| 看黄色毛片网站| 亚洲av中文av极速乱| 久久久久久久久大av| 日韩成人av中文字幕在线观看| 在线免费十八禁| 久久久久久久久中文| 九九在线视频观看精品| 男人舔奶头视频| 视频中文字幕在线观看| 最近2019中文字幕mv第一页| av专区在线播放| 中文天堂在线官网| 美女国产视频在线观看| 少妇裸体淫交视频免费看高清| 国模一区二区三区四区视频| 尤物成人国产欧美一区二区三区| 在线播放无遮挡| 日韩av不卡免费在线播放| 91久久精品国产一区二区三区| 99久国产av精品国产电影| 成人二区视频| 尾随美女入室| av播播在线观看一区| 国产av不卡久久| 日日摸夜夜添夜夜爱| 欧美成人午夜免费资源| 高清视频免费观看一区二区 | 成人午夜高清在线视频| 久久久色成人| 蜜桃久久精品国产亚洲av| 免费看不卡的av| 亚洲真实伦在线观看| 国产黄a三级三级三级人| 在现免费观看毛片| 赤兔流量卡办理| 麻豆久久精品国产亚洲av| 亚洲熟妇中文字幕五十中出| 丝袜美腿在线中文| 国产69精品久久久久777片| 国产成人a区在线观看| 网址你懂的国产日韩在线| 成人国产麻豆网| 两个人视频免费观看高清| 99久国产av精品| 在线a可以看的网站| 国产精品久久久久久久电影| 久久这里只有精品中国| 欧美日韩国产mv在线观看视频 | av福利片在线观看| 自拍偷自拍亚洲精品老妇| 久久久精品欧美日韩精品| 日本欧美国产在线视频| 国产日韩欧美在线精品| 国产高潮美女av| 久久精品综合一区二区三区| 国产人妻一区二区三区在| 高清视频免费观看一区二区 | 日本免费a在线| 亚洲精品国产av成人精品| 日韩视频在线欧美| 高清午夜精品一区二区三区| 精品久久久久久电影网| 丰满少妇做爰视频| 天堂网av新在线| 亚洲电影在线观看av| 男人舔女人下体高潮全视频| 美女cb高潮喷水在线观看| av黄色大香蕉| 夫妻性生交免费视频一级片| 97人妻精品一区二区三区麻豆| 久久6这里有精品| 国产大屁股一区二区在线视频| 国产免费一级a男人的天堂| 亚洲18禁久久av| 成年人午夜在线观看视频 | 久久精品熟女亚洲av麻豆精品 | 成人毛片60女人毛片免费| 97超视频在线观看视频| 成人性生交大片免费视频hd| 久久这里有精品视频免费| 国产乱人视频| 日本猛色少妇xxxxx猛交久久| 亚洲最大成人手机在线| 国产v大片淫在线免费观看| 国产精品不卡视频一区二区| 男女边摸边吃奶| 视频中文字幕在线观看| 亚洲精品中文字幕在线视频 | 久久精品熟女亚洲av麻豆精品 | 亚洲精品视频女| 亚洲欧洲日产国产| 久久久久精品性色| 人妻一区二区av| 亚洲av国产av综合av卡| 嫩草影院精品99| 99久久人妻综合| 亚洲成人av在线免费| 2018国产大陆天天弄谢| 极品教师在线视频| 熟妇人妻久久中文字幕3abv| 亚洲美女视频黄频| 大香蕉97超碰在线| 亚洲在线观看片| av天堂中文字幕网| 国产成人免费观看mmmm| 男人爽女人下面视频在线观看| 人妻一区二区av| 国产永久视频网站| av一本久久久久| a级毛色黄片| 免费看光身美女| 婷婷色综合大香蕉| 久久国内精品自在自线图片| 夜夜爽夜夜爽视频| 国产爱豆传媒在线观看| 午夜亚洲福利在线播放| 欧美极品一区二区三区四区| 夫妻性生交免费视频一级片| 麻豆国产97在线/欧美| 青春草国产在线视频| 国产有黄有色有爽视频| av一本久久久久| 干丝袜人妻中文字幕| 国产成人freesex在线| av又黄又爽大尺度在线免费看| 亚洲成人久久爱视频| 国产精品一区二区三区四区免费观看| 最近最新中文字幕免费大全7| 亚洲av免费高清在线观看| 亚洲av中文字字幕乱码综合| 男人舔女人下体高潮全视频| 欧美xxxx黑人xx丫x性爽| 可以在线观看毛片的网站| av免费观看日本| 午夜福利在线观看吧| 日韩欧美国产在线观看| 日韩精品青青久久久久久| 18禁在线播放成人免费| 国产在线一区二区三区精| 久久久午夜欧美精品| 久久99热这里只频精品6学生| 国产在视频线在精品| 哪个播放器可以免费观看大片| 少妇猛男粗大的猛烈进出视频 | 性色avwww在线观看| 亚洲av中文字字幕乱码综合| 91aial.com中文字幕在线观看| 久久97久久精品| 免费看不卡的av| 久久精品国产自在天天线| 性色avwww在线观看| 在线播放无遮挡| 婷婷六月久久综合丁香| 秋霞在线观看毛片| 国内精品一区二区在线观看| 日韩av免费高清视频| 久久精品久久久久久久性| 黄色一级大片看看| 麻豆国产97在线/欧美| 啦啦啦啦在线视频资源| av天堂中文字幕网| 偷拍熟女少妇极品色| 日本三级黄在线观看| 国产成人精品福利久久| 亚洲美女视频黄频| 亚洲精品影视一区二区三区av| 久久久欧美国产精品| 国产成人91sexporn| 在线免费观看不下载黄p国产| 久久鲁丝午夜福利片| 免费不卡的大黄色大毛片视频在线观看 | 午夜福利视频1000在线观看| 亚洲av二区三区四区| 国产在线一区二区三区精| 国产91av在线免费观看| 久久精品人妻少妇| 丰满少妇做爰视频| 日韩电影二区| 亚洲精品,欧美精品| 乱系列少妇在线播放| 国产黄片美女视频| 18禁在线播放成人免费| 国产激情偷乱视频一区二区| 久久久久久久久久黄片| 九九爱精品视频在线观看| 深夜a级毛片| 三级毛片av免费| 国产成人精品久久久久久| 欧美xxxx黑人xx丫x性爽| 亚洲性久久影院| 夫妻午夜视频| 午夜精品一区二区三区免费看| 亚洲av成人精品一二三区| 国产一级毛片在线| 国产黄色免费在线视频| 一本一本综合久久| 亚洲色图av天堂| 欧美日韩视频高清一区二区三区二| 一级毛片aaaaaa免费看小| 亚洲精品乱久久久久久| 久久精品国产鲁丝片午夜精品| 欧美日韩精品成人综合77777| 成年人午夜在线观看视频 | 日韩精品青青久久久久久| 亚洲内射少妇av| 亚洲国产成人一精品久久久| 人妻夜夜爽99麻豆av| 亚洲欧美清纯卡通| av在线亚洲专区| 色吧在线观看| 别揉我奶头 嗯啊视频| 国产爱豆传媒在线观看| 亚洲四区av| 九草在线视频观看| 九九久久精品国产亚洲av麻豆| 超碰97精品在线观看| 欧美xxⅹ黑人| a级毛片免费高清观看在线播放| 日本一本二区三区精品| 亚洲精品乱码久久久久久按摩| 嫩草影院入口| 十八禁网站网址无遮挡 | 亚洲怡红院男人天堂| 免费黄网站久久成人精品| 一级毛片黄色毛片免费观看视频| 国产综合精华液| 91狼人影院| 亚洲精品色激情综合| 日韩欧美国产在线观看| 婷婷色麻豆天堂久久| 亚洲va在线va天堂va国产| 国产高清不卡午夜福利| 美女xxoo啪啪120秒动态图| 日韩欧美一区视频在线观看 | 亚洲久久久久久中文字幕| 国产一区二区三区av在线| 最后的刺客免费高清国语| 国产精品.久久久| 久久久久久国产a免费观看| 久久久久久久久久久丰满| 26uuu在线亚洲综合色| 日本黄色片子视频| 亚洲av中文av极速乱| 国产成人a∨麻豆精品| av专区在线播放| 亚洲国产欧美在线一区| 免费黄网站久久成人精品| 少妇被粗大猛烈的视频| 最近2019中文字幕mv第一页| 大陆偷拍与自拍| 成人无遮挡网站| 色5月婷婷丁香| 亚洲av电影在线观看一区二区三区 | 一级毛片久久久久久久久女| 国产色爽女视频免费观看| 精品人妻偷拍中文字幕| 亚洲国产精品成人久久小说| 亚洲成色77777| 大片免费播放器 马上看| 亚洲无线观看免费| 精品午夜福利在线看| 精品久久久噜噜| 黄片wwwwww| 欧美成人精品欧美一级黄| 大话2 男鬼变身卡| 日韩成人av中文字幕在线观看| 国产精品不卡视频一区二区| 欧美精品国产亚洲| 久久久久久久亚洲中文字幕| 婷婷色综合大香蕉| 国产人妻一区二区三区在| 1000部很黄的大片| 一区二区三区免费毛片| 国产精品一区二区在线观看99 | 嫩草影院新地址| 青春草视频在线免费观看| 国产国拍精品亚洲av在线观看| 99热这里只有精品一区| 亚洲欧美成人综合另类久久久| 亚洲一区高清亚洲精品| 国产午夜精品久久久久久一区二区三区| 精品一区二区三区视频在线| h日本视频在线播放| 日韩成人av中文字幕在线观看| 国产激情偷乱视频一区二区| 亚洲精品色激情综合| 亚洲人成网站在线播| 亚洲欧美精品专区久久| 啦啦啦啦在线视频资源| 日韩欧美 国产精品| 午夜激情欧美在线| 三级国产精品欧美在线观看| 嫩草影院入口| 成人亚洲欧美一区二区av| 欧美激情在线99| av在线亚洲专区| 搡老妇女老女人老熟妇| 成年女人在线观看亚洲视频 | av在线蜜桃| 国产精品一区www在线观看| 亚洲av二区三区四区| 国产精品精品国产色婷婷| 午夜视频国产福利| 2018国产大陆天天弄谢| 日韩欧美一区视频在线观看 | 国产伦精品一区二区三区视频9| 国产中年淑女户外野战色| 大话2 男鬼变身卡| 国产精品一区www在线观看| 永久免费av网站大全| 免费黄色在线免费观看| 在线a可以看的网站| 亚洲av免费高清在线观看| 色综合色国产| 最后的刺客免费高清国语| 网址你懂的国产日韩在线| av在线天堂中文字幕| 国产伦在线观看视频一区| 舔av片在线| 日韩,欧美,国产一区二区三区| 日本免费a在线| 日韩欧美一区视频在线观看 | 国产男人的电影天堂91| 国产老妇女一区| 日韩三级伦理在线观看| 成人亚洲精品av一区二区| 2022亚洲国产成人精品| 国产69精品久久久久777片| 免费高清在线观看视频在线观看| 婷婷色av中文字幕| 午夜日本视频在线| 又爽又黄无遮挡网站| 只有这里有精品99| 中文乱码字字幕精品一区二区三区 | 身体一侧抽搐| 丝袜美腿在线中文| 看免费成人av毛片| 一区二区三区免费毛片| 欧美区成人在线视频| 国产精品美女特级片免费视频播放器| 成人毛片a级毛片在线播放| 美女高潮的动态| 自拍偷自拍亚洲精品老妇| 欧美 日韩 精品 国产| 成人国产麻豆网| 最近最新中文字幕大全电影3| 我的女老师完整版在线观看|