Compressed sensing has been proposed for several years.1 The technique has been greatly improved and widely used for imaging systems,23.–4 such as the single-pixel camera by Duarte5 and the complicated three-dimensional imaging system. The imaging systems mentioned above fully reflect the advantage of compressed sensing for faster imaging processing time without complicated mechanical scanning structure.
When a system of image reconstruction is based on compressed sensing, the system will inevitably produce various types of errors affecting the image quality. The nonlinear response from the detector is one of the main contributors to the errors.67.–8 To correct the error due to the nonlinear response for compressed sensing, Tao et al. have done a series of studies with the following results:9,10
Assume the measured value can be expressed as
The reconstruction error can be written as9
When the detector has a nonlinear response error (, where is a positive constant), the measured data can be expressed as
The expression below will be approved in this paper11,12 When the detector presents nonlinearity, which is true for any real system, the image quality is impacted. Based on the above considerations, in order to improve the reconstructed image quality and the robustness of the system, a new algorithm for detector nonlinear response compensation is proposed in this paper for compressed sensing. This new algorithm can reduce the response errors and increase the efficiency of the imaging system noticeably based on the theoretical analysis13 and numerical simulation.14
Imaging System Based on Compressed Sensing
The imaging system used for discussion is the single-pixel camera, same as the one described in Ref. 5. The single-pixel camera imaging system is one of the representative imaging systems based on compressed sensing.3 The system block diagram is shown in Fig. 1.
As in Fig. 1, the reflected optical light modulated by the DMD is received by the photo-detector. The measured data from the detector are fed into the computer. If the photo-detector has a nonlinear response, the measured data will be nonlinear, thus the reconstruction result will be impacted. So a nonlinearity compensation algorithm can be used for image reconstruction to improve the quality.
Algorithm of Detector Nonlinear Response Compensation
This section outlines the algorithm in theory in compensating the nonlinear response from the detector for imaging system based on compressed sensing. It also analyzes the reconstruction error affected by the nonlinear response.
Nonlinear Response of the Detector
In the theory of compressed sensing,5 the measured values can be expressed as Eq. (1). In an experiment, the DMD pattern should be changed by M times to obtain different measured values. Giving and the measurement matrix , the reconstructed signal is obtained by
When the detector response is linear and is the response factor, from Eq. (8) it can be derived that9) as
It can be seen from the Eq. (10) that when the detector response is linear, the measured data are a linear transformation. By removing the corresponding response factor, the original measured data are recovered and there is no nonlinear error in the image reconstruction.
When the detector response is nonlinear, note each response coefficient of the measured data as , respectively (excluding ), then the measured data can be expressed as11), by extracting the smallest response factor (assume is the smallest one), it is observed that the measured data are a nonlinear transformation due to the nonlinear response of the detector. When the measured data are changed, it also leads to the decrease of reconstruction efficiency.
Equation (7) can be expressed as
For each measurement12), is a positive constant.
Assuming , then Eq. (13) can be written as13) can be derived. When the measured data are obtained, can be expressed as Fig. 2.
Generally speaking, (as shown in Fig. 2) and the probability of is very small. For example, when the target picture is all black, which means the pixel values are all 0, then . Therefore such a scenario is excluded in discussion. When , , it is evident that & .
Rewrite Eq. (15) as16) becomes
Assume as the original signal. From above, it can be derived that19) that the nonlinear error in the measured data can directly lead to the reconstruction error. When the effect of the nonlinear response becomes noticeable, which means is large, the construction error rate increases, thereby degrading image quality. Therefore, it is very important to compensate the nonlinear response for the detector, which is this paper’s focus.
Let , , then , which can be further expressed as . Therefore,
In the above, is the reconstruction error due to nonlinear response from the detector. As shown in Figs 3Fig. 4–5, the detector nonlinear response apparently affects the quality of the reconstructed images, so needs to be minimized.
Algorithm to Compensate Nonlinear Response from the Detector
It is well known that the photo-detector response model function can be briefly expressed as
In Eq. (23), represents the accuracy error of the linear measurement. is part of the imaging error. represents the incident light intensity. is the detector response function.
Suppose the linear detection range of the detector is . When or , it can be seen that the measured data are in a nonlinear state. Therefore, compensation on the nonlinearity becomes necessary.
The process of the compensation algorithm is as shown in Fig. 6.
Figure 6 depicts detector response curve, in which is the detector response factor. If the detector response factor does not change for a given region, then this region can be treated as a linear detection region. Otherwise, it is nonlinear.
Measured data associated with the nonlinear region in Fig. 6 is rejected in the algorithm. The data corresponding to the linear region is described below.
Assume is the total number of the measured data points while is the number of the measured data points in the nonlinear response region. Define an operation
The dimension of the measured data is after the rejection.
Since both the measured data and the Gaussian random measurement matrix are needed for signal reconstruction and each measured data has its corresponding Gaussian random measurement matrix, therefore after changing the measured data, the Gaussian random measurement matrix should also have the appropriate operation before the final recovery.
Let the Gaussian random measurement matrix be . In the measured data matrix , corresponds to the Gaussian random measurement matrix . When is rejected, the Gaussian random measurement matrix should also be changed accordingly as
Figure 7 is the schematic diagram demonstrating the compensation process of the measured data and the Gaussian random measurement matrix.
In fact, there is no such ideal linear detection. The linear detection region shown in the Fig. 6 is actually an approximation. The response factor also changes in this region but does not change significantly. Such a region is called a near-linear region. The next section is the process to compensate the measured data in the near-linear region in order to get the reconstructed image with improved equality.
As shown in Fig. 8 the near-linear region is divided into equal parts, where is a positive integer. The boundary of the region is noted as and . is the largest measured data in the near-linear region while is the smallest. The interval of the adjacent parts is . The average response factor of each divided part is (). The measured data are assigned for each divided part according to the numerical size and the number of the measured value points that each part has is determined. The part with the largest number of data points is the reference region. The average response factor of the reference region will be the reference response factor and the measured data of the other regions are processed as below
In Eq. (26), is the reference response factor; is the average response factor of the region ; is the measured data in the region , and is the measured data which have been compensated by the response reference factor. Figure 9 illustrates compensation flow for the near-linear region.
Finally, an dimensional measured data matrix is obtained with the improved recovery image for compressed sensing. The overall block diagram of the algorithm proposed in this paper is illustrated in Fig. 10. In short, when obtaining the measured data, data corresponding to the nonlinear region are rejected and the corresponding Gaussian random measurement matrix is modified. Data in the near-linear region are used for processing based on the algorithm proposed here to get the reconstructed image with improved quality for compressed sensing.
Theoretical Analysis of the Algorithm on Error Reduction
Assume is the measured data under the ideal linear detection; is the measured data after processing, and is the measured data that have not been processed. Correspondingly, is the response factor of ; is the reference response factor; and is the average response factor of the region as defined in the previous sections.
The ideal measured data that are linear can be expressed as
The measured data after the compensation can be shown as Eq. (26).
If it can be proved that
Using the above expressions, it can be derived that
The following description consists of two main parts: Part A and Part B. It is outlined previously that the near-linear response curve has a stationary nonlinear factor rate . Generally speaking, we can also assume the nonlinear factor will not change too much in the near-linear region. So in this section it is assumed .The stationary nonlinear factor rate can be written as
is the number of the measured data points in the nonlinear detection region and is the number of the whole measured data points. The nonlinear factor rate shows the changing rate of the nonlinear factor. The larger the nonlinear factor rate is, the more nonlinear the measured data will be.
Figure 11 shows the scenario of . When , for the majority of the measured data, it is a valid assumption that
Please note that the measured data farther from the reference region are very effective when applying the compensation. There are few measured data points that are very close to the reference region. For those data, the compensation effectiveness may not be good, but it will not impact the overall effectiveness of the compensation.
From Eq. (35) it can be obtained that
With the above assumption and Fig. 11, it can be easily derived that for the majority of the measured data, below is a valid assumption
From Eq. (37)
The probability of is small, and it has very little effect on the reconstruction when this event happens. So the case of is not considered here.
From the above discussion, when there is a nonlinear response error , where , is the error of the nonlinear response and is the error of the near-linear response, the measured data can be expressed as
Assume , , after the compensation, can be ignored and is also reduced. Therefore,4)
After the compensation, can be ignored and is reduced, which means the reconstruction error has been minimized with better image quality.
Definition of Parameters
The sampling rate is defined as15 and 16:
The target images to be simulated are the Picture R, Picture Lena, and Picture Camera as shown in Fig. 12. The resolution of Picture R is , and the resolution of Picture Lena and Picture Camera are both .
Figures 3Fig. 4–5 give the comparison between the reconstructed images with and without using the compensation algorithm. The sampling rate is 30%. The nonlinear rate of the five different comparisons is 10%, 20%, 30%, 40%, and 60%, from the left to the right.
The simulation for Picture R is shown in Fig. 3. Figure 3(a)–3(e) shows the reconstruction results without the compensation algorithm. Figure 3(f)–3(j) shows the reconstruction results with the compensation algorithm.
In simulation, the nonlinear factor rate used is 0.025, and the near-linear region of the measured data is divided into six parts. It can be seen that the compensation algorithm significantly improves the quality of image reconstruction.
Figure 4 is the simulation of Picture Lena. In the simulation, the nonlinear factor rate is 0.0025. It is found that the more complex the pictures are, the worse the reconstructed pictures will be due to the nonlinear response. In order to compare the reconstructed pictures with and without using the compensation algorithm, the nonlinear factor rate is reduced to 0.0025. The near-linear region of the measured data is divided into six parts. As shown in Fig. 4, the compensation significantly improves the image quality.
Similar assumption is used for the simulation of Picture Camera and the same observation in improving the image quality is obtained. For simplicity, it will not be repeated here.
From the above analysis, it can be seen that after using the nonlinear compensation algorithm, the PSNR ratio and the reconstruction rate of the images have been significantly improved to mitigate the impact of the nonlinearity from the detector. The algorithm reduces the error caused by the detector nonlinearity and improves the imaging efficiency of the system as well.
In the imaging system based on compressed sensing, the detector nonlinear response incurs reconstruction error and degrades image quality and imaging efficiency. In order to resolve the issue, a compensation algorithm is proposed in this paper. Theoretical analysis and simulation show that the proposed algorithm can effectively mitigate the impact caused by the nonlinear response from the detector. The reconstruction efficiency and PSNR ratio have also been significantly improved after the compensation.
It is unknown whether rejecting the measured data will succeed for the situation of too many measured data in the nonlinear region. This can be future work for continuous improvement.
This project is sponsored by the National Natural Science Foundation of China (NSFC) (Grant Nos. 61101196, 61271332, and 61177091), Weaponry Pre-research Project (Grant No. 40405080401), and National Postdoctoral Foundation (Grant No. 2012M521085).
Y. TsaigD. L. Donoho, “Extensions of compressed sensing,” Technical Report, Department of Statistics, Stanford University (2004).Google Scholar
L. Liet al., “Gated viewing laser imaging with compressive sensing,” Opt. Soc. Am. 51(14), 2706–2712 (2012).Google Scholar
D. Keminget al., “Photon-counting imaging system based on compressive sensing,” Infrared Laser Eng. 41(2), 363–368 (2012).Google Scholar
M. F. Duarteet al., “Single-pixel imaging via compressive sampling,” IEEE Signal Process. Mag. 25(2), 83–91 (2008).Google Scholar
J. N. Laskaet al., “Democracy in action: Quantization, saturation, and compressive sensing,” Applied and Computational Harmonic Anal. 31(3), 429–443 (2011).ACOHE91063-5203http://dx.doi.org/10.1016/j.acha.2011.02.002Google Scholar
T. BlumensathM. E. Davies, “Iterative hard thresholding for compressed sensing,” Appl. Comput. Harmonic Anal. 27(3), 265–274 (2009).Google Scholar
H. T. KungT.-H. LinD. Vlah, “Identifying bad measurements in compressive sensing,” in Proceedings of IEEE Conference on Computer Communications Workshops (INFOCOM WKSHPS), Shanghai, China, pp. 1054–1059 (2011).Google Scholar
E. CandesJ. RombergT. Tao, “Stable signal recovery from incomplete and inaccurate measurements,” Commun. Pure and Appl. Math. 59(8), 1207–1223 (2006).CPMAMV0010-3640http://dx.doi.org/10.1002/(ISSN)1097-0312Google Scholar
E. CandesJ. RombergT. Tao, “Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information,” IEEE Trans. Infom. Theory 52(2), 489–509 (2006).IETTAW0018-9448http://dx.doi.org/10.1109/TIT.2005.862083Google Scholar
W. Lu, “Compressed sensing and sparse signal processing,” Technical Report, Department of Electrical & Computer Engineering, University of Victoria, Canada (2004).Google Scholar
Y. ZhuL. RaoT. Yan, The Analysis, and Calculation of Matrix, 1st ed., Chinese National Defence Industry Press, Beijing, China (2010).Google Scholar
Y. WuP. YeI. O. Mirza, “Experimental demonstration of an optical sectioning compressive sensing microscope (CSM),” Opt. Express 18(24), 24565–24578 (2010).Google Scholar
C. G. LeeI. MoonB. Javidi, “Photon-counting three-dimensional integral imaging with compression of elemental images,” Opt. Soc. Am. 29(6), 854–860 (2012).JOSAAH0030-3941http://dx.doi.org/10.1364/JOSAA.29.000854Google Scholar
Y. LiuR. Zhao, “Matching pursuit algorithm for signal reconstruction based on compressive sensing,” Master’s Thesis Full-text Database (CMFD), Beijing Jiaotong University, China (2010).Google Scholar
Jiayan Zhuang is pursuing his PhD in engineering at the School of Electronic and Optical Engineering at Nanjing University of Science and Technology. His research interests include nonlinear imaging processing, optical applications of compressed sensing, and biomedical imaging.
Qian Chen is a professor at the School of Electronic and Optical Engineering at Nanjing University of Science and Technology. His research interests include photoelectric detection, image processing, and photoelectric signal processing.
Weiji He is a lecturer of the Jiangsu Key Lab of Spectral Imaging and Intelligence Sense (SIIS) at Nanjing University of Science and Technology. His research interests include the mechanism and application of photoelectric detection.