A Novel Sub-Bottom Profiler and Signal Processor

In this paper, we introduce a novel sub-bottom profiler, making good use of the Mills cross configuration of multibeam sonar and synthetic aperture techniques of the synthetic aperture sonar system. The receiver array is mounted along the ship keel, while the transmitter array is mounted perpendicular to the receiver array. With the synthetic aperture technique, the along-track resolution can be greatly improved. The system often suffers from motion error, which severely degrades the imaging performance. To solve this problem, the imaging algorithm with motion compensation (MC) is proposed. With the presented method, the motion error is first estimated based on overlapped elements between successive pulses. Then, the echo data is processed by using the range migration algorithm based on the phase center approximation (PCA) method, which simultaneously performs the MC with the estimated motion error. In order to validate the proposed sub-bottom profiler and data processing method, some simulations and lake trial results are discussed. The processing results of the real data further indicate that the presented configuration has great potential to find buried objects in seabed sediments.


Introduction
Nowadays, the exploration of buried sediment in the sea bottom draws researchers' attention. Using the penetration capability of low frequency signals in water, the sub-bottom profiler [1] is often exploited to detect the structure of soft sediment. Based on matched filtering, the chirp sub-bottom profiler [2,3] uses a chirp signal to obtain high resolution images. The spatial resolution of multibeam sonar [4] is obtained based on the narrow beamwidth. It is often configured with the Mills cross configuration by using a long linear projector aligned with the ship track and a long linear receiver array aligned perpendicular to the receiver array. The transmitter array generates a beam that is narrow in the along-track dimension and wide in the across-track dimension. The receiver array forms multiple beams that are wide in the along-track dimension and narrow in the across-track dimension. By applying different time delays to various receiver signals, the receiver beams can be steered in different directions. The two-way beam patterns are, therefore, narrow in both directions. However, the narrow beamwidth is obtained at the cost of large receiver arrays. The physical length of the receiver array is usually constrained by the sonar platform. Additionally, the resolution decreases with the increase of range. Based on the nonlinear properties of sound propagation, the parametric sub-bottom profiler [5,6] can generate a low frequency via the difference between two transmitted signals. A new signal is used to detect the buried sediment in the sea bottom. However, the efficiency denotes the total number of receivers. The first receiver is close to the transmitter, and the last receiver is far away from the transmitter. The height of the sonar system is supposed to be r. At the beginning, the phase center coordinates of the transmitter array are (0, 0, r), while the coordinates of the i-th receiver are (0, −di, r). When the moving distance of the transmitter array is vt, the coordinates of the i-th receiver are (0, vt − di, r). Here, t denotes the slow time in the along-track dimension. The echo signal of the i-th receiver is the transmitted signal delayed by the two-way time corresponding to the two-way slant range. For the i-th receiver, the two-way slant range from the transmitter to the target and then back to the i-th receiver is given by where i  denotes the accurate propagation time of the transmitted signal [22,23]. Here, represents the moving distance in the along-track dimension during the echo collection;  is the fast time in the across-track dimension.
With the two-way slant range, the backscattered signal of the i-th receiver after demodulation is given by  (2) where () p  denotes the transmitted chirp signal and  represents the wavelength corresponding to the center frequency c f . In order to simplify notations, the beam patterns of the transmitter and receiver are not considered, and constants are also discarded. At the beginning, the phase center coordinates of the transmitter array are (0, 0, r), while the coordinates of the i-th receiver are (0, −d i , r). When the moving distance of the transmitter array is vt, the coordinates of the i-th receiver are (0, vt − d i , r). Here, t denotes the slow time in the along-track dimension. The echo signal of the i-th receiver is the transmitted signal delayed by the two-way time corresponding to the two-way slant range. For the i-th receiver, the two-way slant range from the transmitter to the target and then back to the i-th receiver is given by where τ i denotes the accurate propagation time of the transmitted signal [22,23]. Here, v · τ i represents the moving distance in the along-track dimension during the echo collection; τ is the fast time in the across-track dimension.
With the two-way slant range, the backscattered signal of the i-th receiver after demodulation is given by where p(τ) denotes the transmitted chirp signal and λ represents the wavelength corresponding to the center frequency f c . In order to simplify notations, the beam patterns of the transmitter and receiver are not considered, and constants are also discarded.
Inspecting the configuration shown in Figure 1, it can be found that the presented SAS configuration looks downward along the moving path. This is the major difference compared with traditional side-scan SAS systems [11]. Generally speaking, side-scan SAS systems focus on the wide swath, while the SAS sub-bottom profiler only concentrates on sediment penetration with low frequency signals. That is to say, the narrow swath is the main characteristic of the presented configuration. For this reason, the y-z plane in Figure 1 is considered to be the two-dimensional (2D) imaging geometry for the presented configuration. For simplicity, the medium refraction index is not considered in this paper.

Influence of Motion Error on Imagery
When the sonar platform is towed in water, the system may suffer from the effects of waves, tidal currents, tow boat wake, medium fluctuation, and so on. Due to these reasons, the sonar platform deviates from the nominal path. Under this condition, the SAS system suffers from motion errors, which can severely affect the SAS imagery. In order to easily describe the motion of the sonar platform, a set of three orthogonal axes is placed onto the sonar platform. Figure 2 shows the common terms describing the motion error. Inspecting the configuration shown in Figure 1, it can be found that the presented SAS configuration looks downward along the moving path. This is the major difference compared with traditional side-scan SAS systems [11]. Generally speaking, side-scan SAS systems focus on the wide swath, while the SAS sub-bottom profiler only concentrates on sediment penetration with low frequency signals. That is to say, the narrow swath is the main characteristic of the presented configuration. For this reason, the y-z plane in Figure 1 is considered to be the two-dimensional (2D) imaging geometry for the presented configuration. For simplicity, the medium refraction index is not considered in this paper.

