Magnetic Resonance Imaging Versus Computed Tomography for Three‐Dimensional Bone Imaging of Musculoskeletal Pathologies: A Review

1Magnetic resonance imaging (MRI) is a radiation‐free, noninvasive imaging modality that provides three‐dimensional (3D) visualization of tissues. Its superior soft tissue contrast has made it a preferential diagnostic tool for the imaging of various organ systems, including the musculoskeletal system. Osseous structures are, however, usually visualized using radiography or computed tomography (CT). For imaging complex structures, CT is preferred as it offers high‐resolution 3D images with a radiodensity contrast that highlights bony tissues. Building upon the characteristic high X‐ray attenuation of cortical bone, dedicated (semi‐)automatic bone segmentation tools have been developed for CT images. Resulting 3D bone renderings have proven valuable in the diagnosis and treatment of bone pathologies. Consequently, pathologies affecting both soft and hard tissues, including skull,, spine,,, and joint disorders,, often warrant the acquisition of both MR and CT images. Such a multimodal workflow is logistically complex and induces an adverse radiation burden inherent to CT imaging, especially harmful in young population.

2Recent advances in MR image acquisition and processing, facilitated by the development of new hardware and the increase in computing power, have enabled the improvement of bone contrast on MR images. If reliable, MRI could be a radiation‐free alternative to CT for the diagnosis and treatment planning of certain musculoskeletal pathologies. Transforming a CT‐MR multimodal workflow into a simplified radiation‐free MR‐only workflow, as previously proposed in radiotherapy treatment planning, could lead to less hospital visits, lower costs, allow for the fusion of soft tissue and bone information, and reduce the time under sedation for younger patients., This review will discuss comparative studies of MRI and CT for the diagnosis and treatment planning of bone pathologies in musculoskeletal diseases in multiple anatomical regions, including the skull, the spine, the shoulder, and the pelvis. Four main subjects will be described: MRI‐based techniques for bone imaging, MRI for bone segmentation and 3D reconstruction, MRI for the diagnosis of bone pathologies, and the remaining challenges faced by MRI in the context of bone visualization. Applications in the fields of radiotherapy and positron emission tomography–magnetic resonance (PET‐MR) will not be covered as they have been thoroughly reviewed in the past few years.,,,

3Cortical bone imaging is challenging with MRI due to its low free‐water content. The MR signal that originates from cortical bone is mostly emitted by bound water, causing the signal decay to be rapid. Consequently, in conventional MRI sequences, cortical bone appears as a structure with low signal intensity that is not specific to bone. Although valuable for structural imaging, the poor visualization of cortical bone on conventional sequences has motivated the development of dedicated imaging techniques that facilitate bone visualization and segmentation. The remainder of this section briefly discusses several MR sequences and processing techniques used for bone imaging. The resulting MR images are compiled in Fig. 1 which provides an overview of multiple anatomical regions and in Fig. 2 which displays ankle images of a single patient. For each sequence, Table 1 provides the reported acquisition parameters, and Table 2 summarizes their characteristics and fields of study.

4Within the field of musculoskeletal imaging, T1‐ and proton density (PD)‐weighting is often acquired for structural bone imaging whereas T2‐weighting is acquired for imaging functional and pathophysiological processes. T1‐weighted (T1w) images have been acquired to detect structural lesions using spin‐echo (SE) or gradient‐echo (GRE) sequences. SE images and their derivatives are routinely acquired in musculoskeletal radiology owing to their excellent soft tissue contrast. Compared to GRE, SE sequences are also less prone to susceptibility, chemical shift, and field inhomogeneity artifacts but are nonetheless affected by geometrical distortions, especially at low receiver bandwidth and in regions far from the bore isocenter. On the other hand, GRE sequences are usually faster owing to a shorter minimal repetition time (TR) as shown in the pulse sequence chronograms in Fig. 3. In addition, GRE sequences are more versatile and are increasingly investigated for musculoskeletal radiology using radiofrequency spoiled gradient‐echo (S‐GRE also known as vendor‐specific acronyms FLASH [fast low angle shot], SPGR [spoiled gradient‐recalled], or T1‐FFE [T1 fast field echo]), or volumetric radiofrequency spoiled gradient‐echo (VS‐GRE also known as vendor‐specific acronyms VIBE [volumetric interpolated breath‐hold examination], LAVA [liver acquisition with volume acquisition], or THRIVE [T1‐weighted high‐resolution isotropic volume examination]) that enable post‐acquisition multiplanar reformatting.

5For bone morphology visualization, these sequences have been proposed in combination with fat‐saturation,, or water excitation to suppress the signal from adipose tissues and render bone as structures with a uniform low intensity. Alternatively, water‐only images can be generated by acquiring a GRE sequence with specific echo times to perform a Dixon water–fat separation.,,,,

6To further improve bone visualization, sequences have been developed to enhance bone specificity by providing a uniform soft tissue contrast, or by aiming at a CT‐like contrast.

7The gain in bone specificity was achieved in “black bone” (BB) sequences by applying a low flip angle and short echo time (TE) and TR to GRE‐like sequences, including VS‐GRE and ultrashort echo time (UTE) sequences. With such parameters, cortical bone appears as a low‐intensity structure whereas soft tissues have a uniform intermediate intensity. Originally developed for craniofacial imaging, BB‐MRI has been further applied to image the spine,, and long bones.

8Subsequently, the development of new hardware, which enabled faster transmit/receiving switching coils and more demanding gradients, permitted a drastic lowering of the echo time resulting in UTE sequences. In such sequences, the signal is usually acquired radially, soon after the end of the excitation, before a T2‐induced signal decay and minimal T2* signal decay. An image with a CT‐like contrast containing signal mainly in short T2 components can then be obtained by suppressing the long T2 signal. However, the fast gradient switching between TR's and the acquisition during the gradient ramp up (Fig. 3) renders UTE sequences prone to Eddy currents and susceptible to gradient delays, potentially resulting in imaging artifacts.

