CN113034610A - Astronomical technology-based spatial direction measuring instrument calibration method - Google Patents

Astronomical technology-based spatial direction measuring instrument calibration method Download PDF

Info

Publication number
CN113034610A
CN113034610A CN202110285856.8A CN202110285856A CN113034610A CN 113034610 A CN113034610 A CN 113034610A CN 202110285856 A CN202110285856 A CN 202110285856A CN 113034610 A CN113034610 A CN 113034610A
Authority
CN
China
Prior art keywords
landmark
image
feature point
representing
point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110285856.8A
Other languages
Chinese (zh)
Other versions
CN113034610B (en
Inventor
吴奋陟
王立
郑然�
武延鹏
刘鲁
王苗苗
王艳宝
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Control Engineering
Original Assignee
Beijing Institute of Control Engineering
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Control Engineering filed Critical Beijing Institute of Control Engineering
Priority to CN202110285856.8A priority Critical patent/CN113034610B/en
Publication of CN113034610A publication Critical patent/CN113034610A/en
Application granted granted Critical
Publication of CN113034610B publication Critical patent/CN113034610B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • G06T5/80

Abstract

The invention discloses a space-pointing measuring instrument calibration method based on an astronomical technique, and belongs to the technical field of space-pointing measuring instrument calibration. The method comprises the steps of utilizing terrestrial landmark features on the earth surface to realize on-orbit calibration of a space pointing measuring instrument, firstly extracting and identifying a plurality of terrestrial landmark feature points from a shot earth optical image, and calculating image coordinates of a terrestrial landmark; secondly, carrying out radial distortion, eccentric distortion and image plane distortion correction on the image coordinates of the landmark feature points, and calculating the observation vector of the landmark feature points in the system of the instrument measurement; and then, calculating landmark vectors obtained by matching in the earth landmark library, corresponding the landmark vectors to the observation vectors one by one, and correcting the focal length, the principal point, the radial distortion and the eccentric distortion of the sensor by using a least square method.

Description

Astronomical technology-based spatial direction measuring instrument calibration method
Technical Field
The invention relates to a space-pointing measuring instrument calibration method based on an astronomical technique, and belongs to the technical field of space-pointing measuring instrument calibration.
Background
The space-oriented measuring sensor is used for guiding load aiming targets such as optical cameras and laser weapons on the spacecraft, and has urgent requirements in the fields of low-orbit ultrahigh-scale military surveying and mapping satellites, high-orbit remote sensing satellites, space-based attack defense satellites and the like. Usually, the product is subjected to ground calibration before satellite launching, however, the ground calibration parameters of the product are changed due to severe vibration during launching and changes of the in-orbit running environment, so that the final pointing measurement precision is affected.
In the prior art, a fixed star is mostly adopted as an observation target in a conventional method, and calibration parameters of a space direction measurement sensor during in-orbit work are calibrated by utilizing star-angular distance constraint. However, the fixed star observation quantity which can be obtained by the method is limited, and the requirement of calibrating the instrument parameter in the full view field range cannot be met; on the other hand, the conventional method only considers the focal length and the principal point of the calibrated space pointing measurement sensor, but does not modify the distortion coefficient which affects the pointing measurement precision; although lens distortion is considered in some methods, focal length, principal point correction and distortion correction are carried out separately, and the distortion correction adopts general polynomial fitting, and does not consider the actual rule of the lens distortion, so that the on-orbit calibration precision is difficult to guarantee, and the measurement precision of the space pointing measurement sensor is directly influenced.
Disclosure of Invention
The invention aims to overcome the defects in the prior art, provides a space-oriented measuring sensor calibration method based on astronomical technology, calibrates parameters such as focal length, principal point, radial distortion, eccentric distortion, image plane distortion and the like of a space-oriented measuring instrument by using landmark characteristics on the earth surface, and is suitable for practical engineering application.
The technical solution of the invention is as follows:
a space-pointing measuring instrument calibration method based on astronomical technology comprises the following steps:
step one, collecting an optical image of the earth surface, wherein the specific method comprises the following steps: on the set orbit height H, aligning a space pointing measuring instrument to the earth, so that the earth imaging covers the field range of the instrument as much as possible; setting an integration time to enable the texture of the earth surface to be imaged clearly;
step (II), landmark feature point extraction is carried out, and the specific method comprises the following steps: selecting N landmarks from the optical image, and extracting the image coordinate of the feature as (u)p,i,vp,i) Wherein p ═ 1,2, … N, represents the p-th landmark; (u)p,i,vp,i) Representing the coordinates of the image point of the ith feature of the p-th landmark on the image, up,iRepresenting the coordinates of the image point of the ith feature of the p-th landmark in the image line direction, vp,iThe coordinate of the image point of the ith feature of the pth landmark in the image column direction is represented, i is 1,2, … … k represents the ith feature point for describing the landmark, and k represents the number of the feature points of the pth landmark;
step three, carrying out landmark feature point matching, wherein the specific method comprises the following steps: matching the selected landmark with a terrestrial landmark database, and identifying the selected landmark feature point from the terrestrial landmark database;
step (IV), calculating the landmark feature point observation vector, wherein the specific method comprises the following steps: for each recognized landmark feature point, calculating a feature observation vector corresponding to the image coordinate describing the landmark feature point according to a measurement model of the spatial directional measuring instrument
Figure BDA0002980434310000021
An observation vector representing the ith feature point of the pth landmark;
step five, calculating a landmark feature point constraint vector, wherein the specific method comprises the following steps: for each identified landmark feature point, the right ascension alpha of the landmark in the earth landmark feature database is determinedp,iDeclination deltap,iComputing a feature constraint vector Vp,i;αp,iThe right ascension of the ith feature point of the pth landmark in the landmark feature database is represented, and the declination of the ith feature point of the pth landmark in the landmark feature database is represented, Vp,iRepresenting a feature constraint vector calculated for the ith feature point of the pth landmark;
step six, establishing a constraint equation and solving model parameters, wherein the specific method comprises the following steps: landmark feature point observation vector
Figure BDA0002980434310000022
And a constraint vector Vp,iIf the deviation exists, the parameter when the deviation reaches the minimum is the calibration coefficient of the space direction measuring instrument through an optimization method;
in the calibration method of the spatial direction measuring instrument based on the astronomical technology, in the step (one), the orbit height H should satisfy H ≦ D/Fov, wherein D represents the diameter of the earth, Fov represents the observation field of view of the spatial direction measuring instrument;
in the calibration method of the space-pointing measuring instrument based on the astronomical technology, in the step (II), the selected landmarks comprise artificial landmarks such as roads and the like and natural landmarks such as coastlines and the like, the number N of the landmarks is not less than 15, and the landmarks are uniformly distributed in the full view field range of the space-pointing measuring instrument;
in the above calibration method for a space-pointing measuring instrument based on astronomical techniques, in the step (four), the observation vectors of the landmark feature points
Figure BDA0002980434310000031
Calculated according to the following formula:
Figure BDA0002980434310000032
wherein the content of the first and second substances,
Figure BDA0002980434310000033
representing the image coordinates of the landmark feature point image of the ith feature point describing the p-th landmark after 3-order radial distortion correction, eccentric distortion and image plane distortion correction on the image, and representing:
Figure BDA0002980434310000034
rp,ithe distance of the ith feature point representing the p landmark from the principal point of the spatial orientation measuring instrument on the image is represented as:
rp,i=(xp,i-x0)2+(yp,i-y0)2
f is the focal length of the space direction measuring sensor, and the unit is mm;
(x0,y0) The intersection point of the main optical axis of the spatial direction measuring sensor and the image surface of the image is a main point in mm;
(xp,i,yp,i) Representing the coordinate position of the ith characteristic point of the p landmark imaged on the image plane by the space direction measuring instrument, and the image coordinate (u) of the ith characteristic point of the p landmarkp,i,vp,i) And calculating to obtain:
Figure BDA0002980434310000041
m and n respectively represent the resolution of the image sensor of the space direction measuring sensor in the row direction and the column direction;
ps represents the pixel size of the image sensor of the space direction measuring sensor;
k1represents a first order radial distortion coefficient;
k2representing a second order radial distortion coefficient;
k3represents the third order radial distortion coefficient;
p1representing a first off-center distortion factor;
p2representing a second off-center distortion factor;
ap1representing a first image plane distortion coefficient;
ap2representing a second image plane distortion coefficient;
in the above calibration method for a spatial direction measuring instrument based on astronomical techniques, in the step (V), the constraint vector V of the ith feature point of the pth landmark isp,iCalculated according to the following formula,
Figure BDA0002980434310000042
wherein the content of the first and second substances,
p,ip,i) The right ascension and the declination of the ith characteristic point of the pth landmark are represented;
Figure BDA0002980434310000051
representing an attitude matrix when the space pointing measurement sensor shoots an earth surface image;
Figure BDA0002980434310000052
kappa and omega respectively represent the rotation angles of three axial directions of the measurement body system of the measurement sensor around the space direction, and unit radian;
in the above calibration method for a space-pointing measuring instrument based on astronomical techniques, in the sixth step, the observation vector is observed according to the landmark feature points
Figure BDA0002980434310000053
And a constraint vector Vp,iThe deviation of (d) yields the error equation as follows:
Figure BDA0002980434310000054
wherein the content of the first and second substances,
Figure BDA0002980434310000055
Figure BDA0002980434310000056
Figure BDA0002980434310000057
represents an approximation of Z, Δ Z represents the correction of Z;
Figure BDA0002980434310000058
a partial derivative matrix;
when identifying N landmarks in the field of view, the column
Figure BDA0002980434310000059
Equation, iteratively calculating x according to the least squares principle0,y0,f,k1,k2,k3,p1,p2,ap1,ap2,
Figure BDA00029804343100000510
ω,κ。
Compared with the prior art, the invention has the beneficial effects that:
(1) the on-orbit calibration of the spatial direction measuring instrument is realized by utilizing the abundant landmark features on the earth surface, and the calibration parameters can be ensured to be suitable for the whole imaging range, so that the on-orbit calibration of the spatial direction measuring instrument is ensured;
(2) the calibration model considers the influence of 3-order radial distortion, eccentric distortion and image plane distortion of the lens on the pointing measurement precision, the distortion correction model is consistent with the geometric optics principle and is closer to the imaging rule of the actual lens, and the calibration precision of the spatial pointing measurement instrument can be effectively improved;
(3) in the calibration resolving process, the focal length, the principal point and the 7 distortion coefficients are corrected simultaneously, so that the influence of step-by-step correction on the calibration precision is overcome, and the calibration precision of the spatial direction measuring instrument can be effectively improved.
Drawings
FIG. 1 is a flow chart of the present invention.
Detailed Description
The invention will be further described with reference to the accompanying drawings.
Step 1, collecting an optical image on the surface of the earth; on a certain orbit height, the space pointing measuring instrument is aligned to the earth, so that the earth imaging covers the field range of the instrument as much as possible; setting a proper integration time to enable the texture of the earth surface to be imaged clearly; in the present example, a spatial pointing measurement instrument with a field of view Fov of 20 ° x 20 ° was calibrated: the space direction measuring instrument shoots an optical image of the earth surface at the height H of 36000km of a track.
Step 2, extracting landmark feature points; selecting N landmarks from the optical image, and extracting the image coordinates (u) of each landmark feature pointp,i,vp,i) Where p ═ 1,2, … N denotes the p-th landmark, i ═ 1,2, … … k denotes the i-th feature point describing the landmark; (u)p,i,vp,i) Representing the coordinates of the image point of the ith feature of the p-th landmark on the image, up,iRepresenting the coordinates of the image point of the ith feature of the p-th landmark in the image line direction, vp,iRepresenting the pixel coordinates of the ith feature of the pth landmark in the image column direction; the selected landmarks comprise artificial landmarks such as roads and the like and natural landmarks such as coastlines and the like, the number N of the landmarks is 30, for each landmark, the number of feature points describing the landmark is not less than 10, and the 30 landmarks are uniformly distributed in the field range of 20 degrees multiplied by 20 degrees of a spatial directional measuring instrument.
Step 3, landmark feature point matching; matching the selected 30 landmark feature points with an earth landmark database, and identifying the selected 30 landmark feature points from the landmark database;
step 4, calculating a landmark feature point observation vector; for 30 recognized landmarks, calculating feature observation vectors corresponding to image coordinates describing landmark feature points according to a measurement model of a space-oriented measuring instrument
Figure BDA0002980434310000061
In the embodiment of the present invention, the resolution of the space-pointing measurement sensor in the row and column directions is 2048, m is 2048, and the pixel size ps is 0.0055mm, so that each feature observation vector can be written as:
Figure BDA0002980434310000062
wherein the content of the first and second substances,
Figure BDA0002980434310000071
the coordinate of a landmark feature point image of the ith feature describing the p-th landmark after 3-order radial distortion correction, eccentric distortion and image plane distortion correction on the image can be represented as:
Figure BDA0002980434310000072
rp,ithe distance of the ith feature point representing the pth landmark from the principal point of the spatial pointing instrument on the image can be represented as:
rp,i=(xp,i-x0)2+(yp,i-y0)2
f is the focal length of the space direction measuring sensor, and the unit is mm;
(x0,y0) The intersection point of the main optical axis of the spatial direction measuring sensor and the image surface of the image is a main point in mm;
(xp,i,yp,i) The coordinate position of the ith characteristic point of the pth landmark imaged on the image plane by the space direction measuring instrument is represented, and the image coordinate (u) of the ith characteristic point of the pth landmark can be obtainedp,i,vp,i) And calculating to obtain:
Figure BDA0002980434310000073
k1represents a first order radial distortion coefficient;
k2representing a second order radial distortion coefficient;
k3represents the third order radial distortion coefficient;
p1representing a first off-center distortion factor;
p2representing a second off-center distortion factor;
ap1representing the first image planeA surface distortion coefficient;
ap2representing a second image plane distortion coefficient;
step 5, calculating landmark feature point constraint vectors; for each identified landmark, the right ascension alpha of the landmark feature point in the earth landmark feature database is determinedp,iDeclination deltap,iComputing a feature constraint vector Vp,i
Figure BDA0002980434310000081
Wherein the content of the first and second substances,
p,ip,i) The right ascension and the declination of the ith characteristic point of the pth landmark are represented;
Figure BDA0002980434310000082
representing an attitude matrix when the space pointing measurement sensor shoots an earth surface image;
Figure BDA0002980434310000083
kappa and omega respectively represent the rotation angles of three axial directions of the measurement body system of the measurement sensor around the space direction, and unit radian;
step 6, establishing a constraint equation and solving model parameters; the landmark feature point observation vector and the constraint vector have deviation, and the parameter when the deviation reaches the minimum is the calibration coefficient of the space-pointing measuring instrument by an optimization method; the error equation is as follows:
Figure BDA0002980434310000084
wherein the content of the first and second substances,
Figure BDA0002980434310000085
Figure BDA0002980434310000086
Figure BDA0002980434310000087
indicating the approximate value of Z, Δ Z indicating the correction of Z
Figure BDA0002980434310000088
Is Gp,i(x0,y0,f,k1,k2,k3,p1,p2,ap1,ap2,
Figure BDA0002980434310000089
ω, κ) partial derivative matrix
In the terrestrial landmark image of the earth surface shot by the space-oriented measuring sensor, 30 landmarks are identified, the total number of feature points is 420, 3 × 420 equations can be listed, and the in-orbit calibration parameters of the space-oriented measuring sensor are obtained through iterative calculation according to the least square principle:
Figure BDA0002980434310000091
those skilled in the art will appreciate that the details not described in the present specification are well known.