Influence of Motion Error on Imagery
When the sonar platform is towed in water, the system may suffer from the effects of waves, tidal currents, tow boat wake, medium fluctuation, and so on. Due to these reasons, the sonar platform deviates from the nominal path. Under this condition, the SAS system suffers from motion errors, which can severely affect the SAS imagery. In order to easily describe the motion of the sonar platform, a set of three orthogonal axes is placed onto the sonar platform. Figure 2 shows the common terms describing the motion error. From Figure 2, the sonar system has six possible motion errors [19], which can be further decomposed into translational motion and rotational motion. The translational motion, including the surge, sway, and heave. is defined as variations of a straight, level, and constant moving heights. The surge is a forward or backward shift in the y axis. The sway is a side to side translation, and it is a horizontal shift in the x axis. The heave is a vertical translation in the z axis. The pitch, roll, and yaw belonging to the rotational motion are angular rotations around the x, y, and z axes, respectively. The motion errors can distort the imaging performance. Generally, the motion errors lengthening and reducing the expected time delay would affect the imagery. In this section, simulations are performed to discuss the influence of motion error on the imaging performance. The system parameters are shown in Table 1.  From Figure 2, the sonar system has six possible motion errors [19], which can be further decomposed into translational motion and rotational motion. The translational motion, including the surge, sway, and heave. is defined as variations of a straight, level, and constant moving heights. The surge is a forward or backward shift in the y axis. The sway is a side to side translation, and it is a horizontal shift in the x axis. The heave is a vertical translation in the z axis. The pitch, roll, and yaw belonging to the rotational motion are angular rotations around the x, y, and z axes, respectively. The motion errors can distort the imaging performance. Generally, the motion errors lengthening and reducing the expected time delay would affect the imagery. In this section, simulations are performed to discuss the influence of motion error on the imaging performance. The system parameters are shown in Table 1. Generally, the influence of motion error on processing results cannot be thoroughly presented if a single target is chosen. At this point, several targets are usually set in the imaging scenario. For clarity, the imaging scenario is supposed to consist of five ideal targets, which are depicted in Figure 3a. In order to clearly describe the focusing performance, the targets are marked as T1, T2, T3, T4, and T5. Since the time domain algorithms [23][24][25][26] are often viewed as the most precise methods, we can use the imaging results of time domain methods to evaluate other imaging methods. For comparison, the echo data in the absence of motion error is first processed by using the back projection (BP) method [26], and the imaging result is shown in Figure 3b. From Figure 3b, the five targets are well constructed when the echo data is not contaminated by the motion error. Generally, the influence of motion error on processing results cannot be thoroughly presented if a single target is chosen. At this point, several targets are usually set in the imaging scenario. For clarity, the imaging scenario is supposed to consist of five ideal targets, which are depicted in Figure 3a. In order to clearly describe the focusing performance, the targets are marked as T1, T2, T3, T4, and T5. Since the time domain algorithms [23][24][25][26] are often viewed as the most precise methods, we can use the imaging results of time domain methods to evaluate other imaging methods. For comparison, the echo data in the absence of motion error is first processed by using the back projection (BP) method [26], and the imaging result is shown in Figure 3b. From Figure 3b, the five targets are well constructed when the echo data is not contaminated by the motion error.

Influence of Sway and Yaw
We suppose that the sway deviating from the nominal path is x  . In Figure 2, the negative yaw is defined as a clockwise rotational movement around the z-axis. In Figure 4, the 2D geometry of the sway and yaw in the x-y plane is depicted. For simplicity, the transmitter phase center is defined as the coordinate origin. Here, p  denotes the yaw and the subscript p is the pulse index.  When the transmitter is located at vt in the y-axis, the transmitter coordinates are ( , , ) x vt r  .
When the echo signal is received by the i-th receiver, the receiver coordinates are Here, the moving distance during the signal reception is approximated by 2 i r vv c   . Therefore, the two-way range from the transmitter to ideal target and then back to the i-th receiver is given by

Influence of Sway and Yaw
We suppose that the sway deviating from the nominal path is ∆x. In Figure 2, the negative yaw is defined as a clockwise rotational movement around the z-axis. In Figure 4, the 2D geometry of the sway and yaw in the x-y plane is depicted. For simplicity, the transmitter phase center is defined as the coordinate origin. Here, α p denotes the yaw and the subscript p is the pulse index. Generally, the influence of motion error on processing results cannot be thoroughly presented if a single target is chosen. At this point, several targets are usually set in the imaging scenario. For clarity, the imaging scenario is supposed to consist of five ideal targets, which are depicted in Figure 3a. In order to clearly describe the focusing performance, the targets are marked as T1, T2, T3, T4, and T5. Since the time domain algorithms [23][24][25][26] are often viewed as the most precise methods, we can use the imaging results of time domain methods to evaluate other imaging methods. For comparison, the echo data in the absence of motion error is first processed by using the back projection (BP) method [26], and the imaging result is shown in Figure 3b. From Figure 3b, the five targets are well constructed when the echo data is not contaminated by the motion error.

Influence of Sway and Yaw
We suppose that the sway deviating from the nominal path is x  . In Figure 2, the negative yaw is defined as a clockwise rotational movement around the z-axis. In Figure 4, the 2D geometry of the sway and yaw in the x-y plane is depicted. For simplicity, the transmitter phase center is defined as the coordinate origin. Here, p  denotes the yaw and the subscript p is the pulse index.  When the transmitter is located at vt in the y-axis, the transmitter coordinates are ( , , ) x vt r  .
When the echo signal is received by the i-th receiver, the receiver coordinates are   When the transmitter is located at vt in the y-axis, the transmitter coordinates are (∆x, vt, r). When the echo signal is received by the i-th receiver, the receiver coordinates are (∆x + d i sin α p , vt − d i cos α p + v · 2r c , r). Here, the moving distance during the signal reception is approximated by vτ i ≈ v · 2r c . Therefore, the two-way range from the transmitter to ideal target and then back to the i-th receiver is given by Inspecting Equation (3), the error introduced by the sway and yaw is very slight, as the inequality ∆x 2 +(∆x+d i sin α p ) 2 2r 1 is always met. We conclude that both motion errors do not degrade the imaging performance. In this section, the influence of sway and yaw on the imaging performance is analyzed individually.

Influence of Sway
We first discuss the influence of sway on the imaging performance under the condition of α p = 0. Since the deviation is far less than the target range, the inequality |∆x| r is always met. Here, |∆x| denotes the sway magnitude. Due to this reason, we conclude that the sway nearly has no effect on the SAS imagery. To validate this conclusion, the simulations are carried out.
Here, two canonical examples of synthetic aperture motion errors are discussed. We first introduce the sawtooth perturbation into the echo data, and the motion error is shown in Figure 5. is always met. We conclude that both motion errors do not degrade the imaging performance. In this section, the influence of sway and yaw on the imaging performance is analyzed individually.

