# Patent

## Title
Systems, methods, and media for recording an image using an optical diffuser

## Abstract

Systems, methods, and media for recording an image of a scene ate provided. In accordance with some embodiments, systems for recording an image of a scene are provided. comprising. a diffuser {hat diffuses light representing the scene and that has a scattering function that is independent of aperture coordinates, a sensor that receives diffused light representing the scene and generates data representing an image; and a hardware processes that uses a point spread function to deblur the image.

## Background
CROSS REFERENCE TO RELATED APPLICATIONS
This application claims the benefit of U.S. Provisional Patent Application No. 61/297,667, filed Jan. 22, 2010, which is hereby incorporated by reference herein in its entirety.


STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT
This invention was made with government support under Grant Numbers N00014-08-1-0329 and N00014-08-1-0638 awarded by the Navy/Office of Naval Research. The government has certain rights to the invention.


TECHNICAL FIELD
The disclosed subject matter relates to systems, methods, and media for recording an image using an optical diffuser.
BACKGROUND
For conventional cameras, there is a fundamental trade-off between depth of field (DOF) and noise. Generally, cameras have a single focal plane, and objects that deviate from this plane are blurred due to defocus. The amount of defocus blur depends on the aperture size and the distance from the focal plane. To decrease defocus blur and increase DOF, the aperture size must be decreased, reducing the signal strength of the recorded image as well. In many cases, it is desirable to have a DOF that is as large as possible so that all details in the scene are preserved. This is the case, for instance, in machine vision applications such as object detection and recognition, where it is desirable that all objects of interest be in focus. However, stopping down the lens aperture is not always an option, especially in low light conditions, because it can increase noise, which in turn can materially impact the recorded image.
SUMMARY
Systems, methods, and media for recording an image of a scene are provided. In accordance with some embodiments, systems for recording an image of a scene are provided, comprising: a diffuser that diffuses light representing the scene and that has a scattering function that is independent of aperture coordinates; a sensor that receives diffused light representing the scene and generates data representing an image; and a hardware processor that uses a point spread function to deblur the image.
In accordance with some embodiments, methods for recording an image of a scene are provided, the methods comprising: diffusing light representing the scene using a diffuser that has a scattering function that is independent of aperture coordinates; receiving diffused light representing the scene and generating data representing an image; and using a point spread function to deblur the image



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a diagram of a mechanism for recording an image in accordance with some embodiments.
FIG. 2 is a combination of two images, one not using a diffuser (a) and one using a diffuser (b), in accordance with some embodiments.
FIG. 3 is a diagram of a lens and a sensor in accordance with some embodiments.
FIG. 4 is a diagram of a lens, a diffuser, and a sensor in accordance with some embodiments.
FIG. 5 is a diagram illustrating a light field on a sensor in accordance with some embodiments.
FIG. 6 is a diagram of a ray and the scattering of the ray in accordance with some embodiments.
FIG. 7 is a combination of point-spread-function and modulation-transfer-function graph pairs in accordance with some embodiments.
FIG. 8 is a diagram an optical system including a wedge (a) and a random varying surface in accordance with some embodiments.
FIG. 9 is combination of diagrams of a diffuser profile (a), a diffuser height map (b), a diffuser scatter PDF (c), and a diffuser (d) in accordance with some embodiments.



DETAILED DESCRIPTION
Systems, methods, and media for recording an image using an optical diffuser are provided.
Turning to FIG. 1, an illustration of an image recording mechanism 102 (e.g., a camera, video camera, mobile phone incorporating a camera, and/or any other suitable image recording mechanism) that is being used to capture an image including three objects, A 104, B 106, and C 108, is shown. As can be seen, these objects are at different depths with respect to mechanism 102. Because of limitations in the depth of field of mechanism 102, objects A 104 and C 108 may be out of focus when mechanism 102 is focused on object B 106. For example, these objects may be the toys shown in FIG. 2. As illustrated in FIG. 2(a), when a camera is focused on the center object (which may correspond to object B 106 of FIG. 1), the other objects may be out of focus. By using the mechanisms as described herein, however, an image can be recorded of such objects so that they appear to be in focus as illustrated in FIG. 2(b). This can be referred to as mechanism 102 having an extended depth of field.
In accordance with some embodiments, extended depth of field can be achieved by incorporating a diffuser 110 or 112 into an image recording mechanism 102. Recording an image using a diffuser in the pupil plane of an image recording mechanism can be referred to as diffusion coding. Such a diffuser can be located at any suitable point in the image recording mechanism. For example, a diffuser 110 can be positioned between a light source (e.g., objects 104, 106, and 108) and a lens 114 (e.g., as a lens attachment), a diffuser 112 can be positioned between a lens 114 and a sensor 116 (e.g., as part of a lens or a camera body), etc.
The diffusion coding image can then be detected by sensor 116 and then provided to a hardware processor 118 (incorporated into mechanism 102) and/or a hardware processor 120 (external to mechanism 102) for subsequent processing. This processing can include deblurring the sensed image using a PSF that is matched to the PSF of the diffuser. Any other suitable processing can additionally or alternatively be used. After such processing, an extended depth of field image can be presented on a display 124 (internal to mechanism 102) and/or a display 122 (external to mechanism 102).
In order to illustrate how such an image can be recorded using a diffuser, the optics of some embodiments are now described.
As shown in FIG. 3, a light field L(ū,x) can be used to represent a four-dimensional set of rays propagating from an ideal lens with effective focal length (EFL) f to a sensor. A vector ū=(u,v) can be used to denote the coordinates on the u-v plane, which is coincident with the exit pupil of the lens, and a vector x=(x,y) can be used to denote the coordinates on the x-y plane that is coincident with the sensor. The irradiance E(x) observed on the sensor can be defined as the light field integrated over all ray angles:




  
    
      
        
          
            E
            ⁡
            
              (
              
                x
                _
              
              )
            
          
          =
          
            
              ∫
              
                Ω
                
                  u
                  _
                
              
            
            ⁢
            
              
                L
                ⁡
                
                  (
                  
                    
                      u
                      _
                    
                    ,
                    
                      x
                      _
                    
                  
                  )
                
              
              ⁢
              
                d
                
                  u
                  _
                
              
            
          
        
        ,
      
    
    
      
        (
        1
        )
      
    
  




where, Ωū is the domain of ū. For a scene with smooth depth variation, locally, the captured image E(x) can be modeled as a convolution between a depth-dependent PSF kernel P(x) and an all-in-focus image I(x).

As described further below, in accordance with some embodiments, a camera PSF can be shaped so that an image I(x) can be recovered from the captured image E(x) by deblurring with a single PSF P(x). The depth-dependence of the camera PSF can be analyzed by considering the image produced by a unit energy point source. For example, as shown in FIG. 3, consider a point source whose image comes to focus at a distance d0 from the aperture of the lens. Assuming a rectangular aperture of width A, the light field produced by this point can be represented as:




  
    
      
        
          
            
              L
              δ
            
            ⁡
            
              (
              
                
                  u
                  _
                
                ,
                
                  x
                  _
                
              
              )
            
          
          =
          
            
              1
              
                A
                2
              
            
            ⊓
            
              
                (
                
                  
                    u
                    _
                  
                  A
                
                )
              
              ⁢
              
                δ
                ⁡
                
                  (
                  
                    
                      x
                      _
                    
                    -
                    
                      
                        s
                        0
                      
                      ⁢
                      
                        u
                        _
                      
                    
                  
                  )
                
              
            
          
        
        ,
      
    
    
      
        (
        2
        )
      
    
  




where s0=(d0−f)/d is the defocus slope in light field space, and  is the box function:





  
    
      
        
          ⊓
          
            (
            
              
                x
                _
              
              w
            
            )
          
        
        =
        
          {
          
            
              
                
                  1
                
                
                  
                    
                      
                        if
                        ⁢
                        
                          
                        
                        ⁢
                        
                          
                          
                            
                              x
                              _
                            
                            i
                          
                          
                        
                      
                      <
                      
                        1
                        2
                      
                    
                    ,
                    
                      ∀
                      i
                    
                  
                
              
              
                
                  0
                
                
                  otherwise
                
              
            
            .
          
        
      
    
    
      
        (
        3
        )
      
    
  




The image of this point is the camera PSF at the depth d0, which is a box shaped PSF with defocus blur width s0A:





  
    
      
        
          P
          ⁡
          
            (
            x
            )
          
        
        =
        
          
            1
            
              
                s
                0
                2
              
              ⁢
              
                A
                2
              
            
          
          ⊓
          
            
              (
              
                
                  x
                  _
                
                
                  
                    s
                    0
                  
                  ⁢
                  A
                
              
              )
            
            .
          
        
      
    
    
      
        (
        4
        )
      
    
  




The effect of a general kernel D applied to a light field L, which represents the effect of a diffuser placed in the aperture of a camera lens, can next be analyzed. The kernel can produce a new filtered light field {circumflex over (L)}, from which the modified PSF {circumflex over (P)} can be derived as:




  
    
      
        
          
            
              L
              ^
            
            ⁡
            
              (
              
                
                  u
                  _
                
                ,
                
                  x
                  _
                
              
              )
            
          
          =
          
            
              ∫
              
                Ω
                
                  
                    u
                    _
                  
                  ′
                
              
            
            ⁢
            
              
                ∫
                
                  Ω
                  
                    
                      x
                      _
                    
                    ′
                  
                
              
              ⁢
              
                
                  D
                  ⁡
                  
                    (
                    
                      
                        u
                        _
                      
                      ,
                      
                        
                          u
                          _
                        
                        ′
                      
                      ,
                      
                        x
                        _
                      
                      ,
                      
                        
                          x
                          _
                        
                        ′
                      
                    
                    )
                  
                
                ⁢
                
                  L
                  ⁡
                  
                    (
                    
                      
                        
                          u
                          _
                        
                        ′
                      
                      ,
                      
                        
                          x
                          _
                        
                        ′
                      
                    
                    )
                  
                
                ⁢
                
                  d
                  
                    
                      u
                      _
                    
                    ′
                  
                
                ⁢
                
                  d
                  
                    
                      x
                      _
                    
                    ′
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        5
        )
      
    
  
  
    
      
        
          
            
              P
              ^
            
            ⁡
            
              (
              
                x
                _
              
              )
            
          
          =
          
            
              ∫
              
                Ω
                
                  u
                  _
                
              
            
            ⁢
            
              
                
                  L
                  ^
                
                ⁡
                
                  (
                  
                    
                      u
                      _
                    
                    ,
                    
                      x
                      _
                    
                  
                  )
                
              
              ⁢
              
                d
                
                  u
                  _
                
              
            
          
        
        ,
      
    
    
      
        (
        6
        )
      
    
  




where Ωx is the domain of x. This approach allows a large class of operations applied to a light field to be expressed. For instance, consider a kernel of the form





  
    
      
        
          D
          ⁡
          
            (
            
              
                u
                _
              
              ,
              
                
                  u
                  _
                
                ′
              
              ,
              
                x
                _
              
              ,
              
                
                  x
                  _
                
                ′
              
            
            )
          
        
        =
        
          
            
              1
              
                w
                2
              
            
            ⁢
            
              δ
              ⁡
              
                (
                
                  
                    u
                    _
                  
                  -
                  
                    
                      u
                      _
                    
                    ′
                  
                
                )
              
            
          
          ⊓
          
            
              (
              
                
                  
                    x
                    _
                  
                  -
                  
                    
                      x
                      _
                    
                    ′
                  
                
                w
              
              )
            
            .
          
        
      
    
    
      
        (
        7
        )
      
    
  




Note that here D takes the form of a separable convolution kernel with finite support in the x domain. The geometric meaning of this kernel can be illustrated as shown in FIG. 4. As shown, each ray in the light field is blurred so that, instead of piercing the sensor at a single location, it contributes to a square of width w. In order to understand the effect of the diffuser, an image E captured without the diffuser can be compared to an image Ê captured with it. For this diffuser kernel, substituting Equation 7 into Equations 5 and 6 gives:




  
    
      
        
          
            
              P
              ^
            
            ⁡
            
              (
              
                x
                _
              
              )
            
          
          =
          
            
              1
              
                w
                2
              
            
            ⊓
            
              
                (
                
                  
                    x
                    _
                  
                  w
                
                )
              
              ⊗
              
                P
                ⁡
                
                  (
                  
                    x
                    _
                  
                  )
                
              
            
          
        
        ,
      
    
    
      
        (
        8
        )
      
    
  




where {circle around (x)} denotes convolution. The modified PSF can be the camera PSF blurred with a box function. Therefore, the effect of the diffuser is to blur the image that would be captured were it not present. However, the diffuser given by the kernel in Equation 7 may not be useful for extending depth of field because it does not increase depth independence or preserve high frequencies in the camera PSF.

In general, the kernel for any diffuser that is placed in the aperture can be represented as:

D(ū,ū′,x,x′)=δ(ū−ū′)k(ū,x−x′),  (9)

where k is called the scatter function. As can be seen, the diffuser has no effect in the ū domain, but has the effect of a convolution in the x domain. For the diffuser given by Equation 7, the scatter function can be represented as a two-dimensional box function:





  
    k
    ⁡
    
      (
      
        
          u
          _
        
        ,
        
          x
          _
        
      
      )
    
  
  =
  
    
      1
      
        w
        2
      
    
    ⊓
    
      
        (
        
          
            x
            _
          
          w
        
        )
      
      .
    
  




By changing from rectangular coordinates (u, v, x, y) to polar coordinates (ρ, φ, r, θ) using the relations u=ρ cos φ, v=ρ sin φ, x=r cos θ, and y=r sin θ, a polar system where ρ,r ε (−∞,∞) and θ,φε(0,π) and a circular aperture with diameter A can be considered. In this system, the light field representing a unit-energy point source located at distance d0 can be written as:




  
    
      
        
          
            
              L
              δ
            
            ⁡
            
              (
              
                ρ
                ,
                r
              
              )
            
          
          =
          
            
              4
              
                π
                ⁢
                
                  
                
                ⁢
                
                  A
                  2
                
              
            
            ⊓
            
              
                (
                
                  ρ
                  A
                
                )
              
              ⁢
              
                
                  δ
                  ⁡
                  
                    (
                    
                      r
                      -
                      
                        
                          s
                          0
                        
                        ⁢
                        ρ
                      
                    
                    )
                  
                
                
                  π
                  ⁢
                  
                    
                    r
                    
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        10
        )
      
    
  




