Terahertz ISAR Imaging Joint Compensation Method, Device, Equipment and Medium

Abstract

The present invention relates to a terahertz ISAR imaging joint compensation method, device, equipment and medium. The method includes: performing fractional Fourier transform on the extracted echo data, then performing mutation error elimination and least squares fitting to obtain a rough estimated motion parameter; solving the radial displacement difference between adjacent echoes based on the rough estimated motion parameter; then calculating a rough joint compensation function; performing rough joint compensation on the target echo data through the rough joint compensation function to obtain roughly compensated echo data; solving the synthetic waveform entropy based on the roughly compensated echo data; then using the rough estimated motion parameter as the iterative initial value and the synthetic waveform entropy as the fitness function, and performing iterative optimization solution through the GWO algorithm to obtain finely compensated echo data; performing azimuth Fourier transform on the finely compensated echo data and outputting the final two-dimensional ISAR image. The present invention can effectively eliminate adverse effects such as frequency offset and defocus caused by target motion, and has high compensation accuracy.

Description

Technical Field

[0001] The present invention relates to the technical field of radar signal processing, and particularly to a terahertz ISAR imaging joint compensation method, device, equipment, and medium. Background Art

[0002] Inverse Synthetic Aperture Radar (ISAR) is proposed relative to Synthetic Aperture Radar (SAR), mainly for imaging non-cooperative targets such as ships, aircraft, and satellites, and has characteristics such as all-weather, all-time, and obtaining high-resolution imaging of targets over long distances. The key step in ISAR imaging lies in translational compensation, that is, compensating for the radial motion of the target relative to the Radar Line of Sight (RLOS), so that the target only rotates around the equivalent rotation center relative to the radar.

[0003] Slow targets in traditional scenarios are usually imaged using the walk-stop model, that is, considering that the target is stationary relative to the radar within the time of a single pulse emitted by the radar, and there is only range walk between different pulses emitted by the radar. For high-speed space targets such as satellites, the walk-stop model is no longer applicable. In addition to inter-pulse compensation, intra-pulse compensation of the echo also needs to be considered, otherwise it will greatly affect the ISAR imaging quality.

[0004] Terahertz (THz) waves refer to electromagnetic waves with a frequency range between 0.1 THz and 10 THz, located between microwaves and infrared, and in the transition frequency band from electronics to photonics. Compared with infrared and lasers, terahertz waves have advantages such as good smoke and dust penetration and low sensitivity to aero-optical effects, and can be used for combat in complex environments and space target detection. At the same time, compared with the microwave frequency band, terahertz waves have a shorter wavelength and a larger bandwidth, and have important application value in the application scenario of radar high-resolution imaging. However, at the same time, compared with radars in the microwave frequency band, when a terahertz radar performs ISAR imaging, its higher range resolution also leads to a more serious intra-pulse Doppler phenomenon caused by the high-speed movement of the imaging target. Severe intra-pulse Doppler modulation will cause frequency offset and defocusing in the range image, affecting the subsequent imaging effect. Summary of the Invention

[0005] Based on this, it is necessary to provide a terahertz ISAR imaging joint compensation method, device, equipment, and medium that can solve problems such as frequency offset and defocusing in the range image and low imaging accuracy for the above technical problems.

[0006] A terahertz ISAR imaging joint compensation method, the method includes:

[0007] Obtain the echo data of the target;

[0008] Extract the echo data, then perform fractional Fourier transform on the extracted echo data, and then perform mutation error elimination and least squares fitting to obtain the rough estimated motion parameters;

[0009] Solve the radial displacement difference between adjacent echoes caused by translational motion based on the rough estimated motion parameters; then calculate the rough joint compensation function according to the rough estimated motion parameters and the radial displacement difference;

[0010] Perform rough joint compensation on the target echo data through the rough joint compensation function to obtain the roughly compensated echo data;

[0011] Solve the synthetic waveform entropy based on the roughly compensated echo data; then use the rough estimated motion parameters as the iterative initial value and the synthetic waveform entropy as the fitness function, and perform iterative optimization through the GWO algorithm to obtain the precisely compensated echo data;

[0012] Perform azimuth Fourier transform on the precisely compensated echo data and output the final two-dimensional ISAR image.

[0013] A terahertz ISAR imaging joint compensation device, the device includes:

