Single-Shot Compressed Imaging via Random Phase Modulation

: Compressed sensing (CS) provides an innovative framework for signal sampling, which enables accurate recovery of the sparse or compressible signal from a small set of linear measurements far fewer than the Nyquist rate in traditional signal processing. In compressed sensing, random modulation plays a key role, which can spread out the signal information more or less evenly across all locations. There are many modulation techniques, such as amplitude modulation, frequency modulation, phase modulation, spectrum modulation, and so on. Among these modulation techniques, phase modulation is vital due to the efﬁciency and convenience of modulation. In this paper, we review both the theoretical and application of compressed sensing and several compressed imaging systems using random phase modulation. First, we review the fundamentals of compressed sensing, dividing it into three parts: sparse representation, incoherent measurement, and nonlinear reconstruction algorithm. We then show how phase modulation can be applied to compressed sensing and compressed imaging, where the presentation can be divided into six main parts, corresponding to different aspects of phase modulation applied in compressed sensing or compressed imaging: (1) Fundamentals of compressed sensing. (2) Principles of phase modulation. (3) Single-shot compressed imaging with spatial-domain single random phase mask (CI-SSRPM). (4) Single-shot compressed imaging with a random convolution using a double random phase mask (CI-DRPM). (5) Single-shot compressed imaging with Fourier-domain single random phase mask (CI-FSRPM). (6) Single-shot compressed imaging with double random phase encoding (CI-DRPE).


Introduction
The Nyquist/Shannon theorem has been considered the doctrine of signal acquisition and processing since 1948, which asserts that an arbitrary band-limited signal can be accurately reconstructed from its samples with the sampling rate to be at least twice the bandwidth. One of the major drawbacks is that the Nyquist rate is so high that most of the applications require signal compression before storage or transmission, which takes too much time for signal sampling.
In 2006, compressed sensing (CS) was firstly proposed by Candès [1] and Donoho [2], considered a breakthrough in the signal processing community, aiming at reducing this high sampling rate. Compressed sensing (CS) provides an innovative framework for signal sampling, which enables accurate recovery of the sparse or compressible signal from a small set of linear measurements far fewer than the Nyquist rate in traditional signal processing.
In CS, randomness is exactly crucial because it can spread out the information of the signal more or less evenly across all locations, which is one key for CS to successfully recover from fewer measurements less than the Nyquist rate. In classical CS, the random Gaussian matrix (RGM) is the most popular measurement matrix due to its theoretical guarantees. However, the drawback of RGM is that it is impractical or too expensive for physical implementation, and other random modulation techniques must be introduced to realize random modulation.
CS theory is not only available to traditional imaging applications but also available in more complex cases, e.g., infrared imaging and terahertz imaging. The introduction of CS theory can help enhance the performance of various imaging systems in terms of high-quality reconstruction, low cost, small size and etc. From the aspect of physical realization, various modulation techniques existed, such as amplitude modulation, frequency modulation, phase modulation, spectrum modulation, and so on. Among these modulation techniques, phase modulation plays a vital role in the efficiency and convenience of phase modulation.
In 2009, Romberg [3] proposed a compressed imaging system by random convolution and that the convolution with a specific-designed random pulse followed by random subsampling was considered as a universal compressed sampling operator, which was equivalent to a random phase modulation in the Fourier-domain of the well-known 4-f system. In 2010, Rivenson et al. [4] proposed a single-exposure compressed imaging system by combining the well-known optical scheme of double random phase encoding (DRPE), which could be applied to overcome the resolution limit caused by both diffraction and geometrical limitation. In 2014, Wu et al. proposed a quasi-monochromatic thermo optical phase modulation compressed sensing imaging system, which overcomes the problem of collecting too much data due to high-resolution requirements in the imaging process [5]. In 2017, Liu et al. [6] proposed a single exposure compressed imaging system based on the conventional optical 4-f architecture with an RPM placed in the Fourier plane and a random coded aperture in the image plane for subsampling. In 2021, Li et al. proposed the correlation imaging of moving targets based on phase modulation. The spatial random phase modulator was introduced into the imaging optical path to modulate the threedimensional spectrum data of the target, and then the detector was used to obtain the two-dimensional aliasing signal to realize the reconstruction by single exposure [7].
In this paper, both the theoretical and applications of compressed sensing and several compressed imaging systems using random phase modulation are reviewed. First, the fundamentals of compressed sensing are presented. Second, the principle of phase modulation is described. Finally, several strategies of compressed sensing and compressed imaging using random phase modulation are reviewed. The rest of the paper is organized as follows: In Sections 2 and 3, the fundamentals of compressed sensing and the principles of phase modulation are both described. In Section 4, different schemes of compressed imaging using random phase modulation are presented, which include: (i) singleshot compressed imaging with spatial-domain single random phase mask (CI-SSRPM).
(ii) Single-shot compressed imaging with a random convolution using a double random phase mask (CI-DRPM). (iii) Single-shot compressed imaging with Fourier-domain single random phase mask (CI-FSRPM). (iv) Single-shot compressed imaging with double random phase encoding (CI-DRPE). Section 5 is the conclusion.

