CN102540148A - Inversion algorithm of rainfall three-dimensional distribution parameters - Google Patents

Inversion algorithm of rainfall three-dimensional distribution parameters Download PDF

Info

Publication number
CN102540148A
CN102540148A CN2012100326079A CN201210032607A CN102540148A CN 102540148 A CN102540148 A CN 102540148A CN 2012100326079 A CN2012100326079 A CN 2012100326079A CN 201210032607 A CN201210032607 A CN 201210032607A CN 102540148 A CN102540148 A CN 102540148A
Authority
CN
China
Prior art keywords
rainfall
distribution parameter
horizontal
distribution
radar
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
Application number
CN2012100326079A
Other languages
Chinese (zh)
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.)
University of Shanghai for Science and Technology
Original Assignee
University of Shanghai for Science and Technology
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 University of Shanghai for Science and Technology filed Critical University of Shanghai for Science and Technology
Priority to CN2012100326079A priority Critical patent/CN102540148A/en
Publication of CN102540148A publication Critical patent/CN102540148A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses an inversion algorithm of rainfall three-dimensional distribution parameters; the inversion algorithm comprises the following steps of: one. obtaining total echo data by using space-borne wave band synthesized aperture radar, inverting a starting point of rainfall and the minimal value of total radar echoes, and obtaining an interzone rainfall distribution expression by adopting an VIE inversion algorithm; three. inverting rainfall horizontal distribution parameters according to the interzone rainfall distribution, and obtaining vertical distribution parameters by using a vertical distribution inversion algorithm; three. multiplying the horizontal distribution parameters and the vertical distribution parameters according to the inverted rainfall horizontal distribution parameters and vertical distribution parameter, and obtaining the rainfall three-dimensional distribution parameters. According to the inversion algorithm, when the horizontal distribution parameters are processed, only one length of region needs to be inverted, thereby the computation burden can be substantially reduced; only one line of echo data needs to be processed by adopting a vertical distribution counting manner, thereby the computation burden is substantially reduced, and the accuracy of inverting the rainfall three-dimensional distribution parameters can be kept.

Description

