SE538072C2 - An imaging system parallelizing compressive sensing imaging - Google Patents

Description

25 30 35 538 072 herby incorporated by reference. This has also been demonstrated; please cf.

Howland, G. A., Dixon, P. B., & Howell, J. C. (2011): Photon-counting compressive sensing laser radar for 3D imaging. Applied optics, 50(31), 5917-5920, herby incorporated by reference. ln this technology the 2D detector array in a traditional camera architecture is replaced by a spatial light modulator (SLM), which can e.g. be a digital micro-mirror device (DMD). A pattern applied to the DMD will reflect the light incident on certain pixels towards a lens collecting all the light onto a single detector.

Light incident towards other pixels on the DMD will be directed away from this lens. ln this way a measurement by a single detector will sample a Iinear combination of pixels in the image. A new measurement using a different pattern on the DMD will sample a different linear combination of pixels. lf a number of measurements equal to the number of pixels in the array are performed using patterns that are basis vectors of the space spanned by the array this will produce a linear equation system that can be solved using traditional minimization of the squared error or L2-norm.

The purpose of CS is to reduce the number of measurements that needs to be performed compared to a scanned system. This will produce an underdetermined linear equation system, which has infinitely many solutions. ln CS the fact that most data can be described sparsely in some base is used. A reconstruction base is selected, and the most sparse description, that is the one that could produce the measurement results using the least number of non-zero basis coefficients, is assumed to be correct. The reconstruction basis can be the normal pixel basis or any basis that is suitable for describing the scene in a sparse way, for example different wavelet bases are often suitable for natural scenes in analogy with the jpeg 2000 compression. Different bases should be chosen depending on the type of scene that is imaged. A scene consisting of a few bright points in a dark background, as could happen in thermal imaging, should be described by the pixel basis. A scene consisting of several surfaces with different characteristics should instead be described by a wavelet basis.