which is independent of both θ and φ because the source is isotropic. Note that verifying unit-energy can be carried out trivially by integrating Lδ(ρ,r) in polar coordinates. Comparing the parameterizations for the light field of a point source in Equations 2 and 10, it can be seen than a slice of Lδ(x,y) represents a single ray, while a slice L(ρ,r) represents a 2D set of rays. In the radially symmetric parameterization, a slice of the light field represents a conic surface connecting a circle with radius ρ in the aperture plane to a circle of radius r on the sensor (see FIG. 5).

A radially symmetric diffuser produces a drastically different effect than the diffuser given by Equation 7. When a radially symmetric diffuser is introduced, neither the diffuser nor the lens deflects rays tangentially, and therefore the diffuser kernel and modified light field can be represented using the reduced coordinates (ρ,r). Equations 5 and 6 then become:




  
    
      
        
          
            
              L
              ^
            
            ⁡
            
              (
              
                ρ
                ,
                r
              
              )
            
          
          =
          
            
              π
              2
            
            ⁢
            
              
                ∫
                
                  Ω
                  ρ
                
              
              ⁢
              
                
                  ∫
                  
                    Ω
                    r
                  
                
                ⁢
                
                  
                    D
                    ⁡
                    
                      (
                      
                        ρ
                        ,
                        
                          ρ
                          ′
                        
                        ,
                        r
                        ,
                        
                          r
                          ′
                        
                      
                      )
                    
                  
                  ⁢
                  
                    L
                    ⁡
                    
                      (
                      
                        
                          ρ
                          ′
                        
                        ,
                        r
                      
                      )
                    
                  
                  ⁢
                  
                    
                    
                      ρ
                      ′
                    
                    
                  
                  ⁢
                  
                    d
                    
                      ρ
                      ′
                    
                  
                  ⁢
                  
                    
                    
                      r
                      ′
                    
                    
                  
                  ⁢
                  
                    d
                    
                      r
                      ′
                    
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        11
        )
      
    
  
  
    
      
        
          
            E
            ⁡
            
              (
              r
              )
            
          
          =
          
            π
            ⁢
            
              
                ∫
                
                  Ω
                  ρ
                
              
              ⁢
              
                
                  
                    L
                    ^
                  
                  ⁡
                  
                    (
                    
                      ρ
                      ,
                      r
                    
                    )
                  
                
                ⁢
                
                  
                  ρ
                  
                
                ⁢
                
                  d
                  ρ
                
              
            
          
        
        ,
      
    
    
      
        (
        12
        )
      
    
  




