A method for estimating the skeletal status or
bone quality of a
vertebrate on the basis of two-dimensional image data comprising information relating to the trabecular structure of at least a part of a bone of the
vertebrate, the image data being data obtained by exposing at least the part of the bone to
electromagnetic radiation, such as X-rays, the method comprising subjecting the image data to a
statistical analysis comprising a
background correction procedure in which
low frequency intensity variations not related to the trabecular structure of the bone are reduced relative to image data related to the trabecular structure of the part of the bone, a
feature extraction procedure comprising (a) determining values reflecting the projected trabecular density in the image data, caused by the X-
ray attenuating properties of
cancellous bone in the part of the bone, for each of a number of locations or areas in the image data, (b) deriving one or more features from the variation of the determined PTD-values, preferably in the longitudinal direction of the bone, and an
estimation procedure in which the skeletal status of the
vertebrate is estimated on the basis of the one or more derived features and optionally other features related to the bone of the vertebrate and a predetermined relationship between the features and reference skeletal status parameters. Preferably, a profile describing the projected trabecullar density as a function of the distance along a line substantially at the center of the bone is determined, and information relating to skeletal status is derived from variations, fluctuations or other features of the profile.