Marco Palanca

The use of digital image correlation in the biomechanical area: a review

Abstract This paper offers an overview of the potentialities and limitations of digital image correlation (DIC) as a technique for measuring displacements and strain in biomechanical applications. This review is mainly intended for biomechanists who are not yet familiar with DIC. This review includes over 150 papers and covers different dimensional scales, from the microscopic level (tissue level) up to macroscopic one (organ level). As DIC involves a high degree of computation, and of operator-dependent decisions, reliability of displacement and strain measurements by means of DIC cannot be taken for granted. Methodological problems and existing solutions are summarized and compared, whilst open issues are addressed. Topics addressed include: preparation methods for the speckle pattern on different tissues; software settings; systematic and random error associated with DIC measurement. Applications to hard and soft tissues at different dimensional scales are described and analyzed in terms of strengths and limitations. The potentialities and limitations of DIC are highlighted, also in comparison with other experimental techniques (strain gauges, other optical techniques, digital volume correlation) and numerical methods (finite element analysis), where synergies and complementarities are discussed. In order to provide an overview accessible to different scientists working in the field of biomechanics, this paper intentionally does not report details of the algorithms and codes used in the different studies.

Three-Dimensional Local Measurements of Bone Strain and Displacement: Comparison of Three Digital Volume Correlation Approaches

Different digital volume correlation (DVC) approaches are currently available or under development for bone tissue micromechanics. The aim of this study was to compare accuracy and precision errors of three DVC approaches for a particular three-dimensional (3D) zero-strain condition. Trabecular and cortical bone specimens were repeatedly scanned with a micro-computed tomography (CT). The errors affecting computed displacements and strains were extracted for a known virtual translation, as well as for repeated scans. Three DVC strategies were tested: two local approaches, based on fast-Fourier-transform (DaVis-FFT) or direct-correlation (DaVis-DC), and a global approach based on elastic registration and a finite element (FE) solver (ShIRT-FE). Different computation subvolume sizes were tested. Much larger errors were found for the repeated scans than for the virtual translation test. For each algorithm, errors decreased asymptotically for larger subvolume sizes in the range explored. Considering this particular set of images, ShIRT-FE showed an overall better accuracy and precision (a few hundreds microstrain for a subvolume of 50 voxels). When the largest subvolume (50-52 voxels) was applied to cortical bone, the accuracy error obtained for repeated scans with ShIRT-FE was approximately half of that for the best local approach (DaVis-DC). The difference was lower (250 microstrain) in the case of trabecular bone. In terms of precision, the errors shown by DaVis-DC were closer to the ones computed by ShIRT-FE (differences of 131 microstrain and 157 microstrain for cortical and trabecular bone, respectively). The multipass computation available for DaVis software improved the accuracy and precision only for the DaVis-FFT in the virtual translation, particularly for trabecular bone. The better accuracy and precision of ShIRT-FE, followed by DaVis-DC, were obtained with a higher computational cost when compared to DaVis-FFT. The results underline the importance of performing a quantitative comparison of DVC methods on the same set of samples by using also repeated scans, other than virtual translation tests only. ShIRT-FE provides the most accurate and precise results for this set of images. However, both DaVis approaches show reasonable results for large nodal spacing, particularly for trabecular bone. Finally, this study highlights the importance of using sufficiently large subvolumes, in order to achieve better accuracy and precision.

USE OF DIGITAL IMAGE CORRELATION TO INVESTIGATE THE BIOMECHANICS OF THE VERTEBRA

Digital image correlation (DIC) is being introduced to the biomechanical field. However, as DIC relies on a number of major assumptions, it requires a careful optimization in order to obtain accurate and precise results. The first step was the preparation of the speckle pattern by an airbrush spray gun following a factorial design to explore the different settings: the different speckle patterns created were analyzed to achieve the optimal speckle size, with minimal dispersion of speckle sizes. A benchmark test, with an aluminum specimen prepared with the speckle pattern, was conducted in which the errors affecting the computed strain were measured in a zero-displacement, zero-strain condition. The software parameters (facet size, step, and local regression) were singularly analyzed in order to understand their behavior on the final output. Moreover, the hardware parameters (camera gain, exposure, lens distortion) were analyzed. The output showed that a careful optimization allowed the reducing the systematic and random errors, respectively, from 150 to 10 microstrain and from 600 to 110 microstrain. Finally, the acquired know-how was applied to a biological specimen (human vertebra).

Digital volume correlation can be used to estimate local strains in natural and augmented vertebrae: an organ-level study

Digital Volume Correlation (DVC) has become popular for measuring the strain distribution inside bone structures. A number of methodological questions are still open: the reliability of DVC to investigate augmented bone tissue, the variability of the errors between different specimens of the same type, the distribution of measurement errors inside a bone, and the possible presence of preferential directions. To address these issues, five augmented and five natural porcine vertebrae were subjected to repeated zero-strain micro-CT scan (39μm voxel size). The acquired images were processed with two independent DVC approaches (a local and a global one), considering different computation sub-volume sizes, in order to assess the strain measurement uncertainties. The systematic errors generally ranged within ±100 microstrain and did not depend on the computational sub-volume. The random error was higher than 1000 microstrain for the smallest sub-volume and rapidly decreased: with a sub-volume of 48 voxels the random errors were typically within 200 microstrain for both DVC approaches. While these trends were rather consistent within the sample, two individual specimens had unpredictably larger errors. For this reason, a zero-strain check on each specimen should always be performed before any in-situ micro-CT testing campaign. This study clearly shows that, when sufficient care is dedicated to preliminary methodological work, different DVC computation approaches allow measuring the strain with a reduced overall error (approximately 200 microstrain). Therefore, DVC is a viable technique to investigate strain in the elastic regime in natural and augmented bones.