9With further developments, zero echo time (ZTE) images have been acquired for which the signal is sampled (usually radially) directly after the application of the radiofrequency (RF) pulse. To that end, readout gradients are on during the RF excitation. However, because of the delay in switching from transmit to receive modes, there is a dead time during which the center of the k‐space is not sampled (Fig. 3). Consequently, to reduce the dead time, hard short RF pulses need to be used, which put constraints on the achievable flip angles and bandwidths. ZTE images are acquired only with free induction decay readout, and a CT‐like contrast can be obtained by applying an inverse‐logarithmic rescaling. Since gradient switching is smooth (Fig. 3), ZTE acquisitions are rather silent, and the short achievable TR makes them fast.

10In addition to the advances in image acquisitions, image analysis and processing techniques have been applied to enhance bone visualization, usually aiming to create images with a CT‐like contrast. The simplest processing steps consisted in inverting the intensities or subtracting water intensities from the entire image, thus highlighting low signal in the MR images which is hypothesized to reflect the presence of cortical bone. Such a technique has been applied on standard GRE, VS‐GRE, or Dixon water‐only images., More advanced processing has been searched to convert MR image intensities to CT Hounsfield units (HU), creating so‐called synthetic CT (sCT). The most promising sCT generation models are deep learning‐based and rely on various network architectures, including UNet, generative adversarial network, and their derivatives. The use of sCT images has already been reviewed multiple times for radiotherapy purposes and PET‐MR,,, but their use for orthopedic purposes is rare. sCT generation models for orthopedic care have mainly been developed for the pelvis,,, sacroiliac joint, spine,, and long bones.,

11Three‐dimensional bone renderings are gaining popularity in orthopedic care as they provide an overview of the bone morphology, enable kinematic analyses,, and allow for the patient‐specific design of surgical guides and implants., Hence, 3D bone models facilitate the clinical diagnosis and improve surgical outcomes,,,,,, motivating their use in the treatment management of pathologies in the skull, shoulder, and hip.,,, Therefore, to be a CT surrogate for bone visualization, MRI should provide images on which bone can be segmented within a time and with a level of accuracy similar to or better than what can be achieved on CT. This section describes approaches for bone segmentation on MRI and provides results on segmentation geometrical accuracy and segmentation time with applications related to surgical planning.

12Regardless of the acquired MR contrast, there is a lack of MR‐dedicated, automated software for bone segmentation as exists for CT images. Bone segmentation on MR is mainly manual, or with extensive manual editing,,,,,,,,,,,, although some (semi‐)automated methods based on thresholding,, region growing, or ray casting can be applied. The development of fully automated segmentation approaches is complicated by structures in the vicinity of bones that share the same intensity as (cortical) bone and that can consequently be wrongly included in the bone segmentation. The problematic anatomical areas depend on the acquisition sequence but usually include air,,, and soft tissues like tendons, ligaments, or labrum.,, sCT images are a special case for segmentation as they are quantitative and reproduce HU from CT images. Hence, sCT can benefit from HU‐based segmentation and CT‐dedicated software,, as demonstrated by the segmentations of the knee bones obtained from S‐GRE, Dixon water‐only, sCT, and CT images in Fig. 4.

13In total, MR segmentation of bone lasted from 33 seconds to 5 hours,,,, in the reported literature and were made on standard of care fat‐suppressed MR, GRE‐MR, S‐GRE,,,, VS‐GRE,,, processed GRE‐derived,,,,,, BB,,,,, ZTE,, and sCT,,, images. This duration depends on the anatomy, the user's experience,,,, the segmentation method, and the desired quality of the segmentation, which hinders comparisons between studies. However, compared to CT‐based segmentation within the same study, segmentation on MR images was usually more time‐intensive,,,,, sometimes requiring more than twice the time.,, As an example, Fig. 5 presents timed segmentations of ankle bones obtained from CT and Dixon water‐only images. Nevertheless, when the segmentation was done by experts or companies, no difference was noted in segmenting bone from MR or CT images in terms of processing time., To alleviate the impact of user's experience on bone segmentation, automated methods based on deep learning are being developed, and are becoming commercially available for limited applications (eg, Mimics Innovation Suite 24, Materialize, Leuven, Belgium or CoLumbo, SmartSoft, Varna, Bulgaria).

14Regions prone to motion or magnetic field inhomogeneity could compromise the geometrical accuracy of the bone as seen on MR images. The geometrical integrity of the image can also be altered by nonlinear encoding gradients that may introduce compression or stretching of parts of the image. Because geometrical distortions could alter bone morphology, and consequently MR diagnostic capabilities, the overall geometry of the bone as visualized using MRI has been compared to physical ex vivo specimens and in vivo to CT, which is reviewed below.

15The geometrical accuracy of bone segmentation has been evaluated on long bones and vertebrae in ex vivo studies so that MR segmentation could be compared to the physical bone shape using 3D printing, mechanical contact/optical scanners,,,, or micro‐CT., Bone specimens were processed to remove soft tissues, resulting in a potential shrinkage of the gold standard compared to the bone as scanned using MRI and CT.,,,

16On average, CT segmentation overestimated the actual bone shape, whereas MR segmentation mostly underestimated it,,,, although not consistently., Nevertheless, surface distances between the MR‐based segmentation and the cadaveric specimen were on average submillimeter,,,,, with mean absolute surface distances ranging from 0.23 mm to 0.41 mm for MR‐segmentations and from 0.15 mm to 0.51 mm for CT‐based segmentations.,,, Similarly, root mean square error (RMSE) was mainly submillimeter,,, although it could reach 1.2 mm in the knee for MRI models (vs. 0.5 mm for CT models).

