We propose a finite-difference ghost-point approach for the numerical solution of Cauchy-Navier equations in linear elasticity problems on arbitrary unbounded domains. The technique is based on a smooth coordinate transformation, which maps an unbounded domain into a unit square. Arbitrary geometries are defined by suitable level-set functions. The equations are discretized by classical nine-point stencil on interior points, while boundary conditions and high order reconstructions are used to define the field variables at ghost-points, which are grid nodes external to the domain with a neighbor inside the domain. The linear system arising from such discretization is solved by a multigrid strategy. The approach is then applied to solve elasticity problems in volcanology for computing the displacement caused by pressure sources. The method is suitable to treat problems in which the geometry of the source often changes (explore the effects of different scenarios, or solve inverse problems in which the geometry itself is part of the unknown), since it does not require complex re-meshing when the geometry is modified. Several numerical tests are successfully performed, which asses the effectiveness of the present approach.