Deformable motion compensation in interventional cone-beam CT with a context-aware learned autofocus metric.

PURPOSE: Interventional Cone-Beam CT (CBCT) offers 3D visualization of soft-tissue and vascular anatomy, enabling 3D guidance of abdominal interventions. However, its long acquisition time makes CBCT susceptible to patient motion. Image-based autofocus offers a suitable platform for compensation of deformable motion in CBCT, but it relies on handcrafted motion metrics based on first-order image properties and that lack awareness of the underlying anatomy. This work proposes a data-driven approach to motion quantification via a learned, context-aware, deformable metric, VI F D L ${\bm{VI}}{{\bm{F}}}_{DL}$ , that quantifies the amount of motion degradation as well as the realism of the structural anatomical content in the image. METHODS: The proposed VI F D L ${\bm{VI}}{{\bm{F}}}_{DL}$ was modeled as a deep convolutional neural network (CNN) trained to recreate a reference-based structural similarity metric-visual information fidelity (VIF). The deep CNN acted on motion-corrupted images, providing an estimation of the spatial VIF map that would be obtained against a motion-free reference, capturing motion distortion, and anatomic plausibility. The deep CNN featured a multi-branch architecture with a high-resolution branch for estimation of voxel-wise VIF on a small volume of interest. A second contextual, low-resolution branch provided features associated to anatomical context for disentanglement of motion effects and anatomical appearance. The deep CNN was trained on paired motion-free and motion-corrupted data obtained with a high-fidelity forward projection model for a protocol involving 120 kV and 9.90 mGy. The performance of VI F D L ${\bm{VI}}{{\bm{F}}}_{DL}$ was evaluated via metrics of correlation with ground truth VIF ${\bm{VIF}}$ and with the underlying deformable motion field in simulated data with deformable motion fields with amplitude ranging from 5 to 20 mm and frequency from 2.4 up to 4 cycles/scan. Robustness to variation in tissue contrast and noise levels was assessed in simulation studies with varying beam energy (90-120 kV) and dose (1.19-39.59 mGy). Further validation was obtained on experimental studies with a deformable phantom. Final validation was obtained via integration of VI F D L ${\bm{VI}}{{\bm{F}}}_{DL}$ on an autofocus compensation framework, applied to motion compensation on experimental datasets and evaluated via metric of spatial resolution on soft-tissue boundaries and sharpness of contrast-enhanced vascularity. RESULTS: The magnitude and spatial map of VI F D L ${\bm{VI}}{{\bm{F}}}_{DL}$ showed consistent and high correlation levels with the ground truth in both simulation and real data, yielding average normalized cross correlation (NCC) values of 0.95 and 0.88, respectively. Similarly, VI F D L ${\bm{VI}}{{\bm{F}}}_{DL}$ achieved good correlation values with the underlying motion field, with average NCC of 0.90. In experimental phantom studies, VI F D L ${\bm{VI}}{{\bm{F}}}_{DL}$ properly reflects the change in motion amplitudes and frequencies: voxel-wise averaging of the local VI F D L ${\bm{VI}}{{\bm{F}}}_{DL}$ across the full reconstructed volume yielded an average value of 0.69 for the case with mild motion (2 mm, 12 cycles/scan) and 0.29 for the case with severe motion (12 mm, 6 cycles/scan). Autofocus motion compensation using VI F D L ${\bm{VI}}{{\bm{F}}}_{DL}$ resulted in noticeable mitigation of motion artifacts and improved spatial resolution of soft tissue and high-contrast structures, resulting in reduction of edge spread function width of 8.78% and 9.20%, respectively. Motion compensation also increased the conspicuity of contrast-enhanced vascularity, reflected in an increase of 9.64% in vessel sharpness. CONCLUSION: The proposed VI F D L ${\bm{VI}}{{\bm{F}}}_{DL}$ , featuring a novel context-aware architecture, demonstrated its capacity as a reference-free surrogate of structural similarity to quantify motion-induced degradation of image quality and anatomical plausibility of image content. The validation studies showed robust performance across motion patterns, x-ray techniques, and anatomical instances. The proposed anatomy- and context-aware metric poses a powerful alternative to conventional motion estimation metrics, and a step forward for application of deep autofocus motion compensation for guidance in clinical interventional procedures.

Super-resolution of clinical CT: Revealing microarchitecture in whole bone clinical CT image data.