[0014] An echo data acquisition module, used to acquire the echo data of the target;

[0015] A rough estimated motion parameter acquisition module, used to extract the echo data, then perform fractional Fourier transform on the extracted echo data, and then perform mutation error elimination and least squares fitting to obtain the rough estimated motion parameters;

[0016] A rough joint compensation function acquisition module, used to solve the radial displacement difference between adjacent echoes caused by translational motion based on the rough estimated motion parameters; then calculate the rough joint compensation function according to the rough estimated motion parameters and the radial displacement difference;

[0017] A roughly compensated echo data calculation module, used to perform rough joint compensation on the target echo data through the rough joint compensation function to obtain the roughly compensated echo data;

[0018] Precisely compensated echo data calculation, used to solve the synthetic waveform entropy based on the roughly compensated echo data; then use the rough estimated motion parameters as the iterative initial value and the synthetic waveform entropy as the fitness function, and perform iterative optimization through the GWO algorithm to obtain the precisely compensated echo data;

[0019] A two-dimensional ISAR image generation module, used to perform azimuth Fourier transform on the precisely compensated echo data and output the final two-dimensional ISAR image.

[0020] A computer device includes a memory and a processor. The memory stores a computer program, and when the processor executes the computer program, the steps of the terahertz ISAR imaging joint compensation method are implemented.

[0021] A computer-readable storage medium stores a computer program thereon. When the computer program is executed by a processor, the steps of the terahertz ISAR imaging joint compensation method are implemented.

[0022] For the above-mentioned terahertz ISAR imaging joint compensation method, device, equipment and medium, the echo data of the target is obtained; the echo data is extracted, then the extracted echo data is subjected to fractional Fourier transform, and then mutation error elimination and least square fitting are performed to obtain a rough estimated motion parameter; based on the rough estimated motion parameter, the radial displacement difference between adjacent echoes caused by translational motion is solved; then the rough joint compensation function is calculated according to the rough estimated motion parameter and the radial displacement difference; the target echo data is subjected to rough joint compensation through the rough joint compensation function to obtain roughly compensated echo data; the synthetic waveform entropy is solved based on the roughly compensated echo data; then, taking the rough estimated motion parameter as the iterative initial value and the synthetic waveform entropy as the fitness function, iterative optimization is performed through the GWO algorithm to obtain precisely compensated echo data; the direction Fourier transform is performed on the precisely compensated echo data to output the final two-dimensional ISAR image.

[0023] The beneficial effects of the present invention are as follows: Through the two-step compensation strategy of rough compensation and precise compensation, the adverse effects such as frequency offset and defocus caused by target motion can be effectively eliminated, and the compensation accuracy is higher, so that the target can be more accurately focused during final imaging, improving the image resolution and clarity. Among them, in the rough compensation stage, the extracted echo data is used for fractional Fourier transform, which can reduce the data processing time; mutation error elimination and least square fitting are used, which can effectively remove interference data and obtain more accurate rough estimated motion parameters; the rough joint compensation function is calculated through the rough estimated motion parameter and the radial displacement difference, which can initially correct the echo distortion caused by target motion and make the echo data more conducive to subsequent processing. In the precise compensation stage, taking the rough estimated motion parameter as the iterative initial value and the synthetic waveform entropy as the fitness function can not only greatly improve the accuracy of motion parameter estimation, but also limit the parameter search range in precise compensation, accelerate the convergence of the GWO algorithm, have higher calculation efficiency, and can achieve more accurate motion compensation in a shorter time, with better adaptability and robustness. Description of the Drawings

[0024] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the accompanying drawings required for the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can be obtained based on the structures shown in these drawings.

[0025] Figure 1 It is a schematic flowchart of the terahertz ISAR imaging joint compensation method in an embodiment;

[0026] Figure 2 It is a schematic diagram of the satellite scatter point model in an embodiment;

[0027] Figure 3 It is a schematic diagram of the one-dimensional range profile after coarse compensation in an embodiment, where Figure 3 (a) is a schematic diagram of the coarse-compensated one-dimensional range profile, Figure 3 (b) is a partial enlarged schematic diagram;

[0028] Figure 4 It is a schematic comparison diagram between the method proposed in the present invention and the PSO algorithm in an embodiment;