The inversion algorithm of the three-dimensional distribution parameter of a kind of rainfall
Technical field
The present invention relates to the inversion algorithm of the three-dimensional distribution parameter of a kind of rainfall, belong to the meteorological remote sensing technical field.
Background technology
The comparison of the inversion algorithm of rainfall distribution parameter typical case inversion algorithm is the VIE inversion algorithm at present; For example; Frank Silvio Marzano and James A. Weinman. Inversion of Spaceborne X-Band Synthetic Aperture Radar Measurements for Precipitation Remote Sensing Over Land. IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING; 2008; 46 (11): 3472-3487. this article has proposed a kind of spaceborne
Figure 546844DEST_PATH_IMAGE001
wave band synthetic-aperture radar land rainfall remote sensing survey inversion algorithm; Set forth the VIE inversion algorithm in this article; If
Figure 118771DEST_PATH_IMAGE002
is the piecewise function of rainfall zone horizontal ordinate
Figure 684882DEST_PATH_IMAGE003
, then
Figure 658654DEST_PATH_IMAGE004
(1)
Wherein
Figure 284807DEST_PATH_IMAGE005
is constant; is the volume scattering particle transverse axis length at place, wave beam plane;
Figure 982341DEST_PATH_IMAGE007
is the radar incident angle;
Figure 138515DEST_PATH_IMAGE008
is the rainfall area height;
Figure 607674DEST_PATH_IMAGE009
is , and
Figure 124423DEST_PATH_IMAGE011
is the horizontal ordinate of rainfall area position.
Also set forth the inversion algorithm of a kind of vertical distribution parameter
Figure 135104DEST_PATH_IMAGE012
in this article; That is: in interval
Figure 40743DEST_PATH_IMAGE013
(2)
Wherein,
Figure 848480DEST_PATH_IMAGE015
is the height of rainfall area;
Figure 182509DEST_PATH_IMAGE016
is ground rainfall speed, be the variable of interval
Figure 639215DEST_PATH_IMAGE018
interior rainfall area height.
Set forth the method for distilling of rainfall starting point
Figure 826614DEST_PATH_IMAGE019
in the literary composition, its calculation expression is:
Figure 811887DEST_PATH_IMAGE020
(3)
Wherein and
Figure 926791DEST_PATH_IMAGE022
representes respectively: the mean value and the root mean square of 5 total echo points of radar at before the point; The numerical value that satisfies formula (3) is the starting point that is finally inversed by rainfall,
Figure 503900DEST_PATH_IMAGE024
be the total echo of radar that receives.
Set forth the method for distilling of rectangle rainfall distribution rainfall width
Figure 734024DEST_PATH_IMAGE025
in the literary composition, its calculation expression is:
Figure 325542DEST_PATH_IMAGE026
(4)
Wherein
Figure 285408DEST_PATH_IMAGE027
is the total echo minimum value of radar that is finally inversed by, and
Figure 874653DEST_PATH_IMAGE028
is the starting point of rainfall.
Set forth the inversion method of the three-dimensional distribution parameter of rainfall
Figure 654390DEST_PATH_IMAGE029
in the literary composition: promptly horizontal distribution parameter
Figure 49599DEST_PATH_IMAGE030
and vertical distribution parameter multiply each other, and its expression formula is:
Figure 886766DEST_PATH_IMAGE032
(5)
But; In the above-mentioned VIE inversion algorithm; Because in whole rainfall zone;
Figure 153800DEST_PATH_IMAGE002
is the piecewise function of observed range
Figure 352700DEST_PATH_IMAGE003
; That is: not only every distance at a distance from
Figure 693682DEST_PATH_IMAGE033
need recomputate calculating formula (1); Cause the calculated amount accumulation to increase; And; This algorithm at first need calculate by calculating formula (1) respectively the total echo data of each row, just can obtain the three-dimensional distribution parameter of whole rainfall then, has caused calculated amount significantly to increase like this.
Summary of the invention
To the deficiency that above-mentioned prior art exists, the object of the present invention is to provide the inversion algorithm of the three-dimensional distribution parameter of a kind of rainfall, this method can reduce operand significantly under the prerequisite that guarantees precision.
In order to achieve the above object, the technical scheme that the present invention adopts is: the inversion algorithm of the three-dimensional distribution parameter of a kind of rainfall comprises the following steps:
Step 1: adopt spaceborne wave band synthetic-aperture radar to obtain total echo data
Figure 441376DEST_PATH_IMAGE034
, be finally inversed by starting point
Figure 381650DEST_PATH_IMAGE019
, the total echo minimum value of radar
Figure 905035DEST_PATH_IMAGE035
of rainfall; Adopt the VIE inversion algorithm; Obtain the rainfall distribution and expression formula
Figure 310926DEST_PATH_IMAGE037
of interval
Figure 69300DEST_PATH_IMAGE036
; Wherein,
Figure 54891DEST_PATH_IMAGE038
leaves the horizontal ordinate in rainfall zone for radar beam;
Figure 432782DEST_PATH_IMAGE039
is the radar beam incident angle, and
Figure 33528DEST_PATH_IMAGE040
is the rain belt height;
Step 2: according to the rainfall distribution
Figure 513368DEST_PATH_IMAGE037
of above-mentioned interval
Figure 168974DEST_PATH_IMAGE041
; Be finally inversed by rainfall horizontal distribution parameter
Figure 745766DEST_PATH_IMAGE030
, this rainfall horizontal distribution parameter
Figure 720675DEST_PATH_IMAGE030
confirms that the rainfall horizontal distribution is: the shape of the rainfall horizontal distribution of the whole rain belt of inverting, rainfall area width
Figure 671314DEST_PATH_IMAGE025
, ground rainfall speed , use vertical distribution inversion algorithm obtain vertical distribution parameter
Figure 906303DEST_PATH_IMAGE031
;
Step 3: according to above-mentioned rainfall horizontal distribution parameter
Figure 520955DEST_PATH_IMAGE030
and the vertical distribution parameter
Figure 958890DEST_PATH_IMAGE031
of being finally inversed by; Horizontal distribution parameter
Figure 645086DEST_PATH_IMAGE030
and vertical distribution parameter
Figure 789760DEST_PATH_IMAGE031
are multiplied each other, obtain the three-dimensional distribution parameter
Figure 903209DEST_PATH_IMAGE042
of rainfall.
The inversion algorithm of the three-dimensional distribution parameter of a kind of rainfall of the present invention; Its advantage and effect are: this method only needs the zone of length of inverting for
Figure 828440DEST_PATH_IMAGE033
when processing horizontal distributes; Can significantly reduce operand; Adopt vertical distribution parametric inversion algorithm; Only need to handle delegation's echo data, thereby further significantly reduce operand, and can also keep the precision of inverting rainfall distribution parameter.
Description of drawings
Fig. 1 is the inversion algorithm process flow diagram of the three-dimensional distribution parameter of a kind of rainfall of the present invention
Fig. 2 is total reflectogram that the radar of emulation of the present invention receives
The rainfall horizontal distribution figure that Fig. 3 is finally inversed by for the present invention.
Embodiment
Below in conjunction with accompanying drawing and embodiment to further explain of the present invention.
The inversion algorithm of the three-dimensional distribution parameter of a kind of rainfall of the present invention; The total echo that receives when at first adopting spaceborne
Figure 270658DEST_PATH_IMAGE001
wave band synthetic-aperture radar (TerraSAR-X) through the rainfall zone carries out the inverting that rainfall distributes to total echo then.The incident angle of setting synthetic-aperture radar is 30 o, ground scatter is-5 dB, the three-dimensional distribution parameter of rainfall is: the rainfall starting point is 10 Km, and the rainfall width is 10 Km, and horizontal distribution is shaped as rectangle, and Rainfall height is 7 Km, and rainfall speed in ground is 20 mm/h, wherein, Zone position (Km) falls in the axle expression,
Figure 616506DEST_PATH_IMAGE043
Axle expression rain belt height (Km),
Figure 29033DEST_PATH_IMAGE017
Axle expression rainfall speed (mm/h).
The inversion algorithm of the three-dimensional distribution parameter of a kind of rainfall of the present invention, as shown in Figure 1, may further comprise the steps:
Step 1: adopt spaceborne
Figure 729136DEST_PATH_IMAGE001
wave band synthetic-aperture radar to obtain total echo data
Figure 910718DEST_PATH_IMAGE034
, be finally inversed by starting point
Figure 365970DEST_PATH_IMAGE023
, the total echo minimum value of radar
Figure 937897DEST_PATH_IMAGE035
of rainfall; Adopt the VIE inversion algorithm; Obtain the rainfall distribution and expression formula
Figure 805676DEST_PATH_IMAGE037
of interval
Figure 504008DEST_PATH_IMAGE044
; Wherein, leaves the horizontal ordinate in rainfall zone for radar beam;
Figure 756632DEST_PATH_IMAGE039
is the radar beam incident angle, and
Figure 126433DEST_PATH_IMAGE040
is the rain belt height;
A. inverting rainfall starting point
Figure 282608DEST_PATH_IMAGE045
, calculate following formula:
Figure 751767DEST_PATH_IMAGE046
(3)
Wherein
Figure 360602DEST_PATH_IMAGE047
and
Figure 268516DEST_PATH_IMAGE022
representes respectively: the mean value and the root mean square of 5 total echo points of radar at
Figure 482459DEST_PATH_IMAGE003
before the point, and the numerical value that satisfies formula (3) is the starting point
Figure 450415DEST_PATH_IMAGE048
that is finally inversed by rainfall;
B. the size that compares all numerical value in , wherein minimum value is the total echo minimum value of radar
Figure 930255DEST_PATH_IMAGE050
;
C. the rainfall distribution and expression formula
Figure 934300DEST_PATH_IMAGE037
of computation interval
Figure 795443DEST_PATH_IMAGE051
If
Figure 517728DEST_PATH_IMAGE002
is the piecewise function of rainfall zone horizontal ordinate , then the rainfall distribution and expression formula in whole rainfall zone is:
?(1)
When rainfall zone horizontal ordinate
Figure 159DEST_PATH_IMAGE003
when interval
Figure 271216DEST_PATH_IMAGE054
, substitution (1) formula is got
Figure 836507DEST_PATH_IMAGE037
:
Figure 317167DEST_PATH_IMAGE055
(6)
The total echo data of radar that the radar of emulation receives; As shown in Figure 2; Among the figure;
Figure 875187DEST_PATH_IMAGE003
axle is transverse axis position (Km), rainfall zone; axle is the numerical value of the total echo of radar, begins to the radar return terminating point from observation station, gets 100 points and calculates the total echo of radar respectively; Result of calculation is shown in curve among the figure.Curve can be found out from figure, and the echo that radar receives just reduces rapidly after slowly increasing, and after being reduced to minimum value, fast rise is finally got back to-5 dB places again.The above results occurring is: at first radar beam runs into scatterers such as cloud layer, and radar return is slowly increased; Then, radar beam arrives rainfall area, and the decay of rain belt causes the total echo of radar to reduce rapidly; Then, radar beam leaves the rainfall zone gradually, and the decay of rain belt reduces, and finally is zero, remaining area scattering-5 dB, that is: and after the total echo of radar was reduced to minimum value, fast rise was got back to-5 dB places again;
Step 2: according to the rainfall distribution
Figure 891182DEST_PATH_IMAGE037
of above-mentioned interval ; Be finally inversed by rainfall horizontal distribution parameter
Figure 670919DEST_PATH_IMAGE030
, this rainfall horizontal distribution parameter
Figure 66128DEST_PATH_IMAGE030
confirms that the rainfall horizontal distribution is: the shape of the rainfall horizontal distribution of the whole rain belt of inverting, rainfall area width
Figure 880500DEST_PATH_IMAGE025
, ground rainfall speed
Figure 906225DEST_PATH_IMAGE016
, use vertical distribution inversion algorithm obtain vertical distribution parameter
Figure 173258DEST_PATH_IMAGE031
;
A. calculate
Figure 372159DEST_PATH_IMAGE037
derivative more arbitrarily; Be designated as
Figure 713141DEST_PATH_IMAGE057
; if , then rainfall horizontal distribution is shaped as rectangle;
B. rainfall area width
Figure 460834DEST_PATH_IMAGE025
, its calculating formula is:
Figure 197846DEST_PATH_IMAGE059
(4)
Wherein
Figure 658915DEST_PATH_IMAGE027
is the total echo minimum value of radar that is finally inversed by, and
Figure 88759DEST_PATH_IMAGE060
is the starting point of rainfall;
C. the rainfall level is shaped as rectangle, and the maximal value of function
Figure 330384DEST_PATH_IMAGE061
is ground rainfall speed
Figure 811700DEST_PATH_IMAGE016
;
D. be finally inversed by rainfall horizontal distribution parameter
Figure 455171DEST_PATH_IMAGE030
The rainfall horizontal distribution be shaped as rectangle:
When belonging to interval as
Figure 790337DEST_PATH_IMAGE003
, then rainfall horizontal distribution parameter expression is:
(7)
When not belonging to interval as
Figure 768155DEST_PATH_IMAGE064
, then rainfall horizontal distribution parameter expression is:
Figure 428123DEST_PATH_IMAGE066
(8)
Wherein,
Figure 310628DEST_PATH_IMAGE067
is the rainfall starting point, and
Figure 663112DEST_PATH_IMAGE068
is the rainfall area width;
E. use the vertical distribution inversion algorithm to obtain vertical distribution parameter
Figure 605661DEST_PATH_IMAGE031
In interval
Figure 246857DEST_PATH_IMAGE069
,
Figure 667475DEST_PATH_IMAGE070
(2)
Wherein,
Figure 140044DEST_PATH_IMAGE071
is the height of rainfall area;
Figure 925598DEST_PATH_IMAGE072
is ground rainfall speed, and
Figure 850828DEST_PATH_IMAGE073
is the variable of rainfall area height.
Be illustrated in figure 3 as the horizontal distribution parameter of inverting; Among the figure; Transverse axis is transverse axis position, rain belt (Km), and the longitudinal axis is the rainfall speed (mm/h) of rainfall horizontal distribution.In the radar return data that the radar of this emulation receives; In interval
Figure 136633DEST_PATH_IMAGE074
Km scope; Evenly choose 100 points; Calculate the rainfall speed of these 100 points respectively, connect into figure then, shown in rectangle among the figure.From figure, can draw; The rainfall starting point that is finally inversed by is 9.62 km; The rainfall width that is finally inversed by
Figure 36773DEST_PATH_IMAGE025
is 10.262 km, and the ground rainfall speed of inverting
Figure 64772DEST_PATH_IMAGE075
is 19.25 mm/h; Actual rainfall starting point is 10 Km; The rainfall width is 10 Km; Ground rainfall speed
Figure 918458DEST_PATH_IMAGE076
is 20 mm/h, and then rainfall starting point inversion error is 3.8%; Rainfall width inversion error is 2.6%; Ground rainfall speed inversion error is 3.8%, more exactly inverting rainfall horizontal distribution parameter.
Step 3: according to above-mentioned rainfall horizontal distribution parameter
Figure 108131DEST_PATH_IMAGE030
and the vertical distribution parameter of being finally inversed by; Horizontal distribution parameter
Figure 511748DEST_PATH_IMAGE030
and vertical distribution parameter
Figure 813416DEST_PATH_IMAGE031
are multiplied each other, obtain the three-dimensional distribution parameter of rainfall;
The three-dimensional distribution parameter
Figure 561109DEST_PATH_IMAGE042
of rainfall, its expression formula is:
Figure 131243DEST_PATH_IMAGE077
(5)
Wherein,
Figure 21839DEST_PATH_IMAGE078
is the height of rainfall area;
Figure 818894DEST_PATH_IMAGE079
is the ground rainfall speed of inverting;
Figure 99834DEST_PATH_IMAGE080
is the rainfall area height, and
Figure 273326DEST_PATH_IMAGE081
is the variable of rainfall area height.

