# Ultrasound Imaging Algorithm: Half-Matrix Focusing Method Based on Reciprocity

Authors: Xufei Guo; Yan Han; Pengfei Nie
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

The ultrasonic phased array total focusing method (TFM) has the advantages of high imaging signal-to-noise ratio (SNR) and high defect resolution, but the problem of large amount data capturing and processing limits its practical industrial applications. To reduce the imaging calculation demand of the total focusing method, a half-matrix focusing method (HFM) is proposed based on the acoustic reciprocity. The method simplifies the calculation process of full-matrix data capturing (FMC) and total focus imaging. The experimental results show that the signal obtained by the linear array transceiver sensor is highly consistent, and the imaging resolution and signal-to-noise ratio of the half-matrix focusing method are slightly lower than those of full-matrix focusing method and higher than those of the B-scan imaging. However, compared with TFM, data acquisition and computational efficiency using the HFM have been improved significantly.


---

## Body

## 1. Introduction

Ultrasonic phased array testing adopts multichannel delay transceiver technology to emit each array element so as to achieve the angle deflection of the acoustic beam and focus on depth changes and electronic scanning, which improve the ability to detect defects [1–4]. The ultrasonic phased array imaging method based on postprocessing has been deeply studied by scholars. Holmes et al. [5] first proposed the concepts of full-matrix capture (FMC) and established the full-focus imaging method (TFM) based on FMC. Compared with the conventional ultrasonic phased array focusing algorithm, TFM can realize the synthetic focusing of any point in the detection region and the imaging quality is obviously better [6–9]. This method has therefore been extensively used in the fields of aviation, nuclear power, human organs, and so on [10–12].However, due to the large amount of full-matrix data involved, TFM imaging calculations are time-consuming and have poor real-time processing in real time, which limit their industrial applications. To improve imaging efficiency of TFM, Sutcliffe et al. [13] proposed a method of hardware acceleration of postprocessing using GPU. Bouaziz et al. [14, 15] proposed a real-time method that uses parallel computing with multiple field-programmable gate arrays. But, the cost will increase significantly by improving the hardware architecture. Another research goal is to improve the imaging algorithm while maintaining the image quality. Hunter et al. [16] introduced a generalization of the wave-number algorithm for Fourier-domain full-matrix imaging. Hu et al. [17] proposed a sparse-TFM imaging method based on sparse array optimization and new edge-directed interpolation. For the problem of limited acquisition speed, Moreau et al. [18] proposed an effective aperture method based on far-field approximation. In order to improve the imaging efficiency, Zhao et al. [21] used the upper triangle data to complete the ultrasound focus imaging.In the article, a half-matrix focusing method (HFM) based on transceiver sensor reciprocity is developed to improve inspection efficiency. This method can effectively reduce full-matrix data capturing and TFM imaging time while maintaining the image quality. The remainder of this paper is organized as follows: In Section2, the FMC theory and the imaging methods are introduced. Section 3 presents the experimental methods and procedures for defect detection. Moreover, the experimental results of data processed with conventional TFM and HFM are analyzed. Section 4 compares HFM imaging and B-scan imaging and discusses the subsequent research directions of the HFM imaging. Finally, the conclusions are given in Section 5.


## 2. Methods

### 2.1. FMC and TFM

Figure1 shows the principle of full-matrix ultrasonic data capture [19, 20]. A linear array consisting of N elements excites the ultrasonic wave in sequence. During each excitation, the ultrasonic signal is received by every array element. Finally, an N×N A-scan data matrix is obtained. Here,Aijt is the N×N A-scan signal emitted by element i and received by element j. The format of full-matrix data is shown as follows: (1)A=A11A12…A1j…A1NA21A22…A2j…A2N………………AijAi2…Aij…AiN………………AN1AN2…ANj…ANN.Figure 1
Schematic diagram of full-matrix data capture.Figure2 shows the principle of total focus algorithm [19, 20]. Image reconstruction using full-matrix data is divided into two steps. First, discretization of the inspected region into a grid is performed. Secondly, calculation of the amplitude intensity at a grid point Pxp,zp is calculated by(2)Ipxp,zp=∑i=1N∑j=1NAijtijxp,zp,where tijxp,zp is the time of flight starting from the transmitter element i, passing through the focusing point Pxp,zp, and arriving at the receiver element j. The time of flight is defined by(3)tijxp,zp=xi−xp2+zp2+xj−xp2+zp2c,where c is the propagation velocity of ultrasound in the medium and xi,0 and xj,0 are the Cartesian coordinates of the transmitter and the receiver, respectively.Figure 2
Schematic diagram of the total focus algorithm.


