Uncertainty and Resolution of Speckle Photography on Micro Samples

The efficient development of new materials with defined properties requires fast methods of testing low volumes of material, such as a high-throughput investigation of spherical metallic micro samples with varying compositions and structuring treatments. A classical material testing method for macro samples, the tensile test cannot be employed for analyzing the mechanical properties of spherical samples with diameters below 1 mm since there are currently no methods for holding and stretching them. A combination between the incremental electrohydraulic extrusion as stress actuation unit and the speckle photography as strain measuring method is presented for obtaining the required mechanical characteristics. Positive longitudinal strain is generated at stepwise extrusion through < 1 mm wide forming channels using a liquid punch and the deformation is observed in situ after each forming step at the interface between the micro sample and a transparent window integrated into the forming die. The occurring local strain fields with a lateral extension down to 100 µm and high gradients require displacement measurements with high lateral resolution over a large range of local dislocations between 0.1 and > 10 µm. It is unknown, whether the speckle strain measuring method is able to provide the necessary low uncertainty for the required resolution in the whole measuring range. Therefore, theoretical and experimental investigations of the deformation measurability using the speckle correlation method are presented, showing that local displacements up to 10 µm can be measured with a spatial resolution between 3 and 10 µm depending on the displacement size. The dominant effect influencing the measurement uncertainty for displacements at this high spatial resolution is the speckle noise, resulting into measurement uncertainties of less than 100 nm. Hence, speckle photography is shown to be applicable for tensile test on micro samples.


Introduction
The diversity of technological development requires a constant search for new materials with defined properties. In order to accelerate ample and expansive studies of various metal compositions and structuring treatments, a highthroughput investigation of spherical micro samples is currently considered [1]. The investigation of micro samples results in descriptors, which are connected with the mechanical properties of the macro materials [2]. However, the classical material testing methods cannot be usually applied for spherical micro samples. Major adjustments of existing methods are necessary to accommodate for the dimensions and the geometry of the micro samples or novel methods must be found [2]. Specifically, the tensile test is a typical technique for the material characterization, which gives access to the average mechanical properties through the resulting strain states. The challenge of investigating the mechanical properties of micro samples consists on emulating the conventional testing processes: manipulating and holding the samples, applying and measuring small forces under uniaxial load, measuring the induced uniaxial strain fields.
Small-scale tensile methods are applied for nano and micro structures (microelectromechanical systems, thin 1 3 films, nanowires, etc.) for investigating the deformation mechanisms at these length scales [3]. Depending on the material, dimension and type of structure, specific techniques are developed for every application. For example, a focused ion beam (FIB) system can be used for the production of dogbone tensile samples of less than 10 µm in diameter and as actuation unit for the tensile testing accompanied by FIB/scanning electron microscopy imaging of the deformation [4]. However, there are currently no methods for holding and stretching spherical micro samples.
Changing the form of metallic samples in the micrometer range can be induced through micro forming [5]. Specific challenges of conventional forming are narrow tolerances necessary for the forming dies and samples, a precise guiding between punch and die as well as the requirement of a high stiffness of the mechanical punch. Therefore, the conventional forming methods are not suitable for economical high-throughput techniques. Electrohydraulic forming is an alternative especially for micro production out of sheet metal [6]. The impulse technology eliminates the fracture problem of high-aspect ratio punches and is rendering the setup more flexible to a varying geometry of die and sample. Electrohydraulic forming is already established in the macroscopic range and is used for reducing or expanding tubes and for processing sheet metal. First investigations regarding forming of bulk samples [7] and micro spheres [2] are reported.
Digital image correlation [8] (DIC) is an established optical method for full-field displacement measurements by correlating subimages of a surface after a small deformation with those of the undeformed surface. The ongoing development of the hardware, fast cameras, computers and algorithms allow the use of DIC for real time observation of dynamic processes [9]. The DIC technique is used for micro-scale (100-500 µm) investigations of deformations in thin films [10] or serrated chips by cutting of titanium alloys [11]. However, this method depends on the existence of detectable features on the surface, naturally occurring (for example, texture of titanium alloy) or introduced through mechanical patterning or paint spraying. Even though a contactless optical measuring method, DIC is therefore in most cases invasive. Moreover, the spatial resolution for deformation detection and strain calculation is limited by the surface feature dimensionality. Laser speckle photography [12] transgresses these limitations and is already used for tensile tests [13] as well as an in-process measuring method for the induced deformation during machining [14,15]. Speckle photography is physically limited only by Heisenberg's uncertainty principle and the respective measurement uncertainty limit can be derived with the Cramér-Rao bound [16,17]. It has therefore the potential for precise small-scale measurements. However, applications for deformation/strain investigations at bulk forming of micro samples are not known.
The principle of emulating the tensile test for spherical micro samples was previously introduced as the combination between the incremental electrohydraulic extrusion as actuation unit and the digital image/speckle correlation (DISC) as strain measuring method [2]. Uniaxial loading of the micro samples is realized by extrusion into deep dies with different forming stages that represent defined stress and strain loads. The speckle patterns have to be measured in situ between the forming steps in order to follow the strain history, which requires an optical access to the process. Measuring the occurring displacement fields in situ with DISC limits the extrusion steps to several 10 μm in order to avoid a strong speckle decorrelation [18]. Moreover, the impulse technique leads to strong vibrations in the actuation and measuring system and therefore to additional speckle decorrelation sources such as speckle contrast loss and defocusing. The following questions arise: Can uniaxial tensile strain be generated in the extruded micro sample? Are the induced displacements measurable with a DISC method? Can a setup provide the necessary spatial resolution and measurement uncertainty for displacements?
In this paper we investigate the measurability through DISC of the deformation induced at micro-extrusion. After describing the proposed measurement principle in Sect. 2, the measurement requirements and the actuation/measurement setup are defined in Sect. 3. The measurement uncertainty for displacements is theoretically assessed in Sect. 4 through Monte Carlo simulations of the speckle pattern measurements and the DISC evaluation. The contributions of the photon shot noise, camera noise and speckle noise are presented depending on the resolution of the evaluation. Moreover, other sources of speckle decorrelation due to the measurement setup and the process characteristics are identified and their influence on the measurement uncertainty is estimated. Section 5 comprises experimental validation of the theoretical results and experimental estimation of the measurement uncertainty due to defocusing. The discussion in Sect. 6 aims to define the optimal parameters of the measurement system, followed by conclusions in Sect. 7.