[0029] Figure 5 It is a schematic diagram of the one-dimensional range profile after fine compensation in an embodiment, where Figure 5 (a) is a schematic diagram of the fine-compensated one-dimensional range profile, Figure 5 (b) is a partial enlarged schematic diagram;

[0030] Figure 6 It is a schematic comparison diagram of the two-dimensional ISAR image results in an embodiment, where Figure 6 (a) is a schematic diagram of the image formed by the PSO algorithm, Figure 6 (b) is a schematic diagram of the image formed by the GWO algorithm;

[0031] Figure 7 It is a structural block diagram of the terahertz ISAR imaging joint compensation device in an embodiment;

[0032] Figure 8 It is an internal structure diagram of a computer device in an embodiment.

[0033] The realization of the object, functional features and advantages of the present invention will be further described with reference to the embodiments and the accompanying drawings. Specific embodiments

[0034] Next, in combination with the accompanying drawings in the embodiments of the present invention, the technical solutions in the embodiments of the present invention will be clearly and completely described. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

[0035] It can be understood that the technical solutions between the various embodiments of the present invention can be combined with each other, but it must be based on the ability of those of ordinary skill in the art to implement. When the combination of technical solutions results in contradictions or cannot be implemented, it should be considered that such a combination of technical solutions does not exist and is not within the protection scope required by the present invention.

[0036] Next, in combination with the accompanying drawings in the embodiments of the present invention, the embodiments of the present invention will be described in detail.

[0037] Embodiment 1

[0038] This embodiment discloses a terahertz ISAR imaging joint compensation method. Through a two-step compensation strategy of coarse compensation and fine compensation, it can effectively eliminate the adverse effects such as frequency offset and defocus caused by target movement, and has higher compensation accuracy. As a result, when finally imaging, the target can be focused more accurately, improving the image resolution and clarity. Among them, in the coarse compensation stage, the extracted echo data is used for fractional Fourier transform, which can reduce the data processing time; the mutation error elimination and least squares fitting are used, which can effectively remove the interference data and obtain more accurate rough estimated motion parameters; by calculating the rough joint compensation function through the rough estimated motion parameters and the radial displacement difference, the echo distortion caused by target movement can be initially corrected, making the echo data more conducive to subsequent processing. In the fine compensation stage, taking the rough estimated motion parameters as the iterative initial value and the synthetic waveform entropy as the fitness function, not only can the accuracy of motion parameter estimation be greatly improved, but also the parameter search range in the fine compensation can be restricted, accelerating the convergence of the GWO algorithm, with higher calculation efficiency, and it can achieve more accurate motion compensation in a shorter time, and has better adaptability and robustness.

[0039] Specifically, as Figure 1 shown, the terahertz ISAR imaging joint compensation method provided in this embodiment includes the following steps:

[0040] Step 201, obtain the echo data of the target.

[0041] Step 202, extract the echo data, then perform fractional Fourier transform on the extracted echo data, and then perform mutation error elimination and least squares fitting to obtain the rough estimated motion parameters.

[0042] Step 203: Solve the radial displacement difference between adjacent echoes caused by translational motion based on the roughly estimated motion parameters; then calculate the rough combined compensation function according to the roughly estimated motion parameters and the radial displacement difference.

[0043] Step 204: Perform rough combined compensation on the target echo data through the rough combined compensation function to obtain the roughly compensated echo data.

[0044] Step 205: Solve the synthetic waveform entropy based on the roughly compensated echo data; then use the roughly estimated motion parameters as the initial iteration value and the synthetic waveform entropy as the fitness function, and perform iterative optimization through the GWO algorithm to obtain the precisely compensated echo data.

[0045] Step 206: Perform direction Fourier transform on the precisely compensated echo data and output the final two-dimensional ISAR image.

[0046] In the specific implementation process of Step 201, the Dechirp echo data of the moving target can be expressed as:

[0047] ;

[0048] Let:

[0049] ;

[0050] ;

[0051] ;

[0052] Then the echo data can be simplified to:

[0053] ;

[0054] Furthermore, the second and third terms of the exponential term in the target dechirped received data do not contain fast time , and it is briefly recorded as a constant exponential term . Then the final expression of the echo data is:

[0055] ;

