Case Study of Empirical Beam Hardening Correction Methods for Dimensional X-ray Computed Tomography Using a Dedicated Multi-material Reference Standard

This paper presents a case study of two selected beam hardening correction methods and their effects on dimensional measurements of multi-material objects. The methods under test are empirical cupping correction (ECC) and empirical dual energy calibration (EDEC). These methods were originally developed for medical applications and their potential for the reduction of artefacts is typically only analysed based on grey value images. For testing and benchmarking of the mentioned methods for dimensional metrology, a dedicated multi-material reference standard—a multi-material hole cube—is used. This reference standard was originally developed for acceptance testing of CT systems. This paper shows a second application of this standard. The reference standard has been calibrated by tactile measurements to assess centre–centre distance errors as well as patch-based bidirectional length measurement errors on beam hardening corrected data and on uncorrected data. For the application of the method also to industrial multi-material scenarios, slight modifications of the ECC method are proposed. Practical aspects of both the ECC and the EDEC approaches as well as measurement results are analysed and discussed in detail. ECC was able to significantly improve dimensional measurements and was especially able to reduce extreme errors occurring in particular in multi-material scenarios by a factor of more than 4. EDEC, the dual-energy approach, reduced grey value inhomogeneities caused by artefacts even more. Its performance for dimensional measurements was however a little worse than ECC. EDEC data resulted in a slightly larger total range of residual measurement errors, mainly due to an elevated noise level.


Introduction
The use of X-ray computed tomography (XCT) devices as a coordinate measurement system (CMS) has made the complex task of dimensional measuring multi-material objects possible and also allows the measurement of assemblies in a mounted state. However, the reliable metrology of multimaterial components with XCT is a challenging task, since non-linearities in the projection data, caused by e.g. beam hardening and scattered radiation, lead to severe artefacts (i.e. imaging errors) in the reconstructed XCT images. These artefacts can severely hinder the performance of reliable dimensional measurements.
Several methods to correct the non-linear behaviour of multi-material components have been published. However, the potential for improvement is typically depicted only by grey-value-based image analyses and deducing quality metrics rather than evaluating the influence on material boundaries (surfaces) or errors of feature-based dimensional measurements [1]. This paper presents a comprehensive case study of the metrological performance of multi-material measurements of two selected beam hardening correction (BHC) approaches. Furthermore, this work attempts to perform a critical analysis of these methods from a metrological point of view by applying a calibrated multi-material reference standard for testing and benchmarking the different approaches. This reference standard-a multi-material hole cube (MM-HC) [2]-was originally developed for the acceptance testing of XCT-based CMSs. A second application of this standard is shown in this paper.

