CN111179286A - FCM rotary kiln flame image segmentation method based on energy noise detection - Google Patents

FCM rotary kiln flame image segmentation method based on energy noise detection Download PDF

Info

Publication number
CN111179286A
CN111179286A CN202010000820.6A CN202010000820A CN111179286A CN 111179286 A CN111179286 A CN 111179286A CN 202010000820 A CN202010000820 A CN 202010000820A CN 111179286 A CN111179286 A CN 111179286A
Authority
CN
China
Prior art keywords
image
rotary kiln
flame image
algorithm
energy
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
Application number
CN202010000820.6A
Other languages
Chinese (zh)
Other versions
CN111179286B (en
Inventor
易灵芝
刘宁
刘江永
刘文翰
孙颢一
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xiangtan University
Original Assignee
Xiangtan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xiangtan University filed Critical Xiangtan University
Priority to CN202010000820.6A priority Critical patent/CN111179286B/en
Publication of CN111179286A publication Critical patent/CN111179286A/en
Application granted granted Critical
Publication of CN111179286B publication Critical patent/CN111179286B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10016Video; Image sequence

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

The accurate judgment of the combustion state in the rotary kiln is one of the important steps in the production of the rotary kiln, and the accurate identification of the combustion condition of the flame of the rotary kiln can not only increase the production efficiency of the rotary kiln, save the cost and improve the production quality of products, but also ensure the personal safety of workers and reduce the discharge amount of pollutants in the production process. With the continuous and deep research on computer software technology, the problem of identifying the flame image of the rotary kiln by using the computer technology is one of the main modes of the current industrial production. The method is used for researching the identification of the flame image of the rotary kiln, analyzing the application of image segmentation in the confluent flame image of the rotary kiln and segmenting the flame image of the rotary kiln by using a common algorithm. The method aims at improving poor flame image segmentation effect of the traditional fuzzy clustering algorithm (FCM), the improved algorithm is applied to the flame image, and effectiveness of the improved algorithm and accuracy of flame image segmentation are obtained through experimental analysis and comparison.

Description

FCM rotary kiln flame image segmentation method based on energy noise detection
Technical Field
The identification technology of the rotary kiln has many problems, and how to improve the segmentation identification of the shot flame image by using an intelligent optimization algorithm has very important practical significance.
Background
With the rapid development of the computer industry, image segmentation has become one of the main methods for identifying the flame status during the operation of the rotary kiln. And selecting and extracting characteristic values of the flame image in the kiln by utilizing image segmentation, and then identifying the target flame by utilizing the segmented characteristic values so as to obtain the flame state in the rotary kiln. In recent years, many scholars at home and abroad deeply dig out the application of the algorithm to pattern recognition and recognize and divide the rotary kiln flame image by using image division. The traditional image segmentation method is low in speed and complex in calculation, cannot monitor the flame combustion state in the kiln in real time, and can only segment the acquired flame image at a later stage. However, due to the influence of the combustion environment in the rotary kiln, the state of the flame is easily interfered by dust in the kiln when shooting, and the obtained image is seriously polluted, so that the segmentation is difficult.
Aiming at the defects of the traditional FCM, scholars at home and abroad carry out various improvement researches on the FCM and make a series of progresses. ChenS.H and the like use spatial neighborhood information to constrain an objective function, and a series of algorithms such as FCM-S and the like are provided to improve the segmentation precision, but the algorithm efficiency is low. The Kangjiayin filters the image based on the self-adaptive weighted mean value and corrects the objective function, and provides an improved FCM considering pixel space information, but the algorithm ignores detail information in the image. By analyzing the sample spatial distribution characteristics and the relation thereof, the Xiaochun and the like provide an improved algorithm which gives consideration to image detail information by means of FCM based on spatial correlation, but the algorithm has a large calculation amount. Hocrown English and other algorithms combining space constrained fast FCM and Markov random fields to realize side scan sonar image segmentation, but image detail information is poor to keep. Guo F F and the like propose an adaptive FCM segmentation method based on local noise detection on the basis of documents, and the algorithm distinguishes noise by calculating the variance of neighborhood gray values, but introduces more parameters and easily loses the detailed information of images.
Compared with the image segmentation of other algorithms, the FCM (ENFCM) algorithm based on energy noise detection has certain advantages, has the advantages that an energy curve is introduced to achieve adaptivity of parameters, the target function is updated by using the distance between samples to update the membership degree and the clustering center to improve the spatial correlation of pixels, and the segmentation accuracy and the original information of the image are well reserved.
Disclosure of Invention
The invention provides an FCM algorithm based on energy noise detection on the basis of NDFCM and is used for realizing initial segmentation of a rotary kiln flame image. The NDFCM algorithm obtains a good anti-noise effect by combining with a spatial neighborhood constraint condition, but the method needs to consider and select a proper window to prevent the loss of image detail information when setting parameters. In order to overcome errors caused by a window in the image segmentation process, the FCM algorithm based on energy noise detection introduces data information of an energy curve compressed image by combining spatial information so as to realize accurate and efficient segmentation of a rotary kiln flame image.
The optimization algorithm comprises the following steps:
the method comprises the following steps: let I ═ lijWhere 1. ltoreq. i.ltoreq.m, 1. ltoreq. j.ltoreq.n is an image of size mxn, where lijIs the gray value of the image I at pixel location (I, j). L is the maximum gray value of the image I. Defining the d-order neighborhood system of (i, j) as
Figure BDA0002353314980000021
A second-order neighborhood system is utilized, i.e., (u, v) ∈ { (+ -1, 0), (0, + -1), (1, + -1), (-1, + -1) }. Two-dimensional matrix B generated from I imagel={bij,1≤i≤m,1≤j≤n},C={cijI is more than or equal to 1 and less than or equal to m, j is more than or equal to 1 and less than or equal to n, where cijWhen the energy curve is generated, different positions are adopted according to a d-order neighborhood system, and an energy curve function associated with each gray value in the image is designed by combining spatial context information, and the energy curve function is defined as follows:
Figure BDA0002353314980000022
Figure BDA0002353314980000023
wherein the content of the first and second substances,
Figure BDA0002353314980000024
is a constant term, ensuring an energy value ElGreater than 0. As can be seen from equation (1), the gray scale value for a pixel in image I is [ t ]1,t2]When l is equal to t1Here, a two-dimensional binary matrix B is generatedlMatrix BlThe pixel of the corresponding element in (1) is 1, and when increasing the value of/, only a small number of elements in the matrix change from 1 to-1, so that the energy function will increase, and when increasing to a certain value, the neighbor elements in the matrix change from 1 to-1.
Step two: delineation of pixel second order (N) at spatial locations2) The neighboring pixel (i, j), as shown in fig. 1.
The probability of a certain pixel being noise is calculated based on the energy curve as follows:
Figure BDA0002353314980000025
wherein, the lambda is an adjusting parameter,
Figure BDA0002353314980000026
the objective function of the ENFCM may be defined as follows:
Figure BDA0002353314980000027
wherein:
Figure BDA0002353314980000028
the function of space restriction is realized,
Figure BDA0002353314980000029
(djkrepresents a sample point xkAnd xjEuclidean distance) is a penalty factor representing the sample point xkSurrounding sample points and xkSpatial correlation of the classification results, ωikThe larger the value of (d), the stronger the spatial correlation.
Step three: according toMethod of FCM minimization, order
Figure BDA0002353314980000031
The following can be obtained:
Figure BDA0002353314980000032
Figure BDA0002353314980000033
from the above formula, it can be seen that in the introduction of the energy curve ElAnd ωikThereafter, the cluster centers are readjusted, the sample membership is indirectly changed by the change of the cluster centers in the iterative process, and ElAnd ωikAll taking into account the sample distribution around the sample point, El> 0, and ωikOmega when the sample is evenly distributedikIs a constant value.
As can be seen from the equations (5) and (6), after the energy curve and the parameter term are directly introduced into the objective function, only the cluster center is constrained. In order to improve the spatial correlation of the membership degree, the membership degree shown in formula (5) is redefined as follows by a weighting method:
Figure BDA0002353314980000034
wherein
Figure BDA0002353314980000035
Represents the mean of the degrees of membership.
Drawings
FIG. 1 is a diagram depicting a neighborhood window at a spatial location provided by the present invention.
FIG. 2 is a graph showing the SA variation of the segmentation accuracy of each algorithm in the case of gradual noise enhancement according to the present invention
FIG. 3 is a flow chart of a specific implementation of the ENFCM algorithm provided by the present invention.
Fig. 4 is an abstract drawing provided by the invention.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and embodiments.
In order to test the anti-noise stability of the ENFCM algorithm, noises (salt and pepper noises and Gaussian noises) with different intensities are added into a flame image of the rotary kiln respectively. In view of the effectiveness of the segmentation accuracy SA in measuring the anti-noise aspect of the algorithm, the anti-noise stability of each algorithm is analyzed by adopting the SA. As can be seen from fig. 2: with the gradual increase of the noise intensity, the segmentation accuracy SA of each algorithm is gradually reduced; since spatial correlation between pixels is not considered and noise is sensitive, SA of the conventional FCM rapidly decreases; FGFCM introduces a histogram method to enable gray levels to replace pixel points, and a correction sample method is adopted to accelerate the operation speed of the algorithm, but the algorithm precision is poor; the SCFMFCM algorithm improves the influence of an image on noise by combining MRF while compressing spatial information by using a histogram, and improves the accuracy of the algorithm; the SCMSFCM and the NDFCM combine pixel space information on the basis of a traditional FCM algorithm, and the segmentation accuracy SA is relatively slow to decline; the ENFCM introduces a space energy curve on the basis of the NDFCM, so that the space correlation of the algorithm is increased, and the segmentation quality of the algorithm is further improved. Therefore, although the noise intensity increase has a certain influence on the algorithm of the invention, the influence is relatively small, which shows that the algorithm has good anti-noise stability. Therefore, the method has a good segmentation effect on the rotary kiln image with high noise degree.
The ENFCM is realized by the following steps:
step 1: initializing parameters, including: the cluster number c (c is more than or equal to 2 and less than or equal to n), an update stop threshold epsilon, a fuzzy weighting index m and the like;
step 2: calculating an energy curve function E according to equation (1)l
And step 3: solving a clustering center and a membership degree according to the traditional FCM as the basis of a new algorithm so as to reduce the updating period of the next operation and shorten the operation time;
and 4, step 4: solving for sample x according to traditional FCM algorithmjAnd cluster center xkEuclidean function between, calculating ωik
step 5, calculating alpha by using the formula (3)i *Obtaining a sample xj
Step 6: calculating u according to equation (7)ik *
And 7: calculating the clustering center v of the sample class according to equation (6)ik
And 8: updating the membership and the cluster center when
Figure BDA0002353314980000041
The operation is stopped; otherwise, t is t +1, and the process jumps to Step 6.
The specific implementation flow chart of the ENFCM algorithm is shown in fig. 3.

Claims (1)

1. The invention provides an FCM algorithm based on energy noise detection on the basis of NDFCM and is used for realizing initial segmentation of a rotary kiln flame image. The NDFCM algorithm obtains a good anti-noise effect by combining with a spatial neighborhood constraint condition, but the method needs to consider and select a proper window to prevent the loss of image detail information when setting parameters. In order to overcome errors caused by a window in the image segmentation process, the FCM algorithm based on energy noise detection introduces data information of an energy curve compressed image by combining spatial information so as to realize accurate and efficient segmentation of a rotary kiln flame image.
The optimization algorithm comprises the following steps:
the method comprises the following steps: let I ═ lijWhere 1. ltoreq. i.ltoreq.m, 1. ltoreq. j.ltoreq.n is an image of size mxn, where lijIs the gray value of the image I at pixel location (I, j). L is the maximum gray value of the image I. Defining the d-order neighborhood system of (i, j) as
Figure FDA0002353314970000011
A second-order neighborhood system is utilized, i.e., (u, v) ∈ { (+ -1, 0), (0, + -1), (1, + -1), (-1, + -1) }. Two-dimensional matrix B generated from I imagel={bij,1≤i≤m,1≤j≤n},C={cijI is more than or equal to 1 and less than or equal to m, j is more than or equal to 1 and less than or equal to n, where cijWhen generating the energy curve, different positions are adopted according to a d-order neighborhood system, and the space is combinedThe context information designs an energy curve function associated with each gray value in the image, defined as follows:
Figure FDA0002353314970000012
Figure FDA0002353314970000013
wherein the content of the first and second substances,
Figure FDA0002353314970000014
is a constant term, ensuring an energy value ElGreater than 0. General formula (1)
It can be seen that for image I, the gray scale value for the pixel is at [ t ]1,t2]When l is equal to t1Here, a two-dimensional binary matrix B is generatedlMatrix BlThe pixel of the corresponding element in (1) is 1, and when increasing the value of/, only a small number of elements in the matrix change from 1 to-1, so that the energy function will increase, and when increasing to a certain value, the neighbor elements in the matrix change from 1 to-1.
Step two: delineation of pixel second order (N) at spatial locations2) The neighboring pixel (i, j), as shown in fig. 1.
The probability of a certain pixel being noise is calculated based on the energy curve as follows:
Figure FDA0002353314970000015
wherein, the lambda is an adjusting parameter,
Figure FDA0002353314970000016
the objective function of the ENFCM may be defined as follows:
Figure FDA0002353314970000021
wherein:
Figure FDA0002353314970000022
the function of space restriction is realized,
Figure FDA0002353314970000023
(djkrepresents a sample point xkAnd xjEuclidean distance) is a penalty factor representing the sample point xkSurrounding sample points and xkSpatial correlation of the classification results, ωikThe larger the value of (d), the stronger the spatial correlation.
Step three: a method of minimizing according to FCM, order
Figure FDA0002353314970000024
The following can be obtained:
Figure FDA0002353314970000025
Figure FDA0002353314970000026
from the above formula, it can be seen that in the introduction of the energy curve ElAnd ωikThereafter, the cluster centers are readjusted, the sample membership is indirectly changed by the change of the cluster centers in the iterative process, and ElAnd ωikAll taking into account the sample distribution around the sample point, El> 0, and ωikOmega when the sample is evenly distributedikIs a constant value.
As can be seen from the equations (5) and (6), after the energy curve and the parameter term are directly introduced into the objective function, only the cluster center is constrained. In order to improve the spatial correlation of the membership degree, the membership degree shown in formula (5) is redefined as follows by a weighting method:
Figure FDA0002353314970000027
wherein
Figure FDA0002353314970000028
Represents the mean of the degrees of membership.
CN202010000820.6A 2020-01-02 2020-01-02 FCM rotary kiln flame image segmentation method based on energy noise detection Active CN111179286B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010000820.6A CN111179286B (en) 2020-01-02 2020-01-02 FCM rotary kiln flame image segmentation method based on energy noise detection

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010000820.6A CN111179286B (en) 2020-01-02 2020-01-02 FCM rotary kiln flame image segmentation method based on energy noise detection

