US 20060239336 A1 Abstract A new digital image/video camera that directly acquires random projections of the incident light field without first collecting the pixels/voxels. In one preferred embodiment, the camera employs a digital micromirror array to perform optical calculations of linear projections of an image onto pseudorandom binary patterns. Its hallmarks include the ability to obtain an image with only a single detection element while measuring the image/video fewer times than the number of pixels or voxels—this can significantly reduce the computation required for image/video acquisition/encoding. Since the system features a single photon detector, it can also be adapted to image at wavelengths that are currently impossible with conventional CCD and CMOS imagers.
Claims(29) 1. A method for acquiring an image or video sequence from incoherent projections comprising the steps of:
modulating an incident light field by a series of patterns; optically computing inner products between the light field and said series of patterns; and reconstructing an image or video sequence based upon said inner products. 2. A method for acquiring an image or video sequence from incoherent projections according to 3. A method for acquiring an image or video sequence according to 4. A method for acquiring an image or video sequence according to 5. A method for acquiring an image or video sequence according to 6. A method for acquiring an image or video sequence according to 7. A method for acquiring an image or video sequence according to taking additional measurements; and using said additional measurements for subsequent processing. 8. A method for acquiring an image or video sequence according to 9. An imaging system comprising:
means for modulating an incident light field by a series of patterns; means for optically computing inner products between the light field and a series of patterns; and means for recovering a signal based upon said inner products. 10. An imaging system according to 11. An imaging system according to 12. An imaging system according to a single sensing element; and a processor for applying a reconstruction algorithm to signals received from said single sensing element. 13. An imaging system according to 14. An imaging system according to 15. An imaging system according to a plurality of sensing elements; and a processor for applying a reconstruction algorithm to signals received from said plurality of sensing elements. 16. An imaging system according to 17. An imaging system according to 18. An imaging device comprising:
a lens; a digital micromirror device array; means for modulating mirror orientations of said digital micromirror device array in an incoherent pattern sequence; a pattern generator for supplying incoherent patterns to said means for modulating; a sensor element, wherein each different mirror pattern produces a voltage at the sensor element that corresponds to a measurement; an analog-to-digital converter for quantizing said voltage; and a reconstruction means for reconstructing an image based upon a bitstream received from said analog-to-digital converter. 19. An imaging device according to 20. An imaging device according to 21. An imaging device according to _{1 }reconstruction, l_{0 }reconstruction, and l_{2 }reconstruction. 22. An imaging device according to 23. An imaging device according to 24. An imaging device according to 25. An imaging device according to 26. An imaging device comprising:
a sensing element; and means for focusing light onto said sensing element for measurement; wherein measurements taken by said sensing element correspond to inner products of an incident image with a sequence of incoherent patterns. 27. A method of acquiring a signal comprising the steps of:
measuring and encoding projections of a signal that is compressible in a first set of basis functions onto a second set of basis functions, wherein said second set of basis functions is incoherent with said first set of basis functions; reconstructing said signal based upon said projections. 28. A method for acquiring a signal according to 29. A method for acquiring a signal according to _{1 }reconstruction, l_{0 }reconstruction, and l_{2 }reconstruction.Description The present application claims the benefit of the filing dates of U.S. Provisional Application Ser. No. 60/673,364 entitled “Method and Apparatus for OpticalImage Compression,” and filed on Apr. 21, 2005; U.S. Provisional Application Ser. No. 60/679,237 entitled “Method and Apparatus for Reconstructing Data from Multiple Sources,” and filed on May 10, 2005; U.S. Provisional Application Ser. No. 60/729,983 entitled “Random Filters for Compressive Sampling and Reconstruction,” and filed on Oct. 25, 2005; U.S. Provisional Application Ser. No. 60/732,374 entitled “Method and Apparatus for Compressive Sensing for Analog-to-Information Conversion,” and filed on Nov. 1, 2005; U.S. Provisional Application Ser. No. 60/735,616 entitled “Method and Apparatus for Distributed Compressed Sensing,” and filed on Nov. 10, 2005; and U.S. Provisional Application Ser. No. 60/759,394 entitled “Sudocodes: Efficient Compressive Sampling Algorithms for Sparse Signals,” and filed on Jan. 16, 2006. The above cross-referenced related applications are hereby incorporated by reference herein in their entirety. Supported by NSF grant CCF-0431150, ONR grant N00014-02-1-0353, and AFOSR grant FA9550-04-1-014. 1. Field of the Invention The invention relates to imaging devices such as cameras, video cameras, microscopes, and other visualization techniques, and more particularly, to the acquisition of images and video using fewer measurements than previous techniques. 2. Brief Description of the Related Art The large amount of raw data acquired in a conventional digital image or video often necessitates immediate compression in order to store or transmit that data. This compression typically exploits a priori knowledge, such as the fact that an N-pixel image can be well approximated as a sparse linear combination of K<<N wavelets. These appropriate wavelet coefficients can be efficiently computed from the N pixel values and then easily stored or transmitted along with their locations. Similar procedures are applied to videos containing F frames of P pixels each; we let N=FP denote the number of video “voxels”. This process has two major shortcomings. First, acquiring large amounts of raw image or video data (large N) can be expensive, particularly at wavelengths where CMOS or CCD sensing technology is limited. Second, compressing raw data can be computationally demanding, particularly in the case of video. While there may appear to be no way around this procedure of “sample, process, keep the important information, and throw away the rest,” a new theory known as Compressive Sensing (CS) has emerged that offers hope for directly acquiring a compressed digital representation of a signal without first sampling that signal. See Candès, E., Romberg, J., Tao, T., “Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information,” IEEE Trans. Inform. Theory 52 (2006) 489-509; David Donoho, “Compressed sensing,” IEEE Transactions on Information Theory, Volume 52, Issue 4, April 2006, Pages: 1289-1306; and Candès, E., Tao, T., “Near optimal signal recovery from random projections and universal encoding strategies,” (2004) Preprint. Traditional methods of conserving power in camera monitoring and surveillance applications have either relied upon scheduling sleeping and awake modes, or supplementary sensors such as infrared motion detectors to decide when to power on the camera. In the former case, scheduled power-off periods could result in missing an important event entirely. In the latter case, we require additional hardware that may be costly or undesirable. Moreover, in both cases the system suffers from a “power-on lag,” which delays image or video capture, potentially causing the camera to miss the important event. These problems would be solved by allowing the camera to continuously monitor the scene in a low-power, low-rate mode, and by enabling it to immediately increase its rate when an important or interesting event occurs. This kind of scheme is impossible in the traditional digital camera paradigm, which is an all-or-nothing scheme: either an image/video is captured at full rate, or no image/video is captured at all. Thus a camera that can continuously monitor at low-rate and increase to full rate with no lag-time is not found in the art, but is directly enabled by our unique camera architecture. Other efforts on compressed imaging include Pitsianis, N. P., Brady, D. J., Sun, X.: “Sensor-layer image compression based on the quantized cosine transform,” SPIE Visual Information Processing XIV (2005) and Brady, D. J., Feldman, M., Pitsianis, N., Guo, J. P., Portnoy, A., Fiddy, M., “Compressive optical MONTAGE photography,” SPIE Photonic Devices and Algorithms for Computing VII (2005), which employ optical elements to perform transform coding of multispectral images. The hardware designed for these purposes uses concepts that include optical projections, group testing (see Cormode, G., Muthukrishnan, S., “Towards an algorithmic theory of compressed sensing,” DIMACS Tech. Report 2005-40 (2005)), and signal inference. Two notable previous DMD-driven applications involve confocal microscopy (Lane, P. M., Elliott, R. P., MacAulay, C. E., “Confocal microendoscopy with chromatic sectioning,” Proc. SPIE. Volume 4959 (2003) 23-26) and micro-optoelectromechanical (MOEM) systems (DeVerse, R. A., Coifman, R. R., Coppi, A. C., Fateley, W. G., Geshwind, F., Hammaker, R. M., Valenti, S., Warner, F. J., “Application of spatial light modulators for new modalities in spectrometry and imaging,” Proc. SPIE. Volume 4959 (2003)). The present invention overcomes shortcomings of the prior approaches. Preferred embodiments of the present invention take fewer measurements than prior techniques, enable significant reduction in the resources (power, computation) required for visualization and use only a small number of physical sensors. The reduction in the size of the hardware associated with preferred embodiments of the invention further may significantly reduce costs of visualization systems. The present invention can also acquire and process streaming video data (time-varying images). Finally, the present invention can adjust its data acquisition rate according to the amount of activity in the scene it is imaging. The present invention uses algorithms and hardware to support a new theory of Compressive Imaging (CI). The approach is based on a new digital image/video camera that directly acquires random projections without first collecting the N pixels/voxels. (See Takhar, D., Laska, J. N., Wakin, M., Duarte, M., Baron, D., Sarvotham, S., Kelly, K. K., Baraniuk, R. G., “A new camera architecture based on optical-domain compression,” Proc. IS&T/SPIE Symposium on Electronic Imaging: Computational Imaging. Volume 6065. (2006)). Due to this unique measurement approach, it has the ability to obtain an image with a single detection element while measuring the image far fewer times than the number of pixels/voxels. Note also that additional embodiments using a plurality of detection elements can also be used. The image can be reconstructed, exactly or approximately, from these random projections by using a model, in essence to find the best or most likely image (in some metric) among all possible images that could have given rise to those same measurements. While several preferred embodiments of reconstruction are described below, it should be understood that additional techniques using or incorporating the present invention can also be used. A small number of detectors, even a single detector, can be used. Thus, the camera can be adapted to image at wavelengths of electromagnetic radiation that are currently impossible with conventional CCD and CMOS imagers. This feature is particularly advantageous, because in some cases the usage of many detectors is impossible or impractical, whereas the usage of a small number of detectors, or even a single detector, may become feasible using compressive imaging. A camera in accordance with the present invention can also be used to take streaming measurements of a video signal, which can then be recovered using CS techniques designed for either 2-dimensional (2D) frame-by-frame reconstruction or joint 3D reconstruction. This allows a significant reduction in the computational complexity of the video encoding process. An imaging system in accordance with the present invention enjoys a number of desirable features: -
- Potentially single detector or small number of detectors: By time-multiplexing each detector, we can use a less expensive and yet more sensitive photon detectors. This is particularly important when the detector is expensive, making an N-pixel array prohibitive. A single detector camera can also be adapted to image at wavelengths that are currently impossible with conventional CCD and CMOS imagers.
- Universality: Random and pseudorandom measurement schemes are universal in the sense that they can be paired with any signal model. Therefore, the same encoding strategy can be applied in a variety of different sensing environments; knowledge of the nuances of the environment is needed only at the reconstruction mechanism (decoder). Random measurements are also future-proof: if future research in image processing yields a better signal model then the same set of random measurements can be used to reconstruct an even better quality image or video.
- Encryption: A pseudorandom sequence can be generated using a simple algorithm according to a random seed. Such encoding effectively implements a form of encryption: the randomized measurements will themselves resemble noise and cannot be decoded unless an observer knows the associated seed.
- Robustness and progressivity: Random and pseudorandom measurements are robust in the sense that the measurements have equal priority, unlike the Fourier or wavelet coefficients in current transform coders. Thus they allow a progressively better reconstruction of the data as more measurements are obtained; one or more measurements can also be lost without corrupting the entire reconstruction.
- Scalability: We can adaptively select how many measurements to compute in order to trade off the amount of compression of the acquired image/video versus acquisition time; in contrast, conventional cameras trade off resolution versus the number of pixel sensors.
- Computational asymmetry: compressive imaging (CI) places most of its computational complexity in the decoder, which will often have more substantial computational resources than the encoder/imager. The encoder is very simple; it merely computes incoherent projections and, depending on the specific embodiment, makes few or no decisions.
Still other aspects, features, and advantages of the present invention are readily apparent from the following detailed description, simply by illustrating preferable embodiments and implementations. The present invention is also capable of other and different embodiments and its several details can be modified in various obvious respects, all without departing from the spirit and scope of the present invention. Accordingly, the drawings and descriptions are to be regarded as illustrative in nature, and not as restrictive. Additional objects and advantages of the invention will be set forth in part in the description which follows and in part will be obvious from the description, or may be learned by practice of the invention. For a more complete understanding of the present invention and the advantages thereof, reference is now made to the following description and the accompanying drawings, in which: FIGS. In a first preferred embodiment, a camera architecture of the present invention uses for random measurements a digital micromirror array to spatially modulate an incident image and reflecting the result to a lens, which focuses the light to a single photodiode for measurement. Mathematically, these measurements correspond to inner products of the incident image with a sequence of pseudorandom patterns. For an image model the system assumes sparsity or compressibility; that is, that there exists some basis, frame, or dictionary (possibly unknown at the camera) in which the image has a concise representation. For reconstruction, this system and method uses the above model (sparsity/compressibility) and some recovery algorithm (based on optimization, greedy, iterative, or other algorithms) to find the sparsest or most compressible or most likely image that explains the obtained measurements. The use of sparsity for signal modeling and recovery from incomplete information are the crux of the recent theory of Compressive Sensing (CS), explained below. The camera, however, does not have to rely on reflecting light off a digital micromirror device as in Compressive Sensing Transform Coding Compressive Sensing (CS) builds upon a core tenet of signal processing and information theory: that signals, images, and other data often contain some type of structure that enables intelligent representation and processing. Current state-of-the-art compression algorithms employ a decorrelating transform to compact a correlated signal's energy into just a few essential coefficients. Such transform coders exploit the fact that many signals have a sparse representation in terms of some basis A, meaning that a small number K of adaptively chosen transform coefficients can be transmitted or stored rather than N signal samples, where K<N. Mathematically, we wish to acquire an N-sample signal/image/video x for which a basis or (tight) frame Ψ=[ψ We use the terms “signal” and “image” interchangeably here, since an image is two-dimensional (2D) signal. A video sequence is a sequence of images, or a 3D signal. We use the notations φ The standard procedure for transform coding of sparse signals is to (i) acquire the full N-sample signal x; (ii) compute the complete set {θ(n)} of transform coefficients θ(i)=<x, ψ(i)>, where <·,·> denotes the inner product, θ(i) denotes the i'th coefficient, and ψ(i) denotes the i'th basis vector (i'th column of the matrix Ψ); (iii) locate the K largest, significant coefficients and discard the (many) small coefficients; and (iv) encode the values and locations of the largest coefficients. In cases where N is large and K is small, this procedure is quite inefficient. Much of the output of the analog-to-digital conversion process ends up being discarded (though it is not known a priori which pieces are needed). This raises a simple question: For a given signal, is it possible to directly estimate the set of large coefficients that will not be discarded by the transform coder? While this seems improbable, the recent theory of Compressive Sensing introduced by Candès, Romberg, and Tao and Donoho referenced above demonstrates that a signal that is K-sparse in one basis (call it the sparsity basis) can be recovered from cK nonadaptive linear projections onto a second basis (call it the measurement basis) that is incoherent with the first, where c is a small overmeasuring constant. While the measurement process is linear, the reconstruction process is decidedly nonlinear. Incoherent Projections In CS, we do not measure or encode the K significant θ(n) directly. Rather, we measure and encode M<N projections y(m)=<x, φ The CS theory tells us that when certain conditions hold, namely that the basis cannot sparsely represent the elements of the sparsity-inducing basis (a condition known as incoherence of the two bases) and the number of measurements M is large enough, then it is indeed possible to recover the set of large {θ(n)} (and thus the signal x) from a similarly sized set of measurements {y(m)}. This incoherence property holds for many pairs of bases, including for example, delta spikes and the sine waves of the Fourier basis, or the Fourier basis and wavelets. Significantly, this incoherence also holds with high probability between an arbitrary fixed basis and a randomly generated one (consisting of i.i.d. Gaussian or Bemoulli/Rademacher ±1 vectors). Signals that are sparsely represented in frames or unions of bases can be recovered from incoherent measurements in the same fashion. We call the rows of Φ the measurement basis, the columns of Ψ the sparsity basis or sparsity inducing basis, and the columns of V=ΦΨ=[V Signal/Image Recovery The recovery of the sparse set of significant coefficients {θ(n)} can be achieved using optimization or other algorithms by searching for the signal with l Unfortunately, solving this optimization problem is prohibitively complex and is believed to be NP-hard (see Candès, E., Tao, T., “Error correction via linear programming,” (2005) Preprint). The practical revelation that supports the new CS theory is that it is not necessary to solve the l The optimization problem (2), also known as Basis Pursuit (see Chen, S., Donoho, D., Saunders, M., “Atomic decomposition by basis pursuit,” SIAM J. on Sci. Comp. 20 (1998) 33-61), is significantly more approachable and can be solved with traditional linear programming techniques whose computational complexities are polynomial in N. Although only K+1 measurements are required to recover sparse signals via l We use the notation c to describe the overmeasuring/oversampling constant required in various settings and note the following approximation: The constant c satisfies c≈log 2 (1+N/K). While reconstruction based on linear programming is one preferred embodiment, any reconstruction approach can be used in the present invention. Other examples include the (potentially more efficient) iterative Orthogonal Matching Pursuit (OMP) (see Tropp, J., Gilbert, A. C., “Signal recovery from partial information via orthogonal matching pursuit,” (2005) Preprint), matching pursuit (MP)(see Mallat, S. and Zhang, Z., “Matching Pursuit with Time Frequency Dictionaries”, (1993) IEEE Trans. Signal Processing 41(12): 3397-3415), tree matching pursuit (TMP) (see Duarte, M. F., Wakin, M. B., Baraniuk, R. G., “Fast reconstruction of piecewise smooth signals from random projections,” Proc. SPARS05, Rennes, France (2005)) algorithms, group testing (see Cormode, G., Muthukrishnan, S., “Towards an algorithmic theory of compressed sensing,” DIMACS Tech. Report 2005-40 (2005), Sudocodes (see U.S. Provisional Application Ser. No. 60/759,394 entitled “Sudocodes: Efficient Compressive Sampling Algorithms for Sparse Signals,” and filed on Jan. 16, 2006), or statistical techniques such as Belief Propagation, (see Pearl, J., “Fusion, propagation, and structuring in belief networks”, (1986) Artificial Intelligence, 29(3): 241-288), LASSO (see Tibshirani, R., “Regression shrinkage and selection via the lasso”, (1996) J. Royal. Statist. Soc B., 58(1): 267-288), LARS (see Efron, B., Hastie, T., Johnstone, I., Tibshirani, R., “Least Angle Regression”, (2004) Ann. Statist. 32(2): 407-499), Basis Pursuit with Denoising (see Chen, X., Donoho, D., Saunders, M., “Atomic Decomposition by Basis Pursuit”, (1999), SIAM Journal on Scientific Computing 20(1): 33-61), expectation-maximization (see Dempster, Laird, N., Rubin, D., “Maximum likelihood from incomplete data via the EM algorithm”, (1997) Journal of the Royal Statistical Society, Series B, 39(1): 1-38), and so on. These methods have also been shown to perform well on compressible signals, which are not exactly K-sparse but are well approximated by a K-term representation. Such a model is more realistic in practice. Reconstruction can also be based on other signal models, such as manifolds (see Wakin, M, and Baraniuk, R., “Random Projections of Signal Manifolds” IEEE ICASSP 2006, May 2006, to appear). Manifold models are completely different from sparse or compressible models. Reconstruction algorithms in this case are not necessarily based on sparsity in some basis/frame, yet signals/images can be measured using the systems described here. The systems described here can also be used to acquire a collection of images or video sequences. Each image or video can be viewed as a point in N-dimensional Euclidean space. Therefore, the collection of images/videos forms a point cloud in N dimensional Euclidean space. Incoherent projections as implemented in our systems will keep different images/videos well-separated and preserve the neighborhood relationships among similar signals, even if we never intend to reconstruct these images/videos (see Dasgupta, S., Gupta, A., “An elementary proof of the Johnson-Lindenstrauss lemma,” Tech. Rep. TR-99-006, Berkeley, Calif., 1999). The point cloud approach is useful for posing and solving decision problems with collections of images/videos, such as detection, classification, recognition, tracking, registration, and other problems. The preferred embodiment is to reconstruct an N-pixel image or video sequence from M<N measurements. Additional embodiments using more measurements are possible. For example, if we use M=N or M>N measurements, then the extra measurements can be used for subsequent processing. For example, additional measurements may be used for averaging or filtering when the image is noisy or corrupted in some way. Compressive Imaging The present invention is a new system to support what can be called Compressive Imaging (CI). In one preferred embodiment, the present invention incorporates a microcontrolled mirror array driven by pseudorandom and other measurement bases and a single or multiple photodiode optical sensor. This hardware optically computes incoherent image measurements as dictated by the CS theory; CS reconstruction algorithms are then applied to obtain the acquired images. A camera in accordance with the present invention can also be used to take streaming measurements of a video signal, which can then be recovered using CS techniques designed for either 2D frame-by-frame reconstruction or joint 3D reconstruction. Streaming video can also be supported. Other desirable features of our system include the use of a single detector (potentially enabling imaging at new wavelengths that are currently impossible or infeasible with CCD and CMOS technology), universal measurement bases (incoherent with arbitrary sparse bases), encrypted measurements (tied to a random seed that can be kept secure), and scalable progressive reconstruction (yielding improved quality with more measurements). (See Takhar, D., Laska, J. N., Wakin, M., Duarte, M., Baron, D., Sarvotham, S., Kelly, K. K., Baraniuk, R. G., “A new camera architecture based on optical-domain compression,” Proc. IS&T/SPIE Symposium on Electronic Imaging: Computational Imaging. Volume 6065. (2006).) Camera Hardware One possible hardware realization of the CI concept is a single detector camera; it combines a microcontrolled mirror array displaying a time sequence of M pseudorandom basis images with a single optical sensor to compute incoherent image measurements y as in (1) (see A preferred embodiment of the invention employs a Texas Instruments digital micromirror device (DMD) for generating the random modulation basis patterns. The DMD consists of a 1024×768 array of electrostatically actuated micromirrors where each mirror of the array is suspended above an individual SRAM cell. Each mirror rotates about a hinge and can be positioned in one of two states (+12 degrees and −12 degrees from horizontal); thus light falling on the DMD may be reflected in two directions depending on the orientation of the mirrors. Note that the Texas Instruments DMD is one possible embodiment, but many additional embodiments are possible. Referring again to The photodiode measurements can be interpreted as the inner product of the desired image x with a measurement basis vector φ Equation (3) holds the key for implementing a compressive imaging (CI) system. For a given incident lightfield x, we take M measurements {y(1), y(2), . . . , y(M)} corresponding to mirror configurations {ρ(1), ρ(2), . . . , ρ(M)} and thus measurement basis vectors {φ Mirrors can also be duty-cycled to give the elements of Φ finer precision, for example to approximate Gaussian measurement vectors (see D. Donoho, “Compressed Sensing,” IEEE Transactions on Information Theory, Volume 52, Issue 4, April 2006, Pages: 1289-1306; and Candès, E., Tao, T., “Near optimal signal recovery from random projections and universal encoding strategies,” (2004) Preprint). This duty-cycling technique can be used to emulate inner products with any real-valued vector. Specific embodiments may generate each coefficient of such projection vectors using some continuous probability distribution, but any set of real-valued vector values can be used. This compressive imaging system directly acquires a reduced set of M incoherent projections of an N-pixel image x without first acquiring the N pixel values. Since the camera is “progressive,” better quality images (larger K) can be obtained by taking a larger number of measurements M. Also, since the data measured by the camera is “future-proof,” new reconstruction algorithms based on better sparsifying image transforms can be applied at a later date to obtain even better quality images. Streaming Video Acquisition The CI system and method of the present invention is immediately applicable to video acquisition. As described above, the measurements are taken sequentially in time. Hence, one can view each measurement as a linear projection against a snapshot of the scene at that instant. Viewing the video as a three-dimensional (3D) signal (in which the 2D snapshots are stacked), the measurements vectors {φ(m)} themselves are each localized onto a different 2D snapshot for each m. In order to recover a video sequence from these measurements, some simplifying assumptions are made. Specifically, traditional CS considers an ensemble of measurements taken from a single signal; in the streaming setting, however, each measurement will act on a different snapshot. The present invention overcomes this problem by assuming that the image changes slowly across a group of snapshots, which can then be equated to a single video frame. The number of snapshots assigned to a frame will be determined by the speed of the acquisition system and the desired temporal resolution of the reconstructed video. Under this assumption, the acquired video is represented as a sequence of F frames, each one measured using M/F measurement vectors that we can group as rows of a matrix Ψ Several options exist for reconstructing the video from these measurements. First, one could reconstruct each frame individually using 2D wavelets, performing a total of F CI reconstructions. Each reconstruction would use the same 2D wavelet sparsity basis Ψ but with a different measurement matrix. This process will be referred to herein as frame-by-frame reconstruction. Alternative methods more fully exploit the correlation between frames. One solution is to use 3D wavelets as a sparse representation for the video sequence; i.e., to define the joint measurement matrix
The compressive imaging architecture and method of the present invention can also be extended to acquire full 3D measurements of a video sequence (that is, where each has 3D support). One embodiment of such 3D measurements would combine inner products sampled at different times, but other embodiments are possible. Under this setting with full 3D measurements, the entire video sequence is reconstructed using a single measurement matrix Φ that operates on all of the frames and a suitable 3D sparse basis Ψ such as wavelets. Below, it is demonstrated that such a scheme enjoys better incoherence with the video structure. However, it also increases the complexity of both the measurement and reconstruction processes. Possible solutions to this increased complexity include partitioning the video into blocks, which are then reconstructed separately. For an imaging experiment, we displayed a printout of the letter “R” in front of the camera; To demonstrate the potential for applications in video encoding, we present a series of simulations for video measurement/reconstruction. As shown in Adaptive Compressive Imaging Scheme An alternative adaptive compressive imaging scheme embodiment that takes more and more measurements until image/video reconstruction is possible is described. As shown in Attention Based Processing One way to save power is to use the camera only when it is needed. Blindly turning the power off except for periodic wakeup intervals is easily accommodated in our design, but this is too crude for many modern sensing applications. Another approach is to monitor instances of change. For many applications, the received image light field may not change for long periods of time, or the images may change very slowly, allowing a slower data acquisition rate. For example, a camera used to monitor an infrequently used stairwell does not really need to report the same image of the stairwell over and over. However, if a burglar enters the stairwell, then the camera will need to awake abruptly when the scene changes. To conserve the power consumed by the analog-to-digital converter, the computation determining when the camera should wake up should be as simple as possible, using low-power analog processing and/or low-rate digital processing. Essentially, the camera needs a startle reflex: a quick reaction mode involving little overhead or computation. The CI camera is ideally suited to implement a startle reflex. Simply put, the values or statistics of coefficients at the analog-to-digital converter output, running at a very slow, low-power rate, would be compared with previous values or statistics to determine whether the scene being monitored has changed. Algorithms for achieving this could range from very simple statistical (parametric and non-parametric) tests, for example a test based on an energy detector, a test based on empirical entropies (see Gutman, M., “Asymptotically Optimal Classification for Multiple Tests with Empirically Observed Statistics,” IEEE Trans. Inform. Theory 35 (1989) 401-408), or more sophisticated tests based on detailed models of motion. It is also possible that measurement systems other than pseudo-random, ones more adapted to natural images, may give rise to more predictable distributions of values for natural images; thus, any such algorithms could be enhanced by using a carefully designed basis in the analog-to-digital converter. Moreover, rather than simply detecting a change, different bases could be exploited to provide a more detailed analysis of different kinds of change, allowing a richer set of decision rules to control the camera's rate. A simple embodiment (outlined below) based on the empirical entropy decision rule developed by Gutman (see Gutman above) is described. Experiments have shown that even small changes in images can be detected at an extremely low sampling rate, much lower than would be required to reconstruct the image for viewing. Thus, using only a few (potentially highly quantized) coefficients from the analog-to-digital converter, one can implement an effective startle reflex leading to drastic power savings in applications like monitoring and surveillance. Other statistical tests are also possible beyond the Gutman test. The startle-reflex algorithm may be described as follows: -
- 1. Collect S samples of a sequential measurement vector y, quantized to B bits per measurement. This is the training sample.
- 2. Compute the empirical probability distribution of the training sample.
- 3. Collect the next S samples of y. This is the testing sample.
- 4. Compute the empirical probability distribution of the testing sample.
- 5. Compute the empirical entropy test statistic (see Gutman above).
- 6. If the test statistic exceeds a threshold T, then increase camera to full sampling rate and capture image. Replace training sample with testing sample.
- 7. Otherwise, form a new training sample by concatenating the old training sample with the testing sample.
- 8. Repeat steps 2-8.
The startle-reflex algorithm above can be implemented quite simply and demands very little by way of sensing and computation. Yet, it has proven to be highly effective at detecting changes in sequences of images, even at very low measurement rates. For example, typical values of the parameters examined during simulations were: S=30 samples, 6 coefficients in y, and 3 bit (8 level) uniform scalar quantization of each coefficient. These values yielded good performance in many sequences of test images, correctly rejecting up to 50% of the images when no change had occurred These specific values and the startle-reflex algorithm above are specific embodiments of attention based processing. Additional embodiments may use different parameter values and/or test statistic algorithms differing from that of Gutman (see Gutman above). Attention based processing is not limited to turning a CI camera on and off. Additional embodiments may use the attention information for additional purposes, for example to track different phenomena in space and/or time. Micromirror Modulation for Imaging Current imaging devices rely on CCD or CMOS technology for the optical sensing element. While the scale and cost of this technology is continually reducing, the complexity and power requirements have not similarly scaled. Digital micromirror devices have proven to be a commercially viable MEMs technology for the video/projector display market (see D. Doherty and G. Hewlett, “Phased reset timing for improved digital micromirror device (DMD) brightness,” in Photodiode Sensing Element: By replacing the optical sensor array with a single sensing element (in this case a photodiode), we have greatly reduced the complexity. In shifting the complexity of the camera away from the signal receiving portion and into the signal acquisition/interpretation areas, we are able to work with less expensive and more sensitive photon detectors. The advantages of a photodiode sensing element include low-power, low cost, high photon sensitivity that increases with the diode volume, and very fast response times. Modern photodiodes are routinely operated at hundreds of megahertz and have been extended into the GHz regime. The cost of photodiodes can be as little as $0.50 a chip with the cost increasing with the performance capabilities; still, a very good quality photodiode is around $15. In addition, having one optical receiving element allows us to increase the detection efficiency by increasing its size. Additional applications of this technology occur in the cases where power and cost are not limited. The single photodiode can be replaced with a quadrant photodiode for higher spatial resolution. In specialty applications where greater photon sensitivity may be preferred and cost is not a consideration, the normal photodiode may be replaced with an avalanche photodiode to allow single photon counting. For photodiodes, a variety of semiconductor materials are available, allowing the sensitivity to span the optical, UV, and IR spectrums with equal capabilities. While the initial prototype will be was grayscale, the conversion of this device to a full color camera has been straightforward. For proof of concept we have implemented color imaging in our architecture with RGB filters mounted on a color wheel. However, many other color technologies may also be adapted with our imaging architecture. Our realization of these technologies is due to the broad spectral response of the micromirror and photodiode. Many possible embodiments exist for full-color implementation including a series of prisms to separate the signal between 3 separate photodiodes. In a similar manner we can easily extend the capabilities of our camera for more detailed multispectral or hyperspectral imaging. Any image signal that can be shuttered by an incoherent pattern and then summed at the detector is suitable acquisition and reconstruction by our proposed methods. Color image/video reconstruction can be facilitated by the fact that the color channels share common information (they are correlated or “jointly sparse”). Therefore, the techniques of distributed compressed sensing could be used to lower the number of measurements required for color data acquisition (see Baron, D., Wakin, M. B., Duarte, M. F., Sarvotham, S., Baraniuk, R. G., “Distributed compressed sensing” (2005). There are also alternative embodiments of photodiodes, some of which we describe here. Photodiodes can be implemented in a variety of circuit configurations depending on the application. Output voltage could be set up to be proportional to logarithmic change in the detected light level, the external circuit could optimized to emphasize the converting the light signal to frequency, or an x-ray scintillation detector could be mounted in front of the photodiode for medical or astronomy applications (with the appropriate modifications to the mirror coating). These are specific embodiments; additional specific embodiments of photodiodes are possible. Digital Micromirror Devices: The Texas Instruments Digital Micromirror Device (DMD) is composed of an array of electrostatically actuated micromirrors that has found a great deal of success in the projection screen market (see D. Doherty and G. Hewlett, “Phased reset timing for improved digital micromirror device (DMD) brightness,” in These devices are well above and beyond the necessary requirements of pixel density and contrast ratio for our proposed device. However, any micromirror array, either electrostatically or piezoelectrically driven, is suitable for our camera architecture. Similar driving mechanisms may also be suitable for our camera architecture. Piezoelectric Deformable Mirror: As an alternative to electrostatic manipulation in MEMs devices, piezoelectric materials offer a similar ability to convert electrical signals into mechanical work. At the same time, they are able to actuate at much greater frequencies compared to electrostatic based systems. A piezoelectric transducer can reach its nominal displacement in ⅓ of the period of the resonant frequency. Depending of the final size and scale of the tranducers, these frequencies could correspond to a time on the order of microseconds. Electrostatic structures are still favored in MEMs applications due to the ease of the incorporation with traditional IC manufacturing technologies. The overall goal of our micromirror system is not increased miniaturization and its correspondingly greater pixel density, although a commercial version of the DMD based on piezoelectric has been explored (see Y. J. S. K. H. Hwang and S. G. Kim, “Thin-film micromirror array for high-brightness projection displays,” We envision two possible embodiments between the micromirror 710 and photodiode The second, off-center configuration of the switching mirror illustrated in Another possible embodiment includes a microcontroller that drives the mirror motions in such a manner that the mirror surface structure forms a time-varying 2D smooth surface. Controlling the mirror structure to conform to the desired smooth surface will enable the mirror angle to vary smoothly between spatially close-by locations on the mirror surface. Therefore, in contrast to the measurement via discretely-computed inner products mentioned before, in this system the device will sense the continuous integral of the optical signal modulated by the mirror surface waveform. This capability will enable advanced analog measurement techniques. This integral can be written as an inner product not between two discrete, length-N vectors but between two continuously varying 2D functions. Another possible embodiment is to perform image acquisition using real-space convolution with white-noise as a shuttering mechanism and recorded by an individual detector. Previous embodiments were presented as reflecting light. Additional embodiments can reflect other frequencies and even electron imaging. In another embodiment, an image may be formed using a similar modulation of a (potentially micro-electromechanical) shutter array placed directly over the detector. This would create an essentially flat camera. In fact, the modulation mechanism of our image signal in transmission mode would apply well beyond the optical regime allowing for construction of a camera out of a single sensor in regimes where reflective optical elements do not exist, such as gamma rays. Compressive imaging can be incorporated in distributed systems with multiple imaging devices. This will enable reconstruction of multiple images (e.g., multiple frames of a video sequence) using fewer measurements than before, requiring reduced resource consumption (e.g. power reduction). Alternatively, these techniques could enable better image reconstruction quality. Power reduction can be achieved by minimizing the amount of mirror motion. One way to do this is to specifically design measurement matrices such that adjacent rows are as similar as possible to one another. The foregoing description of the preferred embodiment of the invention has been presented for purposes of illustration and description. It is not intended to be exhaustive or to limit the invention to the precise form disclosed, and modifications and variations are possible in light of the above teachings or may be acquired from practice of the invention. The embodiment was chosen and described in order to explain the principles of the invention and its practical application to enable one skilled in the art to utilize the invention in various embodiments as are suited to the particular use contemplated. It is intended that the scope of the invention be defined by the claims appended hereto, and their equivalents. The entirety of each of the aforementioned documents is incorporated by reference herein. Referenced by
Classifications
Legal Events
Rotate |