[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN109145835B - Rapid calculation method for observing specific region target by on-orbit satellite - Google Patents

Rapid calculation method for observing specific region target by on-orbit satellite Download PDF

Info

Publication number
CN109145835B
CN109145835B CN201810985165.7A CN201810985165A CN109145835B CN 109145835 B CN109145835 B CN 109145835B CN 201810985165 A CN201810985165 A CN 201810985165A CN 109145835 B CN109145835 B CN 109145835B
Authority
CN
China
Prior art keywords
target
area
polygon
satellite
vector
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.)
Expired - Fee Related
Application number
CN201810985165.7A
Other languages
Chinese (zh)
Other versions
CN109145835A (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.)
Unit 63789 Of Pla
Original Assignee
Unit 63789 Of Pla
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 Unit 63789 Of Pla filed Critical Unit 63789 Of Pla
Priority to CN201810985165.7A priority Critical patent/CN109145835B/en
Publication of CN109145835A publication Critical patent/CN109145835A/en
Application granted granted Critical
Publication of CN109145835B publication Critical patent/CN109145835B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00—Scenes; Scene-specific elements
    • G06V20/10—Terrestrial scenes
    • G06V20/13—Satellite images
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/96—Management of image or video recognition tasks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00—Scenes; Scene-specific elements
    • G06V20/10—Terrestrial scenes
    • G06V20/194—Terrestrial scenes using hyperspectral data, i.e. more or other wavelengths than RGB
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V2201/00—Indexing scheme relating to image or video recognition or understanding
    • G06V2201/07—Target detection

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Theoretical Computer Science (AREA)
  • Astronomy & Astrophysics (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention provides a rapid calculation method for observing a target in a specific area by an on-orbit satellite, which is characterized in that the target in the area is regarded as a plane, a polygonal plane closed graph is used for describing the regular or irregular shape of the target in the ground area, three criteria of an Euclidean length method, a vector distance method and a projection method are designed according to the geometric relationship between a polygon and an observation satellite beam, and the target in a satellite coverage area can be rapidly and accurately judged according to the collection of the three criteria. The method has the advantages of quick calculation and high calculation precision, can quickly make a satellite observation plan in real time, shortens the response time of task planning, and better conforms to the practical application of measurement and control engineering.

Description

Rapid calculation method for observing specific region target by on-orbit satellite
Technical Field
The invention belongs to the field of measurement and control management of on-orbit spacecrafts, and relates to a calculation method for observing a ground area target by an on-orbit satellite.
Background
With the development of hyperspectral remote sensing technology, earth observation satellites such as visible light, infrared light and microwave are widely applied to the fields of resource exploration, terrain mapping, disaster dynamic monitoring and the like. In order to realize dense coverage and multi-load observation of ground targets, the ground observation satellite can jointly observe the ground targets with various loads in a multi-satellite networking mode. The networking observation needs to make an observation plan quickly, shorten the response time of task planning, reasonably schedule resources and improve the observation efficiency, and the quick and accurate calculation of the coverage capability of an observation satellite on a ground area target is one of key technologies for improving the planning response. The methods used previously were: firstly, taking a ground area as a point target, calculating the coverage capability of a satellite on a ground observation point, wherein a calculation result has a large error with the actual observation coverage time; secondly, the regional target is regarded as a point set, and the coverage set of the beam pair point set is traversed and solved, so that the method has huge calculation amount and is one of factors for restricting quick response and planning.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention provides a method for rapidly judging the coverage of a polygonal area by an observation satellite, wherein the ground area is described by a polygonal closed graph, so that the calculation of the target observation coverage of the ground area is realized, and the method has the advantages of high calculation speed and high precision.
The technical scheme adopted by the invention for solving the technical problem comprises the following steps:
1) for a given regional target, the regional shape is represented by an N-polygon, wherein N is a positive integer greater than 1; the geodetic coordinates of the N vertices of the N-polygon are known as (L)i,Bi,Hi) 1,2,3, 9, N, wherein L is geodetic longitude, B is geodetic latitude, and H is elevation; the geodetic coordinate of the load beam center point E is (L)e,Be,He) The earth fixation position vector is
Figure BDA0001779534350000011
Calculating the position vector of N vertexes of the N-polygon in the earth fixation system as
Figure BDA0001779534350000012
Wherein,
Figure BDA0001779534350000013
REsemi-major axis, R, representing the equatorial ellipsoid of the earthPRepresents the polar radius of the earth;
Figure BDA0001779534350000021
2) the criterion of the target in the satellite coverage area by adopting the Euclidean length method is
Figure BDA0001779534350000022
Wherein FE is the covering capacity of the load in the ground target direction, and the position vectors of N vertexes of the N-polygon
Figure BDA0001779534350000023
Position vector with point E
Figure BDA0001779534350000024
3) The criterion of the satellite coverage area target by adopting the vector distance method is that ([ U ] d)i) FE ≦ wherein the vector
Figure BDA0001779534350000025
In the vector
Figure BDA0001779534350000026
Projection of
Figure BDA0001779534350000027
When N is more than 2, the position vector of N sides of the N-edge-shaped area in the earth fixation system
Figure BDA0001779534350000028
4) When a projection method is adopted to judge a target in a satellite coverage area, firstly, set primary selection is carried out, and if (| | L |) does not satisfy (| L | computationally branch cells0≤Le≤||L||∞)∩(||B||0≤Be≤||B||∞) The central point E of the load beam is directly discarded, | | L | | non-woven phosphor∞＝max(Li)，||B||∞＝max(Bi)，||L||0＝min(Li)，||B||0＝min(Bi) (ii) a If the load beam center point E is met, further judging that the E point and two adjacent vertexes of the N-polygon form a triangle Se,i,i-1Projected area of
Figure BDA0001779534350000029
Whether or not to satisfy
Figure BDA00017795343500000210
Wherein S is the projection area of the N-shaped polygon, and if the projection area meets the requirement, the N-shaped polygon has covering capability;
5) for a given regional target, one of the criteria in steps 2), 3) and 4) is satisfied, namely that the observation satellite has the coverage capability for the regional target.
The invention has the beneficial effects that: the method has the advantages that the ground area is represented by the polygon, the coverage of the observation satellite on the ground area target is judged through the vector relation between the load beam center point and the polygon, the judgment criterion is concise and clear, the constraint condition is few, the overlarge error between the calculation result of taking the ground area as the point target and the observation time is effectively avoided, and the problems of overlarge calculation amount and poor real-time performance when the point set is used for representing the area target in a traversing mode are solved. The judgment criterion has the advantages of fast calculation and high calculation precision, can quickly make a satellite observation plan in real time, shortens the response time of task planning, and is more in line with the practical application of measurement and control engineering.
Drawings
FIG. 1 is a flow chart of the method calculation.
Fig. 2 is a schematic diagram of vector coverage in a pentagonal region by the euclidean length method.
FIG. 3 is a schematic diagram of vector coverage in a pentagonal region by a vector distance method.
FIG. 4 is a schematic view of vector coverage in a pentagonal region by a projection method.
Detailed Description
The present invention will be further described with reference to the following drawings and examples, which include, but are not limited to, the following examples.
The method comprises the steps of considering the area target as a plane, describing the regular or irregular shape of the ground area target by a polygonal plane closed graph, designing three criteria of an Euclidean length method, a vector distance method and a projection method according to the geometrical relationship between a polygon and an observation satellite beam, and quickly and accurately judging the satellite coverage area target according to the collection of the three criteria.
The specific calculation flow is as follows:
for a given region target S, its region shape is represented by an N (N ═ 2.., N, N being a positive integer) polygon. The geodetic coordinate of the N-vertex geodetic survey of the N-polygon is known as (L)i,Bi,Hi) 1,2,3, 9, N, wherein L is geodetic longitude, B is geodetic latitude, and H is elevation; the load beam center point E is fixed at a position vector of
Figure BDA0001779534350000031
Wherein, the origin of coordinates is earth centroid o, the direction of x is earth centroid pointing to Greenwich mean meridian, z is pointing to the international common origin of north pole, and x, y, z become the right-handed system.
1. Regional target location vector
The position vectors of N vertexes of the N-polygon in the earth's fixation system are
Figure BDA0001779534350000032
Figure BDA0001779534350000033
Figure BDA0001779534350000034
Figure BDA0001779534350000035
Wherein R isEEquator of the earthA major semi-axis of an ellipsoid;
RP-the polar radius of the earth;
2. criterion of Euclidean length method
Position vector of N vertexes of N-polygon
Figure BDA0001779534350000036
Position vector with point E
Figure BDA0001779534350000037
Figure BDA0001779534350000041
When N is 2, the polygonal plane represents a circular area, in which case
Figure BDA0001779534350000042
Radius of
Figure BDA0001779534350000043
Order to
Figure BDA0001779534350000044
Then
Figure BDA0001779534350000045
When N > 2, then
Figure BDA0001779534350000046
The criterion of the Euclidean length method for the satellite coverage area target is as follows:
Figure BDA0001779534350000047
where FE is the coverage capacity of the load in the direction of the ground target.
3. Criterion of vector distance method
When N is more than 2, the position vector of the N sides of the N-edge-shaped area in the earth fixation system is as follows:
Figure BDA0001779534350000048
wherein the value space of i-1 is as follows:
Figure BDA0001779534350000049
then vector
Figure BDA00017795343500000410
In the vector
Figure BDA00017795343500000411
Projection d ofi(i ═ 1,2,3, …, N) is:
Figure BDA00017795343500000412
wherein the value space of i-1 is the same as above.
The criterion of the vector distance method for the satellite coverage area target is as follows:
(∪di)≤FE(i＝1,2,3,…,N)。
4. criterion of projection method
4.1. Initial selection of sets
The known load beam center point E is fixed at a position vector of
Figure BDA00017795343500000413
(wherein, the origin of coordinates is the earth centroid o, the direction of x is the earth centroid pointing to Greenwich mean meridian, z is pointing to the international conventional origin of north pole, x, y, z form a right-hand system), and the geodetic coordinates are (L)e,Be,He)。
Figure BDA0001779534350000051
Figure BDA0001779534350000052
Figure BDA0001779534350000053
Comprises the following steps:
||L||∞＝max(Li)，||B||∞＝max(Bi)；
||L||0＝min(Li)，||B||0＝min(Bi)；
then:
(||L||0≤Le≤||L||∞)∩(||B||0≤Be≤||B||∞)。
and if the load beam center point E does not satisfy the formula, directly discarding the load beam, and if the load beam center point E satisfies the formula, performing the next judgment.
4.2. Criterion of projection method
If the projection area of the N-shaped polygon is S, a triangle formed by two adjacent top points of the N-shaped polygon and the E point meeting the set initial selection condition is Se,i,i-1Then triangle Se,i,i-1Has a projection area of Si：
Figure BDA0001779534350000054
Wherein the value space of i-1 is the same as above.
The criterion of the projection method for the satellite coverage area target is as follows:
Figure BDA0001779534350000055
from the above analysis, it can be seen that for a given target area S, one of the above three criteria conditions is satisfied, i.e., the above-mentioned three criteria are collected to be the coverage of the area target by the observation satellite.
Given that a certain satellite a observes a ground object S, the following takes the observation coverage of a to S on a certain day as an example, and the implementation process is given by taking a specific example.
The coordinate sequence (unit: kilometer) of the center point E of the satellite-borne load beam of a certain satellite A is known:
Figure BDA0001779534350000056
Figure BDA0001779534350000061
the load coverage in the ground target direction is 533.305 Km.
The ground area target S is represented by a pentagonal area, and the geodetic coordinates of its vertices are: j1(55.6,281.5,0), J2(57.5,272.1,0), J3(49.2,276.0,0), J4(53.1,277.9, 0), J5(47.1,276.5,0), these five points constituting the pentagonal region.
1. Vertex vector
Coordinate vectors of a ground fixation system of 5 vertexes of the pentagon are respectively as follows:
J1:(701511.675,-3542892.850,5239409.733)；
J2:(125876.965,-3432851.928,5356022.617)；
J3:(436471.456,-4152748.507,4805121.687)；
J4:(527482.768,-3801363.453,5077233.157)；
J5:(492383.096,-4321590.972,4649339.572)；
2. coverage of regional targets S by satellite A
Start Time(UTC)Stop Time(UTC)
2016-07-09 07:46:41 2016-07-09 07:52:08；
2016-07-09 18:30:47 2016-07-09 18:36:20。