[0056] In the formula, represents the data form after Dechirp processing of a single scatter point and a single echo; represents the data of a single scatter point and a single echo; represents the conjugate form of the reference signal in the Dechirp process; represents the attenuation coefficient of the scatter point echo; represents the amplitude of the transmitted signal; represents the target at The radial velocity of the target relative to the radar at a certain time; Denotes the reference distance during Dechirp processing; Denotes the radar signal wavelength; Denotes the conjugate operation; Denotes the fast time; Denotes the slow time; Denotes the distance from the radar to the target at the moment of each pulse emission; Denotes the pulse repetition period; Denotes the imaginary unit; 、 、 Denotes an intermediate quantity; Denotes the constant exponential term; Denotes the speed of light; Denotes the chirp rate of the transmitted radar signal (LFM signal).

[0057] In the specific implementation process of step 202, in order to improve the operation efficiency, it is only necessary to extract a finite-length LFM signal for fractional Fourier transform (FrFT). Preferably, generally the first half of the echo is extracted for calculation.

[0058] Specifically, a rotation angle is preset; based on the rotation angle the extracted echo data is subjected to fractional Fourier transform along the fast time domain, and the process expression is:

[0059] ;

[0060] Then, according to the relationship between the rotation angle and the chirp rate :

[0061] ;

[0062] ;

[0063] a joint solution is carried out to obtain the radial motion velocity of the target.

[0064] In the formula, denotes the result of fractional Fourier transform; denotes the fractional frequency domain; denotes the normalization factor.

[0065] It can be seen that when a suitable rotation angle is selected for FrFT, the influence of the chirp term is removed, so that a well-focused one-dimensional range image can be obtained.

[0066] Then, to solve the problem of large estimation error of the FrFT method, the method of "eliminating mutation error + least square fitting" is adopted to improve the fitting accuracy of target motion parameters.

[0067] Specifically, by performing mutation error elimination and least square fitting on the radial motion velocity of the target the rough estimated motion parameters are finally obtained, including:

[0068] First, mutation error elimination is performed on the radial motion velocity of the target including:

[0069] For the fractional Fourier transform is performed on the th echo to solve the radial motion velocity of the target at the th echo.

[0070] An empirical value is set. If there exists such that and where then let to complete error elimination.

[0071] Then, least square fitting is performed on the radial motion velocity of the target after error elimination including:

[0072] For the discrete velocity-time observation sequence assuming it follows the parametric model the target minimizes the sum of squared residuals to obtain the rough estimated motion parameters , , .

[0073] In the formula, represents the number of pulses within the coherent integration time (assumed to be even); represents the radial motion velocity of the target at the th echo; represents the radial motion velocity of the target at the th echo; represents the parametric first-order term coefficient with respect to slow time; represents the parametric second-order term coefficient with respect to slow time; represents the estimated velocity at the th echo; represents the rough estimated initial velocity; represents the rough estimated acceleration; represents the rough estimated jerk.

[0074] In the specific implementation process of step 203, based on the rough estimated motion parameters obtained in step 202 , , , let:

[0075] ;

[0076] Based on this, the radial displacement difference between adjacent echoes caused solely by translational motion can be obtained, and the calculation expression is:

[0077] ;

[0078] In the formula, represents the radial displacement difference between adjacent echoes; represents the slow time at the th echo.

[0079] Based on the above analysis, the first compensation function and the second compensation function required for in-pulse compensation of the moving target can be obtained;

[0080] Among them, the expression of the first compensation function is:

[0081] ;

[0082] The expression of the second compensation function is:

[0083] ;

[0084] In the formula, represents the distance from the radar to the target at the moment of each pulse emission; represents the chirp rate of the transmitted radar signal; represents the fast time; represents the radial velocity of the target relative to the radar at ; represents the speed of light;

[0085] Since both in-pulse and inter-pulse compensations rely on the estimation of target motion parameters, therefore, the required joint compensation function is denoted as a whole as , then multiplying the target echo data by the joint compensation function can complete the in-pulse and inter-pulse joint compensation. That is, in the rough compensation stage, based on the rough estimated motion parameters, the radial displacement difference, the first compensation function and the second compensation function, the rough joint compensation function for the target motion parameters is obtained and expressed as .

[0086] In the specific implementation process of step 204, the target echo data is roughly jointly compensated by the rough joint compensation function to obtain the roughly compensated echo data, and the expression is:

[0087] ;

[0088] Among them,

[0089] ;