### 2.2. HFM

Consider that the performance indicators of each channel of the linear ultrasonic array transducer have excellent consistency. That is, the signal that is transmitted from elementi and received on element j is identical to the signal transmitted on element j and received on element i. It is named acoustic reciprocity. According to the principle of acoustic reciprocity, the full-matrix data shown in equation (1) have good symmetry. So the upper triangular or lower triangular matrix signal of the full-matrix data can be used in computational imaging. This method is called the half-matrix focus imaging (HFM), shown as follows (4) [21] or (5):(4)Ipxp,zp=∑i=1N∑j=iNAijtijxp,zp,(5)Ipxp,zp=∑i=jN∑j=1NAijtijxp,zp.During data acquisition and imaging calculation, using the half-matrix focusing method,N×N sets of A-scan data will be reduced to N×N+1/2. Before imaging calculation, it is necessary to normalize the image data. In the article, the amplitude data are normalized to between −128 and 127. It can superimpose and offset the random noise of positive and negative values [21]. The amplitude is normalized by(6)Ipxp,zp=Ipxp,zpImaxxp,zp×127,Ipxp,zp≥0,Ipxp,zp=−Ipxp,zpIminxp,zp×128,Ipxp,zp<0.


## 2.1. FMC and TFM

Figure1 shows the principle of full-matrix ultrasonic data capture [19, 20]. A linear array consisting of N elements excites the ultrasonic wave in sequence. During each excitation, the ultrasonic signal is received by every array element. Finally, an N×N A-scan data matrix is obtained. Here,Aijt is the N×N A-scan signal emitted by element i and received by element j. The format of full-matrix data is shown as follows: (1)A=A11A12…A1j…A1NA21A22…A2j…A2N………………AijAi2…Aij…AiN………………AN1AN2…ANj…ANN.Figure 1
Schematic diagram of full-matrix data capture.Figure2 shows the principle of total focus algorithm [19, 20]. Image reconstruction using full-matrix data is divided into two steps. First, discretization of the inspected region into a grid is performed. Secondly, calculation of the amplitude intensity at a grid point Pxp,zp is calculated by(2)Ipxp,zp=∑i=1N∑j=1NAijtijxp,zp,where tijxp,zp is the time of flight starting from the transmitter element i, passing through the focusing point Pxp,zp, and arriving at the receiver element j. The time of flight is defined by(3)tijxp,zp=xi−xp2+zp2+xj−xp2+zp2c,where c is the propagation velocity of ultrasound in the medium and xi,0 and xj,0 are the Cartesian coordinates of the transmitter and the receiver, respectively.Figure 2
Schematic diagram of the total focus algorithm.


## 2.2. HFM

Consider that the performance indicators of each channel of the linear ultrasonic array transducer have excellent consistency. That is, the signal that is transmitted from elementi and received on element j is identical to the signal transmitted on element j and received on element i. It is named acoustic reciprocity. According to the principle of acoustic reciprocity, the full-matrix data shown in equation (1) have good symmetry. So the upper triangular or lower triangular matrix signal of the full-matrix data can be used in computational imaging. This method is called the half-matrix focus imaging (HFM), shown as follows (4) [21] or (5):(4)Ipxp,zp=∑i=1N∑j=iNAijtijxp,zp,(5)Ipxp,zp=∑i=jN∑j=1NAijtijxp,zp.During data acquisition and imaging calculation, using the half-matrix focusing method,N×N sets of A-scan data will be reduced to N×N+1/2. Before imaging calculation, it is necessary to normalize the image data. In the article, the amplitude data are normalized to between −128 and 127. It can superimpose and offset the random noise of positive and negative values [21]. The amplitude is normalized by(6)Ipxp,zp=Ipxp,zpImaxxp,zp×127,Ipxp,zp≥0,Ipxp,zp=−Ipxp,zpIminxp,zp×128,Ipxp,zp<0.


