CN113763463A - Method for determining position of acquisition equipment based on image data processing - Google Patents

Method for determining position of acquisition equipment based on image data processing Download PDF

Info

Publication number
CN113763463A
CN113763463A CN202111323650.6A CN202111323650A CN113763463A CN 113763463 A CN113763463 A CN 113763463A CN 202111323650 A CN202111323650 A CN 202111323650A CN 113763463 A CN113763463 A CN 113763463A
Authority
CN
China
Prior art keywords
marker
equation
pixel
point
formula
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111323650.6A
Other languages
Chinese (zh)
Other versions
CN113763463B (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.)
Fengmai Changzhong (Nantong) New Energy Co.,Ltd.
Original Assignee
Windmagics Wuhan Co ltd
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 Windmagics Wuhan Co ltd filed Critical Windmagics Wuhan Co ltd
Priority to CN202111323650.6A priority Critical patent/CN113763463B/en
Publication of CN113763463A publication Critical patent/CN113763463A/en
Application granted granted Critical
Publication of CN113763463B publication Critical patent/CN113763463B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Image Analysis (AREA)

Abstract

The invention relates to a method for determining the position of a collecting device based on image data processing, comprising the following steps: the method comprises the following steps of installing markers in a cylinder to ensure that equipment can acquire three markers in the same digital image; identifying pixel position and pixel width information of the marker; calculating the angle of the marker through a relation function of the pixel of the equipment and the angle; establishing a vector relation of the triangle, and performing complex vector operation on the vector relation to obtain a first equation; performing trigonometric function operation according to the imaging relation of the markers to obtain a second equation to form an equation set; and (4) bringing the pixel positions and the pixel width information of the three markers into the three groups of equations to solve the position of the acquisition equipment. According to the invention, the image acquisition is carried out through the equipment, the eccentric direction and the eccentric distance of the image acquisition equipment can be obtained by reading the position and the pixel size of the label in the image, and meanwhile, the physical size of the inner diameter of the cylinder can be obtained without any physical measurement of constructors.

Description