17When the CT‐based segmentation was used as a reference, the MR‐based bone segmentation also showed a submillimeter accuracy. In ex vivo long bones, absolute surface distances ranging from 0.23 mm to 0.61 mm were reported,,, with limits of agreement of the signed surface distance within ±0.72 mm, and RMSE of 1.1 mm. In ex vivo skulls, which can be harder to register and segment, BB‐MRI segmentations deviated on average by ±1.4 mm from CT segmentations.

18Larger differences were generally observed between the reference and MRI models near the joints in the proximal and distal bone ends,, although not always with statistical significance. Such differences resulted from the multiple soft tissues present at these locations (muscle, tendons, cartilage, and ligaments) which induce partial volume effects that hinder bone segmentation and warrant manual editing., Alternatively, errors were observed at the edge of the field of view (FOV), where there is less signal.

19The average submillimeter accuracy of bone segmentations compared to bone specimens shows that MR images have the ability to provide geometrically accurate bone models. To take into account soft tissue, evaluate more complex anatomies, and to make comparisons in an in vivo setting, MR segmentations were compared to CT segmentations. In in vivo hip joints, MR bone models differed on average by 0.4–0.9 mm from CT models,, with average RMSE under 1.8 mm for VS‐GRE Dixon images and under 0.81 mm for sCT images. When considering in vivo knees, there was no difference in the width and volume of the medial tibial plateau, with highly consistent measurements between standard of care PDw MR and CT images.

20This geometrical accuracy was influenced by the MR sequence acquired to perform the segmentation. Compared to other MR sequences, VS‐GRE offered the best correspondence to CT in the knee, with up to 45% differences in surface distance between VS‐GRE and balanced steady‐state GRE or spin‐echo derived images. VS‐GRE sequences had a better soft‐tissue‐to‐bone contrast, offering easier and more reproducible segmentations. Unfortunately, studies did not often report on registration parameters and acquisition parameters, such as the receiver bandwidth (see Table 1) or the built‐in distortion correction that can affect surface distance measurements and geometrical accuracy, preventing further comparisons between sequences and studies.

21Bone models obtained from segmentations can serve as a diagnostic tool in the therapeutic decision‐making, but also for surgical planning by allowing the design of customized surgical guides and implants.

22In the lower arm, saw guides for osteotomy were designed from CT and MR‐based sCT images and placed on cadaveric bones. The average saw guides positioning errors compared to the virtual planning were 2.4 mm and 3.8° for CT‐based guides and 2.8 mm and 4.9° for sCT‐based guides. More specifically, there were no intermodal statistically significant differences in the guides positioning. In addition, the intermodal rotational and translational limits of agreements were within the interobserver limits of agreement, suggesting the interchangeability of CT and sCT for the design of guides for long bone osteotomy. As an example, Fig. 6 shows saw guides positioning differences between CT and sCT.

23In the skull, deviations of ±1.4 mm were reported between BB‐MR and CT segmentations. The surgical guides resulting from the corresponding bone models were positioned on the skull with errors within ±0.6 mm for CT‐guides and ±0.8 mm for MR‐guides relative to their respective virtual planning. Given such differences, the average deviation from planned postoperative craniofacial reconstruction was within ±1.3 mm when using CT‐based guides and ±1.5 mm when using BB‐MRI‐based guides, with no statistical differences between the two modalities.

24Bone visualization and diagnosis on MR images can be hampered by the presence of water–fat interfaces, specific soft tissues like tendons, or air pockets in the vicinity of the bone since they may share the same low signal (Fig. 4). Therefore, the advantages and challenges of employing MRI for diagnosing bone‐related pathologies are anatomy‐specific. Multiple regions, including the skull, spine, shoulder, and pelvis have been assessed in recent years and are discussed in this section. For each anatomical region, the focus has been placed on two aspects: the potential of MRI for 1) detecting structural changes and for 2) measuring morphometric parameters of the bone. Structural changes include the detection of fractures, bone erosion, or sclerosis. On the other hand, morphometric parameters offer a quantitative assessment of bones which provides a standardized discrimination between “normal” and pathological regions and can influence therapeutic decision‐making.

25The development of MR protocols for skull visualization was favored by the routine acquisition of MR images for a wide range of clinical indications. In the standard of care, CT is indicated for trauma patients and for detecting osseous lesions, whereas MR images can be acquired for the detection of intracranial pathologies such as hemorrhage, ischemic changes, tumor, or other neurological disorders.,,,,,

26Skull composition and anatomy vary between stages of life, resulting in an age‐dependent diagnostic power of MRI and CT. Infants under 6 months have a thin skull (~1 mm thick) with high water content. In children under 2 years of age, sutures are wider, and harder to distinguish from trauma‐induced fractures., Lastly, cranial sutures tend to be less conspicuous in adults than children on BB‐MRI and UTE‐MRI as compared to CT because of a lower suture‐to‐skull contrast in adults.

27Nonetheless, the premature fusion of cranial sutures could be evaluated in infants and children with good to excellent inter‐ and intraobserver variability., For detecting skull fractures in children, BB‐MRI had an overall sensitivity of 66.7% and specificity of 87.5%, with errors originating from confusion between linear fractures and sutures, mainly in children under 2 years of age. However, the addition of BB‐MRI‐based 3D skull renderings increased the sensitivity to 83% and the specificity to 100% in a different cohort of patients under 30 months. Skull renderings obtained from MRI as shown in Fig. 7 were stated to be valuable for diagnosis in most patients, which is in line with results obtained in CT images. When UTE‐MR images were used, promising results were reported for the detection of fractures in patients aged from 1 month to 71 years. Sensitivity, specificity, and accuracy were all higher than 90%. The length and depth of the fractures could be measured on UTE‐MR images with no statistical difference compared to measurements made on CT images. For both UTE and BB‐MRI,, good to excellent inter and intraobserver agreement was reported for detecting fractures. Moreover, MRI could detect other pathologies such as edema, axonal injuries, and fractures accompanied by hemorrhages that were not visible on CT.,

