CN114034388A - Whole moon given place time brightness temperature drawing method - Google Patents
Whole moon given place time brightness temperature drawing method Download PDFInfo
- Publication number
- CN114034388A CN114034388A CN202111217891.2A CN202111217891A CN114034388A CN 114034388 A CN114034388 A CN 114034388A CN 202111217891 A CN202111217891 A CN 202111217891A CN 114034388 A CN114034388 A CN 114034388A
- Authority
- CN
- China
- Prior art keywords
- data
- tetrahedron
- brightness temperature
- scatter diagram
- vertex
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 title claims abstract description 27
- 238000010586 diagram Methods 0.000 claims abstract description 64
- 238000007781 pre-processing Methods 0.000 claims abstract description 6
- 239000013598 vector Substances 0.000 claims description 23
- 238000012545 processing Methods 0.000 claims description 19
- 230000002159 abnormal effect Effects 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 239000000126 substance Substances 0.000 claims description 9
- 230000010354 integration Effects 0.000 claims description 6
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 238000003780 insertion Methods 0.000 claims description 3
- 230000037431 insertion Effects 0.000 claims description 3
- 239000004335 litholrubine BK Substances 0.000 claims description 3
- 238000013507 mapping Methods 0.000 abstract description 5
- 239000010936 titanium Substances 0.000 description 9
- 239000000463 material Substances 0.000 description 8
- 239000004973 liquid crystal related substance Substances 0.000 description 6
- 238000012986 modification Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 238000004590 computer program Methods 0.000 description 4
- 239000002689 soil Substances 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- RTAQQCXQSZGOHL-UHFFFAOYSA-N Titanium Chemical compound [Ti] RTAQQCXQSZGOHL-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 229910052719 titanium Inorganic materials 0.000 description 2
- 230000008094 contradictory effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Radiation Pyrometers (AREA)
- Complex Calculations (AREA)
Abstract
The invention relates to a full moon given place time brightness temperature mapping method, which comprises the following steps: s1, preprocessing the pre-acquired 2C-level data to acquire first data; the 2C-level data includes: a plurality of single-track data; each single track data comprises two parts of a data label and a data entity; s2, acquiring a first scatter diagram of the first data based on the first data; s3, acquiring a plurality of tetrahedrons based on the first scatter diagram of the first data; s4, based on a plurality of tetrahedrons and corresponding preset interpolation points, simultaneously carrying out time dimension and space dimension interpolation to obtain a brightness temperature value of each interpolation point; and S5, acquiring a full-month brightness temperature map at a given moment based on the brightness temperature value of each interpolation point.
Description
Technical Field
The invention relates to the technical field of moon brightness temperature mapping, in particular to a full moon brightness temperature mapping method in a given place.
Background
Microwave Radiometer (MRM) carried by ChangE No. 1 and No. 2 lunar exploration satellites obtains microwave brightness temperature data (MRM data) of an on-orbit lunar surface for the first time internationally, the data has the highest global time and spatial resolution, is very sensitive to components and temperature of lunar soil within a certain thickness range, and has important value in lunar scientific research.
However, limited by observation conditions, the current full-moon brightness temperature chart uses brightness temperature observation data of 3 months and above; at 1 month, which corresponds to 27.3 earth hours, the prepared lunar surface temperature map is disadvantageous for accurate identification of lunar soil thermal physical parameters over such a large time span. In the prior art, when the local time span of the selected MRM data of Chang 'e's second number is at least 3 and 2 months respectively when the bright temperature maps of noon and midnight are made, the long time span can cause the bright temperature maps to be interfered by the local time and have stripes and bright temperature abrupt changes along the longitude direction. These streaks and abrupt changes in the light temperature in the longitudinal direction are disadvantageous for the precise identification of the thermophysical parameters of the lunar soil. Alternatively, the prior art does not completely eliminate outliers in MRM data caused by other causes. For abnormal points which cannot be corrected by the current MRM data processing technology, the abnormal points need to be eliminated as dead points.
Disclosure of Invention
Technical problem to be solved
In view of the above disadvantages and shortcomings of the prior art, the present invention provides a method for drawing a full moon at a given time and temperature.
(II) technical scheme
In order to achieve the purpose, the invention adopts the main technical scheme that:
in one aspect, an embodiment of the present invention provides a full moon given place time brightness temperature mapping method, including:
s1, preprocessing the pre-acquired 2C-level data to acquire first data;
the 2C-level data includes: a plurality of single-track data;
each single track data comprises two parts of a data label and a data entity;
the data entity of each single track of data comprises: the observation time, the solar incident angle, the solar azimuth angle, the longitude and the latitude of each observation point on the lunar surface, the brightness temperature of a 3.0GHz channel, the brightness temperature of a 7.8GHz channel, the brightness temperature of a 19.35GHz channel and the brightness temperature of a 37GHz channel, which are collected by each flight track of the Chang' e second microwave radiometer;
s2, acquiring a first scatter diagram of the first data based on the first data;
s3, acquiring a plurality of tetrahedrons based on the first scatter diagram of the first data;
s4, acquiring a brightness temperature value of each interpolation point based on a plurality of tetrahedrons and corresponding preset interpolation points;
and S5, acquiring a full-month brightness temperature map at a given moment based on the brightness temperature value of each interpolation point.
Preferably, the first and second liquid crystal materials are,
the pretreatment comprises the following steps: sequentially carrying out data integration processing, abnormal data rejection processing and time angle conversion processing;
the data integration processing comprises: removing the data label of each single-track data to obtain a data entity in each single-track data, and integrating all the data entities of the single-track data into one file to obtain a first file;
the abnormal data eliminating treatment comprises the following steps: removing the first type data, the second type data and the third type data in the first file;
the first type of data is: data with brightness temperature values exceeding the range of 30K-400K in the first file;
the second type of data is: the brightness temperature value at the same latitude is higher than or lower than the data of the first preset threshold of the brightness temperature value in the adjacent orbit data;
the third type of data is: within the preset longitude and latitude range, the brightness temperature value of the data is different from that of other data by a second preset threshold;
the time angle conversion processing includes:
adopting a formula (1) to perform coordinate system conversion on the data in the first file after the abnormal data elimination processing to obtain first data;
the formula (1) is:
tanh＝sinasini/(cosφcosi+sinφcosasini)；
wherein h is a time angle;
a is the solar azimuth;
i is the sun incident angle;
phi is the latitude;
converting the time angle into the local time of the moon in 24 hours by adopting a formula (2);
the formula (2) is:
t＝12(180+h)/180；
t is the local time of the moon.
Preferably, the first and second liquid crystal materials are,
the first scatter diagram of the first data is a scatter diagram under a first coordinate system;
the first coordinate system is a cartesian coordinate system with longitude on the x-axis, latitude on the y-axis, and the local time of the moon on the z-axis.
Preferably, the S3 specifically includes:
s31, acquiring a second scatter diagram based on the first scatter diagram of the first data;
the second scatter diagram is obtained by adding the first scatter diagram;
the adding treatment comprises the following steps: adding first data of a first preset time period in the first scatter diagram above an upper boundary of the first scatter diagram, adding first data of a second preset time period in the first scatter diagram below a lower boundary of the first scatter diagram, adding first data of a first preset longitude section in the first scatter diagram to the left of the left boundary of the first scatter diagram, and adding first data of a second preset longitude section in the first scatter diagram to the right of the right boundary of the first scatter diagram;
the first preset time period is as follows: local time 0 to 4;
the second preset time period is as follows: local time 8 o 'clock to 24 o' clock;
the third preset longitude segment is: 170 DEG E-180 DEG E;
the fourth preset longitude segment is: -180 ° W-170 ° W;
s32, carrying out Dirony tetrahedron subdivision on the second scatter diagram to obtain a plurality of tetrahedrons.
Preferably, the first and second liquid crystal materials are,
the S32 specifically includes:
and aiming at the second scatter diagram, carrying out Dirony tetrahedron subdivision by adopting a point-by-point insertion algorithm to obtain a plurality of tetrahedrons.
Preferably, the first and second liquid crystal materials are,
the preset interpolation points are as follows: interpolation points are given every 0.5 ° in the longitude and latitude directions and every 1 month in the time direction.
Preferably, the S4 specifically includes:
s41, aiming at each tetrahedron and the interpolation point corresponding to the tetrahedron, obtaining a first proportional coefficient corresponding to each vertex of the tetrahedron;
wherein, the first scale factor of any vertex is the ratio of the vertical distance from the interpolation point corresponding to the tetrahedron to the surface corresponding to the vertex of the tetrahedron to the vertical distance from the vertex to the corresponding surface;
s42, acquiring the lighting temperature value of the interpolation point corresponding to the tetrahedron based on the first scale coefficient corresponding to each vertex of the tetrahedron and the position vector, the lighting temperature value, the longitude, the latitude and the local time corresponding to each vertex of the tetrahedron.
Preferably, the first and second liquid crystal materials are,
in S41, for each tetrahedron and the interpolation point corresponding to the tetrahedron, obtaining a first scaling factor a, b, c, d corresponding to each vertex of the tetrahedron respectively by using formula (3);
formula (3):
wherein, the four vertex position vectors of the tetrahedron are respectively P1, P2, P3 and P4;
wherein the content of the first and second substances,
p1＝(lo1,la1,t1)；p2＝(lo2,la2,t2)；p3＝(lo3,la3,t3)；
p4＝(lo4,la4,t4)。
preferably, the first and second liquid crystal materials are,
in S41, for each tetrahedron and the interpolation point corresponding to the tetrahedron, obtaining a first scaling factor a, b, c, d corresponding to each vertex of the tetrahedron respectively by using formula (4);
a＝D1/D,b＝D2/D,c＝D3/D,d＝D4/D (4)；
wherein the content of the first and second substances,
the four vertex position vectors of the tetrahedron are P1, P2, P3 and P4 respectively;
wherein the content of the first and second substances,
p1＝(lo1,la1,t1)；p2＝(lo2,la2,t2)；p3＝(lo3,la3,t3)；
p4＝(lo4,la4,t4)。
preferably, the S2 specifically includes:
adopting a formula (5) to obtain a brightness temperature value of the interpolation point corresponding to the tetrahedron;
the formula (5) is:
T＝a·T1+b·T2+c·T3+d·T4；
wherein, T1 is the brightness temperature value corresponding to the tetrahedron vertex with the vertex position vector P1;
t2 is the brightness temperature value corresponding to the tetrahedron vertex with the vertex position vector P2;
t3 is the brightness temperature value corresponding to the tetrahedron vertex with the vertex position vector P3;
t4 is the brightness temperature value corresponding to the tetrahedron vertex with the vertex position vector P4.
(III) advantageous effects
The invention has the beneficial effects that: according to the method for drawing the time-temperature and the space-dimension interpolation of the full moon at the given place, the time-dimension and the space-dimension interpolation are simultaneously carried out by adopting the gravity center interpolation method based on the Dirony tetrahedron subdivision, so that the full moon light temperature is obtained.
Drawings
FIG. 1 is a flow chart of a method for plotting brightness and temperature of a given place of a full moon according to the present invention;
FIG. 2 is a scattergram of 37GHz light temperature data for multiple tracks in an embodiment of the invention;
FIG. 3 is a first scatter plot of first data in an embodiment of the present invention;
FIG. 4 is a second scattergram in an embodiment of the present invention;
FIG. 5 is a schematic diagram of any tetrahedral structure in an embodiment of the present invention;
fig. 6(a) is a schematic diagram of an interpolation result of data of a 3.0GHz channel in a three-dimensional euclidean space in the embodiment of the present invention;
fig. 6(b) is a schematic diagram of an interpolation result of data of a 7.8GHz channel in a three-dimensional euclidean space in the embodiment of the present invention;
FIG. 6(c) is a schematic diagram of an interpolation result of data of a 19.35GHz channel in a three-dimensional Euclidean space in the embodiment of the present invention;
FIG. 6(d) is a schematic diagram of an interpolation result of data of a 37GHz channel in a three-dimensional Euclidean space in the embodiment of the present invention;
fig. 7 is a 37.0GHz luminance temperature plot for 24 local times in an embodiment of the present invention.
Detailed Description
For the purpose of better explaining the present invention and to facilitate understanding, the present invention will be described in detail by way of specific embodiments with reference to the accompanying drawings.
In order to better understand the above technical solutions, exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the invention are shown in the drawings, it should be understood that the invention can be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
Referring to fig. 1, the present embodiment provides a method for drawing a brightness temperature at a given place of a full moon, including:
and S1, preprocessing the pre-acquired 2C-level data to acquire first data.
The 2C-level data includes: a plurality of single-track data.
Each single track of data includes two parts, a data tag and a data entity.
The data entity of each single track of data comprises: the observation time, the solar incident angle, the solar azimuth angle, the longitude and the latitude of each observation point on the lunar surface, the brightness temperature of a 3.0GHz channel, the brightness temperature of a 7.8GHz channel, the brightness temperature of a 19.35GHz channel and the brightness temperature of a 37GHz channel are collected by each flight track of the Chang' e second microwave radiometer.
The microwave radiometer Chang' e II observes the lunar surface for about 5000 hours from 10 and 15 days in 2010 to 5 and 20 days in 2011, obtains 2401 orbit data totally, and can completely cover the whole lunar surface for 7 times. The 2C-level data used in this embodiment includes observation time (UTC) of each observation unit on the lunar surface, solar incident angle, solar azimuth, longitude, latitude, satellite altitude, data quality, and microwave brightness temperatures (3.0, 7.8, 19.35, and 37GHz) of 4 channels. On the flight orbit of about 100km of Chang' e second satellite, the spatial resolution of the 3.0GHz channel of the microwave radiometer is about 25km, and the spatial resolution of the remaining three channels is about 17.5 km.
The 2C-level data in this embodiment is a plurality of pieces of single-track data stored in psd (nasaplanetarydatasystem) format, which do not include processed time angle information and have some data anomalies. Therefore, preprocessing is required before three-dimensional interpolation is performed.
S2, acquiring a first scatter diagram of the first data based on the first data.
S3, acquiring a plurality of tetrahedrons based on the first scatter diagram of the first data.
And S4, acquiring the brightness temperature value of each interpolation point based on a plurality of tetrahedrons and corresponding preset interpolation points.
And S5, acquiring a full-month brightness temperature map at a given moment based on the brightness temperature value of each interpolation point.
In practical application of this embodiment, the preprocessing includes: and sequentially performing data integration processing, abnormal data rejection processing and time angle conversion processing.
The data integration processing comprises: and removing the data labels of each single-track data to obtain the data entities in each single-track data, and integrating the data entities of all the single-track data into one file to obtain a first file.
The abnormal data eliminating treatment comprises the following steps: and removing the first type data, the second type data and the third type data in the first file.
The first type of data is: and data with the brightness temperature value exceeding the range of 30K-400K in the first file.
The second type of data is: and the brightness temperature value at the same latitude is higher or lower than the data of the first preset threshold value of the brightness temperature value in the adjacent orbit data.
The third type of data is: and in the preset longitude and latitude range, the brightness temperature value of the data is different from that of other data by a second preset threshold value.
Referring to fig. 2, a box a shows the second type of abnormal data, and a box B shows the third type of abnormal data.
The time angle conversion processing includes:
and (3) performing coordinate system conversion on the data in the first file subjected to the abnormal data elimination processing by adopting a formula (1) to obtain first data.
The formula (1) is:
tanh＝sinasini/(cosφcosi+sinφcosasini)。
wherein h is the hour angle.
a is the solar azimuth.
i is the angle of incidence of the sun.
Phi is the latitude.
In a specific application, the time angle range of the formula (1) in the daytime is-90 degrees, and the time angle range of the formula (1) in the evening is-180-90 degrees and 90-180 degrees.
The time angle is converted into the local time of the moon in 24 hours using equation (2).
The formula (2) is:
t＝12(180+h)/180。
t is the local time of the moon.
In a particular application, the local time of each observation point in the 2C-level data is characterized by the sun's angle of incidence and azimuth in the Horizon Coordinate System (HCS), which describes the position of the sun relative to the observation point. The solar incident angle and azimuth are latitude dependent, however, the local time and latitude of the observation point are independent. To better reflect the position of the sun relative to the observation point, it is more appropriate to convert the horizon coordinate system to the Lunar Equatorial Coordinate System (LECS).
In practical applications of this embodiment, the first scattergram of the first data is a scattergram in a first coordinate system.
The first coordinate system is a cartesian coordinate system with longitude on the x-axis, latitude on the y-axis, and the local time of the moon on the z-axis.
In a specific application, referring to fig. 3, the data in the first scatter plot is evenly distributed in "slices" at each local time, and such "slices" total about 15, i.e., 15 incomplete coverages. Overall, MRM data are uniformly distributed in time and space as a whole, which is advantageous for achieving a better three-dimensional interpolation effect.
In practical application of this embodiment, the S3 specifically includes:
s31, acquiring a second scatter diagram based on the first scatter diagram of the first data.
And the second scatter diagram is obtained by adding the first scatter diagram.
The adding treatment comprises the following steps: adding first data of a first preset time period in the first scatter diagram above an upper boundary of the first scatter diagram, adding first data of a second preset time period in the first scatter diagram below a lower boundary of the first scatter diagram, adding first data of a first preset longitude section in the first scatter diagram to the left of the left boundary of the first scatter diagram, and adding first data of a second preset longitude section in the first scatter diagram to the right of the right boundary of the first scatter diagram.
The first preset time period is as follows: local time 0 to 4.
The second preset time period is as follows: from local time 8 to 24.
The third preset longitude segment is: 170 DEG E-180 DEG E.
The fourth preset longitude segment is: -180 ° W-170 ° W.
In a specific application, since there is less MRM data near the upper boundary (24:00) in fig. 3, this may result in a lack of data points on the boundary that construct a dironi tetrahedron, thereby making interpolation on the boundary incomplete. Considering that the brightness temperature is cyclically changed in a period of 24 hours, we can solve the interpolation problem on the boundary by adopting the following steps: as shown in FIG. 4, we can add data from 0:00 to 4:00 above the upper boundary (24:00) and 20:00 to 24:00 below the lower boundary (0:00), and further, since the MRM data is distributed over this closed geometry of the moon, we add data from 170E to 180E at the left boundary (-180) and data from-180W to 170W at the right boundary (180). A second scatter plot is thus obtained, see fig. 4.
S32, carrying out Dirony tetrahedron subdivision on the second scatter diagram to obtain a plurality of tetrahedrons.
In practical application of this embodiment, the step S32 specifically includes:
and aiming at the second scatter diagram, carrying out Dirony tetrahedron subdivision by adopting a point-by-point insertion algorithm to obtain a plurality of tetrahedrons.
In practical application of this embodiment, the preset interpolation points are: interpolation points are given every 0.5 ° in the longitude and latitude directions and every 1 month in the time direction.
In practical application of this embodiment, the S4 specifically includes:
s41, obtaining a first scale coefficient corresponding to each vertex of each tetrahedron aiming at each tetrahedron and the interpolation point corresponding to the tetrahedron.
The first scale factor of any vertex is the ratio of the vertical distance from the interpolation point corresponding to the tetrahedron to the surface corresponding to the vertex of the tetrahedron to the vertical distance from the vertex to the corresponding surface.
In a specific application, assuming one of the tetrahedrons is shown in FIG. 5, its four vertices are named M1，M2，M3And M4The position vectors of the four vertices are p1＝(lo1,la1,t1)，p2＝(lo2,la2,t2)，p3＝(lo3,la3,t3) And p is4＝(lo4,la4,t4) The brightness temperature values are respectively T1,T2,T3And T4。loi,laiAnd ti(i ═ 1,2,3, and 4) represent the longitude, latitude, and local time of the four vertices, respectively.
S42, acquiring the lighting temperature value of the interpolation point corresponding to the tetrahedron based on the first scale coefficient corresponding to each vertex of the tetrahedron and the position vector, the lighting temperature value, the longitude, the latitude and the local time corresponding to each vertex of the tetrahedron.
In a specific application, the luminance temperature value T of an interpolation point M with an arbitrary vector position p ═ o, la, T satisfies:
T＝a·T1+b·T2+c·T3+d·T4。
where a, b, c, and d in the formula are point M coordinates representing the weight of the influence of the luminance values of the four vertices on the luminance value to be interpolated, and a + b + c + d is 1. The intuitive understanding of the M coordinate is the ratio of the perpendicular distance of the interpolation point to the four faces of the dironi tetrahedron to the perpendicular distance of the corresponding vertex to the four faces. E.g. h in fig. 51And H1The ratio of (a) is the value of a.
In practical application of this embodiment, in S41, for each tetrahedron and the interpolation point corresponding to the tetrahedron, the formula (3) is adopted to obtain the first scaling coefficients a, b, c, and d corresponding to each vertex of the tetrahedron.
Formula (3):
wherein, the four vertex position vectors of the tetrahedron are respectively P1、P2、P3、P4。
Wherein the content of the first and second substances,
p1＝(lo1,la1,t1)；p2＝(lo2,la2,t2)；p3＝(lo3,la3,t3)。
p4＝(lo4,la4,t4)。
in practical application of this embodiment, in S41, for each tetrahedron and the interpolation point corresponding to the tetrahedron, the formula (4) is used to obtain the first scaling coefficients a, b, c, and d corresponding to each vertex of the tetrahedron.
a＝D1/D,b＝D2/D,c＝D3/D,d＝D4/D (4)。
Wherein the content of the first and second substances,
the four vertex position vectors of the tetrahedron are respectively P1、P2、P3、P4。
Wherein the content of the first and second substances,
p1＝(lo1,la1,t1)；p2＝(lo2,la2,t2)；p3＝(lo3,la3,t3)。
p4＝(lo4,la4,t4)。
in practical application of this embodiment, the step S2 specifically includes:
and (5) acquiring the brightness temperature value of the interpolation point corresponding to the tetrahedron by adopting a formula (5).
The formula (5) is:
T＝a·T1+b·T2+c·T3+d·T4。
wherein, T1Is a vertex position vector of P1The brightness temperature value corresponding to the tetrahedron vertex.
T2Is a vertex position vector of P2The brightness temperature value corresponding to the tetrahedron vertex.
T3Is a vertex position vector of P3The brightness temperature value corresponding to the tetrahedron vertex.
T4Is a vertex position vector of P4The brightness temperature value corresponding to the tetrahedron vertex.
In the present embodiment, it should be emphasized that the interpolated brightness temperature obtained in dironi by barycentric interpolation is always located between the maximum brightness temperature value and the minimum brightness temperature value in the four vertices, that is:
min(Ti)≤T≤max(Ti),(i＝1,2,3,4)。
this is because a + b + c + d is 1, if T > max (ti), then we will get from equation (5):
a·T1+b·T2+c·T3+d·T4＞(a+b+c+d)max(Ti)
namely:
a(T1-max(Ti))+b(T2-max(Ti))+c(T3-max(Ti))+d(T4-max(Ti))＞0。
this is clearly contradictory. In the same way, T is more than min (T)i). Therefore, the high and low brightness temperature abnormal values in the brightness temperature map obtained by the barycentric interpolation are caused by the original data, which improves the reliability of the brightness temperature abnormal research.
According to the method for drawing the brightness temperature at the given place of the full moon, due to the fact that the gravity center interpolation method based on the Dirony tetrahedron subdivision is adopted to conduct time dimension and space dimension interpolation at the same time, compared with the prior art, the method can well eliminate the influence on accurate understanding of lunar soil thermal physical parameters caused by the fact that full moon drawing is conducted by mixing data of a plurality of time segments, and meanwhile, the obtained brightness temperature drawing at the given place is more accurate and complete.
Fig. 6(a), 6(b), 6(c), and 6(d) are obtained by a full moon local brightness temperature mapping method, and can be obtained as follows: as the frequency increases, the maximum light temperature increases; the brightness temperature at noon is obviously higher than that at other local time, and the brightness temperature in the early morning is the lowest; the light temperature decreases with increasing latitude.
FIG. 7 shows the 37.0GHz brightness temperature plot for 24 local times (each plot corresponds to a 1-hour 37GHz brightness temperature plot). The 24 local time intensity maps of fig. 7 all have a high degree of sharpness and smoothness, with almost no streaks, compared to the intensity maps obtained by mixing the MRM data of different local times. In addition to the 37GHz brightness temperature plots of 15:00 and 16:00, the 37GHz brightness temperature plots at other time points have little jump in brightness temperature values in longitude. Merle crates, mountains and the moon sea can be clearly identified on the light temperature map at 8 time points from 9:00 to 16: 00. This illustrates three advantages of the method in this embodiment:
(1) the bright temperature graph obtained by the embodiment has higher definition and is smoother, and interference fringes and jumps of bright temperature values in longitude are almost avoided.
(2) The local time of the bright temperature graph obtained by the embodiment is more definite, which is beneficial to comparison and research of bright temperature characteristics in a full month or a large scale.
(3) Although only the interpolation results of 24 local times shown in fig. 6 are shown in this embodiment, the present embodiment can obtain a brightness temperature map of any local time by adjusting the interpolation interval in time.
(4) The special areas, mainly the moon, the basin, the impact pit and the mountain, are continuous over 24 local time light temperature maps, which illustrates the stability of the light temperature maps generated by the method of the present embodiment.
In addition, the quiet sea (marezanquillitatis) with high titanium content exhibits a higher light temperature during the day and a lower light temperature at night than the clear sea (mareserentis) with low titanium content. This is consistent with the simulation results.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions.
It should be noted that in the claims, any reference signs placed between parentheses shall not be construed as limiting the claim. The word "comprising" does not exclude the presence of elements or steps not listed in a claim. The word "a" or "an" preceding an element does not exclude the presence of a plurality of such elements. The invention may be implemented by means of hardware comprising several distinct elements, and by means of a suitably programmed computer. In the claims enumerating several means, several of these means may be embodied by one and the same item of hardware. The use of the terms first, second, third and the like are for convenience only and do not denote any order. These words are to be understood as part of the name of the component.
Furthermore, it should be noted that in the description of the present specification, the description of the term "one embodiment", "some embodiments", "examples", "specific examples" or "some examples", etc., means that a specific feature, structure, material or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present invention. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, various embodiments or examples and features of different embodiments or examples described in this specification can be combined and combined by one skilled in the art without contradiction.
While preferred embodiments of the present invention have been described, additional variations and modifications in those embodiments may occur to those skilled in the art once they learn of the basic inventive concepts. Therefore, the claims should be construed to include preferred embodiments and all changes and modifications that fall within the scope of the invention.
It will be apparent to those skilled in the art that various modifications and variations can be made in the present invention without departing from the spirit or scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention should also include such modifications and variations.
Claims (10)
1. A full moon given place time brightness temperature drawing method is characterized by comprising the following steps:
s1, preprocessing the pre-acquired 2C-level data to acquire first data;
the 2C-level data includes: a plurality of single-track data;
each single track data comprises two parts of a data label and a data entity;
the data entity of each single track of data comprises: the observation time, the solar incident angle, the solar azimuth angle, the longitude and the latitude of each observation point on the lunar surface, the brightness temperature of a 3.0GHz channel, the brightness temperature of a 7.8GHz channel, the brightness temperature of a 19.35GHz channel and the brightness temperature of a 37GHz channel, which are collected by each flight track of the Chang' e second microwave radiometer;
s2, acquiring a first scatter diagram of the first data based on the first data;
s3, acquiring a plurality of tetrahedrons based on the first scatter diagram of the first data;
s4, acquiring a brightness temperature value of each interpolation point based on a plurality of tetrahedrons and corresponding preset interpolation points;
and S5, acquiring a full-month brightness temperature map at a given moment based on the brightness temperature value of each interpolation point.
2. The method of claim 1,
the pretreatment comprises the following steps: sequentially carrying out data integration processing, abnormal data rejection processing and time angle conversion processing;
the data integration processing comprises: removing the data label of each single-track data to obtain a data entity in each single-track data, and integrating all the data entities of the single-track data into one file to obtain a first file;
the abnormal data eliminating treatment comprises the following steps: removing the first type data, the second type data and the third type data in the first file;
the first type of data is: data with brightness temperature values exceeding the range of 30K-400K in the first file;
the second type of data is: the brightness temperature value at the same latitude is higher than or lower than the data of the first preset threshold of the brightness temperature value in the adjacent orbit data;
the third type of data is: within the preset longitude and latitude range, the brightness temperature value of the data is different from that of other data by a second preset threshold;
the time angle conversion processing includes:
adopting a formula (1) to perform coordinate system conversion on the data in the first file after the abnormal data elimination processing to obtain first data;
the formula (1) is:
tanh＝sinasini/(cosφcosi+sinφcosasini)；
wherein h is a time angle;
a is the solar azimuth;
i is the sun incident angle;
phi is the latitude;
converting the time angle into the local time of the moon in 24 hours by adopting a formula (2);
the formula (2) is:
t＝12(180+h)/180；
t is the local time of the moon.
3. The method of claim 2,
the first scatter diagram of the first data is a scatter diagram under a first coordinate system;
the first coordinate system is a cartesian coordinate system with longitude on the x-axis, latitude on the y-axis, and the local time of the moon on the z-axis.
4. The method according to claim 3, wherein the S3 specifically comprises:
s31, acquiring a second scatter diagram based on the first scatter diagram of the first data;
the second scatter diagram is obtained by adding the first scatter diagram;
the adding treatment comprises the following steps: adding first data of a first preset time period in the first scatter diagram above an upper boundary of the first scatter diagram, adding first data of a second preset time period in the first scatter diagram below a lower boundary of the first scatter diagram, adding first data of a first preset longitude section in the first scatter diagram to the left of the left boundary of the first scatter diagram, and adding first data of a second preset longitude section in the first scatter diagram to the right of the right boundary of the first scatter diagram;
the first preset time period is as follows: local time 0 to 4;
the second preset time period is as follows: local time 8 o 'clock to 24 o' clock;
the third preset longitude segment is: 170 DEG E-180 DEG E;
the fourth preset longitude segment is: -180 ° W-170 ° W;
s32, carrying out Dirony tetrahedron subdivision on the second scatter diagram to obtain a plurality of tetrahedrons.
5. The method of claim 4,
the S32 specifically includes:
and aiming at the second scatter diagram, carrying out Dirony tetrahedron subdivision by adopting a point-by-point insertion algorithm to obtain a plurality of tetrahedrons.
6. The method of claim 5,
the preset interpolation points are as follows: interpolation points are given every 0.5 ° in the longitude and latitude directions and every 1 month in the time direction.
7. The method according to claim 6, wherein the S4 specifically includes:
s41, aiming at each tetrahedron and the interpolation point corresponding to the tetrahedron, obtaining a first proportional coefficient corresponding to each vertex of the tetrahedron;
wherein, the first scale factor of any vertex is the ratio of the vertical distance from the interpolation point corresponding to the tetrahedron to the surface corresponding to the vertex of the tetrahedron to the vertical distance from the vertex to the corresponding surface;
s42, acquiring the lighting temperature value of the interpolation point corresponding to the tetrahedron based on the first scale coefficient corresponding to each vertex of the tetrahedron and the position vector, the lighting temperature value, the longitude, the latitude and the local time corresponding to each vertex of the tetrahedron.
8. The method of claim 7,
in S41, for each tetrahedron and the interpolation point corresponding to the tetrahedron, obtaining a first scaling factor a, b, c, d corresponding to each vertex of the tetrahedron respectively by using formula (3);
formula (3):
wherein, the four vertex position vectors of the tetrahedron are respectively P1、P2、P3、P4；
Wherein the content of the first and second substances,
p1＝(lo1,la1,t1)；p2＝(lo2,la2,t2)；p3＝(lo3,la3,t3)；
p4＝(lo4,la4,t4)。
9. the method of claim 7,
in S41, for each tetrahedron and the interpolation point corresponding to the tetrahedron, obtaining a first scaling factor a, b, c, d corresponding to each vertex of the tetrahedron respectively by using formula (4);
a＝D1/D,b＝D2/D,c＝D3/D,d＝D4/D (4)；
wherein the content of the first and second substances,
the four vertex position vectors of the tetrahedron are respectively P1、P2、P3、P4；
Wherein the content of the first and second substances,
p1＝(lo1,la1,t1)；p2＝(lo2,la2,t2)；p3＝(lo3,la3,t3)；
p4＝(lo4,la4,t4)。
10. the method according to claim 8 or 9, wherein S2 specifically includes:
adopting a formula (5) to obtain a brightness temperature value of the interpolation point corresponding to the tetrahedron;
the formula (5) is:
T＝a·T1+b·T2+c·T3+d·T4；
wherein, T1Is a vertex position vector of P1The brightness temperature value corresponding to the tetrahedron vertex;
T2is a vertex position vector of P2The brightness temperature value corresponding to the tetrahedron vertex;
T3is a vertex position vector of P3The brightness temperature value corresponding to the tetrahedron vertex;
T4is a vertex position vector of P4The brightness temperature value corresponding to the tetrahedron vertex.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111217891.2A CN114034388B (en) | 2021-10-19 | 2021-10-19 | Brightness and temperature drawing method for given place of full moon |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111217891.2A CN114034388B (en) | 2021-10-19 | 2021-10-19 | Brightness and temperature drawing method for given place of full moon |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114034388A true CN114034388A (en) | 2022-02-11 |
| CN114034388B CN114034388B (en) | 2024-01-26 |
Family
ID=80135614
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111217891.2A Active CN114034388B (en) | 2021-10-19 | 2021-10-19 | Brightness and temperature drawing method for given place of full moon |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114034388B (en) |
Citations (4)
| 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 |
| CN111006771A (en) * | 2019-12-27 | 2020-04-14 | 广东电网有限责任公司电力科学研究院 | Method and device for judging and identifying fire point based on polar orbit meteorological satellite |
| CN111898687A (en) * | 2020-08-03 | 2020-11-06 | 成都信息工程大学 | Radar reflectivity data fusion method based on Dilongnie triangulation |
| CN112197865A (en) * | 2020-09-02 | 2021-01-08 | 华中科技大学 | Estimation method and system for observation brightness temperature data error of satellite-borne microwave radiometer |
-
2021
- 2021-10-19 CN CN202111217891.2A patent/CN114034388B/en active Active
Patent Citations (4)
| 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 |
| CN111006771A (en) * | 2019-12-27 | 2020-04-14 | 广东电网有限责任公司电力科学研究院 | Method and device for judging and identifying fire point based on polar orbit meteorological satellite |
| CN111898687A (en) * | 2020-08-03 | 2020-11-06 | 成都信息工程大学 | Radar reflectivity data fusion method based on Dilongnie triangulation |
| CN112197865A (en) * | 2020-09-02 | 2021-01-08 | 华中科技大学 | Estimation method and system for observation brightness temperature data error of satellite-borne microwave radiometer |
Non-Patent Citations (2)
| Title |
|---|
| Z.G. MENG.ET.AL: "Mare Deposits Identification and Feature Analysis in Mare Australe Based on CE‐2 CELMS Data", JGR PLANETS * |
| 连懿.ET.AL: "利用嫦娥二号微波辐射计数据的全月亮温制图", 武汉大学学报.信息科学版, vol. 40, no. 6, pages 732 - 737 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114034388B (en) | 2024-01-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Pol et al. | Future prospects for ground-based gravitational-wave detectors: the galactic double neutron star merger rate revisited | |
| Guerrero et al. | What sets the magnetic field strength and cycle period in solar-type stars? | |
| Greybush et al. | The ensemble Mars atmosphere reanalysis system (EMARS) version 1.0 | |
| Newson | Response of a general circulation model of the atmosphere to removal of the Arctic ice-cap | |
| McIntosh et al. | What the sudden death of solar cycles can tell us about the nature of the solar interior | |
| Ren et al. | Gaia EDR3 parallax zero-point offset based on W Ursae Majoris-type eclipsing binaries | |
| Ment et al. | The occurrence rate of terrestrial planets orbiting nearby mid-to-late M dwarfs from TESS Sectors 1–42 | |
| Kodric et al. | Cepheids in M31: The PAndromeda Cepheid Sample | |
| Montero et al. | Solar radiation and shadow modelling with adaptive triangular meshes | |
| Decraemer et al. | Three-dimensional density structure of a solar coronal streamer observed by SOHO/LASCO and STEREO/COR2 in quadrature | |
| Usha et al. | Effect of aerosol-induced snow darkening on the direct radiative effect of aerosols over the Himalayan region | |
| Krista et al. | The longitudinal evolution of equatorial coronal holes | |
| Kleint et al. | Unusual filaments inside the umbra | |
| Daerden et al. | Explaining NOMAD D/H observations by cloud‐induced fractionation of water vapor on Mars | |
| England et al. | The effect of non-migrating tides on the morphology of the equatorial ionospheric anomaly: seasonal variability | |
| Barker et al. | A New View of the Lunar South Pole from the Lunar Orbiter Laser Altimeter (LOLA) | |
| Godsalve | Bi-directional reflectance sampling by ATSR-2: a combined orbit and scan model | |
| Harris et al. | Tracking movement of long-lived equatorial coronal holes from analysis of long-term McIntosh archive data | |
| CN114034388B (en) | Brightness and temperature drawing method for given place of full moon | |
| Stern et al. | Evidence for possible clouds in Pluto’s present-day atmosphere | |
| Jin et al. | Synoptic mapping of chromospheric magnetic flux | |
| Jones et al. | Improving Coronal Magnetic Field Models Using Image Optimization | |
| McNeill et al. | Comparison of the physical properties of the L4 and L5 Trojan asteroids from ATLAS data | |
| Pick et al. | Homologous solar events on 2011 January 27: Build-up and propagation in a complex coronal environment | |
| Sanromá et al. | Reconstructing the photometric light curves of Earth as a planet along its history |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |