CN101545775A - Method for calculating orientation elements of photo and the height of building by utilizing digital map - Google Patents

Method for calculating orientation elements of photo and the height of building by utilizing digital map Download PDF

Info

Publication number
CN101545775A
CN101545775A CN200910026606A CN200910026606A CN101545775A CN 101545775 A CN101545775 A CN 101545775A CN 200910026606 A CN200910026606 A CN 200910026606A CN 200910026606 A CN200910026606 A CN 200910026606A CN 101545775 A CN101545775 A CN 101545775A
Authority
CN
China
Prior art keywords
photo
msub
coordinate system
point
building
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
CN200910026606A
Other languages
Chinese (zh)
Other versions
CN101545775B (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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN2009100266061A priority Critical patent/CN101545775B/en
Publication of CN101545775A publication Critical patent/CN101545775A/en
Application granted granted Critical
Publication of CN101545775B publication Critical patent/CN101545775B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Processing Or Creating Images (AREA)
  • Image Processing (AREA)

Abstract

The invention provides a method for calculating photo orientation elements and building height by utilizing a digital map. The method comprises the following steps: step S1: a digital photo containing a whole building is shot, and a digital map of a building region is obtained, wherein the digital photo is provided with a photo two-dimensional coordinate system, and the digital map is provided with a digital map coordinate system; step S2: one point on the digital map corresponds to a plurality of points on a wall edge of the building in the digital photo, and altitude data is assumed for wall corner points with different heights according to a stratified structure of the building; the wall corner points are photo control points, and coordinates of the photo control points in the photo two-dimensional coordinate system and the digital map coordinate system are respectively obtained; coordinates of a photo principal point in the photo two-dimensional coordinate system are obtained, and the photo principal point is a perpendicular foot point of a shooting lens center on the photo. The method has small field operation workload, low expense and high precision.

Description

Method for calculating photo azimuth element and building height by using digital map
Technical Field
The invention relates to a method for calculating azimuth elements of a photo and the height of a building, in particular to a method for simultaneously calculating the azimuth elements of a close shot and the height of the building by using a digital map. The technology belongs to the field of photogrammetry.
Background
With the increasing and widespread resolution of the photographic equipment, the common digital photos also become important tools for people to carry out photogrammetry. The primary content of the urban three-dimensional modeling technology lies in the acquisition of a large amount of urban basic three-dimensional data, wherein the acquired data comprises the plane position, the ground elevation, the height and the ground texture of a ground object target. The current three-dimensional data acquisition method is to use photogrammetry and remote sensing technology to obtain a Digital ground model of a city, a Digital Line Graphic (DLG), and a Digital camera to obtain texture of ground objects and height data of buildings on the spot to establish a Digital city model.
A digital city three-dimensional geographic space frame (Li Ming, Wang Relay week, Marek. science publishing Co., Ltd., 2008) provides a method for three-dimensional measurement of a building by using a single photo of a measuring camera, and the theory is that orientation elements of the photo are calculated according to three groups of parallel lines of the building, and on the basis, the single photo three-dimensional modeling of the building is realized. However, many buildings do not provide complete parallel lines in three directions and can only be performed with metrology cameras, making this approach very limiting.
Disclosure of Invention
The technical problem is as follows: the invention aims to solve the technical problem of providing a method for calculating azimuth elements and building height by using a digital map, which can calculate the building height from a single digital photo and the digital map.
The technical scheme is as follows: in order to solve the technical problems, the invention adopts the technical scheme that:
step S1: shooting a digital photo of a complete building to obtain a digital map of the area of the building, wherein the digital photo is provided with a photo two-dimensional coordinate system, and the digital map is provided with a digital map coordinate system;
step S2: one point on the digital map corresponds to a plurality of points on the edge of a building wall in the digital photo, and according to the layered structure of the building, elevation data is assumed for wall corner points with different heights, the wall corner points are photo control points, and coordinates of the photo control points in a photo two-dimensional coordinate system and a digital map coordinate system are obtained; acquiring the coordinates of a main point of a photo in a two-dimensional coordinate system of the photo, wherein the main point of the photo is a foot point of the center of a photographic objective lens on the photo;
step S3: establishing a photogrammetric coordinate system D-XDYDZDThe photogrammetry coordinate system takes an approximate photography point position as an origin, a solid east direction as an X axis, a vertical upward direction as a Y axis and a solid south direction as a Z axis, and converts the coordinates of the photo control points in the digital map coordinate system into coordinates in the photogrammetry coordinate system, wherein the conversion process comprises four steps of translation, exchange of the Y axis and the Z axis and plane rotation, and the method for the translation, the exchange of the Y axis and the Z axis comprises the following steps:
Figure A200910026606D00051
the plane rotation method comprises the following steps:
<math> <mrow> <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>X</mi> <mi>t</mi> </msub> <mo>=</mo> <msub> <mi>X</mi> <mi>D</mi> </msub> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Z</mi> <mi>D</mi> </msub> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Y</mi> <mi>t</mi> </msub> <mo>=</mo> <msub> <mi>Y</mi> <mi>D</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Z</mi> <mi>t</mi> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>X</mi> <mi>D</mi> </msub> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Z</mi> <mi>D</mi> </msub> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow></math>
in the formula (I); xD YD ZDTranslating the map coordinate of the representative photo control point to a coordinate with D as an origin;
Xdrawing (A) YDrawing (A) ZDrawing (A)Coordinates of a digital map coordinate system representing photo control points;
Figure A200910026606D0005151132QIETU
map coordinates representing a photogrammetric origin of coordinates D;
Xt Yt Ztcoordinates in a photogrammetric coordinate system representing photographic control points;
α represents a ground azimuth of a photographing direction at the time of photographing;
step S4: calculating orientation elements and building height of the photo according to a collinear equation in photogrammetry; orientation elements of the picture including the location of the photograph (X)S、YS、ZS) Side roll angle of photography
Figure A200910026606D0005151147QIETU
A photographing pitch angle omega and a picture rotation angle kappa, f a focal length of a photographing objective lens and an image coordinate (x) of a picture principal point of the picture0、y0);
Position of photographic point (X)S、YS、ZS) With XOuter straightTo represent; xS、YS、ZSRespectively representing three-dimensional coordinates of the photographing points in the photogrammetric coordinates;
roll angle for photography
Figure A200910026606D0005151147QIETU
The camera pitch angle omega and the image rotation angle kappa in XExternal cornerTo represent;
f is focal length of photographic objective lens, x0、y0Is the image coordinate of the main point of the image, i.e. the point of the center of the objective lens on the image, f, x0、y0With XInner partTo represent;
image distortion coefficient kappa1With XDistortion ofTo represent;
unknown point space coordinate (X)D、YD、ZD) With XUnknown pointTo represent;
the calculation equation is:
Figure A200910026606D00061
in the formula: vx、VyIs the correction number of the error equation;
the classification by unknowns is represented by a matrix:
V=AXouter straight+BYExternal corner+CXInner part+DXDistortion of+EXUnknown point-L (4)
In the formula: v is a correction matrix of the error equation;
l is an error equation constant term matrix;
Xunknown pointThe unknown number of the height of the building is included;
the various unknowns can then be calculated according to the least squares principle.
Has the advantages that: the digital map is used as a condition to integrally calculate the azimuth elements and the building height of the photo, and the digital map is a future development trend for acquiring urban three-dimensional data due to simple equipment, easy operation, small field workload, low cost and high precision, and solves the problem of calculating the building height from a single digital photo and the digital map.
Drawings
FIG. 1 is a flow chart of a method for calculating orientation elements of a photo according to the present invention;
fig. 2 is a digital image of a building. Wherein 1-9 are point numbers (see tables 1 and 2 below)
Detailed Description
The invention is further described below with reference to the accompanying drawings.
The invention can strictly calculate the shooting position, the shooting attitude, the shot focal length and distortion and the building height of the digital photo, can be widely applied to digital city construction and three-dimensional modeling, and belongs to the field of close-range photogrammetry and computer vision.
The invention utilizes the existing large-scale digital map, considers the layered control structure of the building, calculates the orientation elements of the digital photos and the building height, and has the advantages of small field workload, low cost and high precision.
The invention provides a method for strictly calculating the orientation element of a photo and the height of a building according to a collinear equation by utilizing a plane control point provided by a digital map and considering the layered control structure of the building, and the method comprises the following steps as shown in figure 1:
step S11: shooting a digital photo of a complete building to obtain a digital map of the area of the building, wherein the digital photo is provided with a photo two-dimensional coordinate system, and the digital map is provided with a digital map coordinate system;
step S12: one point on the digital map corresponds to a plurality of points on the edge of a building wall in the digital photo, and according to the layered structure of the building, elevation data is assumed for wall corner points with different heights, the wall corner points are photo control points, and coordinates of the photo control points in a photo two-dimensional coordinate system and a digital map coordinate system are obtained;
step S13: establishing a photogrammetric coordinate system D-XDYDZDThe photogrammetry coordinate system takes an approximate photography point position as an origin, a solid east direction as an X axis, a vertical upward direction as a Y axis and a solid south direction as a Z axis, and converts a digital map coordinate system into a photogrammetry coordinate system;
step S14: and calculating the internal and external orientation elements and the building height of the photo by using the coordinates of the photo control points in the photogrammetric coordinate system and the coordinates in the two-dimensional coordinate system of the photo.
In step S12, a photo control point is selected by using the layered structure of the building, and a corner point of the building can select a plurality of photo control points with different heights on the digital photo, and each of the photo control points assumes a height value.
In step S14, the orientation element of the digital photo including the complete building includes the photographing position parameter X of the photographing photoS、YS、ZSPhotographic attitude parametersOmega, kappa, internal orientation elements f, x0、y0(ii) a According to the coordinates of the photo control point in the photogrammetry coordinate system and the coordinates of the photo two-dimensional coordinate system, the collinear equation in photogrammetry is utilized to calculate the orientation element and the photo distortion parameter k of the digital photo1And building height, wherein XS、YS、ZSRespectively representing the three-dimensional coordinates of the photographing points in the photogrammetric coordinates,
Figure A200910026606D00072
omega and kappa respectively represent the roll angle, pitch angle and image rotation angle of the photography, f is the focal length of the objective lens of the photography, x0、y0The image coordinate of the main point of the photo is the point of the center of the objective lens on the photo.
Namely: the first step is as follows: a digital map of the study area is acquired.
The second step is that: the digital photo of the building is shot by a common digital camera. When shooting, the complete building is required to be shot, the high and low sidelines of the building are mainly included, and calculation is facilitated.
The third step: corresponding points on the digital map and digital image are selected. Since one point on the digital map corresponds to a plurality of points on the side line of the building wall, these points are called shot control points by assuming an elevation data for corner points of different heights. And obtaining the coordinates of the photo control points in the digital map coordinate system and the two-dimensional coordinates in the photo coordinate system by programming corresponding software. The photo control points are selected by utilizing the layered structure of the building, and a corner point of the building can select a plurality of photo control points with different heights on the digital photo, and a height value is respectively assumed.
The fourth step: establishing a photogrammetric coordinate system D-XDYDZDA right-handed system having the photographing position as an origin or the approximate photographing position as an origin, the elevation direction as a Y-axis, and the opposite direction of the photographing direction as a Z-axis. The map coordinate system is converted to the photogrammetric coordinate system by programming corresponding software.
The fifth step: and calculating the orientation element and the building height of the photo by using the coordinate of the same name point in photogrammetry and the plane coordinate of the two-dimensional photo.
More specifically, the method comprises the following steps:
(1) a digital map is acquired.
(2) And shooting the digital photo.
(3) And acquiring the map coordinates and the photo coordinates of the photo control points. Since one point on the digital map corresponds to a plurality of points on the side line of the building wall, one elevation data is assumed by giving corner points of different heights.
(4) Establishing a photogrammetric coordinate system D-XDYDZD. And a right-handed system in which the coordinate system uses the approximate photographing position or the photographing position as an origin, the elevation direction is the Y axis, and the reverse direction of the photographing direction is the Z axis. The map coordinate system is converted to the photogrammetric coordinate system by programming corresponding software. Map coordinate system to photogrammetry coordinate system D-XDYDZDThe conversion process comprises four steps of translation, Y-axis and Z-axis exchange and plane rotation, wherein the translation, the Y-axis and the Z-axis exchange are realized by a formula:
Figure A200910026606D00081
the plane rotation formula is:
<math> <mrow> <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>X</mi> <mi>t</mi> </msub> <mo>=</mo> <msub> <mi>X</mi> <mi>D</mi> </msub> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Z</mi> <mi>D</mi> </msub> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Y</mi> <mi>t</mi> </msub> <mo>=</mo> <msub> <mi>Y</mi> <mi>D</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Z</mi> <mi>t</mi> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>X</mi> <mi>D</mi> </msub> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Z</mi> <mi>D</mi> </msub> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow></math>
in the formula (I); xD YD ZDTranslating the map coordinate of the representative photo control point to a coordinate with D as an origin;
Xdrawing (A) YDrawing (A) ZDrawing (A)Coordinates of a digital map coordinate system representing photo control points;
Figure A200910026606D00092
map coordinates representing a photogrammetric origin of coordinates D;
Xt Yt Ztcoordinates in a photogrammetric coordinate system representing photographic control points;
α represents a ground azimuth of a photographing direction at the time of photographing.
(5) And according to the collinear equation, calculating internal and external orientation elements, picture distortion and building height of the picture by using photogrammetric coordinates and picture coordinates of the picture control points. The mathematical model calculated is:
the solution model includes four types of parameters: position, posture, inner orientation element and image distortion coefficient kappa of photographing point1. The solved unknowns include five types of unknowns, which are:
position X of photographic pointS、YS、ZSWith XOuter straightTo represent; xS、YS、ZSRespectively, three-dimensional coordinates of the photographing points in the photogrammetric coordinates.
Pose of camera spot
Figure A200910026606D00093
ω, κ, by XExternal cornerTo represent;
Figure A200910026606D00094
ω and κ represent the roll angle, pitch angle and image rotation angle of the image.
Internal orientation elements f, x0、y0With XInner partTo represent; f is focal length of photographic objective lens, x0、y0The image coordinates of the main point of the photo, namely the drop foot point of the center of the objective lens on the photo are obtained.
Image distortion coefficient kappa1With XDistortion ofTo represent;
unknown point spatial coordinate XD、YD、ZDIn (1) or (2) with XUnknown pointTo indicate.
The error equation of the solution model is:
Figure A200910026606D00095
in the formula: vx、VyIs the correction of the error equation.
The classification by unknowns is represented by a matrix:
V=AXouter straight+BXExternal corner+CXInner part+DXDistortion of+EXUnknown point-L (4)
In the formula: v is a correction matrix of the error equation;
l is an error equation constant term matrix;
Xunknown pointIncluding the building height unknown.
The various unknowns can then be calculated according to the least squares principle.
And (3) programming by adopting Visual C #. NET, and automatically measuring the coordinates of the digital map and the digital photo by utilizing an ArcGIS ENGINE map control.
The coordinates of three room corner points corresponding to 1-9 points in the photo in fig. 2 are shown in table 1.
TABLE 1 two-dimensional plane coordinates on digital map and elevation initial values thereof
Point number x (Pixel) y (Pixel) Map coordinate X (Rice) Map coordinate Y (Rice) Initial value of elevation Z (meter)
1 -1221.2 -338.5 1061.15 1262.77 10
2 -307.6 -273.6 1061.15 1245.10 10
3 1161.4 -292.6 1115.88 1245.10 10
4 -1214.2 119.6 1061.15 1262.77 18
5 -323 360.9 1061.15 1245.10 18
6 1138.5 -39.4 1115.88 1245.10 18
7 -1199.8 559.8 1061.15 1262.77 26
8 -332.5 945.1 1061.15 1245.10 26
9 1119.5 208.4 1115.88 1245.10 26
The coordinates of the established photogrammetric coordinate origin D are as follows:
X S 0 = 1038.892 ( m ) Y S 0 = 1222.861 ( m ) Z S 0 = 9.000 ( m )
wherein,
Figure A200910026606D00102
are approximations. Approximate azimuth angle in photography: alpha is alphaTaking a photographThe inverse of the curve 1.06730270945475 is the Z-axis. Conversion of the map coordinate system in Table 1 to the photogrammetric coordinate system D-XDYDZDThe conversion results are shown in Table 2.
TABLE 2 photogrammetry coordinate system coordinates of control points
Point number x (Pixel) y (Pixel) Map coordinate X (Rice) Map coordinate Y (Rice) Initial value of elevation Z (meter) XD(Rice) YD(Rice) ZD(Rice)
1 -1221.2 -338.5 1061.15 1262.77 10 -24.22 1 -38.75
2 -307.6 -273.6 1061.15 1245.10 10 -8.74 1 -30.23
3 1161.4 -292.6 1115.88 1245.10 10 17.56 1 -78.23
4 -1214.2 119.6 1061.15 1262.77 18 -24.22 9 -38.75
5 -323 360.9 1061.15 1245.10 18 -8.74 9 -30.23
6 1138.5 -39.4 1115.88 1245.10 18 17.56 9 -78.23
7 -1199.8 559.8 1061.15 1262.77 26 -24.22 17 -38.75
8 -332.5 945.1 1061.15 1245.10 26 -8.74 17 -30.23
9 1119.5 208.4 1115.88 1245.10 26 17.56 17 -78.23
The unknowns calculated using the transformed photogrammetric coordinate system coordinates and the photograph coordinates are shown in table 3.
TABLE 3 results of calculation
Calculating parameters Calculation results
Xs(m) -1.483
Ys(m) -0.396
Zs(m) -1.784
Phi (arc) -0.1361996
Omega (radian) 0.1306600
K (arc) -0.039236
f (Pixel) 2901.8
x0 (Pixel) 13.0
y0 (Pixel) -20.6
K1 1.94×10-8
Height (m) 13.13