Osteoporotic fractures, prevalent in the elderly, pose a significant health and economic burden. Current methods for predicting fracture risk, primarily relying on bone mineral density, provide only modest accuracy. If better spatial resolution of trabecular bone in a clinical scan were available, a more complete assessment of fracture risk would be obtained using microarchitectural measures of bone (i.e. trabecular thickness, trabecular spacing, bone volume fraction, etc.). However, increased resolution comes at the cost of increased radiation or can only be applied at small volumes of distal skeletal locations. This study explores super-resolution (SR) technology to enhance clinical CT scans of proximal femurs and better reveal the trabecular microarchitecture of bone. Using a deep-learning-based (i.e. subset of artificial intelligence) SR approach, low-resolution clinical CT images were upscaled to higher resolution and compared to corresponding MicroCT-derived images. SR-derived 2-dimensional microarchitectural measurements, such as degree of anisotropy, bone volume fraction, trabecular spacing, and trabecular thickness were within 16 % error compared to MicroCT data, whereas connectivity density exhibited larger error (as high as 1094 %). SR-derived 3-dimensional microarchitectural metrics exhibited errors <18 %. This work showcases the potential of SR technology to enhance clinical bone imaging and holds promise for improving fracture risk assessments and osteoporosis detection. Further research, including larger datasets and refined techniques, can advance SR's clinical utility, enabling comprehensive microstructural assessment across whole bones, thereby improving fracture risk predictions and patient-specific treatment strategies.

XIOSIS: an X-ray-based intra-operative image-guided platform for oncology smart material delivery.

Image-guided interventional oncology procedures can greatly enhance the outcome of cancer treatment. As an enhancing procedure, oncology smart material delivery can increase cancer therapy's quality, effectiveness, and safety. However, the effectiveness of enhancing procedures highly depends on the accuracy of smart material placement procedures. Inaccurate placement of smart materials can lead to adverse side effects and health hazards. Image guidance can considerably improve the safety and robustness of smart material delivery. In this study, we developed a novel generative deep-learning platform that highly prioritizes clinical practicality and provides the most informative intra-operative feedback for image-guided smart material delivery. XIOSIS generates a patient-specific 3D volumetric computed tomography (CT) from three intraoperative radiographs (X-ray images) acquired by a mobile C-arm during the operation. As the first of its kind, XIOSIS (i) synthesizes the CT from small field-of-view radiographs;(ii) reconstructs the intra-operative spacer distribution; (iii) is robust; and (iv) is equipped with a novel soft-contrast cost function. To demonstrate the effectiveness of XIOSIS in providing intra-operative image guidance, we applied XIOSIS to the duodenal hydrogel spacer placement procedure. We evaluated XIOSIS performance in an image-guided virtual spacer placement and actual spacer placement in two cadaver specimens. XIOSIS showed a clinically acceptable performance, reconstructed the 3D intra-operative hydrogel spacer distribution with an average structural similarity of 0.88 and Dice coefficient of 0.63 and with less than 1 cm difference in spacer location relative to the spinal cord.

Deformable registration of preoperative MR and intraoperative long-length tomosynthesis images for guidance of spine surgery via image synthesis.

PURPOSE: Improved integration and use of preoperative imaging during surgery hold significant potential for enhancing treatment planning and instrument guidance through surgical navigation. Despite its prevalent use in diagnostic settings, MR imaging is rarely used for navigation in spine surgery. This study aims to leverage MR imaging for intraoperative visualization of spine anatomy, particularly in cases where CT imaging is unavailable or when minimizing radiation exposure is essential, such as in pediatric surgery. METHODS: This work presents a method for deformable 3D-2D registration of preoperative MR images with a novel intraoperative long-length tomosynthesis imaging modality (viz., Long-Film [LF]). A conditional generative adversarial network is used to translate MR images to an intermediate bone image suitable for registration, followed by a model-based 3D-2D registration algorithm to deformably map the synthesized images to LF images. The algorithm's performance was evaluated on cadaveric specimens with implanted markers and controlled deformation, and in clinical images of patients undergoing spine surgery as part of a large-scale clinical study on LF imaging. RESULTS: The proposed method yielded a median 2D projection distance error of 2.0 mm (interquartile range [IQR]: 1.1-3.3 mm) and a 3D target registration error of 1.5 mm (IQR: 0.8-2.1 mm) in cadaver studies. Notably, the multi-scale approach exhibited significantly higher accuracy compared to rigid solutions and effectively managed the challenges posed by piecewise rigid spine deformation. The robustness and consistency of the method were evaluated on clinical images, yielding no outliers on vertebrae without surgical instrumentation and 3% outliers on vertebrae with instrumentation. CONCLUSIONS: This work constitutes the first reported approach for deformable MR to LF registration based on deep image synthesis. The proposed framework provides access to the preoperative annotations and planning information during surgery and enables surgical navigation within the context of MR images and/or dual-plane LF images.

Feasibility of bone marrow edema detection using dual-energy cone-beam computed tomography.