Fundamentals of Compressed Sensing
In 2006, Compressed Sensing (CS) was firstly proposed by Donoho and Candès et al. Compared to the traditional Shannon-Nyquist sampling theorem, CS provides a new sampling framework for band-limited signals. As we all know, most signals in the real world are sparse or compressible in some bases or dictionary, which means that they may have a concise representation by selecting the appropriate transform domain or dictionary.
The success of CS relies on two properties: sparsity and incoherence. Sparsity is the inherent property of those signals. The information contained in the signal can be represented only with the help of a few significant components, as compared to the total length of the signal. Incoherence provides a guideline for the design of the measurement matrix, and performance guarantees for successful recovery of the signal from fewer measurements.
CS involves three major foundations: sparse representation, random measurement matrix, and nonlinear recovery algorithm.

Sparse Representation
In compressed sensing, the sparsity of the signal means that the signal can be represented on a suitable basis with fewer significant coefficients than the signal length. Examples of the bases include discrete Fourier transform (DFT), discrete cosine transform (DCT), discrete wavelet transform (DWT), etc.
In mathematical, the sparse representation of the signal can be expressed in the following form where Ψ is a sparsifying matrix, the size is N × N, ψ i is the column of matrix Ψ, Ψ i is the i-th column, i = 1, 2, · · · , N, α is the sparse coefficients whose has only K non-zeros (i.e., sparsity) or significant entries with K << N.

Measurement Matrix
Different from the Shannon sampling scheme, in compressed sensing, the subsampling measurements can be obtained via linear and incoherent projections corrupted by noise as follows where Φ ∈ R M×N is measurement matrix, Θ is the sensing matrix and e is the random additive noise.
Since the measurement process of CS is not adaptive, the measurement matrix is fixed and independent of the signal x, which means that the measurement matrix should be universal for any sparse or compressible signal. The effects of the measurement matrix are reflected in two aspects: On the one hand, a stable measurement matrix Φ can preserve salient information in any K-sparse or compressible signal in a compressed sampling process; On the other hand, the measurement matrix plays a crucial role in the reconstruction process, a 'good' measurement matrix can clearly enhance the recovery performance for recovering x only from the M measurements y.
In order to ensure the successful recovery of a sparse signal from fewer measurements, CS theory indicates that if the Restricted Isometry Property (RIP [8]) is satisfied, i.e., the measurement matrix satisfies the K-order RIP of the parameter δ K ∈ (0, 1), if holds for all sparse vectors α having no more than K nonzero entries.

Nonlinear Recovery Algorithm
Once the measurements vector y and measurement matrix Φ are given, the way to find a signal x in the class of interest so that y = Φx is complete or approximate is a problem that the reconstruction algorithm should solve, termed the sparse recovery algorithms in CS.
In recent years, there have been many sparse recovery algorithms proposed to recover a sparse signal from fewer incoherent measurements. These algorithms can be categorized into three main types: convex relaxation algorithms, non-convex optimization algorithms, and greedy algorithms.

L0-Norm Minimization Recovery
Intuitively, we can solve the problem in Equation (2) by 0 norm minimization using the following optimization.x = argmin Unfortunately, solving Equation (4) relies on an exhaustive search, which has combinatorial computational complexity. Therefore, the goal is to find a computationally feasible algorithm that can successfully recover the sparse signal x from the measurement y.

L1-Norm Minimization Recovery
An alternative is to use 1 norm and defined the 1 norm of a vector x 1 as x 1 = ∑ N n=1 x(n) . Then the optimization problem solved can be rewritten aŝ Due to the convexity of 1 the norm, Equation (5) can be regarded as the convex relaxation version of Equation (4), which can be solved via linear programming, with the computational complexity polynomial over the length of the signal.