Claims (7)

1. A space direction measuring instrument calibration method based on astronomical technology is characterized by comprising the following steps:
step one, collecting an optical image of the surface of the earth;
step (II), landmark feature point extraction is carried out, and the specific method comprises the following steps: selecting N landmarks from the optical image, and extracting the image coordinate of the feature as (u)p,i,vp,i) Wherein p ═ 1,2, … N, represents the p-th landmark; (u)p,i,vp,i) Representing the coordinates of the image point of the ith feature of the p-th landmark on the image, up,iIndicating the ith feature of the pth landmarkCoordinates of image points in the direction of image lines, vp,iThe coordinate of the image point of the ith feature of the pth landmark in the image column direction is represented, i is 1,2, … … k represents the ith feature point for describing the landmark, and k represents the number of the feature points of the pth landmark;
step three, carrying out landmark feature point matching, wherein the specific method comprises the following steps: matching the selected landmark with a terrestrial landmark database, and identifying the selected landmark feature point from the terrestrial landmark database;
step (IV), calculating the landmark feature point observation vector, wherein the specific method comprises the following steps: for each recognized landmark feature point, calculating a feature observation vector corresponding to the image coordinate describing the landmark feature point according to a measurement model of the spatial directional measuring instrument
Figure FDA0002980434300000011
An observation vector representing the ith feature point of the pth landmark;
step five, calculating a landmark feature point constraint vector, wherein the specific method comprises the following steps: for each identified landmark feature point, the right ascension alpha of the landmark in the earth landmark feature database is determinedp,iDeclination deltap,iComputing a feature constraint vector Vp,i;αp,iThe right ascension of the ith feature point of the pth landmark in the landmark feature database is represented, and the declination of the ith feature point of the pth landmark in the landmark feature database is represented, Vp,iRepresenting a feature constraint vector calculated for the ith feature point of the pth landmark;
step six, establishing a constraint equation and solving model parameters, wherein the specific method comprises the following steps: landmark feature point observation vector
Figure FDA0002980434300000012
And a constraint vector Vp,iAnd if the deviation exists, the parameter when the deviation reaches the minimum is the calibration coefficient of the space direction measuring instrument through an optimization method.
2. The astronomical-technology-based spatial direction measurement instrument calibration method according to claim 1, wherein: in the step (one), the specific method for collecting the optical image on the surface of the earth comprises the following steps: on the set orbit height H, aligning a space pointing measuring instrument to the earth, so that the earth imaging covers the field range of the instrument as much as possible; the integration time is set so that the texture of the earth's surface can be imaged sharply.
3. The astronomical-technology-based spatial direction measurement instrument calibration method according to claim 2, wherein: the orbit height H satisfies H ≦ D/Fov, where D represents the diameter of the earth and Fov represents the observation field of view of the spatial pointing instrument.
4. The astronomical-technology-based spatial direction measurement instrument calibration method according to claim 1, wherein: in the step (II), the landmarks comprise artificial landmarks and coastline natural landmarks, the number N of the landmarks is more than or equal to 15, and the landmarks are uniformly distributed in the full-field range of the space-pointing measuring instrument.
5. The astronomical-technology-based spatial direction measurement instrument calibration method according to claim 1, wherein: in the step (IV), the observation vector of the landmark feature point
Figure FDA0002980434300000021
Calculated according to the following formula:
Figure FDA0002980434300000022
wherein the content of the first and second substances,
Figure FDA0002980434300000023
representing the image coordinates of the landmark feature point image of the ith feature point describing the p-th landmark after 3-order radial distortion correction, eccentric distortion and image plane distortion correction on the image, and representing:
Figure FDA0002980434300000024
rp,ithe distance of the ith feature point representing the p landmark from the principal point of the spatial orientation measuring instrument on the image is represented as:
rp,i=(xp,i-x0)2+(yp,i-y0)2
f is the focal length of the space direction measuring sensor, and the unit is mm;
(x0,y0) The intersection point of the main optical axis of the spatial direction measuring sensor and the image surface of the image is a main point in mm;
(xp,i,yp,i) Representing the coordinate position of the ith characteristic point of the p landmark imaged on the image plane by the space direction measuring instrument, and the image coordinate (u) of the ith characteristic point of the p landmarkp,i,vp,i) And calculating to obtain:
Figure FDA0002980434300000031
m and n respectively represent the resolution of the image sensor of the space direction measuring sensor in the row direction and the column direction;
ps represents the pixel size of the image sensor of the space direction measuring sensor;
k1represents a first order radial distortion coefficient;
k2representing a second order radial distortion coefficient;
k3represents the third order radial distortion coefficient;
p1representing a first off-center distortion factor;
p2representing a second off-center distortion factor;
ap1representing a first image plane distortion coefficient;
ap2representing the second image plane distortion coefficient.
6. The method of claim 1The space direction measuring instrument calibration method based on the astronomical technology is characterized by comprising the following steps: in the step (V), the constraint vector V of the ith feature point of the pth landmarkp,iCalculated according to the following formula,
Figure FDA0002980434300000041
wherein the content of the first and second substances,
p,ip,i) The right ascension and the declination of the ith characteristic point of the pth landmark are represented;
Figure FDA0002980434300000042
representing an attitude matrix when the space pointing measurement sensor shoots an earth surface image;
Figure FDA0002980434300000043
and kappa and omega respectively represent the rotation angles of three axial directions of the system measured by the space direction measuring sensor, and the unit radian.
7. The astronomical-technology-based spatial direction measurement instrument calibration method according to claim 1, wherein: in the step (VI), observing vectors according to the landmark feature points
Figure FDA0002980434300000044
And a constraint vector Vp,iThe deviation of (d) yields the error equation as follows:
Figure FDA0002980434300000045
wherein the content of the first and second substances,
Figure FDA0002980434300000046
Figure FDA0002980434300000047
Figure FDA0002980434300000048
represents an approximation of Z, Δ Z represents the correction of Z;
Figure FDA0002980434300000049
is composed of
Figure FDA00029804343000000410
A partial derivative matrix;
when identifying N landmarks in the field of view, the column
Figure FDA0002980434300000051
Equation, iteratively calculating x according to the least squares principle0,y0,f,k1,k2,k3,p1,p2,ap1,ap2,
Figure FDA0002980434300000052
ω,κ。
CN202110285856.8A 2021-03-17 2021-03-17 Astronomical technology-based spatial direction measuring instrument calibration method Active CN113034610B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110285856.8A CN113034610B (en) 2021-03-17 2021-03-17 Astronomical technology-based spatial direction measuring instrument calibration method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110285856.8A CN113034610B (en) 2021-03-17 2021-03-17 Astronomical technology-based spatial direction measuring instrument calibration method