Claims (1)

1. A method for rapidly calculating the observation of an in-orbit satellite on a target in a specific area is characterized by comprising the following steps:
1) for a given regional target, the regional shape is represented by an N-polygon, wherein N is a positive integer greater than 1; the geodetic coordinates of the N vertices of the N-polygon are known as (L)i,Bi,Hi) 1,2,3, 9, N, wherein L is geodetic longitude, B is geodetic latitude, and H is elevation; the geodetic coordinate of the load beam center point E is (L)e,Be,He) The earth fixation position vector is
Figure FDA0001779534340000011
Calculating the position vector of N vertexes of the N-polygon in the earth fixation system as
Figure FDA0001779534340000012
Wherein,
Figure FDA0001779534340000013
REsemi-major axis, R, representing the equatorial ellipsoid of the earthPRepresents the polar radius of the earth;
Figure FDA0001779534340000014
2) the criterion of the target in the satellite coverage area by adopting the Euclidean length method is
Figure FDA0001779534340000015
Wherein FE is the covering capacity of the load in the ground target direction, and the position vectors of N vertexes of the N-polygon
Figure FDA0001779534340000016
Position vector with point E
Figure FDA0001779534340000017
3) The criterion of the satellite coverage area target by adopting the vector distance method is that ([ U ] d)i) FE ≦ wherein the vector
Figure FDA0001779534340000018
In the vector
Figure FDA0001779534340000019
Projection of
Figure FDA00017795343400000110
When N is more than 2, the position vector of N sides of the N-edge-shaped area in the earth fixation system
Figure FDA00017795343400000111
4) When a projection method is adopted to judge a target in a satellite coverage area, firstly, set primary selection is carried out, and if (| | L |) does not satisfy (| L | computationally branch cells0≤Le≤||L||∞)∩(||B||0≤Be≤||B||∞) The central point E of the load beam is directly discarded, | | L | | non-woven phosphor∞＝max(Li)，||B||∞＝max(Bi)，||L||0＝min(Li)，||B||0＝min(Bi) (ii) a If the load beam center point E is met, further judging that the E point and two adjacent vertexes of the N-polygon form a triangle Se,i,i-1Projected area of
Figure FDA00017795343400000112
Whether or not to satisfy
Figure FDA00017795343400000113
Wherein S is the projection area of the N-shaped polygon, and if the projection area meets the requirement, the N-shaped polygon has covering capability;
5) for a given regional target, one of the criteria in steps 2), 3) and 4) is satisfied, namely that the observation satellite has the coverage capability for the regional target.
CN201810985165.7A 2018-08-28 2018-08-28 Rapid calculation method for observing specific region target by on-orbit satellite Expired - Fee Related CN109145835B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810985165.7A CN109145835B (en) 2018-08-28 2018-08-28 Rapid calculation method for observing specific region target by on-orbit satellite

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810985165.7A CN109145835B (en) 2018-08-28 2018-08-28 Rapid calculation method for observing specific region target by on-orbit satellite

