CN110850510B - Method for reconstructing vertical rainfall field based on satellite-ground link scanning - Google Patents
Method for reconstructing vertical rainfall field based on satellite-ground link scanning Download PDFInfo
- Publication number
- CN110850510B CN110850510B CN201911073941.7A CN201911073941A CN110850510B CN 110850510 B CN110850510 B CN 110850510B CN 201911073941 A CN201911073941 A CN 201911073941A CN 110850510 B CN110850510 B CN 110850510B
- Authority
- CN
- China
- Prior art keywords
- satellite
- rainfall
- vertical
- ground link
- attenuation coefficient
- 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
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/14—Rainfall or precipitation gauges
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Abstract
The invention provides a method for scanning and reconstructing a vertical rainfall field based on satellite-to-ground links, which comprises the steps of constructing a plurality of satellite signal receiving ends to receive polar orbit satellite signals to form a plurality of satellite-to-ground links; gridding a vertical area to be reconstructed; the spatial distribution of the attenuation coefficient is obtained by analyzing the relation among the satellite-ground link signal intensity information, the rainfall attenuation coefficient and the satellite-ground link length, the spatial distribution of rainfall is inverted according to the relation model of the attenuation coefficient and the rainfall intensity, rainfall information on a vertical space is generated, and finally the rainfall information and a detection area are combined to reconstruct a vertical rainfall field. The method not only makes up the defects of detecting the vertical rainfall field by the traditional means, but also has the advantages of high space-time resolution, low cost and the like, and provides important data support and guarantee for water resource management, agricultural development, natural disaster early warning and the like.
Description
Technical Field
The invention relates to the field of meteorological detection, in particular to a method for reconstructing a vertical rainfall field based on satellite-ground link scanning.
Background
Monitoring of rainfall vertical distribution plays an important role in water resource management, agricultural development, early warning of natural disasters (debris flow, flood and the like) and the like. Traditional rainfall detection means include rain gauges, weather radars, and meteorological satellites. The rain gauge can only realize the observation of rainfall in a horizontal area through networking, and is difficult to arrange and maintain in regions with rare occurrence, such as deserts, rainforests and the like, so that the operation cost is increased; the weather radar can detect rainfall in a three-dimensional space, the space-time resolution is high, but the rainfall near the ground is difficult to measure under the influence of ground clutter when the elevation angle is low; the meteorological satellite can observe three-dimensional rainfall in the global range, and the rainfall can be monitored only twice a day in the same region, so that the requirement of rainfall real-time monitoring cannot be met due to too low time resolution.
At present, the number of medium and low orbit satellites operating in orbit exceeds 1000, and with the development of satellite constellation plans of companies such as 'one network', 'boeing' and the like, tens of thousands of medium and low orbit satellites enter space to work in the future, and the working frequency bands of the satellites are concentrated in microwave bands such as Ku, Ka and the like. Relevant researches show that rainfall can strongly attenuate signals of a microwave band, and a theoretical basis is provided for calculating rainfall intensity by utilizing microwave rainfall attenuation information. Meanwhile, the satellite antenna in the wave band is small in size, simple to operate and convenient to maintain, and the cost of a single antenna is far lower than that of a weather radar. Therefore, the vertical rainfall field is reconstructed through the satellite-ground link formed by the double (multiple) ground stations and the polar orbit satellite, so that the method has lower operation and maintenance cost and is used for supplementing the existing means, and the method is helpful for understanding the occurrence and development process of rainfall through the analysis of the vertical rainfall field and also has important effects on atmospheric science research, weather forecast and the like.
Disclosure of Invention
The purpose of the invention is as follows: the invention aims to provide a method for scanning and reconstructing a vertical rainfall field based on a satellite-ground link formed by a ground station and a polar satellite.
The technical scheme is as follows: in order to achieve the technical effects, the invention provides the following technical scheme:
a method for reconstructing a vertical rainfall field based on satellite-ground link scanning is characterized by comprising the following steps:
(1) constructing a satellite-ground link: n satellite signal receiving ends are arranged on the ground of an area where a vertical rainfall field needs to be reconstructed, and q satellite-ground links are formed between the n satellite signal receiving ends and m satellites distributed in the air;
(2) dividing a vertical rainfall field two-dimensional plane to be reconstructed into N multiplied by M rectangular grids with equal size;
(3) during rainfall, collecting the radio wave frequency f of each satellite-ground link i at each sampling time point T in a sampling period TiSignal strength information pi(t) and information L of the length of each satellite-ground link i within each cube it crossesi(t); wherein the content of the first and second substances, representing the length of the satellite-ground link i in the kth grid at the moment t;wherein eirp (t) represents the equivalent total radiated power,represents attenuation due to non-rain factors, Ai(t) represents rainfall induced decay; will be provided withRegarded as constant C, to obtain Ai(t)-C=-pi(t);
(4) Constructing a rain attenuation model of the satellite-ground link i:wherein, γkRepresenting the attenuation coefficient in the kth grid; combined standAnd Ai(t)-C=-pi(t), establishing a system of linear equations:
Lγ=P
γ=[γ1,γ2,γ3,…,γN×M,C]T
P=[-p1(1),…,-p1(T),-p2(1),…,-p2(T),…,…,-pq(1),…,-pq(T)]T
solving a linear equation set to obtain an attenuation coefficient vector gamma;
(5) obtaining attenuation coefficient and rainfall intensity R in each grid according to SAM modelkThen substituting the attenuation coefficient vector obtained in the step (4) into the SAM model to calculate the rainfall intensity R in each gridkAnd then obtaining the rainfall intensity distribution in each grid of the vertical rainfall field to be reconstructed.
Specifically, in the step (4), a joint algebraic reconstruction technique is adopted to solve the attenuation coefficient vector, and the specific steps include:
(21) initializing gamma0Setting a maximum number of iterations K for a vector consisting of 0 or a small numbermaxAnd stopping the iteration error threshold epsilon;
(22) according to the linear equation system L gamma is P, the iterative formula of the joint algebraic reconstruction technology is obtained as follows:
wherein s represents the number of iterations, s has an initial value of 0, λsRelaxation coefficient for the s-th iteration; drAnd DcRespectively diagonal matrices related to the row and column vectors of matrix L,wherein the content of the first and second substances,andrespectively representing the ith row vector and the jth column vector of the matrix L;
(23) when the solution of the attenuation coefficient vector satisfies the condition s > KmaxOr gamma ray2Stopping iteration when the epsilon is less than or equal to epsilon, and outputting a result gamma.
wherein the content of the first and second substances,andrespectively representing the incident point coordinates and the emergent point coordinates of the link in the kth grid.
Specifically, the SAM model is: gamma ═ aRbWherein a and b are power law coefficients:
wherein f is the satellite-to-ground link radio frequency.
Has the advantages that: compared with the traditional rainfall field observation method, the method has the following advantages:
the invention not only can realize the detection of rainfall intensity and vertical spatial distribution thereof, but also is beneficial to providing high-precision high-space-time-resolution key data for rainfall disaster early warning. Meanwhile, the invention fully utilizes the advantages of huge number of middle and low orbit satellites, low cost of satellite antennas, easy operation and the like, and completes the real-time monitoring of the vertical distribution and change of rainfall through the data visualization technology. In addition, the method can also be used as a new vertical rainfall field detection means to be put into meteorological observation service.
Drawings
FIG. 1 is a schematic diagram of the construction of the reconstruction of a vertical rainfall field based on satellite-ground link scanning according to the present invention;
FIG. 2 is a block diagram of a system for implementing the method of the present invention;
fig. 3 is a flow chart of the specific steps of the present invention implemented by the system of fig. 2.
Detailed Description
The present invention will be further described with reference to the accompanying drawings and specific embodiments.
First, a system capable of implementing the method of the present invention is provided, which is shown in fig. 2 and includes: the system comprises a signal attenuation information sensing end, a data processing terminal, a database and a vertical rainfall product service end. The signal attenuation information sensing terminal specifically refers to a satellite signal receiving terminal, such as a satellite antenna, for receiving satellite signals and providing the satellite signals and the position information (x) of the antennai,yi) And position information of satelliteTransmitting the data into a data processing terminal and a database; the data processing terminal performs satellite signal processing (extracting frequency information f from satellite signal)iAnd intensity information pi) And performing rainfall information inversion and vertical rainfall field reconstruction, and storing the processing result in a database.
The specific flow for implementing the method of the present invention by the system is shown in fig. 3, and includes the steps of:
the method comprises the following steps: as shown in fig. 1, dividing an area where a vertical rainfall field needs to be reconstructed into grids formed by N × M labeled rectangles with equal sizes, and storing grid coordinate information into a database;
step two: erecting n (n is more than or equal to 1) ground antennas, receiving signals of m (m is more than or equal to 1) satellites, and forming q (q is less than or equal to nxm) satellite-ground links; in rainfall, frequency information f of each link receiving signal is collected at each sampling time point T in a sampling period TiAnd intensity information piAnd together with position information (x) of the antennai,yi) And position information of satelliteThe data are transmitted into a data processing terminal and stored in a database;
step (ii) ofThirdly, the method comprises the following steps: the data processing terminal analyzes the signal intensity information p of the ith link at t moment in real timei(t) and length information L of links within each meshi(t) and storing in a database, specifically:
1. the signal strength information of the ith link at the time t is as follows:
in which eirp (t),and AiAnd (t) respectively represents equivalent full radiation power, attenuation caused by non-rain factors and attenuation caused by rainfall, and the equivalent full radiation power, the attenuation caused by non-rain factors and the attenuation caused by rainfall are constant in one detection time and can be approximately processed as a constant.
The above formula is changed into:
Ai(t)-C=-pi(t) (2)
2. the length of the ith link in the kth grid at the time t is as follows:
in the formulaAndrespectively representing the incident and emergent points of the link within the kth grid.
3. Length information L of ith link at time ti(t) can be expressed as:
step four: the intensity information and the length information in the third step are transmitted from the database to the data processing terminalCombining the formed intensity vector P and the length matrix L to form a linear equation system, and solving the attenuation coefficient gamma in each gridkThe method comprises the following specific steps:
1. the rain fade information for the ith link at time t is shown as:
in the formula of gammakAttenuation coefficient in the kth grid.
2. After time T, the attenuation coefficient gamma to be obtainedkIntensity information pi(t) Length information Li(t) forming an attenuation coefficient vector γ, an intensity vector P and a length matrix L, respectively:
γ=[γ1,γ2,γ3,…,γN×M,C]T (6)
P=[-p1(1),…,-p1(T),-p2(1),…,-p2(T),…,…,-pq(1),…,-pq(T)]T (7)
combining equations (2) and (5) establishes a linear system of equations:
Lγ=P (9)
3. solving the attenuation coefficient vector by a joint algebraic reconstruction technology, specifically comprising the following steps:
in the formula, gammakRepresents the attenuation coefficient vector solution after the kth iteration; lambda [ alpha ]kIs the relaxation coefficient; drAnd DcDiagonal matrices are associated with the row and column vectors of matrix L, respectively:
Stopping iteration when the attenuation coefficient vector solution meets the condition, specifically:
in the formula, epsilon and KmaxRespectively representing an error threshold and a maximum number of iterations.
Step five: computing rainfall R in each grid by using SAM model and attenuation coefficient vector in step fourkThe method specifically comprises the following steps:
1. the attenuation coefficient and the rainfall intensity satisfy a power law relation:
γ=aRb (14)
2. the power law coefficients a and b of the SAM model are respectively as follows:
wherein f is the frequency of the satellite-ground link.
3. From the attenuation coefficient vector γ, the rainfall intensity is calculated:
step six: and (4) transmitting the area grid information from the database to a data processing terminal, fusing the area grid information with the rainfall information intensity in the step five, completing reconstruction of the vertical rainfall field, storing the reconstruction information in the database, and then performing visual and dynamic display through a vertical rainfall product service end, wherein the vertical rainfall product service end comprises an entity terminal capable of performing data display, such as a smart phone, a PC terminal and a related service website.
The above description is only of the preferred embodiments of the present invention, and it should be noted that: it will be apparent to those skilled in the art that various modifications and adaptations can be made without departing from the principles of the invention and these are intended to be within the scope of the invention.
Claims (3)
1. A method for reconstructing a vertical rainfall field based on satellite-ground link scanning is characterized by comprising the following steps:
(1) constructing a satellite-ground link: n satellite signal receiving ends are arranged on the ground of an area where a vertical rainfall field needs to be reconstructed, and q satellite-ground links are formed between the n satellite signal receiving ends and m satellites distributed in the air;
(2) dividing a vertical rainfall field two-dimensional plane to be reconstructed into N multiplied by M rectangular grids with equal size;
(3) during rainfall, collecting the radio wave frequency f of each satellite-ground link i at each sampling time point T in a sampling period TiSignal strength information pi(t) and information L of the length of each satellite-ground link i within each cube it crossesi(t); wherein the content of the first and second substances, representing the length of the satellite-ground link i in the kth grid at the moment t;wherein eirp (t) represents the equivalent total radiated power,represents attenuation due to non-rain factors, Ai(t) represents rainfall induced decay; will be provided withRegarded as constant C, to obtain Ai(t)-C=-pi(t);
(4) Constructing a rain attenuation model of the satellite-ground link i:wherein, γkRepresenting the attenuation coefficient in the kth grid; combined standAnd Ai(t)-C=-pi(t), establishing a system of linear equations:
Lγ=P
γ=[γ1,γ2,γ3,...,γN×M,C]T
P=[-p1(1),...,-p1(T),-p2(1),...,-p2(T),...,...,-pq(1),...,-pq(T)]T
solving a linear equation set to obtain an attenuation coefficient vector gamma;
(5) obtaining attenuation coefficient and rainfall intensity R in each grid according to SAM modelkThen substituting the attenuation coefficient vector obtained in the step (4) into the SAM model to calculate the rainfall intensity R in each gridkThen obtaining the vertical rainfall field to be reconstructed in each gridDistribution of rainfall intensity;
the SAM model is as follows: gamma ═ aRbWherein a and b are power law coefficients:
wherein f is the satellite-to-ground link radio frequency.
2. The method for reconstructing the vertical rainfall field based on the satellite-ground link scanning as claimed in claim 1, wherein the attenuation coefficient vector is solved in the step (4) by using a joint algebraic reconstruction technique, and the specific steps include:
(21) initializing gamma0Setting a maximum number of iterations K for a vector consisting of 0 or a small numbermaxAnd stopping the iteration error threshold epsilon;
(22) according to the linear equation system L gamma is P, the iterative formula of the joint algebraic reconstruction technology is obtained as follows:
wherein s represents the number of iterations, s has an initial value of 0, λsRelaxation coefficient for the s-th iteration; drAnd DcRespectively diagonal matrices related to the row and column vectors of matrix L,wherein the content of the first and second substances,andrespectively representing the ith row vector and the jth column vector of the matrix L;
(23) when the solution of the attenuation coefficient vector satisfies the condition s > KmaxOr gamma ray2Stopping iteration when the epsilon is less than or equal to epsilon, and outputting a result gamma.
3. The method for reconstructing a vertical rainfall field based on satellite-to-ground link scanning according to claim 1, wherein the method comprisesThe calculation formula of (2) is as follows:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911073941.7A CN110850510B (en) | 2019-11-05 | 2019-11-05 | Method for reconstructing vertical rainfall field based on satellite-ground link scanning |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911073941.7A CN110850510B (en) | 2019-11-05 | 2019-11-05 | Method for reconstructing vertical rainfall field based on satellite-ground link scanning |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110850510A CN110850510A (en) | 2020-02-28 |
CN110850510B true CN110850510B (en) | 2021-11-12 |
Family
ID=69598720
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911073941.7A Active CN110850510B (en) | 2019-11-05 | 2019-11-05 | Method for reconstructing vertical rainfall field based on satellite-ground link scanning |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110850510B (en) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111505739B (en) * | 2020-03-30 | 2022-07-22 | 中国人民解放军国防科技大学 | Method for judging weather time by using high-frequency satellite-to-ground link multi-feature quantity |
CN112131989B (en) * | 2020-09-15 | 2021-07-23 | 河海大学 | Millimeter wave rain measurement model parameter obtaining method based on space rainfall data |
CN112712195B (en) * | 2020-12-17 | 2024-01-05 | 杭州鲁尔物联科技有限公司 | Rainfall attenuation coefficient prediction method, rainfall attenuation coefficient prediction device, rainfall attenuation coefficient prediction equipment and storage medium |
CN113009595B (en) * | 2021-02-09 | 2023-05-02 | 海南省气象科学研究所 | Rainfall monitoring method and related device based on vehicle-mounted sensor |
CN113095590B (en) * | 2021-04-29 | 2022-04-29 | 中国人民解放军国防科技大学 | High spatial-temporal resolution reconstruction analysis and short-term prediction method for microwave horizontal rainfall field |
CN113900158A (en) * | 2021-09-30 | 2022-01-07 | 中国人民解放军国防科技大学 | System and method for satellite-ground link chromatography vertical cloud field based on ground antenna array |
CN115225136B (en) * | 2022-06-06 | 2023-06-23 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | Reconstruction method of satellite-ground link propagation attenuation region |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1744185A1 (en) * | 2005-07-13 | 2007-01-17 | CIMA-Università degli Studi di Genova | Method for estimating a variable geophysical field particularly for forecasting atmospheric precipitation |
CN106547036B (en) * | 2016-10-21 | 2019-01-01 | 中国人民解放军理工大学 | A kind of united Regional Precipitation measurement method of multiband microwave link |
CN106324580B (en) * | 2016-11-02 | 2018-11-06 | 中国人民解放军理工大学 | A kind of radar return decaying correction method based on microwave links road network |
KR101875906B1 (en) * | 2017-11-09 | 2018-07-06 | (주)에스이랩 | Apparatus of calculating areal precipitation and method thereof |
CN109697323B (en) * | 2018-12-27 | 2022-08-16 | 武汉大学 | Rainfall observation method integrating satellite remote sensing and mobile communication base station signals |
CN110031916B (en) * | 2019-03-07 | 2021-10-01 | 中国人民解放军国防科技大学 | Rainfall intensity measurement method based on satellite-ground link attenuation effect |
-
2019
- 2019-11-05 CN CN201911073941.7A patent/CN110850510B/en active Active
Non-Patent Citations (1)
Title |
---|
Reconstruction of vertical rainfall fields using satellite communication links;S XI;《Asia-pacific Conference on Communications》;20171231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110850510A (en) | 2020-02-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110850510B (en) | Method for reconstructing vertical rainfall field based on satellite-ground link scanning | |
Fang et al. | Evaluation of the TRMM 3B42 and GPM IMERG products for extreme precipitation analysis over China | |
CN111666656A (en) | Rainfall estimation method and rainfall monitoring system based on microwave rainfall attenuation | |
Xing et al. | Estimation of daily mean land surface temperature at global scale using pairs of daytime and nighttime MODIS instantaneous observations | |
CN110766795B (en) | Method and system for detecting three-dimensional rainfall field by using three-dimensional star-earth link dense-woven net | |
Liang et al. | A composite approach of radar echo extrapolation based on TREC vectors in combination with model-predicted winds | |
CN104636608A (en) | Direct assimilation method of MODIS satellite data | |
Sieglaff et al. | A satellite-based convective cloud object tracking and multipurpose data fusion tool with application to developing convection | |
CN105678085A (en) | PM2.5 concentration estimation method and system | |
Liang et al. | Meteorological applications of precipitable water vapor measurements retrieved by the national GNSS network of China | |
Min et al. | Coverage of China new generation weather radar network | |
Di Paola et al. | Combined MW-IR Precipitation Evolving Technique (PET) of convective rain fields | |
CN115267786A (en) | Resunet-GAN global radar inversion method and device fusing satellite observation brightness and elevation | |
CN112526636B (en) | Near-ground two-dimensional rainfall field reconstruction method and system based on single-station multi-satellite and multi-station networking | |
Juglea et al. | Soil moisture modelling of a SMOS pixel: interest of using the PERSIANN database over the Valencia Anchor Station | |
Chan et al. | Ground-space-sky observing system experiment during tropical cyclone Mulan in August 2022 | |
CN112285709B (en) | Atmospheric ozone remote sensing laser radar data fusion method based on deep learning | |
Féral et al. | Large-scale modeling of rain fields from a rain cell deterministic model | |
CN104360040A (en) | Remote sensing soil moisture content monitoring method based on STARFM fusion technology | |
Braun et al. | Development of a water vapor tomography system using low cost L1 GPS receivers | |
CN116609857A (en) | Cloud vertical structure parameter estimation method based on visible light, infrared and microwave images | |
CN114239379A (en) | Transmission line geological disaster analysis method and system based on deformation detection | |
Takeuchi | An introduction of advanced technology for tropical cyclone observation, analysis and forecast in JMA | |
Xiong et al. | QpefBD: a benchmark dataset applied to machine learning for minute-scale quantitative precipitation estimation and forecasting | |
Zhu et al. | FPCI and SPCI are proposed to distinguish the frontal and saturated precipitation systems based on FY-4A |
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 |