Publications (2)

Publication Number Publication Date
CN113034610A true CN113034610A (en) 2021-06-25
CN113034610B CN113034610B (en) 2022-06-03

Family

ID=76471289

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110285856.8A Active CN113034610B (en) 2021-03-17 2021-03-17 Astronomical technology-based spatial direction measuring instrument calibration method

Country Status (1)

Country Link
CN (1) CN113034610B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103234556A (en) * 2013-04-18 2013-08-07 哈尔滨工业大学 In-orbit calibration method of star sensor lens distortion and satellite attitude determination method based on starlight vector correction
CN105548976A (en) * 2015-12-14 2016-05-04 中国科学院长春光学精密机械与物理研究所 Shipborne radar offshore precision identification method
CN107967700A (en) * 2017-12-01 2018-04-27 北京控制工程研究所 The in-orbit geometric correction of the wide working distance binocular camera of big visual field and precision test method
CN108444449A (en) * 2018-02-02 2018-08-24 中国科学院西安光学精密机械研究所 It is a kind of to the object space attitude measurement method with parallel lines feature
CN108592945A (en) * 2018-03-27 2018-09-28 中国人民解放军国防科技大学 Online calibration method for errors of inertia/astronomical combination system
CN112254743A (en) * 2020-10-15 2021-01-22 长春工业大学 Star sensor on-orbit calibration method based on star angular distance subtraction

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103234556A (en) * 2013-04-18 2013-08-07 哈尔滨工业大学 In-orbit calibration method of star sensor lens distortion and satellite attitude determination method based on starlight vector correction
CN105548976A (en) * 2015-12-14 2016-05-04 中国科学院长春光学精密机械与物理研究所 Shipborne radar offshore precision identification method
CN107967700A (en) * 2017-12-01 2018-04-27 北京控制工程研究所 The in-orbit geometric correction of the wide working distance binocular camera of big visual field and precision test method
CN108444449A (en) * 2018-02-02 2018-08-24 中国科学院西安光学精密机械研究所 It is a kind of to the object space attitude measurement method with parallel lines feature
CN108592945A (en) * 2018-03-27 2018-09-28 中国人民解放军国防科技大学 Online calibration method for errors of inertia/astronomical combination system
CN112254743A (en) * 2020-10-15 2021-01-22 长春工业大学 Star sensor on-orbit calibration method based on star angular distance subtraction

Non-Patent Citations (11)

* Cited by examiner, † Cited by third party
Title
GANG CHEN等: "The Comparison of Two Stereo Vision Sensor Calibration Methods", 《ADVANCED MATERIALS RESEARCH》 *
LI WANG等: "Researches on hazard avoidance cameras calibration of Lunar Rover", 《INTERNATIONAL CONFERENCE ON SPACE OPTICS》 *
PAUL FURGALE等: "Unified Temporal and Spatial Calibration for Multi-Sensor Systems", 《2013 IEEE/RSJ INTERNATIONAL CONFERENCE ON INTELLIGENT ROBOTS AND SYSTEMS (IROS)》 *
刘延芳等: "星上敏感器空间矢量自准直测量方法及标定试验", 《机械工程学报》 *
吴奋陟等: "基于标定场的激光雷达两步标定方法", 《空间控制技术与应用》 *
张成渝等: "小行星探测下降着陆段光学导航技术发展综述", 《红外与激光工程》 *
李建国等: "基于特征匹配的着陆导航敏感器相对校准算法", 《红外与激光工程》 *
牟伟清等: "针对射频相对测量敏感器的快速标定方法", 《北京理工大学学报》 *
王立等: "我国深空探测光学敏感器技术发展与应用", 《红外与激光工程》 *
陈刚等: "空间点的立体视觉传感器标定方法", 《光学精密工程》 *
雷拥军: "星敏感器姿态测量相对基准偏差在轨标方法研究", 《空间控制技术与应用》 *