Publications (2)

Publication Number Publication Date
CN109145835A CN109145835A (en) 2019-01-04
CN109145835B true CN109145835B (en) 2021-06-08

Family

ID=64828310

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810985165.7A Expired - Fee Related CN109145835B (en) 2018-08-28 2018-08-28 Rapid calculation method for observing specific region target by on-orbit satellite

Country Status (1)

Country Link
CN (1) CN109145835B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114665954B (en) * 2022-04-26 2024-07-02 南京熊猫汉达科技有限公司 Satellite beam selection method and device based on communication area beam coverage rate
CN116975504B (en) * 2023-09-22 2023-12-15 中科星图测控技术股份有限公司 Rapid calculation method for satellite reconnaissance coverage area target

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101430379A (en) * 2007-11-07 2009-05-13 中国科学院电子学研究所 Synthetic aperture radar three-dimensional microwave imaging method for circular track of earth synchronization orbit
CN106156472A (en) * 2015-04-20 2016-11-23 中国科学院遥感与数字地球研究所 A kind of remote sensing satellite covering power over the ground analyzes method
CN108334979A (en) * 2018-01-05 2018-07-27 合肥工业大学 More imaging satellite mission planning methods of facing area covering

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8041118B2 (en) * 2007-02-16 2011-10-18 The Boeing Company Pattern recognition filters for digital images

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101430379A (en) * 2007-11-07 2009-05-13 中国科学院电子学研究所 Synthetic aperture radar three-dimensional microwave imaging method for circular track of earth synchronization orbit
CN106156472A (en) * 2015-04-20 2016-11-23 中国科学院遥感与数字地球研究所 A kind of remote sensing satellite covering power over the ground analyzes method
CN108334979A (en) * 2018-01-05 2018-07-27 合肥工业大学 More imaging satellite mission planning methods of facing area covering

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"基于多边形布尔运算的卫星区域覆盖分析算法";汪荣峰;《装备学院学报》;20160430;第27卷(第2期);第83-87页 *
"面向任意几何区域的遥感卫星对地覆盖法";沈夏炯 等;《计算机应用研究》;20160731;第33卷(第7期);第1999-2002、2013页 *

