CN111723335B - Target symmetry axis detection method based on concentric circumference filter - Google Patents
Target symmetry axis detection method based on concentric circumference filter Download PDFInfo
- Publication number
- CN111723335B CN111723335B CN202010434592.3A CN202010434592A CN111723335B CN 111723335 B CN111723335 B CN 111723335B CN 202010434592 A CN202010434592 A CN 202010434592A CN 111723335 B CN111723335 B CN 111723335B
- Authority
- CN
- China
- Prior art keywords
- layer
- relative
- filter
- circumference
- point
- 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.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 15
- 238000005070 sampling Methods 0.000 claims abstract description 86
- 238000006243 chemical reaction Methods 0.000 claims abstract description 4
- 239000010410 layer Substances 0.000 claims description 125
- 230000004044 response Effects 0.000 claims description 57
- 230000006870 function Effects 0.000 claims description 10
- 238000000034 method Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 5
- 230000003287 optical effect Effects 0.000 claims description 3
- 230000000717 retained effect Effects 0.000 claims description 3
- 239000002356 single layer Substances 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 2
- 238000001914 filtration Methods 0.000 claims 1
- 238000013507 mapping Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C1/00—Measuring angles
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/60—Rotation of whole images or parts thereof
- G06T3/604—Rotation of whole images or parts thereof using coordinate rotation digital computer [CORDIC] devices
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Remote Sensing (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Image Analysis (AREA)
Abstract
The invention discloses a target symmetry axis detection method based on a concentric circumference filter. And then taking the overlapped part of the concentric circumference and the target edge as a sampling point, carrying out statistics on angle information under a polar coordinate system through coordinate system conversion, and finally carrying out angle judgment to obtain the direction of the symmetry axis.
Description
Technical Field
The invention relates to a target symmetry axis detection method based on a concentric circumference filter, and belongs to the technical field of digital image processing.
Background
In supervised learning, the classification problem can be generally decomposed into two parts, namely, the image extraction features and the classification method. A strong feature can have the ability to well distinguish different attributes of different targets, so the core part of the classification problem is the feature extraction method. However, many features do not directly have rotation invariance, and therefore, data enhancement is required when the sample direction is not sufficient, and a lot of work is added on training samples.
The circular filter is a filter for detecting a target by utilizing Fourier transform of a circular vector, essentially only amplitude information is utilized, and phase information does not play any role, so that information waste is caused.
Aiming at the situation, polar coordinates are introduced into a circular filter, the collinear situation of the symmetrical targets is solved through relative coordinate clustering, and finally a tangent function is calculated to obtain the direction of the symmetrical axis. After the direction of the symmetry axis is found, the target can be rotated to a specific direction no matter what the original direction of the target is, so that certain rotation invariance can be ensured no matter what features are extracted subsequently.
Disclosure of Invention
The purpose of the invention is as follows: aiming at the problems in the prior art, the invention provides a target symmetry axis detection method based on a concentric circumference filter. According to the method, the concentric circumference filter is constructed, collinear pixel points are accumulated, the collinearity is quantized according to the similarity of relative angles, and the direction of a target symmetry axis can be effectively calculated.
The technical scheme is as follows: in order to realize the purpose of the invention, the technical scheme adopted by the invention is as follows: a target symmetry axis detection method based on a concentric circumference filter comprises the following steps:
(1) Constructing a concentric circular filter: the input target image is taken as an original image img, the original image img is filtered through a circumferential filter, corresponding frequency can be extracted through priori knowledge, a frequency response graph can be obtained, a frequency response threshold value is preset, a series of concentric circles with the same circle center and the same radius and gradually increased can be constructed by taking the frequency position larger than the frequency response threshold value as the center, frequency response is extracted from each concentric circle, and the concentric circles smaller than the frequency response threshold value are discarded, so that a concentric circle filter can be obtained;
(2) Angle statistics: firstly, extracting an edge image edge from an original image img, overlapping a concentric circle filter with the edge image edge, wherein the overlapped part is a sampling point, converting the coordinate of the sampling point from a Cartesian coordinate to a polar coordinate, then calculating the relative angle of the sampling point relative to a reference point by taking the alpha-layer circle sampling point as the reference point, measuring the similarity between the relative angles of the sampling points of different layers, and finally performing collinear statistics by taking the similar angle as a collinear angle;
(3) And (3) angle judgment: and converting the angle into a tangent value, and obtaining the angle in the direction of the symmetry axis.
Further, in the step (1), a specific calculation method for constructing the concentric circular filter is as follows:
(1.1) assuming that a target needing to be subjected to symmetry axis detection is an airplane target, the priori knowledge shows that a general airplane is provided with a left wing, a right wing, a front fuselage and a rear fuselage in front of and behind the wings, and 4 components in total, so that a circumference can be constructed on an input optical image img, a circumference vector is taken as a one-dimensional signal to be subjected to Fourier transform, whether the circumference contains the target or not can be judged according to a transform result, and the coordinates of the target on the img can be judged by updating all coordinates on the img to the circle center of the circumference:
where π represents the circumference ratio, sin and cos represent the sine and cosine functions, respectively, and z k Expressing as the kth pixel value on the circumference, length expresses the total number of circumferential pixels, fourier (p, q) expresses the Fourier transform response corresponding to the circumferential vector with the circle center (p, q), and simultaneously expresses that the response value is placed at the position (p, q) of the frequency response graph, and the frequency response graph fourier can be generated by sequentially updating all coordinates on the original graph img as the circle center for response extraction;
(1.2) constructing a concentric circumference filter by taking the frequency response position larger than the frequency response threshold as the center of a circle from the frequency response graph fourier:
circleFilter=[filter 0 (p,q),...,filter num (p,q)]s.t.fourier(p,q)≥FreqThres
wherein, fourier (p, q) is expressed as frequency response on coordinates (p, q), freqThres is frequency response threshold, filter 0 (p, q) as the layer 0 circumference filter with the circle center (p, q), namely the layer 0 circumference, and the same principle is that the filter num (p, q) is the num layer of circumference filter with the circle center (p, q), namely the outermost layer circumference, num is the total number of the circumferences, and circleFilter represents a concentric circumference filter set;
wherein, the single-layer circumference curve equation with the circle center being (p, q) is as follows:
wherein (x, y) represents the coordinates of a pixel point on the circumference, four k (p, q) and radius k (p, q) are each a filter corresponding to a circular filter k Fourier transform response and radius of (p, q), and filter k (p, q) is the k-th layer circular filter with the center of the circle being (p, q), the frequency response of each circle in the concentric circles is greater than the frequency response threshold value FreqThres, and once the frequency response is less than the threshold value, the circle is determined to be beyond the target boundary.
Further, in the step (2), a specific statistical method of angle statistics is as follows:
(2.1) firstly, extracting an edge image edge from img by using a canny operator, overlapping the edge image edge with a circleFilter, wherein the overlapped part is a sampling point, and a sampling point setWherein->Representing a circular filter j The coordinates of the ith sampling point are (p, q), and edge (p, q) represents the pixel value at the coordinates (p, q) on the edge image edge, and similarly, the filter j (p, q) denotes a circular filter j A pixel value at coordinate (p, q);
(2.2) obtaining the polar coordinate representation of the sampling point by using the following coordinate system conversion formula:
wherein,is the sampling point->Relative angle relative to the origin in a polar coordinate system, andis the sampling point->Expressed in terms of the distance from the origin in a polar coordinate system>Then it is the sampling point pick>Abscissa in the Cartesian coordinate system, for the same reason, in>Then it is the sampling point pick>Ordinate in cartesian coordinate system;
(2.3) taking the sampling points on the alpha layer circumference as reference points, namely calculating the relative angles of the rest sampling points relative to the reference points, wherein the following formula gives calculation by taking the gamma layer and the alpha layer as examplesThe specific mode of (1):
wherein,represents the ith sample point of the most gamma layer, and>represents the jth sampling point of the alpha-most layer, andis the relative angle between the two sampling points;
(2.4) by measuring the relative angle similarity between the sampling point of the other layer and the sampling point of the beta layer, two sampling points with the difference smaller than a similarity threshold value can be considered to be collinear, and the circumferences of the two collinear layers are considered to be in the same weightMarking:
wherein,the method is characterized in that the relative angles of the ith layer q point and the 1 st layer p point are respectively determined, the alpha layer j point is used as a starting point, the similar thresholds are used as similar thresholds, the above formula constraint shows that all sampling points of the ith layer are traversed, and if the sampling points of the ith layer and the/H are available>Is less than the similarity threshold, then it is assumed that the two layers are collinear and that the similarity weight->Indicating whether or not the i-th layer contains collinearly>Is sampled, the similar weights of all layers are summed to obtain->Number of collinear layers in direction:
collinear layer number statistics for representing the relative angle of the jth point of the alpha layer and the pth point of the beta layer, however, each sampling point on the beta layer should be only attributed to one of the components of the front fuselage, the rear fuselage, the left wing and the right wing, therefore, after all the alpha layer points are traversed, the most possible statistical result sW is taken under the condition of all different reference circumference starting points p :
Wherein m represents the total number of alpha layer sampling points, sW p Representing the statistical result with the largest number of collinear layers of the p-th sampling point in the beta-th layer under the condition of all different reference circumference starting points, and obtaining the similarWeight = [ sW ] after the statistical results of all the points on the beta-th layer are all calculated 1 ,...,sW n ]N represents the total number of sampling points of the beta layer, namely a statistical result set;
(2.5) reserving the direction in which the number of layers is greater than the threshold value of the collinear number of layers, and regarding sampling points in the direction smaller than the threshold value as interference:
wherein,represents the statistical result that the p th collinear layer number in the beta layer meets the condition, layerThres is a collinear layer number threshold, pick _ num is the total number of sampling points meeting the condition, and/or is/is selected>And showing that the statistical result set of which the number of collinear layers is smaller than the collinear threshold layerThres is screened out.
Further, in step (3), a specific determination method of the angle determination is as follows:
(3.1) step (2) to obtainIs arranged therein>The corresponding concrete similarity weight is->Thus, by similar weight, the relative angle set retained by the beta layer is combined into ≦ based on the weight of the beta layer>Since the target of the airplane is divided into 4 parts such as the front fuselage, the rear fuselage, the left wing, the right wing and the like, the relative angle set can be divided into 4 parts by the KNN method, and the relative angle can be further expressed by 4 angles, namely [ relative 1 ,...,relative 4 ]Substituting the relative angle value into the tangent function, and determining the direction of the symmetry axis by the following formula:
symmetryAngle=relative j s.t.tan(relative i )-tan(relative j )=tan(relative j )-tan(relative k )
wherein tan represents a tangent function, relative i 、relative k Respectively showing the relative angles of the left wing and the right wing j I.e. the relative angle of the fuselage, and symmertyangle is the direction of the symmetry axis, since the left wing and the right wing are symmetric to the fuselage, the difference between the relative angles should be consistent, and in the actual operation, to avoid the calculation of the negative angle, all the negative angles are added with pi so as to map all the angles to 0, pi]In order to ensure that mapping does not interfere with judgment, the problem can be solved by utilizing the characteristic that tangent values of angles before and after adding pi are consistent, and if the condition of the formula is met, the symmetry angle of the machine body can be judged to be the direction of the symmetry axis.
Has the advantages that: compared with the prior art, the technical scheme of the invention has the following beneficial technical effects:
(1) A target symmetry axis detection method based on a concentric circumference filter is provided. Aiming at the problem that the circumferential filter does not utilize phase information, the collinear condition of the pixel points is counted from the accumulated angle of the relative angle, the relative angle is further judged, the detection task of the symmetry axis is realized, and the rotation invariance of the subsequent characteristics is positively influenced.
Drawings
FIG. 1 is a block diagram of an embodiment of the present invention.
Detailed Description
The present invention is further illustrated by the following examples, which are intended to be purely exemplary and are not intended to limit the scope of the invention, as various equivalent modifications of the invention will occur to those skilled in the art upon reading the present disclosure and fall within the scope of the appended claims.
As shown in fig. 1, a target symmetry axis detection method based on a concentric circular filter includes the following steps:
firstly, an input target image is taken as an original image img, the original image img is filtered through a circumferential filter, corresponding frequencies can be extracted through priori knowledge, a frequency response graph can be obtained, a frequency response threshold value is preset, a series of concentric circles can be constructed by taking the frequency position larger than the frequency response threshold value as the center, frequency response is extracted from each concentric circle, the concentric circles smaller than the frequency response threshold value are discarded, and the concentric circle filter can be obtained.
Firstly, supposing that a target needing to be subjected to symmetry axis detection is an airplane target, the airplane is provided with a left wing, a right wing, a front fuselage and a rear fuselage in front of and behind the wings, and 4 components in total, so that a circumference can be constructed on an input optical image img, the circumference vector is regarded as a one-dimensional signal to be subjected to Fourier transform, whether the circumference contains the target or not can be judged according to a transform result, and all coordinates are updated to the circle center of the circumference on the img so as to judge the coordinates of the target on the img:
where π represents the circumference ratio, sin and cos represent the sine and cosine functions, respectively, and z k The method comprises the steps of representing the k-th pixel value on the circumference, representing the total number of circumferential pixels by length, representing Fourier (p, q) representing Fourier transform response corresponding to a circumferential vector with the circle center being (p, q), simultaneously representing that a response value is placed at the position (p, q) of a frequency response graph, and sequentially updating all coordinates on an original graph img as the circle center to perform response extraction, so that the frequency response graph fourier can be generated.
Then, constructing a concentric circumference filter by taking the frequency response position larger than the frequency response threshold as the center of a circle from the frequency response graph fourier:
circleFilter=[filter 0 (p,q),...,filter num (p,q)]s.t.fourier(p,q)≥FreqThres
wherein, fourier (p, q) is expressed as frequency response on coordinates (p, q), freqThres is frequency response threshold, filter 0 (p, q) as the layer 0 circumference filter with the circle center (p, q), namely the layer 0 circumference, and the same principle is that the filter num (p, q) is the num layer of circumference filter with the circle center of (p, q), namely the outermost circumference, num is the total number of circumferences, and circleFilter represents the concentric circumference filter set.
The single-layer circumference curve equation with the circle center (p, q) is as follows:
wherein (x, y) represents the coordinates of a pixel point on the circumference, fourier k (p, q) and radius k (p, q) are each a filter corresponding to a circular filter k Fourier transform response and radius of (p, q), and filter k (p, q) is the k-th layer circular filter with the center of the circle being (p, q), the frequency response of each circle in the concentric circles is greater than the frequency response threshold value FreqThres, and once the frequency response is less than the threshold value, the circle is determined to be beyond the target boundary.
And secondly, extracting an edge image edge from the original img by using a canny operator, overlapping the concentric circle filter with the edge image edge, wherein the overlapped part is a sampling point, and converting the coordinate of the sampling point from a Cartesian coordinate to a polar coordinate. Then, taking the 0 th layer of circumferential sampling points as reference points, calculating the relative angles of the sampling points relative to the reference points, measuring the similarity among the relative angles of the sampling points on different layers, and finally, regarding the similar angles as collinear angles to carry out collinear statistics.
Firstly, extracting an edge image edge from img by using a canny operator, overlapping the edge image edge with a circleFilter, wherein the overlapped part is a sampling point, and a sampling point setWhereinRepresenting a circular filter j The coordinates of the ith sampling point are (p, q), and edge (p, q) represents the pixel value at the coordinates (p, q) on the edge image edge, and similarly, the filter j (p, q) denotes a circular filter j Pixel value at coordinate (p, q).
Then, the polar coordinate representation of the sampling point can be obtained by using the following coordinate system conversion formula:
wherein,is the sampling point->Relative angle relative to the origin in a polar coordinate system, andis the sampling point->Expressed in terms of the distance from the origin in a polar coordinate system>Then it is the sampling point pick>Abscissa in the Cartesian coordinate system, for the same reason, in>Then it is a sampling point>Ordinate in cartesian coordinate system.
Then, using the sampling points on the circumference of layer 0 as reference points, the relative angles of the rest sampling points to the reference points can be calculated, and the following formula gives the calculation by taking layer 1 and layer 0 as examplesThe specific mode of (1): />
Wherein,represents the jth sampling point of the most 0 layer, <' > or>Represents the ith sample point of the most 1 layer, and->Is the relative angle between the two sample points.
Then, by measuring the relative angle similarity between the sampling points of other layers and the sampling point of the layer 1, two sampling points with the difference smaller than a similar threshold value can be determined to be collinear, and the circumferences of the collinear two layers are determined to be collinear by similar weightsMarking:
wherein,the method is characterized in that the relative angles of the qth point of the ith layer and the pth point of the 1 st layer are respectively determined, the jth point of the 0 th layer is taken as a starting point, the similarThres is taken as a similar threshold, and the above formula constraint shows that all sampling points of the ith layer are traversed, and if the sampling points of the ith layer and the/or the & lt/EN & gt exist>Is less than the similarity threshold, then it is assumed that the two layers are co-linear. And a similar weight->Indicating whether or not the i-th layer contains collinearly>The sampling points of (a). The sum of similar weights for all layers results in->Number of collinear layers in direction:
collinear layer number statistics representing the relative angle of the jth point on the 0 th layer and the pth point on the 1 st layer, however, each sampling point on the 1 st layer should belong to only one of the front fuselage, the rear fuselage, the left wing and the right wing, so after all the 0 th layer points are traversed, the most possible statistical result sW is taken under the condition of all different reference circumference starting points p 。
Wherein m represents the total number of layer 0 samples, sW p Representing the statistical result with the largest number of collinear layers of the p-th sampling point in the 1 st layer in the situation of all different reference circumference starting points, and obtaining the similarWeight = [ sW ] after the statistical results of all the points on the 1 st layer are all calculated 1 ,...,sW n ]And n represents the total number of layer 1 sample points, i.e., the statistical result set.
And finally, reserving the direction in which the number of layers is greater than the threshold value of the collinear number of layers, and regarding sampling points in the direction smaller than the threshold value as interference:
wherein,represents the statistical result that the p th collinear layer number in the layer 1 meets the condition, layerThres is a collinear layer number threshold, pick _ num is the total number of sampling points meeting the condition, and->And showing that the statistical result set of which the number of collinear layers is less than the collinear threshold value layerThres is screened out.
And thirdly, converting the angle into a tangent value, and obtaining the angle in the direction of the symmetry axis.
In the second step canTo obtainIs arranged therein>The corresponding specific similarity weight isTherefore, as can be seen from the similar weights, the set of relative angles retained by layer 1 isSince the target of the airplane is divided into 4 parts such as the front fuselage, the rear fuselage, the left wing, the right wing and the like, the relative angle set can be divided into 4 parts by the KNN method, and the relative angle can be further expressed by 4 angles, namely [ relative 1 ,...,relative 4 ]Substituting the relative angle value into the tangent function, and determining the direction of the symmetry axis by the following formula:
symmetryAngle=relative j s.t.tan(relative i )-tan(relative j )=tan(relative j )-tan(relative k )
wherein tan represents a tangent function, relative i 、relative k Respectively showing the relative angles of the left wing and the right wing j I.e. the relative angle of the fuselage, symmetryAngle is the direction of the symmetry axis, the left wing and the right wing should have the same relative angle difference because they are symmetric to the fuselage, and in the actual operation, to avoid the calculation of the negative number angle, the negative number angle is added by pi so as to map all angles to 0, pi]In order to ensure that mapping does not interfere with judgment, the problem can be solved by utilizing the characteristic that the tangent values of the angle before and after adding pi are consistent, and if the condition of the formula is met, the symmetrical angle of the machine body can be judged to be the direction of the symmetrical axis.
Claims (3)
1. A target symmetry axis detection method based on a concentric circle filter is characterized by comprising the following steps:
(1) Constructing a concentric circular filter: taking an input target image as an original image img, filtering the original image img through a circumferential filter, extracting preset corresponding frequency to obtain a frequency response image, presetting a frequency response threshold, constructing a series of concentric circles with the same circle center and gradually increased radius by taking a frequency position larger than the frequency response threshold as a center, extracting frequency response from each concentric circle, and discarding the concentric circles smaller than the frequency response threshold to obtain a concentric circle filter;
(2) Angle statistics: firstly, extracting an edge image edge from an original image img, overlapping a concentric circle filter with the edge image edge, wherein the overlapped part is a sampling point, converting the coordinate of the sampling point from a Cartesian coordinate to a polar coordinate, then calculating the relative angle of the sampling point relative to a reference point by taking the alpha-layer circle sampling point as the reference point, measuring the similarity between the relative angles of the sampling points of different layers, and finally, regarding the similar angle as a collinear angle to carry out collinear statistics;
(3) And (3) angle judgment: converting the angle into a tangent value to obtain a symmetry axis direction angle;
the angle statistical method in the step (2) is as follows:
(2.1) firstly, extracting an edge image edge from img by using a canny operator, overlapping the edge image edge with a circleFilter, wherein the overlapped part is a sampling point, and the sampling point set comprises:
wherein,representing a circular filter j The coordinates of the ith sampling point are (p, q), and edge (p, q) represents the pixel value at the coordinates (p, q) on the edge image edge, and similarly, the filter j (p, q) denotes a circular filter j A pixel value at coordinate (p, q);
(2.2) obtaining the polar coordinate representation of the sampling point by using the following coordinate system conversion formula:
wherein,as a sampling pointRelative angle relative to the origin in a polar coordinate system, andas a sampling pointThe distance from the origin in a polar coordinate system represents,is the sampling pointThe abscissa in the cartesian coordinate system, similarly,is the sampling pointOrdinate in cartesian coordinate system;
(2.3) taking the sampling points on the alpha layer circumference as reference points, and calculating the relative reference points of the rest sampling pointsThe following formula gives the calculation by taking the gamma layer and the alpha layer as an exampleThe specific mode of (1):
wherein,represents the ith sample point of the gamma-most layer,represents the jth sampling point of the alpha-most layer, andis the relative angle between the two sampling points;
(2.4) determining two sampling points with the difference smaller than a similarity threshold as collinear by measuring the relative angle similarity of the sampling points of other layers and the sampling point of the beta layer, and determining the two sampling points with the difference smaller than the similarity threshold as collinear two-layer circumferences with similar weightsMarking:
wherein,the relative angles of the ith layer of the q point and the beta layer of the p point are respectively shown, the alpha layer of the j point is used as a reference point, the similarThres is used as a similar threshold value, the above formula constraint shows that all sampling points of the ith layer are traversed, if the ith layer of the sampling points and the beta layer of the p point exist, the ith layer of the sampling points and the beta layer of the p point are used as the reference points, the similarity between the sampling points is used as the similarity threshold value, and the formula constraint shows that the ith layer of the p point is used as the similarity threshold valueIs less than the similarity threshold, then the two layers are considered co-linear and the similarity weight is determined to beIndicating whether the ith layer contains collinearity toThe similar weights of all layers are summed to obtainNumber of collinear layers in direction:
collinear layer number statistics for representing the relative angle of the jth point of the alpha layer and the pth point of the beta layer, however, each sampling point on the beta layer should belong to only one of the front fuselage, the rear fuselage, the left wing and the right wing, so after all the alpha layer points are traversed, the statistical result sW with the maximum possibility is taken under the condition of all different reference circumference starting points p :
Wherein m represents the total number of alpha layer sampling points, sW p Representing the statistical result with the maximum number of collinear layers of the p-th sampling point in the beta-th layer under the condition of all different reference circumference starting points, and obtaining the simiarweight = [ sW ] after the statistical results of all the points on the beta-th layer are all calculated 1 ,...,sW n ]N represents the beta layerThe total number of sampling points is the statistical result set;
(2.5) reserving the direction in which the number of layers is greater than the threshold value of the collinear number of layers, and regarding sampling points in the direction smaller than the threshold value as interference:
wherein,showing the statistical result that the p-th collinear layer number in the beta layer meets the condition, layerThres is a collinear layer number threshold value, pick _ num is the total number of sampling points meeting the condition,and showing that the statistical result set of which the number of collinear layers is smaller than the collinear threshold layerThres is screened out.
2. The concentric circular filter-based target symmetry axis detection method according to claim 1, wherein the method for constructing the concentric circular filter in step (1) is as follows:
(1.1) assuming that a target needing to be subjected to symmetry axis detection is an airplane target, the airplane is provided with a left wing, a right wing, a front fuselage and a rear fuselage in front of and behind the wings, and 4 components in total, so that a circumference is constructed on an input optical image img, a circumference vector is regarded as a one-dimensional signal to be subjected to Fourier transform, whether the circumference contains the target or not is judged according to a transform result, all coordinates are updated to the circle center of the circumference on the img, and the coordinates of the target on the img are judged:
where π represents the circumference ratio, sin and cos represent the sine and cosine functions, respectively, and z k Expressed as the kth pixel value on the circumference, and length expresses the circumferential pixel pointThe total number, namely Fourier (p, q) represents Fourier transform response corresponding to a circumferential vector with a circle center (p, q), and simultaneously represents that a response value is placed at the position (p, q) of the frequency response graph Fourier, and all coordinates on the original graph img are sequentially updated as the circle center to perform response extraction to generate the frequency response graph Fourier;
(1.2) constructing a concentric circumference filter by taking the frequency response position larger than the frequency response threshold as the center of a circle from the frequency response graph fourier:
circleFilter=[filter 0 (p,q),...,filter num (p,q)]s.t.fourier(p,q)≥FreqThres
wherein, fourier (p, q) is expressed as frequency response on coordinates (p, q), freqThres is frequency response threshold, filter 0 (p, q) as the layer 0 circumference filter with the circle center (p, q), namely the layer 0 circumference, and the same principle is that the filter num (p, q) is the num layer of circumference filter with the circle center of (p, q), namely the outermost circumference, num is the total number of circumferences, and circleFilter represents a concentric circumference filter set;
wherein, the single-layer circumference curve equation with the circle center being (p, q) is as follows:
wherein (x, y) represents the coordinates of a pixel point on the circumference, four k (p, q) and radius k (p, q) are each a filter corresponding to a circular filter k Fourier transform response and radius of (p, q), and filter k (p, q) is a k-th layer circle filter with the center of the circle being (p, q), the frequency response of each circle in the concentric circles is larger than a frequency response threshold value FreqThres, and if the frequency response is smaller than the threshold value, the circles are determined to be beyond the target boundary.
3. The concentric circular filter-based target symmetry axis detection method according to claim 1, wherein the angle determination in step (3) is as follows:
(3.1) obtaining the product in the step (2)Is provided thereinThe corresponding specific similarity weight isTherefore, from similar weights, the set of relative angles retained by the beta-th layer isAnd because the plane target is divided into 4 parts of the front plane body, the rear plane body, the left wing and the right wing, the relative angle set is divided into 4 parts by the KNN method, and the 4 parts are expressed by 4 angles, namely the relative 1 ,...,relative 4 ]Substituting the relative angle value into a tangent function, and determining the direction of the symmetry axis according to the following formula:
symmetryAngle=relative j s.t.tan(relative i )-tan(relative j )=tan(relative j )-tan(relative k )
wherein tan represents a tangent function, relative i 、relative k Respectively showing the relative angles of the left wing and the right wing j I.e., the relative angle of the fuselage, and symmetryAngle is the direction of the axis of symmetry, since the left wing and the right wing are symmetric to the fuselage, the relative angle difference between them should be consistent, and all negative angles are added by pi so that all angles are mapped to 0, pi]Therefore, if the above formula condition is satisfied:
s.t.tan(relative i )-tan(relative j )=tan(relative j )-tan(relative k )
therefore, the symmetry angle of the fuselage is determined to be the direction of the symmetry axis.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010434592.3A CN111723335B (en) | 2020-05-21 | 2020-05-21 | Target symmetry axis detection method based on concentric circumference filter |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010434592.3A CN111723335B (en) | 2020-05-21 | 2020-05-21 | Target symmetry axis detection method based on concentric circumference filter |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111723335A CN111723335A (en) | 2020-09-29 |
CN111723335B true CN111723335B (en) | 2023-03-24 |
Family
ID=72564854
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010434592.3A Active CN111723335B (en) | 2020-05-21 | 2020-05-21 | Target symmetry axis detection method based on concentric circumference filter |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111723335B (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106778742A (en) * | 2016-12-09 | 2017-05-31 | 东南大学 | A kind of automobile logo detection method suppressed based on Gabor filter background texture |
RU2018109177A (en) * | 2018-03-14 | 2018-04-24 | Федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия войсковой противовоздушной обороны Вооруженных Сил Российской Федерации имени Маршала Советского Союза А.М. Василевского" Министерства обороны Российской Федерации | A method for estimating the spatial size of an air target by the frequency extent of a Doppler portrait |
-
2020
- 2020-05-21 CN CN202010434592.3A patent/CN111723335B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106778742A (en) * | 2016-12-09 | 2017-05-31 | 东南大学 | A kind of automobile logo detection method suppressed based on Gabor filter background texture |
RU2018109177A (en) * | 2018-03-14 | 2018-04-24 | Федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия войсковой противовоздушной обороны Вооруженных Сил Российской Федерации имени Маршала Советского Союза А.М. Василевского" Министерства обороны Российской Федерации | A method for estimating the spatial size of an air target by the frequency extent of a Doppler portrait |
Non-Patent Citations (2)
Title |
---|
"Aircraft Detection of High-Resolution Remote Sensing Image Based on Faster R-CNN Model and SSD Model";Xie Qifang 等;《https://doi.org/10.1145/3191442.3191443》;20180224;全文 * |
"基于稀疏表示的可见光遥感图像飞机检测算法";林煜东 等;《光子学报》;20140513;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111723335A (en) | 2020-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110175615B (en) | Model training method, domain-adaptive visual position identification method and device | |
CN107643181B (en) | A kind of rolling bearing variable working condition method for diagnosing faults based on image recognition | |
CN109635733B (en) | Parking lot and vehicle target detection method based on visual saliency and queue correction | |
CN110018453B (en) | Intelligent model identification method based on aircraft track characteristics | |
CN102298698A (en) | Remote sensing image airplane detection method based on fusion of angle points and edge information | |
CN111340855A (en) | Road moving target detection method based on track prediction | |
CN111369526B (en) | Multi-type old bridge crack identification method based on semi-supervised deep learning | |
Liang et al. | Global context-aware cervical cell detection with soft scale anchor matching | |
CN109558803B (en) | SAR target identification method based on convolutional neural network and NP criterion | |
CN110889399A (en) | High-resolution remote sensing image weak and small target detection method based on deep learning | |
CN108846414A (en) | SAR image subclass classification method based on decision level fusion thought | |
CN111723335B (en) | Target symmetry axis detection method based on concentric circumference filter | |
Yildirim et al. | Ship detection in optical remote sensing images using YOLOv4 and Tiny YOLOv4 | |
CN109145770B (en) | Automatic wheat spider counting method based on combination of multi-scale feature fusion network and positioning model | |
CN111898589B (en) | Unmanned aerial vehicle image rapid registration method based on GPU+feature recognition | |
CN113962324A (en) | Picture detection method and device, storage medium and electronic equipment | |
CN111832463A (en) | Deep learning-based traffic sign detection method | |
CN113298155B (en) | Airport runway underground disease detection method based on SF-SSD algorithm | |
CN111127485B (en) | Method, device and equipment for extracting target area in CT image | |
CN114882375A (en) | Intelligent identification method and device for tailing pond | |
Yang et al. | Edge detection in Cassini astronomy image using Extreme Learning Machine | |
CN113850252A (en) | Method and device for detecting vibration damper bolt and electronic equipment | |
CN109472766B (en) | Bridge bolt area positioning method and terminal equipment | |
CN113553899B (en) | Gun gesture perception recognition method and system | |
Wei et al. | SCV-UNet: Saliency-Combined Complex-Valued U-Net for SAR Ship Target Segmentation |
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 |