Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Two-dimensional time- and space-resolved diagnostic method for an integrated implosion process

Open Access Open Access

Abstract

To precisely measure and evaluate X-ray generation and evolution in a hohlraum during an implosion process, we present a two-dimensional (2D) time- and space-resolved diagnostic method by combining a compressed ultrafast photography (CUP) system and a simplified version of space-resolving flux detector (SSRFD). Numerical experiment results showed that the reconstruction quality of the conventional CUP significantly improved owing to the addition of the external SSRFD, especially when a coded mask with a large pixel size was used in the CUP. Further, the performance of the CUP cooperation with the SSRFD was better than that of adding an external charge-coupled device or streak camera. Compared with existing ultrafast imaging techniques in laser fusion, the proposed method has a prominent advantage of measuring the 2D evolution of implosion by combining high temporal resolution of streak camera and high spatial resolution of SSRFD; moreover, it can provide guidance for designing diagnostic experiments in laser fusion research.

© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

In the study of inertial confinement fusion (ICF), the precise measurement of X-ray generation and evolution plays a crucial role in controlling physical experiment parameters and verifying simulation programs. Two principal approaches are used with lasers to drive an ICF implosion [1]. In the direct-drive approach, the laser directly irradiates the target filled with fuel, while in the indirect-drive approach, the laser is injected into the hohlraum of specific configuration and interacts with the cavity wall to excite X-ray radiation and then irradiates the target. Under the ablation of laser or X-ray, the surface material of the capsule is ablated and sprayed away, and a huge recoil pressure is generated, which compresses the fuel to the extremely high temperature and high-pressure state and achieves the ignition condition. The implosion experiment mainly focuses on the whole process of the implosion compression and the final implosion efficiency. Recording the self-emitting X-ray of the core during the implosion process with high temporal and spatial resolution is essential to evaluate the final implosion symmetry and bang time. With the development of X-ray streak and X-ray framing cameras [25] specifically designed for ICF measurements, these devices can achieve ultrahigh temporal resolution (TR) and spatial resolution (SR) of X-ray diagnosis [69]. However, the high TR of X-ray streak cameras is achieved at the expense of imaging dimension. The streak cameras used to obtain streaked radiographic image of the imploded capsule limit their imaging to only one spatial dimension. Multi-imaging X-ray streak camera can obtain ultrafast two-dimensional (2D) X-ray images [5]. However, its size of the viewing field is only 200 $\mu$m $\times$ 200 $\mu$m [10], and not suitable for imaging targets larger than the field of view due to overlapping. The X-ray framing cameras can simultaneously provide time- and space-resolving X-ray diagnostics; however, owing to its limited frame number, it cannot achieve ultrahigh TR in the implosion process. Up to 8 hot-spot images with sub-30$ps$ temporal resolution were obtained in the experiments [1113]. In 2015, a new 2D space-resolving flux detection technique was developed to measure X-ray flux inside a hohlraum using a time- and space-resolving flux detector (SRFD) [14] and successfully conducted in the China laser facility [15]. Based on the radiation fluxes measured accurately by SRFD, radiation temperature distributions within a hohlraum was first explored in Ref. [15].

In the last few decades, high-speed imaging has made significant progress, especially with the invention of high-performance image sensors based on charge-coupled device (CCD) and complementary metal-oxide semiconductor, high-speed imaging speed has been able to achieve $10^7$ frames per second (fps) [16]. Such frame rates have been unable to meet the increasing demand for ultrafast imaging. Recently, the pump-probe technique has become the dominant method for capturing transient events [17]; however, it requires that the captured events are repeatable. Capturing nonrepeating time evolution events at ultrahigh speed is essential, especially in ICF. Compressed ultrafast photography (CUP) [18], as a receiving-only ultrahigh-speed imaging technology, was proposed in 2014. Such a method realized 2D computational ultrafast imaging on the basis of streak camera and compressed sensing. Related work has been proposed and effective results have been obtained, continuously improving the CUP—from monochromatic to multispectral imaging [1922], from visible to ultraviolet light [23], from three dimensions to higher dimensions [24]; the reported frame rates have been able to achieve $7 \times 10^{13}$ fps [25] and $1.8 \times 10^{14}$ fps [26] by using different techniques. While achieving higher TR, SR is insufficient owing to the existence of encoding and decoding steps. To achieve better quality of the reconstructed dynamic image, except the original used TwIST [27] algorithm , a series of algorithm-based methods [2830] and machine learning methods [3133] have been adopted in the reconstruction process. Using different coded mask designs [34,35] and external CCD [17] also improve the reconstruction results. However, these methods are still not sufficient for the high SR diagnosis of implosion process.

