Biomechanical 3-D finite element modeling of the human breast using MRI data

Construction of multi-component finite element model to predict biomechanical behaviour of breasts during running and quantification of the stiffness impact of internal structure

This study aims to investigate the biomechanical behaviour and the stiffness impact of the breast internal components during running. To achieve this, a novel nonlinear multi-component dynamic finite element method (FEM) has been established, which uses experimental data obtained via 4D scanning technology and a motion capture system. The data are used to construct a geometric model that comprises the rigid body, layers of soft tissues, skin, pectoralis major muscle, fat, ligaments and glandular tissues. The traditional point-to-point method has a relative mean absolute error of less than 7.92% while the latest surface-to-surface method has an average Euclidean distance (d) of 7.05 mm, validating the simulated results. After simulating the motion of the different components of the breasts, the displacement analysis confirms that when the motion reaches the moment of largest displacement, the displacement of the breast components is proportional to their distance from the chest wall. A biomechanical analysis indicates that the stress sustained by the breast components in ascending order is the glandular tissues, pectoralis major muscle, adipose tissues, and ligaments. The ligaments provide the primary support during motion, followed by the pectoralis major muscle. In addition, specific stress points of the breast components are identified. The stiffness impact experiment indicates that compared with ligaments, the change of glandular tissue stiffness had a slightly more obvious effect on the breast surface. The findings serve as a valuable reference for the medical field and sports bra industry to enhance breast protection during motion.

Breast simulation pipeline: From medical imaging to patient-specific simulations

BACKGROUND

Breast-conserving surgery is the most acceptable operation for breast cancer removal from an invasive and psychological point of view. Before the surgical procedure, a preoperative MRI is performed in the prone configuration, while the surgery is achieved in the supine position. This leads to a considerable movement of the breast, including the tumor, between the two poses, complicating the surgeon's task.

METHODS

In this work, a simulation pipeline allowing the computation of patient-specific geometry and the prediction of personalized breast material properties was put forward. Through image segmentation, a finite element model including the subject-specific geometry is established. By first computing an undeformed state of the breast, the geometrico-material model is calibrated by surface acquisition in the intra-operative stance.

FINDINGS

Using an elastic corotational formulation, the patient-specific mechanical properties of the breast and skin were identified to obtain the best estimates of the supine configuration. The final results are a shape-fitting closest point residual of 4.00 mm for the mechanical parameters Ebreast=0.32 kPa and Eskin=22.72 kPa, congruent with the current state-of-the-art. The Covariance Matrix Adaptation Evolution Strategy optimizer converges on average between 5 to 30 min depending on the initial parameters, reaching a simulation speed of 20 s. To our knowledge, our model offers one of the best compromises between accuracy and speed.

INTERPRETATION

Satisfactory results were obtained for the estimation of breast deformation from preoperative to intra-operative configuration. Furthermore, we have demonstrated the clinical feasibility of such applications using a simulation framework that aims at the smallest disturbance of the actual surgical pipeline.

Estimation of the biomechanical mammographic deformation of the breast using machine learning models

BACKGROUND

A typical problem in the registration of MRI and X-ray mammography is the nonlinear deformation applied to the breast during mammography. We have developed a method for virtual deformation of the breast using a biomechanical model automatically constructed from MRI. The virtual deformation is applied in two steps: unloaded state estimation and compression simulation. The finite element method is used to solve the deformation process. However, the extensive computational cost prevents its usage in clinical routine.

METHODS

We propose three machine learning models to overcome this problem: an extremely randomized tree (first model), extreme gradient boosting (second model), and deep learning-based bidirectional long short-term memory with an attention layer (third model) to predict the deformation of a biomechanical model. We evaluated our methods with 516 breasts with realistic compression ratios up to 76%.

FINDINGS

We first applied one-fold validation, in which the second and third models performed better than the first model. We then applied ten-fold validation. For the unloaded state estimation, the median RMSE for the second and third models is 0.8 mm and 1.2 mm, respectively. For the compression, the median RMSE is 3.4 mm for both models. We evaluated correlations between model accuracy and characteristics of the clinical datasets such as compression ratio, breast volume, and tissue types.

INTERPRETATION

Using the proposed models, we achieved accurate results comparable to the finite element model, with a speedup of factor 240 using the extreme gradient boosting model. These proposed models can replace the finite element model simulation, enabling clinically relevant real-time application.

A review of bioengineering techniques applied to breast tissue: Mechanical properties, tissue engineering and finite element analysis

Female breast cancer was the most prevalent cancer worldwide in 2020, according to the Global Cancer Observatory. As a prophylactic measure or as a treatment, mastectomy and lumpectomy are often performed at women. Following these surgeries, women normally do a breast reconstruction to minimize the impact on their physical appearance and, hence, on their mental health, associated with self-image issues. Nowadays, breast reconstruction is based on autologous tissues or implants, which both have disadvantages, such as volume loss over time or capsular contracture, respectively. Tissue engineering and regenerative medicine can bring better solutions and overcome these current limitations. Even though more knowledge needs to be acquired, the combination of biomaterial scaffolds and autologous cells appears to be a promising approach for breast reconstruction. With the growth and improvement of additive manufacturing, three dimensional (3D) printing has been demonstrating a lot of potential to produce complex scaffolds with high resolution. Natural and synthetic materials have been studied in this context and seeded mainly with adipose derived stem cells (ADSCs) since they have a high capability of differentiation. The scaffold must mimic the environment of the extracellular matrix (ECM) of the native tissue, being a structural support for cells to adhere, proliferate and migrate. Hydrogels (e.g., gelatin, alginate, collagen, and fibrin) have been a biomaterial widely studied for this purpose since their matrix resembles the natural ECM of the native tissues. A powerful tool that can be used in parallel with experimental techniques is finite element (FE) modeling, which can aid the measurement of mechanical properties of either breast tissues or scaffolds. FE models may help in the simulation of the whole breast or scaffold under different conditions, predicting what might happen in real life. Therefore, this review gives an overall summary concerning the human breast, specifically its mechanical properties using experimental and FE analysis, and the tissue engineering approaches to regenerate this particular tissue, along with FE models.

Mechanical Model and FEM Simulations for Efforts on Biceps and Triceps Muscles under Vertical Load: Mathematical Formulation of Results