and the general form of the diffuser kernel becomes:





  
    
      
        
          D
          ⁡
          
            (
            
              ρ
              ,
              
                ρ
                ′
              
              ,
              r
              ,
              
                r
                ′
              
            
            )
          
        
        =
        
          
            
              δ
              ⁡
              
                (
                
                  ρ
                  -
                  
                    ρ
                    ′
                  
                
                )
              
            
            
              π
              ⁢
              
                
                
                  ρ
                  ′
                
                
              
            
          
          ⁢
          
            
              k
              ⁡
              
                (
                
                  r
                  -
                  
                    
                      r
                      ′
                    
                    ·
                    ρ
                  
                
                )
              
            
            
              π
              ⁢
              
                
                r
                
              
            
          
        
      
    
    
      
        (
        13
        )
      
    
  




The same box-shaped scattering function as was used for the diffuser kernel in Equation 7 can be used for Equation 13:




  
    
      
        
          k
          ⁡
          
            (
            
              r
              ,
              ρ
            
            )
          
        
        =
        
          
            1
            w
          
          ⊓
          
            
              (
              
                r
                w
              
              )
            
            .
          
        
      
    
    
      
        (
        14
        )
      
    
  




However, the physical interpretation of this diffuser is different than the previous diffuser. For the previous diffuser, each ray in the light field is scattered so that it spreads across a square on the sensor. The effect of the scattering function in Equation 14, however, is as illustrated in FIG. 6. As shown, in the absence of the diffuser, light from an annulus of width dp and radius p in the aperture plane projects to an annulus of width dr and radius r on the sensor. The effect of the scatter function in Equation 14 is to spread the light incident on the sensor so that it produces an annulus of width w instead.