Principles of Phase Modulation
In optical sensing and imaging system, light modulation has the potential to exceed the diffract limit and achieve higher contrast in imaging; among various light modulation techniques, phase modulation has gained more attention than its amplitude modulation due to its notable sensitivity and close relation to spectrum, coherence and polarization of light, and so on. Moreover, phase modulation can be easily implemented using various optical devices, spatial light modulator (SLM), liquid crystal display (LCD) and so on, which makes its use more efficient and convenient.
The phase component carries much more information than the amplitude component, such as the shape, depth, polarization, and so on. Therefore, phase plays a role in many applications, such as system encryption [9] or electro-optic conversion of microwave signals [10]. For example, as illustrated in Figure 1, exchange the Fourier phase of the two images and keep the amplitude unchanged, then the major contexts of both image exchange, which means that the Fourier phase contains the most amount of information about the image.
Phase modulation (PM) is the modulation technique that utilizes various optical devices to change the optical path difference (OPD) of an incident light to produce different phase delays [11]. Phase modulation can be achieved either by changing the refractive index or the thickness of an optical device, as indicated in the following equation where u in is the incident light, and u out is the exit light, ϕ is the modulated phase of the optical device expressed by where d is the thickness of the devices or equivalently the thickness caused by the change of the refractive index, n 0 denotes the refractive index of the passing-through material and λ is the wavelength of the incident light. The phase modulation can vary the direction or the shape of the wavefront of incident light.
Appl. Sci. 2022, 12, x FOR PEER REVIEW 5 of 16 many applications, such as system encryption [9] or electro-optic conversion of microwave signals [10]. For example, as illustrated in Figure 1, exchange the Fourier phase of the two images and keep the amplitude unchanged, then the major contexts of both image exchange, which means that the Fourier phase contains the most amount of information about the image. Phase modulation (PM) is the modulation technique that utilizes various optical devices to change the optical path difference (OPD) of an incident light to produce different phase delays [11]. Phase modulation can be achieved either by changing the refractive index or the thickness of an optical device, as indicated in the following equation where in u is the incident light, and out u is the exit light, j is the modulated phase of the optical device expressed by where d is the thickness of the devices or equivalently the thickness caused by the change of the refractive index, 0 n denotes the refractive index of the passing-through material and l is the wavelength of the incident light. The phase modulation can vary the direction or the shape of the wavefront of incident light.

Compressed Imaging Using Random Phase Modulation
The combination of random phase modulation and compressed sensing has been proposed and applied on many occasions, such as in various types of single-exposure imaging. Using them in imaging can not only overcome the problem of low resolution but also restore the possibility for the reconstruction of complex data or moving targets.
In this section, we briefly review four compressed imaging strategies using random phase modulation: (i) compressed imaging with spatial-domain single random phase mask (CI-SSRPM). (ii) Compressed imaging with double random phase mask (CI-DRPM).

Single-Shot Compressed Imaging with Spatial-Domain Single Random Phase Mask (CI-SSRPM)
In this section, single-shot compressed imaging with spatial-domain single random phase mask (CI-SSRPM [12]) is presented. The optical setup is depicted in Figure 2. As illustrated in Figure 2, the random phase mask is attached to the middle between the two

Compressed Imaging Using Random Phase Modulation
The combination of random phase modulation and compressed sensing has been proposed and applied on many occasions, such as in various types of single-exposure imaging. Using them in imaging can not only overcome the problem of low resolution but also restore the possibility for the reconstruction of complex data or moving targets.
In this section, we briefly review four compressed imaging strategies using random phase modulation: (i) compressed imaging with spatial-domain single random phase mask (CI-SSRPM). (ii) Compressed imaging with double random phase mask (CI-DRPM). (iii) Compressed imaging with Fourier-domain single random phase mask (CI-FSRPM). (iv) Compressed imaging with double random phase encoding (CI-DRPE).

