This contribution deals with the application of the inverse homogenization method to the determination of geometrical properties of cancellous bone. The approach represents a combination of an extended version of the Marquardt-Levenberg method with the multiscale finite element method. The former belongs to the group of gradient-based optimization strategies, while the latter is a numerical homogenization method, suitable for the modeling of materials with a highly heterogeneous microstructure. The extension of the Marquardt-Levenberg method is concerned with the selection strategy for distinguishing the global minimum from the plethora of local minima. Within the numerical examples, the bone is modeled as a biphasic viscoelastic medium and three different representative volume elements are taken into consideration. Different models enable the simulation of the bone either as a purely isotropic or as a transversally anisotropic medium. Main geometrical properties of trabeculae are determined from data on effective shear modulus but alternative schemes are also possible.