CN102663395B - A straight line detection method based on self-adaptation multi-scale fast discrete Beamlet transform - Google Patents

A straight line detection method based on self-adaptation multi-scale fast discrete Beamlet transform Download PDF

Info

Publication number
CN102663395B
CN102663395B CN 201210054953 CN201210054953A CN102663395B CN 102663395 B CN102663395 B CN 102663395B CN 201210054953 CN201210054953 CN 201210054953 CN 201210054953 A CN201210054953 A CN 201210054953A CN 102663395 B CN102663395 B CN 102663395B
Authority
CN
China
Prior art keywords
image
edge
block
straight line
self
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.)
Expired - Fee Related
Application number
CN 201210054953
Other languages
Chinese (zh)
Other versions
CN102663395A (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.)
Xi'an Anmeng Intelligent Technology Co ltd
Original Assignee
Northwestern Polytechnical 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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN 201210054953 priority Critical patent/CN102663395B/en
Publication of CN102663395A publication Critical patent/CN102663395A/en
Application granted granted Critical
Publication of CN102663395B publication Critical patent/CN102663395B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

The present invention belongs to the technical field of image processing, and particularly relates to a straight line detection method based on self-adaptation multi-scale fast discrete Beamlet transform. Technical features of the method rest with following details: an edge detection image E is obtained by carrying out edge detection on an image I using a Canny operator, gradient direction angles of edge pixel points on the edge detection image E are calculated and the edge contour is divided into various linearity regions based on a gradient direction angle of two edge points, and fast discrete Beamlet transform is carried out on a current self-adaptation partition subblock through linearity region screening to detect straight lines. The method reduces calculation amount by a certain degree because the scale of an AMFDBT algorithm initial blockette is in a self-adaptation set by the region where a target point is] in the image. At the same time, the initial Beamlet transform can easily cause straight lines to be truncated while the self-adaptive change in partition sub-block makes result of the detected straight line more complete with no truncation of straight lines.

Description

