CN102998663A - Calibration method for microwave radiometer - Google Patents

Calibration method for microwave radiometer Download PDF

Info

Publication number
CN102998663A
CN102998663A CN2012105830404A CN201210583040A CN102998663A CN 102998663 A CN102998663 A CN 102998663A CN 2012105830404 A CN2012105830404 A CN 2012105830404A CN 201210583040 A CN201210583040 A CN 201210583040A CN 102998663 A CN102998663 A CN 102998663A
Authority
CN
China
Prior art keywords
sky
temperature
value
ref
microwave radiometer
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
Application number
CN2012105830404A
Other languages
Chinese (zh)
Other versions
CN102998663B (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.)
QINGDAO BRANCH, NO.22 INSTITUTE, MINISTRY OF INFORMATION INDUSTRY
Original Assignee
China Research Institute of Radio Wave Propagation CRIRP
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 China Research Institute of Radio Wave Propagation CRIRP filed Critical China Research Institute of Radio Wave Propagation CRIRP
Priority to CN201210583040.4A priority Critical patent/CN102998663B/en
Publication of CN102998663A publication Critical patent/CN102998663A/en
Application granted granted Critical
Publication of CN102998663B publication Critical patent/CN102998663B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Radiation Pyrometers (AREA)

Abstract

The invention discloses a calibration method for a microwave radiometer. The calibration method includes calculating sky brightness temperature that a plurality of zenith angles correspond to according to the brightness temperature transfer function, calculating atmosphere total pad value that the plurality of zenith angles correspond to, conducting linear regression treatment on the secant value of the plurality of zenith angles and the atmosphere total pad value that the plurality fo zenith angles correspond to, utilizing slope of a straight line obtained by the linear regression treatment as the atmosphere total pad value at the time when the zenith angle is 0 degree, calculating the sky brightness temperature when the zenith angle is 0 degree, substitute the sky brightness temperature at the time when the zenith angle is 0 degree into the brightness transfer function to calculate the value of the corresponding calibration coefficient a, updating the value of the calibration coefficient a, repeating the above step to obtain a new a and circulating in this way to enable the restrained a as the calibration value. The calibration method only requires one microwave radiometer and a black body, is low in cost, convenient to operate and high in accuracy, and provides forceful support for detection of atmosphere temperature profile, atmosphere humidity profile and the like.

Description