Claims (1)

1. A method for calculating photo orientation elements and building height using a digital map, comprising: the method comprises the following steps:
step 1: shooting a digital photo of a complete building to obtain a digital map of the area of the building, wherein the digital photo is provided with a photo two-dimensional coordinate system, and the digital map is provided with a digital map coordinate system;
step 2: one point on the digital map corresponds to a plurality of points on the edge of a building wall in the digital photo, and according to the layered structure of the building, elevation data is assumed for wall corner points with different heights, the wall corner points are photo control points, and coordinates of the photo control points in a photo two-dimensional coordinate system and a digital map coordinate system are obtained; acquiring the coordinates of a main point of a photo in a two-dimensional coordinate system of the photo, wherein the main point of the photo is a foot point of the center of a photographic objective lens on the photo;
and step 3: establishing a photogrammetric coordinate system D-XDYDZDThe photogrammetry coordinate system takes an approximate photography point position as an origin, a solid east direction as an X axis, a vertical upward direction as a Y axis and a solid south direction as a Z axis, and converts the coordinates of the photo control points in the digital map coordinate system into coordinates in the photogrammetry coordinate system, wherein the conversion process comprises four steps of translation, exchange of the Y axis and the Z axis and plane rotation, and the method for the translation, the exchange of the Y axis and the Z axis comprises the following steps:
Figure A200910026606C00021
formula 1
The plane rotation method comprises the following steps:
<math> <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>X</mi> <mi>t</mi> </msub> <mo>=</mo> <msub> <mi>X</mi> <mi>D</mi> </msub> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Z</mi> <mi>D</mi> </msub> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Y</mi> <mi>t</mi> </msub> <mo>=</mo> <msub> <mi>Y</mi> <mi>D</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Z</mi> <mi>t</mi> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>X</mi> <mi>D</mi> </msub> <mi>sin</mi> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>Z</mi> <mi>D</mi> </msub> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow></math> formula 2
In the formula (I); xDYDZDTranslating the map coordinate of the representative photo control point to a coordinate with D as an origin;
Xdrawing (A)YDrawing (A)ZDrawing (A)Coordinates of a digital map coordinate system representing photo control points;
Figure A200910026606C00023
map coordinates representing a photogrammetric origin of coordinates D;
XtYtZtcoordinates in a photogrammetric coordinate system representing photographic control points;
α represents a ground azimuth of a photographing direction at the time of photographing;
and 4, step 4: calculating orientation elements and building height of the photo according to a collinear equation in photogrammetry; orientation elements of the picture including the location of the photograph (X)S、YS、ZS) Side roll angle of photography
Figure A200910026606C00031
The shooting pitch angle omega and the image rotation angle kappa, f the focal length of the shooting objective lens and the image coordinate (x) of the image principal point of the image0、y0);
Position of photographic point (X)S、YS、ZS) With XOuter straightTo represent; xS、YS、ZSRespectively representing three-dimensional coordinates of the photographing points in the photogrammetric coordinates;
roll angle for photography
Figure A200910026606C00032
Photographing pitch angle omega and picture rotation angle kappa in XExternal cornerTo represent;
f is focal length of photographic objective lens, x0、y0Is the image coordinate of the main point of the image, i.e. the point of the center of the objective lens on the image, f, x0、y0With XInner partTo represent;
image distortion coefficient kappa1With XDistortion ofTo represent;
unknown point space coordinate (X)D、YD、ZD) With XUnknown pointTo represent;
the calculation equation is:
Figure A200910026606C00033
in the formula: vx、VyIs the correction number of the error equation;
the classification by unknowns is represented by a matrix:
V=AXouter straight+BXExternal corner+CXInner part+DXDistortion of+EXUnknown point-L formula 4
In the formula: v is a correction matrix of the error equation;
l is an error equation constant term matrix;
Xunknown pointThe unknown number of the height of the building is included;
the various unknowns can then be calculated according to the least squares principle.
CN2009100266061A 2009-05-05 2009-05-05 Method for calculating orientation elements of photo and the height of building by utilizing digital map Expired - Fee Related CN101545775B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100266061A CN101545775B (en) 2009-05-05 2009-05-05 Method for calculating orientation elements of photo and the height of building by utilizing digital map

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100266061A CN101545775B (en) 2009-05-05 2009-05-05 Method for calculating orientation elements of photo and the height of building by utilizing digital map