Single-Shot Compressed Imaging with Spatial-Domain Single Random Phase Mask (CI-SSRPM)
In this section, single-shot compressed imaging with spatial-domain single random phase mask (CI-SSRPM [12]) is presented. The optical setup is depicted in Figure 2. As illustrated in Figure 2, the random phase mask is attached to the middle between the two lenses L 1 and L 2 . The object plane is at the position z 1 in front of the lens L 1 , and the right side of the lens is a random phase mask with the phase distribution ϕ(r ϕ ). The right side of the phase mask is closely attached to the lens L 2 , and the image plane is z 2 from the right side of the lens.
The optical field emitted by the object is focused by the lens L 1 and then randomly modulated by a random phase mask; the scattered field is then focused by the lens L 2 . Finally, the complex field of the object reaches the CCD detection array behind the lens. The entire measurement process is achieved in a single exposure.
The advantage of this system is that all the required measurements can be obtained in a single exposure. The forward projection operator (i.e., measurement matrix) between the object point r o and the image point r i can be expressed as where λ denotes the wavelength of the incident light, ϕ(r ϕ ) denotes the random phase mask, f L 1 and f L 2 denotes the focal length of Lens L 1 and L 2 , and C denotes a multiplicative constant. Equation (8) defines the relationship between the object point and image point as integrated overall random phase mask pointing ϕ(r ϕ ). Such relation is spatial random and shift-variant due to the randomness of the phase mask. The relation in Equation (9) can be discretized into the following form ∆ o and ∆ i are the pixel size of the object and the CCD.
As the element of the Gaussian matrix has too much freedom, which cannot be physically achieved or cost too much, some deterministic matrices [13] with fewer independent elements have important practical significance and are easy to implement to satisfy the specific distribution of optical phase mask (such as Toeplitz or Circulant distribution), which can be expressed as follows Equations (10)-(12): The numerical experiments of CI-SSRPM are depicted in Figure 3. Lena image is chosen as the test image with the size 64 × 64. Daubechies 5 wavelet transform is chosen as the sparsifying basis. The subspace pursuit (SP [14]) algorithm is used as the recovery algorithm. The reconstruction quality is quantified using the signal-to-noise ratio (SNR) and shown in Table 1. L and 2 L . The object plane is at the position 1 z in front of the lens 1 L , and the right side of the lens is a random phase mask with the phase distribution ( ) r j j . The right side of the phase mask is closely attached to the lens 2 L , and the image plane is 2 z from the right side of the lens. The optical field emitted by the object is focused by the lens 1 L and then randomly modulated by a random phase mask; the scattered field is then focused by the lens 2 L .
Finally, the complex field of the object reaches the CCD detection array behind the lens. The entire measurement process is achieved in a single exposure. The advantage of this system is that all the required measurements can be obtained in a single exposure. The forward projection operator (i.e., measurement matrix) between the object point o r and the image point i r can be expressed as where l denotes the wavelength of the incident light, ( )

