Assessment of spinal cord injury using ultrasound elastography in a rabbit model in vivo.

The effect of the mechanical micro-environment on spinal cord injury (SCI) and treatment effectiveness remains unclear. Currently, there are limited imaging methods that can directly assess the localized mechanical behavior of spinal cords in vivo. In this study, we apply new ultrasound elastography (USE) techniques to assess SCI in vivo at the site of the injury and at the time of one week post injury, in a rabbit animal model. Eleven rabbits underwent laminectomy procedures. Among them, spinal cords of five rabbits were injured during the procedure. The other six rabbits were used as control. Two neurological statuses were achieved: non-paralysis and paralysis. Ultrasound data were collected one week post-surgery and processed to compute strain ratios. Histologic analysis, mechanical testing, magnetic resonance imaging (MRI), computerized tomography and MRI diffusion tensor imaging (DTI) were performed to validate USE results. Strain ratios computed via USE were found to be significantly different in paralyzed versus non-paralyzed rabbits. The myelomalacia histologic score and spinal cord Young's modulus evaluated in selected animals were in good qualitative agreement with USE assessment. It is feasible to use USE to assess changes in the spinal cord of the presented animal model. In the future, with more experimental data available, USE may provide new quantitative tools for improving SCI diagnosis and prognosis.

A CNN-based method to reconstruct 3-D spine surfaces from US images in vivo.

Three-dimensional (3-D) reconstruction of the spine surface is of strong clinical relevance for the diagnosis and prognosis of spine disorders and intra-operative image guidance. In this paper, we report a new technique to reconstruct lumbar spine surfaces in 3-D from non-invasive ultrasound (US) images acquired in free-hand mode. US images randomly sampled from in vivo scans of 9 rabbits were used to train a U-net convolutional neural network (CNN). More specifically, a late fusion (LF)-based U-net trained jointly on B-mode and shadow-enhanced B-mode images was generated by fusing two individual U-nets and expanding the set of trainable parameters to around twice the capacity of a basic U-net. This U-net was then applied to predict spine surface labels in in vivo images obtained from another rabbit, which were then used for 3-D spine surface reconstruction. The underlying pose of the transducer during the scan was estimated by registering stacks of US images to a geometrical model derived from corresponding CT data and used to align detected surface points. Final performance of the reconstruction method was assessed by computing the mean absolute error (MAE) between pairs of spine surface points detected from US and CT and by counting the total number of surface points detected from US. Comparison was made between the LF-based U-net and a previously developed phase symmetry (PS)-based method. Using the LF-based U-net, the averaged number of US surface points across the lumbar region increased by 21.61% and MAE reduced by 26.28% relative to the PS-based method. The overall MAE (in mm) was 0.24±0.29. Based on these results, we conclude that: 1) the proposed U-net can detect the spine posterior arch with low MAE and large number of US surface points and 2) the newly proposed reconstruction framework may complement and, under certain circumstances, be used without the aid of an external tracking system in intra-operative spine applications.

Identification of ultrasound imaging markers to quantify long bone regeneration in a segmental tibial defect sheep model in vivo.

The healing of large bone defects has been investigated for decades due to its complexity and clinical relevance. Ultrasound (US) methods have shown promise in monitoring bone healing, but no quantitative method to assess regenerated bone morphology in US images has been presented yet. In this study, we investigate new US morphometric parameters to quantify bone regeneration in vivo. A segmental tibial defect was surgically created and stabilized in a sheep animal model. US and computed tomography (CT) imaging data were collected two months post-surgery. New bone was assessed, reconstructed and quantified from the US and CT data using 3 morphometric parameters: the new-bone bulk (NBB), new-bone surface (NBS) and new-bone contact (NBC). The distance (mm) between surface reconstructions from repeated US was [Formula: see text] and from US and CT was [Formula: see text]. In the mid-shaft of the defected tibia, US measurements of NBB, NBS and NBC were significantly higher than the corresponding CT measurements ([Formula: see text]). Based on our results, we conclude that US may complement CT to reconstruct and quantify bone regrowth, especially in its early stages.

Modeling and Analysis of Ultrasound Elastographic Axial Strains for Spine Fracture Identification.

