Next Article in Journal
H-BIM and Artificial Intelligence: Classification of Architectural Heritage for Semi-Automatic Scan-to-BIM Reconstruction
Previous Article in Journal
Key Technologies and Evaluation of a MiniSAR Experimental System for Unmanned Underwater Vehicle Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

X-ray Single Exposure Imaging and Image Processing of Objects with High Absorption Ratio

1
College of Physics, Changchun University of Science and Technology, Changchun 130022, China
2
Electronic Information Engineering College, Changchun University, Changchun 130022, China
*
Authors to whom correspondence should be addressed.
Sensors 2023, 23(5), 2498; https://doi.org/10.3390/s23052498
Submission received: 12 January 2023 / Revised: 19 February 2023 / Accepted: 20 February 2023 / Published: 23 February 2023
(This article belongs to the Section Sensing and Imaging)

Abstract

:
The dynamic range of an X-ray digital imaging system is very important when detecting objects with a high absorption ratio. In this paper, a ray source filter is used to filter the low-energy ray components which have no penetrating power to the high absorptivity object to reduce the X-ray integral intensity. This enables the effective imaging of the high absorptivity objects and avoids the image saturation of low absorptivity objects, thus achieving single exposure imaging of high absorption ratio objects. However, this method will reduce the image contrast and weaken the image structure information. Therefore, this paper proposes a contrast enhancement method for X-ray images based on Retinex. Firstly, based on Retinex theory, the multi-scale residual decomposition network decomposes the image into an illumination component and a reflection component. Then, the contrast of the illumination component is enhanced through the U-Net model with the global–local attention mechanism, and the reflection component is enhanced in detail using the anisotropic diffused residual dense network. Finally, the enhanced illumination component and the reflected component are fused. The results show that the proposed method can effectively enhance the contrast in X-ray single exposure images of the high absorption ratio objects, and can fully display the structure information of images on devices with low dynamic range.

1. Introduction

In the process of digital radiography (DR), objects with a high absorption ratio need higher X-ray photon energy and exposure, while objects with a low absorption ratio need the opposite [1]. This requires an X-ray flat panel detector (FPD) to have a high dynamic range (HDR). The dynamic range of the FPD determines the imaging effect of objects with a high absorption ratio. In fact, the dynamic range of the FPD is limited [2]. When high absorptivity objects can be effectively imaged and resolved, images of low absorptivity objects usually reach saturation and cannot be resolved. Therefore, it is necessary to use an X-ray HDR imaging method to solve the imaging problem of objects with a high absorption ratio [3].
At present, there are two main methods for X-ray HDR imaging [4]. One is that the low-dynamic-range imaging device is used to make multiple exposures of the same scene with variable energy, and then the obtained multiple images are fused into an HDR image [5]. Another method is that the optical path expansion technology is used to expand the light intensity information in the same scene into multiple optical signals, the exposure parameters of imaging devices in different optical paths are set to achieve imaging, and the output images of multiple sensors are fused into an HDR image with complete data information [6].
X-ray HDR imaging of high absorption ratio objects usually uses multiple exposures of different energies. Image fusion keeps the X-ray image information of high absorptivity and low absorptivity objects in one image [7]. In recent years, Chen et al. fused sequences of different voltages at each projection angle, thus extending the dynamic range of the X-ray imaging system [8]. Haidekker et al. proposed a mathematical framework for recovering absorbance from images of objects with high absorption ratios, thereby enhancing the dynamic range of X-ray imaging [9]. Li et al. combined self-adaptive multi-exposure imaging technology and HDR image fusion, and proposed a method for extending the dynamic range of X-ray imaging systems based on linear constraints with variable energy [10].
The above research shows that multi-exposure imaging technology is an effective method to expand the dynamic range of the X-ray imaging system. However, the exposure time and image transmission time of DR are relatively long, resulting in the low detection efficiency of multi-exposure, and the change in the geometric position of the object to be measured will lead to image fusion difficulties [11]. At the same time, the complex imaging system structure of optical path extension technology makes it difficult to determine the exposure conditions [12]. Based on the above reasons, this paper proposes a single exposure imaging technique for objects with a high absorption ratio based on X-ray spectral bandwidth compression. In order to achieve the compression of the X-ray spectral bandwidth, a ray source filter is used to filter the low energy rays that cannot pass through high absorptivity objects, thus reducing the integral intensity of the X-rays that pass through low absorptivity objects and avoiding image saturation. This is equivalent to compressing the dynamic range of the image at the X-ray image end, which will reduce the image contrast. Aiming for the features of such images, in this paper, an X-ray image contrast enhancement method based on Retinex is designed, which realizes that all information can be obtained via the single imaging of objects with a high absorption ratio.
To be specific, a multi-scale residual network is designed to decompose an image into two components (the illumination component and the reflection component). A global–local attention mechanism is proposed for U-Net networks to enhance the contrast in illumination components. An anisotropic diffusion-based residual dense network is designed to deal with the reflection component. Meanwhile, reconstruction loss, illumination smoothness loss, and constant reflectivity loss are introduced to train the decomposition network. Structural similarity loss is introduced to train the enhanced network.
The main ideas and contributions of this paper are summarized as follows:
  • X-ray spectrum bandwidth compression is realized by using a ray source filter, which avoids the image saturation of low absorptivity objects.
  • In view of the low contrast of X-ray images, this paper proposes an X-ray contrast enhancement method based on Retinex. According to the invariant property of the reflection component, the decomposition network is trained using single exposure and multiple exposure X-ray images.
  • Aiming at the phenomenon of uneven global brightness in the illumination components, this paper proposes a U-Net network based on the global–local attention mechanism to extract the illumination features of different scales, and use the complementarity of global and local features to enhance the contrast in the illumination components.
  • In order to avoid ghost artifacts and detail blurring of reflection components, this paper designs a residual dense network based on anisotropic diffusion, which uses hierarchical information to enhance detailed information while suppressing noise and artifacts.