Single-Shot Compressed Imaging with Random Convolution Using Double Random Phase Mask (CI-DRPM)
In 2009, Romberg [4] proposed a compressed imaging system by random convolution, that the convolution with a specific-designed random pulse followed by random sub-sampling was considered as a universal compressed sampling operator, which was equivalent to implementing a random phase modulation in the Fourier-domain of the well-known 4-f system. In Ref. [4], Romberg proposed a novel strategy for compressed sensing which utilizes random convolution as a universal and efficient random measurement operator. Such measurement operator is constructed by firstly convolved with a specific-designed pulse and successively by random sub-sampling. The random convolution for compressed sensing can be algebraically written as where ⊗ denotes convolution operator, F denotes the discrete Fourier matrix with entries (F ij ) = [exp −j2π(i−1)(j−1)/N ] 0≤i,j≤N−1 and F H denotes the inverse, and Σ is the Fourier transform version of pulse h in the following form where Σ is a diagonal matrix with the entries σ i (i = 1, 2, · · · , N) has a unit amplitude and random phases, which can be generated as follows: The random demodulation (RD) method is often used for random subsampling. The image x Σ can be divided into blocks, and successively a measurement vector can be obtained by multiplying the entries in each block by a sequence of random signs and summing. Let Ξ be a diagonal matrix with the entries independently drawn from ±1 Bernoulli distribution with equal probability, and P ∈ R M×N denotes an accumulate-and-dump subsampling operation that can be expressed as follows: For simplicity, a renormalization factor is introduced into the random demodulation process, and finally, the random measurement vector can be algebraically rewritten as where Φ denotes the forward measurement operator (i.e., measurement matrix) in compressed sensing via random convolution followed by random demodulation.
The subsampling strategy of random convolution compressed sensing (RCCS) has several attractive advantages, such as universality, numerical structure, physically realizable, which can provide an excellent strategy to perform universal CS acquisition, fast reconstruction, and a simply equipped setup to record the compressed measurements in a single exposure. The scheme of RCCS is depicted in Figure 4. As illustrated in Figure 4, an SLM coded with the pulse h is placed in the Fourier domain, and another SLM with random sign entries in the image plane. Numerical results will be presented in the next subsection for comparison. Figure 5 is a schematic diagram of a 2D image in a 4-f optical setup with a random phase modulated in the frequency domain. The first lens transforms the standard Lena image and is then followed by a phase modulation that may be performed by an SLM after another lens transforms back into the time domain. Finally, a low-resolution detector array records the measurements. In Refs. [15,16], we proposed a novel compressed imaging using single random phase modulation in the Fourier domain. The optical setup of our proposed is depicted in Figure  6. As depicted in Figure 6, the system is based upon the architecture based on a 4-f optical setup. A phase mask that may be manipulated by an SLM placed in the frequency domain plays a crucial role in maintaining the amplitude of the Fourier transform of the original image unchanged; only the phase of the Fourier transforms of the image makes the type uniformly random delay point-wise. Such modulation can reduce the sampling burden in a spatial domain because the entire information is spread out in the spatial domain. As illustrated in Figure 4, an SLM coded with the pulse h is placed in the Fourier domain, and another SLM with random sign entries in the image plane. Numerical results will be presented in the next subsection for comparison.  As illustrated in Figure 4, an SLM coded with the pulse h is placed in the Fourier domain, and another SLM with random sign entries in the image plane. Numerical results will be presented in the next subsection for comparison. Figure 5 is a schematic diagram of a 2D image in a 4-f optical setup with a random phase modulated in the frequency domain. The first lens transforms the standard Lena image and is then followed by a phase modulation that may be performed by an SLM after another lens transforms back into the time domain. Finally, a low-resolution detector array records the measurements. In Refs. [15,16], we proposed a novel compressed imaging using single random phase modulation in the Fourier domain. The optical setup of our proposed is depicted in Figure  6. As depicted in Figure 6, the system is based upon the architecture based on a 4-f optical setup. A phase mask that may be manipulated by an SLM placed in the frequency domain plays a crucial role in maintaining the amplitude of the Fourier transform of the original image unchanged; only the phase of the Fourier transforms of the image makes the type uniformly random delay point-wise. Such modulation can reduce the sampling burden in a spatial domain because the entire information is spread out in the spatial domain. In Refs. [15,16], we proposed a novel compressed imaging using single random phase modulation in the Fourier domain. The optical setup of our proposed is depicted in Figure 6. As depicted in Figure 6, the system is based upon the architecture based on a 4-f optical setup. A phase mask that may be manipulated by an SLM placed in the frequency domain plays a crucial role in maintaining the amplitude of the Fourier transform of the original image unchanged; only the phase of the Fourier transforms of the image makes the type uniformly random delay point-wise. Such modulation can reduce the sampling burden in a spatial domain because the entire information is spread out in the spatial domain. Appl. Sci. 2022, 12, x FOR PEER REVIEW 11 of 16 For a 2D image, the forward measurement operator of our proposed CI-FSRPM scheme can be expressed as

