Biomechanical properties of soft tissues can provide information regarding the local health status. Often the cells in pathological tissues can be found to form a stiff extracellular environment, which is a sensitive, early diagnostic indicator of disease. Quasi-static ultrasonic elasticity imaging provides a way to image the mechanical properties of tissues. Strain images provide a map of the relative tissue stiffness, but ambiguities and artifacts limit its diagnostic value. Accurately mapping intrinsic mechanical parameters of a region may increase diagnostic specificity. However, the inverse problem, whereby force and displacement estimates are used to estimate a constitutive matrix, is ill conditioned. Our method avoids many of the issues involved with solving the inverse problem, such as unknown boundary conditions and incomplete information about the stress field, by building an empirical model directly from measured data. Surface force and volumetric displacement data gathered during imaging are used in conjunction with the AutoProgressive method to teach artificial neural networks the stress-strain relationship of tissues. The Autoprogressive algorithm has been successfully used in many civil engineering applications and to estimate ocular pressure and corneal stiffness; here, we are expanding its use to any tissues imaged ultrasonically. We show that force-displacement data recorded with an ultrasound probe and displacements estimated at a few points in the imaged region can be used to estimate the full stress and strain vectors throughout an entire model while only assuming conservation laws. We will also demonstrate methods to parameterize the mechanical properties based on the stress-strain response of trained neural networks. This method is a fundamentally new approach to medical elasticity imaging that for the first time provides full stress and strain vectors from one set of observation data.