Mineral Characterization in Human Body: A Dual Energy Approach

: Kidney and uteric stones are a common cause of pain and disturbance in numerous people worldwide, as they tend to reappear. There are several studies investigating the association of urolithiasis and nephrolithiasis with atherosclerosis, as patients suffering from the former diseases were found to have been diagnosed with atherosclerotic plaques. Hydroxyapatite and calcium oxalate are the most common mineral crystals found in both kidney/uteric stones and atherosclerotic plaques’ calciﬁcations. Even though for stones smaller than 5 mm surgery is not recommended, the knowledge of the stone composition is an important tool for the physician in order to provide better treatment for the patient. The mineral crystal characterization of atherosclerotic plaques’ calciﬁcations smaller than 3 mm (spotty calciﬁcations) will assist the physician to limit the possibility of myocardial infraction and stroke, as the presence of hydroxyapatite indicates possible plaque rapture. To this aim, a dual energy (DE) X-ray method was developed in this work. The calcium/phosphorus mass ratio ( m Ca / m P ) was determined through analytical simulations and the results were veriﬁed experimentally. Both monoenergetic and polyenergetic simulation studies were implemented for hydroxyapatite, calcium carbonate and calcium oxalate with thicknesses ranging from 0.50 to 3.00 mm, at 100 µ m increments, to obtain the optimized irradiation conditions. The experimental veriﬁcation of the proposed method was performed using an X-ray tube combined with a high resolution complementary metal-oxide-semiconductor (CMOS) active pixel sensor (APS) detector. The Mann– Whitney U test indicated that statistically signiﬁcant differences were found between the different types of minerals examined for thicknesses of 0.70 mm or higher.