Although isometric contractions in human muscles have been analyzed several times, there are no FEA models that allow us to use the same modeled joint (the elbow under our case) in different conditions. Most elbow joints use 3D elements for meshing. Representing the muscles in the joint is quite useful when the study is focused on the muscle itself, knowing stress distribution on muscle, and checking damage in muscle in a detailed manner (tendon–muscle insertion, for example). However, this technique is not useful for studying muscle behavior at different positions of the joint. This study, based on the mechanical model of the elbow joint, proposes a methodology for modelling muscles that will be studied in different positions by meshing them with 1D elements. Furthermore, the methodology allows us to calculate biceps and triceps efforts under load for different angles of elbow joint aperture. The simulation results have been mathematically modelled to obtain general formulations for these efforts, depending on the load and the aperture angle.

Wound-Healing on Minimally Invasive Systems with Micro-Arc Oxidized (MAO) and Acid-Etched (SLA) Surface Modifications in Finite Element Liver Model

The purpose of this research is to study micro-arc oxidation (MAO) and acid etching (SLA) in the liver after minimally invasive treatment. The three-dimensional (3D) liver model is rebuilt by magnetic resonance imaging (MRI) to simulate the clinical process. The processing layer is an important factor in clinical applications. Several studies have investigated the finite element model of the liver; however, few people have studied the equipment model with treated layers. The results revealed that temperature was significantly decreased when devices using MAO and SLA thin films. In addition, the SLA treatment group showed a relatively low temperature, which was 12.88% lower. It is an effective means to reduce abnormally thermal injury and lateral thermal area. The present study reveals that the novel nanostructured thin film on electrosurgery devices is an effective means of improving the performance of reducing over thermal injury and uniformly distributing temperature in livers.

Experimental characterisation and modelling of breast Cooper's ligaments

The aim of this study was to characterise the mechanical behaviour of Cooper's ligaments. Such ligaments are collagenous breast tissue that create a three-dimensional structure over the entire breast volume. Ten ligaments were extracted from a human cadaver, from which 28 samples were cut and used to perform uniaxial tensile tests. Histological analysis showed that the main direction of the fibres visible to the naked eye corresponds to the orientation of the fibres on a microscopic scale. The specimens were cut according to this orientation, which allowed the sample to be stretched in the main fibre direction. From these experimental stretch/stress curves, an original anisotropic hyperelastic constitutive law is proposed to model the behaviour of Cooper's ligaments and the material parameter validity is discussed.

Improvement of performance and sensitivity of 2D and 3D image reconstruction in EIT using EFG forward model

This paper presents a pure element-free Galerkin method (EFGM) forward model for image reconstruction in 2D and 3D electrical impedance tomography (EIT) using an adaptive current injection method. In EIT systems with the adapting current injection method, both static and dynamic images can be reconstructed; however, determination of electrode contact impedances in the complete electrode model is difficult and the Gap model is used. In this paper, in the EIT forward problem a weak form functional based on the Gap model and a pure EFGM approach are developed, and in the EIT inverse problem, Jacobian matrix is computed by the EFGM, and a fast integration technique is introduced to calculate the entries of the Jacobian matrix within an adequate computation time. The influence of increasing the density of nodes at and near the electrodes with steep electric potential gradients on the accuracy of FEM and EFGM forward solutions is investigated, and the performance of the image reconstruction algorithm with the proposed fast integration technique is examined. The numerical results reveal that the proposed EFGM forward model with the fast integration technique has an efficient performance both in terms of mean relative imaging errors and computational time.

Multimodal Patient-Specific Registration for Breast Imaging Using Biomechanical Modeling with Reference to AI Evaluation of Breast Tumor Change

Background: The strategy to combat the problem associated with large deformations in the breast due to the difference in the medical imaging of patient posture plays a vital role in multimodal medical image registration with artificial intelligence (AI) initiatives. How to build a breast biomechanical model simulating the large-scale deformation of soft tissue remains a challenge but is highly desirable. Methods: This study proposed a hybrid individual-specific registration model of the breast combining finite element analysis, property optimization, and affine transformation to register breast images. During the registration process, the mechanical properties of the breast tissues were individually assigned using an optimization process, which allowed the model to become patient specific. Evaluation and results: The proposed method has been extensively tested on two datasets collected from two independent institutions, one from America and another from Hong Kong. Conclusions: Our method can accurately predict the deformation of breasts from the supine to prone position for both the Hong Kong and American samples, with a small target registration error of lesions.

An anatomically detailed and personalizable head injury model: Significance of brain and white matter tract morphological variability on strain

Finite element head (FE) models are important numerical tools to study head injuries and develop protection systems. The generation of anatomically accurate and subject-specific head models with conforming hexahedral meshes remains a significant challenge. The focus of this study is to present two developmental works: first, an anatomically detailed FE head model with conforming hexahedral meshes that has smooth interfaces between the brain and the cerebrospinal fluid, embedded with white matter (WM) fiber tracts; second, a morphing approach for subject-specific head model generation via a new hierarchical image registration pipeline integrating Demons and Dramms deformable registration algorithms. The performance of the head model is evaluated by comparing model predictions with experimental data of brain-skull relative motion, brain strain, and intracranial pressure. To demonstrate the applicability of the head model and the pipeline, six subject-specific head models of largely varying intracranial volume and shape are generated, incorporated with subject-specific WM fiber tracts. DICE similarity coefficients for cranial, brain mask, local brain regions, and lateral ventricles are calculated to evaluate personalization accuracy, demonstrating the efficiency of the pipeline in generating detailed subject-specific head models achieving satisfactory element quality without further mesh repairing. The six head models are then subjected to the same concussive loading to study the sensitivity of brain strain to inter-subject variability of the brain and WM fiber morphology. The simulation results show significant differences in maximum principal strain and axonal strain in local brain regions (one-way ANOVA test, p < 0.001), as well as their locations also vary among the subjects, demonstrating the need to further investigate the significance of subject-specific models. The techniques developed in this study may contribute to better evaluation of individual brain injury and the development of individualized head protection systems in the future. This study also contains general aspects the research community may find useful: on the use of experimental brain strain close to or at injury level for head model validation; the hierarchical image registration pipeline can be used to morph other head models, such as smoothed-voxel models.