Beam Hardening Correction Methods for Multi-material Parts
The study uses Empirical Cupping Correction (ECC) described by Kachelrieß et al. [3] and Empirical Dual Energy Calibration (EDEC) described by Kachelrieß et al. [4,5] that correct non-linearities in the projection data due to the polychromatic nature of X-rays and partially scattered radiation. Contrary to other correction approaches for beam hardening artefacts [6][7][8][9], ECC and EDEC do not rely on iterative forward projections and filtered-back projections and they do not need any preliminary knowledge about the specimen (geometry, position and materials) and the spectra in use. This is-in principle-a massive advantage for the user of these methods and makes it worth quantifying whether the methods can be proven to provide an advantage for dimensional measurements also for the case of multi-material measurements. ECC and EDEC are practical methods with great potential for metrology too, but so far there have been no qualification studies on these applications. More complex methods with a high calculation burden that rely on iterative reconstructions [10][11][12] or on special data fusion techniques [13] are currently not practical for industrial applications and are therefore out of the scope of this work. The ECC method was originally designed to correct linearity issues for specimens made of one material only. With q 1 (D) being the log-attenuation values, the ECC method models the relation to the line of integration D by the correction function P(q 1 ). The index of q 1 is only of importance in the context of EDEC, the dual-energy approach, where an additional q 2 is used. In [3] and within this work, P(q 1 ) is limited to linear combinations of P n (q 1 ), which leads to corrected attenuation values p ECC approximated by a polynomial function of the order of N p ECC P(q 1 ) N n 0 c n P n (q 1 ) N n 0 c n q n 1 c 1 q 1 + c 2 q 2 1 + c 3 q 3 1 + · · · with c 0 0. (1) The linearity of the Radon transform R is used to calculate N basis images f n (r) R −1 q n 1 (2) that are combined to the reconstructed and corrected XCT image by a vector of coefficients c that minimise the deviation between f (r) and a given template image t(r), with r being the location in the corrected XCT image. This method assumes that the specimen is made of homogeneous material and, for example, does not have any kind of density gradient. In other words, the homogeneity of attenuation coefficients in the corrected XCT image is maximised and as a consequence, artefacts caused by beam hardening and scattered radiation are minimised. Although preliminary knowledge of the specimen is not required, it can be beneficial to use preliminary or expert knowledge for the generation of the template image (for further discussion see Sect. 3.1). The ECC optimisation procedure proposed by Kachelrieß where w(r) is a weight image that can be used to emphasise or exclude regions from the optimisation, e.g. exclude regions that are affected by partial volume effects [3]. Although ECC was initially proposed for medical applications, it can be applied to industrial parts made of one homogeneous material [1]. For industrial parts made of multi-material components, it is proposed to generalise the template t to . . . μ 2 for r ∈ material 2, μ 1 for r ∈ material 1, that contains attenuation coefficients μ for M different materials including air. Furthermore, the weight w(r) is set to zero in regions where the correspondence to a material is uncertain or if a region is not of interest.
w(r) > 0 for r ∈ materials 1 . . . M and air after erosion 0 for r ∈ regions with uncertain material and regions that are not of interest (6) Non-zero values in w(r) can be used to emphasise the importance of any material or air.
It is common for many of today's users of industrial XCT devices to use empirical correction functions that try to compensate for disturbing effects of beam hardening and scattered radiation. This even applies for multi-material parts, e.g. by using predefined correction functions for linearisation that are provided by manufacturers of XCT devices and are, in theory, able to compensate for non-linear effects in monomaterial scenarios. Results of Tan et al. [14] have shown that those predefined correction curves should be used carefully, since they can induce additional measurement errors. This work studies whether the proposed adaption of ECC to multi-materials is adequate for dimensional metrology.
When considering only the linearisation aspect, the basic idea of the second correction technique, EDEC by Kachelrieß et al. [4], is very similar to the ECC method. The main difference is that EDEC needs at least two full XCT scans with different X-ray spectra to perform the homogenisation, while the corrected attenuation values of Eq. (1) are instead calculated by with q 1 being the attenuation value of the low energy (LE) and q 2 that of the high energy (HE) projections. For practical and economic reasons, EDEC is used in this work as a dualenergy approach, although in general, it is a multi-energy approach. Furthermore, P ik (q 1 ,q 2 ) represents the (I + 1)(K + 2)/2 different basis images that are, in contrast to ECC, a polynomial composition of two attenuation values (LE and HE). Note that within this study, only the beam hardening correcting aspect of EDEC is of interest and not the extraction of the effective atomic number and the density of the inspected materials.

Description of the Multi-material Reference Standard Used
Properties of the ECC and EDEC methods are investigated using the MM-HC shown in Fig. 1. It is designed to analyse the multi-material influence on the length measurement error (E-test) in acceptance testing of XCT systems used for dimensional metrology [2]. The cube investigated within this work has a size of 30 mm × 30 mm × 30 mm and consists of two symmetric parts made of aluminium (EN AW 5083) and titanium (Ti6Al4V). This material combination represents a multi-material case where typical multi-material errors can be observed. Conversely, this combination does not reach the domain where the projection data is massively distorted as is the case in more challenging combinations (e.g. steel and plastics, tungsten and plastics, etc.). The cube features 17 holes used as measurement features as well as 12 "V"-shaped grooves, G1-G12, on the outside used to establish reproducible measurement positions, see Fig. 1b and Fig. 1c. Additionally, the cube features a step-like "cut" shape enabling different multi-material ratios along the standard's height, see Fig. 1d. Different ratios between the materials are present at grooves G2-G6, which allows in-and inter-material measurands in a multi-material scenario at a horizontal setting. Measurands in mono-material scenarios are given at levels G1 (100% Al) and G12 (100% Ti).