The rest of this paper is organized as follows: In Section 2, the influence of the filter on the X-ray spectrum is introduced and the experiment is carried out. In Section 3, the network structure of the proposed contrast enhancement algorithm is analyzed, and the effectiveness of the method is verified via simulation experiments. Finally, the conclusion is given in Section 4.

2. X-ray Spectral Bandwidth Compression and Imaging

X-ray imaging equipment is composed of an X-ray tube, high voltage power supply, and X-ray FPD. In this paper, the tungsten target X-ray tube (tube voltage 180 kV, tube current 1 mA, focus size 0.5 mm) was used, the FPD size was 17 inches × 17 inches, the resolution was 3072 × 3072, and the scintillator was CsI crystal. The exposure parameters were tube current of 0.5 mA, exposure time of 2 s, and adjustable tube voltage. The geometry and position of the X-ray tube, imaging sample, and device are shown in Figure 1.

2.1. Effect of Filter on X-ray Spectrum

The absorption of X-rays by substances conforms to the Lambert–Beer law, that is, the intensity of emitted X-rays decays exponentially with the thickness, and the intensity of the X-rays passing through the object is:
I = I 0 e μ m · ρ x
where I 0 is the initial X-ray intensity, μ m is the mass absorption coefficient, and x is the thickness of the object. In the traditional model, the mass absorption coefficient can be expressed as:
μ m = K λ 3 Z 4
In the formula, K is the constant, λ is the X-ray wavelength, and Z is the atomic number. As can be seen from the formula, for a substance with atomic number Z , the mass absorption coefficient rapidly decreases exponentially with the decrease in X-ray wavelength. The nonlinear attenuation of low energy rays and high energy rays can be realized by utilizing the nonlinearity of X-ray absorption by matter, and the purpose of X-ray spectral bandwidth compression can be achieved. Figure 2 shows the relationship curve of the X-ray radiation spectrum with the filter plate at the maximum energy of 180 kV (to simplify the calculation, the characteristic emission peak of the target material and the characteristic absorption peak of the material are ignored). It can be seen from the comparison of the X-ray spectrum curves transmitted by three filters that the larger the atomic number of the filter material, the narrower the X-ray spectrum width of the filter.
By using the filter to realize the X-ray spectral width compression characteristic, the low energy X-ray can be greatly attenuated, while the high energy X-ray can be mostly retained, thus greatly reducing the X-ray integral intensity of penetrating low absorptivity objects and avoiding the saturation of the FPD.

2.2. Analysis of the Influence of Narrowband X-ray on Imaging

2.2.1. Image Gray Distribution Test

The step test block shown in Figure 3a was used for the imaging test. The filter material is silver, and the dynamic range of the FPD image is 16 bits. After imaging, the gray distribution curve of the iron test block and aluminum test block is obtained according to the position shown in Figure 3b.
Figure 4a,b shows the gray distribution curve of iron and aluminum test blocks under 90 kV exposure. The effective penetration depth of the iron test block is 8 mm, and the other six steps cannot be effectively imaged. The aluminum test block can be effectively imaged under various conditions. Figure 4c,d shows the gray distribution curve of the iron test block and aluminum test block under a 180 kV exposure condition. The iron test block can be effectively imaged under various conditions. Under the condition of no filter, the steps of the aluminum test block are saturated and cannot be effectively imaged, that is, the dynamic range of the X-ray transmission intensity exceeds the effective dynamic range of the FPD. With the increase in the thickness of the filter, the saturation phenomenon of the aluminum test block is gradually weakened. When the thickness of the filter reaches 1 mm, the thinnest part of the aluminum test block can also be effectively imaged, achieving a single exposure imaging of 20 mm iron and 1 mm aluminum. The results show that the X-ray image with an HDR can be well compressed into the effective dynamic range of the FPD by using the filter, so as to realize the single exposure imaging of objects with a high absorption ratio.

