CN110793542A - Area array optical remote sensing satellite in-orbit geometric calibration method based on generalized probe element pointing angle - Google Patents
Area array optical remote sensing satellite in-orbit geometric calibration method based on generalized probe element pointing angle Download PDFInfo
- Publication number
- CN110793542A CN110793542A CN201910951317.6A CN201910951317A CN110793542A CN 110793542 A CN110793542 A CN 110793542A CN 201910951317 A CN201910951317 A CN 201910951317A CN 110793542 A CN110793542 A CN 110793542A
- Authority
- CN
- China
- Prior art keywords
- area array
- probe element
- pointing angle
- remote sensing
- coordinate system
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
Abstract
The invention relates to an in-orbit geometric calibration method for an area array optical remote sensing satellite based on a generalized probe element pointing angle, belongs to the technical field of photogrammetry and remote sensing, and particularly relates to in-orbit geometric calibration for an area array optical remote sensing satellite based on a generalized probe element pointing angle. The invention provides an in-orbit geometric calibration model of an area array optical remote sensing satellite based on a generalized probe element pointing angle, which is established by analyzing the influence of each physical parameter of a strict imaging geometric model on positioning accuracy by utilizing a projection ray formed by an image point and a corresponding instantaneous projection center from imaging of the area array optical remote sensing satellite and establishing the generalized CCD probe element pointing angle under an attitude measurement reference coordinate.
Description
Technical Field
The invention relates to an in-orbit geometric calibration method for an area array optical remote sensing satellite based on a generalized probe element pointing angle, belongs to the technical field of photogrammetry and remote sensing, and particularly relates to in-orbit geometric calibration for an area array optical remote sensing satellite based on a generalized probe element pointing angle.
Background
At present, a part of high-resolution optical remote sensing satellites are provided with an area array remote sensing camera, and images of a target area are obtained in a frame format mode. At present, a high-resolution satellite remote sensing earth observation technology becomes one of important means for acquiring earth space information by human beings, and a series of products produced by image data of the high-resolution satellite remote sensing earth observation technology are widely applied to the fields of military reconnaissance, topographic mapping, homeland resource investigation, urban planning and the like. In order to realize the productization of the area array framed satellite image more quickly, better serve the requirements of social sustainable development and create more social and economic benefits, firstly, the problem of accurate target positioning of the area array optical remote sensing satellite image needs to be solved, and the in-orbit geometric calibration is the key for realizing the accurate target positioning of the optical satellite image without ground control.
In order to realize the accurate positioning of the ground target of the area array frame type optical remote sensing satellite image, all imaging parameters (such as camera principal distance, image principal point position, camera installation position and the like) of the remote sensing camera during image acquisition must be accurately obtained, and a strict geometric relation between the area array frame type optical remote sensing satellite image and the ground target is established on the basis of the imaging parameters. Before satellite transmission, a camera is calibrated in a laboratory to determine system parameters such as relative and absolute positions of a CCD (charge coupled device), a principal distance and a principal point position, and the calibrated imaging parameters are provided for a user. The advantages of the geometric calibration of the laboratory are the use of professional equipment, the normalization and standardization of the operation process and high calibration precision. However, due to the influence of excessive acceleration, impulsive force and various disturbance forces during satellite launching, the change of space environments such as temperature, humidity and air pressure during satellite in-orbit running, device loss and aging caused by long-term use of an imaging remote sensor and other factors can cause different differences between imaging parameters and laboratory calibration values before satellite launching, the laboratory calibration values are used for carrying out area array optical remote sensing satellite image to ground target positioning, system errors are inevitably introduced, and the image to ground target positioning precision is reduced. Therefore, how to fully utilize the ground control information to accurately solve the imaging parameters of the optical satellite sensor during in-orbit running and realize the high-precision in-orbit geometric calibration of the satellite remote sensor is a key problem which needs to be solved urgently by photogrammetry students.
The existing on-orbit geometric calibration method generally takes a high-precision manual fixed target of a geometric calibration field or a DEM/DOM product as a ground control condition, the available geometric calibration fields are relatively few and mainly comprise Henan Songshan and Henan Anyang geometric calibration fields, the high-precision manual fixed target has high maintenance cost, the DOM/DEM product also has the problem of difficult updating and needs a large amount of investment of manpower, material resources and financial resources, the existing calibration method is a step-by-step geometric calibration method and has the defects that ① needs to provide a laboratory calibration value, ② internal and external calibration parameters have strong correlation, and a mode of automatically matching DOM and an image to be calibrated needs to be adopted to obtain a large amount of ground control points.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the method completely overcomes the correlation between internal and external parameters caused by step calibration under the condition of no need of laboratory calibration parameters and a small number of ground control points, and realizes the high-precision in-orbit geometric calibration of the area array optical remote sensing satellite.
The technical solution of the invention is as follows:
an in-orbit geometric calibration method for an area array optical remote sensing satellite based on a generalized probe element pointing angle comprises the following steps:
step 1, establishing an in-orbit geometric calibration model of an area array optical remote sensing camera based on a generalized probe element pointing angle under a posture measurement reference coordinate, wherein the in-orbit geometric calibration model is shown as a formula (1);
wherein the content of the first and second substances,
for ground control point under WGS-84 coordinate systemThree-dimensional coordinates of (a);
is a three-dimensional coordinate of the phase center of the GPS antenna under a WGS-84 coordinate system;
m is a scale factor;
the method comprises the following steps of (1) setting a mounting matrix of a camera under a satellite body coordinate system;
the pointing angle of the CCD probe corresponding to the ground control point in the camera coordinate system;
f is the main distance of the camera;
on the basis of the formula (1), the
(x, y, z) describes the coordinates of each probe element of the area array CCD under the attitude measurement coordinate system, and further enables:
(ψ'Y,ψ'X) The pointing angle of each probe element of the area array CCD under the attitude measurement coordinate system is a constant,
dividing the two sides of the medium sign in the formula (2) by-z to obtain
(iv) according to formula (4) (#)'Y,ψ'X) The external calibration parameters and the internal calibration parameters can be described simultaneously; the external calibration parameters include (D)x,Dy,Dz)、(dx,dy,dz)、Andthe internal calibration parameters includeAnd f;
by transforming the formula (1)
Wherein, λ is a scale factor, λ ═ mz;
step 2, establishing a directive angle model of each probe element of the area array CCD under a satellite attitude measurement reference coordinate system as follows,
wherein (a)0,a1,a2,a3,a4,a5,a6,a7,a8,a9,b0,b1,b2,b3,b4,b5,b6,b7,b8,b9) The (l, s) are the row number and the column number of each probe element of the area array CCD;
and 3, solving the pointing angle model coefficient in the step 2 to realize the following steps,
a scale factor lambda is obtained from the equation (5) from each ground control point, and a pointing angle phi of each probe element of the area array CCD corresponding to the ground control point in the attitude measurement coordinate system is obtained from the equation (5)'Y,ψ'X);
Establishing an error equation as shown in the formula (7);
V=AW-L (7)
in the formula (I), the compound is shown in the specification,
W=[a0a1a2a3a4a5a6a7a8a9b0b1b2b3b4b5b6b7b8b9]
then solving an unknown number W according to a least square adjustment principle;
W=(ATA)-1ATL
and 4, solving the pointing angle of each probe element of the area array CCD under the satellite attitude measurement reference coordinate system according to the pointing angle model established in the step 2 and the pointing angle model coefficient solved in the step 3, and finishing the in-orbit geometric calibration of the area array optical remote sensing satellite.
Advantageous effects
(1) The invention provides an in-orbit geometric calibration model of an area array optical remote sensing satellite based on a generalized probe element pointing angle, which is established by analyzing the influence of each physical parameter of a strict imaging geometric model on positioning accuracy by utilizing a projection ray formed by an image point and a corresponding instantaneous projection center from imaging of the area array optical remote sensing satellite and establishing the generalized CCD probe element pointing angle under an attitude measurement reference coordinate.
(2) The method can perform accurate in-orbit geometric calibration on the optical satellite sensor by taking few control points as ground control conditions under the condition that ground control points are difficult to obtain, further eliminate the system error of the positioning result of the area array frame type satellite image to the ground target and obviously improve the positioning accuracy of the image to the ground target without ground control.
Detailed Description
The embodiment provides an in-orbit geometric calibration method of an area array optical remote sensing satellite based on a generalized probe element pointing angle, which comprises the following steps:
step 1, preparing data, acquiring an optical remote sensing satellite image of an area array to be calibrated and a corresponding reference image, selecting a control point on the image to be calibrated and the reference data, wherein the row number and the column number of an area array CCD (charge coupled device) of the control point on the image to be calibrated are (l, s), and the coordinate of the control point under a WGS84 coordinate system is (l, s)Acquiring satellite attitude quaternion (q) at corresponding imaging time according to the imaging time of the image to be calibrated1,q2,q3,q4) Acquiring the coordinates of the corresponding imaging time satellite in a WGS84 coordinate system
Step 2, establishing an in-orbit geometric calibration model of the area array optical remote sensing camera based on the pointing angle of the generalized probe elementWherein:PN (t) is a time and nutation matrix, R (t) is a rotation matrix of the earth, and W (t) is a polar shift matrix.
Step 3, establishing a directive angle model of each probe element of the area array CCD under a satellite attitude measurement reference coordinate system:
step 4, solving the pointing angle model coefficient in the step 3, and utilizing the line number and the column number of the area array CCD of the control point on the image to be calibrated and the coordinates of the control point under the WGS84 coordinate system to calibrate the model according to the in-orbit geometric calibration model of the area array optical remote sensing cameraThe first equation of the equation obtains a scaling factor lambda, and the second equation and the third equation of the equation obtain a pointing angle (psi ') of each probe element of the area array CCD corresponding to the ground control point in an attitude measurement coordinate system'Y,ψ'X) (ii) a Then, an error equation V ═ AW-L is established, which
W=[a0a1a2a3a4a5a6a7a8a9b0b1b2b3b4b5b6b7b8b9],According to the least square adjustment principle, according to the formula W ═ ATA)-1ATAnd L, solving the unknown number W, wherein W is the pointing angle model coefficient.
Step 5, establishing a pointing angle model according to the step 3And the orientation angle model coefficient (a) solved in step 30,a1,a2,a3,a4,a5,a6,a7,a8,a9,b0,b1,b2,b3,b4,b5,b6,b7,b8,b9) And solving the pointing angle of each probe element of the area array CCD under the satellite attitude measurement reference coordinate system to finish the in-orbit geometric calibration of the area array optical remote sensing satellite.
Claims (10)
1. An in-orbit geometric calibration method for an area array optical remote sensing satellite based on a generalized probe element pointing angle is characterized by comprising the following steps:
step 1, establishing an in-orbit geometric calibration model of an area array optical remote sensing camera based on a generalized probe element pointing angle under a posture measurement reference coordinate, wherein the in-orbit geometric calibration model is shown as a formula (1);
wherein the content of the first and second substances,
is a three-dimensional coordinate of the phase center of the GPS antenna under a WGS-84 coordinate system;
m is a scale factor;
a rotation matrix from the star sensor coordinate system to a J2000 coordinate system;
the method comprises the following steps of (1) setting a mounting matrix of a camera under a satellite body coordinate system;
the pointing angle of the CCD probe corresponding to the ground control point in the camera coordinate system;
f is the main distance of the camera;
on the basis of the formula (1), the
(x, y, z) describes the coordinates of each probe element of the area array CCD under the attitude measurement coordinate system, and further enables:
(ψ'Y,ψ'X) Is flourThe pointing angle of each probe element of the array CCD under the attitude measurement coordinate system is a constant;
dividing the two sides of the medium sign in the formula (2) by-z to obtain
(iv) according to formula (4) (#)'Y,ψ'X) Describing an external calibration parameter and an internal calibration parameter at the same time; the external calibration parameters include (D)x,Dy,Dz)、(dx,dy,dz)、Andthe internal calibration parameters includeAnd f;
by transforming the formula (1)
Wherein, λ is a scale factor, λ ═ mz;
step 2, establishing a directive angle model of each probe element of the area array CCD under a satellite attitude measurement reference coordinate system as follows,
wherein (a)0,a1,a2,a3,a4,a5,a6,a7,a8,a9,b0,b1,b2,b3,b4,b5,b6,b7,b8,b9) For the orientation angle model coefficient, (l, s) are the row number and the column number of each probe element of the area array CCD;
And 3, solving the pointing angle model coefficient in the step 2, wherein the method comprises the following steps:
a scale factor lambda is obtained from the equation (5) from each ground control point, and a pointing angle phi of each probe element of the area array CCD corresponding to the ground control point in the attitude measurement coordinate system is obtained from the equation (5)'Y,ψ'X);
Establishing an error equation as shown in the formula (7);
V=AW-L (7)
in the formula (I), the compound is shown in the specification,
W=[a0a1a2a3a4a5a6a7a8a9b0b1b2b3b4b5b6b7b8b9]
then solving an unknown number W according to a least square adjustment principle;
W=(ATA)-1ATL
and 4, solving the pointing angle of each probe element of the area array CCD under the satellite attitude measurement reference coordinate system according to the pointing angle model established in the step 2 and the pointing angle model coefficient solved in the step 3, and finishing the in-orbit geometric calibration of the area array optical remote sensing satellite.
2. The method for geometrically calibrating the in-orbit of the area array optical remote sensing satellite based on the pointing angle of the generalized probe element as claimed in claim 1, wherein: the row number of the area array CCD of the ground control point on the image to be marked is l.
3. The method for geometrically calibrating the in-orbit of the area array optical remote sensing satellite based on the pointing angle of the generalized probe element as claimed in claim 1, wherein: the column number of the area array CCD of the ground control point on the image to be marked is s.
4. The method for geometrically calibrating the in-orbit of the area array optical remote sensing satellite based on the pointing angle of the generalized probe element as claimed in claim 1, wherein: satellite attitude quaternion (q) of imaging time1,q2,q3,q4)。
6. the method for geometrically calibrating the in-orbit of the area array optical remote sensing satellite based on the pointing angle of the generalized probe element as claimed in claim 5, wherein:PN (t) is a time and nutation matrix, R (t) is a rotation matrix of the earth, and W (t) is a polar shift matrix.
7. The method for geometrically calibrating the in-orbit of the area array optical remote sensing satellite based on the pointing angle of the generalized probe element as claimed in claim 1, wherein: the directive angle model of each probe element of the area array CCD under the satellite attitude measurement reference coordinate system is as follows:
8. the method for geometrically calibrating the in-orbit of the area array optical remote sensing satellite based on the pointing angle of the generalized probe element as claimed in claim 7, wherein:
error equation
V=AW-L。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910951317.6A CN110793542A (en) | 2019-10-08 | 2019-10-08 | Area array optical remote sensing satellite in-orbit geometric calibration method based on generalized probe element pointing angle |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910951317.6A CN110793542A (en) | 2019-10-08 | 2019-10-08 | Area array optical remote sensing satellite in-orbit geometric calibration method based on generalized probe element pointing angle |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110793542A true CN110793542A (en) | 2020-02-14 |
Family
ID=69440121
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910951317.6A Pending CN110793542A (en) | 2019-10-08 | 2019-10-08 | Area array optical remote sensing satellite in-orbit geometric calibration method based on generalized probe element pointing angle |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110793542A (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111524196A (en) * | 2020-04-30 | 2020-08-11 | 中国科学院微小卫星创新研究院 | In-orbit geometric calibration method for sweep large-width optical satellite |
CN112802118A (en) * | 2021-01-05 | 2021-05-14 | 湖北工业大学 | On-orbit time-sharing geometric calibration method for optical satellite sensor |
CN113074717A (en) * | 2021-03-24 | 2021-07-06 | 中国科学院空天信息创新研究院 | Method for acquiring scientific satellite observation direction |
CN113313769A (en) * | 2021-06-11 | 2021-08-27 | 湖北工业大学 | Seamless geometric calibration method between optical satellite multi-area array sensor chips |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102901519A (en) * | 2012-11-02 | 2013-01-30 | 武汉大学 | optical push-broom satellite in-orbit stepwise geometric calibration method based on probe element direction angle |
US20140058695A1 (en) * | 2012-08-23 | 2014-02-27 | Raytheon Company | Geometric calibration of a remote sensor |
CN103674063A (en) * | 2013-12-05 | 2014-03-26 | 中国资源卫星应用中心 | On-orbit geometric calibration method of optical remote sensing camera |
CN104764443A (en) * | 2015-04-24 | 2015-07-08 | 国家测绘地理信息局卫星测绘应用中心 | Optical remote sensing satellite rigorous imaging geometrical model building method |
CN109084963A (en) * | 2018-08-02 | 2018-12-25 | 北京空间机电研究所 | A kind of remote sensor on-orbit calibration light source emission system |
-
2019
- 2019-10-08 CN CN201910951317.6A patent/CN110793542A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140058695A1 (en) * | 2012-08-23 | 2014-02-27 | Raytheon Company | Geometric calibration of a remote sensor |
CN102901519A (en) * | 2012-11-02 | 2013-01-30 | 武汉大学 | optical push-broom satellite in-orbit stepwise geometric calibration method based on probe element direction angle |
CN103674063A (en) * | 2013-12-05 | 2014-03-26 | 中国资源卫星应用中心 | On-orbit geometric calibration method of optical remote sensing camera |
CN104764443A (en) * | 2015-04-24 | 2015-07-08 | 国家测绘地理信息局卫星测绘应用中心 | Optical remote sensing satellite rigorous imaging geometrical model building method |
CN109084963A (en) * | 2018-08-02 | 2018-12-25 | 北京空间机电研究所 | A kind of remote sensor on-orbit calibration light source emission system |
Non-Patent Citations (3)
Title |
---|
曹金山 等: ""资源三号卫星成像在轨几何定标的探元指向角法"", 《测绘学报》 * |
隋起胜 等: "《太空安全防御技术》", 30 June 2018, 国防工业出版社 * |
雷星晖 等: ""天绘一号卫星基于密集控制点的定位误差分析和改正"", 《枫林学院》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111524196A (en) * | 2020-04-30 | 2020-08-11 | 中国科学院微小卫星创新研究院 | In-orbit geometric calibration method for sweep large-width optical satellite |
CN111524196B (en) * | 2020-04-30 | 2023-04-07 | 中国科学院微小卫星创新研究院 | In-orbit geometric calibration method for sweep large-width optical satellite |
CN112802118A (en) * | 2021-01-05 | 2021-05-14 | 湖北工业大学 | On-orbit time-sharing geometric calibration method for optical satellite sensor |
CN113074717A (en) * | 2021-03-24 | 2021-07-06 | 中国科学院空天信息创新研究院 | Method for acquiring scientific satellite observation direction |
CN113074717B (en) * | 2021-03-24 | 2023-04-14 | 中国科学院空天信息创新研究院 | Method for acquiring scientific satellite observation direction |
CN113313769A (en) * | 2021-06-11 | 2021-08-27 | 湖北工业大学 | Seamless geometric calibration method between optical satellite multi-area array sensor chips |
CN113313769B (en) * | 2021-06-11 | 2022-08-30 | 湖北工业大学 | Seamless geometric calibration method between optical satellite multi-area array sensor chips |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110793542A (en) | Area array optical remote sensing satellite in-orbit geometric calibration method based on generalized probe element pointing angle | |
CN104897175B (en) | Polyphaser optics, which is pushed away, sweeps the in-orbit geometric calibration method and system of satellite | |
CN108592950B (en) | Calibration method for relative installation angle of monocular camera and inertial measurement unit | |
CN110390697B (en) | Millimeter wave radar and camera combined calibration method based on LM algorithm | |
Cramer et al. | System calibration for direct georeferencing | |
CN107144293A (en) | A kind of geometric calibration method of video satellite area array cameras | |
CN110220491B (en) | Method for estimating installation error angle of optical pod of unmanned aerial vehicle | |
CN109459058B (en) | Calibration method of multi-view-field star sensor based on three-axis turntable | |
CN102346033B (en) | Direct positioning method and system based on satellite observation angle error estimation | |
CN109696182A (en) | A kind of spaceborne push-broom type optical sensor elements of interior orientation calibrating method | |
CN103390102B (en) | A kind of calculation method of satellite image solid intersection angle | |
CN108663043B (en) | Single-camera-assisted distributed POS main node and sub node relative pose measurement method | |
CN106895851A (en) | A kind of sensor calibration method that many CCD polyphasers of Optical remote satellite are uniformly processed | |
CN107316280B (en) | Li Island satellite image RPC model high-precision geometry location method | |
CN110006452B (en) | Relative geometric calibration method and system for high-resolution six-size wide-view-field camera | |
CN111473802A (en) | Optical sensor internal orientation element calibration method based on linear array push-scanning | |
CN105698766A (en) | Satellite image RFM model block adjustment method with orientation parameter precision information taken into consideration | |
CN107798668B (en) | Unmanned aerial vehicle imaging hyperspectral geometric correction method and system based on RGB images | |
CN108447100B (en) | Method for calibrating eccentricity vector and visual axis eccentricity angle of airborne three-linear array CCD camera | |
CN111524196B (en) | In-orbit geometric calibration method for sweep large-width optical satellite | |
CN108154535B (en) | Camera calibration method based on collimator | |
CN114972078B (en) | Method and system for improving uncontrolled geometric quality of domestic optical satellite image by SAR image | |
CN105444778A (en) | Star sensor in-orbit attitude determination error obtaining method based on imaging geometric inversion | |
CN114858133B (en) | Attitude low-frequency error correction method under fixed star observation mode | |
CN113538595A (en) | Method for improving geometric precision of remote sensing stereo image by using laser height measurement data in auxiliary manner |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200214 |