28With regard to the temporomandibular joint and mandible, other structural differences were visible on UTE‐ and ZTE‐MRI with good to excellent inter‐ and intraobserver variability,, benefitting from a good UTE/ZTE‐to‐CT voxelwise intensity correlation in healthy and diseased bones., ZTE‐MRI revealed flattening and osteophytes in the mandibular condyles with near‐perfect agreement to cone beam CT whereas intermodal agreement for detecting medullary sclerosis was excellent on ZTE and moderate on UTE images., However, erosions, osteolysis, and periosteal reactions were more difficult to diagnose with only moderate intermodal agreement. In particular, periosteal reactions could be confused with air pockets on UTE images as both bone and air share similar intensities.

29Overall, all MR images suffered from misdiagnoses at interfaces between bone and air. Air and bone share similar intensities, making the bone/air interface difficult to distinguish.,,, Particularly difficult regions were the mastoid process,,, the paranasal sinuses,, and complex bone/fluid interfaces with high anatomical details like the inner ear. Such interfaces can cause misdiagnoses, and complicate automated processes for segmentation., To facilitate the distinction between tissues and air, phase information could complement the magnitude images, although the processing of phase images is complex and can be error prone.,

30Following the diagnosis, morphometric analysis of the mandible and cranium can be performed to plan craniofacial or maxillofacial surgeries., However, the measurement of such local parameters on MR images was not consistent between studies., On UTE‐based skull segmentations, intermodal differences of up to 2 mm and average deviation to cadaveric measurements of up to 4 mm were reported for eight anatomical parameters. By comparison, on SE‐derived images, performances were deemed statistically equivalent to cone beam CT, with average differences under 0.61 mm and 0.65° for 27 parameters and BB‐based bone segmentations deviated by ±1.4 mm from CT segmentations. Such differences might be due to differences in resolution as the voxel size was twice as small on T1w‐SE and BB‐MRI images as on UTE images (~0.5 mm vs. 1.1 mm).,, The T1w‐SE and BB acquisitions also had high receiver bandwidth (>610 Hz/px—Table 1) to maintain geometrical integrity.

31When imaging the spine, MRI is the modality of choice in many applications as it offers valuable information on the neural structures, the intervertebral discs, bone marrow, and the surrounding soft tissues., CT is typically acquired to assess the osseous involvement of soft tissue pathologies, for assessing bony abnormalities such as fractures, spondylosis, spondylolysis, or for surgical planning,,, owing to the superior cortical bone contrast and to the isotropic resolution of CT that enables multiplanar reformatting.

32When patients are suspected of having a vertebral fracture, CT images are routinely acquired to depict the extent of fracture. In addition, CT imaging is preferred in patients for whom a quick assessment is required, eg, patients who suffered high‐velocity accidents. Additional MRI is sometimes acquired, mainly to rule out occult injuries and to identify spinal cord lesions. In this context, MRI can also aid in distinguishing acute from old fractures, and can help diagnose specific types of fractures, such as stress fractures. Diagnostic performance statistics for detecting acute or stress fractures on MRI were excellent for S‐GRE, VS‐GRE, and UTE images with specificity, sensitivity, and accuracy above 90% when CT was used as ground truth.,, The interobserver agreement was good to excellent,, and was comparable between CT and MRI., For standard of care SE images, the specificity and accuracy for detecting fractures was above 95% while the sensitivity was 75% for incomplete fractures and 91% for complete fractures.

33Overall, GRE images, including S‐GRE and VS‐GRE seemed to outperform standard of care SE images for detecting fractures,, demonstrating a higher sensitivity in the delineation of the fracture line, probably owing to their thinner slices (see Table 1). Regarding S‐GRE and UTE images, differences are more questionable. In a cadaveric study, S‐GRE was reported as the most decisive standard of care sequence for detecting fractures in the pars interarticulares but observers were more confident in their diagnosis and missed fewer fractures when using UTE‐MRI. In particular, UTE imaging demonstrated a better interobserver agreement thanks to its CT‐like characteristics and a better contrast between the bone and fracture gap. On the other hand, in patients with suspected acute vertebral fractures, S‐GRE images outperformed UTE images in terms of the intermodal and interobserver agreements for detecting fractures, which was also the case for sclerosis, osteophytes, and joint degeneration. The difference in diagnostic quality between UTE and S‐GRE images between the cadaveric and in vivo studies might have several sources. In vivo, despite the radial k‐space sampling, the UTE images were reported to be prone to pulsation and motion artifacts that were not present ex vivo. In addition, the presence of multiple tissue types and air in the surrounding of the spine could result in susceptibility artifacts, especially seen in UTE images, more than GRE or ZTE images., A comparison between S‐GRE, UTE, and CT images for the detection of acute fractures and osteophytes is given in Fig. 8.

34SE, S‐GRE, VS‐GRE, and UTE sequences all misdiagnosed fractures in some patients,,, in part because of the misinterpretation of areas of bone sclerosis. These could be confused with subtle fractures or edema‐like changes,, or could mask fractures. However, the addition of a fluid‐sensitive MR sequence like short tau inversion recovery could reveal bone marrow edema and enable the detection of stress reactions that are invisible on CT images,, but can potentially change the patient's clinical management. For diagnosing sclerosis, GRE images seemed superior to UTE images or SE images as seen in the sacroiliac joint.,

