CN109003295A - A kind of unmanned plane aviation image fast matching method - Google Patents

A kind of unmanned plane aviation image fast matching method Download PDF

Info

Publication number
CN109003295A
CN109003295A CN201810811943.0A CN201810811943A CN109003295A CN 109003295 A CN109003295 A CN 109003295A CN 201810811943 A CN201810811943 A CN 201810811943A CN 109003295 A CN109003295 A CN 109003295A
Authority
CN
China
Prior art keywords
unmanned plane
point
image
coordinate system
photo
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
CN201810811943.0A
Other languages
Chinese (zh)
Other versions
CN109003295B (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.)
ZHONGYE SHENKAN ENGINEERING TECHNOLOGY Co Ltd
Original Assignee
ZHONGYE SHENKAN ENGINEERING TECHNOLOGY 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 ZHONGYE SHENKAN ENGINEERING TECHNOLOGY Co Ltd filed Critical ZHONGYE SHENKAN ENGINEERING TECHNOLOGY Co Ltd
Publication of CN109003295A publication Critical patent/CN109003295A/en
Application granted granted Critical
Publication of CN109003295B publication Critical patent/CN109003295B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • G06T7/73Determining position or orientation of objects or cameras using feature-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30181Earth observation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)

Abstract

The present invention provides a kind of unmanned plane aviation image fast matching method, is related to technical field of mapping.This method comprises: computer simulation generates Three Dimensional Ground sampling point set, and point set is transformed under selected coordinate system;It is exposure point by the processing of unmanned plane POS data, and exposure point bit map is arrived under selected coordinate system, then exposure point is added to sampling point set and is measured under coordinate system;The picture plane coordinates that picture point is calculated using central projection imaging equation, determines that sampled point projects to unmanned plane photo to the visibility of photo, and by point set according to as plane coordinates;The matching relationship between image is determined according to the sampled point shared between photo.A kind of unmanned plane aviation image fast matching method provided by the invention, image projecting based on virtual ground sampled point, utilize the visual ground sampled point set of image, the problem of judging the matching relationship between photo, unmanned plane image error hiding caused by traditional treatment method can be thoroughly eliminated during unmanned plane Image Matching.

Description

