A Novel Calibration Method for Signal Leakage between Multiple Channels in Array Imaging

: In this paper, the array imaging process with the presence of signal leakage between channels is derived and analyzed in short-range imaging applications. The echo wave of the system can be divided into two parts: the echo signals we expected

Signal leakage is an important problem which needs to be processed in a radar system.Signal leakage that occurs between transmitter and receiver can cause a strong degradation of the detection sensitivity as a result of the phase noise of the transmit oscillator.Its calibration method has been widely discussed.There are several ways to deal with this leakage.One way is to increase the distance between the Tx and the Rx to reduce the leakage [22].The second way is to add a closed loop leakage canceller which can generate an error vector including leakage amplitude and phase information, and this information can be used to carry out the calibration [23][24][25].
Signal leakage between different channels in multichannel systems is a unique problem, especially in imaging systems with requirements of high resolution and sensitivity.The leakage that occurs between the transmitter and the receiver has been widely studied.There are a lot of discussions about its signal model and processing method [26][27][28].However the impact of the leakage that occurs between transmitters or receivers has rarely been discussed.In this paper, the array imaging process with the presence of signal leakage caused by the lower isolation between transmitters or receivers is researched.First, a common array structure and its signal leakage model between different channels are introduced.Second, the imaging process with the presence of signal leakage is derived based on the model and the influence of the signal leakage is discussed.Signals leaked to adjacent channels can lead to ghosts which can deeply deteriorate the system cross resolution.To eliminate the bad influence of the signal leakages, a calibration method operating in the beam domain is proposed next.Finally, to validate the effectiveness of the proposed method, a simulation and some experiments are carried out.All the derivation and analysis are based on a common array structure applied in short range imaging.
The paper is structured as follows.Section 2 introduces the used array structure and signal leakage model between different channels.The impacts of the channel signal leakages are analyzed in Section 3. Based on the analysis, a calibration method is prosed in Section 4. Finally, in Section 5, we verify the calibration method with a simulation and some experiments.

Array Structure and Signal Leakage Model
Figure 1 is a common antenna array structure of the millimeter wave in short range imaging [29].The architecture of the array consists of two linear arrays.Suppose that the transmitting array is on the above and the receiving array is on the below.The interval between the transmitting or receiving elements is λ.In order to get equivalent phase centers with an interval of λ/2, the receiving array is shifted λ/2 along the x direction with respect to the transmitting array.When one transmitting element is working, its two adjacent receiving elements carry out the signal reception in sequence.The signal is usually regarded as though it is transmitted and received both at the midpoint of the line connecting the corresponding working antennas.Thus, a set of equivalent phase centers spaced λ/2 can be obtained.
Electronics 2020, 9, x FOR PEER REVIEW 2 of 13 transmitters or receivers is researched.First, a common array structure and its signal leakage model between different channels are introduced.Second, the imaging process with the presence of signal leakage is derived based on the model and the influence of the signal leakage is discussed.Signals leaked to adjacent channels can lead to ghosts which can deeply deteriorate the system cross resolution.To eliminate the bad influence of the signal leakages, a calibration method operating in the beam domain is proposed next.Finally, to validate the effectiveness of the proposed method, a simulation and some experiments are carried out.All the derivation and analysis are based on a common array structure applied in short range imaging.The paper is structured as follows.Section 2 introduces the used array structure and signal leakage model between different channels.The impacts of the channel signal leakages are analyzed in Section 3. Based on the analysis, a calibration method is prosed in Section 4. Finally, in Section 5, we verify the calibration method with a simulation and some experiments.