Breast Biomechanics: What Do We Really Know?

Although half the world's population will develop breasts, there is limited research documenting breast structure or motion. Understanding breast structure and motion, however, is imperative for numerous applications, such as breast reconstruction, breast modeling to better diagnose and treat breast pathologies, and designing effective sports bras. To be impactful, future breast biomechanics research needs to fill gaps in our knowledge, particularly related to breast composition and density, and to improve methods to accurately measure the complexities of three-dimensional breast motion. These methods should then be used to investigate breast biomechanics while individuals, who represent the full spectrum of women in the population, participate in a broad range of activities of daily living and recreation.

Breast Tumor Diagnosis Using Finite-Element Modeling Based on Clinical in vivo Elastographic Data

OBJECTIVES

This study exploited finite-element modeling (FEM) to simulate breast tissue multicompression during ultrasound elastography to classify breast tumors based on their nonlinear biomechanical properties.

METHODS

Numeric simulations were first calculated by using 3-dimensional (3D) virtual models with an assumed tumor's geometric dimensions but with actual material properties to test and validate the FEM. Further numeric simulations were used to construct 3D models based on in vivo experimental data to verify our models. The models were designed for each individual in vivo case, emphasizing the geometry, position, and biomechanical properties of the breast tissue. At different compression levels, tissue strains were analyzed between the tumors and the background normal tissues to explore their nonlinearity and classify the tumor type. Tumor classification parameters were deduced by using a power-law relationship between the applied compressive forces and strain differences.

RESULTS

Classification parameters were compared between benign and malignant tumors, for which they were found to be statistically significant in classifying the tumor types (P < .05) by both the validation and verification of FEM. We compared the classification parameters between the in vivo and FEM classifications, for which they were found to be strongly correlated (R = 0.875; P < .001), with no statistical differences between their outcomes (P = .909).

CONCLUSIONS

Good agreement between the model outcomes and the in vivo diagnostics was reported. The implemented models were validated and verified. The introduced 3D modeling method may augment elastographic methods to preliminary classify breast tumors at an early stage.

An anatomically detailed and personalizable head injury model: Significance of brain and white matter tract morphological variability on strain

Finite element head (FE) models are important numerical tools to study head injuries and develop protection systems. The generation of anatomically accurate and subject-specific head models with conforming hexahedral meshes remains a significant challenge. The focus of this study is to present two developmental works: first, an anatomically detailed FE head model with conforming hexahedral meshes that has smooth interfaces between the brain and the cerebrospinal fluid, embedded with white matter (WM) fiber tracts; second, a morphing approach for subject-specific head model generation via a new hierarchical image registration pipeline integrating Demons and Dramms deformable registration algorithms. The performance of the head model is evaluated by comparing model predictions with experimental data of brain-skull relative motion, brain strain, and intracranial pressure. To demonstrate the applicability of the head model and the pipeline, six subject-specific head models of largely varying intracranial volume and shape are generated, incorporated with subject-specific WM fiber tracts. DICE similarity coefficients for cranial, brain mask, local brain regions, and lateral ventricles are calculated to evaluate personalization accuracy, demonstrating the efficiency of the pipeline in generating detailed subject-specific head models achieving satisfactory element quality without further mesh repairing. The six head models are then subjected to the same concussive loading to study the sensitivity of brain strain to inter-subject variability of the brain and WM fiber morphology. The simulation results show significant differences in maximum principal strain and axonal strain in local brain regions (one-way ANOVA test, p < 0.001), as well as their locations also vary among the subjects, demonstrating the need to further investigate the significance of subject-specific models. The techniques developed in this study may contribute to better evaluation of individual brain injury and the development of individualized head protection systems in the future. This study also contains general aspects the research community may find useful: on the use of experimental brain strain close to or at injury level for head model validation; the hierarchical image registration pipeline can be used to morph other head models, such as smoothed-voxel models.

Deformable Mapping Method to Relate Lesions in Dedicated Breast CT Images to Those in Automated Breast Ultrasound and Digital Breast Tomosynthesis Images

This work demonstrates the potential for using a deformable mapping method to register lesions between dedicated breast computed tomography (bCT) and both automated breast ultrasound (ABUS) and digital breast tomosynthesis (DBT) images (craniocaudal [CC] and mediolateral oblique [MLO] views). Two multi-modality breast phantoms with external fiducial markers attached were imaged by the three modalities. The DBT MLO view was excluded for the second phantom. The automated deformable mapping algorithm uses biomechanical modeling to determine corresponding lesions based on distances between their centers of mass (d_COM) in the deformed bCT model and the reference model (DBT or ABUS). For bCT to ABUS, the mean d_COM was 5.2 ± 2.6 mm. For bCT to DBT (CC), the mean d_COM was 5.1 ± 2.4 mm. For bCT to DBT (MLO), the mean d_COM was 4.7 ± 2.5 mm. This application could help improve a radiologist's efficiency and accuracy in breast lesion characterization, using multiple imaging modalities.

Mechanical characterisation of an organic phantom candidate for breast tissue

The aim of this study is to analyse the mechanical properties of agarose, which can be used as a mechanical phantom for breast tissues. In general, tissue mimicking materials may be used to achieve a better understanding of the structure and properties of tissues and organs; this work shares these objectives. The phantom can be used as a tissue surrogate with realistic mechanical behaviour for biomechanical applications. To validate agarose as a suitable mechanical phantom for breast tissues, indentation tests were performed in homogeneous, rectangular agarose blocks. Blocks with different agarose concentrations were moulded and tested. An estimation of the material stiffness was then compared with experimental data on different breast tissues’ types found in literature. The phantom stiffness increased consistently with agarose concentration and stress. The results show that agarose-based mechanical phantoms of stiffer tissues require higher agarose concentrations (0.5% and 0.6%). In contrast, normal tissues can be mimicked with 0.3% and 0.4% of agarose. In addition, it was observed that preconditioning affects the mechanical properties of the gel, as it had already been shown in literature for breast tissues.

In vivo measurement of surface skin strain during human gait to improve the design of rehabilitation devices