Method for determining position of acquisition equipment based on image data processing
Technical Field
The invention belongs to the technical field of image data processing, and particularly relates to a method for determining the position of acquisition equipment based on image data processing.
Background
Measurement in cylindrical structures using ultra wide angle/panoramic digital technology is an emerging measurement method, for example, measurement of foundation sway inside a tower foundation ring of a wind turbine generator system is an example of such measurement. In the measurement of the tower structure by using the panoramic digital technology, the image acquisition device is often required to be placed on the central axis of the tower structure, so that the object distance of the markers on the inner wall of the tower is consistent compared with the center of the tower.
The measurement usually needs to be carried out at the center (circle center) of the foundation ring, but the limitation of field physical conditions is avoided, the requirement is not guaranteed, and other parts are sometimes arranged at the center position, so that the image acquisition equipment cannot be placed; sometimes, the obstacle is arranged near the circle center to block the view, so that the image acquisition information is incomplete, and the position needs to be adjusted; sometimes it is difficult to determine where the center location is specific, simply because of field personnel and condition limitations. Due to the above constraints, the measuring device is often placed at a position offset from the geometric center of the cylinder, and such offset may cause the object distance from the acquisition device to each marker to be inconsistent, resulting in measurement errors.
In addition, the radius of the cylinder structure often comes from project construction and construction data, and measuring personnel often can not obtain at the first site, and sometimes project information only marks external diameter size, can not obtain internal diameter size. Thus, when making panoramic digital measurements within the cylindrical structure, the eccentric position of the acquisition device and the inside diameter of the cylindrical object need to be obtained to be able to determine the precise physical location of other markers within the cylindrical structure.
Disclosure of Invention
The invention aims to provide a method for determining the position of acquisition equipment based on image data processing. According to the scheme of the invention, no physical measurement is required for constructors, the constructors can find a proper position on site to install the image acquisition equipment, and can leave the construction site after confirming that the sight line of the label is not blocked. The input quantity of the method is completely from the digital image, the whole process does not need manual intervention, the requirement of site operators on equipment positioning is reduced, and the technical requirement on site operators is greatly reduced.
In order to achieve the purpose, the invention provides the following technical scheme:
a method for determining the position of an acquisition device based on image data processing, characterized by:
the inner wall of the cylinder is provided with a plurality of markers, so that more than three markers can be acquired by the image acquisition equipment in the same digital image;
identifying pixel position and pixel width information of a marker in a digital image;
calculating the angle of the marker through a relation function between the pixel position of the image acquisition equipment and the incident angle of the light;
establishing a mathematical model, wherein the center of the cylinder, the position of the image acquisition equipment and any marker form a triangle; establishing a vector relation of the triangle, and performing complex vector operation on the vector relation to obtain a first equation;
performing trigonometric function operation according to the imaging relation of the marker to obtain an equation two, and forming an equation set with the equation one;
and (3) at least bringing in the pixel position and pixel width information of the three markers to obtain an equation set consisting of three sets of equations I and II, and solving to obtain the position of the acquisition equipment.
Further, the marker is a rectangular marker, and the pixel position of the marker is represented by a pixel of the transverse span center point or a pixel of the geometric center point of the marker.
Further, according to the optical physical characteristics of the image acquisition device, the relation function between the pixel position and the incident angle of the light ray is known
Figure 426843DEST_PATH_IMAGE001
Passing function
Figure 496430DEST_PATH_IMAGE002
Obtaining pixels in a digital image
Figure 553248DEST_PATH_IMAGE003
Angle value of
Figure 366483DEST_PATH_IMAGE004
Figure 457936DEST_PATH_IMAGE004
The central angle using the image acquisition device as the center of a circle is as follows:
Figure 432845DEST_PATH_IMAGE005
(formula 2)
Angle of view of marker A
Figure 242538DEST_PATH_IMAGE006
Comprises the following steps:
Figure 593885DEST_PATH_IMAGE007
(formula 5)
Wherein,
Figure 539844DEST_PATH_IMAGE008
is the pixel location of marker a.
Further, the pixel width occupied by the marker a is:
Figure 685655DEST_PATH_IMAGE009
(formula 1)
Wherein, the pixel of the start position of the marker A is recorded as
Figure 717065DEST_PATH_IMAGE010
End position pixel is noted
Figure 872102DEST_PATH_IMAGE011
Further, the pass function
Figure 672568DEST_PATH_IMAGE002
The angle value of the initial position in the width direction of the marker A can be obtained
Figure 254859DEST_PATH_IMAGE012
And end position angle value
Figure 773565DEST_PATH_IMAGE013
The method comprises the following steps:
Figure 466715DEST_PATH_IMAGE014
(formula 3)
According to equation 3, the angle of view occupied by marker a is:
Figure 528212DEST_PATH_IMAGE015
(formula 4).
Further, the obtaining of the first equation includes: taking a marker A to establish a mathematical model, wherein the marker A is represented by a point A, the center of the cylinder is a point O, and the position of the equipment is a point P
Figure 406038DEST_PATH_IMAGE016
In a coordinate system with the O point as a center, the vector expression of the marker A is as follows:
Figure 21827DEST_PATH_IMAGE017
wherein
Figure 643301DEST_PATH_IMAGE018
Is a central angle with O point as the center of a circle, and R is a cylinder structureA radius;
in a coordinate system with the P point as a center, the vector expression of the marker A is as follows:
Figure 293725DEST_PATH_IMAGE019
wherein
Figure 342453DEST_PATH_IMAGE020
Is a central angle taking the point P as the center of a circle,
Figure 711117DEST_PATH_IMAGE021
the object distance of the point A compared with the point P is shown;
the position of the point P is expressed by complex plane coordinates (a, b), the origin is the point O, a represents an abscissa in the complex plane, and b represents an ordinate in the complex plane;
according to
Figure 870703DEST_PATH_IMAGE022
The vector relationship of (a) can be found as:
Figure 641213DEST_PATH_IMAGE023
(formula 6)
The complex vector operation is performed on equation 6
Figure 595263DEST_PATH_IMAGE024
(formula 7)
Wherein,
Figure 451223DEST_PATH_IMAGE021
the object distance of point a compared to point P,
Figure 414500DEST_PATH_IMAGE025
is an imaginary unit;
the imaginary part and the real part in equation 7 are respectively given by the following equations:
Figure 39516DEST_PATH_IMAGE026
(formula 8)
Adding the two sides of equation 8 after squaring to eliminate variable
Figure 430046DEST_PATH_IMAGE027
Obtaining:
Figure 507724DEST_PATH_IMAGE028
(formula 9)
Equation 9 is equation one.
Further, the obtaining of the first equation includes:
according to the angle relationship, a triangle can be obtained
Figure 9112DEST_PATH_IMAGE016
In
Figure 488635DEST_PATH_IMAGE029
Comprises the following steps:
Figure 191012DEST_PATH_IMAGE030
(formula 10)
From the imaging relationship of marker a, the equation is derived:
Figure 880619DEST_PATH_IMAGE031
(formula 11)
Wherein W is the actual width of the label A
After deriving from the cosine theorem:
Figure 61065DEST_PATH_IMAGE032
(formula 12)
The formula 10 and the formula 12 are replaced by the formula 11, and the formula is simplified:
Figure 785307DEST_PATH_IMAGE033
(formula 13)
Equation 13 is equation two.
Further, the solving for the location of the acquisition device comprises:
combining equation 9 and equation 13 yields a system of equations:
Figure 393006DEST_PATH_IMAGE034
(formula 14)
Pixel locations brought into three markers A, B, C
Figure 304330DEST_PATH_IMAGE035
Figure 288467DEST_PATH_IMAGE036
Figure 867216DEST_PATH_IMAGE037
And pixel width information
Figure 645816DEST_PATH_IMAGE038
Can be solved out
Figure 44436DEST_PATH_IMAGE039
Six unknowns.
Further, two sides of formula 2 are subjected to derivation to obtain:
Figure 566684DEST_PATH_IMAGE040
(formula 15)
Bringing label a into formula 15, having:
Figure 999940DEST_PATH_IMAGE041
(formula 16)
Substituting equations 3 and 15 into equation 14 can be further simplified as:
Figure 949441DEST_PATH_IMAGE042
(formula 17).
Further, locating other markers in the image includes: the pixel position of a certain marker N in an image is known as
Figure 366516DEST_PATH_IMAGE043
The angle value is:
Figure 426876DEST_PATH_IMAGE044
(formula 18)
Bringing formula 18 into formula 9 can result:
Figure 980217DEST_PATH_IMAGE045
(formula 19)
Solving for the object distance of the available marker N
Figure 100620DEST_PATH_IMAGE046
Compared with the prior art, the invention has the following beneficial effects: according to the invention, the image acquisition is carried out through the equipment, and the eccentric direction and the eccentric distance of the image acquisition equipment can be obtained by reading the position and the pixel size of the label in the image, and meanwhile, the physical size of the inner diameter of the cylinder can also be obtained. According to the scheme of the invention, no physical measurement is required for constructors, the constructors can find a proper position on site to install the image acquisition equipment, and can leave the construction site after confirming that the sight line of the label is not blocked. The input quantity of the method is completely from the digital image, the whole process does not need manual intervention, the requirement of site operators on equipment positioning is reduced, and the technical requirement on site operators is greatly reduced. In addition, the method can obtain the size of the cylinder structure, reduce the investment requirement and simultaneously improve the reliability of data.
Drawings
Fig. 1 is a flowchart of a method of determining a position of an acquisition device based on image data processing in an embodiment.
FIG. 2 is a schematic view of a rectangular tag in an embodiment.
FIG. 3 is a schematic structural relationship diagram of an image capturing device and a tower in an embodiment.
Fig. 4 is a schematic image acquired in the embodiment.
FIG. 5 is a schematic diagram illustrating a relationship between a center of a tower, a position of an acquisition device, and a position of a marker A according to an embodiment.
FIG. 6 is a diagram illustrating the position relationship of an exemplary triangular APO in an embodiment.
FIG. 7 is a schematic view showing the imaging relationship of the marker A in the example.
Fig. 8 is an actual view digital image of the image capturing device according to the first embodiment.
FIG. 9 is a diagram of a geometric model of a digital image according to an embodiment.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to specific embodiments. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In the measurement of the tower structure by using the panoramic digital technology, the image acquisition device is often required to be placed on the central axis of the tower structure, so that the object distance of the markers on the inner wall of the tower is consistent compared with the center of the tower. However, this requirement is often not guaranteed due to the limitations of field operating conditions. Referring to fig. 1, a method for determining the position of the acquisition device based on image data processing is described below, which can obtain the eccentric position of the acquisition device and the internal dimension of the tower structure by means of digital image technology only, and the specific technical route is as follows.
Step one, placing a marker with known width
As shown in fig. 3, at least 3 markers (A, B, C, D) are mounted on the inner wall of the tower structure, so that the image capturing device (such as a camera) can capture more than three markers in the same frame, wherein the markers A, B, C, D are all in the same frame. Only this step requires on-site installation inside the tower;
the shape of the marker is a graphic of known specific geometric dimensions, in the case of the rectangular marker in fig. 2, the width of the marker is denoted as W, a known amount.
Step two, identifying the pixel position and the pixel width information of the marker
In digital images acquired by image acquisition equipment (taking a camera as an example, the equipment for short), performing plane expansion on the acquired digital images to obtain rectangular images inside a tower, and then constructing a geometric model, as shown in fig. 3, extracting pixel position and pixel width information of a marker;
for the pixel position of the marker A, one pixel point can be used
Figure 473832DEST_PATH_IMAGE008
Generally, the lateral span center point, or geometric center point, of the marker is used;
as shown in fig. 4, even for markers of the same physical size, there may be differences in the pixel widths that they represent in the image due to the influence of eccentricity of the image capturing device. In the width direction, the pixel of the start position of the marker A is recorded
Figure 337883DEST_PATH_IMAGE010
(obtained by automatic recognition or manual reading) and the pixels at the end positions are recorded as
Figure 745731DEST_PATH_IMAGE011
The pixel width information is known;
the pixel width occupied by marker a (span) is:
Figure 37035DEST_PATH_IMAGE009
(formula 1).
Step three, calculating the angle of the marker through the lens parameters
From the optical characteristics of the camera lens, the relationship function between the pixel position and the incident angle of the light ray is known
Figure 897543DEST_PATH_IMAGE001
Passing function
Figure 565285DEST_PATH_IMAGE002
The pixels in the digital image can be known
Figure 827639DEST_PATH_IMAGE003
Angle value of
Figure 289845DEST_PATH_IMAGE004
The following formula:
Figure 778595DEST_PATH_IMAGE005
(formula 2)
Angle value of starting position in width direction of marker A
Figure 374661DEST_PATH_IMAGE012
And end position angle value
Figure 366888DEST_PATH_IMAGE013
The method comprises the following steps:
Figure 124628DEST_PATH_IMAGE014
(formula 3)
The angle of view occupied (spanned) by marker a is according to equation 3 (shown in fig. 4):
Figure 835095DEST_PATH_IMAGE015
(formula 4)
Angle of view of marker A
Figure 969274DEST_PATH_IMAGE006
By pixel points
Figure 816007DEST_PATH_IMAGE008
Taking the center pixel of the span of A pixels of the marker as a reference point
Figure 10228DEST_PATH_IMAGE008
The calculation is as follows:
Figure 942412DEST_PATH_IMAGE007
(formula 5)
The central pixel of the marker A can be obtained by automatic identification or manual reading
Figure 880281DEST_PATH_IMAGE008
Step four, complex vector operation
As shown in FIG. 5, a mathematical model is established by taking a marker A on the tower structure, the marker A is represented by point A, the center of the tower is point O, and the equipment position is point P
Figure 847100DEST_PATH_IMAGE016
In a coordinate system with the O point as a center, the vector expression of the marker A is as follows:
Figure 212222DEST_PATH_IMAGE017
wherein
Figure 366123DEST_PATH_IMAGE018
Is a central angle taking the point O as the center of a circle, and R is the radius of the tower barrel structure;
in a coordinate system with the P point as a center, the vector expression of the marker A is as follows:
Figure 107683DEST_PATH_IMAGE019
wherein
Figure 929008DEST_PATH_IMAGE020
Is a circle with P point as the centerThe position of the heart angle, point P, can also be expressed in complex plane coordinates (a, b) (origin is point O);
the object distance of point A compared with point P is
Figure 465032DEST_PATH_IMAGE021
According to
Figure 371808DEST_PATH_IMAGE022
The vector relationship of (a) can be found as:
Figure 651479DEST_PATH_IMAGE023
(formula 6)
The complex vector operation is performed on equation 6
Figure 327311DEST_PATH_IMAGE047
(formula 7)
The imaginary part of (1) in equation 7
Figure 34236DEST_PATH_IMAGE025
In units of imaginary numbers), and the real part are respectively given by the following equations:
Figure 693888DEST_PATH_IMAGE026
(formula 8)
Because of the fact that
Figure 387037DEST_PATH_IMAGE048
For the known quantity obtained by pixel conversion, the two sides of the equation of the formula 8 are squared and added to eliminate the variable
Figure 307589DEST_PATH_IMAGE027
Obtaining:
Figure 795202DEST_PATH_IMAGE049
modifying the above equation yields:
Figure 801204DEST_PATH_IMAGE028
(formula 9)
Equation one (equation 9) is obtained for subsequent calculations.
Step five, imaging relation of marker A
According to the angle relationship, a triangle can be obtained
Figure 298044DEST_PATH_IMAGE016
In
Figure 73102DEST_PATH_IMAGE029
Is (as shown in figure 5):
Figure 997196DEST_PATH_IMAGE030
(formula 10)
Since the digital image of the marker a is acquired by using a reference system with a point P as a center, that is, the projection length of the marker in the AP vertical direction is acquired, but the physical size of the marker a is established by using a reference system with a point O as a center, the two reference systems are connected to obtain an equation:
Figure 224915DEST_PATH_IMAGE031
(formula 11)
Thus, to the left of equation 11 is the marker width multiplied by the cosine of the projection view angle and to the right is the marker distance multiplied by the radian spanned by the label. Although the two methods are different, the apparent width (the width which should be seen by taking P as the center) of the marker is calculated;
in the above actual measurement, the marker is attached to the inner wall of the tower, as shown in fig. 7, the arc M1M2 corresponds to the marker, the length is W, and since the radii OA and W are perpendicular and R is much greater than the width W of the marker, the arc length M1M2 is approximately equal to the line segment M1M2 and is equal to W;
as shown in FIG. 7, two endpoints M of the marker1. The projection points of M2 in the AP perpendicular direction are N1 and N2 respectively. The left side of equation 11 is projected to calculate the length of N1N2, i.e.
Figure 259867DEST_PATH_IMAGE050
(ii) a On the left side of equation 11, since the distance of La is much greater than the length of N1N2 in the triangle N1PN2, N1N2= La × radius (angle N1PN2) can be obtained from the approximate relationship between the arc length and the radian, and since the distance of La is much greater than N1N2, angle M1PM2 is approximately equal to angle N1PN2, and the radian of angle N1PN2 is equal to angle N1PN2
Figure 155011DEST_PATH_IMAGE051
Therefore, there are
Figure 250006DEST_PATH_IMAGE052
In view of the above, it can be seen that,
Figure 965021DEST_PATH_IMAGE053
step six, performing trigonometric function operation
As shown in fig. 6, according to the cosine theorem:
Figure 69243DEST_PATH_IMAGE032
(formula 12)
The formula 10 and the formula 12 are replaced by the formula 11, and the formula is simplified:
Figure 818893DEST_PATH_IMAGE033
(formula 13)
Equation two (equation 13) is obtained for subsequent calculations.
Seventhly, solving global information by equation I and equation II
Combining equation 9 and equation 13 yields a system of equations:
Figure 84790DEST_PATH_IMAGE034
(formula 14)
Read from time to timeGet
Figure 287101DEST_PATH_IMAGE012
Figure 663856DEST_PATH_IMAGE013
Two readings will be in error, and to further reduce the error, the following can be used for optimization:
the derivation is performed on both sides of equation 2 to obtain the corresponding relationship of the angle change caused by the pixel change at a certain point position (the index mark) in the digital image, and the following are:
Figure 268012DEST_PATH_IMAGE040
(formula 15)
Bringing label a into formula 15, having:
Figure 704810DEST_PATH_IMAGE041
(formula 16)
Wherein
Figure 394417DEST_PATH_IMAGE054
Is equal to
Figure 574863DEST_PATH_IMAGE055
Indicating the angle of the viewing angle occupied by the marker a;
Figure 33526DEST_PATH_IMAGE056
the image processing method includes the steps that small changes of pixels of a marker A are shown, and the change of the angle of a visual angle of the marker A can be caused;
substituting equations 3 and 15 into equation 14 can be further simplified as:
Figure 641225DEST_PATH_IMAGE042
(formula 17)
By substituting the pixel position of 3 markers A, B, C
Figure 818128DEST_PATH_IMAGE035
Figure 802265DEST_PATH_IMAGE036
Figure 115434DEST_PATH_IMAGE037
And a pixel width value
Figure 894035DEST_PATH_IMAGE057
Figure 699179DEST_PATH_IMAGE058
Figure 346061DEST_PATH_IMAGE059
The above equation sets are combined to obtain six equations in three sets, i.e. a, b, R,
Figure 654683DEST_PATH_IMAGE060
Figure 752256DEST_PATH_IMAGE061
Figure 779118DEST_PATH_IMAGE062
Six unknowns.
Step eight, solving the positions of other points in the image
After the eccentric positions a and b and the tower radius R are calculated in the above steps, the object distance of any other markers in the image can be derived through vector operation. For example, the pixel position of a certain point N in the image is known to be
Figure 964111DEST_PATH_IMAGE043
Point N is within the view angle around P, the angle values are:
Figure 127239DEST_PATH_IMAGE044
(formula 18)
Bringing formula 18 into formula 9 can result:
Figure 637855DEST_PATH_IMAGE063
(formula 19)
Solving the above quadratic equation to obtain the object distance of point N
Figure 745488DEST_PATH_IMAGE046
The central angle may refer to an angle corresponding to a rotation in a counterclockwise direction with a certain direction fixed as a starting 0 °.
Description of the principle and effects: the method is based on analyzing and calculating the pixel position and the pixel width of the marker in the image, and directly measuring the position of the acquisition equipment. The method has the following innovative effects:
(1) deducing a physical space size by using the relation between the pixel position and the pixel width of the marker in the digital image, determining the geometric relation between the acquisition equipment and the tower drum, and determining the geometric position of the image acquisition equipment in the tower drum;
(2) the position of the image acquisition equipment can be obtained only by the physical position information of 3 or more markers, and the radius information of the tower can be obtained at the same time;
(3) and based on the digital image information, uniquely determining the specific position of the image acquisition equipment in the tower drum structure, and describing the position of the equipment by adopting the deviation direction angle and the deviation amount.
Example one
As shown in fig. 8 and 9, the panoramic camera position is denoted as P, and the pixels of the center points of the 3 markers A, B, C on the inner wall of the tower in the digital image are obtained, and the pixel positions and pixel widths of the marker A, B, C are shown in the following table.
Pixel position and pixel width information table for marker A, B, C
Figure 609539DEST_PATH_IMAGE064
1. The tower radius R, the eccentricity position (a, b) is solved.
From A, B, C pixel position and pixel width parameter pairs for three points:
Figure 17387DEST_PATH_IMAGE065
Figure 308691DEST_PATH_IMAGE066
Figure 310145DEST_PATH_IMAGE067
the radius R =2110 mm of the tower can be obtained by respectively substituting the equation 17, the eccentric position a = 310mm and the eccentric position b = 100mm of the device, and the converted angle and amplitude are
Figure 836941DEST_PATH_IMAGE068
Object distances of A, B, C are 1996mm, 1785mm and 2242 mm.
2. And solving the position of a new marker N on the tower structure.
Given a pixel position of 1879p and a pixel width of 16.83pixels for point N, the solution is given by substituting equation 19, and the object distance of N is 1845 mm.
As in fig. 9, all position parameters for each marker table:
parameter list of marker A, B, C
Figure 99295DEST_PATH_IMAGE069
Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.

