CN102521802A - Mathematical morphology and LoG operator combined edge detection algorithm - Google Patents

Mathematical morphology and LoG operator combined edge detection algorithm Download PDF

Info

Publication number
CN102521802A
CN102521802A CN2011103868671A CN201110386867A CN102521802A CN 102521802 A CN102521802 A CN 102521802A CN 2011103868671 A CN2011103868671 A CN 2011103868671A CN 201110386867 A CN201110386867 A CN 201110386867A CN 102521802 A CN102521802 A CN 102521802A
Authority
CN
China
Prior art keywords
image
mathematical morphology
edge detection
log operator
log
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.)
Pending
Application number
CN2011103868671A
Other languages
Chinese (zh)
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.)
AUTOMATION ENGINEERING R&M CENTER GUANGDONG ACADEMY OF SCIENCES
Original Assignee
AUTOMATION ENGINEERING R&M CENTER GUANGDONG ACADEMY OF SCIENCES
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 AUTOMATION ENGINEERING R&M CENTER GUANGDONG ACADEMY OF SCIENCES filed Critical AUTOMATION ENGINEERING R&M CENTER GUANGDONG ACADEMY OF SCIENCES
Priority to CN2011103868671A priority Critical patent/CN102521802A/en
Publication of CN102521802A publication Critical patent/CN102521802A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

The invention discloses a mathematical morphology and LoG operator combined edge detection algorithm. The method comprises the following steps of: smoothing an image through a mathematical morphology method; and using the LoG operator to perform edge detection on the image treated by the mathematical morphology method. Since the mathematical morphology has the advantages of simplifying image data, keeping basic shape characteristics of the image, removing unrelated structures, being easy for realizing hardware and the like, when executing edge detection on the image with noise, by combining the mathematical morphology method with the LoG operator, the noise is removed better without influence on the detection efficiency of the LoG operator. The mathematical morphology and LoG operator combined edge detection algorithm is suitable for the edge detection field.

Description