Array Structure and Signal Leakage Model
Figure 1 is a common antenna array structure of the millimeter wave in short range imaging [29].The architecture of the array consists of two linear arrays.Suppose that the transmitting array is on the above and the receiving array is on the below.The interval between the transmitting or receiving elements is λ.In order to get equivalent phase centers with an interval of λ/2, the receiving array is shifted λ/2 along the x direction with respect to the transmitting array.When one transmitting element is working, its two adjacent receiving elements carry out the signal reception in sequence.The signal is usually regarded as though it is transmitted and received both at the midpoint of the line connecting the corresponding working antennas.Thus, a set of equivalent phase centers spaced λ/2 can be obtained.In order to implement the signal's independent emission and reception of the array units described in Figure 1, the imaging system usually adopts a multistage switch matrix.Figure 2 is a common switch network structure.For the convenience of description, a three-stage switch matrix is used in Figure 2. The first, second, and third stage of the transmitting network adopt SPTK, SPTM, and SPTN switch chips, respectively, as shown in Figure 2. The structure of the receiving switch network is the same as the transmitting.Therefore, this architecture can implement a 1 × K × M × N switch network.In order to implement the signal's independent emission and reception of the array units described in Figure 1, the imaging system usually adopts a multistage switch matrix.Figure 2 is a common switch network structure.For the convenience of description, a three-stage switch matrix is used in Figure 2. The first, second, and third stage of the transmitting network adopt SPTK, SPTM, and SPTN switch chips, respectively, as shown in Figure 2. The structure of the receiving switch network is the same as the transmitting.Therefore, this architecture can implement a 1 × K × M × N switch network.
Electronics 2020, 9, x FOR PEER REVIEW 2 of 13 transmitters or receivers is researched.First, a common array structure and its signal leakage model between different channels are introduced.Second, the imaging process with the presence of signal leakage is derived based on the model and the influence of the signal leakage is discussed.Signals leaked to adjacent channels can lead to ghosts which can deeply deteriorate the system cross resolution.To eliminate the bad influence of the signal leakages, a calibration method operating in the beam domain is proposed next.Finally, to validate the effectiveness of the proposed method, a simulation and some experiments are carried out.All the derivation and analysis are based on a common array structure applied in short range imaging.The paper is structured as follows.Section 2 introduces the used array structure and signal leakage model between different channels.The impacts of the channel signal leakages are analyzed in Section 3. Based on the analysis, a calibration method is prosed in Section 4. Finally, in Section 5, we verify the calibration method with a simulation and some experiments.