Publications (2)

Publication Number Publication Date
CN101545775A true CN101545775A (en) 2009-09-30
CN101545775B CN101545775B (en) 2011-06-08

Family

ID=41193045

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100266061A Expired - Fee Related CN101545775B (en) 2009-05-05 2009-05-05 Method for calculating orientation elements of photo and the height of building by utilizing digital map

Country Status (1)

Country Link
CN (1) CN101545775B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101825452A (en) * 2010-04-20 2010-09-08 中国第一冶金建设有限责任公司 Method for performing inter-conversion between survey coordinates and construction coordinates
CN103471583A (en) * 2013-09-09 2013-12-25 西南交通大学 Method for determining evacuating route of persons in tunnel based on video image
CN106323240A (en) * 2015-06-27 2017-01-11 腾讯科技(深圳)有限公司 Image object measuring method and apparatus
CN106871864A (en) * 2017-02-13 2017-06-20 中科宇图科技股份有限公司 A kind of method that depth of building is automatically extracted based on three-dimensional satellite image
CN108067939A (en) * 2016-11-18 2018-05-25 成都飞机工业(集团)有限责任公司 A kind of point position in space Measuring datum error compensation method
CN105279372B (en) * 2015-09-29 2018-12-14 百度在线网络技术(北京)有限公司 A kind of method and apparatus of determining depth of building
US11195324B1 (en) 2018-08-14 2021-12-07 Certainteed Llc Systems and methods for visualization of building structures
CN115359114A (en) * 2022-08-16 2022-11-18 中建一局集团第五建筑有限公司 Positioning method, positioning device, electronic equipment and computer readable storage medium

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1664884A (en) * 2005-03-02 2005-09-07 西安四维航测遥感中心 Method for making metropolitan area three-dimensional visible and measurable hypsometric maps
CN100533486C (en) * 2007-08-07 2009-08-26 深圳先进技术研究院 Digital city full-automatic generating method
CN100580371C (en) * 2008-05-16 2010-01-13 南京师范大学 Building single image three-dimensional measurement method based on space right-angle relationship
CN101325724B (en) * 2008-07-23 2010-06-02 四川虹微技术有限公司 Method for correcting polar line of stereoscopic picture pair

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101825452A (en) * 2010-04-20 2010-09-08 中国第一冶金建设有限责任公司 Method for performing inter-conversion between survey coordinates and construction coordinates
CN103471583A (en) * 2013-09-09 2013-12-25 西南交通大学 Method for determining evacuating route of persons in tunnel based on video image
CN103471583B (en) * 2013-09-09 2015-10-28 西南交通大学 A kind of tunnel evacuating personnel gesture determination method based on video image
CN106323240A (en) * 2015-06-27 2017-01-11 腾讯科技(深圳)有限公司 Image object measuring method and apparatus
CN106323240B (en) * 2015-06-27 2019-03-01 腾讯科技(深圳)有限公司 A kind of image object measurement method and device
CN105279372B (en) * 2015-09-29 2018-12-14 百度在线网络技术(北京)有限公司 A kind of method and apparatus of determining depth of building
CN108067939A (en) * 2016-11-18 2018-05-25 成都飞机工业(集团)有限责任公司 A kind of point position in space Measuring datum error compensation method
CN108067939B (en) * 2016-11-18 2021-08-03 成都飞机工业(集团)有限责任公司 Space point location measurement reference error compensation method
CN106871864A (en) * 2017-02-13 2017-06-20 中科宇图科技股份有限公司 A kind of method that depth of building is automatically extracted based on three-dimensional satellite image
US11195324B1 (en) 2018-08-14 2021-12-07 Certainteed Llc Systems and methods for visualization of building structures
US11704866B2 (en) 2018-08-14 2023-07-18 Certainteed Llc Systems and methods for visualization of building structures
CN115359114A (en) * 2022-08-16 2022-11-18 中建一局集团第五建筑有限公司 Positioning method, positioning device, electronic equipment and computer readable storage medium