2.2.2. Luminance Compression Ratio

The ratio of the gray value of the image without the filter to the gray value with the filter is defined as the luminance compression ratio. The thickness of the test block that can be imaged normally under different exposure conditions is counted, and the results are shown in Figure 5. It can be seen from the curve that the luminance compression ratio increases with the increase in the thickness of the filter, and the smaller the absorbance is, the greater the compression ratio is. This indicates that the proportion of the remaining high energy rays in the filtered X-ray spectrum increases, the proportion of the integral intensity of the rays transmitted by high absorptivity objects increases, and the proportion of the integral intensity of the X-ray transmitted by low absorptivity objects decreases.

2.2.3. Image Contrast Feature

The use of a filter can achieve single exposure imaging of objects with a high absorption ratio, but with the compression of the dynamic range of X-ray images, the image contrast will be reduced. Under different exposure conditions, there are impenetrable phenomena of iron test blocks and image saturation phenomena of aluminum test blocks. Therefore, it is impossible to calculate the contrast generated by a gray scale of 20 mm iron and 1 mm aluminum. To this end, steps that can be distinguished under different exposure conditions are used for contrast statistics. In this paper, 4 mm iron and 6 mm iron are used. As shown in Figure 6, the image contrast decreases with the increase in the filter thickness.
In order to verify the change rule of the global maximum contrast of the image, take the gray value of the image at the thickest part of iron and the gray value of the image at the thinnest part of aluminum for the calculation. Under 150 kV and 180 kV exposure conditions, 20 mm iron can effectively be imaged, so only these two sets of data are taken. The “---” symbol in the table represents image saturation and cannot be calculated, so only partial data can be obtained in global contrast statistics. The rule in Figure 6 can also be obtained from Table 1, that is, under the same exposure voltage, the contrast decreases with the increase in filter thickness and decreases with the increase in working voltage. Contrast is one of the important parameters of the image, so it is necessary to enhance it to improve the human eye recognition ability.

3. X-ray Image Contrast Enhancement Method Based on Retinex

Inspired by the Retinex theory [13], a deep Retinex network was designed, and the overall architecture is presented in Figure 7. The network consists of three subnetworks: illumination and reflection decomposition network, illumination adjustment network, and reflection enhancement network.
The original image S was first decomposed into illuminance I ˜ and reflectance R ˜ :
I ˜ , R ˜ = f d e c o m S
where f d e c o m   signifies the image decomposition process. Illumination I represents the various lightness on objects. Reflectance R describes the intrinsic property of captured objects, which is considered to be consistent under any lightness conditions.
Then, the decomposed I ˜ and R ˜ were enhanced to obtain the enhanced illumination component I ^ and the restored reflection component R ^ :
I ^ = f e n h a n c e I ˜
R ^ = f r e s t o r e R ˜
Finally, the reconstructed image S ^ was obtained via element-wise multiplication:
S ^ = I ^ · R ^

3.1. Image Decomposition Network

At the decomposition stage, a multi-scale residual decomposition network (Decom-Net) was designed to generate the illumination component and the reflection component, as shown in Figure 8. Decom-Net uses a 3 × 3 convolution layer to generate the characteristic of the input image. Then, it uses six convolution layers with the activation function rectified linear unit (ReLU) to change the size of the feature map and learn the characteristics of the illumination part and the reflection part. After the second convolution, multi-scale feature extraction was added. Finally, convolution and Sigmoid function were used to map the learned image features onto illuminated images and reflected images, and then output them.
In the training process, paired single/multiple exposure images were taken as the inputs of Decom-Net. Since single and multiple exposure images of the same scene have the same reflectance, when they are fed into two Decom-Nets, the network weights are shared between them. In the testing process, only the single exposure image needs to be inputted.

3.2. Illumination Enhancement Network

