AN ARITHMETIC METHOD OF SEGMENTING MOVING AEOLIAN SAND PARTICLES＇ IMAGES FROM BACKGROUND INFORMATION OF DIGITAL HIGH-SPEED PHOTOGRAPHY IMAGES

为了深入理解多因素驱动下风沙颗粒起动的动态演化规律,需要准确地获得沙质床面附近沙粒群起跳的方式、速度和运动轨迹.以连续强激光源为照明的数字高速摄影技术是研究这类问题的有效手段,但由于风沙运动的高速摄影图像具有运动沙粒和静止床面的对比度小、相邻两帧图像相似性小等特点,原始图像叠加算法难以有效实现目标与背景的分割.该文提出了基于相邻的风沙运动图像灰度差值变化原理的图像分割算法.实例显示,只要选择合适的相邻图像灰度差值阈值和自适应二值化处理方法就能实现图像中运动沙粒与床面分割.当起跳沙粒浓度较低情况下,基于MATLAB平台的最小距离匹配的粒子追踪算法（PTV算法）能较为准确地恢复床面附近沙粒的运动轨迹.     

一种改进的高浓度风沙图像的动态灰度阈值分割算法

风沙跃移是塑造干旱区地貌的最主要的动力,其伴生的粉尘的释放、输送和沉降不但严重地影响了大气环境质量,也引起了全球气候 系统和海洋系统的变化. 单个沙粒轨迹形成是理解粒--床碰撞过程和风沙两相流耦合过程的重要纽带. 高浓度风沙图像处理算法是深入理解沙粒轨迹形成机制的关键技术. 为了实现对高浓度风沙图像精细处理, 本文提出了不依赖于经验参数的动态灰度阈值分割算法,它包括背景模板去噪、绿光通道灰度化处理、图像微分、灰度方差阈值目标检测和 最大类间方差灰度阈值分割等,其中背景模板去噪和灰度方差阈值目标检测等是新算法的主要亮点. 高浓度风沙图像分割实验显示, 扣减背景模板去除了条纹状和斑点状的稳定噪声;图像微分和灰度方差阈值目标检测显著地提高了暗沙粒识别的数量并有效地去除了随机噪声; 改进算法的沙粒有效识别个数(人 机判读坐标一致的沙粒数目)、查全率(计算机提取沙粒数目与人工判读的实际沙粒数比值)和查准率(有效沙粒数目与实际沙粒数目的比值) 分别为461, 71%和86%,显著地高于传统算法对应的85, 13%和82%,这表明新算法对高浓度风沙图像的分割效果良好. 避免单粒子分割和表观重叠问题是进一步完善高浓度风沙图像的分割算法的可能途径.      

一种改进的高浓度风沙图像的动态灰度阈值分割算法

风沙跃移是塑造干旱区地貌的最主要的动力,其伴生的粉尘的释放、输送和沉降不但严重地影响了大气环境质量,也引起了全球气候系统和海洋系统的变化.单个沙粒轨迹形成是理解粒-床碰撞过程和风沙两相流耦合过程的重要纽带.高浓度风沙图像处理算法是深入理解沙粒轨迹形成机制的关键技术.为了实现对高浓度风沙图像精细处理,本文提出了不依赖于经验参数的动态灰度阈值分割算法,它包括背景模板去噪、绿光通道灰度化处理、图像微分、灰度方差阈值目标检测和最大类间方差灰度阈值分割等,其中背景模板去噪和灰度方差阈值目标检测等是新算法的主要亮点.高浓度风沙图像分割实验显示,扣减背景模板去除了条纹状和斑点状的稳定噪声;图像微分和灰度方差阈值目标检测显著地提高了暗沙粒识别的数量并有效地去除了随机噪声;改进算法的沙粒有效识别个数(人机判读坐标一致的沙粒数目)、查全率(计算机提取沙粒数目与人工判读的实际沙粒数比值)和查准率(有效沙粒数目与实际沙粒数目的比值)分别为461,71%和86%,显著地高于传统算法对应的85,13%和82%,这表明新算法对高浓度风沙图像的分割效果良好.避免单粒子分割和表观重叠问题是进一步完善高浓度风沙图像的分割算法的可能途径.     

非标识PIV图像的模糊聚类识别方法

文中将模糊分析的原理应用于分析粒子在一幅ＰＩＶ图像中的运动情况。聚类由几个相邻的粒子开始,聚类的原则是流体运动的基本规律,随着流体的运动,粒子簇在第二次曝光时改变其形状及位置,通过模糊聚类可获得粒子运动的规律。应用聚类分析方法可改善非标识ＰＩＶ图像的处理效率,基于这种方法的粒子运动轨迹追踪算法被提出并应用于ＰＩＶ图像分析。     