BACKGROUND: Dual-energy (DE) detection of bone marrow edema (BME) would be a valuable new diagnostic capability for the emerging orthopedic cone-beam computed tomography (CBCT) systems. However, this imaging task is inherently challenging because of the narrow energy separation between water (edematous fluid) and fat (health yellow marrow), requiring precise artifact correction and dedicated material decomposition approaches. PURPOSE: We investigate the feasibility of BME assessment using kV-switching DE CBCT with a comprehensive CBCT artifact correction framework and a two-stage projection- and image-domain three-material decomposition algorithm. METHODS: DE CBCT projections of quantitative BME phantoms (water containers 100-165 mm in size with inserts presenting various degrees of edema) and an animal cadaver model of BME were acquired on a CBCT test bench emulating the standard wrist imaging configuration of a Multitom Rax twin robotic x-ray system. The slow kV-switching scan protocol involved a 60 kV low energy (LE) beam and a 120 kV high energy (HE) beam switched every 0.5° over a 200° angular span. The DE CBCT data preprocessing and artifact correction framework consisted of (i) projection interpolation onto matched LE and HE projections views, (ii) lag and glare deconvolutions, and (iii) efficient Monte Carlo (MC)-based scatter correction. Virtual non-calcium (VNCa) images for BME detection were then generated by projection-domain decomposition into an Aluminium (Al) and polyethylene basis set (to remove beam hardening) followed by three-material image-domain decomposition into water, Ca, and fat. Feasibility of BME detection was quantified in terms of VNCa image contrast and receiver operating characteristic (ROC) curves. Robustness to object size, position in the field of view (FOV) and beam collimation (varied 20-160 mm) was investigated. RESULTS: The MC-based scatter correction delivered > 69% reduction of cupping artifacts for moderate to wide collimations (> 80 mm beam width), which was essential to achieve accurate DE material decomposition. In a forearm-sized object, a 20% increase in water concentration (edema) of a trabecular bone-mimicking mixture presented as ∼15 HU VNCa contrast using 80-160 mm beam collimations. The variability with respect to object position in the FOV was modest (< 15% coefficient of variation). The areas under the ROC curve were > 0.9. A femur-sized object presented a somewhat more challenging task, resulting in increased sensitivity to object positioning at 160 mm collimation. In animal cadaver specimens, areas of VNCa enhancement consistent with BME were observed in DE CBCT images in regions of MRI-confirmed edema. CONCLUSION: Our results indicate that the proposed artifact correction and material decomposition pipeline can overcome the challenges of scatter and limited spectral separation to achieve relatively accurate and sensitive BME detection in DE CBCT. This study provides an important baseline for clinical translation of musculoskeletal DE CBCT to quantitative, point-of-care bone health assessment.

Macroscopic differences in adult human femora are linked to body mass index.

Bone functional adaptation is routinely invoked to interpret skeletal morphology despite ongoing debate regarding the limits of the bone response to mechanical stimuli. The wide variation in human body mass presents an opportunity to explore the relationship between mechanical load and skeletal response in weight-bearing elements. Here, we examine variation in femoral macroscopic morphology as a function of body mass index (BMI), which is used as a metric of load history. A sample of 80 femora (40 female; 40 male) from recent modern humans was selected from the Texas State University Donated Skeletal Collection. Femora were imaged using x-ray computed tomography (voxel size ~0.5 mm), and segmented to produce surface models. Landmark-based geometric morphometric analyses based on the Coherent Point Drift algorithm were conducted to quantify shape. Principal components analyses were used to summarize shape variation, and component scores were regressed on BMI. Within the male sample, increased BMI was associated with a mediolaterally expanded femoral shaft, as well as increased neck-shaft angle and decreased femoral neck anteversion angle. No statistically significant relationships between shape and BMI were found in the female sample. While mechanical stimulus has traditionally been applied to changes in long bong diaphyseal shape it appears that bone functional adaptation may also result in fundamental changes in the shape of skeletal elements.

Improving model-data mismatch for photon-counting detector model using global and local model parameters.

BACKGROUND: An energy-discriminating capability of a photon counting detector (PCD) can provide many clinical advantages, but several factors, such as charge sharing (CS) and pulse pileup (PP), degrade the capability by distorting the measured x-ray spectrum. To fully exploit the merits of PCDs, it is important to characterize the output of PCDs. Previously proposed PCD output models showed decent agreement with physical PCDs; however, there were still scopes to be improved: a global model-data mismatch and pixel-to-pixel variations. PURPOSES: In this study, we improve a PCD model by using count-rate-dependent model parameters to address the issues and evaluate agreement against physical PCDs. METHODS: The proposed model is based on the cascaded model, and we made model parameters condition-dependent and pixel-specific to deal with the global model-data mismatch and the pixel-to-pixel variation. The parameters are determined by a procedure for model parameter estimation with data acquired from different thicknesses of water or aluminum at different x-ray tube currents. To analyze the effects of having proposed model parameters, we compared three setups of our model: a model with default parameters, a model with global parameters, and a model with global-and-local parameters. For experimental validation, we used CdZnTe-based PCDs, and assessed the performance of the models by calculating the mean absolute percentage errors (MAPEs) between the model outputs and the actual measurements from low count-rates to high count-rates, which have deadtime losses of up to 24%. RESULTS: The outputs of the proposed model visually matched well with the PCD measurements for all test data. For the test data, the MAPEs averaged over all the bins were 49.2-51.1% for a model with default parameters, 8.0-9.8% for a model with the global parameters, and 1.2-2.7% for a model with the global-and-local parameters. CONCLUSION: The proposed model can estimate the outputs of physical PCDs with high accuracy from low to high count-rates. We expect that our model will be actively utilized in applications where the pixel-by-pixel accuracy of a PCD model is important.