Influence of Sway
We first discuss the influence of sway on the imaging performance under the condition of 0 p  = . Since the deviation is far less than the target range, the inequality xr  is always met.
Here, x  denotes the sway magnitude. Due to this reason, we conclude that the sway nearly has no effect on the SAS imagery. To validate this conclusion, the simulations are carried out.
Here, two canonical examples of synthetic aperture motion errors are discussed. We first introduce the sawtooth perturbation into the echo data, and the motion error is shown in Figure 5.  Figure 6a shows the imaging results. To visually compare the imaging performance, Figure 6bf shows the along-track slices of the focused target. In Figure 6, the ideal case denotes the along-track slices without motion error by using the back projection (BP) algorithm [26]. Inspecting Figure 6, the imaging performance with the sawtooth perturbation is nearly identical to that without the motion error. That is to say, the sawtooth perturbation does not affect the imaging performance.
The next experiment concentrates on the case of sinusoidal perturbation, which is shown in Figure 7.  Figure 6a shows the imaging results. To visually compare the imaging performance, Figure 6b-f shows the along-track slices of the focused target. In Figure 6, the ideal case denotes the along-track slices without motion error by using the back projection (BP) algorithm [26]. Inspecting Figure 6, the imaging performance with the sawtooth perturbation is nearly identical to that without the motion error. That is to say, the sawtooth perturbation does not affect the imaging performance.
The next experiment concentrates on the case of sinusoidal perturbation, which is shown in Figure 7.       After processing the echo data, the results are shown in Figure 8. From Figure 8, the along-track slices with sinusoidal sway mostly agree with slices in the absence of motion error, and we nearly obtain the same imaging result with the ideal case. After processing the echo data, the results are shown in Figure 8. From Figure 8, the along-track slices with sinusoidal sway mostly agree with slices in the absence of motion error, and we nearly obtain the same imaging result with the ideal case. The peak sidelobe level ratio (PSLR), the integral sidelobe level ratio (ISLR) and deviation of along-track coordinates (DAC) are exploited to quantitatively evaluate the imaging performance.  The peak sidelobe level ratio (PSLR), the integral sidelobe level ratio (ISLR) and deviation of along-track coordinates (DAC) are exploited to quantitatively evaluate the imaging performance. Table 2 lists the quality parameters for the two aforementioned experiments. From Table 2, the PSLR and ISLR with sway error are mostly identical to the corresponding values without the motion error. The along-track coordinates of the ideal target are perfectly reconstructed, even if the echo data is contaminated by the sway. Based on the two aforementioned experiments, we draw a conclusion that the sway does not noticeably affect the imaging performance. Therefore, we can neglect the sway and focus our attention on other motion errors.

Influence of Yaw
As ∆x = 0 was set in Equation (3), the yaw influence on imaging performance is discussed now. It should be noted that the large rotational motion errors are directly obtained from INS. Due to this reason, we mainly concentrate on the microrotation motion error, which is usually beyond INS accuracy. Figure 9 shows the yaw perturbation, which is attributed to the Gaussian processes.   Table 2, the PSLR and ISLR with sway error are mostly identical to the corresponding values without the motion error. The along-track coordinates of the ideal target are perfectly reconstructed, even if the echo data is contaminated by the sway. Based on the two aforementioned experiments, we draw a conclusion that the sway does not noticeably affect the imaging performance. Therefore, we can neglect the sway and focus our attention on other motion errors.

Influence of Yaw
As 0 x = was set in Equation (3), the yaw influence on imaging performance is discussed now. It should be noted that the large rotational motion errors are directly obtained from INS. Due to this reason, we mainly concentrate on the microrotation motion error, which is usually beyond INS accuracy. Figure 9 shows the yaw perturbation, which is attributed to the Gaussian processes. The data contaminated by the yaw is directly processed, and the resultant images are shown in Figure 10. The data contaminated by the yaw is directly processed, and the resultant images are shown in Figure 10.   Table 2, the PSLR and ISLR with sway error are mostly identical to the corresponding values without the motion error. The along-track coordinates of the ideal target are perfectly reconstructed, even if the echo data is contaminated by the sway. Based on the two aforementioned experiments, we draw a conclusion that the sway does not noticeably affect the imaging performance. Therefore, we can neglect the sway and focus our attention on other motion errors.

Influence of Yaw
As 0 x = was set in Equation (3), the yaw influence on imaging performance is discussed now. It should be noted that the large rotational motion errors are directly obtained from INS. Due to this reason, we mainly concentrate on the microrotation motion error, which is usually beyond INS accuracy. Figure 9 shows the yaw perturbation, which is attributed to the Gaussian processes. The data contaminated by the yaw is directly processed, and the resultant images are shown in Figure 10. From Figure 10, we find that the yaw perturbation has little influence on the imaging performance. To some degree, it can be found from Equation (3) that the yaw influence is equivalent to that of sway on SAS imaging performance.
Based on simulations, we draw a conclusion that both sway and yaw essentially do not affect the imaging performance. Therefore, both motion errors can be neglected.

Influence of Heave and Pitch
In Figure 2, the positive pitch is defined as a clockwise rotational movement around the x-axis. In Figure 11, we depict the 2D geometry of heave and yaw in the y-z plane. Here, r  in Figure 11 is supposed to be the heave. For simplicity, the phase center of the transmitter is still defined as the coordinate origin, p  denotes the yaw, and subscript p is the pulse index. When the transmitter is located at vt in the y-axis, the transmitter coordinates are (0, , ) vt r r + . The coordinates of the i-th receiver are 2 (0, cos , sin ) when the echo signal is received by this receiver. Therefore, the two-way range from the transmitter to the ideal target and then back to the i-th receiver is given by From Figure 10, we find that the yaw perturbation has little influence on the imaging performance. To some degree, it can be found from Equation (3) that the yaw influence is equivalent to that of sway on SAS imaging performance.
Based on simulations, we draw a conclusion that both sway and yaw essentially do not affect the imaging performance. Therefore, both motion errors can be neglected.

Influence of Heave and Pitch
In Figure 2, the positive pitch is defined as a clockwise rotational movement around the x-axis. In Figure 11, we depict the 2D geometry of heave and yaw in the y-z plane. Here, ∆r in Figure 11 is supposed to be the heave. For simplicity, the phase center of the transmitter is still defined as the coordinate origin, β p denotes the yaw, and subscript p is the pulse index. From Figure 10, we find that the yaw perturbation has little influence on the imaging performance. To some degree, it can be found from Equation (3) that the yaw influence is equivalent to that of sway on SAS imaging performance.
Based on simulations, we draw a conclusion that both sway and yaw essentially do not affect the imaging performance. Therefore, both motion errors can be neglected.

Influence of Heave and Pitch
In Figure 2, the positive pitch is defined as a clockwise rotational movement around the x-axis. In Figure 11, we depict the 2D geometry of heave and yaw in the y-z plane. Here, r  in Figure 11 is supposed to be the heave. For simplicity, the phase center of the transmitter is still defined as the coordinate origin, p  denotes the yaw, and subscript p is the pulse index. When the transmitter is located at vt in the y-axis, the transmitter coordinates are (0, , ) vt r r + . The coordinates of the i-th receiver are 2 (0, cos , sin ) when the echo signal is received by this receiver. Therefore, the two-way range from the transmitter to the ideal target and then back to the i-th receiver is given by  When the transmitter is located at vt in the y-axis, the transmitter coordinates are (0, vt, r + ∆r). The coordinates of the i-th receiver are (0, vt − d i cos β p + v · 2r c , r + ∆r + d i sin β p ) when the echo signal is received by this receiver. Therefore, the two-way range from the transmitter to the ideal target and then back to the i-th receiver is given by From Equation (4), it can be seen that the heave and pitch lengthen or reduce the expected time delay. Therefore, the imaging performance must be sensitive to both motion errors. Here, the influence of pith and heave on the imaging performance is also discussed individually.