Claims (4)

1. the inversion algorithm of the three-dimensional distribution parameter of rainfall comprises the following steps:
Step 1: adopt spaceborne wave band synthetic-aperture radar to obtain total echo data
Figure 773271DEST_PATH_IMAGE002
, be finally inversed by starting point , the total echo minimum value of radar
Figure 802407DEST_PATH_IMAGE004
of rainfall; Adopt the VIE inversion algorithm; Obtain the rainfall distribution and expression formula of interval ; Wherein, leaves the horizontal ordinate in rainfall zone for radar beam;
Figure 603321DEST_PATH_IMAGE008
is the radar beam incident angle, and
Figure DEST_PATH_IMAGE009
is the rain belt height;
Step 2: according to the rainfall distribution of above-mentioned interval
Figure 548143DEST_PATH_IMAGE010
; Be finally inversed by rainfall horizontal distribution parameter
Figure DEST_PATH_IMAGE011
, this rainfall horizontal distribution parameter
Figure 381856DEST_PATH_IMAGE011
confirms that the rainfall horizontal distribution is: the shape of the rainfall horizontal distribution of the whole rain belt of inverting, rainfall area width
Figure 819791DEST_PATH_IMAGE012
, ground rainfall speed
Figure DEST_PATH_IMAGE013
, use vertical distribution inversion algorithm obtain vertical distribution parameter
Figure 646932DEST_PATH_IMAGE014
;
Step 3: according to above-mentioned rainfall horizontal distribution parameter
Figure 853923DEST_PATH_IMAGE011
and the vertical distribution parameter
Figure 764110DEST_PATH_IMAGE014
of being finally inversed by; Horizontal distribution parameter
Figure 689341DEST_PATH_IMAGE011
and vertical distribution parameter
Figure 913649DEST_PATH_IMAGE014
are multiplied each other, obtain the three-dimensional distribution parameter of rainfall.
2. the inversion algorithm of the three-dimensional distribution parameter of a kind of rainfall according to claim 1; It is characterized in that: above-mentioned steps one described employing spaceborne
Figure 319353DEST_PATH_IMAGE001
wave band synthetic-aperture radar is obtained total echo data
Figure 338125DEST_PATH_IMAGE016
, is finally inversed by starting point
Figure DEST_PATH_IMAGE017
, the total echo minimum value of radar
Figure 547389DEST_PATH_IMAGE018
of rainfall; Adopt the VIE inversion algorithm; Obtain the rainfall distribution and expression formula
Figure 575388DEST_PATH_IMAGE020
of interval
Figure DEST_PATH_IMAGE019
; Wherein,
Figure 835599DEST_PATH_IMAGE007
leaves the horizontal ordinate in rainfall zone for radar beam;
Figure 290851DEST_PATH_IMAGE008
is the radar beam incident angle, and
Figure 925095DEST_PATH_IMAGE009
is the rain belt height;
A. inverting rainfall starting point
Figure DEST_PATH_IMAGE021
, calculate following formula:
Figure 553523DEST_PATH_IMAGE022
(3)
Wherein
Figure DEST_PATH_IMAGE023
and
Figure 855191DEST_PATH_IMAGE024
representes respectively: the mean value and the root mean square of 5 total echo points of radar at before the point, and the numerical value that satisfies formula (3) is the starting point
Figure 25885DEST_PATH_IMAGE026
that is finally inversed by rainfall;
B. the size that compares all numerical value in
Figure DEST_PATH_IMAGE027
, wherein minimum value is the total echo minimum value of radar
Figure 147425DEST_PATH_IMAGE028
;
C. the rainfall distribution and expression formula
Figure 1297DEST_PATH_IMAGE006
of computation interval
Figure 845122DEST_PATH_IMAGE019
If
Figure DEST_PATH_IMAGE029
is the piecewise function of rainfall zone horizontal ordinate
Figure 611401DEST_PATH_IMAGE025
, then the rainfall distribution and expression formula in whole rainfall zone is:
(1)
When rainfall zone horizontal ordinate
Figure 3516DEST_PATH_IMAGE025
when interval
Figure DEST_PATH_IMAGE031
, substitution (1) formula is got
Figure 778891DEST_PATH_IMAGE006
:
Figure 875023DEST_PATH_IMAGE032
(6)
Wherein,
Figure DEST_PATH_IMAGE033
axle is transverse axis position (Km), rainfall zone, and
Figure 665256DEST_PATH_IMAGE034
axle is the numerical value of the total echo of radar.
3. the inversion algorithm of the three-dimensional distribution parameter of a kind of rainfall according to claim 1; It is characterized in that: above-mentioned steps two described rainfall distributions
Figure 530444DEST_PATH_IMAGE006
according to above-mentioned interval
Figure DEST_PATH_IMAGE035
; Be finally inversed by rainfall horizontal distribution parameter , this rainfall horizontal distribution parameter
Figure 315046DEST_PATH_IMAGE011
confirms that the rainfall horizontal distribution is: the shape of the rainfall horizontal distribution of the whole rain belt of inverting, rainfall area width
Figure 830341DEST_PATH_IMAGE012
, ground rainfall speed
Figure 550035DEST_PATH_IMAGE013
, use vertical distribution inversion algorithm obtain vertical distribution parameter
Figure 410194DEST_PATH_IMAGE036
;
A. calculate derivative more arbitrarily; Be designated as
Figure DEST_PATH_IMAGE037
; if
Figure 534324DEST_PATH_IMAGE038
, then rainfall horizontal distribution is shaped as rectangle;
B. rainfall area width
Figure 702001DEST_PATH_IMAGE012
, its calculating formula is:
(4)
Wherein
Figure 182660DEST_PATH_IMAGE040
is the total echo minimum value of radar that is finally inversed by, and
Figure DEST_PATH_IMAGE041
is the starting point of rainfall;
C. the rainfall level is shaped as rectangle, and the maximal value of function
Figure 553730DEST_PATH_IMAGE006
is ground rainfall speed ;
D. be finally inversed by rainfall horizontal distribution parameter
Figure 839535DEST_PATH_IMAGE011
The rainfall horizontal distribution be shaped as rectangle:
When belonging to interval as
Figure 553413DEST_PATH_IMAGE025
, then rainfall horizontal distribution parameter expression is:
(7)
When
Figure 806988DEST_PATH_IMAGE025
does not belong to interval
Figure 621360DEST_PATH_IMAGE044
when the level of rainfall distribution parameter expression:
Figure DEST_PATH_IMAGE045
(8)
Wherein,
Figure 443823DEST_PATH_IMAGE046
is the rainfall starting point, and
Figure 507594DEST_PATH_IMAGE012
is the rainfall area width;
E. use the vertical distribution inversion algorithm to obtain vertical distribution parameter
In interval
Figure DEST_PATH_IMAGE047
Figure 454001DEST_PATH_IMAGE048
(2)
Wherein,
Figure DEST_PATH_IMAGE049
is the height of rainfall area;
Figure 447365DEST_PATH_IMAGE013
is ground rainfall speed, and
Figure 201694DEST_PATH_IMAGE050
is the variable of rainfall area height.
4. the inversion algorithm of the three-dimensional distribution parameter of a kind of rainfall according to claim 1; It is characterized in that: above-mentioned steps three is described according to above-mentioned rainfall horizontal distribution parameter
Figure 1023DEST_PATH_IMAGE011
and the vertical distribution parameter
Figure 524408DEST_PATH_IMAGE036
of being finally inversed by; Horizontal distribution parameter
Figure 954253DEST_PATH_IMAGE011
and vertical distribution parameter
Figure 740418DEST_PATH_IMAGE036
are multiplied each other, obtain the three-dimensional distribution parameter
Figure DEST_PATH_IMAGE051
of rainfall;
The three-dimensional distribution parameter
Figure 15542DEST_PATH_IMAGE051
of rainfall, its expression formula is:
Figure 659013DEST_PATH_IMAGE052
(5)
Wherein,
Figure 322075DEST_PATH_IMAGE049
is the height of rainfall area;
Figure 50997DEST_PATH_IMAGE013
is the ground rainfall speed of inverting;
Figure DEST_PATH_IMAGE053
is the rainfall area height, and
Figure 942861DEST_PATH_IMAGE050
is the variable of rainfall area height.
CN2012100326079A 2012-02-15 2012-02-15 Inversion algorithm of rainfall three-dimensional distribution parameters Pending CN102540148A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2012100326079A CN102540148A (en) 2012-02-15 2012-02-15 Inversion algorithm of rainfall three-dimensional distribution parameters

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012100326079A CN102540148A (en) 2012-02-15 2012-02-15 Inversion algorithm of rainfall three-dimensional distribution parameters