[0090] In the formula, represents the one-dimensional range profile after combined coarse compensation in intra-pulse and inter-pulse; represents the one-dimensional range profile after intra-pulse coarse compensation; represents the displacement caused by target translational motion; represents the frequency domain; represents the range walk caused by target translational motion; represents the range walk caused by target rotational motion.

[0091] In the specific implementation process of step 205, first determine the fine combined compensation function for target motion parameters.

[0092] According to the coarsely estimated motion parameters , , , obtain the displacement caused by target translational motion, and the expression is:

[0093] ;

[0094] According to the displacement caused by target translational motion, perform range Fourier transform on the coarsely compensated echo data obtained in step 204, which is expressed as:

[0095] ;

[0096] Perform discretization processing on the above formula to obtain:

[0097] ;

[0098] For the convenience of expression, let:

[0099] ;

[0100] Then the simplified expression is:

[0101] ;

[0102] In the formula, represents the single scatter point and single echo data after combined compensation; represents the discretized ; represents the range sampling rate; represents the number of range sampling points; Denote an arbitrary sampling point in the range direction; Denote the pulse frequency of the radar transmitted signal; Denote the sampling pulse in the azimuth direction.

[0103] is the echo data format after intra-pulse and inter-pulse compensation for a single scatterer and single echo. For scatterers, the storage format of the echo data is expressed as:

[0104] ;

[0105] In the formula, Denote scatterers, the echo data within the coherent integration time; Denote the number of received pulses.

[0106] Based on the above analysis, during the fine compensation process, use to denote the one-dimensional range image after translational compensation, where Denote taking the absolute value.

[0107] Express the motion parameters to be estimated during the fine compensation process as:

[0108] ;

[0109] Then during the fine compensation process, the fine joint compensation function can be further written as , and the one-dimensional range image after fine joint compensation can be expressed as .

[0110] Drawing on the concept of waveform entropy, sum the one-dimensional range images of each pulse by rows after synthesis (in actual processing, generally the column direction is the range direction and the row direction is the azimuth direction), so as to obtain the entropy value of the synthesized waveform. If during intra-pulse compensation, the broadening of the main lobe is well corrected, the peak is more prominent, and during inter-pulse compensation, the envelope alignment is well completed, then the peaks will be superimposed on each other, and the synthesized waveform will be more "sharpened", that is, the waveform entropy value is lower; conversely, if the intra-pulse and inter-pulse compensation is not well completed, it will cause the peaks to be staggered from each other, and the envelope of the synthesized waveform will be more "flat", that is, the waveform entropy value is higher. Denote the synthesized waveform as , then should be a column vector, and for an arbitrary sampling point in the range direction, the synthesized waveform has the following relationship:

[0111] ;

[0112] Assume the synthesized range vector , normalize the values of this range vector to , that is, the value of each unit after normalization is , and there is , at this time, it has the property of similar probability, and the synthetic waveform entropy can be expressed as:

[0113] ;

[0114] When solving for any element in it, in order to achieve the best estimation, it is necessary to make:

[0115] ;

[0116] In the formula, represents any sampling point in the range direction; represents the joint compensation function; takes values of , , .

[0117] When using the synthetic waveform entropy as the fitness function, the initial waveform entropy is solved based on the coarsely compensated echo data, and the joint compensation function is the coarse joint compensation function , and the initial waveform entropy is used as the initial fitness function for iteration; in the iterative optimization process, the synthetic waveform entropy is solved based on the finely compensated echo data, and the joint compensation function is the fine joint compensation function .

[0118] After obtaining the coarsely estimated motion parameters and the synthetic waveform entropy, using the coarsely estimated motion parameters as the initial value of iteration and the synthetic waveform entropy as the fitness function, iterative optimization is carried out through the GWO algorithm to obtain the finely compensated echo data.

[0119] It can be understood that in the GWO algorithm, in order to mathematically simulate the social hierarchy of wolves, the optimal solution is defined as for simulating the position of the leading wolf, the sub-optimal solution and the second sub-optimal solution are respectively named and , for simulating the positions of the second-order wolf and the third-order wolf. The remaining candidate solutions are assumed to be , for simulating the position of the subordinate wolf. In the GWO algorithm, hunting (solving) is guided by wolves, wolves and wolves, and wolves follow these three wolves to change their own positions. In this embodiment, the coarsely estimated motion parameters in the coarse compensation stage limit the search range of parameters in the iterative process, accelerating the convergence speed of GWO and having higher compensation efficiency. In addition, iterative optimization through the GWO algorithm is a conventional technical means and will not be elaborated here.

