Abstract:
An image cytometer for capturing and analyzing the image of a sample in the momentum domain. The cytometer is provided with a light source for illuminating a sample with a light beam, an optical transforming system positioned behind the sample in the beam propagation direction for generating the Fourier transform in the space plane, a light sensor array and a spatially selective filter positioned with respect to the optical system such that the Fourier transform is imaged onto the light sensor array.

Description:
FIELD OF THE INVENTION 
       [0001]    The present invention relates to optical devices, more particularly, an image cytometer. 
       BACKGROUND OF THE INVENTION 
       [0002]    Optical technologies can be very powerful when it comes to detecting minute quantities of particulate in biological and pathogenic samples. More specifically, optical microscopy is also widely used in combination with fluorescence labeling for high resolution imaging of particulates (cells, microorganisms, etc.) in large laboratories and clinics. Fluorescent microscopes on the market are bulky and expensive due to their cumbersome assembly and high cost of optical elements forming them. Over the past few years, efforts have been made to find more compact and inexpensive imaging solutions to meet market needs by making use of low cost imaging technologies based on charge-coupled devices (CCD) or complementary metal-oxide-semiconductor (CMOS) sensor arrays. 
         [0003]    Flow cytometry is a well-known laser based fluorescence technique that has experienced significant growth and innovation in recent years. This analytical laboratory technique can rapidly and very reliably measure different parameters on single cells and particles. Advances in the use of charge-coupled devices (CCD) have great potential in lowering prices by substituting expensive laser sources with much less expensive, up to two orders of magnitude, light-emitting diodes and sophisticated microscopes with simpler and more economic proximity detection schemes. These devices are known as image cytometers (I-CYTS) and can be easily operated through direct cell imaging on a computer screen. Unlike flow cytometers, I-CYTs do not work by illuminating the cells with a laser one by one but rather by imaging and analyzing thousands of cells in a single picture. For these reasons I-CYTs are gradually entering the market as they offer similar characteristics and benefits as conventional flow cytometers do but at a lower cost. U.S. Pat. No. 8,866,063 and U.S. Pat. No. 7,872,796 disclose microscope systems that use an image sensor, as is also the case for an image cytometer, that detect an image which is the replica of the sample in the space domain (known as the real or coordinate space). However, detecting the image in the space domain implies a tradeoff between spatial resolution, field of view (FOV) and depth of field (DOF). Resolution in optical microscope systems is limited by the diffraction limit (of the order of the wavelength of the light source, λ). In U.S. Pat. No. 7,872,796, a lenslet array is used to increase DOF of the device. However, FOV is still limited by the microscope objective in the system; therefore a λ/2 resolution can potentially be achieved with the device therein disclosed, but for reduced FOV. U.S. Pat. No. 8,866,063 describes a device with improved FOV by using an illumination source configured to scan the sample in at least two dimensions, capturing images in a plurality of scan locations. However, DOF is limited and the proposed solution to increase FOV implies lengthy computation to reconstruct the sample from the plurality of images captured; it also requires a cumbersome mechanical adaptation to provide the light source scanning. 
       BRIEF DESCRIPTION OF THE INVENTION 
       [0004]    The present invention captures and analyzes the image of the sample in the spatial frequency domain (also known as momentum domain or Fourier domain), in a plane that contains spatial frequency information. The real (space) image can then be reconstructed through Fourier analysis. The term “frequency” instead of “spatial frequency” will be used hereinafter. 
         [0005]    For this purpose, the present invention is provided with a light source for illuminating a sample with a light beam, an optical transforming system positioned behind the sample in the beam propagation direction for generating the Fourier transform in the space plane, a light sensor array and a spatially selective filter positioned with respect to the optical transforming system such that the Fourier transform is imaged onto the light sensor array. 
         [0006]    Given that sample information is retrieved through data in the frequency domain instead of the space domain, the invention is capable of analyzing large volumes (from one to ten milliliters) in a single capture as a result of a combination of a large FOV and DOF, automatically and accurately estimating the concentration of the particulate and differentiating population in terms of size and complexity (absorbance). Furthermore, with the disclosed system, spatial resolution can reach sub-micron resolution for a wide field of view choosing the proper light sensor array, doing so without the need for any mechanical adaptation or complex computation. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0007]    A set of drawings is attached to complete the description and to provide for better understanding of the invention. Said drawings illustrate a preferred embodiment of the invention, which should not be interpreted as restricting the scope of the invention, but merely as an example of how the invention can be embodied. 
           [0008]      FIG. 1  shows a general implementation of the invention. 
           [0009]      FIG. 2  shows one embodiment with an optical lens as the OTS. 
           [0010]      FIG. 3  shows a microlens array used in the invention as the spatially selective filter. 
           [0011]      FIG. 4  shows an absorptive polymer mask as the spatially selective filter. 
           [0012]      FIG. 5  shows one embodiment with a fluorescent filter. 
           [0013]      FIG. 6  shows another embodiment of a fluorescent cytometer according to the invention, with a reflective filter. 
           [0014]      FIG. 7  shows one embodiment where the source is placed at an angle in order to avoid the use of filters. 
           [0015]      FIG. 8  shows one embodiment with dual-band fluorescent bandpass filter. 
           [0016]      FIG. 9  schematically shows how data is processed. 
           [0017]      FIG. 10  shows a microlens array used in the invention. 
           [0018]      FIG. 11  shows an absorptive polymer mask used in the invention. 
           [0019]      FIG. 12  shows raw capture from the embodiment of  FIG. 3  proving the extended device FOV. 
           [0020]      FIG. 13  shows how the same data is retrieved from the sample regardless of distance to the OTS, thus proving the large device DOF capability (≈1 mm). 
           [0021]      FIG. 14  shows a graph acquired from processing sample volumes of two different microorganisms captured with the embodiment of  FIG. 3 . 
           [0022]      FIG. 15  shows a graph acquired from processing a water sample with a mixed population of microorganisms captured with the embodiment of  FIG. 3 . 
           [0023]      FIG. 16  shows a graph acquired from processing the water sample of  FIG. 14 , labeled with an anti-E. coil antibody conjugated to FITC and captured with the embodiment of  FIG. 5 . 
       
    
    
     DETAILED DESCRIPTION 
       [0024]    The invention works as follows. The sample is illuminated by a light source. For a given sample volume, interaction of said volume with the light source results in the emission of a broad band of frequencies. The emission process is incoherent; thus the intensities sum up and do not give rise to interference effects. The sample is placed in front of and in proximity to an optical transforming system (OTS). The OTS is an essential part, as it allows retrieving the Fourier transform of the beam generated by interaction between the sample volume and light beam. Examples of possible OTSs are a lens, a curved mirror with a radius of curvature equal to twice the focal length of a lens or a pinhole. 
         [0025]    Finally, a spatially selective filter is placed after the OTS at a distance suitable for achieving the Fourier transform of the beam interacting with the sample. The incoming beam is divided into regions in the XY plane through sampling with the sub-structures of the spatially selective filter, thus creating a set of sub-images; in other words, frequency information is selectively contained in the sub-images. Spatial sampling consists of the sub-images produced by the spatially selective filter and detected by the corresponding light sensors. The sensor array is located right after and in proximity to the spatially selective filter. The resulting image will be a composition of all the sub-images, each one containing a part of the Fourier transform of the beam after interaction with the sample volume. A measure on the resulting intensity distribution of these sub-images implies knowledge of the sample&#39;s energy spectrum. 
         [0026]    A preferred embodiment of the OTS is a converging lens. Such a lens inherently performs a two-dimensional Fourier transform, as will be explained. Consider the general geometry where a sample, located in front of the lens, is illuminated by a normally incident plane wave of amplitude A. Amplitude transmittance of the sample is represented by t A . In this case, the beam leaving the sample and incident on the lens can be described as: 
         [0000]        U   L ( x,y )= A·t   A ( x,y ) 
         [0027]    The effects of the lens on the incident optical beam due to the optic can be described by: 
         [0000]        t   I ( x,y )=exp[− j ·( k/ 2· f )·( x   2   +y   2 )]
 
         [0028]    Amplitude distribution behind the lens therefore becomes 
         [0000]        U   L ′( x,y )= U   L ( x,y )·exp[− j ·( k/ 2· f )˜( x   2   +y   2 )]
 
         [0029]    Using the Fresnel diffraction formula, one can find the distribution U f (u,v) at a distance z from the lens. 
         [0000]        U   f ( u,v )={exp[ j ·( k /(2· z ))·( u   2   +v   2 )]·(1/( j·λ·z ))·∫∫ U   L ( x,y )·exp[− j ·( k/ 2· f )·( x   2   +y   2 )]·exp[ j ·( k /(2· z ))·( x   2   +y   2 )]·exp[− j ·(2π/(λ· z ))·( x·u+y·v )] dxdy  
 
         [0030]    Field distribution U f (u,v) is therefore proportional to the two-dimensional Fourier transform of the incident field subtended by lens aperture with a quadratic phase factor, and the amplitude and phase of the light at coordinates (u, v) are related to the amplitude and phase of the input spectrum at frequencies (u/λf), v/(λf)) 
         [0031]    Let the converging lens have a focal length “f”. The sample volume is assumed to be uniformly illuminated by a normally incident beam from the light source. After the lens, complex field distribution is proportional to the two-dimensional Fourier transform within the lens aperture. The amplitude and phase of the light field that has traversed the sample volume in the spatial domain can then be reconstructed from the amplitude and phase of the measured Fourier components at frequencies (f′=u/(λ·f)). Therefore, it is possible to retrieve sample information from the measured power (energy) distribution after the OTS. Consider the following more specific geometry: the input sample volume, located at a distance “d”, where “d” is much smaller than “f”, in front of (in proximity to) the converging lens, is illuminated by a normally incident collimated beam. The Fourier transform of the beam after the sample volume is imaged at a distance “f” after the converging lens. 
         [0032]    For a pinhole, transformation from the space domain into the frequency domain works as follows: 
         [0033]    A plane wave of light incident upon an aperture [A(x)] will produce the Fourier transform of A(x) on the image plane. In particular, consider a diffracting aperture that is circular rather than rectangular, and let the radius of the aperture be w. Therefore, if q is a radius coordinate in the plane of the aperture, then t A (q)=circ(q/w). The circular symmetry of the problem suggests that the Fourier transform be rewritten as a Fourier-Bessel transform. Amplitude distribution in the Fraunhofer diffraction pattern is 
         [0000]        U ( r )= A ·exp[ j·k·z ]exp[ j·k·r ̂2/(2 z )]·(1/( j ·λ))·[2· J   1 ( k·w·r/z )/( k·w·r/z )],
 
         [0000]    where J 1  is the Bessel function of the first kind. 
         [0034]    A curved mirror is equivalent to a lens when the radius of curvature of the mirror is equal to twice the focal length (f) of the lens. 
         [0035]    Once the optical signal has been transformed by means of the OTS, an image is selected by means of a spatial filter. Examples of such a filter are apertures, a microlens array and an absorptive polymer mask, i.e., a two-dimensional array-like structure built of several apertures capable of spatially filtering an incoming beam. Each aperture selects a portion of the incoming beam preserving information about both its amplitude and phase. Sample information contained in the optical signal detected after the spatially selective filter can thereby be retrieved after suitable data processing of the image detected by the light sensor array. Multiple two-dimensional array-like structures are known to comply with the spatially selective filtering requirements to implement the disclosed invention; among them, a microlens array and an absorptive polymer mask made up of a transmitting well array are considered preferred solutions for this invention. The pitch of the lens array, polymer mask or any other structure complying with the specifications, will define the periodicity of the sub-images captured by the light sensor array. The aperture of the microlenses, wells or equivalent will define spatial filtering capabilities of the spatially selective filter. 
         [0036]    In the case of a microlens array configuration, filtering capabilities are defined by the angular aperture of each microlens (N m ): 
         [0000]    
       
         
           
             
               N 
               m 
             
             = 
             
               
                 f 
                 m 
               
               
                 D 
                 m 
               
             
           
         
       
     
         [0037]    In the equation, “f m ” represents the focal length of each microlens and “D m ” its diameter. 
         [0038]    A microlens is a small lens, generally with a diameter less than one millimeter (mm) and often as small as 10 micrometers (μm). A typical microlens may be a single element with one plane surface and one spherical convex surface to refract the light. Since microlenses are so small, the substrate that supports them is usually thicker than the lens. More sophisticated lenses may use aspheric surfaces, and others may use several layers of optical material to achieve the designed performance. Microlens arrays contain multiple microlenses formed in a one-dimensional or two-dimensional array on a supporting substrate. If the individual lenses have circular apertures and are not allowed to overlap, they may be placed in a hexagonal array to obtain maximum coverage of the substrate. However, there will still be gaps between lenses which can only be reduced by making microlenses with non-circular apertures. With optical sensor arrays, very small lens systems serve to focus and concentrate the light onto the photodiode surface instead of allowing it to fall on non-photosensitive areas of the pixel device. Fill-factor is the ratio of the active refracting area, i.e., that area which directs light to the detecting surface, to the total contiguous area occupied by the microlens array. The microlens array works as an array of point scanning microscope objectives over the optical Fourier transform of the volumetric sample. 
         [0039]    In the case of an absorptive polymer mask, the filtering effect can be described as follows: 
         [0000]    
       
         
           
             
               tan 
                
               
                   
               
                
               
                 ϕ 
                 m 
               
             
             = 
             
               
                 d 
                 m 
               
               
                 L 
                 m 
               
             
           
         
       
     
         [0040]    In the equation, “d m ” represents the diameter of the mask aperture element, “L m ” the height of the mask and the angle phi, “(φ m ”, represents the acceptance angle with respect to the normal which is inherently linked to the spatial frequencies of the incoming beam. An absorptive polymer mask is a well array with a pitch equivalent to that of the microlens array but apertures (1% to 30%) smaller than the diameter of the lenses, thus providing a robust structure. The height of the mask defines filtering capabilities. Said height should be within a range that avoids both aliasing and under-sampling; it should preferably result in an aspect ratio between the height of the mask and the well aperture (AR=L m /d m ) from 1 to 10, which is equivalent to an acceptance angle from 6 to 45 degrees. An acceptance angle below 6 degrees may result in under-sampling the optical signal and therefore not enough information for sample recovery, while an acceptance angle higher than 45 degrees may result in serious detrimental aliasing. The thicker the mask (increase in height), the more selective the structure. However, there is a tradeoff between high selectivity and under-sampling of the signal which would result in loss of information. 
         [0041]    Finally, the resulting image is detected by means of a light sensor array. Examples of such an array are CODs, cameras in conventional mobile phones or any other portable devices, etc. 
         [0042]    Data extraction is based on Fourier optics principles; the detected pattern intensity is given by the following equation: 
         [0000]    
       
         
           
             
               I 
               
                 ( 
                 
                   x 
                   , 
                   y 
                 
                 ) 
               
             
             ≈ 
             
               
                 
                   
                     [ 
                     
                       A 
                       
                         2 
                          
                         
                             
                         
                          
                         λ 
                          
                         
                             
                         
                          
                         
                           Z 
                           i 
                         
                       
                     
                     ] 
                   
                   2 
                 
                 · 
                 
                   
                     sync 
                     2 
                   
                    
                   
                     ( 
                     
                       
                         2 
                          
                         wy 
                       
                       
                         λ 
                          
                         
                             
                         
                          
                         
                           Z 
                           i 
                         
                       
                     
                     ) 
                   
                 
               
                
               
                 { 
                 
                   
                     
                       sync 
                       2 
                     
                      
                     
                       ( 
                       
                         
                           2 
                            
                           wx 
                         
                         
                           λ 
                            
                           
                               
                           
                            
                           
                             Z 
                             i 
                           
                         
                       
                       ) 
                     
                   
                   + 
                   
                     
                       
                         m 
                         2 
                       
                       4 
                     
                      
                     
                       
                         sync 
                         2 
                       
                        
                       
                         [ 
                         
                           
                             
                               2 
                                
                               w 
                             
                             
                               λ 
                                
                               
                                   
                               
                                
                               
                                 Z 
                                 i 
                               
                             
                           
                            
                           
                             ( 
                             
                               x 
                               + 
                               
                                 
                                   f 
                                   0 
                                 
                                  
                                 λ 
                                  
                                 
                                     
                                 
                                  
                                 
                                   Z 
                                   i 
                                 
                               
                             
                             ) 
                           
                         
                         ] 
                       
                     
                   
                   + 
                   
                     
                       
                         m 
                         2 
                       
                       4 
                     
                      
                     
                       
                         sync 
                         2 
                       
                        
                       
                         [ 
                         
                           
                             
                               2 
                                
                               w 
                             
                             
                               λ 
                                
                               
                                   
                               
                                
                               
                                 Z 
                                 i 
                               
                             
                           
                            
                           
                             ( 
                             
                               x 
                               - 
                               
                                 
                                   f 
                                   0 
                                 
                                  
                                 λ 
                                  
                                 
                                     
                                 
                                  
                                 
                                   Z 
                                   i 
                                 
                               
                             
                             ) 
                           
                         
                         ] 
                       
                     
                   
                 
                 } 
               
             
           
         
       
     
         [0043]    In the case of two-dimensional geometry, this represents spatial frequency distribution. The center of the image is zero spatial frequency and each of the other intensity points (nλZ i f 0 ) represents a harmonic of the captured signal. 
         [0044]    Each of the sub-images contains information about statistical parameters of the sample volume, such as particulate size and complexity, also allowing counting of particles. Complexity and size information can be paired into a dispersion graph to present the results to the user. Said dispersion graphs will allow differentiating multiple particulates within a single sample volume. The complexity parameter refers to a measure of the absorbance of the sample, i.e., how much of the incident beam is absorbed by the sample. 
         [0045]    Different embodiments are shown below in reference to the drawings. 
         [0046]      FIG. 1  shows a general schematic of the invention, comprising: a wavelength band-limited light source (A 1 ), an optical transforming system (OTS) (B 1 ), a spatially selective filter (B 2 ), and a light sensor array (A 2 ). The figure also indicates the position of the sample (C 1 ) within the apparatus. The sample is to be placed between the light source and the OTS. 
         [0047]    In  FIG. 2 , an optical lens is used as the OTS (B 11 ). The spatially selective filter (B 2 ) is placed at the lens&#39;s angular frequency plane and the sample volume (C 1 ) is placed in proximity to the OTS lens. 
         [0048]      FIG. 3  shows an embodiment using a microlens array as a spatially selective filter and an optical lens to focus the light source onto the sample for improved illumination. The focusing optical lens (e.g. focal length of 50 mm) (B 3 ) focuses the wavelength band-limited LED light source (A 1 ) onto an area of the sample volume; the resulting beam is collected by the OTS (B 1 ). The beam leaving the OTS is incident on the microlens array (B 21 ), which decomposes the spatial Fourier transform components generated by the OTS (B 1 ). Said array can have an area of 10 mm×10 mm array with a pitch of 300 μm for a total of 1111 microlenses and with a thickness of 1.2 mm. The light sensor array (A 2 ) can be a commercially available CMOS image sensor with Bayer Filter technology, a size of 5.70 mm×4.28 mm and a pitch of 2.2 μm. 
         [0049]      FIG. 4  shows another embodiment of the invention that uses an absorptive polymer mask instead of a microlens array in order to prevent or at least reduce aliasing effects in the detected signal. The rest of the components in the apparatus are the same as in the embodiment of  FIG. 2 . The absorptive polymer mask can be a structure measuring 10 mm×10 mm×6 mm (XYZ) made of a two-dimensional aperture array 250 μm in diameter and 300 μm in pitch. The resulting overall structure is made of 1111 different apertures each of which independently samples the incoming beam onto the light sensor array. 
         [0050]    One embodiment using a fluorescent filter is shown in  FIG. 5 . An appropriate fluorescent filter attenuates (absorbs) the pump light and lets the fluorescent signal be transmitted. This introduces an additional degree of specificity into the system, especially when size and/or complexity are not enough to differentiate between particulates. The system as described in  FIG. 1  has been proven to differentiate particulates by size down to a 3 um difference; however, many biological samples to analyze might have several microorganisms of equal size. In this scenario, labeling and adding the fluorescent filter of  FIG. 5  allows distinguishing the target microorganism. The focusing optical lens (B 3 ) (e.g. focal length of 50 mm) focuses the wavelength band-limited LED light source (A 1 ) onto an area of the sample volume (C 2 ). This sample volume generates a fluorescent beam (C 3 ) which is collected by an OTS (B 1 ). The remaining pump signal is attenuated (filtered) by the fluorescent filter (B 4 ), thereby leaving only fluorescent emission going through the spatially selective filter (B 2 ) and detected by the light sensor array (A 2 ). 
         [0051]      FIG. 6  shows another embodiment using a reflective filter (B 5 ) to suppress the pump beam remaining after the sample volume. Rejection of the pump is especially important in the case of fluorescence detection. If the pump is rejected more than the fluorescent signal, a fluorescent filter might not be needed. Fluorescent emission is omnidirectional; therefore, even though some fluorescent signals will be reflected by the reflective filter, enough fluorescent intensity will reach the light sensor array (A 2 ) so as to retrieve sample information from the capture. 
         [0052]    In the embodiment of  FIG. 7 , the incident pump beam enters the sample volume at an angle. Use of filters can thereby be avoided. 
         [0053]      FIG. 8  shows a multi-fluorescent (dual-fluorescent) I-CYT configuration. The light source (A 1 ) represents one or two wavelength band-limited LED light sources. A dual-band fluorescent bandpass filter (B 6 ) is used, replacing the fluorescent filter in  FIG. 5 . More specifically, this system could be applicable to analyze a sample volume fluorescently marked with FITC (excitation and emission wavelengths of 495 nm and 519 nm, respectively) and PerCP (excitation and emission wavelengths of 470 nm and 670 nm, respectively) using a single LED light source centered within the excitation peaks and a 524 nm/676 nm dual-band fluorescent bandpass filter. Information can be singularized at the sensor by using a color image sensor and retrieving data from the corresponding color channel for each emitted wavelength. This introduces further specificity into the system, since particulates can now be sorted by complexity, size and fluorescence emission. 
         [0054]    What has been said for fluorescence signals can be applied to autofluorescence geometries. This means in this case that the pump induces fluorescence directly from the particulate to detect without the need for fluorescent labeling. This is because the fluorescence from labels is attached to particulates and autofluorescence from the particulates themselves comes from the same sample volume regions, so the signal can be processed in the same way and the invention can thus be applied. 
         [0055]    The operability of the invention allows to use standard optical chambers, cuvettes and fluidics (devices with no moving parts), to allocate the volume sample. An optical chamber or cuvette is a small tube of circular or square cross section, sealed at one end, made of plastic, glass, or fused quartz (for UV light) and designed to hold samples for experiments. Disposable plastic cuvettes are often used in fast assays, where speed is more important than high accuracy. 
         [0056]    Furthermore, the raw sample can be filtered and purified while being concentrated using standard concentrator filters. Particulate concentrators are disposable, single use only ultrafiltration devices with polymer membranes for the concentration and/or purification of biological samples. The optical signal response of the invention to the number of particulates is linear, hence the invention requires a single value correction factor when using concentrator devices to better measure the real particulate concentration. 
         [0057]      FIG. 9  shows optical data processing. If the sample volume (C 1 ) is placed in proximity to the OTS (B 1 ), the principles of Fourier optics establish that at a distance (B 8 ) from (B 1 ), the position known as the spatial frequency plane or Fourier plane, the received signal (C 4 ) will be equivalent to the Fourier transform of the input beam collected by the OTS (B 1 ). By placing a spatially selective filter (B 2 ) at said point (B 8 ), several sub-images of the spatial Fourier transform will be generated, each of them corresponding to one sub-structure of the filter. A complete image made up of the aforementioned sub-images is captured at the sensor (A 2 ). 
         [0058]    A schematic of a microlens array (B 21 ) used in one embodiment of the invention can be seen in  FIG. 10 . (B 211 ) refers to the full size of the array in one direction, (B 212 ) to the size in the orthogonal direction, (B 213 ) indicates the thickness of the array and (B 214 ) is the pitch between lenses. For the preferred embodiment of  FIG. 2 , the microlens array has a length of 10 mm along (B 211 ), 10 mm along ( 212 ), 1.2 mm of thickness (B 213 ) and a pitch (B 214 ) of 300 μm, resulting in an array with a total of 1111 microlenses. 
         [0059]    A schematic of the absorptive polymer mask (B 22 ) used in one embodiment of the invention can be seen in  FIG. 11 . (B 221 ) refers to the full size of the array in directions X and Y, (B 222 ) refers to the size in the orthogonal Z direction, (B 223 ) indicates the pitch between the apertures of the mask and (B 224 ) the diameter of said apertures. For the preferred embodiment of  FIG. 3 , the absorptive polymer mask has a length of 10 mm along (B 221 ), 6 mm along (B 222 ), a pitch (B 223 ) of 300 μm and an aperture size (B 224 ) of 250 μm, resulting in a structure with a total of 1111 apertures. 
         [0060]      FIG. 12  shows a raw capture from the embodiment of  FIG. 3 , proving the extended device FOV. In the case of prior art for imaging applications, including light-field microscopy, FOV is dependent on the microscope objective used in the device. For example, for a 16× objective lens FOV is 1.3 mm 2 . In the current invention, the sample is placed behind and very close to the lens capturing the Fourier transform. Hence, the magnification factor in the system is approximately 1. Therefore, the area of sample capture is equal to that imaged at the sensor, i.e. 5.70 mm×4.28 mm=24.39 mm 2 , achieving a large FOV=24.39 mm 2 . In contrast, in the prior art capturing images in real space means that the sample must be placed apart from the lens, and therefore a magnification factor greater than 1 reduces the effective field-of-view. 
         [0061]    In  FIG. 13  a sample of 5 μm particles was captured at several distances from the OTS and processed to retrieve sample size and complexity information. The resulting graph shown indicates how same data is retrieved from the sample regardless of its distance to the OTS, thus proving the large device DOF capability ≈1 mm). Detection of the particulate is repetitive in terms of the retrieved size and complexity information. 
         [0062]      FIG. 14  shows a dispersion graph acquired from processing two sample volumes captured with the embodiment of  FIG. 3 , in which the OTS (B 1 ) is an optical lens (B 11 ) as shown in  FIG. 2 . A sample of  Escherichia coli  ( E. coli ) at a concentration of 10 4  CFU/ml and a sample of  Saccharomyces cerevisiae  ( S. cerevisiae ) having a similar concentration, both diluted in phosphate buffered saline (PBS), were captured and processed.  E. coli  microorganisms are known to a have an average size of 2 μm, whereas  S. cerevisiae  have an average size of 5 μm. The graph displays complexity (horizontal axis) and size (vertical axis) of the particulate within the analyzed volumetric samples. From the graph it is evident that different microorganisms can be detected, characterized and differentiated in terms of their size and complexity parameters retrieved from the captured signal 
         [0063]      FIG. 15  depicts a dispersion graph acquired from processing a contaminated water sample captured with the embodiment of  FIG. 3 . The sample initially had an unknown variety of microorganisms. To ensure the presence of  E. coli,  1 ml of a 10 4  CFU/ml  E. coli  concentration, diluted in PBS with anti- E. coli  antibody conjugated to fluorescein isothiocyanate (FITC) fluorescent marker was added to 2 ml of a contaminated water sample. This process ensures fluorescently labeling the added 10 4  CFU/ml concentration of  E. coli . The resulting graph shows information corresponding to the complete microbiological load of the sample. 
         [0064]      FIG. 16  shows the dispersion graph corresponding to the same sample, this time captured with the embodiment of  FIG. 5 . In this case, only data corresponding to the FITC emitting  E. coli  microorganisms is retrieved. 
         [0065]    In this text, the term “comprises” and its derivatives (such as “comprising”, etc.) should not be understood in an excluding sense, that is, these terms should not be interpreted as excluding the possibility that what is described and defined may include further elements, steps, etc. 
         [0066]    On the other hand, the invention is obviously not limited to the specific embodiment(s) described herein, but also encompasses any variations that may be considered by any person skilled in the art (for example, as regards the choice of materials, dimensions, components, configuration, etc.), within the general scope of the invention as defined in the claims.