Principle of the Micro-Extrusion Strain Measurement
Electrohydraulic extrusion is adapted as a micro-forming method which generates regions with uniaxial tensile and compressive strain. These strain fields are measured via an inspection window with DISC as an in situ measurement method, see

Electrohydraulic Incremental Micro Forming
In the electrohydraulic incremental forming the mechanical punch is replaced by a shock wave which allows the extrusion of the micro samples through high aspect ratio channels ( Fig. 1, left side). One approach to generate a shock wave is to explode a wire in a water filled pressure chamber. The required energy is provided by the electrical discharge of a current pulse generator and the resulting reactions in the arising plasma. For the explosion a minimum energy is required, determined by the material, diameter and length of the wire. The resulting shock waves transmitting through water cause a pressure increase up to several giga pascal and convey the punching force within a few microseconds. The shock wave necessary for the extrusion is reached by increasing the loading energy above the vaporization energy of the wire. However, an incremental forming with small consecutive steps is required for the DISC measurements. The size of the extrusion steps is controlled through the shock wave energy by changing the loading energy of the pulse generator, adapted to the material, size and geometry of the sample and of the forming channel. Micro samples covering a wide formability range can be therefore analyzed without changing the extrusion tool by changing the exploding wire characteristics or the provided energy.

Speckle Photography
The functional principle of the image processing (the DISC evaluation) is shown in Fig. 2. The elastic deformations are determined by comparing speckle images of loaded and unloaded workpiece states. This comparison is performed locally by calculating the position of the maximum cross correlation in an evaluation window of size W eval in the field of view ( FOV ). The local surface deformation in the image plane (for example the in-plane displacement in x-direction) is determined from the difference between the position of the cross correlation's maximum and the evaluation window center. In order to obtain the global deformation field, the evaluation window scans with half-window wide steps across the whole field of view, whereby the calculated local displacements x shift , y shift are stored for each scan position (i,j). Due to the almost Gaussian profile of the cross-correlation function, a two-dimensional Gaussian curve is fitted to the correlation data to determine the maximum position with subpixel resolution [19,20].

Spatial Resolution of DISC
For speckle photography the spatial resolution for the displacement distribution is defined as half of the evaluation window width W eval set for the cross-correlation calculations projected on the investigated surface.
The evaluation window width, expressed in pixels in the image plane, is multiplied with the physical width of a detector pixel s pixel and divided by the absolute value of the magnification M of the measurement optical system. The best achievable spatial resolution is however limited by the optical resolution, which is 0.61 ∕NA for a system with the numerical aperture NA and illumination with the laser wavelength . With other words, the minimal reasonable evaluation window is approximately 2 times the size of a fully developed speckle s speckle [12] in the image plane: (2) W eval,min = 2 0.61 * NA |M| s pixel ≈ 2s speckle . Note that by the definition of Goodman [21] for speckles observed through a circular aperture, the speckle size is ~ 92% of the Airy disk projected on the detector for systems with NA ≪ |M| , with NA � = NA∕|M| the numerical aperture on the image side of the optical system. Improving the spatial resolution beyond 1 μm means a numerical aperture NA ≥ 0.247 for = 0.405 μm. Additional improvement of the spatial resolution through a decrease of the physical pixel size of the detector or increase of the imaging system magnification, as suggested by Eq. (1), occurs only as long as the required resolution is larger than the diffraction limit. Otherwise this leads to an increase of the speckle size and of the minimum evaluation window size W eval,min , respectively. Note finally that the magnification is ultimately limited by the minimum required field of view FOV reported to the sensor size s sensor , i.e. M max = s sensor ∕FOV.

