CN109118535A - A method of accurately calculating low orbit satellite windward area - Google Patents
A method of accurately calculating low orbit satellite windward area Download PDFInfo
- Publication number
- CN109118535A CN109118535A CN201810622993.4A CN201810622993A CN109118535A CN 109118535 A CN109118535 A CN 109118535A CN 201810622993 A CN201810622993 A CN 201810622993A CN 109118535 A CN109118535 A CN 109118535A
- Authority
- CN
- China
- Prior art keywords
- satellite
- area
- windward
- referring
- windward side
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/90—Determination of colour characteristics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10024—Color image
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Image Processing (AREA)
Abstract
The present invention provides a kind of methods for accurately calculating low orbit satellite windward area, include the following steps: step 1, establish satellite threedimensional model and carry out color mark to satellite threedimensional model and background；Step 2, obtain satellite windward side various in-orbit state diagrams and referring to figure；Step 3 using colouring information identification satellite region and generates binary segmentation figure；Step 4 calculates the elemental area of satellite windward side according to binary segmentation figure and referring to the elemental area in figure Satellite region；Step 5 utilizes the real area that satellite windward side is calculated referring to figure information.The present invention can accurately calculate the satellite windward area of complex contour, and can accurately analyze influence of the celestial body surface pivots component to satellite windward area, can provide important evidence for the design of the in-orbit management system of satellite orbital altitude.
Description
Technical field
The present invention relates to technical field of satellite overall design, and in particular to one kind accurately calculates low orbit satellite windward area
Method.
Background technique
Low orbit satellite in orbit when, unavoidably influenced by atmospheric drag, satellite orbital altitude can constantly decline,
The satellite transit orbital period is caused gradually to shorten, Track of Sub-Satellite Point is gradually deviated from nominal position, pair of Satellite Payloads
Ground is covering the whole world or coverage hole occurs in region picture mosaic, is unable to satisfy the regular covering performance requirement of satellite on a surface target.
The atmospheric drag that satellite is subject to mainly with the track speed of service, windward area and the Atmospheric Density Models of satellite
It is related.Wherein, the track speed of service of satellite is stablized near the same orbit altitude, can accurately calculate；Atmospheric density mould
Type also comparative maturity, if it is possible to accurately calculate the windward area of satellite, so that it may accurately calculate the in-orbit atmospheric drag of satellite
Influence situation, so as to realize that the control of satellite Autonomous orbit altitude or autonomous Orbit height such as adjust at the autonomous in-orbit pipe in a small amount
Reason measure changes and infuses the mode that number carries out orbit altitude control operation by ground at present, and ground is greatly reduced and carries out satellite
The workload of maintenance.Meanwhile the maintenance of satellite Autonomous orbit altitude can not also interrupt the imaging job process of load, improve satellite
Effective utilization.
Since there are the solar battery array of one or two large area, traditional windward area meters for most of low orbit satellites
Satellite is mostly reduced to the combination of several regular voxels (such as cube, cylindrical body, cone) by calculation method using method of geometry,
The satellite capsule plate area in satellite flight direction is simply added the area of solar battery array as the windward area of satellite.But
The conventional method can not accurately calculate the satellite windward area of complex contour, while can not calculate celestial body surface and load rotation section
The in-orbit situation of change of satellite windward area of part (such as solar battery array, antenna, payload hood) so that satellite without
Method realizes the measures such as autonomous Orbit height controls or autonomous Orbit height adjusts in a small amount.It is therefore proposed that one kind accurately calculates low rail
The method of satellite windward area is one of the main path solved the problems, such as.
Currently without the explanation or report for finding technology similar to the present invention, it is also not yet collected into money similar both at home and abroad
Material.
Summary of the invention
The present invention is to overcome the above deficiencies in the existing technologies, provides one kind and accurately calculates low orbit satellite windward side
The method of area.
To achieve the above object, the present invention is realized especially by following technical scheme:
A kind of method accurately calculating low orbit satellite windward area includes the following steps:
Step 1 establishes satellite threedimensional model and carries out color mark to satellite threedimensional model and background；
Step 2, obtain satellite windward side various in-orbit state diagrams and referring to figure；
Step 3 using colouring information identification satellite region and generates binary segmentation figure；
Step 4 calculates the elemental area of satellite windward side according to binary segmentation figure and referring to the pixel in figure Satellite region
Area；
Step 5 utilizes the real area that satellite windward side is calculated referring to figure information.
Preferably, the step 1 Satellite threedimensional model and background need to be labeled as different colors.
Preferably, the step 2 is specifically includes the following steps: firstly, by the fixed step size traversal satellite orbit complete period
The various working conditions for influencing windward area observe satellite against satellite flight direction, obtain the throwing of satellite using sciagraphy
Then satellite is reduced to known to outer dimension parameter cube by shadow figure, the as various in-orbit state diagrams of satellite windward side
Body still observes satellite against satellite flight direction, obtains the perspective view of satellite using sciagraphy, as referring to figure.
Preferably, the step 3 includes the following steps: the difference that using low orbit satellite color be black, background is white,
It identifies all windward side perspective views and referring to the satellite region in figure and generates binary segmentation figure, wherein satellite region
Pixel value is 1, and background area pixels value is 0.
Preferably, the step 3 specifically comprises the following steps: that reading each width windward side using MATLAB software programming throws
Shadow figure and referring to figure, traversal calculates the color of each pixel and the difference of black in every width picture, if being less than threshold value, in phase
Label is in the binary segmentation figure answered, and is otherwise labeled as 0.
Preferably, the elemental area in the step 4 is by counting the pixel number that pixel value is 1 in binary segmentation figure
It measures to obtain.
Preferably, the real area of the satellite windward side in the step 5 is the pixel faces according to itself and satellite windward side
Product, the real area referring to figure Satellite region have the feature of equal proportion relationship with the elemental area referring to figure Satellite region
It is calculated.
Compared with the prior art, the present invention have it is following the utility model has the advantages that
1, the satellite windward area of complex contour can be accurately calculated；
2, influence of the celestial body surface pivots component to satellite windward area can accurately be analyzed；
3, important evidence can be provided for the design of the in-orbit management system of satellite orbital altitude.
Detailed description of the invention
Upon reading the detailed description of non-limiting embodiments with reference to the following drawings, other feature of the invention and
Purpose will become more apparent upon:
Fig. 1 is a kind of method flow diagram for accurately calculating low orbit satellite windward area of the present invention；
Fig. 2 is certain low orbit satellite appearance schematic diagram；
In figure: 1 is satellite body, and 2 be solar battery array.
Fig. 3 is for certain low orbit satellite windward side perspective view and referring to figure；
In figure: satellite windward side perspective view when a) for solar battery array corner equal to 0 degree, b) it is solar battery array corner
Satellite windward side perspective view when equal to 72 degree, c) be solar battery array corner be equal to 144 degree when satellite windward side perspective view,
D) for solar battery array corner be equal to 216 degree when satellite windward side perspective view, e) be solar battery array corner be equal to 288 degree when
Satellite windward side perspective view, f) be reference figure when satellite body is simplified to the rib long regular hexahedron for being equal to 2 meters.
Fig. 4 is change curve of certain low orbit satellite windward area with solar battery array corner.
Specific embodiment
Elaborate below to the embodiment of the present invention: the present embodiment carries out under the premise of the technical scheme of the present invention
Implement, the detailed implementation method and specific operation process are given.It should be pointed out that those skilled in the art
For, without departing from the inventive concept of the premise, various modifications and improvements can be made, these belong to guarantor of the invention
Protect range.
As shown in Figure 1, a kind of method for accurately calculating low orbit satellite windward area provided in this embodiment, including it is as follows
Step:
Step 1 establishes certain low orbit satellite threedimensional model using Pro/E software, and sets black, background for satellite color
Color is set as white, as shown in Fig. 2, the low orbit satellite is made of ontology and solar battery array two parts, in orbit when, too
Positive cell array will carry out periodic rotary according to the variation at track illumination angle, therefore need to consider too when carrying out windward area calculating
Positive cell array rotation bring influences；
Step 2, firstly, solar cell when traversing the low orbit satellite in-orbit flight with 2 ° by solar battery array corner for step-length
The various working conditions of battle array rotation observe satellite against satellite flight direction, obtain the projection of the low orbit satellite using sciagraphy
Figure, the as various in-orbit state diagrams of satellite windward side, wherein windward side perspective view when 5 specific solar battery array corners is such as
Shown in Fig. 3, then, satellite body is simplified to the long regular hexahedron for being equal to 2 meters of rib, is still defended against the observation of satellite flight direction
Star obtains the perspective view of satellite using sciagraphy, as referring to figure, as shown in Figure 3, it is notable that described all to meet
It wind face perspective view and is finally saved referring to figure by Pro/E software with the picture of BMP format.
Step 3 using low orbit satellite color is black, the difference that background is white, by all windward side perspective views and ginseng
It is identified according to the satellite region in figure and generates binary segmentation figure, wherein satellite region pixel value is 1, background area pixels
Value is 0.Method particularly includes: each width windward side perspective view is read using MATLAB software programming and referring to figure, traversal calculates every
The difference of the color of each pixel and black in width picture, if be less than threshold value, in corresponding binary segmentation figure label for
1, otherwise it is labeled as 0；
Step 4, statistics calculate the pixel quantity that pixel value is 1 in a width binary segmentation figure, as figure centre halfback's constellation
The elemental area S of the low orbit satellite windward side can be calculated in the elemental area in domain in the methodiCentre halfback's constellation is schemed with reference
The elemental area S in domainref；
Step 5, due to when being generated referring to figure satellite be reduced to the long regular hexahedron for being equal to 2 meters of rib, therefore referring to figure Satellite
Real area be Aref=2m × 2m=4m2, the real area A of satellite windward sideiIt can be calculated according to the following formula:
Finally, the low orbit satellite windward area is as shown in Figure 4 with the change curve of solar battery array corner, wherein most
Big value is 25.3719m2, minimum value 4.9056m2。
Specific embodiments of the present invention are described above.It is to be appreciated that the invention is not limited to above-mentioned
Particular implementation, those skilled in the art can make various deformations or amendments within the scope of the claims, this not shadow
Ring substantive content of the invention.
Claims (7)
1. a kind of method for accurately calculating low orbit satellite windward area, which comprises the steps of:
Step 1 establishes satellite threedimensional model and carries out color mark to satellite threedimensional model and background；
Step 2, obtain satellite windward side various in-orbit state diagrams and referring to figure；
Step 3 using colouring information identification satellite region and generates binary segmentation figure；
Step 4 calculates the elemental area of satellite windward side according to binary segmentation figure and referring to the elemental area in figure Satellite region；
Step 5 utilizes the real area that satellite windward side is calculated referring to figure information.
2. a kind of method for accurately calculating satellite windward area according to claim 1, which is characterized in that the step
1 Satellite threedimensional model and background need to be labeled as different colors.
3. a kind of method for accurately calculating satellite windward area according to claim 1, which is characterized in that the step
2 specifically includes the following steps: firstly, by the various work for influencing windward area in the fixed step size traversal satellite orbit complete period
State observes satellite against satellite flight direction, obtains the perspective view of satellite using sciagraphy, as satellite windward side is various
Then satellite is reduced to cube known to outer dimension parameter by in-orbit state diagram, still observe against satellite flight direction
Satellite obtains the perspective view of satellite using sciagraphy, as referring to figure.
4. a kind of method for accurately calculating satellite windward area according to claim 1, which is characterized in that the step
3 include the following steps: the difference that using low orbit satellite color be black, background is white, by all windward side perspective views and ginseng
It is identified according to the satellite region in figure and generates binary segmentation figure, wherein satellite region pixel value is 1, background area pixels
Value is 0.
5. a kind of method for accurately calculating satellite windward area according to claim 1, which is characterized in that the step
3 are specifically comprised the following steps: to be read each width windward side perspective view using MATLAB software programming and referring to figure, traverse and calculate often
The difference of the color of each pixel and black in width picture, if be less than threshold value, in corresponding binary segmentation figure label for
1, otherwise it is labeled as 0.
6. a kind of method for accurately calculating satellite windward area according to claim 1, which is characterized in that the step
Elemental area in 4 is obtained by counting the pixel quantity that pixel value is 1 in binary segmentation figure.
7. a kind of method for accurately calculating satellite windward area according to claim 1, which is characterized in that the step
The real area of satellite windward side in 5 is according to it with the elemental area of satellite windward side, referring to the reality in figure Satellite region
Interphase product with the elemental area referring to figure Satellite region there is the feature calculation of equal proportion relationship to obtain.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810622993.4A CN109118535A (en) | 2018-06-15 | 2018-06-15 | A method of accurately calculating low orbit satellite windward area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810622993.4A CN109118535A (en) | 2018-06-15 | 2018-06-15 | A method of accurately calculating low orbit satellite windward area |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109118535A true CN109118535A (en) | 2019-01-01 |
Family
ID=64822416
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810622993.4A Pending CN109118535A (en) | 2018-06-15 | 2018-06-15 | A method of accurately calculating low orbit satellite windward area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109118535A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111125832A (en) * | 2019-12-11 | 2020-05-08 | 上海卫星工程研究所 | Method for acquiring position and area of pressure center of windward side of spacecraft |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150178926A1 (en) * | 2013-12-19 | 2015-06-25 | Google Inc. | Obtaining geographic-location related information based on shadow characteristics |
| CN106600434A (en) * | 2016-10-18 | 2017-04-26 | 河南省农业科学院农业经济与信息研究所 | Remote crop growth status monitoring method based on crop model and assimilation technology |
| CN106780646A (en) * | 2016-12-01 | 2017-05-31 | 中科唯实科技（北京）有限公司 | A kind of printenv background modeling method for being applicable many scenes |
| CN107590104A (en) * | 2017-09-07 | 2018-01-16 | 四川九洲北斗导航与位置服务有限公司 | Vehicular load computational methods and device |
-
2018
- 2018-06-15 CN CN201810622993.4A patent/CN109118535A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150178926A1 (en) * | 2013-12-19 | 2015-06-25 | Google Inc. | Obtaining geographic-location related information based on shadow characteristics |
| CN106600434A (en) * | 2016-10-18 | 2017-04-26 | 河南省农业科学院农业经济与信息研究所 | Remote crop growth status monitoring method based on crop model and assimilation technology |
| CN106780646A (en) * | 2016-12-01 | 2017-05-31 | 中科唯实科技（北京）有限公司 | A kind of printenv background modeling method for being applicable many scenes |
| CN107590104A (en) * | 2017-09-07 | 2018-01-16 | 四川九洲北斗导航与位置服务有限公司 | Vehicular load computational methods and device |
Non-Patent Citations (2)
| Title |
|---|
| 吴娱 等: "《数字图像处理》", 31 October 2017 * |
| 杨成 等: "应用阴影图的航天器迎风面积计算方法", 《计算机辅助设计与图形学学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111125832A (en) * | 2019-12-11 | 2020-05-08 | 上海卫星工程研究所 | Method for acquiring position and area of pressure center of windward side of spacecraft |
| CN111125832B (en) * | 2019-12-11 | 2023-11-10 | 上海卫星工程研究所 | Method for acquiring windward side core pressing position and area of spacecraft |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110554704B (en) | Unmanned aerial vehicle-based fan blade autonomous inspection method | |
| CN110244766B (en) | Planning method and system for unmanned aerial vehicle routing inspection route of photovoltaic power station | |
| CN110879601B (en) | Unmanned aerial vehicle inspection method for unknown fan structure | |
| CN103413272B (en) | Low spatial resolution multi-source Remote Sensing Images Space Consistency bearing calibration | |
| Li et al. | Current issues in high-resolution earth observation technology | |
| CN107966112A (en) | A kind of large scale rotor movement parameter measurement method | |
| CN104463166B (en) | The extracting method and system of target water body | |
| CN109002048B (en) | Multi-rotor unmanned aerial vehicle large-scale centralized photovoltaic power station image data acquisition method | |
| CN110006407A (en) | Based on rotor wing unmanned aerial vehicle close to photogrammetric survey method | |
| CN111311650A (en) | Point cloud data registration method and device and storage medium | |
| CN110727903B (en) | Satellite task planning method for realizing maximum observation area by limited coverage resources | |
| CN109118535A (en) | A method of accurately calculating low orbit satellite windward area | |
| CN108716919A (en) | Plant protection drone path planning method based on arbitrary polygon clear area | |
| CN105551032A (en) | Pole image collection system and method based on visual servo | |
| US20220383755A1 (en) | Unmanned aerial vehicle positioning method based on millimeter-wave radar | |
| CN106331684A (en) | Three-dimensional image obtaining method based on small unmanned aerial vehicle video recording in engineering geological survey | |
| CN110389114A (en) | A kind of medium and small water quality kind identification method based on unmanned plane imaging spectral | |
| CN109840360A (en) | A kind of satellite faces constellation smallest size design method under the detection mode of side | |
| CN113776504B (en) | High-precision photographing and quality control method for water engineering unmanned aerial vehicle with complex structure | |
| CN114926739A (en) | Unmanned collaborative acquisition and processing method for underwater and overwater geographic spatial information of inland waterway | |
| CN109178345A (en) | A kind of holder direction and celestial body posture cooperative control method for aerial tracking of maneuvering target | |
| CN114463521B (en) | Building target point cloud rapid generation method for air-ground image data fusion | |
| Zhang et al. | Task-driven on-Board Real-Time Panchromatic Multispectral Fusion Processing Approach for High-Resolution Optical Remote Sensing Satellite | |
| CN206709846U (en) | A kind of aerial survey of unmanned aerial vehicle photo control point surface mark | |
| CN109684593A (en) | Mountainous region photovoltaic project array pitch calculation method and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190101 |