In this work we focus on the development of a numerical algorithm for the inverse elastography problem. The goal is to perform an efficient material parameter identification knowing the elastic displacement field induced by a mechanical load. We propose to define the inverse problem through a quadratic optimization program which uses the direct problem formulation to define the objective function. In this way, we end up with a convex minimization problem which attains its minimum at the solution of a linear system. The effectiveness of our method is illustrated through numeral examples.