[0120] In the specific implementation process of step 206, before performing the azimuth Fourier transform on the fine-compensated echo data, it is also necessary to perform phase fine calibration PGA on the fine-compensated echo data to correct the phase error and achieve image autofocus.

[0121] In one of the embodiments, a simulation experiment is carried out to verify the method provided by the present invention. As Figure 2 shown, the satellite model used in the experiment consists of 544 scatterers.

[0122] The simulation parameters used are shown in Table 1:

[0123] Table 1 Radar and target simulation parameters

[0124]

[0125] In the experiment, 512 echoes are taken for processing. The one-dimensional range profile of the target after the first-step coarse compensation is as Figure 3 shown.

[0126] It is found that the envelope alignment effect is poor and cannot meet the subsequent imaging requirements. Set the number of gray wolf populations to 30 and the number of iterations to 50 times. Based on the same roughly estimated motion parameters , , Compare the particle swarm optimization algorithm (PSO) with the method proposed by the present invention. In the particle swarm algorithm, the number of particles is also set to 30 and the number of iterations is 50 times, and the search range is controlled the same. The comparison chart of the iterative optimization results of the two algorithms is as Figure 4 shown. It can be seen from the curves of the method proposed by the present invention and PSO that when iterating 5 times, there is a large change in the entropy value of the synthesized waveform of the one-dimensional range profile, but neither of them reaches the optimal value. When iterating to about the 8th time, the PSO algorithm tends to be flat but fails to reach the optimal solution; the method proposed by the present invention begins to tend to be flat at the 20th iteration and approaches the optimal value. It shows that compared with the two algorithms, the method proposed by the present invention can converge faster and has a better final convergence effect.

[0127] Use the fine joint compensation function after 50 iterations to complete the intra-pulse and inter-pulse joint compensation of the target. The one-dimensional range profile is as Figure 5 shown.

[0128] The two-dimensional ISAR imaging results of the echo data after the intra-pulse and inter-pulse compensation of the two algorithms are as Figure 6 shown. At the same time, record the operation time of the two algorithms iterating 50 times and the contrast of the final imaging results as shown in Table 2.

[0129] Table 2 Comparison of algorithm effects

[0130]

[0131] It can be seen that the method proposed by the present invention takes less time, has higher calculation efficiency, and has higher contrast in the two-dimensional ISAR image and better compensation accuracy.

[0132] Although each step in this embodiment Figure 1 is shown in sequence according to the indication of the arrow, these steps do not necessarily need to be executed in the order indicated by the arrow. Unless there is a clear description in the present invention, the execution of these steps has no strict order limit, and these steps can be executed in other orders. Moreover, Figure 1 at least a part of the steps in

[0133] Embodiment 2

[0134] Based on the terahertz ISAR imaging joint compensation method in Embodiment 1, this embodiment discloses a terahertz ISAR imaging joint compensation device. As Figure 7 shown, the terahertz ISAR imaging joint compensation device includes: an echo data acquisition module 401, a rough estimation motion parameter acquisition module 402, a rough joint compensation function acquisition module 403, a rough compensation echo data calculation module 404, a fine compensation echo data calculation 405, and a two-dimensional ISAR image generation module 406, where:

[0135] The echo data acquisition module 401 is used to acquire the echo data of the target.

[0136] The rough estimation motion parameter acquisition module 402 is used to extract the echo data, then perform a fractional Fourier transform on the extracted echo data, and then perform mutation error elimination and least squares fitting to obtain the rough estimation motion parameters.

[0137] The rough joint compensation function acquisition module 403 is used to solve the radial displacement difference between adjacent echoes caused by translational motion based on the rough estimation motion parameters; then calculate the rough joint compensation function according to the rough estimation motion parameters and the radial displacement difference.

[0138] The rough compensation echo data calculation module 404 is used to perform rough joint compensation on the target echo data through the rough joint compensation function to obtain the rough compensation echo data.

[0139] The fine compensation echo data calculation 405 is used to solve the synthetic waveform entropy based on the coarse compensation echo data; then, taking the coarsely estimated motion parameters as the initial iteration values and the synthetic waveform entropy as the fitness function, iterative optimization is performed through the GWO algorithm to obtain the fine compensation echo data.