Publications (2)

Publication Number Publication Date
CN111179286A true CN111179286A (en) 2020-05-19
CN111179286B CN111179286B (en) 2022-09-30

Family

ID=70649236

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010000820.6A Active CN111179286B (en) 2020-01-02 2020-01-02 FCM rotary kiln flame image segmentation method based on energy noise detection

Country Status (1)

Country Link
CN (1) CN111179286B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105678766A (en) * 2016-01-06 2016-06-15 福州大学 Fuzzy c-means image segmentation method based on local neighborhood and global information
CN109509196A (en) * 2018-12-24 2019-03-22 广东工业大学 A kind of lingual diagnosis image partition method of the fuzzy clustering based on improved ant group algorithm

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105678766A (en) * 2016-01-06 2016-06-15 福州大学 Fuzzy c-means image segmentation method based on local neighborhood and global information
CN109509196A (en) * 2018-12-24 2019-03-22 广东工业大学 A kind of lingual diagnosis image partition method of the fuzzy clustering based on improved ant group algorithm

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LIANG, F.ET.AL: "Recognition algorithm based on improved FCM and rough sets for meibomian gland morphology", 《APPLIED SCIENCES》 *
R. M. PRAKASH,ET.AL: "Segmentation of thermal infrared breast images using K-means", 《IN PROC. INT. CONF. INNOV. INF., EMBEDDED COMMUN. SYST. (ICIIECS)》 *
王秀秀: "基于自适应模糊C一均值与机会约束支持向量机的图像分割算法", 《中国优秀硕士学位论文全文数据库》 *