As illustrated by volume 602 in FIG. 6, in polar coordinates, a ray can be a small annular section that travels from the aperture plane to the sensor plane. The effect of the diffuser, which is to scatter a ray along a radial line of width w, can be as illustrated by volume 604.
A box-shaped scatter function can be used here for notational convenience, but a Gaussian scattering function (e.g., as illustrated in FIG. 9(c)) can be superior for extended DOF imaging. The light field of a point source filtered by this diffuser kernel and PSF can be shown to be:




  
    
      
        
          
            
              L
              ^
            
            ⁡
            
              (
              
                ρ
                ,
                r
              
              )
            
          
          =
          
            
              4
              
                π
                ⁢
                
                  
                
                ⁢
                
                  A
                  2
                
              
            
            ⊓
            
              
                (
                
                  ρ
                  A
                
                )
              
              ⁢
              
                
                  ⊓
                  
                    (
                    
                      
                        r
                        -
                        
                          
                            s
                            0
                          
                          ⁢
                          ρ
                        
                      
                      w
                    
                    )
                  
                
                
                  π
                  ⁢
                  
                    
                  
                  ⁢
                  w
                  ⁢
                  
                    
                    r
                    
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        15
        )
      
    
  
  
    
      
        
          
            P
            ^
          
          ⁡
          
            (
            r
            )
          
        
        =
        
          
            4
            
              π
              ⁢
              
                
              
              ⁢
              
                s
                0
                4
              
              ⁢
              
                A
                2
              
            
          
          ⁢
          
            
              
                1
                
                  w
                  ⁢
                  
                    
                    r
                    
                  
                
              
              ⁡
              
                [
                
                  ⊓
                  
                    
                      (
                      
                        r
                        w
                      
                      )
                    
                    ⊗
                    
                      (
                      
                        ⊓
                        
                          
                            (
                            
                              r
                              
                                
                                  s
                                  0
                                
                                ⁢
                                A
                              
                            
                            )
                          
                          ·
                          
                            
                            r
                            
                          
                        
                      
                      )
                    
                  
                
                ]
              
            
            .
          
        
      
    
    
      
        (
        16
        )
      
    
  




