We apply the hybrid method for determining the shape of a bounded object in the elastostatic half plane from given Cauchy data on the boundary. The identifiability of the shape is investigated. For the integral representation of the function and the traction on the boundaries, the Green's function approach based on Kelvin's fundamental solution is used. The approximation of the integral operators with various singularities is made by trigonometrical and sinc quadratures. The presented numerical experiments exhibit the feasibility of the hybrid method for the system of differential equations and its stability in the case of noisy data.