Line detection method based on self-adapting multi-dimension fast discrete Beamlet conversion
Technical field
The invention belongs to image and process, relate in particular to a kind of line detection method based on self-adapting multi-dimension fast discrete Beamlet conversion.
Background technology
The extraction of line feature is a study hotspot in image processing and area of pattern recognition always, although traditional line detection algorithm can both detect straight line quickly, length information and the starting point of line segment and accurate location that endpoint information can't realize line segment because line segment is not provided all.
The Beamlet conversion is one of effective tool of multi-scale geometric analysis, be used for recovering straight line from noisy image, curve and boxed area, it is easily realized the accurate location of line segment and the multiple dimensioned of line segment is similar to, on noisy linear feature detection and fitting a straight line, obvious advantage and potentiality is arranged.But the calculated amount of traditional B eamlet conversion is very large, and many improved fast algorithms are all wasted a lot of calculated amount because will carry out a large amount of repeated works.And advance square because the blockette size of traditional Beamlet conversion is two, if it is relatively loaded down with trivial details to be applied to general rectangular image.Simultaneously, because algorithm is to advance on the basis of recurrence to realize two, the situation that easily causes straight line to be truncated.
In sum, there are the deficiencies such as calculated amount is large, testing result is imperfect in traditional Beamlet conversion.
Summary of the invention
The technical matters that solves
For fear of the deficiencies in the prior art part, the present invention proposes a kind of line detection method based on self-adapting multi-dimension fast discrete Beamlet conversion.
Technical scheme
A kind of line detection method based on self-adapting multi-dimension fast discrete Beamlet conversion is characterized in that step is as follows:
Step 1: utilize the Canny operator to carry out rim detection to image I, obtain edge-detected image E;
Step 2: the deflection of edge calculation detected image E coboundary pixel gradient:
Figure BDA0000140613440000021
When the gradient direction angle of two marginal points differs by more than predetermined threshold value threshold, in this this marginal point place's disconnection, edge contour is divided into a plurality of ranges of linearity and gives mark; In formula, G xAnd G yThe horizontal and vertical gradient component that represents respectively marginal point;
Step 3: range of linearity inward flange is counted is less than the district's deletion of 20 line, obtain through the image E ' after line district's screening, and with the marginal point after rim detection as impact point;
Step 4: take edge image E ' as initial subregion sub-block, seek the Minimum Area of impact point existence as the adaptive partition sub-block; If this sub-block is more than or equal to defined threshold t ' time execution in step 5; Otherwise execution in step 6;
Step 5: current adaptive partition sub-block is carried out fast discrete Beamlet conversion, detection of straight lines.If find the beamlet straight line, record end points and the slope information of straight line, then wipe this straight line under current sub-block; Otherwise, with the current blockette quartern, and as the initial size of next yardstick subregion sub-block, repeated execution of steps 4;
Step 6: if the quantity of beamlet straight line greater than 0, according to the slope information of all straight lines, find straight slope poor 1 the degree with all interior straight lines as parallel lines.
Beneficial effect
A kind of line detection method based on self-adapting multi-dimension fast discrete Beamlet conversion that the present invention proposes, because the yardstick of AMFDBT algorithm primary partition piece is to be set by the region adaptivity that impact point in figure exists, in the situation that the most extreme impact point just can spread all over full figure, this makes calculated amount obtain certain reducing.Simultaneously, the situation that initial Beamlet conversion easily causes straight line to be truncated, and the adaptive change of subregion sub-block makes the result of final detection of straight lines more complete, phenomenon is blocked in difficult generation.
Description of drawings
Fig. 1: based on the parallel lines overhaul flow chart of AMFDBT
Embodiment
Now in conjunction with the embodiments, the invention will be further described for accompanying drawing:
The present invention has defined the self-adaptation tile and has returned subregion (Adaptive Tile Recursive Dyadic Partitioning is called for short ATRDP) algorithm on the basis that the recursive partitioning algorithm is advanced in traditional discrete Beamlet conversion two.Here use S I, j, lRepresent each subregion sub-block, wherein (i, j) represents x coordinate and the y coordinate of this sub-block upper left angle point in original image, and l represents the length of side of sub-block, and k represents k subregion sub-block, and J represents to carry out for the J time the self-adaptation tile and returns subregion.The ATRDP algorithm is defined as follows:
1. ATP JAn ATRDP.
If 2. ATP J = { S i 0 , j 0 , l 0 , . . . . S i k , j k , l k , . . . , S i n , j n , l n } An ATRDP, and If satisfy the adaptive partition condition (refer to find beamlet b line segment in this blockette and eliminate from figure the contained impact point of b line segment can change blockette in the distributed areas of impact point),
Figure BDA0000140613440000033
Be adjusted into by the region adaptivity that exists by its intrapartition destination punctuate
Figure BDA0000140613440000034
ATP J + 1 = { S i 0 , j 0 , l 0 , . . . , S i k , j k , l k , S ( adap ) i k , j k , l k , . . . , S i n , j n , l n }
Also ATRDP.
3. If do not satisfy the adaptive partition condition, and
Figure BDA0000140613440000037
Can be subdivided into four equal subregion sub-blocks
Figure BDA0000140613440000038
Its yardstick is readjusted in the zone that exists by impact point in these four sub-blockettes respectively, obtains
Figure BDA0000140613440000039
Figure BDA00001406134400000310
ATP J + 1 = { S i 0 , j 0 , l 0 , . . . , S i k , j k , l k , . . . , S ( adap ) i k , j k , l k / 2 , S ( adap ) i k + l k / 2 , j k , l k / 2 , S ( adap ) i k , j k + l k / 2 , l k / 2 , S ( adap ) i k + l k / 2 , j k + l k / 2 , l k / 2 . . . , S i n , j n , l n } Also ATRDP, wherein
Figure BDA00001406134400000312
Be
Figure BDA00001406134400000313
Father's piece.
With ATRDP partition method and fast discrete Beamlet conversion (Fast Discrete Beamlet Transform is called for short FDBT) combination, definition self-adapting multi-dimension fast discrete Beamlet is transformed to:
Figure BDA00001406134400000315
Wherein
Figure BDA0000140613440000041
Be fast discrete Beamlet conversion, be defined as follows:
Figure BDA0000140613440000042
G ( f ( x , y ) ) = 1 f ( x , y ) &GreaterEqual; t 0 f ( x , y ) < t - - - ( 3 )
In formula (1), beamlet b is a line segment in subregion sub-block S, and S is ATRDP subregion sub-block.The cut set of AMFDBT has comprised all energy functions
Figure BDA0000140613440000044
(wherein || b|| is the element number on beamlet b) over assign thresholds t *Beamlet b, will no longer consider in the fine dimension sub-block that next stage decomposes at beamlet b that sub-block under the last layer large scale has existed, the decomposition of the beamlet b under large scale can be expressed as b=U jb jSo the cut set of AMFDBT is:
B t * = { b | E ( b ) > t * , | b &Element; S , S &Element; ATP J , b &NotElement; D ~ t * } - - - ( 4 )
D ~ t * = { b j | b = U j b j , E ( b ) > t * , b &Element; P } - - - ( 5 )
Wherein P is father's piece of S.
(2) parallel lines based on AMFDBT detect
The present invention is applied to that with AMFDBT image is carried out parallel lines and detects.
At first will carry out series of preprocessing to image before parallel detection, its concrete steps are as follows:
1. utilize the Canny operator to carry out rim detection to image I, obtain edge-detected image E.
2. range of linearity screening is rejected curve or than short lines.Due in image I, exist a large amount of complex backgrounds toward contact except target, after rim detection, their main manifestations go out a little irregular or crooked edge features, and these edge features can cause larger impact in the straight-line detection process.Therefore before straight-line detection, first edge image E carries out the edge contour screening, only stays the more range of linearity of marginal point in the same way, removes some crooked irregular or too short edge features.
The present invention is mainly the Gradient Phase information of utilizing edge pixel point, realizes the screening of the range of linearity.At first the deflection of edge calculation pixel gradient:
Figure BDA0000140613440000047
In formula, G xAnd G yThe horizontal and vertical gradient component that represents respectively marginal point.When the gradient direction angle of two marginal points differed by more than predetermined threshold value threshold, expression edge contour direction transformation was larger, in this marginal point place's disconnection.Edge contour can be divided into a lot of ranges of linearity like this, mark these ranges of linearity, statistics range of linearity inward flange is counted, and detect with further minimizing and disturb in deletion marginal point less line district.We will be designated as E ' through the image after line district's screening.
Utilize the AMFDBT parallel detection on the basis of edge image E ', specific algorithm is described below:
1. traversing operation subgraph (initial operation subgraph is edge image E '), the zone of finding impact point in figure (being that image is through the marginal point after rim detection) to exist.If this area size is carried out next step greater than defined threshold t ' (t ' can adjust according to the image size is generally picture traverse/100); Otherwise execution in step 4..
2. the impact point zone of 1. finding with step is carried out fast discrete Beamlet conversion as the subregion sub-block to this blockette, finds out all the beamlet straight lines in current blockette, calculates its slope according to the rectilinear end dot information.
If 3. find the beamlet straight line in step in 2., extract straight line information and also wipe this straight line; Otherwise, with the current blockette quartern, and as the initial size of next yardstick subregion sub-block.Execution in step 1..
If 4. detect the beamlet straight line, according to the slope information of straight line, find out all parallel lines.
Embodiment of the present invention reference flow sheet (1), as follows based on the parallel lines detecting step of self-adapting multi-dimension Beamlet conversion:
The present invention is applied to that with AMFDBT image is carried out parallel lines and detects, and its concrete steps are as follows:
(1) image pre-service
1. utilize the Canny operator to carry out rim detection to image I, obtain edge-detected image E.
2. range of linearity screening is rejected curve or than short lines.Due in image I, exist a large amount of complex backgrounds toward contact except target, after rim detection, their main manifestations go out a little irregular or crooked edge features, and these edge features can cause larger impact in the straight-line detection process.Therefore before straight-line detection, first edge image E carries out the edge contour screening, only stays the more range of linearity of marginal point in the same way, removes some crooked irregular or too short edge features.
The present invention is mainly the Gradient Phase information of utilizing edge pixel point, realizes the screening of the range of linearity.At first the deflection of edge calculation pixel gradient:
Figure BDA0000140613440000061
In formula, G xAnd G yThe horizontal and vertical gradient component that represents respectively marginal point.When the gradient direction angle of two marginal points differed by more than predetermined threshold value threshold (threshold=20 herein), expression edge contour direction transformation was larger, in this marginal point place's disconnection.Edge contour can be divided into a lot of ranges of linearity like this, mark these ranges of linearity, statistics range of linearity inward flange is counted, and detect with further minimizing and disturb in deletion marginal point less line district.We will be designated as E ' through the image after line district's screening.
(2) parallel lines based on AMFDBT detect
1. traversing operation subgraph (initial operation subgraph is edge image E '), the zone of finding impact point in figure (being that image is through the marginal point after rim detection) to exist.If this area size is carried out next step greater than defined threshold t ' (t ' can adjust according to the image size is generally picture traverse/100); Otherwise execution in step 4..
2. the impact point zone of 1. finding with step is carried out fast discrete Beamlet conversion as the subregion sub-block to this blockette, finds out all the beamlet straight lines in current blockette, calculates its slope according to the rectilinear end dot information.
If 3. find the beamlet straight line in step in 2., extract straight line information (end points and slope) and also wipe this straight line; Otherwise, with the current blockette quartern, and as the initial size of next yardstick subregion sub-block.Execution in step 1..
If 4. detect the beamlet straight line, according to the slope information of straight line, find out all parallel lines.

Claims (1)

1. line detection method based on self-adapting multi-dimension fast discrete Beamlet conversion is characterized in that step is as follows:
Step 1: utilize the Canny operator to carry out rim detection to image I, obtain edge-detected image E;
Step 2: the deflection of edge calculation detected image E up contour point gradient:
Figure FDA00003383106600011
When the gradient direction angle of two marginal points differs by more than predetermined threshold value threshold, expression edge contour direction transformation is larger, judge that these two marginal points do not belong to the same range of linearity, are divided into several ranges of linearity and separate marking with this with all marginal points; In formula, G xAnd G yThe horizontal and vertical gradient component that represents respectively marginal point;
Step 3: all are contained the edge count and be less than the deletion of 20 the range of linearity, obtain through the image E ' after the screening of the range of linearity, and with edge image E ' as subregion sub-block at first, with all marginal points in E ' as impact point;
Step 4: seek Minimum Area that impact point exists as the adaptive partition sub-block under current subregion sub-block, if the size of this adaptive partition sub-block during more than or equal to predetermined threshold value t', execution in step 5; Otherwise execution in step 6;
Step 5: current adaptive partition sub-block is carried out fast discrete Beamlet conversion, detection of straight lines; If find the beamlet straight line, record end points and the slope information of straight line, then wipe this straight line under current subregion sub-block; Otherwise, with the current subregion sub-block quartern, and as the initial size of next yardstick subregion sub-block, turn back to step 4;
Step 6: if the quantity of beamlet straight line greater than 0, according to the slope information of all straight lines, find straight slope poor 1 the degree with all interior straight lines as parallel lines.
CN 201210054953 2012-03-05 2012-03-05 A straight line detection method based on self-adaptation multi-scale fast discrete Beamlet transform Expired - Fee Related CN102663395B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201210054953 CN102663395B (en) 2012-03-05 2012-03-05 A straight line detection method based on self-adaptation multi-scale fast discrete Beamlet transform

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201210054953 CN102663395B (en) 2012-03-05 2012-03-05 A straight line detection method based on self-adaptation multi-scale fast discrete Beamlet transform

Publications (2)

Publication Number Publication Date
CN102663395A CN102663395A (en) 2012-09-12
CN102663395B true CN102663395B (en) 2013-11-06

Family

ID=46772879

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201210054953 Expired - Fee Related CN102663395B (en) 2012-03-05 2012-03-05 A straight line detection method based on self-adaptation multi-scale fast discrete Beamlet transform

Country Status (1)

Country Link
CN (1) CN102663395B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103914830B (en) 2014-02-22 2017-02-01 小米科技有限责任公司 Straight line detection method and device
CN110298845A (en) * 2019-06-17 2019-10-01 中国计量大学 It transmits electricity under a kind of complex background based on image procossing line detecting method
CN111047615B (en) * 2019-12-09 2024-02-02 Oppo广东移动通信有限公司 Image-based straight line detection method and device and electronic equipment
CN111445510A (en) * 2020-03-24 2020-07-24 杭州东信北邮信息技术有限公司 Method for detecting straight line in image

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ATE525705T1 (en) * 2007-04-13 2011-10-15 Pasteur Institut CHARACTERISTIC BEAMLET TRANSFORMATION APPARATUS AND RELATED METHODOLOGY FOR DETECTING CURVILINE OBJECTS OF AN IMAGE
CN101159009B (en) * 2007-11-09 2010-04-21 西北工业大学 Method for detecting bridge from remote sense image
CN101915764B (en) * 2010-08-10 2012-02-01 武汉武大卓越科技有限责任公司 Road surface crack detection method based on dynamic programming
CN102156882A (en) * 2011-04-14 2011-08-17 西北工业大学 Method for detecting airport target based on high-resolution remote sensing image

Also Published As

Publication number Publication date
CN102663395A (en) 2012-09-12

Similar Documents

Publication Publication Date Title
EP3082066B1 (en) Road surface gradient detection device
CN102663395B (en) A straight line detection method based on self-adaptation multi-scale fast discrete Beamlet transform
CN107833206B (en) Accurate power line extraction method under complex background
CN102289810B (en) Quick rectangle detection method of images high resolution and high order of magnitude
CN104331689B (en) The recognition methods of a kind of cooperation mark and how intelligent individual identity and pose
CN103258201A (en) Form line extraction method integrating global information and local information
CN104408445A (en) Automatic real-time human body detecting method
CN105678760A (en) Method for recognizing insulator image on the basis of Canny edge detection algorithm
CN107832674B (en) Lane line detection method
CN109001757A (en) A kind of parking space intelligent detection method based on 2D laser radar
CN103870834A (en) Method for searching for sliding window based on layered segmentation
CN105023265A (en) Checkerboard angular point automatic detection method under fish-eye lens
CN104008542A (en) Fast angle point matching method for specific plane figure
CN106504247B (en) A kind of quick-fried automatic testing method of the insulator chain based on visible images and device
CN104182728A (en) Vehicle logo automatic location and recognition method based on pattern recognition
CN106875430B (en) Single moving target tracking method and device based on fixed form under dynamic background
CN102054275B (en) Real-time detection method for sea sky line of gray level image
CN104036514A (en) Circle detection method based on histogram peak value search
CN104156722A (en) Airport target detection method based on high-resolution remote sensing image
CN105354575A (en) Image binarization threshold determination method in sea-sky-line detection
CN110084779A (en) A kind of extraction of aircraft thickness covering end surface features point and denoising method based on laser scanning
CN106897989B (en) A kind of fingerprint image dividing method calculated based on line density
CN110490887B (en) 3D vision-based method for quickly identifying and positioning edges of rectangular packages
CN108898080B (en) Ridge line neighborhood evaluation model-based crack connection method
CN111738094B (en) Automatic extraction method and device for remote sensing image road

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20190226

Address after: 710075 Oak Constellation 1 Building 2 Unit 20301, North of Science and Technology Fifth Road, Xi'an High-tech Zone, Shaanxi Province

Patentee after: Xi'an Xu Yao Network Technology Co.,Ltd.

Address before: 710072 No. 127 Youyi West Road, Shaanxi, Xi'an

Patentee before: Northwestern Polytechnical University

TR01 Transfer of patent right

Effective date of registration: 20190425

Address after: 710000 Unit 2, Unit 3, Building 5, No. 2, Hongpu South Road, Yanta District, Xi'an City, Shaanxi Province

Patentee after: Zhao Feng

Address before: 710075 Oak Constellation 1 Building 2 Unit 20301, North of Science and Technology Fifth Road, Xi'an High-tech Zone, Shaanxi Province

Patentee before: Xi'an Xu Yao Network Technology Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20190917

Address after: 710075 Room 1005, Tongda International Building, 41 Tangyan Road, Xi'an High-tech Zone, Shaanxi Province

Patentee after: Xi'an Anmeng Intelligent Technology Co.,Ltd.

Address before: 710000 Unit 2, Unit 3, Building 5, No. 2, Hongpu South Road, Yanta District, Xi'an City, Shaanxi Province

Patentee before: Zhao Feng

TR01 Transfer of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20131106

CF01 Termination of patent right due to non-payment of annual fee