Experimental XCT Setup
XCT scans of the MM-HC are performed with PTB's Nikon MCT 225 system (using Nikon Metrology X-Inspect 3.1.9 acquisition software) with a voxel size of approximately (32 µm) 3 . It should be noted that the standard is scanned in horizontal orientation with a vertical rotation axis as shown in Fig. 1d to allow conclusions for different well-defined material ratios along the standard. Due to the circular scanning trajectory with a middle plane approximately at the height of groove G4, cone beam artefacts-especially at the top and bottom planar surfaces of the cube-would hinder measurements. Within this work however these plane surfaces are not of interest. Reconstructions are performed by filtered-back projection [15] without using the manufacturer's reconstruction software. Details on the evaluation procedure can be found in Sect. 2.4. The goal of this study is to compare measurement results achieved with uncorrected and beam hardening corrected scans of the MM-HC for a total given scan time of 200 min. The scan parameters of all XCT scans performed on the MM-HC are shown in Table 1, where ME refers to medium energy, LE to low energy and HE to high energy. Varying energy levels of the spectra are achieved by using different combinations of X-ray tube voltages as well as prefilter materials and prefilter thicknesses. The ME parameters represent a typical parameter selection by an experienced XCT operator for a single energy scan of the MM-HC. The ECC correction scheme is applied to the ME scan data. Furthermore, the LE and HE scans are used for the dual-energy-based beam hardening correction (EDEC). They have a combined total scan time of 200 min, i.e. the same time which is used for the single ME scan. The integration time of the HE scan is selected to be slightly longer to compensate for the intensity loss caused by the high amount of pre-filtering to harden the spectrum.
The X-ray tube voltages and prefilters of the LE and HE scan are selected to assure that they result in different artefacts in regions that should be represented by homogeneous grey values in the XCT images. The parameters are a trade-off between spectral separation and noise. X-ray tube power and thus the target current are selected with respect to the voxel size in order to reduce the geometrical blur of the focal spot.
Small residual scaling errors of the XCT scans are corrected by a multi-sphere scale correction to further improve the accuracy of the measurements additionally to the BHC methods and to allow an evaluation of the correction techniques free from other error sources. Each scale correction consists of a bracket of two scans of a multi-sphere reference standard. The scale correction is applied as the average Parameters mentioned at the bottom are constant of these two scans. Figure 2 shows flowcharts for the two series of measurements, one for the single energy ME scan and one for the dual-energy approach (LE + HE). For the ME and LE scan, the flat-field corrections (FFCs) have been performed at the beginning of both scan series. There have been no movements of axes during one series to allow a proper scaling correction. Scans for scale correction are done with scan parameters of the following or previous MM-HC scans, except that the number of projections is reduced to 2000 to speed up the process. The average of the scale correction values determined by the multi-sphere scans before and after scans of the MM-HC scans are used as correction values. The thermal expansion of the cube materials has been neglected since the temperature in the XCT cabinet was recorded to be (20.5 ± 0.3)°C. Thus, thermal expansion could lead to errors smaller than 0.3 µm only. Temperatureinduced focal spot movements due to different operational settings (ME, LE and HE) of the X-ray tube are avoided by appropriate stabilisation times of at least 1.5 h.