Influence of Heave
When the condition β p = 0 in Equation (4) is met, we can analyze the heave influence on imaging performance. Inspecting Equation (4), the last term denotes the motion error, which is introduced by heave. With the same magnitude of sway and heave (i.e., |∆r| = |∆x|), the inequality ∆x 2 r 2|∆r| is always met. Slight heave will severely influence the SAS imagery. Using the SAS parameters shown in Table 1 and the BP algorithm [26], we perform simulations to validate this conclusion.
We first focus on the processing of echo data when the heave is caused by sawtooth perturbation. The motion error is shown in Figure 5. The processing results and along-track slices are shown in Figure 12. Inspecting Figure 12, all targets are defocused, as the motion error has already lengthened or reduced the expected time delay. The motion error would lead to the deviation of along-track coordinates. Additionally, the sidelobe level is increasing. The focusing results of T1, T3, and T5 enhance this conclusion. Considering T2 and T4, the mainlobe is dramatically influenced by the sidelobe, and imaging results cannot provide useful information about targets. Therefore, the heave severely distorts the imaging performance. That is to say, the heave should be simultaneously compensated when targets are focused.
From Equation (4), it can be seen that the heave and pitch lengthen or reduce the expected time delay. Therefore, the imaging performance must be sensitive to both motion errors. Here, the influence of pith and heave on the imaging performance is also discussed individually.

Influence of Heave
When the condition 0 p  = in Equation (4) is met, we can analyze the heave influence on imaging performance. Inspecting Equation (4), the last term denotes the motion error, which is introduced by heave. With the same magnitude of sway and heave (i.e., rx  =  ), the inequality  is always met. Slight heave will severely influence the SAS imagery. Using the SAS parameters shown in Table 1 and the BP algorithm [26], we perform simulations to validate this conclusion.
We first focus on the processing of echo data when the heave is caused by sawtooth perturbation. The motion error is shown in Figure 5. The processing results and along-track slices are shown in Figure 12. Inspecting Figure 12, all targets are defocused, as the motion error has already lengthened or reduced the expected time delay. The motion error would lead to the deviation of along-track coordinates. Additionally, the sidelobe level is increasing. The focusing results of T1, T3, and T5 enhance this conclusion. Considering T2 and T4, the mainlobe is dramatically influenced by the sidelobe, and imaging results cannot provide useful information about targets. Therefore, the heave severely distorts the imaging performance. That is to say, the heave should be simultaneously compensated when targets are focused. Considering the heave caused by sinusoidal perturbation shown in Figure 7, we show the processing results in Figure 13. Inspecting Figure 13, we nearly obtain the same conclusions drawn from Figure 12, because the sinusoidal perturbation also lengthens or reduces the expected time delay. Table 3 lists the PSLR, ISLR, and DAC when the echo data is contaminated by the heave. Based on Table 3, it can be found that the quality parameters with the motion error are significantly lowered. Besides, the along-track position of the focused target dramatically deviates from that shown in Figure 3. Therefore, the microheave should be estimated and compensated when the SAS imagery is carried out. Considering the heave caused by sinusoidal perturbation shown in Figure 7, we show the processing results in Figure 13. Inspecting Figure 13, we nearly obtain the same conclusions drawn from Figure 12, because the sinusoidal perturbation also lengthens or reduces the expected time delay. Table 3 lists the PSLR, ISLR, and DAC when the echo data is contaminated by the heave. Based on Table 3, it can be found that the quality parameters with the motion error are significantly lowered. Besides, the along-track position of the focused target dramatically deviates from that shown in Figure 3. Therefore, the microheave should be estimated and compensated when the SAS imagery is carried out. Considering the heave caused by sinusoidal perturbation shown in Figure 7, we show the processing results in Figure 13. Inspecting Figure 13, we nearly obtain the same conclusions drawn from Figure 12, because the sinusoidal perturbation also lengthens or reduces the expected time delay. Table 3 lists the PSLR, ISLR, and DAC when the echo data is contaminated by the heave. Based on Table 3, it can be found that the quality parameters with the motion error are significantly lowered. Besides, the along-track position of the focused target dramatically deviates from that shown in Figure 3. Therefore, the microheave should be estimated and compensated when the SAS imagery is carried out.     (3), we now discuss the pitch influence on imaging performance. In Figure 9, the rotational motion error is also considered to be the pitch perturbation. After processing the data contaminated by the pitch error, the resulting images are depicted in Figure 14.  As ∆r = 0 was set in Equation (3), we now discuss the pitch influence on imaging performance. In Figure 9, the rotational motion error is also considered to be the pitch perturbation. After processing the data contaminated by the pitch error, the resulting images are depicted in Figure 14.   (3), we now discuss the pitch influence on imaging performance. In Figure 9, the rotational motion error is also considered to be the pitch perturbation. After processing the data contaminated by the pitch error, the resulting images are depicted in Figure 14. As seen in Figure 14, the pitch does affect the imaging performance. Fortunately, the influence of micropitch is tolerable. Figure 14a further enhances this conclusion. In practice, the large pitch is easily obtained from the INS. However, it is hard to obtain the micropitch, which is beyond the accuracy of INS. Inspecting Equation (4), the influence of pitch on SAS imaging performance is mostly equivalent to that of heave on the SAS focusing performance. At this point, we estimate the average heave for the SAS system in the case of 0 p  = . This assumption is reasonable, as the receiver array is very short.

Influence of Surge
The waves, tidal currents, tow boat wake, and medium fluctuation make the platform velocity unstable. For these reasons, the system would suffer from surge, which is a forward or backward shift in the along-track dimension. Under this condition, the sampled data in the along-track dimension is nonuniform [27,28]. Figure 15 shows the platform velocity. Based on sonar parameters shown in Table 1, the ideal velocity satisfying the uniform sampling in the along-track dimension should be 2 m/s. In practice, the platform velocity usually vibrates around the ideal velocity.  As seen in Figure 14, the pitch does affect the imaging performance. Fortunately, the influence of micropitch is tolerable. Figure 14a further enhances this conclusion. In practice, the large pitch is easily obtained from the INS. However, it is hard to obtain the micropitch, which is beyond the accuracy of INS. Inspecting Equation (4), the influence of pitch on SAS imaging performance is mostly equivalent to that of heave on the SAS focusing performance. At this point, we estimate the average heave for the SAS system in the case of β p = 0. This assumption is reasonable, as the receiver array is very short.

Influence of Surge
The waves, tidal currents, tow boat wake, and medium fluctuation make the platform velocity unstable. For these reasons, the system would suffer from surge, which is a forward or backward shift in the along-track dimension. Under this condition, the sampled data in the along-track dimension is nonuniform [27,28]. Figure 15 shows the platform velocity. Based on sonar parameters shown in Table 1, the ideal velocity satisfying the uniform sampling in the along-track dimension should be 2 m/s. In practice, the platform velocity usually vibrates around the ideal velocity. As seen in Figure 14, the pitch does affect the imaging performance. Fortunately, the influence of micropitch is tolerable. Figure 14a further enhances this conclusion. In practice, the large pitch is easily obtained from the INS. However, it is hard to obtain the micropitch, which is beyond the accuracy of INS. Inspecting Equation (4), the influence of pitch on SAS imaging performance is mostly equivalent to that of heave on the SAS focusing performance. At this point, we estimate the average heave for the SAS system in the case of 0 p  = . This assumption is reasonable, as the receiver array is very short.