Three-dimensional analysis of ductular reactions and their correlation with liver regeneration and fibrosis.

The liver has multiple regeneration modes, including hepatocellular hypertrophy and self-renewal of hepatocytes. When hepatocyte proliferation is impaired, hepatic progenitor cells may proliferate through ductular reaction (DR), differentiate into hepatocytes, and contribute to fibrosis. However, the three-dimensional spatial relationship between DR and regenerating hepatocytes and dynamic changes in DR associated with fibrosis remain poorly understood. Here, we performed three-dimensional (3D) imaging of cleared 42 liver explants with chronic and acute liver diseases and 4 normal livers to visualize DR. In chronic hepatic liver diseases, such as viral hepatitis, steatohepatitis, autoimmune hepatitis, and cryptogenic cirrhosis, the total length and number of branches of DR showed a significant positive correlation. We studied the spatial relationship between DR and GS-expressing cells using glutamine synthetase (GS) and cytokeratin 19 (CK19) as markers of liver regeneration and DR, respectively. The percentage of CK19-positive cells that co-expressed GS was less than 10% in chronic liver diseases. In contrast, nearly one-third of CK19-positive cells co-expressed GS in acute liver diseases, and chronic cholestatic liver diseases, including primary biliary cholangitis and primary sclerosing cholangitis, showed no co-expression. We also found that DR was longer and had more branching in livers with progressive fibrosis compared to those with regressive fibrosis. Our results suggest that DR displays varying degrees of spatial complexity and contribution to liver regeneration. DR may serve as hepatobiliary junctions that maintain continuity between hepatocytes and bile ducts rather than hepatocyte regeneration in chronic liver diseases.

Geometric Mismatch Promotes Anatomic Repair in Periorbital Bony Defects in Skeletally Mature Yucatan Minipigs.

Porous tissue-engineered 3D-printed scaffolds are a compelling alternative to autografts for the treatment of large periorbital bone defects. Matching the defect-specific geometry has long been considered an optimal strategy to restore pre-injury anatomy. However, studies in large animal models have revealed that biomaterial-induced bone formation largely occurs around the scaffold periphery. Such ectopic bone formation in the periorbital region can affect vision and cause disfigurement. To enhance anatomic reconstruction, geometric mismatches are introduced in the scaffolds used to treat full thickness zygomatic defects created bilaterally in adult Yucatan minipigs. 3D-printed, anatomically-mirrored scaffolds are used in combination with autologous stromal vascular fraction of cells (SVF) for treatment. An advanced image-registration workflow is developed to quantify the post-surgical geometric mismatch and correlate it with the spatial pattern of the regenerating bone. Osteoconductive bone growth on the dorsal and ventral aspect of the defect enhances scaffold integration with the native bone while medio-lateral bone growth leads to failure of the scaffolds to integrate. A strong positive correlation is found between geometric mismatch and orthotopic bone deposition at the defect site. The data suggest that strategic mismatch >20% could improve bone scaffold design to promote enhanced regeneration, osseointegration, and long-term scaffold survivability.

Musculoskeletal CT Imaging: State-of-the-Art Advancements and Future Directions.

CT is one of the most widely used modalities for musculoskeletal imaging. Recent advancements in the field include the introduction of four-dimensional CT, which captures a CT image during motion; cone-beam CT, which uses flat-panel detectors to capture the lower extremities in weight-bearing mode; and dual-energy CT, which operates at two different x-ray potentials to improve the contrast resolution to facilitate the assessment of tissue material compositions such as tophaceous gout deposits and bone marrow edema. Most recently, photon-counting CT (PCCT) has been introduced. PCCT is a technique that uses photon-counting detectors to produce an image with higher spatial and contrast resolution than conventional multidetector CT systems. In addition, postprocessing techniques such as three-dimensional printing and cinematic rendering have used CT data to improve the generation of both physical and digital anatomic models. Last, advancements in the application of artificial intelligence to CT imaging have enabled the automatic evaluation of musculoskeletal pathologies. In this review, the authors discuss the current state of the above CT technologies, their respective advantages and disadvantages, and their projected future directions for various musculoskeletal applications.