A kind of unmanned plane aviation image fast matching method
Technical field
The present invention relates to technical field of mapping, and in particular to a kind of unmanned plane aviation image fast matching method.
Background technique
In photogrammetric data processing, Image Matching provides identical point coordinates observation for image orientation, and precision is very The precision of photogrammetric end result is determined in big degree, and the speed of Image Matching is to influence photogrammetric data processing effect An important factor for rate, traditional aerial survey using manned aircraft as flying platform, good flight control can guarantee lesser scale, Visual angle and corner variation, therefore the point of same place is it is expected that and unmanned aerial vehicle platform is larger by weather and the influence of topography, image Scale, visual angle and corner change greatly, the steady matching of image is difficult to realize using traditional aerial survey image matching method.It is existing Image feature (SIFT) operator the characteristic point of a large amount of scales and invariable rotary can be extracted from image, and characteristic point is to affine Transformation and illumination variation have good adaptability, therefore are widely used in photogrammetric Image Matching, however, based on nearest The Image Matching of neighbour's search usually there is error hiding as a result, the image of error hiding, which is introduced compensating computation, will cause empty three precision Decline, even results in image global orientation and does not restrain.
Summary of the invention
In view of the problems of the existing technology, the present invention provides a kind of unmanned plane aviation image fast matching method, in terms of Based on calculation machine simulates Three Dimensional Ground sampling point set, projected using central projection imaging equation virtual image, according between image Shared sampled point determines the matching relationship between image, eliminates image error hiding problem caused by traditional treatment method.
To achieve the goals above, a kind of unmanned plane aviation image fast matching method, comprising the following steps:
Step 1: according to the landform and range of tested region, generating ground three using computer simulation in different ways Sampling point set is tieed up, it then will be under the coordinate transform of sampling point set to selected measurement coordinate system;Wherein, the area flat for landform Domain carries out three-dimensional coordinate sampling in an elevation plane, for region with a varied topography, using global digital elevation model, Generate three-dimensional sample point set;
Step 2: the position that unmanned plane positioning and orientation system is obtained and attitude data processing transform to expose point Under selected measurement coordinate system, the exposure point of unmanned plane and Three Dimensional Ground sampling point set are added to and measured under coordinate system;
Step 3: picture plane coordinates (x, y) of picture point is calculated using central projection imaging equation, it is flat according to the picture of sampled point Definitely ground sampling point set and is projected to unmanned plane photo to the visibility of photo to areal coordinate by surface sample point;
Step 4: the matching relationship between image being determined according to the sampled point shared between photo, if for any two photographs Piece concludes that there are matching relationships for this two photos if there are intersections for visible ground sampled point.
Further, the formula of central projection imaging equation is as follows in the step 3:
Wherein, f is photogrammetric focal length, and (X, Y, Z) is the measurement coordinate system coordinate of ground sampled point, (XS, YS, ZS) For the measurement coordinate system coordinate of projection centre, a1、a2、a3、b1、b2、b3、c1、c2、c3The cosine in respectively 9 directions;9 directions Cosine calculation formula it is as follows:
Wherein, ω is corner of the photogrammetric image space coordinate system around X-axis,For photogrammetric image space coordinate system Around the corner of Y-axis, κ is the corner of photogrammetric image space coordinate system about the z axis.
Beneficial effects of the present invention:
The present invention proposes a kind of unmanned plane aviation image fast matching method, utilizes computer simulation ground sampled point, knot Unmanned plane positioning and orientation (POS) data are closed, the image projecting of virtual ground sampled point is calculated based on central projection imaging equation, On the basis of this, friendship is asked using the visual ground sampled point set of image, the matching relationship between photo is judged, in unmanned plane image The problem of unmanned plane image error hiding caused by traditional treatment method can be thoroughly eliminated in matching process.
Detailed description of the invention
Fig. 1 is the flow chart of unmanned plane aviation image fast matching method in the embodiment of the present invention;
Fig. 2 is the Three Dimensional Ground sampling point set for utilizing computer simulation to generate in tested region in the embodiment of the present invention;
Fig. 3 is unmanned plane camera exposure point in the embodiment of the present invention;
Fig. 4 is that the exposure point of unmanned plane in the embodiment of the present invention is superimposed point with Three Dimensional Ground sampling point set;
Fig. 5 is the ground sampled point of single photo covering in the embodiment of the present invention;
Fig. 6 is the three-dimensional exposure point of unmanned plane photo in the embodiment of the present invention.
Specific embodiment
It is right in the following with reference to the drawings and specific embodiments in order to be more clear the purpose of the present invention, technical solution and advantage The present invention is described in further details.Described herein specific examples are only used to explain the present invention, is not used to limit this Invention.
A kind of unmanned plane aviation image fast matching method, process is as shown in Figure 1, that the specific method is as follows is described:
Step 1: according to the landform and range of tested region, generating ground three using computer simulation in different ways Sampling point set is tieed up, it then will be under the coordinate transform of sampling point set to selected measurement coordinate system;Wherein, the area flat for landform Domain carries out three-dimensional coordinate sampling in an elevation plane, for region with a varied topography, using global digital elevation model, Generate three-dimensional sample point set.
It is as shown in Figure 2 in the Three Dimensional Ground sampling point set that tested region is generated using computer simulation in the present embodiment.
Step 2: the position that unmanned plane positioning and orientation system is obtained and attitude data processing transform to expose point Under selected measurement coordinate system, the exposure point of unmanned plane and Three Dimensional Ground sampling point set are added to and measured under coordinate system.
In the present embodiment, unmanned plane camera exposure point is as shown in figure 3, the exposure point and Three Dimensional Ground of unmanned plane sample It is as shown in Figure 4 that point set is superimposed point.
Step 3: picture plane coordinates (x, y) of picture point is calculated using central projection imaging equation, it is flat according to the picture of sampled point Definitely ground sampling point set and is projected to unmanned plane photo to the visibility of photo to areal coordinate by surface sample point.
Shown in the formula such as formula (1) of the central projection imaging equation:
Wherein, f is photogrammetric focal length, and (X, Y, Z) is the measurement coordinate system coordinate of ground sampled point, (XS, YS, ZS) For the measurement coordinate system coordinate of projection centre, a1、a2、a3、b1、b2、b3、c1、c2、c3The cosine in respectively 9 directions;9 directions Cosine calculation formula such as formula (2) shown in:
Wherein, ω is corner of the photogrammetric image space coordinate system around X-axis,For photogrammetric image space coordinate system Around the corner of Y-axis, κ is the corner of photogrammetric image space coordinate system about the z axis.
In the present embodiment, the ground sampled point of single photo covering is as shown in figure 5, the ground in image coverage area is adopted Sampling point is visible to the photo.
Step 4: the matching relationship between image being determined according to the sampled point shared between photo, for any two photographs Piece concludes that there are matching relationships for this two photos if there are intersections for visible ground sampled point.
Image Matching operation is also carried out using the image matching method based on SIFT operator in the prior art, SIFT operator is The classic algorithm of Image Matching can extract the characteristic point of a large amount of scales and invariable rotary from image, therefore be widely used in Photogrammetric and 3D computer vision.Neighborhood where SIFT operator can extract position, scale and the characteristic point of characteristic point Feature description.SIFT feature description is with the expression of 128 dimensional vectors, for matching to characteristic point.Ordinary matches method utilizes core The characteristic point of most of error hiding is rejected in line constraint, but cannot thoroughly eliminate the error hiding problem of image.It is based on SIFT below The image matching method of operator and unmanned plane aviation image fast matching method proposed by the present invention compare test.
To 4 course lines, SIFT is matched two-by-two for totally 120 unmanned plane photos progress, and Fig. 6 shows 120 unmanned plane photos The sparse cloud that three-dimensional exposure point and sky three are rebuild.Table 1 is the unmanned plane photo matching result based on SIFT operator, in conjunction with figure 6, it is known that there are error hidings for SIFT matching result.Number is under every photo and No. 1 course line under current course line number in table 1 Photo 1 matches obtained points of the same name, and number, which is greater than 1, indicates two photos there are matching relationship, the matching result underlined For the image pair of error hiding.
Unmanned plane photo matching result of the table 1 based on SIFT operator
When the photo matching relationship being calculated using method proposed by the present invention is as shown in table 2.Matching result in table 2 For the matching relationship of the photo 1 under each photo and No. 1 course line, wherein 1 indicates two photos, and there are matching relationships, 0 to indicate Matching relationship is not present in photo, and the matching result underlined is the image pair of error hiding in table 1.By comparison it can be found that originally The problem of method that invention proposes completely eliminates unmanned plane image error hiding caused by SIFT image matching method.
The matching result of 2 unmanned plane aviation image fast matching method of table
In terms of computational efficiency, match time is 8 minutes to matched 120 photos of SIFT two-by-two, and proposed by the present invention The calculating time of photo matching relationship calculation method is 5 seconds, hence it is evident that is better than SIFT image matching method.
Finally, it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although Present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that;It still may be used To modify to technical solution documented by previous embodiment, or some or all of the technical features are equal Replacement;Thus these are modified or replaceed, defined by the claims in the present invention that it does not separate the essence of the corresponding technical solution Range.