Influence of Surge
The waves, tidal currents, tow boat wake, and medium fluctuation make the platform velocity unstable. For these reasons, the system would suffer from surge, which is a forward or backward shift in the along-track dimension. Under this condition, the sampled data in the along-track dimension is nonuniform [27,28]. Figure 15 shows the platform velocity. Based on sonar parameters shown in Table 1, the ideal velocity satisfying the uniform sampling in the along-track dimension should be 2 m/s. In practice, the platform velocity usually vibrates around the ideal velocity.  After processing the echo data contaminated by surge, the resultant images are shown in Figure 16. It should be noted that the average velocity (i.e., 2 m/s) is used by the imaging algorithm. From Figure 16, the surge has a slight influence on the imaging results. Generally, the imaging algorithm is directly performed when the platform velocity slightly deviates from the ideal one. If the platform velocity largely deviates from the ideal one, the nonuniform sampled data should be coerced into the uniform data before the imaging [27,28]. After processing the echo data contaminated by surge, the resultant images are shown in Figure 16. It should be noted that the average velocity (i.e., 2 m/s) is used by the imaging algorithm. From Figure 16, the surge has a slight influence on the imaging results. Generally, the imaging algorithm is directly performed when the platform velocity slightly deviates from the ideal one. If the platform velocity largely deviates from the ideal one, the nonuniform sampled data should be coerced into the uniform data before the imaging [27,28]. Inspecting Equations (3) and (4), the yaw and pitch also lead to the surge in the along-track dimension. However, the slight surge does not noticeably affect the imaging performance. At this point, the echo data suffering from the microsurge is still considered to be uniform data. Therefore, we can also neglect the influence of surge on SAS imagery.  Inspecting Equations (3) and (4), the yaw and pitch also lead to the surge in the along-track dimension. However, the slight surge does not noticeably affect the imaging performance. At this point, the echo data suffering from the microsurge is still considered to be uniform data. Therefore, we can also neglect the influence of surge on SAS imagery.

Motion Estimation and Imagery
The sonar system often suffers from motion error, which seriously affects the imaging performance. Therefore, we must estimate and compensate for the motion error. In this paper, the microheave beyond the INS accuracy is directly estimated from the sonar data. With the estimated motion error, we focus the sonar data using the range migration algorithm [29,30] based on the phase center approximation (PCA) method [16], because this algorithm is very efficient compared with the BP algorithm [26].

Estimation of Motion Error
The sonar platform is assumed to travel in a straight line with a constant speed. Due to unknown platform movement, deviations from this straight path often happen in real environments. We first start with the estimation of motion error.
Since the echo signal received by the receiver highly depends on the time delay of the transmitted signal, any change of time delay would degrade the focusing performance. With the SAS configuration shown in Figure 1, heave is the main motion error that directly affects the along-track resolution, which will be discussed in detail. The heave causes the platform to be either lifted up toward the ocean surface or pushed down toward the ocean floor. It has the effect of shortening or lengthening the overall time delay from the moment a pulse is transmitted to the echo from the received target. Since heave equally affects all of the receivers in the case of a multiple receiver system, the extra time delay is mostly identical for each receiver.
From Figure 17, the transmitter and each receiver can be modelled by the transducer located at the midpoint based on the DPCA. K = N − M equivalent elements overlap from pulse to pulse by operating at a speed of v = Md 2PRI . Here, PRI represents the pulse repetition interval. The receiver array includes N uniformly spaced receivers. The echo of M (M < N) receivers is used to focus the target. In fact, the SAS parameters are chosen in such a way that there is always overlapping of equivalent receiving arrays between two successive pulses in the opposite direction of a subset of equivalent receiving arrays, and there are some equivalent monostatic channels that receive exactly the same signal during two successive acquisitions in the absence of motion errors, regardless of the spatial distribution of elementary scatterers. For simplicity, we assume that the noise has been disregarded and that scatterer geometry and propagation medium have not changed.

Motion Estimation and Imagery
The sonar system often suffers from motion error, which seriously affects the imaging performance. Therefore, we must estimate and compensate for the motion error. In this paper, the microheave beyond the INS accuracy is directly estimated from the sonar data. With the estimated motion error, we focus the sonar data using the range migration algorithm [29,30] based on the phase center approximation (PCA) method [16], because this algorithm is very efficient compared with the BP algorithm [26].