The edge detection algorithm that a kind of mathematical morphology and LoG operator combine
Technical field
The invention belongs to the rim detection field, relate in particular to the edge detection algorithm that a kind of mathematical morphology and LoG operator combine.
Background technology
The edge is the reflection of gradation of image uncontinuity, and it has comprised the bulk information of image, has reflected the essential characteristic of object.Rim detection plays an important role in application such as computer vision, image segmentation, feature extraction and images match, is the important step of graphical analysis and identification.
The method of rim detection has a lot, can be divided into two types substantially: the spatial domain is detected and transform domain detects.Classical edge detection algorithm is based on the spatial domain mostly and detects, like Laplace operator, Sobel operator, Canny operator and LoG operator etc.Various operators are each targeted and characteristic all, and wherein, LoG operator (Gauss-Laplace operator) method that is proposed by Marr and Hildreth has the high characteristics of rim detection efficient.Its basic thought is: carry out smoothing processing with the Gauss function earlier and suppress noise, detect the edge with the image after the processing smoothly of Laplace operator then.Because the threshold value of LoG operator can be calculated by computing machine fully, calculated amount is also few relatively, thereby makes that the detection speed of LoG operator is faster.But because these operators mainly act on is high-pass filtering, all relates to directivity, generally responsive to noise ratio, therefore is difficult to detect complex edge.And can suppress noise effectively based on the edge detection method of transform domain, but this type algorithm computation amount is bigger, is difficult to requirement of real time in a lot of occasions.
In recent years, in order to improve the real-time of rim detection, improve the performance of rim detection, the application of mathematical morphology is more and more paid attention to.Mathematical morphology is to be based upon on random collection opinion and the integral geometry basis, and basic thought is to survey and extract image as basic tool with a structural element.It comprises that mainly expansion, burn into are opened, closed these four kinds of computings, and makes up and derive other algorithm by these four kinds of computings, thereby realizes the rim detection to detected image, accomplishes the pictures different analysis.
Summary of the invention
The purpose of this invention is to provide to having noise image, not only processing speed is fast, and the good again a kind of mathematical morphology of effect and the edge detection algorithm of LoG operator combination.
The edge detection algorithm that a kind of mathematical morphology of the present invention and LoG operator combine may further comprise the steps:
Adopt Mathematical Morphology Method that image is carried out smoothing processing;
Use the LoG operator, carry out rim detection adopting the image after Mathematical Morphology Method is handled.
Further, said employing Mathematical Morphology Method is carried out smoothing processing to image and is specially image is carried out closure operation.
The invention has the beneficial effects as follows: since mathematical morphology have the simplified image data, keep image the grown form characteristic, remove irrelevant structure, be easy to advantage such as hardware realization; When the image that contains noise is carried out rim detection; Through combining Mathematical Morphology Method with the LoG operator; In the detection efficiency that does not influence the LoG operator, can well remove noise again.
Description of drawings
Be described further below in conjunction with the accompanying drawing specific embodiments of the invention:
Fig. 1 is the inventive method flow chart of steps;
Fig. 2 is the glass fragmentation image that contains salt-pepper noise;
Fig. 3 uses the LoG operator image shown in Figure 1 is carried out the image after the rim detection;
To be the edge detection algorithm that adopts mathematical morphology of the present invention and the combination of LoG operator carry out the image after the rim detection to image shown in Figure 1 to Fig. 4.
Embodiment
LoG operator, just Gauss-Laplace operator.Laplace operator is a kind of second derivative operator, and its expression formula is:
▿ 2 f ( x , y ) = ∂ 2 f ( x , y ) ∂ x 2 + ∂ 2 f ( x , y ) ∂ y 2
When being applied to digital picture, then can with difference approximation be expressed as:
▿ 2 f ( x , y ) = f ( x + 1 , y ) + f ( x - 1 , y ) + f ( x , y + 1 ) + f ( x , y - 1 ) - 4 f ( x , y )
Because Laplace operator is a scalar, so it only needs a template just passable when calculating.Because Laplace operator is a second derivative operator, thereby it is very sensitive to noise, so when carrying out Flame Image Process, generally will carry out smothing filtering earlier and carry out second-order differential again.The most frequently used smooth function is exactly a Gaussian function, because Gauss's smoothing filter to the noise effects of removing Normal Distribution clearly.Two-dimensional Gaussian function and it one, second-order partial differential coefficient is as follows:
h ( x , y ) = 1 2 πσ 2 e - x 2 + y 2 2 σ 2
∂ h ( x , y ) ∂ x = - x 2 πσ 4 e - x 2 + y 2 2 σ 2 , ∂ h ( x , y ) ∂ y = - y 2 πσ 4 e - x 2 + y 2 2 σ 2
∂ 2 h ( x , y ) ∂ x 2 = 1 2 πσ 4 [ x 2 σ 2 - 1 ] e - x 2 + y 2 2 σ 2 , ∂ 2 h ( x , y ) ∂ y 2 = 1 2 πσ 4 [ y 2 σ 2 - 1 ] e - x 2 + y 2 2 σ 2
Wherein σ is the standard variance of Gaussian distribution, and it has determined the width of Gaussian filter, uses this function to the result that image carries out behind the smothing filtering to be:
g ( x , y ) = h ( x , y ) ⊗ f ( x , y )
Wherein
Figure BDA0000113346040000042
is the convolution symbol; Image is used Laplace operator again after level and smooth, the result is:
▿ 2 g ( x , y ) = ▿ 2 ( h ( x , y ) ⊗ f ( x , y ) )
In linear system, the order of differential and convolution can exchange, so have:
▿ 2 ( h ( x , y ) ⊗ f ( x , y ) ) = ▿ 2 h ( x , y ) ⊗ f ( x , y ) = 1 πσ 4 [ x 2 + y 2 2 σ 2 - 1 ] e - x 2 + y 2 2 σ 2 ⊗ f ( x , y )
Operator after wherein level and smooth and differential merges
▿ 2 h ( x , y ) = 1 πσ 4 [ x 2 + y 2 2 σ 2 - 1 ] e - x 2 + y 2 2 σ 2
Be exactly Gauss-Laplace operator (Laplacian of Gaussian, LoG).
LoG can remove the noise in the certain size effectively, just more can suppress interference of noise.Because LoG is isotropic, so can save many calculated amount.And its bearing accuracy is high, and detected edge continuity is better.But LoG has some very important shortcomings, because it is to utilize the zero cross point of the second derivative of image to obtain the edge, so very sensitive to noise.
Mathematical morphology is recent emerging a kind of image processing method, and it is based upon on random collection opinion and the integral geometry basis, and basic thought is to survey and extract image as basic tool with a structural element.It comprises that mainly expansion, burn into are opened, closed these four kinds of computings, and makes up and derive other algorithm by these four kinds of computings, thereby realizes the rim detection to detected image, accomplishes the pictures different analysis.
Expanding and corroding is operation the most basic in the mathematical morphology.If A and B are the set among the Z, wherein A is an original image, and B is a structural element, and Z is an integer space.Then have B to the dilation operation note work
Figure BDA0000113346040000051
of A:
Figure BDA0000113346040000052
Wherein
Figure BDA0000113346040000053
representes empty set;
Figure BDA0000113346040000054
is the reflection collection of set B:
Figure BDA0000113346040000055
can find out according to following formula; B is exactly in fact a set of being made up of all translational movement z to the expansion of A; These translational movements z satisfies: when the reflection collection translation of B after the z, with the common factor of set A be sky.The image pixel more shared than original image is many after overexpansion.
Equally, note A Θ B is the erosion operation of B to A, then:
AΘB = { z | ( B ) z ⊆ A }
B is the set of translational movement z to the corrosion of A, and these translational movements still belong to set A after satisfying set B translation z.Result images after the corrosion shrinks with respect to original image to some extent, and the image after the corrosion is a sub-set of original image.
Be expansion and corrode the compound operation that combines and open with closure.Open operation can play the effect of smoothed image, removes burr outstanding on the profile, blocks narrow mountain valley.Closed procedure also has smoothing effect to image outline, but the result is opposite, and it can remove the aperture in the zone, fills and leads up the profile breach.
The unlatching note of A is done A ο B to B, then has:
Figure BDA0000113346040000057
Just with structural element B image A is corroded earlier, with B Corrosion results is done expansive working then.
Similarly, the closed procedure note of A is done AB to B, then:
Figure BDA0000113346040000058
Opposite with open operation, it is earlier original image to be done dilation operation, again expansion results is corroded operation.
With reference to Fig. 1, the edge detection algorithm that a kind of mathematical morphology of the present invention and LoG operator combine may further comprise the steps:
Adopt Mathematical Morphology Method that image is carried out smoothing processing;
Use the LoG operator, carry out rim detection adopting the image after Mathematical Morphology Method is handled.
As the step of further optimization, when adopting Mathematical Morphology Method to handle, specifically adopt closure operation that image is carried out smoothing processing.
In order to prove the superiority that adopts the more simple LoG operator of edge detection algorithm provided by the present invention to handle, on the Matlab7.0 platform same image that has noise carried out handling relatively now.In the test, the threshold value of the LoG operator that both adopt is the same.
In test, be object with the glass fragmentation image, and former figure has been added certain salt-pepper noise, as shown in Figure 2, to have the image of noise in the simulation actual production processing procedure.
After with the LoG operator image shown in Figure 2 being carried out rim detection, result is as shown in Figure 3.
Adopt technical scheme provided by the present invention, earlier image is carried out closure operation, and then carry out edge extracting with the LoG operator, it is as shown in Figure 4 to handle the back result.
Through relatively finding out, for the image that contains noise, when independent use LoG operator carries out edge extracting, noise and graph outline have been hard to tell, and effect is unsatisfactory.And adopt technical scheme of the present invention, after image was carried out closure operation earlier, the LoG operator then can clearly detect the profile of figure.
Therefore, for noisy image, through adopting mathematical morphology earlier it is handled, in smooth noise, kept the elementary contour characteristic of image effectively, thereby helped the rim detection of LoG operator, it is also better to detect effect.
More than be that preferable enforcement of the present invention is specified; But the invention is not limited to said embodiment; Those of ordinary skill in the art make all equivalent variations or replacement under the prerequisite of spirit of the present invention, also can doing, and distortion that these are equal to or replacement all are included in the application's claim institute restricted portion.

