CN107220946A - A kind of real-time eliminating method of bad lumpiness image on rock transportation band - Google Patents
A kind of real-time eliminating method of bad lumpiness image on rock transportation band Download PDFInfo
- Publication number
- CN107220946A CN107220946A CN201710353593.3A CN201710353593A CN107220946A CN 107220946 A CN107220946 A CN 107220946A CN 201710353593 A CN201710353593 A CN 201710353593A CN 107220946 A CN107220946 A CN 107220946A
- Authority
- CN
- China
- Prior art keywords
- image
- mrow
- mtd
- gradient
- msub
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 31
- ZPUCINDJVBIVPJ-LJISPDSOSA-N cocaine Chemical compound O([C@H]1C[C@@H]2CC[C@@H](N2C)[C@H]1C(=O)OC)C(=O)C1=CC=CC=C1 ZPUCINDJVBIVPJ-LJISPDSOSA-N 0.000 claims abstract description 30
- 238000009499 grossing Methods 0.000 claims abstract description 14
- 238000001914 filtration Methods 0.000 claims abstract description 11
- 230000018199 S phase Effects 0.000 claims description 5
- 230000004069 differentiation Effects 0.000 claims description 2
- 230000008030 elimination Effects 0.000 abstract description 4
- 238000003379 elimination reaction Methods 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 description 9
- 239000004575 stone Substances 0.000 description 6
- 239000002245 particle Substances 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000005065 mining Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000013467 fragmentation Methods 0.000 description 2
- 238000006062 fragmentation reaction Methods 0.000 description 2
- 238000010191 image analysis Methods 0.000 description 2
- 238000003709 image segmentation Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000005422 blasting Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005188 flotation Methods 0.000 description 1
- 238000003706 image smoothing Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 238000001303 quality assessment method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
- G06T7/0008—Industrial image inspection checking presence/absence
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10024—Color image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30108—Industrial image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30168—Image quality inspection
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Processing (AREA)
Abstract
本发明属于图像处理技术领域,公开了一种岩石运输带上不良块度图像的实时剔除方法,包括:获取一幅RGB块度图像;计算对应缩小后的灰度图像进行平滑滤波后的灰度图像的平均灰度值和相对方差;设定第一正常平均灰度阈值,剔除一部分不良岩石块度图像;设定第一正常相对方差阈值,剔除一部分不良岩石块度图像;设定第二正常平均灰度阈值和第二正常相对方差阈值,剔除一部分不良岩石块度图像;得到对应的梯度图像,计算梯度平均值和梯度相对方差;设定平均梯度阈值和梯度相对方差阈值,剔除一部分不良岩石块度图像;获取下一幅块度图像,并重复上述步骤,能够快速准确地检测出运输带上的不良岩石块度图像。
The invention belongs to the technical field of image processing, and discloses a real-time elimination method of a bad blockiness image on a rock conveyor belt, comprising: acquiring an RGB blockiness image; calculating the grayscale corresponding to the reduced grayscale image after smoothing and filtering The average gray value and relative variance of the image; set the first normal average gray threshold to remove some bad rock blockiness images; set the first normal relative variance threshold to reject some bad rock blocky images; set the second normal The average gray threshold and the second normal relative variance threshold remove some bad rock blockiness images; get the corresponding gradient image, calculate the gradient average value and gradient relative variance; set the average gradient threshold and gradient relative variance threshold to remove some bad rocks Looseness image: Obtain the next piece of lumpiness image and repeat the above steps to quickly and accurately detect bad rock lumpiness images on the conveyor belt.
Description
技术领域technical field
本发明属于图像处理技术领域,尤其涉及一种岩石运输带上不良块度图像的实时剔除方法,适用于运输带上运动岩石块度图像的在线检测和分析。The invention belongs to the technical field of image processing, and in particular relates to a real-time elimination method for bad lumpiness images on a rock conveyor belt, which is suitable for online detection and analysis of moving rock lumpiness images on the conveyor belt.
背景技术Background technique
在采石和矿业工程中,对岩石块度的尺寸分布的测量是非常重要的。石料就是自然岩块和爆破及机械破碎的岩块的混合体,主要是用于建筑,公路,铁路和大坝等。为了判断石料的质量,对石料颗粒的尺寸和形状参数进行估算是必要的。石料的尺寸分布不但是用来评估产品质量的一个数据,而且还是调整破碎机或爆破生产的重要信息,例如:在采石生产中,调节破碎机间隙及在矿业工程中调节打孔的孔径等。破碎机通常被设定用来生产某个严格指定的相对较窄尺寸范围内的石料,比如从16mm到30mm。通常破碎机操作的一个主要指标就是平均尺寸。在自动破碎控制系统中,从实时系统发回的包括平均石料尺寸的反馈信号,就显示了流水线上破碎过程的实际进展。在实际应用中,从破碎机出来的破碎颗粒在一条传送带上传输,在其上方放置一个CCD(charge coupled device,电荷耦合器件)摄像头向下拍摄,然后用图像处理、分割和分析对获取的图像中的颗粒进行测量。In quarrying and mining engineering, the measurement of the size distribution of rock lumpiness is very important. Stone is a mixture of natural rock blocks and blasted and mechanically broken rock blocks, which are mainly used for construction, roads, railways and dams. In order to judge the quality of the stone, it is necessary to estimate the size and shape parameters of the stone particles. The size distribution of stone is not only a data used to evaluate product quality, but also important information for adjusting crusher or blasting production, for example: in quarrying production, adjusting the gap of crusher and adjusting the hole diameter of drilling in mining engineering, etc. . Crusher is usually set to produce stone within a strictly specified relatively narrow size range, such as from 16mm to 30mm. Usually a major indicator of crusher operation is the average size. In the automatic crushing control system, the feedback signal from the real-time system, including the average stone size, shows the actual progress of the crushing process on the line. In practical applications, the broken particles from the crusher are transported on a conveyor belt, and a CCD (charge coupled device, charge-coupled device) camera is placed above it to shoot downwards, and then image processing, segmentation and analysis are used to analyze the acquired images. Particles in are measured.
现有技术主要是用图像处理和分析及计算机视觉技术对复杂岩石块度图像进行高速和高精度的处理和分析,为提高采矿和选矿的生产流水线上的自动监测和控制水平奠定新的应用基础。矿岩块度图像是最为复杂的多物体图像,因为矿岩块度的颜色、粒度尺寸、形状、粗糙度、三维结构等特性使图像处理、分析和描述远远难于其它的粒度物体图像,所以在这方面进行模式识别、图像分析及机器视觉是有重要意义的。问题是:图像变化太大太频繁致使图像分割产生错误,从而导致错误的测量和分析结果,为了克服这个问题,对处理的图像应该是有选择性的,尽量在图像处理前,去掉那些质量低下的图像,例如:无岩石块度图像,雨水运输带反光照成的白板,由运动抖动照成的模糊图像等等。The existing technology mainly uses image processing and analysis and computer vision technology to process and analyze complex rock mass images at high speed and high precision, laying a new application foundation for improving the automatic monitoring and control level of mining and beneficiation production lines . The ore lumpiness image is the most complex multi-object image, because the color, grain size, shape, roughness, three-dimensional structure and other characteristics of the ore lumpiness make image processing, analysis and description far more difficult than other granular object images, so In this regard, pattern recognition, image analysis and machine vision are of great significance. The problem is: the image changes too much and too frequently, resulting in errors in image segmentation, resulting in wrong measurement and analysis results. In order to overcome this problem, the processed images should be selective, and try to remove those low-quality images before image processing For example: images without rock massiness, whiteboards made of reflections from rainwater conveyor belts, blurred images made of motion jitter, etc.
发明内容Contents of the invention
针对上述问题,本发明的目的在于提供一种实时剔除运输带上不良岩石块度图像的方法,能够快速准确地检测出运输带上的不良岩石块度图像,以解决图像质量评定和剔除的问题,从而能够服务于矿山及采石场上的工业流水生产线上的监测和控制。In view of the above problems, the purpose of the present invention is to provide a real-time method for removing bad rock lumpiness images on the conveyor belt, which can quickly and accurately detect bad rock lumpiness images on the conveyor belt, so as to solve the problems of image quality assessment and elimination , so that it can serve the monitoring and control of the industrial assembly line in mines and quarries.
为达到上述目的,本发明采用如下技术方案予以实现。In order to achieve the above object, the present invention adopts the following technical solutions to achieve.
一种岩石运输带上不良块度图像的实时剔除方法,所述方法包括如下步骤:A real-time method for removing bad lumpiness images on a rock conveyor belt, said method comprising the steps of:
步骤1,获取运输带上运动岩石的一幅块度图像,所述一幅块度图像为RGB图像;Step 1, obtaining a blockiness image of moving rocks on the conveyor belt, the blockiness image is an RGB image;
步骤2,将所述RGB图像转换为对应的灰度图像,并对所述灰度图像按照预先设定的比例因子进行缩小,得一幅缩小后的灰度图像;Step 2, converting the RGB image into a corresponding grayscale image, and reducing the grayscale image according to a preset scale factor to obtain a reduced grayscale image;
步骤3,对所述一幅缩小后的灰度图像进行平滑滤波,得到平滑滤波后的灰度图像;Step 3, smoothing and filtering the reduced grayscale image to obtain a smoothed and filtered grayscale image;
步骤4,计算所述平滑滤波后的灰度图像的平均灰度值和相对方差;Step 4, calculating the average gray value and relative variance of the gray image after smoothing and filtering;
步骤5,设定第一正常平均灰度阈值,若所述平滑滤波后的灰度图像的平均灰度值小于或者等于所述第一正常平均灰度阈值,则该平滑滤波后的灰度图像对应的运输带上运动岩石的该幅块度图像为不良岩石块度图像,将该不良岩石块度图像进行剔除;Step 5, setting the first normal average grayscale threshold, if the average grayscale value of the smooth-filtered grayscale image is less than or equal to the first normal average grayscale threshold, then the smooth-filtered grayscale image The piece size image of the moving rock on the corresponding conveyor belt is a bad rock piece size image, and the bad rock piece size image is removed;
若所述平滑滤波后的灰度图像的平均灰度值大于所述第一正常平均灰度阈值,则继续执行步骤6;If the average gray value of the smooth-filtered gray image is greater than the first normal average gray threshold, proceed to step 6;
步骤6,设定第一正常相对方差阈值,若所述平滑滤波后的灰度图像的相对方差小于或者等于所述第一正常相对方差阈值,则该平滑滤波后的灰度图像对应的运输带上运动岩石的该幅块度图像为不良岩石块度图像,将该不良岩石块度图像进行剔除;Step 6, setting the first normal relative variance threshold, if the relative variance of the smooth-filtered grayscale image is less than or equal to the first normal relative variance threshold, then the conveyor belt corresponding to the smooth-filtered grayscale image The piece size image of the moving rock is a bad rock piece size image, and the bad rock piece size image is removed;
若所述平滑滤波后的灰度图像的相对方差大于所述第一正常相对方差阈值,则继续执行步骤7;If the relative variance of the smooth-filtered grayscale image is greater than the first normal relative variance threshold, proceed to step 7;
步骤7,设定第二正常平均灰度阈值和第二正常相对方差阈值,若所述平滑滤波后的灰度图像的平均灰度值小于或者等于所述第二正常平均灰度阈值,且所述平滑滤波后的灰度图像的相对方差小于或者等于所述第二正常相对方差阈值,则该平滑滤波后的灰度图像对应的运输带上运动岩石的该幅块度图像为不良岩石块度图像,将该不良岩石块度图像进行剔除;否则,继续执行步骤8;Step 7, setting a second normal average grayscale threshold and a second normal relative variance threshold, if the average grayscale value of the smooth-filtered grayscale image is less than or equal to the second normal average grayscale threshold, and the If the relative variance of the smooth-filtered grayscale image is less than or equal to the second normal relative variance threshold, the blockiness image of the moving rock on the conveyor belt corresponding to the smooth-filtered grayscale image is bad rock blockiness image, remove the bad rock blockiness image; otherwise, proceed to step 8;
步骤8,根据所述平滑滤波后的灰度图像,得到对应的梯度图像,计算所述梯度图像的梯度平均值和梯度相对方差;Step 8, obtaining a corresponding gradient image according to the grayscale image after smoothing and filtering, and calculating the average gradient value and the relative gradient variance of the gradient image;
步骤9,设定平均梯度阈值和梯度相对方差阈值,若所述梯度图像的梯度平均值小于或者等于所述梯度阈值,且所述梯度图像的梯度相对方差小于或者等于所述梯度相对方差阈值,则该梯度图像对应的运输带上运动岩石的该幅块度图像为不良岩石块度图像,将该不良岩石块度图像进行剔除;Step 9, setting the average gradient threshold and the gradient relative variance threshold, if the average gradient value of the gradient image is less than or equal to the gradient threshold, and the gradient relative variance of the gradient image is less than or equal to the gradient relative variance threshold, Then the piece size image of the moving rock on the conveyor belt corresponding to the gradient image is a bad rock piece size image, and the bad rock piece size image is removed;
步骤10,获取运输带上运动岩石的下一幅块度图像,并依次重复执行步骤2至步骤9,从而实时剔除运输带上运动岩石的不良块度图像。Step 10, acquire the next blockiness image of the moving rock on the conveyor belt, and repeat step 2 to step 9 in sequence, so as to remove the bad blockiness image of the moving rock on the conveyor belt in real time.
本发明的目的是进行动态图像的质量评价,能够快速去除不良质量的岩石块度图像,从而保证下一步要进行的处理与分割分析的图像有好的质量。若没有这一过程,不良质量的图像会给后续的图像处理造成困难及分析误差,从而使实时在线准确检测结果不能得到保证。该方法为了适应实时处理的需要,分几个步骤来进行图像的质量分析,步骤与步骤之间避免了重复的计算,分析和计算方法以求尽可能简单有效,避免了运用复杂的计算和分析。非常适用于矿岩块度生产现场的应用,也很容易扩展到其他类似的在线颗粒检测,如:动态的浮选气泡,运输带上的木质碎片,运动的粮食颗粒及水果等图像的质量分析与检测。The purpose of the present invention is to evaluate the quality of dynamic images, and can quickly remove bad-quality rock blockiness images, so as to ensure that the images to be processed and segmented and analyzed in the next step have good quality. Without this process, poor-quality images will cause difficulties and analysis errors in subsequent image processing, so that real-time online accurate detection results cannot be guaranteed. In order to meet the needs of real-time processing, this method analyzes the image quality in several steps, avoiding repeated calculations between steps, and the analysis and calculation methods are as simple and effective as possible, avoiding the use of complex calculations and analysis . It is very suitable for the application of ore lumpiness production site, and it is also easy to expand to other similar online particle detection, such as: dynamic flotation air bubbles, woody debris on the conveyor belt, moving grain particles and fruit image quality analysis with detection.
附图说明Description of drawings
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings that need to be used in the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only These are some embodiments of the present invention. Those skilled in the art can also obtain other drawings based on these drawings without creative work.
图1为本发明实施例提供的一种岩石运输带上不良块度图像的实时剔除方法的流程示意图。Fig. 1 is a schematic flowchart of a real-time removal method for bad blockiness images on a rock conveyor belt provided by an embodiment of the present invention.
具体实施方式detailed description
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only some, not all, embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without making creative efforts belong to the protection scope of the present invention.
本发明实施例提供一种岩石运输带上不良块度图像的实时剔除方法,如图1所示,所述方法包括如下步骤:An embodiment of the present invention provides a method for real-time removal of bad lumpiness images on a rock conveyor belt, as shown in FIG. 1 , the method includes the following steps:
步骤1,获取运输带上运动岩石的一幅块度图像,所述一幅块度图像为RGB图像。Step 1, acquiring a blockiness image of moving rocks on the conveyor belt, where the blockiness image is an RGB image.
一般来说,为了长期实时地获取运输带上运动岩石的块度图像,运输带的上方需架设一台普通的CCD相机获取运输带上运动的岩石的块度图像,为了避免外来光照的不均及不稳定,也为了防尘防雨雪等,需要架设一个上面及侧面封闭的光照箱,箱内安装均匀的光源(灯)及CCD摄像头,为了清晰地抓取运动(一般运动速度为2-3米/秒)块度图像,CCD相机应该可以设置取像参数,既常说的快门和光圈,获取的图像可以通过连接的图像板传输到电脑进行后续图像处理。Generally speaking, in order to obtain long-term real-time fragmentation images of moving rocks on the conveyor belt, an ordinary CCD camera needs to be installed above the conveyor belt to obtain fragmentation images of rocks moving on the conveyor belt. In order to avoid uneven external light And instability, and in order to prevent dust, rain and snow, etc., it is necessary to set up a light box with a closed top and side, and install a uniform light source (light) and a CCD camera in the box, in order to clearly capture the movement (generally the movement speed is 2- 3 m/s) block size image, the CCD camera should be able to set the image acquisition parameters, which are often referred to as the shutter and aperture, and the acquired image can be transmitted to the computer through the connected image board for subsequent image processing.
步骤2,将所述RGB图像转换为对应的灰度图像,并对所述灰度图像按照预先设定的比例因子进行缩小,得一幅缩小后的灰度图像。Step 2, converting the RGB image into a corresponding grayscale image, and reducing the grayscale image according to a preset scale factor to obtain a reduced grayscale image.
步骤2具体包括:Step 2 specifically includes:
由于岩石块度彩色不明显,为了减少计算量,将所述RGB图像转换为对应的灰度图像:Since the color of the rock mass is not obvious, in order to reduce the amount of calculation, the RGB image is converted into a corresponding grayscale image:
f(x,y)灰=(f(x,y)R+f(x,y)G+f(x,y)B)/3f(x,y) = (f(x,y) R +f(x,y) G +f(x,y) B )/3
其中,f(x,y)R、f(x,y)G、f(x,y)B分别表示RGB图像中位于(x,y)处的像素的红、绿、蓝像素值,f(x,y)灰表示灰度图像中位于(x,y)处的像素的灰度值;x∈(1,...,2×M),y∈(1,...,2×N),2×M为灰度图像宽度维的总像素个数,2×N为灰度图像高度维的总像素个数;Among them, f(x, y) R , f(x, y) G , f(x, y) B respectively represent the red, green and blue pixel values of the pixel located at (x, y) in the RGB image, f( x, y) gray represents the gray value of the pixel located at (x, y) in the gray image; x∈(1,...,2×M), y∈(1,...,2×N ), 2×M is the total number of pixels in the width dimension of the grayscale image, and 2×N is the total number of pixels in the height dimension of the grayscale image;
为了消除黑点和亮点噪声及进一步减少后续处理的工作量,将灰度图像进行缩小;预先设定的比例因子为1/4,则按照所述比例因子对所述灰度图像进行缩小,每四个相邻的像素取灰度平均值作为缩小后的灰度图像在对应位置处的灰度值,从而缩小后的灰度图像在宽度维的总像素个数为M,在高度维的总像素个数为N。In order to eliminate the noise of black spots and bright spots and further reduce the workload of subsequent processing, the grayscale image is reduced; the preset scale factor is 1/4, then the grayscale image is reduced according to the scale factor, every The average gray value of four adjacent pixels is taken as the gray value of the reduced gray image at the corresponding position, so that the total number of pixels of the reduced gray image in the width dimension is M, and the total number of pixels in the height dimension is The number of pixels is N.
步骤3,对所述一幅缩小后的灰度图像进行平滑滤波,得到平滑滤波后的灰度图像。Step 3, smoothing and filtering the reduced grayscale image to obtain a smoothed and filtered grayscale image.
如果按常规的图像平滑算法(如邻域平均法或高斯平滑法等)对缩小后的灰度图像进行平滑,可能会把微弱的边界平滑掉,故此,采取分数阶积分平滑方法来去除所述缩小后的灰度图像中的噪声,但为了既保持图像的平滑又能保持块度间的边界不被丢失,对传统的滤波器模板进行了改进:If the reduced grayscale image is smoothed by a conventional image smoothing algorithm (such as the neighborhood average method or Gaussian smoothing method, etc.), the weak boundary may be smoothed away. Therefore, the fractional integral smoothing method is used to remove the The noise in the reduced grayscale image, but in order to keep the smoothness of the image and keep the boundary between blocks from being lost, the traditional filter template is improved:
采用分数阶积分平滑方法对所述缩小后的灰度图像进行平滑滤波,所使用的5×5模板的滤波器系数h为:The fractional order integral smoothing method is used to smooth and filter the reduced grayscale image, and the filter coefficient h of the 5×5 template used is:
则平滑滤波后的灰度图像在(x1,y1)处的灰度值F(x1,y1)为:F(x1,y1)=(f(x1,y1)*h)/8,且平滑滤波后的灰度图像的左边两列、右边两列、上边两行以及下边两行像素的灰度值与缩小后的灰度图像在对应位置处的灰度值相同;Then the gray-scale value F(x 1 , y 1 ) of the smooth-filtered gray-scale image at (x 1 , y 1 ) is: F(x 1 , y 1 )=(f(x 1 , y 1 )* h)/8, and the gray value of pixels in the left two columns, right two columns, upper two rows and lower two rows of the smoothed gray image is the same as the gray value of the corresponding position of the reduced gray image ;
需要说明的是,平滑滤波后的灰度图像在(x1,y1)处的灰度值F(x1,y1)为:F(x1,y1)=(f(x1,y1)*h)/8,该式所表示的滤波运算是指,获取平滑滤波后的灰度图像以(x1,y1)处的像素为中心的25个像素,从而这25个像素与滤波器系数分别卷积,从而得到(x1,y1)处的灰度值F(x1,y1)。It should be noted that the gray value F(x 1 , y 1 ) of the smoothed gray image at (x 1 , y 1 ) is: F(x 1 , y 1 )=(f(x 1 , y 1 )*h)/8, the filtering operation represented by this formula refers to obtaining 25 pixels centered on the pixel at (x 1 , y 1 ) of the grayscale image after smoothing and filtering, so that the 25 pixels Convolute with the filter coefficients respectively, so as to obtain the gray value F(x 1 , y 1 ) at (x 1 , y 1 ).
其中,f(x1,y1)表示缩小后的灰度图像的在(x1,y1)处的灰度值,且x1∈(0,...,M),y1∈(0,...,N)。Among them, f(x 1 , y 1 ) represents the gray value at (x 1 , y 1 ) of the reduced gray-scale image, and x 1 ∈ (0, ..., M), y 1 ∈ ( 0,...,N).
步骤4,计算所述平滑滤波后的灰度图像的平均灰度值和相对方差。Step 4, calculating the average gray value and relative variance of the smoothed and filtered gray image.
步骤4具体包括:Step 4 specifically includes:
计算所述平滑滤波后的灰度图像的平均灰度值y和相对方差S相:Calculate the average gray value y and the relative variance S phase of the smooth filtered gray image:
S相=(s/v)×100S phase =(s/v)×100
其中,F(x1,y1)表示平滑滤波后的灰度图像在(x1,y1)处的灰度值,x1∈(0,...,M),y1∈(0,...,N),平滑滤波后的灰度图像在宽度维的总像素个数为M,在高度维的总像素个数为N,S表示平滑滤波后的灰度图像的方差。Among them, F(x 1 , y 1 ) represents the gray value of the smooth-filtered gray-scale image at (x 1 , y 1 ), x 1 ∈ (0,..., M), y 1 ∈ (0 ,..., N), the total number of pixels in the width dimension of the smoothed and filtered grayscale image is M, the total number of pixels in the height dimension is N, and S represents the variance of the smoothed and filtered grayscale image.
步骤5,设定第一正常平均灰度阈值,若所述平滑滤波后的灰度图像的平均灰度值小于或者等于所述第一正常平均灰度阈值,则该平滑滤波后的灰度图像对应的运输带上运动岩石的该幅块度图像为不良岩石块度图像,将该不良岩石块度图像进行剔除;Step 5, setting the first normal average grayscale threshold, if the average grayscale value of the smooth-filtered grayscale image is less than or equal to the first normal average grayscale threshold, then the smooth-filtered grayscale image The piece size image of the moving rock on the corresponding conveyor belt is a bad rock piece size image, and the bad rock piece size image is removed;
若所述平滑滤波后的灰度图像的平均灰度值大于所述第一正常平均灰度阈值,则继续执行步骤6。If the average gray value of the smooth-filtered gray image is greater than the first normal average gray threshold, proceed to step 6.
一般来说,由于受光照及运动速度的影响,运动岩石的块度图像质量会低下,而当块度图像平均灰度值很低时,图像中的岩石块度会模糊一片或暗灰色甚至黑色一片(如:无块度图像),无法进行后续的图像处理和分割,所以为了避免选取这种图像,需要根据现场情况(岩石色度,尺寸大小等)设定第一正常平均灰度阈值。Generally speaking, due to the influence of light and motion speed, the image quality of the blockiness of moving rocks will be low, and when the average gray value of the blockiness image is very low, the blockiness of rocks in the image will be blurred or dark gray or even black One piece (such as: non-blockiness image) cannot be followed by image processing and segmentation, so in order to avoid selecting such an image, it is necessary to set the first normal average gray threshold according to the site conditions (rock color, size, etc.).
设定第一正常平均灰度阈值具体为:人为选取运动岩石的一幅优质图像,得到该优质图像的平均灰度值,设定该优质图像的平均灰度值的30%作为第一正常平均灰度阈值。Setting the first normal average grayscale threshold is specifically: artificially select a high-quality image of moving rocks, obtain the average grayscale value of the high-quality image, and set 30% of the average grayscale value of the high-quality image as the first normal average Gray threshold.
步骤6,设定第一正常相对方差阈值,若所述平滑滤波后的灰度图像的相对方差小于或者等于所述第一正常相对方差阈值,则该平滑滤波后的灰度图像对应的运输带上运动岩石的该幅块度图像为不良岩石块度图像,将该不良岩石块度图像进行剔除;Step 6, setting the first normal relative variance threshold, if the relative variance of the smooth-filtered grayscale image is less than or equal to the first normal relative variance threshold, then the conveyor belt corresponding to the smooth-filtered grayscale image The piece size image of the moving rock is a bad rock piece size image, and the bad rock piece size image is removed;
若所述平滑滤波后的灰度图像的相对方差大于所述第一正常相对方差阈值,则继续执行步骤7。If the relative variance of the smooth-filtered grayscale image is greater than the first normal relative variance threshold, proceed to step 7.
尽管在步骤5中能够剔除一部分不良质量的图像,被选取的图像虽然平均灰度值较高,但也还有些图像是低质量图像,例如:由运动速度引起的模糊的图像,或是由突然的强光引起的白板图像等,这些图像的方差都会很低,因此可以用图像的方差来判断图像中是否有块度:如果方差值高,证明图像灰度差异大,也就是块度较多。Although some bad-quality images can be eliminated in step 5, although the average gray value of the selected images is relatively high, there are still some low-quality images, such as: blurred images caused by motion speed, or blurred images caused by sudden Whiteboard images caused by strong light, etc., the variance of these images will be very low, so the variance of the image can be used to judge whether there is blockiness in the image: if the variance value is high, it proves that the gray level of the image has a large difference, that is, the blockiness is relatively large many.
问题是:在不同的光照条件下,会有不同的图像平灰度均值,而此时,同样的块度数目及大小也会导致差异较大的方差。为了避免这种难以统一判断的问题,本技术方案引用相对误差来判断图像的质量。The problem is: under different lighting conditions, there will be different mean values of image flat gray levels, and at this time, the same number and size of blocks will also lead to large differences in variance. In order to avoid such a problem that it is difficult to judge uniformly, the technical solution uses relative error to judge the image quality.
步骤6中,设定第一正常相对方差阈值具体为:人为选取运动岩石的一幅优质图像,得到该优质图像的相对方差,设定该优质图像的相对方差的40%作为第一正常相对方差阈值。In step 6, setting the first normal relative variance threshold is specifically: artificially select a high-quality image of the moving rock, obtain the relative variance of the high-quality image, and set 40% of the relative variance of the high-quality image as the first normal relative variance threshold.
步骤7,设定第二正常平均灰度阈值和第二正常相对方差阈值,若所述平滑滤波后的灰度图像的平均灰度值小于或者等于所述第二正常平均灰度阈值,且所述平滑滤波后的灰度图像的相对方差小于或者等于所述第二正常相对方差阈值,则该平滑滤波后的灰度图像对应的运输带上运动岩石的该幅块度图像为不良岩石块度图像,将该不良岩石块度图像进行剔除;否则,继续执行步骤8。Step 7, setting a second normal average grayscale threshold and a second normal relative variance threshold, if the average grayscale value of the smooth-filtered grayscale image is less than or equal to the second normal average grayscale threshold, and the If the relative variance of the smooth-filtered grayscale image is less than or equal to the second normal relative variance threshold, the blockiness image of the moving rock on the conveyor belt corresponding to the smooth-filtered grayscale image is bad rock blockiness image, remove the bad rock blockiness image; otherwise, proceed to step 8.
进一步的,步骤5和步骤6均为以单一指标来判断图像的质量,但有些情况需要结合两种参数来判断才能可靠。Further, step 5 and step 6 both use a single index to judge the quality of the image, but in some cases it is necessary to combine two parameters to make the judgment reliable.
步骤7中,设定第二正常平均灰度阈值和第二正常相对方差阈值具体为:In step 7, setting the second normal average gray threshold and the second normal relative variance threshold are specifically:
人为选取运动岩石的一幅优质图像,得到该优质图像的平均灰度值和相对方差,设定该优质图像的平均灰度值的45%作为第二正常平均灰度阈值,设定该优质图像的相对方差的60%作为第二正常相对方差阈值。Artificially select a high-quality image of moving rocks, obtain the average gray value and relative variance of the high-quality image, set 45% of the average gray value of the high-quality image as the second normal average gray threshold, and set the high-quality image 60% of the relative variance of is used as the second normal relative variance threshold.
步骤8,根据所述平滑滤波后的灰度图像,得到对应的梯度图像,计算所述梯度图像的梯度平均值和梯度相对方差。Step 8: Obtain the corresponding gradient image according to the smoothed and filtered grayscale image, and calculate the gradient average value and the gradient relative variance of the gradient image.
需要说明的是,通过上述步骤5、步骤6和步骤7的剔除,70%-80%左右的不良质量图像会被筛选掉,剩下的20%-30%不良质量的图像可以通过步骤8和步骤9中的准则来剔除。It should be noted that through the elimination of the above steps 5, 6 and 7, about 70%-80% of poor quality images will be screened out, and the remaining 20%-30% of poor quality images can be passed through steps 8 and 7. Criteria in step 9 to eliminate.
步骤8具体包括:Step 8 specifically includes:
对所述平滑滤波后的灰度图像进行一阶微分得到对应的梯度图像;计算所述梯度图像的梯度平均值V1和梯度相对方差S1相:Carry out first-order differentiation to the grayscale image after the smooth filter to obtain the corresponding gradient image; calculate the gradient mean value V of the gradient image and the gradient relative variance S phase :
S1相=(s1/v1)×100S 1 phase = (s 1 /v 1 )×100
其中,G(x2,y2)表示梯度图像在(x2,y2)处的梯度值,x2∈(0,...,M),y2∈(0,...,N),梯度图像在宽度维的总像素个数为M,在高度维的总像素个数为N,S1表示梯度图像的方差。Among them, G(x 2 , y 2 ) represents the gradient value of the gradient image at (x 2 , y 2 ), x 2 ∈ (0, ..., M), y 2 ∈ (0, ..., N ), the total number of pixels of the gradient image in the width dimension is M, the total number of pixels in the height dimension is N, and S 1 represents the variance of the gradient image.
步骤9,设定平均梯度阈值和梯度相对方差阈值,若所述梯度图像的梯度平均值小于或者等于所述梯度阈值,且所述梯度图像的梯度相对方差小于或者等于所述梯度相对方差阈值,则该梯度图像对应的运输带上运动岩石的该幅块度图像为不良岩石块度图像,将该不良岩石块度图像进行剔除。Step 9, setting the average gradient threshold and the gradient relative variance threshold, if the average gradient value of the gradient image is less than or equal to the gradient threshold, and the gradient relative variance of the gradient image is less than or equal to the gradient relative variance threshold, Then the blockiness image of the moving rock on the conveyor belt corresponding to the gradient image is a bad rock blockiness image, and the bad rock blockiness image is removed.
步骤9中,设定平均梯度阈值和梯度相对方差阈值具体为:In step 9, the average gradient threshold and the gradient relative variance threshold are set as follows:
人为选取运动岩石的一幅优质图像,得到该优质图像对应的梯度图像的平均梯度阈值和梯度相对方差阈值,设定该优质图像对应的梯度图像的平均梯度阈值的50%作为平均梯度阈值,设定该优质图像对应的梯度图像的梯度相对方差阈值的60%作为梯度相对方差阈值。Artificially select a high-quality image of moving rocks, obtain the average gradient threshold and gradient relative variance threshold of the gradient image corresponding to the high-quality image, set 50% of the average gradient threshold of the gradient image corresponding to the high-quality image as the average gradient threshold, set 60% of the gradient relative variance threshold of the gradient image corresponding to the high-quality image is determined as the gradient relative variance threshold.
步骤10,获取运输带上运动岩石的下一幅块度图像,并依次重复执行步骤2至步骤9,从而实时剔除运输带上运动岩石的不良块度图像。Step 10, acquire the next blockiness image of the moving rock on the conveyor belt, and repeat step 2 to step 9 in sequence, so as to remove the bad blockiness image of the moving rock on the conveyor belt in real time.
最终剩下的块度图像是质量较好的图像,只要后续的处理算法适合于图像的类别,就不会对后续的处理产生大的误差,从而可以进行后续的图像分割和分析。The final remaining blocky image is an image with better quality. As long as the subsequent processing algorithm is suitable for the category of the image, there will be no large error in the subsequent processing, so that subsequent image segmentation and analysis can be performed.
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。Those of ordinary skill in the art can understand that all or part of the steps to realize the above method embodiments can be completed by hardware related to program instructions, and the aforementioned programs can be stored in computer-readable storage media. When the program is executed, the execution includes The steps of the above-mentioned method embodiments; and the aforementioned storage medium includes: ROM, RAM, magnetic disk or optical disk and other various media that can store program codes.
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。The above is only a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Anyone skilled in the art can easily think of changes or substitutions within the technical scope disclosed in the present invention. Should be covered within the protection scope of the present invention. Therefore, the protection scope of the present invention should be determined by the protection scope of the claims.
Claims (9)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710353593.3A CN107220946B (en) | 2017-05-18 | 2017-05-18 | Real-time removing method for bad block degree images on rock conveyer belt |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710353593.3A CN107220946B (en) | 2017-05-18 | 2017-05-18 | Real-time removing method for bad block degree images on rock conveyer belt |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107220946A true CN107220946A (en) | 2017-09-29 |
CN107220946B CN107220946B (en) | 2020-06-23 |
Family
ID=59944259
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710353593.3A Expired - Fee Related CN107220946B (en) | 2017-05-18 | 2017-05-18 | Real-time removing method for bad block degree images on rock conveyer belt |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107220946B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109241867A (en) * | 2018-08-16 | 2019-01-18 | 清能艾科(深圳)能源技术有限公司 | Using the method and device of intelligent algorithm identification digital cores image |
CN110689550A (en) * | 2019-10-12 | 2020-01-14 | 嘉应学院 | Efficient and automatic screening system and method for lumbar vertebra sagittal plane CT (computed tomography) images |
WO2023279556A1 (en) * | 2021-07-09 | 2023-01-12 | 福建南方路面机械股份有限公司 | Aggregate particle size sampling and monitoring method |
CN116687442A (en) * | 2023-08-08 | 2023-09-05 | 汕头市超声仪器研究所股份有限公司 | A fetal face imaging method based on three-dimensional volume data |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101567044A (en) * | 2009-05-22 | 2009-10-28 | 北京大学 | Method for detecting quality of human face image |
CN102567996A (en) * | 2012-01-05 | 2012-07-11 | 中国地质大学(武汉) | Remote sensing satellite image quality evaluation method |
CN104881878A (en) * | 2015-06-12 | 2015-09-02 | 武汉大学 | Image quality evaluation method based on multi-scale gradient difference entropies |
-
2017
- 2017-05-18 CN CN201710353593.3A patent/CN107220946B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101567044A (en) * | 2009-05-22 | 2009-10-28 | 北京大学 | Method for detecting quality of human face image |
CN102567996A (en) * | 2012-01-05 | 2012-07-11 | 中国地质大学(武汉) | Remote sensing satellite image quality evaluation method |
CN104881878A (en) * | 2015-06-12 | 2015-09-02 | 武汉大学 | Image quality evaluation method based on multi-scale gradient difference entropies |
Non-Patent Citations (2)
Title |
---|
王卫星: "岩石颗粒计算机视觉系统", 《有色金属(矿山部分)》 * |
王志明: "无参考图像质量评价综述", 《自动化学报 》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109241867A (en) * | 2018-08-16 | 2019-01-18 | 清能艾科(深圳)能源技术有限公司 | Using the method and device of intelligent algorithm identification digital cores image |
CN110689550A (en) * | 2019-10-12 | 2020-01-14 | 嘉应学院 | Efficient and automatic screening system and method for lumbar vertebra sagittal plane CT (computed tomography) images |
CN110689550B (en) * | 2019-10-12 | 2022-06-21 | 嘉应学院 | A system and method for efficient automatic screening of lumbar sagittal CT images |
WO2023279556A1 (en) * | 2021-07-09 | 2023-01-12 | 福建南方路面机械股份有限公司 | Aggregate particle size sampling and monitoring method |
CN116687442A (en) * | 2023-08-08 | 2023-09-05 | 汕头市超声仪器研究所股份有限公司 | A fetal face imaging method based on three-dimensional volume data |
Also Published As
Publication number | Publication date |
---|---|
CN107220946B (en) | 2020-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110286126B (en) | Wafer surface defect regional detection method based on visual image | |
CN103051872B (en) | Based on the belt deflection detection method of Edge extraction | |
CN107220946B (en) | Real-time removing method for bad block degree images on rock conveyer belt | |
CN103413290B (en) | Multiple features and the multi-level ore grain size image partition method combined | |
CN104504388B (en) | A kind of pavement crack identification and feature extraction algorithm and system | |
CN107314957B (en) | A method for measuring the size distribution of rock fragmentation | |
CN109598715B (en) | Material granularity online detection method based on machine vision | |
Schmugge et al. | Detection of cracks in nuclear power plant using spatial-temporal grouping of local patches | |
CN102663726B (en) | Method and device for material segmentation of material accumulation images | |
CN115345885A (en) | Method for detecting appearance quality of metal fitness equipment | |
CN104867144B (en) | IC element welding point defect detection methods based on mixed Gauss model | |
CN103886290B (en) | Coal stream video detecting method under a kind of complex background | |
CN110414341A (en) | A measurement and statistics method of open-pit mine explosion pile based on unmanned aerial vehicle technology | |
WO2020228111A1 (en) | X-ray image-based spasm detection method for steel cable core conveyor belt connector | |
CN103824304A (en) | Method for performing fault diagnosis on ores on conveying belt during conveying process | |
CN103226106A (en) | Machine vision based bearing dust cap quality monitoring system | |
CN115535525B (en) | System and method for detecting longitudinal tearing of conveyor belt based on image matching | |
CN104063866A (en) | Method for detecting granularity in ore transmission process | |
CN108711149A (en) | Ore-rock particle size detection method based on image procossing | |
CN107784646B (en) | Adaptive detection method for road aggregate | |
CN108470173B (en) | Ore particle segmentation method | |
CN106447673A (en) | Chip pin extraction method under non-uniform illumination condition | |
CN114212452A (en) | Coal flow detection method and energy-saving control system based on laser assistance and image processing | |
CN103927553A (en) | Coal and rock recognition method based on multi-scale micro-lamination and contrast ratio joint distribution | |
CN103914829B (en) | Method for detecting edge of noisy image |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200623 |
|
CF01 | Termination of patent right due to non-payment of annual fee |