CN106525651B - The method for covering day observation inverting near space atmospheric density based on X-ray - Google Patents

The method for covering day observation inverting near space atmospheric density based on X-ray Download PDF

Info

Publication number
CN106525651B
CN106525651B CN201610948414.6A CN201610948414A CN106525651B CN 106525651 B CN106525651 B CN 106525651B CN 201610948414 A CN201610948414 A CN 201610948414A CN 106525651 B CN106525651 B CN 106525651B
Authority
CN
China
Prior art keywords
ray
point
atmospheric density
atmospheric
energy
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
Application number
CN201610948414.6A
Other languages
Chinese (zh)
Other versions
CN106525651A (en
Inventor
李保权
李海涛
桑鹏
高露露
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National Space Science Center of CAS
Original Assignee
National Space Science Center of CAS
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 National Space Science Center of CAS filed Critical National Space Science Center of CAS
Priority to CN201610948414.6A priority Critical patent/CN106525651B/en
Publication of CN106525651A publication Critical patent/CN106525651A/en
Application granted granted Critical
Publication of CN106525651B publication Critical patent/CN106525651B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N9/00Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity
    • G01N9/24Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity by observing the transmission of wave or particle radiation through the material

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

The present invention relates to a kind of methods for covering day observation inverting near space atmospheric density based on X-ray, it include: to pass through the solar X-ray intensity of X-ray while measurement without decaying and cover the X-ray intensity after decaying on day path, and then the atmospheric density being fitted in 50-200km kilometers of altitude ranges.Method of the invention has the advantages such as detection altitude range is big, technology complexity is low, easy to implement, measurement accuracy is high, and this method energy real-time monitoring near space atmospheric density changes the response to solar activity, the foundation for facilitating near space atmospheric density three-dimensional chromatography dynamic model meets the needs of near space vehicle exploitation, verifying, control etc. are to atmospheric parameter.

Description

The method for covering day observation inverting near space atmospheric density based on X-ray
Technical field
It is the present invention relates to space exploration field, in particular to a kind of that day observation inverting near space atmosphere is covered based on X-ray The method of density.
Background technique
The atmospheric density and its variation characteristic for detecting near space height (50-200km) set near space vehicle Meter, control etc. are of great significance.Also lack the measurement side to 50-200km height airspace global atmosphere density in the world at present Method or Detection Techniques.
Summary of the invention
It is an object of the invention to overcome the prior art to lack the measurement to 50-200km height airspace global atmosphere density The defect of method, to provide a kind of atmospheric density measurement method easy to implement, measurement accuracy is high.
To achieve the goals above, the present invention provides one kind based on X-ray cover day observation inverting near space it is airtight greatly The method of degree, comprising:
The X-ray for measuring the solar X-ray intensity without decaying simultaneously by X-ray and covering after decaying on day path is strong Degree, and then the atmospheric density being fitted in 50-200km kilometers of altitude ranges.
In above-mentioned technical proposal, this method is specifically included:
Step 1) acquires a series of X-ray spectral lines covered on day path Jing Guo atmospheric attenuation using spaceborne X-ray energy spectrometer Intensity data, X-ray spectral line intensity data and background data without atmospheric attenuation;
Step 2, the X-ray spectral line intensity data obtained to step 1) pre-process, and the pretreatment changes including background Just, exceptional data point, the X-ray spectral line intensity data after being pre-processed are removed;
Step 3) calculates moonscope moment atmospheric density corresponding with point of contact height using Atmospheric models;
Step 4), the atmospheric density obtained using step 3) obtain step 2 using Beer law pre- as initial value Unbated X-ray intensity data after processing carry out numerical reduction to get model of fit is arrived;
Step 5), by the decaying X-ray spectrum after the model of fit that step 4) obtains and the pretreatment that step 2 obtains Line intensity data carries out nonlinear least square fitting, the atmospheric density parameter being fitted;
Step 6) repeats step 4) using the atmospheric density parameter that step 5) obtains, and then repeats step 5), by anti- Multiple iteration, until reaching stopping criterion for iteration, to obtain the big airtight of the optimal fitting in 50-200km kilometers of altitude ranges Degree.
In above-mentioned technical proposal, in step 1), the energy range of X-ray energy spectrometer measurement is 0.1-100keV or should A certain energy section in range.
In above-mentioned technical proposal, the step 3) includes: that the atmospheric density at moonscope moment is calculated using Atmospheric models, The atmospheric density value being calculated is interpolated into atmospheric density corresponding with point of contact height;Wherein, point of contact height referred to the earth's core A series of height above sea level of intersection points of the straight line and sight vertical with sight, these height above sea level angle value correspond to sampling number.
In above-mentioned technical proposal, the step 4) includes:
The X-ray radiation intensity S that spaceborne X-ray detector receivesmIt is described by Beer law, obtains being fitted mould as follows Type:
And
Wherein, R=RjlFor detector response matrix, subscript l represents high-resolution energy lattice point number;For spaceborne X-ray energy The collected the intensity of spectral line without atmospheric attenuation of spectrometer,For ambient noise;Opticalthicknessτ is defined as two-dimensional matrix τ (El, Zk), wherein ElIndicate high-resolution energy spectral coverage lattice point and ZkIndicate point of contact height;SIndicate the position coordinates of the sun;S0It represents The position coordinates of satellite;ng(Z (s)) is the number density of a certain height component gas g of direction of visual lines, which passes through atmosphere mould Type is calculated;σgFor the absorption cross-section of component gas;Normalization factor α and complex scalar factor-betagFor freely joining in model Number;
Tangential number density ng(Z (s)) and radial number density n0(Z) relationship is provided by following formula integral:
Wherein, n0(Z) the point of contact height atmospheric density obtained by atmospheric model, i.e., radial number density;Z represents point of contact A series of values of height;Z(s the point of contact height of observation moment direction of visual lines) is represented.
In above-mentioned technical proposal, in the step 5), maximum-likelihood estimator is used as using C statistics:
Wherein, diFor the X-ray spectral line intensity data point Jing Guo atmospheric attenuation, miFor i-th point of model theory value, it may be assumed that
mi=Sm(Ej,Zk) (4)
SmIt is ENERGY EjWith point of contact height ZkFunction, subscript j, k respectively corresponds detection channels and point of contact height, by covering day Sight determines;Subscript i is defined as:
I ≡ j+Nk (5)
Wherein, N is X-ray energy spectrometer energy measurement road number;
Then using Levenberg-Marquardt algorithm as the solver of nonlinear least square fitting, by making C The minimum inverting to realize point of contact height atmospheric density of statistics solves;It is high to all energy channels and point of contact simultaneously in solution procedure The day observation data point of covering of degree is fitted, and disposably obtains the atmospheric density value of all energy channels and point of contact height.
In above-mentioned technical proposal, in step 6), the condition of the iteration ends is the atmospheric density value that back calculates It is less than a certain threshold value with the atmospheric density value residual error that latter step calculates;The setting of this threshold value is needed depending on specific observation mission It asks.
The present invention has the advantages that
It is excellent that method of the invention has that detection altitude range is big, technology complexity is low, easy to implement, measurement accuracy is high etc. Gesture, and this method energy real-time monitoring near space atmospheric density changes the response to solar activity, and it is big to facilitate near space Air tightness three-dimensional chromatographs the foundation of dynamic model, meets near space vehicle exploitation, verifying, control etc. to atmospheric parameter Demand.
Detailed description of the invention
Fig. 1 is the schematic diagram that X-ray covers day observation inverting atmospheric density moonscope geometry;
Fig. 2 is the flow chart of the method for the present invention.
Specific embodiment
Now in conjunction with attached drawing, the invention will be further described.
Signified near space refers to the airspace apart from earth sea level 50-200km in the present invention.Method of the invention is wanted Realize the measurement near space atmospheric density.
Fig. 1 is that X-ray covers the schematic diagram for day observing inverting atmospheric density moonscope geometry.Under normal circumstances, satellite Orbit altitude is higher than earth atmosphere top, and spaceborne X-ray detector receives the X-ray radiation from the sun.Satellite sight (is defended The line of star and the sun) height above sea level be referred to as point of contact height Z.When point of contact height is higher than the height on atmosphere top, satellite The solar X-ray radiation received is not by atmospheric attenuation, and X-ray radiation intensity at this time is denoted as I0;When point of contact height When degree is lower than the height on atmosphere top, satellite reception to solar X-ray radiation be the X-ray at this time by atmospheric attenuation Radiation intensity is denoted as I.
With reference to Fig. 2, the method for the invention for covering day observation inverting near space atmospheric density based on X-ray includes following step It is rapid:
Step 1), acquired using spaceborne X-ray energy spectrometer a series of X-ray spectral line intensity datas by atmospheric attenuation and Without the X-ray spectral line intensity data and background data of atmospheric attenuation;
Step 2, the X-ray spectral line intensity data obtained to step 1) pre-process, and the pretreatment changes including background Just, exceptional data point etc., the X-ray spectral line intensity data after being pre-processed are removed;
Step 3) calculates moonscope moment atmospheric density corresponding with point of contact height using Atmospheric models;Wherein, described Existing Atmospheric models, such as USSA-1976 in the prior art, NRL-MSISE-00 etc. can be used in Atmospheric models.
Step 4), the atmospheric density obtained using step 3) obtain step 2 using Beer law pre- as initial value Unbated X-ray intensity data after processing carry out numerical reduction to get model of fit is arrived;
Step 5), by the decaying X-ray spectrum after the model of fit that step 4) obtains and the pretreatment that step 2 obtains Line intensity data carries out nonlinear least square fitting, the atmospheric density parameter being fitted;
Step 6) repeats step 4) using the atmospheric density parameter that step 5) obtains, and then repeats step 5), by anti- Multiple iteration, until reaching stopping criterion for iteration, to obtain the atmospheric density of optimal fitting.
Each step in the method for the present invention is described further below.
In order to meet the measurement of 50-200km altitude range, the energy range that X-ray energy spectrometer measures in step 1) is 0.1-100keV or a certain energy section within the scope of this, wherein energy channel number and sampling number it is as needed measurement essence Degree is to determine.
In step 3), the atmospheric density at moonscope moment is calculated using Atmospheric models, it is big airtight by what is be calculated Angle value is interpolated into atmospheric density corresponding with point of contact height.Wherein point of contact height referred to the earth's core straight line vertical with sight and view A series of height above sea level of intersection points of line, these height above sea level angle value correspond to sampling number.
In step 4), the X-ray radiation intensity that spaceborne X-ray detector receives is also referred to as the X of X-ray detector Ray counts, and S can be usedmIt indicates, SmIt is a channel energy E about pulsed height analyzerjWith point of contact height ZkLetter Number, subscript j represents the energy channel of pulsed height analyzer, k represents point of contact height, is determined by covering a day sight;Define Sm(Ej, Zk) it is two-dimensional matrix, then x-ray count can be described as follows by Beer law:
And
Wherein, R=RjlFor detector response matrix, subscript l represents high-resolution energy lattice point number.For spaceborne X-ray energy The collected the intensity of spectral line without atmospheric attenuation of spectrometer,For ambient noise.Opticalthicknessτ is defined as two-dimensional matrix τ (El, Zk), wherein ElIndicate high-resolution energy spectral coverage lattice point, ZkIndicate point of contact height.SIndicate the position coordinates of infinity, the present invention The position coordinates of literary middle finger x-ray radiation source, the i.e. position coordinates of the sun;Since the solar radiation that satellite reception arrives is parallel Light, it is believed that the sun is located at infinite point;Therefore the position coordinates of the sun are indicated with subfix " ∞ ".S0Represent the position of satellite Set coordinate.ng(Z (s)) is the number density of a certain height component gas g of direction of visual lines, which is calculated by Atmospheric models It arrives.σgFor the absorption cross-section of component gas.Free parameter in model is normalization factor α and complex scalar factor-betag
Tangential number density ng(Z (s)) and radial number density n0(Z) relationship is provided by following formula integral:
Wherein, n0(Z) the point of contact height atmospheric density obtained by atmospheric model, i.e., radial number density;Z is variable, Represent a series of values of point of contact height;Z(s) it is constant, represents a fixed value of point of contact height, is i.e. observation moment sight side To point of contact height.
From mathematical form, first item is that spaceborne X-ray energy spectrometer is collected without atmospheric attenuation in formula (1) The intensity of spectral line decayed by e index, this prolongs the decaying on the path of direction of visual lines in the physical sense with solar X-ray It is identical, therefore formula (1) is the measurement model that X-ray covers day observation, that is, model of fit described in step 4), This model is fitted with collected X-ray spectral line intensity data, finally obtains the value of point of contact height atmospheric density.
In order to simplify expression, m is definediThe theoretical value for day observing measurement model is covered for i-th point of X-ray, then X-ray The expression formula for covering the measurement model of day observation can transform to:
mi=Sm(Ej,Zk) (4)
Wherein, SmIt is the x-ray count of estimation, is the channel energy E of pulsed height analyzerjWith point of contact height Zk's Function, subscript j represents the energy channel of pulsed height analyzer, k represents point of contact height, is determined by covering a day sight.Subscript i definition For
I ≡ j+Nk (5)
Wherein, N is X-ray energy spectrometer energy measurement road number.
In step 5), it is contemplated that the Poisson distribution property of collected X-ray spectral line intensity data, we are united using C It is counted as maximum-likelihood estimator:
Wherein, diFor the X-ray spectral line intensity data point Jing Guo atmospheric attenuation, miFor i-th point of model theory value, benefit Use Levenberg-Marquardt (L-M) algorithm as the solver of nonlinear least square fitting, by keeping C statistics minimum To realize that the inverting of point of contact height atmospheric density solves.Solution procedure simultaneously sees the day of covering of all energy channels and point of contact height Measured data point is fitted, and disposably obtains the atmospheric density value of all energy channels and point of contact height.This data processing Journey can maximally utilise effective signal-to-noise ratio and detection data.
In step 6), a day inversion problem is covered for X-ray proposed by the present invention, we use iteration inverse approach, that is, exist It needs to carry out successive ignition to C statistical fit in fit procedure, this avoid simplify theoretical model it is assumed that sufficiently examining Actual physical process is considered, solving result is more accurate.
In step 6), the conditions of the iteration ends is the atmospheric density value that back calculates and latter step calculate it is big Airtight angle value residual error is less than a certain threshold value.The setting of this threshold value depends on specific observation mission demand.
It should be noted last that the above examples are only used to illustrate the technical scheme of the present invention and are not limiting.Although ginseng It is described the invention in detail according to embodiment, those skilled in the art should understand that, to technical side of the invention Case is modified or replaced equivalently, and without departure from the spirit and scope of technical solution of the present invention, should all be covered in the present invention Scope of the claims in.