Estimation of Motion Error
The sonar platform is assumed to travel in a straight line with a constant speed. Due to unknown platform movement, deviations from this straight path often happen in real environments. We first start with the estimation of motion error.
Since the echo signal received by the receiver highly depends on the time delay of the transmitted signal, any change of time delay would degrade the focusing performance. With the SAS configuration shown in Figure 1, heave is the main motion error that directly affects the along-track resolution, which will be discussed in detail. The heave causes the platform to be either lifted up toward the ocean surface or pushed down toward the ocean floor. It has the effect of shortening or lengthening the overall time delay from the moment a pulse is transmitted to the echo from the received target. Since heave equally affects all of the receivers in the case of a multiple receiver system, the extra time delay is mostly identical for each receiver.
From Figure 17, the transmitter and each receiver can be modelled by the transducer located at the midpoint based on the DPCA. K N M =− equivalent elements overlap from pulse to pulse by operating at a speed of 2 Here, PRI represents the pulse repetition interval. The receiver array includes N uniformly spaced receivers. The echo of M (M < N) receivers is used to focus the target. In fact, the SAS parameters are chosen in such a way that there is always overlapping of equivalent receiving arrays between two successive pulses in the opposite direction of a subset of equivalent receiving arrays, and there are some equivalent monostatic channels that receive exactly the same signal during two successive acquisitions in the absence of motion errors, regardless of the spatial distribution of elementary scatterers. For simplicity, we assume that the noise has been disregarded and that scatterer geometry and propagation medium have not changed. In the presence of motion errors, two equivalent receiving arrays between two successive pulses will no longer overlap but form a synthetic interferometric baseline. The change in round-trip travel time to each scatterer will be determined by the projection of the baseline in the radial direction to the scatterer. If the angular spread of this direction within a given sonar resolution cell or a set of resolution cells within a given time window remains small compared to the angular resolution of the interferometric baseline, then the changes in travel time will be the same for all targets in the resolution cell. For this time window, two received signals will, therefore, be nearly identical in shape, but will be received with slightly different delays after transmission. This difference provides an estimate of the line-of-sight motion error. In the presence of motion errors, two equivalent receiving arrays between two successive pulses will no longer overlap but form a synthetic interferometric baseline. The change in round-trip travel time to each scatterer will be determined by the projection of the baseline in the radial direction to the scatterer. If the angular spread of this direction within a given sonar resolution cell or a set of resolution cells within a given time window remains small compared to the angular resolution of the interferometric baseline, then the changes in travel time will be the same for all targets in the resolution cell. For this time window, two received signals will, therefore, be nearly identical in shape, but will be received with slightly different delays after transmission. This difference provides an estimate of the line-of-sight motion error.
Considering the influence of heave on the echo signal, the received echo of the k-th element is where k(k ∈ [1, K]) is the index of overlapped elements. Here, p(p ∈ [1, P]) is the pulse index, s k (τ, p) = ss(τ, Np − N + M + k) represents the signal magnitude, and ϕ p (τ, p) denotes the signal phase. We suppose that the sonar receivers suffer from the common motion error ∆r p for all targets. After transmitting the (p + 1)-th pulse, the received signal of the overlapped element is Since the computer can only obtain the wrapped phase distribution in (−π, π], the measured phase is written as where mod(·) denotes the modulus operator. After integrating the motion error, the unwrapped phase is obtained. We now concentrate on the estimation steps of motion error in detail. Inspecting the first term in Equation (4), the i-th receiver and transmitter can be approximated by a single transducer located midway between them, as the along-track coordinates of the i-th receiver, transmitter, and equivalent transducer are vt, vt − d i + v · 2r c , and vt − d i 2 + r v c , respectively. This approximation is equivalent to the PCA method [16]. With this approximation, the two-way path from the transmitter to the target and then back to the i-th receiver can be approximated by the two-way path between the target and equivalent transducer. However, this introduces range error, which is called the PCA error. The second term in Equation (4) is the PCA error. The last term in Equation (4) is introduced by the heave. Applying a 2D Fourier transformation to Equation (4), we obtain the point target reference spectrum (PTRS) based on the stationary phase of the method [31]. After neglecting the spectrum magnitude, the PTRS is given by with − 4π( f c + f τ )∆r (9) where P( f τ ) is the spectrum of the transmitted signal, f τ is the instantaneous frequency corresponding to the fast time in the across-track dimension, and f t is the Doppler frequency corresponding to the slow time in the along-track dimension. Inspecting Equation (9), the first term is similar to the PTRS phase of the monostatic SAS system. The second term is introduced by the stop-and-hop approximation. The target would suffer from the deviation of the along-track coordinates if it were not compensated. The third term denotes the spatial sampling in the along-track dimension. The fourth and last terms are introduced by the PCA error and motion error, respectively.
The heave is estimated by using the data of the overlapped elements between two successive pulses [20,32]. Inspecting Equation (4), the monostatic conversion from multireceiver data to monostatic equivalent data should first be performed. The PCA error that is the second term in Equation (4) should be compensated. Generally, the PCA error is decomposed into the Doppler phase error and microrange migration, which is caused by the spatial displacement of the transmitter and receiver. In order to discriminate this range migration from that introduced by the movement of the sonar system, we name it the bistatic range migration. We start with the compensation of the Doppler phase error introduced by the PCA error. Inspecting Equation (4), the filtering function used for the compensation of Doppler phase error is expressed as Using the spectrum of the transmitted signal, the pulse compression in the across-track dimension is performed. Examining the second term in Equation (4), the bistatic range migration is a range variant. Considering the wide swath SAS system, the interpolation should be used to deal with the space variance of the bistatic range migration. The presented system is mainly used for sub-bottom detection, and the mapping swath is very narrow. At this point, the compensation of bistatic range migration can also be carried out in the instantaneous frequency domain. The term is given by where * denotes the conjugate operation and r ref represents the reference range, which is often considered to be the height of the sonar platform.
After the monostatic conversion, the echo signal of all overlapped elements within the p-th pulse is superposed: The superposition of the echo signal gives us two advantages. One improves the robustness of motion error estimation, and the other improves the signal-to-noise ratio. Based on the same method, we superpose the echo signal of all overlapped elements within the (p + 1)-th pulse: Using s(τ, p) and s(τ, p + 1) between two successive pulses, we perform the correlation operation with lag τ = 0. This is given by where Σ denotes the summation in the across-track dimension. Based on Equation (14), the average phase difference for the p-th pulse is expressed as where arctan(·) denotes the operation of arctangent. Im(·) and Re(·) are the imaginary and real parts of the complex signal, respectively. Due to the wrapped characteristic, the phase error shown as Equation (15) is constrained within (−π, π]. The unwrapped phase error is further obtained by integrating the average phase error. This operation is given by With the center frequency f c , the time delay error introduced by the platform displacement is expressed as With the estimated time delay, the motion error is simultaneously compensated when the SAS imagery is carried out.

Imagery with the MC
Here, the MC is integrated into the SAS imaging algorithm. After the range compression and cancelation of bistatic range migration, we extract the M receiver data used for the SAS imagery. Then, we compensate the motion error in the across-track frequency domain. The filtering function is (18) It should be noted that the variance of time delay in the along-track dimension is ignored due to the narrow swath of the SAS sub-bottom profiler.
After this operation, the M × P receiver data collected with P pulses is arranged in order. The monostatic conversation from the multireceiver data to the monostatic equivalent data is also carried out. With the resultant data, the range migration algorithm [27,28] is used. Applying along-track FT to the data yields the spectrum in the 2D frequency domain. Based on Equation (19), the bulk focusing is performed.
With this step, the targets in the reference range are correctly focused. However, targets away from the reference range are only partially focused. Therefore, extra steps should be used to handle this issue. The Stolt mapping, which is also called the interpolation in the 2D frequency domain, can deal with this issue. This step remaps the instantaneous frequency axis based on Equation (20).
After this operation, targets away from the reference range are all focused. Then, the inverse Fourier transformation (IFT) in the across-track dimension is performed. After this processing, the along-track offset is corrected by using Equation (21): The last operation performs the IFT in the along-track dimension. After this step, the high resolution image is obtained.
Based on processing steps, the block diagram of the presented imaging method is shown in Figure 18.

Simulation and Real Data Processing
In this section, the simulated data and real data are used to verify the presented sub-bottom profiler and signal processor.

Simulation Results
This section concentrates on the processing of simulated data. The system parameters can be found in Table 1. The first experiment focuses on the sawtooth heave. With the presented method in Section 4.1, the motion errors shown in Figures 5 and 7 can be estimated. The original perturbation, the estimated error, and the difference between both curves are shown in Figure 19. From Figure 19, it can be seen that the presented method can estimate the motion error.

Simulation and Real Data Processing
In this section, the simulated data and real data are used to verify the presented sub-bottom profiler and signal processor.