Array Structure and Signal Leakage Model
Figure 1 is a common antenna array structure of the millimeter wave in short range imaging [29].The architecture of the array consists of two linear arrays.Suppose that the transmitting array is on the above and the receiving array is on the below.The interval between the transmitting or receiving elements is λ.In order to get equivalent phase centers with an interval of λ/2, the receiving array is shifted λ/2 along the x direction with respect to the transmitting array.When one transmitting element is working, its two adjacent receiving elements carry out the signal reception in sequence.The signal is usually regarded as though it is transmitted and received both at the midpoint of the line connecting the corresponding working antennas.Thus, a set of equivalent phase centers spaced λ/2 can be obtained.In order to implement the signal's independent emission and reception of the array units described in Figure 1, the imaging system usually adopts a multistage switch matrix.Figure 2 is a common switch network structure.For the convenience of description, a three-stage switch matrix is used in Figure 2. The first, second, and third stage of the transmitting network adopt SPTK, SPTM, and SPTN switch chips, respectively, as shown in Figure 2. The structure of the receiving switch network is the same as the transmitting.Therefore, this architecture can implement a 1 × K × M × N switch network.Based on the structure, we can make a simple analysis about the signal leakages occurred in the switch network due to the lower isolation of the switch chips when the signal is transmitted and received.Suppose that the red path in Figure 2 is the signal path we expected; therefore, the corresponding switch channels are opening while the other channels on the same chip are closed at the same time.As the isolation of the switch is limited, some signals are leaked to the closed channels, indicated by the green lines in Figure 2. The signal leakages occur in every switch stage, but those which do not occur at the last stage do not need to be considered because they have been attenuated at least twice when they arrive at the antenna and their impacts to the imaging result are not significant.Thus, only the signal leakages that occurred at the same switch chip with the opened channel at the last stage need to be considered.That means only the N − 1 channels' signal leakages need to be considered on the transmitting switch array or the receiving array at a certain moment.
Because the signal leaked by the transmitting array and the receiving array in the received signal at a certain time is linearly superimposed, the analysis method in this paper, such as Fourier transform, has linear properties.Therefore, in order to simplify the analysis, we can independently analyze the signal leakage's influence of the transmitting array and the receiving array.Because the two analysis processes are similar, only the signal leakage of the transmitting array is considered, and an assumption is made that signal leakage does not exist between the receiving switch matrixes.Under this assumption, for the signal transmitted from the expected transmitting antenna and the signals leaked from the adjacent antenna, there are two sets of paths to the receiving antenna, as Figure 3 shows.In one set of paths (indicated by green and red lines in Figure 3), the signals illuminate the targets and then are reflected to the receiving antenna.In the other set of paths (indicated by yellow line in Figure 3), the signal is directly leaked to the receiving antenna without the reflection of the targets, which could be called the Tx/Rx coupling interference.The Tx/Rx coupling interference has a huge impact on the imaging quality.It usually has a higher power level compared to the target echo wave such that the phase noise occurred with it sometimes may raise up the system noise floor.It will deteriorate the dynamic range of the imaging system.However, the noise could be cancelled mostly by using the empty background noise.As such, the Tx/Rx coupling interference would not be concerned in this paper.Only the signals (original signal or the leaked signals) reflected by the targets are concerned.
Electronics 2020, 9, x FOR PEER REVIEW 3 of 13 Based on the structure, we can make a simple analysis about the signal leakages occurred in the switch network due to the lower isolation of the switch chips when the signal is transmitted and received.Suppose that the red path in Figure 2 is the signal path we expected; therefore, the corresponding switch channels are opening while the other channels on the same chip are closed at the same time.As the isolation of the switch is limited, some signals are leaked to the closed channels, indicated by the green lines in Figure 2. The signal leakages occur in every switch stage, but those which do not occur at the last stage do not need to be considered because they have been attenuated at least twice when they arrive at the antenna and their impacts to the imaging result are not significant.Thus, only the signal leakages that occurred at the same switch chip with the opened channel at the last stage need to be considered.That means only the N − 1 channels' signal leakages need to be considered on the transmitting switch array or the receiving array at a certain moment.
Because the signal leaked by the transmitting array and the receiving array in the received signal at a certain time is linearly superimposed, the analysis method in this paper, such as Fourier transform, has linear properties.Therefore, in order to simplify the analysis, we can independently analyze the signal leakage's influence of the transmitting array and the receiving array.Because the two analysis processes are similar, only the signal leakage of the transmitting array is considered, and an assumption is made that signal leakage does not exist between the receiving switch matrixes.Under this assumption, for the signal transmitted from the expected transmitting antenna and the signals leaked from the adjacent antenna, there are two sets of paths to the receiving antenna, as Figure 3 shows.In one set of paths (indicated by green and red lines in Figure 3), the signals illuminate the targets and then are reflected to the receiving antenna.In the other set of paths (indicated by yellow line in Figure 3), the signal is directly leaked to the receiving antenna without the reflection of the targets, which could be called the Tx/Rx coupling interference.The Tx/Rx coupling interference has a huge impact on the imaging quality.It usually has a higher power level compared to the target echo wave such that the phase noise occurred with it sometimes may raise up the system noise floor.It will deteriorate the dynamic range of the imaging system.However, the noise could be cancelled mostly by using the empty background noise.As such, the Tx/Rx coupling interference would not be concerned in this paper.Only the signals (original signal or the leaked signals) reflected by the targets are concerned.Suppose that the nth transmitting element's position is (xn Tx , 0, za) where nϵ{1, 2, …, K × M × N}.The mth receiving element's position is (xm Rx , 0, −za) where mϵ{1, 2, …, K × M × N}.The target position is (x, y, 0).The flight trip from the nth transmitting element to the target is: and from the target to the mth, the receiving element is: and from the target to the mth, the receiving element is: Then, the echo wave when the nth transmitting element and the mth receiving antenna is working can be expressed as: where: • N i f rem(n, N) 0 which represents the transmitting element' position set of the leaked signals.
[A] represents the nearest integer less than or equal to A. Rem(n,N) represents the remainder after division of n by N. λ is the wavelength.k is the wavenumber 2π/λ.η is the attenuation coefficient of the leaked signal compared to the original signal (suppose that the original signal amplitude is 1).σ is the radar scattering coefficient (RCS) of the target.The echo wave consists of two parts.The first term represents the echo signal we expected, and the second term represents the echo signals that result from the leaked signals.
It is difficult to make a fast imaging calculation according to Equation (1).If the distance between the transmitting element and the receiving element is very close, then some phase center approximation can be made according to the monostatic measurement method developed by Kell [30] and we can approximate that the signal is transmitted and received at the same place at the equivalent phase center located at the midpoint of the transmitting and receiving antenna connection.The echo wave can be expressed as: According to Equation ( 1), the elements of the set Φ are integers in the range of n − N + 1 to n + N − 1 according to the configuration of the array described earlier.However, there are only N − 1 elements in the set Φ, not 2N − 1.As such, the Φ's elements cannot traverse all of these numbers from n − N + 1 to n + N − 1 at a specific transmitting and receiving position.Its elements change with the transmitting element's position.Therefore, the second term of echo Equation (2) varies with the transmitting element position.It is difficult to make a direct and clear analysis by using it.The first term can be regarded as a special case of the second term when I = n and η = 1, so we can merge the first and the second items.Based on the above reasons, we rewrite the echo wave Expression (2) as follows: where Set that (x n Tx + x m Rx )/2 = x a , considering that the transmitting interval is λ according to the system model shown in Figure 1, then we can get ( Using this relation, the echo Expression (3) can be rewritten to: In this section, a common array structure applied in short range imaging is introduced.Its signal model with the presence of signal leakage is constructed as Expression (4).Based on the expression, we can do further analysis about the impacts of the signal leakage on array imaging.

The Analysis of the Signal Leakage Impact on Array Imaging
Image reconstruction can be made based on Expression (4).Firstly, one-dimensional Fourier transform is implemented along the variable x a .η i and σ would not change with n, so we can obtain the wavenumber distribution of the echo data as follows: a periodic rectangular pulse function of receiving number n due to the special switch matrix configuration of the array and its value is between 0 and 1.Its spectrum U i is shown in Figure 4, which can be expressed as: where j = 0, ±1, ±2, . . ., Sa(x) = sin(x)/x, τ is the pulse width of u i (n), T is the period of the u i (n).U i can be seen as some δ-function combination with the envelope of Sa(jπτ/T).
Electronics 2020, 9, x FOR PEER REVIEW 5 of 13 In this section, a common array structure applied in short range imaging is introduced.Its signal model with the presence of signal leakage is constructed as Expression (4).Based on the expression, we can do further analysis about the impacts of the signal leakage on array imaging.

The Analysis of the Signal Leakage Impact on Array Imaging
Image reconstruction can be made based on Expression (4).Firstly, one-dimensional Fourier transform is implemented along the variable xa.ηi and σ would not change with n, so we can obtain the wavenumber distribution of the echo data as follows: Ui is the Fourier transform of ui(n).When i = n, un = 1, ui(n) is a constant, its spectrum is a δ-function.When i ≠ n, ui(n) is a periodic rectangular pulse function of receiving number n due to the special switch matrix configuration of the array and its value is between 0 and 1.Its spectrum Ui is shown in Figure 4, which can be expressed as:  where j = 0, ±1, ±2, …, Sa(x) = sin(x)/x, τ is the pulse width of ui(n), T is the period of the ui(n).Ui can be seen as some δ-function combination with the envelope of Sa(jπτ/T).We can see in Figure 4 that the power of other δ-function (δ(kx − Δkj), Δkj ≠ 0) is weaker compared to the δ(kx) because the value of Sa(kx) decreases sharply as kx increases.It would not have an important influence on the image.δ-functions (δ(kx − Δkj), Δkj ≠ 0) will be ignored in the following derivation, namely Ui ≈ δ(kx).γi is the ratio of the amplitude of δ(kx) in Ui to the δ(kx) in Un.Obviously, γn = 1.The convolution result of δ(kx − Δkj) with A(kx) is still A(kx − Δkj).Then, the wavenumber distribution (5) can be simplified as follows: 6) can be simplified as: It is easy to find that the phase of Sif(ky,kx) is a linear combination of variable kx and ky.Thus, applying a 2D-FFT or 2D-IFFT on it, the image result will be obtained.One thing to notice is that the map from k to ky is not a linear process, therefore, the ky directly achieved form the uniformly sampled k is not equally spaced.Before the 2D-IFFT, a process called Stolt interpolation [31,32] needs We can see in Figure 4 that the power of other δ-function (δ(k x − ∆k j ), ∆k j 0) is weaker compared to the δ(k x ) because the value of Sa(k x ) decreases sharply as k x increases.It would not have an important influence on the image.δ-functions (δ(k x − ∆k j ), ∆k j 0) will be ignored in the following derivation, namely U i ≈ δ(k x ).γ i is the ratio of the amplitude of δ(k x ) in U i to the δ(k x ) in U n .Obviously, γ n = 1.The convolution result of δ(k x − ∆k j ) with A(k x ) is still A(k x − ∆k j ).Then, the wavenumber distribution (5) can be simplified as follows: Set k y = 4k 2 − k x 2 , Expression (6) can be simplified as: It is easy to find that the phase of S if (k y ,k x ) is a linear combination of variable k x and k y .Thus, applying a 2D-FFT or 2D-IFFT on it, the image result will be obtained.One thing to notice is that the map from k to k y is not a linear process, therefore, the k y directly achieved form the uniformly sampled k is not equally spaced.Before the 2D-IFFT, a process called Stolt interpolation [31,32] needs to be carried out which will not be discussed in this paper.Suppose that the band of the k x and k y is B x and B y , respectively.Then, the final image result is obtained as follows after 2D-IFFT: Set that which can be seen as the ideal point spread function of the system.The P original (x,y) can be rewritten as: It is obvious that the original point spread function P original (x,y) influenced by the leaked signals is a combination of the P(x,y) shifted (i − n)λ/2 in the x direction.This can deteriorate the system cross resolution.It is easy to imagine that the superposition of several pulses is usually wider than a single pulse.The larger the coefficient γ i , the more serious the deterioration of the cross resolution is.On the contrary, the leaked signals have little impact on the range resolution according to the P original (x,y) expression.
In this section, the imaging process is derived.The final imaging result is conveyed as Expression (10).The influence of the signal leakage proposed above is analyzed.It can deteriorate the system cross resolution but have little impact on the range resolution.