Subchondral bone in knee osteoarthritis: bystander or treatment target?

The subchondral bone is an important structural component of the knee joint relevant for osteoarthritis (OA) incidence and progression once disease is established. Experimental studies have demonstrated that subchondral bone changes are not simply the result of altered biomechanics, i.e., pathologic loading. In fact, subchondral bone alterations have an impact on joint homeostasis leading to articular cartilage loss already early in the disease process. This narrative review aims to summarize the available and emerging imaging techniques used to evaluate knee OA-related subchondral bone changes and their potential role in clinical trials of disease-modifying OA drugs (DMOADs). Radiographic fractal signature analysis has been used to quantify OA-associated changes in subchondral texture and integrity. Cross-sectional modalities such as cone-beam computed tomography (CT), contrast-enhanced cone beam CT, and micro-CT can also provide high-resolution imaging of the subchondral trabecular morphometry. Magnetic resonance imaging (MRI) has been the most commonly used advanced imaging modality to evaluate OA-related subchondral bone changes such as bone marrow lesions and altered trabecular bone texture. Dual-energy X-ray absorptiometry can provide insight into OA-related changes in periarticular subchondral bone mineral density. Positron emission tomography, using physiological biomarkers of subchondral bone regeneration, has provided additional insight into OA pathogenesis. Finally, artificial intelligence algorithms have been developed to automate some of the above subchondral bone measurements. This paper will particularly focus on semiquantitative methods for assessing bone marrow lesions and their utility in identifying subjects at risk of symptomatic and structural OA progression, and evaluating treatment responses in DMOAD clinical trials.

Combining physics-based models with deep learning image synthesis and uncertainty in intraoperative cone-beam CT of the brain.

BACKGROUND: Image-guided neurosurgery requires high localization and registration accuracy to enable effective treatment and avoid complications. However, accurate neuronavigation based on preoperative magnetic resonance (MR) or computed tomography (CT) images is challenged by brain deformation occurring during the surgical intervention. PURPOSE: To facilitate intraoperative visualization of brain tissues and deformable registration with preoperative images, a 3D deep learning (DL) reconstruction framework (termed DL-Recon) was proposed for improved intraoperative cone-beam CT (CBCT) image quality. METHODS: The DL-Recon framework combines physics-based models with deep learning CT synthesis and leverages uncertainty information to promote robustness to unseen features. A 3D generative adversarial network (GAN) with a conditional loss function modulated by aleatoric uncertainty was developed for CBCT-to-CT synthesis. Epistemic uncertainty of the synthesis model was estimated via Monte Carlo (MC) dropout. Using spatially varying weights derived from epistemic uncertainty, the DL-Recon image combines the synthetic CT with an artifact-corrected filtered back-projection (FBP) reconstruction. In regions of high epistemic uncertainty, DL-Recon includes greater contribution from the FBP image. Twenty paired real CT and simulated CBCT images of the head were used for network training and validation, and experiments evaluated the performance of DL-Recon on CBCT images containing simulated and real brain lesions not present in the training data. Performance among learning- and physics-based methods was quantified in terms of structural similarity (SSIM) of the resulting image to diagnostic CT and Dice similarity metric (DSC) in lesion segmentation compared to ground truth. A pilot study was conducted involving seven subjects with CBCT images acquired during neurosurgery to assess the feasibility of DL-Recon in clinical data. RESULTS: CBCT images reconstructed via FBP with physics-based corrections exhibited the usual challenges to soft-tissue contrast resolution due to image non-uniformity, noise, and residual artifacts. GAN synthesis improved image uniformity and soft-tissue visibility but was subject to error in the shape and contrast of simulated lesions that were unseen in training. Incorporation of aleatoric uncertainty in synthesis loss improved estimation of epistemic uncertainty, with variable brain structures and unseen lesions exhibiting higher epistemic uncertainty. The DL-Recon approach mitigated synthesis errors while maintaining improvement in image quality, yielding 15%-22% increase in SSIM (image appearance compared to diagnostic CT) and up to 25% increase in DSC in lesion segmentation compared to FBP. Clear gains in visual image quality were also observed in real brain lesions and in clinical CBCT images. CONCLUSIONS: DL-Recon leveraged uncertainty estimation to combine the strengths of DL and physics-based reconstruction and demonstrated substantial improvements in the accuracy and quality of intraoperative CBCT. The improved soft-tissue contrast resolution could facilitate visualization of brain structures and support deformable registration with preoperative images, further extending the utility of intraoperative CBCT in image-guided neurosurgery.

Computed Tomography: State-of-the-Art Advancements in Musculoskeletal Imaging.