The constructed illumination enhancement network (Enhance-Net) adopts the U-Net structure, which is composed of an encoder and decoder, as shown in Figure 9. In the encoder part, the input image is first encoded into a lower dimension and then entered into the global–local attention mechanism. Finally, the original resolution image is obtained via input to the decoder. The encoder–decoder structure is composed of four convolution blocks. The difference is that the encoder is composed of two convolution layers and a pooled down-sampling layer, and the decoder is composed of two convolution layers and a deconvolution up-sampling layer.
The global and local relationships in the illumination component have an important impact on the training of the network. In order to enhance the modeling ability of the network, the global–local attention mechanism (GLAM) is proposed, and its structure is shown in Figure 10. Enhance-Net integrates GLAM to filter the features inputted into the network, enhancing useful information and suppressing invalid information.
Global attention mechanism: Because the illumination distribution is not uniform and consistent, the designed global attention mechanism builds the relationship between each channel, making the illumination component uniform on the overall pixel. By calculating the covariance between channels in the feature map, the global correlation and the covariance matrix are obtained. Then, the matrix is row convolved and normalized to retain the structure information, and the attention weight of each channel is obtained. Finally, the original feature map is multiplied by the attention weight.
Local attention mechanism: Since different objects have different features in the illumination component, the local information of the illumination component is used to enhance contrast. Firstly, the feature graph is divided into several feature blocks of the same size, and the mean and standard deviation of the feature blocks are used to improve the ability of capturing feature information. The mean value reflects the brightness of the image. The standard deviation reflects the dispersion degree between the image pixels and the mean value, and the clarity degree of the edge of the illumination component. The features after the mean and standard deviation are convolved and normalized to obtain the weight map of the feature block, and then multiplied with the original feature block to obtain local features.

3.3. Reflectance Restore Network

As for the reflectance, it is usually contaminated by the blurring of details; the reflectance restoration network (Restore-Net) is designed to remove noise and enhance details, as shown in Figure 11.
In Restore-Net, reflectance passes through three residual dense blocks, each focused on a task without interruption by other groups. The central convolution layer in the residual dense block is set into two groups, which avoids influence from the previous layer and makes effective use of the hierarchical information. Image edge information is enhanced by anisotropic diffusion filtering (AD). Dense connections can effectively boost performance and preserve semantic information. To avoid the disappearance of network gradients, a residual skip connection to connect all of these dense blocks as well as the successive concatenation and convolution layer is added.

3.4. Loss Function

3.4.1. Decomposition Process

Motivated by Retinex-Net [14], the loss functions L D e c o m of Decom-Net include reconstruction loss L r e c o n , illumination smoothness loss L i s , and invariable reflectance loss L i r .
L D e c o m = L r e c o n + λ i s L i s + λ i r L i r
where λ i s and λ i r denote the coefficients of illumination smoothness and balance reflectance consistency, respectively.
Reconstruction loss: In Decom-Net, the similarity between the reconstructed results of the illumination and reflection components of the image’s own decomposition and the original image not only need attention. It is also necessary to pay attention to the similarity between the combined reconstruction results of the illuminance component and the reflection component of the paired image decomposition and the single/multiple exposure image. Reconstruction loss is shown in Equation (8):
L r e c o n = i = s e , m e j = s e , m e λ i j R i I j S j 1
where S is the original image, R is the reflection component, I is the illumination component, and λ i j is the weighting factor.
Illumination smoothness loss: The single/multiple exposure image decomposition does not have the same illumination components but should highlight the structure of the image and local details while maintaining overall smoothness. The traditional total variation (TV) minimization as a loss function will cause over-smoothing.
L i s = i = s e , m e I i exp λ g R i
where represents the gradient (including horizontal and vertical gradient), λ g represents the gradient-aware balance coefficient, and exp λ g R i relaxes the smoothing constraint at locations with more complex image structures and illumination discontinuities.
Invariable reflectance loss: According to the Retinex theory, it is known that the reflection components are invariant by the nature of the object itself, so the reflection components of single and multiple exposure images are similar. Furthermore, the invariable reflectance loss is used as the constraint of the training stage.
L i r = R s e R m e 1
where R s e is the reflection component of the single exposure image, R m e is the reflection component of the multiple exposure image, and 1 indicates the L1 norm operation.

3.4.2. Illumination Enhancement Process

The illumination enhancement loss L I E maximizes the structural similarity between the illumination components before enhancement I and after enhancement I ^ SSIM:
L I E = S S I M I ^ , I

3.4.3. Reflectance Restore Process

Reflectance component enhancement loss L R E constrains the enhanced reflection component R ^ . There is similarity between this and the initial reflection component R :
L R E = S S I M R ^ , R

3.4.4. Fusion Process

Fusion loss L F u s i o n limits the fusion image I ^ R ^ . There is similarity with the original image S :
L F u s i o n = S S I M I ^ R ^ , S

3.5. Experiment and Analysis

In order to verify the effectiveness of the X-ray image contrast enhancement algorithm proposed in this paper, the proposed method was compared with other contrast enhancement algorithms. These include X-ray contrast enhancement methods, CLAHE [15], LCM-CLAHE [16], methods of contrast enhancement based on deep learning, Retinex-Net [14] and Zero-DCE [17].