沙粒跃移运动的数值模拟

沙粒在风场作用下形成的风沙运动,直接导致一系列环境问题.现有的数值模拟并没有揭示风沙流的二维性质.本文提出一种合理有效的模拟方法来揭示风沙流时空的变化.气相运动采用大涡模型,而固相(沙粒)处理为高散体系(DEM),并利用已有的沙粒与床面碰撞的激溅函数,以获得整个风沙系统运动特性.结果表明:其发展过程与以往的试验结果相比较符合:风沙流达到稳定后,不同截面输沙率沿高度分布与平均输沙率相比差别很大,并大约在9cm高处,平均输沙率出现极值,单宽输沙率沿风向分布同样具有非均匀性.     

风沙两相流PIV测量算法研究

在风沙两相流图像特征的基础上，提出了一种基于模式识别动态聚类方法中$K$-均值 算法的数字面具(Digital Mask)自动生成算法来求解风沙两相流动. 简化了传统生成Digital Mask过程中手动设置参数的操作，减少了人为误差，为批量处理风沙两相流PIV图像提供 了一种安全快捷的方法. 并将该算法应用于风沙两相流的PIV实际测量，分别得到了不同流 动状态下的气流、沙粒以及风沙气固两相流的速度场.     

NUMERICAL SIMULATION OF WIND-BLOWN SAND MOVEMENT BASED ON SPH

运用光滑粒子流体动力学(smoothed particles hydrodynamics, SPH)方法对沙粒和气流的相互耦合运动特性进行了分析,研究提出了风沙流的SPH数值方法并进行了数值模拟. 首先提出了风沙流的SPH建模方法和基本理论,建立了风沙流动的SPH数值模拟平台. 其次通过建立风沙流的SPH 模型并施加边界条件,对自然风作用下沙粒的运动情况进行了数值模拟,详细分析了沙粒运动轨迹及特性,最后通过与相关研究成果对比分析,验证了完善后的SPH方法有效性. 通过考虑气流场的可变性,在风沙流SPH计算模型中引入了加载(起风)和卸载(停风)方式,观察并对比分析了沙粒的运动轨迹和特性. 为进一步研究风沙流的实时动态非线性行为提供了SPH理论基础和数值分析方法.     

风沙运动的沙粒带电机理及其影响的研究进展

土壤风蚀以及沙尘暴是人类面临的重大自然环境灾害之一,它直接关系着人类的生存环境和生活质量.许多国家和政府以及科学家们对此一直十分关注并开展了大量研究.由于大气和地表条件的多样性以及地表沙粒的运动与风场作用的相互耦合,导致风沙运动问题研究的复杂性.除此之外,伴随着对风沙运动机理的研究和对沙尘暴的观测,人们发现在风沙跃移运动和沙尘暴中的沙粒带有电荷并形成风沙电场.沙粒的运动不仅导致沙粒带电和风沙电场形成,而且也直接受到沙粒带电和风沙电场的影响.为了弄清风沙电的起因以及揭示其影响规律,研究人员对此进行了野外观测、实验室测量和理论推测等.本文将就这一问题重点介绍有关风沙运动中沙粒带电的机理、风沙电场分布规律、沙粒带电对风沙运动的微宏观物理量以及沙尘暴中沙粒带电对电磁波散射的影响等方面的实验和理论研究的基本概况和进展以及作者等人在此方面的研究工作.主要包含内容有:沙粒带电及其电场的基本实验测量,粒带电对风沙运动的影响,沙尘暴对电磁波衰减的影响研究以及有关风沙起电的机理研究等. 同时,对这一研究领域中目前关注与存在的主要问题给予评述.图6参81     

基于可见光与红外图像特征融合的目标跟踪

针对单一图像源下目标跟踪精度不高的问题，利用跟踪状态下的目标存在于可见光与红外图像中的特征对连续自适应均值移动跟踪算法做出改进。首先选取可见光图像的“颜色梯度背投影”作为改进的目标模型，选取红外图像的“灰度梯度背投影”作为改进的目标模型；然后根据可见光序列图像和红外序列图像各自进行连续自适应均值移动跟踪算法得到的对应的qi系数判定两种图像跟踪的效果，对两种图像的权重进行自适应调整，得到这两种图像的特征级融合图像和跟踪结果。实验结果表明，对于320像素×240像素的可见光和红外图像，基于可见光与红外图像特征融合的目标跟踪算法在复杂背景下能够较准确的跟踪目标，目标跟踪精度为0.5像素，跟踪速度为30~32 ms/帧。     