## 3. Simulation and Results

### 3.1. Experimental Conditions

The finite element method (FEM) was used to validate our proposed imaging algorithm, and ABAQUS 2016 (Dassault,Velizy-Villacoublay, France) software was used to simulate the full-matrix detection of defects in a Q235 steel plate. The finite element model is shown in Figure3. The 32-element linear array transducer had a central frequency f = 5 MHz with an interelement spacing of 0.5 mm and an element width of 0.4 mm. Full-array probe parameters are listed in Table 1.Figure 3
Finite element inspection model (the enlarged part is array chip): (a) round through-hole; (b) oval through-hole; (c) rectangular through-hole.
(a)(b)(c)Table 1
Experimental array parameters.
ParameterValueNumber of elements32Element width0.4 mmElement pitch0.1 mmSampling time interval2e-8sCenter frequency5 MHzWave velocity5800 m/sThe defect positions in the steel plates are shown in Figure4. The TFM and HFM were tested on three finite element models of isotropic steel plates (100 mm × 30 mm). To suppress the interference of reflected waves, each model has infinite border with a thickness of 5 mm. In model 4(a), there are three round through-hole defects (Nos. 1, 2, and 3) with diameters of 2 mm. Between each defect, the lateral spacing was 5 mm. The distance between defect No. 1 and the upper edge of the steel plate model was 14 mm. In model 4(b), there is an oval through-hole (No. 4) with a long axis of 4 mm and a short axis of 2 mm. The distance between the defect and the upper edge of the steel plate model was 19 mm. In model 4(c), there is a rectangular through-hole (No. 5) defect with a length of 3 mm and width of 1 mm. The distance between the defect and the upper edge of the steel plate model was 20 mm.Figure 4
Schematic of defects in: (a) round through-hole; (b) oval through-hole; (c) rectangular through-hole.
(a)(b)(c)In this experiment, multiple analysis steps were set. Each element was excited by a 4-cycle sinusoidal signal modulated by a Hanning window in sequence at each analysis step. At the same time, every array element received ultrasonic signal. Finally, data for32×32 full matrix were captured. The signals were processed on MATLAB R2016a (MathWorks, Natick, MA, USA).


### 3.2. Results and Analysis

#### 3.2.1. Symmetry of Full-Matrix Data

Figure5 presents the amplitude data for 32×32 full matrix of round through-hole, oval through-hole, and rectangular through-hole, respectively. It can be seen that the data for full matrix were captured to be symmetrical. The result indirectly shows that only the upper or lower triangle data can be used for focus imaging, that is, half-matrix focus imaging. Therefore, the resolution of half-matrix imaging is very close to that of full-matrix focus imaging in theory.Figure 5
Full-matrix data: (a) round through-hole; (b) oval through-hole; (c) rectangular through-hole.
(a)(b)(c)


#### 3.2.2. Consistency of the Captured Signal

According to the reciprocity of the ultrasonic transceiver sensor, the acoustic signal captured by each transceiver channel remains consistent. The data of transceiver array elements exchange were extracted from full-matrix data. At the same time, the data were analyzed and compared. Here, six sets of A-scan data were randomly selected from full-matrix data. Signals received by different channels are shown in Figures6–8. Figure 6 presents A-scan signal of defect Nos. 1, 2, and 3 detection. Figure 7 presents A-scan signal of defect No. 4 detection. Figure 8 presents A-scan signal of defect No. 5 detection. It can be seen that the transceiver element is exchanged, and the collected signal basically coincides. This fully demonstrates the reciprocity of the ultrasonic transceiver channel. This further illustrates that half-matrix data can be used in efficient focus imaging.Figure 6
A-scan data of round through-hole: (a)A24−19 and A19−24; (b) A15−4 and A4−15.
(a)(b)Figure 7
A-scan data of oval through-hole: (a)A18−2 and A2−18; (b) A32−1 and A1−32
(a)(b)Figure 8
A-scan data of rectangular through-hole: (a)A22−21 and A21−22; (b) A16−13 and A13−16.
(a)(b)