A kind of calibrating method for microwave radiometer
Technical field
The present invention relates to the microwave radiometer technical field, particularly relate to a kind of calibrating method for microwave radiometer.
Background technology
Microwave radiometer is a kind of passive type ground microwave remote sensing equipment, the calibrating method of microwave radiometer is divided into two kinds substantially: a kind of is the substep scaling method, carry out respectively for the calibration of receiver and the calibration of antenna, need the intermediate parameters of measurement more, therefore the source of error of substep scaling method is many, and calibration precision is low; Another kind is whole scaling method, once finish for the calibration that is scaled to final output from antenna, without intermediate link, thus calibration precision the substep scaling method want high.In whole scaling method, usually need to use two noise sources stable, different temperatures, in the reality often with black matrix as the high temperature noise source, with liquid nitrogen as the low temperature noise source, wherein, although use liquid nitrogen to finish calibration accurately, shortcoming is the liquid nitrogen price, and transport and place all inconveniences, it is higher to cause calibrating cost.
Summary of the invention
In order to solve the high problem of calibration cost of microwave radiometer in the prior art, the invention provides a kind of calibrating method for microwave radiometer.
Calibrating method for microwave radiometer of the present invention comprises: step 1, according to bright temperature transport function, bright temperature and corresponding voltage thereof based on reference blackbody, and the initial value of calibration coefficient a, calculate sky temperature corresponding to a plurality of zenith angles, the initial value of wherein said calibration coefficient a is predefined; Step 2 based on the sky temperature under described a plurality of zenith angles, is calculated each self-corresponding atmosphere overall attenuation value of described a plurality of zenith angle; Step 3 for secant value and each self-corresponding atmosphere overall attenuation value of described a plurality of zenith angle of described a plurality of zenith angles, is carried out linear regression processing, and take the slope of linear regression processing gained straight line as zenith angle the atmosphere overall attenuation value during as 0 °; Step 4, the atmosphere overall attenuation value when being 0 ° based on described zenith angle, the sky temperature when the calculating zenith angle is 0 °; Step 5, the described bright temperature transport function of the sky temperature substitution when described zenith angle is 0 ° solves the value of corresponding calibration coefficient a; Step 6, as the described initial value in the step 1, repeating step one to five obtains the value of new a with the value of the calibration coefficient a of described correspondence; So circulation, constantly upgrading the value of a, with the value of a after the convergence scaled values as a, calibration is finished.
Further, described bright temperature transport function is T Sky=T Ref+ G (V Sky-V Ref) f w, wherein, T RefThe bright temperature of reference blackbody, V RefThe voltage of corresponding catoptron when pointing to reference blackbody; T SksSky temperature, V SkyThe voltage of corresponding catoptron when pointing into the sky, f wBe the constant relevant with the emissivity of the material that covers the microwave radiometer minute surface, G is the calibration gain,
Figure BDA00002663565600021
Wherein, G = T nd V ref + nd - V ref , Then a = T ref - T nd · f w · V ref V ref + nd - V ref ; Wherein, V Ref+ndThe voltage of microwave radiometer catoptron when pointing to reference blackbody and adding noise diode, T NdBe noise temperature, noise temperature is with the environmental change change; Calibration is equivalent to T to a NdCalibration.
Further, described bright temperature transport function is T Sky=f (a, V, σ), wherein, V represents parameter relevant with voltage in the microwave radiometer, and σ represents other parameter in the microwave radiometer.
Further, in the step 2, according to sky temperature, the bright temperature of average radiation and atmosphere overall attenuation three's funtcional relationship, calculate each self-corresponding atmosphere overall attenuation value of described a plurality of zenith angle.
Further, in the step 4, the atmosphere overall attenuation value when being 0 ° according to described zenith angle, the bright gentle sky temperature three's of the average radiation when zenith angle is 0 ° funtcional relationship, the sky temperature when calculating described zenith angle and being 0 °.
Wherein, any two in described a plurality of zenith angles differ 15 ° at least, and arbitrary zenith angle is no more than 70 °.
Beneficial effect of the present invention is as follows:
The present invention with black matrix as high temperature source, according to bright temperature transport function, the bright temperature of zenith when utilizing the sky temperature of microwave radiometer under different zenith angles to calculate zenith angle to be 0 °, substituted the effect as the liquid nitrogen of cold temperature source, saved input cost, and then obtain the value of calibration coefficient by cycle calculations, be applied to the temperature coefficient timing signal to the microwave radiometer noise diode, the temperature that can be good at calibrating noise diode.The present invention only needs a microwave radiometer and a black matrix, and cost is low, and is easy and simple to handle, and precision is high, but and all weather operations.Bright temperature transport function wherein is out of shape, this calibrating method can be extended to linear between the output voltage of receiver and the input noise temperature or be in the Ground-Based Microwave Radiometer of nonlinear relationship, for the profiles such as atmospheric envelope temperature, humidity in real time, efficient, the low-cost detection provide strong support.The present invention can be applicable to the inverting of atmospheric envelope temperature profile, target detection and affirmation, supervision and drawing etc.
Description of drawings
Fig. 1 is the calibrating method process flow diagram that the present invention is used for microwave radiometer.
Fig. 2 is the result schematic diagram that circulates for the first time in the embodiment of the invention calibration process.
Fig. 3 is the result schematic diagram of convergence after Fig. 2 embodiment repeatedly circulates.
Embodiment
Below in conjunction with accompanying drawing and embodiment, the present invention is further elaborated.Should be appreciated that specific embodiment described herein only in order to explain the present invention, does not limit the present invention.
Fig. 1 is the calibrating method process flow diagram that the present invention is used for microwave radiometer, and the method may further comprise the steps:
Step 1, according to bright temperature transport function, based on bright temperature and the corresponding voltage thereof of reference blackbody, and the initial value of calibration coefficient a, calculating sky temperature corresponding to a plurality of zenith angles, the initial value of wherein said calibration coefficient a is predefined;
Step 2 based on the sky temperature under described a plurality of zenith angles, is calculated each self-corresponding atmosphere overall attenuation value of described a plurality of zenith angle;
Step 3 for secant value and each self-corresponding atmosphere overall attenuation value of described a plurality of zenith angle of described a plurality of zenith angles, is carried out linear regression processing, and take the slope of linear regression processing gained straight line as zenith angle the atmosphere overall attenuation value during as 0 °;
Step 4, the atmosphere overall attenuation value when being 0 ° based on described zenith angle, the sky temperature when the calculating zenith angle is 0 °;
Step 5, the described bright temperature transport function of the sky temperature substitution when described zenith angle is 0 ° solves the value of corresponding calibration coefficient a;
Step 6, as the described initial value in the step 1, repeating step one to five obtains the value of new a with the value of the calibration coefficient a of described correspondence; So circulation, constantly upgrading the value of a, with the value of a after the convergence scaled values as a, calibration is finished.
The method with black matrix as high temperature source, according to bright temperature transport function, the bright temperature of zenith when utilizing the sky temperature of microwave radiometer under different zenith angles to calculate zenith angle to be 0 °, and then calculate corresponding calibration coefficient obtains the value of final calibration coefficient through cycle calculations repeatedly.
Specifically, 1. above-mentioned calibrating method is applied in the Ground-Based Microwave Radiometer linear between the receiver output voltage and input noise temperature, slightly different from the situation in the Ground-Based Microwave Radiometer that is applied to 2. to be nonlinear relationship, below both of these case is described respectively.
1. for Ground-Based Microwave Radiometer linear between receiver output voltage and the input noise temperature, described bright temperature transport function is:
T sky=T ref+G(V sky-V ref)f w
=T ref-Gf wV ref+Gf wV sky
=a+bV sky
In the formula, T RefThe bright temperature of reference blackbody, V RefThe catoptron of the microwave radiometer voltage when pointing to reference blackbody; T SkySky temperature, V SkyThe voltage of catoptron when pointing into the sky, f wBe the constant relevant with the emissivity of the material that covers the microwave radiometer minute surface, G is the calibration gain; Further, a=T Ref-Gf wV Ref, b=Gf w, also be
Figure BDA00002663565600041
Under the prerequisite of voltage corresponding to the bright gentleness of known reference black matrix, the relation of a and b: a=T Ref-bV Ref, b=(T Ref-a)/V Ref, a or b are regarded calibration coefficient, a, b can change mutually, and present embodiment is take a as calibration coefficient.
According to bright temperature transport function, the initial value of a given a, then sky temperature value corresponding to different zenith angle θ:
T sky ( θ ) = a + bV sky ( θ ) = a + ( T ref - a ) V ref V sky ( θ )
During zenith angle θ=0 °, the corresponding bright temperature T of zenith Sky(0 °)=a+bV Sky(0 °), T Ref=a+bV Ref, can be by the bright temperature T of zenith Sky(0 °) solves a new a value:
Figure BDA00002663565600051
This a value back substitution is circulated in bright temperature transport function next time, constantly upgrade a until its convergence, as calibration coefficient accurately, calibration is finished the value of a that obtains.
Again further, the below describes technical scheme of the present invention so that the temperature of microwave radiometer noise diode is calibrated as example.
Under the prerequisite of the bright gentle corresponding voltage of known reference black matrix, ARM(Atmosphere RadiationMeasure, atmosphere radiation is measured) bright temperature computing formula that microwave radiometer provides, namely bright temperature transport function is:
T sky=T ref+G(V sky-V ref)f w (1)
In the formula
G=T nd/(V ref+nd-V ref) (2)
In the formula, V Ref+ndThe voltage of radiometer catoptron when pointing to reference blackbody and adding noise diode, T NdThe temperature of noise diode, although T NdRelatively stable, but its meeting is along with environmental evolution produces very little change, so T NdTo need definite calibration coefficient.Here, relate to the calibration coefficient a of an embodiment, have
Figure BDA00002663565600052
To a calibrate namely be equivalent to present embodiment to T NdCalibrate.
So, adopt calibrating method of the present invention, in known first noise source---the bright temperature T of reference blackbody Ref---prerequisite under, need calculate again second noise source---the bright temperature T of zenith Sky(0 °) just can solve accurately calibration coefficient T Nd, concrete computation process is as follows.
Give the calibration coefficient T NdAn initial value is set, and substitution (1) formula is obtained the bright temperature under the different zenith angles (can choose several bright temperature and differ larger zenith angle):
T sky ( θ ) = T ref + T nd ( V ref + nd - V ref ) ( V sky ( θ ) - V ref ) f w - - - ( 3 )
With the bright temperature of bright gentle average radiation under each zenith angle that obtains (the bright temperature of average radiation can obtain by itself and the temperature on ground and the statistical relationship of relative humidity) substitution following formula, can obtain the atmosphere overall attenuation under each zenith angle:
τ ( 0 , ∞ ) ( θ ) = ln [ T m ( θ ) - 2.73 T m ( θ ) - T sky ( θ ) ] - - - ( 4 )
Herein, formula (4) is zenith angle θ and atmosphere overall attenuation, the bright temperature three's of the bright gentle average radiation of zenith functional relation, the bright temperature T of average radiation mCan be by the temperature on ground and the statistical relationship of relative humidity be obtained.
Below briefly introduce the derivation of formula (4): this area, the fundamental formular of atmosphere overall attenuation and sky temperature is:
T sky 2.73 e - τ ( 0 , ∞ ) + ∫ 0 ∞ k e ( r ) T ( r ) e - τ ( 0 , r ) dr - - - ( 5 )
In the formula, r is the distance apart from ground,
Figure BDA00002663565600063
Be the atmosphere overall attenuation (or claiming atmosphere optical thickness, the Np of unit) between 0~r, k eThe body extinction coefficient, the Npm of unit -1
The bright temperature T of the average radiation of atmosphere mBe defined as
Figure BDA00002663565600064
Wherein denominator is variable is changed to:
∫ 0 ∞ k e ( r ) e - τ ( 0 , r ) dr
= ∫ 0 τ ( 0 , ∞ ) e - τ ( 0 , r ) dτ ( 0 , r )
= 1 - e - τ ( 0 , ∞ )
Then the bright temperature of zenith can be expressed as the function of atmosphere overall attenuation and average radiation brightness, for:
T sky = 2.73 e - τ ( 0 , ∞ ) + T m ∫ 0 ∞ k e ( r ) e - τ ( 0 , r ) dr
= 2.73 e - τ ( 0 , ∞ ) + T m [ 1 - e - τ ( 0 , ∞ ) ]
So the atmosphere overall attenuation can be expressed as the function of the bright temperature of the bright gentle average radiation of zenith, i.e. formula (4):
Figure BDA000026635656000610
Bright temperature when being 0 ° in order to obtain zenith angle, carry out following calculating:
According to the relation of zenith angle and atmosphere overall attenuation, vertical height is z, and the atmosphere overall attenuation when zenith angle is 0 ° is
Figure BDA00002663565600071
And for the uniform atmosphere of level, vertical height is z, and the wrong path electrical path length is r, and the atmosphere overall attenuation when zenith angle is θ is:
Figure BDA00002663565600072
Figure BDA00002663565600073
Figure BDA00002663565600074
Following formula r gets just infinite, obtains:
τ (0,∞)(θ)=τ (0,∞)(0°)sec(θ) (6)
Can be found out that by formula (6) the secant value sec (θ) of atmosphere overall attenuation τ (θ) and zenith angle is linear relationship, crosses initial point by the straight line that linear regression processing obtains, and the atmosphere overall attenuation τ of the slope of this straight line when equaling zenith angle θ and being 0 ° (0, ∞)(0 °) value.
In the reality, the secant value of atmosphere overall attenuation and zenith angle is carried out linear regression, if related coefficient, thinks that linear regression is effective greater than 0.998, then make the atmosphere overall attenuation τ of zenith direction (0, ∞)(0 °) equals its slope, otherwise the weather condition of explanation this moment is relatively poor or be affected by other factors and be not suitable for self calibration.
Atmosphere overall attenuation τ with zenith direction (0, ∞)In (0 °) and the bright temperature substitution of the zenith direction average radiation formula (5), calculate the bright temperature of zenith:
With the bright temperature T of the zenith that obtains Sky(0 °) substitution formula (1) solves a new T NdValue:
Figure BDA00002663565600076
Then, with this new T NdValue is to calibration coefficient T NdUpgrade, repeat said process, constantly upgrade T Nd, through repeatedly circulating until it is restrained, the T that obtains at last NdAs calibration coefficient accurately, calibration is finished.
More than, describe the present invention by embodiment and be applied to 1. situation in the Ground-Based Microwave Radiometer linear between the receiver output voltage and input noise temperature, below be the process that is applied to 2. to be in the Ground-Based Microwave Radiometer of nonlinear relationship.
2. for the Ground-Based Microwave Radiometer that is nonlinear relationship between receiver output voltage and the input noise temperature, described bright temperature transport function is:
T sky=f(a,V,σ) (8)
Wherein, V represents parameter relevant with voltage in the microwave radiometer, may be a plurality of, and σ represents other parameter in the microwave radiometer, also may be a plurality of, and then the bright temperature value of different zenith angles is:
T sky(θ)=f(a,V(θ),σ)
At calibrating method according to the present invention, when calculating new calibration coefficient by the bright temperature of zenith, with T SkyBe independent variable, a is dependent variable, and its inverse function is:
a=f -1(T sky(θ),V,σ)
The bright temperature substitution of zenith following formula can be obtained new calibration coefficient a=f -1(T Sky(0 °), V (0 °), σ).Follow-up value with new a renewal calibration coefficient repeats said process, constantly upgrades a, is the calibration coefficient value that needs after the convergence.
In addition, calibrating method of the present invention is applied in the dual-frequency microwave radiometer of certain independent development, this radiometer belongs to real-time calibration microwave radiometer, the output equation of radiometer is written as the form of T=a+bV, that is:
T = T 1 ′ - T 1 ′ - T 2 ′ V 1 - V 2 V 1 + T 1 ′ - T 2 ′ V 1 - V 2 V = a + bV
In the formula
Figure BDA00002663565600082
Figure BDA00002663565600083
T and V are respectively the bright gentle corresponding voltage that this antenna for radiometer receives, and it can utilize calibrating method of the present invention that calibration coefficient a or b are calibrated calibration.
The bright temperature data of emulation that below will survey by experiment bright temperature data and sounding contrast, and have verified the validity of this calibrating method.
Certain dual-frequency microwave radiometer is placed tall building roof field all around, the zenith angle that (positive north) is chosen at 0 ° of position angle during calibration is 0 °, 45 °, 60 °, for the compensation apparatus lateral error, two equal zenith angles have been increased at 180 ° of position angles again, 45 °, 60 °; The frequency of present embodiment is 23.8GHz.Aspect the selection of zenith angle, consider that the side lobe effect at the low elevation angle is larger on the impact of bright temperature measurement, should guarantee that any two differ 15 ° at least in the selected zenith angle, maximum zenith angle is no more than 70 °.
Fig. 2 is the result of circulation for the first time in the calibration process, and horizontal ordinate is the secant value of zenith angle, the atmosphere overall attenuation value that the circle representative is calculated by initial five bright temperature of zenith angle, and straight line is the figure line after the linear regression, figure line initial point only among Fig. 1.
Fig. 3 is the afterwards result of convergence that repeatedly circulates, and figure line is near initial point, and calibration finishes.The related coefficient of final regression curve is 0.9983, satisfies the requirement (related coefficient is greater than 0.998) of effective self calibration.
Table 1 is the comparing data of surveying the bright temperature of emulation of bright temperature and Qingdao sounding before and after this dual-frequency microwave radiometer calibration on Dec 10th, 2011.Relatively calibrate the difference of front and back bright temperature that radiometer is surveyed, can find out the validity of calibrating method of the present invention.
The comparison of the table 1 bright temperature of actual measurement and the bright temperature of emulation
Figure BDA00002663565600091
In sum, calibrating method provided by the invention only needs a microwave radiometer and a black matrix, be applicable to the calibration of general Ground-Based Microwave Radiometer, cost is low, easy and simple to handle, precision is high, but and all weather operations, can be applicable to the fields such as the inverting of atmospheric envelope temperature profile, target detection and affirmation, supervision, drawing.
Although be the example purpose, the preferred embodiments of the present invention are disclosed, it also is possible those skilled in the art will recognize various improvement, increase and replacement, therefore, scope of the present invention should be not limited to above-described embodiment.