Claims (2)

1. the edge detection algorithm that combines of mathematical morphology and LoG operator is characterized in that: may further comprise the steps:
Adopt Mathematical Morphology Method that image is carried out smoothing processing;
Use the LoG operator, carry out rim detection adopting the image after Mathematical Morphology Method is handled.
2. the edge detection algorithm that a kind of mathematical morphology according to claim 1 and LoG operator combine is characterized in that: said employing Mathematical Morphology Method is carried out smoothing processing to image and is specially image is carried out closure operation.
CN2011103868671A 2011-11-28 2011-11-28 Mathematical morphology and LoG operator combined edge detection algorithm Pending CN102521802A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011103868671A CN102521802A (en) 2011-11-28 2011-11-28 Mathematical morphology and LoG operator combined edge detection algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011103868671A CN102521802A (en) 2011-11-28 2011-11-28 Mathematical morphology and LoG operator combined edge detection algorithm

Publications (1)

Publication Number Publication Date
CN102521802A true CN102521802A (en) 2012-06-27

Family

ID=46292707

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011103868671A Pending CN102521802A (en) 2011-11-28 2011-11-28 Mathematical morphology and LoG operator combined edge detection algorithm

Country Status (1)

Country Link
CN (1) CN102521802A (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103136753A (en) * 2013-02-25 2013-06-05 哈尔滨工业大学 Depth image segmentation method based on mean shift algorithm and mathematical morphology
CN104778710A (en) * 2015-04-24 2015-07-15 大连理工大学 Morphological image edge detecting method based on quantum theory
CN105930847A (en) * 2016-03-31 2016-09-07 中国人民解放军空军航空大学 Combined edge detection-base SAR image linear feature extraction method
CN107169982A (en) * 2017-05-17 2017-09-15 重庆邮电大学 A kind of quantum LoG edge detection methods
CN107633525A (en) * 2017-08-23 2018-01-26 天津理工大学 A kind of noise reduction edge detection method based on FPGA
CN108830873A (en) * 2018-06-29 2018-11-16 京东方科技集团股份有限公司 Depth image object edge extracting method, device, medium and computer equipment
CN109766892A (en) * 2018-12-21 2019-05-17 西安交通大学 A kind of ray detection image tagged information character dividing method based on edge detection
CN111738990A (en) * 2020-06-03 2020-10-02 东北林业大学 LOG algorithm-based damaged fruit temperature field detection method
CN114187267A (en) * 2021-12-13 2022-03-15 沭阳县苏鑫冲压件有限公司 Stamping part defect detection method based on machine vision
CN115346048A (en) * 2022-08-25 2022-11-15 珠江水利委员会珠江水利科学研究院 Remote sensing image semantic segmentation method based on boundary point selection algorithm

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101519981A (en) * 2009-03-19 2009-09-02 重庆大学 Mine locomotive anti-collision early warning system based on monocular vision and early warning method thereof
CN101691994A (en) * 2009-09-30 2010-04-07 浙江大学 Method for automatically positioning and detecting maximum width of tunnel crack
CN101718870A (en) * 2009-11-13 2010-06-02 西安电子科技大学 High-speed weak target flight path detection method of image field
CN101763513A (en) * 2010-02-26 2010-06-30 成都三泰电子实业股份有限公司 Foreground extraction method for removing light effect

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101519981A (en) * 2009-03-19 2009-09-02 重庆大学 Mine locomotive anti-collision early warning system based on monocular vision and early warning method thereof
CN101691994A (en) * 2009-09-30 2010-04-07 浙江大学 Method for automatically positioning and detecting maximum width of tunnel crack
CN101718870A (en) * 2009-11-13 2010-06-02 西安电子科技大学 High-speed weak target flight path detection method of image field
CN101763513A (en) * 2010-02-26 2010-06-30 成都三泰电子实业股份有限公司 Foreground extraction method for removing light effect

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
杨雨迎 等: "结合数学形态学的LoG算法在嵌入式射击模拟训练中的应用研究", 《2008系统仿真技术及其应用学术会议论文集》, vol. 10, 31 December 2008 (2008-12-31), pages 451 - 4 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103136753B (en) * 2013-02-25 2016-02-17 哈尔滨工业大学 A kind of depth image segmentation method based on mean shift algorithm and mathematical morphology
CN103136753A (en) * 2013-02-25 2013-06-05 哈尔滨工业大学 Depth image segmentation method based on mean shift algorithm and mathematical morphology
CN104778710A (en) * 2015-04-24 2015-07-15 大连理工大学 Morphological image edge detecting method based on quantum theory
CN104778710B (en) * 2015-04-24 2018-01-09 大连理工大学 A kind of morphological images edge detection method based on quantum theory
CN105930847A (en) * 2016-03-31 2016-09-07 中国人民解放军空军航空大学 Combined edge detection-base SAR image linear feature extraction method
CN107169982A (en) * 2017-05-17 2017-09-15 重庆邮电大学 A kind of quantum LoG edge detection methods
CN107633525A (en) * 2017-08-23 2018-01-26 天津理工大学 A kind of noise reduction edge detection method based on FPGA
CN107633525B (en) * 2017-08-23 2021-03-02 天津理工大学 FPGA-based noise reduction edge detection method
CN108830873B (en) * 2018-06-29 2022-02-01 京东方科技集团股份有限公司 Depth image object edge extraction method, device, medium and computer equipment
CN108830873A (en) * 2018-06-29 2018-11-16 京东方科技集团股份有限公司 Depth image object edge extracting method, device, medium and computer equipment
US11379988B2 (en) 2018-06-29 2022-07-05 Boe Technology Group Co., Ltd. Method and apparatus for extracting edge of object in depth image and computer readable storage medium
CN109766892A (en) * 2018-12-21 2019-05-17 西安交通大学 A kind of ray detection image tagged information character dividing method based on edge detection
CN111738990A (en) * 2020-06-03 2020-10-02 东北林业大学 LOG algorithm-based damaged fruit temperature field detection method
CN114187267A (en) * 2021-12-13 2022-03-15 沭阳县苏鑫冲压件有限公司 Stamping part defect detection method based on machine vision
CN114187267B (en) * 2021-12-13 2023-07-21 沭阳县苏鑫冲压件有限公司 Stamping part defect detection method based on machine vision
CN115346048A (en) * 2022-08-25 2022-11-15 珠江水利委员会珠江水利科学研究院 Remote sensing image semantic segmentation method based on boundary point selection algorithm

Similar Documents

Publication Publication Date Title
CN102521802A (en) Mathematical morphology and LoG operator combined edge detection algorithm
US11620735B2 (en) Method for restoring video data of pipe based on computer vision
Hwang et al. Gaussian filtering detection based on features of residuals in image forensics
CN104794685A (en) Image denoising realization method and device
CN104331869A (en) Image smoothing method with combination of gradient and curvature
CN102169581A (en) Feature vector-based fast and high-precision robustness matching method
CN102096915B (en) Camera lens cleaning method based on precise image splicing
CN105447828B (en) The one-view image deblurring method of one-dimensional deconvolution is carried out along motion blur path
Gao et al. Based on soft-threshold wavelet de-noising combining with Prewitt operator edge detection algorithm
CN104992403A (en) Hybrid operator image redirection method based on visual similarity measurement
Wei et al. Detection of lane line based on Robert operator
CN110599509B (en) Edge detection method based on eight-direction fractional order differential operator
WO2020114134A1 (en) Visual processing method for identifying emery particles
CN102509265B (en) Digital image denoising method based on gray value difference and local energy
CN106600598B (en) Color image tampering detection method based on local grid matching
Modi et al. Skew correction for vehicle license plates using principal component of Harris corner feature
Zhang et al. Detection of composite forged image
CN100365665C (en) Three-D model characteristic line pick-up method based on sharpening filtering
CN107292859B (en) Chaotic medium polarization image acquisition method based on optical correlator
Liu et al. Spatial‐temporal fusion for flotation froth image denoising based on BLS‐GSM method in curvelet domain
Pizano Extracting line features from images of business forms and tables
Lin et al. Abrasive Segmentation of Multiple Diamond Images Based on Secondary Morphological Reconstruction
Li et al. Research on Motion Blur Image of Infrared Target Deblurring Based on Wavelet Transform
CN109840896A (en) A kind of image de-noising method based on gradient and adaptive curvature feature
Chen et al. Research on contour extraction and edge detection of sugarcane image based on MATLAB

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20120627