Abstract:We propose a reliable direct imaging method based on the reverse time migration for finding extended obstacles with phaseless total field data. We prove that the imaging resolution of the method is essentially the same as the imaging results using the scattering data with full phase information. The imaginary part of the cross-correlation imaging functional always peaks on the boundary of the obstacle. Numerical experiments are included to illustrate the powerful imaging quality.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: to reliably reconstruct the position and shape of an unknown obstacle in two - dimensional space by using the total - field data with missing phase (i.e., only amplitude information). Specifically, the author proposes a direct imaging method based on Reverse Time Migration (RTM) to find extended obstacles from full - field data containing only intensity measurements. The key of this method lies in proving that the imaging resolution using phase - missing data is basically the same as that obtained using data with complete phase information, and the imaging function always reaches its peak at the obstacle boundary.
### Problem Background
In diffraction - optical imaging and radar imaging systems, it is easier to measure the intensity of the full - field than to measure its phase information. Therefore, it is very necessary to develop reliable numerical methods for reconstructing obstacles using only phase - missing data (i.e., amplitude information). In recent years, there have been many studies in the literature on using phase - missing data to solve inverse scattering problems, but most methods rely on iteration or require prior information.
### Paper Contributions
This paper proposes a non - iterative direct imaging algorithm based on the Reverse Time Migration method, which can reconstruct obstacles using only intensity measurement data \(|u|\) under point - source excitation. The author proves that the imaging resolution of this algorithm is comparable to that obtained using data with complete phase information, and the imaging function always reaches its peak at the obstacle boundary. Moreover, this method does not require any prior information about the physical properties of the obstacle, such as whether it is penetrable and the type of boundary conditions for impenetrable obstacles.
### Mathematical Model
Consider a soft acoustic obstacle occupying a bounded Lipschitz domain \(D\subset\mathbb{R}^2\), where \(\nu\) is the unit outer normal vector of its boundary. The incident wave is \(u_i\), and the total field \(u = u_i+u_s\), where \(u_s\) is the solution of the following acoustic scattering problem:
\[
\begin{cases}
\Delta u_s + k^2u_s = 0 & \text{in }\mathbb{R}^2\setminus\overline{D},\\
u_s=-u_i & \text{on }\Gamma_D,\\
\sqrt{r}\left(\frac{\partial u_s}{\partial r}-ik u_s\right)\to 0 & \text{as }r = |x|\to+\infty.
\end{cases}
\]
where \(k > 0\) is the wave number, and the last condition is the outgoing Sommerfeld radiation condition, which ensures the uniqueness of the solution.
### Imaging Algorithm
The RTM algorithm includes two main steps:
1. **Back - propagation**: Calculate the modified data \(\Delta(x_r,x_s)=|u(x_r,x_s)|^2-|u_i(x_r,x_s)|^2/u_i(x_r,x_s)\), and then back - propagate it into the domain using the fundamental solution \(\Phi(x_r,z)\).
2. **Cross - correlation**: Calculate the imaginary part cross - correlation between the incident wave field \(u_i(z,x_s)\) and the back - propagated field to obtain the final imaging profile.
Through the above method, the author proves that the imaging function always reaches its peak at the obstacle boundary, and the imaging resolution is comparable to that obtained using data with complete phase information.
### Summary
The main contribution of this paper is to propose a direct imaging method based on Reverse Time Migration, which can effectively reconstruct the position and shape of an obstacle using only total - field data with missing phase. This method not only avoids prior assumptions about the physical properties of the obstacle, but also has a high imaging resolution.