Whether hydroxyapatite (HA) orientation in fossilised bone samples can be non-destructively retrieved and used to determine the arrangement of the bone matrix and the location of muscle attachments (entheses), is a question of high relevance to palaeontology, as it facilitates a detailed understanding of the (micro-)anatomy of extinct species with no damage to the precious fossil specimens. Here, we report studies of two fossil bone samples, specifically the tibia of a 300-million-year-old tetrapod, Discosauriscus austriacus, and the humerus of a 370-million-year-old lobe-finned fish, Eusthenopteron foordi, using XRD-CT - a combination of X-ray diffraction (XRD) and computed tomography (CT). Reconstructed 3D images showing the spatial mineral distributions and the local orientation of HA were obtained. For Discosauriscus austriacus, details of the muscle attachments could be discerned. For Eusthenopteron foordi, the gross details of the preferred orientation of HA were deduced using three tomographic datasets obtained with orthogonally oriented rotation axes. For both samples, the HA in the bone matrix exhibited preferred orientation, with the unit cell c-axis of the HA crystallites tending to be parallel with the bone surface. In summary, we have demonstrated that XRD-CT combined with an intuitive reconstruction procedure is becoming a powerful tool for studying palaeontological samples.