Claims (8)

1. a calibrating method that is used for microwave radiometer is characterized in that, comprising:
Step 1, according to bright temperature transport function, based on bright temperature and the corresponding voltage thereof of reference blackbody, and the initial value of calibration coefficient a, calculating sky temperature corresponding to a plurality of zenith angles, the initial value of wherein said calibration coefficient a is predefined;
Step 2 based on the sky temperature under described a plurality of zenith angles, is calculated each self-corresponding atmosphere overall attenuation value of described a plurality of zenith angle;
Step 3 for secant value and each self-corresponding atmosphere overall attenuation value of described a plurality of zenith angle of described a plurality of zenith angles, is carried out linear regression processing, and take the slope of linear regression processing gained straight line as zenith angle the atmosphere overall attenuation value during as 0 °;
Step 4, the atmosphere overall attenuation value when being 0 ° based on described zenith angle, the sky temperature when the calculating zenith angle is 0 °;
Step 5, the described bright temperature transport function of the sky temperature substitution when described zenith angle is 0 ° solves the value of corresponding calibration coefficient a;
Step 6, as the described initial value in the step 1, repeating step one to five obtains the value of new a with the value of the calibration coefficient a of described correspondence; So circulation, constantly upgrading the value of a, with the value of a after the convergence scaled values as a, calibration is finished.
2. the calibrating method for microwave radiometer as claimed in claim 1 is characterized in that, described bright temperature transport function is T Sky=T Ref+ G (V Sky-V Ref) f w, wherein, T RefThe bright temperature of reference blackbody, V RefThe catoptron of the microwave radiometer voltage when pointing to reference blackbody; T SkySky temperature, V SkyThe voltage of catoptron when pointing into the sky, f wBe the constant relevant with the emissivity of the material that covers the microwave radiometer minute surface, G is the calibration gain, G = T ref - a f w V ref .
3. the calibrating method for microwave radiometer as claimed in claim 2 is characterized in that, wherein, G = T nd V ref + nd - V ref , Then a = T ref - T nd · f w · V ref V ref + nd - V ref ; Wherein, V Ref+ndThe voltage of microwave radiometer catoptron when pointing to reference blackbody and adding noise diode, T NdThe temperature of noise diode, the temperature T of noise diode NdWith the environmental change change, a calibrated be equivalent to T NdCalibrate.
4. the calibrating method for microwave radiometer as claimed in claim 1 is characterized in that, described bright temperature transport function is T Sky=f (a, V, σ), wherein, V represents parameter relevant with voltage in the microwave radiometer, and σ represents other parameter in the microwave radiometer.
5. the calibrating method for microwave radiometer as claimed in claim 1, it is characterized in that, in the step 2, according to the funtcional relationship of zenith angle, sky temperature, the bright temperature of average radiation and atmosphere overall attenuation, calculate each self-corresponding atmosphere overall attenuation value of described a plurality of zenith angle.
6. the calibrating method for microwave radiometer as claimed in claim 1, it is characterized in that, in the step 4, the funtcional relationship of the atmosphere overall attenuation value when being 0 ° according to sky temperature and described zenith angle, the bright temperature of average radiation when zenith angle is 0 °, the sky temperature when calculating described zenith angle and being 0 °.
7. the calibrating method for microwave radiometer as claimed in claim 1 is characterized in that, any two in described a plurality of zenith angles differ 15 ° at least.
8. the calibrating method for microwave radiometer as claimed in claim 1 is characterized in that, any is no more than 70 ° in described a plurality of zenith angles.
CN201210583040.4A 2012-12-27 2012-12-27 Calibration method for microwave radiometer Active CN102998663B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210583040.4A CN102998663B (en) 2012-12-27 2012-12-27 Calibration method for microwave radiometer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210583040.4A CN102998663B (en) 2012-12-27 2012-12-27 Calibration method for microwave radiometer

Publications (2)

Publication Number Publication Date
CN102998663A true CN102998663A (en) 2013-03-27
CN102998663B CN102998663B (en) 2014-10-29

Family

ID=47927472

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210583040.4A Active CN102998663B (en) 2012-12-27 2012-12-27 Calibration method for microwave radiometer

Country Status (1)

Country Link
CN (1) CN102998663B (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104181507A (en) * 2013-05-23 2014-12-03 中国科学院空间科学与应用研究中心 Data processing and controlling apparatus used for ground-based microwave radiometer
CN104181511A (en) * 2013-05-23 2014-12-03 中国科学院空间科学与应用研究中心 Scaling method of foundation microwave radiometer
CN104504256A (en) * 2014-12-12 2015-04-08 中国电子科技集团公司第二十二研究所 Estimation algorithm for accurately inverting boundary layer temperature profile
CN104573816A (en) * 2014-12-18 2015-04-29 中国电子科技集团公司第二十二研究所 Neural network clustering method for microwave radiometer remote-sensing atmospheric parameter
CN105372610A (en) * 2015-10-27 2016-03-02 中国科学院国家空间科学中心 Four-point scaling device and method for microwave radiometer
CN106054147A (en) * 2016-05-13 2016-10-26 国家卫星气象中心 Microwave radiometer brightness temperature non-linear correction method and system
CN108267739A (en) * 2016-12-30 2018-07-10 中国科学院国家空间科学中心 A kind of Ground-Based Microwave Radiometer and its calibrating method, Atmospheric Survey method
CN109521405A (en) * 2018-12-05 2019-03-26 国家卫星气象中心 A kind of unified calibrating method suitable for spaceborne large aperture antenna microwave radiometer
CN110261934A (en) * 2019-07-18 2019-09-20 常州感通威视微波科技有限公司 Ground-Based Microwave Radiometer system and its calibrating method
CN113484917A (en) * 2021-06-30 2021-10-08 中国电波传播研究所(中国电子科技集团公司第二十二研究所) Automatic calibration method for microwave radiometer
US11567015B2 (en) 2020-12-30 2023-01-31 Boulder Environmental Sciences and Technology Systems for passive microwave remote sensing and their calibration methods

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003050286A (en) * 2001-08-06 2003-02-21 Communication Research Laboratory Rainfall intensity estimation method and system in mountainous region
CN102243294A (en) * 2010-05-14 2011-11-16 中国科学院空间科学与应用研究中心 Nonlinear calibrating method and device for foundation microwave radiometer
CN102032949B (en) * 2009-09-28 2012-06-27 中国科学院空间科学与应用研究中心 External calibration device of geosynchronous orbit atmospheric temperature detecting instrument and external calibration method thereof

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003050286A (en) * 2001-08-06 2003-02-21 Communication Research Laboratory Rainfall intensity estimation method and system in mountainous region
CN102032949B (en) * 2009-09-28 2012-06-27 中国科学院空间科学与应用研究中心 External calibration device of geosynchronous orbit atmospheric temperature detecting instrument and external calibration method thereof
CN102243294A (en) * 2010-05-14 2011-11-16 中国科学院空间科学与应用研究中心 Nonlinear calibrating method and device for foundation microwave radiometer

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李靖等: "微波辐射计定标", 《遥感技术与应用》, vol. 14, no. 1, 31 March 1999 (1999-03-31), pages 1 - 4 *
肖志辉等: "地基、 空基、 星基微波辐射计定标技术概览", 《遥感技术与应用》, vol. 15, no. 2, 30 June 2000 (2000-06-30), pages 113 - 120 *

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104181511A (en) * 2013-05-23 2014-12-03 中国科学院空间科学与应用研究中心 Scaling method of foundation microwave radiometer
CN104181507A (en) * 2013-05-23 2014-12-03 中国科学院空间科学与应用研究中心 Data processing and controlling apparatus used for ground-based microwave radiometer
CN104504256A (en) * 2014-12-12 2015-04-08 中国电子科技集团公司第二十二研究所 Estimation algorithm for accurately inverting boundary layer temperature profile
CN104573816A (en) * 2014-12-18 2015-04-29 中国电子科技集团公司第二十二研究所 Neural network clustering method for microwave radiometer remote-sensing atmospheric parameter
CN104573816B (en) * 2014-12-18 2017-12-26 中国电子科技集团公司第二十二研究所 The neural network clustering method of microwave radiometer remote sensing atmosphere parameter
CN105372610A (en) * 2015-10-27 2016-03-02 中国科学院国家空间科学中心 Four-point scaling device and method for microwave radiometer
CN105372610B (en) * 2015-10-27 2017-12-29 中国科学院国家空间科学中心 A kind of device and method for 4 points of calibrations of microwave radiometer
CN106054147B (en) * 2016-05-13 2019-01-25 国家卫星气象中心 A kind of method and system of the bright temperature gamma correction of microwave radiometer
CN106054147A (en) * 2016-05-13 2016-10-26 国家卫星气象中心 Microwave radiometer brightness temperature non-linear correction method and system
CN108267739A (en) * 2016-12-30 2018-07-10 中国科学院国家空间科学中心 A kind of Ground-Based Microwave Radiometer and its calibrating method, Atmospheric Survey method
CN109521405A (en) * 2018-12-05 2019-03-26 国家卫星气象中心 A kind of unified calibrating method suitable for spaceborne large aperture antenna microwave radiometer
CN109521405B (en) * 2018-12-05 2020-05-08 国家卫星气象中心 Full-aperture calibration method suitable for satellite-borne large-aperture antenna microwave radiometer
CN110261934A (en) * 2019-07-18 2019-09-20 常州感通威视微波科技有限公司 Ground-Based Microwave Radiometer system and its calibrating method
CN110261934B (en) * 2019-07-18 2022-03-08 常州感通威视微波科技有限公司 Foundation microwave radiometer system and calibration method thereof
US11567015B2 (en) 2020-12-30 2023-01-31 Boulder Environmental Sciences and Technology Systems for passive microwave remote sensing and their calibration methods
CN113484917A (en) * 2021-06-30 2021-10-08 中国电波传播研究所(中国电子科技集团公司第二十二研究所) Automatic calibration method for microwave radiometer
CN113484917B (en) * 2021-06-30 2023-08-11 中国电波传播研究所(中国电子科技集团公司第二十二研究所) Automatic calibration method for microwave radiometer

Also Published As

Publication number Publication date
CN102998663B (en) 2014-10-29

Similar Documents

Publication Publication Date Title
CN102998663B (en) Calibration method for microwave radiometer
CN1882854B (en) Correction of humidity measurement results of a radiosonde
CN106569186B (en) A kind of deviation correction method of satellite-borne microwave radiometer
CN102243294B (en) Nonlinear calibrating method and device for foundation microwave radiometer
CN104181511A (en) Scaling method of foundation microwave radiometer
US20180329022A1 (en) Method, apparatus and system for locating an object using cluster-type magnetic field
Sherwood et al. Atmospheric changes through 2012 as shown by iteratively homogenized radiosonde temperature and wind data (IUKv2)
CN101915914B (en) Lookup table based pixel-by-pixel atmospheric correction method of remote sensing images
US20110057840A1 (en) Method of Positioning RFID Tags
US20210048557A1 (en) Method and device for simulating discharge, and computer device
CN106324580A (en) Method for correcting radar echo attenuation based on microwave link network
CN111596145B (en) Phase control array emission calibration device and method
CN108318871A (en) A kind of cold empty scaled correction method of satellite-borne microwave radiometer
CN104865274A (en) Soil water content inversion method and system based on single polarization radar data
CN110764153B (en) System and method for correcting on-orbit error of hot mirror back lobe of satellite-borne microwave imager
CN116297068A (en) Aerosol optical thickness inversion method and system based on earth surface reflectivity optimization
CN102830448B (en) Method, device and system for microwave hyper-spectral clear air scaling
CN104483646B (en) Real-time calibration device and method for ground-based microwave radiometer
CN108983172B (en) SAR radiation calibration method based on random evolution constraint
CN109684712B (en) BRDF (bidirectional reflectance distribution function) model creating method and device, electronic equipment and storage medium
CN103092815B (en) To the method transferring function by calibration in monitoring device
CN109241565A (en) Sampling location error calibration method and device in antenna spherical surface near field measurement
CN111337065A (en) Small-size environmental simulation calibration test device of microwave radiometer
KR20180086373A (en) Meteorological radiation cooling temperature measuring system and method of thereof
US20190346498A1 (en) Field probe isotropic compensation using orthogonal scalar field components

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20170523

Address after: Chengyang District of Shandong city in Qingdao Province mountain road 266107 No. 36

Patentee after: QINGDAO BRANCH, NO.22 INSTITUTE, MINISTRY OF INFORMATION INDUSTRY

Address before: Chengyang District of Shandong city in Qingdao Province mountain road 266107 No. 36

Patentee before: No.22 Institute, China Electronic Science &. Technology Group Corp.

TR01 Transfer of patent right