Measurement Setup and Requirements
Both the forming process, in particular the geometry of the forming channel, and the optical measurement system must be designed in such a way that tensile strain can be generated and measured in a micro sample. Both areas are subject to individual restrictions and influence each other so that a common optimum must be found.

Measurement Requirements
Spherical metallic micro samples with diameters between 600 µm and 1000 µm are deformed in micro extrusion channels. Constant stress and strain conditions are desired independent of the sample diameter and sample material. Furthermore, manufacturing the micro channels with the required dimensions is a challenge. Finally, an optical access to the deformed sample is required for the in situ DISC measurements. For these reasons simulations of the process are performed to support the design of the extrusion channels which lead to uniaxial tensile strain.

Forming Simulations
A channel geometry which simplifies manufacturing and causes tensile strain fields over large areas independent of spheres diameter is a V-groove with a constant reduction of width and depth along the extrusion path, see Fig. 3a. The flat inspection window acts as a part of the channel. Figure 3b depicts exemplary results of 3D FEM simulations of the extrusion of an aluminum microsphere in a V-groove channel with a 3% width reduction. A bilinear hardening model is assumed for the calculations with the software ANSYS. Positive longitudinal strain in the extrusion direction y results after an extrusion step. Investigation of strain fields with a lateral extension down to 100-200 μm-depending on sample dimensions and position in the extrusion channel-with relatively high gradients require a spatial resolution Δx of less than 10 μm. Due to the incremental procedure the temporal resolution Δt is not influencing the deformation measurability. While in the extrusion direction global displacements up to 100 μm are induced, local displacements in the extrusion and in the lateral direction of only several μm or even less than 1 μm have to be detected. Consequently, a measurement uncertainty for displacements of less than 100 nm is needed.

Extrusion Setup
An oblique 1 mm wide cylindrical channel through the steel die leads the micro sample without modifications to the forming and measuring region. The extrusion channel is realized as a V-groove with the width decreasing from 1 mm at the beginning of the inspection window to 0.5 mm over a length of 7 mm. The optical access is realized by replacing a half of the channel wall with a transparent material which is not interfering with the forming process or impeding on the DISC measurement of the displacement fields at the sample surface. As low friction between the window surface and sample as well as a high dynamic stiffness is required, sapphire is used as window material. Forming tests performed on micro samples with different properties (for example, aluminum Al99.5, steel 100Cr6) showed the ability of 10 mm thick sapphire to withstand even high forming energies, in contrast to borosilicate glass and pre-stressed glass. Less wear and no cracks were observed on the inspection window after the forming steps.

Optical Measurement
Considering the symmetry of strain fields, minimum fields of view FOV slightly larger than 500 μm are necessary in order to detect the strain history over the entire length of the extrusion channel. The optical system is translated to follow the deformed micro sample across the inspection window without changing the illumination conditions. The 10 mm diameter inspection window consists of 10 mm thick sapphire with a refraction index of 1.76. The accessible measurement region is maximized by orthogonally illuminating the sample surface. Therefore, a telecentric objective with inline illumination is used.