Claims (10)

1. A method for determining the position of an acquisition device based on image data processing, characterized by:
the inner wall of the cylinder is provided with a plurality of markers, so that more than three markers can be acquired by the image acquisition equipment in the same digital image;
identifying pixel position and pixel width information of a marker in a digital image;
calculating the angle of the marker through a relation function between the pixel position of the image acquisition equipment and the incident angle of the light;
establishing a mathematical model, wherein the center of the cylinder, the position of the image acquisition equipment and any marker form a triangle; establishing a vector relation of the triangle, and performing complex vector operation on the vector relation to obtain a first equation;
performing trigonometric function operation according to the imaging relation of the marker to obtain an equation two, and forming an equation set with the equation one;
and (3) at least bringing in the pixel position and pixel width information of the three markers to obtain an equation set consisting of three sets of equations I and II, and solving to obtain the position of the acquisition equipment.
2. The method for determining a position of an acquisition device based on image data processing of claim 1, wherein: the marker is a rectangular marker, and the pixel position of the marker is represented by a pixel of the transverse span center point or a pixel of the geometric center point of the marker.
3. The method for determining a position of an acquisition device based on image data processing of claim 1, wherein: according to the optical physical characteristics of the image acquisition equipment, the relation function between the position of the pixel and the incident angle of the light ray is known
Figure 8657DEST_PATH_IMAGE001
Passing function
Figure 929209DEST_PATH_IMAGE001
Obtaining pixels in a digital image
Figure 682401DEST_PATH_IMAGE002
Angle value of
Figure 422824DEST_PATH_IMAGE003
Figure 919665DEST_PATH_IMAGE003
The central angle using the image acquisition device as the center of a circle is as follows:
Figure 694723DEST_PATH_IMAGE004
angle of view of marker A
Figure 618816DEST_PATH_IMAGE005
Comprises the following steps:
Figure 112114DEST_PATH_IMAGE006
wherein,
Figure 147067DEST_PATH_IMAGE007
is the pixel location of marker a.
4. The method of determining a position of an acquisition device based on image data processing of claim 3, wherein: the pixel width occupied by marker a is:
Figure 776631DEST_PATH_IMAGE008
wherein the initial position of the marker APixel is expressed as
Figure 871626DEST_PATH_IMAGE009
End position pixel is noted
Figure 852220DEST_PATH_IMAGE010
5. The method of determining a position of an acquisition device based on image data processing of claim 4, wherein: passing function
Figure 425284DEST_PATH_IMAGE001
The angle value of the initial position in the width direction of the marker A can be obtained
Figure 174934DEST_PATH_IMAGE011
And end position angle value
Figure 440831DEST_PATH_IMAGE012
The method comprises the following steps:
Figure 643142DEST_PATH_IMAGE013
according to equation 3, the angle of view occupied by marker a is:
Figure 285476DEST_PATH_IMAGE014
6. the method of determining a position of an acquisition device based on image data processing of claim 5, wherein: the obtaining of the first equation comprises: taking a marker A to establish a mathematical model, wherein the marker A is represented by a point A, the center of the cylinder is a point O, and the position of the equipment is a point P
Figure 889632DEST_PATH_IMAGE015
In a coordinate system with the O point as a center, the vector expression of the marker A is as follows:
Figure 592009DEST_PATH_IMAGE016
wherein
Figure 156983DEST_PATH_IMAGE017
Is a central angle taking the point O as the center of a circle, and R is the radius of the cylinder structure;
in a coordinate system with the P point as a center, the vector expression of the marker A is as follows:
Figure 196483DEST_PATH_IMAGE018
wherein
Figure 796091DEST_PATH_IMAGE019
Is a central angle taking the point P as the center of a circle,
Figure 528424DEST_PATH_IMAGE020
the object distance of the point A compared with the point P is shown;
the position of the point P is expressed by complex plane coordinates (a, b), the origin is the point O, a represents an abscissa in the complex plane, and b represents an ordinate in the complex plane;
according to
Figure 580694DEST_PATH_IMAGE021
The vector relationship of (a) can be found as:
Figure 423885DEST_PATH_IMAGE022
the complex vector operation is performed on equation 6
Figure 143579DEST_PATH_IMAGE023
Wherein,
Figure 781234DEST_PATH_IMAGE024
is an imaginary unit;
the imaginary part and the real part in equation 7 are respectively given by the following equations:
Figure 586379DEST_PATH_IMAGE025
adding the two sides of equation 8 after squaring to eliminate variable
Figure 233261DEST_PATH_IMAGE026
Obtaining:
Figure 541882DEST_PATH_IMAGE027
equation 9 is equation one.
7. The method of determining a position of an acquisition device based on image data processing of claim 6, wherein: the obtaining of the first equation comprises:
according to the angle relationship, a triangle can be obtained
Figure 616018DEST_PATH_IMAGE015
In
Figure 377300DEST_PATH_IMAGE028
Comprises the following steps:
Figure 827873DEST_PATH_IMAGE029
from the imaging relationship of marker a, the equation is derived:
Figure 991001DEST_PATH_IMAGE030
wherein W is the actual width of the label A
After deriving from the cosine theorem:
Figure 236038DEST_PATH_IMAGE031
the formula 10 and the formula 12 are replaced by the formula 11, and the formula is simplified:
Figure 750196DEST_PATH_IMAGE032
equation 13 is equation two.
8. The method for determining a position of an acquisition device based on image data processing of claim 7, wherein: solving for the location of the acquisition device comprises:
combining equation 9 and equation 13 yields a system of equations:
Figure 762318DEST_PATH_IMAGE033
pixel locations brought into three markers A, B, C
Figure 779952DEST_PATH_IMAGE034
Figure 461470DEST_PATH_IMAGE035
Figure 462924DEST_PATH_IMAGE036
And pixel width information
Figure 865086DEST_PATH_IMAGE037
Can be solved out
Figure 127440DEST_PATH_IMAGE038
Six unknowns.
9. The method for determining a position of an acquisition device based on image data processing of claim 8, wherein: and (3) performing derivation on two sides of the formula 2 to obtain:
Figure 855225DEST_PATH_IMAGE039
bringing label a into formula 15, having:
Figure 937450DEST_PATH_IMAGE040
substituting equations 3 and 15 into equation 14 can be further simplified as:
Figure 408883DEST_PATH_IMAGE041
10. the method for determining a position of an acquisition device based on image data processing of claim 9, wherein: locating other markers in the image, comprising: the pixel position of a certain marker N in an image is known as
Figure 525743DEST_PATH_IMAGE042
The angle value is:
Figure 158850DEST_PATH_IMAGE043
bringing formula 18 into formula 9 can result:
Figure 993951DEST_PATH_IMAGE044
solving for the object distance of the available marker N
Figure 3495DEST_PATH_IMAGE045
CN202111323650.6A 2021-11-10 2021-11-10 Method for determining position of acquisition equipment based on image data processing Active CN113763463B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111323650.6A CN113763463B (en) 2021-11-10 2021-11-10 Method for determining position of acquisition equipment based on image data processing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111323650.6A CN113763463B (en) 2021-11-10 2021-11-10 Method for determining position of acquisition equipment based on image data processing

