1. Introduction
Magnetic Induction Tomography (MIT) is a non-invasive and contactless imaging technique that is used to display the images describing passive electromagnetic properties (PEP), i.e., conductivity, permeability and permittivity [
1,
2]. The fundamental principle of MIT can be explained by the mutual inductance and eddy current theory: injecting an alternating current through excitation coils to generate a primary magnetic field, which interact with a conductive sample to induce a secondary magnetic field that can be detected by the sensor coil [
3,
4]. MIT has typically focused on imaging conductivity distributions of materials with high electrical conductivity by modelling the eddy currents in the forward model and then acquiring the conductivity distribution by solving the inverse problem [
5]. The inverse problem in MIT is a major challenge, but this has been solved conventionally by The Tikhonov [
3] and Total Variation (TV) [
6] regularization algorithms. However, previous works have demonstrated that the Tikhonov regularization method may produce overly smoothed reconstructed images with blurred boundaries between different materials. The quality of the reconstructed images can be improved by an enhanced inverse solver, such as the Total Variation regularization technique [
7,
8].
In recent years, MIT has been widely proposed as an appropriate imaging technique in both industrial and medical applications, such as brain stroke detection [
5], molten metal flow monitoring [
9,
10,
11], pipelines inspection [
12,
13], multi-phase flow imaging [
14] and non-destructive testing (NDT) for material characterization [
15]. It’s well known that steel production has been attracting an increasing level of interest for potential use in various applications that require a contactless method for remotely monitoring the casting process. Many previous works [
9,
10,
11] have sought to monitor the molten steel flow by reconstructing the conductivity distributions of the internal area of interest. However, MIT presents difficulties in detecting and imaging external surface defects and deformations due to its low spatial resolution. The inspection of the outer surface of pipelines using MIT was analyzed in [
12] through a traditional MIT pixel-based reconstruction method (PBRM) and a narrowband pass filtering method (NPFM). Metal samples in a cylindrical shape with different types of damage at the outer surface have been investigated as defected pipelines using an 8-channel MIT system, which demonstrated that PBRM can retrieve information with a detectability of 10%, while NPFM can achieve a resolution of 2%.
This paper focuses on detecting the defects or deformations on the external surface of metallic targets using a Total variation algorithm-based MIT. This novel algorithm is described in [
16]. Here, several simulation and experimental tests were carried out to investigate the ability of the novel algorithm to identify and image deformations and defects in metal. The proposed TV-based threshold-differencing algorithm was compared to the TV reconstruction method, and significant improvements were observed, which indicates that the novel algorithm for MIT may provide a suitable and valuable method for monitoring surface defects, deformations and displacements in hot steel continuous casting processes.
2. Methodology
The forward problem in MIT for surface defect application is based on the laws of induction and the eddy currents which are induced in the magnetic field with an alternating current [
3]. The governing equation can be written as Equation (1):
where
is the permeability,
is the angular frequency,
is the total magnetic vector potential as a result of the current source
and the effect of eddy current induced by the electrical conductivity
. The current density can be determined by the magnetic vector potential according to the Biot-Savart Law. Equation (1) is solved by approximating the system as a combination of linear equations in small elements with appropriate boundary conditions using the Galerkin’s approximation [
4,
7]:
where
is the linear combination of edge shape functions,
is the current source region (excitation coil),
is the entire region (current source region and eddy current region),
is the electric vector potential, which is defined as:
Then, after applying the volume integration equation, the induced voltage in the measuring coil can be calculated:
where
is the unit current density passing through the coil.
The Jacobian matrix
can be expressed by the relationship between the induced voltage in the sensing coil and the conductivity:
where
is the conductivity of pixel
,
is the volume of the perturbation,
is the forward solver of excitation coil
excited by
and
is the forward solver of sensor coil excited by unit current.
The inverse problem in MIT is defined as the retrieval of the unknown conductivity distributions
of targets from the measured voltage
, expressed by the linear equation:
where
,
,
F is the forward operator,
means the initial estimate voltage obtained from forward problem,
is the initial estimate conductivity,
is the Jacobian matrix obtained from forward problem. The resolution of the reconstructed images can be improved by increasing the size of Jacobian matrix.
The inverse problem can then be represented by solving the least-square problem:
The theory and application of both the Tikhinov and TV inverse solvers have been investigated and explained detailed in [
7]. The total variation method can be defined by adding a penalty term to Equation (7):
where
is the regularization parameter,
is the gradient and
is the
. The anisotropic version of the discrete TV functional [
6,
8] is adopted in this paper:
It should be mentioned that in this study, .
Then, the difficulty is to solve the constrained optimization problem:
The above constrained optimization problem can be controlled by applying the Split Bregman iteration, which is an iterative method based on Bregman distance [
8]. Then, the problem may be transferred to solve the Split Bregman equations [
7]:
where
is the
iteration and
is an auxiliary variable.
Based on [
16], the TV algorithm seeks to achieve the conductivity distribution of a perfect target
which retains as reference data. Consequently, the conductivity distribution of the target contain a deformation,
, which can be obtained using the TV method as well. Subtracting the reference conductivity distribution from the target conductivity distribution gives the deformed target conductivity as:
The image of the deformed target conductivity distribution was expected to indicate the location and the deformation information. And it’s obvious that the negative differential conductivity values represent the depression deformation part, and that the positive differential conductivity values represent the bulging. However, the reconstructed image of the deformed target is inadequate. For instance, the differential reconstructed image for a metal square with a depression deformation on one side may also show a bulging deformation on the opposite side due to the ill-posed nature of the inverse problem. This mirror effect can be tackled by comparing the absolute value of the maximum differential conductivity of deformed target |max{
}| to the absolute value of the minimum differential conductivity |min{
}|, and then eliminating the minimum to obtain the threshold-differential conductivity distribution, defined as [
16]
In fact, the quality and accuracy of deformed image can be further improved by forcing values below or above a certain threshold
to zero, such that:
where
. Introducing the
will remove some amount of the noise data to make the results more stable and precise. After applying different values of
during the research, we discovered that 0.3 ≤
γ ≤ 0.5 gave the optimal reconstructed images for the scenarios considered here. Appropriate values for
give threshold-differential images that display the position of deformation individually. The complete images with deformation of the target can be obtained simply by adding the reference conductivity distribution to the threshold-differential conductivity distribution, such as:
where
is the total conductivity distribution.
3. Simulation Results
In this section, we displayed the results in [
16], which investigated the ability of the TV-based threshold-differential reconstruction method by reconstructing a 74 × 74 mm
2 square sample with surface deformations using simulation data based on the laboratory prototype (
Figure 1A). The metal square sample was located at the center of a 16 coil array MIT system, as illustrated in
Figure 1B, with an operational frequency of
.
Figure 2 shows the conductivity distribution obtained for a perfect square metal without any deformation using the TV reconstruction algorithm, which indicates that the algorithm can produce reconstructed images representing the square metal with reasonable accuracy. Then, bulging and depression deformation were considered.
Figure 3 illustrates the simulation scenarios and the corrosponding reconstructed conductivity distribution images using TV algorithm for (A) a metal with a surface bulge on top, and (B) a metal with a surface depression on top, where the black lines represent the location of perfect square, the red line represents the boundaries of simulation setting of sample with bulging and the blue line represents the boundaries of sample with depression. The two scenarios were deformed by 9.25 mm at the midsection of one external surface. It’s obvious that although the bulge deformation in
Figure 3A can be detected, the image may also be interpreted to include a depression at the bottom, as shown in
Figure 3B.
Figure 4 illustrates the image obtained for the deformed target conductivity from applying Equation (14) for (A) the scenario with a bulge deformation on top, and (B) the scenario with a depression deformation on top, which therefore show only the effects of the deformation. It is obvious that a bulge (red) represented by positive numbers on the colorbar. A depression (blue) represented by negative numbers on the colorbar in both images can also be observed. However, the maximum and minimum values of
in the colorbar of the images shown in
Figure 4 give a clue about the actual deformation in the considered scenario. In
Figure 4A, the absolute value of the maximum deformed target conductivity, which represents a bulge, is greater than the absolute value of the minimum deformed target conductivity that represents a depression. Similarly,
Figure 4B shows that the the absolute value of the minimum deformed target conductivity, which represents a depression, is greater than the absolute value of the maximum deformed target conductivity that represents a bulge.
The threshold-differential conductivity distribution can be plotted from Equation (15), obtained for the scenarios with deformation to get a more accurate representation of the particular deformation.
Figure 5 shows the threshold-differential conductivity distribution obtained for (A) the metal with a bulge, and (B) the metal with a depression. Hence, a bulge represented by positive numbers on the colorbar in
Figure 5A and a depression represented by negative numbers on the colorbar in
Figure 5B can be observed, which indicates that we can obtain a better representation of the deformation in the scenario from
than from
.
By applying an appropriate value of
to the Equation (16), the threshold-differential conductivity distribution images illustrate the location of the deformation alone.
Figure 6 shows the threshold-differential conductivity distribution images obtained for (A) the metal with a bulge, and (B) the metal with a depression deformation when
. The threshold-differential conductivity distribution images located the deformation more clearly by applying
in
Figure 6 than in
Figure 5.
The TV-based threshold-differential image reconstruction algorithm with reference sample was analyzed for the outer shape reconstruction with surface deformations (bulging and surface depression defects), showing potential. Therefore, further research was conducted based on this method to analyze the image reconstruction of shapes with surface deformations. Supplementary simulation data were produced and analyzed: bulging or depression with different depth of the defect; bulging and depression at the same time on different surfaces and on the same surface.
Figure 7 shows the reconstructed image obtained for a square metal (
without any deformation using the TV reconstruction algorithm and larger Jacobian matrix (
). The resolution in this dissertation can be defined by a pixel resolution in pixel per
. The size of the square in the reconstructed image is
, which means that
; the resolution of this simulation scenario is
.
Samples with a bulging or depression defect on the top surface with a length of 74 mm (i.e., in the whole side) and with the 5 mm deformation depths were analyzed.
Table 1 indicates the threshold differential reconstructed images and their contour images obtained for the different forms of defects. The second column is the threshold differential image obtained for the sample with a bulging deformation on top with different defect depths. The fourth column is the threshold differential image obtained for the sample with depression deformation on top. The third and fifth columns are the corresponding contour images of reconstructed images. The blue line is the sample without any deformation, while the black lines are deformation parts. The threshold factor
was set to 0.5. The iteration number used for solving TV function equals 500, in reference to Equations (7)–(9), which can provide reconstructed images with sharp boundaries.
The deformation detectability can be defined as:
So, the deformation detectability of the proposed simulation system is
. In order to evaluate the accuracy of the proposed method applied to detect the surface defect, the normalized mean square error can be defined as [
17]:
where
is the reconstructed defect depth value and
is the actual one. For the case in which bulging and depression deformations occur at the same time:
By referring to these quantitative metrics and the above results, the conclusion can be drawn that the value of the defect depth can be detected and acquired accurately with an error ( and ) through the proposed method.
Table 2 illustrates the simulation scenario with two deformations at the same time (one bulging on the top surface and one depression at the bottom surface of the sample, both with a depth of 9.25 mm), as well as the differential reconstructed images.
The results show that the TV-based differencing algorithm has the ability to detect two different types of deformation at the same time. By referring to Equations (19) and (20), the quantitative metrics and the reconstructed defect depth can be acquired with an average error level of .
Further simulation work was conducted on bulging and depressions at the same surface, which was achieved by transferring the square sample (
) to a rhomboid one. The rhomboid samples can be defined by the parameter
.
Table 3 illustrates the simulation scenario and the results of the rhomboid problem with
and
.
The positive values on the color bar represent bulging as well as the negative value mean depression. The conductivity values also give a clue of the amount of deformation. The results shows that the TV-based differencing algorithm has the ability to detect two deformations on the same surface.
4. Experimental Tests
Regarding the inspection of pipelines, according to [
12], simillar results can be obtained when testing the same damaged cylinder samples using the proposed algorithms. In this section, three sets of MIT machine experimental tests were conducted to investigate and display the ability of the threshold-differential algorithm to utilise the Total Variation (TV) reconstruction method to image a 50 × 50 mm
2 metal cube with surface deformations located at the center of an 8-coil array, as illustrated in
Figure 8. The MIT system in this research consists of (i) a host computer (ii) an National Instrument data acquisation system, and (iii) an equally-spaced eight x 50 turns coils with 40 mm diameter circular-array forming a 110 mm diameter imaging-region. The working frequency is
. One metal cuboid represents the surface deformation. The demonstrated experiment system detectability is
. In this case:
Then, Equation (14) turns out to be:
Table 4 indicates the differential and threshold differential reconstructed images obtained for the sample with a bulging deformation at the center of different surfaces with different defect depths. The metal cuboid was placed at different positions which were close enough to the center of the metal cube surface. The two side lengths of the metal cuboid represent different defect depths (6 mm and 9 mm) according to different orientations. The threshold factor
was set to 0.5. It can be determined that
, so the resolution of this experimental scenario is
.
The reconstructed images obtained by differential and threshold differential TV algorithms both precisely indicate the position of the defect. The threshold differential reconstructed images and the corresponding contour images shown in Column 3 and Column 4 demonstrate better image performance than the differential images displayed in Column 2. When compared to the results of Scenario 1 and 2, the threshold differential image shown in Row 2 has a longer reconstructed length of the bulging defect than that of the image shown in Row 1, which is consistent with the real experimental scenario. So, the results displayed above prove that the TV-based differencing algorithm has the ability to locate and reconstruct the bulging defect.
A second stage threshold can be adopted to refine the accuracy of position locating and depth detection by applying a larger number of
to focus on the values around the maximum or minimum values and force the other values to zero.
was set at 0.8 in this scenario. The metal cuboid sample was located at a different corner of the metal cube with a different orientation to achieve different positions and variable defect depths (6 mm and 9 mm).
Table 5 displays the second stage threshold differential reconstructed images obtained for the sample with a bulging deformation in different corners with variable defect depths.
It’s obviously that the conductivity distribution images from the further threshold differential algorithm shown in Column 3 can locate and image the bulge more accurately than those obtained directly from the differential algorithm shown in Column 2. The defect depth can be acquired from the contour images displayed in Column 4 with an acceptable degree of error (, , and ).
Based on the simulation scenario with two deformation defects at the same time displayed in
Table 2, the corresponding experimental tests were implemented by the displacement of the square sample. In this case, Equation (14) turns out to be:
Table 6 indicates the differential reconstructed images obtained for the sample with bulging and depression deformations at same time, which were achieved by moving the metal cube sample and then applying the abstraction. The red color in the differential reconstructed images, formed by positive numbers in the colorbar, and blue color parts, produced by negative numbers in the colorbar, represent the bulging and depression defect, respectively.
The differential reconstructed images in Column 2 prove that the TV-based differencing algorithm can locate bulging and depression defects at the same time, as predicted from the simulation results shown in
Table 2. It can be seen from the colorbar that the bulging amount is approximately the same as the depression amount, as it should be. Meanwhile, the reconstructed depth can be acquired by the corresponding contour images shown in Column 3, with a low level of error (
and
). Moreover, the dispalcements of the sample can be detected by analyzing the differential resonstructed images.
Based on the simulation scenario of the rhomboid problem displayed in
Table 3, the corresponding experimental tests were implemented by rotating the square sample with 45°.
Table 7 indicates the differential reconstructed images obtained for the sample with bulging and depression deformations on the same surface, which were achieved by rotating the metal cube sample and then applying the subtraction.
The red and blue color parts in differential reconstructed images represent the bulging and depression defects respectively. The results indicates that the TV-based differencing algorithm has the potential to reconstruct the two deformations at the same surface, but further research should be conducted to obtain more accurate results.
5. Validation and Quantitative Evaluation of the Method
The TV-based threshold-differential algorithm in MIT was adopted into real factory hot test data to track the movement of the strand during a continuous casting process. In this case:
where is
is the conductivity distribution of the strand at one fixed time and
is any random time selected during the casting process.
Figure 9 shows the test set up for the 16 channel MTI system in the secondary stage of continuous casting. The system model is the same as coil geometry for
Figure 1.
Figure 10 illustrates the movement distance of the strand in X- and Y-directions with its corrosponding differential reconstructed images. The positive number in the colorbar form the red colour part of the diffrential conductivity distribution, which indicated bulging and also the moving direction of the strand.
Moreover, two factors were defined to evaluate and validate the performance of the proposed method for monitoring the movement of the strand during a continuous casting process:
Figure 11 displays the relationship between the sum value of bulging and depression with respect to the reference, and the movement distance of the strand to the reference point.
It is clear that the trend of these two plots is nearly the same, which indicates that the movement pattern of the strand is related to the sum conductivity distribution values of the differential reconstructed images. Therefore, the TV-based threshold-differential algorithm can be treated as a suitable method for monitoring surface defects, deformations and displacements.