Measurement Uncertainty
Simulations of the measurement and the evaluation process of DISC are employed to investigate the measurement uncertainty for displacements of the proposed system. Physical limitation of the optical measurement is the photon shot noise, while camera noise was identified in a previous work [14] as a dominant contribution to the measurement uncertainty. The uncertainty of the DISC evaluation is further influenced by speckle noise, i.e. the influence of the different speckle patterns for the observed surfaces on the cross correlation result. In order to analyze the three contributions to the measurement uncertainty separately, simulations are performed at first for zero displacement.
The shot noise and camera noise induce speckle decorrelation in the simulated speckle images, which is evaluated as non-zero displacement. Speckle noise, however, cannot be estimated in the absence of displacements and other sources of noise since it does not change the speckle patterns when evaluating the cross-correlation of an image pair. Therefore, global shifts up to 3 pixel, in subpixel steps, are introduced in the simulated images before the DISC evaluation. For each type of noise, the measurement uncertainty for displacements is determined as the standard deviation of the calculated displacements.
Speckle patterns with the size of 40 × 40 evaluation windows with statistically varying speckle contrast were numerically generated for speckle sizes up to s speckle = 10 pixel. The pseudo-random normally distributed phasor matrices which are used for the simulations lead to patterns of fully developed speckle [21]. Averaging over the large number of windows reduces the influence of the speckle noise on the calculated measurement uncertainty due to shot noise and camera noise. The speckle images were normed to the full well capacity FWC = 12300 e − of the camera CP70, which was also used in the measurement setup (Sect. 3.3). With other words, the speckle patterns are normed to full illumination independent of illumination time. For zero displacement, each of these patterns was overlaid 40 times with Poisson distributed photon shot noise and camera noise, respectively. Camera noise is simulated as constant white Gaussian noise. The resulting differences between the induced and the calculated average displacements are mostly below 10 −6 pixel or at least 2 orders of magnitude below the respective standard deviation, meaning that no systematical contribution was introduced to the measurement uncertainty through the simulation and evaluation procedure. This way, up to 40 × 40 × 40 = 64000 measurement points for each speckle/evaluation window size are generated. Note that averaging over the ensemble of simulated images means averaging over the speckle contrast. The limits of measurability and the measurement uncertainty budget presented in the next two subsections are evaluated for an average speckle contrast and for zero or 1D-displacements y shift .
The displacement measurements are influenced also by the optical access and the process characteristics. Scattering at the inspection window interface to the material, volume scattering in sapphire or traces of water which was pressed by the shock wave between the window and the sample are phenomena leading to an overall contrast loss. The influence of the speckle contrast on the total measurement uncertainty is discussed in Sect. 4.3. As a result of extrusion, 2D strain fields are generated in the measurement plane (x, y) . Therefore, the influence of a displacement component x shift which is perpendicular to the evaluation direction on the measurement uncertainty of y shift is additionally theoretically estimated in Sect. 4.3. An important decorrelation source for this application is the defocusing ( z shift ≠ 0) in the optical system subjected to strong excitation of the setup during the extrusion. This effect cannot be simulated and will be treated experimentally in Sect. 5 after the validation of the theoretical considerations of Sect. 4.

Limits of Measurability due to Photon Shot Noise
For fully developed speckle patterns, the fundamental limit of the measurement uncertainty due to photon shot noise is expressed by the square root of the Cramér-Rao bound 1 3 where N total represents the total number of photons involved in the process [16]. The measurement uncertainty for an ideal speckle photography system is 0 = √ CRB , which is therefore proportional to the speckle size s speckle . However, Eq. (4) is only valid for speckles larger than a pixel, since it does not consider pixel quantization effects. Moreover, it is valid for evaluation windows with a large number of speckles. With this conditions, the total number of photons N total is expected to vary proportional to the area of the evaluation window W 2 eval or, for a given speckle size, to the average number of speckles evaluated as N speckle = W eval ∕s speckle 2 .
This behavior is illustrated in Fig. 4 for speckles with s speckle = 5 pixel, where N total is calculated as the average number of photons illuminating each of the simulated speckle images which were overlaid only with Poisson distributed noise and considering a quantum efficiency of 1.
The measurement uncertainty 0 varies linearly with N speckle in a double logarithmic representation, with a slope of − 1/2 (dashed line) for N speckle > 10. The obtained deviation of 0 from this trend for N speckle < 10 is fully in agreement with the findings in [16]. Additionally, the measurement uncertainty due to the photon shot noise shot is calculated from the speckle decorrelation of the simulated speckle images, as described in the 2 nd paragraph of Sect. 4. The difference between the empirical shot and the analytical 0 , as seen in Fig. 4, is due to averaging over many patterns with statistically varying speckle contrast. The minimum value of the measurement uncertainty, an average over the lowest 10 values of the measurement uncertainty due to the photon shot noise, is min,shot ≈ 0 as expected since the CRB describes the lower limit of measurability. Note that the CRB is calculated using average values for N total , which explains the values of min,shot < 0 .
For N speckle > 4, shot is also proportional to √ N speckle . The deviation from this trend of shot at N speckle < 4 roots in the cross-correlation algorithm. The cross-correlation peak loses in amplitude and sharpness and finding the position of the maximum is subject to a larger uncertainty. Therefore, an increase at small N speckle is expected for all components of the measurement uncertainty. Moreover, with only a few speckles in the window, the influence of the speckle noise is enhanced. Note, however, that reducing the size of the evaluation window below 2s speckle [cf. Eq. (2)], which means N speckle < 4, does not lead to an increase in spatial resolution and, thus, is not applied here.

Measurement Uncertainty due to Speckle Noise
Increasing the spatial resolution for the displacement measurements by choosing a small evaluation window leads to a reduced number of speckles that contribute to the correlation. The influence of the average number of speckles N speckle on the measurement uncertainty speckle (due to speckle noise) is shown in Fig. 5 over different displacements for the case of s speckle = 5 pixel. The measurement uncertainty speckle is determined as the standard deviation of the calculated displacements for 100 × 100 evaluation windows for each N speckle between 2 2 and 10 2 . For full-pixel Measurement uncertainty speckle as a function of the displacement y shift for an average speckle size s speckle = 5 pixel. The variation of the average number of speckles N speckle in the evaluation window W eval is indicated displacements, the cross-correlation algorithm reproduces exactly the displacement due to the absence of other sources of noise and the measurement uncertainty is zero. For subpixel displacements, the maximum of the cross-correlation peak is found through a 3 pixel × 3 pixel Gaussian fit around the pixel with the highest value. Therefore the results of the subpixel interpolation are subject to a larger variance as the shift of the peak approaches half a pixel, which causes the highest asymmetry of the 3 × 3 matrix around the cross-correlation maximum. The effect of this asymmetry is strongly enhanced for the broad cross-correlation peaks which occur at N speckle ≤ 4, as can be observed in Fig. 5. Note that this behavior of speckle noise repeats itself for each subsequent full-pixel interval of displacements (not shown here). Figure 6 shows the influence of N speckle on the measurement uncertainty due to the shot, camera and speckle noise at its maximum, for half a pixel displacements, cf. . For displacements in the range of half a pixel ( y shift = 0.208 μm, for a camera with s pixel = 6.45 μm and M = 15.51), speckle dominates the measurement uncertainty over the whole simulated range of N speckle . For displacements < 0.3 pixel and N speckle > 4, however, the total measurement uncertainty is dominated by camera (not shown in Fig. 6). At very small numbers of speckle in the evaluation window ( N speckle < 4) speckle is again the major component of the measurement uncertainty, independent of the displacement size. Figure 7 shows the influence of the speckle size on the total measurement uncertainty for different average N speckle (grey curves). For s speckle ≳ 4 pixel, total depends mainly on N speckle . Below this value, quantization effects diminish the efficiency of the correlation algorithm. However, for small evaluation windows, an optimum between the effects of s speckle and N speckle must be found. For illustration, Fig. 7 also depicts the isocurves for W eval = 30 pixel, 20 pixel and 10 pixel, respectively. The measurement uncertainty displays an increase for s speckle below 2 pixel due to undersampling. Note that the cross-correlation algorithm of DISC does not work for s speckle < 1 pixel. The minimum uncertainty is reached between 2 and 3 pixel and increases with increasing s speckle , which agrees with the findings in [14].

Decorrelation
An important aspect influencing the measurement uncertainty is the speckle contrast, which is usually defined as C = I ∕I , the ratio between the standard deviation and the average value of the intensity in the image [21]. This definition of contrast is however valid only for large images [22] and is therefore not applicable for small evaluation windows. A more relevant parameter for DISC is the cross-correlation coefficient , which is also influenced by the deviation of intensity from the expected (average) value [21]. A decorrelation factor according to [23] is defined as Fig. 6 Measurement uncertainty as a function of the average number of speckles N speckle in the evaluation window W eval for an average speckle size s speckle = 5 pixel. speckle is evaluated for a displacement y shift = 0.5 pixel The influence of speckle decorrelation on the measurement uncertainty total is shown exemplarily in Fig. 8 for W eval = 10 pixel and s speckle = 4 pixel, which means an average N speckle = 6.25. For these results, speckle patterns with large variations of contrast were generated and global shifts y shift up to 2 pixel were induced. The images consisted of 100 × 100 evaluation windows, each image with different W eval and different s speckle . Both the original and the shifted patterns were overlaid with shot and camera noise. For zero and full-pixel displacements (Fig. 8a) total depends linearly on with a slope which is not affected by the displacement value. For s speckle ≥ 4 pixel the slope decreases with increasing N speckle , meaning that for evaluation over larger numbers of speckle, the measurement uncertainty is less influenced by decorrelation (not shown here). At = 0 the measurement uncertainty approaches 0 , since this implies best correlation ( = 1 ) and best contrast.
In the presence of subpixel displacements (Fig. 8b), the linear behavior of the measurement uncertainty with is only valid for very small displacements y shift ≲ 0.1 pixel. For subpixel displacements 0.2 pixel < y shift ≤ 0.5 pixel the effect of finite pixel size dominates the measurement uncertainty, which no longer depends on in the simulated range ( ≤ 0.4 or ≥ 0.86). This means that an improvement of the correlation beyond a certain limit does not lead to a significant decrease of the measurement uncertainty. A slight increase of the measurement uncertainty and average with full-pixel value of y shift can be observed in Fig. 8c (red filled markers), a trend which is followed also for the corresponding subpixel y shift . The variation of total versus average , respective versus y shift , mirrors the variation of speckle (Fig. 5) combined with a displacement-proportional increase.
The extrusion process leads to large in-plane displacements at the interface to the inspection window. Therefore, the influence of large displacements y shift of several µm and of non-zero perpendicular components x shift on the speckle decorrelation and the measurement uncertainty for displacements has to be estimated. Figure 9a shows the overall increase of average over a large range of displacements y shift considering also simultaneous perpendicular displacements x shift . The decorrelation induced by the finite pixel size is enhanced in the presence of arbitrary oriented in-plane displacements. The increase of total with y shift is stronger than that of average , even for x shift = 0 (not shown here), because for some of the 10 × 10 pixel 2 evaluation windows Fig. 8 Measurement uncertainty total as a function of the decorrelation factor for different displacements y shift , indicated in pixel by the numbers in the graphs. An average speckle size s speckle = 4 pixel and W eval = 10 pixel ( N speckle = 6.25) is considered. a full-pixel displacements up to 2 pixel; b subpixel displacements up to 0.5 pixel; c total versus the average for displacements up to 2 pixel in 0.1 pixel steps (full-pixel displacements are marked with red, the subpixel fraction is specified on top of the graph) Fig. 9 a Average decorrelation factor average versus the y-component y shift of a displacement with the x-component x shift = 0, y shift ∕10 and y shift ∕5 ; b average and measurement uncertainty total versus y shift for x shift = y shift ∕5 . An average speckle size s speckle = 4 pixel and W eval = 10 pixel ( N speckle = 6.25) is considered the correlation is completely lost ( ≈ 0). Moreover, the measurement uncertainty total follows only tendentially average for complex displacements, as shown in Fig. 9b for the measurement uncertainty of the y-component y shift of a displacement with the x-component x shift = y shift ∕5 . Note that total is considered as the standard deviation of the calculated y shift over all the 10000 evaluation windows with a large contrast variation. However, an optimized speckle contrast leads to decreased measurement uncertainties. Selecting only the windows with the correlation coefficient > 0.8 ( < 0.5) for the evaluation, not only that average and total decrease, but they are also no longer significantly influenced by the full pixel value of y shift . For evaluation windows with N speckle > 6.25, the resulting measurement uncertainty allows for displacement measurements with even lower correlation than 0.8, particularly for displacements of several pixels and larger.
For the example of a camera with a pixel width s pixel = 6.45 μm and an optical system with a magnification M = 15.51, as also used in the experimental section, an image pixel corresponds to 0.416 µm. A 1 µm displacement in y-direction can be therefore measured with a measurement uncertainty between 30 nm and 90 nm at a spatial resolution Δx = 2.1 µm, depending on the contrast. For comparison, for a camera with s pixel = 3.45 μm (1 pixel = 0.222 µm) the measurement uncertainty for the same displacement lies between 16 nm and 50 nm at the same Δx as above. The theoretical measurement uncertainty lies below the required 100 nm. For displacements of several µm total may reach larger values (up to 300 nm). However, the relative measurement uncertainty remains below 5%, which would be sufficient for the displacement measurements.
Constant illumination conditions are assumed for the simulated displacements. This assumption is not necessarily valid for larger deformations, which could lead to higher speckle decorrelation than estimated. Another important decorrelation source for this application is the defocusing in the optical system subjected to strong excitation of the setup during the extrusion. These effects cannot be accessed by simulations with reasonable efforts and will be treated experimentally in the next section, after validation measurements for the theoretical considerations concerning the measurement uncertainty. Figure 10 depicts the measurement uncertainty for displacements as resulting from DISC simulations total and from validation measurements measurement for zero displacements. For imaging, the camera CP70 (FWC = 12300 e − , parameter also used for simulating the camera noise) and an imaging lens (Apo Rodagon) with adjustable aperture (f/4 to f/22) were employed. A laser with = 0.405 µm illuminated the surface of a reference sample with a Gaussian beam profile at an incidence angle of about 30°. Speckle patterns with s speckle between 2 and 10 pixel were generated by varying the aperture. Both the simulated and the measured speckle patterns were evaluated with W eval = 30 pixel in order to increase the range of accessible speckle sizes. The measurement results correspond very well with the simulated data, for average decorrelation (considering all evaluation windows) and for lowest decorrelation (considering only ten windows with the highest value of , not shown in Fig. 10) as well. Note that the simulated data are also presented in Fig. 7. For completion, the measurement uncertainty 0 due to CRB and shot as functions of the speckle size are shown.

Validation Measurements
The reference measurements were repeated with the camera DX4 -285 FW (FWC = 23000 e) and the telecentric objective with NA = 0.23 and inline illumination ( = 0.405 µm). A magnification M = 15.51 was achieved by using a 4 cm extension tube. The DISC evaluation of the measurements for both optical systems was realized with an evaluation window of 10 pixel. Figure 11 shows that for s speckle of about 4 pixel the linear behavior of the measurement uncertainty with the decorrelation is independent on the measurement optical system. Moreover, by comparing with the simulation results, an only slightly larger slope can be calculated for the measurements (0.18 pixel −1 ) than for simulations (0.16 pixel −1 ).
Thus, the estimated measurement uncertainty and its linear behavior with the speckle decorrelation is validated for zero displacement. Subsequently, the influence of large and 2D displacements on and total is experimentally investigated and compared with the theoretical results.

Displacement Measurements
On a reference sample global in-plane displacements were induced using a P-611.3S NanoCube xyz-nanopositioning system (PI GmbH & Co. KG). The measuring system consists of the DX4 -285 FW camera and the telecentric objective with M = 15.51. Displacements in the y-direction up to 10 µm in 0.1 µm steps, meaning 0.24 pixel steps, were realized for three different x positions at a 0.2 µm interval.
The DISC evaluation showed that the displacements had an x-component x shift ≈ y shift ∕12 due to a slight misalignment between the y-axis of the NanoCube and that of the camera. This result is illustrated by the modulation of the average decorrelation factor (Fig. 12a), which is a consequence of the finite pixel size. For W eval = 10 pixel (spatial resolution Δx = 2.1 µm), average for the measurements is slightly higher than that for simulations, indicating regions of the measured speckle patterns with lower correlation (and lower contrast) than realized by simulations. The influence of the speckle decorrelation on the measurement uncertainty of y shift is mirrored in Fig. 12b by the strongly increased standard deviation especially for half-pixel x shift Since N speckle = 6.6 (for s speckle = 3.9 pixel), the behavior of the speckle noise (Fig. 5) is therefore reproduced. Selecting only the windows with the correlation coefficient > 0.8 for the evaluation, the uncertainty is strongly reduced while still showing the same behavior (not shown here), rendering the measurement suitable for the application in the investigated displacement range. An alternative for reducing the measurement uncertainty is the resolution decrease. For comparison, the evaluation with W eval = 40 pixel (Fig. 12c) is no longer influenced by the finite pixel size and the speckle noise. The observed increase of the standard deviation is explained by the constant change of the optical way for each surface element contributing to the speckle formation, leading to a progressive change of the speckle pattern. For displacements ≳ 9 µm, the relative measurement uncertainty is larger than 10% and therefore no longer suitable for our application. Note that the presented standard deviations mirror the uncertainty of a single measurement.
The influence of speckle decorrelation at high resolution is particularly relevant for subpixel displacements. Figure 13 shows the measurement uncertainty as a function of the decorrelation factor for different measured displacements y shift in comparison to the results of simulations (Fig. 8). The subpixel value of y shift determines the minimum achievable decorrelation. For y shift = 0, the expected linear behavior is observed, however with a slightly larger slope than resulting from simulations. Up to a threshold which is varying with the non-zero displacement value, the measurement uncertainty is not significantly influenced by the change of . Beyond this threshold, total follows qualitatively the slope shown for y shift = 0. Therefore, the relative measurement uncertainty of ≲ 15% cannot be reduced only by improving the speckle correlation and a different approach for compensating the speckle noise must be found.

Defocusing
With the same measurement setup as described in the previous subsection, out-of-plane global shifts z shift of ± 10 µm, with 0.1 µm steps, were also realized in order to investigate the influence of defocusing on the in-plane displacement measurement. Figure 14a shows the decorrelation effect of defocusing, at first for y shift = 0. The rapid increase of average for z shift ≠ 0 is accompanied by a strong increase of the standard deviation for the evaluation of y shift (Fig. 14b) with N speckle = 6.6 ( W eval = 10 pixel, s speckle = 3.9 pixel). Moreover, a systematic effect of defocusing can be observed as the average measured y shift deviates significantly from zero beyond z shift = ± 4 µm ( y shift = 0.5 µm at z shift = ± 6 µm). Speckle displacements can be expressed through a term directly proportional to object inplane displacements and a term proportional to the defocusing [24], which explains the systematic error. Note however that measurement uncertainties below 100 nm can only be reached in the range of negligible systematic error, region marked in the figure by the green rectangle. Considering only the windows with > 0.8 for the cross-correlation evaluation diminishes the systematic error, but the measurement uncertainties do not decrease relevantly and the measurability region is therefore not enlarged (not shown here). The systematic influence of defocusing is reduced at lower spatial resolution, as shown in Fig. 14c for W eval = 40 pixel ( N speckle = 105.2). Additionally, the measurement uncertainties are also reduced leading to a larger measurability region. The so defined measurability region is determined by the measurement uncertainties due to defocusing in the presence of noise (shot, camera and speckle noise) at zero in-plane displacement. For non-zero (subpixel) displacements, the measurability region is smaller. Defocusing with 1 µm changes only slightly the variation of the measurement uncertainty with decorrelation, as illustrated in Fig. 15 for an in-plane displacement y shift = 0.1 µm (0.24 pixel). For z shift = 2 µm, however, a relative measurement uncertainty up to 30% results, independent of the decorrelation factor. The measurability region is reduced by a factor of 2 as compared to Fig. 14b. With other words, at the spatial resolution Δx = 2.1 µm the in-plane displacement measurement is very sensitive to defocusing.

Discussion
Simulations of the forming process showed that longitudinal tensile strain can be generated in micro samples at the incremental electrohydraulic extrusion depending on the shape of Fig. 13 Measurement uncertainty as a function of the decorrelation factor for different displacements y shift . Measurements with the DX4 -285 FW camera ( s pixel = 6.45 μm) and s speckle = 3.9 pixel and simulations with s speckle = 4 pixel; for both measurements and the simulation W eval = 10 pixel Fig. 14 a Average decorrelation factor average versus the induced out-of-plane displacements z shift for y shift = 0. Average of the measured y-displacements y shift and standard deviation for evaluation with b W eval = 10 pixel ( N speckle = 6.6) and c W eval = 40 pixel ( N speckle = 105.2). Green rectangles mark the measurability region at zero displacement the extrusion channel. Furthermore, the displacement measurement requirements could be determined with the simulations (Sect. 3.1). In-plane displacements of from less than a μm to several μm have to be detected with a measurement uncertainty < 100 nm and a spatial resolution Δx < 10 µm. Theoretical and experimental investigations of the measurability of the resulting deformations with the DISC method, respectively its limitations, are presented in the Sects. 4 and 5.
The questions left open are whether the proposed optical setup fulfils the described measurement requirements and which restrictions for the extrusion process result from the limitations of the measurement method.
As shown in the previous sections, the achievable spatial resolution is typically limited at Δx ≈ 3 µm for defocusing below 1.5 µm and in-plane maximal displacements of less than 10 µm. As the spatial resolution is not an issue for the sample translation at extrusion, global displacements of several tens of micrometers can be accurately measured from the speckle images and subtracted from the displacement field of the micro sample surface. However, the resulting global speckle decorrelation adds to that induced by the local displacements, cannot be compensated by the algorithm and impedes therefore the evaluation with high spatial resolution of the deformation. Measurements of the deformation with a spatial resolution Δx ≈ 10 µm are less sensitive to the global displacements and can be realized with the present setup also for images with lower speckle correlation. A measurement strategy which is currently tested involves adjusting locally the spatial resolution of the displacement evaluation to the calculated correlation coefficient.
Measurement uncertainties below 100 nm at high spatial resolution can be reached with the chosen setup. Particularly for small displacements ≤ 200 nm, however, the achievable spatial resolution has to be limited to Δx ≥ 6 µm in order to provide relative measurement uncertainties ≤ 10%. Therefore, the achievable spatial resolution of the displacement measurement is determined by the smallest detected displacements. Alternatively, the spatial resolution can be locally adjusted during the evaluation depending on the measured displacement.

Conclusions
The realization of the tensile test for spherical micro samples by the combination of the incremental electrohydraulic extrusion as actuation unit and the digital image/speckle correlation as deformation measuring method was investigated theoretically (through simulations of the forming process and of the speckle measurement) and experimentally. Particularly the measurability of the local displacements with a high spatial resolution was studied in detail to identify important sources for the measurement uncertainty. A dominant effect influencing the measurement uncertainty at high spatial resolutions Δx ≤ 6 µm is the speckle noise. Regarding the required spatial resolution, the theoretical and experimental results demonstrate the feasibility of the deformation measurement approach. With the presented experimental setup, local displacements of less than 10 µm induced by extrusion steps of ≈ 10 µm (global displacements) can be measured with a spatial resolution Δx ≥ 3 µm. For the measurement of small displacements (≤ 200 nm) a decrease of the measurement uncertainty is needed, which can be realized by minimizing the influence of the speckle noise through averaging over different speckle patterns. As a next step, the characterized displacement measurement system for micro samples is applied for investigating the strain fields which are generated at extrusion depending on the mechanical properties of the micro samples.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.

Fig. 15
Measurement uncertainty as a function of the decorrelation factor for different displacements y shift and different out-of-plane displacements z shift . Measurements with the DX4 -285 FW camera, s speckle = 3.9 pixel and W eval = 10 pixel