3.5.1. X-ray Images of Materials with High Absorption Ratio under Different Voltages

X-ray images of cables under different voltages (90 kV, 120 kV, 150 kV, and 180 kV) were selected. It can be seen from Figure 12 that the overall brightness of the original image is low, the contrast is not obvious, and the cable details are not prominent enough. With the increase in voltage, the internal structure of the cable gradually appears, but the overall brightness and contrast are still very low.
Qualitative analysis
Five methods were used to enhance the contrast in the original image. CLAHE and LCM-CLAHE improved the contrast of the image, but the brightness did not improve significantly. In the Retinex-Net and Zero-DCE methods, the image brightness is enhanced, but the overall over-exposure is not suitable for practical applications. The method proposed in this paper can reasonably improve the brightness and contrast of the image. At the same time, the detailed structure of the cable is also clearly visible, which is suitable for X-ray images with different voltages.
Quantitative analysis
The purpose of X-ray image enhancement is to enhance brightness and contrast and highlight details. Therefore, five evaluation indicators were selected to evaluate the performance of different methods. These included structure definition (NRSS), information entropy (entropy), spatial frequency (SF), average gradient (AG), and image standard deviation (STD). The higher their values, the better the image enhancement effect.
  • NRSS: reflects the clarity of the image;
  • AG: reflects the detail change rate of the image, which can be used to characterize the clarity of the image;
  • SF: reflects the overall spatial activity of the image;
  • STD: reflects the dispersion between the gray value and the average value of each pixel in the image, which can be used to evaluate the contrast of the image;
  • Entropy: represents the average information of the image, which is used to measure the information richness of the image.
Table 2 lists the quality indicators of each method for processing the above images, and the best value of these methods is highlighted in bold. It can be seen that the proposed method has the best evaluation index for all voltage values. Compared with the other four methods, the proposed method has higher clarity, contrast, and information richness.

3.5.2. Qualitative Analysis of X-ray Images of Objects with Different Absorption Ratios

Qualitative analysis
Figure 13 and Figure 14 are X-ray images of the brake pump and water pump at 90 kV. It can be seen that the structures of the brake pump and water pump are complex. After CLAHE and LCM-CLAHE processing, the image does not improve significantly, and local details are not prominent enough. The images processed using Retinex-Net and Zero-DCE have clear details, but there are artifacts and the overall contrast is not enhanced. The proposed method improves the contrast and brightness of the image, and enhances the detail structure to a certain extent. The method is suitable for X-ray images with complex details.
Figure 15 and Figure 16 are X-ray images of an iron pipe at 120 kV and pressure doubling boxes at 100 kV, respectively. Among them, the iron pipe and pressure doubling box are defective, but the original X-ray image cannot be well displayed. After CLAHE and LCM-CLAHE, the defect is still dark and cannot be clearly displayed. Retinex-Net and Zero-DCE will blur the image, and the defect still has low contrast. In terms of the proposed method, in addition to the increase in global brightness, the contrast in defects is better than the other methods. Therefore, the proposed method can effectively reproduce the structure and defects in the dark area.
Quantitative analysis
Table 3 objectively reflects that the proposed method is very effective in improving the overall visual effect and the detail effect of the image. NRSS, AG, and SF reflect the sharpness of the image. Table 3 shows that the proposed method is clearer than the original image and other methods. STD reflects the contrast of the image, and the proposed method is far superior to other algorithms in contrast. Entropy reflects the richness in image information. The proposed method is superior to the other four methods in improving information entropy, which indicates that the method can retain more information. In conclusion, compared with the other methods, the proposed method achieved the best results in general.

4. Conclusions

The ray source filter is used to filter the low energy X-ray, realizing the compression of the X-ray radiation spectrum bandwidth, effectively reducing the ray integral intensity, penetrating the low absorptivity object, avoiding the saturation of the low absorptivity object image, and preserving all of the information of the high absorption ratio object in the single exposure image. In order to solve the problems of low contrast and structural information weakening caused by the dynamic compression of X-ray images, this paper proposes a Retinex global–local contrast enhancement model. In the decomposition phase, Decom-Net is constructed to decompose the original X-ray image into an illumination component and reflection components. In Enhance-Net, the global–local attention mechanism is used to adjust the contrast in the illumination components. In Restore-Net, anisotropic diffusion filtering is added to the residual dense network to enhance the image details and suppress noise and artifacts. Finally, the enhanced X-ray image is obtained by multiplying the enhanced illumination component and reflection component. The experimental results show that this method can effectively improve the visual effect and contrast, enhance detail, suppress noise, and avoid artifacts. The single exposure imaging method of objects with high absorption ratios used in this paper has higher detection efficiency and environmental adaptability than the traditional multiple exposure image fusion method, and has important practical value in portable X-ray detection equipment.