In this paper, we present a 2D time- and space-resolved diagnostic method for the integrated implosion process by combining a CUP system and a simplified version of SRFD (SSRFD). Figure 1 shows the schematic diagram of the proposed method. Measurement via the CUP system is imaged through a pinhole on the left side of the hohlraum, whereas the SSRFD measurement is obtained through a pinhole on the right side of the hohlraum. Finally, the dynamic scenes of the integrated implosion process are reconstructed from these measurements by adopting a compressive sensing algorithm. To the best of our knowledge, no 2D time- and space-resolved diagnostic method with high resolution in the entire implosion process exists. Numerical experiment results show that the proposed method achieves better reconstruction quality than the conventional CUP. Compared with the CUP, the high TR of the streak camera was retained, and insufficient SR of CUP significantly improved by adding the SSRFD.

 figure: Fig. 1.

Fig. 1. Schematic of the proposed method, where FXRD is the flat-response X-ray detector.

Download Full Size | PDF

2. Method

2.1 CUP

In the conventional CUP system, dynamic scenes $I\left (x,y,t\right )$ are imaged via an objective lens and a 4f imaging system, spatially modulated by a digital micromirror device (DMD), and then measured using a streak camera. In the simulation, the pixel size of DMD needs to match with CCD. However, in the experiment, since the pixel size of DMD is not equal to that of CCD, the two need to bin pixels to match each other. The measured image $E\left (m,n\right )$ using the conventional CUP system can be mathematically expressed as follows [18]:

$$E\left(m,n\right)=TSCI\left(x,y,t\right),$$
where $C$ is the spatial-encoding operator which means the dynamic screens are encoded by the coded mask, $S$ is the temporal-shearing operator which means the dynamic screens are sheared along one spatial direction by streak camera with time, and $T$ is the spatiotemporal-integrating operator which means the temporally integrating of the sheared dynamic screens. Given the information on the coded mask, the inverse problem of Eq. (1) is solved to reconstruct the dynamic scenes from the CUP measurement. This process can be formulated as follows:
$$\mathop{\arg\min}_{I}{\left\{\frac{1}{2}\left\|E-TSCI\right\|_{2}^{2}+\lambda TV\left(I\right)\right\}},$$
where $TV\left (\right )$ denotes the total variation (TV) regularizer, The three-dimensional TV operator operates on the temporal and spatial dimensions equally and $\lambda$ is the regularization parameter. Eq. (2) can be solved by adopting a compressed sensing algorithm.

2.2 SSRFD

The concept of SRFD was first proposed in Ref. [14]. The SSRFD comprises an X-ray CCD at the periphery and a flat-response X-ray detector (FXRD) [36] at the center (Fig. 1). The FXRD is an absolutely calibrated flux detector for the X-ray energy range of 0.1–4 keV [14]. The size of the SSRFD is the same as that of CCD in the streak camera which is total 100$\times$100 pixels with pixel size of 10 $\mu$m $\times$ 10 $\mu$m. The peripheral CCD provides the time-integrated image of dynamic scenes with high SR, whereas the internal FXRD provides the X-ray flux information with TR. The diameter of internal FXRD is one-tenth that of the SSRFD. The internal FXRD can be customized and placed in an area where the dynamic scenes change the most with time. In this study, the internal FXRD was placed at the center of the SSRFD and and the temporal resolutions of FXRD were chosen range from 25$ps$ to 2$ns$.

2.3 Combination of the CUP and SSRFD

To reconstruct a 2D dynamic scene with high TR and SR, the measurement system is composed of a CUP and an SSRFD (Fig. 1). The data $E$ using the proposed method can be mathematically expressed as follows:

