CN110031841A - The method and system of InSAR atmospheric delay correction based on ECMWF - Google Patents
The method and system of InSAR atmospheric delay correction based on ECMWF Download PDFInfo
- Publication number
- CN110031841A CN110031841A CN201910257685.0A CN201910257685A CN110031841A CN 110031841 A CN110031841 A CN 110031841A CN 201910257685 A CN201910257685 A CN 201910257685A CN 110031841 A CN110031841 A CN 110031841A
- Authority
- CN
- China
- Prior art keywords
- delay
- zenith
- atmosphere
- value
- point
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000012937 correction Methods 0.000 title claims abstract description 52
- 238000011160 research Methods 0.000 claims abstract description 83
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000010276 construction Methods 0.000 claims description 6
- 239000007788 liquid Substances 0.000 claims description 6
- 230000002035 prolonged effect Effects 0.000 claims description 3
- 239000005436 troposphere Substances 0.000 description 13
- 230000000694 effects Effects 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 241000208340 Araliaceae Species 0.000 description 3
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 3
- 235000003140 Panax quinquefolius Nutrition 0.000 description 3
- 235000008434 ginseng Nutrition 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000002706 hydrostatic effect Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- CCEKAJIANROZEO-UHFFFAOYSA-N sulfluramid Chemical group CCNS(=O)(=O)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)F CCEKAJIANROZEO-UHFFFAOYSA-N 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Abstract
The method and system of the invention discloses a kind of InSAR atmospheric delay correction based on ECMWF.Method includes: to calculate the tropospheric zenith atmosphere dry delay and zenith wet stack emission in research area according to the dem data of ECMWF data of weather forecast and research area；Total zenith delay that research area corresponds to the thick scale of each SAR image is obtained with zenith wet stack emission according to the dry delay of the zenith in research area；Layering interpolation model is constructed according to total zenith delay of each SAR image；Obtain the atmosphere delay phase value of every bit on interference SAR image；Atmosphere delay phase value is twined in figure from interference SAR solution and is removed, the interference pattern after obtaining atmospheric delay correction.The present invention improves atmospheric delay correction resolution ratio and calculated result is accurate.
Description
Technical field
The present invention relates to InSAR technical fields, more particularly, to a kind of InSAR atmospheric delay correction based on ECMWF
Method and system.
Background technique
InSAR (Interferometric Synthetic Aperture Radar, interfering synthetic aperture radar), is one
Kind is applied to the Radar Technology of mapping and remote sensing.It is the two width complex values shadows observed using synthetic aperture radar areal
As (image that existing amplitude has phase) data progress Coherent processing again, to obtain the technology of earth's surface elevation information.InSAR technology
There is very widely research and application in many fields, as topographic survey, earthquake, landslide, volcano, large area earth's surface are heavy
Drop, marine charting and disaster monitoring etc..
InSAR technology is by electromagnetic wave return path signal in the application, and electromagnetic wave is during atmospheric propagation, will be by right
The influence of fluid layer atmosphere especially vapor scattering, causes Electromagnetic Wave Propagation path to bend, makes radar echo signal
Phase delay, referred to as atmosphere delay.For InSAR, imaging time and visual angle are different from radar twice, and atmospheric conditions
Temporally and spatially there is variation in (relative humidity, temperature, air pressure etc.), Electromagnetic Wave Propagation is experienced big when being imaged twice
Gas bar part is just different, and atmospheric phase delay can have an impact interference and differential interferometry phase, and influence the essence of InSAR technology
Degree.
Currently, the atmospheric delay correction method in InSAR have based on MERIS (a kind of Moderate Imaging Spectroradiomete),
A kind of method of MODIS (Moderate Imaging Spectroradiomete) steam product, method based on GPS data are based on Mesoscale
The method etc. of forecasting model WRF (Weather Research and Forecast Model).Wherein, MERIS, MODIS are obtained
Moisture content data are influenced by cloud, therefore in the region for having cloud, the moisture content inaccuracy of the data, while this method is counted
According to the limitation of temporal resolution, the receiving time of data be it is fixed, the time of data time and practical interference pattern has deviation,
Influence calculated result；Using GPS data carry out atmospheric delay correction precision it is higher, but the GPS station data density in China compared with
It is low, therefore obtained atmospheric correction image resolution ratio is relatively low.And the data resolution of WRF model is also relatively low.
Therefore it provides a kind of method and system of the InSAR atmospheric delay correction based on ECMWF, is that this field urgently solves
Certainly the technical issues of.
Summary of the invention
In view of this, the method and system of the present invention provides a kind of InSAR atmospheric delay correction based on ECMWF, mentions
High atmospheric delay correction resolution ratio and calculated result is accurate.
The side of in order to solve the above technical problem, the present invention provides a kind of InSAR atmospheric delay correction based on ECMWF
Method characterized by comprising
According to the dem data of ECMWF data of weather forecast and research area, the tropospheric zenith atmosphere in calculating research area is dry to be prolonged
Slow and zenith wet stack emission；
Obtaining studying area according to the dry delay of the zenith in research area and zenith wet stack emission, to correspond to each SAR image thick
Total zenith delay of scale；
Layering interpolation model is constructed according to total zenith delay of each SAR image, is layered interpolation model are as follows: Δ L=S (h)+
T (x)+w, wherein △ L is total zenith delay difference of two scape SAR images, and S is the hierarchical delay part of atmosphere delay, and T is big
The turbulent flow decay part of gas delay, w are residual retardance value, and x is the certain point in SAR image, and h is the elevation at corresponding x point
Value；
The hierarchical delay value S of each point in research area's internal interference SAR image is calculated using interpolation method, comprising:
Iteration function S (h)=qexp {-b (h-hmin)/(hmax-hmin) find the optimum parameter value of q and b, q and b be to
The index parameters of interpolation area, h are the height value of interpolation point, hmaxFor the elevation maximum value in SAR image region, hminFor SAR
The elevation minimum value of image-region；
The hierarchical delay value S of each point in research area is obtained using interpolation method by the optimum parameter value of q and b；
The hierarchical delay value S for studying each point in area is separated according to layering interpolation model, is inserted using anti-distance weighting
Value method obtains the turbulent flow length of delay T and residual retardance value w of each point in research area；
According to the hierarchical delay value S of each point and turbulent flow length of delay T of each point in research area, obtain every on interference SAR image
The atmosphere delay phase value of a bit；
Atmosphere delay phase value is twined in figure from interference SAR solution and is removed, the interference pattern after obtaining atmospheric delay correction.
Optionally, according to the dem data of ECMWF data of weather forecast and research area, the tropospheric zenith in research area is calculated
Atmosphere dry delay and zenith wet stack emission, comprising:
Using formula
Calculate the dry delay of zenith, wherein ZHD is the dry delay of zenith, and unit is mm；Ps is earth's surface air pressure, unit hpa；Ground point dimension is represented, unit is degree；H represents ground point height, unit km；
Using formulaCalculate zenith wet stack emission, wherein ZWD is
Zenith wet stack emission, k' and k are refraction constant, ρwIt is liquid water density, TMIt is tropospheric weighted mean, RvIt is special
Fixed atmosphere constant, PWV are precipitable water vapor content.
Optionally, according to the dem data of ECMWF data of weather forecast and research area, the tropospheric zenith in research area is calculated
Atmosphere dry delay and zenith wet stack emission, comprising:
Using formula
Calculate the dry delay of zenith, wherein ZHD is the dry delay of zenith, and unit is mm；Ps is earth's surface air pressure, unit hpa；Ground point dimension is represented, unit is degree；H represents ground point height, unit km；
Zenith wet stack emission is calculated using formula ZWD=Π PWV, wherein ZWD is zenith wet stack emission, and PWV is
Precipitable water vapor content, Π are scale factor.
Optionally, 6.0≤Π≤6.5.
Optionally, each SAR in research area is obtained according to the dry delay of zenith in research area and zenith wet stack emission
Total zenith delay of image, comprising: ZTD=ZHD+ZWD, ZTD are total zenith delay.
The system of the present invention also provides a kind of InSAR atmospheric delay correction based on ECMWF characterized by comprising number
Prolong according to module, total zenith delay computing module, interpolation model computing module, atmosphere delay phase value generation module and atmosphere is obtained
Correction module late；Wherein,
Data acquisition module is connected with total zenith delay computing module, for obtaining ECMWF data of weather forecast and grinding
Study carefully the dem data in area, and the dem data of ECMWF data of weather forecast and research area is sent to total zenith delay computing module；
Total zenith delay computing module is connected with difference model computing module, for according to ECMWF data of weather forecast
With the dem data in research area, the tropospheric zenith atmosphere dry delay and zenith wet stack emission in research area are calculated, and according to grinding
The dry delay of zenith and zenith wet stack emission for studying carefully area, which obtain studying area and correspond to total zenith of the thick scale of each SAR image, prolongs
Late, and by the total zenith delay for studying each SAR image in area it is sent to interpolation model computing module；
Interpolation model computing module is connected with atmosphere delay phase value generation module, for according to each SAR image
Total zenith delay building layering interpolation model, and the hierarchical delay value S for studying each point in area is calculated according to layering interpolation model
With the turbulent flow length of delay T of each point, comprising: model construction unit, hierarchical delay value S computing unit and turbulent flow length of delay T calculate single
Member, wherein
Model construction unit is layered interpolation model are as follows: Δ L=S (h)+T (x)+w, △ L for constructing layering interpolation model
For total zenith delay difference of two scape SAR images, S is the hierarchical delay part of atmosphere delay, and T is that the turbulent flow of atmosphere delay postpones
Part, w are residual retardance value, and x is the certain point in SAR image, and h is the height value at corresponding x point；
Hierarchical delay value S computing unit is used to calculate point of each point in research area's internal interference SAR image using interpolation method
Layer length of delay S, first according to iteration function S (h)=qexp {-b (h-hmin)/(hmax-hmin) find q and b optimal parameter
Value, q and b are the index parameters in interpolation region, and h is the height value of interpolation point, hmaxIt is maximum for the elevation in SAR image region
Value, hminFor the elevation minimum value in SAR image region；Then it is studied by the optimum parameter value of q and b using interpolation method
The hierarchical delay value S of each point in area, and the hierarchical delay value S for studying each point in area is sent to atmosphere delay phase value and generates mould
Block；
Turbulent flow length of delay T computing unit is used to be divided the hierarchical delay value S for studying each point in area according to layering interpolation model
It separates out and, the turbulent flow length of delay T and residual retardance value w of each point in research area are obtained using anti-distance weighting interpolation method, and will
The turbulent flow length of delay T of each point is sent to atmosphere delay phase value generation module in research area；
Atmosphere delay phase value generation module is connected with atmospheric delay correction module, for according to each point in research area
The turbulent flow length of delay T of hierarchical delay value S and each point, obtain the atmosphere delay phase value of every bit on interference SAR image, and send
Give atmospheric delay correction module；
Atmospheric delay correction module is removed for twining in figure atmosphere delay phase value from interference SAR solution, obtains atmosphere and prolongs
Interference pattern after correction late.
Optionally, total zenith delay computing module includes atmosphere dry delay computing unit and atmospheric moisture delay calculation unit,
Wherein,
Atmosphere dry delay computing unit is used for the dem data according to ECMWF data of weather forecast and research area, using formulaIt is dry to calculate zenith
Delay, wherein ZHD is the dry delay of zenith, and unit is mm；Ps is earth's surface air pressure, unit hpa；Represent ground point dimension
Degree, unit are degree；H represents ground point height, unit km；
Atmospheric moisture delay calculation unit is used for the dem data according to ECMWF data of weather forecast and research area, using formulaCalculating zenith wet stack emission, wherein ZWD is zenith wet stack emission,
K' and k is refraction constant, ρwIt is liquid water density, TMIt is tropospheric weighted mean, RvIt is specific atmosphere constant, PWV
For precipitable water vapor content.
Optionally, total zenith delay computing module includes atmosphere dry delay computing unit and atmospheric moisture delay calculation unit,
Wherein,
Atmosphere dry delay computing unit is used for the dem data according to ECMWF data of weather forecast and research area, using formulaIt is dry to calculate zenith
Delay, wherein ZHD is the dry delay of zenith, and unit is mm；Ps is earth's surface air pressure, unit hpa；Represent ground point dimension
Degree, unit are degree；H represents ground point height, unit km；
Atmospheric moisture delay calculation unit is used for the dem data according to ECMWF data of weather forecast and research area, using formula
ZWD=Π PWV calculates zenith wet stack emission, wherein ZWD is zenith wet stack emission, and PWV is precipitable water vapor content, Π
For scale factor.
Optionally, total zenith delay computing module further includes total delay computing module；Total delay computing module is used for basis
The dry delay of the zenith in research area and zenith wet stack emission obtain total zenith delay of each SAR image in research area,
In, ZTD=ZHD+ZWD, ZTD are total zenith delay.
Compared with prior art, the method and system of the InSAR atmospheric delay correction provided by the invention based on ECMWF,
At least realize it is following the utility model has the advantages that
The method and system of troposphere atmospheric delay correction provided by the invention carries out InSAR atmosphere using ECMWF data
Delay computing, ECMWF data used, forecast data time interval are 3 hours, are updated once every 6 hours, and timeliness is strong,
Spatial resolution reaches as high as 0.125 °, is the weather forecast data of current highest resolution.The present invention utilizes effective in real time
Meteorologic parameter calculates accurate high-precision atmosphere delay.In addition, the present invention considers turbulence effect and the layering of atmosphere simultaneously
The regularity of distribution of effect carries out interpolation calculation to delay in spatial altitude distribution using layering interpolation model, improves resolution
Rate, obtained atmosphere delay calculated result accuracy are higher.
By referring to the drawings to the detailed description of exemplary embodiment of the present invention, other feature of the invention and its
Advantage will become apparent.
Detailed description of the invention
It is combined in the description and the attached drawing for constituting part of specification shows the embodiment of the present invention, and even
With its explanation together principle for explaining the present invention.
Fig. 1 is the method flow diagram of the InSAR troposphere atmospheric delay correction provided by the invention based on ECMWF；
Fig. 2 is the system block diagram of the InSAR troposphere atmospheric delay correction provided in an embodiment of the present invention based on ECMWF.
Specific embodiment
Carry out the various exemplary embodiments of detailed description of the present invention now with reference to attached drawing.It should also be noted that unless in addition having
Body explanation, the unlimited system of component and the positioned opposite of step, numerical expression and the numerical value otherwise illustrated in these embodiments is originally
The range of invention.
Be to the description only actually of at least one exemplary embodiment below it is illustrative, never as to the present invention
And its application or any restrictions used.
Technology, method and apparatus known to person of ordinary skill in the relevant may be not discussed in detail, but suitable
In the case of, technology, method and apparatus should be considered as part of specification.
It is shown here and discuss all examples in, any occurrence should be construed as merely illustratively, without
It is as limitation.Therefore, other examples of exemplary embodiment can have different values.
It should also be noted that similar label and letter indicate similar terms in following attached drawing, therefore, once a certain Xiang Yi
It is defined in a attached drawing, then in subsequent attached drawing does not need that it is further discussed.
European Center for Medium Weather Forecasting (European Centre for Medium-Range Weather
Forecasts abbreviation ECMWF) it is the international organizations including 24 EU member countries, the current whole world is taken the course of its own
International weather forecast research and operational agency.The current real-time update of third generation ERA-Interim of ECMWF analysis of data again,
Using four-dimensional variation analysis, which includes 60 layers in vertical direction, and top arrival 0.1hPa is provided four times per day
Assimilation data include the various atmospheric parameters such as temperature, air pressure, humidity.Inventor considers that the data based on ECMWF utilize it
Higher spatial resolution and temporal resolution correct InSAR troposphere atmosphere delay, have improved atmospheric delay correction
Resolution ratio and calculated result accuracy.
Fig. 1 is the method flow diagram of the InSAR troposphere atmospheric delay correction provided by the invention based on ECMWF.Such as Fig. 1
It is shown, the method include that
Step S101: according to the dem data of ECMWF data of weather forecast and research area, the troposphere day in research area is calculated
Push up atmosphere dry delay and zenith wet stack emission；Wherein, ECMWF data of weather forecast includes air pressure, temperature, humidity, longitude and latitude
Etc. parameters.Optionally, method of the invention selects ECMWF data of weather forecast intermediate-resolution for 0.125 ° of weather in application
Forecast data.
The dry delay (i.e. zenith hydrostatic delay, Zenith Hydrostatic Delay, ZHD) of zenith, can be by ground
The data such as point latitude, elevation (height) and earth's surface air pressure are calculated, and optionally, can use formulaCalculate that zenith is dry prolongs
Late, wherein ZHD is the dry delay of zenith, and unit is mm；Ps is earth's surface air pressure, unit hpa；Ground point dimension is represented,
Unit is degree；H represents ground point height, unit km.
Zenith wet stack emission (Zenith Wet Delay, ZWD), by perceptible water vapor can calculate
It arrives, wherein perceptible water vapor is the physical quantity of water-air regime in atmosphere, and it is straight can to monitor system by existing space base
It obtains.
Optionally, formula can be usedCalculate that zenith is wet prolongs
Late, wherein ZWD is zenith wet stack emission and k is refraction constant, is liquid water density, is tropospheric weighted average temperature
Degree, is specific atmosphere constant, and PWV is precipitable water vapor content.
Optionally, zenith wet stack emission can be calculated using formula ZWD=Π PWV, wherein ZWD is zenith
Wet stack emission, PWV are precipitable water vapor content, and Π is scale factor.Further, 6.0≤Π≤6.5.For between ZWD and PWV
When being converted roughly, it is 6.2 that Π, which can choose an average conversion factor,.
Step S102: according to the dry delay of the zenith in research area and zenith wet stack emission obtain studying area correspond to it is each
Total zenith delay of the thick size of SAR image；Troposphere atmosphere delay is mainly made of the dry delay of zenith and Zenith wet delay.We
In method, definition ZTD is total zenith delay, wherein ZTD=ZHD+ZWD.
Step S103: layering interpolation model (Stratified is constructed according to total zenith delay of each SAR image
Decomposition Interpolated model, SDI model), it is layered interpolation model are as follows: Δ L=S (h)+T (x)+w,
In, △ L is total zenith delay difference of two scape SAR images, i.e. atmosphere delay value on interference SAR image, and S is atmosphere delay
Hierarchical delay part, T are the turbulent flow decay part of atmosphere delay, and w is residual retardance value, and x is the certain point in SAR image, and h is
Height value at corresponding x point.
Inventor considers that the resolution ratio of 0.125 ° of data of weather forecast in ECMWF is higher, closer to the resolution of SAR image
Rate, but the resolution ratio of SAR image is generally several meters or tens meters, the resolution ratio for individually ECMWF data being leaned on to calculate is also
It is not nearly enough, therefore further construct layering interpolation model.In view of changing due to weather condition and strong steam, will lead to
The turbulent fluctuation of temperature and humidity, vertical wind shear power or the strong convection effect generated due to the thin turbosphere in cumulus are one
A random variable quantity makes air index show spatially heterogeneous, so as to cause local phase gradient, turbulent flow is caused to be imitated
It answers.In addition, earth atmosphere section shows different degrees of layering i.e. layering effect in extensive vertical dimension, in consideration office
When portion's scale (such as interference pattern), this layering may be by the severe jamming of troposphere turbulence.Therefore inventor is prolonging
Slow result considers layering effect and turbulence effect when calculating, and above-mentioned layering interpolation model is constructed, to improve atmospheric delay correction
Accuracy.
Step S104: the hierarchical delay value S of each point in research area's internal interference SAR image is calculated using interpolation method, comprising:
Iteration function S (h)=qexp {-b (h-hmin)/(hmax-hmin) find the optimum parameter value of q and b, q and b be to
The index parameters of interpolation area, h are the height value of interpolation point, hmaxFor the elevation maximum value in SAR image region, hminFor SAR
The elevation minimum value of image-region；
Step S105: the hierarchical delay value of each point in research area is obtained using interpolation method by the optimum parameter value of q and b
S；In this step, optimal parameter q and b are found by iteration function, it is possible thereby to by the layering of different height in survey region
Delay component is separated, i.e., obtains the hierarchical delay value S of each point in survey region using interpolation method.
Step S106: will study the hierarchical delay value S of each point in area according to layering interpolation model and separate, using instead away from
The turbulent flow length of delay T and residual retardance value w of each point in research area are obtained from Power Interpolation method；Carrying out turbulent flow length of delay T meter
When calculation, inventor is suitable for the case where gradually changing in wide area in view of batten or bilinear interpolation method, but uncomfortable
Conjunction varies widely in shorter horizontal distance, this makes above two method be not suitable for extreme weather condition, because
This is closed present invention employs the turbulent flow length of delay T that the interpolation method of IDW (anti-distance weighting) calculates each point in research area with realizing
Reason calculates turbulent flow length of delay T.
Step S107: according to the hierarchical delay value S of each point and turbulent flow length of delay T of each point in research area, interference SAR is obtained
The atmosphere delay phase value of every bit on image；
Step S108: atmosphere delay phase value being twined in figure from interference SAR solution and is removed, dry after obtaining atmospheric delay correction
Relate to figure.
The method of troposphere atmospheric delay correction provided by the invention carries out InSAR atmosphere delay meter using ECMWF data
It calculates, ECMWF data used, forecast data time interval is 3 hours, is updated once every 6 hours, and timeliness is strong, space point
Resolution reaches as high as 0.125 °, is the weather forecast data of current highest resolution.The present invention utilizes meteorological ginseng effective in real time
Number calculates accurate high-precision atmosphere delay.In addition, the present invention considers the turbulence effect and layering effect of atmosphere simultaneously
The regularity of distribution carries out interpolation calculation to delay in spatial altitude distribution using layering interpolation model, improves resolution ratio, obtain
Atmosphere delay calculated result accuracy it is higher.
The system of the present invention also provides a kind of InSAR troposphere atmospheric delay correction based on ECMWF, Fig. 2 are the present invention
The system block diagram for the InSAR troposphere atmospheric delay correction based on ECMWF that embodiment provides.As shown in Fig. 2, the present invention provides
System 200 include: data acquisition module 201, total zenith delay computing module 202, interpolation model computing module 203, atmosphere
Postpone phase value generation module 204 and atmospheric delay correction module 205；Wherein,
Data acquisition module 201 is connected with total zenith delay computing module 202, for obtaining ECMWF weather forecast number
Total zenith delay meter is sent to according to the dem data with research area, and by the dem data of ECMWF data of weather forecast and research area
Calculate module 202；
Total zenith delay computing module 202 is connected with interpolation model computing module 203, for pre- according to ECMWF weather
The dem data of count off evidence and research area calculates the tropospheric zenith atmosphere dry delay and zenith wet stack emission in research area, and
Research area, which is obtained, with zenith wet stack emission according to the dry delay of the zenith in research area corresponds to the total of the thick scale of each SAR image
Zenith delay, and the total zenith delay for studying each SAR image in area is sent to interpolation model computing module 203；
Optionally, total zenith delay computing module 202 includes atmosphere dry delay computing unit, atmospheric moisture delay calculation unit
With total delay computing module；Wherein,
Atmosphere dry delay computing unit is used for the dem data according to ECMWF data of weather forecast and research area, using formulaIt is dry to calculate zenith
Delay, wherein ZHD is the dry delay of zenith, and unit is mm；Ps is earth's surface air pressure, unit hpa；Represent ground point dimension
Degree, unit are degree；H represents ground point height, unit km；
Atmospheric moisture delay calculation unit is used for the dem data according to ECMWF data of weather forecast and research area.
In one embodiment, formula can be usedCalculate zenith
Atmosphere wet stack emission, wherein ZWD is zenith wet stack emission, and k' and k are refraction constant, ρwIt is liquid water density, TMIt is troposphere
Weighted mean, RvIt is specific atmosphere constant, PWV is precipitable water vapor content.
In another embodiment, zenith wet stack emission can be calculated using formula ZWD=Π PWV, wherein
ZWD is zenith wet stack emission, and PWV is precipitable water vapor content, and Π is scale factor.
Total delay computing module is used to be studied according to the dry delay of zenith in research area and zenith wet stack emission
Total zenith delay of each SAR image in area, wherein ZTD=ZHD+ZWD, ZTD are total zenith delay.
Interpolation model computing module 203 is connected with atmosphere delay phase value generation module 204, for according to each SAR
Total zenith delay building layering interpolation model of image, and prolonged according to the layering that layering interpolation model calculates each point in research area
The turbulent flow length of delay T of slow value S and each point, comprising: model construction unit 2031, hierarchical delay value S computing unit 2032 and turbulent flow
Length of delay T computing unit 2033, wherein
Model construction unit 2031 is layered interpolation model for constructing layering interpolation model are as follows: Δ L=S (h)+T (x)+w,
△ L is total zenith delay difference of two scape SAR images, and S is the hierarchical delay part of atmosphere delay, and T is the turbulent flow of atmosphere delay
Decay part, w are residual retardance value, and x is the certain point in SAR image, and h is the height value at corresponding x point；
Hierarchical delay value S computing unit 2032 is used to calculate each point in research area's internal interference SAR image using interpolation method
Hierarchical delay value S, first according to iteration function S (h)=qexp {-b (h-hmin)/(hmax-hmin) find q and b best ginseng
Numerical value, q and b are the index parameters in interpolation region, and h is the height value of interpolation point, hmaxFor SAR image region elevation most
Big value, hminFor the elevation minimum value in SAR image region；Then it is ground by the optimum parameter value of q and b using interpolation method
Study carefully the hierarchical delay value S of each point in area, and the hierarchical delay value S for studying each point in area is sent to atmosphere delay phase value and is generated
Module 204；
Turbulent flow length of delay T computing unit 2033 is used to study the hierarchical delay value of each point in area according to layering interpolation model
S is separated, and obtains the turbulent flow length of delay T and residual retardance value w of each point in research area using anti-distance weighting interpolation method, and
The turbulent flow length of delay T for studying each point in area is sent to atmosphere delay phase value generation module 204；
Atmosphere delay phase value generation module 204 is connected with atmospheric delay correction module 205, for according in research area
The hierarchical delay value S of each point and turbulent flow length of delay T of each point, obtains the atmosphere delay phase value of every bit on interference SAR image,
And it is sent to atmospheric delay correction module 205；
Atmospheric delay correction module 205 is removed for twining in figure atmosphere delay phase value from interference SAR solution, obtains big
Interference pattern after gas delay correction.
Through the foregoing embodiment it is found that the method for the InSAR atmospheric delay correction provided by the invention based on ECMWF and being
System, at least realize it is following the utility model has the advantages that
The method of troposphere atmospheric delay correction provided by the invention carries out InSAR atmosphere delay meter using ECMWF data
It calculates, ECMWF data used, forecast data time interval is 3 hours, is updated once every 6 hours, and timeliness is strong, space point
Resolution reaches as high as 0.125 °, is the weather forecast data of current highest resolution.The present invention utilizes meteorological ginseng effective in real time
Number calculates accurate high-precision atmosphere delay.In addition, the present invention considers the turbulence effect and layering effect of atmosphere simultaneously
The regularity of distribution carries out interpolation calculation to delay in spatial altitude distribution using layering interpolation model, improves resolution ratio, obtain
Atmosphere delay calculated result accuracy it is higher.
Although some specific embodiments of the invention are described in detail by example, the skill of this field
Art personnel it should be understood that example above merely to being illustrated, the range being not intended to be limiting of the invention.The skill of this field
Art personnel are it should be understood that can without departing from the scope and spirit of the present invention modify to above embodiments.This hair
Bright range is defined by the following claims.
Claims (9)
1. a kind of method of the InSAR atmospheric delay correction based on ECMWF characterized by comprising
According to ECMWF data of weather forecast and research area dem data, calculate research area tropospheric zenith atmosphere dry delay and
Zenith wet stack emission；
Research area, which is obtained, with zenith wet stack emission according to the dry delay of the zenith in research area corresponds to the thick scale of each SAR image
Total zenith delay；
Layering interpolation model is constructed according to total zenith delay of each SAR image, is layered interpolation model are as follows: Δ L=S (h)+T (x)
+ w, wherein △ L is total zenith delay difference of two scape SAR images, and S is the hierarchical delay part of atmosphere delay, and T prolongs for atmosphere
Slow turbulent flow decay part, w are residual retardance value, and x is the certain point in SAR image, and h is the height value at corresponding x point；
The hierarchical delay value S of each point in research area's internal interference SAR image is calculated using interpolation method, comprising:
Iteration function S (h)=qexp {-b (h-hmin)/(hmax-hmin) optimum parameter value of q and b is found, q and b are interpolation
The index parameters in region, h are the height value of interpolation point, hmaxFor the elevation maximum value in SAR image region, hminFor SAR image
The elevation minimum value in region；
The hierarchical delay value S of each point in research area is obtained using interpolation method by the optimum parameter value of q and b；
The hierarchical delay value S for studying each point in area is separated according to layering interpolation model, using anti-distance weighting interpolation side
Method obtains the turbulent flow length of delay T and residual retardance value w of each point in research area；
According to the hierarchical delay value S of each point and turbulent flow length of delay T of each point in research area, every bit on interference SAR image is obtained
Atmosphere delay phase value；
Atmosphere delay phase value is twined in figure from interference SAR solution and is removed, the interference pattern after obtaining atmospheric delay correction.
2. the method for the InSAR atmospheric delay correction according to claim 1 based on ECMWF, which is characterized in that according to ECMWF
The dem data of data of weather forecast and research area calculates the tropospheric zenith atmosphere dry delay in research area and zenith is wet prolongs
Late, comprising:
Using formulaMeter
Calculate the dry delay of zenith, wherein ZHD is the dry delay of zenith, and unit is mm；Ps is earth's surface air pressure, unit hpa；
Ground point dimension is represented, unit is degree；H represents ground point height, unit km；
Using formulaCalculate zenith wet stack emission, wherein ZWD is zenith
Atmosphere wet stack emission, k' and k are refraction constant, ρwIt is liquid water density, TMIt is tropospheric weighted mean, RvIt is specific
Atmosphere constant, PWV are precipitable water vapor content.
3. the method for the InSAR atmospheric delay correction according to claim 1 based on ECMWF, which is characterized in that according to ECMWF
The dem data of data of weather forecast and research area calculates the tropospheric zenith atmosphere dry delay in research area and zenith is wet prolongs
Late, comprising:
Using formulaMeter
Calculate the dry delay of zenith, wherein ZHD is the dry delay of zenith, and unit is mm；Ps is earth's surface air pressure, unit hpa；
Ground point dimension is represented, unit is degree；H represents ground point height, unit km；
Zenith wet stack emission is calculated using formula ZWD=Π PWV, wherein ZWD is zenith wet stack emission, and PWV is that can drop
Moisture content, Π are scale factor.
4. the method for the InSAR atmospheric delay correction according to claim 3 based on ECMWF, which is characterized in that 6.0≤Π≤
6.5。
5. according to the method for the InSAR atmospheric delay correction based on ECMWF of Claims 2 or 3, which is characterized in that according to grinding
The dry delay of zenith and zenith wet stack emission of studying carefully area obtain total zenith delay of each SAR image in research area, comprising:
ZTD=ZHD+ZWD, ZTD are total zenith delay.
6. a kind of system of the InSAR atmospheric delay correction based on ECMWF characterized by comprising data acquisition module, total
Zenith delay computing module, interpolation model computing module, atmosphere delay phase value generation module and atmospheric delay correction module；Its
In,
Data acquisition module is connected with total zenith delay computing module, for obtaining ECMWF data of weather forecast and research area
Dem data, and by ECMWF data of weather forecast and research area dem data be sent to total zenith delay computing module；
Total zenith delay computing module is connected with difference model computing module, for according to ECMWF data of weather forecast with grind
The dem data for studying carefully area calculates the tropospheric zenith atmosphere dry delay and zenith wet stack emission in research area, and according to research area
The dry delay of zenith and zenith wet stack emission obtain studying total zenith delay that area corresponds to the thick scale of each SAR image,
And the total zenith delay for studying each SAR image in area is sent to interpolation model computing module；
Interpolation model computing module is connected with atmosphere delay phase value generation module, for total day according to each SAR image
Zenith delay building layering interpolation model, and the hierarchical delay value S of each point and each is calculated in research area according to layering interpolation model
The turbulent flow length of delay T of point, comprising: model construction unit, hierarchical delay value S computing unit and turbulent flow length of delay T computing unit,
In,
Model construction unit is layered interpolation model for constructing layering interpolation model are as follows: Δ L=S (h)+T (x)+w, △ L is two
Total zenith delay difference of scape SAR image, S are the hierarchical delay part of atmosphere delay, and T is the turbulent flow delay portion of atmosphere delay
Point, w is residual retardance value, and x is the certain point in SAR image, and h is the height value at corresponding x point；
The layering that hierarchical delay value S computing unit is used to calculate each point in research area's internal interference SAR image using interpolation method is prolonged
Slow value S, first according to iteration function S (h)=qexp {-b (h-hmin)/(hmax-hmin) find q and b optimum parameter value, q and
B is the index parameters in interpolation region, and h is the height value of interpolation point, hmaxFor the elevation maximum value in SAR image region, hmin
For the elevation minimum value in SAR image region；Then it is obtained by the optimum parameter value of q and b using interpolation method each in research area
The hierarchical delay value S of point, and the hierarchical delay value S for studying each point in area is sent to atmosphere delay phase value generation module；
Turbulent flow length of delay T computing unit is used to be isolated the hierarchical delay value S for studying each point in area according to layering interpolation model
Come, the turbulent flow length of delay T and residual retardance value w of each point in research area is obtained using anti-distance weighting interpolation method, and will research
The turbulent flow length of delay T of each point is sent to atmosphere delay phase value generation module in area；
Atmosphere delay phase value generation module is connected with atmospheric delay correction module, for the layering according to each point in research area
The turbulent flow length of delay T of length of delay S and each point, obtain the atmosphere delay phase value of every bit on interference SAR image, and are sent to big
Gas delay correction module；
Atmospheric delay correction module is removed for twining in figure atmosphere delay phase value from interference SAR solution, obtains atmosphere delay and changes
Interference pattern after just.
7. the system of the InSAR atmospheric delay correction according to claim 6 based on ECMWF, which is characterized in that total zenith delay
Computing module includes atmosphere dry delay computing unit and atmospheric moisture delay calculation unit, wherein
Atmosphere dry delay computing unit is used for the dem data according to ECMWF data of weather forecast and research area, using formulaIt is dry to calculate zenith
Delay, wherein ZHD is the dry delay of zenith, and unit is mm；Ps is earth's surface air pressure, unit hpa；Represent ground point dimension
Degree, unit are degree；H represents ground point height, unit km；
Atmospheric moisture delay calculation unit is used for the dem data according to ECMWF data of weather forecast and research area, using formulaCalculating zenith wet stack emission, wherein ZWD is zenith wet stack emission,
K' and k is refraction constant, ρwIt is liquid water density, TMIt is tropospheric weighted mean, RvIt is specific atmosphere constant, PWV
For precipitable water vapor content.
8. the system of the InSAR atmospheric delay correction according to claim 6 based on ECMWF, which is characterized in that total zenith delay
Computing module includes atmosphere dry delay computing unit and atmospheric moisture delay calculation unit, wherein
Atmosphere dry delay computing unit is used for the dem data according to ECMWF data of weather forecast and research area, using formulaIt is dry to calculate zenith
Delay, wherein ZHD is the dry delay of zenith, and unit is mm；Ps is earth's surface air pressure, unit hpa；Represent ground point dimension
Degree, unit are degree；H represents ground point height, unit km；
Atmospheric moisture delay calculation unit is used for the dem data according to ECMWF data of weather forecast and research area, using formula ZWD
=Π PWV calculate zenith wet stack emission, wherein ZWD be zenith wet stack emission, PWV be precipitable water vapor content, Π be than
The example factor.
9. according to the system of the InSAR atmospheric delay correction based on ECMWF of claim 7 or 8, which is characterized in that total zenith
Delay computing module further includes total delay computing module；Total delay computing module is used for the dry delay of zenith according to research area
Total zenith delay of each SAR image in research area is obtained with zenith wet stack emission, wherein ZTD=ZHD+ZWD, ZTD are total
Zenith delay.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910257685.0A CN110031841B (en) | 2019-04-01 | 2019-04-01 | ECMWF-based InSAR (interferometric synthetic aperture radar) atmospheric delay correction method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910257685.0A CN110031841B (en) | 2019-04-01 | 2019-04-01 | ECMWF-based InSAR (interferometric synthetic aperture radar) atmospheric delay correction method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110031841A true CN110031841A (en) | 2019-07-19 |
| CN110031841B CN110031841B (en) | 2021-07-23 |
Family
ID=67237134
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910257685.0A Expired - Fee Related CN110031841B (en) | 2019-04-01 | 2019-04-01 | ECMWF-based InSAR (interferometric synthetic aperture radar) atmospheric delay correction method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110031841B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112711022A (en) * | 2020-12-18 | 2021-04-27 | 中国矿业大学 | GNSS chromatography-assisted InSAR (interferometric synthetic aperture radar) atmospheric delay correction method |
| CN112816983A (en) * | 2021-01-06 | 2021-05-18 | 中南大学 | Time sequence InSAR turbulence atmospheric delay correction method based on optimized interferogram set |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020138194A1 (en) * | 1999-12-22 | 2002-09-26 | Lorraine Flynn | System and method for estimating aircraft flight delay |
| CN103885046A (en) * | 2012-12-20 | 2014-06-25 | 河南省电力勘测设计院 | InSAR atmosphere delay correction method based on GPS |
| CN105629263A (en) * | 2015-12-21 | 2016-06-01 | 广州中海达卫星导航技术股份有限公司 | Troposphere atmosphere delay error correction method and correction system |
| CN105842692A (en) * | 2016-03-17 | 2016-08-10 | 中国科学院遥感与数字地球研究所 | Atmospheric correction method during INSAR measurement |
| CN106526577A (en) * | 2016-10-09 | 2017-03-22 | 中国船舶重工集团公司第七〇五研究所 | Array shape estimation method using cooperation sound source information |
| CN107367716A (en) * | 2017-07-04 | 2017-11-21 | 武汉大学 | A kind of high-precision satellite-borne SAR geometric calibration method |
| CN107389029A (en) * | 2017-08-24 | 2017-11-24 | 北京市水文地质工程地质大队 | A kind of surface subsidence integrated monitor method based on the fusion of multi-source monitoring technology |
| CN107839581A (en) * | 2017-09-05 | 2018-03-27 | 鲁东大学 | A kind of automatic control feedback compartment and method based on temperature field three-dimensional imaging |
| JP2018159668A (en) * | 2017-03-23 | 2018-10-11 | 株式会社デンソー | On-vehicle device |
-
2019
- 2019-04-01 CN CN201910257685.0A patent/CN110031841B/en not_active Expired - Fee Related
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020138194A1 (en) * | 1999-12-22 | 2002-09-26 | Lorraine Flynn | System and method for estimating aircraft flight delay |
| CN103885046A (en) * | 2012-12-20 | 2014-06-25 | 河南省电力勘测设计院 | InSAR atmosphere delay correction method based on GPS |
| CN105629263A (en) * | 2015-12-21 | 2016-06-01 | 广州中海达卫星导航技术股份有限公司 | Troposphere atmosphere delay error correction method and correction system |
| CN105842692A (en) * | 2016-03-17 | 2016-08-10 | 中国科学院遥感与数字地球研究所 | Atmospheric correction method during INSAR measurement |
| CN106526577A (en) * | 2016-10-09 | 2017-03-22 | 中国船舶重工集团公司第七〇五研究所 | Array shape estimation method using cooperation sound source information |
| JP2018159668A (en) * | 2017-03-23 | 2018-10-11 | 株式会社デンソー | On-vehicle device |
| CN107367716A (en) * | 2017-07-04 | 2017-11-21 | 武汉大学 | A kind of high-precision satellite-borne SAR geometric calibration method |
| CN107389029A (en) * | 2017-08-24 | 2017-11-24 | 北京市水文地质工程地质大队 | A kind of surface subsidence integrated monitor method based on the fusion of multi-source monitoring technology |
| CN107839581A (en) * | 2017-09-05 | 2018-03-27 | 鲁东大学 | A kind of automatic control feedback compartment and method based on temperature field three-dimensional imaging |
Non-Patent Citations (4)
| Title |
|---|
| PRATS P等: ""Interpolation-free Coregistration and Phase-Correction of Airborne SAR Interferograms"", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
| 占文俊 等: ""一种INSAR大气相位建模与估计新方法"", 《地球物理学报》 * |
| 周淼 等: ""综合时间序列与高程的天顶对流层延迟模型研究"", 《城市勘测》 * |
| 李博峰 等: ""不同全球对流层天顶延迟产品在中国区域的比较"", 《同济大学学报（自 然科学版）》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112711022A (en) * | 2020-12-18 | 2021-04-27 | 中国矿业大学 | GNSS chromatography-assisted InSAR (interferometric synthetic aperture radar) atmospheric delay correction method |
| CN112711022B (en) * | 2020-12-18 | 2022-08-30 | 中国矿业大学 | GNSS chromatography-assisted InSAR (interferometric synthetic aperture radar) atmospheric delay correction method |
| CN112816983A (en) * | 2021-01-06 | 2021-05-18 | 中南大学 | Time sequence InSAR turbulence atmospheric delay correction method based on optimized interferogram set |
| CN112816983B (en) * | 2021-01-06 | 2023-09-19 | 中南大学 | Time sequence InSAR turbulence atmosphere delay correction method based on optimized interference atlas |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110031841B (en) | 2021-07-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ohlmann et al. | Interpretation of coastal HF radar–derived surface currents with high-resolution drifter data | |
| Breivik et al. | Real time assimilation of HF radar currents into a coastal ocean model | |
| CN107037439B (en) | Atmospheric delay ranging error correction method for laser altimeter aiming at land target | |
| Miniscloux et al. | Geostatistical analysis of orographic rainbands | |
| CN109100717A (en) | A kind of multi-source microwave remote sensing Ocean Wind-field data fusion method and its device | |
| Mateus et al. | Three-dimensional variational assimilation of InSAR PWV using the WRFDA model | |
| CN105842692B (en) | A kind of atmospheric correction method in INSAR measurements | |
| CN105069295B (en) | Satellite and surface precipitation measured value assimilation method based on Kalman filtering | |
| CN107609713B (en) | Land-air coupling real-time forecasting method through rainfall and runoff double-factor correction | |
| CN110031841A (en) | The method and system of InSAR atmospheric delay correction based on ECMWF | |
| Sik Kim et al. | Flood simulation using the gauge‐adjusted radar rainfall and physics‐based distributed hydrologic model | |
| Wang et al. | Performance of ERA5 data in retrieving precipitable water vapor over Hong Kong | |
| Ralph et al. | A study of mountain lee waves using clear‐air radar | |
| CN113009531A (en) | Small-scale high-precision low-altitude troposphere atmosphere refractive index model | |
| Pany et al. | Atmospheric GPS slant path delays and ray tracing through numerical weather models, a comparison | |
| CN112666634A (en) | Precipitation prediction method and system | |
| Mitra et al. | Daily rainfall analysis for Indian summer monsoon region | |
| Mestas-Nuñez et al. | An evaluation of ECMWF-based climatological wind stress fields | |
| Bakış et al. | Analysis and comparison of spatial rainfall distribution applying different interpolation methods in Porsuk river basin, Turkey | |
| Manjula et al. | Diurnal variation of ducts observed over a tropical station, Gadanki, using high-resolution GPS radiosonde observations | |
| Weissman et al. | Calibrating the QuikSCAT/SeaWinds radar for measuring rainrate over the oceans | |
| Jiang et al. | Assessment of reprocessed sea surface height measurements derived from HY-2A radar altimeter and its application to the observation of 2015–2016 El Niño | |
| CN115755103A (en) | Robust self-adaptive GNSS (Global navigation satellite System) water vapor chromatography method | |
| Gregow et al. | Improving the precipitation accumulation analysis using lightning measurements and different integration periods | |
| CN114252875B (en) | High-precision meshing method for imaging altitude data |
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 |
Granted publication date: 20210723 |