ABSTRACT: Although musculoskeletal magnetic resonance imaging (MRI) plays a dominant role in characterizing abnormalities, novel computed tomography (CT) techniques have found an emerging niche in several scenarios such as trauma, gout, and the characterization of pathologic biomechanical states during motion and weight-bearing. Recent developments and advancements in the field of musculoskeletal CT include 4-dimensional, cone-beam (CB), and dual-energy (DE) CT. Four-dimensional CT has the potential to quantify biomechanical derangements of peripheral joints in different joint positions to diagnose and characterize patellofemoral instability, scapholunate ligamentous injuries, and syndesmotic injuries. Cone-beam CT provides an opportunity to image peripheral joints during weight-bearing, augmenting the diagnosis and characterization of disease processes. Emerging CBCT technologies improved spatial resolution for osseous microstructures in the quantitative analysis of osteoarthritis-related subchondral bone changes, trauma, and fracture healing. Dual-energy CT-based material decomposition visualizes and quantifies monosodium urate crystals in gout, bone marrow edema in traumatic and nontraumatic fractures, and neoplastic disease. Recently, DE techniques have been applied to CBCT, contributing to increased image quality in contrast-enhanced arthrography, bone densitometry, and bone marrow imaging. This review describes 4-dimensional CT, CBCT, and DECT advances, current logistical limitations, and prospects for each technique.

Quantitative Dual-Energy Imaging of Bone Marrow Edema Using Multisource Cone-Beam CT with Model-Based Decomposition.

Purpose: We investigated the feasibility of dual-energy (DE) detection of bone marrow edema (BME) using a dedicated extremity cone-beam CT (CBCT) with a unique three-source x-ray unit. The sources can be operated at different energies to enable single-scan DE acquisitions. However, they are arranged parallel to the axis of rotation, resulting in incomplete sampling and precluding the application of DE projection-domain decompositions (PDD) for beam-hardening reduction. Therefore, we propose a novel combination of a model-based "one-step" DE two-material decomposition followed by a constrained image-domain change-of-basis to obtain virtual non-calcium (VNCa) images for BME detection. Methods: DE projections were obtained using an "alternating-kV" protocol by operating the peripheral two sources of the CBCT system at low-energy (60 kV, 0.105 mAs/frame) and the central source at high-energy (100 kV, 0.028 mAs/frame), for a total of 600 frames over 216° of gantry rotation. Projections were processed with detector lag, glare and fast Monte Carlo (MC)-based iterative scatter corrections. Model-based material decomposition (MBMD) was then implemented to obtain aluminum (Al) and polyethylene (PE) volume fraction images with minimal beam-hardening. Statistical ray weights in MBMD were modified to account for regions with highly oblique sampling by the peripheral sources. To generate the VNCa maps, image-domain decomposition (IDD) constrained by the volume conservation principle (VCP) was performed to convert the Al and PE MBMD images into volume fractions of water, fat and cortical bone. Accuracy of BME detection was evaluated using physical phantom data acquired on the multi-source extremity CBCT scanner. Results: The proposed framework estimated the volume of BME with ~10% error. The MC-based scatter corrections and the modified MBMD ray weights were essential to achieve such performance - the error without MC scatter corrections was >30%, whereas the uniformity of estimated VNCa images was 3x improved using the modified weights compared to the conventional weights. Conclusions: The proposed DE decomposition framework was able to overcome challenges of high scatter and incomplete sampling to achieve BME detection on a CBCT system with axially-distributed x-ray sources.

Quantification of manipulation forces needed for robot-assisted reduction of the ankle syndesmosis: an initial cadaveric study.

PURPOSE: Manual surgical manipulation of the tibia and fibula is necessary to properly align and reduce the space in ankle fractures involving sprain of the distal tibiofibular syndesmosis. However, manual reduction is highly variable and can result in malreduction in about half of the cases. Therefore, we are developing an image-guided robotic assistant to improve reduction accuracy. The purpose of this study is to quantify the forces associated with reduction of the ankle syndesmosis to define the requirements for our robot design. METHODS: Using a cadaveric specimen, we designed a fixture jig to fix the tibia securely on the operating table. We also designed a custom fibula grasping plate to which a force-torque measuring device is attached. The surgeon manually reduced the fibula utilizing this construct while translational and rotational forces along with displacement were being measured. This was first performed on an intact ankle without ligament injury and after the syndesmosis ligaments were cut. RESULTS: Six manipulation techniques were performed on the three principal directions of reduction at the cadaveric ankle. The results demonstrated the maximum force applied to the lateral direction to be 96.0 N with maximum displacement of 8.5 mm, applied to the anterior-posterior direction to be 71.6 N with maximum displacement of 10.7 mm, and the maximum torque applied to external-internal rotation to be 2.5 Nm with maximum rotation of 24.6°. CONCLUSIONS: The specific forces needed to perform the distal tibiofibular syndesmosis manipulation are not well understood. This study quantified these manipulation forces needed along with their displacement for accurate reduction of ankle syndesmosis. This is a necessary first step to help us define the design requirements of our robotic assistance from the aspects of forces and displacements.