Simulation Results
This section concentrates on the processing of simulated data. The system parameters can be found in Table 1. The first experiment focuses on the sawtooth heave. With the presented method in Section 4.1, the motion errors shown in Figures 5 and 7 can be estimated. The original perturbation, the estimated error, and the difference between both curves are shown in Figure 19. From Figure 19, it can be seen that the presented method can estimate the motion error. In order to compare the imaging performance, Figure 20 depicts the along-track slices of the focusing target. When we process the echo data not contaminated by the motion error based on the method in Section 4.2, we call it the ideal case. From Figure 20, the imaging performance is dramatically degraded if the echo data contaminated by the sawtooth perturbation are directly processed by the imaging algorithm. In Figure 20, slice sidelobes before MC are increasing. The imaging results shown in Figure 20b,c,f enhance this conclusion. Besides, the positions of focused targets deviate from the target positions shown in Figure 3a. The results shown in Figure 20c,d strengthen this conclusion. Furthermore, the ghost targets shown in Figure 20e are introduced by the motion error. Fortunately, the slices after MC mostly agree well with the slices of the ideal case. The presented method nearly gives the same performance as that in the absence of motion error. Therefore, we draw the conclusion that the imaging performance is greatly improved by using the presented method. In order to compare the imaging performance, Figure 20 depicts the along-track slices of the focusing target. When we process the echo data not contaminated by the motion error based on the method in Section 4.2, we call it the ideal case. From Figure 20, the imaging performance is dramatically degraded if the echo data contaminated by the sawtooth perturbation are directly processed by the imaging algorithm. In Figure 20, slice sidelobes before MC are increasing. The imaging results shown in Figure 20b,c,f enhance this conclusion. Besides, the positions of focused targets deviate from the target positions shown in Figure 3a. The results shown in Figure 20c,d strengthen this conclusion. Furthermore, the ghost targets shown in Figure 20e are introduced by the motion error. Fortunately, the slices after MC mostly agree well with the slices of the ideal case. The presented method nearly gives the same performance as that in the absence of motion error. Therefore, we draw the conclusion that the imaging performance is greatly improved by using the presented method. In order to compare the imaging performance, Figure 20 depicts the along-track slices of the focusing target. When we process the echo data not contaminated by the motion error based on the method in Section 4.2, we call it the ideal case. From Figure 20, the imaging performance is dramatically degraded if the echo data contaminated by the sawtooth perturbation are directly processed by the imaging algorithm. In Figure 20, slice sidelobes before MC are increasing. The imaging results shown in Figure 20b,c,f enhance this conclusion. Besides, the positions of focused targets deviate from the target positions shown in Figure 3a. The results shown in Figure 20c,d strengthen this conclusion. Furthermore, the ghost targets shown in Figure 20e are introduced by the motion error. Fortunately, the slices after MC mostly agree well with the slices of the ideal case. The presented method nearly gives the same performance as that in the absence of motion error. Therefore, we draw the conclusion that the imaging performance is greatly improved by using the presented method.  Table 4 lists the quality parameters. Based on quality parameters of T4 in Table 4, the PSLR and ISLR before MC differ from the values of the ideal case by about 15 and 24.7 dB, respectively. The DAC is 0.44 m. This conclusion is consistent with the focused results in Figure 20e, as the ghost target is introduced by the motion error. Inspecting quality parameters of other target in Table 4, the PSLR and ISLR before MC are also lowered compared with the corresponding values in the ideal case. The motion error leads to different DACs for various targets. That is to say, the distributed targets may suffer from distortion. Fortunately, the imaging performance is greatly improved by using the presented method. The PSLR and ISLR after MC mostly agree with the corresponding values in the ideal case. With the presented method, all target coordinates are reconstructed.  Figure 21 shows the processing results when the echo data is contaminated by the sinusoidal perturbation. In general, the motion error leads to four negative results. The first one is that the reconstructed targets before MC are no longer located at the ideal positions shown in Figure 3a. All focused targets before MC in Figure 21 enhance this conclusion. Since each scatterer owns various DACs, the distributed targets would be distorted. The second result is that the sidelobes are increasing. The slices in Figure 21b,d also strengthens this conclusion. Under this condition, it is very difficult to detect targets with weak reflectivity. The third result is that the ghost targets are introduced. Figure 21c demonstrates this conclusion. The last result is that the mainlobes of slices are broadening. In this circumstance, the resolution is lowered, and it is impossible for the system to distinguish between two neighboring scatterers. With the presented method, the aforementioned four issues can be solved. Inspecting Figure 21a, the presented method can reconstruct the targets, and the focusing performance is nearly identical to that without motion error.  Table 4 lists the quality parameters. Based on quality parameters of T4 in Table 4, the PSLR and ISLR before MC differ from the values of the ideal case by about 15 and 24.7 dB, respectively. The DAC is 0.44 m. This conclusion is consistent with the focused results in Figure 20e, as the ghost target is introduced by the motion error. Inspecting quality parameters of other target in Table 4, the PSLR and ISLR before MC are also lowered compared with the corresponding values in the ideal case. The motion error leads to different DACs for various targets. That is to say, the distributed targets may suffer from distortion. Fortunately, the imaging performance is greatly improved by using the presented method. The PSLR and ISLR after MC mostly agree with the corresponding values in the ideal case. With the presented method, all target coordinates are reconstructed.  Figure 21 shows the processing results when the echo data is contaminated by the sinusoidal perturbation. In general, the motion error leads to four negative results. The first one is that the reconstructed targets before MC are no longer located at the ideal positions shown in Figure 3a. All focused targets before MC in Figure 21 enhance this conclusion. Since each scatterer owns various DACs, the distributed targets would be distorted. The second result is that the sidelobes are increasing. The slices in Figure 21b,d also strengthens this conclusion. Under this condition, it is very difficult to detect targets with weak reflectivity. The third result is that the ghost targets are introduced. Figure 21c demonstrates this conclusion. The last result is that the mainlobes of slices are broadening. In this circumstance, the resolution is lowered, and it is impossible for the system to distinguish between two neighboring scatterers. With the presented method, the aforementioned four issues can be solved. Inspecting Figure 21a, the presented method can reconstruct the targets, and the focusing performance is nearly identical to that without motion error. With this experiment, we also present the quality parameters listed in Table 5.  Figure 21c. Inspecting the performance of the other target in Table 5  With this experiment, we also present the quality parameters listed in Table 5.  Figure 21c. Inspecting the performance of the other target in Table 5, the imaging performance before MC is seriously degraded compared with that of the ideal case. With the presented method, the motion error is well estimated and compensated, and the quality parameters with the presented method are mostly close to the corresponding values in the ideal case.
Based on simulations, we conclude that the presented method can estimate the motion error well. Additionally, the presented imaging method with the MC dramatically improves the focusing performance.

Real Data Processing
In this section, the real data is used to test the presented method. The test was performed in the Zhanghe reservoir. The transmitter array and receiver array have eight elements each. The diameter of each element is 0.08 m. The chirp signal is used by the system. The center frequency is 12 kHz, while the signal bandwidth is 7 kHz. The platform velocity is about 2 m/s. The pulse repetition interval is 0.06 s. Without the MC, the echo data is directly processed by the imaging processor presented in Section 4.2. The resulting image is shown in Figure 22. From Figure 22, we find that the targets circled by the red line are defocused. with that of the ideal case. With the presented method, the motion error is well estimated and compensated, and the quality parameters with the presented method are mostly close to the corresponding values in the ideal case. Based on simulations, we conclude that the presented method can estimate the motion error well. Additionally, the presented imaging method with the MC dramatically improves the focusing performance.