#### 3.2.3. Comparison of Imaging Results

In the section, the difference between the imaging results of TFM and HFM was compared. According to the principle of TFM and HFM, we performed imaging experiments on round through-hole defects, oval through-hole defect, and rectangular through-hole defect, respectively. The imaging results are shown in Figures9–11.Figure 9
Focus imaging of round through-hole: (a) TFM; (b) HFM.
(a)(b)Figure 10
Focus imaging of oval through-hole: (a) TFM; (b) HFM.
(a)(b)Figure 11
Focus imaging of rectangular through-hole: (a) TFM; (b) HFM.
(a)(b)Comparing the two-dimensional imaging performance of the two algorithms, the five defects can be clearly displayed by TFM and HFM. However, it can be seen there is scattered noise around defects (mainly in the horizontal direction) by using the half-matrix imaging method. And, the shape and size of the defect have changed slightly. This is due to the small amount of data processed by the half-matrix imaging algorithm and poor signal average effect.


#### 3.2.4. Image Synthesis A-Scan Signal and SNR

In order to quantify the imaging performance of the two methods, image synthesis A-scan signal is used as shown in Figure12. It can be seen that image synthesis signals of Nos. 1–5 by TFM and HFM coincide. The amplitude and phase of the signal are basically equal.Figure 12
Image synthesis amplitude signal of defects: (a) round through-hole; (b) oval through-hole; (c) rectangular through-hole.
(a)(b)(c)For quantitatively analyzing the performance of the two imaging methods, the signal-to-noise ratio (SNR) is used to indicate the relationship between the defect signal and noise. The definition of SNR for a defect can be expressed [21–23] by(7)SNR=20×log10AmaxAave,where Amax is the maximum value (peak-to-peak) of the defect signal in the surrounding region and Aave is the average of noise signal. Except for the defect signal, the rest of the data is regarded as a noise signal. Here, we analyzed the SNR of five defects in the three steel plate models for TFM and HFM. The results are shown in Figure 13.Figure 13
Comparison of signal-to-noise ratio of TFM and HFM.The SNR for different defects by using TFM is better than the SNR by using HFM. But, the range of difference in SNR of TFM and HFM is very small. It can be seen that the difference range is between 1.09 dB and 1.87 dB. Since more data are superimposed during total focus imaging, its waveform average effect is more obvious and the SNR is slightly higher.In terms of imaging calculation time, the conditions are that the imaging pixel is set to 0.1 mm× 0.1 mm and the imaging point is set to 1024000. The CPU model of computer is Intel Core i5-9300H (2.4 GHz), and RAM is 8 GB. The full-matrix focus imaging takes 8.35 seconds, while the half-matrix focus imaging takes 4.1 seconds. Compared with TFM, the imaging time by using HFM is approximately doubled and imaging efficiency greatly improved.


## 3.1. Experimental Conditions