Model-based three-material decomposition in dual-energy CT using the volume conservation constraint.

Objective. We develop a model-based optimization algorithm for 'one-step' dual-energy (DE) CT decomposition of three materials directly from projection measurements.Approach.Since the three-material problem is inherently undetermined, we incorporate the volume conservation principle (VCP) as a pair of equality and nonnegativity constraints into the objective function of the recently reported model-based material decomposition (MBMD). An optimization algorithm (constrained MBMD, CMBMD) is derived that utilizes voxel-wise separability to partition the volume into a VCP-constrained region solved using interior-point iterations, and an unconstrained region (air surrounding the object, where VCP is violated) solved with conventional two-material MBMD. Constrained MBMD (CMBMD) is validated in simulations and experiments in application to bone composition measurements in the presence of metal hardware using DE cone-beam CT (CBCT). A kV-switching protocol with non-coinciding low- and high-energy (LE and HE) projections was assumed. CMBMD with decomposed base materials of cortical bone, fat, and metal (titanium, Ti) is compared to MBMD with (i) fat-bone and (ii) fat-Ti bases.Main results.Three-material CMBMD exhibits a substantial reduction in metal artifacts relative to the two-material MBMD implementations. The accuracies of cortical bone volume fraction estimates are markedly improved using CMBMD, with ∼5-10× lower normalized root mean squared error in simulations with anthropomorphic knee phantoms (depending on the complexity of the metal component) and ∼2-2.5× lower in an experimental test-bench study.Significance.In conclusion, we demonstrated one-step three-material decomposition of DE CT using volume conservation as an optimization constraint. The proposed method might be applicable to DE applications such as bone marrow edema imaging (fat-bone-water decomposition) or multi-contrast imaging, especially on CT/CBCT systems that do not provide coinciding LE and HE ray paths required for conventional projection-domain DE decomposition.

Technical assessment of 2D and 3D imaging performance of an IGZO-based flat-panel X-ray detector.

BACKGROUND: Indirect detection flat-panel detectors (FPDs) consisting of hydrogenated amorphous silicon (a-Si:H) thin-film transistors (TFTs) are a prevalent technology for digital x-ray imaging. However, their performance is challenged in applications requiring low exposure levels, high spatial resolution, and high frame rate. Emerging FPD designs using metal oxide TFTs may offer potential performance improvements compared to FPDs based on a-Si:H TFTs. PURPOSE: This work investigates the imaging performance of a new indium gallium zinc oxide (IGZO) TFT-based detector in 2D fluoroscopy and 3D cone-beam CT (CBCT). METHODS: The new FPD consists of a sensor array combining IGZO TFTs with a-Si:H photodiodes and a 0.7-mm thick CsI:Tl scintillator. The FPD was implemented on an x-ray imaging bench with system geometry emulating intraoperative CBCT. A conventional FPD with a-Si:H TFTs and a 0.6-mm thick CsI:Tl scintillator was similarly implemented as a basis of comparison. 2D imaging performance was characterized in terms of electronic noise, sensitivity, linearity, lag, spatial resolution (modulation transfer function, MTF), image noise (noise-power spectrum, NPS), and detective quantum efficiency (DQE) with entrance air kerma (EAK) ranging from 0.3 to 1.2 µGy. 3D imaging performance was evaluated in terms of the 3D MTF and noise-equivalent quanta (NEQ), soft-tissue contrast-to-noise ratio (CNR), and image quality evident in anthropomorphic phantoms for a range of anatomical sites and dose, with weighted air kerma, K w ${K_w}$ , ranging from 0.8 to 4.9 mGy. RESULTS: The 2D imaging performance of the IGZO-based FPD exhibited up to ∼1.7× lower electronic noise than the a-Si:H FPD at matched pixel pitch. Furthermore, the IGZO FPD exhibited ∼27% increase in mid-frequency DQE (1 mm-1 ) at matched pixel size and dose (EAK ≈ 1.0 µGy) and ∼11% increase after adjusting for differences in scintillator thickness. 2D spatial resolution was limited by the scintillator for each FPD. The IGZO-based FPD demonstrated improved 3D NEQ at all spatial frequencies in both head (≥25% increase for all dose levels) and body (≥10% increase for K w ${K_w}$ ≤2 mGy) imaging scenarios. These characteristics translated to improved low-contrast visualization in anthropomorphic phantoms, demonstrating ≥10% improvement in CNR and extension of the low-dose range for which the detector is input-quantum limited. CONCLUSION: The IGZO-based FPD demonstrated improvements in electronic noise, image lag, and NEQ that translated to measurable improvements in 2D and 3D imaging performance compared to a conventional FPD based on a-Si:H TFTs. The improvements are most beneficial for 2D or 3D imaging scenarios involving low-dose and/or high-frame rate.