Calibration Method of the Signal Leakage
To eliminate the bad influences of the signal leakages, a calibration method will be proposed next.The shift in the time domain can be compensated in the frequency domain.Doing 1-D FFT along the x direction on both side of the Expression (10), the following expression can be obtained: Then, the ideal PSF can be obtained: The denominator in the IFFT x parenthesis can be regarded as a calibration factor in k x domain.We call it array calibration factor.It makes the ghost caused by signal leakages shift to the original position.The factor can be obtained according to the configuration of the switch array.When the calibration factor is ready, the calibration can be implemented according to Expression (11) in the frequency domain.Since the calibration factor is determined by the structure of the array, this calibration method can be named array calibration factor method (ACFM).
The calibration process can be carried out according to the following steps: Step 1: Measure the attenuation coefficient η of the antenna array, then the coefficient η i can be obtained.
Step 2: Calculate the ratio γ i according to Expression (6).It is usually difficult to obtain its analytical form; however, we can get its numerical solution through some simulation according to the application scenario of the antenna array.
Step 3: Calculate the array calibration factor Step 2.
Step 4: Calculate the calibration result P(x,y) by using P original (x,y) according to Expression (11).
Because the structure and isolation of the antenna array usually do not change, coefficient η and ratio γ i are constants.The calibration factor in Step 3 only needs to be calculated once.The computational complexity of the calibration method proposed above is only related to Step 4. If the dimension of P(x,y) is N x × N y , the computational complexity of Step 4 is O(N x, N y ).It has no significant effect on the imaging speed.
In this section, a calibration method operating in the beam domain is proposed.Its calibration steps and computational complexity are introduced.It has no significant effect on imaging speed.

