CN103247047A - Image edge detection method based on fractional order partial differential - Google Patents
Image edge detection method based on fractional order partial differential Download PDFInfo
- Publication number
- CN103247047A CN103247047A CN2013101588035A CN201310158803A CN103247047A CN 103247047 A CN103247047 A CN 103247047A CN 2013101588035 A CN2013101588035 A CN 2013101588035A CN 201310158803 A CN201310158803 A CN 201310158803A CN 103247047 A CN103247047 A CN 103247047A
- Authority
- CN
- China
- Prior art keywords
- differential
- fractional order
- image edge
- edge detection
- fractional
- 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
Links
Images
Landscapes
- Image Analysis (AREA)
Abstract
The invention discloses an image edge detection method based on fractional order partial differential, which comprises the following steps: the Sobel operator and a 3*3 pixel neighborhood with the gray-scale function of F (x, y) are subjected to convolution; the convolution sum is subjected to central difference to obtain an integral order differential expression; the integral order differential expression is replaced by a v order differential expression; based on definition of fractional differential, first three items of the fractional differential expression of a digital image substitute values in the v order differential expression; and different results can be obtained through regulation parameter, namely the differential order v. The method can satisfactorily extract edge contour information, has an excellent effect on texture detail detection, is superior to conventional edge detection methods, and satisfactorily achieves the purpose of image edge detection.
Description
Technical field
The present invention relates to the Image Edge-Detection field, relate in particular to a kind of method for detecting image edge based on the fractional order partial differential.
Background technology
Image Edge-Detection is one of important topic of Digital Image Processing research field, and existing method for detecting image edge is more, and has relative merits separately.
Carrying out rim detection with the Sobel operator is a kind of of partial differential method rim detection, has not been a kind of new edge detection method.Its basic detection principle is directly to utilize horizontal gradient operator and vertical gradient operator to calculate this pixel value, choosing suitable threshold carries out emulation experiment and just can obtain the result, this method can realize soon on Matlab software, and ready-made program code is arranged.Be with the code typing in this method implementation procedure, even copy to paste and just can finish.But the effect that this method realizes is unsatisfactory.
As a kind of popularization of integer rank differential theory, the fractional order differential theory had applied to the image processing rapidly in recent years, and obtained certain achievement.Possessing on the differential basis, integer rank, and definition and the character of fractional order differential have been grasped, and the fractional order differential image is handled relevant document, important references documents such as " based on the rim detection of fractional order differential ", " based on the figure image intensifying of fractional order differential ", " based on the image denoising new model of fractional order partial differential " particularly, after understanding the fractional order differential image processing method, we combine fractional order differential and derive and make new advances with the Sobel operator rim detection model.
We know that the anti-noise ability of Sobel operator is better, can understand the fractional order differential computing can keep image as much as possible when carrying out Image Edge-Detection grain details by read documents.
But the unsatisfactory problem of detection effect that existing method for detecting image edge exists, particularly detected image edge and grain details exist can be improved aspect.
Summary of the invention
The object of the present invention is to provide a kind of method for detecting image edge, this method is effective, can not only quite good detecting go out edge of image, can also detect a large amount of grain details.
For reaching this purpose, the present invention by the following technical solutions:
A kind of method for detecting image edge based on the fractional order partial differential may further comprise the steps:
A, with Sobel operator and gray scale function be F (x, 3 * 3 neighborhood of pixels y) are done convolution;
B, obtain integer rank differential expressions with convolution with central difference again;
C, integer rank differential is replaced with v rank differential;
D, from the definition of fractional order differential, in the first three items substitution v rank differential expressions with the fractional order differential expression formula of digital picture, and can obtain different results by regulating parameter differential exponent number v.
Preferably, integer rank differential is replaced with v rank differential among the described step C, wherein v is greater than 0 and less than 1.
Preferably, integer rank differential is replaced with v rank differential among the described step C, wherein v equals at 0.7 o'clock, and effect is best.
Beneficial effect of the present invention is: this kind extracts edge contour information well based on the method for detecting image edge of fractional order partial differential, detection effect to grain details is also fine, be better than existing other several edge detection methods, reach the purpose of Image Edge-Detection well.
Description of drawings
Fig. 1 is new fractional order row gradient and the row gradient former of the embodiment of the invention;
Fig. 2 is the corresponding differential mask of the embodiment of the invention;
Fig. 3 is the rim detection rose figure of the different differential exponent numbers of new model of the embodiment of the invention;
Fig. 4 is new model and other edge detection methods experiment contrast palace figure of the embodiment of the invention;
Fig. 5 be the embodiment of the invention to the girl figure lab diagram that superposes;
Fig. 6 is the palace figure rim detection objective evaluation index of the embodiment of the invention.
Embodiment
Further specify technical scheme of the present invention below in conjunction with accompanying drawing and by embodiment.
Be that (x, 3 * 3 neighborhood of pixels y) are done convolution to F, and then obtain integer rank differential expressions with convolution with central difference with Sobel operator and gray scale function at first.Again integer rank differential is replaced (0<v<1) with v rank differential, at last from the definition of fractional order differential, in the first three items substitution v rank differential expressions with the fractional order differential expression formula of digital picture, and can obtain different experimental results by regulating parameter differential exponent number v.This has just obtained the Image Edge-Detection new model based on the fractional order partial differential, new fractional order row gradient and row gradient former such as Fig. 1, corresponding differential mask such as Fig. 2.
At first experimentize by choosing different differential exponent number v, as shown in Figure 3.(a) being original image, (b)~(f) is that this paper model is chosen the edge detection results that different fractional order parameters obtains.When 0<v<1, experiment effect improves gradually with the increase of parameter v; Be not difficult to find out that when v=0.7, experiment effect is best, can not only quite good detecting go out edge of image, can also detect a large amount of grain details, shown in Fig. 3 (d).
As can be seen from Figure 4, Sobel operator, Canny operator, LOG operator edge detection method shortcoming separately are not easy to overcome, and the rim detection effect is relatively poor; Existing [8] are that the definition with the fractional order differential method defines the mask operator and carries out Image Edge-Detection as template, though the fractional order differential method can better detect grain details, but the lifting to some edge is not enough, therefore easy lost part edge.New model can extract edge contour information well, and is also fine to the detection effect of grain details, is better than existing other several edge detection methods, reaches the purpose of Image Edge-Detection well.
Below will continue to do the validity of proving new method, be the accuracy that outline map after example will detect and original graph superpose to verify the new method bearing accuracy with girl figure.From the result of Fig. 5 as can be seen, new method is accurate to the edge detection and location.
Below also will come review new initiative to the superiority of rim detection from objective experimental data.At the objective evaluation of rim detection, objective evaluation standard commonly used has linear connection degree L, fallout ratio N and loss F.The connection degree L of efficient frontier is more high, and the edge evaluation is more high; Fallout ratio N is more little, and the edge evaluation is more high; Loss F is more little, and the edge evaluation is more high.3 indexs of objective evaluation are combined, define new edge evaluation tolerance M
e, M
eWeighted mean by 3 indexs obtains, and is shown below.
M
e=α L+ β (1-N)+γ (1-F) is wherein: alpha+beta+γ=1, α, β, γ three numbers represent that the different evaluation index is to estimating the influence degree of tolerance.
We are the example contrast that experimentizes with palace figure, and we might as well get
Result of calculation as shown in Figure 6.Can draw objectively from Fig. 6,3 evaluation indexes of new model rim detection all are better than all the other SOME METHODS, and loss F is 0, can well detect the edge of image grain details; Compare with additive method, the linearity connection degree L of this method model is higher, and fallout ratio N is lower, on the basis of fractional order differential certain improvement is arranged, and obviously is better than existing integer rank differential algorithm.
The above; only for the preferable embodiment of the present invention, but protection scope of the present invention is not limited thereto, and anyly is familiar with the people of this technology in the disclosed technical scope of the present invention; the variation that can expect easily or replacement all should be encompassed within protection scope of the present invention.Therefore, protection scope of the present invention should be as the criterion with the protection domain of claim.
Claims (3)
1. the method for detecting image edge based on the fractional order partial differential is characterized in that, may further comprise the steps:
A, with Sobel operator and gray scale function be F (x, 3 * 3 neighborhood of pixels y) are done convolution;
B, obtain integer rank differential expressions with convolution with central difference again;
C, integer rank differential is replaced with v rank differential;
D, from the definition of fractional order differential, in the first three items substitution v rank differential expressions with the fractional order differential expression formula of digital picture, and can obtain different results by regulating parameter differential exponent number v.
2. a kind of method for detecting image edge based on the fractional order partial differential as claimed in claim 1 is characterized in that, integer rank differential is replaced with v rank differential among the described step C, and wherein v is greater than 0 and less than 1.
3. a kind of method for detecting image edge based on the fractional order partial differential as claimed in claim 1 is characterized in that, integer rank differential is replaced with v rank differential among the described step C, and wherein v equals 0.7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013101588035A CN103247047A (en) | 2013-04-23 | 2013-04-23 | Image edge detection method based on fractional order partial differential |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013101588035A CN103247047A (en) | 2013-04-23 | 2013-04-23 | Image edge detection method based on fractional order partial differential |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103247047A true CN103247047A (en) | 2013-08-14 |
Family
ID=48926552
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2013101588035A Pending CN103247047A (en) | 2013-04-23 | 2013-04-23 | Image edge detection method based on fractional order partial differential |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103247047A (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103927725A (en) * | 2014-05-07 | 2014-07-16 | 哈尔滨工业大学 | Movie nuclear magnetic resonance image sequence motion field estimation method based on fractional order differential |
CN104680548A (en) * | 2015-03-20 | 2015-06-03 | 福州大学 | Omnidirectional M-type kinetocardiogram detection method |
CN107223266A (en) * | 2017-04-27 | 2017-09-29 | 香港应用科技研究院有限公司 | Kernel approximation for the fractional order differential operator of rim detection |
CN107424121A (en) * | 2017-06-30 | 2017-12-01 | 中原智慧城市设计研究院有限公司 | A kind of blurred picture ultra-resolution ratio reconstructing method based on fractional order differential |
WO2018196018A1 (en) * | 2017-04-27 | 2018-11-01 | Hong Kong Applied Science and Technology Research Institute Company Limited | Kernel approximation on fractional differential operator for edge detection |
CN110599509A (en) * | 2019-08-02 | 2019-12-20 | 西安理工大学 | Edge detection method based on eight-direction fractional order differential operator |
CN111402204A (en) * | 2020-02-26 | 2020-07-10 | 哈尔滨工业大学 | Chip appearance defect detection method based on multi-order fractional order wavelet packet transformation |
CN111724318A (en) * | 2020-06-15 | 2020-09-29 | 石家庄铁道大学 | Image denoising method based on mixed high-order partial differential equation model |
CN112179813A (en) * | 2020-08-26 | 2021-01-05 | 清华大学 | Liquid contact angle on-line measurement method based on experimental image |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1917576A (en) * | 2006-08-30 | 2007-02-21 | 蒲亦非 | Fractional order differential filter for digital image |
US20080101716A1 (en) * | 2006-10-27 | 2008-05-01 | Quanta Computer Inc. | Image sharpening apparatus and method thereof |
CN101800847A (en) * | 2010-04-14 | 2010-08-11 | 蒲亦非 | Fractional order differential filter with 1-2 orders of digital image defined based on Riemann-Liouville |
-
2013
- 2013-04-23 CN CN2013101588035A patent/CN103247047A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1917576A (en) * | 2006-08-30 | 2007-02-21 | 蒲亦非 | Fractional order differential filter for digital image |
US20080101716A1 (en) * | 2006-10-27 | 2008-05-01 | Quanta Computer Inc. | Image sharpening apparatus and method thereof |
CN101800847A (en) * | 2010-04-14 | 2010-08-11 | 蒲亦非 | Fractional order differential filter with 1-2 orders of digital image defined based on Riemann-Liouville |
Non-Patent Citations (1)
Title |
---|
蒋伟 等: "基于分数阶微分和Sobel算子的边缘检测新模型", 《计算机工程与应用》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103927725A (en) * | 2014-05-07 | 2014-07-16 | 哈尔滨工业大学 | Movie nuclear magnetic resonance image sequence motion field estimation method based on fractional order differential |
CN103927725B (en) * | 2014-05-07 | 2017-04-26 | 哈尔滨工业大学 | Movie nuclear magnetic resonance image sequence motion field estimation method based on fractional order differential |
CN104680548A (en) * | 2015-03-20 | 2015-06-03 | 福州大学 | Omnidirectional M-type kinetocardiogram detection method |
US10210616B2 (en) | 2017-04-27 | 2019-02-19 | Hong Kong Applied Science And Technology Research Institute Co., Ltd. | Kernal approximation on fractional differential operator for edge detection |
WO2018196018A1 (en) * | 2017-04-27 | 2018-11-01 | Hong Kong Applied Science and Technology Research Institute Company Limited | Kernel approximation on fractional differential operator for edge detection |
CN107223266A (en) * | 2017-04-27 | 2017-09-29 | 香港应用科技研究院有限公司 | Kernel approximation for the fractional order differential operator of rim detection |
CN107223266B (en) * | 2017-04-27 | 2020-07-07 | 香港应用科技研究院有限公司 | Kernel approximation of fractional order differential operators for edge detection |
CN107424121A (en) * | 2017-06-30 | 2017-12-01 | 中原智慧城市设计研究院有限公司 | A kind of blurred picture ultra-resolution ratio reconstructing method based on fractional order differential |
CN110599509A (en) * | 2019-08-02 | 2019-12-20 | 西安理工大学 | Edge detection method based on eight-direction fractional order differential operator |
CN110599509B (en) * | 2019-08-02 | 2021-10-08 | 西安理工大学 | Edge detection method based on eight-direction fractional order differential operator |
CN111402204A (en) * | 2020-02-26 | 2020-07-10 | 哈尔滨工业大学 | Chip appearance defect detection method based on multi-order fractional order wavelet packet transformation |
CN111402204B (en) * | 2020-02-26 | 2021-07-06 | 哈尔滨工业大学 | Chip appearance defect detection method based on multi-order fractional order wavelet packet transformation |
CN111724318A (en) * | 2020-06-15 | 2020-09-29 | 石家庄铁道大学 | Image denoising method based on mixed high-order partial differential equation model |
CN111724318B (en) * | 2020-06-15 | 2022-04-08 | 石家庄铁道大学 | Image denoising method based on mixed high-order partial differential equation model |
CN112179813A (en) * | 2020-08-26 | 2021-01-05 | 清华大学 | Liquid contact angle on-line measurement method based on experimental image |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103247047A (en) | Image edge detection method based on fractional order partial differential | |
CN110428433B (en) | Canny edge detection algorithm based on local threshold | |
CN106709928B (en) | fast two-dimensional maximum inter-class variance threshold method for noisy images | |
CN103942797B (en) | Scene image text detection method and system based on histogram and super-pixels | |
CN109035274B (en) | Document image binarization method based on background estimation and U-shaped convolution neural network | |
CN101710425B (en) | Self-adaptive pre-segmentation method based on gray scale and gradient of image and gray scale statistic histogram | |
CN103455991A (en) | Multi-focus image fusion method | |
CN109671035B (en) | Infrared image enhancement method based on histogram | |
CN110766708B (en) | Image comparison method based on contour similarity | |
CN103116879A (en) | Neighborhood windowing based non-local mean value CT (Computed Tomography) imaging de-noising method | |
CN114418904B (en) | Infrared image enhancement method based on improved histogram equalization and enhanced high-pass filtering | |
CN103337073B (en) | A kind of two dimensional image threshold segmentation method based on three-dimensional entropy | |
CN111429372A (en) | Method for enhancing edge detection effect of low-contrast image | |
CN111445386A (en) | Image correction method based on four-point detection of text content | |
CN104021528A (en) | Dictionary learning algorithm based on sparse model analysis | |
CN107729863B (en) | Human finger vein recognition method | |
CN117893550A (en) | Moving object segmentation method under complex background based on scene simulation | |
CN106023184A (en) | Depth significance detection method based on anisotropy center-surround difference | |
González-Hidalgo et al. | Impulsive noise removal with an adaptive weighted arithmetic mean operator for any noise density | |
Guanglin et al. | A simple and fast method of recognizing license plate number | |
CN106446904A (en) | Image recognition method based on global binarization | |
Sanmiguel et al. | Pixel-based colour contrast for abandoned and stolen object discrimination in video surveillance | |
CN116563140A (en) | Image enhancement method based on target protection | |
CN102855612A (en) | Adaptive enhancement algorithm of CCD (charged coupled device) image based on linear gray array | |
CN108932520A (en) | In conjunction with the SAR image water body probability drafting method of prior probably estimation |
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: 20130814 |