35Other structural anomalies, including degenerative changes in the craniocervical junction and in the cervical spine that can cause neck pain were also investigated. In these cases, MRI is suitable for detecting ligamentous or intervertebral disc pathologies whereas CT can detect stenosis of the cervical spinal canal or neuroforamina. In both the craniocervical junction and cervical spine, degenerative changes were graded with good intermodal agreement,, similar to the interobserver agreement on CT. This was facilitated by multiplanar reformatting possible on isotropic ZTE images. Overall, good to excellent inter‐ and intraobserver agreement was reported with MRI., However, when using MR images with inverted intensities, care needs to be taken not to misinterpret the apparent high signal intensity of ligaments, or of the gas accumulation in the intervertebral discs as calcifications as seen in Fig. 8.

36The diagnosis of degenerative changes can also be made quantitatively by measuring morphometric parameters including vertebral body and intervertebral disc parameters. Despite a good to excellent intermodal agreement in measuring vertebral body height on UTE, S‐GRE, and sCT images, and an excellent interobserver agreement,, the accuracy of morphological vertebral assessment was highly dependent on the MR acquisition. Average differences in the vertebral height of 0.26 mm were reported in the sCT with limits of agreement within ±2 mm for S‐GRE images and within 6–10 mm for UTE images. For intervertebral disc heights, limits of agreement of ±2–3 mm were reported between S‐GRE and CT images and of ±4 mm between UTE and CT images. Similarly, CT/UTE intermodal limits of agreements were within ±1 mm for the distance between the cranium and C1 and within ±2–4 mm between the cranium and C2. However, these intermodal differences in the distances between the cranium and cervical spine were not significant and may partly originate from the differences in resolution between MR (0.8 mm × 0.8 mm × 1.2 mm) and CT (<0.6 mm × 0.6 mm × 0.6 mm).

37In the shoulder, MR examinations are commonly performed to examine the ligaments, the rotator cuff, the labrum, and the joint capsule, eg, after shoulder dislocation.,, However, standard of care T1w‐SE images have similar low intensity for cortical bone and labrum, warranting a CT examination to assess the glenohumeral bone architecture and review bone changes. In particular, the amount of glenoid bone loss, and to a lesser extent humeral deformity, often associated with shoulder dislocation, determines the clinical management plan.

38On MR images acquired with, or without intra‐articular contrast injection, a strong correlation (r > 0.8) was found between MR and CT for the glenoid width and percentage bone loss. The mean glenoid bone loss error was under 2.5% for both modalities, with intermodal differences not statistically different., Furthermore, MR and CT measurements had good correlation with arthroscopy as percentage bone loss differences under 3% were reported between MR and arthroscopy and under 1% between CT and arthroscopy.

39When MR‐based 3D bone renderings were compared to CT, bone defects were also equivalently visible on the bone reconstructions with no statistically significant differences between MR‐ and CT‐based renderings,,, good to excellent intermodal correlation,, excellent intermodal agreement,, and submillimeter/<1% average defect size difference.,, Although small on average, some intermodal differences could reach up to 3 mm/10% difference in glenoid bone loss,, potentially influencing clinical management for a minority of patients. However, such large differences were not systematically reported, with some maximal differences within ±7.5%.,,

40In addition to diagnosing bone loss, MR was used to detect fractures in the humerus and scapula with good intermodal agreement, and excellent sensitivity and specificity (>90%), using ZTE‐MR and VS‐GRE images. Fracture extent was measured equivalently on MR and CT. ZTE‐MR images were able to reveal bone depression, bone resorption, and bone fragments better than standard of care PDw images in most patients, along with a good ZTE‐to‐CT intermodal agreement for detecting bone fragments and osteoarthritis. Moreover, ZTE‐MRI surpassed CT in revealing cortical bone and intraosseous lesions within a single image. In particular, bone marrow edema and cysts that remained undetected in CT were visible on ZTE images. As cysts indicate regions of lower bone quality, it is important to accurately detect them, particularly when the images are used for guidance of surgical planning.

41Treatment planning of shoulder instability might include the measurement of morphometric parameters of the glenoid on 2D images or 3D renderings. The glenoid morphometric accuracy on MRI was comparable to CT as demonstrated by an excellent intermodal agreement with no statistical difference in the measurement on glenoid vault and glenoid version angle,,, using ZTE or standard of care MR images., Intermodal agreement was good when comparing certain shoulder‐specific parameters, with limits of agreements within 6 mm for measuring glenoid vault depth and within 5° for the version angle, for most patients. However, for some patients, these measurements could differ drastically due to blurring and reduced FOVs on the MR images., The intermodal limits of agreements were within the interobserver limits of agreement.

42In a similar way, the geometrical accuracy of MR‐based 3D bone renderings was compared to CT. Measuring glenoid/humeral width, height, and surface areas was equivalent between CT and MRI, although some statistically significant differences could be found but not systematically., In particular, the average intermodal differences in glenoid and humeral surfaces were within ±10%,, and in glenoid/humeral width and height were within ±1 mm.,

43Based on bone rendering and radiography, kinematic analysis of joints can also be performed to quantify changes to the joint position and contact points., Hence, using S‐GRE or CT‐like VS‐GRE images, digitally reconstructed radiographs (DRR) have been generated and registered to radiography images for shoulder and knee kinematic analyses., Registration errors and kinematics measurements errors were larger for MR‐based DRRs than CT‐based DRR and showed an RMSE under 2.2 mm and 2.6° (vs. 1.6 mm and 2.2° for CT)., Because CT and radiography share the same contrast mechanism, intensity‐based metrics can be used for the CT‐DRR to radiograph registration whereas only edge‐information was used for the MRI‐DRR to radiograph registration, potentially explaining the larger registration errors of MRI‐based DRRs.