Publications (1)

Publication Number Publication Date
CN102540148A true CN102540148A (en) 2012-07-04

Family

ID=46347484

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012100326079A Pending CN102540148A (en) 2012-02-15 2012-02-15 Inversion algorithm of rainfall three-dimensional distribution parameters

Country Status (1)

Country Link
CN (1) CN102540148A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103499813A (en) * 2013-10-11 2014-01-08 邓勇 Method for distinguishing cloud and precipitation by radar precipitation factor
CN111474606A (en) * 2020-03-30 2020-07-31 中国人民解放军国防科技大学 Rainfall and water vapor comprehensive measurement device and method based on dual-frequency dual-polarization microwave link
CN112131989A (en) * 2020-09-15 2020-12-25 河海大学 Millimeter wave rain measurement model parameter obtaining method based on space rainfall data
CN113325418A (en) * 2021-07-08 2021-08-31 上海海洋大学 Rainfall inversion method based on synthetic aperture radar
CN113672864A (en) * 2021-07-15 2021-11-19 中国电波传播研究所（中国电子科技集团公司第二十二研究所） Weighted algorithm applied to rain attenuation statistical prediction for annual average rain top height

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101609172A (en) * 2009-07-14 2009-12-23 国家***第三海洋研究所 A kind of water surface rainfall intensity measuring method and device thereof based on patter of rain spectrum

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101609172A (en) * 2009-07-14 2009-12-23 国家***第三海洋研究所 A kind of water surface rainfall intensity measuring method and device thereof based on patter of rain spectrum

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
MARZANO, F.S.: "Inversion of Spaceborne X-Band Synthetic Aperture Radar Measurements for Precipitation Remote Sensing Over Land", 《GEOSCIENCE AND REMOTE SENSING, IEEE TRANSACTIONS ON》 *
谢亚楠等: "Study on the Algorithm to Retrieve Precipitation with X-Band Synthetic Aperture Radar", 《ACTA METEOROLOGICA SINICA》 *
陈捷等: "用合成孔径雷达测量降雨信杂比的方法研究", 《科学技术与工程》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103499813A (en) * 2013-10-11 2014-01-08 邓勇 Method for distinguishing cloud and precipitation by radar precipitation factor
CN111474606A (en) * 2020-03-30 2020-07-31 中国人民解放军国防科技大学 Rainfall and water vapor comprehensive measurement device and method based on dual-frequency dual-polarization microwave link
CN111474606B (en) * 2020-03-30 2022-04-05 中国人民解放军国防科技大学 Rainfall and water vapor comprehensive measurement device and method based on dual-frequency dual-polarization microwave link
CN112131989A (en) * 2020-09-15 2020-12-25 河海大学 Millimeter wave rain measurement model parameter obtaining method based on space rainfall data
CN112131989B (en) * 2020-09-15 2021-07-23 河海大学 Millimeter wave rain measurement model parameter obtaining method based on space rainfall data
CN113325418A (en) * 2021-07-08 2021-08-31 上海海洋大学 Rainfall inversion method based on synthetic aperture radar
CN113325418B (en) * 2021-07-08 2023-07-25 上海海洋大学 Rainfall inversion method based on synthetic aperture radar
CN113672864A (en) * 2021-07-15 2021-11-19 中国电波传播研究所（中国电子科技集团公司第二十二研究所） Weighted algorithm applied to rain attenuation statistical prediction for annual average rain top height
CN113672864B (en) * 2021-07-15 2023-08-15 中国电波传播研究所（中国电子科技集团公司第二十二研究所） Annual average roof height weighting algorithm applied to statistical prediction of rain attenuation

