A computer-implemented method of processing a medical image, the method comprising: receiving medical image data, the medical image data representing a medical image of at least a portion of a vertebral column; processing the medical image data to determine a plurality of positions within the image, each of the plurality of positions comprising a centre coordinate within the medical image, each centre coordinate representing a central point of a respective vertebral bone within the vertebral column, wherein the plurality of positions comprises positions representing a plurality of regions of interest, each region of interest bounding a respective vertebral bone; and processing data representing the plurality of positions, to determine a degree of deformity by quantifying and comparing vertebral compression fractures of at least one vertebral bone within the vertebral column, wherein the method comprises: determining, for each of the plurality of positions, a corresponding anatomical feature; assigning, to each of the plurality of regions of interest, a label representing the respective corresponding anatomical feature; generating a mask comprising the plurality of regions of interest and representing the imaged vertebral column on the basis of the corresponding labels and centre coordinates; segmenting each of the plurality of regions of interest of the image; processing each of the plurality of regions of interest to determine one or more sagittal points for each of the corresponding vertebral bones; and determining a Mahalanobis distance between each of the one or more sagittal points and the corresponding centre coordinates. receiving medical image data, the medical image data representing a medical image of at least a portion of a vertebral column; processing the medical image data to determine a plurality of positions within the image, each of the plurality of positions comprising a centre coordinate within the medical image, each centre coordinate representing a central point of a respective vertebral bone within the vertebral column, wherein the plurality of positions comprises positions representing a plurality of regions of interest, each region of interest bounding a respective vertebral bone; and processing data representing the plurality of positions, to determine a degree of deformity by quantifying and comparing vertebral compression fractures of at least one vertebral bone within the vertebral column, wherein the method comprises: determining, for each of the plurality of positions, a corresponding anatomical feature; assigning, to each of the plurality of regions of interest, a label representing the respective corresponding anatomical feature; generating a mask comprising the plurality of regions of interest and representing the imaged vertebral column on the basis of the corresponding labels and centre coordinates; segmenting each of the plurality of regions of interest of the image; processing each of the plurality of regions of interest to determine one or more sagittal points for each of the corresponding vertebral bones; and determining a Mahalanobis distance between each of the one or more sagittal points and the corresponding centre coordinates. determining, for each of the plurality of positions, a corresponding anatomical feature; assigning, to each of the plurality of regions of interest, a label representing the respective corresponding anatomical feature; generating a mask comprising the plurality of regions of interest and representing the imaged vertebral column on the basis of the corresponding labels and centre coordinates; segmenting each of the plurality of regions of interest of the image; processing each of the plurality of regions of interest to determine one or more sagittal points for each of the corresponding vertebral bones; and determining a Mahalanobis distance between each of the one or more sagittal points and the corresponding centre coordinates.