CN113960555B - Target terahertz time-domain echo processing method, device, equipment and storage medium - Google Patents
Target terahertz time-domain echo processing method, device, equipment and storage medium Download PDFInfo
- Publication number
- CN113960555B CN113960555B CN202111224460.9A CN202111224460A CN113960555B CN 113960555 B CN113960555 B CN 113960555B CN 202111224460 A CN202111224460 A CN 202111224460A CN 113960555 B CN113960555 B CN 113960555B
- Authority
- CN
- China
- Prior art keywords
- pulse peak
- echo
- processed
- determining
- sub
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 35
- 239000012634 fragment Substances 0.000 claims description 34
- 238000012545 processing Methods 0.000 claims description 26
- 238000001228 spectrum Methods 0.000 claims description 23
- 238000010606 normalization Methods 0.000 claims description 15
- 238000004458 analytical method Methods 0.000 claims description 13
- 238000004590 computer program Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 abstract description 14
- 238000005259 measurement Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 6
- 230000005855 radiation Effects 0.000 description 6
- 238000005070 sampling Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000002592 echocardiography Methods 0.000 description 3
- 230000036544 posture Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/411—Identification of targets based on measurements of radar reflectivity
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention provides a target terahertz time-domain echo processing method, device, equipment and storage medium, wherein the method comprises the following steps: acquiring a reference echo of a calibration body and acquiring a target echo to be processed; determining a main pulse peak in the reference echo; identifying whether each pulse peak contained in the target echo is a real pulse peak or not according to the main pulse peak; and obtaining the processed target echo according to the identification result. According to the scheme, whether each pulse peak in the target echo is a real pulse peak or not can be automatically identified by utilizing the main pulse peak in the reference echo, and the identification process is an automatic identification process without human participation, so that the identification efficiency and the accuracy can be improved.
Description
Technical Field
The embodiment of the invention relates to the technical field of signal processing, in particular to a target terahertz time-domain echo processing method, device and equipment and a storage medium.
Background
In terahertz time-domain spectrum target characteristic measurement, the measured echo data are utilized to obtain a target echo in the time domain, and the real echo pulse of each scattering center of the target is reflected through the target echo. However, due to the influence of noise and other factors, the actual pulse of the target echo in the echo data is prone to change in value and shape, for example, a plurality of pulse peaks often exist in the actual target echo, and a signal with a magnitude far greater than that of the noise is carried behind the pulse peak included in the target echo obtained by using the echo data, and the signal is similar to the pulse peak shape of the target echo. Therefore, there may be a case where these signals are misjudged as true pulse peaks.
Currently, it is manually identified whether these pulse peaks are true pulse peaks or not according to the shape and amplitude of these pulse peaks. The manual identification mode is poor in accuracy and efficiency.
Disclosure of Invention
Based on the problems of poor accuracy and low efficiency of an artificial identification mode, the embodiment of the invention provides a target terahertz echo processing method, device, equipment and storage medium, which can quickly and accurately identify a real pulse peak in a target echo.
In a first aspect, an embodiment of the present invention provides a target terahertz time-domain echo processing method, including:
Acquiring a reference echo of a calibration body and acquiring a target echo to be processed;
Determining a main pulse peak in the reference echo;
Identifying whether each pulse peak contained in the target echo is a real pulse peak or not according to the main pulse peak;
And obtaining the processed target echo according to the identification result.
Preferably, the identifying whether each pulse peak included in the target echo is a true pulse peak according to the main pulse peak includes:
Determining at least one pulse peak to be processed according to each pulse peak included in the target echo;
For each pulse peak to be processed, determining the pulse peak to be processed and an adjacent pulse peak positioned before the pulse peak to be processed as a first sub-echo;
determining a second sub-echo from the reference echo, wherein the starting point of the second sub-echo is the main pulse peak, and the length of the second sub-echo is equal to that of the first sub-echo;
And judging whether the pulse peak to be processed is a real pulse peak or not according to the amplitude of the second sub-echo and the amplitude of the pulse peak to be processed.
Preferably, determining at least one pulse peak to be processed according to each pulse peak included in the target echo includes:
dividing the target echo into a plurality of echo fragments according to the length of the main pulse peak, wherein the length of each echo fragment is equal to the length of the main pulse peak, and each echo fragment comprises a pulse peak;
For each of the echo segments: comparing the main pulse peak with the current echo fragment, and if the comparison result meets the set first condition and/or the set second condition, determining the pulse peak included in the current echo fragment as the pulse peak to be processed.
Preferably, said comparing said main pulse peak with the current echo segment comprises:
Determining a corresponding amplitude range according to the main pulse peak;
And carrying out frequency domain analysis on the current echo fragment, determining whether the maximum amplitude of the spectrum power of the current echo fragment is in the amplitude range after the frequency domain analysis, and if so, determining that the comparison result meets the first condition.
Preferably, said comparing said main pulse peak with the current echo segment comprises:
and calculating the similarity between the main pulse peak and the current echo fragment, and if the similarity is larger than a set similarity threshold, determining that the comparison result meets the second condition.
Preferably, the determining whether the pulse peak to be processed is a real pulse peak according to the amplitude of the second sub-echo and the amplitude of the pulse peak to be processed includes:
determining a first distance length and a second distance length according to the distance between the pulse peak to be processed and the adjacent pulse peak;
Determining a third sub-echo between a first position and a second position on the second sub-echo, the first position being located after the main pulse peak and at a distance from the main pulse peak of the first distance length, the second position being located after the main pulse peak and at a distance from the main pulse peak of the second distance length;
And carrying out normalization processing on the pulse peak to be processed and the third sub-echo, judging whether the peak value of the pulse peak to be processed after normalization processing is larger than or equal to the maximum amplitude value of the third sub-echo after normalization processing, and if so, determining that the pulse peak to be processed is a real pulse peak.
Preferably, the determining the corresponding amplitude range according to the main pulse peak includes:
Carrying out frequency domain analysis on the main pulse peak to obtain the maximum amplitude of the spectrum power of the main pulse peak, and determining the maximum amplitude of the spectrum power of the main pulse peak as the maximum value of the amplitude range; and determining the difference value between the maximum amplitude value of the spectrum power of the main pulse front and a preset value as the minimum value of the amplitude range.
In a second aspect, an embodiment of the present invention further provides a target echo processing device, including:
the echo acquisition unit is used for acquiring a reference echo of the calibration body and acquiring a target echo to be processed;
a main pulse peak determining unit for determining a main pulse peak in the reference echo;
The identification unit is used for identifying whether each pulse peak contained in the target echo is a real pulse peak or not according to the main pulse peak;
And the echo generating unit is used for obtaining the processed target echo according to the identification result.
In a third aspect, an embodiment of the present invention further provides a computing device, including a memory and a processor, where the memory stores a computer program, and the processor implements a method according to any embodiment of the present specification when executing the computer program.
In a fourth aspect, embodiments of the present invention also provide a computer-readable storage medium having stored thereon a computer program which, when executed in a computer, causes the computer to perform a method according to any of the embodiments of the present specification.
The embodiment of the invention provides a target terahertz time-domain echo processing method, device, equipment and storage medium, wherein a reference echo of a calibration body is used as a reference for comparison with a target echo, so that a main pulse peak in the reference echo can be utilized to automatically identify whether each pulse peak in the target echo is a real pulse peak, and the identification process is an automatic identification process without human participation, so that the identification efficiency and accuracy can be improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 is a flowchart of a target terahertz time-domain echo processing method according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of reference echoes in the time domain of a calibration body according to an embodiment of the present invention;
FIG. 3 is a time domain diagram of a pre-processed target echo in an embodiment of the invention;
FIG. 4 is a flow chart of a method for identifying real pulse peaks in an embodiment of the present invention;
FIG. 5 is a time domain diagram of a target echo obtained after processing in an embodiment of the present invention;
FIG. 6 is a hardware architecture diagram of a computing device according to one embodiment of the invention;
Fig. 7 is a block diagram of a target terahertz time-domain echo processing device according to an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments, and all other embodiments obtained by those skilled in the art without making any inventive effort based on the embodiments of the present invention are within the scope of protection of the present invention.
As described above, in the conventional target echo processing method, it is often manually determined whether the target echo is a real pulse peak or not simply according to the shape and amplitude of the target echo, and the process is complicated and easy to miss and misjudge. In order to solve the problem in the process of human participation identification, an automatic identification mode can be considered for identification. In order to enable automatic recognition of the target echo, it is conceivable to use the reference echo of the calibration body as a reference for comparison of the target echo in order to recognize whether each pulse peak included in the target echo is a true pulse peak.
Specific implementations of the above concepts are described below.
Referring to fig. 1, an embodiment of the present invention provides a target terahertz time-domain echo processing method, which includes:
Step 100, obtaining a reference echo of the calibration body and obtaining a target echo to be processed.
Step 102, determining a main pulse peak in the reference echo.
And step 104, identifying whether each pulse peak contained in the target echo is a real pulse peak according to the main pulse peak.
And 106, obtaining the processed target echo according to the identification result.
In the target terahertz time-domain echo processing method shown in fig. 1, by taking the reference echo of the calibration body as a reference for comparison with the target echo, whether each pulse peak in the target echo is a real pulse peak or not can be automatically identified by utilizing the main pulse peak in the reference echo, and the identification process is an automatic identification process without human participation, so that the identification efficiency and the accuracy can be improved.
The manner in which the individual steps shown in fig. 1 are performed is described below.
First, for step 100, a reference echo of the calibration volume is acquired, and a target echo to be processed is acquired.
In the embodiment of the invention, the reference echo is a terahertz echo obtained by carrying out terahertz radiation measurement on the calibration body. The reference echo is a reference for comparison with the target echo, and the accuracy of identifying each pulse peak in the target echo is higher as the noise in the reference echo is smaller. When the number of reflecting surfaces of the calibration body is two or more, noise in the reference echo will be larger, and thus, when the calibration body is selected, a calibration body having a single reflecting surface, such as a metal ball, a flat plate, or the like, may be selected.
In addition, in order to be able to acquire the scattering characteristics of the complete target, it is necessary to ensure that the time domain range measured when the terahertz radiation is performed on the target is greater than or equal to the length corresponding to the actual time domain range of the target, and in order to ensure that the reference echo can be used as a reference for comparison with the target echo, it is also necessary to ensure that the reference echo for performing the terahertz radiation measurement on the calibration body should be ensured, and the time domain range should be able to cover the time domain range for performing the terahertz radiation measurement on the target to the target. For example, the length corresponding to the time domain range of the target echo is 300ps, and the length corresponding to the time domain range of the reference echo measured on the calibration body should be not less than 300ps.
In this step 100, the target echo may be acquired first, then the time domain range that needs to be satisfied by the reference echo is determined based on the time domain range of the target echo, and then the reference echo is acquired according to the determined time domain range. In this step 100, the reference echo may be acquired first, and then the target echo may be acquired, or the reference echo and the target echo may be acquired simultaneously, which may be understood that when the reference echo is acquired in this manner, the reference echo having a larger time domain range may be acquired, so as to ensure that the time domain range of the target echo can be covered.
For example, a metal sphere is selected as a calibration body, and the calibration body is measured to obtain a reference echo s_cali (N) shown in fig. 2, where n=1..n, where s_cali (N) is the amplitude corresponding to the nth sampling point, and N is the total number of sampling points.
In the embodiment of the present invention, in order to ensure that the reference echo can be used as a reference for comparison with the target echo, in this step 100, when the target is subjected to terahertz radiation to measure the target echo to be processed, and the calibration body is subjected to terahertz radiation to measure the reference echo, measurement parameters in the measurement process of the reference echo and the target echo need to be kept the same. For example, the measurement parameters of bias voltage, humidity and the like in the measurement device, and the measurement parameters of sampling interval, phase-locked integration time and the like all need to be kept at the same setting.
In order to obtain scattering characteristics of the target in different postures, target echoes of the target in different postures can be measured, and the target echoes in each posture can be processed by using the embodiment.
For example, an aircraft model may be selected as a target, and a target echo of s_tar (n) of the aircraft model in a certain attitude is measured, as shown in fig. 3, where A, B, C, D, E, F in fig. 3 are respectively a part of pulse peaks in the target echo.
Then, for step 102, a main pulse peak in the reference echo is determined.
Taking the reference echo shown in fig. 2 as an example, a section of echo intercepted by a rectangular frame in fig. 2 is taken as a complete main pulse front in the reference echo, denoted as s_cali_std (K), where k=1..k, where s_cali_std (K) represents an amplitude corresponding to a kth sampling point, and K is a total number of sampling points included in a main pulse peak.
Next, for step 104, it is identified from the main pulse peak whether each pulse peak included in the target echo is a true pulse peak.
In an embodiment of the present invention, referring to fig. 4, at least the following steps 400-406 may be used for identification:
And 400, determining at least one pulse peak to be processed according to each pulse peak included in the target echo.
In the embodiment of the present invention, the manner of determining the pulse peak to be processed may at least include the following two ways:
Mode one: all pulse peaks included in the target echo are determined as pulse peaks to be processed.
Mode two: and setting screening conditions to screen each pulse peak included in the target echo, and determining the screened pulse peak as a pulse peak to be processed.
Specifically, when determining the pulse peak to be processed by the second mode, S1-S2 may be included:
S1: dividing the target echo into a plurality of echo fragments according to the length of the main pulse peak, wherein the length of each echo fragment is equal to the length of the main pulse peak, and each echo fragment comprises a pulse peak.
In the embodiment of the invention, the main pulse peak in the reference echo is used as a comparison standard to identify whether each pulse peak contained in the target echo is a real pulse peak, and in order to better identify, the target echo needs to be divided into a plurality of echo fragments with the same length as the main pulse peak.
Wherein the target echo may be divided using a hamming window, a rectangular window, etc. In order to ensure that the divided echo fragments are equal to the length of the main pulse peak, the width of the hamming window and the rectangular window can be set to be the length of the main pulse peak.
For example, since the total number of sampling points of the main pulse peak in the reference echo is K, a hamming window with a window length of K is set to divide the target echo into a plurality of echo segments in turn, and serial numbers are respectively performed on each echo segment, denoted as s_tra_i (z), where z represents the serial number of the echo segment in the target echo.
S2: for each of the echo segments: comparing the main pulse peak with the current echo fragment, and if the comparison result meets the set first condition and/or the set second condition, determining the pulse peak included in the current echo fragment as the pulse peak to be processed.
In determining whether each echo segment is a pulse peak to be processed, since the determination manner is the same, in the embodiment of the present invention, one of the echo segments s_tra_i (z) may be used as a current echo segment, and the determination of the current echo segment is illustrated as an example.
In one embodiment of the present invention, the first condition (hereinafter, described as condition one) and/or the second condition (hereinafter, described as condition two) may be preset to determine the comparison result of the main pulse peak and the current echo segment. The following describes whether the comparison result satisfies the first condition or the second condition.
First, it is explained whether the comparison result satisfies the first condition.
In the step S2, comparing the main pulse peak with the current echo segment may include: determining a corresponding amplitude range according to the main pulse peak; and carrying out frequency domain analysis on the current echo fragment, determining whether the maximum amplitude of the spectrum power of the current echo fragment is in the amplitude range after the frequency domain analysis, and if so, determining that the comparison result meets a first condition.
It should be noted that, the method for determining the corresponding amplitude range according to the main pulse peak may adopt a plurality of determination methods, and in one embodiment of the present invention, the amplitude range may be determined at least according to the maximum amplitude of the spectrum power of the main pulse peak, which specifically may be the following steps: carrying out frequency domain analysis on the main pulse peak to obtain the maximum amplitude of the spectrum power of the main pulse peak, and determining the maximum amplitude of the spectrum power of the main pulse peak as the maximum value of the amplitude range; and determining the difference value between the maximum amplitude value of the spectrum power of the main pulse front and a preset value as the minimum value of the amplitude range.
In the embodiment of the present invention, the preset value may be set according to the recognition accuracy, for example, may be set to 5dB. When the main pulse peak is subjected to frequency domain analysis, the main pulse peak can be converted into the spectrum power amplitude by adopting Fourier transformation and the like, so that the spectrum power maximum amplitude of the main pulse peak is obtained, for example, the spectrum power maximum amplitude of the main pulse peak is-79 dB, and the amplitude range can be [ -84, -79dB ].
In another embodiment of the present invention, the amplitude range may also be determined using the following method: and determining the frequency corresponding to the available frequency band in the power spectrogram of the pulse peak to be processed, determining the amplitude of the main pulse peak corresponding to the frequency as the maximum value of the amplitude range, and determining the interpolation of the amplitude corresponding to the frequency and the set value as the minimum value of the amplitude range.
For example, a frequency band of 20dB higher than the spectrum power of the background noise in the power spectrum of the pulse peak to be processed is determined as an available frequency band, the frequency corresponding to the available frequency band is 1.5THz, and then the amplitude corresponding to the main pulse peak is determined as the maximum value of the amplitude range.
Then, whether the comparison result satisfies the second condition is described.
When determining whether the comparison result satisfies the second condition, in the step S2, comparing the main pulse peak with the current echo segment may include: and calculating the similarity between the main pulse peak and the current echo fragment, and if the similarity is larger than a set similarity threshold, determining that the comparison result meets a second condition.
In the embodiment of the present invention, the manner of calculating the similarity may include the following manner of calculating the similarity. For example, if the similarity threshold is set to 0.8, when the calculated similarity is greater than 0.8, it is determined that the comparison result satisfies the second condition.
In using similarity calculations, the formula may be utilizedCalculation is performed, wherein Cov (X, Y) is the covariance of X and Y, var [ X ] is the variance of X, var [ Y ] is the variance of Y, X, Y are the current echo segment and the main pulse peak, respectively.
It should be noted that, each echo segment needs to be compared with the main pulse peak, if the comparison result of the current echo segment meets the first condition, the current echo segment can be determined as the pulse peak to be processed; if the comparison result of the current echo fragment meets the second condition, the current echo fragment can be determined to be a pulse peak to be processed; if the comparison result of the current echo fragment meets the first condition and the second condition at the same time, the current echo fragment can be determined as a pulse peak to be processed.
Step 402, determining, for each of the to-be-processed pulse peaks, the to-be-processed pulse peak and an adjacent pulse peak located before the to-be-processed pulse peak as a first sub-echo.
For example, the pulse peak to be processed is L1, and the adjacent pulse peak before L1 is L0, then the adjacent pulse peak L0 and the pulse peak to be processed L1 form a first sub-echo. The length of the first sub-echo is the time domain distance between the start point of the adjacent pulse peak L0 and the end point of the pulse peak L1 to be processed.
And step 404, determining a second sub-echo from the reference echo, wherein the starting point of the second sub-echo is the main pulse peak, and the length of the second sub-echo is equal to that of the first sub-echo.
For example, the second sub-echo starts from the start point of the main pulse peak L, and the main pulse peak L and a section of echo following the main pulse peak L are intercepted from the reference echo, so that the length of the obtained second sub-echo is the same as the length of the first sub-echo.
And step 406, judging whether the pulse peak to be processed is a real pulse peak or not according to the amplitude of the second sub-echo and the amplitude of the pulse peak to be processed.
In one embodiment of the present invention, in step 406, when determining whether the pulse peak to be processed is a true pulse peak, steps H1-H3 may be specifically included:
H1: and determining a first distance length and a second distance length according to the distance between the pulse peak to be processed and the adjacent pulse peak.
And if the time domain distance between the peak point of the pulse peak L1 to be processed in the first sub-echo and the peak point of the adjacent pulse peak L0 is d and the length of the main pulse peak is K, determining d-K as a first distance length and d+K as a second distance length.
H2: and determining a third sub-echo between a first position and a second position on the second sub-echo, wherein the first position is positioned behind the main pulse peak and is at a distance of the first distance length from the main pulse peak, and the second position is positioned behind the main pulse peak and is at a distance of the second distance length from the main pulse peak.
And H3: and carrying out normalization processing on the pulse peak to be processed and the third sub-echo, judging whether the peak value of the pulse peak to be processed after normalization processing is larger than or equal to the maximum amplitude value of the third sub-echo after normalization processing, and if so, determining that the pulse peak to be processed is a real pulse peak.
For comparison, the peak point of the main pulse peak L in the second sub-echo and the peak point of the adjacent pulse peak L0 in the first sub-echo may be aligned, and then whether the peak value of the normalized pulse peak to be processed is greater than or equal to the maximum amplitude of the normalized third sub-echo is compared. If yes, determining the pulse peak L1 to be processed as a real pulse peak.
Steps 402-406 may be determined as condition three, with which each pulse peak to be processed is identified to identify whether it is a true pulse peak.
In the embodiment of the invention, all pulse peaks in the target echo are determined to be processed in the first mode, namely, each pulse peak in the target echo is identified in the third mode, and all pulse peaks can be identified in parallel due to the same identification process, so that the identification speed can be improved.
In the embodiment of the invention, the first pulse peak and/or the second pulse peak are/is firstly identified by the second utilization mode, then the second identification is carried out by the third utilization mode on the first pulse peak and/or the second pulse peak which meet the first condition and/or the second condition, and the pulse peak which meets the conditions in the second identification is identified as the real pulse peak, so that the accuracy of the identification result can be ensured.
Finally, for step 106, the processed target echo is obtained according to the recognition result.
Taking the target echo shown in fig. 3 as an example, a section of echo taken by a rectangular frame in fig. 3 includes A, B, C, D, E, F pulse peaks. After the automatic recognition in step 104 of this embodiment, a target echo containing a real pulse peak as shown in fig. 5 is obtained, and the real pulse peaks in the pulse peaks included in the section of echo intercepted by the rectangular frame in fig. 5 are A, B, C, E, F pulse peaks, so that it can be seen that the pulse peak D is a non-real pulse peak, and the waveform formed by the real pulse peak is determined as the processed target echo.
As shown in fig. 6 and 7, the embodiment of the invention provides a target terahertz time-domain echo processing device. The apparatus embodiments may be implemented by software, or may be implemented by hardware or a combination of hardware and software. In terms of hardware, as shown in fig. 6, a hardware architecture diagram of a computing device where a target terahertz time-domain echo processing apparatus provided by an embodiment of the present invention is located, where in addition to a processor, a memory, a network interface, and a nonvolatile memory shown in fig. 6, the computing device where the embodiment is located may generally include other hardware, such as a forwarding chip responsible for processing a packet, and so on. Taking a software implementation as an example, as shown in fig. 7, as a device in a logic sense, the device is formed by reading a corresponding computer program in a nonvolatile memory into a memory by a CPU of a computing device where the device is located. The target terahertz time-domain echo processing device provided in this embodiment includes:
An echo acquisition unit 701, configured to acquire a reference echo of a calibration body and acquire a target echo to be processed;
A main pulse peak determining unit 702 for determining a main pulse peak in the reference echo;
a recognizing unit 703, configured to recognize whether each pulse peak included in the target echo is a real pulse peak according to the main pulse peak;
and the echo generating unit 704 is configured to obtain a processed target echo according to the identification result.
In one embodiment of the present invention, the identifying unit 703 is specifically configured to determine at least one pulse peak to be processed according to each pulse peak included in the target echo; for each pulse peak to be processed, determining the pulse peak to be processed and an adjacent pulse peak positioned before the pulse peak to be processed as a first sub-echo; determining a second sub-echo from the reference echo, wherein the starting point of the second sub-echo is the main pulse peak, and the length of the second sub-echo is equal to that of the first sub-echo; and judging whether the pulse peak to be processed is a real pulse peak or not according to the amplitude of the second sub-echo and the amplitude of the pulse peak to be processed.
In one embodiment of the present invention, when determining at least one pulse peak to be processed according to each pulse peak included in the target echo, the identifying unit 703 is specifically configured to divide the target echo into a plurality of echo segments according to the length of the main pulse peak, where the length of each echo segment is equal to the length of the main pulse peak, and each echo segment includes one pulse peak; for each of the echo segments: comparing the main pulse peak with the current echo fragment, and if the comparison result meets the set first condition and/or the set second condition, determining the pulse peak included in the current echo fragment as the pulse peak to be processed.
In one embodiment of the present invention, the identifying unit 703 is specifically configured to determine a corresponding amplitude range according to the main pulse peak when performing the comparison between the main pulse peak and the current echo segment; and carrying out frequency domain analysis on the current echo fragment, determining whether the maximum amplitude of the spectrum power of the current echo fragment is in the amplitude range after the frequency domain analysis, and if so, determining that the comparison result meets the first condition.
In one embodiment of the present invention, the identifying unit 703 is specifically configured to calculate a similarity between the main pulse peak and the current echo segment when performing the comparison between the main pulse peak and the current echo segment, and determine that the comparison result meets the second condition if the similarity is greater than a set similarity threshold.
In one embodiment of the present invention, the identifying unit 703 is specifically configured to determine a first distance length and a second distance length according to a distance between the pulse peak to be processed and the adjacent pulse peak when determining whether the pulse peak to be processed is a real pulse peak according to the amplitude of the second sub-echo and the amplitude of the pulse peak to be processed; determining a third sub-echo between a first position and a second position on the second sub-echo, the first position being located after the main pulse peak and at a distance from the main pulse peak of the first distance length, the second position being located after the main pulse peak and at a distance from the main pulse peak of the second distance length; and carrying out normalization processing on the pulse peak to be processed and the third sub-echo, judging whether the peak value of the pulse peak to be processed after normalization processing is larger than or equal to the maximum amplitude value of the third sub-echo after normalization processing, and if so, determining that the pulse peak to be processed is a real pulse peak.
In one embodiment of the present invention, the identifying unit 703 is specifically configured to perform frequency domain analysis on the main pulse peak when the determining of the corresponding amplitude range according to the main pulse peak is performed, so as to obtain a maximum amplitude of the spectrum power of the main pulse peak, and determine the maximum amplitude of the spectrum power of the main pulse peak as a maximum value of the amplitude range; and determining the difference value between the maximum amplitude value of the spectrum power of the main pulse front and a preset value as the minimum value of the amplitude range.
It will be appreciated that the structure illustrated in the embodiments of the present invention does not constitute a specific limitation on a target terahertz time-domain echo processing apparatus. In other embodiments of the present invention, a target terahertz time-domain echo processing device may include more or fewer components than shown, or certain components may be combined, certain components may be split, or different component arrangements. The illustrated components may be implemented in hardware, software, or a combination of software and hardware.
The content of information interaction and execution process between the modules in the device is based on the same conception as the embodiment of the method of the present invention, and specific content can be referred to the description in the embodiment of the method of the present invention, which is not repeated here.
The embodiment of the invention also provides a computing device, which comprises a memory and a processor, wherein the memory stores a computer program, and the processor realizes the target terahertz time-domain echo processing method in any embodiment of the invention when executing the computer program.
The embodiment of the invention also provides a computer readable storage medium, and the computer readable storage medium stores a computer program, and when the computer program is executed by a processor, the computer program causes the processor to execute the target terahertz time-domain echo processing method in any embodiment of the invention.
Specifically, a system or apparatus provided with a storage medium on which a software program code realizing the functions of any of the above embodiments is stored, and a computer (or CPU or MPU) of the system or apparatus may be caused to read out and execute the program code stored in the storage medium.
In this case, the program code itself read from the storage medium may realize the functions of any of the above-described embodiments, and thus the program code and the storage medium storing the program code form part of the present invention.
Examples of storage media for providing program code include floppy disks, hard disks, magneto-optical disks, optical disks (e.g., CD-ROMs, CD-R, CD-RWs, DVD-ROMs, DVD-RAMs, DVD-RWs, DVD+RWs), magnetic tapes, nonvolatile memory cards, and ROMs. Alternatively, the program code may be downloaded from a server computer by a communication network.
Further, it should be apparent that the functions of any of the above-described embodiments may be implemented not only by executing the program code read out by the computer, but also by causing an operating system or the like operating on the computer to perform part or all of the actual operations based on the instructions of the program code.
Further, it is understood that the program code read out by the storage medium is written into a memory provided in an expansion board inserted into a computer or into a memory provided in an expansion module connected to the computer, and then a CPU or the like mounted on the expansion board or the expansion module is caused to perform part and all of actual operations based on instructions of the program code, thereby realizing the functions of any of the above embodiments.
It is noted that relational terms such as first and second, and the like, are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
Those of ordinary skill in the art will appreciate that: all or part of the steps for implementing the above method embodiments may be implemented by hardware related to program instructions, and the foregoing program may be stored in a computer readable storage medium, where the program, when executed, performs steps including the above method embodiments; and the aforementioned storage medium includes: various media in which program code may be stored, such as ROM, RAM, magnetic or optical disks.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and are not limiting; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.
Claims (8)
1. The target terahertz time-domain echo processing method is characterized by comprising the following steps of:
Acquiring a reference echo of a calibration body and acquiring a target echo to be processed;
Determining a main pulse peak in the reference echo;
Identifying whether each pulse peak contained in the target echo is a real pulse peak or not according to the main pulse peak;
obtaining a processed target echo according to the identification result;
The identifying whether each pulse peak included in the target echo is a true pulse peak according to the main pulse peak includes:
Determining at least one pulse peak to be processed according to each pulse peak included in the target echo;
For each pulse peak to be processed, determining the pulse peak to be processed and an adjacent pulse peak positioned before the pulse peak to be processed as a first sub-echo;
determining a second sub-echo from the reference echo, wherein the starting point of the second sub-echo is the main pulse peak, and the length of the second sub-echo is equal to that of the first sub-echo;
Judging whether the pulse peak to be processed is a real pulse peak or not according to the amplitude of the second sub-echo and the amplitude of the pulse peak to be processed;
And judging whether the pulse peak to be processed is a real pulse peak according to the amplitude of the second sub-echo and the amplitude of the pulse peak to be processed, wherein the method comprises the following steps:
determining a first distance length and a second distance length according to the distance between the pulse peak to be processed and the adjacent pulse peak;
Determining a third sub-echo between a first position and a second position on the second sub-echo, the first position being located after the main pulse peak and at a distance from the main pulse peak of the first distance length, the second position being located after the main pulse peak and at a distance from the main pulse peak of the second distance length;
And carrying out normalization processing on the pulse peak to be processed and the third sub-echo, judging whether the peak value of the pulse peak to be processed after normalization processing is larger than or equal to the maximum amplitude value of the third sub-echo after normalization processing, and if so, determining that the pulse peak to be processed is a real pulse peak.
2. The method of claim 1, wherein determining at least one pulse peak to be processed from each pulse peak included in the target echo comprises:
dividing the target echo into a plurality of echo fragments according to the length of the main pulse peak, wherein the length of each echo fragment is equal to the length of the main pulse peak, and each echo fragment comprises a pulse peak;
For each of the echo segments: comparing the main pulse peak with the current echo fragment, and if the comparison result meets the set first condition and/or the set second condition, determining the pulse peak included in the current echo fragment as the pulse peak to be processed.
3. The method of claim 2, wherein said comparing said main pulse peak with said current echo segment comprises:
Determining a corresponding amplitude range according to the main pulse peak;
And carrying out frequency domain analysis on the current echo fragment, determining whether the maximum amplitude of the spectrum power of the current echo fragment is in the amplitude range after the frequency domain analysis, and if so, determining that the comparison result meets the first condition.
4. The method of claim 2, wherein said comparing said main pulse peak with said current echo segment comprises:
and calculating the similarity between the main pulse peak and the current echo fragment, and if the similarity is larger than a set similarity threshold, determining that the comparison result meets the second condition.
5. A method according to claim 3, wherein said determining a corresponding amplitude range from said main pulse peak comprises:
Carrying out frequency domain analysis on the main pulse peak to obtain the maximum amplitude of the spectrum power of the main pulse peak, and determining the maximum amplitude of the spectrum power of the main pulse peak as the maximum value of the amplitude range;
And determining the difference value between the maximum amplitude value of the spectrum power of the main pulse front and a preset value as the minimum value of the amplitude range.
6. A target terahertz time-domain echo processing apparatus, characterized by comprising:
the echo acquisition unit is used for acquiring a reference echo of the calibration body and acquiring a target echo to be processed;
a main pulse peak determining unit for determining a main pulse peak in the reference echo;
The identification unit is used for identifying whether each pulse peak contained in the target echo is a real pulse peak or not according to the main pulse peak;
The echo generating unit is used for obtaining a processed target echo according to the identification result;
The identification unit is specifically configured to determine at least one pulse peak to be processed according to each pulse peak included in the target echo; for each pulse peak to be processed, determining the pulse peak to be processed and an adjacent pulse peak positioned before the pulse peak to be processed as a first sub-echo; determining a second sub-echo from the reference echo, wherein the starting point of the second sub-echo is the main pulse peak, and the length of the second sub-echo is equal to that of the first sub-echo; judging whether the pulse peak to be processed is a real pulse peak or not according to the amplitude of the second sub-echo and the amplitude of the pulse peak to be processed;
The identifying unit is specifically configured to determine a first distance length and a second distance length according to a distance between the pulse peak to be processed and the adjacent pulse peak when determining whether the pulse peak to be processed is a real pulse peak according to the amplitude of the second sub-echo and the amplitude of the pulse peak to be processed; determining a third sub-echo between a first position and a second position on the second sub-echo, the first position being located after the main pulse peak and at a distance from the main pulse peak of the first distance length, the second position being located after the main pulse peak and at a distance from the main pulse peak of the second distance length; and carrying out normalization processing on the pulse peak to be processed and the third sub-echo, judging whether the peak value of the pulse peak to be processed after normalization processing is larger than or equal to the maximum amplitude value of the third sub-echo after normalization processing, and if so, determining that the pulse peak to be processed is a real pulse peak.
7. A computing device comprising a memory and a processor, the memory having stored therein a computer program, the processor implementing the method of any of claims 1-5 when the computer program is executed.
8. A computer readable storage medium having stored thereon a computer program which, when executed in a computer, causes the computer to perform the method of any of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111224460.9A CN113960555B (en) | 2021-10-19 | 2021-10-19 | Target terahertz time-domain echo processing method, device, equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111224460.9A CN113960555B (en) | 2021-10-19 | 2021-10-19 | Target terahertz time-domain echo processing method, device, equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113960555A CN113960555A (en) | 2022-01-21 |
| CN113960555B true CN113960555B (en) | 2024-06-07 |
Family
ID=79465208
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111224460.9A Active CN113960555B (en) | 2021-10-19 | 2021-10-19 | Target terahertz time-domain echo processing method, device, equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113960555B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0841548A2 (en) * | 1996-11-06 | 1998-05-13 | Lucent Technologies Inc. | Systems and methods for processing and analyzing terahertz waveforms |
| CA2692942A1 (en) * | 2007-07-12 | 2009-01-15 | Picometrix, Llc | System and method to measure the transit time position(s) of pulses in a time domain data |
| CN211856905U (en) * | 2020-03-15 | 2020-11-03 | 上海铭剑电子科技有限公司 | Multifunctional active terahertz imaging system device |
-
2021
- 2021-10-19 CN CN202111224460.9A patent/CN113960555B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0841548A2 (en) * | 1996-11-06 | 1998-05-13 | Lucent Technologies Inc. | Systems and methods for processing and analyzing terahertz waveforms |
| CA2692942A1 (en) * | 2007-07-12 | 2009-01-15 | Picometrix, Llc | System and method to measure the transit time position(s) of pulses in a time domain data |
| CN211856905U (en) * | 2020-03-15 | 2020-11-03 | 上海铭剑电子科技有限公司 | Multifunctional active terahertz imaging system device |
Non-Patent Citations (1)
| Title |
|---|
| 太赫兹时域光谱分析中的回波滤除方法;王凌辉;王迎新;赵自然;红外与毫米波学报;20131231;第32卷(第3期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113960555A (en) | 2022-01-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2041667B1 (en) | Method and apparatus for target discrimination within return signals | |
| CN109564270B (en) | Positioning information determining method and device | |
| CN113567604A (en) | Detection and analysis method of chromatographic spectrogram and electronic equipment | |
| CN113567603A (en) | Detection and analysis method of chromatographic spectrogram and electronic equipment | |
| CN110988834B (en) | Pulse arrival time measuring method based on self-adaptive threshold value double thresholds | |
| US20170296081A1 (en) | Frame based spike detection module | |
| CN111239704A (en) | Atmospheric detection radar target echo identification processing method, device, equipment and medium | |
| CN113050071A (en) | Laser radar data processing method, device, equipment and storage medium | |
| US10378937B2 (en) | Dynamic clustering of transient signals | |
| CN111999726B (en) | Personnel positioning method based on millimeter wave radar | |
| CN113960555B (en) | Target terahertz time-domain echo processing method, device, equipment and storage medium | |
| CN109752633B (en) | Method and system for positioning partial discharge signal of transformer substation | |
| CN115166650B (en) | Radar signal identification and parameter estimation method and system | |
| JP2010127833A (en) | Radio source identification device | |
| RU2694271C2 (en) | Device for classifying noisy objects | |
| CN112436863B (en) | Signal rapid detection search method based on frequency spectrum dispersion analysis | |
| CN110632563A (en) | Intra-pulse frequency coding signal parameter measuring method based on short-time Fourier transform | |
| Shirer et al. | Estimating the correlation dimension of atmospheric time series | |
| JP2014032082A (en) | Passive sonar device, azimuth concentration processing method, and passive sonar signal processing program | |
| US20100153063A1 (en) | Automatic procedure for merging tracks and estimating harmonic combs | |
| CN115541021A (en) | Method for locating characteristic peak of Raman spectrum, electronic device and storage medium | |
| Aouada et al. | Application of the bootstrap to source detection in nonuniform noise | |
| CN116403080B (en) | Face clustering evaluation method, system, computer and readable storage medium | |
| CN115598615B (en) | Power spectrum information geometric radar target detection method and device based on sub-band filtering | |
| CN109581423A (en) | Linear frequency sweep interference suppression method and system using subspace projection |
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 |