Simulation Results
Suppose that the configuration of the switch array is K = 5, M = 4, and N = 4.The operating frequency of the system is 27-33 GHz and the isolation of the switch is 10 dB.Then, the switch array can match with a 1 × 80 Tx's antenna array and a 1 × 80 Rx's antenna array.The point spread function of the array is simulated.The imaging results are shown in Figures 5 and 6. Figure 5a,b are the PSF when the signal leakages between different channels don't exist and exist, respectively.Figure 6 shows the comparison of the cross range PSF and the comparison of the range PSF in the two different cases, respectively.It is obvious that the signal leakages between different channels make the PSF broaden in cross range direction but it has little effect on the PSF in range direction.According to Expression (10), the PSF P original (x,y) influenced by leaked signals is a combination of the P(x,y) shifted (i − n)λ/2 in the x direction but not shift in the y direction.The result of the superposition of several pulses with misaligned centers will be wider than the pulse before superposition, thus, the system cross resolution along the x axis is deteriorated and the range resolution along y axis is not affected.We can see that in Figure 6, the cross resolution of signal leakage is worse, and the range resolution does not change compared to the results without signal leakage.It is consistent with the analysis above.
Electronics 2020, 9, x FOR PEER REVIEW 7 of 13 Step 3: Calculate the array calibration factor  by using result of Step 1 and Step 2.
Step 4: Calculate the calibration result P(x,y) by using Poriginal(x,y) according to Expression (11).
Because the structure and isolation of the antenna array usually do not change, coefficient η and ratio γi are constants.The calibration factor in Step 3 only needs to be calculated once.The computational complexity of the calibration method proposed above is only related to Step 4. If the dimension of P(x,y) is Nx × Ny, the computational complexity of Step 4 is O(Nx,Ny).It has no significant effect on the imaging speed.
In this section, a calibration method operating in the beam domain is proposed.Its calibration steps and computational complexity are introduced.It has no significant effect on imaging speed.

