CN111141260B - Satellite resolution detection method and system based on three-dimensional spherical target - Google Patents
Satellite resolution detection method and system based on three-dimensional spherical target Download PDFInfo
- Publication number
- CN111141260B CN111141260B CN201811308362.1A CN201811308362A CN111141260B CN 111141260 B CN111141260 B CN 111141260B CN 201811308362 A CN201811308362 A CN 201811308362A CN 111141260 B CN111141260 B CN 111141260B
- Authority
- CN
- China
- Prior art keywords
- points
- coordinates
- spherical target
- dimensional spherical
- meridian
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Image Analysis (AREA)
Abstract
The invention provides a satellite resolution detection method and a satellite resolution detection system based on a three-dimensional spherical target, wherein the detection method comprises the following steps: determining the geometric center coordinates of the projection curve of each meridian on the three-dimensional spherical target on the satellite image; at least three identification points are taken on each warp of the three-dimensional spherical target, and the coordinates of projection points of the identification points on the satellite image are determined; determining demarcation points on each projection curve; according to the cross ratio invariant property in the photographic theory, calculating to obtain the three-dimensional coordinates of the corresponding points of the demarcation points on the three-dimensional spherical target; and calculating the distance between corresponding points of the demarcation points on the adjacent projection curves on the three-dimensional spherical target, and taking the distance as the geometric resolution of the satellite image at the corresponding position. According to the technical scheme provided by the invention, the three-dimensional spherical target is used for replacing the traditional two-dimensional plane three-line target and the plane radial target, so that the problem of inaccurate detection result caused by the shooting angle during the detection of the satellite resolution can be solved.
Description
Technical Field
The invention belongs to the technical field of photogrammetry and remote sensing mapping, and particularly relates to a satellite resolution detection method and system based on a three-dimensional spherical target.
Background
The geometric resolution of the high-resolution remote sensing satellite image is the primary evaluation index of the imaging quality of the satellite sensor. At present, the optical imaging system in the field of aerospace remote sensing is mainly a photoelectric imaging system, and the geometric resolution generally refers to the ground resolution (GSD) corresponding to a single pixel on a sensor. The GSD is the size of a pixel projected onto the ground via an optical path, and for a particular on-orbit sensor, the size is related to the spacecraft orbit height, the sensor side view angle, and other factors.
In order to detect the geometric resolution of the on-orbit running sensor, a three-line target of a two-dimensional plane is mainly used in the early stage, has a simple structure, is easy to manufacture, is convenient to arrange, is visual and reliable, and is widely used for detecting the ground resolution of a film optical imaging satellite by countries around the world. After the three-line target is applied to the photoelectric imaging system, due to the fact that the number of lines is small, the length is short, the moire effect, the aliasing effect, the edge effect and the like cause distortion of an obtained image, even serious distortion, and therefore the pixel resolution of the photoelectric imaging system can not be obtained or correct results can not be obtained by using the three-line target to check. In order to avoid the error of the three-line target detection result, the pixel resolution of a radial target detection photoelectric imaging system is widely adopted in various countries at the present stage. The number of the radial target lines is large (more than tens of radial target lines), so that the accidental detection result with small number of the three-line target lines is effectively avoided. And because the line width becomes the trapezoidal continuous gradual change, can reduce the influence of factors such as phase difference effectively, can improve the accuracy that detects. In addition, the radiation angle is generally larger than 90 degrees, the full range from the CCD linear array direction to the flight direction is covered, the image quality in different directions can be detected, and the interference of the possibility of pseudo resolution is eliminated.
The patent document of the invention with the publication number of CN104764465B discloses a method for measuring the ground resolution of a remote sensing satellite, wherein a sector radial target, a rectangular target and a black substrate are laid, then a satellite to be measured is used for imaging the sector radial target and the rectangular target area to obtain the target edge position of the target image, the circle center coordinate of the sector radial target is calculated, the coordinate of the separable position is determined, and the ground resolution is obtained through calculation.
The two-dimensional radial target has more advantages, but because the space sensor generally has a certain inclination angle when observing the ground, the two-dimensional plane target image generates affine deformation, and therefore, certain errors still exist by utilizing the geometric resolution of the two-dimensional radial target for detecting.
Disclosure of Invention
The invention aims to provide a satellite resolution detection method based on a three-dimensional spherical target, which is used for solving the problem that the detection result is inaccurate due to the existence of a shooting angle when the satellite resolution is detected in the prior art; correspondingly, the invention also provides a satellite resolution detection system based on the three-dimensional spherical target, which is used for solving the problem that the detection result is inaccurate due to the existence of a shooting angle when the satellite resolution is detected in the prior art.
In order to achieve the purpose, the technical scheme provided by the invention is as follows:
a satellite resolution detection method based on a three-dimensional spherical target comprises the following steps:
determining the geometric center coordinates of the projection curve of each meridian on the three-dimensional spherical target on the satellite image;
at least three identification points are taken on each warp of the three-dimensional spherical target, and the coordinates of projection points of the identification points on the satellite image are determined;
determining demarcation points on each projection curve; the demarcation point is used for distinguishing a fuzzy region and a clear region of the corresponding projection curve;
calculating to obtain the three-dimensional coordinates of the corresponding points of the demarcation points on the three-dimensional spherical target according to the coordinates of the identification points, the coordinates of the geometric center, the coordinates of the projection points, the coordinates of the demarcation points and the coordinates of the sphere center of the three-dimensional spherical target in combination with the constant cross ratio property in the photography theory;
and calculating the distance between corresponding points of the demarcation points on the adjacent projection curves on the three-dimensional spherical target, and taking the distance as the geometric resolution of the satellite image at the corresponding position.
According to the technical scheme provided by the invention, the three-dimensional spherical target is used for replacing the traditional two-dimensional plane three-line target and the plane radial target, so that the sight vectors of the satellite sensor pixels at different side visual angles can be vertically intersected with the three-dimensional spherical target, the error of the plane target in detecting the geometric resolution of the satellite sensor is effectively reduced, and the problem of inaccurate detection result caused by the existence of a shooting angle in the prior art when the satellite resolution is detected is solved.
Furthermore, the average value of the included angle between the photographing light and the spherical surface of the three-dimensional spherical target is set as the inclination angle corresponding to the geometric resolution of the satellite image, a change relation curve between the geometric resolution of the satellite image and the corresponding inclination angle is firstly drawn, and then the curve is fitted by utilizing a polynomial to obtain the relation between the geometric resolution of the satellite image and the inclination angle corresponding to the geometric resolution of the satellite image.
The invention establishes the relation between the satellite resolution and the shooting angle, and is beneficial to directly determining the resolution of the image according to the shooting angle, thereby quickly obtaining the actual size of the shot object and improving the working efficiency of obtaining the geometric resolution of the satellite image.
Further, the equation of the meridian projection curve is set as:
Ax2+By2+Cxy+Dx+Ey+F＝0
setting the coordinate of the ith pixel point on the meridian projection curve as (x)i，yi) Then, n sampling points are taken from the meridian projection curve and substituted into the above formula to obtain the following linear equation set:
solving the values of the parameters A, B, C, D, E and F by using a least square method;
let the geometric center of the meridian projection curve be (x)o，yo) And then:
according to the invention, the projection curve equation is calculated by adopting the least square method according to the coordinates of the points on the projection curve, the geometric center point of the projection curve can be accurately obtained, and a foundation is laid for subsequent calculation.
Further, three identification points on one meridian on the three-dimensional spherical target are respectively Pa，PbAnd PcThe three marksThe projection points of the identification points on the satellite images are respectively pa，pbAnd pc(ii) a Setting the dividing point of the projection curve corresponding to the meridian as plThe corresponding point of the dividing point on the meridian line is PLAccording to the formula:
calculating to obtain LLAnd LbAngle L betweenLOLbThen according to the coordinates of the centre of sphere, PaCoordinates of points, PbCoordinates of points and PcCalculating the coordinates of the points to obtain the three-dimensional coordinates of the demarcation points; l isL，La，Lb，LcFor three-dimensional spherical target sphere center and PL，Pa，Pb，PcDetermined straight line, /)l，la，lb，lcDistributed as the geometric center and point p of the corresponding projection curve of the meridianl，pa，pbAnd pcAnd (4) determining a straight line.
According to the invention, the three-dimensional coordinates of the corresponding points of the boundary points on the target are calculated and obtained by combining the boundary points on the satellite images according to the nature of the constant cross ratio in the photographic theory, thereby laying a foundation for calculating the satellite resolution.
Further, the three-dimensional coordinate of the corresponding point of the demarcation point on the three-dimensional spherical target on the projection curve corresponding to the ith meridian is set as (X)i，Yi，Zi) Then, the geometric resolution of the satellite image at the corresponding position is:
the method takes the distance between the corresponding points of the demarcation points on the two adjacent projection curves on the satellite image on the target as the satellite resolution at the corresponding position of the target, fully considers the influence of the fuzzy area and the clear area on the satellite image, and ensures that the calculation result of the satellite resolution is more accurate.
A satellite resolution detection system based on a three-dimensional spherical target, comprising a processor and a memory, the memory having stored thereon a computer program for execution on the processor; wherein the processor implements the following steps when executing the computer program:
determining the geometric center coordinates of the projection curve of each meridian on the three-dimensional spherical target on the satellite image;
at least three identification points are taken on each warp of the three-dimensional spherical target, and the coordinates of projection points of the identification points on the satellite image are determined;
determining demarcation points on each projection curve; the demarcation point is used for distinguishing a fuzzy region and a clear region of the corresponding projection curve;
calculating to obtain the three-dimensional coordinates of the corresponding points of the demarcation points on the three-dimensional spherical target according to the coordinates of the identification points, the coordinates of the geometric center, the coordinates of the projection points, the coordinates of the demarcation points and the coordinates of the sphere center of the three-dimensional spherical target in combination with the constant cross ratio property in the photography theory;
and calculating the distance between corresponding points of the demarcation points on the adjacent projection curves on the three-dimensional spherical target, and taking the distance as the geometric resolution of the satellite image at the corresponding position.
According to the technical scheme provided by the invention, the three-dimensional spherical target is used for replacing the traditional two-dimensional plane three-line target and the plane radial target, so that the sight vectors of the satellite sensor pixels at different side visual angles can be vertically intersected with the three-dimensional spherical target, the error of the plane target in detecting the geometric resolution of the satellite sensor is effectively reduced, and the problem of inaccurate detection result caused by the existence of a shooting angle in the prior art when the satellite resolution is detected is solved.
Furthermore, the average value of the included angle between the photographing light and the spherical surface of the three-dimensional spherical target is set as the inclination angle corresponding to the geometric resolution of the satellite image, a change relation curve between the geometric resolution of the satellite image and the corresponding inclination angle is firstly drawn, and then the curve is fitted by utilizing a polynomial to obtain the relation between the geometric resolution of the satellite image and the inclination angle corresponding to the geometric resolution of the satellite image.
The invention establishes the relation between the satellite resolution and the shooting angle, and is beneficial to directly determining the resolution of the image according to the shooting angle, thereby quickly obtaining the actual size of the shot object and improving the working efficiency of obtaining the geometric resolution of the satellite image.
Further, the equation of the meridian projection curve is set as:
Ax2+By2+Cxy+Dx+Ey+F＝0
setting the coordinate of the ith pixel point on the meridian projection curve as (x)i，yi) Then, n sampling points are taken from the meridian projection curve and substituted into the above formula to obtain the following linear equation set:
solving the values of the parameters A, B, C, D, E and F by using a least square method;
let the geometric center of the meridian projection curve be (x)o，yo) And then:
according to the invention, the projection curve equation is calculated by adopting the least square method according to the coordinates of the points on the projection curve, the geometric center point of the projection curve can be accurately obtained, and a foundation is laid for subsequent calculation.
Further, three identification points on one meridian on the three-dimensional spherical target are respectively Pa，PbAnd PcThe projection points of the three identification points on the satellite image are respectively pa，pbAnd pc(ii) a Setting the dividing point of the projection curve corresponding to the meridian as plThe corresponding point of the dividing point on the meridian line is PLAccording to the formula:
calculating to obtain LLAnd LbAngle L betweenLOLbThen according to the coordinates of the centre of sphere, PaCoordinates of points, PbCoordinates of points and PcCalculating the coordinates of the points to obtain the three-dimensional coordinates of the demarcation points; l isL，La，Lb，LcFor three-dimensional spherical target sphere center and PL，Pa，Pb，PcDetermined straight line, /)l，la，lb，lcDistributed as the geometric center and point p of the corresponding projection curve of the meridianl，pa，pbAnd pcAnd (4) determining a straight line.
According to the invention, the three-dimensional coordinates of the corresponding points of the boundary points on the target are calculated and obtained by combining the boundary points on the satellite images according to the nature of the constant cross ratio in the photographic theory, thereby laying a foundation for calculating the satellite resolution.
Further, the three-dimensional coordinate of the corresponding point of the demarcation point on the three-dimensional spherical target on the projection curve corresponding to the ith meridian is set as (X)i，Yi，Zi) Then, the geometric resolution of the satellite image at the corresponding position is:
the method takes the distance between the corresponding points of the demarcation points on the two adjacent projection curves on the satellite image on the target as the satellite resolution at the corresponding position of the target, fully considers the influence of the fuzzy area and the clear area on the satellite image, and ensures that the calculation result of the satellite resolution is more accurate.
Drawings
FIG. 1 is a top view of a three-dimensional spherical target in an embodiment of a method of the invention;
FIG. 2 is a front view of a three-dimensional spherical target in an embodiment of a method of the invention;
FIG. 3 is a side view of a three-dimensional spherical target in an embodiment of a method of the invention;
FIG. 4 is a schematic diagram of a straight line determined by the center of a three-dimensional spherical target sphere and corresponding points of each identification point and dividing point on a meridian in the embodiment of the method of the invention;
FIG. 5 is a schematic diagram of a straight line defined by the geometric center of the meridian projection, the demarcation point on the meridian projection curve and the corresponding point of each identification point in the embodiment of the method.
Detailed Description
The invention aims to provide a satellite resolution detection method based on a three-dimensional spherical target, which is used for solving the problem that the detection result is inaccurate due to the existence of a shooting angle when the satellite resolution is detected in the prior art; correspondingly, the invention also provides a satellite resolution detection system based on the three-dimensional spherical target, which is used for solving the problem that the detection result is inaccurate due to the existence of a shooting angle when the satellite resolution is detected in the prior art.
In order to achieve the purpose, the technical scheme provided by the invention is as follows:
a satellite resolution detection method based on a three-dimensional spherical target comprises the following steps:
determining the geometric center coordinates of the projection curve of each meridian on the three-dimensional spherical target on the satellite image;
at least three identification points are taken on each warp of the three-dimensional spherical target, and the coordinates of projection points of the identification points on the satellite image are determined;
determining demarcation points on each projection curve; the demarcation point is used for distinguishing a fuzzy region and a clear region of the corresponding projection curve;
calculating to obtain the three-dimensional coordinates of the corresponding points of the demarcation points on the three-dimensional spherical target according to the coordinates of the identification points, the coordinates of the geometric center, the coordinates of the projection points, the coordinates of the demarcation points and the coordinates of the sphere center of the three-dimensional spherical target in combination with the constant cross ratio property in the photography theory;
and calculating the distance between corresponding points of the demarcation points on the adjacent projection curves on the three-dimensional spherical target, and taking the distance as the geometric resolution of the satellite image at the corresponding position.
The following further describes embodiments of the present invention with reference to the drawings.
The method comprises the following steps:
the embodiment provides a satellite resolution detection method based on a three-dimensional spherical target, which is used for detecting the resolution of a satellite by adopting the three-dimensional spherical target, so that the accuracy of satellite resolution detection is improved.
The satellite resolution detection method based on the three-dimensional spherical target provided by the embodiment comprises the following steps:
determining the geometric center coordinates of the projection curve of each meridian on the three-dimensional spherical target on the satellite image;
at least three identification points are taken on each warp of the three-dimensional spherical target, and the coordinates of projection points of the identification points on the satellite image are determined;
determining demarcation points on each projection curve; the demarcation point is used for distinguishing a fuzzy region and a clear region of the corresponding projection curve;
calculating to obtain the three-dimensional coordinates of the corresponding points of the demarcation points on the three-dimensional spherical target according to the coordinates of the identification points, the coordinates of the geometric center, the coordinates of the projection points, the coordinates of the demarcation points and the coordinates of the sphere center of the three-dimensional spherical target in combination with the constant cross ratio property in the photography theory;
and calculating the distance between corresponding points of the demarcation points on the adjacent projection curves on the three-dimensional spherical target, and taking the distance as the geometric resolution of the satellite image at the corresponding position.
In this embodiment, the method for establishing the three-dimensional spherical target includes:
measuring and recording coordinates of bottom net edge point positions of the three-dimensional spherical target by using an RTK-GPS coordinate measuring instrument, a total station and the like, making point location identification by using an obvious mark, identifying the boundary of the target layout by using a nylon rope with small white elasticity, and uniformly arranging a white area and a black area with regular shapes on the surface of the three-dimensional spherical target; the three-dimensional spherical target in this embodiment is shown in fig. 1 in a top view, fig. 2 in a front view, and fig. 3 in a side view. On the three-dimensional spherical target, the vertical grid boundary lines of the spherical grid are called longitude lines, and the horizontal grid boundary lines are called latitude lines.
The method for determining the geometric center coordinates of the projection curve corresponding to the meridian projection on the three-dimensional spherical target in the embodiment comprises the following steps:
an equation of one of the meridian projection curves is set as
Ax2+By2+Cxy+Dx+Ey+F＝0
Setting the coordinate of the ith pixel point on the meridian projection curve as (x)i，yi) Then, n sampling points are taken from the meridian and substituted into the above formula to obtain the following linear equation set:
solving the values of the parameters A, B, C, D, E and F by using a least square method;
let the geometric center of the meridian projection curve be (x)o，yo) And then:
the geometric center of each meridian projection curve is solved by the method.
On a satellite image after a three-dimensional spherical target is shot, the closer to the center of the image, the easier the shot image is blurred due to insufficient resolution, so that a demarcation point for distinguishing a blurred area and a clear area needs to be determined on a meridian projection curve, the demarcation point can be determined manually, and the demarcation point can also be automatically identified through an image processing means; when the image processing means is adopted to automatically identify the demarcation point, a set resolution is firstly set, and when the resolution of a certain point is identified to be the set resolution, the point can be judged to be the demarcation point.
After determining a demarcation point used for distinguishing a fuzzy area from a clear area on a projection curve, calculating to obtain a three-dimensional coordinate corresponding to the demarcation point on the three-dimensional spherical target according to coordinates of each identification point, coordinates of a geometric center, coordinates of a projection point, coordinates of the demarcation point and coordinates of a sphere center of the three-dimensional spherical target and by combining the nature of invariant cross ratio in the photographic theory, taking one meridian as an example, the method specifically comprises the following steps:
the division point on the meridian projection curve is plPoints for distinguishing fuzzy and clear areas on the meridian projection curve, demarcation points p on the meridian projection curvelThe corresponding point on the longitude line is PLSetting the point as the dividing point of the meridian; three mark points P are taken on the longitude linea，PbAnd PcThe corresponding positions of the three identification points on the meridian projection curve are pa，pbAnd pcThen the sphere center of the three-dimensional spherical target and the point P on the meridian lineL、Pa、Pb、PcThe determined straight lines are respectively straight lines LL、La、LbAnd LcAs shown in fig. 4, the geometric center of the meridian projection curve and the point pl，pa，pbAnd pcThe determined straight lines are respectively straight lines ll，la，lb，lcAs shown in fig. 5;
straight line ll，la，lb，lcThe cross ratio of (A) to (B) is as follows:
straight line LL，La，Lb，LcThe cross ratio of (A) to (B) is as follows:
sin (A, B) represents the sine value of the included angle between the straight line A and the straight line B; based on the nature of the constant ratio between the common point and the straight line in the imaging geometry
From this, a straight line L can be calculatedLAnd a straight line LbBetweenAngle of less than LLOLbThen according to the sphere center coordinate P of the three-dimensional spherical targetaCoordinates of points, PbCoordinate sum P ofcCalculating the coordinates of the points to obtain the meridian demarcation point PLThree-dimensional coordinates of (a).
Solving the dividing point of each meridian according to the method, and setting the three-dimensional coordinate of the dividing point of the jth meridian as (X)j，Yj，Zj) Then the geometric resolution of the jth meridian image is
In this embodiment, the method for detecting a satellite resolution based on a three-dimensional spherical target further includes a step of determining a relationship between a geometric resolution of a satellite image and a corresponding tilt angle, specifically:
according to the three-dimensional position coordinates (X) of the satellite sensor at the right distinguishable boundary image point photographing time obtained by the position detection systemS，YS，ZS) And calculating the coordinates (X) of the center position of the satellite camera at the time of photographing by using the eccentricity vectors (dX, dY, dZ) between the photographing center of the camera and the center of the positioning system calibrated in the laboratoryC，YC，ZC) Comprises the following steps:
[XC，YC，ZC]T＝(XS，YS，ZS)T+(dX，dY，dZ)T
the coordinate (X) of the shooting center position of the satellite camera is determined according to the three-dimensional coordinates of each partitionable pointC，YC，ZC) And the sphere center coordinates (X) of the three-dimensional spherical targetO，YO，ZO) Calculating the included angle alpha between the photographic light and the spherical surface at each boundary pointj：
Taking included angle between photographic light at adjacent meridian demarcation point and three-dimensional spherical target sphereThe average value is the inclination angle corresponding to the geometric resolution of the image, and the inclination angle is set as thetaj-1And then:
θj-1＝(αj+αj-1)/2
geometric resolution res for drawing satellite imagej-1To the angle of inclination thetaj-1The mathematical expression between the geometric resolution of the satellite image and the corresponding inclination angle is obtained by utilizing a polynomial fitting curve.
The embodiment of the system is as follows:
the embodiment provides a satellite resolution detection system based on a three-dimensional spherical target, which comprises a processor and a memory, wherein the memory is used for storing a computer program executed on the processor, and when the processor executes the computer program, the satellite resolution detection method based on the three-dimensional spherical target provided in the above method embodiments is realized.
Claims (10)
1. A satellite resolution detection method based on a three-dimensional spherical target is characterized by comprising the following steps:
determining the geometric center coordinates of the projection curve of each meridian on the three-dimensional spherical target on the satellite image;
at least three identification points are taken on each warp of the three-dimensional spherical target, and the coordinates of projection points of the identification points on the satellite image are determined;
determining demarcation points on each projection curve; the demarcation point is used for distinguishing a fuzzy region and a clear region of the corresponding projection curve;
calculating to obtain the three-dimensional coordinates of the corresponding points of the demarcation points on the three-dimensional spherical target according to the coordinates of the identification points, the coordinates of the geometric center, the coordinates of the projection points, the coordinates of the demarcation points and the coordinates of the sphere center of the three-dimensional spherical target in combination with the constant cross ratio property in the photography theory;
and calculating the distance between corresponding points of the demarcation points on the adjacent projection curves on the three-dimensional spherical target, and taking the distance as the geometric resolution of the satellite image at the corresponding position.
2. The method as claimed in claim 1, wherein the average value of the angles between the photographing light and the spherical surface of the three-dimensional spherical target is defined as the tilt angle corresponding to the geometric resolution of the satellite image, and a curve of the variation relationship between the geometric resolution of the satellite image and the corresponding tilt angle is first drawn, and then a polynomial is used to fit the curve to obtain the relationship between the geometric resolution of the satellite image and the tilt angle corresponding to the geometric resolution of the satellite image.
3. The method for detecting the resolution of the satellite based on the three-dimensional spherical target according to claim 1, wherein the equation of the meridian projection curve is:
Ax2+By2+Cxy+Dx+Ey+F＝0
and (3) setting the coordinate of the ith pixel point on the meridian projection curve as (xi, yi), substituting n sampling points on the meridian projection curve into the formula to obtain the following linear equation set:
solving the values of the parameters A, B, C, D, E and F by using a least square method;
if the geometric center of the meridian projection curve is (xo, yo), then:
4. the method for detecting the resolution of a satellite based on a three-dimensional spherical target according to claim 1, wherein the three identification points on one meridian of the three-dimensional spherical target are respectively designated as Pa，PbAnd PcThe projection points of the three identification points on the satellite image are respectively pa，pbAnd pc(ii) a Setting the dividing point of the projection curve corresponding to the meridian as plThe dividing point corresponds to the point on the meridian linePLAccording to the formula:
calculating to obtain LLAnd LbAngle L betweenLOLbThen according to the coordinates of the centre of sphere, PaCoordinates of points, PbCoordinates of points and PcCalculating the coordinates of the points to obtain the three-dimensional coordinates of the demarcation points; l isL，La，Lb，LcFor three-dimensional spherical target sphere center and PL，Pa，Pb，PcDetermined straight line, /)l，la，lb，lcDistributed as the geometric center and point p of the corresponding projection curve of the meridianl，pa，pbAnd pcAnd (4) determining a straight line.
5. The method for detecting the resolution of a satellite based on a three-dimensional spherical target according to claim 1, wherein the three-dimensional coordinates of the corresponding point on the three-dimensional spherical target of the demarcation point on the projection curve corresponding to the ith meridian are (Xi, Yi, Zi), and the geometric resolution of the satellite image at the corresponding position is:
6. a satellite resolution detection system based on a three-dimensional spherical target, comprising a processor and a memory, the memory having stored thereon a computer program for execution on the processor; wherein the processor implements the following steps when executing the computer program:
determining the geometric center coordinates of the projection curve of each meridian on the three-dimensional spherical target on the satellite image;
at least three identification points are taken on each warp of the three-dimensional spherical target, and the coordinates of projection points of the identification points on the satellite image are determined;
determining demarcation points on each projection curve; the demarcation point is used for distinguishing a fuzzy region and a clear region of the corresponding projection curve;
calculating to obtain the three-dimensional coordinates of the corresponding points of the demarcation points on the three-dimensional spherical target according to the coordinates of the identification points, the coordinates of the geometric center, the coordinates of the projection points, the coordinates of the demarcation points and the coordinates of the sphere center of the three-dimensional spherical target in combination with the constant cross ratio property in the photography theory;
and calculating the distance between corresponding points of the demarcation points on the adjacent projection curves on the three-dimensional spherical target, and taking the distance as the geometric resolution of the satellite image at the corresponding position.
7. The system of claim 6, wherein an average value of an included angle between the photographing light and the spherical surface of the three-dimensional spherical target is defined as an inclination angle corresponding to the geometric resolution of the satellite image, and a curve of a variation relationship between the geometric resolution of the satellite image and the corresponding inclination angle is first drawn, and then a polynomial is used to fit the curve to obtain a relationship between the geometric resolution of the satellite image and the inclination angle corresponding to the geometric resolution of the satellite image.
8. The system for detecting the resolution of the satellite based on the three-dimensional spherical target according to claim 6, wherein the equation of the meridian projection curve is:
Ax2+By2+Cxy+Dx+Ey+F＝0
and (3) setting the coordinate of the ith pixel point on the meridian projection curve as (xi, yi), substituting n sampling points on the meridian projection curve into the formula to obtain the following linear equation set:
solving the values of the parameters A, B, C, D, E and F by using a least square method;
if the geometric center of the meridian projection curve is (xo, yo), then:
9. the system according to claim 6, wherein the three identification points on one meridian of the three-dimensional spherical target are respectively designated as Pa，PbAnd PcThe projection points of the three identification points on the satellite image are respectively pa，pbAnd pc(ii) a Setting the dividing point of the projection curve corresponding to the meridian as plThe corresponding point of the dividing point on the meridian line is PLAccording to the formula:
calculating to obtain LLAnd LbAngle L betweenLOLbThen according to the coordinates of the centre of sphere, PaCoordinates of points, PbCoordinates of points and PcCalculating the coordinates of the points to obtain the three-dimensional coordinates of the demarcation points; l isL，La，Lb，LcFor three-dimensional spherical target sphere center and PL，Pa，Pb，PcDetermined straight line, /)l，la，lb，lcDistributed as the geometric center and point p of the corresponding projection curve of the meridianl，pa，pbAnd pcAnd (4) determining a straight line.
10. The system according to claim 6, wherein the three-dimensional coordinates of the corresponding point on the three-dimensional spherical target at the demarcation point on the projection curve corresponding to the ith meridian are (Xi, Yi, Zi), and the geometric resolution of the satellite image at the corresponding position is:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811308362.1A CN111141260B (en) | 2018-11-05 | 2018-11-05 | Satellite resolution detection method and system based on three-dimensional spherical target |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811308362.1A CN111141260B (en) | 2018-11-05 | 2018-11-05 | Satellite resolution detection method and system based on three-dimensional spherical target |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111141260A CN111141260A (en) | 2020-05-12 |
| CN111141260B true CN111141260B (en) | 2021-03-26 |
Family
ID=70515645
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811308362.1A Active CN111141260B (en) | 2018-11-05 | 2018-11-05 | Satellite resolution detection method and system based on three-dimensional spherical target |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111141260B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN204612710U (en) * | 2015-01-22 | 2015-09-02 | 中国测绘科学研究院 | Be applicable to the removable control target ball device of aerial survey |
-
2018
- 2018-11-05 CN CN201811308362.1A patent/CN111141260B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN204612710U (en) * | 2015-01-22 | 2015-09-02 | 中国测绘科学研究院 | Be applicable to the removable control target ball device of aerial survey |
Non-Patent Citations (2)
| Title |
|---|
| 《基于辐射状靶标的天绘一号卫星CCD相机分辨率在轨检测》;彭宇;《测绘与空间地理信息》;20170731;第36卷(第7期);全文 * |
| 《航天遥感成像系统像元分辨率在轨检测方法研究》;高正清;《航天返回与遥感》;20090930;第30卷(第3期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111141260A (en) | 2020-05-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106441138B (en) | The deformation monitoring method of view-based access control model measurement | |
| CN110057295B (en) | Monocular vision plane distance measuring method without image control | |
| US10060739B2 (en) | Method for determining a position and orientation offset of a geodetic surveying device and such a surveying device | |
| CN109655040B (en) | Side slope displacement monitoring method based on unmanned aerial vehicle targeting technology | |
| US8970853B2 (en) | Three-dimensional measurement apparatus, three-dimensional measurement method, and storage medium | |
| CN108362712B (en) | Substrate mother board and detection method thereof | |
| CN102768762B (en) | Digital camera calibration method targeted to shield tunnel defect digital radiography detection and device thereof | |
| CN105823416A (en) | Method for measuring object through multiple cameras and device thereof | |
| US11346666B2 (en) | System and method for measuring a displacement of a mobile platform | |
| CN106871787A (en) | Large space line scanning imagery method for three-dimensional measurement | |
| US20130113897A1 (en) | Process and arrangement for determining the position of a measuring point in geometrical space | |
| US20210190484A1 (en) | Method and apparatus for single camera optical measurements | |
| CN106643567A (en) | Lane deviation system production line calibration board verification method and system | |
| CN110736447A (en) | method for checking horizontal position of integrated image acquisition equipment in vertical direction | |
| CN104930976A (en) | Portable crack length-measuring apparatus and method | |
| CN113465573A (en) | Monocular distance measuring method and device and intelligent device | |
| CN111141260B (en) | Satellite resolution detection method and system based on three-dimensional spherical target | |
| CN105758339B (en) | Optical axis and object plane measuring for verticality method based on geometric error correction technique | |
| CN116743973A (en) | Automatic correction method for noninductive projection image | |
| CN108286960B (en) | Focusing type light tube array device and photographic detection method | |
| CN113421300B (en) | Method and device for determining actual position of object in fisheye camera image | |
| CN113610782B (en) | Building deformation monitoring method, equipment and storage medium | |
| CN112734857B (en) | Calibration method for camera internal reference and camera relative laser radar external reference and electronic equipment | |
| US11640680B2 (en) | Imaging system and a method of calibrating an image system | |
| CN114663486A (en) | Building height measurement method and system based on binocular vision |
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 |