Claims (2)

1. a kind of unmanned plane aviation image fast matching method, which comprises the following steps:
Step 1: according to the landform and range of tested region, generating Three Dimensional Ground using computer simulation in different ways and adopt Sampling point collection, then will be under the coordinate transform of sampling point set to selected measurement coordinate system;Wherein, the region flat for landform, Three-dimensional coordinate sampling is carried out in an elevation plane, and region with a varied topography is generated using global digital elevation model Three-dimensional sample point set;
Step 2: the position that unmanned plane positioning and orientation system is obtained and attitude data processing transform to selected to expose point Measurement coordinate system under, the exposure point of unmanned plane and Three Dimensional Ground sampling point set are added to and measured under coordinate system;
Step 3: calculating picture plane coordinates (x, y) of picture point using central projection imaging equation, sat according to the picture plane of sampled point Definitely ground sampling point set and is projected to unmanned plane photo to the visibility of photo to mark by surface sample point;
Step 4: the matching relationship between image being determined according to the sampled point shared between photo, for any two photos, such as There are intersections for the visible ground sampled point of fruit, that is, conclude that there are matching relationships for this two photos.
2. unmanned plane aviation image fast matching method according to claim 1, which is characterized in that in the step 3 The formula that the heart projects imaging equation is as follows:
Wherein, f is photogrammetric focal length, and (X, Y, Z) is the measurement coordinate system coordinate of ground sampled point, (XS, YS, ZS) it is to throw The measurement coordinate system coordinate at shadow center, a1、a2、a3、b1、b2、b3、c1、c2、c3The cosine in respectively 9 directions;9 directions it is remaining String calculation formula is as follows:
Wherein, ω is corner of the photogrammetric image space coordinate system around X-axis,It is photogrammetric image space coordinate system around Y-axis Corner, κ be the corner of photogrammetric image space coordinate system about the z axis.
CN201810811943.0A 2018-04-11 2018-07-23 Rapid matching method for aerial images of unmanned aerial vehicle Active CN109003295B (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN2018103205317 2018-04-11
CN201810320531 2018-04-11

Publications (2)

Publication Number Publication Date
CN109003295A true CN109003295A (en) 2018-12-14
CN109003295B CN109003295B (en) 2021-07-23

Family

ID=64596156

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810811943.0A Active CN109003295B (en) 2018-04-11 2018-07-23 Rapid matching method for aerial images of unmanned aerial vehicle

Country Status (1)

Country Link
CN (1) CN109003295B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112002007A (en) * 2020-08-31 2020-11-27 胡翰 Model obtaining method and device based on air-ground image, equipment and storage medium
CN113485424A (en) * 2021-07-19 2021-10-08 武汉中测晟图遥感技术有限公司 Design method of pole tower inspection route

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101114022A (en) * 2007-09-04 2008-01-30 国家海洋局第二海洋研究所 Navigation multiple spectrum scanner geometric approximate correction method under non gesture information condition
CN101127078A (en) * 2007-09-13 2008-02-20 北京航空航天大学 Unmanned machine vision image matching method based on ant colony intelligence
CN101464149A (en) * 2008-12-31 2009-06-24 武汉大学 POS auxiliary aviation image matching method
CN102136136A (en) * 2011-03-17 2011-07-27 南京航空航天大学 Luminosity insensitivity stereo matching method based on self-adapting Census conversion
CN102607534A (en) * 2012-03-13 2012-07-25 上海交通大学 Satellite relative attitude measuring method based on structure from motion
US8315794B1 (en) * 2006-09-05 2012-11-20 Honeywell International Inc. Method and system for GPS-denied navigation of unmanned aerial vehicles
CN104766302A (en) * 2015-02-05 2015-07-08 武汉大势智慧科技有限公司 Method and system for optimizing laser scanning point cloud data by means of unmanned aerial vehicle images
CN106940181A (en) * 2017-03-10 2017-07-11 中国电建集团昆明勘测设计研究院有限公司 A kind of unmanned plane image picture control distribution of net is built and the optional commensurate in scope method of aerophotograph

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8315794B1 (en) * 2006-09-05 2012-11-20 Honeywell International Inc. Method and system for GPS-denied navigation of unmanned aerial vehicles
CN101114022A (en) * 2007-09-04 2008-01-30 国家海洋局第二海洋研究所 Navigation multiple spectrum scanner geometric approximate correction method under non gesture information condition
CN101127078A (en) * 2007-09-13 2008-02-20 北京航空航天大学 Unmanned machine vision image matching method based on ant colony intelligence
CN101464149A (en) * 2008-12-31 2009-06-24 武汉大学 POS auxiliary aviation image matching method
CN102136136A (en) * 2011-03-17 2011-07-27 南京航空航天大学 Luminosity insensitivity stereo matching method based on self-adapting Census conversion
CN102607534A (en) * 2012-03-13 2012-07-25 上海交通大学 Satellite relative attitude measuring method based on structure from motion
CN104766302A (en) * 2015-02-05 2015-07-08 武汉大势智慧科技有限公司 Method and system for optimizing laser scanning point cloud data by means of unmanned aerial vehicle images
CN106940181A (en) * 2017-03-10 2017-07-11 中国电建集团昆明勘测设计研究院有限公司 A kind of unmanned plane image picture control distribution of net is built and the optional commensurate in scope method of aerophotograph

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112002007A (en) * 2020-08-31 2020-11-27 胡翰 Model obtaining method and device based on air-ground image, equipment and storage medium
CN112002007B (en) * 2020-08-31 2024-01-19 胡翰 Model acquisition method and device based on air-ground image, equipment and storage medium
CN113485424A (en) * 2021-07-19 2021-10-08 武汉中测晟图遥感技术有限公司 Design method of pole tower inspection route

Also Published As

Publication number Publication date
CN109003295B (en) 2021-07-23

Similar Documents

Publication Publication Date Title
JP7413321B2 (en) Daily scene restoration engine
CN109658461B (en) Unmanned aerial vehicle positioning method based on cooperation two-dimensional code of virtual simulation environment
US9799139B2 (en) Accurate image alignment to a 3D model
US20090154793A1 (en) Digital photogrammetric method and apparatus using intergrated modeling of different types of sensors
CN106780729A (en) A kind of unmanned plane sequential images batch processing three-dimensional rebuilding method
CN110033489A (en) A kind of appraisal procedure, device and the equipment of vehicle location accuracy
CN106485690A (en) Cloud data based on a feature and the autoregistration fusion method of optical image
IL214151A (en) Method and apparatus for three-dimensional image reconstruction
CN107831515B (en) Underwater Navigation method and system
CN102506867B (en) SINS (strap-down inertia navigation system)/SMANS (scene matching auxiliary navigation system) combined navigation method based on Harris comer matching and combined navigation system
CN109708649A (en) A kind of attitude determination method and system of remote sensing satellite
JP2012118666A (en) Three-dimensional map automatic generation device
CN108917753A (en) Method is determined based on the position of aircraft of structure from motion
CN111710040B (en) High-precision map construction method, system, terminal and storage medium
CN109003295A (en) A kind of unmanned plane aviation image fast matching method
Fu-Sheng et al. Batch reconstruction from UAV images with prior information
CN111612829B (en) High-precision map construction method, system, terminal and storage medium
CN107063191B (en) A kind of method of photogrammetric regional network entirety relative orientation
JP6761388B2 (en) Estimator and program
CN113129422A (en) Three-dimensional model construction method and device, storage medium and computer equipment
Kupervasser et al. Robust positioning of drones for land use monitoring in strong terrain relief using vision-based navigation
CN109946682A (en) GF3 data baseline estimation method based on ICESat/GLAS
CN109341685B (en) Fixed wing aircraft vision auxiliary landing navigation method based on homography transformation
Geva et al. Estimating camera pose using bundle adjustment and digital terrain model constraints
CN104331882B (en) Method for measuring speed of aircraft

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