44Overall, MR and CT showed good to excellent inter‐ and intraobserver variability in diagnosing bone pathologies or performing morphometric measurements in the shoulder,,,,,, with good to excellent intermodal agreement.,,, Figure 9 compares the CT image of a shoulder to the corresponding ZTE images acquired at different resolutions.

45The pelvic bone connects the upper body to the lower limbs through the sacroiliac and hip joints. Both joints can be subjected to degenerative osteoarthritic changes, affecting the bone and the surrounding soft tissues. In the sacroiliac joint, spondyloarthritis induces bone marrow edema and inflammatory lesions that can be detected with MRI, and structural lesions such as erosions, sclerosis, or ankylosis that may be detected with MRI,, but are better defined on CT.,, In the hip joint, hip dysplasia and femoroacetabular impingement are morphological hip conditions that affect bone, and soft tissues including but not limited to cartilage and labrum. For these conditions, clinical care usually includes morphometric assessment of the joint made on radiograph, with optional addition of CT or MRI for diagnosis, a soft tissue evaluation with MRI and a bone rendering based on CT for surgical planning. To limit adverse ionizing radiation, MR has been investigated as a diagnostic tool for detecting structural bone lesions and performing bone morphometric assessments and 3D renderings.

46In the sacroiliac joint, SE, GRE, and sCT images were used to assess structural changes. T1w‐SE images were shown to strongly correlate with low‐dose CT for detecting erosions and were able to reveal 88% of erosions. However, standard MRI missed some cases of axial spondyloarthritis when used alone. On more dedicated images to visualize bone, including VS‐GRE or sCT images, higher diagnostic accuracy and diagnostic confidence were achieved for detecting erosions,,, especially when sclerosis was present. In particular, the sensitivity and specificity for detecting erosions increased between standard and dedicated images, reaching a sensitivity above 70%,,, a specificity around 90%,,, and an accuracy above 90%. A qualitative comparison between T1w, VS‐GRE, and CT images for diagnosing erosions is presented in Fig. 10.

47Standard MR images were also able to detect 92% of joint space alterations, and to a lesser extent sclerosis, while sCT could diagnose sclerosis and ankylosis with accuracies higher than 90%. Dedicated MR imaging was also reliable owing to a good to excellent interobserver agreement, and repeatable with good intraobserver agreement, in accordance with CT imaging., Moreover, observers were more confident when scoring VS‐GRE images than low‐dose CT images because of the noise of low‐dose CT and equally confident when scoring sCT compared to CT images. Some erosions were only visible on dedicated MR images and not on CT images,, especially in young patients and for small lesions. As no erosions were detected in healthy controls, this suggests that the observed destructive changes were no artifact and that MR was superior in revealing those erosions.

48In addition, MR images have been acquired to describe fractures in the hip joint. MRI has been shown to perform better than CT for the detection of fractures in the hip in elderly patients,, but also in adolescents and children.,, All fractures detected on CT were also detected on MR, and some fractures were visible on MR but not CT, or misdiagnosed on CT resulting in changes in the clinical management plan, especially regarding instructions for weight bearing. In patients under the age of 13 years, the posterior acetabular wall is not fully ossified and the MRI findings of traumatic hip dislocations with acetabular fractures were better correlated with intraoperative findings than CT findings, which did not always directly detect acetabular fractures. Some soft tissue defects, oblivious to CT, were also identified on MRI. These included entrapment of labra and posterior acetabular cartilage fractures. The detection of entrapment of labra, in particular, had an influence on the patient clinical management.

49The measurement of morphometric parameters is especially important in the hip joint for the diagnosis of hip dysplasia and femoroacetabular impingement and has been investigated on VS‐GRE,, ZTE, sCT, or intermediate‐weighted images. All imaging techniques found good to excellent intermodal agreement for measuring the acetabular version,,, with excellent inter‐ and intraobserver agreements,, and statistically equivalent measures given an acceptable error below 4.3°., The reported limits of agreements were, however, mixed. On VS‐GRE and sCT images, intermodal limits of agreements were in line with the intraobserver variability, within ±4.2°., On ZTE images, on the other hand, intermodal limits of agreements of acetabular version reached 11.3°, higher than the 8° obtained for the interobserver variability on CT. These differences might originate from the fact that pelvic tilt was not standardized in the ZTE study. Other parameters that were compared include the lateral center edge and alpha angles. The intermodal agreement was good to excellent, with intermodal limits of agreements roughly within 12°,, and bounded by the interobserver limits of agreement achieved on CT.

50Femoral parameters such as the femoral anteversion were also measured and compared between CT and MRI. In two studies using standard clinical sequences,, a strong correlation between CT and MR measurements was reported with correlation coefficients of 0.77 and 0.80 between the two modalities. However, the intermodal absolute agreement was poor with biases ranging from 5° to 10°, probably because of interscan positioning differences. The MR examination being long (30–45 minutes), patients might be given knee wedges to bend the knee, or can relax into greater external rotation of the hip, possibly explaining such differences. When measured in infants with developmental dysplasia of the hip on CT and T1w‐SE images, the intermodal, intraobserver, and interobserver agreements for the femoral version were all excellent (intraclass correlation coefficient >0.9), demonstrating the reliability and reproducibility of the methods. In such a young population, MRI had the advantage over CT that it was able to visualize the not fully ossified femoral condyle in infants under the age of 6 months. In these cases, the condylar plane could be defined more accurately on MR than CT.