$$E=\left\{ \begin{aligned} E_{CUP} & =TSCI\left(x,y,t\right), \\ E_{SSFRD} & =\left\{ \begin{aligned} E\left(m,n\right) & =TI\left(x,y,t\right), & & m,n\in CCD, \\ E\left(i\right) & =\sum_{m,n}^{} \int_{\left(i-1\right)\cdot t_0}^{i\cdot t_0}I\left(x,y,t\right)dt, & & m,n\in FXRD, \end{aligned} \right. \end{aligned} \right.$$
where $i$ is the $i$-th time interval, and $t_0$ is the TR of the SSRFD. The measured data E contain two parts: the measurement of the CUP and SSRFD. The distances between the measured image and the target of CUP and SSRFD are equal to each other. The pixel size of CCDs used in CUP and SSRFD is also the same. Based on the above settings, we assume that the measured data of CUP and SSRFD correspond to the same dynamic image. Measurements from CUP and SSRFD are combined to reconstruct the dynamic screens. But they are not of equal contributions. Due to its poor temporal resolution, measurements from SSRFD can partially contribute to the reconstructed dynamic image in the backward model. Similar to the conventional CUP, the dynamic scenes can be reconstructed from the measurement by adopting a compressed sensing algorithm. In this study, the Generalized Alternating Projection (GAP) algorithm [37] with a TV regularizer is adopted to reconstruct the dynamic scenes. Eq. (2) can be solved as follow [37]:
$$\left\{ \begin{aligned} x^{(t+1)} & =\theta^{(t)}+{\Phi}^{\top}{\left(\Phi{\Phi}^{\top}\right)}^{{-}1}\left(E-\Phi\theta^{(t)}\right), \\ \theta^{(t+1)} & =TVdenoise(x^{(t+1)}),\\ \end{aligned} \right.$$
where $\Phi$ is the sensing matrix which is the whole measurement process of the proposed system. $\theta$ is an auxiliary variable. $TVdenoise()$ is the TV denoising step which can be solved by the iterative clipping algorithm [38] and a BM3D [39] denoiser is used after every 50 iterations to obtain better reconstruction results. The parameter of the contribution of SSRFD and the regularization parameter used in compressed sensing algorithm need manual adjustment for better reconstruction quality. In this simulation, they are 0.2 and 0.5, respectively.

2.4 Evaluation metrics

To quantitatively assess the quality of recovery, the peak signal-to-noise ratio (PSNR) is defined. The PSNR between the source $x$ and the restored $y$ is defined as follows:

$$PSNR\left(x,y\right)=10\log \frac{peakval^{2}}{MSE\left(x,y\right)},$$
where $MSE(x,y)$ is the mean square error between $x$ and $y$. the $peakval$ is the maximum value of the image data type. The larger the value of PSNR, the higher is the reconstruction quality. Moreover, to evaluate the similarity of the reconstructed dynamic scene, the structural similarity (SSIM) [40] is defined as follows:
$$SSIM\left(x,y\right)=\frac{\left(2\mu_x \mu_y+C_1\right)\left(2\sigma_{xy}+C_2\right)}{\left(\mu^2_x +\mu^2_y+C_1\right)\left(\sigma^2_x+\sigma^2_y+C_2\right)},$$
where $\mu _x$ is the mean of $x$, $\sigma _x$ is the variance of $x$, $\sigma _{xy}$ is the covariance of $x$ and $y$, and $C_1$, $C_2$ are constants. The value range of SSIM is $[0, 1]$. The larger the value of SSIM, the higher the SSIM.

3. Numerical experiments setup

X-ray self-emission images of hot-spot have been obtained experimentally in Ref. [12] and Ref. [13] from different laser facilities. Regardless of the inhomogeneity of the image intensity in some directions, the distributions of targets could be approximated as a Gaussian distribution. To make the simulation closer to the real data, the physical parameters of the hohlraum and capsule are derived from experiments at laser facilities in China. So, in these simulations, a 100 $\times$ 100 image with two Gaussian distributions was used as a base image, corresponding to the physical size of 1 mm $\times$ 1 mm, which means the CCD’s pixel size is 10 $\mu$m $\times$ 10 $\mu$m. We consider a dynamic scene set with 80 frames and a total time of 2 $ns$ with a TR of 25 $ps$. Similar to the result in Ref. [41], the maximum intensity of each frame with respect to time is shown in Fig. 2(a). Every pixel of each frame follows the relative change over time. The means of the Gaussian distribution along the x and y-axis are both 0 and the variances of the Gaussian distribution along the x and y-axis are both 0.03. Under the consideration of inhomogeneity in each frame, One of the Gaussian distributions with half intensity is shifted 40 pixels from the center. The correlation coefficient between the x and y-axis of each frame is slightly different ranging from −0.01 to 0.01 so that the difference in the radiation field distribution among frames is not only the intensity but also the outline. (Figs. 2(b)-2(d)). In ICF experiments, the laser interacts with the target to produce strong electromagnetic and ionizing radiation, including neutrons, gamma rays, X-rays and electromagnetic pulses. The sources of noise in experimental measurement are various and the radiation flux is usually very large, so we simply choose Gaussian noise as the noise model. In these simulations, the mean of noise was 0, and the variance of noise was 1$\%$ of the mean measured values of CUP and SSRFD, respectively.

 figure: Fig. 2.

Fig. 2. Numerical experiments setup of intensity: (a) maximum intensity of each frame; (b-d) different intensity distributions of three frames.

Download Full Size | PDF

3.1 Validation test

In the conventional CUP, the coded mask is usually a DMD displayed at a pseudorandom binary pattern. In laser fusion, DMD is unsuitable as a coded mask for X-ray. We consider using a transmission mask with a large pixel size because of manufacturing difficulties. In this test, we use a mask of 10 times CCD unit pixel, that is, a 10 $\times$ 10 mask with a pixel size of 0.1 mm $\times$ 0.1 mm. The TR of SSRFD is 100 $ps$.

3.2 Comparison tests

To illustrate the advantages of our method, we set up four groups of numerical experiments, the first was the conventional CUP (marked as CUP), the second was two CUP systems using complementary masks (marked as Two_SC), the third was CUP with an external CCD camera placed at the right of a hohlraum (marked as SC_CCD), and the last was the proposed method with different TR for the SSRFD. The contribution to reconstruction quality from the pixel size of coded mask was also investigated in three cases: 10 $\times$ 10, 20 $\times$ 20, and 100 $\times$ 100 coded masks (pixel sizes of 0.1 mm $\times$ 0.1 mm, 50 $\mu$m $\times$ 50 $\mu$m, and 10 $\mu$m $\times$ 10 $\mu$m, respectively).

4. Results and discussion

4.1 Validation test

Figure 3 shows three frames of the ground truth and the reconstructed intensity profile; the general shape and intensity of distributions are all maintained using the proposed method. Focusing on the absolute pixel-wise error in the figure, the intensities of the proposed method were closer to the ground truth. The reconstructed intensity distribution of the conventional CUP distorted owing to noise and the large pixel size of the mask, but the proposed method does not. The proposed method achieves better smoothness.

 figure: Fig. 3.

Fig. 3. Comparison of three frames: (a) simulated ground truth, reconstructed frames of the conventional CUP system and the proposed method. (b) Absolute pixel-wise error of the conventional CUP system and the proposed method, respectively.

Download Full Size | PDF

The PSNR between ground truth and reconstructed intensity distribution was calculated for the conventional CUP and the proposed method as 30.518 and 36.633, respectively. The PSNRs of all frames were calculated to evaluate the reconstruction difference during the implosion (Fig. 4(a)). The proposed method was superior to the conventional CUP at each point. Further, the SSIM of the proposed method was 0.9513, whereas that of the conventional CUP was 0.8958. Figure 4(b) shows the SSIMs of all frames, which has the same trend as Fig. 4(a). That is, the proposed method was superior to the conventional CUP. In indirectly driven laser fusion, the measurement of the radiation flux crucial in the diagnostic process. To compare the reconstruction quality of the two methods across frames, we plotted the normalized total intensity against the frame index (Fig. 4(c)). The reconstruction of the proposed method got closer to the ground truth than the conventional CUP, demonstrating a better reconstruction performance in the time domain.

 figure: Fig. 4.

Fig. 4. Comparisons of the conventional CUP with the proposed method for different frames: (a) PSNRs of different frames; (b) SSIMs of different frames. (c) Normalized total intensity of different frames.

Download Full Size | PDF

4.2 Comparison tests

To quantitatively evaluate the reconstruction quality of the four methods, we calculate PSNR and SSIM in different cases, and the results are shown in Fig. 5. The PSNRs of the first three methods were 30.518, 32.640, and 32.389, respectively. From Fig. 5(a), when the TR of SSRFD exceeded 500 $ps$, the PSNR of the proposed method was higher than that of the first three methods. The SSIMs of the first three methods were 0.8958, 0.9439, and 0.9296, respectively. From Fig. 5(b), when the TR of SSRFD exceeded 200 $ps$, the SSIM of the proposed method was higher than that of the first three methods. The above result implies that when the TR of SSRFD was more than 200 $ps$, the proposed method was better than the other three methods in terms of SSIM and PSNR. This requirement for the TR of SSRFD can be easily achieved.

 figure: Fig. 5.

Fig. 5. Comparison of different SSRFD settings: (a) PSNRs of the proposed method under different SSRFD settings; (b) SSIMs of the proposed method under different SSRFD settings.

Download Full Size | PDF

The mask size used in the above experimental results was 10 times the pixel size of CCD. To evaluate the effect of mask unit size on reconstruction quality, we tested two other mask sizes: one and five times the pixel size of CCD. We calculated the SSIMs and PSNRs of the four methods in different mask settings. The results of our method also showed the performance of the SSRFD with different TRs.

As shown in Table 1, the PSNRs and SSIMs of the reconstructed dynamic scenes improved with the decrease in the minimum pixel size of the mask. The improvements of the first three methods were larger than the last, implying that the first three methods were more dependent on the pixel size of the mask than the proposed method. With the increase in the pixel size of coded masks, the advantages of the proposed method were more. In fusion experiment diagnosis, the pixel size of coded masks is unlikely to be as fine as that of the conventional CUP system for visible light. When the mask size was small, the proposed method was better than the first three methods in terms of both SSIM and PSNR. Moreover, the reconstructed results improved by elevating the TR of the SSRFD.

Tables Icon

Table 1. PSNRs and SSIMs of the reconstructed dynamic scenes using different methods

To illustrate the relationship between the number of FXRD and the reconstruction of dynamic images, a case was studied where one FXRD was located in the center of the CCD and there were some FXRDs in the adjacent regions. In this case, we did not consider the effect of noise, and the temporal resolution of FXRD was fixed at 100$ps$. The PSNR and SSIM of the reconstructed dynamic image were calculated under different mask sizes as shown in Fig. 6. With the increased number of FXRD, the reconstructed quality of dynamic images showed a trend of decreasing at first and then increased slightly. The reconstructed results of multiple FXRD were not as good as those of single FXRD. One possible reason was that FXRD had a poor spatial resolution, and more FXRD would reduce the high spatial resolution characteristics of SSRFD. In addition, the results were affected by the size of FXRD which was not considered yet and could be investigated in future work.

 figure: Fig. 6.

Fig. 6. Comparison of different number of FXRD: (a) PSNRs of the proposed method under different number of FXRD; (b) SSIMs of the proposed method under number of FXRD.

Download Full Size | PDF

5. Conclusion

In summary, based on the SSRFD and CUP, a 2D space- and time-resolved diagnostic method for integrated implosion targets was developed. Compared with the CUP, the high TR of the streak camera was retained, and insufficient SR of CUP was significantly improved by adding an SSRFD. Situations of coded masks with different pixel size were also investigated. Further, the performance of CUP cooperation with the SSFRD is better than that of adding an external CCD or streak camera. Numerical experiments showed that when the TR of the entire system was limited by a TR of the streak camera, a smaller pixel size of coded masks or an SSFRD with higher TR significantly improved the SR of reconstruction. In future studies, the proposed method can be further combined with pinhole array imaging, which has potential applications in three-dimensional space- and time-resolved diagnosis.

Funding

National Natural Science Foundation of China (11675014, 11675157, 11805180).

Acknowledgments

The authors thank Prof. Molin Ge for the valuable discussions.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. J. D. Lindl, P. Amendt, R. L. Berger, S. G. Glendinning, S. H. Glenzer, S. W. Haan, R. L. Kauffman, O. L. Landen, and L. J. Suter, “The physics basis for ignition using indirect-drive targets on the national ignition facility,” Phys. Plasmas 11(2), 339–491 (2004). [CrossRef]  

2. F. Ze, R. Kauffman, J. Kilkenny, J. Wielwald, P. Bell, R. Hanks, J. Stewart, D. Dean, J. Bower, and R. Wallace, “A new multichannel soft x-ray framing camera for fusion experiments,” Rev. Sci. Instrum. 63(10), 5124–5126 (1992). [CrossRef]  

3. P. Gallant, P. Forget, F. Dorchies, Z. Jiang, J. Kieffer, P. Jaanimagi, J. Rebuffie, C. Goulmy, J. Pelletier, and M. Sutton, “Characterization of a subpicosecond x-ray streak camera for ultrashort laser-produced plasmas experiments,” Rev. Sci. Instrum. 71(10), 3627–3633 (2000). [CrossRef]  

4. M. Heya, S. Fujioka, H. Shiraga, N. Miyanaga, and T. Yamanaka, “Development of wide-field, multi-imaging x-ray streak camera technique with increased image-sampling arrays,” Rev. Sci. Instrum. 72(1), 755–758 (2001). [CrossRef]  

5. H. Shiraga, S. Fujioka, P. Jaanimagi, C. Stoeckl, R. Stephens, H. Nagatomo, K. Tanaka, R. Kodama, and H. Azechi, “Multi-imaging x-ray streak camera for ultrahigh-speed two-dimensional x-ray imaging of imploded core plasmas,” Rev. Sci. Instrum. 75(10), 3921–3925 (2004). [CrossRef]  

6. J. Kimbrough, P. Bell, G. Christianson, F. Lee, D. Kalantar, T. Perry, N. Sewall, and A. Wootton, “National ignition facility core x-ray streak camera,” Rev. Sci. Instrum. 72(1), 748–750 (2001). [CrossRef]  

7. J. A. Oertel, R. Aragonez, T. Archuleta, C. Barnes, L. Casper, V. Fatherley, T. Heinrichs, R. King, D. Landers, F. Lopez, P. Sanchez, G. Sandoval, L. Schrank, P. Walsh, P. Bell, M. Brown, R. Costa, J. Holder, S. Montelongo, and N. Pederson, “Gated x-ray detector for the national ignition facility,” Rev. Sci. Instrum. 77(10), 10E308 (2006). [CrossRef]  

8. J. Kimbrough, P. Bell, D. Bradley, J. Holder, D. Kalantar, A. MacPhee, and S. Telford, “Standard design for national ignition facility x-ray streak and framing cameras,” Rev. Sci. Instrum. 81(10), 10E530 (2010). [CrossRef]  

9. C. Trosseille, D. Aubert, L. Auger, S. Bazzoli, T. Beck, P. Brunel, M. Burillo, C. Chollet, J. Gazave, S. Jasmin, P. Maruenda, I. Moreau, G. Oudot, J. Raimbourg, G. Soullié, P. Stemmler, and C. Zuber, “Overview of the argos x-ray framing camera for laser megajoule,” Rev. Sci. Instrum. 85(11), 11D620 (2014). [CrossRef]  

10. H. Shiraga, H. Nagatomo, W. Theobald, A. Solodov, and M. Tabak, “Fast ignition integrated experiments and high-gain point design,” Nucl. Fusion 54(5), 054005 (2014). [CrossRef]  

11. K. Engelhorn, T. Hilsabeck, J. Kilkenny, D. Morris, T. Chung, A. Dymoke-Bradshaw, J. Hares, P. Bell, D. Bradley, A. Carpenter, M. Dayton, S. R. Nagel, L. Claus, J. Porter, G. Rochau, M. Sanchez, S. Ivancic, C. Sorce, and W. Theobald, “Sub-nanosecond single line-of-sight (slos) x-ray imagers,” Rev. Sci. Instrum. 89(10), 10G123 (2018). [CrossRef]  

12. W. Theobald, C. Sorce, M. Bedzyk, S. Ivancic, F. Marshall, C. Stoeckl, R. Shah, M. Lawrie, S. Regan, T. Sangster, E. M. Campbell, T. J. Hilsabeck, K. Englehorn, J. D. Kilkenny, D. Morris, T. M. Chung, J. D. Hares, A. K. L. Dymoke-Bradshaw, P. Bell, J. Celeste, A. C. Carpenter, M. Dayton, D. K. Bradley, M. C. Jackson, L. Pickworth, S. R. Nagel, G. Rochau, J. Porter, M. Sanchez, L. Claus, G. Robertson, and Q. Looker, “The single-line-of-sight, time-resolved x-ray imager diagnostic on omega,” Rev. Sci. Instrum. 89(10), 10G117 (2018). [CrossRef]  

13. Y. Pu, T. Huang, F. Ge, S. Zou, F. Wang, J. Yang, S. Jiang, and Y. Ding, “First integrated implosion experiments on the sg-iii laser facility,” Plasma Phys. Controlled Fusion 60(8), 085017 (2018). [CrossRef]  

14. K. Ren, S. Liu, L. Hou, H. Du, G. Ren, W. Huo, L. Jing, Y. Zhao, Z. Yang, M. Wei, K. Deng, L. Yao, Z. Li, D. Yang, C. Zhang, J. Yan, G. Yang, S. Li, S. Jiang, Y. Ding, J. Liu, and K. Lan, “Direct measurement of x-ray flux for a pre-specified highly-resolved region in hohlraum,” Opt. Express 23(19), A1072–A1080 (2015). [CrossRef]  

15. K. Ren, S. Liu, X. Xie, H. Du, L. Hou, L. Jing, D. Yang, Y. Zhao, J. Yan, Z. Yang, Z. Li, J. Dong, G. Yang, S. Li, Z. Cao, K. Lan, W. Huo, J. Liu, G. Ren, Y. Ding, and S. Jiang, “First exploration of radiation temperatures of the laser spot, re-emitting wall and entire hohlraum drive source,” Sci. Rep. 9(1), 5050 (2019). [CrossRef]  

16. Y. Kondo, K. Takubo, H. Tominaga, R. Hirose, N. Tokuoka, Y. Kawaguchi, Y. Takaie, A. Ozaki, S. Nakaya, F. Yano, and T. Daigen, “Development of hypervision hpv-x high-speed video camera,” Shimadzu Rev. 69, 285–291 (2012).

17. L. Zhu, Y. Chen, J. Liang, Q. Xu, L. Gao, C. Ma, and L. V. Wang, “Space-and intensity-constrained reconstruction for compressed ultrafast photography,” Optica 3(7), 694–697 (2016). [CrossRef]  

18. L. Gao, J. Liang, C. Li, and L. V. Wang, “Single-shot compressed ultrafast photography at one hundred billion frames per second,” Nature 516(7529), 74–77 (2014). [CrossRef]  

19. J. Liang, L. Zhu, and L. V. Wang, “Single-shot real-time femtosecond imaging of temporal focusing,” Light: Sci. Appl. 7(1), 42 (2018). [CrossRef]  

20. Y. Lu, T. T. Wong, F. Chen, and L. Wang, “Compressed ultrafast spectral-temporal photography,” Phys. Rev. Lett. 122(19), 193904 (2019). [CrossRef]  

21. J. Yao, D. Qi, C. Yang, F. Cao, Y. He, P. Ding, C. Jin, Y. Yao, T. Jia, Z. Sun, and S. Zhang, “Multichannel-coupled compressed ultrafast photography,” J. Opt. 22(8), 085701 (2020). [CrossRef]  

22. C. Yang, F. Cao, D. Qi, Y. He, P. Ding, J. Yao, T. Jia, Z. Sun, and S. Zhang, “Hyperspectrally compressed ultrafast photography,” Phys. Rev. Lett. 124(2), 023902 (2020). [CrossRef]  

23. Y. Lai, Y. Xue, C.-Y. Côté, X. Liu, A. Laramée, N. Jaouen, F. Légaré, L. Tian, and J. Liang, “Single-shot ultraviolet compressed ultrafast photography,” Laser Photonics Rev. 14(10), 2000122 (2020). [CrossRef]  

24. J. Liang, P. Wang, L. Zhu, and L. V. Wang, “Single-shot stereo-polarimetric compressed ultrafast photography for light-speed observation of high-dimensional optical transients with picosecond resolution,” Nat. Commun. 11(1), 5252 (2020). [CrossRef]  

25. P. Wang, J. Liang, and L. V. Wang, “Single-shot ultrafast imaging attaining 70 trillion frames per second,” Nat. Commun. 11(1), 2091 (2020). [CrossRef]  

26. D. Qi, F. Cao, S. Xu, Y. Yao, Y. He, J. Yao, P. Ding, C. Jin, L. Deng, T. Jia, J. Liang, Z. Sun, and S. Zhang, “100-trillion-frame-per-second single-shot compressed ultrafast photography via molecular alignment,” Phys. Rev. Appl. 15(2), 024051 (2021). [CrossRef]  

27. J. M. Bioucas-Dias and M. A. Figueiredo, “A new twist: Two-step iterative shrinkage/thresholding algorithms for image restoration,” IEEE Trans. Image Process. 16(12), 2992–3004 (2007). [CrossRef]  

28. C. Yang, D. Qi, F. Cao, Y. He, X. Wang, W. Wen, J. Tian, T. Jia, Z. Sun, and S. Zhang, “Improving the image reconstruction quality of compressed ultrafast photography via an augmented lagrangian algorithm,” J. Opt. 21(3), 035703 (2019). [CrossRef]  

29. M. Tsumori, S. Nagai, R. Harakawa, T. Sasaki, and M. Iwahashi, “Contour-adaptive image restoration based on compressed sensing for ultrafast phenomena,” Appl. Phys. Lett. 117, 044102 (2020). [CrossRef]  

30. J. Yao, D. Qi, Y. Yao, F. Cao, Y. He, P. Ding, C. Jin, T. Jia, J. Liang, L. Deng, Z. Sun, and S. Zhang, “Total variation and block-matching 3d filtering-based image reconstruction for single-shot compressed ultrafast photography,” Opt. Lasers Eng. 139, 106475 (2021). [CrossRef]  

31. Y. Ma, X. Feng, and L. Gao, “Deep-learning-based image reconstruction for compressed ultrafast photography,” Opt. Lett. 45(16), 4400–4403 (2020). [CrossRef]  

32. A. Zhang, J. Wu, J. Suo, L. Fang, H. Qiao, D. D.-U. Li, S. Zhang, J. Fan, D. Qi, Q. Dai, and C. Pei, “Single-shot compressed ultrafast photography based on u-net network,” Opt. Express 28(26), 39299–39310 (2020). [CrossRef]  

33. C. Yang, Y. Yao, C. Jin, D. Qi, F. Cao, Y. He, J. Yao, P. Ding, L. Gao, T. Jia, J. Liang, Z. Sun, and S. Zhang, “High-fidelity image reconstruction for compressed ultrafast photography via an augmented-lagrangian and deep-learning hybrid algorithm,” Photonics Res. 9(2), B30–B37 (2021). [CrossRef]  

34. C. Yang, D. Qi, J. Liang, X. Wang, F. Cao, Y. He, X. Ouyang, B. Zhu, W. Wen, T. Jia, J. Tian, L. Gao, Z. Sun, S. Zhang, and L. V. Wang, “Compressed ultrafast photography by multi-encoding imaging,” Laser Phys. Lett. 15(11), 116202 (2018). [CrossRef]  

35. C. Yang, D. Qi, X. Wang, F. Cao, Y. He, W. Wen, T. Jia, J. Tian, Z. Sun, L. Gao, S. Zhang, and L. V. Wang, “Optimizing codes for compressed ultrafast photography by the genetic algorithm,” Optica 5(2), 147–151 (2018). [CrossRef]  

36. Z. Li, X. Jiang, S. Liu, T. Huang, J. Zheng, J. Yang, S. Li, L. Guo, X. Zhao, H. Du, T. Song, R. Yi, Y. Liu, S. Jiang, and Y. Ding, “A novel flat-response x-ray detector in the photon energy range of 0.1–4 kev,” Rev. Sci. Instrum. 81(7), 073504 (2010). [CrossRef]  

37. X. Liao, H. Li, and L. Carin, “Generalized alternating projection for weighted-ℓ2, 1 minimization with applications to model-based compressive sensing,” SIAM J. Imaging Sci. 7(2), 797–823 (2014). [CrossRef]  

38. A. Beck and M. Teboulle, “Fast gradient-based algorithms for constrained total variation image denoising and deblurring problems,” IEEE Trans. Image Process. 18(11), 2419–2434 (2009). [CrossRef]  

39. K. Dabov, A. Foi, V. Katkovnik, and K. Egiazarian, “Image denoising by sparse 3-d transform-domain collaborative filtering,” IEEE Trans. Image Process. 16(8), 2080–2095 (2007). [CrossRef]  

40. Z. Wang, A. C. Bovik, H. R. Sheikh, and E. P. Simoncelli, “Image quality assessment: from error visibility to structural similarity,” IEEE Trans. Image Process. 13(4), 600–612 (2004). [CrossRef]  

41. Y. Pu, X. Luo, L. Zhang, C. Sun, Z. Hu, G. Shen, X. Wang, Q. Tang, Z. Yuan, F. Wang, D. Yang, J. Yang, S. Jiang, Y. Ding, and J. Wang, “Chunk mixing implosion experiments using deuterated foam capsules with gold dopant,” Phys. Rev. E 102(2), 023204 (2020). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Fig. 1.
Fig. 1. Schematic of the proposed method, where FXRD is the flat-response X-ray detector.
Fig. 2.
Fig. 2. Numerical experiments setup of intensity: (a) maximum intensity of each frame; (b-d) different intensity distributions of three frames.
Fig. 3.
Fig. 3. Comparison of three frames: (a) simulated ground truth, reconstructed frames of the conventional CUP system and the proposed method. (b) Absolute pixel-wise error of the conventional CUP system and the proposed method, respectively.
Fig. 4.
Fig. 4. Comparisons of the conventional CUP with the proposed method for different frames: (a) PSNRs of different frames; (b) SSIMs of different frames. (c) Normalized total intensity of different frames.
Fig. 5.
Fig. 5. Comparison of different SSRFD settings: (a) PSNRs of the proposed method under different SSRFD settings; (b) SSIMs of the proposed method under different SSRFD settings.
Fig. 6.
Fig. 6. Comparison of different number of FXRD: (a) PSNRs of the proposed method under different number of FXRD; (b) SSIMs of the proposed method under number of FXRD.

Tables (1)

Tables Icon

Table 1. PSNRs and SSIMs of the reconstructed dynamic scenes using different methods

Equations (6)

Equations on this page are rendered with MathJax. Learn more.

E ( m , n ) = T S C I ( x , y , t ) ,
arg min I { 1 2 E T S C I 2 2 + λ T V ( I ) } ,
E = { E C U P = T S C I ( x , y , t ) , E S S F R D = { E ( m , n ) = T I ( x , y , t ) , m , n C C D , E ( i ) = m , n ( i 1 ) t 0 i t 0 I ( x , y , t ) d t , m , n F X R D ,
{ x ( t + 1 ) = θ ( t ) + Φ ( Φ Φ ) 1 ( E Φ θ ( t ) ) , θ ( t + 1 ) = T V d e n o i s e ( x ( t + 1 ) ) ,
P S N R ( x , y ) = 10 log p e a k v a l 2 M S E ( x , y ) ,
S S I M ( x , y ) = ( 2 μ x μ y + C 1 ) ( 2 σ x y + C 2 ) ( μ x 2 + μ y 2 + C 1 ) ( σ x 2 + σ y 2 + C 2 ) ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.