Local displacement and strain uncertainties in different bone types by digital volume correlation of synchrotron microtomograms

Understanding bone mechanics at different hierarchical levels is fundamental to improve preclinical and clinical assessments of bone strength. Digital Volume Correlation (DVC) is the only experimental measurement technique used for measuring local displacements and calculating local strains within bones. To date, its combination with laboratory source micro-computed tomography (LS-microCT) data typically leads to high uncertainties, which limit its application. Here, the benefits of synchrotron radiation micro-computed tomography (SR-microCT) for DVC are reported. Specimens of cortical and trabecular bovine bone and murine tibiae, were each scanned under zero-strain conditions with an effective voxel size of 1.6μm. In order to consider the effect of the voxel size, analyses were also performed on downsampled images with voxel size of 8μm. To evaluate displacement and strain uncertainties, each pair of tomograms was correlated using a global DVC algorithm (ShIRT-FE). Displacement random errors for original SR-microCT ranged from 0.024 to 0.226μm, depending on DVC nodal spacing. Standard deviation of strain errors was below 200 microstrain (ca. 1/10 of the strain associated with physiological loads) for correlations performed with a measurement spatial resolution better than 40μm for cortical bovine bone (240μm for downsampled images), 80μm for trabecular bovine bone (320μm for downsampled images) and murine tibiae (120μm for downsampled images). This study shows that the uncertainties of SR-microCT-based DVC, estimated from repeated scans, are lower than those obtained from LS-microCT-based DVC on similar specimens and low enough to measure accurately the local deformation at the tissue level.

Precision of Digital Volume Correlation Approaches for Strain Analysis in Bone Imaged with Micro-Computed Tomography at Different Dimensional Levels

Accurate measurement of local strain in heterogeneous and anisotropic bone tissue is fundamental to understand the pathophysiology of musculoskeletal diseases, to evaluate the effect of interventions from preclinical studies, and to optimize the design and delivery of biomaterials. Digital volume correlation (DVC) can be used to measure the three-dimensional displacement and strain fields from micro-Computed Tomography (µCT) images of loaded specimens. However, this approach is affected by the quality of the input images, by the morphology and density of the tissue under investigation, by the correlation scheme, and by the operational parameters used in the computation. Therefore, for each application the precision of the method should be evaluated. In this paper we present the results collected from datasets analyzed in previous studies as well as new data from a recent experimental campaign for characterizing the relationship between the precision of two different DVC approaches and the spatial resolution of the outputs. Different bone structures scanned with laboratory source µCT or Synchrotron light µCT (SRµCT) were processed in zero-strain tests to evaluate the precision of the DVC methods as a function of the subvolume size that ranged from 8 to 2500 micrometers. The results confirmed that for every microstructure the precision of DVC improves for larger subvolume size, following power laws. However, for the first time large differences in the precision of both local and global DVC approaches have been highlighted when SRµCT or in vivo µCT images were used instead of conventional ex vivo µCT. These findings suggest that in situ mechanical testing protocols applied in SRµCT facilities should be optimized in order to allow DVC analyses of localized strain measurements. Moreover, for in vivo µCT applications DVC analyses should be performed only with relatively course spatial resolution for achieving a reasonable precision of the method. In conclusion, we have extensively shown that the precision of both tested DVC approaches is affected by different bone structures, different input image resolution and different subvolume sizes. Before each specific application DVC users should always apply a similar approach to find the best compromise between precision and spatial resolution of the measurements.

The Size of Simulated Lytic Metastases Affects the Strain Distribution on the Anterior Surface of the Vertebra

Metastatic lesions of the vertebra are associated with risk of fracture, which can be disabling and life-threatening. In the literature, attempts are found to identify the parameters that reduce the strength of a metastatic vertebra leading to spine instability. However, a number of controversial issues remain. Our aim was to quantify how the strain distribution in the vertebral body is affected by the presence and by the size of a simulated metastatic defect. Five cadaveric thoracic spine segments were subjected to non-destructive presso-flexion while intact, and after simulation of metastases of increasing size. For the largest defect, the specimens were eventually tested to failure. The full-field strain distribution in the elastic range was measured with digital image correlation (DIC) on the anterior surface of the vertebral body. The mean strain in the vertebra remained similar to the intact when the defects were smaller than 30% of the vertebral volume. The mean strains became significantly larger than in the intact for larger defects. The map of strain and its statistical distribution indicated a rather uniform condition in the intact vertebra and with defects smaller than 30%. Conversely, the strain distribution became significantly different from the intact for defects larger than 30%. A strain peak appeared in the region of the simulated metastasis, where fracture initiated during the final destructive test. This is a first step in understanding how the features of metastasis influence the vertebral strain, and for the construction of a mechanistic model to predicted fracture.