Simulation Results
Suppose that the configuration of the switch array is K = 5, M = 4, and N = 4.The operating frequency of the system is 27-33 GHz and the isolation of the switch is 10 dB.Then, the switch array can match with a 1 × 80 Tx's antenna array and a 1 × 80 Rx's antenna array.The point spread function of the array is simulated.The imaging results are shown in Figures 5 and 6. Figure 5a,b are the PSF when the signal leakages between different channels don't exist and exist, respectively.Figure 6 shows the comparison of the cross range PSF and the comparison of the range PSF in the two different cases, respectively.It is obvious that the signal leakages between different channels make the PSF broaden in cross range direction but it has little effect on the PSF in range direction.According to Expression (10), the PSF Poriginal(x,y) influenced by leaked signals is a combination of the P(x,y) shifted (i − n)λ/2 in the x direction but not shift in the y direction.The result of the superposition of several pulses with misaligned centers will be wider than the pulse before superposition, thus, the system cross resolution along the x axis is deteriorated and the range resolution along y axis is not affected.We can see that in Figure 6, the cross resolution of signal leakage is worse, and the range resolution does not change compared to the results without signal leakage.It is consistent with the analysis above.To eliminate the influence of the signal leakage, a calibration was made according to Expression (11).The array calibration factor is shown in Figure 7a and the calibration results are shown in Figure 7b and Figure 8.The imaging result after calibration is basically the same as the result without the signal leakages.It indicates that the calibration method can eliminate the effects of leakage signals.On the other hand, it can be found that the calibration factor is similar to a window function.Thus, the calibration process can be seen as an inverse process of windowing in the frequency domain.We can use the window function to replace the array calibration factor.We call this method To eliminate the influence of the signal leakage, a calibration was made according to Expression (11).The array calibration factor is shown in Figure 7a and the calibration results are shown in Figures 7b  and 8.The imaging result after calibration is basically the same as the result without the signal leakages.It indicates that the calibration method can eliminate the effects of leakage signals.To eliminate the influence of the signal leakage, a calibration was made according to Expression (11).The array calibration factor is shown in Figure 7a and the calibration results are shown in Figure 7b and Figure 8.The imaging result after calibration is basically the same as the result without the signal leakages.It indicates that the calibration method can eliminate the effects of leakage signals.On the other hand, it can be found that the calibration factor is similar to a window function.Thus, the calibration process can be seen as an inverse process of windowing in the frequency domain.We can use the window function to replace the array calibration factor.We call this method  To eliminate the influence of the signal leakage, a calibration was made according to Expression (11).The array calibration factor is shown in Figure 7a and the calibration results are shown in Figure 7b and Figure 8.The imaging result after calibration is basically the same as the result without the signal leakages.It indicates that the calibration method can eliminate the effects of leakage signals.On the other hand, it can be found that the calibration factor is similar to a window function.Thus, the calibration process can be seen as an inverse process of windowing in the frequency domain.We can use the window function to replace the array calibration factor.We call this method  On the other hand, it can be found that the calibration factor is similar to a window function.Thus, the calibration process can be seen as an inverse process of windowing in the frequency domain.We can use the window function to replace the array calibration factor.We call this method the window function calibration method (WFCM).A Kaiser window with beta of 2.5 created by MATLAB shown in Figure 9a is used and the calibration result is shown in Figures 9b and 10.
the window function calibration method (WFCM).A Kaiser window with beta of 2.5 created by MATLAB shown in Figure 9a is used and the calibration result is shown in Figure 9b and Figure 10.It can be seen that both methods can eliminate the broadening of the main lobe in the cross range.However, the side lobes raise evidently when using the WFCM method at the cross direction compared to the ACFM method described in Expression (11).It indicates that the array calibration factor created by Expression (11) can improve the image quality compared to the general window function.Impulse response width (IRW) and peak side lobe ratio (PSLR) [33] can be used to evaluate the calibration performance.The comparison of two calibration method performance is shown in Table 1.IRW is normalized with its theoretical resolution, so its unit is 1.The measurement results of PSLR are also marked in Figure 10.