Real Data Processing
In this section, the real data is used to test the presented method. The test was performed in the Zhanghe reservoir. The transmitter array and receiver array have eight elements each. The diameter of each element is 0.08 m. The chirp signal is used by the system. The center frequency is 12 kHz, while the signal bandwidth is 7 kHz. The platform velocity is about 2 m/s. The pulse repetition interval is 0.06 s. Without the MC, the echo data is directly processed by the imaging processor presented in Section 4.2. The resulting image is shown in Figure 22. From Figure 22, we find that the targets circled by the red line are defocused. Inspecting Figure 22, the targets circled by the red circle are distorted by the uncompensated motion error. Based on the method in Section 4.1, we can estimate the motion error, which is shown in Figure 23.
With the estimated motion error and imaging method in Section 4.2, the real data is processed. The resultant image is shown in Figure 24. Inspecting Figure 24, the defocused targets in Figure 22 are well reconstructed by using the presented method. The imaging performance is greatly improved by using the presented method in this paper. Additionally, the sediment penetration with the low-frequency signal is successful. That is to say, the targets in the sediment can be detected by exploiting the presented SAS sub-bottom profiler.

Along-track(m)
Across-track Inspecting Figure 22, the targets circled by the red circle are distorted by the uncompensated motion error. Based on the method in Section 4.1, we can estimate the motion error, which is shown in Figure 23.
With the estimated motion error and imaging method in Section 4.2, the real data is processed. The resultant image is shown in Figure 24. Inspecting Figure 24, the defocused targets in Figure 22 are well reconstructed by using the presented method. The imaging performance is greatly improved by using the presented method in this paper. Additionally, the sediment penetration with the low-frequency signal is successful. That is to say, the targets in the sediment can be detected by exploiting the presented SAS sub-bottom profiler.
With the estimated motion error and imaging method in Section 4.2, the real data is processed. The resultant image is shown in Figure 24. Inspecting Figure 24, the defocused targets in Figure 22 are well reconstructed by using the presented method. The imaging performance is greatly improved by using the presented method in this paper. Additionally, the sediment penetration with the low-frequency signal is successful. That is to say, the targets in the sediment can be detected by exploiting the presented SAS sub-bottom profiler.

Conclusions
This paper first proposes a novel sub-bottom profiler based on the Mills cross configuration of multibeam sonar and the synthetic aperture technique of the SAS system, and then presents a signal processor improving the imaging performance of the presented sub-bottom profiler. On the one hand, the presented sub-bottom profiler exploits the Mills cross configuration of multibeam sonar. The receiver array is mounted along the ship's moving track, while the transmitter array is mounted perpendicular to the receiver array. This is the major difference between the presented system and multibeam sonar. On the other hand, the synthetic aperture technique of the SAS system is used by the presented system. The sub-bottom profiler presented in this paper looks downward along the moving path, while the traditional SAS systems use the side-scan configuration. Compared with traditional side-scan SAS systems, the transmitter array generates a beam that is narrow in the across-track dimension and wide in the along-track dimension. With the downward-looking configuration, the seafloor strongly reflects the energy in the incoming direction at the normal incidence rate. Furthermore, the synthetic aperture technique allows us to circumvent the trade-off between the sonar aperture size and along-track resolution, as a virtual array that is much longer than its physical length can be synthesized based on the motion of the sonar platform and coherently signal processing. Compared with traditional sub-bottom profilers, the presented system can generate high resolution images in the along-track dimension based on the synthetic aperture technique. Additionally, the signal energy can be efficiently used for sediment penetration with the presented system. At this point, the presented system is very suitable for sub-bottom detection.
The presented sub-bottom system often suffers from motion errors, which play an important role in the imaging performance. In general, motion errors lengthening and reducing the expected time delay would affect the imagery. This paper discusses the influence of motion error on the imaging performance, and indicates that except for heave, the motion errors have nearly no effect on the processing results. Considering the serious influence of heave, an error estimation method based on overlapped elements between successive pulses is presented. Then, a range migration algorithm integrated with MC is discussed. To validate the presented sub-bottom profiler and signal processor, the simulated data and real data based on the presented sub-bottom profiler are

Conclusions
This paper first proposes a novel sub-bottom profiler based on the Mills cross configuration of multibeam sonar and the synthetic aperture technique of the SAS system, and then presents a signal processor improving the imaging performance of the presented sub-bottom profiler. On the one hand, the presented sub-bottom profiler exploits the Mills cross configuration of multibeam sonar. The receiver array is mounted along the ship's moving track, while the transmitter array is mounted perpendicular to the receiver array. This is the major difference between the presented system and multibeam sonar. On the other hand, the synthetic aperture technique of the SAS system is used by the presented system. The sub-bottom profiler presented in this paper looks downward along the moving path, while the traditional SAS systems use the side-scan configuration. Compared with traditional side-scan SAS systems, the transmitter array generates a beam that is narrow in the across-track dimension and wide in the along-track dimension. With the downward-looking configuration, the seafloor strongly reflects the energy in the incoming direction at the normal incidence rate. Furthermore, the synthetic aperture technique allows us to circumvent the trade-off between the sonar aperture size and along-track resolution, as a virtual array that is much longer than its physical length can be synthesized based on the motion of the sonar platform and coherently signal processing. Compared with traditional sub-bottom profilers, the presented system can generate high resolution images in the along-track dimension based on the synthetic aperture technique. Additionally, the signal energy can be efficiently used for sediment penetration with the presented system. At this point, the presented system is very suitable for sub-bottom detection.
The presented sub-bottom system often suffers from motion errors, which play an important role in the imaging performance. In general, motion errors lengthening and reducing the expected time delay would affect the imagery. This paper discusses the influence of motion error on the imaging performance, and indicates that except for heave, the motion errors have nearly no effect on the processing results. Considering the serious influence of heave, an error estimation method based on overlapped elements between successive pulses is presented. Then, a range migration algorithm integrated with MC is discussed. To validate the presented sub-bottom profiler and signal processor, the simulated data and real data based on the presented sub-bottom profiler are processed. The processing results of the simulated data show that the heave can be efficiently estimated by using the presented motion estimation method. Besides, the processing results are greatly improved based on the presented signal processor. The processing results of the real data show that the presented system does penetrate the sub-bottom. The results further indicate that the presented sub-bottom profiler has great potential in sub-bottom detection. Moreover, the results after MC are improved compared with the results before MC. This further validates the signal processor discussed in this paper.
Our next task is to consider the medium refraction index when it comes to the imagery. In addition, other robust motion error estimation methods should be taken into account to improve the performance of the signal processor.
Author Contributions: C.T. analyzed the results. X.Z. checked the paper. P.Y. and M.S. conducted experiments and wrote the paper.