Also Published As

Publication number Publication date
CN109145835A (en) 2019-01-04

Similar Documents

Publication Publication Date Title
CN104931022B (en) Satellite image stereoblock adjustment method based on spaceborne laser altimeter system data
CN111680354B (en) Method for calculating self-intersection point of orbit of near-earth regression orbit satellite subsatellite point and photographing point
WO2020134856A1 (en) Remote sensing satellite system
JP6290539B2 (en) An adaptive method for estimating the number of electrons in the ionosphere.
CN108253987B (en) Unmanned aerial vehicle trajectory planning method, device and storage device based on A-x algorithm
CN106871932B (en) Satellite-borne laser on-orbit pointing calibration method based on pyramid search terrain matching
CN105698762A (en) Rapid target positioning method based on observation points at different time on single airplane flight path
CN109145835B (en) Rapid calculation method for observing specific region target by on-orbit satellite
CN104913780A (en) GNSS-CCD-integrated zenith telescope high-precision vertical deflection fast measurement method
CN113343777B (en) Satellite rapid planning method and system for natural disaster emergency response
CN111369436A (en) Airborne LiDAR point cloud rarefying method considering multi-terrain features
CN114490815B (en) Unmanned aerial vehicle task load reconnaissance coverage rate calculation method based on topography
CN109934868A (en) One kind is based on three-dimensional point cloud and the matched vehicle positioning method of satellite mapping
CN111121785A (en) Non-road path planning method based on graph search
US11568638B2 (en) Image targeting via targetable 3D data
CN117389305A (en) Unmanned aerial vehicle inspection path planning method, system, equipment and medium
CN110456319A (en) A kind of radar intervisibility calculation method based on SRTM
Tan et al. An adaptive crack inspection method for building surface based on BIM, UAV and edge computing
CN114741907B (en) Earth center angle-based method for rapidly forecasting satellite transit in ground circular area
CN110160503B (en) Unmanned aerial vehicle landscape matching positioning method considering elevation
Hassani et al. A new point-cloud-based lidar/imu localization method with uncertainty evaluation
CN114754741A (en) Slope section gradient measuring method
Zongjian et al. Accuracy analysis of low altitude photogrammetry with wide-angle camera
CN109444915A (en) A kind of danger zone pre-judging method based on laser radar data
CN110310370B (en) Method for point-plane fusion of GPS (Global positioning System) and SRTM (short Range TM)

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210608