基于U-net神经网络算法和改进的细化算法的水坝混凝土裂缝测量

以水坝裂缝的视觉识别与计算为目标,融合U-net神经网络算法和改进图像细化算法,提出大坝裂缝识别和宽度计算的方法.首先,采用U-net处理1500张水坝裂缝图像样本,得到裂缝预分割模型;接着,采用改进的细化算法提取裂缝骨架分割模型;最后,结合预分割结果和骨架分割模型得到水坝裂缝的准确信息.结果表明,使用U-net神经网...     

A THREE-VALUE-SEGMENTATION METHOD OF COAL CONTAINING INCLUSION BASED ON GRAY DISTRIBUTION OF COMPUTED TOMOGRAPHY IMAGE1)

煤体CT图像阈值分割是其三维真实结构重构的前提. 基于单轴压缩状态下煤体的CT图像,提出了一种基于图像灰度分布的孔隙度计算方法, 结合灰度直方图,利用实测孔隙度反推法对煤体CT图像进行了裂隙--煤基质--煤杂质三值化,同时与 最大类间方差法和最大熵法的三值化结果进行了对比分析. 结果表明, 最大熵法对于煤基质和煤杂质不能较好区分; 最大类间方差法对于图像微小区域灰度的突变(裂隙的产生)敏感性低,无法有效区分裂隙与煤基质; 本文改进的灰度直方图法避免了最大类间方差法和最大熵法的劣势,具有较好的图像三值化效果,且煤基质与煤杂质间的阈值基 于实验结果得到,具有更高的可信度和准确度.     

两相流PIV粒子图像处理方法的研究

本文在单相PIV技术的基础上研究了两相流动PIV图像处理方法，采用摸板匹配法和灰度加权标定法对两相粒子进行了识别、区分和标定，采用灰度互相关法对区分后的单相粒子图像进行了处理，应用基于以上方法编制的Windows应用软件，首先对由美国Minnesota大学复杂流动实验室提供的两相流动粒子图片进行了处理，通过对比分析可见，应用本文所采用的方法能对两相粒子进行有效的识别和区分，然后以搅拌槽内液固两相流场为例对此方法进行了应用。     

Sand particle dislodgement in windblown sand

The incipient motion of sand particle from sand bed plays a very important role in the prediction of windblown sand.In this paper,we proposed a new method for predicting the incipient motion of sand particle based on wind speed fluctuation as follows,when the wind speed is larger than the critical wind speed,if the total impulse on sand particle is larger than the critical impulse,incipient motion of sand particle would take place,otherwise if not.Furthermore,from the analysis of entrainment in the rolling and lifting modes,we come to the following conclusion.When the average wind speed is smaller than the critical wind speed,if the average wind speed is used to judge the incipient motion of sand particle,one will underestimate the number of sand particles jumping from the bed,if the instantaneous wind speed is used to judge incipient motion of sand particle,one will overestimate the number of sand particles jumping from the bed;When the average wind speed is larger than the critical wind speed,either the average or the instantaneous wind speeds is used to judge the incipientmotion of sand particles,one will overestimate the number of sand particles jumping from the bed.     

Rotation speed measurement of moving particles in a CFB riser

This paper presents a study on measuring rotation speed of moving glass beads with an average diameter of 500 p.m in a pilot-scale circulating fluidized bed （CFB） riser with a high-speed digital imaging system. Two methods have been developed to calculate particle rotation speed from the particle images. The first method consists of a fully automated algorithm based on cross-correlation of gray distribution of particle images for particles whose rotation axes are （nearly） perpendicular to the imaging plane, and the second method calculates the speed of particle rotation by identifying its rotation axis using two or more characteristic points on its surface. The reliability of the two methods is verified by using a small sphere with known speed and direction of rotation. The first method is shown to be capable of measuring accurately the rotation speed for the particle with a rotation axis （nearly） perpendicular to the imaging plane and filtering off other particles using an appropriate threshold of correlation coefficient. The second method is shown to be capable of yielding both the speed and direction of particle rotation, with a measurement error of less than 10%. Results of both methods on real glass beads in a CFB riser are compared against each other.     

Determination of the particle velocity,size, and distribution function along the diameters in a two-phase jet