The finite element method (FEM) was used to validate our proposed imaging algorithm, and ABAQUS 2016 (Dassault,Velizy-Villacoublay, France) software was used to simulate the full-matrix detection of defects in a Q235 steel plate. The finite element model is shown in Figure3. The 32-element linear array transducer had a central frequency f = 5 MHz with an interelement spacing of 0.5 mm and an element width of 0.4 mm. Full-array probe parameters are listed in Table 1.Figure 3
Finite element inspection model (the enlarged part is array chip): (a) round through-hole; (b) oval through-hole; (c) rectangular through-hole.
(a)(b)(c)Table 1
Experimental array parameters.
ParameterValueNumber of elements32Element width0.4 mmElement pitch0.1 mmSampling time interval2e-8sCenter frequency5 MHzWave velocity5800 m/sThe defect positions in the steel plates are shown in Figure4. The TFM and HFM were tested on three finite element models of isotropic steel plates (100 mm × 30 mm). To suppress the interference of reflected waves, each model has infinite border with a thickness of 5 mm. In model 4(a), there are three round through-hole defects (Nos. 1, 2, and 3) with diameters of 2 mm. Between each defect, the lateral spacing was 5 mm. The distance between defect No. 1 and the upper edge of the steel plate model was 14 mm. In model 4(b), there is an oval through-hole (No. 4) with a long axis of 4 mm and a short axis of 2 mm. The distance between the defect and the upper edge of the steel plate model was 19 mm. In model 4(c), there is a rectangular through-hole (No. 5) defect with a length of 3 mm and width of 1 mm. The distance between the defect and the upper edge of the steel plate model was 20 mm.Figure 4
Schematic of defects in: (a) round through-hole; (b) oval through-hole; (c) rectangular through-hole.
(a)(b)(c)In this experiment, multiple analysis steps were set. Each element was excited by a 4-cycle sinusoidal signal modulated by a Hanning window in sequence at each analysis step. At the same time, every array element received ultrasonic signal. Finally, data for32×32 full matrix were captured. The signals were processed on MATLAB R2016a (MathWorks, Natick, MA, USA).


## 3.2. Results and Analysis

### 3.2.1. Symmetry of Full-Matrix Data

Figure5 presents the amplitude data for 32×32 full matrix of round through-hole, oval through-hole, and rectangular through-hole, respectively. It can be seen that the data for full matrix were captured to be symmetrical. The result indirectly shows that only the upper or lower triangle data can be used for focus imaging, that is, half-matrix focus imaging. Therefore, the resolution of half-matrix imaging is very close to that of full-matrix focus imaging in theory.Figure 5
Full-matrix data: (a) round through-hole; (b) oval through-hole; (c) rectangular through-hole.
(a)(b)(c)


### 3.2.2. Consistency of the Captured Signal

According to the reciprocity of the ultrasonic transceiver sensor, the acoustic signal captured by each transceiver channel remains consistent. The data of transceiver array elements exchange were extracted from full-matrix data. At the same time, the data were analyzed and compared. Here, six sets of A-scan data were randomly selected from full-matrix data. Signals received by different channels are shown in Figures6–8. Figure 6 presents A-scan signal of defect Nos. 1, 2, and 3 detection. Figure 7 presents A-scan signal of defect No. 4 detection. Figure 8 presents A-scan signal of defect No. 5 detection. It can be seen that the transceiver element is exchanged, and the collected signal basically coincides. This fully demonstrates the reciprocity of the ultrasonic transceiver channel. This further illustrates that half-matrix data can be used in efficient focus imaging.Figure 6
A-scan data of round through-hole: (a)A24−19 and A19−24; (b) A15−4 and A4−15.
(a)(b)Figure 7
A-scan data of oval through-hole: (a)A18−2 and A2−18; (b) A32−1 and A1−32
(a)(b)Figure 8
A-scan data of rectangular through-hole: (a)A22−21 and A21−22; (b) A16−13 and A13−16.
(a)(b)


### 3.2.3. Comparison of Imaging Results

In the section, the difference between the imaging results of TFM and HFM was compared. According to the principle of TFM and HFM, we performed imaging experiments on round through-hole defects, oval through-hole defect, and rectangular through-hole defect, respectively. The imaging results are shown in Figures9–11.Figure 9
Focus imaging of round through-hole: (a) TFM; (b) HFM.
(a)(b)Figure 10
Focus imaging of oval through-hole: (a) TFM; (b) HFM.
(a)(b)Figure 11
Focus imaging of rectangular through-hole: (a) TFM; (b) HFM.
(a)(b)Comparing the two-dimensional imaging performance of the two algorithms, the five defects can be clearly displayed by TFM and HFM. However, it can be seen there is scattered noise around defects (mainly in the horizontal direction) by using the half-matrix imaging method. And, the shape and size of the defect have changed slightly. This is due to the small amount of data processed by the half-matrix imaging algorithm and poor signal average effect.


### 3.2.4. Image Synthesis A-Scan Signal and SNR