When designing any rehabilitation, sportswear or exoskeleton device the mechanical behaviour of the body segment must be known, specifically the skin, because an excessive tissue strain may lead to ulceration and bedsores. To date, it is not known if the kinematic variability between subjects have an effect on the skin strain field, and therefore, in the design and manufacturing of rehabilitation products, such as orthoses. Several studies have analysed the skin deformation during human motion, nevertheless, the comparison between the skin strain field in different subjects during normal or pathological gait has not been reported yet. This work presents a comparison of skin strain analysis for different gait patterns to study the differences between people and, specifically, if it is possible to standardize the orthotic design between subjects with the same gait disorder. Moreover, the areas with relatively minimum strain during the ankle-foot motion are compared to improve the design of structural parts of rehabilitation devices. In this case, a validated 3D digital image correlation system has been used for this purpose combined with strain ellipse theory. The results demonstrate variations in the skin strain field between subjects with the same pathology and similarities between subjects with normal gait. However, more studies and experiments are necessaries to validate this hypothesis and also to test it between different gait pathologies.

Deformable mapping technique to correlate lesions in digital breast tomosynthesis and automated breast ultrasound images

PURPOSE

To develop a deformable mapping technique to match corresponding lesions between digital breast tomosynthesis (DBT) and automated breast ultrasound (ABUS) images.

METHODS

External fiducial markers were attached to the surface of two CIRS multi-modality compressible breast phantoms (A and B) containing multiple simulated lesions. Both phantoms were imaged with DBT (upright positioning with cranial-caudal compression) and ABUS (supine positioning with anterior-to-chest wall compression). The lesions and markers were manually segmented by three different readers. Reader segmentation similarity and reader reproducibility were assessed using Dice similarity coefficients (DSC) and distances between centers of mass (d_COM ). For deformable mapping between the modalities each reader's segmented dataset was processed with an automated deformable mapping algorithm as follows: First, Morfeus, a finite element (FE) based multi-organ deformable image registration platform, converted segmentations into triangular surface meshes. Second, Altair HyperMesh, a FE pre-processor, created base FE models for the ABUS and DBT data sets. All deformation is performed on the DBT image data; the ABUS image sets remain fixed throughout the process. Deformation was performed on the external skin contour (DBT image set) to match the external skin contour on the ABUS set, and the locations of the external markers were used to morph the skin contours to be within a user-defined distance. Third, the base DBT-FE model was deformed with the FE analysis solver, Optistruct. Deformed DBT lesions were correlated with matching lesions in the base ABUS FE model. Performance (lesion correlation) was assessed with d_COM for all corresponding lesions and lesion overlap. Analysis was performed to determine the minimum number of external fiducial markers needed to create the desired correlation and the improvement of correlation with the use of external markers.

RESULTS

Average DSC for reader similarity ranged from 0.88 to 0.91 (ABUS) and 0.57 to 0.83 (DBT). Corresponding d_COM ranged from 0.20 to 0.36 mm (ABUS) and 0.11 to 1.16 mm (DBT). Lesion correlation is maximized when all corresponding markers are within a maximum distance of 5 mm. For deformable mapping of phantom A, without the use of external markers, only two of six correlated lesions showed overlap with an average lesion d_COM of 6.8 ± 2.8 mm. With use of three external fiducial markers, five of six lesions overlapped and average d_COM improved to 4.9 ± 2.4 mm. For deformable mapping of Phantom B without external markers analysis, four lesions were correlated of seven with overlap between only one of seven lesions, and an average lesion d_COM of 9.7 ± 3.5 mm. With three external markers, all seven possible lesions were correlated with overlap between four of seven lesions. The average d_COM was 8.5 ± 4.0 mm.

CONCLUSION

This work demonstrates the potential for a deformable mapping technique to relate corresponding lesions in DBT and ABUS images by showing improved lesion correspondence and reduced lesion registration errors with the use of external fiducial markers. The technique should improve radiologists' characterization of breast lesions which can reduce patient callbacks, misdiagnoses and unnecessary biopsies.

Analytical derivation of elasticity in breast phantoms for deformation tracking

PURPOSE

Patient-specific biomedical modeling of the breast is of interest for medical applications such as image registration, image guided procedures and the alignment for biopsy or surgery purposes. The computation of elastic properties is essential to simulate deformations in a realistic way. This study presents an innovative analytical method to compute the elastic modulus and evaluate the elasticity of a breast using magnetic resonance (MRI) images of breast phantoms.

METHODS

An analytical method for elasticity computation was developed and subsequently validated on a series of geometric shapes, and on four physical breast phantoms that are supported by a planar frame. This method can compute the elasticity of a shape directly from a set of MRI scans. For comparison, elasticity values were also computed numerically using two different simulation software packages.

RESULTS

Application of the different methods on the geometric shapes shows that the analytically derived elongation differs from simulated elongation by less than 9% for cylindrical shapes, and up to 18% for other shapes that are also substantially vertically supported by a planar base. For the four physical breast phantoms, the analytically derived elasticity differs from numeric elasticity by 18% on average, which is in accordance with the difference in elongation estimation for the geometric shapes. The analytic method has shown to be multiple orders of magnitude faster than the numerical methods.

CONCLUSION

It can be concluded that the analytical elasticity computation method has good potential to supplement or replace numerical elasticity simulations in gravity-induced deformations, for shapes that are substantially supported by a planar base perpendicular to the gravitational field. The error is manageable, while the calculation procedure takes less than one second as opposed to multiple minutes with numerical methods. The results will be used in the MRI and Ultrasound Robotic Assisted Biopsy (MURAB) project.

Generation of polychromatic projection for dedicated breast computed tomography simulation using anthropomorphic numerical phantom

Numerical simulations are fundamental to the development of medical imaging systems because they can save time and effort in research and development. In this study, we developed a method of creating the virtual projection images that are necessary to study dedicated breast computed tomography (BCT) systems. Anthropomorphic software breast phantoms of the conventional compression type were synthesized and redesigned to meet the requirements of dedicated BCT systems. The internal structure of the breast was randomly constructed to develop the proposed phantom, enabling the internal structure of a naturally distributed real breast to be simulated. When using the existing monochromatic photon incidence assumption for projection-image generation, it is not possible to simulate various artifacts caused by the X-ray spectrum, such as the beam hardening effect. Consequently, the system performance could be overestimated. Therefore, we considered the polychromatic spectrum in the projection image generation process and verified the results. The proposed method is expected to be useful for the development and optimization of BCT systems.