A method for simultaneous determination of the particle size and velocity in a supersonic, two-phase luminescent jet consisting of a gas phase and the particles is considered. The jet is subjected to external transverse blowing by a hypersonic gas stream in two modes under whose effect it separates into the gas phase and a particle stream. The particle trajectories were determined by photography. A comparison between the theoretically computed and the experimental trajectories permits determination of the particle diameter and velocity. The particle distribution function over the diameters was found by means of the light intensity emitted by the stream of particles.     

Experimental and numerical investigation of the shock-induced fluidization of a particles bed

An experimental study on unsteady two phase flow is conducted in a vertical shock tube. Shock Mach numbers range from 1.3 to 1.5 in 1 bar. The particles are initially positioned in horizontal beds of various thicknesses. Our research covers a large domain of void fraction from 1 (single particles) to 0.35 (compact beds). The experiments provide shadowgraph images for the recording of particle trajectories (effect of the gas on the particles) and side-wall pressures (action of the particles on the gas). A dense two phase flow model has been elaborated and numerically solved using a finite difference scheme with pseudoviscosity. The simulated shock-induced fluidization of a 2 cm thick bed of 1.5 mm diameter glass particles is compared to the experiment. Received 10 September 1996 / Accepted 4 January 1997     

风沙运动力学机理研究的历史、进展与趋势

回顾了风沙运动力学机理研究历史,阐述了中国在风沙运动力学机理研究方 面所取得的最新进展，提出了目前存在的问题和今后发展的趋势.     

Background extraction from double-frame PIV images

This study presents a simple image pre-processing scheme to extract background information from double-frame particle image velocimetry (PIV) images. Everything that stays stagnant in the image (e.g., image background and light reflections from stationary objects) is assumed to be a source of disturbance and is removed by subtracting the second frame of the image pair from the first frame. This yields a single frame difference image, which is transferred back to a background extracted double-frame image. After the procedure the background in the image does not correlate with itself anymore and therefore a bias error in PIV analysis towards zero displacement is avoided. The simulations show that the procedure conserves the profile of tracer particle images when a displacement between the frames is larger than a particle image size. The performance of this procedure is emphasized with various examples, and extensions of the procedure are introduced. The extended procedure extracts background objects that move between the image frames, e.g. dispersed phase particles in a two-phase flow or laser light sheet reflections from moving objects.     

基于颗粒离散元的沙纹演变大涡模拟研究

应用大涡模拟、"点球" 浸入边界法和基于"事件驱动" 模型和"弹簧-阻尼" 模型的颗粒离散元法, 数值模拟了明渠湍流中沙纹的演变. 通过对不同谢尔兹数下无沙纹床面的推移质输沙率进行计算, 并与经典输沙率公式进行对比, 验证了模型的精度和可靠性. 随后, 对明渠湍流中沙纹床面的演变过程进行了数值模拟, 计算了有沙纹床面的推移质输沙率、沙纹长度和高度、等效床面高度的最大值、最小值和平均值、沙纹形状阻力、体积流速随时间的变化曲线. 研究发现:初始平整的床面在较短的时间内(tU_b=h≈100) 发展出数条沙纹, 随后沙纹逐渐发展, 在tU_b=h为1 600~2 000 时, 沙纹发生合并. 在沙纹数量不变的条件下, 沙纹高度随时间近似呈线性增长, 而沙纹的长度的平均值却保持恒定. 随着沙纹高度的增大, 输沙率和体积流速逐渐降低, 沙纹形状阻力则逐渐增大;当沙纹发生合并时, 沙纹高度快速增加, 输沙率、体积流速和沙纹形状阻力也出现了大幅跳跃. 在同等的水流强度条件下, 有沙纹床面的输沙率小于平整床面的输沙率.     

In-line imaging measurements of particle size,velocity and concentration in a particulate two-phase flow

A novel method is developed for in-line measurements of particle size, velocity and concentration in a dilute, particulate two-phase flow based on trajectory image processing. The measurement system consists of a common industrial CCD camera, an inexpensive LED light and a telecentric lens. In this work, the image pre-processing steps include stitching, illumination correction, binarization, denoising, and the elimination of unreal and defocused particles. A top-hat transformation is found to be very effective for the binarization of images with non-uniform background illumination. Particle trajectories measured within a certain exposure time are used to directly obtain particle size and velocity. The particle concentration is calculated by using the statistics of recognized particles within the field of view. We validate our method by analyzing experiments in a gas-droplet cyclone separator. This in-line image processing method can significantly reduce the measurement cost and avoid the data inversion process involved in the light scattering method.