This study reports the first use of ultrasound (US) elastography for imaging spinal fractures by assessing the mechanical response of the soft tissue at the posterior vertebra boundary to a uniaxial compression in rabbit ex vivo samples. Three-dimensional finite-element (FE) models of the vertebra-soft tissue complex in rabbit samples are generated and analyzed to evaluate the distribution of the axial normal and shear strains at the vertebra-soft tissue interface. Experiments on the same samples are performed to corroborate simulation findings. Results of this study indicate that the distribution of the axial strains manifests as distinct patterns around intact and fractured vertebrae. Numerical characteristics of the axial strain's spatial distribution are further used to construct two shape descriptors to make inferences on spinal abnormalities: 1) axial normal strain asymmetry for assessing the presence of fractures and 2) principal orientation of axial shear strain concentration regions (shear zones) for measurement of spinous process dislocation. This study demonstrates that axial normal strain and axial shear strain maps obtained via US elastography can provide a new means to detect spine fractures and abnormalities in the selected ex vivo animal models. Spinal fracture detection is important for the assessment of spinal cord injuries and stability. However, identification of spinal fractures using US is currently challenging. Our results show that features resulting from strain elastograms can serve as a useful adjunct to B-mode images in identifying spine fractures in the selected animal samples, and this information could be helpful in clinical settings.

Assessment of the long bone inter-fragmentary gap size in ultrasound strain elastograms.

The inter-fragmentary gap size (IFGS) is a critical factor affecting the propensity of bone healing. In this paper, we present a study to analyze ultrasound strain elastographic numerical features in samples with distinct IFGS using both simulations and experiments. Six fractured rabbit hind leg samples in total were used in this study with controlled IFGS of 1 mm, 5 mm and 1 cm. For the simulation, computed tomography (CT) scans of all six samples were used to create solid models. Finite element analysis (FEA) and subsequent elastography simulations were performed on the 3D models to produce tensorial strain field data. Features of bony fragment separation were defined on different strain components and computed for strains segmented at varying thresholds to evaluate their performance in estimating the IFGS. A threshold for each strain component was then determined, based on which extra 3D features of interest were defined and extracted from the segmented strain data. Then, all 3D features were compared statistically among the three nominal groups. Additional simulations and experiments of axial shear strain elastography (ASSE) on the median coronal plane of the same samples were also performed. Our results indicate that coronal plane axial shear (CPAS) strain elastography produces a separation feature which is statistically correlated with the IFGS, and that our elastography simulation module is effective in predicting the CPAS elastographic strain behavior for different IFGS.

A Model-Based Approach to Investigate the Effect of a Long Bone Fracture on Ultrasound Strain Elastography.

The mechanical behavior of long bones and fractures has been under investigation for many decades due to its complexity and clinical relevance. In this paper, we report a new subject-specific methodology to predict and analyze the mechanical behavior of the soft tissue at a bone interface with the intent of identifying the presence and location of bone abnormalities with high accuracy, spatial resolution, and contrast. The proposed methodology was tested on both intact and fractured rabbit femur samples with finite element-based 3-D simulations, created from actual femur computed tomography data, and ultrasound elastography experiments. The results included in this study demonstrate that elastographic strains at the bone/soft tissue interface can be used to differentiate fractured femurs from the intact ones on a distribution level. These results also demonstrate that coronal plane axial shear strain creates a unique contrast mechanism that can be used to reliably detect fractures (both complete and incomplete) in long bones. Kruskal-Wallis test further demonstrates that the contrast measure for the fracture group (simulation: 2.1286±0.2206; experiment: 2.7034 ± 1.0672) is significantly different from that for the intact group (simulation: 0 ± 0; experiment: 1.1540±0.6909) when using coronal plane axial shear strain elastography ( < 0.01). We conclude that: 1) elastography techniques can be used to accurately identify the presence and location of fractures in a long bone and 2) the proposed model-based approach can be used to predict and analyze strains at a bone fracture site and to better interpret experimental elastographic data.

Spine surface detection from local phase-symmetry enhanced ridges in ultrasound images.