51Using bone renderings, the measurement of local morphometric parameters, including the center‐edge angle and acetabular version was similar between CT and MR, with average intermodal differences under 4°.,, Hip range of motion measurements were also compared between CT and MR with average differences under 4° for all rotations and with limits of agreements within ±6°. All measurements had excellent intermodal correlation, intermodal agreement, and interobserver agreement., Correspondingly, such models were able to diagnose femoroacetabular impingement or hip dysplasia with 100% agreement reported between MR and CT for the presence and location of cam deformity, and good to excellent intermodal and interobserver agreements. Figure 11 presents 3D bone reconstructions as obtained from CT and sCT.

52Overall, the use of MRI as a radiation‐free alternative to CT for bone visualization has received a lot of attention in the last decade and has been valued by multiple editorials,, and overviews., Although promising, MRI does suffer from challenges related to data acquisition and accessibility to novel technologies. This section describes these challenges and discusses how they might affect the adoption of MRI for bone imaging in clinical practice.

53When planning an MRI scan, a trade‐off has to be made between the FOV, the resolution, the signal‐to‐noise ratio (SNR), and the acquisition time. The MRI sequences used for bone visualization pose different constraints to this trade‐off which may limit their applicability in specific situations. These constraints are related to several external factors, including the size of the region to be imaged, or also the intrinsic tissue‐specific factors like magnetic susceptibility. In this section, we address these aspects in relation to their use for bone visualization.

54In general, a limited FOV is chosen when planning MRI acquisitions to reduce scan time. However, some sequences like UTE or ZTE have spatially nonselective excitations which induce large fields of view and reduce their flexibility regarding other acquisition parameters (eg, spatial resolution) to keep a reasonable scan time.

55For other sequences, the freedom to reduce the FOV has two constraints. First, the FOV must be large enough to make a proper diagnosis. This includes the visualization of landmarks for post‐acquisition image standardization (eg, correction of the anterior pelvic plane for hip imaging) and of a sufficiently large region for the measurement of morphometric parameters (eg, scapula for measuring the glenoid version in the shoulder and femoral shaft and condyles for measuring the femoral neck shaft angle, ). Second, care needs to be taken to avoid the edge of the FOV where the lower signal and field inhomogeneity may compromise the measurements. As a solution, MR images can be acquired in multiple blocks, overlapping or not, to obtain the necessary information.,,,, Such multi‐station acquisitions are, however, susceptible to slight changes in position between individual acquisitions that can compromise the geometric integrity of the bone.

56Image resolution was often lower on MRI than on CT, with voxel sizes usually ranging between 0.6 mm and 1 mm in the reviewed literature (see Table 1). Note that some of the reported resolutions are reconstructed resolutions and not acquired resolutions. Low‐resolution images induce more partial volume effects that can mask or, on the contrary, enlarge structures of interest, potentially resulting in the under‐ or over‐segmentation of bone on MR images. Furthermore, for 3D bone modeling, low resolutions result in high interpolation uncertainty and can cause stair‐step artifacts.,, However, increasing the resolution is not always beneficial as it is accompanied by a decrease in SNR or an increase in acquisition time without necessarily improving the diagnostic capabilities of the images.,,

57MR images usually had an (almost) isotropic resolution in the literature assessed for this review (Table 1). The voxel isotropy enables multiplanar reformatting of the images for an improved visualization of the vertebrae, of the glenoid, or of the femoral neck for the measurement of morphometric parameters in these regions. In addition, for 3D bone modeling, voxel isotropy makes the interpolation uncertainty equal in all directions, facilitating bone models interpretation and the subsequent modeling of surgical tools.

58Long acquisition times are problematic as they induce higher costs and potential motion artifacts. In children, in particular, motion artifacts could compromise the diagnostic quality of the MR images., Voluntary motion in the youngest patients can be avoided, by either using immobilization, or sedation which includes the feed‐and‐sleep method or general anesthesia., Although also sometimes required for CT acquisition, deeper sedation is usually needed during MRI acquisition because of the longer acquisition time. The use of anesthesia is however not risk‐free, especially on repeated occurrences. In addition, uncontrolled motion was seen in the spine, and the jaw, and is common in clinical care of the shoulder, weakening the diagnostic power of MRI. Nonetheless, motion artifacts can be reduced by using motion insensitive acquisition methods, including breath holds, interleaved scanning, increased parallel imaging with higher signal averaging, or radial sampling of the k‐space.,,,

59MRI can also be impaired by magnetic susceptibility‐related distortions due to the shape of the body, or the presence of air or of implanted devices. Areas of concern for such artifacts are the spine, where the bone is surrounded by multiple magnetically differing soft tissues, air and/or metal instrumentation, the jaw which can contain orthodontic devices,, but also long bones with screw fixations or the skull with ventriculoperitoneal shunts. When expected, susceptibility artifacts can be partly mitigated by choosing the adequate MR sequence and acquisition parameters. At equivalent acquisition parameters, GRE sequences are more prone to susceptibility artifacts than SE sequences, and sequences such as UTE are more prone to field inhomogeneity artifacts than S‐GRE or ZTE. The geometrical distortions induced by susceptibility artifacts can be mitigated by increasing the receiver bandwidth at the cost of SNR, by applying the scanner's built‐in distortion correction, or by limiting the FOV around the scanner's isocenter, but they are never completely removed. In addition, although 3 T acquisitions are usually equivalent or better than 1.5 T acquisitions for bone visualization and segmentation,, lower field acquisitions should be favored when inhomogeneity artifacts are expected., Low‐field MRI (<0.5 T) in particular could be acquired to diagnose pathologies associated with orthopedic hardware, given the assumption that low‐field MRI is not overly impacted by susceptibility artifacts and is able to image soft tissues in the vicinity of the implant. Other advantages of low‐field MRI include its low cost (purchase and maintenance), and, when considering musculoskeletal radiology, the ability to scan in weight‐bearing position. This, however, comes at the expanse of SNR and resolution.