Publications (2)

Publication Number Publication Date
CN113763463A true CN113763463A (en) 2021-12-07
CN113763463B CN113763463B (en) 2022-02-11

Family

ID=78784900

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111323650.6A Active CN113763463B (en) 2021-11-10 2021-11-10 Method for determining position of acquisition equipment based on image data processing

Country Status (1)

Country Link
CN (1) CN113763463B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120327220A1 (en) * 2011-05-31 2012-12-27 Canon Kabushiki Kaisha Multi-view alignment based on fixed-scale ground plane rectification
US20130291749A1 (en) * 2010-05-07 2013-11-07 Manfred Loddenkoetter Method and apparatus for determining and setting an optimized operating distance between at least two cylinders involved in a printing process
CN105371816A (en) * 2014-08-29 2016-03-02 同济大学 Wireless inclination angle measurement system applied to wind driven generator tower cylinder detection
CN109559354A (en) * 2018-11-30 2019-04-02 北京金风科创风电设备有限公司 Method for measuring tower clearance, calibration method and device of binocular system
CN110687583A (en) * 2019-11-05 2020-01-14 中国计量科学研究院 Position energy time testing system and device based on CZT detector
CN111721238A (en) * 2020-07-22 2020-09-29 上海图漾信息科技有限公司 Depth data measuring apparatus and target object data collecting method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130291749A1 (en) * 2010-05-07 2013-11-07 Manfred Loddenkoetter Method and apparatus for determining and setting an optimized operating distance between at least two cylinders involved in a printing process
US20120327220A1 (en) * 2011-05-31 2012-12-27 Canon Kabushiki Kaisha Multi-view alignment based on fixed-scale ground plane rectification
CN105371816A (en) * 2014-08-29 2016-03-02 同济大学 Wireless inclination angle measurement system applied to wind driven generator tower cylinder detection
CN109559354A (en) * 2018-11-30 2019-04-02 北京金风科创风电设备有限公司 Method for measuring tower clearance, calibration method and device of binocular system
CN110687583A (en) * 2019-11-05 2020-01-14 中国计量科学研究院 Position energy time testing system and device based on CZT detector
CN111721238A (en) * 2020-07-22 2020-09-29 上海图漾信息科技有限公司 Depth data measuring apparatus and target object data collecting method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
REZA MAALEK.ET.: "Robust detection of non-overlapping ellipses from points with applications to circular target extraction in images and cylinder detection in point clouds", 《ISPRS JOURNAL OF PHOTOGRAMMETRY AND REMOTE SENSING》 *
李思亮等: "基于SCADA数据的风机功率特性变化", 《第五届中国风电后市场专题研讨会论文集》 *

Also Published As

Publication number Publication date
CN113763463B (en) 2022-02-11

Similar Documents

Publication Publication Date Title
CN108921901B (en) Large-view-field camera calibration method based on precise two-axis turntable and laser tracker
CN102364299B (en) Calibration technology for multiple structured light projected three-dimensional profile measuring heads
CN107167169B (en) Readings of pointer type meters based on NI Vision Builder for Automated Inspection identifies measurement method
CN108828554B (en) Coordinate transformation-based measuring method, system and device without laser drop point
JP3728900B2 (en) Calibration method and apparatus, and calibration data generation method
CN104200086A (en) Wide-baseline visible light camera pose estimation method
US20130113897A1 (en) Process and arrangement for determining the position of a measuring point in geometrical space
CN106197292B (en) A kind of building displacement monitoring method
CN102376089A (en) Target correction method and system
CN108429908B (en) Camera module testing method, device, equipment and medium
CN112254663B (en) Plane deformation monitoring and measuring method and system based on image recognition
CN104964673A (en) Close-shot photography measurement system capable of realizing positioning and attitude determination and close-shot photography measurement method capable of realizing positioning and attitude determination
CN106705962B (en) A kind of method and system obtaining navigation data
CN110146030A (en) Side slope surface DEFORMATION MONITORING SYSTEM and method based on gridiron pattern notation
CN107421473A (en) The two beam laser coaxial degree detection methods based on image procossing
CN105488808A (en) Device and method of field measurement of folding pipes on the basis of three-dimensional space shooting positioning technology
CN115143887B (en) Method for correcting measurement result of visual monitoring equipment and visual monitoring system
CN108958256A (en) A kind of vision navigation method of mobile robot based on SSD object detection model
CN110398199A (en) A kind of track clearance detection method
CN203100724U (en) Device for ranging by utilizing photo taken by camera
CN102223483B (en) Method and device for adjusting position of road traffic monitoring equipment
Huang et al. Measurement method and recent progress of vision-based deflection measurement of bridges: A technical review
CN113658279B (en) Camera internal reference and external reference estimation method, device, computer equipment and storage medium
CN113763463B (en) Method for determining position of acquisition equipment based on image data processing
CN113781546B (en) Method for determining position of image acquisition equipment based on tower drum outer dimension

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
TR01 Transfer of patent right

Effective date of registration: 20240823

Address after: 226300 South facing first office building on the west side of the fourth floor of Building 10, Jufeng Science and Technology Innovation Industrial Park, Xingyuan West Road South and Jinding Road West, Nantong High tech Industrial Development Zone, Jiangsu Province, China

Patentee after: Fengmai Changzhong (Nantong) New Energy Co.,Ltd.

Country or region after: China

Address before: 430000 999 new high road, East Lake Development Zone, Wuhan, Hubei

Patentee before: WINDMAGICS (WUHAN) CO.,LTD.

Country or region before: China

TR01 Transfer of patent right