CN101976297B - 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
- CN101976297B CN101976297B CN 201010297361 CN201010297361A CN101976297B CN 101976297 B CN101976297 B CN 101976297B CN 201010297361 CN201010297361 CN 201010297361 CN 201010297361 A CN201010297361 A CN 201010297361A CN 101976297 B CN101976297 B CN 101976297B
- 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.)
- Active
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 on moonscape, being covered with 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 is different, and data processing method is also different.Under the condition that the antenna automaticity is not high in early days, adopted is grid scanning more, and promptly antenna is observed the position one by one, the antenna temperature that record antenna direction each day empty position obtains.Draw the brightness temperature distribution plan in the 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 with handling after the data recording, obtain that brightness temperature distributes in the observation area.And extanded source of single antenna observation can not obtain actual extanded source Luminance Distribution, distributes and just obtain restrainting mean flow rate.How from bundle mean flow rate distributions of observation, recover to obtain the actual brightness temperature distribution of extanded source, 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; φ) be the observed ray of antenna; The useful area of antenna is
on this direction is the actual brightness temperature distribution of the moon;
is antenna radiation pattern; The antenna temperature that
obtains for observation, λ is the observation wavelength.It is thus clear that 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 expectation has a kind of single antenna moon brightness temperature observation data disposal route, makes it possible to obtain real moon brightness temperature distribution plan.
Summary of the invention
The technical matters that (one) 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; Said observation data is calibrated to obtain dirty figure; Utilize said 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 said antenna radiation pattern; According to said 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), said position is row and the position of row in the image array; 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 said 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 that obtain; 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 object of the invention, technical scheme and advantage 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 to qualification of the present 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 through 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:
Step 101 obtains the moon and radio source observation raw data.
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 ease through 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-representes to observe scan mode, and e is pitching scanning, and a is an azimuth scan; Aa representes the antenna bearingt actual value of this observation station; Ea representes the antenna pitching actual value of this observation station; Ar representes the antenna bearingt theoretical value of this observation station; Er representes the antenna pitching theoretical value of this observation station; Ad representes the poor of this observation station antenna bearingt theoretical value and actual value; Ed representes the poor of this observation station antenna pitching theoretical value and actual value; Level value when Rv representes this observation station of antenna alignment.
Step 102 is removed the bad data points in the lunar observation raw data.
The data that this step is directed against are the Rv in the observation 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.
Step 103, the background match.The data that step 103 is directed against are through removing bad data points lunar observation data afterwards.The purpose of background match is in order to carry out the calibration of observation background, removes the factors such as background drift that background causes owing to the reasons such as drift of receiver.There is certain instability in observation receiving equipments such as receiver, make observation data occur wafing, wafing down, even the phenomenon of periodic vibration, therefore before carrying out follow-up data processing, need remove the gradual composition in the data.Adopt fitting of a polynomial to carry out the background match.
When carrying out the background match, at first confirm 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 each item 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)
Can know that by Fig. 3 B the residual error of a polynomial fitting curve representes that this order polynomial can reflect the funtcional relationship that curve is inherent more exactly in (0.02,0.02) scope.For the system drifting of complicacy, also can select for use quadratic polynomial to carry out match.
Step 104 is carried out pointing accuracy to observation data and is proofreaied and correct.The data that this step is directed against are data relevant with aerial position, for example observe Ar and Er in the raw data, i.e. the orientation of antenna and pitching coordinate.
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, and this deviate is observation position poor of observation position and the Theoretical Calculation of the actual aligning of antenna.We deduct the resulting deviation of model with the position of actual observation, obtain the observation position of Theoretical Calculation.Linear correction model is following:
Δ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.
Step 105 after carrying out the pointing accuracy correction, carried out the conversion of coordinate system.The position of antenna is to write down with orientation under the horizontal system of coordinates and pitching coordinate in the raw data file; Need to convert this coordinate system into hour angle system of coordinate earlier; Finally convert right ascension and declination coordinate under the equatorial system of coordinates into, thereby obtain the brightness temperature data of menology diverse location.
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 convert equatorial coordinate again into 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
is the geographic latitude of observation station, and LST is local sidereal time.
Step 106 after carrying out coordinate transform, is carried out the elimination of moon motion.In the scanning observation of the moon, generally to spend 2～4 hours time, during this period, the motion of the moon on celestial sphere can not ignore.
In data processing, the coordinate of observation all is transformed into the 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.
Step 107, the calibration of data.
After the background match, data are calibrated.Calibration is the process that measured value is converted to the astrophysics amount, in the method, is through being converted to moon brightness temperature T with measured voltage value Rv
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 are similar with 103, repeat no more at this.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.In the scan-data of radio source, seeking the maximal value of radio source scanning, 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, the promptly actual bright temperature distribution plan of the moon through 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 use the set of δ pulse with goal description to be:
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).
CLEAN deconvolution process practical implementation step is following:
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), said position is row and the position of row in the image array;
3, the maximal value of finding 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 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;
The residual figure that 7, 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; Through 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; The object of the invention, technical scheme and beneficial effect have been carried out further explain, and institute it should be understood that the above is merely specific embodiment of the present invention; Be not limited to the present invention; All within spirit of the present invention and principle, 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 said ground single antenna that the known radio source of flux density is observed, and utilizes observed result through said moon brightness temperature observation data is calibrated processing, to obtain the dirty figure of the moon;
Step 3 utilizes said 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 said antenna radiation pattern;
Step 4 according to said 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, said moon brightness temperature observation data comprise antenna observation scan mode, observation station corresponding antenna orientation actual value and theoretical value, observation station corresponding antenna pitching actual value 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, converts the level value in the observation data into moon brightness temperature value through calibration, utilizes said moon right ascension declination coordinate to become figure with said 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 said dirty bundle, obtains clean bundle.
9. method according to claim 6, the reference source that said 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, 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), said position is row and the position of row in the image array;
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 said 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 that obtain;
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 CN101976297A (en) | 2011-02-16 |
| CN101976297B true 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) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103512673B (en) * | 2012-06-28 | 2016-01-20 | 中国科学院国家天文台 | A kind of Microwave sounder observation moon brightness temperature data processing method |
| CN103353922B (en) * | 2013-06-21 | 2016-09-21 | 中国科学院紫金山天文台 | A kind of OTF observes scan method |
| CN106405558B (en) * | 2016-11-09 | 2019-05-21 | 武汉华梦科技有限公司 | Ground meteorology microwave radiation meter systems and measurement method based on dual-band antenna |
| CN106908662B (en) * | 2017-02-24 | 2019-03-08 | 中国科学院上海天文台 | A kind of high frequency continuous spectrum scanning observation method based on dualbeam difference |
| CN108665426A (en) * | 2018-04-02 | 2018-10-16 | 上海应用技术大学 | Radio astronomy image Fast Restoration method |
| CN109859135B (en) * | 2019-01-31 | 2021-05-07 | 北京邮电大学 | Image enhancement processing method applied to associated imaging |
| CN112162158B (en) * | 2020-09-03 | 2022-09-20 | 国家卫星气象中心(国家空间天气监测预警中心) | Method and device for evaluating pointing mismatch of antenna of on-orbit terahertz detector |
| CN115690244B (en) * | 2022-09-21 | 2023-07-21 | 贵州大学 | High dynamic range reconstruction method for sparse interference array |
| 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 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101976297A (en) | 2011-02-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101976297B (en) | Processing method of moon brightness temperature data observed by foundation single antenna | |
| Carlstrom et al. | Interferometric imaging of the sunyaev-zeldovich effect at 30 ghz | |
| Price et al. | Midcourse space experiment survey of the galactic plane | |
| Li et al. | The global image of the Moon obtained by the Chang’E-1: Data processing and lunar cartography | |
| Chang et al. | Beam calibration of radio telescopes with drones | |
| Ovodenko et al. | Investigating range error compensation in UHF radar through IRI-2007 real-time updating: Preliminary results | |
| CN104101297A (en) | Space object dimension acquisition method based on photoelectric observation | |
| Yang et al. | Special Sensor Microwave Imager (SSM/I) intersensor calibration using a simultaneous conical overpass technique | |
| Peng et al. | Soil moisture active/passive (SMAP) L-band microwave radiometer post-launch calibration upgrade | |
| Nieke et al. | Ocean and Land Color Imager on Sentinel‐3 | |
| Kahan et al. | Mars precession rate determined from radiometric tracking of the InSight lander | |
| 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 | |
| Anterrieu et al. | Preliminary system studies on a high-resolution SMOS follow-on: SMOS-HR | |
| de Bernardis et al. | A short-wavelength measurement of the cosmic background radiation anisotropy | |
| Lowe et al. | 30 GHz flux density measurements of the Caltech-Jodrell flat-spectrum sources with OCRA-p | |
| Lang et al. | Very large array observations of solar active regions | |
| Yang et al. | ATMS radiance data products’ calibration and evaluation | |
| Takaku et al. | High resolution DSM generation from ALOS PRISM—Status updates on over three year operations | |
| Hu et al. | Scan planning optimization for 2-D beam scanning using a future geostationary microwave radiometer | |
| Geldzahler | Coherent uplink arraying techniques for next generation space communications and planetary radar systems | |
| Dearie et al. | The first image of the Milky Way's central black hole and the unique enhancement Africa could offer future tests of gravity | |
| Tan et al. | Improving IPATS channel-to-channel registration assessment | |
| Abraham et al. | Background emissions during cold sky calibration of Aquarius | |
| Keller et al. | Lunar reconnaissance orbiter: Instrument suite and objectives |
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 |