Abstract A numerical reconstruction method is proposed, which is applied to image identification of defects detected in elastic solid samples, in the case when a circular Ultrasonic scanning provides a measurement of the scattering pattern over full interval of the incident polar angle. The problem is first formulated as a system of respective boundary integral equations whose solution is used to calculate the far-field scattering diagram. Then the inverse reconstruction problem is reduced to a minimization of a certain strongly nonlinear functional. The proposed numerical algorithm is tested on some examples of volumetric flaw. It is also evaluated the influence of the error in the input data on precision of the reconstruction.