Fast compressive measurements acquisition using optimized binary sensing matrices for low-light-level imaging

: Compressive measurements beneﬁt low-light-level imaging ( L 3 -imaging) due to the signiﬁcantly improved measurement signal-to-noise ratio (SNR). However, as with other compressive imaging (CI) systems, compressive L 3 -imaging is slow. To accelerate the data acquisition, we develop an algorithm to compute the optimal binary sensing matrix that can minimize the image reconstruction error. First, we make use of the measurement SNR and the reconstruction mean square error (MSE) to deﬁne the optimal gray-value sensing matrix. Then, we construct an equality-constrained optimization problem to solve for a binary sensing matrix. From several experimental results, we show that the latter delivers a similar reconstruction performance as the former, while having a smaller dynamic range requirement to system sensors.


Introduction
The idea of compressive imaging (CI) has been used when a large amount of data need to be acquired [1], while most of the collected data will be discarded [2]. The method is also useful when high-resolution sensor arrays are not available due to limitations of manufacturing technology, for example, in LIDAR 3D imaging [3,4] or terahertz imaging [5]. CI is also suitable for applications where the data acquisition device is expensive, such as low-light-level imaging (L 3 -imaging) [6][7][8]. Using CI for L 3 -imaging takes advantage of better measurement signalto-noise ratio (SNR). In such a system, the measurements are linear combinations of multiple pixels, consequently the measurement power is much larger than that of a single pixel.
Nevertheless, inherent challenges in CI make it difficult for implementation in many fields for real-world applications. One difficulty is its sensitivity to system calibration [9]. System calibration error can also be amplified by a system's data processing. Another difficulty is the size of the signal space [9]. This is a challenge for computational memory allocation and computation power. Besides these difficulties, CI system's speed also limits its applications, especially for real-time imaging. Several groups, including us, have investigated the use of GPU or FPGA to accelerate the data processing [8,[10][11][12][13][14]. However, few have considered accelerating the data acquisition. In this paper, we develop a methodology to compute the optimized binary sensing matrices to improve the speed in the system's measurement collection.
In compressive imaging, sensing matrix is always a crucial factor. Initially, random sensing matrices, such as matrices with elements selected from Bernoulli or Gaussian distribution, or matrices with row vectors randomly selected from a Fourier transformation matrix, have been studied. A CI system using these matrices can successfully reconstruct sparse signals, as they satisfy the restricted isometry property (RIP) [15,16]. However, as noise increases, the performances degrades fast [17,18]. Even without noise, given a sparse representation dictionary, searching an optimal sensing matrix is an open question for CI. Therefore, several research groups work on sensing matrix design for different applications [19,20]. A major trend focuses on the mutual coherence of the multiplication of the sensing matrix and the sparse representation dictionary [21][22][23][24]. When this coherence is small, a wide set of signals can be reconstructed successfully using the sensing matrix. Different from using mutual coherence, by analyzing the oracle estimator MSE as Gaussian noise present in system measurements, Chen et al. [25,26] find that a good sensing matrix ought to be close to a Parseval tight frame. Based on the work in literature [21,26], Yang et al. [27] further optimize the sensing matrix combining with a matrix decomposition, thus to have the optimal reconstruction and the optimal mutual coherence performance simultaneously. Different from the works discussed above, in which the sparse representation dictionary is predefined, Duarte-Carvajalino and Sapiro [22] design the dictionary and the sensing matrix jointed from a set of training images to obtain better reconstructions. Besides the mutual coherence, several other groups use information-theoretic framework to design a sensing matrix [28][29][30][31][32].
Note that the design methods discussed above are for gray value sensing matrices, which has their limitation in CI system implementation. In many CI systems investigated in literature, a DMD is the spatial light modulator for measurement collection. A DMD consists of thousands of micromirrors. These mirrors flip between two angles at a predefined frequency. During one measurement collection time T 0 , if the mirrors use half of the time to reflect light to output, then the object light is modulated by a factor 0.5. Therefore, light modulation is based on time division in DMD. Due to this modulation methodology, for a fixed DMD working frequency f , collecting a measurement using a binary projection needs one period 1/ f . However, for gray values, such as values which require 8 bit binary numbers for representation, DMD needs 256 periods to realize the light modulation. Therefore, the measurement needs much longer time( 1/ f ). This slows down a CI system's measurement collection and imaging speed. Besides this major disadvantage, using gray value projections, CI systems require more complicated system calibration process, thus potentially have larger systematical error compared to using binary projections. Because of these issues, we study binary sensing matrix design in this work.
In CI, we cannot design a sensing matrix without considering a reconstruction method. In the classical definition of CI, an object is reconstructed by solving an optimization problem leveraging object sparsity. Besides such nonlinear methods, in our previous works we also investigated the linear reconstruction Wiener operator [18,33]. The operator minimizes reconstruction mean square error (MSE). It can be written as where R x is the object autocorrelation matrix, F is the sensing matrix, and σ 2 is the noise variance. This operator does not use object sparsity for reconstruction. Its performance is slightly worse than nonlinear methods [34,35] for sparse signals when noise is low. However, the performance difference is small [18]. When random features are collected for objects of size (80 × 80) with dynamic range [0, 255] and sparsity K = 400 in the wavelet domain, the reconstruction RMSE (root mean square error) values are 9.21 and 6.54 using Wiener operator and a nonlinear Nowak method [36], respectively [18]. The noise standard deviation is 0.21. Same level RMSE difference was observed using PCA (Principal Component Analysis) features. Although Wiener operator does not use object sparsity, it is defined based on an object's second order statistical parameter, the object's autocorrelation matrix R x . This statistical parameter is not used in the nonlinear reconstruction methods of CI. Additionally, because noise is directly used to define Wiener operator, the operator is more tolerant to noise. Thus, when noise is high or object light level is low, to obtain a better reconstruction, object sparsity is not the only object prior which can be used for reconstruction. Wiener operator presents even better reconstruction in terms of smaller RMSE. Besides these advantages, more importantly to the work discussed in this paper, the computational speed using the Wiener operator is much faster than nonlinear iterative methods, especially with GPU implementation [8,11,12]. Therefore we use the Wiener operator for sensing matrix design.
Note that, Wiener-type operators recently have been used in geophysics, and medical imaging specifically in diffuse optical tomographic imaging [37][38][39]. In a near-surface geophysical problem, observed data is emulated using a linear layered-earth model by representing the data as a multiplication between a matrix and a model vector [37]. To solve the geophysical problem is to understand how well the calculated data match with the observed data, and then select those well predicted data. A Wiener-type matrix, data-resolution matrix is used to do the selection. Then another Wiener-type matrix, model-resolution matrix is used to evaluate how well the layered-earth model works with the selected data. In diffuse optical tomographic imaging, the forward problem is different from CI or geophysical problem [38,39]. However, the natural logarithm of the amplitude of experimental data can still be written as a function of biological tissue's absorption coefficients. Using the Jacobian of this function and a regularization parameter, data-resolution matrix and model-resolution matrix for this problem can be defined. Then based on data-resolution matrix, independent measurements are chosen to reduce data collection time without compromising reconstructed image quality much [38]. Model-resolution matrix is used with basic pursuit deconvolution algorithm to improve reconstructed image quality [39].
With the focuses on binary sensing matrices and Wiener operator for reconstruction, the paper is organized as follows. Before studying binary projection design, we first use maximizing measurement SNR and minimizing reconstruction MSE as the criteria to search the optimal gray value sensing matrixes in Section 2. Then based on the solutions, we design binary matrices in Section 3. Experimental results are presented in Section 4 to demonstrate the idea. Conclusions are drawn in Section 5.

Sensing matrix design based on maximizing measurement SNR
We first define the mathematic model for CI. In compressive imaging, system measurements can be represented as where y(M × 1) is the measurement vector, F = f 1 f 2 · · · f M T is the sensing matrix of size (M × N), x(N × 1) and n(M × 1) are the object and noise vectors, respectively. The noise is assumed as additive white Gaussian noise N (0, σ 2 ). Then the CI system measurement SNR becomes In Eq. (2), the object autocorrelation matrix R x = E{ x x T 2 } can be decomposed as QDQ T , where Q is an orthonormal matrix, and D is diagonal. The i th column of Q, q i , is the eigenvector of matrix R. The diagonal elements of D are sorted in descending order, such as d 1 ≥ d 2 ≥ · · · ≥ d N . d i is the eigenvalue of R corresponding to q i . Note for natural images, all d i are positive .
To obtain the optimal projection in terms of maximizing SNR, we solve the following optimization problem, Using the decomposition of R x , we rewrite Tr{FR x F T } as Tr{(FQ)D(FQ) T }. If we redefine FQ as H, then problem P 1 becomes The matrix H can be further written as H = H 1 H 2 , where H 1 and H 2 have sizes (M × M) and (M × (N − M)). One special case for P 1 or P 1a is H 1 H T 1 = I and H 2 = 0. The cost function for this case is then this special case can be written as We modify F slightly as with 0 ≤ α ≤ 1, then the cost function of P 1 becomes Thus, we can observe that as F is driven away from the space spanned by vectors q 1 q 2 · · · q M , Tr{FR x F T } will be reduced. Therefore, d i is the largest cost function value of P 1 , and its solution is

Sensing matrix design based on minimizing the reconstruction MSE
In compressive imaging, the final goal is to obtain high quality object reconstructions. Instead of improving the intermediate measurement SNR, in this subsection we study the matrix design to minimize the reconstruction error. As discussed in Sec. 1, we use Wiener operator for object reconstruction. The Wiener operator is defined as Using W, the object is reconstructed as x est = W y. Then the reconstruction MSE becomes [40] The first term in Eq. (7) is independent of F. To minimize the reconstruction error MSE CS is equivalent to maximizing the second term. Thus, we define the design problem as Due to the term (FR x F T + σ 2 I) −1 , solving P 2 with a close form solution, such as using KKT theory [41], is hard. Note that, for then we have Thus, the maximizer for P 2 is F = AQ T PCA . Note the multiple optimal solutions to P 2 provide more options for CI system implementation. A practical issue in CI is the high feature measurement dynamic range. To keep a good reconstruction performance, CI system measurement SNR needs to be large. Then, for projections such as PCA and Hadamard, the first feature collected is the summation of all object pixel values. This feature is significantly larger than others. Thus detectors having large dynamic range are required. However, as what we find in section 2.1 and 2.2, any matrix F = AQ T PCA presents the optimized measurement SNR or reconstruction MSE. Therefore, we can choose one among all of these matrices, which has the most relaxed dynamic range requirement to system detectors.

Optimized binary sensing matrices
In Sec. 2.2, we conclude that, F opt = AQ T PCA presents the minimal MSE for any orthonomal matrix A. Note the row vectors of F opt are M orthonomal vectors in a space Q, which is spanned by { q 1 , q 2 , · · · , q M }. A matrix F presents a small MSE, if its distance from the space Q, defined as F − AQ T PCA F , is small. Hence we can consider searching the optimal binary projection to minimize MSE as to search a binary matrix sgn{AQ T PCA } close to Q. sgn() is the signum function, which has value {1, 0, −1} for positive, zero, or negative input, respectively. Note that, if all elements in sgn{AQ T PCA } are assumed nonzero, then the L 2 norm of each sensing vector is N. To make both of the binary and the non-binary vectors having same L 2 norms, we multiply √ N to AQ T PCA . Then designing binary matrices becomes solving the following problem, where · F represents the Frobenius norm of a matrix. We redefine the constraint as Using these definitions, the Lagrangian for P 3 can be written as where [ · ] i j represents the (i, j) th element of a matrix, λ is a symmetric matrix of size (M × M), and . * represents the element multiplication between two matrices. P 3 is an equality-constrained optimization problem. To solve P 3 , we use the Newton-Raphson method [41], which is specifically designed for nonlinear equality-constrained optimization problem. The algorithm can be represented using the following equations, and In Eq. (13) and (14), ∇ 2 AA L is the Hessian matrix of the Lagrangian, J is the Jacobi of the constraint for P 3 , and ∇ε is the gradient vector of the cost function. The di- , and (M 2 × 1), respectively. Different from Eq. (12), the variables in Eq. (13) and (14) are vectors a and λ , which are obtained by lexicographically ordering the elements of A and the upper triangle part of λ . The vector a is defined as a In this section, the definitions of ∇ 2 AA L , J, ∇ε are presented directly. More details about the calculation of these parameters can be found in the Appendix section.
The elements of the Hessian matrix ∇ 2 AA L ( a (ν) , λ (ν) ) are defined as Function δ (i) is the delta function. For the Jacobi matrix of G, note that G is symmetric. Therefore only the upper triangle part of G is used to find J. We have The third parameter in the Newton-Raphson method is ∇ε( a) with To ensure the step direction Δ a (ν) is a descent direction for ∇ε( a (ν) ) + J( a (ν) ) T λ (ν) , the Hessian matrix ∇ 2 AA L is modified to become positive definite [41,42]. In addition, the stepsize α in Eq. (14) is set to be small, and the number of the iterations ν is set to be large, so that ∇ε will converge to a cost value without oscillation. As discussed in section 1, collecting features or linear combinations of object pixels will benefit system measurement SNR. Therefore, CI is well suitable for low-light-level imaging (L 3 -imaging) [6,7]. Figure 1 is a block-based system diagram for L 3 -imaging. In such a system, objects are defined into blocks. Then, for each block, multiple compressive measurements are collected using a detector. The advantage of using block-based CI (BCI) instead of CI with a single detector is to obtain high resolution object reconstructions [33]. Figure 2(a) and 2(b) present the two kinds of objects used for the experimental study [43,44]. They are low-light-level images and are rather dark, so for illustrations only we enhanced To make this SNR range having more physical meaning, we use Hamamatsu SI PIN photodiode S1223 as an example to calculate the input light power. The detector S1223 has NEP (noise equivalent power) 9.4 × 10 −15 W/Hz 1/2 . If the system bandwidth is assumed as 30MHz, then when SNR is 0 dB, the input light power or the noise power is 51.486pW. In our experimental study, the lowest SNR is −13.02 dB. At this SNR level, the input light power is 2.568pW. In other words, for such level of light power, directly using PIN photodiode S1223 for measurement collection, the input signal will be overwhelmed by noise. However, using a BCI system as shown in Fig. 1, we can still make the measurements, and get an object reconstructed. In the first experiment, we studied seven kinds of projections, such as three Hadamard projections, two DCT projections, PCA projection Q T PCA , and the projection AQ T PCA . In the first Hadamard projection, H freq , we sort the sensing vectors based on the variation of elements' signs. For example, the first vector has all one values. The element's sign varies 0 time. In the second vector, the first half elements are +1, and then the rest are −1. The sign varies once. In the third vector, the element values change from +1 to −1, and then back to +1. The sign varies twice. Note that, using this sorting method, the system collects the DC value, the low frequency components, and the high frequency components of the object x sequentially. In the second Hadamard projection, H mean , we use a training sample average to sort the sensing vectors. The sensing vectors are in the same order as the Hadamard coefficients of the sample average. In the third Hadamard projection, H snr , we sort the vectors using measurement SNR [7]. The two DCT for this experiment are the regular DCT, F DCT , and the DCT with sensing vectors sorted using measurement SNR, F DCT snr . Table 1 presents the cost function ε(F) = Tr{FR x F T } for problem P 1 using these projections. It can be observed, the features obtained using PCA and AQ T PCA have the largest ε(F) as we discussed in Section 2.1. In Fig. 3(a) and 3(b), we present the normalized reconstruction root mean square error (RMSE) vs. noise σ for the seven kinds of projections. The RMSE is defined as x est − x / x . As expected, for both kinds of objects, Q T PCA and AQ T PCA present the same reconstruction error performance at different noise levels. These errors are smaller than the errors obtained using other kinds of projections. This is consistent with the conclusion obtained in section 2.2, that the projection AQ T PCA for any orthonomal matrix A is optimal in terms of minimizing RMSE. Note that, comparing Fig. 3(a) and 3(b), we find the RMSE values for the large size object are smaller. The reason is that, for a large size object, pixel values in object blocks vary slowly, or blocks are sparse in frequency domain. Therefore, using a few feature measurements can reconstruct the block with small error.

Numerical results for fast compressive
In the second experiment, we use the iterative algorithm discussed in section 3 to search for the optimal binary projection for the object shown in Fig. 2(b). Note that, the solution to this algorithm is a local minimizer. To find a global solution, 10M 2 initializations with 4000 iterations in each trail are used, where M is the number of features, and M 2 is the number of unknown variables in problem P 3 . For example, if M = 8, the number of initializations used to search the optimal matrix A of size (8 × 8) is 640. In addition, to maximize the object signal in a feature measurement, we normalize the maximum absolute value of each sensing vector to 1.
In Fig. 4(a), the 8 Hadamard vectors in H snr , which are selected based on measurement SNR [7], are presented. Figure 4(b) presents the optimization result, F rmse , for M = 8. It can be observed, except the all-one vector, the solution for P 3 is not directly related to Hadamard projection.
Besides the optimal solution, we also compare the 640 local minimizers obtained by the 10M 2 trials. The RMSE values obtained using these local solutions are calculated. We find 217 Fig. 4. (a) The Hadamard vectors sorted using measurement SNR, (b) the optimal binary vector set or the optimal solution for problem P 3 , (c) the suboptimal solution for P 3 which has the minimal feature measurement dynamic range, when M is 8 for 32 × 32 blocks. local minimizers including the optimal one presents very similar RMSE. The RMSE difference between these local minimizers and the optimal solution is less than 10 −4 . Then among these sub-optimizers, the solution having the smallest feature dynamic range, labeled as F dyn , is presented in Fig. 4(c). The all-one vector is not among the result. Table 2 summarizes the minima, the maxima, and the dynamic ranges of features obtained using H mean , H snr , F DCT snr , PCA, F rmse , and F dyn . It is clear that, F dyn needs about 30% less dynamic range than the other five projections. Figure 5(a) presents the reconstruction errors using H mean , H snr , F DCT snr , Q T PCA , F rmse , and F dyn vs. noise σ . It can be observed that F rmse presents even smaller RMSE than PCA, although it is binary. The reason is, after normalizing each vector having maximum 1, a system using F rmse can collect more object light into a measurement than using Q T PCA . Thus system measurement SNR and reconstruction RMSE are improved. From the figure we can also observe that F dyn presents similar RMSE values as F rmse . As a comparison, we also present the reconstruction RMSE for the standard Hadamard and random binary projections in Fig. 5(b). It is clear that, both projections have larger reconstruction errors. To investigate the reason that the optimized binary projections work well, we calculate the Fourier transform of each projection vector. Figure 6 shows the normalized FFT of the vectors in the standard Hadamard, H snr , and F dyn matrices. The FFT for F rmse is similar to F dyn . From the figure, it is clear that, F dyn will collect object information at more frequencies. H snr will collect at fewer frequencies, while the standard Hadamard collects at the fewest. Thus the reconstruction using F dyn has the smallest RMSE.
The experiment is repeated for M = 16. Sensing matrices F rmse and F dyn are shown in Fig. 7  as functions of noise in Fig. 7(c). Once again, the optimal and the suboptimal sensing matrices present very similar reconstruction performance. We also calculate the feature dynamic range values. We find the feature dynamic range is reduced 28% from 7.47 × 10 4 for H snr to 5.39 × 10 4 for F dyn . Figure 7(d) and 7(e) presents the normalized FFT of the vectors in H snr and F dyn matrices. Once again, we observe that F dyn will collect more object information at more frequencies.
Thus the system using this matrix has better reconstruction performance. To make the reconstruction more visually pleasant, we increase the intensity 3 times. It is clear, all reconstructions using CI have smaller RMSE compared with the conventional result. The RMSE values for the conventional imaging, H mean , H snr , and F dyn are 0.7312 , 0.0733, 0.0645 and 0.0635, respectively. Note that, the reconstructions using H snr , and F dyn are very similar to each other, and are close to the original.
To have a better observation, Fig. 9 shows two magnified detail parts of both reconstructions.   Figure 9(c) and 9(d) are for F dyn . It can be observed that, the reconstruction using F dyn has slightly smoother reconstruction with less block effect compared with the reconstruction using H snr .
In the next experiment, we add a bias to each projection to make it nonnegative. Then the maximum of the obtained projection is normalize to 1. For the PCA and DCT projections, the vector values are quantized to 128 levels. As discussed in the beginning of the paper, CI systems using gray value projections, require longer time to make measurements. In the time to collect a measurement with a 8-bit projection, a system can make 128 measurements using binary vectors. Therefore, binary projections benefit a CI system imaging speed. Figure 10 presents the reconstruction RMSE with these projections. We get the same conclusions as in Fig. 7(c).
The RMSE values for F rmse and F dyn are almost same. They are smaller than the RMSE values using other projections. Among the rest four projections, H snr has the smallest reconstruction at almost all noise levels.  In the last set of experiments, we studied the performance of F rmse and F dyn with sparse objects. The object blocks in Fig. 2(b) are transformed into the wavelet domain using the 4-level wavelet transformation matrix. Except the largest 64 coefficients in each block, the rest are set to zero. Using this sparse object, we repeat the experiment as in Fig. 10. Besides the six matrices we discussed in last several experiments, wavelet and curvelet projections are also studied. The wavelet projection matrix is for 4 level transformation. Both kinds of projection vectors are sorted using measurement SNR. These vectors are also added with a bias to make them nonnegative, and then quantized to 128 discrete levels. Figure 11 presents the reconstruction RMSE vs. noise σ . We observe similar trend as in Fig. 10, except that the RMSE values are smaller because the object sparsity benefits the reconstruction. For wavelet and curvelet transforms, the reconstruction error are larger than other projections at most noise levels. The reason is that, both transformation vectors are spatial localized. They are suitable to extract object's local information and make the transformation sparse. However, this property is not helpful in compressive low-light-level imaging, in which more light from an object are required to surpass detector noise. Between wavelet and curvelet transforms, the performance of the latter is worse. This is caused by the non-orthogonality among the transformation vectors.
In the last experiment, we make the object further sparse with K = 4 for each block. The sparse object is shown in Fig. 12(a). Then we use a nonlinear method, fast linearized Bregman iterative algorithm, to reconstruct the object. In Fig. 12 to make the feature measurements. In Fig. 12(c), F dyn is used. The noise level in the experiment is σ = 148.67. We can observe that the reconstruction using F dyn is visually better. The RMSE for Fig. 12(b) is 0.0623, while the RMSE for 12(c) is 0.0268. The designed binary sensing matrix works well for the sparse object. Because it is more structured in Fourier domain, it works better than random projections. Additionally, we also present the reconstruction using Wiener operator in Fig. 12(d) when F dyn is used. The performance is even better than using the linearized Bregman iterative algorithm. As we discussed in the beginning of the paper, although Wiener operator does not use sparsity for reconstruction, it incorporates object second order statistical parameter, object autocorrelation matrix R x . It is more suitable when system measurements have much noise. The RMSE value obtained using Wiener operator is 0.0131 for the sparse object. Table 4 summarizes the measurement collection time and the object reconstruction time using different kinds of sensing matrices and linear/nonlinear reconstruction methods. We assume TI Discovery 4100 development tool set is used to make measurements. Then the DMD working frequency is set as 32kHz. With this working frequency, for 16 binary features, we need 0.5ms to make the measurements. For 128-level gray sensing matrix, the time becomes 64ms. Using the linear reconstruction method, Wiener operator, 12ms is required to obtain the reconstruction. Using the nonlinear method, fast linearized Bregman iterative algorithm, the reconstruction time is more than 100s. Thus, to make the system imaging speed faster, binary sensing matrix and Wiener operator should be used.

Conclusions
In this paper, we study to design binary sensing matrices for fast compressive L 3 -imaging. Two optimization problems, maximizing measurement SNR and minimizing RMSE, are solved first. We find the solutions to both problems are matrices AQ T PCA , where A is orthonomal. Based on the observation, a binary matrix design problem is defined as to search a matrix sgn(AQ T PCA ). Then from the experimental study, we find the optimal solution presents similar reconstruction as PCA, or even better than PCA when the maximum absolute values of projection vectors are normalized to 1. More interestingly, many suboptimal binary matrices are found to have almost the same reconstructions as the optimal. Among them, we can choose the matrix presenting the least measurement dynamic range to relax the requirement to CI sensors.
We focus on Gaussian noise in this work. When the light from object is extremely low, the noise becomes Poisson. For this case, the compressive imaging idea will still help the system's reconstruction performance, because system measurement SNR will be improved using BCI. However, the linear Wiener operator and the nonlinear reconstruction algorithms used for CI need to be modified for Poisson noise. Besides Gaussian noise, in a physical system there are other kinds of detector noise, long-term drift error, and systematic errors from stray-light. Different kinds of detectors or imaging elements have different noise. For photomultipliers, photoconductive and photovoltaic detectors, shot noise, generation-recombination noise or g-r noise, thermal noise, and 1/ f noise are the four main factors. Shot noise matches Poisson noise probability density distribution. When the signal frequency is not very high, g-r noise can be modeled as white noise. Thermal noise generally is white Gaussian noise. 1/ f noise can be ignored in most cases when a signal frequency is higher than 1kHz. For different kinds of noise, using BCI, we can always have larger detector signal. When the signal power is larger, Poisson noise becomes Gaussian noise. For long-term drift, except calibrating the imaging system, using the optimized binary sensing matrix will shorten the measurement collection time. Thus it will help the system to deal with such kind of error. Stray-light will cause detector collecting undesired background or scattered light. For active imaging, modulating the illumination source and then preprocessing the detected signal using a bandpass filter will suppress stray-light in a system.
We study fast compressive L 3 -imaging in this paper. Thus, we use binary sensing matrices to make measurements and linear Wiener operator for object reconstruction. However, for applications in which imaging speed is not restricted and object light is not low or noise is not high, gray-value sensing matrices and nonlinear reconstruction methods have their advantages. To further improve the reconstruction performance of fast compressive L 3 -imaging, in the future we will investigate linear methods including object sparsity as an prior.
To find ∇ 2 AA L ( a (ν) , λ (ν) ), we first calculate ∇ε( a). Note that the cost function of P 3 , ε( a), can be written as ε( a) = Taking the derivative of ε( a) respects to A i j , we have We assume AQ T PCA nonzero. Therefore, ∂ ε/∂ A i j becomes where Q T PCA Q PCA = I. To define the Hessian matrix, we need to find ∂ 2 L /∂ A i j ∂ A i j . If we represent the second term of Eq. (12) as ε λ , then ∂ 2 L /∂ A i j ∂ A i j = ∂ 2 ε/∂ A i j ∂ A i j + ∂ 2 ε λ /∂ A i j ∂ A i j . Using the