Also Published As

Publication number Publication date
CN113034610B (en) 2022-06-03

Similar Documents

Publication Publication Date Title
CN106871927B (en) Unmanned aerial vehicle photoelectric pod installation error calibration method
CN101246590B (en) Star loaded camera spacing deformation image geometric correction method
CN111798523B (en) Satellite camera on-orbit calibration attitude determination and remote sensing image geometric positioning method and system
CN107564057B (en) High-orbit planar array optical satellite in-orbit geometric calibration method considering atmospheric refraction correction
CN109269512B (en) Relative navigation method for fusing planet landing image and distance measurement
CN111486864B (en) Multi-source sensor combined calibration method based on three-dimensional regular octagon structure
CN109613583A (en) The passive object localization method of the time difference is surveyed based on the direction finding of Dan Xingyu earth station and joint
CN112435301A (en) Remote sensing camera on-orbit geometric calibration method based on star locus
CN113029132A (en) Spacecraft navigation method combining ground image and astrolabe measurement
CN110660099B (en) Rational function model fitting method for remote sensing image processing based on neural network
CN113034610B (en) Astronomical technology-based spatial direction measuring instrument calibration method
CN111508028A (en) Autonomous in-orbit geometric calibration method and system for optical stereo mapping satellite camera
CN114972078A (en) Method and system for improving uncontrolled geometric quality of domestic optical satellite image by applying SAR image
CN113514055B (en) Joint estimation method for atmospheric refraction and ground attitude of foundation star sensor
CN111044076B (en) Geometric calibration method for high-resolution first-number B satellite based on reference base map
CN115388890A (en) Visual sense-based multi-unmanned aerial vehicle cooperative ground target positioning method
CN115015931A (en) Real-time differential stereo SAR geometric positioning method and system without external error correction
CN109143295B (en) Internal orientation element calibration method combining digitized geometric calibration field and GCP
CN108106634A (en) Star sensor internal parameter calibration method for direct star observation
Zhan et al. Adaptive celestial positioning for the stationary mars rover based on a self-calibration model for the star sensor
CN113049006A (en) Star light vector gas masking difference correction method based on model parameter estimation
CN112013874A (en) Satellite-borne laser altimeter on-orbit calibration method based on laser footprint prior coordinate
CN113900125B (en) Satellite-ground combined linear array imaging remote sensing satellite full-autonomous geometric calibration method and system
CN114858186B (en) On-satellite geometric calibration method for linear array camera under fixed star observation mode
CN116091546B (en) Observation construction method under push-broom mode of optical camera

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant