# Design of an Automatic Defect Identification Method Based ECPT for Pneumatic Pressure Equipment

Authors: Bo Zhang; YuHua Cheng; Chun Yin; Xuegang Huang; Sara Dadras; Hadi Malek
Journal: Complexity
 (2018)

---

## Abstract

In this paper, in order to achieve automatic defect identification for pneumatic pressure equipment, an improved feature extraction algorithm eddy current pulsed thermography (ECPT) is presented. The presented feature extraction algorithm contains four elements: data block selection; variable step search; relation value classification; and between-class distance decision function. The data block selection and variable step search are integrated to decrease the redundant computations in the automatic defect identification. The goal of the classification and between-class distance calculation is to select the typical features of thermographic sequence. The main image information can be extracted by the method precisely and efficiently. Experimental results are provided to demonstrate the capabilities and benefits (i.e., reducing the processing time) of the proposed algorithm in automatic defect identification.


---

## Body

## 1. Introduction

As one of the most significant energy storage and transmission kind of equipment, the pneumatic pressure equipment has been widely used in modern aerospace industries, such as manned spacecraft and ground wind tunnel testing in Figure1. Moreover, given its wide range of applications, the equipment possesses a variety of specifications that make it flexible. The flexibility of their specifications adds complexity to problem of inspecting and testing this important equipment to identify defects. Researchers have studied methods to characterize pressure vessel properties in recent years. J. Lee et al. proposed a method of estimating the design pressure in a prismatic storage vessel for liquefied natural gas [1]. H. Al-Gahtaniet et al. studied the local pressure testing of spherical vessels with nozzles [2]. J. Proczka et al. reported guidelines for the efficient design and sizing of Small-Scale Compressed Air Energy Storage (SS-CAES) pressure vessels [3]. P. Blanc studied the residual burst pressure of composite cylinders after mechanical impacts on their cylindrical part [4]. These studies, however, do not explore the impact of defects related to the on-going use of pressure vessels made of diverse materials in a wide range of operating environments.Figure 1
Pneumatic pressure equipment.The complicated usage environments and diversity of materials cause the pneumatic pressure equipment to generate fatigue cracks, corrosion pits, and other defects. Therefore it is important to detect the potential defects for the pressure equipment. Nondestructive testing and evaluation (NDT&E) [5–10] are important to ensure the safety of manufacturing environments and equipment operation. K. Schabowicz et al. presented some nondestructive methods for the testing of concrete members and fibre cement boards [11]. Furthermore, the criteria for special specimen inspection are distinct, so various detection methods are needed including surface nondestructive testing, appearance inspection, and physicochemical test. The surface nondestructive testing mainly includes magnetic particle testing and penetrant testing. Due to the advantages of surface nondestructive testing, like high sensitivity and easy identification, this method been broadly adopted. However, there are many problems that need to be resolved in its practical application, like high labor intensity, low efficiency, pollution of the environment, poor safety, and so on. In recent years, ECPT has been developing rapidly. In addition to not damaging the body, the method is fast, is efficient, and can solve the issues of regular methods. So in this paper, the eddy current pulsed thermography is implemented in the detection of defects of pneumatic pressure equipment. Researchers also have made valuable contributions in processing data [12–15], which can be applied in (NDT&E). P. Zhu et al. presented an improved feature extraction algorithm for automatic defect identification based on eddy current pulsed thermography [16]. H. Yu utilized the normalized fuzzy weighting functions to dispose the records of experiment [17].Nevertheless, the above-mentioned researches possess low efficiency and limitations in processing data from low-alloy steel material. Besides, there are many methods, such as Independent Component Analysis (ICA); the processed feature information will be distorted after restructuring, because of data normalization, lending to no visual effect in detection. In recent years, adaptive technique has gained a lot of interest to enhance the efficiency of processing data [18–21]. To make the detection system be more useful for detecting the defect, it may own the ability of adapting property. A novel characteristic identification algorithm that uses the similarity of the typical transient thermal responses (TTRs) with the mixing vectors (the vectors of the pseudo-inverse matrix of the demixing matrix) in ICA is proposed. It solves the poor accuracy and low efficiency of the ICA efficiently by using the choice of a known message. The proposed algorithm has four steps. First the thermal imager collects the thermal image sequence of the pneumatic pressure equipment. Second, the TTRs are separated into parts using threshold values with a variable interval search. The variable interval search in an infrared image sequence can reduce the repetitive computation and hold the typical feature of the pneumatic pressure equipment. Third, the membership matrix and the clustering center are calculated to classify the acquired TTRs. Moreover, the rule of the largest distance of two class is used to identify the typical TTRs. Fourth, the typical TTRs are transformed into a two-dimension matrix in linear time. The main features of the infrared image sequence in the pneumatic pressure equipment can be extracted by the typical TTRs. Experimental results indicate the algorithm can select the typical feature more precise than ones of the ICA. Furthermore, the processing time of the presented algorithm is shorter than ones of the ICA. Therefore, the new algorithm can extract the main image information of the pneumatic pressure equipment accurately and efficiently in comparison to the ICA.


## 2. Background: Pneumatic Pressure Equipment and Nondestructive Testing

### 2.1. Pneumatic Pressure Equipment

The manufacture of high pressure vessels is complex, involving the coordination of subject knowledge and professional technology in many industries, including metallurgy, corrosion and protection, mechanical processing, safety protection, chemical engineering, and testing. The recent advances in technologies across various industries, including nondestructive testing, have enabled considerable progress in pressure vessel manufacturing.Pneumatic pressure equipment includes high, medium, and low pressure storage containers and inlet/outlet gas pipelines, as shown in Figure2. The pneumatic pressure equipment stores and conveys compressed air. The main features of these pressure vessels and the significance of their defect inspection are reflected in the following aspects.Figure 2
Pneumatic pressure equipment in China Aerodynamics Research & Development Center.( 1 ) Varieties and specifications: the types of material can be mainly divided into carbon steel, low-alloy steel, and a small amount of stainless steel. Moreover, the main specification of external diameter has dozens of kinds.( 2 ) Complex usage environment: pressure vessels and pipes often suffer from high pressure, some special aerodynamic erosion, and environmental corrosion in the open air, which usually lead to cracks, perforation, fatigue damage, and corrosion of the equipment.( 3 ) Devastating consequences of defects: accidents not only endanger the safety of personnel, equipment, and the plant, but also lead to delay of test task, failures, and so on.


### 2.2. Nondestructive Testing

Surface nondestructive testing is the key test method of pressure equipment, such as pressure piping and pressure vessels. Surface nondestructive testing basically includes magnetic particle testing and penetrant testing. The method possesses numerous advantages such as high defect detection rate and high sensitivity. But there are also many problems, including high labor intensity, long overhaul cycle, low efficiency, pollution of the environment, and poor safety. So there exist a great many obstructors for development of detection skill. In flammable explosive environment, because of the safety risks, the traditional types of surface nondestructive testing methods (magnetic particle, penetrant) cannot be implemented in the field.In recent years, surface detection by eddy current has been (developing) rapid. In addition to not damaging the body, the method is fast and efficient and can achieve large area rapid detection reducing manpower and material resources. The surface detection by eddy current is a method that utilizes thermal infrared image analysis.The temperature field distribution is based on a vortex phenomenon that can be measured using infrared detection. The device is a thermal infrared imager with high speed and high resolution. Moreover, through the analysis of the infrared thermal image sequence, we could detect the changes of electromagnetic and thermal characteristics of the structure and material defects. It can detect the surface and near surface by using the eddy current effect; it also can detect deeper defects. In the light of Joule’s law, the eddy current is turned into Joule heat in the testing piece generating high temperature and low temperature zones. The infrared thermal imager collects these temperature changes in the specimen. The collected data are analyzed to compute anticorrelation results and detect the presence of defects.Eddy density can be indicated asJe=σ×(∂A/∂t), in which A expresses the magnetic vector bit and σ indicates the electrical conductivity of materials. In the depth direction, the eddy current density is Je(z)=Je(0)·e-zπμσf, in which z is the depth, μ is the permeability of testing material, and f is the frequency of alternating current flowing through an exciting coil. The depth of eddy density has an attenuation to the 1/e of the strength surface, which is called as skin depth δ=1/πμσf. According to Joule’s law, the heat power is Pw=(1/σ)Je2=(1/σ)σE2 and E is electric field intensity. The Joule heat Q is transmitted inside the material, and the heat propagation equation is ρCp(∂T/∂t)-Δ(σT∇T)=Q, in which ρ is material density, Cp is specific heat capacity of materials, σT is thermal conductivity of materials, and T is temperature of materials. The crack depth is larger than eddy current skin depth in this specimen, so(1)Je=2βHcosh⁡2βx-cos⁡2βxcosh⁡βb+cos⁡βb, in which b is height, H is magnetic field intensity, and(2)β=πfσμ=1δ.The thermal power in the unit volume of the two sides of the crack isPw=0.5(1/σ)∫‍Je2dx, when the eddy skin depth is smaller than the plate thickness (δ≪b), the thermal power Pw=βH2/(2σ). Recognizing the defects change the heat distribution in the parts, or zones of a specimen, the principal component analysis enhances the zones where the thermal response and the eigenvector has similar trends. So, in the enhancement of defects in a principal component, the amplitude of other regions will be inhibited and distributed in a very small range. In this paper, a new algorithm in view of principal component selection is presented to automatically extract the principal components containing defects.


## 2.1. Pneumatic Pressure Equipment

The manufacture of high pressure vessels is complex, involving the coordination of subject knowledge and professional technology in many industries, including metallurgy, corrosion and protection, mechanical processing, safety protection, chemical engineering, and testing. The recent advances in technologies across various industries, including nondestructive testing, have enabled considerable progress in pressure vessel manufacturing.Pneumatic pressure equipment includes high, medium, and low pressure storage containers and inlet/outlet gas pipelines, as shown in Figure2. The pneumatic pressure equipment stores and conveys compressed air. The main features of these pressure vessels and the significance of their defect inspection are reflected in the following aspects.Figure 2
Pneumatic pressure equipment in China Aerodynamics Research & Development Center.( 1 ) Varieties and specifications: the types of material can be mainly divided into carbon steel, low-alloy steel, and a small amount of stainless steel. Moreover, the main specification of external diameter has dozens of kinds.( 2 ) Complex usage environment: pressure vessels and pipes often suffer from high pressure, some special aerodynamic erosion, and environmental corrosion in the open air, which usually lead to cracks, perforation, fatigue damage, and corrosion of the equipment.( 3 ) Devastating consequences of defects: accidents not only endanger the safety of personnel, equipment, and the plant, but also lead to delay of test task, failures, and so on.


## 2.2. Nondestructive Testing

Surface nondestructive testing is the key test method of pressure equipment, such as pressure piping and pressure vessels. Surface nondestructive testing basically includes magnetic particle testing and penetrant testing. The method possesses numerous advantages such as high defect detection rate and high sensitivity. But there are also many problems, including high labor intensity, long overhaul cycle, low efficiency, pollution of the environment, and poor safety. So there exist a great many obstructors for development of detection skill. In flammable explosive environment, because of the safety risks, the traditional types of surface nondestructive testing methods (magnetic particle, penetrant) cannot be implemented in the field.In recent years, surface detection by eddy current has been (developing) rapid. In addition to not damaging the body, the method is fast and efficient and can achieve large area rapid detection reducing manpower and material resources. The surface detection by eddy current is a method that utilizes thermal infrared image analysis.The temperature field distribution is based on a vortex phenomenon that can be measured using infrared detection. The device is a thermal infrared imager with high speed and high resolution. Moreover, through the analysis of the infrared thermal image sequence, we could detect the changes of electromagnetic and thermal characteristics of the structure and material defects. It can detect the surface and near surface by using the eddy current effect; it also can detect deeper defects. In the light of Joule’s law, the eddy current is turned into Joule heat in the testing piece generating high temperature and low temperature zones. The infrared thermal imager collects these temperature changes in the specimen. The collected data are analyzed to compute anticorrelation results and detect the presence of defects.Eddy density can be indicated asJe=σ×(∂A/∂t), in which A expresses the magnetic vector bit and σ indicates the electrical conductivity of materials. In the depth direction, the eddy current density is Je(z)=Je(0)·e-zπμσf, in which z is the depth, μ is the permeability of testing material, and f is the frequency of alternating current flowing through an exciting coil. The depth of eddy density has an attenuation to the 1/e of the strength surface, which is called as skin depth δ=1/πμσf. According to Joule’s law, the heat power is Pw=(1/σ)Je2=(1/σ)σE2 and E is electric field intensity. The Joule heat Q is transmitted inside the material, and the heat propagation equation is ρCp(∂T/∂t)-Δ(σT∇T)=Q, in which ρ is material density, Cp is specific heat capacity of materials, σT is thermal conductivity of materials, and T is temperature of materials. The crack depth is larger than eddy current skin depth in this specimen, so(1)Je=2βHcosh⁡2βx-cos⁡2βxcosh⁡βb+cos⁡βb, in which b is height, H is magnetic field intensity, and(2)β=πfσμ=1δ.The thermal power in the unit volume of the two sides of the crack isPw=0.5(1/σ)∫‍Je2dx, when the eddy skin depth is smaller than the plate thickness (δ≪b), the thermal power Pw=βH2/(2σ). Recognizing the defects change the heat distribution in the parts, or zones of a specimen, the principal component analysis enhances the zones where the thermal response and the eigenvector has similar trends. So, in the enhancement of defects in a principal component, the amplitude of other regions will be inhibited and distributed in a very small range. In this paper, a new algorithm in view of principal component selection is presented to automatically extract the principal components containing defects.


## 3. Introduction of Proposed Algorithm in ECPT

First of all, some definitions are introduced in the Appendix to make the algorithm clear. Then, the proposed algorithm is explicit as follows.Step 1.
3D matrixD saves the initial image sequence. Every pixels value of thermal images is saved in the matrix D. The third-dimension of D is the time axis t.
Step 1.1. Search for D(Rmm,Cmm,Tmm). To find the length of the interval in the column coordinate, this step sets (K=1,2,⋯k⋯) temperature thresholds T(m),(m=1,2,⋯K) from large to small; moreover the row which includes the highest temperature point D(Rmm,Cmm,Tmm) is divided into K+1 parts.
Step 1.2. Set temperature thresholds T_CLk, for each part, and the highest temperature point Dk(Rmm,Ck,:), (k=1,2,⋯,K) will be found in every part.
Step 1.3. Set P(P=1,2,⋯p⋯) temperature thresholds T(n),(n=1,2,⋯P) from large to small, at the same time, the column which includes the highest temperature point D(Rmm,Cmm,Tmm) is divided into P+1 parts. After the row is divided, the temperature value of the column also should be cut apart rationally in order to extract the characteristics of the defect more accurately.
Step 1.4. Set temperature thresholds T_RLP, for each part, and the highest temperature point Dp(Rmm,Cp,:), (p=1,2,⋯,P), will be found in every part. The number of each block is regarded as lenp, (p=1,2,⋯). Moreover, the same method as that of column is employed to find the row variable interval RLp. The specific procedure of searching for interval is in the following.
(a) Set(K=1,2,⋯k⋯) temperature thresholds to divide the TTR into K+1 parts.
(b) Set temperature thresholdsT_CLk for each part.
(c) Search for the highest temperature pointDk(Rmm,Ck,:).
(d) Calculate the PCC of(Dk(Rmm,Ck,:) with Dk(Rmm,c,:))
(e) until their PCC is less thanT_CLk.
(f) The number of(Dk(Rmm,Ck,:) whose PCCs with Dk(Rmm,c,:)) are more than T_CLk is regarded as CLk.
(g) Set(P=1,2,⋯p⋯) temperature thresholds to divide the TTR into P+1 parts.
(h) The same process as the searching for column interval to find the row interval.
(i) The row interval is regarded asRLp.
Step 1.5. What is more, set the threshold DD, initialize p=1,n=1,z=1, and carry out the following procedures:
(a) Setr=1, save the Dkp(r,c,:) into the X(z,:), r=r+RLp
(b) Ifr≤lenp, Compute the PCC of Dkp(r,c,:) with X(z,:), else go to step (d)
(c) IfPCC<DD, z=z+1, X(z,:)=Dkp(r,c,:), r=r+RLk, return step (b), else r=r+RLk, back to step (b)
(d) Ifr<M, p=p+1, z=z+1 back to step (a), else c=c+CLk
(e) Ifc<lenk, back to step (a); else go to step (f)
(f) Ifc<N, k=k+1, c=1, back to step (a), else the steps are finished
The specific calculation process is shown as Figure3.Figure 3
The specific process of Step1.Remark 1.
For the intervalCLk, the suitable method to seek out is searching for the length of region with largest temperature variation. What is more, the TTR with largest peak value is always contained in the region with the largest temperature variation. As analysis of Step 1, the length of region with largest temperature variation can be found by the coordinate value of TCL. For RLp, the homologous setting rule is analogous to CLk. The proposed algorithm can include the total typical temperature variations. And the representative temperature variations express the feature in the corresponding image sequence. characteristic pick-up algorithm is put forward to handle data. In the proposed algorithm, the appropriate interval values of column and row CIk, (k=1,2,⋯K), RIp, (p=1,2,⋯P) are set to reduce the repeated calculation. Moreover, the variable interval also can save the significant features. At the same time the proposed algorithm possesses less redundant computation than ICA.Step 2.
DivideX(:,z) into L parts. The specific calculation process is shown as Figure 4.
(a) Set the cluster number L. And, at the same time, initialize the cluster centermo; moreover set the number of iterations c, and set the weighting coefficient b; finally set terminating iterative threshold ε.
(b) Utilize formula to calculate membership functionuj(xi)=dji-2/(b-1)/∑s=1L‍dsi-2/(b-1), in which i=1,2,⋯,g, djk=Xk-mj, j=1,2,⋯,L; dji expresses the Euclidean distance between the ith information weight value and the jth cluster center. Xk is the information weight value. b is weighting coefficient. uj(xi) indicates the degree of the information weight value xi attached to the jth part.
(c) Update cluster centermj=∑i=1g‍[ujxi]bxi/∑i=1g‍[uj(xi)]b, in which mj presents the jth cluster center.
(d) Judge whether the absolute value of the objective function difference is smaller than the threshold value. IfJf(c)-Jf(c-1)≥ε, and i<g, i=i+1, go back to step (b). If i>g and j<L, j=j+1, return to step (b). If j≥L, stop. If Jf(c)-Jf(c-1)<ε, stop too.
(e) The membership maximization rule is used to deblur all characteristic information weight and get the category of each information weight.Mk=argmax(uj(xi)).Figure 4
The specific process of Step2.Remark 2.
Making use of the method of computing membership function and updating the clustering center in algorithm can extract feature information accurately and effectively. Taking advantage of the COV to classify data is called hard division. Hard division divides every object into a certain category strictly, and each class is unrelated to each other. However, the actual defect information objective exists intermediation in form and category. Moreover there is no definite boundary conditions to distinguish the class. Therefore, classifying the categories according to the membership and cluster centers of each class can be more accurate. The feature information weights in each category represent the optimal weight of the feature information, making the extracted features more accurate and more reliable.Step 3.
The mean value of theith classification is MEANi=(1/Ri)∑j=1,2,⋯,R(i)‍CLASji,i=1,2,⋯,L, in which R(i) presents the total number of TTRs in the ith classification. NPCCjti,(j=1,2,⋯,R(t)) is the noncorrelation value of MEANi and CLASjt, t=1,2,⋯,L and i=1,2,⋯,t-1,t+1,⋯,L (the jth  TTR does not belong to the ith classification). Let DISjti=NPCCjtip, where DISjti expresses the distance of two classes of TTRs  j in the tth classification with other classes. In tth (t=1,2,⋯,L) classification, FTt=DISjti∞, in which FTt is the final representation TTR of tth classification. Save FTt,(t=1,2,⋯,L) into Y (i.e., Y(:,t),t=1,2,⋯,L saves FTt,t=1,2,⋯,L). The specific calculation process is shown in Figure 5.Figure 5
The specific process of Step3.Remark 3.
The noncorrelation value ofMEANi and CLASjt, t=1,2,⋯,L is expressed as NPCCjti=1-PCC(MEANi,CLASjt).Remark 4.
∗ p is P-norm, ∗∞ is infinity-norm, and ∗:X→R satisfies the following:
( 1 ) positivity: x≥0, and x=0⇔x=0
( 2 ) positive homogeneity: cx=cx
( 3 ) trigonometric inequality: x+y≤x+y
x p = ( x 1 p + x 1 p + ⋯ + x 1 p ) 1 / p, and x∞=max(x1,x2,⋯,xn).Remark 5.
The purpose of the proposed algorithm in Step1 to Step 3 is to select the typical thermal responses. These typical responses have less intimately connection with each other in ECPT. The noncorrelation is larger and the thermal responses are more representative. Moreover, the selected accurate data is classified by Steps 2 and 3; the cluster center and membership function are carried out to classify the defect feature. The classification simplifies the data of image processing and makes the postprocessing more precise.Step 4.
The3d initial image sequence matrix D will be transformed into 2d matrix U. The elements in one row of U are taken columnwise from D(:,:,p),p=1,2,⋯,P. Calculate Y^ and solve this linear transformation S=Y^∗U, in which Y^ expresses the pseudo-inverse matrix of Y and R represents the result of the proposed algorithm. It includes the features of the initial image sequence processed by new algorithm. Moreover, continue to utilize the following way to achieve the extraction of defect feature O=(1/N)∑(x,y)∈S‍f(x,y) in which O expresses the mean of sum of the pixel values, f(x,y) is the pixel values, and N is the number of pixel values. Then we can obtain the uniform measurement of S; i.e., maxfx,y-m(x,y)∈S<V, in which V is a threshold. Finally, we can pick up the defect feature from the image.Remark 6.
The final purpose of the proposed algorithm is to extract the defect feature and ignore redundant information. The mean value of selected region is obtained by the similarity criterion of the method. Moreover the calculation of the uniform measurement divides the region into defect feature and other areas.Remark 7.
In order to highlight the efficiency and reasonability of the proposed method, the comparison of the ICA and the proposed algorithm is shown in Figure6. In the first part “Data whitening”, the data in the ICA is dealt with by normalization and then is computed to obtain the whitening vectors according to the corresponding feature vectors and feature matrix. These processes of calculation take a lot of time to dispose data and reduce efficiency of processing. However, the new algorithm does not have these computational procedures, to avoid the redundancy calculation. Hence, the amount of data in the proposed algorithm after processing of pretreatment and whitening is much less than that in the ICA, so, in the postprocessing, the proposed algorithm avoids redundancy calculation and possesses higher efficiency than the ICA. 2) In the second part “Data analysis and comparison”, compared with the extracted defect feature information graph, the amplitude of the graph extracted by the new algorithm contains the physical characteristics of the defect, and the content is closer to the defect feature information. Thus the proposed algorithm is more accurate and reasonable than the ICA.Figure 6
The comparison between the proposed method and ICA.Remark 8.
Step1 realizes the segmentation of data block and the calculation of column interval and row interval. Step 2 and Figure 3 show the process of the variable interval search. Step 3 is the correlation value classification. Step 4 selects the L typical TTRs based on the distance of two classes and extracts the typical features.