Evaluation Procedure
Dimensional length measurement error evaluations, i.e. deviation from the length reference values of tactile CMS The reference measurements are based on the measurements of all 17 holes at 7 heights indicated by the respective groove position. For each groove height, seven circumferen-tial lines at different heights (circumferential lines distance 25 µm) are measured inside the holes, see Fig. 1c. In total, about 500 000 tactile CMS points have been assessed. The tactile data serves as references for the XCT scans performed on PTB's Nikon MCT 225 system evaluated with a measurement strategy similar to the calibration.
For multi-material assemblies with a large difference in the X-rays attenuation, the standard procedure with a single threshold value for the method "advanced" surface determination of both materials does not provide acceptable results in the current version of the analysis software used (VG Studio Max 3.0.3). This is due to the fact that the starting value for the surface determination process does not consider two or more material peaks in the grey value histogram. When the peaks of the two materials are significantly apart from each other, as visible in Fig. 3a, there is a large offset for the initial threshold that is typically not compensated for by local threshold refinements of VG Studio Max 3.0.3. This is why an improved two-step surface determination approach was used for the evaluation of the MM-HC data. In the first step, the threshold value is customised for aluminium and in the second step for titanium, see Fig. 3. This approach is due to the large difference in the X-ray attenuation between Al and Ti. Slices of the resulting XCT volume (uncorrected ME) are shown in Fig. 3b and c, where the white contour represents the surface determined. The workflow of the analysis is described in [2] and below: The pre-processing approach is performed to remove wrong surface regions. A CAD model of one half of the MM-HC is registered to the surface of one material in the XCT volume to create a region of interest (ROI) in the XCT volume. Then a dilation process is applied to this ROI to ensure that the entire surface along the holes is not cut and a small region of air around the material remains in this ROI. This pre-processing approach is applied to Al and Ti, respectively. Although the CAD model of the HC is used to create a CAD-based ROI in the XCT data, it has no metrological influence other than reducing the time of the fitting process of the tactile points to the XCT surface.
The dimensional analysis of the correction methods for the MM-HC is based on centre-centre length measurements as well as on bidirectional inner and outer length measurements. Centre-centre measurements are one special example of unidirectional (averaged) measurements and are evaluated as the distances between centre points of cylinders which have each been created by least-square fitting of a short cylinder (with~175 µm height) at a specific height, see Fig. 4a. Bidirectional length measurements are evaluated as inner and outer length measurements. The distance is defined by two representative points created through dedicated patch operators, see Fig. 4b (outer measurement) and 4c (inner measurement). The patch operator used in this analysis decreases the noise of point measurements and allows a better analysis of the different beam hardening correction procedures. For the measurement analysed in this work, the patch operator covers a surface area of approx. 175 µm 2 and contains typically 300 calibration/measurement points.
For the centre-centre lengths, all possible distances between the holes are measured in the plane perpendicular to the hole axes. However, for the bidirectional measurements, at least six lengths (not independent) in seven main directions of the MM-HC are measured, see Fig. 4d. In total, 956 unidirectional centre-to-centre lengths and 2 × 378 bidirectional patch-based inner and outer lengths at different heights (i.e. different material ratios) were measured. The difference in the total number of analysed lengths is due to the calculation time which, for the patch-based analysis, is much longer than the one for centre-to-centre measurements. To provide a better view of the material influence depending on the material ratios, the average, the range and the standard deviation of length measurement errors per groove are evaluated for both centre-centre and bidirectional inner and outer length measurements.