A robust, autonomous, volumetric quality assurance method for 3D printed porous scaffolds.

Bone tissue engineering strategies aimed at treating critical-sized craniofacial defects often utilize novel biomaterials and scaffolding. Rapid manufacturing of defect-matching geometries using 3D-printing strategies is a promising strategy to treat craniofacial bone loss to improve aesthetic and regenerative outcomes. To validate manufacturing quality, a robust, three-dimensional quality assurance pipeline is needed to provide an objective, quantitative metric of print quality if porous scaffolds are to be translated from laboratory to clinical settings. Previously published methods of assessing scaffold print quality utilized one- and two-dimensional measurements (e.g., strut widths, pore widths, and pore area) or, in some cases, the print quality of a single phantom is assumed to be representative of the quality of all subsequent prints. More robust volume correlation between anatomic shapes has been accomplished; however, it requires manual user correction in challenging cases such as porous objects like bone scaffolds. Here, we designed porous, anatomically-shaped scaffolds with homogenous or heterogenous porous structures. We 3D-printed the designs with acrylonitrile butadiene styrene (ABS) and used cone-beam computed tomography (CBCT) to obtain 3D image reconstructions. We applied the iterative closest point algorithm to superimpose the computational scaffold designs with the CBCT images to obtain a 3D volumetric overlap. In order to avoid false convergences while using an autonomous workflow for volumetric correlation, we developed an independent iterative closest point (I-ICP10) algorithm using MATLAB®, which applied ten initial conditions for the spatial orientation of the CBCT images relative to the original design. Following successful correlation, scaffold quality can be quantified and visualized on a sub-voxel scale for any part of the volume.

Point-of-care treatment of geometrically complex midfacial critical-sized bone defects with 3D-Printed scaffolds and autologous stromal vascular fraction.

Critical-sized midfacial bone defects present a unique clinical challenge due to their complex three-dimensional shapes and intimate associations with sensory organs. To address this challenge, a point-of-care treatment strategy for functional, long-term regeneration of 2 cm full-thickness segmental defects in the zygomatic arches of Yucatan minipigs is evaluated. A digital workflow is used to 3D-print anatomically precise, porous, biodegradable scaffolds from clinical-grade poly-ε-caprolactone and decellularized bone composites. The autologous stromal vascular fraction of cells (SVF) is isolated from adipose tissue extracts and infused into the scaffolds that are implanted into the zygomatic ostectomies. Bone regeneration is assessed up to 52 weeks post-operatively in acellular (AC) and SVF groups (BV/DV = 0.64 ± 0.10 and 0.65 ± 0.10 respectively). In both treated groups, bone grows from the adjacent tissues and restores the native anatomy. Significantly higher torque is required to fracture the bone-scaffold interface in the SVF (7.11 ± 2.31 N m) compared to AC groups (2.83 ± 0.23 N m). Three-dimensional microcomputed tomography analysis reveals two distinct regenerative patterns: osteoconduction along the periphery of scaffolds to form dense lamellar bone and small islands of woven bone deposits growing along the struts in the scaffold interior. Overall, this study validates the efficacy of using 3D-printed bioactive scaffolds with autologous SVF to restore geometrically complex midfacial bone defects of clinically relevant sizes while also highlighting remaining challenges to be addressed prior to clinical translation.

Sensitivity of the stress field of the proximal femur predicted by CT-based FE analysis to modeling uncertainties.

Proximal femur anatomy and bone mineral density vary widely among individuals, precluding the use of one predefined finite element (FE) model to determine the stress field for all proximal femurs. This variability poses a challenge in current prosthetic hip design approach. Given the numerous options for generating computed tomography (CT)-based FE models, selecting the best methods for defining the mechanical behavior of the proximal femur is difficult. In this study, a combination of computational and experimental approaches was used to explore the susceptibility of the predicted stress field of the proximal femur to different combinations of density-elasticity relationships, element type, element size, and calibration error. Our results suggest that FE models with first-order voxelized elements generated by the Keyak and Falkinstein density-elasticity relationship or quadratic tetrahedral elements generated by the Morgan density-elasticity relationship lead to accurate estimations of the mechanical behavior of human femurs. Other combinations of element size, element type, and mathematical relationships produce less accurate results, especially in the cortical bone of the femoral neck and calcar region. The voxelized model was more susceptible to variation of element size and density-elasticity relationships than FE models with quadratic tetrahedral elements. Regardless of element type, the stress fields predicted by the Keyak and Falkinstein and the Morgan relationships were the most robust to calibration error when deriving material density from CT-generated Hounsfield data. These results provide insight into the implementation of a robust platform for designing patient-specific implants capable of maintaining or modifying the stress in bones.