Towards an in-plane methodology to track breast lesions using mammograms and patient-specific finite-element simulations

In breast cancer screening or diagnosis, it is usual to combine different images in order to locate a lesion as accurately as possible. These images are generated using a single or several imaging techniques. As x-ray-based mammography is widely used, a breast lesion is located in the same plane of the image (mammogram), but tracking it across mammograms corresponding to different views is a challenging task for medical physicians. Accordingly, simulation tools and methodologies that use patient-specific numerical models can facilitate the task of fusing information from different images. Additionally, these tools need to be as straightforward as possible to facilitate their translation to the clinical area. This paper presents a patient-specific, finite-element-based and semi-automated simulation methodology to track breast lesions across mammograms. A realistic three-dimensional computer model of a patient's breast was generated from magnetic resonance imaging to simulate mammographic compressions in cranio-caudal (CC, head-to-toe) and medio-lateral oblique (MLO, shoulder-to-opposite hip) directions. For each compression being simulated, a virtual mammogram was obtained and posteriorly superimposed to the corresponding real mammogram, by sharing the nipple as a common feature. Two-dimensional rigid-body transformations were applied, and the error distance measured between the centroids of the tumors previously located on each image was 3.84 mm and 2.41 mm for CC and MLO compression, respectively. Considering that the scope of this work is to conceive a methodology translatable to clinical practice, the results indicate that it could be helpful in supporting the tracking of breast lesions.

Automated modification and fusion of voxel models to construct body phantoms with heterogeneous breast tissue: Application to MRI simulations

BACKGROUND

Human voxel models incorporating detailed anatomical features are vital tools for the computational evaluation of electromagnetic (EM) fields within the body. Besides whole-body human voxel models, phantoms representing smaller heterogeneous anatomical features are often employed; for example, localized breast voxel models incorporating fatty and fibroglandular tissues have been developed for a variety of EM applications including mammography simulation and dosimetry, magnetic resonance imaging (MRI), and ultra-wideband microwave imaging. However, considering wavelength effects, electromagnetic modeling of the breast at sub-microwave frequencies necessitates detailed breast phantoms in conjunction with whole-body voxel models.

METHODS

Heterogeneous breast phantoms are sized to fit within radiofrequency coil hardware, modified by voxel-wise extrusion, and fused to whole-body models using voxel-wise, tissue-dependent logical operators. To illustrate the utility of this method, finite-difference time-domain simulations are performed using a whole-body model integrated with a variety of available breast phantoms spanning the standard four tissue density classifications representing the majority of the population.

RESULTS

The software library uses a combination of voxel operations to seamlessly size, modify, and fuse eleven breast phantoms to whole-body voxel models. The software is publicly available on GitHub and is linked to the file exchange at MATLAB^® Central. Simulations confirm the proportions of fatty and fibroglandular tissues in breast phantoms have significant yet predictable implications on projected power deposition in tissue.

CONCLUSIONS

Breast phantoms may be modified and fused to whole-body voxel models using the software presented in this work; user considerations for the open-source software and resultant phantoms are discussed. Furthermore, results indicate simulating breast models as predominantly fatty tissue can considerably underestimate the potential for tissue heating in women with substantial fibroglandular tissue.

Patient-specific Deformation Modelling via Elastography: Application to Image-guided Prostate Interventions

Image-guided prostate interventions often require the registration of preoperative magnetic resonance (MR) images to real-time transrectal ultrasound (TRUS) images to provide high-quality guidance. One of the main challenges for registering MR images to TRUS images is how to estimate the TRUS-probe-induced prostate deformation that occurs during TRUS imaging. The combined statistical and biomechanical modeling approach shows promise for the adequate estimation of prostate deformation. However, the right setting of the biomechanical parameters is very crucial for realistic deformation modeling. We propose a patient-specific deformation model equipped with personalized biomechanical parameters obtained from shear wave elastography to reliably predict the prostate deformation during image-guided interventions. Using data acquired from a prostate phantom and twelve patients with suspected prostate cancer, we compared the prostate deformation model with and without patient-specific biomechanical parameters in terms of deformation estimation accuracy. The results show that the patient-specific deformation model possesses favorable model ability, and outperforms the model without patient-specific biomechanical parameters. The employment of the patient-specific biomechanical parameters obtained from elastography for deformation modeling shows promise for providing more precise deformation estimation in applications that use computer-assisted image-guided intervention systems.

Building a virtual simulation platform for quasistatic breast ultrasound elastography using open source software: A preliminary investigation

PURPOSE

Quasistatic ultrasound elastography (QUE) is being used to augment in vivo characterization of breast lesions. Results from early clinical trials indicated that there was a lack of confidence in image interpretation. Such confidence can only be gained through rigorous imaging tests using complex, heterogeneous but known media. The objective of this study is to build a virtual breast QUE simulation platform in the public domain that can be used not only for innovative QUE research but also for rigorous imaging tests.

METHODS

The main thrust of this work is to streamline biomedical ultrasound simulations by leveraging existing open source software packages including Field II (ultrasound simulator), VTK (geometrical visualization and processing), FEBio [finite element (FE) analysis], and Tetgen (mesh generator). However, integration of these open source packages is nontrivial and requires interdisciplinary knowledge. In the first step, a virtual breast model containing complex anatomical geometries was created through a novel combination of image-based landmark structures and randomly distributed (small) structures. Image-based landmark structures were based on data from the NIH Visible Human Project. Subsequently, an unstructured FE-mesh was created by Tetgen. In the second step, randomly positioned point scatterers were placed within the meshed breast model through an octree-based algorithm to make a virtual breast ultrasound phantom. In the third step, an ultrasound simulator (Field II) was used to interrogate the virtual breast phantom to obtain simulated ultrasound echo data. Of note, tissue deformation generated using a FE-simulator (FEBio) was the basis of deforming the original virtual breast phantom in order to obtain the postdeformation breast phantom for subsequent ultrasound simulations. Using the procedures described above, a full cycle of QUE simulations involving complex and highly heterogeneous virtual breast phantoms can be accomplished for the first time.

RESULTS