Claims (5)

1. a kind of method for covering day observation inverting near space atmospheric density based on X-ray, comprising:
It measures the solar X-ray intensity without decaying simultaneously by X-ray and covers the X-ray intensity after decaying on day path, And then it is fitted the atmospheric density in 50-200 kilometers of altitude ranges;
This method specifically includes:
Step 1) acquires a series of X-ray the intensity of spectral line covered on day path Jing Guo atmospheric attenuation using spaceborne X-ray energy spectrometer Data, X-ray spectral line intensity data and background data without atmospheric attenuation;
Step 2), to step 1) obtain X-ray spectral line intensity data pre-process, it is described pretreatment include background correction, Remove exceptional data point, the X-ray spectral line intensity data after being pre-processed;
Step 3) calculates height corresponding atmospheric density in moonscope moment point of contact using Atmospheric models;
The step 3) includes: that the atmospheric density at moonscope moment is calculated using Atmospheric models, big airtight by what is be calculated Angle value is interpolated into atmospheric density corresponding with point of contact height;Wherein, point of contact height referred to the earth's core straight line vertical with sight and A series of height above sea level of intersection points of sight, these height above sea level angle value correspond to sampling number;
Step 4), the atmospheric density obtained using step 3) are as initial value, the pretreatment obtained using Beer law to step 2) Unbated X-ray intensity data later carry out numerical reduction to get model of fit is arrived;
It is step 5), the X-ray spectral line of decaying after the model of fit that step 4) obtains and the pretreatment that step 2) obtains is strong Degree is according to progress nonlinear least square fitting, the atmospheric density parameter being fitted;
Step 6) repeats step 4) using the atmospheric density parameter that step 5) obtains, and then repeats step 5), by changing repeatedly Generation, until reaching stopping criterion for iteration, to obtain the atmospheric density of the optimal fitting in 50-200km kilometers of altitude ranges.
2. the method according to claim 1 for covering day observation inverting near space atmospheric density based on X-ray, feature exist In in step 1), the energy range of X-ray energy spectrometer measurement is 0.1-100keV.
3. the method according to claim 1 for covering day observation inverting near space atmospheric density based on X-ray, feature exist In the step 4) includes:
The X-ray radiation intensity S that spaceborne X-ray detector receivesmIt is described by Beer law, obtains following model of fit:
And
Wherein, R=RjlFor detector response matrix, subscript l represents high-resolution energy lattice point number;For spaceborne X-ray energy spectrum The collected the intensity of spectral line without atmospheric attenuation of instrument,For ambient noise;Opticalthicknessτ is defined as two-dimensional matrix τ (El, Zk), wherein ElIndicate high-resolution energy spectral coverage lattice point and ZkIndicate point of contact height;SIndicate the position coordinates of the sun;S0It represents The position coordinates of satellite;ng(Z (s)) is the number density of a certain height component gas g of direction of visual lines, which passes through atmosphere mould Type is calculated;σgFor the absorption cross-section of component gas;Normalization factor α and complex scalar factor-betagFor freely joining in model Number;
Tangential number density ng(Z (s)) and radial number density n0(Z) relationship is provided by following formula integral:
Wherein, n0(Z) the point of contact height atmospheric density obtained by atmospheric model, i.e., radial number density;Z represents point of contact height A series of values;Z (s) represents the point of contact height of observation moment direction of visual lines.
4. the method according to claim 3 for covering day observation inverting near space atmospheric density based on X-ray, feature exist In in the step 5), using C statistics as maximum-likelihood estimator:
Wherein, diFor the X-ray spectral line intensity data point Jing Guo atmospheric attenuation, miFor i-th point of model theory value, it may be assumed that
mi=Sm(Ej,Zk) (4)
SmIt is ENERGY EjWith point of contact height ZkFunction, subscript j, k respectively corresponds detection channels and point of contact height, by covering a day sight It determines;Subscript i is defined as:
i≡j+Nk (5)
Wherein, N is X-ray energy spectrometer energy measurement road number;
Then using Levenberg-Marquardt algorithm as the solver of nonlinear least square fitting, by counting C Minimum come realize point of contact height atmospheric density inverting solve;Simultaneously to all energy channels and point of contact height in solution procedure It covers day observation data point to be fitted, disposably obtains the atmospheric density value of all energy channels and point of contact height.
5. the method according to claim 1 for covering day observation inverting near space atmospheric density based on X-ray, feature exist In, in step 6), the conditions of the iteration ends is the atmospheric density value that back calculates and latter step calculate it is big airtight Angle value residual error is less than a certain threshold value;The setting of this threshold value depends on specific observation mission demand.
CN201610948414.6A 2016-10-24 2016-10-24 The method for covering day observation inverting near space atmospheric density based on X-ray Active CN106525651B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610948414.6A CN106525651B (en) 2016-10-24 2016-10-24 The method for covering day observation inverting near space atmospheric density based on X-ray

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610948414.6A CN106525651B (en) 2016-10-24 2016-10-24 The method for covering day observation inverting near space atmospheric density based on X-ray