Practical Aspects of the ECC and EDEC Approaches
In practice, it is useful to perform the optimisation procedure of ECC and EDEC only on some selected slices, instead of using the reconstructed volume of the complete specimen. This reduces the calculation effort of the optimisation drastically. In the presented case, seven slice images, located at the grooves G1 (Al mono-material layer), G2- G6 (Al/Ti multi-material layers) and G12 (Ti mono-material layer), respectively, are used for the determination of coefficients c. Therefore, also seven template t(r) and weighting images w(r) are created that define regions that should be homogenised. The generation of templates was done in three steps. First, global thresholding (ISO-50%) was used on uncorrected ME data for ECC templates and on HE data for EDEC. The thresholds from air to aluminium and air to titanium were calculated based on cumulated histograms that contained all seven slices. This initial binarisation led to over-and undersegmented regions where severe artefacts were present. Therefore, in the second step, the initial binarisations were manually corrected by using regional growing methods and manual filling in combination with morphological operations (erosion, dilation, closing, opening) to refine the binary images. The sequence and parameters of the applied methods were adapted to the local appearance of the artefacts. It is of great importance to create binarisations and subsequently a weighted template w(r)t(r) that excludes uncertain regions, where the correspondence to a material or air is unclear, e.g. voxels affected by partial volume effects. On the contrary, regions close to those voxels (e.g. at outer surfaces of the specimen) should be included to get rid of beam hardening where it is often most significant. Simple thresholding methods without visual checks and corrections would especially cause template errors in regions that are most relevant for homogenisation. As an example, Fig. 5 shows a set of three binarisations (M = 2: air, Al, Ti) for one ME slice image at groove G4. It can be seen from Fig. 5b that local region weights of 0 were used to exclude non-relevant air regions far from the outer surface of the MM-HC as well as four fixing screws made of polyether ether ketone (PEEK) (cf. Figure 4d) from the optimisation.
In the third step, we switched to linear indexing of the slice images, while consecutively using the three binary images (e.g. Fig. 5b, c and d) to mask zero values and concatenate the voxel values of all slices into one large one-dimensional vector f (r). Three images each filled with the material's template values are masked and concatenate similarly to generate the template vector t(r). Furthermore, three images each filled with 1/V M material are masked and concatenate to generate the weighting vector w(r) with V M being the total number of voxels that correspond to material M in the M binary images. The applied weighting makes Eq. (4) a minimisation of summed mean squared errors between template values and the corrected XCT image.
Within this work, only polynomial functions as basis functions for the optimisation were used. Another common basis function type are exponential functions, which at first seem to be best suitable due to the exponential terms in Beer-Lampert's law in general and for polychromatic X-rays in particular. Nevertheless, real X-ray acquisitions are affected by additional effects besides beam hardening, e.g. scattered radiation or geometrical blur, whereby polynomial functions typically result in smaller values of the optimisation metric E with the same number of coefficients in use (Eqs. (1) and (7)) for most applications and especially for the MM-HC.
Preliminary tests for the polynomial degrees have shown that the most appropriate degree is 3, whereby 3 leads to slightly smaller E values in comparison to a degree of 2, but to similar values when compared to 4 th order polynomials. Polynomials of degree 4 seem to over-parametrise the optimisation problem and are not advisable for the MM-HC on the XCT device in use in this study, since they lead to slight overcompensation of beam hardening artefacts especially at outer surfaces of the MM-HC. Conclusions were similar for the ECC and the EDEC correction techniques.  The minimisation problem of Eq. (4) within this work is solved by a Nelder-Mead Simplex approach [16], which is an iterative heuristic method that is not based on gradients and is practical for the numerical optimisation of multi-dimensional problems. Note that this method does not guarantee results to be optimal or perfect. Different start values for the iteration have been tested to initialise c for the ECC and EDEC optimisation. The coefficients for ECC (3 rd order polynomial) shown in Table 2 have been determined reproducible for various start values, without actual proof for optimality. In contrast, the EDEC minimisation ended at significantly different coefficients for most of the tested initial c values, which might be caused by the large number of coefficients. This means that there are multiple solutions for the chosen stop criteria of the iteration. Two termination criteria had to be satisfied to stop the iteration: (i) a lower bound on the step size for the coefficients of 10 −6 and a lower bound on the step size for the optimisation metric E of 10 −7 . The highest grey value changes in the EDEC image by varying initial coefficients could be seen along the two outer lateral edges of titanium with approximately ± 5% of titanium's template value. In internal regions where the cube has measurement features, the changes were below 1%. Therefore, this possible influence on dimensional measurements has not been further investigated for the EDEC procedure. Table 2 shows the initial coefficients and initial search steps used to determine the final coefficients. The ECC search started with p ECC q 1 and ended after 87 iterations. EDEC's search started at p EDEC= 0.5q 1 + 0.5q 2 and was terminated after 1389 iterations. The final 3rd order polynomials of ECC and EDEC that satisfied the above-stated termination criterion are further evaluated regarding metrology in the following Sect. 3.3.
It can be noted that the coefficients of the beam hardening correction "preset 2" (e.g. used in [2]) provided by the system manufacturer (which are not used in this study) are close to the ECC result for a 2 nd order polynomial.
The ECC and EDEC methods are implemented as MAT-LAB/Octave scripts with additional use of VG Studio Max 3.0.5 for the generation and analysis of the weighted template. Figure 6 shows slice images for the comparison of the uncorrected ME scan at the top, the ECC corrected ME scan in the middle and the EDEC corrected (LE + HE) scan at the bottom for the groove heights G2, G4 and G6 from left to right. Both BHC methods are able to homogenise regions that are strongly affected by artefacts caused by beam hardening and scattered radiation, at all groove heights, which means for all Al/Ti ratios, when compared to the uncorrected ME scan. In principle, the homogenisation is quantified by the optimisation metric, but the metric cannot be used to compare ECC and EDEC results. Therefore, contrast-to-noise ratio (CNR) values between air and aluminium as well as air and titanium are provided in Table 3. CNR values have been calculated using Eqs. (9) and (10) for all seven representative slices used for the optimisation. Template regions are utilised to calculate mean grey values μ and standard deviations σ . Note that  For ECC, q 1 are attenuation values of ME data. For EDEC, q 1 are attenuation values of LE and q 2 of HE data Template regions define regions of air, Al and Ti. The calculation was performed on a common volume of all seven slices σ Al and σ Ti were used in the denominator to quantify the homogeneity and noise of the two materials [17]. All seven slices and template regions are treated as one volume.