Representative examples were used to demonstrate capabilities of this virtual simulation platform. In the first set of three ultrasound simulation examples, three heterogeneous volumes of interest were selected from a virtual breast ultrasound phantom to perform sophisticated ultrasound simulations. These resultant B-mode images realistically represented the underlying complex but known media. In the second set of three QUE examples, advanced applications in QUE were simulated. The first QUE example was to show breast tumors with complex shapes and/or compositions. The resultant strain images showed complex patterns that were normally seen in freehand clinical ultrasound data. The second and third QUE examples demonstrated (deformation-dependent) nonlinear strain imaging and time-dependent strain imaging, respectively.

CONCLUSIONS

The proposed virtual QUE platform was implemented and successfully tested in this study. Through show-case examples, the proposed work has demonstrated its capabilities of creating sophisticated QUE data in a way that cannot be done through the manufacture of physical tissue-mimicking phantoms and other software. This open software architecture will soon be made available in the public domain and can be readily adapted to meet specific needs of different research groups to drive innovations in QUE.

A review of biomechanically informed breast image registration

Breast radiology encompasses the full range of imaging modalities from routine imaging via x-ray mammography, magnetic resonance imaging and ultrasound (both two- and three-dimensional), to more recent technologies such as digital breast tomosynthesis, and dedicated breast imaging systems for positron emission mammography and ultrasound tomography. In addition new and experimental modalities, such as Photoacoustics, Near Infrared Spectroscopy and Electrical Impedance Tomography etc, are emerging. The breast is a highly deformable structure however, and this greatly complicates visual comparison of imaging modalities for the purposes of breast screening, cancer diagnosis (including image guided biopsy), tumour staging, treatment monitoring, surgical planning and simulation of the effects of surgery and wound healing etc. Due primarily to the challenges posed by these gross, non-rigid deformations, development of automated methods which enable registration, and hence fusion, of information within and across breast imaging modalities, and between the images and the physical space of the breast during interventions, remains an active research field which has yet to translate suitable methods into clinical practice. This review describes current research in the field of breast biomechanical modelling and identifies relevant publications where the resulting models have been incorporated into breast image registration and simulation algorithms. Despite these developments there remain a number of issues that limit clinical application of biomechanical modelling. These include the accuracy of constitutive modelling, implementation of representative boundary conditions, failure to meet clinically acceptable levels of computational cost, challenges associated with automating patient-specific model generation (i.e. robust image segmentation and mesh generation) and the complexity of applying biomechanical modelling methods in routine clinical practice.

The Three-Dimensional Techniques in the Objective Measurement of Breast Aesthetics

Breast symmetry, size, and shape are key components of aesthetic outcomes of augmentation mammoplasty, reduction, and reconstruction. Many have claimed that the 3D scanning technique, which measures breast volumes directly and assesses the asymmetry of the chest and breast on a 3D model, is superior to anthropometric measuring in accuracy, precision, and reproducibility. The documented methods of 3D body surface imaging include laser scanning, stereo photography and so on. To achieve ideal aesthetic results, individualized surgery planning based on a reliable virtual model of the prospective surgery outcome could be of considerable value in decision making and assisting in guidance for the surgery procedure. Additionally, the 3D scanning technique is applicable in postoperative monitoring of morphological change, notably, in a dynamic way. Another distinguishing feature is that it enables virtual division of breast volume, thus surgeons could virtually divide the breast volumes into portions using 3D scanning during the programming and evaluation of surgery plans. However, because 3D surface scanning cannot look through the breast substances and reach the interspace between the chest and posterior border of the breast/dorsal limit of the breast, the inframammary fold in larger breasts cannot be correctly imaged, leaving the preoperative inframammary fold reference lacking. Therefore, 3D scanning is thought to be inaccurate in large and/or ptotic breasts. Another fact that prevents 3D scanning from wide application is its high cost and lack of access.

LEVEL OF EVIDENCE IV

This journal requires that authors assign a level of evidence to each article. For a full description of these Evidence-Based Medicine ratings, please refer to the Table of Contents or the online Instructions to Authors www.springer.com/00266.

Symmetric Biomechanically Guided Prone-to-Supine Breast Image Registration

Prone-to-supine breast image registration has potential application in the fields of surgical and radiotherapy planning, image guided interventions, and multi-modal cancer diagnosis, staging, and therapy response prediction. However, breast image registration of three dimensional images acquired in different patient positions is a challenging problem, due to large deformations induced to the soft breast tissue caused by the change in gravity loading. We present a symmetric, biomechanical simulation based registration framework which aligns the images in a central, virtually unloaded configuration. The breast tissue is modelled as a neo-Hookean material and gravity is considered as the main source of deformation in the original images. In addition to gravity, our framework successively applies image derived forces directly into the unloading simulation in place of a subsequent image registration step. This results in a biomechanically constrained deformation. Using a finite difference scheme avoids an explicit meshing step and enables simulations to be performed directly in the image space. The explicit time integration scheme allows the motion at the interface between chest and breast to be constrained along the chest wall. The feasibility and accuracy of the approach presented here was assessed by measuring the target registration error (TRE) using a numerical phantom with known ground truth deformations, nine clinical prone MRI and supine CT image pairs, one clinical prone-supine CT image pair and four prone-supine MRI image pairs. The registration reduced the mean TRE for the numerical phantom experiment from initially 19.3 to 0.9 mm and the combined mean TRE for all fourteen clinical data sets from 69.7 to 5.6 mm.

A class of kernel based real-time elastography algorithms

In this paper, a novel real-time kernel-based and gradient-based Phase Root Seeking (PRS) algorithm for ultrasound elastography is proposed. The signal-to-noise ratio of the strain image resulting from this method is improved by minimizing the cross-correlation discrepancy between the pre- and post-compression radio frequency signals with an adaptive temporal stretching method and employing built-in smoothing through an exponentially weighted neighborhood kernel in the displacement calculation. Unlike conventional PRS algorithms, displacement due to tissue compression is estimated from the root of the weighted average of the zero-lag cross-correlation phases of the pair of corresponding analytic pre- and post-compression windows in the neighborhood kernel. In addition to the proposed one, the other time- and frequency-domain elastography algorithms (Ara et al., 2013; Hussain et al., 2012; Hasan et al., 2012) proposed by our group are also implemented in real-time using Java where the computations are serially executed or parallely executed in multiple processors with efficient memory management. Simulation results using finite element modeling simulation phantom show that the proposed method significantly improves the strain image quality in terms of elastographic signal-to-noise ratio (SNRe), elastographic contrast-to-noise ratio (CNRe) and mean structural similarity (MSSIM) for strains as high as 4% as compared to other reported techniques in the literature. Strain images obtained for the experimental phantom as well as in vivo breast data of malignant or benign masses also show the efficacy of our proposed method over the other reported techniques in the literature.