Mathematically the process can be described as follows for 2D-imaging. Let f be the scene as it would be seen by a normal camera in the position of the DMD. Let the DMD and the imaginary camera have NZ pixels. The DMD could also be rectangular, but a quadratic array is assumed here for illustrational purposes. Randomly selected patterns for the mirrors can be written as NZ long vectors of zeroes and ones, 10 15 20 25 30 35 538 072 placed as rows of the matrix (D. Conducting M measurements with different patterns for the DMD can then be written as b = Cbf, where b is an M-long column vector containing the measurement results. The scene can also be described as f = Wx, where ll' is a basis matrix containing all basis vectors for the reconstruction basis. lf the pixel basis is used for reconstruction W is the identity matrix. lt is important that Gb and ll' are uncorrelated to each other. This is valid for all reconstruction bases when using randomly generated patterns for the DMD. The Nz-element vectorx is the description of the scene in the reconstruction basis. For CS to be of use x should be a sparse vector with only a small number of non-zero values.

The problem to solve can then be written as b = GDWx = Ax, where A is an M> linear equation system can according to the theory of CS be found by minimization of the L1-norm, which is the sum of the absolute values of all coefficients in x, while keeping the equality. Methods for this and extensions to handle noise in measure- ments include basis pursuit and other similar methods. Functions to perform this minimization are available e.g. in the SPGL1-(Spectral Projected Gradient for L1 minimization) package at http://www.cs.ubc.ca/labs/scl/spql1/2.

One problem of CS is that for high definition imaging the number of measurements needed are not small, thus the sequential measurements using different patterns on the DMD take time. ln addition the reconstruction also becomes very computationally demanding when the equation system becomes large. Kelly et al. have suggested reducing this problem by directing sub-images to different discrete detectors; please cf. Kelly, K. F., Baraniuk, R. G., Mcmackin, L., Bridge, R. F., Chatterjee, S., & Weston, T. H. (2012): Decreasing image acquisition time for compressive imaging devices, hereby incorporated by reference. Baraniuk et al. have further discussed the use of re-imaging optics between the DMD and a smaller detector array to multiply the resolution of the detector; please cf. Baraniuk, R. G., Kelly, K. F., & Woods, G. (2013): Number of pixels in detector arrays using compressive sensing, hereby incorporated by reference. 10 15 20 25 30 35 538 072 The present invention solves the problem of long measurement times in compressed sensing by parallelizing the measurement using astigmatic optics and a linear de- tector array in the way that is evident from the following independent claim. The remaining claims concern advantageous embodiments of the invention.

The invention will in the following be described with reference to the accompanying drawings, in which: Fig. 1 is an illustration of an embodiment of the invention where the scene is imaged onto a spatial light modulator (SLM) using standard imaging optics. The SLM imposes a line pattern mask onto the image. Each row of SLM pixels is then re-imaged onto one pixel of a linear detector array using astigmatic optics and Fig. 2 is an illustration of an embodiment of the invention where the pattern is created by the illumination source and an astigmatic camera lens images the scene onto a linear array detector. ln many more complex imaging systems fabrication of large array detectors is a problem. lt may be simply a problem of manufacturing technology where large detectors would have low yield and very high cost, as for e.g. infrared imaging. lt may also be a problem of complex electronics necessary for every pixel, as in 3D laser radar detectors. ln a linear detector array the electronics can expand to the sides without increasing the pixel pitch along the array dimension. This is of course not possible in a 2D detector. Another situation where 2D detector arrays are difficult is hyper-spectral imaging where the spectrum needs to be resolved in addition to the two spatial dimensions. Here it is common to use a 2D detector for the spectral and one spatial dimension and scan the second spatial dimension. CS using astigmatic optics could improve the efficiency of this setup, by removing the need to scan the slit-shaped field of view.

Current DMD technology allows 1920><108O pixels with 23148 Hz frame rate and 10.8 pm pixel pitch (Texas Instruments chipset 0.95 1080p). The size of DMD arrays is expected to continue to increase. lf the full DMD is used for a single CS measure- ment the number of dimensions will be very high (2073600), causing the need for many measurements and hence slow frame rates. By using a linear detector array 10 15 20 25 30 35 538 072 with 1><1080 pixels and astigmatic optics this is reduced to 1080 CS measurements, each with 1920 dimensions. This is a very reasonable problem size where each reconstructed frame can be collected with fifty to a few hundred DMD patterns, using integration times of 10-200 us for each mirror pattern, and hence a frame rate of around 100 Hz can be achieved for low information content scenes and good illumi- nation conditions. For lower illumination levels longer integration times for each mirror pattern can be used to acquire the signal at the cost of lower frame rates.

There is basically no limit to what integration times can be used, it only depends on the dynamic range of the detector and the light conditions. For an active illumination system multiple laser pulses can be used for the same mirror pattern and the signals added to improve the signal to noise ratio. For moderately complex scenes the compressed sensing algorithm will need a larger number of mirror patterns, but the method may be of advantage compared to classical scanning up to over 50 % of the number of dimensions.

The smaller pixel pitch of the DMD makes long focal length imaging lenses unnecessary, potentially reducing the overall size of the imaging system even with the increased complexity of the CS setup compared to a normal camera. ln a preferred embodiment, illustrated in Fig. 1, suitable for passive imaging, e.g. infrared imaging, but also for active 3D imaging with pulsed laser illumination, the invention is an imaging detector where the varying pattern used for the compressed sensing (CS) processing is applied in the detection system. The imaging system consists of a lens system imaging the scene onto a spatial light modulator (SLM) comprising N> direct the radiation into a further re-imaging system or block the radiation depending on the pixel values in the pattern applied to the SLM. ln a preferred embodiment all P rows would use the same patterns, but different patterns for different rows are also possible. The re-imaging system comprises astigmatic optical elements so that the radiation from each row of N pixels of the SLM is collected onto different pixels in a P pixel Iinear detector array. ln this way P simultaneous measurements are performed for each pattern on the SLM and M patterns will produce data to solve P different underdetermined Iinear equation systems with a M> equation system. 10 15 20 25 30 35 538 072 ln one preferred embodiment the SLM is a digital micro-mirror device (DMD). Other possibilities for the SLM include pixelated liquid crystal cells.

The illustration in Fig. 1 shows an imaging system that studies a field of view 101.

The scene inside the field of view could be illuminated by a light source included in the system, be illuminated by ambient light from e.g. the sun, or the thermal radiation from the objects in the scene can be used as light source. lf a dedicated light source is included this could be e.g. a pulsed laser for 3D-imaging or a super-continuum laser for hyper-spectral imaging. This scene is imaged by optics 102 onto an SLM 103. The optics 102 could be a standard camera lens or a telescope suitable for the wavelength of interest. The optics images a small area 104 onto one position 105 on the SLM and other areas 106 onto other positions of the SLM 107, just like regions of the scene are imaged onto pixels of a CCD detector in a standard camera. A second astigmatic optical system 108 images the radiation reflected from or transmitted by the SLM 103 onto a linear detector array 109. The SLM is used to create patterns of vertical lines 110 on the SLM 103 where all or none of the radiation is directed towards the linear detector array 109 based on if that line on the SLM is assigned 1 or 0 in the pattern mask. The astigmatic optical system 108 images slit like regions, e.g. 111 and 113 of the SLM, that are crossing the stripe pattern 110, onto different pixels, 112 and 114 respectively, on the linear detector array 109. Different patterns 110 are used sequentially with one detector reading taken for each pattern to produce a dataset than can be used in compressed sensing reconstruction of the scene. The data from each pixel in the linear detector array produces the image of one line in the scene and these linear images are then stacked together to form a 2D image. ln one embodiment the astigmatic part of the re-imaging system consists of one or more cylinder lenses. ln another embodiment the re-imaging system consists solely of mirrors, where a cylindrical or toroidal mirror provides the astigmatism. ln one preferred embodiment for 3D imaging applications an off axis cylindrical mirror is used as the astigmatic re-imaging optics in such a way as to keep the time delay between SLM and detector equal for all pixels on the SLM. ln one preferred embodiment the scene is illuminated by a pulsed laser and each pixel in the linear detector array comprises a temporally resolved detector circuit to provide 3D information about the scene through the time-of-flight laser radar 10 15 20 25 30 35 538 072 principle. ln one embodiment this temporally resolved detector circuit is a photodiode and a sampling circuit comprising a number of memory registers to provide a dense temporal sampling of the received radiation intensity. The Iinear architecture of the detector array allows dense packing of the detectors along the line at the same time as there is ample space for electronics for the sampling. ln another embodiment the detector array consists of a row of single photon avalanche diode (SPAD) detectors, each with separate electronics for collecting histograms of photon arrival times. This detector system comprlses a time-correlated single-photon counting (TCSPC) laser radar system. The Iinear detector array for a TCSPC-system may also consist of other photon counting detectors, e.g. superconducting nanowire single photon detectors. ln one embodiment the Iinear detector is the slit of a streak camera, allowing very high temporal resolution. ln one preferred embodiment the astigmatic re-imaging system also includes a dispersive element to re-image the NXP pixels of the SLM onto a Q> array, where each row of N pixels is redirected onto one row of Q pixels so that one wavelength component arrives at each of the Q pixels to produce a hyper-spectral imaging system. Every column of the QXP pixel array is then a sensor of the type described in the monochromatic implementations of this invention. The hyper- spectral sensor can be implemented either by placing the dispersive element in front of the focus of the astigmatic re-imaging system, or in the focus with a second re- imaging system directing the light to the detector array. ln one embodiment the dispersive element is a prism. ln another embodiment the dispersive element is a grating.

A simpler multispectral embodiment uses one or more chromatic beam splitters to direct the light to two or more discrete Iinear detector arrays. ln one embodiment the two mirror positions of the DMD reflect radiation into two different but identical astigmatic optical system and Iinear detector array systems, that by subtraction of the measurement data produce a random sampling matrix (<1>) consisting of values -1 and 1 instead of 0 and 1. This is used to improve numerical stability in the reconstruction process and hence reduce the number of measure- ments necessary, following the results of Sale et al.; please cf. Sale, D., Rozell, C.

J., Romberg, J. K., & Lanterman, A. D. (2012): Compressive ladar in realistic 10 15 20 25 30 35 538 072 environments. ln 2012 /EEE Statistical Signal Processing Workshop (pp. 720-723), hereby incorporated by reference. ln one preferred embodiment illustrated in Fig. 2 the patterns for compressed sensing processing are applied in the illumination source. A spatial light modulator projects a pattern of illuminated lines on the scene. A detector system comprising an astigmatic imaging system and a Iinear detector array is used so that the field of view of each detector is a stripe perpendicular to the illuminated lines on the target. ln one embodiment the illumination source is a pulsed laser to provide 3D information about the scene.

The illustration in Fig. 2 shows an imaging system where the light source 201 illumi- nates the whole field of view 202 in a pattern of vertical stripes 203. The light source includes a spatial light modulator to produce a changing set of vertical stripes. The spatial light modulator may be a DMD, and the full light source may be a standard computer projector. Light sources based on pulsed lasers, but otherwise similar to a projector, are suitable for longer ranges and 3D-imaging. The receiver subsystem consists of a Iinear detector array 211 and an astigmatic optical system 212. ln the simplest implementation the astigmatic optical system is a cylindrical lens. More complex systems consisting of multiple lens elements or cylindrical or toroidal mirrors to improve the light collection capacity of the detector subsystem are possible. A single pixel 213 of the Iinear detector array will have a horizontal slit like field of view 214 crossing the stripes produced by the light source. A different pixel 215 will have a similar field of view 216 at a different vertical position in the total field of view 202.

By performing a number of measurements with different patterns of vertical light stripes each detector element in the Iinear detector array will produce a set of collected data, which together with applied patterns of light stripes can be used to reconstruct the scene inside the horizontal slit seen by that detector element using compressive sensing reconstruction where the solution to a underdetermined Iinear equation system that maximizes the spasity of the scene is found. By adding these slit like scenes as lines in an image a two-dimensional image can be built.

A number of other concrete embodiments of the invention are possible and obvious within the inventive concept to the skilled man implementing the invention.

Claims (7)

Patentkrav:

1. Avbildande anordning innefattande en matrisdetektor (109,211) och en spatial ljusmodulator (103), vilket avbildningssystem upploser ett tvadimensionellt omrade (101,202) med anvandning av compressive sensing, k a n n et eck n ad a v att detektorn ar en linjar matrisdetektor som upplOser bildinformation tangs sin utstrackning genom att fokusera infallande straining pa detektorpixlar med hjalp av astigmatisk optik (108,202) och som upplOser bildinformation vinkelrat mot detektorns utstrackning med anvandning av ett flertal spatiala mOnster pa den spatiala ljusmodulatorn och compressive sensing-berakning, varigenom skapas ett antal compressive sensing-rekonstruktionsproblem som till antalet är lika manga som pixlarna i linjardetektorn och vardera av en matematisk dimension som motsvarar antalet element i den spatiala ljusmodulatorns mOnster vinkelratt mot matrisdetektorns utstrackning.

2. Avbildande system enligt patentkravet 1, kann ete ck n a t a v att namnda ljusmodulator (103) skapar ett monster av linjer (110,203) parallellt med riktningen for den linjara matrisdetektorn (109,211).

3. Avbildande system enligt patentkravet 1 eller 2, k ë n netecknat a v att ett system (201) som belyser scenen som ska avbildas innefattar den spatiala ljusmodulatorn.

4. Avbildande system enligt patentkravet 1 eller 2, k an n eteck nat a v att det innefattar standardoptik for avbildning som skapar en bild av scenen (101) som ska avbildas pa den spatiala ljusmodulatorn (103) och att av ljusmodulatorn transmitterat eller reflekterat ljus är ater-avbildat pa den linjara matrisdetektorn (109) av den astigmatiska optiken.

5. Avbildande system enligt patentkravet 4, k a n n et e ck n at a v att den spatiala ljusmodulatorn (103) är en digital mikrospegelanordning och att avbildningssystemet innefattar tva uppsattningar av en linjar matrisdetektor (109) och dess astigmatiska ater-avbildningsoptik (108), att ljus som reflekteras i tva riktningar fran den digitala mikrospeglingsanordningen uppsamlas av respektive linjara matris- detektor, och att det ena detektorresultatet fran matrisdetektorerna subtraheras fran det andra fOr att Oka den numeriska stabiliteten. 538 072

6. Avbildande system enligt nagot av de tidigare patentkraven, k a n net eck n at a v att den linjara matrisdetektorn (109,211) bestar av hyperspektrala detektorer implementerade som ett dispersivt element och en tva- dimensionell matrisdetektor.

7. Avbildande system enligt nagot av de tidigare patentkraven, k an neteck n at a v att det innefattar en pulsad ljuskalla (201) som belyser scenen som ska avbildas och att den linjara matrisdetektorn (211) bestar av tidsupplOsta detektorer for att skapa en 3D-bild av scenen. 538 072 M MIME Ia Mil • NIFAMEM II II k••••••In...wpme.gmempeopmtimi Iwo .■ 0 538 072 .4zr