60Compared to CT, MRI suffers from a multitude of contraindications that make it unavailable for some patients. For trauma patients, the access to an MRI can be limited by obstacles related to diagnostic speed, transport of the patient to the MRI, MRI incompatibility with life‐support or monitoring equipment, and patient implants. Metallic MR‐compatible devices are problematic when in the vicinity of the region of interest as they can generate susceptibility artifacts hampering the diagnosis. Devices that are not MR‐compatible, including some pacemakers and cochlear implants, preclude any MR acquisition. In addition, claustrophobic patients or patients unable to stay motionless might require sedation to undergo MRI, complicating the workflow, potentially causing adverse effects, and hindering compliance with the breathing instructions required for some sequences. Overall, in an emergency department, more than a quarter of the elderly patients coming after a trauma could have at least one contraindication for MRI.

61Another issue of MRI is its availability for acquisition. CT being faster, it is more accessible, especially in cases of emergency., When an MR system is acquired, the choice of the MR sequence might be driven by the available hardware. The sequences described in this review are not commonly present on all scanners. Dixon reconstruction is now usually built in the scanner, but sequences like UTE and ZTE might require modern hardware or specific chargeable licensing. They tend to be increasingly available and offered as standard sequences but as an example, all ZTE images presented in this review were obtained only on GE scanners. Tools for sCT generation from GRE‐derived images are also becoming commercially available.

62CT is considered the modality of choice for visualizing cortical bone in 3D. However, its adverse radiation burden has motivated the research into alternative modalities with lower radiation doses, including radiography‐based, low‐dose CT‐based, or MR‐based methods. In this competition, despite some challenges in the acquisition, MRI has favorable properties including its superior soft tissue contrast that can be exploited to concurrently assess the soft tissue involvement of musculoskeletal pathologies without the need for image registration, and the complete absence of ionizing radiation.

63The utility of MRI as an imaging modality for visualizing bone has been shown in many areas of the human body. Overall, the CT‐to‐MR intermodal agreement for the diagnosis of osseous pathologies and for the measurement of anatomical parameters was good to excellent with multiple reports of statistical equivalence.,,, In addition, MRI could provide 3D bone renderings, critical in the clinical care for the skull, shoulder, or hip, with a submillimeter accuracy compared to CT, although in general representing an underestimation of the actual bone size.

64MRI presented several advantages compared to CT in the diagnosis of musculoskeletal pathologies. First, immature bone as seen in the femur and pelvis of young children was better visualized on MR than on CT images., Second, MRI can acquire soft tissue and bone information in a single examination. Sequences like VIBE, UTE, or multi‐echo steady state (MESS) can provide bone structural information while providing complementary information on other tissue, including cartilage.,,, This can promote joint biomechanical and kinematic modeling by limiting the need for registration.,, Some of the dedicated images also revealed fractures and lesions, like cysts or edema, that can improve patients' clinical management but which were not visible on CT. Furthermore, MRI can be used to generate simulated radiographs with a diagnostic quality similar to CT for imaging bone tumors, while providing additional information on tumor architecture and soft tissue extension. Lastly, MR sequences for imaging bone can also be combined with other sequences for specific imaging such as venous or fluid‐sensitive, imaging, magnetic resonance angiography, or quantitative susceptibility mapping for a more comprehensive diagnosis within a single modality.

65The use of an MR‐based bone visualization could in the future be extended to facilitate clinical care motivated by the benefits of CT/MR fusion. Such fusions could be useful for the design of patient‐specific implants by combining bone and joint capsule information and have proven their potential for diagnostic and treatment purposes, by easing the diagnosis for junior readers, and by facilitating treatment guidance, and surgical navigation.,,, However, fusing MR and CT requires an intermodal registration that is not necessarily straightforward. With MR providing a CT‐like visualization of bone, visualization of soft tissues and bones can be obtained in one scanning session with similar body geometry, offering new perspectives for diagnosis, treatment planning, and guidance. Fusion of MRI with radiography has also been performed for kinematic analysis.,

66Is one sequence better than the other? VS‐GRE and ZTE sequences seemed to stand out with validation in multiple anatomies, owing to their rather fast acquisition offering isotropic images with good cortical bone‐to‐bone marrow and cortical bone‐to‐muscles contrasts that facilitate bone segmentation. VS‐GRE sequences are also robust to respiratory motion through breath holds while ZTE acquisitions are robust to motion in general thanks to their radial k‐space sampling. In parallel, sCT is gaining interest for orthopedics,, building upon its CT‐like HU, although care still needs to be taken in interpreting such artificial intelligence‐based images. In addition, the validity and robustness of single sCT generation models need to be carefully assessed across multiple MR vendors and sites. In general, awareness of the possible artifacts and MR image specificities, especially regarding air, ligaments, tendons, or water/fat interfaces, is required for all anatomical regions and MR sequences.,,,,,, However, getting acquainted with the use of MR images for measurements and diagnosis might be easier and faster on images with a CT‐like contrast, like ZTE or sCT, that have a high correspondence to CT images.,,, As of now, only a few studies compared multiple MR sequences,,,,,,, with equivalent acquisition parameters to CT, complicating definite conclusions, which might be specific to an anatomical region.

67To conclude, MRI is a promising radiation‐free alternative to CT for the diagnosis and treatment planning of bone pathologies. The recent advances in hardware and software provide MR images with a spatial resolution and contrast that are similar to CT images for the detection of structural and degenerative bone changes. MRI will probably not replace CT for all its applications in the near future, especially not in emergency settings. However, for clinical indications where both bone and soft tissue information are required, these new approaches open new perspectives for comprehensive protocols that facilitate bone and soft tissue visualization and fusion, for diagnosis, treatment planning, and surgical guidance.