The imaging function \(I_{DSM}(z)\)  is expected to peak when \(z\in \partial D\)  and decays like the Bessel functions for sampling points away from the boundary. We refer to [1]} for the corresponding theoretical analysis and its connections with the other sampling methods.
