CN109684317B - Data editing method for lunar sounding radar detection data - Google Patents

Data editing method for lunar sounding radar detection data Download PDF

Info

Publication number
CN109684317B
CN109684317B CN201811578788.9A CN201811578788A CN109684317B CN 109684317 B CN109684317 B CN 109684317B CN 201811578788 A CN201811578788 A CN 201811578788A CN 109684317 B CN109684317 B CN 109684317B
Authority
CN
China
Prior art keywords
lunar
data
detection data
radar detection
sounding
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811578788.9A
Other languages
Chinese (zh)
Other versions
CN109684317A (en
Inventor
苏彦
李春来
刘建军
张洪波
邢树果
薛喜平
戴舜
孔德庆
朱新颖
李俊铎
李臣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National Astronomical Observatories of CAS
Original Assignee
National Astronomical Observatories of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by National Astronomical Observatories of CAS filed Critical National Astronomical Observatories of CAS
Priority to CN201811578788.9A priority Critical patent/CN109684317B/en
Publication of CN109684317A publication Critical patent/CN109684317A/en
Application granted granted Critical
Publication of CN109684317B publication Critical patent/CN109684317B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

The invention discloses a data editing method of lunar sounding radar detection data, which comprises the following steps: step 1, preprocessing lunar radar detection data and auxiliary telemetering data; step 2, inserting the telemetering data into the lunar sounding radar detection data according to time, and then selecting the lunar sounding radar detection data in a walking state according to position information in the lunar sounding radar detection data; and 3, carrying out data alignment and distance scale consistency processing on the selected lunar sounding radar detection data to obtain a final processing result. The data editing method of the lunar sounding radar detection data provided by the invention is used for processing and editing the lunar sounding radar detection data by combining the telemetering data to obtain the lunar sounding radar detection data processing result with consistent distance scale, and is beneficial to effectively identifying the target and the horizon by combining the lunar sounding radar detection data at the later stage.

Description

Data editing method for lunar sounding radar detection data
Technical Field
The invention relates to a data editing method of lunar sounding radar detection data, which is used for realizing data editing work of the lunar sounding radar detection data by combining telemetering data.
Background
The Chang' e No. three satellite is successfully launched in the West Chang satellite launching center in 2013, 12 months and 2 days, lands in the northern region of the lunar sea in 12 months and 14 days, and the patrol instrument lunar vehicle is successfully separated from the lander, so that detection work is respectively carried out. The lunar ranging radar is used as an important scientific load on the Chang 'e' III 'rabbit' lunar rover, and the total effective detection mileage is larger than 100m during lunar surface detection.
The lunar radar works on the lunar surface in three stages, the first stage is detection work in the first lunar day time, a lunar vehicle travels from a navigation point N0101 to a navigation point N0108, the lunar radar is mainly in an on-orbit testing stage, a test is carried out on several important parameters of two channels of the lunar radar, such as a time window, accumulation times, repetition frequency of pulse emission, a receiver signal gain mode and attenuation setting, and finally parameters of a better working state of the radar are determined. The second phase is a detection operation in the second month day, the lunar vehicle travels from the navigation point N0201 to the navigation point N209, and the radar basically operates under a fixed parameter to perform detection. The third stage is the working stage of the third month day and later, and the lunar radar carries out fixed point detection at the navigation point N0209. The center frequency of the first channel of the lunar ranging radar is 60MHz, the corresponding bandwidth is 40MHz-80MHz, the center frequency of the second channel of the lunar ranging radar is 500MHz, and the corresponding bandwidth is 250MHz-750 MHz.
The lunar radar has the following characteristics in the working of the lunar surface: firstly, the lunar rover is not continuously detected on the lunar surface, and is stopped when the lunar rover moves, so that a large amount of fixed-point repeated data exist; secondly, the traveling speeds of different areas are different, so that the problem of inconsistent scale exists after the radar data are spliced; thirdly, when the lunar rover is started up every time, due to the fact that the space environment and the lunar rover are different in posture, the initial sampling points of the data are different, and the phenomenon of misalignment exists after the data are spliced; fourthly, the position information of the track head in the radar data is too little, and is updated about every minute, so that the usability is not high no matter in early data screening or later data processing and summarizing. In order to better interpret the lunar radar detection data later, data editing needs to be carried out on the original radar data by combining the telemetering data.
At present, the existing data editing work mostly adopts a subjective data editing method, radar data is checked manually, and operations such as deleting, supplementing and the like are carried out on original radar detection data through data characteristics, so that an effective reference basis is lacked. The methods have the disadvantages of no accuracy and reference, and the common processing result is multi-solution. Therefore, an effective data editing method for lunar sounding radar sounding data is expected.
Disclosure of Invention
Technical problem to be solved
In order to overcome the defects in the prior art, the invention provides a data editing method for lunar radar detection data by combining telemetering data.
(II) technical scheme
The invention provides a data editing method of lunar sounding radar detection data, which comprises the following steps: step 1: preprocessing the lunar radar detection data and the auxiliary telemetering data; step 2: inserting the telemetering data into the lunar sounding radar detection data according to time, and then selecting the lunar sounding radar detection data in a walking state according to position information in the lunar sounding radar detection data; and step 3: and carrying out data alignment and distance scale consistency processing on the selected lunar sounding radar detection data.
Preferably, the preprocessing of the lunar radar detection data and the auxiliary telemetry data in step 1 comprises: and carrying out time homogenization correction processing on the lunar radar detection data, and carrying out interpolation processing on the telemetering data.
Preferably, the time-uniformization correction processing of the lunar sounding radar detection data includes: and calculating the total working time of the radar according to the starting time and the ending time of the detection data of the lunar sounding radar in one detection period time in two adjacent lunar sounding radar navigation points, and then averagely inserting the time information into the detection data of the lunar sounding radar in each channel by using the number of radar channels acquired in the area.
Preferably, the interpolation processing of the telemetry data includes: the telemetry data for the original 4 second update position is interpolated to an update position every 0.05 seconds using a linear interpolation method.
Preferably, the step 2 of inserting the telemetry data into the lunar sounding radar detection data according to time comprises: and finding two moments before and after the closest time point in the telemetering data according to the time by using the time of the detection data of each lunar ranging radar, and averaging the corresponding position coordinates of the two moments to obtain the position coordinates of the detection data of the lunar ranging radar.
Preferably, the selecting the lunar sounding radar detection data in the walking state according to the position information in the lunar sounding radar detection data in the step 2 includes: and averaging and storing the lunar sounding radar detection data with the same coordinate (x, y, z) value into a lunar sounding radar detection data by taking the position information (x, y, z) of the lunar sounding radar detection data of each channel as a judgment basis.
Preferably, the data alignment of the selected lunar sounding radar detection data in step 3 includes: and taking the first trough sampling point of the first lunar rover detection data as a calibration point, aligning the first trough point of all the lunar rover detection data with the calibration point, cutting the excessive effective data points, and filling zero in the insufficient effective data points. The number of effective data points of the detection data of the first channel of the lunar sounding radar is 4096, and the number of effective data points of the detection data of the second channel is 2048.
Preferably, the distance scale consistency processing on the selected lunar sounding radar detection data includes: firstly, performing 100-time linear interpolation on lunar ranging radar detection data containing position information; secondly, calculating the path difference between the next path and the previous path from the first path in the newly generated lunar ranging radar detection data, accumulating the path difference, and recording the path difference as the first lunar ranging radar detection data when the accumulated path is larger than a preset scale value; and thirdly, starting from the last month measuring radar detection data, continuing accumulation operation, and repeating the operation in the same way, traversing the whole month measuring radar detection data, and thus obtaining the month measuring radar detection data with the horizontal distance interval as the scale preset value. The scale threshold value of the detection data of the first channel of the lunar rover is 7.053cm, and the scale threshold value of the detection data of the second channel of the lunar rover is 0.045 cm.
(III) advantageous effects
The data editing method of the lunar sounding radar detection data provided by the invention is combined with the lunar sounding radar remote measurement data, realizes effective arrangement and editing work of radar data, and has important significance for carrying out target extraction, horizon identification and data interpretation by utilizing the radar data subsequently.
Drawings
FIG. 1 is a flow chart of a data editing method for lunar sounding radar detection data according to an embodiment of the present invention;
FIG. 2 is a comparison of lunar radar survey data before and after time normalization correction processing according to an embodiment of the present invention;
FIG. 3 is a comparison of lunar radar telemetry data before and after interpolation, in accordance with an embodiment of the present invention;
FIG. 4 is a comparison of lunar sounding radar data before and after travel state data selection according to an embodiment of the present invention;
FIG. 5 is a comparison of lunar radar data before and after data alignment according to an embodiment of the present invention;
fig. 6 is a comparison diagram before and after distance scale consistency processing is performed on lunar radar detection data according to an embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to specific embodiments and the accompanying drawings.
Fig. 1 is a flowchart of a data editing method of lunar sounding radar detection data according to an embodiment of the present invention, the method including: step 1: preprocessing the lunar radar detection data and the auxiliary telemetering data, wherein the preprocessing mainly comprises the steps of carrying out time homogenization correction processing on the lunar radar detection data and carrying out interpolation processing on the telemetering data; step 2: inserting the telemetering data into the lunar sounding radar detection data according to time, and then selecting the lunar sounding radar detection data in a walking state according to position information in the lunar sounding radar detection data; and step 3: and carrying out data alignment and distance scale consistency processing on the selected lunar sounding radar detection data so as to obtain a final processing result.
In the embodiment of the present invention, the performing time homogenization correction processing on the lunar radar detection data in step 1 includes: and calculating the total working time of the radar according to the starting time and the ending time of the detection data of the lunar sounding radar in one detection period time in two adjacent lunar sounding radar navigation points, and then averagely inserting the time information into the detection data of the lunar sounding radar in each channel by using the number of radar channels acquired in the area.
Fig. 2 is a comparison of lunar radar detection data before and after time-normalized correction processing according to an embodiment of the present invention. Taking the data of the two channels of the lunar ranging radar as an example, the time difference between the two adjacent channels of data is obtained through theoretical calculation through sampling frequency and accumulated times and is 0.8192 s. In the figure, the left figure is before correction, the time difference of two data before correction is found by comparison to be about 0.8192s, and fluctuation exists; the right graph shows that after correction, the time difference of the two corrected data is consistent at 0.8192s, and the correction result is consistent with the theoretical calculation value.
In this embodiment of the present invention, the interpolating the telemetry data in step 1 includes: the telemetry data for the original 4 second update position is interpolated to an update position every 0.05 seconds using a linear interpolation method.
FIG. 3 is a comparison of lunar radar telemetry data before and after interpolation, in accordance with an embodiment of the present invention. Since the y and z coordinates are in line with the x coordinate trend, only the x coordinate variation curve is listed in the figure. The original position information of the telemetering data is shown in the left graph, the telemetering data is updated every 4 seconds, and it cannot be guaranteed that the telemetering data of the lunar sounding radar (the data interval of the telemetering data of the lunar sounding radar is about 1 s) corresponds to one position, so that the telemetering data is interpolated to be updated every 0.05s, and the interpolated result is shown in the right graph.
In the embodiment of the present invention, the inserting the telemetry data into the lunar sounding radar detection data according to time in step 2 includes: and finding two moments before and after the closest time point in the telemetering data according to the time by using the time of the detection data of each lunar ranging radar, and averaging the corresponding position coordinates of the two moments to obtain the position coordinates of the detection data of the lunar ranging radar.
In the embodiment of the present invention, the selecting the lunar sounding radar detection data in the walking state according to the position information in the lunar sounding radar detection data in step 2 includes: and averaging and storing the lunar sounding radar detection data with the same coordinate (x, y, z) value into a lunar sounding radar detection data by taking the position information (x, y, z) of the lunar sounding radar detection data of each channel as a judgment basis.
Fig. 4 is a comparison graph of lunar sounding radar data before and after walking state data selection according to an embodiment of the invention. Before selecting data, 19934 original lunar sounding radar detection data are shared in the left image; after the detection data of the lunar ranging radar in the walking state are selected and shown in the right graph, 2699 tracks of detection data of the lunar ranging radar are walked after repeated detection data of the lunar ranging radar are removed.
In the embodiment of the present invention, the performing data alignment on the selected lunar sounding radar detection data in step 3 includes: and taking the first trough sampling point of the first lunar rover detection data as a calibration point, aligning the first trough point of all the lunar rover detection data with the calibration point, cutting the excessive effective data points, and filling zero in the insufficient effective data points.
Fig. 5 is a comparison of before and after data alignment of lunar sounding radar data in accordance with an embodiment of the present invention. It is clear that the right image after alignment is significantly more realistic than the left image before alignment, in particular when the surface echoes are at the same horizontal line.
In the embodiment of the present invention, the distance scale consistency processing on the selected lunar radar detection data in step 3 includes: firstly, performing 100-time linear interpolation on lunar ranging radar detection data containing position information; secondly, calculating the path difference between the next path and the previous path from the first path in the newly generated lunar ranging radar detection data, accumulating the path difference, and recording the path difference as the first lunar ranging radar detection data when the accumulated path is larger than a preset scale value; and thirdly, starting from the last month measuring radar detection data, continuing accumulation operation, and repeating the operation in the same way, traversing the whole month measuring radar detection data, and thus obtaining the month measuring radar detection data with the horizontal distance interval as the scale preset value.
Fig. 6 is a comparison diagram before and after distance scale consistency processing is performed on lunar radar detection data according to an embodiment of the invention. In the figure, the radar result of the area marked by the arrow on the left figure is more stretched due to the slower speed of the corresponding lunar vehicle; in the image after the right image correction, the arrow mark area is greatly improved through distance scale consistency processing.
In the embodiment of the invention, firstly, preprocessing the lunar sounding radar detection data and the auxiliary telemetering data, and mainly comprises the steps of carrying out time homogenization and telemetering data interpolation on the lunar sounding radar detection data; secondly, inserting the telemetering data into the lunar sounding radar detection data according to time, and then selecting the lunar sounding radar detection data in a walking state according to position information in the lunar sounding radar detection data; finally, carrying out data alignment and distance scale consistency processing on the selected lunar sounding radar detection data to obtain a final processing result; the work has important significance for carrying out target extraction, horizon identification and data interpretation by utilizing the lunar radar detection data subsequently.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the present invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (9)

1. A data editing method for lunar sounding radar detection data is characterized by comprising the following steps:
step 1: preprocessing the lunar radar detection data and the auxiliary telemetering data;
step 2: inserting the telemetering data into the lunar sounding radar detection data according to time, and then selecting the lunar sounding radar detection data in a walking state according to position information in the lunar sounding radar detection data;
and step 3: carrying out data alignment and distance scale consistency processing on the selected lunar sounding radar detection data;
wherein, the step 3 of performing distance scale consistency processing on the selected lunar sounding radar detection data comprises the following steps:
firstly, performing 100-time linear interpolation on lunar ranging radar detection data containing position information;
secondly, calculating the path difference between the next path and the previous path from the first path in the newly generated lunar ranging radar detection data, accumulating the path difference, and recording the path difference as the first lunar ranging radar detection data when the accumulated path is larger than a preset scale value;
and thirdly, starting from the last month measuring radar detection data, continuing accumulation operation, and repeating the operation in the same way, traversing the whole month measuring radar detection data, and obtaining the month measuring radar detection data with the horizontal distance interval as the scale preset value.
2. The method for editing data of lunar radar detection data as claimed in claim 1, wherein the preprocessing of the lunar radar detection data and the auxiliary telemetry data in step 1 comprises:
and carrying out time homogenization correction processing on the lunar radar detection data, and carrying out interpolation processing on the telemetering data.
3. The method for editing data of lunar sounding radar detection data according to claim 2, wherein the performing time-normalization correction processing on the lunar sounding radar detection data includes:
and calculating the total working time of the radar according to the starting time and the ending time of the detection data of the lunar sounding radar in a detection period time in two adjacent lunar sounding radar navigation points, and then averagely inserting the time information into the detection data of the lunar sounding radar in each channel by using the number of radar channels acquired in the detection period time.
4. The method for editing data of sounding data of lunar sounding radar according to claim 2, wherein the interpolating the telemetry data comprises:
the telemetry data for the original 4 second update position is interpolated to an update position every 0.05 seconds using a linear interpolation method.
5. The method for editing data of lunar sounding radar detection data according to claim 1, wherein the inserting telemetry data into the lunar sounding radar detection data according to time in step 2 comprises:
and finding two moments before and after the closest time point in the telemetering data according to the time by using the time of the detection data of each lunar ranging radar, and averaging the corresponding position coordinates of the two moments to obtain the position coordinates of the detection data of the lunar ranging radar.
6. The method for editing data of lunar ranging radar detection data according to claim 1, wherein the step 2 of selecting lunar ranging radar detection data in a walking state according to position information in the lunar ranging radar detection data comprises:
and averaging and storing the lunar sounding radar detection data with the same coordinate (x, y, z) value into a lunar sounding radar detection data by taking the position information (x, y, z) of the lunar sounding radar detection data of each channel as a judgment basis.
7. The method for editing data of lunar sounding radar detection data according to claim 1, wherein the step 3 of performing data alignment on the selected lunar sounding radar detection data includes:
and taking the first trough sampling point of the first lunar rover detection data as a calibration point, aligning the first trough point of all the lunar rover detection data with the calibration point, cutting the part exceeding the number of the effective data points, and filling zero in the part lacking the number of the effective data points.
8. The method for editing data of lunar sounding radar detection data as claimed in claim 7, wherein the number of valid data points of lunar sounding radar for one channel detection data is 4096 and the number of valid data points of two channel detection data is 2048.
9. The method for editing data of lunar radar survey data as claimed in claim 1, wherein the preset value of the first channel scale of the lunar radar is 7.053cm, and the preset value of the second channel scale of the lunar radar is 0.045 cm.
CN201811578788.9A 2018-12-21 2018-12-21 Data editing method for lunar sounding radar detection data Active CN109684317B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811578788.9A CN109684317B (en) 2018-12-21 2018-12-21 Data editing method for lunar sounding radar detection data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811578788.9A CN109684317B (en) 2018-12-21 2018-12-21 Data editing method for lunar sounding radar detection data

Publications (2)

Publication Number Publication Date
CN109684317A CN109684317A (en) 2019-04-26
CN109684317B true CN109684317B (en) 2020-10-09

Family

ID=66188922

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811578788.9A Active CN109684317B (en) 2018-12-21 2018-12-21 Data editing method for lunar sounding radar detection data

Country Status (1)

Country Link
CN (1) CN109684317B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412285A (en) * 2013-08-30 2013-11-27 中国科学院国家天文台 Lunar surface inspection tour device moon detection radar data preprocessing method
CN104636608A (en) * 2015-01-30 2015-05-20 国家电网公司 Direct assimilation method of MODIS satellite data
WO2016108904A1 (en) * 2014-12-31 2016-07-07 Halliburton Energy Services, Inc. Electromagnetic telemetry for sensor systems deployed in a borehole environment

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105044696B (en) * 2015-07-30 2017-07-04 中国科学院国家天文台 A kind of method that moon sight radar detection depth is calculated based on correlation analysis method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412285A (en) * 2013-08-30 2013-11-27 中国科学院国家天文台 Lunar surface inspection tour device moon detection radar data preprocessing method
WO2016108904A1 (en) * 2014-12-31 2016-07-07 Halliburton Energy Services, Inc. Electromagnetic telemetry for sensor systems deployed in a borehole environment
CN104636608A (en) * 2015-01-30 2015-05-20 国家电网公司 Direct assimilation method of MODIS satellite data

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
《嫦娥三号测月雷达第一通道数据处理与分析》;高云泽 等;《雷达学报》;20151031;第4卷(第5期);第518-526页 *
《月球数据预处理工作流模型的构建及应用》;李春来 等;《吉林大学学报(工学版)》;20151130;第45卷(第6期);第2007-2013页 *
《美军某卫星系统抗干扰性能仿真分析》;李臣 等;《中国通信学会青年工作委员会会议论文集》;20070831;第1416-1423页 *
《雷达探测技术在探月中的应用》;苏彦 等;《天文研究与技术》;20150430(第2期);第228-242页 *

Also Published As

Publication number Publication date
CN109684317A (en) 2019-04-26

Similar Documents

Publication Publication Date Title
US10215868B2 (en) Automated horizon auto-picking on multiple volumes
CN106990401B (en) full-waveform airborne laser radar data-based two-class elevation error correction method
CN110031898B (en) Data optimization method and integral method prestack depth migration method
CN112180433B (en) Method and device for picking up first arrival wave of earthquake
CN101551467A (en) Automatic first break picking method based on edge detection
US6980482B2 (en) Seismic survey system
CN109684317B (en) Data editing method for lunar sounding radar detection data
CN113075732B (en) Method for eliminating high-resolution small multi-channel seismic stratum abnormal fluctuation
CN112305595B (en) Method for analyzing geologic body structure based on refraction wave and storage medium
CN110568499B (en) First-arrival time difference correction method and device for VSP seismic data
CN114545495A (en) OVT trace gather processing method and device for PS wave seismic trace data and electronic equipment
CN110824559B (en) Method and device for generating co-coordinate point gather and computer-readable storage medium
CN111257969A (en) High-precision speed modeling method under fault control and processing terminal
CN108260359A (en) A kind of seismic data processing methods, devices and systems
CN107356970B (en) High-precision well seismic data matching method
CN112379412A (en) Quality monitoring method and device for collecting seismic data
CN112014816A (en) Double-pass travel time calculation method based on improved horizon tracking algorithm
CN111239822A (en) Fault lower speed modeling method for eliminating well seismic closure difference and processing terminal
CN112904422B (en) Method and system for picking up first arrival time jump point of explosive source seismic data
DE102019215440A1 (en) Recognition of traffic signs
CN111198400B (en) Speed modeling method capable of eliminating well seismic closure difference and processing terminal
CN114859414B (en) Method and device for automatically extracting stratum dip angle information from seismic data
US2858069A (en) Seismic cross-section plotting calculator
CN109752760B (en) Method and device for determining reliability of seismic attribute
Lines et al. Hybrid GPR layer picking method using average square difference function

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