The analytic solution for this PSF is a piecewise function due to the contribution from the term in brackets, which is a convolution between the two rect functions (one weighted by |r|). Note that as the scattering width w is reduced to zero, the first rect (combined with 1/w) approaches a delta function and the result is a pillbox-shaped defocus PSF. Also note that if a different diffuser with different scattering function is used, the first rect is simply replaced with the new scattering function. However, the convolution term is far less significant than the 1/|r| term, whose effect dominates, resulting in a PSF which can be strongly depth independent while still maintaining a strong peak and preserving high frequencies.

As illustrated in FIG. 6, light incident on a small annular region of width δr and radius r emanates from an annulus in the aperture, and its energy can be proportional to ρ or equivalently to r/s0. This explains the presence of the |r| multiplier within the term in brackets of Equation 16. This term in brackets states that the energy in a pillbox defocus PSF annulus is spread uniformly along radial lines of width w by the diffuser, as shown on the right hand side of FIG. 6. The 1/|r| term in Equation 16 can be attributed to the fact that the energy density becomes larger for light that is scattered closer to the center of the PSF.
FIG. 7 shows several PSF 702 and Modulation Transfer Function (MTF) 704 graph pairs for a camera with (714, 716, 718, 720, 722, 724, 726, and 728) and without (715, 717, 719, 721, 723, 725, 727, and 729) the diffuser given by Equation 16. The defocus blur diameter s0A changes between 0 pixels 706, 25 pixels 708, 50 pixels 710, and 100 pixels 712. The scatter function of Equation 14 is a Gaussian instead of a box function, and the diffuser parameter w (the variance of the Gaussian) is chosen so that w=100 pixels. Note that when the diffuser is present, there is little variation with depth for either the PSF or MTF. Introducing the diffuser also eliminates the zero crossings in the MTF. For smaller defocus values, the diffuser suppresses high frequencies in the MTF. However, because the diffuser MTF does not vary significantly with depth, high frequencies can be recovered via deconvolution.
In accordance with some embodiments, diffusers of the “kinoform” type (as described in Caufield, H. J., “Kinoform Diffusers,” In SPIE Conference Series, vol. 25, p. 111, 1971, which is hereby incorporated by reference herein in its entirety) where the scattering effect is caused entirely by roughness variations across a surface can be used. Such a diffuser can be considered to be a random phase screen, and according to statistical optics, for a camera with effective focal length f, and center wavelength λ, the effect of placing this screen in the aperture of the camera can result in the following:




  
    
      
        
          
            
              P
              ^
            
            ⁡
            
              (
              
                x
                ,
                y
              
              )
            
          
          ∝
          
            
              p
              
                
                  φ
                  u
                
                ,
                
                  φ
                  v
                
              
            
            ⁡
            
              (
              
                
                  x
                  
                    
                      λ
                      _
                    
                    k
                  
                
                ,
                
                  y
                  
                    
                      λ
                      _
                    
                    k
                  
                
              
              )
            
          
        
        ,
      
    
    
      
        (
        18
        )
      
    
  




