CN101976297A - Processing method of moon brightness temperature data observed by foundation single antenna - Google Patents
Processing method of moon brightness temperature data observed by foundation single antenna Download PDFInfo
- Publication number
- CN101976297A CN101976297A CN2010102973619A CN201010297361A CN101976297A CN 101976297 A CN101976297 A CN 101976297A CN 2010102973619 A CN2010102973619 A CN 2010102973619A CN 201010297361 A CN201010297361 A CN 201010297361A CN 101976297 A CN101976297 A CN 101976297A
- Authority
- CN
- China
- Prior art keywords
- dirty
- moon
- antenna
- clean
- brightness temperature
- 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
Images
Landscapes
- Radiation Pyrometers (AREA)
Abstract
The invention discloses a processing method of moon brightness temperature data observed by a foundation single antenna, comprising the following steps: observing by the foundation single antenna to obtain moon brightness temperature observation data; calibrating the observation data to obtain a dirty diagram; observing a radio source with known flow density by the foundation single antenna to obtain an antenna directional diagram; according to the antenna directional diagram, obtaining normalized dirty beam and clean beam; according to the dirty diagram, the dirty beam and the clean beam, cleaning the dirty diagram by a modified CLEAN algorithm; and obtaining a clean diagram. In the processing method of moon brightness temperature data observed by foundation single antenna of the invention, the effective deconvolution CLEAN algorithm is introduced so as to eliminate influence on the antenna directional diagram to certain degree, thus obtaining a result which more approaches to the practical brightness distribution of the moon.
Description
Technical field
The present invention relates to a kind of data processing method, relate in particular to a kind of ground single antenna observation moon brightness temperature data processing method.
Background technology
The moon is the celestial body nearest apart from the earth, the ground moon center mean distance be 384400km; Its diameter is about 3476km, is about 1/3.67 of the earth; Quality is about 7.350 * 1025g, is about 1/81.301 of the earth; The visual angle directly is 33 ' 31 to the maximum ", minimum is 29 ' 22 ", average visual angle directly is 31 ' 5 "; The surface gravity acceleration is about 1/6 of the earth.
The moonscape information of utilizing telemetry approaches to obtain all comes from electromagnetic radiation, the electromagnetic radiation field can be divided into fields such as gamma-rays, X ray, ultraviolet ray, visible light, infrared ray, radio wave according to the difference of wavelength, in these fields moon characteristic is all had research.The ground remote measurement utilizes visible light, infrared ray, radio means that moon characteristic is studied more.Galileo has write down the moonscape image that utilizes optical telescope observation to obtain first.Utilize the menology image to obtain the menology albedo in the optics frequency range, can infer the menology mineral composition, grain size etc.
No matter be which kind of the means research moon, understand lunar material and form and to study lunar soil.Lunar soil is that one deck of being covered with on moonscape is by material potpourris that form, loosely organized such as rock chip, powder, dust, bump melten glass.Lunar soil is to obtain the important source that lunar material is formed information, and the telescopical observed object of space probe and ground is also mostly to be lunar soil.In the radio frequency range moon is surveyed heat distribution, the thermal characteristics that can obtain the lunar interior, thereby obtain the physical characteristics of lunar soil.
Brightness temperature is the important indicator of reflected microwave radiation.If celestial body equates that with the absolute black body of a certain temperature T B then the temperature T B of this absolute black body is called the brightness temperature of this object in this frequency in the spectrum radiancy of a certain frequency in the radiation intensity of same frequency.Research lunar soil characteristic is absolutely necessary in microwave frequency band research brightness temperature.In addition, the present normally used frequency range of satellite communication is S-band, X-band and Ka wave band.The lunar exploration satellite when the ground receiving antenna points to lunar exploration satellite reception detection and measuring and control data, just must be introduced the influence of moon radio radiation so around moonwork.Therefore must consider of the influence of moon radio radiation, promptly to receiving the increase of system noise temperature to antenna reception satellite downlink data signal.The research of brightness temperature also has important effect for the design of satellite communication link.
At present, the ground radio telescope is mainly adopted in the research of moon brightness temperature, and the radar that microwave radiometer that ground radar, circumlunar satellite carry and circumlunar satellite carry carries out.Wherein the ground radio telescope mainly contains the brightness temperature mean value of the whole moon of observation or the brightness temperature value of lunar surface central point to moon temperature observation research, and moon brightness temperature space distribution.In the Study on spatial distribution to moon brightness temperature, utilize the observation of single antenna radio telescope is main approach always.
Utilize single antenna that extanded source is scanned observation, can obtain observing the radio map picture of extanded source.Antenna is in pitching or azimuth direction scanning, and scan mode difference, data processing method are also different.Under the condition that the antenna automaticity is not high, adopted is grid scanning more in early days, and promptly antenna is observed the position one by one, the antenna temperature that record antenna direction each day empty position obtains.Draw the interior brightness temperature distribution plan in zone of surveying by these grid point temperature values through difference.To the scanning observation procedure of the moon, except grid scanning observation, be exactly OTF (on-the-sky) method.In the OTF method, radio telescope in the orientation or pitch orientation open certain angle partially, with the inswept a slice of certain speed observation area, and will handle after the data recording, obtain in the observation area brightness temperature and distribute.And extanded source of single antenna observation can not obtain actual extanded source Luminance Distribution, distributes and just obtain restrainting mean flow rate.How to distribute from the bundle mean flow rate of observing, recover to obtain the brightness temperature distribution of extanded source reality, this always is the problem of main research in radio astronomy.
The moon is exactly an extanded source, and the footpath, visual angle is bigger relatively.Fig. 1 is the principle schematic of the single antenna observation moon.A is that the actual brightness temperature of the moon distributes, and b is an antenna radiation pattern, and c distributes for the antenna temperature that lunar observation is obtained.In prior art, generally adopt direct one-tenth figure method to obtain moon brightness temperature space distribution to moon ground based observa tion.The moon brightness temperature figure (c) that obtains with the direct method of measurement is the convolution of true brightness temperature figure of the moon (a) and antenna radiation pattern (b), that is:
Wherein (θ φ) is the observed ray of antenna, and the useful area of antenna is on this direction
Be the brightness temperature distribution of moon reality,
Be antenna radiation pattern,
Be the antenna temperature that observation obtains, λ is the observation wavelength.As seen, the moon brightness temperature figure that is obtained not is the distribution of real moon brightness temperature.
Prior art is not consider the influence of antenna radiation pattern to cause the bright temperature precision of images of the moon that is become low to the shortcoming of the disposal route of single antenna moon brightness temperature observation data, can not reflect the real change of moon brightness temperature.So expect to have a kind of single antenna moon brightness temperature observation data disposal route, make it possible to obtain real moon brightness temperature distribution plan.
Summary of the invention
(1) technical matters that will solve
In order to overcome the above-mentioned defective of prior art, the invention provides a kind of ground single antenna observation moon brightness temperature processing method of data.
(2) technical scheme
Ground single antenna observation moon brightness temperature processing method of data provided by the invention comprises step: utilize the ground single antenna to observe to obtain moon brightness temperature observation data; Described observation data is calibrated to obtain dirty figure; Utilize described ground single antenna that the known radio source of flux density is observed to obtain antenna radiation pattern, obtain normalized dirty bundle and clean bundle according to described antenna radiation pattern; According to described dirty figure, dirty bundle and clean bundle, utilize the clean CLEAN algorithm of revising that dirty figure is carried out clean, obtain clean figure.
Wherein the clean algorithm further comprises: step 401, and (wherein (x y) is M the set that intensity is the point source of J to O, and wherein M is a natural number for x, the y) initial value of Jacobian matrix as O with the equal-sized empty matrix of dirty figure in foundation; Step 402 finds the intensity I of brightness temperature maximum point in the dirty figure
MaxAnd position (X
Max, Y
Max), described position is the row in the image array and the position of row; Step 403 multiply by ring gain g with the maximal value of finding, with result of product g*I
MaxAdd to δ Jacobian matrix relevant position (Xmax, Ymax), g＜1 wherein; Step 404 multiply by dirty bundle B with described result of product
d, and with g*I as a result
Max* the Bd integral translation is to (the X of maximal value place of step 402 discovery
Max, Y
Max), again the gained result is cut from dirty figure, obtain residual figure R, return step 402, the residual figure that utilizes step 404 to obtain replaces dirty figure, promptly new dirty figure I
New(x, y)=(x y), so repeats to the iterations N of setting to R always, and perhaps the variance of residual figure is less than setting value, and wherein N is a natural number; Step 405, (x is y) with clean bundle convolution with a series of δ function O of obtaining; Step 406, the result with the residual figure that obtains at last is added to step 405 forms clean figure at last.
(3) beneficial effect
Use ground single antenna observation moon brightness temperature processing method of data of the present invention, wherein introduced effective deconvolution CLEAN algorithm, the feasible influence of having eliminated antenna radiation pattern to a certain extent, thus obtain and the more approaching result of the actual bright temperature distribution of the moon.
Description of drawings
Fig. 1 be in the prior art single antenna to the principle schematic of lunar observation;
Fig. 2 is the process flow diagram of ground single antenna moon sight brightness temperature processing method of data of the present invention;
Fig. 3 A illustration the synoptic diagram that the lunar observation data are carried out the background match; And
Fig. 3 B illustration the residual plot of one order polynomial background match.
Embodiment
For making the purpose, technical solutions and advantages of the present invention clear more understandable, the specific embodiment of the invention is described below in conjunction with accompanying drawing.Be used to explain the present invention in this illustrative example of the present invention, but not as a limitation of the invention.
With reference to Fig. 2, Fig. 2 is the process flow diagram of ground single antenna moon sight brightness temperature processing method of data of the present invention.Realization principle of the present invention is to utilizing single antenna to observe acquisition moon brightness temperature observation raw data handle, and the calibration of passing through this raw data is handled, obtain dirty figure, dirty bundle and clean bundle, utilize the dirty figure, dirty bundle and the clean bundle that obtain, obtain the actual bright temperature distribution plan of the moon by the Clean algorithm.
At first, describe the distribute procurement process of dirty figure of moon brightness temperature in detail with reference to Fig. 2, comprising following steps:
Utilize single antenna to observe and obtain moon brightness temperature observation raw data.Moon brightness temperature observation raw data is that antenna carries out the data that pitching scanning or azimuth scan obtain.Radio telescope is made up of antenna system, receiver system and terminal system.Wherein antenna system is the topmost ingredient of radio telescope, and its effect is the radiation from the sky observation area of collecting selectively; The effect of receiver system is the radiation signal in the limited frequency band scope that amplifies from antenna output; The effect of computer system is to show, write down, process and storage radiometer output data, and has the function of control telescope each several part co-ordination.
Understand for convenience by the single antenna data layout that observation obtains to the moon, provide with the form of giving an example below and utilize Miyun 50 meters single antenna in Beijing to carry out observation station to get observation data:
| Time | Scan mode | Aa | Ea | Ar | Er | Ad | Ed | Rv |
| 4:30:41 | e- | 0.9956 | 42.8223 | 209.2484 | 67.0085 | 208.2528 | -24.1852 | 2.249 |
| 4:30:42 | e- | 1.0001 | 42.8233 | 209.2446 | 67.0101 | 208.2445 | 24.1868 | 2.084 |
The implication of each parameter is in the form: e-a-represents to observe scan mode, and e is pitching scanning, and a is an azimuth scan; Aa represents the antenna bearingt actual value of this observation station; Ea represents the antenna pitching actual value of this observation station; Ar represents the antenna bearingt theoretical value of this observation station; Er represents the antenna pitching theoretical value of this observation station; Ad represents the poor of this observation station antenna bearingt theoretical value and actual value; Ed represents the poor of this observation station antenna pitching theoretical value and actual value; Level value when Rv represents this observation station of antenna alignment.
This step at data are observation Rv in the raw data.Bad data points is meant that those significantly depart from the data of range of normal value.The appearance of these data can make moon brightness temperature distribution plan big deviation occur.The method of removing bad data points is: the data of selecting certain-length, calculate the root mean square (rms) of this segment data, then this segment data is carried out fitting of a polynomial, obtain the rms of residual error data after the match, check the residual error of every bit in this segment data then, if its absolute value greater than 5 times rms (being threshold value above-mentioned), then removes this point or replaces with match value.
When carrying out the background match, at first determine the exponent number of fitting of a polynomial, select the data area of match then.After selected data carried out fitting of a polynomial, fitting data is deducted from observation data, obtain the data behind the background.
With reference to Fig. 3 A and 3B,, the detailed process of background match is described with the data instance of certain observation lunar surface scanning.Thin curve representation among Fig. 3 A promptly carries out bad data points to observation data and removes data afterwards the scanning curve of target source, and thick dotted line is represented once polynomial background matched curve.Shown in Fig. 3 A, select the data outside 26 to the 79 time point scopes to carry out match.Then select the polynomial fitting exponent number.
In MATLAB software can be in the hope of the function expression of a polynomial fitting curve and the fiducial interval of every coefficient:
y＝f(x)＝p
1x+p
2 (3.2)
Wherein, P1 and P2 are a multinomial coefficient (95% fiducial interval):
p1＝0.0002959(0.0002457，0.0003461)
p2＝2.152(2.149，2.155)
By Fig. 3 B as can be known, the residual error of a polynomial fitting curve represents that this order polynomial can reflect the funtcional relationship of curve inherence more exactly in (0.02,0.02) scope.For the system drifting of complexity, also can select for use quadratic polynomial to carry out match.
In observation process, caused the error in pointing of antenna owing to many-sided reason, in order to eliminate the error in pointing of antenna, we adopt linear method to set up the single antenna correction model, calculate the deviate of antenna bearingt and pitching, this deviate is observation position poor of the observation position of the actual aligning of antenna and Theoretical Calculation.We deduct the resulting deviation of model with the position of actual observation, obtain the observation position of Theoretical Calculation.Linear correction model is as follows:
ΔA＝C
1+C
3/cos(E
r)+C
4cos(A
r)tan(E
r)+C
5sin(A
r)tan(E
r)+C
6tan(E
r)
ΔE＝C
2-C
4sin(A
r)+C
5cos(A
r)+C
7cos(E
r)+C
8cos
3(E
r)
In the formula, Δ A and Δ E are respectively the linear model correction functions of orientation and pitching, and E and A are respectively the pitching and the orientation of actual observation, Δ A-azimuth angle error, Δ E-angle of pitch error, C
1-azimuth angular encoders zero-error, C
2-angle of pitch encoder zero-crossing error, C
3The non-orthogonal error of-optical axis or electric axis and pitch axis, C
4The error that-azimuth axis East and West direction deflection causes, C
5The error that-azimuth axis north-south deflection causes, C
6The non-orthogonal error in-pitch axis and axle center, orientation, C
7The gravity deformation error of-electric axis is item once, C
8The high-order term of-electric axis gravity deformation error.
If the horizontal coordinate of target source be (A, z), wherein A is the position angle, z is the angle of pitch.Time angular coordinate after the conversion be (t, δ), wherein t is a hour angle, δ is a declination.Further be converted to equatorial coordinate again for (α, δ), wherein α is a right ascension.
cosδsin?t＝sin?z?sin?A
cosδcos?t＝sinδsin?z?cos?A+cos?z?cosδ
α＝LST-t
Other parameter
Be the geographic latitude of observation station, LST is local sidereal time.
In data processing, the coordinate of observation all is transformed into moon right ascension declination position when beginning to scan.Concrete method is the motion on right ascension declination direction when deducting moon ball center and beginning to observe relatively of the right ascension declination coordinate with record.Right ascension declination conversion process is formulated as follows:
RADif＝RA-RASta
RACor＝RAObs-RADif
REDif＝RE-RESta
RECor＝REObs-REDif
Wherein RASta is a moon center right ascension coordinate when beginning to scan, and RA is month real-time right ascension coordinate of ball center, and RAObs is the right ascension coordinate of each analyzing spot.RADif is the right ascension position deviation that brings owing to moon motion, and RACor is through the right ascension behind the motion correction.RESta is a moon center declination coordinate when beginning to scan, and RE is month real-time declination coordinate of ball center, and REObs is the declination coordinate of each analyzing spot.REDif is the declination position deviation that brings owing to moon motion, and RECor is through the declination behind the motion correction.
After the background match, data are calibrated.Calibration is the process that measured value is converted to the astrophysics amount, in the method, is that the magnitude of voltage Rv that will measure is by being converted to moon brightness temperature T
MoonIn astronomical sight, employed reference source is the fine and close source of known mass flow density.
When calibrating, earlier known mass flow density radio source is scanned the record value RV when obtaining antenna and departing from radio source
Off_calAnd the record value RV when aiming at radio source
On_calIn step 201 and 202, the radio source observation data is removed bad data points and background process of fitting treatment, concrete processing procedure and step 102 and 103 similar do not repeat them here.Behind the scanning moon, obtain departing from moon record value RV
Off_moonWith aligning lunar surface record value RV
On_moon, have
(RV
on_cal-RV
off_cal)C＝S
cal
(RV
on_moon-RV
off_moon)C＝S
moon
S wherein
CalBe the flux density of reference source, S
MoonBe the flux density of the moon, C is a constant.Because reference source is the known radio source of flux density, thereby obtains:
Because antenna beam is less than moon angular diameter, Δ Ω is the antenna beam solid angle, T
MoonShu Pingjun brightness temperature for observation station.
In step 108, the lunar surface observation coordinate (right ascension and declination) that step 105 is obtained becomes figure with the lunar surface brightness temperature value that step 107 obtains, and obtains dirty figure.The dirty figure that obtains is the distribution of coordinate system under the line of moon brightness temperature.And scan pattern is orientation or pitching scanning, and this has just caused observed reading is equal interval sampling in the horizontal system of coordinates, then is not equal interval sampling under the line in the coordinate system.We equidistantly carry out interpolation to observation brightness temperature data with the right ascension declination, obtain dirty figure at last.
Continuation is described the acquisition process of dirty bundle and clean bundle with reference to Fig. 2.
In step 203, the antenna radiation pattern that obtains is carried out normalization obtain dirty bundle.
Antenna radiation pattern mainly comprises main lobe, secondary lobe, the back lobe in addition that has.Antenna is exactly the directional diagram of this antenna to the response of radio point source.Antenna is at the response basically identical of azimuth plane and pitching face, and symmetry is all relatively good.For simplicity, do rotation, obtain antenna three-dimensional figure with the normalized response of azimuth plane or pitching face.Seeking the maximal value of radio source scanning in the scan-data of radio source, is the center with the right ascension or the declination position of this point, calculates the angle at radio source scanned offset center, and carries out volume normalization, obtains dirty bundle.
In step 204, the main lobe of dirty bundle is carried out Gauss curve fitting, obtain clean bundle.The purpose that obtains clean bundle mainly is to consider the effect of main lobe, the influence of deduction secondary lobe (even back lobe).
After obtaining dirty figure, dirty bundle and clean bundle, in step 106, obtain clean figure, i.e. the bright temperature distribution plan of the Shi Ji moon by Clean clean algorithm.
Because the brightness temperature distribution plan that the direct method of measurement obtains is the result of true brightness temperature distribution and antenna radiation pattern convolution,, adopt the CLEAN algorithm that it is optimized in order to eliminate the influence of antenna radiation pattern.In step 109, adopt the CLEAN algorithm of revising that dirty figure is carried out clean.
The CLEAN algorithm be with target O (therefore x, the set that y) to regard M intensity as be the point source of J can with goal description be with the set of δ pulse:
If dirty bundle is normalized antenna radiation pattern is B
d(x y), uses B
d(x, y) convolution O (x, y), then obtain dirty figure D (x y), can be expressed as:
Equally, can be with the set goal reconstructed object I
Clean(x y) is expressed as:
I wherein
Clean(x y) is called clean figure, B
c(x y) is called normalized clean bundle, R
Final(x y) is final residual figure.I
Clean(x is as a true estimation to target image y).
The concrete implementation step of CLEAN deconvolution process is as follows:
1, the equal-sized empty matrix of foundation and dirty figure is as O (x, y) initial value of Jacobian matrix;
2, find the intensity I of brightness temperature maximum point in the dirty figure
MaxAnd position (X
Max, Y
Max), described position is the row in the image array and the position of row;
3, the maximal value of finding be multiply by ring gain g (g＜1), with result of product g*I
MaxAdd δ Jacobian matrix relevant position (X to
Max, Y
Max);
4, this result of product be multiply by normalized dirty bundle B
d, and with g*I as a result
Max* B
d(the X of maximal value place that integral translation is found to step 2
Max, Y
Max), again the gained result is cut from dirty figure, obtain residual figure R, promptly
R(x，y)＝I(x，y)-g*I
max(X
max，Y
max)*B
d(x，y)；
5, return step 2, change dirty figure into the residual figure that step 4 obtains, i.e. I
New(x, y)=(x y), so repeats to the iterations N of setting to R always, and perhaps the variance of residual figure is less than setting value;
(x is y) with the clean bundle of normalization B for a series of δ function O that 6, will obtain
c(x, y) convolution, promptly
I wherein
Mod el(x y) is illustraton of model;
7, the residual figure that will obtain at last is added to the result of step 4, forms clean figure at last, thereby the brightness temperature that obtains the moon distributes, i.e. I
Clean(x, y)=I
Mod el(x, y)+R
Final(x, y).
For the CLEAN cleaning process of single antenna, be a process that iterates, set different loop limit values, can obtain different results.Set 2 parameters altogether, i.e. the ratio S of gain factor g, and residual figure maximal value and residual figure variance.The theoretical implication of g is the number percent that contains actual signal in the dirty figure maximal value, handles for the face source data, and be 0.001～0.01 according to the experience of processing value.The theoretical implication of S is the signal to noise ratio (S/N ratio) of processing signals, it is generally acknowledged that the signal of S＞3 just contains actual signal, when S＜3, thinks that the CLEAN algorithm finishes.
In the present invention, by introducing CLEAN clean algorithm ground single antenna observation moon brightness temperature data are handled, removed the influence of antenna radiation pattern effectively to the bright temperature distribution plan of the moon, thereby obtain the bright temperature distribution plan of the comparatively desirable moon, and then can study microwave property of menology etc. more objectively.
Above-described specific embodiment; purpose of the present invention, technical scheme and beneficial effect are further described; institute is understood that; the above only is specific embodiments of the invention; be not limited to the present invention; within the spirit and principles in the present invention all, any modification of being made, be equal to replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (10)
1. a ground single antenna is observed moon brightness temperature processing method of data, and this method may further comprise the steps:
Step 1 utilizes the ground single antenna to observe to obtain moon brightness temperature observation data;
Step 2 utilizes described ground single antenna that the known radio source of flux density is observed, and utilizes observed result by described moon brightness temperature observation data is calibrated processing, to obtain the dirty figure of the moon;
Step 3 utilizes described ground single antenna that the known radio source of flux density is observed to obtain antenna radiation pattern, obtains dirty bundle of normalization and clean bundle according to described antenna radiation pattern;
Step 4 according to described dirty figure, dirty bundle and clean bundle, utilizes the clean CLEAN algorithm of revising that dirty figure is carried out clean, obtains clean figure.
2. the level value when method according to claim 1, described moon brightness temperature observation data comprise the antenna pitching actual value of the antenna bearingt actual value of antenna observation scan mode, observation station correspondence and theoretical value, observation station correspondence and theoretical value, antenna alignment observation station.
3. method according to claim 2 comprises also after step 1 that wherein the level value in the observation data is carried out bad data points to be removed and the background process of fitting treatment.
4. method according to claim 2 wherein also comprises after step 1 data relevant with aerial position in the observation data is carried out antenna-point accuracy correction, the horizontal system of coordinates and equatorial system of coordinates conversion and eliminated the moon motion processing.
5. method according to claim 4, wherein the conversion of the horizontal system of coordinates and the equatorial system of coordinates further comprises and all converts the coordinate in the observation data to when beginning to scan moon right ascension declination coordinate.
6. method according to claim 5, wherein step 2 further comprises, by calibration the level value in the observation data is converted to moon brightness temperature value, utilizes described moon right ascension declination coordinate to become figure with described moon brightness temperature value, and then obtains dirty figure.
7. method according to claim 6 finds the maximal value of radio source scanning in the scan-data of radio source, be the center with the right ascension or the declination position of this point, calculates the angle at radio source scanned offset center, and carries out volume normalization, obtains dirty bundle.
8. method according to claim 7 is carried out Gauss curve fitting to the main lobe of described dirty bundle, obtains clean bundle.
9. method according to claim 6, the reference source that described calibration is adopted are the fine and close sources of known mass flow density.
10. according to each described method of claim 1-9, wherein step 4 further comprises:
Step 401, (wherein (x y) is M the set that intensity is the point source of J to O for x, the y) initial value of Jacobian matrix as O with the equal-sized empty matrix of dirty figure in foundation;
Step 402 finds the intensity I of brightness temperature maximum point in the dirty figure
MaxAnd position (X
Max, Y
Max), described position is the row in the image array and the position of row;
Step 403 multiply by ring gain g with the maximal value of finding, with result of product g*I
MaxAdd δ Jacobian matrix relevant position (X to
Max, Y
Max), g＜1 wherein;
Step 404 multiply by dirty bundle B with described result of product
d, and with g*I as a result
Max* B
d(the X of maximal value place that integral translation is found to step 402
Max, Y
Max), again the gained result is cut from dirty figure, obtain residual figure R, return step 402, the residual figure that utilizes step 404 to obtain replaces dirty figure, promptly obtains new dirty figure I
New(x, y)=(x y), so repeats to the iterations N of setting to R always, and perhaps the variance of residual figure is less than setting value, and wherein N is a natural number;
Step 405, (x is y) with clean bundle convolution with a series of δ function O of obtaining;
Step 406, the result with the residual figure that obtains at last is added to step 405 forms clean figure at last.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010297361 CN101976297B (en) | 2010-09-30 | 2010-09-30 | Processing method of moon brightness temperature data observed by foundation single antenna |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010297361 CN101976297B (en) | 2010-09-30 | 2010-09-30 | Processing method of moon brightness temperature data observed by foundation single antenna |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101976297A true CN101976297A (en) | 2011-02-16 |
| CN101976297B CN101976297B (en) | 2012-09-26 |
Family
ID=43576182
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 201010297361 Active CN101976297B (en) | 2010-09-30 | 2010-09-30 | Processing method of moon brightness temperature data observed by foundation single antenna |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101976297B (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103353922A (en) * | 2013-06-21 | 2013-10-16 | 中国科学院紫金山天文台 | OTF (optical transfer function) observation scanning method |
| CN103512673A (en) * | 2012-06-28 | 2014-01-15 | 中国科学院国家天文台 | Method for processing moon brightness and temperature data observed by microwave sounder |
| CN106405558A (en) * | 2016-11-09 | 2017-02-15 | 武汉华梦科技有限公司 | Ground-based meteorological microwave radiometer system based on dual-frequency antenna and measurement method |
| CN106908662A (en) * | 2017-02-24 | 2017-06-30 | 中国科学院上海天文台 | A kind of high frequency continuous spectrum scanning observation procedure based on dualbeam difference |
| CN108665426A (en) * | 2018-04-02 | 2018-10-16 | 上海应用技术大学 | Radio astronomy image Fast Restoration method |
| CN109859135A (en) * | 2019-01-31 | 2019-06-07 | 北京邮电大学 | A kind of image enhancement processing method applied to relevance imaging |
| CN112162158A (en) * | 2020-09-03 | 2021-01-01 | 国家卫星气象中心(国家空间天气监测预警中心) | Method and device for evaluating pointing mismatch of antenna of on-orbit terahertz detector |
| CN115690244A (en) * | 2022-09-21 | 2023-02-03 | 贵州大学 | High dynamic range reconstruction method for sparse interference array |
| CN117647316A (en) * | 2023-11-08 | 2024-03-05 | 中山大学 | Method and system for calculating moon observation brightness temperature of satellite-borne microwave radiometer |
| CN117647316B (en) * | 2023-11-08 | 2024-06-11 | 中山大学 | Method and system for calculating moon observation brightness temperature of satellite-borne microwave radiometer |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| ZA200303429B (en) * | 2000-11-24 | 2004-05-05 | Breath Ltd | Sterilisation of pharmaceuticals. |
| CN101424741A (en) * | 2008-12-08 | 2009-05-06 | 中国海洋大学 | Real time extracting method for satellite remote sensing sea fog characteristic quantity |
| CN101452078A (en) * | 2008-12-30 | 2009-06-10 | 国家卫星气象中心 | Daytime and nighttime sea fog detecting method based on polarorbiting meteorological satellite remote sense |
-
2010
- 2010-09-30 CN CN 201010297361 patent/CN101976297B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| ZA200303429B (en) * | 2000-11-24 | 2004-05-05 | Breath Ltd | Sterilisation of pharmaceuticals. |
| CN101424741A (en) * | 2008-12-08 | 2009-05-06 | 中国海洋大学 | Real time extracting method for satellite remote sensing sea fog characteristic quantity |
| CN101452078A (en) * | 2008-12-30 | 2009-06-10 | 国家卫星气象中心 | Daytime and nighttime sea fog detecting method based on polarorbiting meteorological satellite remote sense |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103512673A (en) * | 2012-06-28 | 2014-01-15 | 中国科学院国家天文台 | Method for processing moon brightness and temperature data observed by microwave sounder |
| CN103512673B (en) * | 2012-06-28 | 2016-01-20 | 中国科学院国家天文台 | A kind of Microwave sounder observation moon brightness temperature data processing method |
| CN103353922A (en) * | 2013-06-21 | 2013-10-16 | 中国科学院紫金山天文台 | OTF (optical transfer function) observation scanning method |
| CN103353922B (en) * | 2013-06-21 | 2016-09-21 | 中国科学院紫金山天文台 | A kind of OTF observes scan method |
| CN106405558A (en) * | 2016-11-09 | 2017-02-15 | 武汉华梦科技有限公司 | Ground-based meteorological microwave radiometer system based on dual-frequency antenna and measurement method |
| CN106908662A (en) * | 2017-02-24 | 2017-06-30 | 中国科学院上海天文台 | A kind of high frequency continuous spectrum scanning observation procedure based on dualbeam difference |
| CN108665426A (en) * | 2018-04-02 | 2018-10-16 | 上海应用技术大学 | Radio astronomy image Fast Restoration method |
| CN109859135A (en) * | 2019-01-31 | 2019-06-07 | 北京邮电大学 | A kind of image enhancement processing method applied to relevance imaging |
| CN109859135B (en) * | 2019-01-31 | 2021-05-07 | 北京邮电大学 | Image enhancement processing method applied to associated imaging |
| CN112162158A (en) * | 2020-09-03 | 2021-01-01 | 国家卫星气象中心(国家空间天气监测预警中心) | Method and device for evaluating pointing mismatch of antenna of on-orbit terahertz detector |
| CN115690244A (en) * | 2022-09-21 | 2023-02-03 | 贵州大学 | High dynamic range reconstruction method for sparse interference array |
| CN117647316A (en) * | 2023-11-08 | 2024-03-05 | 中山大学 | Method and system for calculating moon observation brightness temperature of satellite-borne microwave radiometer |
| CN117647316B (en) * | 2023-11-08 | 2024-06-11 | 中山大学 | Method and system for calculating moon observation brightness temperature of satellite-borne microwave radiometer |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101976297B (en) | 2012-09-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101976297B (en) | Processing method of moon brightness temperature data observed by foundation single antenna | |
| Kalberla et al. | The leiden/argentine/bonn (lab) survey of galactic hi-final data release of the combined lds and iar surveys with improved stray-radiation corrections | |
| Hajj et al. | COSMIC GPS ionospheric sensing and space weather | |
| Ushio et al. | A Kalman filter approach to the Global Satellite Mapping of Precipitation (GSMaP) from combined passive microwave and infrared radiometric data | |
| Carlstrom et al. | Interferometric imaging of the sunyaev-zeldovich effect at 30 ghz | |
| Chang et al. | Beam calibration of radio telescopes with drones | |
| Xie et al. | Design and data processing of China's first spaceborne laser altimeter system for earth observation: GaoFen-7 | |
| Ovodenko et al. | Investigating range error compensation in UHF radar through IRI-2007 real-time updating: Preliminary results | |
| Le Vine et al. | The Aquarius simulator and cold-sky calibration | |
| CN104101297A (en) | Space object dimension acquisition method based on photoelectric observation | |
| Spider Collaboration | A Constraint on Primordial B-modes from the First Flight of the Spider Balloon-borne Telescope | |
| Yang et al. | Special Sensor Microwave Imager (SSM/I) intersensor calibration using a simultaneous conical overpass technique | |
| Reul et al. | Earth-viewing L-band radiometer sensing of sea surface scattered celestial sky radiation—Part II: Application to SMOS | |
| Ellis | Galactic Radio Emission Below 16· 5 MHz and the Galactic Emission Measure | |
| CN109782280B (en) | Method for measuring distance direction directional diagram of medium and high orbit SAR antenna based on calibration satellite | |
| Li et al. | Intercalibration of advanced Himawari imager's infrared channels with IASI/Metop-B 1C data | |
| de Bernardis et al. | A short-wavelength measurement of the cosmic background radiation anisotropy | |
| Piles Guillem | Multiscale soil moisture retrievals from microwave remote sensing observations | |
| Gao et al. | Radiometric cross-calibration of GF-4/VNIR sensor with Landsat8/OLI, Sentinel-2/MSI, and Terra/MODIS for monitoring its degradation | |
| Lang et al. | Very large array observations of solar active regions | |
| Yang et al. | ATMS Radiance Data Products’ Calibration and Evaluation | |
| Hu et al. | Scan planning optimization for 2-D beam scanning using a future geostationary microwave radiometer | |
| Lukashin et al. | CLARREO reflected solar spectrometer: restrictions for instrument sensitivity to polarization | |
| Feng et al. | Data processing and error analysis for the CE-1 Lunar microwave radiometer | |
| Tan et al. | Improving IPATS channel-to-channel registration assessment |
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 |