Author Contributions

Conceptualization, Y.L. (Ye Li); methodology, Y.L. (Yanxiu Liu); software, K.L.; validation, D.D.; formal analysis, P.Z.; resources, Y.L. (Ye Li); writing—original draft preparation, Y.L. (Yanxiu Liu); writing—review and editing, K.L. and D.D.; visualization, Y.L. (Yanxiu Liu); All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sukovic, P.; Clinthorne, N.H. A method for extending the dynamic range of flat panel imagers for use in cone beam computed tomography. In Proceedings of the 2001 IEEE Nuclear Science Symposium Conference Record (Cat. No. 01CH37310), San Diego, CA, USA, 4–10 November 2001; pp. 1647–1651. [Google Scholar]
  2. Kuntz, K.D.; Atz, E.; Collier, M.R.; Collado-Vega, Y.; Connor, H.K.; Porter, F.S.; Sibeck, D.G.; Walsh, B.M. Making the invisible visible: X-ray imaging. In Magnetospheric Imaging; Elsevier: Amsterdam, The Netherlands, 2022; pp. 59–99. [Google Scholar]
  3. Wilbur, S.; Anastopoulos, C.; Angelmahr, M.; Asfis, G.; Koch, J.; Lindblom, M.; Lohwasser, K.; Margulis, W.; on behalf of the FleX-RAY Collaboration. Flexible X-ray imaging detectors using scintillating fibers. J. Instrum. 2022, 17, C10013. [Google Scholar] [CrossRef]
  4. Keerativittayanun, S.; Kondo, T.; Kotani, K.; Phatrapornnant, T.; Karnjana, J. Two-layer pyramid-based blending method for exposure fusion. Mach. Vis. Appl. 2021, 32, 48. [Google Scholar] [CrossRef]
  5. Nayar, S.K.; Mitsunaga, T. High dynamic range imaging: Spatially varying pixel exposures. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition. CVPR 2000 (Cat. No. PR00662), Hilton Head, SC, USA, 15 June 2000; pp. 472–479. [Google Scholar]
  6. Jin, X.; Chen, Z.B.; Ni, X.X.; Lu, Z.K. Research of HDRI system based on PID controller. Opt. Tech. 2008, 34, 547–551. [Google Scholar]
  7. Zhao, T.; Zhang, S.X. X-ray Image Enhancement Based on Nonsubsampled Shearlet Transform and Gradient Domain Guided Filtering. Sensors 2022, 22, 4074. [Google Scholar] [CrossRef] [PubMed]
  8. Chen, P.; Han, Y.; Pan, J. High-dynamic-range CT reconstruction based on varying tube-voltage imaging. PLoS ONE 2015, 10, e0141789. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Haidekker, M.A.; Morrison, L.D.; Sharma, A.; Burke, E. Enhanced dynamic range X-ray imaging. Comput. Biol. Med. 2017, 82, 40–48. [Google Scholar] [CrossRef] [PubMed]
  10. Li, Y.; Han, Y.; Chen, P. X-ray energy self-adaption high dynamic range (HDR) imaging based on linear constraints with variable energy. IEEE Photonics J. 2017, 10, 3400114. [Google Scholar] [CrossRef]
  11. Chen, P.; Li, L.; Zhang, X.; Wang, S.; Tan, A. Blind quality index for tone-mapped images based on luminance partition. Pattern Recognit. 2019, 89, 108–118. [Google Scholar] [CrossRef]
  12. Ma, K.; Zeng, K.; Wang, Z. Perceptual quality assessment for multi-exposure image fusion. IEEE Trans. Image Process. 2015, 24, 3345–3356. [Google Scholar] [CrossRef] [PubMed]
  13. Parihar, A.S.; Singh, K. A study on Retinex based method for image enhancement. In Proceedings of the 2018 2nd International Conference on Inventive Systems and Control (ICISC), Coimbatore, India, 19–20 January 2018; pp. 619–624. [Google Scholar]
  14. Li, X.; Shang, J.; Song, W.; Chen, J.; Zhang, G.; Pan, J. Low-Light Image Enhancement Based on Constraint Low-Rank Approximation Retinex Model. Sensors 2022, 22, 6126. [Google Scholar] [CrossRef] [PubMed]
  15. Reza, A.M. Realization of the contrast limited adaptive histogram equalization (CLAHE) for real-time image enhancement. J. VLSI Signal Process. Syst. Signal Image Video Technol. 2004, 38, 35–44. [Google Scholar] [CrossRef]
  16. Mohan, S.; Ravishankar, M. Modified contrast limited adaptive histogram equalization based on local contrast enhancement for mammogram images. In International Conference on Advances in Information Technology and Mobile Communication; Springer: Berlin/Heidelberg, Germany, 2012; pp. 397–403. [Google Scholar]
  17. Guo, C.; Li, C.; Guo, J.; Loy, C.C.; Hou, J.; Kwong, S.; Cong, R. Zero-reference deep curve estimation for low-light image enhancement. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition, Seattle, WA, USA, 13–19 June 2020; pp. 1780–1789. [Google Scholar]
