New Four Points Initialization for Digital Image Correlation in Metal-Sheet Strain Measurements

Nowadays, the deformation measurement in metal sheets is important for industries such as the automotive and aerospace industries during its mechanical stamping processes. In this sense, Digital Image Correlation (DIC) has become the most relevant measurement technique in the field of experimental mechanics. This is mainly due to its versatility and low-cost compared with other techniques. However, traditionally, DIC global image registration implemented in software, such as MATLAB 2018, did not find the complete perspective transformation needed successfully and with high precision, because those algorithms use an image registration of the type “afine” or “similarity”, based on a 2D information. Therefore, in this paper, a DIC initialization method is presented to estimate the surface deformation of metal sheets used in the bodywork automotive industry. The method starts with the 3D points reconstruction from a stereoscopic digital camera system. Due to the problem complexity, it is first proposed that the user indicates four points, belonging to reference marks of a “Circle grid”. Following this, an automatic search is performed among the nearby marks, as far as one desires to reconstruct it. After this, the local DIC is used to verify that those are the correct marks. The results show reliability by reason of the high coincidence of marks in experimental cases. We also consider that the quality of mark stamping, lighting, and the initial conditions also contribute to trustworthy effects.


Introduction
The quality control of mechanical stamping is a need in modern industries.Engineers are always looking to improve the mechanical work optimization.Specifically, in the mechanical industry of metal-sheets manufacture, it is necessary to know their mechanical properties to use them in the most suitable way [1,2].In this industry, the material properties such as Young's modulus, Poisson's ratio, anisotropic plastic ratio parameters and others are required for improving the design and manufacture processes by finite element analyses.On the other hand, the most commonly used method for finding the material properties is a tensile test with a strain-gauge type extensometer.However, the results from the extensometer are not applicable to measure the strain at a local point and the onset of diffuse necking.In contrast, the Digital Image Correlation (DIC) method is a state-of-the-art technique that can be used for an accurate strain measurement of material properties [3].DIC has become a very popular technique in mechanics, particularly for measurement materials deformation without contact.It is an economic and versatile method which also offers a big amount of experimental data [4].As digital cameras technology constantly improves in terms of resolution, lenses and frames per second (FPS), this technique allows covering a wide range of scale in space and time.DIC in stereo systems even reach non-planar cases, very common in the industry, meaning the three-dimensional reconstruction and access to volumetric data about the material and structure behavior through time.Nonetheless, the use of DIC data to validate models in a quantitative way or to identify with precision several constitutive parameters, remains an open problem.One of the reasons for this is the complicate compromise between the measuring resolution and the large space for the measuring to be done.A second reason is the state of the frontiers.A third reason is that the measured displacements are not directly compared with common simulations.Of course, there are several efforts presenting promising results, as well as [5], a work analyzing the damage in a super pression aerostatic balloon (SPB) used by the Japanese Aerospace Exploration Agency (JAXA) during the insufflate prove.The expansion, when it is not uniform, generates a concentration of stress and makes the balloon explode.They proved a new method called the Simplified Digital Image Correlation Method (SiDIC), in a rubber balloon, to verify if the SPB could be measured.The SiDIC identified in a correct way the non-deformed region even though the deformation was not precise.
Below, a brief historical background of DIC is given.An open code 3D-DIC Toolbox for MATLAB, the one implemented algorithm for multiple cameras, inspired by biomechanics, was reported in [6].They were verified using the Multiple Digital Image Correlation Method (MultiDIC) and low-cost hardware.In [7], a quick method to correlate digital images (3D-DIC) was proposed to implement real-time measuring.Two improvements came up: the development of an efficient algorithm of the inverse composition of Gauss-Newton (IC-GN) with parallel computing to avoid redundant calculation, and an efficient IC-GN algorithm which was used to reach a speed of 10 frames/second with a resolution of 5000 points per frame.For validation, the displacement field of a four-point bending beam was determined by the real-time 3D-DIC.The experimental results were verified in traditional Chinese medicine.In the publication of [8], a procedure was developed to generate a sequence of intermediate synthetic images for a gradual following of the transformation of the pixel's subset between the two external configurations.An adequate distortion function of image was defined in the whole image through adopting an algorithm based in characteristics and followed by a scheme of interpolation based in non-uniform rational B-spline (NURBS).This allowed a fast and trustworthy estimation of the initial deformation parameters for the refinement phase after the DIC analysis, and tests were made on aluminum plates.Document [9] proposed a method of deconvolution to recover the real fields of displacement and tension of their counterparts provided by local DIC.The proposed algorithm could be considered an extension of the Van Cittert deconvolution, based on the small tension assumption.The authors claimed an improvement on the fine details in displacement and tension maps, as well as in space resolution.The guidelines to evaluate the precision and resolution of the 2D analysis under the patronage of the Society of Experimental Mechanics were reported in [10].
On the other hand, algorithms to map relative displacements of the deformed material points in opposition to non-deformed material by using DIC were reported in [11,12].The current possibilities for the DIC scale allows the study of deformation on different levels, from meters to the nanoscale [13,14], with the condition of recording it correctly, following a known pattern like a circle grid or square grid.Additionally, studies have looked at how different systems behave, like biological materials [15][16][17][18], metallic alloys [19][20][21], memory-shape alloys [22,23], porous metals [24][25][26], polymers [27] and polymeric foams [28].A key step in the process of DIC image tracking is the definition of the initial assumption, for the non-linear optimization routine aimed at finding the parameters that describe the transformation of the subset of pixels.This initialization can be very challenging and possibly fail when it comes to pairs of highly deformed images, such as those obtained from two angled views of non-planar objects.For example, in Figure 1, it can be seen how the global DIC cannot exactly match the circle mark to measure the strain.This was performed with the MATLAB function called "imregcorr", which was reported in [29].As can be observed from the preceding discussion, there still is the question of whether it is possible to improve the DIC accurateness.In brief, these previous works highlight the need for developing a systematic approach for achieving this goal, which to the author's awareness has not yet been stated.Therefore, in this paper we propose the use of a mark-search as an initial condition and of DIC as the system to verify the point-pairing.We propose the use of a mark-search as an initial condition and the use of DIC as the system to verify the point-pairing, because the number of mechanical part images are plentiful, the search for features to measure is very broad, as a result of which providing the initial conditions is necessary to fully automate the process gradually.The lack of a correct initial condition generates slowness and inaccurateness in current implementations of DIC.An appropriate initial condition is especially useful to know the material properties, specifically, its deformation curve and characterization zones such as: elastic, inelastic and rupture [30].The proposed method in this work is applied to galvanized steel sheet, used in the outside bodywork of a pickup truck at 620 µm-thick and 96.5 µ hardness.Figure 2 shows the surface view of the material obtained from a metallographic microscope model:AX @ 200X, Company: Carl Zeiss AG, Jena, Turinga, Germany.It is necessary to mention that the experiments presented in this paper were all made for the automotive industry sector.As can be observed from the preceding discussion, there still is the question of whether it is possible to improve the DIC accurateness.In brief, these previous works highlight the need for developing a systematic approach for achieving this goal, which to the author's awareness has not yet been stated.Therefore, in this paper we propose the use of a mark-search as an initial condition and of DIC as the system to verify the point-pairing.We propose the use of a mark-search as an initial condition and the use of DIC as the system to verify the point-pairing, because the number of mechanical part images are plentiful, the search for features to measure is very broad, as a result of which providing the initial conditions is necessary to fully automate the process gradually.The lack of a correct initial condition generates slowness and inaccurateness in current implementations of DIC.An appropriate initial condition is especially useful to know the material properties, specifically, its deformation curve and characterization zones such as: elastic, inelastic and rupture [30].The proposed method in this work is applied to galvanized steel sheet, used in the outside bodywork of a pickup truck at 620 µm-thick and 96.5 µ hardness.Figure 2 shows the surface view of the material obtained from a metallographic microscope model:AX @ 200X, Company: Carl Zeiss AG, Jena, Turinga, Germany.It is necessary to mention that the experiments presented in this paper were all made for the automotive industry sector.As can be observed from the preceding discussion, there still is the question of whether it is possible to improve the DIC accurateness.In brief, these previous works highlight the need for developing a systematic approach for achieving this goal, which to the author's awareness has not yet been stated.Therefore, in this paper we propose the use of a mark-search as an initial condition and of DIC as the system to verify the point-pairing.We propose the use of a mark-search as an initial condition and the use of DIC as the system to verify the point-pairing, because the number of mechanical part images are plentiful, the search for features to measure is very broad, as a result of which providing the initial conditions is necessary to fully automate the process gradually.The lack of a correct initial condition generates slowness and inaccurateness in current implementations of DIC.An appropriate initial condition is especially useful to know the material properties, specifically, its deformation curve and characterization zones such as: elastic, inelastic and rupture [30].The proposed method in this work is applied to galvanized steel sheet, used in the outside bodywork of a pickup truck at 620 µm-thick and 96.5 µ hardness.Figure 2 shows the surface view of the material obtained from a metallographic microscope model:AX @ 200X, Company: Carl Zeiss AG, Jena, Turinga, Germany.It is necessary to mention that the experiments presented in this paper were all made for the automotive industry sector.

Theoretical Bases
The dimensional estimation using stereoscopic vision systems was described in detail by Tsai in [31], who quantitatively concluded that the accuracy and precision relied on the distance between the camera and the object to be measured, as well as the resolution of the cameras used.Stereoscopic vision is a technique frequently used to locate points in three dimensions (3D) based on points in two or more 2D images [6,32].To achieve a stereoscopic point triangulation, it is necessary to calibrate both cameras.In this paper, we used the PinHole model, which is described by Equation (1) [33]: since s is the number that defines the scale of the objects with respect to their real size in the image, and [R|t] is the matrix of extrinsic parameters of the camera that describe the rigid transformation (rotation and translation) between the coordinate system of the camera and the coordinate system of an object outside the camera.A is the matrix of intrinsic parameters that describes the position of the center of the image in pixels (u 0 , v 0 ), and the ratio of the size of the pixel ( f x, f y) has units expressed in pixels meters in the axes x and y regarding the focal distance between the entrance hole of the light and the matrix of the light sensors of the camera.Parameter c describes the asymmetry of the two axes of a pixel where a zero expresses an angle of 90 degrees.On the other hand, the entry ) in homogeneous coordinates of the scene or object expressed in meters in its own coordinate system, and m = is the corresponding 2D point in the image expressed in pixels (u, v).Finally, the correction of the lens distortion is made with a polynomial based on the idea that the distortion changes as a circumference, as in Equations ( 2) and (3): Subsequently, to achieve a stereoscopic calibration, the translation vector that joins each camera reference system is calculated [34,35].The point position in three dimensions can be estimated from the coordinates in two dimensions and the Equation ( 4) [33,[36][37][38]: where (u, v) and (u , v ) are the coordinates of the paired points of the left and right cameras that correspond to the 3D point to be reconstructed.The scalars m ij are obtained by multiplying the intrinsic and extrinsic parameters of the left camera A[R|t]; in an analogous way, the scalars m ij are obtained by multiplying the intrinsic and extrinsic parameters of the right chamber A [R |t ].

Normalized Cross Correlation (NCC)
Points matching, which consists of locating a point (photographed from different positions) in the left photograph and in the right image, is necessary to carry out the triangulation process.2D-DIC is used to correlate a given set of points in the two stereo views of the reference configuration and match these points along the sequence of images.The correlated image points that are set are used to reconstruct and track the 3D position of the material points of the ROI (region of interest) through time [6].In the literature, there are several works that report using 2D-DIC to accomplish this process, because this technique is so meaningful for 3D point reconstruction based on 2D images.The first advantage is that cross correlation is equally simple to calculate.Once a coincidence for a patch in a typical position inside an image is achieved, the Fourier methods can be used to calculate the correlation in a fast way.The second advantage is that cross correlation is independent from translations and invariant in the scaling dominance of the intensity.Equation ( 5) presents the definition of a normalized cross correlation: Here, f and g are each the grayscale functions of the windows of the current image at a specific location (x, y).The functions f and g correspond to the gray scale mean of the reference image and the current subset.

Superficial Strain Estimation
To estimate the level of superficial strain, the engravings of speckle, circular or square marks on the sheet are used, which involves an electrochemical process in which an electrolyte is applied as a reagent on the surface of the sheet to engrave a thin pattern so that when the sheet is painted, the marking is imperceptible.The use of a circles pattern was chosen because it is considered a non-uniform mark deformation.In this proposed method, the facet sizes (window sizes) for the correlation function are defined by the distance between the circles' centroids.The advantages and disadvantages of a circle pattern in contrast with a speckle pattern are show in Table 1.
Table 1.The advantages and disadvantages of circle pattern versus speckle pattern.

Characteristics Speckle Pattern Circle Pattern
The normalized correlation matches well with the patterns.YES NO Mark deformation is assumed as irregular.
NO YES A pixel can represent a mark.Some cases YES Independence of the distance between the cameras and the specimen.

Some cases YES
After the sheet is subjected to the sausage process, the circles are distorted in the form of ellipses and the distances between their centroids are modified; the latter are used in the digital images to determine the deformation states in the testing metal sheet.Figure 3a,b shows the changes of the circular grid with the deformation of the metal sheet.A simple way to estimate the deformation in the centroid of each ellipse is by averaging the distance with its four neighboring centroids Ĉk−1,l , Ĉk,l−1 , Ĉk+1,l and Ĉk,l+1 , as expressed by Equation ( 6): where d Appl.Sci.

Materials
To capture the images, a pair of Posilica GT 2750 (Allied vision, Exton, PA, USA) cameras with a sensitivity in the range of 400 nm to 670 nm as used.It has ethernet communication, a maximum capture rate of 33 fps, and resolution of 2750 to 2200 pixels at 230-250 lux.They were mounted on a metal structure that allows lighting control.Figure 4 shows the arrangement of the cameras, the light source and the metal test piece used to test the proposed method.LED lighting was chosen due to the contrast with the circles marked on the sheet: a blue LED of 640 nm, which coincides with the work [39].To compare the dimensional measurements, a digital microscope (Jiusion 6-06814-24289-8) with a 100 µm scale was used.

Materials
To capture the images, a pair of Posilica GT 2750 (Allied vision, Exton, PA, USA) cameras with a sensitivity in the range of 400 nm to 670 nm as used.It has ethernet communication, a maximum capture rate of 33 fps, and resolution of 2750 to 2200 pixels at 230-250 lux.They were mounted on a metal structure that allows lighting control.Figure 3 shows the arrangement of the cameras, the light source and the metal test piece used to test the proposed method.LED lighting was chosen due to the contrast with the circles marked on the sheet: a blue LED of 640 nm, which coincides with the work [39].To compare the dimensional measurements, a digital microscope (Jiusion 6-06814-24289-8) with a 100 µm scale was used.

Method
The proposed method has the purpose of measuring the surface deformation in a metal sheet specimen used for truck exteriors, and it was divided into ten stages from (a) to (j).The peculiarity of this procedure consists of the incorporation of an initial condition based on points normalization to perform the matching using 2D-DIC.The general procedure is as follows:

Method
The proposed method has the purpose of measuring the surface deformation in a metal sheet specimen used for truck exteriors, and it was divided into ten stages from (a) to (j).The peculiarity of this procedure consists of the incorporation of an initial condition based on points normalization to perform the matching using 2D-DIC.The general procedure is as follows: (j) Strain estimation from averaging the centroids' differences with the four neighbors using Equation (6).
The process includes the stamping of the known circle grid on the unformed metal sheet and the deformation of the sheet through the mechanical stamping process, because it is important to delimit the application of the proposed method.
The individual camera calibration can be performed with the "calib" function, while the stereoscopic calibration can be performed with the function "stereo_gui".Both are from the library "Camera Calibration Toolbox for Matlab" published by Jean-Yves Bouguet [40].The illumination can be performed through a blue LED set.It is important that the image capture is done with the calibrated camera system.

Four-Points Initialization
The construction of the 3D centroid mesh CM = C i,j of the deformed circles pattern begins with the process of obtaining the binarized image, which is detailed in Appendix A. After this, two clicks given on two consecutive ellipses of the left binarized image and two on the same circles of the right binarized image (Figure 5), allows for a knowledge of the distance and slope between two centroids of pseudo-ellipses.

Appl. Sci. 2019, 9, x FOR PEER REVIEW 7 of 17
The individual camera calibration can be performed with the "calib" function, while the stereoscopic calibration can be performed with the function "stereo_gui".Both are from the library "Camera Calibration Toolbox for Matlab" published by Jean-Yves Bouguet [40].The illumination can be performed through a blue LED set.It is important that the image capture is done with the calibrated camera system.

Four-Points Initialization
The construction of the 3D centroid mesh  =  , of the deformed circles pattern begins with the process of obtaining the binarized image, which is detailed in Appendix A. After this, two clicks given on two consecutive ellipses of the left binarized image and two on the same circles of the right binarized image (Figure 4), allows for a knowledge of the distance and slope between two centroids of pseudo-ellipses.
The clicks do not necessarily have to be done on the centroids; it is only necessary that they are done on the two consecutives pseudo-ellipses.The selected points are referred to as  , ,  , and ( , ,  , ), while in the right image the points are  , ,  , and ( , ,  , ).To work with discrete values, the values of each point, in each image, are rounded up because the computational functions give fractions of the pixels where the mouse is clicked via the computer.From the rounded coordinates provided by the user, the blob number  (considering a blob as a group of pixels with a connected connection equal to four) is retrieved to provide the coordinates of the centroid, also rounded; as an index of the tagged image (7).
where  can be the left or right side, and  (, ) is the labeled image as explained in Appendix A.
Next, the centroid of each blob  , is calculated using the statistical moments for 2D (8) with the help of Equations ( 9) and (10): (, ) = 1   (, ) =  0 ℎ (10) i (considering a blob as a group of pixels with a connected connection equal to four) is retrieved to provide the coordinates of the centroid, also rounded; as an index of the tagged image (7).
where s can be the left or right side, and L s (x, y) is the labeled image as explained in Appendix A.
Next, the centroid of each blob C s i,j is calculated using the statistical moments for 2D (8) with the help of Equations ( 9) and (10): x p y q BO s (x, y) The distance and the slope between the two centroids in Equations ( 11) and ( 12) are calculated to initiate an iterative automatic process of estimation of the 2D position of blobs.The estimation of the next blob position from the information of the two clicks on each image yc s i+1, j , yc s i,j , xc s i+1, j and xc s i,j : where g and h are the index of the distances and slopes matrixes.From Equations ( 11) and ( 12), the pseudo-row of marks is estimated using Equations ( 13) and ( 14): Meanwhile, to find the next row, Equations ( 15) and ( 16) are used: The normalized cross-correlation ( 5) is used to ensure a good pairing of centroids.This drastically decreases the number of operations performed for the matching between the interest points and at the same time takes advantage of the reliability granted by the NCC.The Matlab code can be downloaded from https://sites.google.com/itcelaya.edu.mx/dic/.

Results
The results of a typical initialization of the DIC are displayed in Figure 6 to show an initialization comparison.The working conditions were: one facet of 21 × 21, the search for the maximum normalized correlation, one-pixel step, and 100 marks that were pre-selected in the left image.
To observe the behavior of the normalized correlation internally, the correlation between a left image ellipse against all the possible 21 × 21 facets of the right image was plotted to qualitatively locate the maximum peak of the correlation (Figure 7).can be downloaded from Supplementary Materials.

Results
The results of a typical initialization of the DIC are displayed in Figure 5 to show an initialization comparison.The working conditions were: one facet of 21 × 21, the search for the maximum normalized correlation, one-pixel step, and 100 marks that were pre-selected in the left image.To observe the behavior of the normalized correlation internally, the correlation between a left image ellipse against all the possible 21 × 21 facets of the right image was plotted to qualitatively locate the maximum peak of the correlation (Figure 6).The engraving of light reference circles on the metal sheet resulted from an equidistance of 1.55 mm between the centroids when the metal sheet had not been deformed.The measurement made with a microscope is shown in Figure 7, where the before (Figure 7a) and after (Figure 7b) deformation is appreciated.On the other hand, the calibration parameters of the right and left cameras, as well as the translation vector between both cameras are the following: Intrinsic parameters of the left camera: The engraving of light reference circles on the metal sheet resulted from an equidistance of 1.55 mm between the centroids when the metal sheet had not been deformed.The measurement made with a microscope is shown in Figure 8, where the before (Figure 8a) and after (Figure 8b) deformation is appreciated.The engraving of light reference circles on the metal sheet resulted from an equidistance of 1.55 mm between the centroids when the metal sheet had not been deformed.The measurement made with a microscope is shown in Figure 7, where the before (Figure 7a) and after (Figure 7b) deformation is appreciated.The results of the image processing step are shown in Figure 9.In the top row, the pair of stereo images before the area filter is presented, while in the row below, the images after the area filter application are reported.The results of the image processing step are shown in Figure 8.In the top row, the pair of stereo images before the area filter is presented, while in the row below, the images after the area filter application are reported.

Left image
Right image Binarized image after pre-process of image.

Results
of filtering regions with less than 80 pixels areas.We selected one hundred points to show the method´s accuracy in a space of 100 × 100 marks.To verify if each point matching was right, the local DIC was calculated for the proposed points estimated in the initialization stage.The best correspondence was in the same position projected by the four-point method.Figure 9 shows a case of the correlation using a window of 31 × 31 pixels on a zone of 21 × 21 pixels.The MATLAB function used for this aim was "corr2".We selected one hundred points to show the method´s accuracy in a space of 100 × 100 marks.To verify if each point matching was right, the local DIC was calculated for the proposed points estimated in the initialization stage.The best correspondence was in the same position projected by the four-point method.Figure 10 shows a case of the correlation using a window of 31 × 31 pixels on a zone of 21 × 21 pixels.The MATLAB function used for this aim was "corr2".We selected one hundred points to show the method´s accuracy in a space of 100 × 100 marks.To verify if each point matching was right, the local DIC was calculated for the proposed points estimated in the initialization stage.The best correspondence was in the same position projected by the four-point method.Figure 9 shows a case of the correlation using a window of 31 × 31 pixels on a zone of 21 × 21 pixels.The MATLAB function used for this aim was "corr2".Figure 11 shows the points matching of an area of 10 × 10 deformed circular marks.Those are the same taken to perform the 3D reconstruction.On one side, the detail of the points matching is shown with a zoom.Figure 10 shows the points matching of an area of 10 × 10 deformed circular marks.Those are the same taken to perform the 3D reconstruction.On one side, the detail of the points matching is shown with a zoom.Figure 11 shows the three-dimensional reconstruction of the marks found and verified for the point matching stage.From the 3D reconstruction the strain was calculated using Equation (6), averaging the distances between the four neighboring marks of the mesh of circles of a radius of 0.5 mm and a separation between the centers of 2.55 mm.The strain estimated by the proposed method is shown in Figure 12. Figure 12 shows the three-dimensional reconstruction of the marks found and verified for the point matching stage.From the 3D reconstruction the strain was calculated using Equation (6), averaging the distances between the four neighboring marks of the mesh of circles of a radius of 0.5 mm and a separation between the centers of 2.55 mm.The strain estimated by the proposed method is shown in Figure 13.
Figure 11 shows the three-dimensional reconstruction of the marks found and verified for the point matching stage.From the 3D reconstruction the strain was calculated using Equation (6), averaging the distances between the four neighboring marks of the mesh of circles of a radius of 0.5 mm and a separation between the centers of 2.55 mm.The strain estimated by the proposed method is shown in Figure 12.The strains estimated in points having their four neighbors from a matrix of 100 × 100 were reduced to a matrix of 8 × 8.The values in Tables 2 and 3 were obtained with the stereoscopic vision system.From these results, it can be calculated that the average error in the experiment was, in this case,  = 1.290878 % (0.019995 mm.).This takes into account a distance between the centroids of the marks of 1.55 mm, without deformation.
Table 2.Estimated strain (%) using the proposed method.Table 3.Estimated strain (%) using a digital microscope.The strains estimated in points having their four neighbors from a matrix of 100 × 100 were reduced to a matrix of 8 × 8.The values in Tables 2 and 3 were obtained with the stereoscopic vision system.From these results, it can be calculated that the average error in the experiment was, in this case, = 1.290878% (0.019995 mm.).This takes into account a distance between the centroids of the marks of 1.55 mm, without deformation.

Discussion
The principal contribution of this work consists in the correct points-pairing between images (100%), in the specific case of almost-identical marks, to measure the deformation of a metal sheet with a small error.The method is semi-automatic due to the high complexity of the shown problem.This is evident when using the MATLAB function called "imregcorr" [41], and it is solved when the user indicates only two pairs of coinciding points.This works as an initial condition to begin a search of correspondent marks among the images.Basically, we noticed that the initial conditions of the search increased, in a notable way, the efficiency of the point-pairing system (with efficiency referring to the amount of paired marks, divided between the total of marks in the space of analysis).We assumed that the quality of the mark stamping, lighting, and the initial conditions, also contribute to trustworthy effects.According to the results, it can be observed that the correlation applied to the detection of deformed circular patterns (ellipses) is not ordered.This is because the objective of the construction of the pattern is to make the marks identical and that, in opposition to this, the correlation uses the differences between marks to emit a factor of similarity.The system was also compared with a manual deformation measurement using a microscope, which achieved a low sensitivity of the camera system and a displacement.This was due to the resolution of the camera system.In terms of processing time, it is well known that the correlation has a computational complexity of m 2 n 2 for an nxn image and an mxm template.Meanwhile, the proposed method depends on the number of points on the analyzed sheet; which is many times less than the number of pixels of an image, and it also depends on the distance between the camera and the printed sheet that is to be studied.

Conclusions
The proposed methodology solves the corresponding points-pairing for cases where multiple marks tend to be identical.This is one of the most complex and difficult problems DIC experiences in stereoscopic vision.The philosophical basis where most of points-pairing algorithms reside is the search for differences between the objects inside images.This goes against the goal of the Circle-grid of obtaining identical circles.The proposal in this work takes advantage of the characteristics of the problem to reduce its fallibility and to reach very important results, as well as the knowledge of an equal-points grid and a superficial deformation of less than 50%; which is of interest.The results are compared to the predetermined MATLAB function to show the advantages of this work´s proposal, which may be of use with a high degree of trustworthiness, thanks to the large amount of accurate correspondences.The processing time for the search in the initial conditions is relatively low because only distances, slopes and cosine functions are calculated.This way of solving the problem could be used in systems requiring patterns or points-grids, such as squares, hexagons and others.As a future work, we suggest to automatize in a total way the points´correspondence without losing the gained benefits.common to have small non-representative stains in that set of blobs  (, ), and this is why an interconnected pixels selection is applied with areas less than 80 pixels (A6):  (, ) =  (, ) ∀   (, )  (A6) The overall previous process is shown in Figure A1 and dot matrix position indexes Figure A2.

Figure 1 .
Figure 1.Global correlation with an error of image registration transform.

Figure 2 .
Figure 2. The surface view of a metal sheet from a metallographic microscope AX ZEISS @ 200X.

Figure 1 .
Figure 1.Global correlation with an error of image registration transform.

17 Figure 1 .
Figure 1.Global correlation with an error of image registration transform.

Figure 2 .
Figure 2. The surface view of a metal sheet from a metallographic microscope AX ZEISS @ 200X.

Figure 2 .
Figure 2. The surface view of a metal sheet from a metallographic microscope AX ZEISS @ 200X.

Figure 2 .
Figure 2. (a) Circles on metal sheet without deformation; and (b) circles on formed metal sheet.Figure 3. (a) Circles on metal sheet without deformation; and (b) circles on formed metal sheet.

Figure 3 .
Figure 2. (a) Circles on metal sheet without deformation; and (b) circles on formed metal sheet.Figure 3. (a) Circles on metal sheet without deformation; and (b) circles on formed metal sheet.

Figure 3 .
Figure 3. Setting of stereo cameras with illumination and metal sheet.
(a) Stamping of known circle grid on the unformed metal sheet.(b) Deformation of the metal-sheet through the mechanical stamping process.(c) Calibration of cameras.(d) Illumination of the piece with LED blue light for measuring.(e) Capture of stereo images.(f) Selection of four landmarks.(g) Search for neighbor's centroid.(h) Calculation of the NCC in the proposed neighborhood.(i) Triangulation of the points to obtain their position in 3D space.(j) Strain estimation from averaging the centroids' differences with the four neighbors using

Figure 4 .
Figure 4. Setting of stereo cameras with illumination and metal sheet.
(a) Stamping of known circle grid on the unformed metal sheet.(b) Deformation of the metal-sheet through the mechanical stamping process.(c) Calibration of cameras.(d) Illumination of the piece with LED blue light for measuring.(e) Capture of stereo images.(f) Selection of four landmarks.(g) Search for neighbor's centroid.(h) Calculation of the NCC in the proposed neighborhood.(i) Triangulation of the points to obtain their position in 3D space.

Figure 4 .
Figure 4.The user with blue (+) and the centroids of each blob in red (*).

Figure 5 .
Figure 5.The user with blue (+) and the centroids of each blob in red (*).The clicks do not necessarily have to be done on the centroids; it is only necessary that they are done on the two consecutives pseudo-ellipses.The selected points are referred to as x le f t 0,0 , y le f t 0,0 and x le f t 1,0 , y le f t 1,0 , while in the right image the points are x right 0,0 , y right 0,0 and x right 1,0 , y right 1,0 .To work with discrete values, the values of each point, in each image, are rounded up because the computational functions give fractions of the pixels where the mouse is clicked via the computer.From the rounded coordinates provided by the user, the blob number nb si (considering a blob as a group of pixels with a connected connection equal to four) is retrieved to provide the coordinates of the centroid, also rounded; as an index of the tagged image(7).

Figure 5 .
Figure 5. Pairing of points with typical initialization of the Digital Image Correlation.

Figure 6 .
Figure 6.Pairing of points with typical initialization of the Digital Image Correlation.Appl.Sci.2019, 9, x FOR PEER REVIEW 9 of 17

Figure 6 .
Figure 6.Correlation between a left image ellipse versus the overall possible 21x21 facets of the right image.

Figure 7 .
Figure 7. (a) Image captured before the forming mechanical process of the metal sheet; and (b) image captured after the forming mechanical process of the metal sheet.

Figure 7 .
Figure 7. Correlation between a left image ellipse versus the overall possible 21 × 21 facets of the right image.

Figure 6 .
Figure 6.Correlation between a left image ellipse versus the overall possible 21x21 facets of the right image.

Figure 7 .
Figure 7. (a) Image captured before the forming mechanical process of the metal sheet; and (b) image captured after the forming mechanical process of the metal sheet.

Figure 8 .
Figure 8.(a) Image captured before the forming mechanical process of the metal sheet; and (b) image captured after the forming mechanical process of the metal sheet.

Figure 8 .
Figure 8. Image segmentation to obtain referential ellipses for the deformation measurement.

Figure 9 .
Figure 9. Image segmentation to obtain referential ellipses for the deformation measurement.

Figure 8 .
Figure 8. Image segmentation to obtain referential ellipses for the deformation measurement.

Figure 9 .
Figure 9. Correlation on one matching projected by a four-point initialization using a window of 31 × 31 pixels on a zone of 21 × 21 pixels.

Figure 10 .
Figure 10.Correlation on one matching projected by a four-point initialization using a window of 31 × 31 pixels on a zone of 21 × 21 pixels.

Figure 10 .
Figure 10.Points matching using the initial condition proposed and verified by Normalized Cross Correlation.

Figure 11 .
Figure 11.Points matching using the initial condition proposed and verified by Normalized Cross Correlation.

Figure 11 .
Figure 11.Three-dimensional reconstruction of marks on the metal-sheet surface from point matching.

Figure 12 . 17 Figure 12 .
Figure 12.Three-dimensional reconstruction of marks on the metal-sheet surface from point matching.Appl.Sci.2019, 9, x FOR PEER REVIEW 12 of 17

Figure 13 .
Figure 13.Estimated strain in the selected zone.

17 Figure A1 .
Figure A1.Pair of images of the stereoscopic system showing the result of each stage of the preprocessing of images.The upper row presents the left-side images, and the lower row presents those from the right-side.

Figure A2 .
Figure A2.Binarized left-side photography to obtain a candidate's marks.

Figure A1 . 17 Figure A1 .
Figure A1.Pair of images of the stereoscopic system showing the result of each stage of the pre-processing of images.The upper row presents the left-side images, and the lower row presents those from the right-side.

Figure A2 .
Figure A2.Binarized left-side photography to obtain a candidate's marks.