Image-Based Evaluation
CNR Al,air μ Al − μ air σ Al (9) CNR Ti,air μ Ti − μ air σ Ti (10) Table 3 depicts an increase in CNR from uncorrected to ECC corrected data, whereby CNR Ti,air might be slightly higher than CNR Al,air since, during the ECC optimisation, the squared errors for the optimisation metric E (see Eq. (4)) in titanium regions are larger and thus more relevant than those of aluminium. This might put more weight on the optimisation of titanium. Hence, ECC reduced artefacts caused by any kind of non-linear behaviour, which is confirmed by CNR values in Table 3 and images of Fig. 6. CNR values of EDEC corrected data are overall the lowest (i.e. worst). This is due to slightly smaller contrast terms in the enumerator of EDEC data in combination with denominators σ Ti and σ Al that are at least twice as high as standard deviations of uncorrected and ECC data. Furthermore, visual compar- Fig. 6 Example of slice image of top row the uncorrected ME scan, middle row ECC corrected ME scan and bottom row EDEC corrected LE plus HE scan for the groove heights left column G2, middle column G4 and right column G6. Slices of uncorrected ME and ECC corrected ME were generated using a grey value windowing of [− 0.1, 0.09] and the EDEC images with [− 0.1, 0.075] isons in Fig. 6 depict that EDEC was able to compensate for artefacts even more than ECC, but is prone to noise and ring artefacts, which severely affected the CNR values. Possible influences of these findings on material transitions and thus dimensional measurements are discussed in the next section.