[0140] The two-dimensional ISAR image generation module 406 is used to perform direction-position Fourier transform on the fine compensation echo data and output the final two-dimensional ISAR image.

[0141] In this embodiment, the specific working processes and working principles of the echo data acquisition module 401, the coarse joint compensation function acquisition module 403, the coarse compensation echo data calculation module 404, the fine compensation echo data calculation 405, and the two-dimensional ISAR image generation module 406 are the same as those in the method of Embodiment 1. Therefore, they will not be elaborated herein. Each of these unit modules can be implemented in whole or in part by software, hardware, and their combination. Each unit module can be embedded in the processor of the computer device in hardware form or be independent of it, or can be stored in the memory of the computer device in software form, so that the processor can call and execute the operations corresponding to each of these unit modules.

[0142] Embodiment 3

[0143] As Figure 8 shown, a terminal device disclosed in this embodiment includes a transmitter, a receiver, a memory, and a processor. Among them, the transmitter is used to send instructions and data, the receiver is used to receive instructions and data, the memory is used to store computer execution instructions, and the processor is used to execute the computer execution instructions stored in the memory to implement the method in Embodiment 1 above.

[0144] It should be noted that the above-mentioned memory can be either independent or integrated with the processor. When the memory is set independently, the terminal device further includes a bus for connecting the memory and the processor.

[0145] Embodiment 4

[0146] This embodiment discloses a computer-readable storage medium in which computer execution instructions are stored. When the processor executes the computer execution instructions, the method in Embodiment 1 above is implemented.

[0147] Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be completed by instructing relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above methods. Among them, any reference to a memory, storage, database, or other medium used in the embodiments provided in the present application can include non-volatile and / or volatile memories. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDR SDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), Rambus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and Rambus dynamic RAM (RDRAM), etc.

[0148] The technical features of the above embodiments can be combined arbitrarily. For the sake of brevity of description, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, it should be considered as the scope described in this specification.

[0149] The above-described embodiments merely represent several implementation manners of the present invention. The description thereof is relatively specific and detailed, but it should not be construed as a limitation on the scope of the invention. It should be noted that for those of ordinary skill in the art, without departing from the concept of the present invention, several modifications and improvements can still be made, and these all belong to the protection scope of the present invention. Therefore, the protection scope of the present invention should be subject to the appended claims.

Claims

1. A terahertz ISAR imaging joint compensation method, characterized in that: The method comprises: Get the target's echo data; Extracting the echo data, then performing fractional Fourier transform on the extracted echo data, and then performing mutation error elimination and least square fitting to obtain a rough estimated motion parameter; Solving the radial displacement difference between adjacent echoes caused by translation based on the coarse estimated motion parameters; then calculating the coarse joint compensation function according to the coarse estimated motion parameters and the radial displacement difference; Performing coarse joint compensation on the target echo data by using the coarse joint compensation function to obtain coarse compensated echo data; Solving the synthetic waveform entropy based on the coarse compensation echo data; then taking the coarse estimated motion parameter as the iteration initial value and the synthetic waveform entropy as the fitness function, performing iterative optimization through the GWO algorithm to obtain the fine compensation echo data; The precision compensated echo data is subjected to directional Fourier transform to output a final two-dimensional ISAR image.

2. The terahertz ISAR imaging joint compensation method according to claim 1, characterized in that: The extracted echo data is subjected to fractional Fourier transform, and then mutation error elimination and least squares fitting are performed to obtain rough estimated motion parameters, including: Preset rotation angle ; Based on the rotation angle The extracted echo data is subjected to fractional Fourier transform in the fast time domain. The process expression is: ； According to the rotation angle Frequency The relationship between: ； ； Perform joint solution to obtain the radial motion speed of the target ; The radial velocity of the target Perform mutation error elimination and least squares fitting to finally obtain rough estimated motion parameters; In the formula, Represents the fractional Fourier transform result; represents the fractional frequency domain; Indicates slow time; Indicates the rotation angle; represents the normalization factor; Indicates the attenuation coefficient of the echo from the scattering point; Indicates the amplitude of the transmitted signal; represents the pulse repetition period; represents an imaginary unit; , , Indicates an intermediate quantity; Indicates the frequency modulation; represents a constant exponential term; Represents the speed of light.