Interrogating a multifactorial model of breast conserving therapy with clinical data

Most women with early stage breast cancer do not require removal of the entire breast to treat their cancer; instead, up to 70% of women can be effectively and safely treated by breast conserving therapy (BCT) with surgical removal of the tumor only (lumpectomy) followed by radiation treatment of the remaining breast tissue. Unfortunately, the final contour and cosmesis of the treated breast is suboptimal in approximately 30% of patients. The ability to accurately predict breast contour after BCT for breast cancer could significantly improve patient decision-making regarding the choice of surgery for breast cancer. Our overall hypothesis is that the complex interplay among mechanical forces due to gravity, breast tissue constitutive law distribution, inflammation induced by radiotherapy and internal stress generated by the healing process play a dominant role in determining the success or failure of lumpectomy in preserving the breast contour and cosmesis. We have shown here from a first patient study that even in the idealistic situation of excellent cosmetic outcome this problem requires multiscale modeling. We propose a method to decide which component of the model works best for each phase of healing and what parameters should be considered dominant and patient specific. This patient study is part of a clinical trial registered on ClinicalTrial.gov, identifier NCT02310711.

Multi-axis dose accumulation of noninvasive image-guided breast brachytherapy through biomechanical modeling of tissue deformation using the finite element method

PURPOSE

Noninvasive image-guided breast brachytherapy delivers conformal HDR (192)Ir brachytherapy treatments with the breast compressed, and treated in the cranial-caudal and medial-lateral directions. This technique subjects breast tissue to extreme deformations not observed for other disease sites. Given that, commercially-available software for deformable image registration cannot accurately co-register image sets obtained in these two states, a finite element analysis based on a biomechanical model was developed to deform dose distributions for each compression circumstance for dose summation.

MATERIAL AND METHODS

The model assumed the breast was under planar stress with values of 30 kPa for Young's modulus and 0.3 for Poisson's ratio. Dose distributions from round and skin-dose optimized applicators in cranial-caudal and medial-lateral compressions were deformed using 0.1 cm planar resolution. Dose distributions, skin doses, and dose-volume histograms were generated. Results were examined as a function of breast thickness, applicator size, target size, and offset distance from the center.

RESULTS

Over the range of examined thicknesses, target size increased several millimeters as compression thickness decreased. This trend increased with increasing offset distances. Applicator size minimally affected target coverage, until applicator size was less than the compressed target size. In all cases, with an applicator larger or equal to the compressed target size, > 90% of the target covered by > 90% of the prescription dose. In all cases, dose coverage became less uniform as offset distance increased and average dose increased. This effect was more pronounced for smaller target-applicator combinations.

CONCLUSIONS

The model exhibited skin dose trends that matched MC-generated benchmarking results within 2% and clinical observations over a similar range of breast thicknesses and target sizes. The model provided quantitative insight on dosimetric treatment variables over a range of clinical circumstances. These findings highlight the need for careful target localization and accurate identification of compression thickness and target offset.

In-vivo quantification of human breast deformation associated with the position change from supine to upright

Stereophotographic imaging and digital image correlation are used to determine the variation of breast skin deformation as the subject orientation is altered from supine to upright. A change in subject's position from supine to upright can result in significant stretches in some parts of the breast skin. The maximum of the major principal stretch ratio of the skin is different in different subjects and varies in the range of 1.25-1.60. It is also found that the boundaries of the breast move significantly relative to the skeletal structure and other fixed points such as the sternal notch. Such measurements are crucial since they provide basic data for validation of biomechanical breast models based on finite element formulations.

Probabilistic estimation of mechanical properties of biomaterials using atomic force microscopy

Nanoindentation using contact-mode atomic force microscopy (AFM) has emerged as a powerful tool for effective material characterization of a wide variety of biomaterials across multiple length scales. However, the interpretation of force-indentation experimental data from AFM is subject to some debate. Uncertainties in AFM data analysis stems from two primary sources: The exact point of contact between the AFM probe and the biological specimen and the variability in the spring constant of the AFM probe. While a lot of attention has been directed toward addressing the contact-point uncertainty, the effect of variability in the probe spring constant has not received sufficient attention. In this paper, we report on an error-in-variables-based Bayesian change-point approach to quantify the elastic modulus of human breast tissue samples after accounting for variability in both contact point and the probe spring constant. We also discuss the efficacy of our approach to a wide range of hyperparameter values using a sensitivity analysis.

Phase-based direct average strain estimation for elastography

In this paper, a phase-based direct average strain estimation method is developed. A mathematical model is presented to calculate axial strain directly from the phase of the zero-lag cross-correlation function between the windowed precompression and stretched post-compression analytic signals. Unlike phase-based conventional strain estimators, for which strain is computed from the displacement field, strain in this paper is computed in one step using the secant algorithm by exploiting the direct phase-strain relationship. To maintain strain continuity, instead of using the instantaneous phase of the interrogative window alone, an average phase function is defined using the phases of the neighboring windows with the assumption that the strain is essentially similar in a close physical proximity to the interrogative window. This method accounts for the effect of lateral shift but without requiring a prior estimate of the applied strain. Moreover, the strain can be computed both in the compression and relaxation phases of the applied pressure. The performance of the proposed strain estimator is analyzed in terms of the quality metrics elastographic signal-to-noise ratio (SNRe), elastographic contrast-to-noise ratio (CNRe), and mean structural similarity (MSSIM), using a finite element modeling simulation phantom. The results reveal that the proposed method performs satisfactorily in terms of all the three indices for up to 2.5% applied strain. Comparative results using simulation and experimental phantom data, and in vivo breast data of benign and malignant masses also demonstrate that the strain image quality of our method is better than the other reported techniques.