Metrology Results
Computed tomography datasets of the multi-material hole cube standards were obtained with the PTB's Nikon MCT 225 with different measurement parameters in order to compare the metrological performance of the selected BHC methods. An ME scan representing a typical set of scan parameters used by an experienced XCT user as well as an LE scan and an HE scan providing enough difference in the artefacts in the XCT image to apply the dual-energy-based EDEC correction were carried out.
Unidirectional (averaged) centre-centre measurements as well as bidirectional inner and outer measurements (cf.

Al-Ti Ti
Increasing amount of Ti Fig. 7 Boxplots of performed bidirectional inner and outer measurements as well as centre-centre measurements (cf. Figure 4) for left column uncorrected, middle column ECC and right column EDEC data; voxel size ≈ 31.738 µm age, standard deviation and maximum and minimum values per groove are given.
In general, a significant improvement of the results for both BHC methods could be observed, whereby the ECC method seems to give somewhat better results than the EDEC approach. Improvements are noticeable as smaller mean, minimum and maximum values at all grooves, but are especially strong at grooves G5 and G6 featuring a high Ti percentage.
For inner bidirectional measurements, a total error range of 75 µm (i.e. approx. 2.3 × voxel size) was observed for the uncorrected data, see For outer bidirectional measurements, a total range of 98 µm (3.1× voxel size) was observed for the uncorrected data, see Fig. 7d. Again, better results were observed with the ECC methods, the total range of measurement errors was around 15 µm (i.e. 0.4 × voxel size), while for the EDEC approach, the total range was bigger with 38 µm (i.e. 1.2 × voxel size), see Fig. 7e and f.
An even better improvement of the measurement results was also observed for the centre-centre unidirectional measurements. The uncorrected data presented a 39 µm (i.e. 1.2 × voxel size) range of the measurement error, see Fig. 7g. On the contrary, ECC decreased the total range of the measurement error to 9 µm (i.e. 0.3 × voxel size), see Fig. 7h. Applying the EDEC method, the total range of the measurement range to 17 µm (0.5× voxel size) was reduced, see Fig. 7i. However, only the extreme cases were successfully corrected by EDEC (i.e. at G5 and G6) and the measurement results for some of the additional grooves were slightly degraded using the EDEC method, although the remaining scaling errors seem to be very low.
For the centre-centre measurements, a detailed analysis has been performed on groove 5, since it shows the largest maximum and minimum values in the box plots of Fig. 7. The detailed analysis in Fig. 8 is based on a colour plot of all centre-centre measurements carried out at G5, where the colour code represents the measurement error in µm. The vertical and horizontal axes of the plot indicate the hole where the length measurand starts and ends, respectively. This plot gives a quick overview of all measurements at G5 and allows the identification of holes where the highest deviations occur. The two holes H10 and H11 represent measurements at groove G5 height where a single hole consists of two materials. For this reason and due to a potential non-concentricity of the two halves, the holes H10 and H11 have been excluded from the analysis in Fig. 8. With this colour plot, it becomes clear that the lengths contributed by the two holes H12 and H14 represent the highest deviations for the uncorrected ME scan data. A significant improvement of the results is noticeable in the colour plot of the ECC and EDEC methods. Figure 8 further shows slice images at G5 and grey value profiles through hole H12 at height G5. Additionally, the CNR air,Al was calculated in the H12 region based on Eq. (11). The analysis reveals that holes H12 and H14 in the aluminium part at groove G5 are strongly influenced by artefacts due to the long penetration length containing titanium and the nearby presence of titanium. Similarly, holes 15, 16 and 17 in aluminium are affected at G6 and cause high measurement errors (see Fig. 7a, d and g).

CNR air,Al
μ Al − μ air σ air (11) CNR analysis in Table 3 and the grey value profiles in Fig. 8g indicate that measurements using EDEC data might suffer from an elevated noise level and slightly more ring artefacts, when compared to ECC. This might be the reason why ECC performs slightly better than EDEC especially for the total range of measurement errors. The mean errors are comparable between the ECC and EDEC approaches, which indicates that the material influence is reduced for both BHC methods.
Although, minimising temperature-related geometrical misalignments and scaling errors between the LE and HE scans was tried, some influences were still noticed. Prior to applying the scaling correction, the relative stability between the scaling correction scans before and after the actual scans was less than 0.3 µm for the ME scan series and 3.0 µm for the LE + HE scan series over a length of 40 mm. This indicates that the LE and HE scans had slightly different magnifications. One additional cause worth noting is the Xray spectra-dependent magnifications that have been reported by Illemann et al. [18].

Conclusions
The two beam hardening correction techniques, EDEC and ECC, were successfully applied to an industrial multimaterial scenario, with materials that show a significant difference in their absorption properties. The test scenario was created using a calibrated multi-material hole cube made of aluminium and titanium originally designed for developing future multi-material acceptance testing. For ECC, originally designed for mono-material cases, modifications were proposed and implemented to enable the method to also be used for multi-material scenarios. Even though the procedures under study do not need any preliminary knowledge of the specimen (geometry, position and materials), a manual refinement of the template images for both procedures was necessary. This procedure, mainly working with morphological operations, ensures optimised conditions for the reduction of artefacts caused by beam hardening and partially scattered radiation, leading to more reliable measurements by minimising measurement errors. Table 4 shows a summary of the achieved improvements in terms of total ranges of measurement errors for all measurements on uncorrected, ECC and EDEC corrected data.
Both BHC approaches show significant improvements for the dimensional measurands under study (centre-centre length measurement errors and as well bidirectional inner and outer length measurement errors), especially for the highest errors which occur at low absorbing aluminium areas in mea-surement slices with a high content of the strong absorbing titanium. The improvement in the total range of measurement errors reaches a factor between 2.1 (for the EDEC correction of bidirectional inner lengths) and 4.3 (for the ECC correction of centre-centre lengths).
Only on rare occasions-slices with 100% titanium and analysing centre-centre lengths only-does the EDEC correction result in larger outliers compared to the uncorrected state mainly due to noise, even though the remaining inhomogeneities related to artefacts and mean measurement errors are smaller compared to ECC.
For the experimental setup used to scan the multi-material hole cube, the use of an additional spectrum does not seem to be beneficial for a constant total scan time of 200 min. ECC made more efficient use of the total exposure time and led to even better corrections. It is also easier for operators of XCT devices to handle on a daily basis, since no special care is needed to assure minimal misalignments between the LE and the HE scans and only one FFC needs to be assessed for the ME scan. It has to be mentioned that the potential for reducing beam hardening artefacts with ECC potentially worsens, the more different the attenuation behaviour of scanned materials is, since this linearisation method is actually only capable of dealing with mono-material scenarios. ECC may not be able to handle even larger absorption differences properly and thus EDEC is expected to outperform ECC in more challenging cases with extreme variations in the attenuation behaviour of investigated materials (steel and plastics, tungsten and plastics etc.).
All in all, both correction techniques may help to improve the measurement capabilities of recent XCT systems when measuring multi-materials on a level that today's user expects from mono-material measurements. At the moment both applied procedures require manual adjustment-a step which the authors believe could be automated using rule-based refinement procedures and preliminary knowledge.