Also Published As

Publication number Publication date
CN101545775B (en) 2011-06-08

Similar Documents

Publication Publication Date Title
CN107492069B (en) Image fusion method based on multi-lens sensor
CN101545775A (en) Method for calculating orientation elements of photo and the height of building by utilizing digital map
CN109115186B (en) 360-degree measurable panoramic image generation method for vehicle-mounted mobile measurement system
AU2011312140C1 (en) Rapid 3D modeling
US7944547B2 (en) Method and system of generating 3D images with airborne oblique/vertical imagery, GPS/IMU data, and LIDAR elevation data
CN104330074B (en) Intelligent surveying and mapping platform and realizing method thereof
CN106228538B (en) Binocular vision indoor orientation method based on logo
WO2021004416A1 (en) Method and apparatus for establishing beacon map on basis of visual beacons
CN108168521A (en) One kind realizes landscape three-dimensional visualization method based on unmanned plane
CN101545776B (en) Method for obtaining digital photo orientation elements based on digital map
US10467726B2 (en) Post capture imagery processing and deployment systems
CN101763643A (en) Automatic calibration method for structured light three-dimensional scanner system
CN109272574A (en) Linear array rotary scanning camera imaging model building method and scaling method based on projective transformation
CN110782498B (en) Rapid universal calibration method for visual sensing network
CN109900205A (en) A kind of quick calibrating method of high-precision single line laser device and optical camera
CN105374067A (en) Three-dimensional reconstruction method based on PAL cameras and reconstruction system thereof
CN104655106B (en) Autonomous positioning based on GPS RTK and full-view image orients plotting method
CN104732557A (en) Color point cloud generating method of ground laser scanner
CN117274499B (en) Unmanned aerial vehicle oblique photography-based steel structure processing and mounting method
CN106023207B (en) It is a kind of to be enjoyed a double blessing the Municipal Component acquisition method of scape based on traverse measurement system
CN115511956A (en) Unmanned aerial vehicle imaging positioning method
Liu et al. A new approach to fast mosaic UAV images
CN115423863A (en) Camera pose estimation method and device and computer readable storage medium
CN111724446A (en) Zoom camera external parameter calibration method for building three-dimensional reconstruction
CN112857328B (en) Calibration-free photogrammetry method

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
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110608

Termination date: 20140505