CN109684317A - A kind of data editing method of moon sight GPR Detection Data - Google Patents

A kind of data editing method of moon sight GPR Detection Data Download PDF

Info

Publication number
CN109684317A
CN109684317A CN201811578788.9A CN201811578788A CN109684317A CN 109684317 A CN109684317 A CN 109684317A CN 201811578788 A CN201811578788 A CN 201811578788A CN 109684317 A CN109684317 A CN 109684317A
Authority
CN
China
Prior art keywords
detection data
moon sight
gpr detection
data
gpr
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201811578788.9A
Other languages
Chinese (zh)
Other versions
CN109684317B (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

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 kind of data editing methods of moon sight GPR Detection Data, comprising: step 1, pre-processes to the telemetry of moon sight GPR Detection Data and auxiliary;Step 2, telemetry is inserted into moon sight GPR Detection Data according to the time, the moon sight GPR Detection Data under walking states is then chosen according to the location information in moon sight GPR Detection Data;Step 3, alignment of data and distance scale consistency treatment are carried out to the moon sight GPR Detection Data of selection, and then obtains final processing result.The data editing method of moon sight GPR Detection Data provided by the invention, processing and edit operation are carried out to moon sight GPR Detection Data in conjunction with telemetry, the consistent moon sight GPR Detection Data processing result of distance scale is obtained, is conducive to the later period in conjunction with effective identification of moon sight GPR Detection Data progress target and layer position.

Description

A kind of data editing method of moon sight GPR Detection Data
Technical field
The present invention relates to a kind of data editing methods of moon sight GPR Detection Data, in conjunction with telemetry for realizing moon sight The data edition of GPR Detection Data works.
Background technique
No. three satellites of the goddess in the moon are succeeded in sending up on December 2nd, 2013 in Xichang Launching Site, landing on December 14 moon Ball Mare Imbrium northern territory, rover lunar rover are successfully separated with lander, respectively carry out detection operations.Moon sight radar is as the goddess in the moon The upper important scientific load of No. three " Jade Hare number " lunar rovers, during lunar surface detection, the effective mileage detected in total is greater than 100m。
Moon sight radar is divided into three phases in the work of lunar surface, and the first stage is the detection work in the time in first month daytime Make, lunar rover advances to navigation spots N0108 point from navigation spots N0101, which is mainly the Orbital detection stage, to moon sight radar Several important parameters in two channels for example when window, accumulative frequency, the repetition rate of impulse ejection, receiver signal gain mode and Decaying setting, is tested, and has finally been determined that radar works the parameter of preferable state.Second stage is the time in second month daytime Interior detection operations, lunar rover advance to navigation spots N209 point from navigation spots N0201, and the stage radar groundwork is fixed one Parameter under be unfolded detection.Three phases are third month daytime and later working stage, and moon sight radar is in navigation spots N0209 The fixed point of expansion detects.One channel central frequency of moon sight radar is 60MHz, and corresponding bandwidth is 40MHz-80MHz, moon sight radar two Channel central frequency is 500MHz, and corresponding bandwidth is 250MHz-750MHz.
There are following features in the work of lunar surface for moon sight radar: first, detection of the lunar rover on lunar surface is not to connect Continuous, it loiters, there is a large amount of fixed point repeated data;Second, there is also differences for the speed that different zones are advanced, and will lead to There is a problem of that scale is inconsistent after radar data is spliced;Third, every time when booting, due to space environment and Rover attitude Difference, data starting sample point can have differences, and can have the phenomenon that misalignment after data splicing;4th, road in radar data The location information of head is very little, about updates per minute, can no matter when Primary Stage Data screening or later data processing summarize It is not high with property.For subsequent better interpretation moon sight GPR Detection Data, need combine telemetry to original radar data into Row data edition.
Existing data edition work mostly uses subjective data edit methods at present, by manually checking radar data, leads to Cross data characteristic original GPR Detection Data is deleted, the operation such as polishing, lack effective reference frame.These methods Disadvantage is do not have accuracy, referential, and common processing result solves more.It is desirable to have a kind of effective moon sight thunder Up to the data editing method of detection data.
Summary of the invention
(1) technical problems to be solved
In order to overcome the drawbacks described above of the prior art, the present invention provides a kind of combination telemetries to moon sight radar detection The data editing method of data.
(2) technical solution
The data editing method of a kind of moon sight GPR Detection Data provided by the invention, this method comprises: step 1: to survey The telemetry of month GPR Detection Data and auxiliary is pre-processed;Step 2: telemetry is inserted into moon sight thunder according to the time Up in detection data, the moon sight radar under walking states is then chosen according to the location information in moon sight GPR Detection Data Detection data;Step 3: alignment of data and distance scale consistency treatment are carried out to the moon sight GPR Detection Data of selection.
Preferably, carrying out pretreatment to the telemetry to moon sight GPR Detection Data and auxiliary described in step 1 includes: Time homogenization correction process is carried out to moon sight GPR Detection Data, and interpolation processing is carried out to telemetry.
Preferably, described that time homogenization correction process is carried out to moon sight GPR Detection Data, comprising: according to two neighboring Thunder is calculated with the end time at the beginning of moon sight GPR Detection Data in a detection cycle time in moon sight radar navigation point Up to net cycle time, then using the radar road number obtained in the area, temporal information is averagely inserted into per pass moon sight thunder Up in detection data.
Preferably, described that interpolation processing is carried out to telemetry, comprising: to utilize linear interpolation method, original 4 seconds are updated The telemetry of position is interpolated into the position of update in every 0.05 second.
Preferably, telemetry is inserted into moon sight GPR Detection Data according to the time described in step 2, comprising: benefit With the time of per pass moon sight GPR Detection Data, gone for according to the time near before and after this time point two in telemetry A moment takes the average position coordinates as the road moon sight GPR Detection Data to the two moment corresponding position coordinates.
Preferably, it is chosen under walking states described in step 2 according to the location information in moon sight GPR Detection Data Moon sight GPR Detection Data, comprising: with the location information (x, y, z) of per pass moon sight GPR Detection Data be judgment basis, will The same moon sight GPR Detection Data of all coordinate (x, y, z) values takes average preservation at one of moon sight GPR Detection Data.
Preferably, alignment of data is carried out to the moon sight GPR Detection Data of selection described in step 3, comprising: with first First trough sampled point of moon sight GPR Detection Data is calibration point, by first trough of all moon sight GPR Detection Datas Point is aligned with calibration point, and the valid data points part exceeded is cut, the part zero filling of insufficient valid data points Processing.Wherein, one channel detection data valid data of moon sight radar points are 4096, two channel detection data valid data points It is 2048.
Preferably, the moon sight GPR Detection Data to selection carries out distance scale consistency treatment, comprising: first Step carries out 100 times of linear interpolation to the moon sight GPR Detection Data containing location information;Second step, in newly-generated moon sight thunder Up in detection data, since the first track data, lower one and preceding path length difference together are calculated, and add up the path length difference, when tired When the distance of meter is greater than scale preset value, it is recorded as first of moon sight GPR Detection Data;Third step, terminate from previous step that Road moon sight GPR Detection Data starts, and continues accumulative operation, and so on, traverse entire moon sight GPR Detection Data, so that it may Obtain the moon sight GPR Detection Data that scale preset value is divided between horizontal distance.Wherein, one channel detection data ruler of moon sight radar Degree threshold value is 7.053cm, and two channel detection data scale threshold value of moon sight radar is 0.045cm.
(3) beneficial effect
Using the data editing method of moon sight GPR Detection Data provided by the invention, moon sight radar telemetering number is combined According to, realize radar data it is effective arrangement and editing, to later use radar data carry out Objective extraction, layer position identify And data interpretation is significant.
Detailed description of the invention
Fig. 1 is the data editing method flow chart of moon sight GPR Detection Data according to an embodiment of the present invention;
Fig. 2 be according to an embodiment of the present invention to moon sight GPR Detection Data carry out the time homogenization correction process before and after Comparison diagram;
Fig. 3 is the comparison diagram carried out before and after interpolation to moon sight radar telemetry according to an embodiment of the present invention;
Fig. 4 is comparison of the moon sight GPR Detection Data before and after walking status data is chosen according to an embodiment of the present invention Figure;
Fig. 5 is the comparison diagram carried out before and after alignment of data to moon sight radar data according to an embodiment of the present invention;
Fig. 6 be according to an embodiment of the present invention to moon sight GPR Detection Data carry out distance scale consistency treatment before and after Comparison diagram.
Specific embodiment
To make the objectives, technical solutions, and advantages of the present invention clearer, below in conjunction with specific embodiment, and reference Attached drawing, the present invention is described in more detail.
Fig. 1 is the data editing method flow chart of moon sight GPR Detection Data according to an embodiment of the present invention, this method packet It includes: step 1: the telemetry of moon sight GPR Detection Data and auxiliary being pre-processed, mainly includes to moon sight radar detection Data carry out time homogenization correction process, and carry out interpolation processing to telemetry;Step 2: by telemetry according to the time It is inserted into moon sight GPR Detection Data, is then chosen under walking states according to location information in moon sight GPR Detection Data Moon sight GPR Detection Data;Step 3: alignment of data and distance scale consistency are carried out to the moon sight GPR Detection Data of selection Processing, and then obtain final processing result.
In embodiments of the present invention, time homogenization correction process is carried out to moon sight GPR Detection Data described in step 1, Include: according to moon sight GPR Detection Data in the detection cycle time in two neighboring moon sight radar navigation point at the beginning of Radar net cycle time is calculated with the end time temporal information is averaged then using the radar road number obtained in the area It is inserted into per pass moon sight GPR Detection Data.
Fig. 2 be according to an embodiment of the present invention to moon sight GPR Detection Data carry out the time homogenization correction process before and after Comparison diagram.By taking two channel data of moon sight radar as an example, by sample frequency, cumulative number, theoretical calculation learns adjacent twice number According to time difference should be 0.8192s.In figure, left figure is before correcting, and the time difference that comparison discovery corrects preceding two track data exists Near 0.8192s, and there is fluctuation;Right figure is after correcting, and the time difference of two track datas is unanimously in 0.8192s, correction after correction As a result consistent with calculated value.
In embodiments of the present invention, interpolation processing is carried out to telemetry described in step 1, comprising: utilize linear interpolation The telemetry that original 4 seconds update a position is interpolated into the position of update in every 0.05 second by method.
Fig. 3 is the comparison diagram carried out before and after interpolation to moon sight radar telemetry according to an embodiment of the present invention.Because y and Z coordinate is consistent with x coordinate trend, and the change curve of x coordinate is only listed in figure.The original position-information of telemetry such as left figure institute Show, update once within every 4 seconds, it cannot be guaranteed that (data break of moon sight radar telemetry is 1s to per pass moon sight radar telemetry Left and right) a corresponding position, therefore telemetry is interpolated into every 0.05s location updating is primary, the result such as right figure after interpolation It is shown.
In embodiments of the present invention, telemetry is inserted into moon sight GPR Detection Data according to the time described in step 2 In, comprising: using the time of per pass moon sight GPR Detection Data, gone for according to the time near this time in telemetry At former and later two moment of point, the average position as the road moon sight GPR Detection Data is taken to the two moment corresponding position coordinates Coordinate.
In embodiments of the present invention, it is in described in step 2 according to the location information selection in moon sight GPR Detection Data Moon sight GPR Detection Data under walking states, comprising: with the location information (x, y, z) of per pass moon sight GPR Detection Data be sentence The same moon sight GPR Detection Data of all coordinates (x, y, z) value is taken average preservation at one of moon sight radar detection by disconnected foundation Data.
Fig. 4 is comparison of the moon sight GPR Detection Data before and after walking status data is chosen according to an embodiment of the present invention Figure.Before choosing data, one shares 19934 original moon sight GPR Detection Datas in left figure;By the moon sight radar under walking states Detection data is as shown on the right after choosing, and rejects after repeating moon sight GPR Detection Data, the moon sight GPR Detection Data of walking has 2699.
In embodiments of the present invention, alignment of data, packet are carried out to the moon sight GPR Detection Data of selection described in step 3 It includes: using first trough sampled point of first of moon sight GPR Detection Data as calibration point, by all moon sight GPR Detection Datas First trough point be aligned with calibration point, the valid data that exceed points part is cut, insufficient significant figure strong point Several part zero filling processing.
Fig. 5 is the comparison diagram carried out before and after alignment of data to moon sight radar data according to an embodiment of the present invention.Obviously may be used To find out, compared with the left figure before being aligned, the right figure after alignment obviously more meets practical situation, is embodied at surface echo In same horizontal line.
In embodiments of the present invention, consistent to the moon sight GPR Detection Data progress distance scale of selection described in step 3 Property processing, comprising: the first step carries out 100 times of linear interpolation to the moon sight GPR Detection Data containing location information;Second step, In newly-generated moon sight GPR Detection Data, since the first track data, lower one and preceding path length difference together are calculated, and tire out The path length difference is counted, when accumulative distance is greater than scale preset value, is recorded as first of moon sight GPR Detection Data;Third step, That road moon sight GPR Detection Data terminated since previous step, continue accumulative operation, and so on, traverse entire moon sight radar Detection data, so that it may obtain the moon sight GPR Detection Data that scale preset value is divided between horizontal distance.
Fig. 6 be according to an embodiment of the present invention to moon sight GPR Detection Data carry out distance scale consistency treatment before and after Comparison diagram.In figure, the band radar result of arrow mark is stretched more serious in left figure, is the speed advanced due to corresponding lunar rover More slowly cause;And in the image after right figure correction, by distance scale consistency treatment, arrow mark region has to be changed greatly very much It is kind.
In embodiments of the present invention, the telemetry of moon sight GPR Detection Data and auxiliary is pre-processed first, it is main Including carrying out time homogenization, telemetry interpolation to moon sight GPR Detection Data;Secondly by telemetry according to when interleave Enter into moon sight GPR Detection Data, is then chosen under walking states according to location information in moon sight GPR Detection Data Moon sight GPR Detection Data;Finally the moon sight GPR Detection Data of selection is carried out at alignment of data and distance scale consistency Reason, and then obtain final processing result;The work carries out Objective extraction to later use moon sight GPR Detection Data, layer position is known Not and data interpretation is significant.
Particular embodiments described above has carried out further in detail the purpose of the present invention, technical scheme and beneficial effects It describes in detail bright, it should be understood that the above is only a specific embodiment of the present invention, is not intended to restrict the invention, it is all Within the spirit and principles in the present invention, any modification, equivalent substitution, improvement and etc. done should be included in guarantor of the invention Within the scope of shield.

Claims (10)

1. a kind of data editing method of moon sight GPR Detection Data, which is characterized in that this method comprises:
Step 1: the telemetry of moon sight GPR Detection Data and auxiliary is pre-processed;
Step 2: telemetry being inserted into moon sight GPR Detection Data according to the time, then according to moon sight GPR Detection Data In location information choose the moon sight GPR Detection Data under the walking states;
Step 3: alignment of data and distance scale consistency treatment are carried out to the moon sight GPR Detection Data of selection.
2. the data editing method of moon sight GPR Detection Data according to claim 1, which is characterized in that institute in step 1 It states that the telemetry of moon sight GPR Detection Data and auxiliary pre-process and includes:
Time homogenization correction process is carried out to moon sight GPR Detection Data, and interpolation processing is carried out to telemetry.
3. the data editing method of moon sight GPR Detection Data according to claim 2, which is characterized in that described to moon sight GPR Detection Data carries out time homogenization correction process, comprising:
At the beginning of moon sight GPR Detection Data in the detection cycle time in two neighboring moon sight radar navigation point Radar net cycle time is calculated with the end time temporal information is averaged then using the radar road number obtained in the area It is inserted into per pass moon sight GPR Detection Data.
4. the data editing method of moon sight GPR Detection Data according to claim 2, which is characterized in that described to telemetering Data carry out interpolation processing, comprising:
Using linear interpolation method, the telemetry that original 4 seconds update a position is interpolated into every 0.05 second and updates primary position It sets.
5. the data editing method of moon sight GPR Detection Data according to claim 1, which is characterized in that institute in step 2 It states and is inserted into telemetry in moon sight GPR Detection Data according to the time, comprising:
Using the time of per pass moon sight GPR Detection Data, gone for according to the time near this time point in telemetry Former and later two moment take the average position as the road moon sight GPR Detection Data to sit on the two moment corresponding position coordinates Mark.
6. the data editing method of moon sight GPR Detection Data according to claim 1, which is characterized in that institute in step 2 State the moon sight GPR Detection Data chosen under walking states according to the location information in moon sight GPR Detection Data, comprising:
It is with the location information (x, y, z) of per pass moon sight GPR Detection Data for judgment basis, all coordinates (x, y, z) value is the same Moon sight GPR Detection Data, take it is average preservation at one of moon sight GPR Detection Data.
7. the data editing method of moon sight GPR Detection Data according to claim 1, which is characterized in that institute in step 3 It states and alignment of data is carried out to the moon sight GPR Detection Data of selection, comprising:
Using first trough sampled point of first of moon sight GPR Detection Data as calibration point, by all moon sight GPR Detection Datas First trough point be aligned with calibration point, the valid data that exceed points part is cut, insufficient significant figure strong point Several part zero filling processing.
8. the data editing method of moon sight GPR Detection Data according to claim 7, which is characterized in that moon sight radar one Channel detection data valid data points are 4096, and two channel detection data valid data points are 2048.
9. the data editing method of moon sight GPR Detection Data according to claim 1, which is characterized in that institute in step 3 It states and distance scale consistency treatment is carried out to the moon sight GPR Detection Data of selection, comprising:
The first step carries out 100 times of linear interpolation to the moon sight GPR Detection Data containing location information;
Second step, in newly-generated moon sight GPR Detection Data, since the first track data, calculate it is lower together with it is preceding together Path length difference, and add up the path length difference, when accumulative distance is greater than scale preset value, it is recorded as first of moon sight radar detection number According to;
Third step that road moon sight GPR Detection Data terminated since previous step, continues accumulative operation, and so on, traversal Entire moon sight GPR Detection Data obtains the moon sight GPR Detection Data that scale preset value is divided between horizontal distance.
10. the data editing method of moon sight GPR Detection Data according to claim 9, which is characterized in that moon sight radar One channel detection data scale threshold value is 7.053cm, and two channel detection data scale threshold value of moon sight radar is 0.045cm.
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 true CN109684317A (en) 2019-04-26
CN109684317B 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 (4)

* 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
CN105044696A (en) * 2015-07-30 2015-11-11 中国科学院国家天文台 Method for calculating probing depth of moon-probing radar based on correlation coefficient analysis method
WO2016108904A1 (en) * 2014-12-31 2016-07-07 Halliburton Energy Services, Inc. Electromagnetic telemetry for sensor systems deployed in a borehole environment

Patent Citations (4)

* 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
CN105044696A (en) * 2015-07-30 2015-11-11 中国科学院国家天文台 Method for calculating probing depth of moon-probing radar based on correlation coefficient analysis method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
李春来 等: "《月球数据预处理工作流模型的构建及应用》", 《吉林大学学报(工学版)》 *
李臣 等: "《美军某卫星系统抗干扰性能仿真分析》", 《中国通信学会青年工作委员会会议论文集》 *
苏彦 等: "《雷达探测技术在探月中的应用》", 《天文研究与技术》 *
高云泽 等: "《嫦娥三号测月雷达第一通道数据处理与分析》", 《雷达学报》 *

Also Published As

Publication number Publication date
CN109684317B (en) 2020-10-09

Similar Documents

Publication Publication Date Title
CN111833333B (en) Method and system for measuring pose of cantilever tunneling equipment based on binocular vision
US10634500B2 (en) Aircraft and obstacle avoidance method and system thereof
CN108868805B (en) Shield method for correcting error based on statistical analysis in conjunction with XGboost
CN108931795B (en) Positioning equipment track optimization and boundary extraction method and device
CN104049636B (en) Navigation altitude obtaining method combining relative altitude and absolute altitude
Sagnotti Paleomagnetic evidence for a Pleistocene counterclockwise rotation of the Sant'Arcangelo basin, southern Italy
CN103079059A (en) Single-gun-camera-multi-dome-camera linkage method based on grid combination interpolation
CN104537671A (en) Cigarette filter online counting and quality detecting method based on machine vision
CN110197318A (en) Surface mine Grade Model update method
DE60217218T2 (en) Method and device for GPS, navigation system and computer program
CN116224382A (en) Post-processing method for weakening GPS (Global positioning System) related LOD (LOD) system errors
CN107966729B (en) A kind of three-dimensional TTI medium ray-tracing procedure and system
CN109684317A (en) A kind of data editing method of moon sight GPR Detection Data
CN109655890A (en) A kind of shallow mid-deep strata joint chromatography inversion speed modeling method of Depth Domain and system
DE102018008402A1 (en) METHOD AND SYSTEM FOR DETERMINING A MOTION DIRECTION OF AN OBJECT
EP1408344B1 (en) Survey apparatus and survey method using a laserscanner
CN113158427B (en) Shield tunneling machine soil body relative displacement calculation method based on guide data
CN113177270B (en) Real-time calculation method for stratum deformation caused by shield attitude
CN105592294A (en) VSP excited cannon group monitoring system
CN106354333B (en) A method of improving touch-control system tracking precision
CN109725354A (en) Anisotropic velocity modeling method and system
CN105044789B (en) A kind of the tunnel effect removing method and device of Mine transient electromagnetic forward probe
CN110568499B (en) First-arrival time difference correction method and device for VSP seismic data
CN109143327B (en) Q value initial modeling method and system based on block constraint
CN111060970B (en) Reservoir constrained deconvolution method and application

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