Publications (2)

Publication Number Publication Date
CN106525651A CN106525651A (en) 2017-03-22
CN106525651B true CN106525651B (en) 2019-04-02

Family

ID=58293406

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610948414.6A Active CN106525651B (en) 2016-10-24 2016-10-24 The method for covering day observation inverting near space atmospheric density based on X-ray

Country Status (1)

Country Link
CN (1) CN106525651B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107132155A (en) * 2017-05-02 2017-09-05 中国西安卫星测控中心 Atmospheric models modification method in short-term based on actual measurement atmospheric density
CN109827870B (en) * 2019-03-11 2021-10-19 中国科学院国家空间科学中心 Ground atmospheric density measurement method based on X-ray absorption
CN109827871B (en) * 2019-03-11 2021-08-31 中国科学院国家空间科学中心 Atmospheric density measurement system based on X-ray absorption
CN113203764A (en) * 2021-05-07 2021-08-03 北京科技大学 Material tissue quantitative analysis system using scanning electron microscope and energy spectrometer
CN113657014A (en) * 2021-08-10 2021-11-16 南京师范大学 Near space atmospheric state simulation method and device based on quantum migration

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5983705A (en) * 1996-09-20 1999-11-16 Fujitsu Limited Method for measuring atmospheric gas and system for measuring atmospheric gas
CN101866021A (en) * 2010-05-20 2010-10-20 武汉大学 Data treatment method of parallelization Abel transformation atmospheric parameters
CN103684628A (en) * 2013-11-18 2014-03-26 航天东方红卫星有限公司 Method of simultaneously acquiring target track and sub-satellite point track
CN104391299A (en) * 2014-11-19 2015-03-04 北京卫星环境工程研究所 Middle and high layer atmosphere density separation type detector
CN105372720A (en) * 2015-10-27 2016-03-02 合肥工业大学 SPSA meteorological parameter solution method based on assimilation of navigation satellite occultation data
CN105891856A (en) * 2016-04-07 2016-08-24 中国科学院国家空间科学中心 GNSS radio occultation detector

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5983705A (en) * 1996-09-20 1999-11-16 Fujitsu Limited Method for measuring atmospheric gas and system for measuring atmospheric gas
CN101866021A (en) * 2010-05-20 2010-10-20 武汉大学 Data treatment method of parallelization Abel transformation atmospheric parameters
CN103684628A (en) * 2013-11-18 2014-03-26 航天东方红卫星有限公司 Method of simultaneously acquiring target track and sub-satellite point track
CN104391299A (en) * 2014-11-19 2015-03-04 北京卫星环境工程研究所 Middle and high layer atmosphere density separation type detector
CN105372720A (en) * 2015-10-27 2016-03-02 合肥工业大学 SPSA meteorological parameter solution method based on assimilation of navigation satellite occultation data
CN105891856A (en) * 2016-04-07 2016-08-24 中国科学院国家空间科学中心 GNSS radio occultation detector