In order to quantify the imaging performance of the two methods, image synthesis A-scan signal is used as shown in Figure12. It can be seen that image synthesis signals of Nos. 1–5 by TFM and HFM coincide. The amplitude and phase of the signal are basically equal.Figure 12
Image synthesis amplitude signal of defects: (a) round through-hole; (b) oval through-hole; (c) rectangular through-hole.
(a)(b)(c)For quantitatively analyzing the performance of the two imaging methods, the signal-to-noise ratio (SNR) is used to indicate the relationship between the defect signal and noise. The definition of SNR for a defect can be expressed [21–23] by(7)SNR=20×log10AmaxAave,where Amax is the maximum value (peak-to-peak) of the defect signal in the surrounding region and Aave is the average of noise signal. Except for the defect signal, the rest of the data is regarded as a noise signal. Here, we analyzed the SNR of five defects in the three steel plate models for TFM and HFM. The results are shown in Figure 13.Figure 13
Comparison of signal-to-noise ratio of TFM and HFM.The SNR for different defects by using TFM is better than the SNR by using HFM. But, the range of difference in SNR of TFM and HFM is very small. It can be seen that the difference range is between 1.09 dB and 1.87 dB. Since more data are superimposed during total focus imaging, its waveform average effect is more obvious and the SNR is slightly higher.In terms of imaging calculation time, the conditions are that the imaging pixel is set to 0.1 mm× 0.1 mm and the imaging point is set to 1024000. The CPU model of computer is Intel Core i5-9300H (2.4 GHz), and RAM is 8 GB. The full-matrix focus imaging takes 8.35 seconds, while the half-matrix focus imaging takes 4.1 seconds. Compared with TFM, the imaging time by using HFM is approximately doubled and imaging efficiency greatly improved.


## 3.2.1. Symmetry of Full-Matrix Data

Figure5 presents the amplitude data for 32×32 full matrix of round through-hole, oval through-hole, and rectangular through-hole, respectively. It can be seen that the data for full matrix were captured to be symmetrical. The result indirectly shows that only the upper or lower triangle data can be used for focus imaging, that is, half-matrix focus imaging. Therefore, the resolution of half-matrix imaging is very close to that of full-matrix focus imaging in theory.Figure 5
Full-matrix data: (a) round through-hole; (b) oval through-hole; (c) rectangular through-hole.
(a)(b)(c)


## 3.2.2. Consistency of the Captured Signal

According to the reciprocity of the ultrasonic transceiver sensor, the acoustic signal captured by each transceiver channel remains consistent. The data of transceiver array elements exchange were extracted from full-matrix data. At the same time, the data were analyzed and compared. Here, six sets of A-scan data were randomly selected from full-matrix data. Signals received by different channels are shown in Figures6–8. Figure 6 presents A-scan signal of defect Nos. 1, 2, and 3 detection. Figure 7 presents A-scan signal of defect No. 4 detection. Figure 8 presents A-scan signal of defect No. 5 detection. It can be seen that the transceiver element is exchanged, and the collected signal basically coincides. This fully demonstrates the reciprocity of the ultrasonic transceiver channel. This further illustrates that half-matrix data can be used in efficient focus imaging.Figure 6
A-scan data of round through-hole: (a)A24−19 and A19−24; (b) A15−4 and A4−15.
(a)(b)Figure 7
A-scan data of oval through-hole: (a)A18−2 and A2−18; (b) A32−1 and A1−32
(a)(b)Figure 8
A-scan data of rectangular through-hole: (a)A22−21 and A21−22; (b) A16−13 and A13−16.
(a)(b)


## 3.2.3. Comparison of Imaging Results

In the section, the difference between the imaging results of TFM and HFM was compared. According to the principle of TFM and HFM, we performed imaging experiments on round through-hole defects, oval through-hole defect, and rectangular through-hole defect, respectively. The imaging results are shown in Figures9–11.Figure 9
Focus imaging of round through-hole: (a) TFM; (b) HFM.
(a)(b)Figure 10
Focus imaging of oval through-hole: (a) TFM; (b) HFM.
(a)(b)Figure 11
Focus imaging of rectangular through-hole: (a) TFM; (b) HFM.
(a)(b)Comparing the two-dimensional imaging performance of the two algorithms, the five defects can be clearly displayed by TFM and HFM. However, it can be seen there is scattered noise around defects (mainly in the horizontal direction) by using the half-matrix imaging method. And, the shape and size of the defect have changed slightly. This is due to the small amount of data processed by the half-matrix imaging algorithm and poor signal average effect.