Stiffness, compliance, resilience, and creep deformation: understanding implant-soft tissue dynamics in the augmented breast: fundamentals based on materials science

BACKGROUND

Postoperative tissue stretch deformities are among the possible complications in breast augmentation. These deformities are responsible for many potential risks such as bottoming-out deformity, breakdown of the inframammary fold, permanent tissue atrophy, sensory loss, and breast distortion (visible implant edges and traction rippling), among others. Although the elastic properties of the breast are a major concern for plastic surgeons, concepts such as stiffness, compliance, elasticity, and resilience have not been sufficiently defined or explored in the plastic surgery literature.

METHODS

Similar to any other material, living tissues are subject to the fundamentals of the mechanics of materials. Based on their experience with more than 5,000 breast augmentations, the authors explored the basic fundamentals of the mechanics of materials in search of a rational explanation for long-term results in breast augmentation and augmentation-mastopexy.

RESULTS

A basic law of the mechanics of materials determines that when a material (e.g., breast) is loaded with a force (e.g., implant), it produces a stress that causes the material to deform (e.g., breast augmentation), and this behavior might be graphed in a theoretical material's stress-stress curve. This deformation will increase with time although the load (implant) remains constant, a concept termed "creep deformation." Because the breast, like all human tissues, is a viscoelastic material, the application of concepts such as elastic and plastic deformation, stiffness, compliance, resilience, and creep deformation can and should be applied to breast augmentation surgery.

CONCLUSIONS

The authors have found that the principles of the mechanics of materials can provide plastic surgeons with some clues for a predictable, long-lasting good result in breast augmentation and augmentation-mastopexy. Future studies are needed to develop these concepts and evaluate how they might individually determine the mid- and long-term outcomes of augmented breasts.

Deformable medical image registration: a survey

Deformable image registration is a fundamental task in medical image processing. Among its most important applications, one may cite: 1) multi-modality fusion, where information acquired by different imaging devices or protocols is fused to facilitate diagnosis and treatment planning; 2) longitudinal studies, where temporal structural or anatomical changes are investigated; and 3) population modeling and statistical atlases used to study normal anatomical variability. In this paper, we attempt to give an overview of deformable registration methods, putting emphasis on the most recent advances in the domain. Additional emphasis has been given to techniques applied to medical images. In order to study image registration methods in depth, their main components are identified and studied independently. The most recent techniques are presented in a systematic fashion. The contribution of this paper is to provide an extensive account of registration techniques in a systematic manner.

Breast image analysis for risk assessment, detection, diagnosis, and treatment of cancer

The role of breast image analysis in radiologists' interpretation tasks in cancer risk assessment, detection, diagnosis, and treatment continues to expand. Breast image analysis methods include segmentation, feature extraction techniques, classifier design, biomechanical modeling, image registration, motion correction, and rigorous methods of evaluation. We present a review of the current status of these task-based image analysis methods, which are being developed for the various image acquisition modalities of mammography, tomosynthesis, computed tomography, ultrasound, and magnetic resonance imaging. Depending on the task, image-based biomarkers from such quantitative image analysis may include morphological, textural, and kinetic characteristics and may depend on accurate modeling and registration of the breast images. We conclude with a discussion of future directions.

Morbid obesity may increase dislocation in total hip patients: a biomechanical analysis

BACKGROUND

Obesity has reached epidemic proportions in the United States. Recently, obesity, especially morbid obesity, has been linked to increased rates of dislocation after THA. The reasons are unclear. Soft tissue engagement caused by increased thigh girth has been proposed as a possible mechanism for decreased joint stability.

QUESTIONS/PURPOSES

We asked (1) whether thigh soft tissue impingement could decrease THA stability, and if so, at what level of BMI this effect might become evident; and (2) how THA construct factors (eg, head size, neck offset, cup abduction) might affect stability in the morbidly obese.

METHODS

The obesity effect was explored by augmenting a physically validated finite element model of a total hip construct previously comprising just implant hardware and periarticular (capsular) soft tissue. The model augmentation involved using anatomic and anthropometric data to include graded levels of increased thigh girth. Parametric computations were run to assess joint stability for two head sizes (28 and 36 mm), for normal versus high neck offset, and for multiple cup abduction angles.

RESULTS

Thigh soft tissue impingement lowered the resistance to dislocation for BMIs of 40 or greater. Dislocation risk increased monotonically above this threshold as a function of cup abduction angle, independent of hardware impingement events. Increased head diameter did not substantially improve joint stability. High-offset necks decreased the dislocation risk.

CONCLUSIONS

Excessive obesity creates conditions that compromise stability of THAs. Given such conditions, our model suggests reduced cup abduction, high neck offset, and full-cup coverage would reduce the risks of dislocation events.

Data-driven breast decompression and lesion mapping from digital breast tomosynthesis

Digital breast tomosynthesis (DBT) emerges as a new 3D modality for breast cancer screening and diagnosis. Like in conventional 2D mammography the breast is scanned in a compressed state. For orientation during surgical planning, e.g., during presurgical ultrasound-guided anchor-wire marking, as well as for improving communication between radiologists and surgeons it is desirable to estimate an uncompressed model of the acquired breast along with a spatial mapping that allows localizing lesions marked in DBT in the uncompressed model. We therefore propose a method for 3D breast decompression and associated lesion mapping from 3D DBT data. The method is entirely data-driven and employs machine learning methods to predict the shape of the uncompressed breast from a DBT input volume. For this purpose a shape space has been constructed from manually annotated uncompressed breast surfaces and shape parameters are predicted by multiple multi-variate Random Forest regression. By exploiting point correspondences between the compressed and uncompressed breasts, lesions identified in DBT can be mapped to approximately corresponding locations in the uncompressed breast model. To this end, a thin-plate spline mapping is employed. Our method features a novel completely data-driven approach to breast shape prediction that does not necessitate prior knowledge about biomechanical properties and parameters of the breast tissue. Instead, a particular deformation behavior (decompression) is learned from annotated shape pairs, compressed and uncompressed, which are obtained from DBT and magnetic resonance image volumes, respectively. On average, shape prediction takes 26s and achieves a surface distance of 15.80 +/- 4.70 mm. The mean localization error for lesion mapping is 22.48 +/- 8.67 mm.