Single-Shot Compressed Imaging with Fourier-Domain Single Random Phase Mask (CI-FSRPM)
And the adjoint operator can be constructed as For a 2D image, the forward measurement operator of our proposed CI-FSRPM scheme can be expressed as where denotes Hadamard product. F 2D and F −1 2D denote 2D Fourier transform and 2D inverse Fourier transform operator and matrix Σ is a N × N matrix whose entries are the (renormalized) 2D Fourier transform coefficients of h, i.e., Σ = F 2D (h). Nevertheless, we first generate Σ at random by taking where the entries σ ij , i = 1, · · · , N, j = 1, · · · , N are unit magnitude complex numbers with random phases which means pure phase mask. Note that the entries of Σ must satisfy conjugate symmetric property for simplicity so that real object x can still be real after performing the operator Φ. Specifically, the entries Σ are generated as follows: For N even, For N odd, And the adjoint operator can be constructed as where Σ H is the conjugate of Σ, D ↑ is the integral upsampling. Note that one advantage of the construction is that the modulation scheme M is orthogonal, i.e., where I is the matrix with all entries equal to one. In this section, numerical experiments are presented to verify the effectiveness of CI-DRPM and CI-FSRPM. The test environment is set with a 32-bit Windows7 operating system, the Intel (R) Core (TM) i5-2320 and 4-core processor, a basic frequency 3.00 GHz, an effective memory of 3.0 GB, and the test software Matlab2010a. Figure 7 shows the potential applications of random phase encoding in frequency domain simulation based on 4-f optical architecture. We create 1536 × 1536 original images using Lena and Phantom and take them using Figure 6 to achieve a random phase encoding and integral subsampling.
where H S is the conjugate of S , ↑   is the integral upsampling. Note that one advantage of the construction is that the modulation scheme  is orthogonal, i.e., where  is the matrix with all entries equal to one.
In this section, numerical experiments are presented to verify the effectiveness of CI-DRPM and CI-FSRPM. The test environment is set with a 32-bit Windows7 operating system, the Intel (R) Core (TM) i5-2320 and 4-core processor, a basic frequency 3.00 GHz, an effective memory of 3.0 GB, and the test software Matlab2010a. Figure 7 shows the potential applications of random phase encoding in frequency domain simulation based on 4-f optical architecture. We create 1536 × 1536 original images using Lena and Phantom and take them using Figure 6 to achieve a random phase encoding and integral subsampling.    [19]. The wavelet scale is the same as 3.
In another numerical experiment, we compare our proposed method with Random Demodulation (RD) method proposed by Justin Romberg [16] and the RecPC method proposed by W. Yin [20]. The test images are Lena, Cameraman, and Phantom image; the image size is chosen into seven levels of 256 × 256, 512 × 512, 768 × 768, 1024 × 1024, 1280 × 1280, 1536 × 1536, and 1792 × 1792, the total downsampling rate is 1/8 (the reason for this rate selection is due to the hybrid measurements of RD method, which mean two lowresolution measurements, one is the measurement with random demodulation, another is a direct low-resolution version of the original image, both sampling rates are 1/16, with   Figure 7c,f are compressed measurements of Lena and Phantom images, respectively. The measurement size is 384 × 384; that is to say, the sampling rate is four in the horizontal and vertical directions. Figure 7d,e are the reconstruction results of Lena with the measurements shown in Figure 7c using iterative soft thresholding (IST [17]) algorithm and gradient projection sparse reconstruction (GPSR [18]) algorithm. The SNR is 29.84 dB and 29.68 dB. Figure 7g,h are the reconstruction results from measurements shown in Figure 7f using the IST and GPSR algorithm. The SNR is 39.87 dB and 39.62 dB. The sparse operators Ψ have chosen the Daubechies 10 wavelet of Rice University (rice wavelet toolbox, RWT) [19]. The wavelet scale is the same as 3.
In another numerical experiment, we compare our proposed method with Random Demodulation (RD) method proposed by Justin Romberg [16] and the RecPC method proposed by W. Yin [20]. The test images are Lena, Cameraman, and Phantom image; the image size is chosen into seven levels of 256 × 256, 512 × 512, 768 × 768, 1024 × 1024, 1280 × 1280, 1536 × 1536, and 1792 × 1792, the total downsampling rate is 1/8 (the reason for this rate selection is due to the hybrid measurements of RD method, which mean two low-resolution measurements, one is the measurement with random demodulation, another is a direct low-resolution version of the original image, both sampling rates are 1/16, with the amount to 1/8. For a fair comparison, while RecPC method uses random location downsampling, that is, in all positions, the uniform random decimation measurement value of the position. In this experiment, the sample rate for horizontal and vertical are 1/2 and 1/4, respectively. The test environment is a 32-bit window7 operating system, the processor is Intel (R) Core (TM) i5-2320 quad-core, the frequency is 3.00 GHz, available memory is 3.0 G, and we used Matlab2010a software. Numerical results are listed in Table 2. As depicted in Table 2, the CI-DRPM method proposed by Romberg with a random demodulator (RD) for subsampling using tvqc_logbarrier solver can only deal with 1280 × 1280 pixels; larger is impossible. Our proposed CI-FSRPM is superior to CI-DRPM in terms of image size, reconstruction time, and reconstruction quality. The reason is that our proposed CI-FSRPM performs in 2D form, which is more effective than CI-DRPM performs in 1D form. Although the RecPC method has the highest SNR (especially sparse gradient images, such as the Phantom images), the reason is the subsampling type is by random sampling with a random position rather than random demodulation, which severely limits its practical application.

Single-Shot Compressed Imaging with Double Random Phase Encoding (CI-DRPE)
In 2010, Rivenson et al. [15] proposed a single-exposure compressed imaging system by combining the well-known optical scheme of double random phase encoding (DRPE), which could be applied to overcome the resolution limit caused by both diffraction and geometrical limitation. The schematic of compressed imaging using double random phase encoding (CI-DRPE) is depicted in Figure 8. As shown in Figure 8, the image is randomly modulated with the two random phase masks (RPMs) placed in the input plane and the Fourier plane of the 4-f optical system. random sampling with a random position rather than random demodulation, which severely limits its practical application.

Single-Shot Compressed Imaging with Double Random Phase Encoding (CI-DRPE)
In 2010, Rivenson et al. [15] proposed a single-exposure compressed imaging system by combining the well-known optical scheme of double random phase encoding (DRPE), which could be applied to overcome the resolution limit caused by both diffraction and geometrical limitation. The schematic of compressed imaging using double random phase encoding (CI-DRPE) is depicted in Figure 8. As shown in Figure 8, the image is randomly modulated with the two random phase masks (RPMs) placed in the input plane and the Fourier plane of the 4-f optical system. The modulation scheme of CI-DRPE can be expressed as follows   The modulation scheme of CI-DRPE can be expressed as follows where denotes Hadamard product. Σ RPM1 = exp(jP 1 ) and Σ RPM2 = exp(jP 2 ) denotes two random phase masks, both the size of P 1 and P 2 are the same as the image x, D ↓ denotes a random subsampling operator. The Equation (21) can be reformulated as a universal measurement operator in matrixvector form as follows where D ↓ is the random subsampling operation, F denote the discrete Fourier matrix (DFT) with size N × N (N = n 2 ) and F H denotes its inverse, Θ RPM1 = diag(Σ RPM1 ) and Θ RPM2 = diag(Σ RPM2 ), i.e., Θ RPM1 take the elements of random phase mask Σ RPM1 on the diagonal, so as the Θ RPM1 . x denotes the vector rearranged the input image lexicographically, and y denotes the compressed measurements. Note that the operation FΘ RPM1 in Equation (25) denotes random scrambling in the Fourier domain in the matrix form as: where W N = exp(j2π/N). The entries {p i } N i=1 are drawn independently from a uniform distribution [0, 1), which leads to inter-column statistical independence for FΘ RPM1 and equivalence to the random Gaussian measurement operator. Numerical experiments are shown in Figure 9. In Figure 9, the parameters are set as follows: Cameramen image is selected as the test image, the size is 512 × 512, the subsampling rate is 0.5, the measurements are complex, the type of subsampling is random location. The sparsifying basis is RWT with a wavelet scale is 3. The performance of CI-DRPE under different subsampling rates is tested; the results are depicted in Table 3. As shown in Table 3, we can find that the quality of the recovered image can be improved as the subsampling rate increases. equivalence to the random Gaussian measurement operator. Numerical experiments are shown in Figure 9. In Figure 9, the parameters are set as follows: Cameramen image is selected as the test image, the size is 512 × 512, the subsampling rate is 0.5, the measurements are complex, the type of subsampling is random location. The sparsifying basis is RWT with a wavelet scale is 3. The performance of CI-DRPE under different subsampling rates is tested; the results are depicted in Table 3. As shown in Table 3, we can find that the quality of the recovered image can be improved as the subsampling rate increases.

Conclusions
In this paper, we have demonstrated the effectiveness of several different compressive imaging methods using random phase modulation. The principles of compressive sensing theory and random phase modulation are presented. A brief review of the implementation of compressive imaging with different optical architectures is described. Numerical experiments demonstrate that compressive imaging together with random phase modulation can help in recovering the image from real-valued measurements or complexvalued measurements.
In summary, we have shown that random phase modulation can effectively enhance the performance of the compressive imaging method. We believe that random phase modulation combined with compressive sensing can further enhance the performance of three-

Conclusions
In this paper, we have demonstrated the effectiveness of several different compressive imaging methods using random phase modulation. The principles of compressive sensing theory and random phase modulation are presented. A brief review of the implementation of compressive imaging with different optical architectures is described. Numerical experiments demonstrate that compressive imaging together with random phase modulation can help in recovering the image from real-valued measurements or complex-valued measurements.
In summary, we have shown that random phase modulation can effectively enhance the performance of the compressive imaging method. We believe that random phase modulation combined with compressive sensing can further enhance the performance of three-dimensional imaging, video sensing, and phase retrieval, which will lead to new designs of optical setups and enable new applications.