Experiment Results
An experiment was conducted to verify the effectiveness of the method mentioned above.The real array, targets, and photo of measurements in the experiments are shown in Figure 11.The configuration of the switch array is K = 5, M = 4, and N = 4.The operating frequency of the system is 27-33 GHz.The isolation of the array between different channels is about 25 dB.A two-dimension (2D) imaging was carried out with a ball shown in Figure 11, which can be seen as a point target.The the window function calibration method (WFCM).A Kaiser window with beta of 2.5 created by MATLAB shown in Figure 9a is used and the calibration result is shown in Figure 9b and Figure 10. ( It can be seen that both methods can eliminate the broadening of the main lobe in the cross range.However, the side lobes raise evidently when using the WFCM method at the cross direction compared to the ACFM method described in Expression (11).It indicates that the array calibration factor created by Expression (11) can improve the image quality compared to the general window function.Impulse response width (IRW) and peak side lobe ratio (PSLR) [33] can be used to evaluate the calibration performance.The comparison of two calibration method performance is shown in Table 1.IRW is normalized with its theoretical resolution, so its unit is 1.The measurement results of PSLR are also marked in Figure 10.

Experiment Results
An experiment was conducted to verify the effectiveness of the method mentioned above.The real array, targets, and photo of measurements in the experiments are shown in Figure 11.The configuration of the switch array is K = 5, M = 4, and N = 4.The operating frequency of the system is 27-33 GHz.The isolation of the array between different channels is about 25 dB.A two-dimension (2D) imaging was carried out with a ball shown in Figure 11, which can be seen as a point target.The It can be seen that both methods can eliminate the broadening of the main lobe in the cross range.However, the side lobes raise evidently when using the WFCM method at the cross direction compared to the ACFM method described in Expression (11).It indicates that the array calibration factor created by Expression (11) can improve the image quality compared to the general window function.Impulse response width (IRW) and peak side lobe ratio (PSLR) [33] can be used to evaluate the calibration performance.The comparison of two calibration method performance is shown in Table 1.IRW is normalized with its theoretical resolution, so its unit is 1.The measurement results of PSLR are also marked in Figure 10.

Experiment Results
An experiment was conducted to verify the effectiveness of the method mentioned above.The real array, targets, and photo of measurements in the experiments are shown in Figure 11.The configuration of the switch array is K = 5, M = 4, and N = 4.The operating frequency of the system is 27-33 GHz.The isolation of the array between different channels is about 25 dB.A two-dimension (2D) imaging was carried out with a ball shown in Figure 11, which can be seen as a point target.The results are shown in Figure 12. Figure 12a is the 2D imaging result without calibration.Figure 12b is the imaging result after calibration.Figure 13a,b are the PSF comparisons between the results before and after calibration.Obviously, it can be seen that the calibration process has an improvement on the cross-range resolution and has no impact on the range resolution.One point to note is that the calibration process improves the cross resolution and also raises the sidelobe levels.However, a set of calibration factors which is designed according to the structure of the array can bring as little side-lobe increases as possible compared to other window functions.
Electronics 2020, 9, x FOR PEER REVIEW 10 of 13 results are shown in Figure 12. Figure 12a is the 2D imaging result without calibration.Figure 12b is the imaging result after calibration.Figure 13a,b are the PSF comparisons between the results before and after calibration.Obviously, it can be seen that the calibration process has an improvement on the cross-range resolution and has no impact on the range resolution.One point to note is that the calibration process improves the cross resolution and also raises the sidelobe levels.However, a set of calibration factors which is designed according to the structure of the array can bring as little side-lobe increases as possible compared to other window functions.A simple resolution test board shown in Figure 11 was tested with three-dimension (3D) imaging.The maximum projection results are shown in Figure 14. Figure 14a shows the result results are shown in Figure 12. Figure 12a is the 2D imaging result without calibration.Figure 12b is the imaging result after calibration.Figure 13a,b are the PSF comparisons between the results before and after calibration.Obviously, it can be seen that the calibration process has an improvement on the cross-range resolution and has no impact on the range resolution.One point to note is that the calibration process improves the cross resolution and also raises the sidelobe levels.However, a set of calibration factors which is designed according to the structure of the array can bring as little side-lobe increases as possible compared to other window functions.A simple resolution test board shown in Figure 11 was tested with three-dimension (3D) imaging.The maximum projection results are shown in Figure 14. Figure 14a shows the result results are shown in Figure 12. Figure 12a is the 2D imaging result without calibration.Figure 12b is the imaging result after calibration.Figure 13a,b are the PSF comparisons between the results before and after calibration.Obviously, it can be seen that the calibration process has an improvement on the cross-range resolution and has no impact on the range resolution.One point to note is that the calibration process improves the cross resolution and also raises the sidelobe levels.However, a set of calibration factors which is designed according to the structure of the array can bring as little side-lobe increases as possible compared to other window functions.A simple resolution test board shown in Figure 11 was tested with three-dimension (3D) imaging.The maximum projection results are shown in Figure 14. Figure 14a shows the result A simple resolution test board shown in Figure 11 was tested with three-dimension (3D) imaging.The maximum projection results are shown in Figure 14. Figure 14a shows the result without calibration.
Figure 14b shows the result calibrated with the method mentioned above.We can see that the gap between the aluminum bars in Figure 14b is clearer than the gap in Figure 14a.This shows that the calibration method can eliminate most of the impacts of the signal leakage and improve the cross resolution of the system.
Electronics 2020, 9, x FOR PEER REVIEW 11 of 13 without calibration.Figure 14b shows the result calibrated with the method mentioned above.We can see that the gap between the aluminum bars in Figure 14b is clearer than the gap in Figure 14a.This shows that the calibration method can eliminate most of the impacts of the signal leakage and improve the cross resolution of the system.In this section, a simulation and some experiments are carried out to validate the effectiveness of the proposed method.The simulation and the experiment results indicate that this calibration method can mostly cancel the influence of the signal leakages.

Conclusions
This paper discussed the influence of the signal leakages between different channels in the application of array imaging.Signal leakage between different channels can deteriorate the system cross resolution but has no impact on the range resolution.When the channel isolation is not good enough, it will lead to the appearance of grating lobes, which we tried to eliminate in this paper.A calibration method was presented to eliminate the influences of signal leakages and was validated with the simulation and some experiments.The calibration factor proposed in this paper can better improve the image quality compared to the used window function.However, we have not yet obtained an optimal window function selection method which may obtain better calibration results.Some optimization methods can be tried in related directions.It will be carried out in future work.The effect of noise on the calibration process proposed in the paper also can be a research direction in the next work.The same analysis method and calibration method can be extended to complex array imaging application analyses.In this section, a simulation and some experiments are carried out to validate the effectiveness of the proposed method.The simulation and the experiment results indicate that this calibration method can mostly cancel the influence of the signal leakages.

Conclusions
This paper discussed the influence of the signal leakages between different channels in the application of array imaging.Signal leakage between different channels can deteriorate the system cross resolution but has no impact on the range resolution.When the channel isolation is not good enough, it will lead to the appearance of grating lobes, which we tried to eliminate in this paper.A calibration method was presented to eliminate the influences of signal leakages and was validated with the simulation and some experiments.The calibration factor proposed in this paper can better improve the image quality compared to the used window function.However, we have not yet obtained an optimal window function selection method which may obtain better calibration results.Some optimization methods can be tried in related directions.It will be carried out in future work.The effect of noise on the calibration process proposed in the paper can be a research direction in the next work.The same analysis method and calibration method can be extended to complex array imaging application analyses.

Figure 1 .
Figure 1.Array configuration of the transmitting and receiving elements.

Figure 1 .
Figure 1.Array configuration of the transmitting and receiving elements.

Figure 1 .
Figure 1.Array configuration of the transmitting and receiving elements.

Figure 3 .
Figure 3. Signal path form Tx to Rx.

Figure 4 .
Figure 4.The Fourier transform of u i (n).

Figure 5 .Figure 5 .Figure 6 .
Figure 5. Simulated PSF in two different cases, (a) is the PSF when the signal leakages do not exist and (b) is the PSF when leakages exist.

Figure 7 .Figure 8 .
Figure 7.The array calibration factor and the imaging result after calibration, (a) is the array calibration factor and (b) is the imaging result after calibration.

Figure 6 .
Figure 6.The comparison of PSF in two different dimensions, (a) is the cross range PSF and (b) is the range PSF.

Figure 6 .
Figure 6.The comparison of PSF in two different dimensions, (a) is the cross range PSF and (b) is the range PSF.

Figure 7 .Figure 8 .
Figure 7.The array calibration factor and the imaging result after calibration, (a) is the array calibration factor and (b) is the imaging result after calibration.

Figure 7 .Figure 6 .
Figure 7.The array calibration factor and the imaging result after calibration, (a) is the array calibration factor and (b) is the imaging result after calibration.

Figure 7 .Figure 8 .
Figure 7.The array calibration factor and the imaging result after calibration, (a) is the array calibration factor and (b) is the imaging result after calibration.

Figure 8 .
Figure 8.The PSF comparison of three different cases, (a) is the cross range PSF and (b) is the range PSF.

Figure 9 .Figure 10 .
Figure 9.The Kaiser window with beta of 2.5 and its calibration result, (a) is Kaiser window and (b) is the calibration result.

Figure 9 .
Figure 9.The Kaiser window with beta of 2.5 and its calibration result, (a) is Kaiser window and (b) is the calibration result.

Figure 9 .
Figure 9.The Kaiser window with beta of 2.5 and its calibration result, (a) is Kaiser window and (b) is the calibration result.

Figure 10 .
Figure 10.The PSF comparison of different cases in cross range and range, (a) is the cross range PSF and (b) is the range PSF.

Figure 10 .
Figure 10.The PSF comparison of different cases in cross range and range, (a) is the cross range PSF and (b) is the range PSF.

Figure 11 .Figure 12 .Figure 13 .
Figure 11.The array, targets and photo of measurements in experiments, (a) is array and targets; (b) is the experiment scene.

Figure 11 .
Figure 11.The array, targets and photo of measurements in experiments, (a) is array and targets; (b) is the experiment scene.

Figure 11 .Figure 12 .Figure 13 .
Figure 11.The array, targets and photo of measurements in experiments, (a) is array and targets; (b) is the experiment scene.

Figure 12 .
Figure 12. 2D imaging result of the ball in experiments, (a) is the imaging result before calibration and (b) is the imaging result after calibration.

Figure 11 .Figure 12 .Figure 13 .
Figure 11.The array, targets and photo of measurements in experiments, (a) is array and targets; (b) is the experiment scene.

Figure 13 .
Figure 13.The PSF comparisons of the ball in experiments before and after calibration, (a) is the cross range PSF and (b) is the range PSF.

Figure 14 .
Figure 14.Imaging results of the resolution board in experiments before and after calibration, (a) is the imaging result before calibration and (b) is the imaging result after calibration.

Figure 14 .
Figure 14.Imaging results of the resolution board in experiments before and after calibration, (a) is the imaging result before calibration and (b) is the imaging result after calibration.

Table 1 .
The comparison of two calibration method performance.

Table 1 .
The comparison of two calibration method performance.

Table 1 .
The comparison of two calibration method performance.