## 3.2.4. Image Synthesis A-Scan Signal and SNR

In order to quantify the imaging performance of the two methods, image synthesis A-scan signal is used as shown in Figure12. It can be seen that image synthesis signals of Nos. 1–5 by TFM and HFM coincide. The amplitude and phase of the signal are basically equal.Figure 12
Image synthesis amplitude signal of defects: (a) round through-hole; (b) oval through-hole; (c) rectangular through-hole.
(a)(b)(c)For quantitatively analyzing the performance of the two imaging methods, the signal-to-noise ratio (SNR) is used to indicate the relationship between the defect signal and noise. The definition of SNR for a defect can be expressed [21–23] by(7)SNR=20×log10AmaxAave,where Amax is the maximum value (peak-to-peak) of the defect signal in the surrounding region and Aave is the average of noise signal. Except for the defect signal, the rest of the data is regarded as a noise signal. Here, we analyzed the SNR of five defects in the three steel plate models for TFM and HFM. The results are shown in Figure 13.Figure 13
Comparison of signal-to-noise ratio of TFM and HFM.The SNR for different defects by using TFM is better than the SNR by using HFM. But, the range of difference in SNR of TFM and HFM is very small. It can be seen that the difference range is between 1.09 dB and 1.87 dB. Since more data are superimposed during total focus imaging, its waveform average effect is more obvious and the SNR is slightly higher.In terms of imaging calculation time, the conditions are that the imaging pixel is set to 0.1 mm× 0.1 mm and the imaging point is set to 1024000. The CPU model of computer is Intel Core i5-9300H (2.4 GHz), and RAM is 8 GB. The full-matrix focus imaging takes 8.35 seconds, while the half-matrix focus imaging takes 4.1 seconds. Compared with TFM, the imaging time by using HFM is approximately doubled and imaging efficiency greatly improved.


## 4. Discussion

The method used in the article is similar to that described in literature [21], but the experimental results are different. It can be seen from the experimental results that, compared with the TFM, the disadvantage of the HFM is its low signal-to-noise ratio. However, compared with traditional ultrasound B-scan imaging, HFM imaging has a higher signal-to-noise ratio. In the article, under the same experimental conditions, B-scan imaging is performed on three kinds of defects, and the image obtained is shown in Figure14. It can be seen that there is background noise in the image.Figure 14
B-scan imaging of 3 different defects: (a) round through-hole; (b) oval through-hole; (c) rectangular through-hole.
(a)(b)(c)In the literature [24, 25], a microseismic waveform extraction method using unsupervised machine learning technique is proposed to improve the resolution of microseismic imaging. With reference to the method of seismic imaging, machine learning and waveform feature extraction methods can also be used in ultrasound half-matrix imaging to increase the signal-to-noise ratio and improve imaging resolution. For imaging efficiency, the literature [13] mentioned the use of GPU parallel acceleration to reduce the time of ultrasound TFM imaging. Therefore, the ultrasound HFM imaging can be combined with GPU parallel computing to further reduce the imaging time. The two methods discussed above deserve further research.


## 5. Conclusions

In this paper, a half-matrix focusing method based on transceiver sensor reciprocity is developed to focus imaging. The imaging quality and computational efficiency are analyzed. In terms of image performance, the synthesized signals of the half-matrix focus imaging algorithm and the full-matrix focus imaging algorithm have good consistency. But, the imaging resolution and SNR of TFM are slightly higher than those of HFM. However, the HFM can simplify the data acquisition and imaging operation process of the full-matrix focusing algorithm. The imaging time was approximately doubled. The HFM provide an important reference for the industrial application.



---

*Source: 8888469-2021-02-02.xml*