## 4. Experimental Design

The setup of the experiment is shown in Figure7. It includes five functional parts: induction heater; coil; PC; timing trigger; and IR camera. The induction heater produces high frequency alternating current for coil excitation. A rectangular coil applies directional excitation to heat the sample. The timing trigger controls the time for heating the sample. The IR camera records the thermal image sequence of the sample. The analysis of experimental data is represented below.Figure 7
The experimental schematic diagram.Figure8 shows the damaged pneumatic pressure equipment. All the test materials in this paper are provided by the China Aerodynamics Research & Development Center, which broadly supports research on aerodynamic equipment. In view of the devices volume and degree of damage, specimen 1 and specimen 2, shown in Figure 9, are extracted from the pneumatic pressure equipment to analyze.Figure 8
The damaged pneumatic pressure device.Figure 9
The specimen 1 and specimen 2.The defects have been detected using the method of ECPT and the algorithm presented in Section3. The two specimens are analyzed as follows. For sample 1, the parameter index of specimen 1 is presented in Table 1.Table 1
Parameter index of specimen 1.
SampleTime T e m p e r a t u r e T h r e s h o l d s THRE_CL REFR L 17s T ( 1 ) = 40 , T ( 2 ) = 80 0.98 R E F R 1 = 0.98 , R E F R 2 = 0.96 , R E F R 3 = 0.94 2The process of Steps3 and 4 is represented in Figure 10. In Step 3, set p=1, b=1, c=100, ε=10-5. There are 46, 13, 35, and 2 TTRs in the corresponding parts. Through the algorithm of Step 3, in the 1st part, compared with other TTRs, the 2nd possesses the maximum noncorrelation, DIS21i=1.3124. The same as in the 2nd part, the 11(th) possesses the maximum noncorrelation, DIS112i=1.3876. In the 3rd part, the 35th possesses the maximum noncorrelation, DIS353i=1.3097. In the 4th, the 1st possesses the maximum noncorrelation, DIS14i=1.3951. TTR1, TTR2, TTR3, and TTR4 have been extracted by Steps 3 and 4. The extraction results of proposed algorithm and the ICA result are shown in Figure 11. It is obvious that the features of the selected TTRs are similar to the ICs. In Figure 12, the difference of TTRs is shown. Comparing the proposed algorithm result with the result of ICA in Figure 13, the trends of the red curves are similar to the blue curves, respectively. Moreover, Figure 14 illustrates the precise and accurate defect feature processed by the algorithm. When comparing the proposed algorithm result with the result of the ICA in Figure 15, the final result of the algorithm ignores more redundant information than the ICA result. At the same time, the defect feature extracted by the proposed algorithm is more precise and accurate than the ICA result. Consequently, the algorithm has selected the typical thermal responses and extracted main features successfully. The proposed method can not only extract the main features like the ICA, but also reduce the processed time substantially. Its biggest advantage is the efficiency. Figure 16 shows the processing time of the ICA and the proposed algorithm. It is obvious that the proposed algorithm needs less time to complete the feature extraction process.Figure 10
Steps3 and 4 in proposed algorithm of sample 1.Figure 11
The proposed algorithm and ICA result of sample 1.Figure 12
The comparison ofTTR1, TTR2, TTR3, and TTR4 of sample 1.Figure 13
The normalized mixing vector 1, mixing vector 2, mixing vector 3, mixing vector 4, andTTR1, TTR2, TTR3, and TTR4.Figure 14
The defect feature of final result.Figure 15
The specimen defect, the ICA result, and the final result.Figure 16
The time comparison for specimen 1.For sample 2, the parameter index of specimen 2 is presented in Table2.Table 2
Parameter index of specimen 2.
SampleTime T e m p e r a t u r e T h r e s h o l d s THRE_CL REFR L 18s T ( 1 ) = 29 , T ( 2 ) = 40 0.97 R E F R 1 = 0.97 , R E F R 2 = 0.95 , R E F R 3 = 0.93 2The process of Steps3 and 4 is represented in Figure 17. In Step 3, set p=1, b=1, c=100, ε=10-5. There are 30, 30, 24, and 19 TTRs in the corresponding parts. Through the algorithm of Step 3, in the 1st part, compared with other TTRs, the 12th possesses the maximum noncorrelation, DIS121i=1.6693. The same as in the 2nd part, the 27th possesses the maximum noncorrelation, DIS272i=1.8238. In the 3rd part, the 6th possesses the maximum noncorrelation, DIS63i=1.3921. In the 4th, the 18th possesses the maximum noncorrelation, DIS184i=1.1842. TTR1, TTR2, TTR3, and TTR4 have been extracted by Step 3. The extraction results of proposed algorithm and the ICA result are shown in Figure 18. The Pearson correlation coefficients between mixing vectors 1, 2, 3, and 4 and TTR1, TTR2, TTR3, and TTR4 are 0.9786, 0.9435, 0.9678, and 0.9878, respectively. It is obvious that the features of the selected TTRs are similar to the ICs. In Figure 19, the difference of TTRs is shown. Comparing the proposed algorithm result with the result of ICA in Figure 20, the trends of the red curves are similar to the blue curves, respectively. Moreover, Figure 21 illustrates the precise and accurate defect feature processed by the algorithm. When comparing the proposed algorithm result with the result of the ICA in Figure 22, the final result of the algorithm ignores more redundant information than the ICA result. At the same time, the defect feature extracted by the proposed algorithm is more precise and accurate than the ICA result. Consequently, the algorithm has selected the typical thermal responses and extracted main features successfully. The proposed method can not only extract the main features like the ICA, but also reduce the processed time substantially. Its biggest advantage is the efficiency. Figure 23 shows the processing time of the ICA and the proposed algorithm. It is obvious that the proposed algorithm needs less time to complete the feature extraction process.Figure 17
Steps3 and 4 in proposed algorithm of sample 2.Figure 18
The proposed algorithm and ICA result of sample 2.Figure 19
The comparison ofTTR1, TTR2, TTR3, and TTR4 of sample 2.Figure 20
The normalized mixing vector 1, mixing vector 2, mixing vector 3, mixing vector 4, andTTR1, TTR2, TTR3, and TTR4.Figure 21
The defect feature of final result.Figure 22
The specimen defect, the ICA result, and the final result.Figure 23
The time comparison for specimen 2.


## 5. Conclusion and Future Work

In this paper, an accurate and more efficient algorithm in ECPT is proposed. The validity and efficiency of the proposed method are demonstrated with experimental results. The physical meaning of ECPT and the mathematical foundation of ICA are integrated in this proposed method. The results contribute to advancing the use of ECPT to detect defects as follows.( 1 ) The primary features of the thermal image sequences can be extracted by the proposed approach. Meanwhile the main features can be utilized to detect defects.( 2 ) Experimental results show that whitening preprocedure is time-consuming. The proposed algorithm has no data whitening procedure. So the proposed algorithm is more effective.( 3 ) After processing of pretreatment and whitening, the amount of data in the proposed algorithm is much less than that of in ICA. Hence, in the postprocessing, the proposed algorithm avoids redundancy calculation and possesses higher efficiency than the ICA. So the proposed algorithm is more accurate and reasonable than the ICA.Future work will pay more attention to how to enhance the efficiency and accuracy of the novel algorithm. In the experiment, some external interference will bring plenty of noise to influence the precision. Moreover, the change of threshold values in the algorithm is also anticipated to impact the efficiency of the detection. Therefore, it is worth researching how to reduce the external interference and to adjust the threshold value.



---

*Source: 5423924-2018-10-22.xml*