Similar Documents

Publication Publication Date Title
CN102540148A (en) Inversion algorithm of rainfall three-dimensional distribution parameters
CN110488393B (en) X-band dual-polarization weather radar quantitative measurement precipitation method and system
Rosenow et al. Vertical velocity and physical structure of generating cells and convection in the comma head region of continental winter cyclones
CN108681652A (en) A kind of vegetation-covered area Soil Moisture Inversion method based on three number of high score
Mateus et al. Can spaceborne SAR interferometry be used to study the temporal evolution of PWV?
CN103198314A (en) Remote sensing image radiation correction method
Fries et al. Catchment precipitation processes in the San Francisco valley in southern Ecuador: Combined approach using high-resolution radar images and in situ observations
CN104316930A (en) Raininess estimation method based on dual polarization Doppler weather radar detection
CN103926635A (en) Method for monitoring rain area distribution by utilization of microwave link network
KR20160057911A (en) System for deducting agricultural drought index and method therefor
CN104833974A (en) SAR imaging quick backward projection method based on image spectrum compression
CN103777186A (en) Method for calculating near-field radar echo characteristics of moving object based on beam decomposition and local irradiation
CN102707336A (en) Novel method of using A-Train series satellite data for synergetic inversion of cloud phase states and cloud parameters
CN106569226B (en) A method of have and utilizes laser radar Data Inversion Boundary Layer Height in the case of cloud
KR101174466B1 (en) System and method of all-sky cloud liquid water content production
DE112016007092T5 (en) RADAR DEVICE
Levin et al. Observation impacts on the Mid-Atlantic Bight front and cross-shelf transport in 4D-Var ocean state estimates: Part I—Multiplatform analysis
Li et al. Improving a spectral bin microphysical scheme using TRMM satellite observations
Vogelzang et al. NWP model error structure functions obtained from scatterometer winds
KR100922128B1 (en) The correction technique of weather radar reflectivity
Mitchell et al. A new snow growth model with application to radar precipitation estimates
CN113064130B (en) Method, device, storage medium and program product for determining particle spectral distribution
Shrestha et al. Spatial Variations in the diurnal pattern of precipitation over Nepal Himalayas
Holleman Bias adjustment of radar-based 3-hour precipitation accumulations
Van Dinther et al. Using the time-lag correlation function of dual-aperture scintillometer measurements to obtain the crosswind

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20120704