CN106485752B - The satellite image picture combined based on forward and reverse ground coordinate analogy method - Google Patents
The satellite image picture combined based on forward and reverse ground coordinate analogy method Download PDFInfo
- Publication number
- CN106485752B CN106485752B CN201610569289.8A CN201610569289A CN106485752B CN 106485752 B CN106485752 B CN 106485752B CN 201610569289 A CN201610569289 A CN 201610569289A CN 106485752 B CN106485752 B CN 106485752B
- Authority
- CN
- China
- Prior art keywords
- coordinate
- image
- picpointed
- elements
- forward sight
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
Abstract
A kind of the invention discloses satellite image picture combined based on forward and reverse coordinate analogy method leads to the technical issues of the practicability is poor because being restricted by data source for solving existing satellite image picture coordinate analogy method.Technical solution is to establish scan coordinate system, generates forward sight and backsight image picpointed coordinate, and the elements of exterior orientation at photography moment is calculated according to the sample frequency of image when photography；After picpointed coordinate as unit of pixel is converted to the coordinate in photo coordinate system, topocentric coordinates are calculated using forward intersection formula；According to calculated ground coordinate and internal and external orientation, the picpointed coordinate on collinearity equation calculating forward sight, backsight image is utilized；Satellite image picture point and topocentric coordinates simulation are realized by loop iteration.The present invention is supported without digital elevation model, by way of space intersection and iterative calculation picpointed coordinate, solves the problems, such as that positive simulation coordinate value on vertical flight direction has differences, practical.
Description
Technical field
A kind of the present invention relates to satellite image picture coordinate analogy method is based on forward and reverse phase more particularly to one kind
In conjunction with satellite image picture ground coordinate analogy method.
Background technique
The simulation of optical survey satellite picture point metric data, it is triangular right exactly to establish picture point, projection centre and ground point
It should be related to, be to carry out the master data of subsequent Photogrammetric Processing, and objective evaluation satellite image positioning accuracy has in advance
Effect support.Satellite develop initial stage, meet its state of the art real satellite image data lack, be unable to fully proof theory and
The correctness of algorithm.Generally mould can be strictly imaged according to satellite optical system in satellite image picture point and ground point digital simulation
Type carries out, and main method includes being fitted to calculate picture point by topocentric coordinates by the direct die that picpointed coordinate calculates topocentric coordinates
The inverse analog of coordinate.
" three_line array CCD image satellite photogrammetry principle, Mapping Press, 2006, p14 " discloses positive simulation to document
With the specific steps of inverse analog, but the problem of two methods are not individually present, is analysed in depth.Forward direction simulation basis
The principle of two light forward intersection, directly calculating topocentric coordinates.This method is simple and easy, and iteration is not necessarily in calculating.But by
The influence of elements of exterior orientation measurement accuracy and satellite attitude stabilization degree takes the photograph station and the right side by a left side and takes the photograph on the basis of standing and calculate separately ground point
After coordinate, there are larger differences for coordinate value on vertical flight direction.It has been investigated that the difference magnitude is suitable, symbol on the contrary,
It can be offset after taking number.Three-dimensional intersection analyzing and positioning precision is not influenced using this method, but it is used for bundle adjustment, monolithic
When positioning, theoretical validation and precision analysis result will affect.Inverse analog calculates picture point by topocentric coordinates loop iteration
Coordinate, can avoid the problem of positive simulation has differences, but precondition be must be based on existing digit elevation model data, and
Digital elevation model range is almost the same with elements of exterior orientation range, is restricted by data source.
Summary of the invention
In order to which with overcoming existing satellite image picture coordinate analogy method leads to because being restricted by data source that the practicability is poor not
Foot, a kind of present invention satellite image picture combined based on forward and reverse is provided coordinate analogy method.This method establish with
The image upper left corner is the scan coordinate system of coordinate origin, generates the picpointed coordinate of forward sight and backsight image as unit of pixel, and
The elements of exterior orientation at photography moment is calculated according to the sample frequency of image when photography；By the picpointed coordinate conversion as unit of pixel
After the coordinate in photo coordinate system, topocentric coordinates are calculated using forward intersection formula；It is sat according to calculated ground
Mark and internal and external orientation utilize the picpointed coordinate on collinearity equation calculating forward sight, backsight image；Nothing is realized by loop iteration
Under conditions of digital elevation model, satellite image picture point and topocentric coordinates simulation.The characteristics of due to according to optical satellite image,
In conjunction with existing picture point and topocentric coordinates analogy method, using the method combined based on forward and reverse to satellite image picture
Point and topocentric coordinates are simulated, and in the case where supporting without digital elevation model, pass through space intersection and iterative calculation picture
The mode of point coordinate solves the problems, such as that positive simulation coordinate value on vertical flight direction has differences, practical.
The technical solution adopted by the present invention to solve the technical problems is：A kind of satellite combined based on forward and reverse
Image picture ground coordinate analogy method, its main feature is that including the following steps：
Step 1: generating forward sight image picpointed coordinate in scan coordinate systemWith backsight image picpointed coordinateCoordinate unit is pixel, and origin is located at the upper left corner of forward sight, backsight image.T is consistent with heading, and y hangs down
Directly in heading.
Step 2: being utilized according to the sample frequency of image when photographyWithIt is corresponding to calculate the forward sight camera shooting moment
Elements of exterior orientationAnd rear view camera photography moment corresponding elements of exterior orientation
Step 3: utilizing the picpointed coordinate generated in step 1WithCalculate the coordinate in corresponding photo coordinate systemWithCalculation formula is：
In formula, yol,yorRespectively to the pixel number of principal point, pixel is for forward sight and rear view camera CCD linear array upper extreme point
Pixel size.
Step 4: being calculated according to picpointed coordinate, elements of interior orientation and elements of exterior orientation using space intersection formula
To coordinate (X of the ground point in earth coordinatesA,YA,ZA)；
Step 5: being calculated according to topocentric coordinates, elements of interior orientation and elements of exterior orientation data using collinearity condition equation
Coordinate of the new forward sight, backsight picture point in photo coordinate system
The beneficial effects of the invention are as follows：This method is established using the image upper left corner as the scan coordinate system of coordinate origin, is generated
The picpointed coordinate of forward sight and backsight image as unit of pixel, and the photography moment is calculated according to the sample frequency of image when photography
Elements of exterior orientation；After picpointed coordinate as unit of pixel is converted to the coordinate in photo coordinate system, forward intersection is utilized
Formula calculates topocentric coordinates；According to calculated ground coordinate and internal and external orientation, using collinearity equation calculate forward sight,
Picpointed coordinate on backsight image；Under conditions of realizing nil elevation model by loop iteration, satellite image picture point and ground
The simulation of millet cake coordinate.The characteristics of due to according to optical satellite image, in conjunction with existing picture point and topocentric coordinates analogy method, benefit
Satellite image picture point and topocentric coordinates are simulated with the method combined based on forward and reverse, are being not necessarily to digital elevation
Under model supports, by way of space intersection and iterative calculation picpointed coordinate, solves positive simulation in vertical flight
The problem of coordinate value has differences on direction, it is practical.
It elaborates With reference to embodiment to the present invention.
Specific embodiment
The satellite image picture combined the present invention is based on forward and reverse coordinate analogy method specific step is as follows：
No.1 satellite technology state parameter is drawn with day, carries out digital simulation.In example data, forward sight camera focus fl'=
717.0mm, rear view camera focal length fr'=717.0mm, camera intersection angle αl=25 °, αr=-25 ° of high 500km of satellite orbit, as
First size pixel=0.0065mm, 12000 pixel of linear array length, the satellite flight time about 100 seconds.
Step 1: generating forward sight image picpointed coordinate in scan coordinate systemWith backsight image picpointed coordinateCoordinate unit is pixel, and origin is located at the upper left corner of forward sight, backsight image.T is consistent with heading, and y hangs down
Directly in heading.
Step 2: being utilized according to the sample frequency of image when photographyWithIt is corresponding to calculate the forward sight camera shooting moment
Elements of exterior orientationAnd rear view camera photography moment corresponding elements of exterior orientation
Step 3: utilizing the picpointed coordinate generated in step 1WithCalculate the coordinate in corresponding photo coordinate systemWithCalculation formula is：
In formula, yol,yorRespectively to the pixel number of principal point, pixel is for forward sight and rear view camera CCD linear array upper extreme point
Pixel size.
Step 4: being calculated according to picpointed coordinate, elements of interior orientation and elements of exterior orientation using space intersection formula
To coordinate (X of the ground point in earth coordinatesA,YA,ZA)；
Step 5: being calculated according to topocentric coordinates, elements of interior orientation and elements of exterior orientation data using collinearity condition equation
Coordinate of the new forward sight, backsight picture point in photo coordinate system
Claims (1)
1. a kind of satellite image picture combined based on forward and reverse ground coordinate analogy method, it is characterised in that including following step
Suddenly：
Step 1: generating forward sight image picpointed coordinate in scan coordinate systemWith backsight image picpointed coordinateCoordinate unit is pixel, and origin is located at the upper left corner of forward sight, backsight image；T is consistent with heading, and y hangs down
Directly in heading；
Step 2: being utilized according to the sample frequency of image when photographyWithCalculate forward sight camera shooting moment corresponding foreign side
Bit elementAnd rear view camera photography moment corresponding elements of exterior orientation
Step 3: utilizing the picpointed coordinate generated in step 1WithCalculate the coordinate in corresponding photo coordinate systemWithCalculation formula is：
In formula, yol,yorRespectively to the pixel number of principal point, pixel is that pixel is big for forward sight and rear view camera CCD linear array upper extreme point
It is small；
Step 4: being calculated ground using space intersection formula according to picpointed coordinate, elements of interior orientation and elements of exterior orientation
Coordinate (X of the millet cake in earth coordinatesA,YA,ZA)；
Step 5: being calculated newly according to topocentric coordinates, elements of interior orientation and elements of exterior orientation data using collinearity condition equation
The coordinate of forward sight, backsight picture point in photo coordinate system
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610569289.8A CN106485752B (en) | 2016-07-19 | 2016-07-19 | The satellite image picture combined based on forward and reverse ground coordinate analogy method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610569289.8A CN106485752B (en) | 2016-07-19 | 2016-07-19 | The satellite image picture combined based on forward and reverse ground coordinate analogy method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106485752A CN106485752A (en) | 2017-03-08 |
| CN106485752B true CN106485752B (en) | 2018-11-27 |
Family
ID=58273125
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610569289.8A Active CN106485752B (en) | 2016-07-19 | 2016-07-19 | The satellite image picture combined based on forward and reverse ground coordinate analogy method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106485752B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101509784A (en) * | 2009-03-20 | 2009-08-19 | 西安煤航信息产业有限公司 | GPS//INS data direct directing precision assessment method |
| CN103604417A (en) * | 2013-11-15 | 2014-02-26 | 南京师范大学 | Multi-view image bidirectional matching strategy with constrained object information |
| CN104197898A (en) * | 2014-09-15 | 2014-12-10 | 中国人民解放军总参谋部测绘研究所 | Epipolar ray image generating method of linear array satellite remote sensing image based on projection track method |
-
2016
- 2016-07-19 CN CN201610569289.8A patent/CN106485752B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101509784A (en) * | 2009-03-20 | 2009-08-19 | 西安煤航信息产业有限公司 | GPS//INS data direct directing precision assessment method |
| CN103604417A (en) * | 2013-11-15 | 2014-02-26 | 南京师范大学 | Multi-view image bidirectional matching strategy with constrained object information |
| CN104197898A (en) * | 2014-09-15 | 2014-12-10 | 中国人民解放军总参谋部测绘研究所 | Epipolar ray image generating method of linear array satellite remote sensing image based on projection track method |
Non-Patent Citations (4)
| Title |
|---|
| Accuracy Evaluation of Ground Points from IKONOS High-Resolution Satellite imagery;Guoqlng Zhou 等;《PHOTOGRAMMETRIC ENGINEERING & REMOTE SENSING》;20000930;第66卷(第9期);1103-1112 * |
| POS数据辅助的航空影像变化检测方法研究;袁修孝 等;《武汉大学学报· 信息科学版》;20070430;第32卷(第4期);283-286 * |
| 三线阵CCD影像直接前方交会精度估算;王建荣 等;《测绘科学》;20090731;第34卷(第7期);9-10,17 * |
| 天绘一号无地面控制点摄影测量;王任享 等;《测绘学报》;20130228;第42卷(第1期);1-5 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106485752A (en) | 2017-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110057295B (en) | Monocular vision plane distance measuring method without image control | |
| CN108510551B (en) | Method and system for calibrating camera parameters under long-distance large-field-of-view condition | |
| CN105716542B (en) | A kind of three-dimensional data joining method based on flexible characteristic point | |
| CN110675450B (en) | Method and system for generating orthoimage in real time based on SLAM technology | |
| CN104240262B (en) | Calibration device and calibration method for outer parameters of camera for photogrammetry | |
| CN109919911B (en) | Mobile three-dimensional reconstruction method based on multi-view photometric stereo | |
| CN102661717A (en) | Monocular vision measuring method for iron tower | |
| CN105627991A (en) | Real-time panoramic stitching method and system for unmanned aerial vehicle images | |
| CN106780321B (en) | CBERS-02 satellite HR sensor image overall tight orientation and correction splicing method | |
| KR101308744B1 (en) | System for drawing digital map | |
| CN104677277B (en) | A kind of method and system for measuring object geometric attribute or distance | |
| CN105931200A (en) | Quick geometric precise correction method for small area array spaceborne TDI CCD camera | |
| CN106885585B (en) | Integrated calibration method of satellite-borne photogrammetry system based on light beam adjustment | |
| CN108594255B (en) | Laser ranging auxiliary optical image joint adjustment method and system | |
| Zhao et al. | Direct georeferencing of oblique and vertical imagery in different coordinate systems | |
| CN108154535B (en) | Camera calibration method based on collimator | |
| CN110986888A (en) | Aerial photography integrated method | |
| CN107798668B (en) | Unmanned aerial vehicle imaging hyperspectral geometric correction method and system based on RGB images | |
| WO2024093635A1 (en) | Camera pose estimation method and apparatus, and computer-readable storage medium | |
| CN107063191B (en) | A kind of method of photogrammetric regional network entirety relative orientation | |
| Eiken et al. | Photogrammetric methods applied to Svalbard glaciers: accuracies and challenges | |
| CN105571598A (en) | Satellite laser altimeter footprint camera pose measuring method | |
| JP2013160602A (en) | Photographic surveying apparatus | |
| CN106485752B (en) | The satellite image picture combined based on forward and reverse ground coordinate analogy method | |
| JP2016223934A (en) | Position correction system, position correcting method, and position correction program |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 |