3. The terahertz ISAR imaging joint compensation method according to claim 2, characterized in that: The radial velocity of the target Perform mutation error elimination and least squares fitting to finally obtain rough estimated motion parameters, including: First, the radial velocity of the target Perform mutation error elimination, including: for , for Perform fractional Fourier transform on the echo and solve the The radial velocity of the target at the second echo ; Setting Experience Points , if exists , so that and ,in , then let , complete error elimination; Then, the radial motion speed of the target after error elimination is Perform a least squares fit, including: For discrete speed-time observation sequences , assuming that it obeys the parameterized model , then the objective is to minimize the residual sum of squares , get the rough estimated motion parameters , , ; In the formula, Indicates the number of pulses within the coherent accumulation time; Indicates The radial velocity of the target at the time of the second echo; Indicates The radial velocity of the target at the time of the second echo; represents the coefficient of the first-order term parameterized with respect to slow time; represents the coefficient of the parameterized quadratic term with respect to slow time; Indicates the estimated The speed of the secondary echo; represents a rough estimate of the initial velocity; represents a rough estimate of acceleration; Represents a rough estimate of jerk.

4. The terahertz ISAR imaging joint compensation method according to claim 3, characterized in that: Solving the radial displacement difference between adjacent echoes caused by translation based on the coarse estimated motion parameters includes: set up: ； Based on this solution, the radial displacement difference between adjacent echoes caused by translation is calculated as follows: ； In the formula, Indicates the radial displacement difference between adjacent echoes; Indicates Slow time during the second echo.

5. The terahertz ISAR imaging joint compensation method according to any one of claims 1 to 4, characterized in that: Calculating a rough joint compensation function according to the rough estimated motion parameter and the radial displacement difference includes: First, construct a first-order compensation function and a second-order compensation function; The primary compensation function expression is: ； The quadratic compensation function expression is: ； In the formula, Indicates the distance from the radar to the target at each pulse emission moment; Indicates the modulation rate of the transmitted radar signal; Indicates fast time; Indicates the target The radial velocity relative to the radar at time; represents the speed of light; Then, based on the rough estimated motion parameters, the radial displacement difference, the first compensation function and the second compensation function, a rough joint compensation function for the target motion parameters is obtained. .

6. The terahertz ISAR imaging joint compensation method according to claim 5, characterized in that: The target echo data is coarsely combined and compensated by the coarse joint compensation function to obtain coarse compensated echo data, which is expressed as: ； In the formula, It represents the one-dimensional distance image after the intra-pulse and inter-pulse joint coarse compensation; It represents the one-dimensional range image after intra-pulse coarseness compensation; Represents the displacement caused by the target translation; Indicates the wavelength of radar signal; represents the radial distance movement caused by the target translation; Indicates the radial distance movement caused by the target rotation; Indicates the reference distance during Dechirp processing.

7. The terahertz ISAR imaging joint compensation method according to claim 6, characterized in that: The calculation expression of the synthetic waveform entropy is: ； In the formula, Indicates the distance to any sampling location; represents the joint compensation function.

8. A terahertz ISAR imaging joint compensation device, characterized in that: The device comprises: An echo data acquisition module, used to acquire echo data of a target; A rough estimation motion parameter acquisition module is used to extract the echo data, then perform fractional Fourier transform on the extracted echo data, and then perform mutation error elimination and least square fitting to obtain rough estimation motion parameters; A coarse joint compensation function acquisition module is used to solve the radial displacement difference between adjacent echoes caused by translation based on the coarse estimated motion parameters; and then calculate the coarse joint compensation function according to the coarse estimated motion parameters and the radial displacement difference; A coarse compensation echo data calculation module, used for performing coarse joint compensation on the target echo data by using the coarse joint compensation function to obtain coarse compensation echo data; Fine compensation echo data calculation is used to solve the synthetic waveform entropy based on the coarse compensation echo data; then the coarse estimated motion parameter is used as the iteration initial value, the synthetic waveform entropy is used as the fitness function, and the GWO algorithm is used to perform iterative optimization and solve to obtain fine compensation echo data; The two-dimensional ISAR image generation module is used to perform directional Fourier transform on the precision compensated echo data and output a final two-dimensional ISAR image.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the terahertz ISAR imaging joint compensation method according to any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the terahertz ISAR imaging joint compensation method according to any one of claims 1 to 7 are implemented.

Images