Introduction
Urolithiasis and nephrolithiasis are frequent causes of intense pain affecting higher percentages of the population and increasing the expenses for health systems [1][2][3]. Based on previous works, the majority of ureteric and kidney stones will pass spontaneously, especially when they are smaller than 5 mm [4][5][6][7][8]. In this case, no external assistance is needed. On the other hand, when a stone cannot pass spontaneously, elective surgery will be applied [3]. It is not only the presence of stones but also the repetition of their presence in the patients that urge the physicians in deeper learning of stones' pathophysiology through investigation over their composition [9,10]. The known stone composition will lead to patients' better treatment since the appropriate strategy will be followed for each patient and unneeded operation may be prevented [3,[9][10][11][12][13][14]. The earlier the stone is found, the lower the probability of an operation and reappearance [10,[15][16][17]. Calcium oxalate minerals are most commonly found in kidney/ureteric stones, followed by hydroxyapatite [9,10,[18][19][20][21][22][23][24][25][26][27][28]. Calcium carbonate can also be found in stones; however, its presence is rare enough [23,29,30]. Computed Tomography (CT) is the gold standard method for calcified plaques were more accurately quantified with their technique [68]. Henzler et al., by using human specimens, explored the feasibility of atherosclerotic plaques' characterization on a density basis. The results of this study were compared to histopathology findings. Their results indicated that they could discriminate atheroma, fibroatheroma and calcific plaques [69]. Nakajima et al. proposed the use of the effective atomic number as an indicator for differentiation between soft and fibrous non-calcified coronary plaques. Compared to the Hounsfield unit-based differentiation obtained when CT is used, effective atomic number-based differentiation with the use of DECT resulted in more accurate classification of the plaques [70]. The effective atomic number was also used by Nishizawa et al. in an in vivo study in an attempt to analyze the calcific lesions' composition. They concluded that although hydroxyapatite is mostly found in plaques' calcifications, calcium oxalate was found in the majority of atherosclerotic plaques of dialysis patients [63]. Ding et al. conducted a simulation study to investigate if three and four material decomposition with the use of DECT could lead to composition analysis of atherosclerotic plaques. Three and four low atomic number materials were used in three and four material decomposition, respectively. The results showed that the proposed technique has the potential to quantify the atherosclerotic plaques' composition [71].
Even though kidney stones are not fatal, they may provoke cardiovascular diseases that are among the diseases with the highest mortality rates [72,73]. Several studies investigated the association between urolithiasis/nephrolithiasis and cardiovascular diseases including atherosclerotic plaques [72][73][74][75][76][77][78]. Aydin et al. investigated data from two hundred patients with kidney stones composed of calcium oxalate and two hundred controls to estimate the ten-year risk of cardiovascular disease and mortality. The results indicated higher risk for both cardiovascular disease and mortality for the patients that had a history of kidney stones compared to controls [75]. Davarajan published a review paper presenting the studies correlating kidney stones to cardiovascular diseases, and especially atherosclerosis. This review paper clarifies how kidney stones lead to atherosclerosis and atherosclerosis to kidney stones. Based on all the preclinical and clinical studies included in this research, the association between these two diseases is a two-way process [73]. Huang et al. worked on the assumption that the presence of atherosclerotic plaques is related to the existence of kidney stones in a patient, and especially stones composed of calcium. They studied the relation of serum total cholesterol, high-density lipoprotein cholesterol, and low-density lipoprotein cholesterol on twenty-four-hour urine chemical components and calcium composed stones. The evaluation of the association between calcium composed stones, serum cholesterol levels and the extent of atherosclerotic plaques was accomplished by comparing the results with the controls. The results indicated that there was significant correlation between calcium composed kidney stones and atherosclerotic plaques [76]. Luo et al. in their review paper collected and presented the results of numerous studies investigating the association between urolithiasis and atherosclerotic plaques. More than 70,000 samples were examined in all these studies. The assessment of all the studies' results revealed that the possibility of atherosclerotic plaque formation when ureteric stones were found in the samples was higher [79]. Although the association between ureteric/kidney stones and atherosclerotic plaques requires more investigation, the screening for cardiovascular diseases in patients with ureteric/kidney stones shall be taken into consideration [77,78].
Our research team has worked on the discrimination between the different mineral crystals found in breast cancer [79,80]. Calcium oxalate is the mineral crystal found in Type I breast calcifications and is associated with benignancy [81,82]. Type II malignant breast calcifications are composed of hydroxyapatite, while the presence of calcium carbonate in Type II indicates benignancy [83][84][85]. The differentiation was accomplished through the estimation of the effective calcium/phosphorus mass ratio (m Ca /m P ) with the dual energy (DE) X-ray method that we developed. Since phosphorus is present only in hydroxyapatite, an effective m Ca /m P was estimated for calcium carbonate and calcium oxalate [79]. The method was experimentally evaluated with the use of photon counting energy dispersive and energy integrating detectors [80,86,87]. There was a statistically significant difference between the minerals examined for calcifications 300 µm thick or higher.
In the present study, a dual energy X-ray method was developed in order to discriminate between calcium composed kidney/ureteric stones. The proposed method was also used to differentiate the calcium composed atherosclerotic plaques. In both cases, m Ca /m P was estimated. Simulation monoenergetic and polyenergetic studies were conducted based on a previously validated modified analytical model to obtain the optimized irradiation conditions. The method was experimentally verified with the use of an X-ray tube combined with a high resolution complementary metal-oxide-semiconductor (CMOS) active pixel sensor (APS) detector. The results of the study indicate that the differentiation of hydroxyapatite, calcium carbonate and calcium oxalate is feasible when m Ca /m P is used as a differentiation parameter.

Monoenergetic Beams
The analytical model of the study calculated the attenuated intensities for the low-/high-energy beams as follows. Assume that along the monoenergetic X-ray path a torso of thickness T consists of: (i) surrounding tissue t s and (ii) kidney/uteric stone or calcification in atherosclerotic plaque t c (Figure 1). The surrounding tissue was simulated using attenuation coefficients of water. Since the phantoms of kidney/uteric stones and calcifications in atherosclerotic plaque were composed of a mixture of mineral and epoxy resin, the same quantity of epoxy resin was also used in the surrounding tissue. Thus, when only surrounding tissue was in the path of the X-ray beam, the low-/high-energy intensities, I s,El and I s,Eh , were calculated as: (1) For both kidney/uteric stones and calcifications in atherosclerotic plaque, three different types of mineral crystals were investigated: (i) hydroxyapatite (Ca10(PO4)6(OH)2), HAp, with a density of 3.18 g cm -3 [88], (ii) calcium carbonate (CaCO3) with a density of 2.93 g cm -3 [89], and (iii) calcium oxalate (CaC2O4) with a density of 2.20 g cm −3 [90]. The epoxy resin (C19H23CIO4) used in the simulation study had a density of 1.06 g cm −3 (Resoltech, Rousset, France) [91]. The mineral crystals are referred to as HAp, CaCO3 and CaC2O4.
The unattenuated intensities, Io,El and Io,Eh, for the low and high energy, respectively, were calculated through the entrance surface dose (Ka) (Equation (3)) for all photon ener- When both surrounding tissue and kidney/uteric stone or calcification in atherosclerotic plaque were in the path of the X-ray beam, the low-/high-energy intensities, I c,El and I c,Eh , were calculated as: where I o,El and I o,Eh are the unattenuated intensities for the low and high energy, respectively. The energy-dependent linear attenuation coefficients (cm −1 ) and thickness (cm) for the surrounding tissue and the mineral crystals were given by µ s,Ei , µ mineral,Ei and t s , t c , respectively. For both kidney/uteric stones and calcifications in atherosclerotic plaque, three different types of mineral crystals were investigated: (i) hydroxyapatite (Ca 10 (PO 4 ) 6 (OH) 2 ), HAp, with a density of 3.18 g cm -3 [88], (ii) calcium carbonate (CaCO 3 ) with a density of 2.93 g cm -3 [89], and (iii) calcium oxalate (CaC 2 O 4 ) with a density of 2.20 g cm −3 [90]. The epoxy resin (C 19 H 23 CIO 4 ) used in the simulation study had a density of 1.06 g cm −3 (Resoltech, Rousset, France) [91]. The mineral crystals are referred to as HAp, CaCO 3 and CaC 2 O 4 .
The unattenuated intensities, I o,El and I o,Eh , for the low and high energy, respectively, were calculated through the entrance surface dose (K a ) (Equation (3)) for all photon energies between 40 and 140 keV at 1 keV increments, according to Equation (4).
where I o,Ei is the unattenuated intensity (photons/mm 2 ) at energy E i . The X-ray mass energy absorption coefficient of air, (µ en /ρ) air , was obtained from the literature [92]. K a was set at 8 mGy and was evenly split between the low-/high-energy. The attenuated intensities, I s,Ei and I c,Ei , were calculated with the use of energydependent linear attenuation coefficients, which were obtained by multiplying the mass attenuation coefficients from published data (NIST) [92] and the corresponding densities. In order to differentiate the mineral crystals, the attenuated intensities after passing through the torso were used. For both simulation and experimental studies, for every HAp thickness the equivalent thickness of CaCO 3 and CaC 2 O 4 was defined, leading to equivalent attenuation. The examined HAp kidney/uteric stone and calcification in atherosclerotic plaque thicknesses ranged from 0.5 mm to 3 mm, at 100 µm increments. The thicknesses of calcium carbonate and calcium oxalate, which resulted in an equal number of photons after the attenuation, were calculated and averaged in order to obtain the corresponding equivalent thicknesses. The standard deviation (SD) of equivalent thicknesses, for both CaCO 3 and CaC 2 O 4 , was lower than 24.42 µm. The HAp thicknesses and the equivalent thicknesses of calcium carbonate and calcium oxalate are presented in Table 1 for a torso thickness of 20.8 cm [93-96]. The effective calcium/phosphorus mass ratio (m Ca /m P ) was calculated for all mineral crystal types using Equation (5) [79,97].
for the low-and high-energy, respectively.
for the low-and high-energy, respectively. µ Ca , µ PO4 and µ s are the energy-depended linear attenuation coefficients (1 cm −1 ) for calcium, phosphate and soft tissue, respectively. For all calcification types, the m Ca /m P was estimated using the equation above. The modyfied analytical model of the present study calculates the m Ca /m P of HAp where both calcium and phosphorus are found in its molecule. However, phosphorus is not present in CaCO 3 and CaC 2 O 4 molecules. Since the aim of the present study is to characterize an unknown mineral type based on attenuation intensity measurements, in the m Ca /m P calculation the linear attenuation coefficients of PO 4 are used, regardless of the calcification type. The m Ca /m P is mentioned as "effective" as the amount of CO 3 and C 2 O 4 corresponds to a smaller amount of PO 4 in order to preserve equal photon beam attenuation. This is attributed to the fact that CO 3 and C 2 O 4 have lower linear attenuation coefficients compared to PO 4 [79]. Hence, the calculated effective m Ca /m P for these mineral crystals is expected to be lower.
In an incident beam with a number of I photons (considering Poisson distribution), the coefficient of variation (CV) can be calculated as CV = √ I/I. The use of the error theory and the performance of derivations [79] resulted in Equation (6) for the simulation study with analytical methods. The coefficient of variation of the effective m Ca /m P (CV m Ca /m P (%)) was calculated according to Equation (6) for all mineral crystals [79]. The minimization of CV m Ca /m P (%) was the criterion to select the optimum energy pair. The CV m Ca /m P (%) of each mineral crystal is referred to as CV m Ca /m P ,H Ap (%), CV m Ca /m P ,CaCO 3 (%) and CV m Ca /m P ,CaC 2 O 4 (%) for hydroxyapatite, calcium carbonate and calcium oxalate, respectively.

Polyenergetic X-rays
Unfiltered spectra from TASMIP spectral models generated from a tungsten anode were used to obtain low-/high-energy spectra. The mean energies of the spectra should be similar to those indicated by the monoenergetic study [98]. The peak tube voltages were 50 kVp for the low energy and 90 kVp for the high energy. The examined filter materials for the low energy were selected based on their K-edges. In addition to these filters, aluminum was also examined. Although aluminum does not have a strong K-edge that could modify a spectrum into a quasi-monoenergetic spectrum, it is a common filter used in X-ray systems as it attenuates the low energies in spectra. In lanthanide filters, which are those with the most effective K-edges, there is no filter with a K-edge higher than 64 keV [79]. Hence, the selection of high energy filter materials was based on the filters' density. In Table 2, the filter materials with their K-edges and densities are presented for the low-and high-energy, respectively. All filter materials were examined for thicknesses ranging from 0.01 to 3 cm at 0.01 cm increments. The filters for the low-/high-energy were applied to unfiltered spectra obtained at 200 and 400 mAs, respectively. The Radcal 2026C (Radcal Corporation, Monrovia, CA, USA) ionization chamber [99] was placed at 66 cm from the tube output in order to measure the entrance surface doses for the low-/high-energy. The measured doses were 2.31 mGy and 30.50 mGy for the low-/high-energy, respectively.
A nonparametric statistical analysis was conducted with the use of the attenuated intensities (I s,Ei and I c,Ei ). Assuming Poisson distribution for the I s,Ei and I c,Ei , 5000 random values were generated for each one, leading to 5000 effective m Ca /m P values. The calculations were accomplished for HAp, CaCO 3 and CaC 2 O 4 in order to determine the distribution of m Ca /m P of each mineral crystal. Based on a previous work of our team, the normal Kernel distribution describes the random variable, effective m Ca /m P [79]. For all mineral crystal types and thicknesses, the Kernel probability functions were obtained and the false positive and false negative values were calculated.

Irradiation Process
The radiographic system Del Medical Eureka (DEL MEDICAL, Harrison, NY and Bloomingdale, IL, USA) with a tungsten (W) anode was used in the experimental evaluation of the method [100]. The low energy images were acquired at 50 kVp and 200 mAs with the inherent filtration of 3 mm aluminum, while the high energy images were acquired at 90 kVp, 400 mAs and added filtration of 3 mm 99.99% pure copper (Cu) foil. The detection system was a terbium-doped gadolinium oxysulfide (Gd 2 O 2 S:Tb) phosphor screen (Min-R 2190 with a mass thickness of 33.91 mg cm −2 ) coupled to an optical readout device including a complementary metal-oxide-semiconductor (CMOS) Remote RadEye HR (Rad-icon Imaging Corp, Santa Clara, CA, USA) photodiode pixel array. The CMOS photodiode array has a format of 1200 × 1600 pixels. The active area is 27 × 36 mm 2 , and the pixel pitch is 22.5 µm. The Gd 2 O 2 S:Tb screen was directly overlaid onto the active area of the CMOS [87]. The source-to-detector distance (SDD) was set at 66 cm. A schematic representation of the experimental configuration is illustrated in Figure 2.
The analytical model requires the total number of photons of the attenuated intensities for the determination of effective m Ca /m P . Thus, the pixel values of the low-and highenergy images should be converted to photons. To this aim, the detector response curve was measured and used in the conversion of pixel values to dose and finally to photons, since the I s,Ei and I c,Ei are known. Hence, Equation (3) is rewritten as [79,101]: R 2190 with a mass thickness of 33.91 mg cm −2 ) coupled to an optical readout device including a complementary metal-oxide-semiconductor (CMOS) Remote RadEye HR (Radicon Imaging Corp, Santa Clara, CA, USA) photodiode pixel array. The CMOS photodiode array has a format of 1200 × 1600 pixels. The active area is 27 × 36 mm 2 , and the pixel pitch is 22.5 μm. The Gd2O2S:Tb screen was directly overlaid onto the active area of the CMOS [87]. The source-to-detector distance (SDD) was set at 66 cm. A schematic representation of the experimental configuration is illustrated in Figure 2.
For the measurement of the detector response curve, the Radcal 2026C ionization chamber was placed at 66 cm from the output of the X-ray tube. The entrance surface air kerma was measured for several tube current-time products. The low-and high-energy images were acquired in 16 bit and regions of interest (ROIs) of 500 × 500 pixels were For the measurement of the detector response curve, the Radcal 2026C ionization chamber was placed at 66 cm from the output of the X-ray tube. The entrance surface air kerma was measured for several tube current-time products. The low-and high-energy images were acquired in 16 bit and regions of interest (ROIs) of 500 × 500 pixels were measured from each low-and high-energy image. The mean pixel value (MPV) and the standard deviation (SD) of each ROI were calculated. The relationship between MPV and detector entrance dose was determined using linear regression, and the corresponding graph is presented in Figure 3.
Crystals 2021, 11, 345 9 of 21 measured from each low-and high-energy image. The mean pixel value (MPV) and the standard deviation (SD) of each ROI were calculated. The relationship between MPV and detector entrance dose was determined using linear regression, and the corresponding graph is presented in Figure 3.

Phantoms
As was stated earlier, the spotty calcifications found in atherosclerotic plaques are 3 mm or smaller [43,[58][59][60], and although kidney/uteric stones of the same size do not need surgery (pass spontaneously), the knowledge of their composition is important. The 3 mm was set as the highest examined mineral crystal thickness based on the fact that 3 mm is the upper limit for a calcification to be considered as spotty calcification. The discrimination between HAp and either of the two mineral crystal types resulted in lower false positive (%), false negative (%) and overlap area (%) values compared with the discrimintation between CaCO3 and CaC2O4. For the experimental evaluation of the method, hydroxyapatite phantoms of 0.7 mm, 1 mm, 1.5 mm and 3 mm thickness were constructed. Phantoms with the equivalent thicknesses of CaCO3 and CaC2O4 were also constructed,

Phantoms
As was stated earlier, the spotty calcifications found in atherosclerotic plaques are 3 mm or smaller [43,[58][59][60], and although kidney/uteric stones of the same size do not need surgery (pass spontaneously), the knowledge of their composition is important. The 3 mm was set as the highest examined mineral crystal thickness based on the fact that 3 mm is the upper limit for a calcification to be considered as spotty calcification. The discrimination between HAp and either of the two mineral crystal types resulted in  [85]. Each mineral crystal was mixed with the same quantity of epoxy resin and placed into silicon tubes (22.31 mm inner diameter). An epoxy resin phantom, with the same quantity of epoxy used in the mineral mixtures, was constructed and positioned in the beam path when measurements without mineral crystal phantoms were implemented. A vacuum tube VACIOTEM-T (SELECTA, Barcelona, Spain) was used to diminish the remaining air in the phantoms. Each phantom was immersed in a PMMA tank with dimensions of 20 × 20 × 25 cm 3 filled with double distilled water (20.8 cm thickness) to simulate the torso. Figure 4 illustrates the constructed mineral crystal phantoms for the experimental verification.  Ten low and ten high-energy images were acquired for each phantom for the determination of the effective ⁄ and the CV ( ). After every low-and high-energy image combination acquisition, the phantoms were repositioned in order to take into consideration the possible inhomogeneity of the phantoms. Twenty ROIs of 44 × 44 pixels, corresponding to 0.99 × 0.99 mm 2 , were obtained in each image as this area is close to the dimensions of kidney/uteric stones or atherosclerotic plaques' calcifications and achieving an adequate number of photons [79]. The MPVs of the ROIs were converted into photon values for the calculation of the effective ⁄ . Figure 5 illustrates indicative low energy images of the 3 mm thick phantoms and the profiles of each phantom image. The graph with the profiles shows that the calcium oxalate and calcium carbonate phantoms with equivalent thicknesses of 3 mm resulted in almost equal photon beam attenuation to the 3 mm hydrohyapatite phantom. Ten low and ten high-energy images were acquired for each phantom for the determination of the effective m Ca /m P and the CV (CV m Ca /m P ,experimental (%) = SD m Ca /m P /mean m Ca /m P ). After every low-and high-energy image combination acquisition, the phantoms were repositioned in order to take into consideration the possible inhomogeneity of the phantoms. Twenty ROIs of 44 × 44 pixels, corresponding to 0.99 × 0.99 mm 2 , were obtained in each image as this area is close to the dimensions of kidney/uteric stones or atherosclerotic plaques' calcifications and achieving an adequate number of photons [79]. The MPVs of the ROIs were converted into photon values for the calculation of the effective m Ca /m P . Figure 5 illustrates indicative low energy images of the 3 mm thick phantoms and the profiles of each phantom image. The graph with the profiles shows that the calcium oxalate and calcium carbonate phantoms with equivalent thicknesses of 3 mm resulted in almost equal photon beam attenuation to the 3 mm hydrohyapatite phantom. mm 2 , were obtained in each image as this area is close to the dimensions of kidney/uteric stones or atherosclerotic plaques' calcifications and achieving an adequate number of photons [79]. The MPVs of the ROIs were converted into photon values for the calculation of the effective ⁄ . Figure 5 illustrates indicative low energy images of the 3 mm thick phantoms and the profiles of each phantom image. The graph with the profiles shows that the calcium oxalate and calcium carbonate phantoms with equivalent thicknesses of 3 mm resulted in almost equal photon beam attenuation to the 3 mm hydrohyapatite phantom. The Cadmium Telluride AMPTEK XR-100T (Amptek, Inc., Bedford, MA, USA) was used to measure the filtered spectra for the low-and high-energy, respectively, and thus the mean pixel values were converted into the number of photons [79,87,102,103]. The The Cadmium Telluride AMPTEK XR-100T (Amptek, Inc., Bedford, MA, USA) was used to measure the filtered spectra for the low-and high-energy, respectively, and thus the mean pixel values were converted into the number of photons [79,87,102,103]. The spectrophotometer was calibrated for energy scales, linearity checks and energy resolution with the use of 125 I and 99m Tc γ-ray calibration sources [101]. The total entrance dose was 8.38 mGy.
In order to assess the experimental results, a statistical analysis was conducted. As the effective m Ca /m P probability distribution is not Gaussian, the Mann-Whitney U non-parametric statistical test was applied to the measured effective m Ca /m P values [79]. Through the two-tailed Mann-Whitney U non-parametric statistical test, we investigated if there was a statistically significant difference between the mineral crystal types tested [79,104]. For all tests, the significance level was set at 5%.  Table 1. In all cases, for low energies ranging from 40 to 50 keV combined with all high energies, the CV m Ca /m P (%) ranged from the minimum CV m Ca /m P (%) value up to +20% of the minimum. spectrophotometer was calibrated for energy scales, linearity checks and energy resolution with the use of 125 I and 99m Tc γ-ray calibration sources [101]. The total entrance dose was 8.38 mGy.

Results and Discussion
In order to assess the experimental results, a statistical analysis was conducted. As the effective ⁄ probability distribution is not Gaussian, the Mann-Whitney U nonparametric statistical test was applied to the measured effective ⁄ values [79]. Through the two-tailed Mann-Whitney U non-parametric statistical test, we investigated if there was a statistically significant difference between the mineral crystal types tested [79,104]. For all tests, the significance level was set at 5%.  Table 1. In all cases, for low energies ranging from 40 to 50 keV combined with all high energies, the ⁄ (%) ranged from the minimum ⁄ (%) value up to +20% of the minimum.

The effective
⁄ is calculated with the use of a three-equation system. As the low-and high-energy converge, they lead to linear dependency. This results in the increase in the ⁄ (%), which is depicted in Figure 6. Figure 7 shows the mean energies and photons/mm 2 of the examined filters, as a function of the surface density (g cm −2 ). Figure 7b illustrates the results for the low energy filters. Although aluminum resulted in lower mean energies compared to the other filters, it provided the highest number of photons through the whole surface density range, which is an important factor for the analytical model. Additionally, aluminum is a lowcost filter commonly used for beam filtration in X-ray systems. Figure 7c,d show that in both mean energy and the number of photons, the high energy filters appeared to be grouped. As in mean energy (Figure 7c), antimony, tin and bismuth resulted in higher mean energies, while chromium and vanadium resulted in lower mean energies. When the aforementioned filters were examined for the number of photons, the results of the The effective m Ca /m P is calculated with the use of a three-equation system. As the low-and high-energy converge, they lead to linear dependency. This results in the increase in the CV m Ca /m P (%), which is depicted in Figure 6. Figure 7 shows the mean energies and photons/mm 2 of the examined filters, as a function of the surface density (g cm −2 ). Figure 7b illustrates the results for the low energy filters. Although aluminum resulted in lower mean energies compared to the other filters, it provided the highest number of photons through the whole surface density range, which is an important factor for the analytical model. Additionally, aluminum is a low-cost filter commonly used for beam filtration in X-ray systems. Figure 7c,d show that in both mean energy and the number of photons, the high energy filters appeared to be grouped. As in mean energy (Figure 7c), antimony, tin and bismuth resulted in higher mean energies, while chromium and vanadium resulted in lower mean energies. When the aforementioned filters were examined for the number of photons, the results of the groups mentioned above were reversed. Antimony, tin and bismuth resulted in a lower number of photons, while chromium and vanadium resulted in a higher number of photons. Regardless, the examined factors gallium and copper resulted in intermediate values. More specifically, the mean energies of copper and gallium were on average 7.98% lower than chromium and vanadium, while the mean energies of antimony, tin and bismuth were on average 15.55% lower than chromium and vanadium. As for the number of photons, chromium and vanadium resulted in 5.40 times (on average) higher number of photons than copper and vanadium and 38.52 times (on average) higher than antimony, tin and bismuth. Thus, the use of copper or gallium is a good compromise among the factors examined. Since vanadium cannot be found as foil, copper was selected as the high energy filter. groups mentioned above were reversed. Antimony, tin and bismuth resulted in a lower number of photons, while chromium and vanadium resulted in a higher number of photons. Regardless, the examined factors gallium and copper resulted in intermediate values. More specifically, the mean energies of copper and gallium were on average 7.98% lower than chromium and vanadium, while the mean energies of antimony, tin and bismuth were on average 15.55% lower than chromium and vanadium. As for the number of photons, chromium and vanadium resulted in 5.40 times (on average) higher number of photons than copper and vanadium and 38.52 times (on average) higher than antimony, tin and bismuth. Thus, the use of copper or gallium is a good compromise among the factors examined. Since vanadium cannot be found as foil, copper was selected as the high energy filter. A study, based on the nonparametric statistical analysis, was conducted in order to find the proper thicknesses for both aluminum and copper used to modify the X-ray spectra. From all the examined filter thickness combinations, ranging from 0.01 to 3 cm at 0.01 cm increments, the false negative (%), false positive (%) and the overlap area (%) were estimated. Based on the results of this study, the irradiation conditions were the following: (i) for the low energy (50 kVp) images, a 3 mm aluminum filter was applied (inherent filtration), and 3 mm of copper for the high energy (90 kVp) images. This filter thickness combination led to lower false negative (%) and false positive (%) values, and as a result a lower overlap area (%) value. This can be attributed to the fact that both low-and highenergy spectra resulted in an adequate number of photons. Additionally, the aforementioned filter thickness led to spectral separation between the low-and high-energy spectra, which is a key factor for mineral characterization. The mean energies were 34 and 72 keV for the low-and high-energy, respectively. Table 3 shows the results of the nonparametric statistical analysis for the above irradiation conditions. The false negative (%), false positive (%) and the overlap area (%) are A study, based on the nonparametric statistical analysis, was conducted in order to find the proper thicknesses for both aluminum and copper used to modify the X-ray spectra. From all the examined filter thickness combinations, ranging from 0.01 to 3 cm at 0.01 cm increments, the false negative (%), false positive (%) and the overlap area (%) were estimated. Based on the results of this study, the irradiation conditions were the following: (i) for the low energy (50 kVp) images, a 3 mm aluminum filter was applied (inherent filtration), and 3 mm of copper for the high energy (90 kVp) images. This filter thickness combination led to lower false negative (%) and false positive (%) values, and as a result a lower overlap area (%) value. This can be attributed to the fact that both lowand high-energy spectra resulted in an adequate number of photons. Additionally, the aforementioned filter thickness led to spectral separation between the low-and high-energy spectra, which is a key factor for mineral characterization. The mean energies were 34 and 72 keV for the low-and high-energy, respectively. Table 3 shows the results of the nonparametric statistical analysis for the above irradiation conditions. The false negative (%), false positive (%) and the overlap area (%) are presented for mineral crystal thicknesses ranging from 0.50 to 3.00 mm, at 0.  Regardless of the examined case, the overlap areas corresponding to smaller thicknesses of the mineral crystals were increased compared to higher thicknesses. Additionally, when the case of discrimination between CaCO 3 and CaC 2 O 4 was investigated, the overlap areas were increased compared to the corresponding HAp; CaCO 3 and HAp; CaC 2 O 4 cases for all thicknesses. Table 4 shows the averaged measured effective m Ca /m P , CV m Ca /m P ,experimental (%) values and the results of the Mann-Whitney U test for the examined thicknesses. In all cases, the critical U value was 127 [104] since in each phantom image 20 ROIs were selected. Based on the results shown in Table 4, there are statistically significant differences (p < 0.05) between the mineral crystals examined for thicknesses of 0.70 mm or higher for HAp; CaCO 3 and HAp; CaC 2 O 4 cases. The CV m Ca /m P ,experimental (%) values of all mineral crystals were found to be higher than the ones calculated in the simulation study. This can be attributed to the fact that the effective mass attenuation coefficients were used instead of the mass attenuation coefficients, due to the fact that the neighboring energy pairs that result in linear dependence of the equation system are not taken into consideration. Additionally, in the determination of the effective mass attenuation coefficients, the effect of beam hardening is not considered [74]. Finally, the instability of both the electronics and the X-ray tube was an additional factor to the above that contributed to the CV m Ca /m P ,experimental (%) value increase.

Experimental Verification
The proposed dual energy method characterizes mineral crystals of 0.70 mm thick or higher, found in kidney/uteric stones and calcifications in atherosclerotic plaques. The statistical noise that has an impact on the number of photons finally detected was the main reason for the increase in the CV m Ca /m P ,experimental (%) value and the higher U values in the case of CaCO 3 ; CaC 2 O 4 . The latter case was proved to be the most difficult among the three cases examined. Since CaCO 3 is rarely found in kidney/uteric stones [23,29,30] and CaC 2 O 4 and hydroxyapatite are the main calcium mineral crystals found in atherosclerotic plaques [42,43,[45][46][47][48][49][50][51][52][53][54][55], the results indicate that the mainly found mineral crystals' characterization can be accomplished by the proposed method.
The discrimination of calcifications smaller than 3.00 mm (spotty calcifications) is of great importance. Apart from the fact that calcifications are clinical indicator for plaques' instability [43,44,49,[58][59][60][61][62][63][64][65][66], the presence of hydroxyapatite increases the possibility of rapture. The results of this work are also important for patients diagnosed with urolithiasis or nephrolithiasis. Despite the fact that kidney/uteric stones smaller than 3.00 mm, such as the ones examined in this work, do not need surgery, the knowledge of their composition will assist the physician to plan a better treatment for the patient and may prevent surgery [3,[9][10][11][12][13][14][15][16][17].
The results presented above, for both simulation study and experimental verifications, can be improved by increasing the entrance dose without risking the patients' health. The low energy irradiation could guide the physician to identify and locate the kidney/uteric stones or atherosclerotic plaques' calcifications, while the high energy irradiation can be restricted to a smaller area including only the stones or calcifications. As it has already been mentioned, the total entrance dose was 8.38 mGy, resulting in 0.22 mSv and 3.84 mSv effective equivalent doses for urolithiasis/nephrolithiasis and atherosclerotic plaques, respectively. Previous studies for discrimination between different mineral crystal types found in kidney/uteric stones using CT resulted in an entrance dose of 10 mGy and an effective dose of 8 up to 10 mSv [105][106][107][108][109]. When DECT was used in phantom studies, the entrance dose ranged from 6 to 26.2 mGy and 8.2 to 14.6 mGy when patients were involved [10,[110][111][112][113]. Low dose (LDCT) and ultra-low dose CT (ULDCT) were also used, resulting in an absorbed dose of 3 mSv or less [105,[107][108][109]114]. In all cases, there were difficulties in discrimination between mineral crystals when the stones were 3.00 mm or smaller, and especially when LDCT and ULDCT were used [105,[107][108][109]114]. As for atherosclerotic plaques, the effective dose of CT ranges from 5.4 to 8 mSv, while the typical entrance dose for DECT is 19 mGy and 23 mGy for the low-and high-energy, respectively [71,[115][116][117][118].
The radiation field that was used in order to obtain the low-and high-energy images was almost equal to the dimensions of the CMOS detector (27 × 36 mm 2 ) of the DE system. Although the radiation field's dimensions are small, the scatter radiation still exists but in a smaller amount compared to that of usual radiation fields used in radiography. When the kidney/uteric stones or atherosclerotic plaques' calcifications are identified in a radiographic image, the field will be even smaller so only the kidney/uteric stones or atherosclerotic plaques' calcifications will be irradiated.
In order to achieve an adequate photon flux, especially in the high energy where a 3 mm Cu filter is applied, the time-current product was increased to 400 mAs. The results of the proposed method will be further improved if a specialized X-ray system is used to preserve a higher photon flux. As it is already mentioned, apart from the adequate number of photons, spectral separation is also an important parameter which improves mineral characterization. An X-ray system with higher photon flux and thicker filters will further improve the spectral separation. In future work, a more realistic torso phantom will be constructed in order to further evaluate the method, and Monte Carlo methods or convolution kernel-based approaches will be adopted for scatter correction [118][119][120][121].

Conclusions
In this work, a dual energy method was proposed for the discrimination of mineral crystals found in kidney/uteric stones and calcifications in atherosclerotic plaques through the effective m Ca /m P . The potential of the method to differentiate HAp, CaCO 3 and CaC 2 O 4 was investigated using both simulations (monoenergetic and polyenergetic X-ray beams) and experimental verification. The experimental results revealed that with the proposed method, the mineral crystals can be characterized for thicknesses of 0.70 mm or higher for HAp; CaCO 3 and HAp; CaC 2 O 4 cases. Further studies will show if the proposed method can provide in vivo information for the mineral crystals.