Also Published As

Publication number Publication date
CN106525651A (en) 2017-03-22

Similar Documents

Publication Publication Date Title
CN106525651B (en) The method for covering day observation inverting near space atmospheric density based on X-ray
Antoine et al. Assessment of uncertainty in the ocean reflectance determined by three satellite ocean color sensors (MERIS, SeaWiFS and MODIS‐A) at an offshore site in the Mediterranean Sea (BOUSSOLE project)
EP3262440B1 (en) Iterative ray-tracing for autoscaling of oblique ionograms
CN105067534A (en) Pollutant transport flux measurement and calculation method based on ground-based MAX-DOAS
CN113219496A (en) Spaceborne laser radar detection system for ocean atmospheric parameters
CN103792538A (en) Atmosphere profile inversion method based on foundation hyperspectral microwave radiometer
CN110909449B (en) Multi-source data ionization layer region reporting method
CN106405612A (en) Sea water body [Gamma] radiation in-site detector and detection data processing method
CN103743679B (en) A kind of method accelerating greenhouse gases vertical column retrieving concentration speed
Song et al. Tomographic reconstruction of atmospheric gravity wave parameters from airglow observations
CN102707336A (en) Novel method of using A-Train series satellite data for synergetic inversion of cloud phase states and cloud parameters
CN103197303B (en) Earth surface two-direction reflection characteristic retrieval method and earth surface two-direction reflection characteristic retrieval system based on multiple sensors
CN103728612A (en) Passive distance measuring method based on target infrared radiation spectrum and band model
Foelsche et al. Errors in GNSS radio occultation data: Relevance of the measurement geometry and obliquity of profiles
CN113253233A (en) Analysis processing method and system based on all-sky meteor radar signals
CN103809171B (en) High speed passive ranging method based on oxygen absorption and multiple regression
US9335258B2 (en) System and method of retrieving mass density distributions and thermal profiles from the atmosphere to identify molecular constituents that may absorb spectral energy
Wang et al. Correcting negatively biased refractivity below ducts in GNSS radio occultation: an optimal estimation approach towards improving planetary boundary layer (PBL) characterization
Liu et al. Retrieval of atmospheric thermodynamic state from synergistic use of radio occultation and hyperspectral infrared radiances observations
CN115480032A (en) Point source discharge intensity prediction method based on ground remote sensing measurement
CN107908838A (en) The passage choosing method and device of EO-1 hyperion satellite data inverting nitrous oxide profile
CN109827870A (en) A kind of surface air density measuring method based on X-ray absorption
CN103808413A (en) Method and device for temperature-emissivity separation based on noise separation
Frank et al. A Bayesian model to correct underestimated 3-D wind speeds from sonic anemometers increases turbulent components of the surface energy balance
CN106296630A (en) A kind of website LAI observes spatial representative evaluation methodology in Remote Sensing Products grid cell size

Legal Events

Date Code Title Description
C06 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