Also Published As

Publication number Publication date
CN111179286B (en) 2022-09-30

Similar Documents

Publication Publication Date Title
CN110211126B (en) Image segmentation method based on intuitive fuzzy C-means clustering
CN111340824B (en) Image feature segmentation method based on data mining
CN110322445B (en) Semantic segmentation method based on maximum prediction and inter-label correlation loss function
CN113870235A (en) Method for detecting defects of circular stamping part based on quantum firework arc edge extraction
CN101923637B (en) A kind of mobile terminal and method for detecting human face thereof and device
CN109801305B (en) SAR image change detection method based on deep capsule network
CN107590427B (en) Method for detecting abnormal events of surveillance video based on space-time interest point noise reduction
CN111242026B (en) Remote sensing image target detection method based on spatial hierarchy perception module and metric learning
CN111310768B (en) Saliency target detection method based on robustness background prior and global information
CN114139631B (en) Multi-target training object-oriented selectable gray box countermeasure sample generation method
CN116071339A (en) Product defect identification method based on improved whale algorithm optimization SVM
CN117522864B (en) European pine plate surface flaw detection method based on machine vision
CN117314940B (en) Laser cutting part contour rapid segmentation method based on artificial intelligence
CN117237902B (en) Robot character recognition system based on deep learning
CN111179286B (en) FCM rotary kiln flame image segmentation method based on energy noise detection
CN109829511B (en) Texture classification-based method for detecting cloud layer area in downward-looking infrared image
CN105512682B (en) A kind of security level identification recognition methods based on Krawtchouk square and KNN-SMO classifier
CN113177563B (en) Post-chip anomaly detection method integrating CMA-ES algorithm and sequential extreme learning machine
CN111723808B (en) Infrared weak target extraction method based on maximum value background prediction
CN114862786A (en) Retinex image enhancement and Ostu threshold segmentation based isolated zone detection method and system
CN110517326B (en) Colorimetric sensor array optimization method based on weight dragonfly algorithm
CN113838066A (en) Color image segmentation method based on improved fuzzy c-means clustering algorithm
CN110910392B (en) Image segmentation method based on fuzzy clustering of multiple spatial information weighting
JP2003216953A (en) Method and device for extracting object from image
CN113920391A (en) Target counting method based on generated scale self-adaptive true value graph

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