Figure 1. (a) X-ray imaging equipment, and (b) geometric location diagram of the equipment.
Figure 1. (a) X-ray imaging equipment, and (b) geometric location diagram of the equipment.
Sensors 23 02498 g001
Figure 2. Influence curve of filter on X-ray spectrum. (a) Aluminum filter; (b) iron filter; (c) silver filter; and (d) comparison of filtering characteristics of three filters.
Figure 2. Influence curve of filter on X-ray spectrum. (a) Aluminum filter; (b) iron filter; (c) silver filter; and (d) comparison of filtering characteristics of three filters.
Sensors 23 02498 g002
Figure 3. (a) Step test block used in imaging experiment and (b) obtaining position of gray curve of X-ray image of test block.
Figure 3. (a) Step test block used in imaging experiment and (b) obtaining position of gray curve of X-ray image of test block.
Sensors 23 02498 g003
Figure 4. (a) Gray scale curve of iron image under 90 kV exposure; (b) gray scale curve of aluminum image under 90 kV exposure; (c) gray distribution curve of iron image under 180 kV exposure; and (d) gray distribution curve of aluminum image under 180 kV exposure.
Figure 4. (a) Gray scale curve of iron image under 90 kV exposure; (b) gray scale curve of aluminum image under 90 kV exposure; (c) gray distribution curve of iron image under 180 kV exposure; and (d) gray distribution curve of aluminum image under 180 kV exposure.
Sensors 23 02498 g004
Figure 5. Relation curve between image brightness compression ratio and filter thickness: (a) 180 kV iron brightness compression ratio change curve and (b) 90 kV aluminum brightness compression ratio curve.
Figure 5. Relation curve between image brightness compression ratio and filter thickness: (a) 180 kV iron brightness compression ratio change curve and (b) 90 kV aluminum brightness compression ratio curve.
Sensors 23 02498 g005
Figure 6. Relation curve between image contrast and filter thickness.
Figure 6. Relation curve between image contrast and filter thickness.
Sensors 23 02498 g006
Figure 7. The network architecture of the proposed method.
Figure 7. The network architecture of the proposed method.
Sensors 23 02498 g007
Figure 8. The structure of Decom-Net.
Figure 8. The structure of Decom-Net.
Sensors 23 02498 g008
Figure 9. The structure of Enhance-Net.
Figure 9. The structure of Enhance-Net.
Sensors 23 02498 g009
Figure 10. The detail architecture of GLAM.
Figure 10. The detail architecture of GLAM.
Sensors 23 02498 g010
Figure 11. The structure of Restore-Net.
Figure 11. The structure of Restore-Net.
Sensors 23 02498 g011
Figure 12. The contrast enhancement results of different methods (the voltage is 90 kV, 120 kV, 150 kV, and 180 kV from top to bottom). (a) The raw images; (b) CLAHE; (c) LCM-CLAHE; (d) Retinex-Net; (e) Zero-DCE; and (f) our method.
Figure 12. The contrast enhancement results of different methods (the voltage is 90 kV, 120 kV, 150 kV, and 180 kV from top to bottom). (a) The raw images; (b) CLAHE; (c) LCM-CLAHE; (d) Retinex-Net; (e) Zero-DCE; and (f) our method.
Sensors 23 02498 g012
Figure 13. The contrast enhancement results of different methods. Top: raw image and the enhancement results; bottom: zoomed-in views on the highlighted area. (a) The raw images; (b) CLAHE; (c) LCM-CLAHE; (d) Retinex-Net; (e) Zero-DCE; and (f) our method.
Figure 13. The contrast enhancement results of different methods. Top: raw image and the enhancement results; bottom: zoomed-in views on the highlighted area. (a) The raw images; (b) CLAHE; (c) LCM-CLAHE; (d) Retinex-Net; (e) Zero-DCE; and (f) our method.
Sensors 23 02498 g013
Figure 14. The contrast enhancement results of different methods. Top: raw image and the enhancement results; bottom: zoomed-in views on the highlighted area. (a) The raw images; (b) CLAHE; (c) LCM-CLAHE; (d) Retinex-Net; (e) Zero-DCE; and (f) our method.
Figure 14. The contrast enhancement results of different methods. Top: raw image and the enhancement results; bottom: zoomed-in views on the highlighted area. (a) The raw images; (b) CLAHE; (c) LCM-CLAHE; (d) Retinex-Net; (e) Zero-DCE; and (f) our method.
Sensors 23 02498 g014
Figure 15. The contrast enhancement results of different methods. Top: raw image and the enhancement results; bottom: zoomed-in views on the highlighted area. (a) The raw images; (b) CLAHE; (c) LCM-CLAHE; (d) Retinex-Net; (e) Zero-DCE; and (f) our method.
Figure 15. The contrast enhancement results of different methods. Top: raw image and the enhancement results; bottom: zoomed-in views on the highlighted area. (a) The raw images; (b) CLAHE; (c) LCM-CLAHE; (d) Retinex-Net; (e) Zero-DCE; and (f) our method.
Sensors 23 02498 g015
Figure 16. The contrast enhancement results of different methods. Top: raw image and the enhancement results; bottom: zoomed-in views on the highlighted area. (a) The raw images; (b) CLAHE; (c) LCM-CLAHE; (d) Retinex-Net; (e) Zero-DCE; and (f) our method.
Figure 16. The contrast enhancement results of different methods. Top: raw image and the enhancement results; bottom: zoomed-in views on the highlighted area. (a) The raw images; (b) CLAHE; (c) LCM-CLAHE; (d) Retinex-Net; (e) Zero-DCE; and (f) our method.
Sensors 23 02498 g016
Table 1. Image contrast of Fe (20 mm thick) and Al (1 mm thick).
Table 1. Image contrast of Fe (20 mm thick) and Al (1 mm thick).
Filter Thickness (mm)LminLmaxContrast Ratio
150 kV0449------
0.521944,475203
0.7518034,598192
116128,064174
180 kV01135------
0.5523------
0.75425------
134541,229120
Table 2. Evaluation indicators of different methods.
Table 2. Evaluation indicators of different methods.
Original ImageCLAHELCM_CLAHERetinexNetZero_DCEProposed Method
90 kVNRSS0.77910.85990.77910.76150.79690.8668
AG0.6140.8140.6140.2150.2560.815
SF4.74.724.74.394.625.85
STD51.2250.4751.2227.1744.473.17
Entropy6.55496.56776.55495.91346.83497.013
120 kVNRSS0.94570.91890.94570.91940.94570.9838
AG0.2850.3490.2850.5620.6160.623
SF2.492.532.493.984.545.4
STD40.5240.0740.5222.6837.7287.19
Entropy4.5954.59854.5954.64595.21045.7642
150 kVNRSS0.96630.91860.96630.93170.93410.9881
AG0.2080.2570.2080.3030.3870.566
SF1.962.011.964.834.515.26
STD31.2131.2331.2119.0131.3882.33
Entropy3.53943.58423.53943.95374.04674.1662
180 kVNRSS0.94280.97740.94280.89960.80860.9856
AG0.2030.2480.2030.5070.5120.578
SF1.71.771.75.125.76.03
STD25.2825.2825.2817.3226.3175.05
Entropy3.39533.57363.39533.93324.16284.8621
Table 3. Evaluation indicators of different methods.
Table 3. Evaluation indicators of different methods.
Original ImageCLAHELCM_CLAHERetinexNetZero_DCEProposed Method
Average valueNRSS0.88820.91210.90730.94790.93650.9723
AG0.3490.4060.3490.78130.79260.8786
SF1.861.9261.863.96673.68664.6766
STD24.45624.16324.52422.01330.3763.636
Entropy4.99665.0014.99964.35615.33916.2685
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, Y.; Li, K.; Ding, D.; Li, Y.; Zhao, P. X-ray Single Exposure Imaging and Image Processing of Objects with High Absorption Ratio. Sensors 2023, 23, 2498. https://doi.org/10.3390/s23052498

AMA Style

Liu Y, Li K, Ding D, Li Y, Zhao P. X-ray Single Exposure Imaging and Image Processing of Objects with High Absorption Ratio. Sensors. 2023; 23(5):2498. https://doi.org/10.3390/s23052498

Chicago/Turabian Style

Liu, Yanxiu, Kaitai Li, Dan Ding, Ye Li, and Peng Zhao. 2023. "X-ray Single Exposure Imaging and Image Processing of Objects with High Absorption Ratio" Sensors 23, no. 5: 2498. https://doi.org/10.3390/s23052498

APA Style

Liu, Y., Li, K., Ding, D., Li, Y., & Zhao, P. (2023). X-ray Single Exposure Imaging and Image Processing of Objects with High Absorption Ratio. Sensors, 23(5), 2498. https://doi.org/10.3390/s23052498

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop