CN102508260A - Geometric imaging construction method for side-looking medium resolution ratio satellite - Google Patents
Geometric imaging construction method for side-looking medium resolution ratio satellite Download PDFInfo
- Publication number
- CN102508260A CN102508260A CN2011103889150A CN201110388915A CN102508260A CN 102508260 A CN102508260 A CN 102508260A CN 2011103889150 A CN2011103889150 A CN 2011103889150A CN 201110388915 A CN201110388915 A CN 201110388915A CN 102508260 A CN102508260 A CN 102508260A
- Authority
- CN
- China
- Prior art keywords
- satellite
- reference mark
- parameter
- model
- equation
- 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
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
The invention relates to a geometric imaging construction method for a side-looking medium resolution ratio satellite. The method comprises the following step: (1) importing satellite image metadata and satellite almanac data to determine a satellite orbit and an attitude fitting equation, parameter initial value of a datalink node and other parameter initial values; (2) automatically primarily acquiring an image control point initial value according to an image metadata containing projection information of a reference; (3) manually setting at least three control points, automatically matching image control points by gray matching, automatically providing a control point missing area, and modifying distribution of control points in the control point missing area; (4) modifying distribution of control points for which the step (3) is completed, transforming coordinates, and constructing a geometric imaging model of the satellite; (5) modifying, iterating and calculating model parameters according to the acquired parameter initial values of step (1) and the geometric imaging model of step (4); and (6) outputting the corrected model result and accuracy index.
Description
Technical field
The invention belongs to the image processing field, relate to a kind of how much imaging construction methods towards side-looking intermediate-resolution satellite.
Background technology
The when and where that disaster takes place always has great randomness; Progress along with modern humans society; Timely the physical environment disaster being responded fast and corresponding assistance is provided is the effective means that reduces casualty loss, thereby the monitoring of regional disaster is become particularly important.The core mode of setting up monitoring system is to set up cover reaction time remote sensing satellite system rapidly; Become today of mainstream development direction at lower-cost moonlet constellation, the intermediate-resolution moonlet constellation that uses several big field angle and be observed imaging mode with side-looking becomes as the information obtain manner of environment mitigation monitoring system and not only guarantees response speed but also cost-effective optimal selection.And remotely-sensed data is supplying need to give certain Geographic Reference information before the use, promptly through specific mode it is carried out geometric manipulations, the foundation of how much imaging models of the core of geometric manipulations.
Be set up in the image that the line array sensor of the intermediate-resolution wide visual field angle side-looking observation on the moonlet collects following some principal character is arranged: sweep limit is big, receives the earth curvature image serious, owing to adopt this influence of side-looking observation more serious; Field angle is big, on single scan line, appears with the central projection mode fully, and differs bigger with parallel projection; Satellite platform stability is general, for lower-cost moonlet, because the various factors satellite attitude jitter is more serious, can't carry out match with parallel projection.And present stage centering low resolution image the geometric manipulations means mainly rely on classical fitting of a polynomial mode, promptly carry out picture point the single mode of picture point carried out registration.But adopt the intermediate-resolution satellite of the big field angle of line array sensor, because area coverage is big, landform and earth curvature with the non-constant of single fitting of a polynomial precision, can not satisfy the drawing demand of most applications to the influence of distortion distortion.Its also inapplicable affined transformation is simulated its imaging process simultaneously, and this application of model condition restriction is neglected the high-resolution satellite of rink corner in long-focus; Because the factor of satellite platform own, the attitude jitter pattern is comparatively complicated again, and satellite to start at data precision limited, can't use direct linear transformation's model (DLT) or rational polynominal model (RFM) to describe imaging process as geometric model.Describe the wide visual field angle imaging and generally can only lean on traditional collinearity equation model; But describe of the variation clearly definition of the collinearity equation of line array sensor, and the description that parameter is changed and answer and separate and related to a series of data demands such as reference mark distributions to the outer orientation parameter of sensor.
Summary of the invention
The present invention solves the existing in prior technology technical matters; Provide a kind of satellite metadata that makes full use of to consider earth curvature simultaneously; Topographic relief and sensor deformation; The distortion that composite factors such as satellite platform shake cause; Improved the precision of geometric manipulations, and proposed a kind of how much imaging construction methods of the integrality of respective algorithms assurance resolved data towards side-looking intermediate-resolution satellite.
Above-mentioned technical matters of the present invention mainly is able to solve through following technical proposals:
A kind of how much imaging construction methods towards side-looking intermediate-resolution satellite is characterized in that, may further comprise the steps:
Step 1 imports satellite image metadata and satellite almanac data, confirms parameter initial value and other parameter initial values of satellite orbit and attitude fit equation and data chainning node;
Step 2, the image metadata that has comprised projection information according to benchmark tentatively obtains image reference mark initial value automatically;
Step 3, at least three reference mark of manual setting are through the automatic matching image of the method reference mark of gray scale coupling; The required condition of having set according to the step of setting 1 then of parameter initial value, it is regional also to propose the reference mark disappearance automatically according to satellite orbit and the distribution of the automatic detected image of attitude fitting data node location reference mark, carries out changing of reference mark distribution to lacking the zone;
Step 4; Changing of reference mark distribution with completing steps 3; The row-coordinate conversion of going forward side by side; Promptly by projection coordinate's solid the earth's core three-dimensional cartesian coordinate system with forwarding to; Combine the satellite orbit fitting data to be planned under the approximate satellite body coordinate system
again, and make up how much imaging models of satellite;
Step 5 is carried out changing of model parameter iterative calculation according to how much imaging models of the satellite in parameter initial value that has obtained in the step 1 and the step 4;
Step 6, model result and the precision index of output after correcting.
At above-mentioned a kind of how much imaging construction methods towards side-looking intermediate-resolution satellite, in the described step 1, equation and data chainning and other parameter Determination of Initial Values may further comprise the steps:
Step 1.1, three the quadratic polynomial position vector variations of match satellite under the solid geocentric coordinate system in ground respectively, wherein parameters
X 0 , X 1 , X 2 , Y 0 , Y 1 , Y 2 , Z 0 , Z 1 , Z 2 Through the back end in the metadata according to least-squares calculation; Step 1.2, three the quadratic polynomial velocity variations of match satellite under the solid geocentric coordinate system in ground respectively, wherein parameters
V X0 , V X1 , V X2 , V Y0 , V Y1 , V Y2 , V Z0 , V Z1 , V Z2 Through the back end in the metadata according to least-squares calculation; Step 1.3, three data chainnings, each node comprise satellite this time be engraved in attitude angle and this orientation attitude angle rate of change under the body coordinate system, every data chainning comprises 7 ~ 9 data nodes, the swing angle back end of three directions can be expressed as respectively:
, to t constantly the interpolation of certain direction attitude angle according to two 3 Hermite interpolation calculation, for
The angle at t interpolation calculation formula constantly is:
Step 1.4 is laid the initial value of inclination angle as constant offset with the sensor of satellite nominal, and calculates the equivalent focal length initial value of digital sensor according to the satellite flight height h of resolution r and nominal under the star of nominal:
At above-mentioned a kind of how much imaging construction methods towards side-looking intermediate-resolution satellite, in the said step 2, manual setting section reference mark is promptly manual chooses 3-6 reference mark, lays respectively to wait to correct near central and four angles of image.
At above-mentioned a kind of how much imaging construction methods towards side-looking intermediate-resolution satellite, in the described step 3, to carry out reference mark the changing that distribute and promptly carry out autotelicly the reference mark is carried out in disappearance zone, reference mark in the step 3 adding, the concrete operations step is following:
Step 3.1 is set up bipartite graph
, C is the reference mark node set, and each reference mark accounts for two nodes, and P is parameter sets to be asked, and each parameter accounts for a node, if reference mark corresponding node c
iWith parameter p
jRelevant, c then
iWith p
jBetween a limit e is arranged
k, otherwise boundless; Wherein total reference mark number n, number of parameters m,
Then the point among its P only and in the corresponding time period has the limit between the point among the C in the imaging scope among the bigraph (bipartite graph) B; Step 3.2 is according to the maximum match of Hungary algorithm search bigraph (bipartite graph) B; Step 3.3; Generate the result of maximum match in the markers step 3.2; Wherein P concentrates the parameter point do not find coupling to need corresponding increase reference mark, and mark is the range of control of such point, preferentially reference mark of interpolation in by the zone of repeating label repeatedly; Refresh the result who generates maximum match in the weight of structure retry row labels step 3.2 of bipartite graph, up to the complete compatible coupling that finds P to C.
At above-mentioned a kind of how much imaging construction methods towards side-looking intermediate-resolution satellite; In the said step 4; Carry out the control point coordinate transform; How much imaging model principles that make up satellite can be described to: the satellite body coordinate system points to known topocentric vector down through going to picture side's coordinate system after the similarity transformation of space; Itself and picture side's coordinate system point to the vectorial conllinear of known point corresponding diagram picture point down; Its equation is following:
Wherein: angle value
is obtained by the interpolation of data chainning described in the 2-c in
formula,
be the distortion mould
Type, this step adopt five order polynomials to be used for factors such as distortion of fit line array sensor and refractive power and cause
A kind of how much imaging construction methods towards side-looking intermediate-resolution satellite above-mentioned, in the described step 5, carry out changing of model parameter iterative calculation and comprise:
Step 6.1 is set up error equation, makes up linear system and is used for iterative;
Step 6.2, regulation satellite geometry model parameter resolve the order, loop iteration is found the solution;
Step 6.3 is carried out the reference mark according to model and parameter and is just calculated.
At above-mentioned a kind of how much imaging construction methods towards side-looking intermediate-resolution satellite, the concrete operations step of described step 6.1 is following:
Step 6.11; According to imaging model described in the step 4; The initial error equation of setting up the correction model parameter is:
Definition
; Then
is the first order derivative of
, and promptly the coefficient of
equation is respectively:
;
;
;
;
;
;
;
;
;
;
;
;
;
;
;
is easy for expressing, and makes
Then error equation can be expressed as:
, then directly use least square to resolve its normal equation of correction to be:
, in the formula
For
Valuation; Step 6.12;, use single general correction solution to strengthen the reliability of resolving in the parameter calculation when finding the solution a plurality of parameter; Can't restrain or converge to improper value with what avoid that normal equation morbid state causes; Promptly single iteration is introduced one deck internal layer iteration again in order to resolve normal equation in the correction iteration, and the internal layer iteration is rewritten normal equation and is:
, in the formula
Be
iThe correction of inferior calculating,
IBe unit matrix, this solution is that the nothing of least square is estimated partially.
At above-mentioned a kind of how much imaging construction methods towards side-looking intermediate-resolution satellite, the concrete operations step of described step 6.2 is following: step 6.21, resolve equivalent focal length
Coarse value; Step 6.22 is resolved the angular deflection constant
,
,
Coarse value; Step 6.23 is resolved the distortion model coefficient
Coarse value, the constant of changing angular deflection simultaneously
,
,
Step 6.24, repeated execution of steps 6.21 to step 6.23 to reference mark residual error reduces speed less than 10
-2Step 6.25 is revised all angle-data chain node datas
Step 6.26 is refined and is calculated the distortion model coefficient
Residual error no longer diminishes to the some position.
At above-mentioned a kind of how much imaging construction methods towards side-looking intermediate-resolution satellite, the concrete operations step of described step 6.3 is following:
Step 6.31; Equation is resolved as shown in the formula
in the skeleton diagram image position of using quadratic polynomial universal model calculation level position just calculating; Step 6.32; With the corresponding time of sweep trace coordinate of general location;
(
is the sweep speed of sweep trace) is according to inserting instantaneous elements of exterior orientation
in the data chainning in the model parameter; Step 6.33; With initial instantaneous elements of exterior orientation substitution imaging model; Changing of residual values iteration scan line position according to the line coordinate; Computing formula is following:
;
;
,
is respectively two coordinate components in the topocentric coordinates that changes the satellite body coordinate system over to; Step 6.34; With the scanning instantaneous elements of exterior orientation substitution image deformation model constantly of refining; Iteration is obtained sampling point range coordinate; Thereby accomplishing prototype accurately just calculates; Photogrammetric distortion model correction computing formula is following:
,
.
Therefore; The present invention has following advantage: make full use of the satellite metadata and consider earth curvature simultaneously, topographic relief and sensor deformation, the distortion that composite factors such as satellite platform shake cause; Improve the precision of geometric manipulations, and proposed the integrality of respective algorithms assurance resolved data.
Description of drawings
Fig. 1 is a schematic flow sheet of the present invention.
Fig. 2 is the single imaging synoptic diagram of wide cut satellite.
Fig. 3 is an imaged image reference mark initial distribution synoptic diagram.
Fig. 4 is the reference mark residual plot (wherein residual error is amplified 10 times) of indication model output of the present invention.
Fig. 5 is the reference mark residual plot (wherein residual error is amplified 10 times) of multinomial model output.
Embodiment
Pass through embodiment below, and combine accompanying drawing, do further bright specifically technical scheme of the present invention.
Embodiment:
This instance through the wide cut side-looking satellite image that has the satellite partial orbit and decide appearance and locating information with and with known geographic with reference to image; Lay the accurate parameters of trying to achieve the geometric model that the present invention mentions under the unknown prerequisite of matrix at satellite sensor, thereby mapping relations accurately are provided for the geometric exact correction correction.Its flow process is seen Fig. 1.
This instance is chosen 30m resolution (12000 * 12000) environment mitigation moonlet 1 grade data and (is only done overshoot and correct, do not make any geometric manipulations.Its imaging mode be typical wide cut side-looking to sweeping imaging, as shown in Figure 2, wherein S is a projection centre; I is the sweep trace that holds shadow, and I ' is a sweep trace upper thread accordingly, and O is the principal point on the single line; O ' is a principal point millet cake accordingly, T be satellite to arriving, T ' is a path); And with reference to the dem data of the Landsat7 image and the corresponding region SRTM form that have utm projection information as a reference, model building method and the general multinomial model of using the present invention to mention compare test.
Step 1 imports satellite image metadata and satellite almanac data, confirms parameter initial value and other parameter initial values of satellite orbit and attitude fit equation and data chainning node;
(a) xml form meta data file wherein comprises the initial termination time that data are obtained, and calculates sweep speed according to line number;
(b) read the GPS node data of EPH form; Change and rate variation with quadratic polynomial difference match satellite position, fit equation is expressed as:
and:
(c) read the attitude angle node data of ATT form, data are directly increased in the data chainning of attitude match.
(d) calculate the equivalent focal length initial value according to nominal resolution and nominal satellite flying height.
Step 2, the image metadata that has comprised projection information according to benchmark tentatively obtains image reference mark initial value automatically;
Step 3, at least three reference mark of manual setting are through the automatic match map of method of gray scale coupling
The picture reference mark; The required condition of having set according to the step of setting 1 then of parameter initial value; Distribute and the automatic disappearance zone, reference mark that proposes according to satellite orbit and the automatic detected image of attitude fitting data node location reference mark; Carry out changing of reference mark distribution to lacking the zone, in the present embodiment, through the automatic matching image of the method reference mark of gray scale coupling; Matching image reference mark result is totally 358 points; Then according to the required condition of setting of parameter, set up the bipartite graph that expression parameter and position, reference mark concern, and carry out separating of maximum match according to bipartite graph and provide position range and the number that needs are set up the reference mark; And a small amount of reference mark of manual increase, revise the distribution of reference mark, back and show like Fig. 3.This reference mark can be used for general quadratic polynomial model resolves, and it is as shown in Figure 4 to resolve result images point position residual error.
Step 4, with the changing of reference mark distribution of completing steps 3, the row-coordinate conversion of going forward side by side is promptly by projection
The solid the earth's core of coordinate three-dimensional cartesian coordinate system, it is approximate to combine the satellite orbit fitting data to be planned for again with forwarding to
Under the satellite body coordinate system
, and how much imaging models of structure satellite; It should be noted that
Be: this step adopts five order polynomials to be used for factors such as fit line array sensor distortion and refractive power and makes
Wherein, this distortion model is the match of pure mathematics angle, and each parameter can be interpreted as respectively: 0 time item is main
Relevant with sensor translation equal error, 1 item is main, and lay factor such as inclination with sensor relevant, and 2,4
Factors such as inferior main and lens distortion are relevant, and 3,5 items mainly and sensor self deformation and foldings
Factors such as light are relevant.
In the present embodiment,
(a) according to the radiation transformation parameter that provides in the Landsat7 data, coupling and manual corrected image reference mark are changed in the subsidiary Geographic Reference system of reference data, promptly transfer UTM to and divide in the plane right-angle coordinate of being with;
(b) change the some position under the UTM coordinate system over to the earth latitude and longitude coordinates system (this instance adopts the geocentric coordinate system of WGS84 ellipsoid) by the inverse transformation of projective transformation, and from dem data, read the height value of this point with the bilinear interpolation mode;
(c) insert height value in thinking and be a some position geodetic height, the earth latitude and longitude coordinates the earth's core three-dimensional cartesian coordinate system admittedly with changing over to;
(d) the picture line row position calculation imaging time that belongs to according to the reference mark; In insert out the outer orientation element of satellite of instantaneous moment; Position earth rotation linear velocity in conjunction with the satellite place obtains the transition matrix that the satellite body coordinate is consolidated the earth's core three-dimensional cartesian coordinate system with being tied to, further changes ground point over to the satellite body coordinate system.Transition matrix is:
; Wherein:
;
;
; Indication contains time variable in the following formula parameter is like
;
;
;
;
;
, wherein constant
is the angular velocity of earth rotation:
Step 5 is carried out changing of model parameter iterative calculation according to how much imaging models of the satellite in parameter initial value that has obtained in the step 1 and the step 4;
Step 6, model result and the precision index of output after correcting: result of calculation output map picture point position residual error distribution plan is as shown in Figure 5, and the error quantization comparison is as shown in the table in the picture point residual error that itself and general multinomial model are resolved
? | The column direction maximum residul difference | The line direction maximum residul difference | Maximum total residual error | Error in the residual error |
Quadratic polynomial | 2.581551 | 2.728294 | 2.973319 | 0.993866287 |
Model among the present invention | 35.010054 | 2.994242 | 36.230516 | 10.87620294 |
Specific embodiment described herein only is that the present invention's spirit is illustrated.Person of ordinary skill in the field of the present invention can make various modifications or replenishes or adopt similar mode to substitute described specific embodiment, but can't depart from spirit of the present invention or surmount the defined scope of appended claims.
Claims (9)
1. how much imaging construction methods towards side-looking intermediate-resolution satellite is characterized in that, may further comprise the steps:
Step 1 imports satellite image metadata and satellite almanac data, confirms parameter initial value and other parameter initial values of satellite orbit and attitude fit equation and data chainning node;
Step 2, the image metadata that has comprised projection information according to benchmark tentatively obtains image reference mark initial value automatically;
Step 3, at least three reference mark of manual setting are through the automatic matching image of the method reference mark of gray scale coupling; The required condition of having set according to the step of setting 1 then of parameter initial value, it is regional also to propose the reference mark disappearance automatically according to satellite orbit and the distribution of the automatic detected image of attitude fitting data node location reference mark, carries out changing of reference mark distribution to lacking the zone;
Step 4; Changing of reference mark distribution with completing steps 3; The row-coordinate conversion of going forward side by side; Promptly by projection coordinate's solid the earth's core three-dimensional cartesian coordinate system with forwarding to; Combine the satellite orbit fitting data to be planned under the approximate satellite body coordinate system
again, and make up how much imaging models of satellite;
Step 5 is carried out changing of model parameter iterative calculation according to how much imaging models of the satellite in parameter initial value that has obtained in the step 1 and the step 4;
Step 6, model result and the precision index of output after correcting.
2. a kind of how much imaging construction methods towards side-looking intermediate-resolution satellite according to claim 1 is characterized in that in the described step 1, equation and data chainning and other parameter Determination of Initial Values may further comprise the steps:
Step 1.1, three the quadratic polynomial position vector variations of match satellite under the solid geocentric coordinate system in ground respectively, wherein parameters
X 0 , X 1 , X 2 , Y 0 , Y 1 , Y 2 , Z 0 , Z 1 , Z 2 Through the back end in the metadata according to least-squares calculation; Step 1.2, three the quadratic polynomial velocity variations of match satellite under the solid geocentric coordinate system in ground respectively, wherein parameters
V X0 , V X1 , V X2 , V Y0 , V Y1 , V Y2 , V Z0 , V Z1 , V Z2 Through the back end in the metadata according to least-squares calculation; Step 1.3, three data chainnings, each node comprise satellite this time be engraved in attitude angle and this orientation attitude angle rate of change under the body coordinate system, every data chainning comprises 7 ~ 9 data nodes, the swing angle back end of three directions can be expressed as respectively:
, to t constantly the interpolation of certain direction attitude angle according to two 3 Hermite interpolation calculation, for
The angle at t interpolation calculation formula constantly is:
Step 1.4 is laid the initial value of inclination angle as constant offset with the sensor of satellite nominal, and calculates the equivalent focal length initial value of digital sensor according to the satellite flight height h of resolution r and nominal under the star of nominal:
3. a kind of how much imaging construction methods according to claim 1 towards side-looking intermediate-resolution satellite; It is characterized in that; In the said step 2, manual setting section reference mark is promptly manual chooses 3-6 reference mark, lays respectively to wait to correct near central and four angles of image.
4. a kind of how much imaging construction methods according to claim 1 towards side-looking intermediate-resolution satellite; It is characterized in that; In the described step 3; Carry out reference mark the changing that distribute and promptly carry out autotelicly the reference mark is carried out in disappearance zone, reference mark in the step 3 adding, the concrete operations step is following:
Step 3.1 is set up bipartite graph
, C is the reference mark node set, and each reference mark accounts for two nodes, and P is parameter sets to be asked, and each parameter accounts for a node, if reference mark corresponding node c
iWith parameter p
jRelevant, c then
iWith p
jBetween a limit e is arranged
k, otherwise boundless; Wherein total reference mark number n, number of parameters m,
Then the point among its P only and in the corresponding time period has the limit between the point among the C in the imaging scope among the bigraph (bipartite graph) B; Step 3.2 is according to the maximum match of Hungary algorithm search bigraph (bipartite graph) B; Step 3.3; Generate the result of maximum match in the markers step 3.2; Wherein P concentrates the parameter point do not find coupling to need corresponding increase reference mark, and mark is the range of control of such point, preferentially reference mark of interpolation in by the zone of repeating label repeatedly; Refresh the result who generates maximum match in the weight of structure retry row labels step 3.2 of bipartite graph, up to the complete compatible coupling that finds P to C.
5. a kind of how much imaging construction methods according to claim 1 towards side-looking intermediate-resolution satellite; It is characterized in that; In the said step 4; Carry out the reference mark coordinate transform; How much imaging model principles that make up satellite can be described to: the satellite body coordinate system points to known topocentric vector down through going to picture side's coordinate system after the similarity transformation of space; Itself and picture side's coordinate system point to the vectorial conllinear of known point corresponding diagram picture point down, and its equation is following:
Wherein: angle value
is obtained by the interpolation of data chainning described in the 2-c in
formula,
be the distortion mould
Type, this step adopt five order polynomials to be used for factors such as distortion of fit line array sensor and refractive power and cause
6. a kind of how much imaging construction methods towards side-looking intermediate-resolution satellite according to claim 1 is characterized in that, in the described step 5, carry out changing of model parameter iterative calculation and comprise:
Step 6.1 is set up error equation, makes up linear system and is used for iterative;
Step 6.2, regulation satellite geometry model parameter resolve the order, loop iteration is found the solution;
Step 6.3 is carried out the reference mark according to model and parameter and is just calculated.
7. according to claim 6, it is characterized in that the concrete operations step of described step 6.1 is following:
Step 6.11; According to imaging model described in the step 4; The initial error equation of setting up the correction model parameter is:
Definition
; Then
is the first order derivative of
, and promptly the coefficient of
equation is respectively:
;
;
;
;
;
;
;
;
;
;
;
;
;
;
;
Easy for expressing, order
Then error equation can be expressed as:
, then directly use least square to resolve its normal equation of correction to be:
, in the formula
For
Valuation; Step 6.12;, use single general correction solution to strengthen the reliability of resolving in the parameter calculation when finding the solution a plurality of parameter; Can't restrain or converge to improper value with what avoid that normal equation morbid state causes; Promptly single iteration is introduced one deck internal layer iteration again in order to resolve normal equation in the correction iteration, and the internal layer iteration is rewritten normal equation and is:
, in the formula
Be
iThe correction of inferior calculating,
IBe unit matrix, this solution is that the nothing of least square is estimated partially.
8. a kind of how much imaging construction methods towards side-looking intermediate-resolution satellite according to claim 1 is characterized in that the concrete operations step of described step 6.2 is following: step 6.21, resolve equivalent focal length
Coarse value; Step 6.22 is resolved the angular deflection constant
,
,
Coarse value; Step 6.23 is resolved the distortion model coefficient
Coarse value, the constant of changing angular deflection simultaneously
,
,
Step 6.24, repeated execution of steps 6.21 to step 6.23 to reference mark residual error reduces speed less than 10
-2Step 6.25 is revised all angle-data chain node datas
Step 6.26 is refined and is calculated the distortion model coefficient
Residual error no longer diminishes to the some position.
9. a kind of how much imaging construction methods towards side-looking intermediate-resolution satellite according to claim 1 is characterized in that the concrete operations step of described step 6.3 is following:
Step 6.31; Equation is resolved as shown in the formula
in the skeleton diagram image position of using quadratic polynomial universal model calculation level position just calculating; Step 6.32; With the corresponding time of sweep trace coordinate of general location;
(
is the sweep speed of sweep trace) is according to inserting instantaneous elements of exterior orientation
in the data chainning in the model parameter; Step 6.33; With initial instantaneous elements of exterior orientation substitution imaging model; Changing of residual values iteration scan line position according to the line coordinate; Computing formula is following:
;
;
,
is respectively two coordinate components in the topocentric coordinates that changes the satellite body coordinate system over to; Step 6.34; With the scanning instantaneous elements of exterior orientation substitution image deformation model constantly of refining; Iteration is obtained sampling point range coordinate; Thereby accomplishing prototype accurately just calculates; Photogrammetric distortion model correction computing formula is following:
,
.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110388915 CN102508260B (en) | 2011-11-30 | 2011-11-30 | Geometric imaging construction method for side-looking medium resolution ratio satellite |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110388915 CN102508260B (en) | 2011-11-30 | 2011-11-30 | Geometric imaging construction method for side-looking medium resolution ratio satellite |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102508260A true CN102508260A (en) | 2012-06-20 |
CN102508260B CN102508260B (en) | 2013-10-23 |
Family
ID=46220366
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110388915 Expired - Fee Related CN102508260B (en) | 2011-11-30 | 2011-11-30 | Geometric imaging construction method for side-looking medium resolution ratio satellite |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102508260B (en) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102760301A (en) * | 2012-06-07 | 2012-10-31 | 何潇 | Image control point interactive editing system |
CN102901519A (en) * | 2012-11-02 | 2013-01-30 | 武汉大学 | optical push-broom satellite in-orbit stepwise geometric calibration method based on probe element direction angle |
CN103115627A (en) * | 2013-01-21 | 2013-05-22 | 武汉大学 | Multi-track combined on-track geometrical checking method of remote sensing satellite linear array sensor |
CN106846385A (en) * | 2016-12-30 | 2017-06-13 | 广州地理研究所 | Many sensing Remote Sensing Images Matching Methods, device and system based on unmanned plane |
CN107085856A (en) * | 2017-04-06 | 2017-08-22 | 上海航天测控通信研究所 | A kind of in-orbit high-precision real-time location method based on optical image |
CN107330136A (en) * | 2017-05-16 | 2017-11-07 | 上海卫星工程研究所 | It is adapted to bilateral and regards work remote sensing satellite navigation antenna layout method |
CN107689064A (en) * | 2017-08-08 | 2018-02-13 | 武汉大学 | Take the strict geometry imaging model construction method of satellite optical of aberration correction into account |
CN108593556A (en) * | 2017-12-26 | 2018-09-28 | 中国科学院电子学研究所 | Satellite imagery geometry based on vector characteristic is refined the construction method of model |
CN108761453A (en) * | 2018-03-29 | 2018-11-06 | 航天东方红卫星有限公司 | A kind of imaged viewing angle optimization method of optical satellite and SAR Satellite Images Fusions |
CN111486850A (en) * | 2020-04-20 | 2020-08-04 | 中国科学院微小卫星创新研究院 | In-orbit frame frequency adjustment strategy for earth observation satellite |
CN113066132A (en) * | 2020-03-16 | 2021-07-02 | 天目爱视(北京)科技有限公司 | 3D modeling calibration method based on multi-device acquisition |
CN113393571A (en) * | 2021-06-15 | 2021-09-14 | 北京市遥感信息研究所 | Cloud-free satellite image generation method and device |
CN114882105A (en) * | 2022-05-18 | 2022-08-09 | 武汉大学 | Improved satellite ground track intersection point extraction method |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110500995B (en) * | 2019-07-12 | 2020-06-23 | 武汉大学 | Method for establishing high-resolution satellite image equivalent geometric imaging model by using RPC parameters |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001091649A (en) * | 1999-09-20 | 2001-04-06 | Mitsubishi Electric Corp | Ground control point device for synthetic aperture radar image precise geometrical correction |
US20030044085A1 (en) * | 2001-05-01 | 2003-03-06 | Dial Oliver Eugene | Apparatuses and methods for mapping image coordinates to ground coordinates |
CN101598785A (en) * | 2009-06-30 | 2009-12-09 | 武汉大学 | The rational function imaging model generation method of satellite image products at different levels |
CN101839713A (en) * | 2010-04-20 | 2010-09-22 | 武汉大学 | Satellite image system error correction method based on bias matrix with time factor |
-
2011
- 2011-11-30 CN CN 201110388915 patent/CN102508260B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001091649A (en) * | 1999-09-20 | 2001-04-06 | Mitsubishi Electric Corp | Ground control point device for synthetic aperture radar image precise geometrical correction |
US20030044085A1 (en) * | 2001-05-01 | 2003-03-06 | Dial Oliver Eugene | Apparatuses and methods for mapping image coordinates to ground coordinates |
CN101598785A (en) * | 2009-06-30 | 2009-12-09 | 武汉大学 | The rational function imaging model generation method of satellite image products at different levels |
CN101839713A (en) * | 2010-04-20 | 2010-09-22 | 武汉大学 | Satellite image system error correction method based on bias matrix with time factor |
Non-Patent Citations (1)
Title |
---|
谢仁伟等: "《基于控制点库的SIFT多源影像自动配准方法》", 《测绘科学》, vol. 36, no. 4, 31 July 2011 (2011-07-31) * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102760301A (en) * | 2012-06-07 | 2012-10-31 | 何潇 | Image control point interactive editing system |
CN102760301B (en) * | 2012-06-07 | 2015-09-30 | 武汉天目经纬科技有限公司 | image control point interactive editing system |
CN102901519A (en) * | 2012-11-02 | 2013-01-30 | 武汉大学 | optical push-broom satellite in-orbit stepwise geometric calibration method based on probe element direction angle |
CN102901519B (en) * | 2012-11-02 | 2015-04-29 | 武汉大学 | optical push-broom satellite in-orbit stepwise geometric calibration method based on probe element direction angle |
CN103115627A (en) * | 2013-01-21 | 2013-05-22 | 武汉大学 | Multi-track combined on-track geometrical checking method of remote sensing satellite linear array sensor |
CN106846385A (en) * | 2016-12-30 | 2017-06-13 | 广州地理研究所 | Many sensing Remote Sensing Images Matching Methods, device and system based on unmanned plane |
CN106846385B (en) * | 2016-12-30 | 2020-11-24 | 广州地理研究所 | Multi-sensing remote sensing image matching method, device and system based on unmanned aerial vehicle |
CN107085856A (en) * | 2017-04-06 | 2017-08-22 | 上海航天测控通信研究所 | A kind of in-orbit high-precision real-time location method based on optical image |
CN107085856B (en) * | 2017-04-06 | 2020-01-31 | 上海航天测控通信研究所 | on-orbit high-precision real-time positioning method based on optical images |
CN107330136A (en) * | 2017-05-16 | 2017-11-07 | 上海卫星工程研究所 | It is adapted to bilateral and regards work remote sensing satellite navigation antenna layout method |
CN107330136B (en) * | 2017-05-16 | 2020-07-14 | 上海卫星工程研究所 | Layout method of remote sensing satellite navigation antenna suitable for double-side-view work |
CN107689064A (en) * | 2017-08-08 | 2018-02-13 | 武汉大学 | Take the strict geometry imaging model construction method of satellite optical of aberration correction into account |
CN108593556A (en) * | 2017-12-26 | 2018-09-28 | 中国科学院电子学研究所 | Satellite imagery geometry based on vector characteristic is refined the construction method of model |
CN108761453A (en) * | 2018-03-29 | 2018-11-06 | 航天东方红卫星有限公司 | A kind of imaged viewing angle optimization method of optical satellite and SAR Satellite Images Fusions |
CN108761453B (en) * | 2018-03-29 | 2021-02-05 | 航天东方红卫星有限公司 | Imaging view angle optimization method for image fusion of optical satellite and SAR (synthetic aperture radar) satellite |
CN113066132A (en) * | 2020-03-16 | 2021-07-02 | 天目爱视(北京)科技有限公司 | 3D modeling calibration method based on multi-device acquisition |
CN111486850A (en) * | 2020-04-20 | 2020-08-04 | 中国科学院微小卫星创新研究院 | In-orbit frame frequency adjustment strategy for earth observation satellite |
CN111486850B (en) * | 2020-04-20 | 2021-10-26 | 中国科学院微小卫星创新研究院 | In-orbit frame frequency adjustment strategy for earth observation satellite |
CN113393571A (en) * | 2021-06-15 | 2021-09-14 | 北京市遥感信息研究所 | Cloud-free satellite image generation method and device |
CN113393571B (en) * | 2021-06-15 | 2023-06-27 | 北京市遥感信息研究所 | Cloud-free satellite image generation method and device |
CN114882105A (en) * | 2022-05-18 | 2022-08-09 | 武汉大学 | Improved satellite ground track intersection point extraction method |
CN114882105B (en) * | 2022-05-18 | 2024-09-20 | 武汉大学 | Improved satellite ground track intersection point extraction method |
Also Published As
Publication number | Publication date |
---|---|
CN102508260B (en) | 2013-10-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102508260A (en) | Geometric imaging construction method for side-looking medium resolution ratio satellite | |
CN105627991B (en) | A kind of unmanned plane image real time panoramic joining method and system | |
CN103822615B (en) | A kind of multi-control point extracts and the unmanned aerial vehicle target real-time location method be polymerized automatically | |
CN102506824B (en) | Method for generating digital orthophoto map (DOM) by urban low altitude unmanned aerial vehicle | |
CN104848860B (en) | A kind of agile satellite imagery process attitude maneuver planing method | |
CN103345737B (en) | A kind of UAV high resolution image geometric correction method based on error compensation | |
CN102645209B (en) | Joint positioning method for spatial points by means of onboard LiDAR point cloud and high resolution images | |
CN101852623B (en) | On-track calibration method for internal element of satellite optical remote sensing camera | |
CN106780321B (en) | CBERS-02 satellite HR sensor image overall tight orientation and correction splicing method | |
CN107144293A (en) | A kind of geometric calibration method of video satellite area array cameras | |
CN102735216B (en) | CCD stereoscopic camera three-line imagery data adjustment processing method | |
CN108226982B (en) | Single linear array satellite laser combined high-precision positioning processing method | |
CN104820984B (en) | A kind of satellite remote sensing three line scanner stereopsis processing system and method | |
CN107341778A (en) | SAR image ortho-rectification methods based on satellite control point storehouse and DEM | |
CN103129752A (en) | Dynamic compensation method for attitude angle errors of optical remote sensing satellite based on ground navigation | |
CN102322863B (en) | Remote sensing satellite multi-satellite combined converse orbit and attitude determination method | |
CN102721410B (en) | Island-air triangular measuring method based on GPS/IMU positioning and orientating technology | |
CN104978761B (en) | A kind of RPC models correction coefficient acquisition methods | |
CN112698664A (en) | Sight line sector dynamic estimation method for unmanned cluster collaborative navigation optimization | |
CN106443676A (en) | Scarce control point space-borne synthetic aperture radar image ground positioning method | |
CN107316280A (en) | Li Island satellite image RPC models high accuracy geometry location method | |
CN101477682A (en) | Method for remotely sensing image geometric correction by weighted polynomial model | |
CN110632636B (en) | Carrier attitude estimation method based on Elman neural network | |
CN102147249B (en) | Method for precisely correcting satellite-borne optical linear array image based on linear characteristic | |
CN117194869B (en) | Attitude-considered low-orbit satellite antenna phase center forecasting and fitting 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20131023 Termination date: 20191130 |