where φu and φv are the u and v derivatives of the phase shift induced by the surface, and pφx,φy is the joint probability of these derivatives. The result of Equation 18 is that a diffuser can be implemented by creating an optical element with thickness t(u,v), where the gradient of this surface ∇t(u,v) is sampled from a probability distribution which is also a desired PSF. Intuitively, this equation can be understood as follows: pφu,φv denotes the fraction of the surface t(u,v) with slope (φu, φv). For small angles, all incoming rays incident on this fraction of the surface will be deflected at the same angle, since the slope is constant over this region. Thus, the quantity pφu,φv also reflects the portion of light that will be deflected by the slope (φx,φy).

A Kinoform diffuser has a randomly varying surface with a general probability distribution of slopes as illustrated in FIG. 8(b). Kinoform diffusers can be thought of as generalized phase plates. For example, a regular deterministic phase plate with thickness t(u)=αλu, as shown in FIG. 8(a), can also be thought of as having a slope drawn from a probability function p(φu) which is a delta function. The result of placing this phase plate in the pupil plane of a camera is to shift the PSF, which can be thought of as convolving p(φu) with the PSF.
To implement the diffuser defined in Equation 14, the diffuser surface can be implemented as a sequence of quadratic elements whose diameter and sag is drawn from a random distribution as described in Sales, T. R. M., “Structured microlens arrays for beam shaping,” Optical Engineering 42, 11, pp. 3084-3085, 2003, which is hereby incorporated by reference herein in its entirety. The scatter function of the diffuser can be designed to be roughly Gaussian with 0.5 mm variance (corresponding to w=1 mm in Equation 16) as shown in FIG. 9(c). To create a radially symmetric diffuser, a one-dimensional random profile can be created and then a polar transformation applied to create a two-dimensional surface (see, e.g., FIGS. 9(a) and 9(b)).
In some embodiments, a diffuser can be made using laser etching.
In some embodiments, the maximum height of the diffuser surface can be 3 μm, and the diffuser can be fabricated using a laser machining technology which has a minimum spot size of about 10 μm. To ensure that each quadratic element in the diffuser is fabricated with high accuracy, the minimum diameter of a single element can be chosen to be 200 μm, resulting in a diffuser with 42 different annular sections.
Any suitable hardware can be used to implement a mechanism 102 in accordance with some embodiments. For example, a Canon EOS 450D sensor from Canon U.S.A., Inc. can be used as sensor 116, a 22 mm diameter diffuser (e.g., as illustrated in FIG. 9(d)) that is laser etched in a piece of suitable optical glass by RPC Photonics of Rochester, N.Y. can be used as diffuser 110 or 112, and a 50 mm f/1.8 lens from Canon U.S.A., Inc. can be used as lens 114. As another example, lens 114 can have any focal length and consist of refractive optics, reflective optics, or both. For instance, a 3048 mm focal length Meade LX200 telescope (available from) can be used in some embodiments.
In accordance with some embodiments, any suitable processing can be performed to deblur the image hitting a camera sensor after passing through a lens and diffuser (in either order). For example, the Wiener deconvolution with the PSF at the center depth can be used to deblur the sensed images. Any suitable additional or alternative processing on the images can be used. For example, additional deblurring of diffusion coded images can performed using the BM3D deblurring algorithm as described in Dabov, K., Foi, A., Katkovnik, V., and Egiazarian, K., “Image restoration by sparse 3D transform-domain collaborative filtering,” In SPIE Conference Series, vol. 6812, 681207, 2008, which is hereby incorporated by reference herein in its entirety. In some embodiments, the BM3D deblurring algorithm enforces a piecewise smoothness prior that suppresses the noise amplified by the deblurring process.
Any suitable hardware processor, such as a microprocessor, digital signal processor, special purpose computer (which can include a microprocessor, digital signal processor, a controller, etc., memory, communication interfaces, display controllers, input devices, etc.), general purpose computer suitably programmed (which can include a microprocessor, digital signal processor, a controller, etc., memory, communication interfaces, display controllers, input devices, etc.), server, programmable gate array, etc. can be used to deblur the image captured by the sensor. Any suitable hardware can by used to transfer the image from the sensor to the processor. Any suitable display, storage device, or printer can then be used to display, store, or print the deblurred image.
In some embodiments, any suitable computer readable media can be used for storing instructions for performing the processes described herein. For example, in some embodiments, computer readable media can be transitory or non-transitory. For example, non-transitory computer readable media can include media such as magnetic media (such as hard disks, floppy disks, etc.), optical media (such as compact discs, digital video discs, Blu-ray discs, etc.), semiconductor media (such as flash memory, electrically programmable read only memory (EPROM), electrically erasable programmable read only memory (EEPROM), etc.), any suitable media that is not fleeting or devoid of any semblance of permanence during transmission, and/or any suitable tangible media. As another example, transitory computer readable media can include signals on networks, in wires, conductors, optical fibers, circuits, any suitable media that is fleeting and devoid of any semblance of permanence during transmission, and/or any suitable intangible media.
Although the invention has been described and illustrated in the foregoing illustrative embodiments, it is understood that the present disclosure has been made only by way of example, and that numerous changes in the details of implementation of the invention can be made without departing from the spirit and scope of the invention, which is only limited by the claims which follow. Features of the disclosed embodiments can be combined and rearranged in various ways.