PURPOSE: Improper administration of epidural anesthesia can result in nerve complications. This problem is exacerbated for obese patients whose vertebrae cannot be palpated. Ultrasound (US) has recently emerged as an attractive imaging modality for accurate epidural placement. However, anesthesiologists untrained in US have difficulty interpreting the anatomy in noisy spinal US images. Furthermore, the complex geometry in spinal US images is characterized by a discontinuous intensity profile because the transducer is often not perpendicularly oriented to spine surface regions such as laminae, articular and transverse processes. This makes the interpretation of spinal images more challenging than typical long bone surface images. In this article, we propose a new method to segment the spine anatomy in US images obtained in both the transverse and paramedian planes. METHODS: A set of 108 B-mode images were randomly chosen from 35 cine loops obtained from scanning the lumbar and thoracic vertebrae of 17 healthy volunteers with a BMI ranging from 19.5 to 27.9. A local phase-symmetry technique was applied to the B-mode images for enhancement of bone-like ridges, and the spine blobs were subsequently classified. The segmented spine surface from the blobs was compared against a radiologist's manual delineation of the spine surface. RESULTS: For the performance of the spine blob classifier, we obtain a Matthews Correlation Coefficient (MCC) of 0.77 and a geometric mean (G-mean) of 0.96. The mean absolute error between the manual delineation of the laminae by the radiologist and the automatic laminae segmentation is found to be 0.26 mm with a maximum possible absolute error of 2.01 mm for spinal US images of 70 mm depth. CONCLUSIONS: Our proposed technique successfully performs automatic segmentation of the spine surface - specifically the laminae, ligamentum flava, spinous, transverse, and articular processes - and can be extended to any bone anatomy present in an US image. This has implications for 3D visualization of bone surfaces and, without loss of generality, the vertebral column.

Characterization of ventral incisional hernia and repair using shear wave elastography.

BACKGROUND: To assess the integrity of hernia repair, imaging modalities such as computed tomography or ultrasound (US) are commonly used. Neither modality has currently the capacity to simultaneously image the mesh and quantify a prosthetic and surrounding tissue stiffness. In this pilot study, we hypothesize that US shear wave elastography (SWE) can be used to identify a polyester mesh and a biologic graft and to assess their stiffness noninvasively in a rat model of bridging hernia repair. METHODS: Lewis rats underwent hernia creation and repair with Parietex or Strattice at 30 d. After 3 mo, the animals were euthanized, and the Young's Modulus was measured using SWE. Three-dimensional reconstructions of the hernia pre- and post-repair were performed using in-house image processing algorithms. RESULTS: SWE was capable of accurate and real-time assessment and diagnosis of the hernia defects in vivo. Young's Modulus of Parietex meshes and Strattice grafts as estimated from the shear wave elastograms were found to be statistically different from each other (P < 0.05). Accurate three-dimensional reconstructions of the hernia defects pre- and post-repair were generated. CONCLUSIONS: In this study, we demonstrate the feasibility of using US SWE to detect ventral hernias and evaluate mesh repair in vivo. Our results indicate that the presence of a hernia and repair can be reliably visualized by SWE and three dimensionally reconstructed. Thus, this technique may provide both structural and functional information regarding the hernia and the repair.

Ultrasound elastography assessment of bone/soft tissue interface.

We report on the use of elastographic imaging techniques to assess the bone/soft tissue interface, a region that has not been previously investigated but may provide important information about fracture and bone healing. The performance of axial strain elastograms and axial shear strain elastograms at the bone/soft tissue interface was studied ex vivo on intact and fractured canine and ovine tibias. Selected ex vivo results were corroborated on intact sheep tibias in vivo. The elastography results were statistically analyzed using elastographic image quality tools. The results of this study demonstrate distinct patterns in the distribution of the normalized local axial strains and axial shear strains at the bone/soft tissue interface with respect to the background soft tissue. They also show that the relative strength and distribution of the elastographic parameters change in the presence of a fracture and depend on the degree of misalignment between the fracture fragments. Thus, elastographic imaging modalities might be used in the future to obtain information regarding the integrity of bones and to assess the severity of fractures, alignment of bone fragments as well as to follow bone healing.