## Claims

1. A system for recording an image of a scene, comprising:
an image sensor that receives light representing the scene and generates image data representing a first image based on light representing the scene;
a radially symmetric diffuser having a first side that receives light representing the scene and a second side and a one-dimensional thickness profile of the radially symmetric diffuser along a radial direction is random such that each ray of light passing through the radially symmetric diffuser from the first side to the second side forms a section of an annulus on the image sensor having a width w that is the difference between an inner radius and an outer radius of the annulus, such that w is larger than a width of a section of annulus that the ray of light would form in the absence of the radially symmetric diffuser, wherein width w is substantially constant regardless of where along the radial direction the ray of light passes through the radially symmetric diffuser; and
a hardware processor that uses a point spread function corresponding to width w to generate a second image based on the image data representing the first image and that corresponds to a deblurred version of the first image.

2. The system of claim 1, wherein the diffuser is made by laser etching.

3. The system of claim 1, wherein the diffuser is of the kinoform type.

4. The system of claim 1, wherein a scattering function of the diffuser is a symmetric distribution function.

5. The system of claim 4, wherein the symmetric distribution function of the diffuser is perturbed Gaussian.

6. The system of claim 1, further comprising a lens that passes the light representing the scene before it is incident on the diffuser.

7. The system of claim 1, further comprising a lens that is positioned between the diffuser and the sensor.

8. The system of claim 1, further comprising a display that displays the deblurred image.

9. A method for recording an image of a scene, comprising:
diffusing light representing the scene using a radially symmetric diffuser having a first side that receives light representing the scene and a second side and a one-dimensional thickness profile of the radially symmetric diffuser along a radial direction is random such that each ray of light passing through the radially symmetric diffuser from the first side to the second side forms a section of an annulus on the image sensor having a width w that is the difference between an inner radius and an outer radius of the annulus, such that w is larger than a width of a section of annulus that the ray of light would form in the absence of the radially symmetric diffuser, wherein width w is substantially constant regardless of where along the radial direction the ray of light passes through the radially symmetric diffuser; and
receiving light representing the scene using an image sensor;
generating image data representing a first image based on the light representing the scene received by the image sensor;
using a point spread function corresponding to width w to generate a second image based on the image data representing the first image that corresponds to a deblurred version of the first image.

10. The method of claim 9, wherein the diffuser is made by laser etching.

11. The method of claim 9, wherein the diffuser is of the kinoform type.

12. The method of claim 9, wherein a scattering function of the diffuser is a symmetric distribution function.

13. The method of claim 12, wherein the symmetric distribution function of the diffuser is perturbed Gaussian.

14. The method of claim 9, further comprising positioning a lens so that the light representing the scene passes through the lens before it is incident on the diffuser.

15. The method of claim 9, further comprising positioning a lens between the diffuser and the sensor.

16. The method of claim 9, further comprising displaying the deblurred image.