WO2007054938A2 - Optics for an extended depth of field - Google Patents

Optics for an extended depth of field Download PDF

Info

Publication number
WO2007054938A2
WO2007054938A2 PCT/IL2006/001294 IL2006001294W WO2007054938A2 WO 2007054938 A2 WO2007054938 A2 WO 2007054938A2 IL 2006001294 W IL2006001294 W IL 2006001294W WO 2007054938 A2 WO2007054938 A2 WO 2007054938A2
Authority
WO
WIPO (PCT)
Prior art keywords
optical
assembly
image
lenses
design
Prior art date
Application number
PCT/IL2006/001294
Other languages
French (fr)
Other versions
WO2007054938A3 (en
Inventor
Alex Alon
Irina Alon
Haim Bezdin
Original Assignee
D-Blur Technologies Ltd.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from US11/278,255 external-priority patent/US20070236573A1/en
Application filed by D-Blur Technologies Ltd. filed Critical D-Blur Technologies Ltd.
Priority to US12/093,394 priority Critical patent/US7773316B2/en
Priority to TW096108966A priority patent/TWI427554B/en
Publication of WO2007054938A2 publication Critical patent/WO2007054938A2/en
Publication of WO2007054938A3 publication Critical patent/WO2007054938A3/en
Priority to US12/847,651 priority patent/US8014084B2/en
Priority to US13/196,051 priority patent/US8294999B2/en
Priority to US13/625,390 priority patent/US8611030B2/en

Links

Classifications

    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/0075Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00 with means for altering, e.g. increasing, the depth of field or depth of focus

Definitions

  • the present invention relates generally to digital imaging, and specifically to optics and methods for designing digital cameras with an enhanced depth of field.
  • the objective optics used in digital cameras are typically designed so as to minimize the optical point spread function (PSF) and maximize the modulation transfer function (MTF), subject to the limitations of size, cost, aperture size, and other factors imposed by the camera manufacturer.
  • the PSF of the resulting optical system may still vary from the ideal due to focal variations and aberrations.
  • a number of methods are known in the art for measuring and compensating for such PSF deviations by digital image processing.
  • U.S. Patent 6,154,574 whose disclosure is incorporated herein by reference, describes a method for digitally focusing an out-of-focus image in an image processing system.
  • a mean step response is obtained by dividing a defocused image into sub-images, and calculating step responses with respect to the edge direction in each sub-image.
  • the mean step response is used in calculating PSF coefficients, which are applied in turn to determine an image restoration transfer function.
  • An in-focus image is obtained by multiplying this function by the out-of-focus image in the frequency domain.
  • U.S. Patent 6,567,570 whose disclosure is incorporated herein by reference, describes an image scanner, which uses targets within the scanner to make internal measurements of the PSF. These measurements are used in computing convolution kernels, which are applied to images captured by the scanner in order to partially compensate for imperfections of the scanner lens system.
  • an optical imaging assembly that may be used in a digital camera is constructed so as to generate a distorted image at an imaging plane of the assembly.
  • the digital camera incorporates a sensor at the imaging plane, and the distorted image formed on the sensor is corrected by a deconvolution engine so as to produce an undistorted image.
  • the optical imaging assembly is configured to produce a high defocus aberration coefficient, greater than approximately 0.1, across the distorted image when the image is optimally focused to a focal plane of the assembly.
  • the high defocus aberration coefficient causes a modulation transfer function (MTF) for the assembly to have generally equal low values for all objects in a large field, typically from infinity to approximately 10 cm from the assembly.
  • MTF modulation transfer function
  • the deconvolution engine may be configured to improve the MTF of all the images at the different object distances.
  • the final images thus produced by the assembly and the deconvolution engine are substantially free of aberrations for all objects within the field, so that the digital camera has a depth of field which is significantly higher than depths of field of prior art cameras.
  • the imaging assembly is designed by assuming that a virtual phase plate is added to the assembly, the virtual phase plate having characteristics that generate aberration coefficients which are the negatives of those required by the deconvolution engine.
  • the assembly with the phase plate is configured to give an undistorted image.
  • an optical imaging assembly having cylindrical symmetry, including a plurality of lenses having surfaces with curvatures and spacings between the surfaces, such that an optical image formed by the plurality of lenses has a defocus aberration coefficient greater than 0.1 at a focal plane of the assembly.
  • the defocus aberration coefficient may be greater than or equal to 0.6, and less than or equal to 0.9.
  • a fifth order spherical aberration is typically between -0.2 and -0.3, a seventh order spherical aberration may be -0.5, and a ninth order spherical abei ⁇ ation may be 0.
  • the curvatures are non-logarithmic.
  • the optical image is of one or more objects within a field of view subtending up to 33° to an optic axis of the assembly.
  • the optical image is of one or more objects which are located at distances greater than or equal to 10 cm from the assembly.
  • the assembly is configured so that a modulation transfer function generated by the one or more objects at the focal plane is less than or equal to 0.35 for spatial frequencies greater than 50 cycles per mm.
  • the optical image has an extent of a point spread function equal to or greater than 14 ⁇ m at the focal plane.
  • a method for optical imaging including: formulating an optical specification of an optical imaging assembly including optical elements having surfaces, the optical specification including one or more aberration coefficients; generating a virtual phase element having phase coefficients which are negatives of the aberration coefficients; generating a design of the optical assembly including the optical elements and the virtual phase element so as to determine curvatures and spacings of the surfaces that will focus an image of an object on a focal plane of the assembly in the presence of the virtual phase element; and fabricating the optical assembly in accordance with the curvatures and spacings determined by the design.
  • the method includes generating a digital filter to compensate for the one or more aberration coefficients, and coupling the digital filter to the optical assembly so that a filtered image formed by the digital filter from the image on the focal plane is free of aberrations.
  • the filtered image is free of aberrations for the object being located at a distance greater than 10 cm from the optical assembly.
  • an optical imaging assembly having cylindrical symmetry, including a plurality of lenses having surfaces with curvatures and spacings between the surfaces, such that an optical image formed by the plurality of lenses at a focal plane of the assembly has a modulation transfer function (MTF) for spatial frequencies greater than 50 cycles per millimeter that is less than or equal to 0.35.
  • MTF modulation transfer function
  • a method for forming an optical imaging assembly including: providing a plurality of lenses having cylindrical symmetry and surfaces with curvatures and spacings between the surfaces; and arranging the lenses so that an optical image formed by the plurality of lenses has a defocus aberration coefficient greater than 0.1 at a focal plane of the assembly.
  • apparatus for optical imaging including: an optical imaging assembly including optical elements having surfaces; and a processor which is configured to: formulate an optical specification of the optical imaging assembly, the optical specification including one or more aberration coefficients, generate a virtual phase element having phase coefficients which are negatives of the aberration coefficients; generate a design of the optical assembly including the optical elements and the virtual phase element so as to determine curvatures and spacings of the surfaces that will focus an image of an object on a focal plane of the assembly in the presence of the virtual phase element; and fabricate the optical assembly in accordance with the curvatures and spacings determined by the design.
  • a method for forming an optical imaging assembly including: providing a plurality of lenses having cylindrical symmetry and surfaces with curvatures and spacings between the surfaces; and arranging the lenses so that an optical image formed by the plurality of lenses at a focal plane of the assembly has a modulation transfer function (MTF) for spatial frequencies greater than 50 cycles per millimeter that is less than or equal to 0.35.
  • an optical imaging assembly including: an aperture having a semi-diameter of 7.5-10"! mm; and a plurality of lenses having surfaces with curvatures and spacings between the surfaces, wherein respective surfaces of the plurality of lenses are defined by an equation
  • z is a sag of the surface
  • r is a radius measured to an optic axis of the assembly
  • c, k are respective curvature and conic constants for the surface
  • ⁇ 2 , ... ⁇ g are parameters for the surface
  • Fig. 1 is a block diagram that schematically illustrates a digital camera, in accordance with an embodiment of the present invention
  • Fig. 2 is a schematic, pictorial illustration of a system for designing a digital camera, in accordance with an embodiment of the present invention
  • Fig. 3 A is a schematic, pictorial illustration showing conceptual elements of a digital camera used in a design process, in accordance with an embodiment of the present invention
  • Fig. 3B is a plot of modulation transfer functions (MTF) for a digital camera with and without application of a deconvolution filter, in accordance with an embodiment of the present invention
  • Fig. 4 is a flow chart that schematically illustrates a method for designing a digital camera, in accordance with an embodiment of the present invention
  • Figs. 5A-5C are schematic, perspective plots of DCF kernels used in a digital camera, in accordance with an embodiment of the present invention
  • Fig. 6 is an image that simulates the output of an image sensor using objective optics with specifications that have been relaxed in accordance with an embodiment of the present invention
  • Fig. 7 is an image that simulates the effect of application of an appropriate DCF to the image of Fig. 6, in accordance with an embodiment of the present invention
  • Fig. 8 is a schematic pictorial representation of a camera imaging a large depth of field, in accordance with an embodiment of the present invention
  • FIG. 9 schematically shows detail of a camera sensor, in accordance with an embodiment of the present invention.
  • Fig. 10 is a flowchart schematically describing procedures performed in designing a digital camera, in accordance with an embodiment of the present invention.
  • Figs. 1 IA-I IE schematically show results generated in the procedures of Fig. 10, in accordance with an embodiment of the present invention
  • Figs. 12A-12G show Modulation Transfer Function curves for an exemplary optical imaging assembly, in accordance with an embodiment of the present invention.
  • Fig. 13 is a schematic diagram of the exemplary optical imaging assembly, in accordance with an embodiment of the present invention.
  • Pitch of a detector array refers to the center-to-center distance between elements of the array.
  • Cylindrical symmetry describes a structure, such as a simple or compound lens, which has an optical axis such that the structure is invariant under rotation about the optical axis for any and all angles of rotation.
  • Point spread function is the impulse response of an optical system in the spatial domain, i.e., the image formed by the system of a bright point object against a dark background.
  • Extent of the PSF is the full width at half maximum (FWHM) of the PSF.
  • OTF Optical transfer function
  • Modulation transfer function is the modulus of the OTF.
  • Optical radiation refers to electromagnetic radiation in any of the visible, infrared and ultraviolet regions of the spectrum.
  • SYSTEM OVERVIEW Fig. 1 is a block diagram that schematically illustrates a digital camera 20, in accordance with an embodiment of the present invention.
  • the camera comprises refractive objective optics 22, which focus an image onto an image sensor 24, so that a surface of the sensor acts as a focal plane of the optics.
  • Optics 22 are also referred to herein as optical imaging assembly 22.
  • Optics 22 are designed in an iterative process together with a deconvolution engine 26 that operates on image data that are output by image sensor 24.
  • the deconvolution engine applies one or more digital filters, typically comprising at least one deconvolution filter (DCF), to the image data.
  • DCF deconvolution filter
  • the design process and method of filtering are described in detail hereinbelow.
  • the DCF kernel is typically chosen so as to correct for blur in the image formed by optics 22.
  • the image data are processed by an image signal processor (ISP) 28, which performs standard functions such as color balance and format conversion and outputs the resulting image.
  • ISP
  • Fig. 1 The optical and digital processing schemes illustrated in Fig. 1 are shown here solely for the sake of example, as an aid to understanding the techniques and tools that are described hereinbelow.
  • the principles of the present invention may be applied in conjunction with a wide variety of electronic imaging systems, using substantially any sort of optical design and substantially any type of image sensor, including both two- dimensional detector matrices and linear detector arrays, as are known in the art.
  • Deconvolution engine 26 and ISP 28 may be implemented as separate devices or as a single integrated circuit component. In either case, the deconvolution engine and ISP are typically combined with other I/O and processing elements, as are known in the art.
  • the term "digital camera” should therefore be understood as referring to any and all sorts of electronic imaging systems that comprise an image sensor, objective optics for focusing optical radiation onto the image sensor, and electronic circuits for processing the sensor output.
  • Fig. 2 is a schematic, pictorial illustration showing a system 30 for designing a digital camera, in accordance with an embodiment of the present invention.
  • System comprises a digital processing design station 32 and an optical design station 34.
  • Processing design station 32 receives a camera specification as input, from a camera manufacturer, for example, specifying the key dimensions, sensor type and desired optical characteristics (referred to hereinafter as the target optical specification) of the camera.
  • the specified optical characteristics may include, for example, the number of optical elements, materials, tolerances, focal length, magnification, aperture (F-number), depth of field, and resolution performance.
  • the optical resolution performance is typically defined in. terms of the MTF, but it may alternatively be specified in terms of PSF, wavefront quality, aberrations, and/or other measures of optical and image quality that are known in the art.
  • Processing design station 32 analyzes and modifies the target optical specification, taking into account the expected operation of engine 26, in order to provide a modified optical specification to the optical design station.
  • both the original camera specification and the modified optical specification use cylindrically-symmetrical optical elements. Specialized phase plates or other elements that break the cylindrical symmetry of the optics are generally undesirable, due to their added cost, and engine 26 is able to correct the aberrations of optics 22 without requiring the use of such elements.
  • processing design station 32 may compute and provide to optical design station 34 a merit function, indicating target values of the aberrations of optics 22 or scoring coefficients to be used in weighting the aberrations in the course of optimizing the optical design.
  • the aberrations express deviations of the optical wavefront created by optics 22 from the ideal, and may be expressed, for example, in terms of Zernike polynomials or any other convenient mathematical representation of the wavefront that is known in the art.
  • Optical design station 34 is typically operated by a lens designer, in order to produce a lens design according to the modified optical specification provided by processing design station 32.
  • the processing design station determines the optimal DCF (and possibly other filters) to be used in engine 26 in conjunction with this lens design.
  • the DCF computation is tied to the specific lens design in question so that the filter coefficients reflect the "true" PSF of the actual optical system with which the DCF is to be used.
  • the processing design station then evaluates the optical design together with the DCF in order to assess the combined result of the expected optical quality of optics 22 and the enhancement expected from engine 26, and to compare the result to the target optical specification.
  • the assessment may take the form of mathematical analysis, resulting in a quality score.
  • a quality scoring schemes that may be used in this context is described hereinbelow.
  • other quality scoring schemes may be used, such as that described, for example, in the above-mentioned PCT publication WO 2004/063989 A2.
  • station 32 may generate and display a simulated image 36, which visually demonstrates the output image to be expected from the camera under design based on the current choice of optical specifications and DCF.
  • the processing design station may perform further design iterations internally, or it may generate a further modified optical specification, which it passes to optical design station 34 for generation of a modified optical design. This process may continue iteratively until a suitable optical design and DCF are found. Details of this process are described hereinbelow with reference to Fig. 4.
  • stations 32 and 34 comprise general-purpose computers running suitable software to carry out the functions described herein.
  • the software may be downloaded to the computers in electronic form, over a network, for example, or it may alternatively be furnished on tangible media, such as optical, magnetic, or electronic memory media.
  • some of the functions of stations 32 and/or 34 may be implemented using dedicated or programmable hardware components.
  • the functions of optical design station 34 may be carried out using off-shelf optical design software, such as ZEMAX® (produced by ZEMAX Development Corp., San Diego, California).
  • ZEMAX® produced by ZEMAX Development Corp., San Diego, California
  • Fig. 3A is a schematic, pictorial illustration showing conceptual elements of camera 20, as they are applied in the design process used in system 30, in accordance with an embodiment of the present invention.
  • System 30 takes engine 26 into account in the design of optics 22, as explained hereinabove, and thus relates to the DCF as a sort of "virtual lens" 40.
  • the design constraints on the actual objective optics are relaxed by the use of this virtual lens, as though the optical designer had an additional optical element to incorporate in the design for purposes of aberration correction.
  • the virtual lens that is implemented in engine 26 is chosen, in conjunction with the actual optical lenses, to give an image output that meets the manufacturer's camera specifications.
  • Fig. 3B is a plot showing the MTF of a camera designed using system 30, in accordance with an embodiment of the present invention.
  • the plot includes an uncorrected curve 44, corresponding to the modified optical specification generated by station 32 for use in designing optics 22 on station 34.
  • the low MTF permitted by curve 44 is indicative of the expected improvement in MTF that can be achieved by use of DCF 26.
  • a corrected curve 46 shows the net MTF of the camera that is achieved by applying the DCF to the image sensor output. These curves show the MTF at the center of the optical field, with the object at a certain distance from the camera.
  • the MTF may be specified at multiple different focal depths and field angles.
  • Figs. 3 A and 3 B permits the camera manufacturer to achieve the desired level of optical performance with fewer, smaller and/or simpler optical components than would be required to achieve the same result by optical means alone. Additionally or alternatively, the camera may be designed for enhanced performance, such as reduced aberrations, reduced F-number, wide angle, macro operation, or increased depth of field.
  • Fig. 4 is a flow chart that schematically illustrates a method for designing a digital camera, in accordance with an embodiment of the present invention. The method will be described hereinbelow, for the sake of clarity, with reference to camera 20 and system 30, although the principles of this method may be applied generally to other cameras and using other design systems.
  • Processing design station 32 translates the target optical specification of the camera into a modified optical specification, at a specification translation step 50.
  • station 32 uses an estimate of the DCF to be implemented in the camera.
  • the image enhancement to be expected due to this DCF is then applied to the optical specification in order to estimate how far the optical design parameters, such as the MTF, can be relaxed.
  • Image enhancement by the DCF tends to amplify noise in the output of image sensor 24.
  • the noise gain NG is proportional to the norm of the DCF (-vD fc D , wherein D is the DCF kernel and the superscript t indicates the He ⁇ nitian transpose).
  • the processing design station uses the maximum permissible noise gain as a limiting condition.
  • engine 26 may also comprise a noise filter.
  • the limit placed on the DCF coefficients by the noise gain may thus be mitigated by the noise reduction that is expected due to the noise filter.
  • the norm of the DCF kernel is approximately given by the product of the maximum permissible noise gain with the expected noise reduction factor (i.e., the ratio of image noise following the noise filter to image noise without noise filtering).
  • a more accurate estimate of the overall noise gain may be obtained by taking the norm of the product of the noise filter multiplied by the DCF in the frequency domain.
  • the OTF may be assumed, at first approximation, to be linear as a function of spatial frequency q, which is normalized to the Nyquist frequency of image sensor 24:
  • OTF 1 - ⁇ q ⁇ q ⁇ 1 / ⁇
  • the PSF may be determined analytically from the OTF of equation (1). Because of the zeroes in the OTF, the frequency-domain representation of the DCF to be used in the camera may be estimated as:
  • the maximal permissible noise gain may be determined by comparing the expected noise characteristic of image sensor 24 to the noise specification. As noted above, digital smoothing of the noise in the output image may also be taken into account in order to permit the constraint on noise gain in the DCF to be relaxed.
  • noise removal methods may be used in engine 26.
  • a morphological operation may be used to identify edges in the image, followed by low-pass filtering of non-edge areas.
  • the choice of noise removal method to be used in engine 26, however, is beyond the scope of the present invention.
  • Equations (4) and (5) may be used in estimating how far the MTF of optics 22 may be reduced relative to the original target specification, subject to a given noise gain limit. This reduction factor may be applied, for example, to the MTF required by the original camera specification at a benchmark frequency, such as half the Nyquist frequency. In the example shown in Fig. 3B, the target MTF has been reduced to about 1/3 of its original specified value. The MTF will be restored in the output image from camera 20 by operation of the DCF in engine 26.
  • processing design station 32 may also generate a merit function at step 50 for use by the optical designer.
  • the merit function may take the form of aberration scores, which are assigned to each significant aberration that may characterize optics 22.
  • the aberrations may be expressed, for example, in terms of the Zernike polynomials, for each of the colors red, green and blue individually.
  • Standard software packages for optical design such as ZEMAX, are capable of computing the Zernike polynomial coefficients for substantially any design that they generate.
  • Values of the merit functions may be provided in tabular form. Generation of these values is described in detail in the above-mentioned PCT Publication WO 2004/063989 A2.
  • processing design station 32 may generate target wavefront characteristics that the optical design should achieve in the image plane (i.e., the plane of sensor 24). These wavefront characteristics may conveniently be expressed in terms of values of the aberrations of optics 22, such as Zernike coefficient values. Typically, aberrations that can be corrected satisfactorily by deconvolution engine 26 may have high values in the optical design, whereas aberrations that are difficult to correct should have low values. In other words, an aberration that would have a high score in the merit function will have a low target value, and vice versa.
  • the target aberration values can be seen as the inverse of the wavefront corrections that can be achieved by "virtual lens" 40.
  • the target aberration values may also include aberrations that reduce the sensitivity of the optics to various undesirable parameters, such as manufacturing deviations and defocus.
  • An optical designer working on station 34 uses the specification, along with the merit function and/or aberration target values provided at step 50, in generating an initial design of optics 22, at an optical design step 52.
  • the designer may use the merit function in determining a design score, which indicates how to trade off one aberration against another in order to generate an initial design that maximizes the total of the merit scores subject to the optical specification.
  • the optical designer may insert a dummy optical element, with fixed phase characteristics given by the target aberration values as an additional element in the optical design. This dummy optical element expresses the wavefront correction that is expected to be achieved using engine 26 and thus facilitates convergence of the calculations made by the optical design software on station 34 to the desired design of the elements of optics 22.
  • Control of the design process now passes to processing design station 32, in a design optimization stage 53.
  • the processing design station analyzes the optical design, at a design analysis step 54.
  • the analysis at this step may include the effect of virtual lens
  • station 32 typically computes the optical performance of the optics as a function of wavelength and of location in the image plane. For example, station 32 may perform an accurate ray trace computation based on the initial optical design in or to calculate a phase model at the image plane, which may be expressed in terms of Zernike polynomial coefficients.
  • the total aberration - and hence the PSF - at any point in the image plane may be obtained from the total, wavefront aberration, which is calculated by summing the values of the Zernike polynomials.
  • Station 32 determines a design quality score, at a scoring step 55.
  • this score combines the effects of the PSF on image resolution and on artifacts in the image, and reflects the ability of engine 26 to compensate for these effects.
  • the score measures the extent to which the current optical design, taken together with filtering by engine 26, will satisfy the camera specification that was originally provided as input to station 32 as input to step 50.
  • the score computed at step 55 is based on the camera specification and on a set of weights assigned to each parameter in the camera specification.
  • the camera specification is expressed in a list of desired parameter values at various image plane locations and wavelengths, such as:
  • the weight assigned to each parameter is typically determined by its scaling, subjective importance, and likelihood of satisfying the desired parameter value relative to other parameters.
  • the overall score is computed by summing the weighted contributions of all the relevant parameters. In this embodiment, if a given parameter is within the specified range, it makes no contribution to the score. If the value is outside the specified range, the score is decreased by the square difference between the parameter value and the closest permissible value within the specified range, multiplied by the appropriate weight. A design that fully complies with the camera specification will thus yield a zero score, while non-compliance will yield negative values. Alternatively, other parameters and other methods may be used in computing numerical values representing how well the current design satisfies the camera specification.
  • the score computed at step 55 is assessed to determine whether it indicates that the current design is acceptable, at a quantitative assessment step 56. If the design does not meet the specification, station 32 modifies the optical design parameters at an optimization step 58. For this purpose, the station may estimate the effects of small changes in the aberrations on the PSF. This operation gives a multi-dimensional gradient, which is used in computing a change to be made in the optical design parameters by linear approximation. The DCF parameters may be adjusted accordingly. A method for computing and using gradients of this sort is described, for example, in . the above- mentioned PCT Publication WO 2004/063989 A2. The results of step 58 are input to step 54 for recomputation of the optical performance analysis.
  • the process continues iteratively through steps 55 and 56 until the design quality score reaches a satisfactory result.
  • the design parameters are presented by processing design station 32 to the system operator, at a design checking step 60.
  • the system operator reviews the optical design (as modified by station 32 in step 58, if necessary), along with the results of the design analysis performed at step 54. Additionally or alternatively, the optical design and DCF may be used at this point in generating a simulated output image, representing the expected performance of the camera in imaging a known scene or test pattern. (Exemplary simulated images of this sort are shown below in Figs. 6 and 7.)
  • the system operator reviews the design in order to verify that the results are indeed satisfactory for use in manufacture of camera 20.
  • the operator may change certain parameters, such as specification parameters and/or scoring weights, and return to stage 53.
  • certain parameters such as specification parameters and/or scoring weights
  • the operator may initiate changes to the • original camera specification and return the process to step 50. This sort of operator involvement may also be called for if stage 53 fails to converge to an acceptable score at step 56.
  • processing design station 32 uses processing design station 32 to generate tables of values to be used in camera 20, at a DCF creation step 62.
  • the DCF tables vary according to location in the image plane.
  • a different DCF kernel is computed for each region of 50 x 50 pixels in image sensor 24.
  • sensor 24 is a color image sensor
  • different kernels are computed for the different color planes of sensor 24.
  • common mosaic image sensors may use a Bayer pattern of red, green, and blue pixels 42.
  • the output of the image sensor is an interleaved stream of sub-images, comprising pixel samples belonging to different, respective colors.
  • DCF 26 applies different kernels in alternation, so that the pixels of each color are filtered using values of other nearby pixels of the same color. Appropriate kernel arrangements for performing this sort of filtering are described in U.S. Provisional Patent Application 60/735,519, filed November 10, 2005, which is assigned to the assignee of the present patent application and is incorporated herein by reference.
  • Figs. 5A, 5B and 5C are schematic, perspective plots of DCF kernels 70, 72 and 74 for red, green, and blue pixels, respectively, which are computed in accordance with an embodiment of the present invention.
  • Each kernel extends over 15 x 15 pixels, but contains non-zero values only at pixels of the appropriate color.
  • red kernel 70 for example, in each square of four pixels, only one - the red pixel - has a nonzero value.
  • Blue kernel 74 is similarly constructed, while green kernel 72 contains two non-zero values in each four-pixel square, corresponding to the greater density of green pixels in the Bayer matrix.
  • the central pixel has a large positive value, while surrounding values are lower and may include negative values 76.
  • the DCF values are chosen so that the norm does not exceed the permitted noise gain.
  • design station 32 uses the DCF tables from step 62 and the optical design output from stage 53 in simulating the performance of camera 20, at a simulation step 64.
  • the simulation may also use characteristics, such as noise figures, of image sensor 24 that is to be installed in the camera, as well as other factors, such as manufacturing tolerances to be applied in producing the camera and/or operation of ISP 28.
  • the results of this step may include simulated images, like image 36 (Fig. 2), which enable the system operator to visualize the expected camera performance.
  • Figs. 6 and 7 are images that simulate the expected output of camera 20, as may be generated at step 64, in accordance with an embodiment of the present invention.
  • Fig. 6 shows a standard test pattern as it would be imaged by optics and captured by image sensor 24, without the use of DCF 26.
  • the image of the test pattern is blurred, especially at higher spatial frequencies, due to the low MTF of camera 20. (The MTF is given roughly by uncorrected curve 44 in Fig. 3B.)
  • the image pixels are decimated due to the use of a color mosaic sensor, and random noise is added to the image corresponding to the expected noise characteristics of the image sensor.
  • Fig. 7 shows the image of Fig. 6 after simulation of processing by DCF 26, including noise removal as described hereinbelow.
  • the MTF of this image is given roughly by curve 46 in Fig. 3B.
  • the aliasing apparent in the images of the high- frequency test patterns is the result of a true simulation of the performance of a low- resolution image sensor following DCF processing.
  • the system operator, viewing this image, is able to ascertain visually whether the camera performance will meet the original camera specifications that were provided at step 50.
  • system operator's visual assessment is combined with the numerical results of the design analysis, in order to determine whether the overall performance of the design is acceptable, at an acceptance step 66. If there are still flaws in the simulated image or in other design quality measures, the design iteration through stage is repeated, as described above. Alternatively, in case of serious flaws, the camera specification may be modified, and the process may return to step 50. Otherwise, system 30 outputs the final optical design and DCF tables, together with other aspects of the hardware circuit implementation of the camera (such as a netlist of engine 26), and the design process is thus complete.
  • the DCF tables may be tested and modified in a testbench calibration procedure. Such a procedure may be desirable in order to correct the DCF for deviations between the actual performance of the optics and the simulated performance that was used in the design process of Fig. 4. A calibration procedure that may be used for this purpose is described in the above- mentioned provisional application.
  • optic aberrations may be quantified in terms of coefficients of the Zernike polynomials.
  • the Zernike polynomials are described, for example, by Born & Wolf in Principles of Optics, 7 th edition (Cambridge University Press, 2002), in section 9.2, pages 523-527, and in Appendix VII, pages 905-910, which are incorporated herein by reference. Table I below lists different Zeraike polynomials:
  • ZA and Z9 contribute to spherical aberrations.
  • Z16, Z25, Z36, and Z37 also contribute to spherical aberrations.
  • Zl, Z2, and Z3 do not contribute to loss of resolution of the imaging system.
  • Fig. 8 is a schematic pictorial representation of camera 20 imaging a large depth of field, in accordance with an embodiment of the present invention.
  • the iterative detailed design process described above with reference to Fig. 4, and the design procedures described below with reference to Fig. 10, are followed.
  • the specifications of Tables II, III, IV, V, and VI are exemplary, and it will be understood that embodiments of the present invention provide numerous camera assemblies having a large depth of field.
  • step 50 (Fig. 4), is assumed to include cylindrically symmetric optics 22 having an F-number of 2.8.
  • the optics are to be produced using spherical or non-logarithmic aspherical surfaces.
  • the specification also includes a depth of field (DOF), for objects in a scene 100, from approximately 10 cm to infinity.
  • DOF is the range of object distances from camera 20 for which the image generated by the camera has an extent of an equivalent PSF, i.e., a PSF of an image produced by ISP 28 (Fig. 1), that is less than or equal to approximately 5 ⁇ m.
  • a more complete specification provided to step 50 is given in Table II below.
  • the depth of field of a camera is related to the F-number of the camera lens so that, for example, an F/2.8 lens has a largest depth of field when focused at the hyperfocal distance of the lens.
  • a depth of field of a 4.3 mm focal length lens is from infinity to of the order of 50 cm from the camera.
  • Fig. 9 shows detail of sensor 24, in accordance with an embodiment of the present invention.
  • Sensor 24 comprises an array of 1200 x 1600 pixels, which deconvolution engine 26 processes as 192 sub-arrays of 100 x 100 pixels.
  • deconvolution engine 26 applies a different 15 x 15 pixel DCF kernel (Figs 5A, 5B, 5C) to each sub-array.
  • the same DCF kernel may be applied to groups of four or eight sub-arrays of sensor 24.
  • the four sub- arrays 120 are all symmetric with respect to the center of sensor 24, and so engine 26 may apply the same DCF kernel to the four.
  • each DCF kernel applies to a group of sub-arrays, and each sub-array in a given group subtends the same angle to the optic axis of the camera.
  • the DCF kernels are iteratively created in step 62, typically using iterations of the optical design as described herein.
  • sensor 24 is considered as 768 sub-arrays of 50 x 50 pixels.
  • deconvolution engine applies 768 respective DCF kernels to the sub-arrays.
  • An example of a DCF kernel for the alternative embodiment, for a corner 50 x 50 array, is given in the matrix below.
  • the matrix above corresponds to a sub-array which subtends 33 degrees to the optic axis of the camera.
  • DCF kernels for sub-arrays subtending other angles to the optic axis by the iterative methods described herein.
  • application of the matrix above over other sub-arrays of sensor 24 yielded satisfactory results.
  • the initial optical design uses a virtual or dummy phase element, as described above with reference to step 52, and as exemplified below with reference to the flowchart of Fig. 10.
  • the virtual phase element is formed to have Zeraike coefficients with the same absolute value as those of Table III, but with opposite sign.
  • Fig. 10 is a flowchart describing procedures performed in steps 52 and 53 (Fig. 4), and Figs. HA, HB, ... , HE, show results generated in the procedures, in accordance with an embodiment of the present invention.
  • a modeling procedure 132 to evaluate the feasibility of the optical design, the optical designer may initially generate results for an ideal lens with a focal length corresponding the specification of Table II, together with a phase element having opposite Zernike coefficients to those of Table III.
  • Table IV the specification is assumed that the specification is as given below for Table IV, with Zernike coefficients as given below in Table V.
  • Fig. HA shows graphs of optical path differences (OPDs) for an ideal lens of focal length 3.3 mm combined with a phase element, in accordance with an embodiment of the present invention.
  • OPD graphs plot the number of wavefronts (W) vs. PX, PY ( ⁇ m) at wavelengths 486 nm and 656 nm. Inspection of the graphs shows that they are relatively smooth, so that it is likely that a final optical design corresponding to the specifications of Tables IV and V will be successful. In the event that the graphs are not smooth, the final design may not be successful, so that a change in the specifications may be advisable.
  • the optical designer in a classical design procedure 134, the optical designer generates a best version of the lens assembly according to the initial specifications, but ignoring the aberration specifications.
  • the number of lenses may be varied, although typically the fewer the lenses the lower the cost.
  • the parameters associated with the lens manufacture such as a maximum angle made by a surface of the lens relative to the optic axis, may be varied, although typically this may not be more than 55°. From these assemblies and with regard to the constraints, the designer selects a best version, and simulates results produced by the version.
  • the results typically include MTF curves, spot diagrams, and Zernike polynomial coefficients.
  • the MTF values are typically high, corresponding to the overall good image quality generated for a best version, whereas the spot diagrams and/or coefficient distributions may be poor.
  • Typical MTF results of the classical design procedure are shown in Fig. 1 IB.
  • a modified classical procedure 136 the designer attempts to improve the spot diagrams and/or coefficient distributions, typically at the cost of reducing the MTF values.
  • the improvement is usually achieved by applying merit functions to the Zernike coefficients, to attempt to reduce the coefficients to zero over the complete field of view of the assembly.
  • merit functions is described above. Comparison of the results shown in Fig. HC with those of Fig. 11B shows that zeroing the Zernike coefficients in procedure 136 reduces the MTF values.
  • a Zernike design procedure 138 a phase element with Zernike coefficients having the same values as those of Table V is applied to the version produced in procedure 136.
  • a number of designs are produced incorporating the phase element, and each of the designs is adjusted to attempt to zero the final Zernike coefficients. Zeroing the final Zernike coefficients normally reduces the MTF of the design.
  • a best design is chosen from the designs found in this stage.
  • Fig. HD shows examples of different MTF curves that may be produced by procedure 138.
  • a refinement procedure 140 the phase element is removed from the best design chosen in procedure 138.
  • the resulting design gives results shown in Fig. HE and, as expected, the Zernike coefficient values are close to those of the initial request.
  • the resulting design may be further refined to attempt to bring the coefficients closer to the initial request.
  • Fig. HE shows examples of MTF curves and values of Zernike coefficients before and after the refinement procedure.
  • the design is adjusted so that the MTF curves are all approximately the same shape. The conditions on the MTF curves, that the design meets, are described in more detail with respect to Figs. 12A-12G below.
  • the inventors have found that application of the procedures given in Fig. 10, together with the steps described in Fig. 4, gives a camera having a depth of field from approximately 10 cm to infinity.
  • the resulting value of Z4 is typically between approximately 0.6 and approximately 0.9, but the inventors have found that Z4 > approximately 0.1 gives good results, and that deconvolution engine 26 may be configured to satisfactorily correct the aberrations produced corresponding to these values of Z4, and to other non-zero values of the Zernike coefficients.
  • the coefficients of Z9, the third order spherical aberration are between approximately 0 and approximately -0.2
  • the coefficients of Zl 6, the fifth order spherical aberration are between approximately -0.1 and approximately -0.3
  • the coefficients of Z25, the seventh order spherical aberration are between approximately -0.4 and approximately -0.6
  • the coefficients of Z36, the ninth order spherical aberration are between approximately 0 and approximately 0.2.
  • the aberrations give an extent of the PSF that is typically greater than or equal to approximately 5 pixels, i.e., 14 ⁇ m, on sensor 24.
  • Figs. 12A-12G shows MTF curves for an exemplary optical imaging assembly 152, in accordance with an embodiment of the present invention.
  • Exemplary assembly 152 is described in more detail with reference to Fig. 13 below, and is one example of assembly 22.
  • Assembly 152 is designed according to the steps and procedures of Figs. 4 and 10.
  • the MTF curves may be produced during, and/or at the conclusion of, the design.
  • Fig. 12A is for an object 100 mm from camera 20;
  • Fig. 12B is for an object 150 mm from the camera;
  • Fig. 12C is for an object 200 mm from the camera;
  • Fig. 12D is for an object 300 mm from the camera;
  • Fig. 12E is for an object 500 mm from the camera;
  • Fig. 12F is for an object 1000 mm from the camera; and Fig. 12G is for an object at infinity.
  • Each figure comprises four fields of view 0°, 12°, 25°, and 33°.
  • the sagittal (S) and tangential (T) orientation curves are shown for each field of view.
  • Fig. 13 is a schematic diagram of optical imaging assembly 152, in accordance with an embodiment of the present invention.
  • Assembly 152 comprises five components: an aperture 151, three lenses 154, 156, 158 and an infra-red filter 160.
  • Assembly 152 forms its image on a front surface 162 of sensor 24, so that the front surface acts as a focal plane of the assembly.
  • Each lens has two aspheric surfaces, each surface being defined by expression (6):
  • z is a surface sag
  • r is a radius measured to an optic axis
  • c are curvature and conic constants for the surface
  • ⁇ X2, ... ⁇ g are parameters for the surface.
  • Tables VI, VII, and VTH below show dimensions in mm and parameters of each of the ten surfaces of the assembly. In the tables, the surfaces of each component are distinguished by a suffix A or B. Thus, aperture 151 has surfaces 15 IA and 15 IB.
  • the thickness gives the distance, measured along the optic axis, to the next surface, and thus gives spacings between the components.
  • aperture 151 has a thickness of 5.0-10 ⁇ 2 mm, the distance from surface 151B to the following surface, 154A, is l.l-10"2, and the distance from the final surface of the assembly, surface 160B of filter 160, to the sensor, is 8.61- 10"! J 31111-
  • Plastic E48R is produced by Zeon Corporation of Tokyo, Japan and has a refractive index approximately equal to 1.53.
  • Material OKP4 is produced by Osaka Gas Chemicals Ltd., of Japan and has a refractive index approximately equal to 1.61.
  • Glass BK7 is produced by Schott AG of Mainz, Germany and has a refractive index approximately equal to 1.52.
  • Table IX gives Zernike fringe coefficient values for assembly 152.
  • Table IX Inspection of Table IX shows that for all fields of view the values of Z4 are between approximately 0.6 and approximately 0.9, and the values of Z16 are between approximately -0.2 and approximately -0.3. Furthermore, the values of Zl, Z4, Z9, Z16, Z25, Z36, and Z37 in Table IX are all close to the target values in Table III.
  • DCF deconvolution filter

Abstract

An optical imaging assembly (22) having cylindrical symmetry, comprising a plurality of lenses having surfaces with curvatures and spacings between the surfaces, such that an optical image formed by the plurality of lenses has a defocus aberration coefficient greater than 0.1 at a focal plane of the assembly.

Description

OPTICS FOR AN EXTENDED DEPTH QF FIELD
CROSS-REFERENCE TO RELATED APPLICATIONS
This application claims the benefit of U.S. Provisional Patent Application 60/735,441, filed November 10, 2005. This application is a continuation-in-part of U.S. Patent Application 10/541,967, which was filed July 8, 2005, in the national phase of PCT Patent Application PCT/IL2004/000040, claiming the benefit of U.S. Provisional Patent Application 60/440,561, filed January 16, 2003, and is also a continuation-in-part of U.S. Patent Application 11/278,255, filed 31 March, 2006. All of these related applications are incorporated herein by reference. FIELD OF THE INVENTION
The present invention relates generally to digital imaging, and specifically to optics and methods for designing digital cameras with an enhanced depth of field.
BACKGROUND OF THE INVENTION
The objective optics used in digital cameras are typically designed so as to minimize the optical point spread function (PSF) and maximize the modulation transfer function (MTF), subject to the limitations of size, cost, aperture size, and other factors imposed by the camera manufacturer. The PSF of the resulting optical system may still vary from the ideal due to focal variations and aberrations. A number of methods are known in the art for measuring and compensating for such PSF deviations by digital image processing. For example, U.S. Patent 6,154,574, whose disclosure is incorporated herein by reference, describes a method for digitally focusing an out-of-focus image in an image processing system. A mean step response is obtained by dividing a defocused image into sub-images, and calculating step responses with respect to the edge direction in each sub-image. The mean step response is used in calculating PSF coefficients, which are applied in turn to determine an image restoration transfer function. An in-focus image is obtained by multiplying this function by the out-of-focus image in the frequency domain.
As another example, U.S. Patent 6,567,570, whose disclosure is incorporated herein by reference, describes an image scanner, which uses targets within the scanner to make internal measurements of the PSF. These measurements are used in computing convolution kernels, which are applied to images captured by the scanner in order to partially compensate for imperfections of the scanner lens system.
It is also possible to add a special-purpose blur to an image so as to create invariance to certain optical aberrations. Signal processing is then used to remove the blur. A technique of this sort is described by Kubala et al., in "Reducing Complexity in Computational Imaging Systems," Optics Express 11 (2003), pages 2102-2108, which is incorporated herein by reference. The authors refer to this technique as "Wavefront Coding." A special aspheric optical element is used to create the blur in the image. This optical element may be a separate stand-alone element, or it may be integrated into one or more of the lenses in the optical system. Optical designs and methods of image processing based on Wavefront Coding of this sort are described, for example, in U.S. Patent 5,748,371 and in U.S. Patent Application Publications US 2002/0118457 Al, US 2003/0057353 Al and US 2003/0169944 Al, whose disclosures are incorporated herein by reference.
SUMMARY OF THE INVENTION
In embodiments of the present invention, an optical imaging assembly that may be used in a digital camera is constructed so as to generate a distorted image at an imaging plane of the assembly. The digital camera incorporates a sensor at the imaging plane, and the distorted image formed on the sensor is corrected by a deconvolution engine so as to produce an undistorted image. The optical imaging assembly is configured to produce a high defocus aberration coefficient, greater than approximately 0.1, across the distorted image when the image is optimally focused to a focal plane of the assembly.
The high defocus aberration coefficient causes a modulation transfer function (MTF) for the assembly to have generally equal low values for all objects in a large field, typically from infinity to approximately 10 cm from the assembly. Because of the equal values of the MTF at the different object distances within the field, the deconvolution engine may be configured to improve the MTF of all the images at the different object distances. The final images thus produced by the assembly and the deconvolution engine are substantially free of aberrations for all objects within the field, so that the digital camera has a depth of field which is significantly higher than depths of field of prior art cameras. In some embodiments the imaging assembly is designed by assuming that a virtual phase plate is added to the assembly, the virtual phase plate having characteristics that generate aberration coefficients which are the negatives of those required by the deconvolution engine. The assembly with the phase plate is configured to give an undistorted image.
There is therefore provided, according to an embodiment of the present invention, an optical imaging assembly having cylindrical symmetry, including a plurality of lenses having surfaces with curvatures and spacings between the surfaces, such that an optical image formed by the plurality of lenses has a defocus aberration coefficient greater than 0.1 at a focal plane of the assembly. The defocus aberration coefficient may be greater than or equal to 0.6, and less than or equal to 0.9. A fifth order spherical aberration is typically between -0.2 and -0.3, a seventh order spherical aberration may be -0.5, and a ninth order spherical abeiτation may be 0.
In some embodiments, the curvatures are non-logarithmic. Typically, the optical image is of one or more objects within a field of view subtending up to 33° to an optic axis of the assembly.
In an embodiment the optical image is of one or more objects which are located at distances greater than or equal to 10 cm from the assembly. Typically, the assembly is configured so that a modulation transfer function generated by the one or more objects at the focal plane is less than or equal to 0.35 for spatial frequencies greater than 50 cycles per mm.
In a disclosed embodiment the optical image has an extent of a point spread function equal to or greater than 14 μm at the focal plane.
There is further provided, according to an embodiment of the present invention, a method for optical imaging, including: formulating an optical specification of an optical imaging assembly including optical elements having surfaces, the optical specification including one or more aberration coefficients; generating a virtual phase element having phase coefficients which are negatives of the aberration coefficients; generating a design of the optical assembly including the optical elements and the virtual phase element so as to determine curvatures and spacings of the surfaces that will focus an image of an object on a focal plane of the assembly in the presence of the virtual phase element; and fabricating the optical assembly in accordance with the curvatures and spacings determined by the design. In a disclosed embodiment, the method includes generating a digital filter to compensate for the one or more aberration coefficients, and coupling the digital filter to the optical assembly so that a filtered image formed by the digital filter from the image on the focal plane is free of aberrations. Typically, the filtered image is free of aberrations for the object being located at a distance greater than 10 cm from the optical assembly. There is further provided, according to an embodiment of the present invention, an optical imaging assembly having cylindrical symmetry, including a plurality of lenses having surfaces with curvatures and spacings between the surfaces, such that an optical image formed by the plurality of lenses at a focal plane of the assembly has a modulation transfer function (MTF) for spatial frequencies greater than 50 cycles per millimeter that is less than or equal to 0.35.
There is further provided, according to an embodiment of the present invention, a method for forming an optical imaging assembly, including: providing a plurality of lenses having cylindrical symmetry and surfaces with curvatures and spacings between the surfaces; and arranging the lenses so that an optical image formed by the plurality of lenses has a defocus aberration coefficient greater than 0.1 at a focal plane of the assembly.
There is further provided, according to an embodiment of the present invention, apparatus for optical imaging, including: an optical imaging assembly including optical elements having surfaces; and a processor which is configured to: formulate an optical specification of the optical imaging assembly, the optical specification including one or more aberration coefficients, generate a virtual phase element having phase coefficients which are negatives of the aberration coefficients; generate a design of the optical assembly including the optical elements and the virtual phase element so as to determine curvatures and spacings of the surfaces that will focus an image of an object on a focal plane of the assembly in the presence of the virtual phase element; and fabricate the optical assembly in accordance with the curvatures and spacings determined by the design.
There is further provided, according to an embodiment of the present invention, a method for forming an optical imaging assembly, including: providing a plurality of lenses having cylindrical symmetry and surfaces with curvatures and spacings between the surfaces; and arranging the lenses so that an optical image formed by the plurality of lenses at a focal plane of the assembly has a modulation transfer function (MTF) for spatial frequencies greater than 50 cycles per millimeter that is less than or equal to 0.35. There is further provided, according to an embodiment of the present invention, an optical imaging assembly, including: an aperture having a semi-diameter of 7.5-10"! mm; and a plurality of lenses having surfaces with curvatures and spacings between the surfaces, wherein respective surfaces of the plurality of lenses are defined by an equation
Figure imgf000006_0001
where z is a sag of the surface, r is a radius measured to an optic axis of the assembly, c, k are respective curvature and conic constants for the surface, and α2, ... αg are parameters for the surface, the lenses including: a first aspheric lens following the aperture, having a refractive index 1.53, a first surface wherein c = 9.57-10"2, k = 9.74-101, α2 = -1.38-10"2, 0,3 = -9.0-10"1, 014 = 2.83,
015 = -2.26, αg = -2.96, αγ = 2.69, a$ = 1.61, and a second surface wherein c = -4.46- 10"*, k = 0, α2 = -9.83-10-2, α3 = 1.24-10"1, α4 = -4.25 10"1, α5 = 7.62-10-1, α6 = -8.04-10"1, α7 = 4.73-10-1, α8 = -1.28-10-1; a second aspheric lens following the first aspheric lens, having a refractive index 1.61, a third surface wherein c = -7.60-1O-I, k = -3.37, and absolute values of (X2, 013, 04, C15, αg, aη, αg less than 1, and a fourth surface wherein c = -2.42-1O-I, k = 2.38, and absolute values of α2, α3, 04, α5, αg, aq, αg less than 1; and a third aspheric lens following the second aspheric lens, having a refractive index .
1.53, a fifth surface wherein c = 8.1-1O-I, ^ = -4.03, and absolute values of α2, α3, α4, 0.5, αg, aη, αg less than 1, and a sixth surface wherein c = 5.17-10"!, k = -4.07, and absolute values of ot2, 013, (X4, CC5, αg, αy, αg less than 1.
The present invention will be more fully understood from the following detailed description of the embodiments thereof, taken together with the drawings in which: BRIEF DESCRIPTION OF THE DRAWINGS
Fig. 1 is a block diagram that schematically illustrates a digital camera, in accordance with an embodiment of the present invention;
Fig. 2 is a schematic, pictorial illustration of a system for designing a digital camera, in accordance with an embodiment of the present invention; Fig. 3 A is a schematic, pictorial illustration showing conceptual elements of a digital camera used in a design process, in accordance with an embodiment of the present invention;
Fig. 3B is a plot of modulation transfer functions (MTF) for a digital camera with and without application of a deconvolution filter, in accordance with an embodiment of the present invention;
Fig. 4 is a flow chart that schematically illustrates a method for designing a digital camera, in accordance with an embodiment of the present invention;
Figs. 5A-5C are schematic, perspective plots of DCF kernels used in a digital camera, in accordance with an embodiment of the present invention; Fig. 6 is an image that simulates the output of an image sensor using objective optics with specifications that have been relaxed in accordance with an embodiment of the present invention;
Fig. 7 is an image that simulates the effect of application of an appropriate DCF to the image of Fig. 6, in accordance with an embodiment of the present invention; Fig. 8 is a schematic pictorial representation of a camera imaging a large depth of field, in accordance with an embodiment of the present invention;
Fig. 9 schematically shows detail of a camera sensor, in accordance with an embodiment of the present invention;
Fig. 10 is a flowchart schematically describing procedures performed in designing a digital camera, in accordance with an embodiment of the present invention;
Figs. 1 IA-I IE schematically show results generated in the procedures of Fig. 10, in accordance with an embodiment of the present invention; Figs. 12A-12G show Modulation Transfer Function curves for an exemplary optical imaging assembly, in accordance with an embodiment of the present invention; and
Fig. 13 is a schematic diagram of the exemplary optical imaging assembly, in accordance with an embodiment of the present invention.
DETAILED DESCRIPTION OF EMBODIMENTS
DEFINITIONS
The following is a non-exhaustive list of technical terms that are used in the present patent application and in the claims. Although these terms are used herein in accordance with the plain meaning accorded the terms in the art, they are listed below for the convenience of the reader in understanding the following description and the claims.
• Pitch of a detector array refers to the center-to-center distance between elements of the array. • Cylindrical symmetry describes a structure, such as a simple or compound lens, which has an optical axis such that the structure is invariant under rotation about the optical axis for any and all angles of rotation.
• Point spread function (PSF) is the impulse response of an optical system in the spatial domain, i.e., the image formed by the system of a bright point object against a dark background.
• Extent of the PSF is the full width at half maximum (FWHM) of the PSF.
• Optical transfer function (OTF) is the two-dimensional Fourier transform of the PSF to the frequency domain. Because of the ease with which a PSF may be transformed into an OTF, and vice versa, computation of the OTF is considered to be equivalent to computation of the PSF for the purposes of the present invention.
• Modulation transfer function (MTF) is the modulus of the OTF.
• Optical radiation refers to electromagnetic radiation in any of the visible, infrared and ultraviolet regions of the spectrum.
SYSTEM OVERVIEW Fig. 1 is a block diagram that schematically illustrates a digital camera 20, in accordance with an embodiment of the present invention. The camera comprises refractive objective optics 22, which focus an image onto an image sensor 24, so that a surface of the sensor acts as a focal plane of the optics. Optics 22 are also referred to herein as optical imaging assembly 22. Optics 22 are designed in an iterative process together with a deconvolution engine 26 that operates on image data that are output by image sensor 24. The deconvolution engine applies one or more digital filters, typically comprising at least one deconvolution filter (DCF), to the image data. The design process and method of filtering are described in detail hereinbelow. The DCF kernel is typically chosen so as to correct for blur in the image formed by optics 22. After filtering, the image data are processed by an image signal processor (ISP) 28, which performs standard functions such as color balance and format conversion and outputs the resulting image.
The optical and digital processing schemes illustrated in Fig. 1 are shown here solely for the sake of example, as an aid to understanding the techniques and tools that are described hereinbelow. In practice, the principles of the present invention may be applied in conjunction with a wide variety of electronic imaging systems, using substantially any sort of optical design and substantially any type of image sensor, including both two- dimensional detector matrices and linear detector arrays, as are known in the art. Deconvolution engine 26 and ISP 28 may be implemented as separate devices or as a single integrated circuit component. In either case, the deconvolution engine and ISP are typically combined with other I/O and processing elements, as are known in the art. In the context of the present patent application, the term "digital camera" should therefore be understood as referring to any and all sorts of electronic imaging systems that comprise an image sensor, objective optics for focusing optical radiation onto the image sensor, and electronic circuits for processing the sensor output.
Fig. 2 is a schematic, pictorial illustration showing a system 30 for designing a digital camera, in accordance with an embodiment of the present invention. System comprises a digital processing design station 32 and an optical design station 34. Processing design station 32 receives a camera specification as input, from a camera manufacturer, for example, specifying the key dimensions, sensor type and desired optical characteristics (referred to hereinafter as the target optical specification) of the camera. The specified optical characteristics may include, for example, the number of optical elements, materials, tolerances, focal length, magnification, aperture (F-number), depth of field, and resolution performance. The optical resolution performance is typically defined in. terms of the MTF, but it may alternatively be specified in terms of PSF, wavefront quality, aberrations, and/or other measures of optical and image quality that are known in the art.
Processing design station 32 analyzes and modifies the target optical specification, taking into account the expected operation of engine 26, in order to provide a modified optical specification to the optical design station. Typically, both the original camera specification and the modified optical specification use cylindrically-symmetrical optical elements. Specialized phase plates or other elements that break the cylindrical symmetry of the optics are generally undesirable, due to their added cost, and engine 26 is able to correct the aberrations of optics 22 without requiring the use of such elements. In addition, processing design station 32 may compute and provide to optical design station 34 a merit function, indicating target values of the aberrations of optics 22 or scoring coefficients to be used in weighting the aberrations in the course of optimizing the optical design. The aberrations express deviations of the optical wavefront created by optics 22 from the ideal, and may be expressed, for example, in terms of Zernike polynomials or any other convenient mathematical representation of the wavefront that is known in the art.
Optical design station 34 is typically operated by a lens designer, in order to produce a lens design according to the modified optical specification provided by processing design station 32. The processing design station determines the optimal DCF (and possibly other filters) to be used in engine 26 in conjunction with this lens design. The DCF computation is tied to the specific lens design in question so that the filter coefficients reflect the "true" PSF of the actual optical system with which the DCF is to be used.
The processing design station then evaluates the optical design together with the DCF in order to assess the combined result of the expected optical quality of optics 22 and the enhancement expected from engine 26, and to compare the result to the target optical specification. The assessment may take the form of mathematical analysis, resulting in a quality score. A quality scoring schemes that may be used in this context is described hereinbelow. Alternatively, other quality scoring schemes may be used, such as that described, for example, in the above-mentioned PCT publication WO 2004/063989 A2. Alternatively or additionally, station 32 may generate and display a simulated image 36, which visually demonstrates the output image to be expected from the camera under design based on the current choice of optical specifications and DCF. If the result of the analysis by station 32 indicates that the combined optical and DCF design will meet the target specifications, then the complete camera design, including optics and DCF, is output for production. Otherwise, the processing design station may perform further design iterations internally, or it may generate a further modified optical specification, which it passes to optical design station 34 for generation of a modified optical design. This process may continue iteratively until a suitable optical design and DCF are found. Details of this process are described hereinbelow with reference to Fig. 4.
Typically, stations 32 and 34 comprise general-purpose computers running suitable software to carry out the functions described herein. The software may be downloaded to the computers in electronic form, over a network, for example, or it may alternatively be furnished on tangible media, such as optical, magnetic, or electronic memory media. Alternatively, some of the functions of stations 32 and/or 34 may be implemented using dedicated or programmable hardware components. The functions of optical design station 34 may be carried out using off-shelf optical design software, such as ZEMAX® (produced by ZEMAX Development Corp., San Diego, California). Although stations 32 and 34 are shown and described, for the sake of conceptual clarity, as separate computer workstations, the functions of these stations may alternatively be combined in a single physical machine, running software processes for both optical design and digital processing design.
Fig. 3A is a schematic, pictorial illustration showing conceptual elements of camera 20, as they are applied in the design process used in system 30, in accordance with an embodiment of the present invention. System 30 takes engine 26 into account in the design of optics 22, as explained hereinabove, and thus relates to the DCF as a sort of "virtual lens" 40. In other words, the design constraints on the actual objective optics are relaxed by the use of this virtual lens, as though the optical designer had an additional optical element to incorporate in the design for purposes of aberration correction. The virtual lens that is implemented in engine 26 is chosen, in conjunction with the actual optical lenses, to give an image output that meets the manufacturer's camera specifications.
Fig. 3B is a plot showing the MTF of a camera designed using system 30, in accordance with an embodiment of the present invention. The plot includes an uncorrected curve 44, corresponding to the modified optical specification generated by station 32 for use in designing optics 22 on station 34. The low MTF permitted by curve 44 is indicative of the expected improvement in MTF that can be achieved by use of DCF 26. A corrected curve 46 shows the net MTF of the camera that is achieved by applying the DCF to the image sensor output. These curves show the MTF at the center of the optical field, with the object at a certain distance from the camera. In practice, the MTF may be specified at multiple different focal depths and field angles.
The design concept exemplified by Figs. 3 A and 3 B permits the camera manufacturer to achieve the desired level of optical performance with fewer, smaller and/or simpler optical components than would be required to achieve the same result by optical means alone. Additionally or alternatively, the camera may be designed for enhanced performance, such as reduced aberrations, reduced F-number, wide angle, macro operation, or increased depth of field.
DETAILED DESIGN PROCESS
Fig. 4 is a flow chart that schematically illustrates a method for designing a digital camera, in accordance with an embodiment of the present invention. The method will be described hereinbelow, for the sake of clarity, with reference to camera 20 and system 30, although the principles of this method may be applied generally to other cameras and using other design systems.
The point of departure of the design is the camera specification, as noted above. Processing design station 32 translates the target optical specification of the camera into a modified optical specification, at a specification translation step 50. For this purpose, station 32 uses an estimate of the DCF to be implemented in the camera. The image enhancement to be expected due to this DCF is then applied to the optical specification in order to estimate how far the optical design parameters, such as the MTF, can be relaxed. Image enhancement by the DCF, however, tends to amplify noise in the output of image sensor 24. Generally speaking, the noise gain NG is proportional to the norm of the DCF (-vDfcD , wherein D is the DCF kernel and the superscript t indicates the Heπnitian transpose). Therefore, in estimating the DCF, and hence in estimating the degree to which the optical design parameters can be relaxed, the processing design station uses the maximum permissible noise gain as a limiting condition. Typically, engine 26 may also comprise a noise filter. The limit placed on the DCF coefficients by the noise gain may thus be mitigated by the noise reduction that is expected due to the noise filter. In other words, the norm of the DCF kernel is approximately given by the product of the maximum permissible noise gain with the expected noise reduction factor (i.e., the ratio of image noise following the noise filter to image noise without noise filtering). Alternatively, a more accurate estimate of the overall noise gain may be obtained by taking the norm of the product of the noise filter multiplied by the DCF in the frequency domain.
In order to determine the noise gain and permissible MTF reduction, the OTF may be assumed, at first approximation, to be linear as a function of spatial frequency q, which is normalized to the Nyquist frequency of image sensor 24:
OTF = 1 - λq } q ≤ 1 / λ
OTF = 0 } q > 1 I λ
The PSF may be determined analytically from the OTF of equation (1). Because of the zeroes in the OTF, the frequency-domain representation of the DCF to be used in the camera may be estimated as:
OTF
DCF = (2)
OTF^ + (f
wherein a is a small number that keeps the DCF from exploding for small PSF. The noise gain NG due to the DCF of equation (2) depends on the two parameters λ, a :
, z s, arctan (l/a) f l
+ 1 (3) λ a \ a
These parameters are chosen so that the noise gain does not exceed a target bound, for example, 300%. If the original camera specifications include a noise figure, the maximal permissible noise gain may be determined by comparing the expected noise characteristic of image sensor 24 to the noise specification. As noted above, digital smoothing of the noise in the output image may also be taken into account in order to permit the constraint on noise gain in the DCF to be relaxed.
Various noise removal methods, as are known in the art, may be used in engine 26. For example, a morphological operation may be used to identify edges in the image, followed by low-pass filtering of non-edge areas. The choice of noise removal method to be used in engine 26, however, is beyond the scope of the present invention.
Having chosen appropriate values of the parameters, the average MTF over the normalized frequency range [0,1] is given by:
MTFavg = — (1 - α * arctan (1 / α) ) (4)
The formulas given above in equations (3) and (4) apply for λ > 1 , which will be the case in most simple camera designs. Alternative estimates may be developed for high- resolution cameras in which λ < 1. For α «\, the noise gain may be expressed as a polynomial series in α or in the form:
Figure imgf000014_0001
Other representations will be apparent to those skilled in the art. Equations (4) and (5) may be used in estimating how far the MTF of optics 22 may be reduced relative to the original target specification, subject to a given noise gain limit. This reduction factor may be applied, for example, to the MTF required by the original camera specification at a benchmark frequency, such as half the Nyquist frequency. In the example shown in Fig. 3B, the target MTF has been reduced to about 1/3 of its original specified value. The MTF will be restored in the output image from camera 20 by operation of the DCF in engine 26.
Referring back now to Fig. 4, processing design station 32 may also generate a merit function at step 50 for use by the optical designer. The merit function may take the form of aberration scores, which are assigned to each significant aberration that may characterize optics 22. For this purpose, the aberrations may be expressed, for example, in terms of the Zernike polynomials, for each of the colors red, green and blue individually. Standard software packages for optical design, such as ZEMAX, are capable of computing the Zernike polynomial coefficients for substantially any design that they generate. Values of the merit functions may be provided in tabular form. Generation of these values is described in detail in the above-mentioned PCT Publication WO 2004/063989 A2.
Alternatively or additionally, processing design station 32 may generate target wavefront characteristics that the optical design should achieve in the image plane (i.e., the plane of sensor 24). These wavefront characteristics may conveniently be expressed in terms of values of the aberrations of optics 22, such as Zernike coefficient values. Typically, aberrations that can be corrected satisfactorily by deconvolution engine 26 may have high values in the optical design, whereas aberrations that are difficult to correct should have low values. In other words, an aberration that would have a high score in the merit function will have a low target value, and vice versa. The target aberration values can be seen as the inverse of the wavefront corrections that can be achieved by "virtual lens" 40. The target aberration values may also include aberrations that reduce the sensitivity of the optics to various undesirable parameters, such as manufacturing deviations and defocus.
An optical designer working on station 34 uses the specification, along with the merit function and/or aberration target values provided at step 50, in generating an initial design of optics 22, at an optical design step 52. The designer may use the merit function in determining a design score, which indicates how to trade off one aberration against another in order to generate an initial design that maximizes the total of the merit scores subject to the optical specification. Additionally or alternatively, the optical designer may insert a dummy optical element, with fixed phase characteristics given by the target aberration values as an additional element in the optical design. This dummy optical element expresses the wavefront correction that is expected to be achieved using engine 26 and thus facilitates convergence of the calculations made by the optical design software on station 34 to the desired design of the elements of optics 22. Control of the design process now passes to processing design station 32, in a design optimization stage 53. The processing design station analyzes the optical design, at a design analysis step 54. The analysis at this step may include the effect of virtual lens
40. At step 54, station 32 typically computes the optical performance of the optics as a function of wavelength and of location in the image plane. For example, station 32 may perform an accurate ray trace computation based on the initial optical design in or to calculate a phase model at the image plane, which may be expressed in terms of Zernike polynomial coefficients. The total aberration - and hence the PSF - at any point in the image plane may be obtained from the total, wavefront aberration, which is calculated by summing the values of the Zernike polynomials.
Station 32 determines a design quality score, at a scoring step 55. Typically, this score combines the effects of the PSF on image resolution and on artifacts in the image, and reflects the ability of engine 26 to compensate for these effects. The score measures the extent to which the current optical design, taken together with filtering by engine 26, will satisfy the camera specification that was originally provided as input to station 32 as input to step 50.
In an exemplary embodiment, the score computed at step 55 is based on the camera specification and on a set of weights assigned to each parameter in the camera specification. The camera specification is expressed in a list of desired parameter values at various image plane locations and wavelengths, such as:
• MTF
• Geometrical distortion
• Field of view • Chromatic aberrations
• Chief ray angle
• F-number
• Relative illumination
• Artifact level • Glare
• Back focal length
• Manufacturing tolerances
• Depth of field
• Noise level • Total length of optics. The weight assigned to each parameter is typically determined by its scaling, subjective importance, and likelihood of satisfying the desired parameter value relative to other parameters.
The overall score is computed by summing the weighted contributions of all the relevant parameters. In this embodiment, if a given parameter is within the specified range, it makes no contribution to the score. If the value is outside the specified range, the score is decreased by the square difference between the parameter value and the closest permissible value within the specified range, multiplied by the appropriate weight. A design that fully complies with the camera specification will thus yield a zero score, while non-compliance will yield negative values. Alternatively, other parameters and other methods may be used in computing numerical values representing how well the current design satisfies the camera specification.
The score computed at step 55 is assessed to determine whether it indicates that the current design is acceptable, at a quantitative assessment step 56. If the design does not meet the specification, station 32 modifies the optical design parameters at an optimization step 58. For this purpose, the station may estimate the effects of small changes in the aberrations on the PSF. This operation gives a multi-dimensional gradient, which is used in computing a change to be made in the optical design parameters by linear approximation. The DCF parameters may be adjusted accordingly. A method for computing and using gradients of this sort is described, for example, in . the above- mentioned PCT Publication WO 2004/063989 A2. The results of step 58 are input to step 54 for recomputation of the optical performance analysis. The process continues iteratively through steps 55 and 56 until the design quality score reaches a satisfactory result. Once the design has converged, the design parameters are presented by processing design station 32 to the system operator, at a design checking step 60. Typically, the system operator reviews the optical design (as modified by station 32 in step 58, if necessary), along with the results of the design analysis performed at step 54. Additionally or alternatively, the optical design and DCF may be used at this point in generating a simulated output image, representing the expected performance of the camera in imaging a known scene or test pattern. (Exemplary simulated images of this sort are shown below in Figs. 6 and 7.) The system operator reviews the design in order to verify that the results are indeed satisfactory for use in manufacture of camera 20. If not, the operator may change certain parameters, such as specification parameters and/or scoring weights, and return to stage 53. Alternatively, if it appears that there are serious problems with the design, the operator may initiate changes to the original camera specification and return the process to step 50. This sort of operator involvement may also be called for if stage 53 fails to converge to an acceptable score at step 56.
Once the design is found to be acceptable, processing design station 32 generates tables of values to be used in camera 20, at a DCF creation step 62. Typically, because of the non-uniform performance of optics 22, the DCF tables vary according to location in the image plane. In an exemplary embodiment, a different DCF kernel is computed for each region of 50 x 50 pixels in image sensor 24.
Furthermore, when sensor 24 is a color image sensor, different kernels are computed for the different color planes of sensor 24. For example, referring back to Fig. 3 A, common mosaic image sensors may use a Bayer pattern of red, green, and blue pixels 42. In this case, the output of the image sensor is an interleaved stream of sub-images, comprising pixel samples belonging to different, respective colors. DCF 26 applies different kernels in alternation, so that the pixels of each color are filtered using values of other nearby pixels of the same color. Appropriate kernel arrangements for performing this sort of filtering are described in U.S. Provisional Patent Application 60/735,519, filed November 10, 2005, which is assigned to the assignee of the present patent application and is incorporated herein by reference.
Figs. 5A, 5B and 5C are schematic, perspective plots of DCF kernels 70, 72 and 74 for red, green, and blue pixels, respectively, which are computed in accordance with an embodiment of the present invention. Each kernel extends over 15 x 15 pixels, but contains non-zero values only at pixels of the appropriate color. In other words, in red kernel 70, for example, in each square of four pixels, only one - the red pixel - has a nonzero value. Blue kernel 74 is similarly constructed, while green kernel 72 contains two non-zero values in each four-pixel square, corresponding to the greater density of green pixels in the Bayer matrix. In each kernel, the central pixel has a large positive value, while surrounding values are lower and may include negative values 76. As explained above, the DCF values are chosen so that the norm does not exceed the permitted noise gain.
Referring back to Fig. 4, design station 32 uses the DCF tables from step 62 and the optical design output from stage 53 in simulating the performance of camera 20, at a simulation step 64. The simulation may also use characteristics, such as noise figures, of image sensor 24 that is to be installed in the camera, as well as other factors, such as manufacturing tolerances to be applied in producing the camera and/or operation of ISP 28. The results of this step may include simulated images, like image 36 (Fig. 2), which enable the system operator to visualize the expected camera performance.
Figs. 6 and 7 are images that simulate the expected output of camera 20, as may be generated at step 64, in accordance with an embodiment of the present invention. Fig. 6 shows a standard test pattern as it would be imaged by optics and captured by image sensor 24, without the use of DCF 26. The image of the test pattern is blurred, especially at higher spatial frequencies, due to the low MTF of camera 20. (The MTF is given roughly by uncorrected curve 44 in Fig. 3B.) In addition, the image pixels are decimated due to the use of a color mosaic sensor, and random noise is added to the image corresponding to the expected noise characteristics of the image sensor.
Fig. 7 shows the image of Fig. 6 after simulation of processing by DCF 26, including noise removal as described hereinbelow. The MTF of this image is given roughly by curve 46 in Fig. 3B. (The aliasing apparent in the images of the high- frequency test patterns is the result of a true simulation of the performance of a low- resolution image sensor following DCF processing.) The system operator, viewing this image, is able to ascertain visually whether the camera performance will meet the original camera specifications that were provided at step 50.
The system operator's visual assessment is combined with the numerical results of the design analysis, in order to determine whether the overall performance of the design is acceptable, at an acceptance step 66. If there are still flaws in the simulated image or in other design quality measures, the design iteration through stage is repeated, as described above. Alternatively, in case of serious flaws, the camera specification may be modified, and the process may return to step 50. Otherwise, system 30 outputs the final optical design and DCF tables, together with other aspects of the hardware circuit implementation of the camera (such as a netlist of engine 26), and the design process is thus complete. Optionally, after prototypes of optics 22 have been fabricated, the DCF tables may be tested and modified in a testbench calibration procedure. Such a procedure may be desirable in order to correct the DCF for deviations between the actual performance of the optics and the simulated performance that was used in the design process of Fig. 4. A calibration procedure that may be used for this purpose is described in the above- mentioned provisional application.
ZERMKE POLYNOMIALS
As stated above, optic aberrations may be quantified in terms of coefficients of the Zernike polynomials. The Zernike polynomials are described, for example, by Born & Wolf in Principles of Optics, 7th edition (Cambridge University Press, 2002), in section 9.2, pages 523-527, and in Appendix VII, pages 905-910, which are incorporated herein by reference. Table I below lists different Zeraike polynomials:
Figure imgf000020_0001
Table I
As shown in Table I, ZA and Z9 contribute to spherical aberrations. In addition, Z16, Z25, Z36, and Z37 also contribute to spherical aberrations. Zl, Z2, and Z3 do not contribute to loss of resolution of the imaging system. EXTENDED DEPTH OF FIELD
Fig. 8 is a schematic pictorial representation of camera 20 imaging a large depth of field, in accordance with an embodiment of the present invention. In the following description, unless otherwise stated, it is assumed that the iterative detailed design process described above with reference to Fig. 4, and the design procedures described below with reference to Fig. 10, are followed. Also in the following description, the specifications of Tables II, III, IV, V, and VI are exemplary, and it will be understood that embodiments of the present invention provide numerous camera assemblies having a large depth of field. A camera specification, corresponding to that provided to design station 32 (Fig.
1) at step 50 (Fig. 4), is assumed to include cylindrically symmetric optics 22 having an F-number of 2.8. The optics are to be produced using spherical or non-logarithmic aspherical surfaces. The specification also includes a depth of field (DOF), for objects in a scene 100, from approximately 10 cm to infinity. The DOF is the range of object distances from camera 20 for which the image generated by the camera has an extent of an equivalent PSF, i.e., a PSF of an image produced by ISP 28 (Fig. 1), that is less than or equal to approximately 5μm. A more complete specification provided to step 50 is given in Table II below.
Figure imgf000021_0001
Table II In prior art systems, the depth of field of a camera is related to the F-number of the camera lens so that, for example, an F/2.8 lens has a largest depth of field when focused at the hyperfocal distance of the lens. In this case a depth of field of a 4.3 mm focal length lens is from infinity to of the order of 50 cm from the camera.
In addition to the specification generated in step 50, Zernike polynomial coefficient values shown in Table III are provided for continuing to optical design step
52. The description hereinbelow assumes that the coefficient values are for Zernike fringe coefficients, Those skilled in the art will be able to convert the values for the fringe coefficients to Zernike standard coefficients, without undue experimentation.
Figure imgf000022_0001
Table III
Fig. 9 shows detail of sensor 24, in accordance with an embodiment of the present invention. Sensor 24 comprises an array of 1200 x 1600 pixels, which deconvolution engine 26 processes as 192 sub-arrays of 100 x 100 pixels. Because of the non-uniform performance of optics 22 over sensor 24, as described above, deconvolution engine 26 applies a different 15 x 15 pixel DCF kernel (Figs 5A, 5B, 5C) to each sub-array. However, since optics 22 are cylindrically symmetric, the same DCF kernel may be applied to groups of four or eight sub-arrays of sensor 24. For example, the four sub- arrays 120 are all symmetric with respect to the center of sensor 24, and so engine 26 may apply the same DCF kernel to the four. Similarly, the eight sub-arrays 122 are also symmetric with respect to the center, so that engine 26 may apply the same DCF kernel to the eight. Inspection of the sub-arrays of sensor 24 shows that 33 different DCF kernels suffice for engine 26 to process the complete sensor. Each DCF kernel applies to a group of sub-arrays, and each sub-array in a given group subtends the same angle to the optic axis of the camera.
The DCF kernels are iteratively created in step 62, typically using iterations of the optical design as described herein.
In an alternative embodiment of the present invention, sensor 24 is considered as 768 sub-arrays of 50 x 50 pixels. In the alternative embodiment, deconvolution engine applies 768 respective DCF kernels to the sub-arrays. An example of a DCF kernel for the alternative embodiment, for a corner 50 x 50 array, is given in the matrix below.
2 8 1 -1 2 14 1 10 1 19 2 11 0 6 0
0 -2 10 0 -13 1 23 0 23 -I 10 0 0 -1 -5
3 21 -2 -21 -3 -25 -1 -4 1 -3 0 -11 -1 2 1
0 2 31 -5 -23 -13 -151 -11 -165 -4 -49 0 33 0 -6
3 7 0 -15 -8 33 -34 -182 -29 -163 -3 -76 1 13 2
0 -4 -41 -5 -6 25 161 -38 330 -17 -125 3 31 -1 9
3 14 -7 -16 -8 -123 93 401 70 470 -12 -18 2 -8 3
0 0 14 -15 -199 -52 304 358 479 -37 -176 -1 -7 -1 15
1 26 0 14 -25 -171 62 136 -42 177 -22 -144 1 -9 3
0 -1 33 -1 -90 -26 -35 -19 -81 -29 -129 -3 -3 -2 17
0 12 0 3 -8 -55 -9 4 -20 -45 -3 -58 -5 18 4
0 -1 0 0 1 1 -7 0 9 1 -11 -3 -13 -4 13
0 7 -1 -1 0 3 1 -8 1 -8 1 7 -5 -15 3
0 -1 -2 0 6 0 9 1 13 1 8 -1 1 -1 10
0 0 1 0 2 0 3 0 3 0 2 0 1 0 3
The matrix above corresponds to a sub-array which subtends 33 degrees to the optic axis of the camera. Those having ordinary skill in the art will be able to produce
DCF kernels for sub-arrays subtending other angles to the optic axis, by the iterative methods described herein. In one embodiment, application of the matrix above over other sub-arrays of sensor 24 yielded satisfactory results.
In one embodiment of the present invention, the initial optical design uses a virtual or dummy phase element, as described above with reference to step 52, and as exemplified below with reference to the flowchart of Fig. 10. The virtual phase element is formed to have Zeraike coefficients with the same absolute value as those of Table III, but with opposite sign.
Fig. 10 is a flowchart describing procedures performed in steps 52 and 53 (Fig. 4), and Figs. HA, HB, ... , HE, show results generated in the procedures, in accordance with an embodiment of the present invention. In a modeling procedure 132, to evaluate the feasibility of the optical design, the optical designer may initially generate results for an ideal lens with a focal length corresponding the specification of Table II, together with a phase element having opposite Zernike coefficients to those of Table III. In the following description, it is assumed that the specification is as given below for Table IV, with Zernike coefficients as given below in Table V.
Figure imgf000024_0001
Table IV
Figure imgf000024_0002
Table V
An example of results generated in procedure 132 is given in Fig. 1 IA. Fig. HA shows graphs of optical path differences (OPDs) for an ideal lens of focal length 3.3 mm combined with a phase element, in accordance with an embodiment of the present invention. The OPD graphs plot the number of wavefronts (W) vs. PX, PY (μm) at wavelengths 486 nm and 656 nm. Inspection of the graphs shows that they are relatively smooth, so that it is likely that a final optical design corresponding to the specifications of Tables IV and V will be successful. In the event that the graphs are not smooth, the final design may not be successful, so that a change in the specifications may be advisable.
Returning to Fig. 10, in a classical design procedure 134, the optical designer generates a best version of the lens assembly according to the initial specifications, but ignoring the aberration specifications. It will be understood that for any given set of specifications, there are many possible assemblies that may be considered, and that there are also many practical constraints. For example, the number of lenses may be varied, although typically the fewer the lenses the lower the cost. Alternatively or additionally, the parameters associated with the lens manufacture, such as a maximum angle made by a surface of the lens relative to the optic axis, may be varied, although typically this may not be more than 55°. From these assemblies and with regard to the constraints, the designer selects a best version, and simulates results produced by the version. The results typically include MTF curves, spot diagrams, and Zernike polynomial coefficients. The MTF values are typically high, corresponding to the overall good image quality generated for a best version, whereas the spot diagrams and/or coefficient distributions may be poor. Typical MTF results of the classical design procedure are shown in Fig. 1 IB.
In a modified classical procedure 136, the designer attempts to improve the spot diagrams and/or coefficient distributions, typically at the cost of reducing the MTF values. The improvement is usually achieved by applying merit functions to the Zernike coefficients, to attempt to reduce the coefficients to zero over the complete field of view of the assembly. The application of merit functions is described above. Comparison of the results shown in Fig. HC with those of Fig. 11B shows that zeroing the Zernike coefficients in procedure 136 reduces the MTF values. In a Zernike design procedure 138, a phase element with Zernike coefficients having the same values as those of Table V is applied to the version produced in procedure 136. Typically, a number of designs are produced incorporating the phase element, and each of the designs is adjusted to attempt to zero the final Zernike coefficients. Zeroing the final Zernike coefficients normally reduces the MTF of the design. A best design is chosen from the designs found in this stage. Fig. HD shows examples of different MTF curves that may be produced by procedure 138.
In a refinement procedure 140, the phase element is removed from the best design chosen in procedure 138. The resulting design gives results shown in Fig. HE and, as expected, the Zernike coefficient values are close to those of the initial request. The resulting design may be further refined to attempt to bring the coefficients closer to the initial request. Fig. HE shows examples of MTF curves and values of Zernike coefficients before and after the refinement procedure. In both the design and the refinement procedures, in addition to bringing the Zernike coefficients close to the initial request, the design is adjusted so that the MTF curves are all approximately the same shape. The conditions on the MTF curves, that the design meets, are described in more detail with respect to Figs. 12A-12G below.
The inventors have found that application of the procedures given in Fig. 10, together with the steps described in Fig. 4, gives a camera having a depth of field from approximately 10 cm to infinity. The resulting value of Z4 is typically between approximately 0.6 and approximately 0.9, but the inventors have found that Z4 > approximately 0.1 gives good results, and that deconvolution engine 26 may be configured to satisfactorily correct the aberrations produced corresponding to these values of Z4, and to other non-zero values of the Zernike coefficients. In addition, typically the coefficients of Z9, the third order spherical aberration, are between approximately 0 and approximately -0.2, the coefficients of Zl 6, the fifth order spherical aberration, are between approximately -0.1 and approximately -0.3, the coefficients of Z25, the seventh order spherical aberration, are between approximately -0.4 and approximately -0.6, and the coefficients of Z36, the ninth order spherical aberration, are between approximately 0 and approximately 0.2. The aberrations give an extent of the PSF that is typically greater than or equal to approximately 5 pixels, i.e., 14 μm, on sensor 24.
Figs. 12A-12G shows MTF curves for an exemplary optical imaging assembly 152, in accordance with an embodiment of the present invention. Exemplary assembly 152 is described in more detail with reference to Fig. 13 below, and is one example of assembly 22. Assembly 152 is designed according to the steps and procedures of Figs. 4 and 10. The MTF curves may be produced during, and/or at the conclusion of, the design. Fig. 12A is for an object 100 mm from camera 20; Fig. 12B is for an object 150 mm from the camera; Fig. 12C is for an object 200 mm from the camera; Fig. 12D is for an object 300 mm from the camera; Fig. 12E is for an object 500 mm from the camera;
Fig. 12F is for an object 1000 mm from the camera; and Fig. 12G is for an object at infinity.
Each figure comprises four fields of view 0°, 12°, 25°, and 33°. The sagittal (S) and tangential (T) orientation curves are shown for each field of view.
Inspection of the MTF curves shows that all MTF values for spatial frequencies above approximately 50 cycles per mm are low, equaling approximately 0.35 or less. The low MTF values are a consequence of the distorted image on sensor 24. However, because the MTF values are all approximately equal, the corrections applied by engine 26 generate substantially the same improvement to all the images, so as to produce approximately equal and high MTFs for all the corrected images.
It will be understood that the above MTF properties are not just for assembly 152, but apply to a range of assemblies 22 generated according to the steps and procedures of Figs. 4 and 10, generating a camera with a depth of field from approximately 10 cm to infinity.
Fig. 13 is a schematic diagram of optical imaging assembly 152, in accordance with an embodiment of the present invention. Assembly 152 comprises five components: an aperture 151, three lenses 154, 156, 158 and an infra-red filter 160. Assembly 152 forms its image on a front surface 162 of sensor 24, so that the front surface acts as a focal plane of the assembly. Each lens has two aspheric surfaces, each surface being defined by expression (6):
Figure imgf000027_0001
where z is a surface sag, r is a radius measured to an optic axis, c, k are curvature and conic constants for the surface, and <X2, ... αg are parameters for the surface. Tables VI, VII, and VTH below show dimensions in mm and parameters of each of the ten surfaces of the assembly. In the tables, the surfaces of each component are distinguished by a suffix A or B. Thus, aperture 151 has surfaces 15 IA and 15 IB. In Table VI, the thickness gives the distance, measured along the optic axis, to the next surface, and thus gives spacings between the components. Thus, aperture 151 has a thickness of 5.0-10~2 mm, the distance from surface 151B to the following surface, 154A, is l.l-10"2, and the distance from the final surface of the assembly, surface 160B of filter 160, to the sensor, is 8.61- 10"! J31111-
Figure imgf000028_0001
Table VI
Figure imgf000028_0002
Table VII
Figure imgf000029_0001
Table VIII
Plastic E48R is produced by Zeon Corporation of Tokyo, Japan and has a refractive index approximately equal to 1.53. Material OKP4 is produced by Osaka Gas Chemicals Ltd., of Japan and has a refractive index approximately equal to 1.61. Glass BK7 is produced by Schott AG of Mainz, Germany and has a refractive index approximately equal to 1.52.
Table IX below gives Zernike fringe coefficient values for assembly 152.
Figure imgf000029_0002
Figure imgf000030_0001
Table IX
Inspection of Table IX shows that for all fields of view the values of Z4 are between approximately 0.6 and approximately 0.9, and the values of Z16 are between approximately -0.2 and approximately -0.3. Furthermore, the values of Zl, Z4, Z9, Z16, Z25, Z36, and Z37 in Table IX are all close to the target values in Table III.
Although the embodiments described above refer to certain specific digital filters, and particularly to a deconvolution filter (DCF), the principles of the present invention may similarly be applied in electronic cameras that use other types of digital image filters, as are known in the art.
It will thus be appreciated that the embodiments described above are cited by way of example, and that the present invention is not limited to what has been particularly shown and described hereinabove. Rather, the scope of the present invention includes both combinations and subcombinations of the various features described hereinabove, as well as variations and modifications thereof which would occur to persons skilled in the art upon reading the foregoing description and which are not disclosed in the prior art.

Claims

1. An optical imaging assembly having cylindrical symmetry, comprising a plurality of lenses having surfaces with curvatures and spacings between the surfaces, such that an optical image formed by the plurality of lenses has a defocus aberration coefficient greater than 0.1 at a focal plane of the assembly.
2. The assembly according to claim I3 wherein the defocus aberration coefficient is greater than or equal to 0.6.
3. The assembly according to claim 2, wherein the defocus aberration coefficient is less than or equal to 0.9.
4. The assembly according to claim 3, wherein a fifth order spherical aberration is between -0.2 and -0.3.
5. The assembly according to claim 4, wherein a seventh order spherical aberration is -0.5.
6. The assembly according to claim 5, wherein a ninth order spherical aberration is 0.
7. The assembly according to claim 1, wherein the curvatures are non-logarithmic.
8. The assembly according to claim 1, wherein the optical image is of one or more objects within a field of view subtending up to 33° to an optic axis of the assembly.
9. The assembly according to claim 1, wherein the optical image is of one or more objects which are located at distances greater than or equal to 10 cm from the assembly.
10. The assembly according to claim 9, wherein the assembly is configured so that a modulation transfer function generated by the one or more objects at the focal plane is less than or equal to 0.35 for spatial frequencies greater than 50 cycles per mm.
11. The assembly according to claim 1, wherein the optical image has an extent of a point spread function equal to or greater than 14 μm at the focal plane.
12. A method for optical imaging, comprising: formulating an optical specification of an optical imaging assembly comprising optical elements having surfaces, the optical specification comprising one or more aberration coefficients; generating a virtual phase element having phase coefficients which are negatives of the aberration coefficients; generating a design of the optical assembly comprising the optical elements and the virtual phase element so as to determine curvatures and spacings of the surfaces that will focus an image of an object on a focal plane of the assembly in the presence of the virtual phase element; and fabricating the optical assembly in accordance with the curvatures and spacings determined by the design.
13. The method according to claim 12, wherein the aberration coefficients comprise a de focus aberration coefficient greater than 0.1.
14. The method according to claim 13, wherein the defocus aberration coefficient is greater than or equal to 0.6.
15. The method according to claim 14, wherein the defocus aberration coefficient is less than or equal to 0.9.
16. The method according to claim 15, wherein the one or more aberrations comprise a fifth order spherical aberration that is between -0.1 and -0.3.
17. The method according to claim 12, wherein the curvatures are non-logarithmic.
18. The method according to claim 12, wherein the optical specification comprises a field of view of the object subtending up to 33° to an optic axis of the assembly.
19. The method according to claim 12, wherein the optical specification comprises a distance of the object from the assembly being greater than or equal to 10 cm from the assembly.
20. The method according to claim 19, wherein the optical specification comprises a modulation transfer function generated by the object at the focal plane being less than or ' equal to 0.35 for spatial frequencies greater than 50 cycles per mm.
21. The method according to claim 12, wherein the image has an extent of a point spread function equal to or greater than 14 urn at the focal plane.
22. The method according to claim 12, and comprising generating a digital filter to compensate for the one or more aberration coefficients, and coupling the digital filter to the optical assembly so that a filtered image formed by the digital filter from the image on the focal plane is free of aberrations.
23. The method according to claim 22, wherein the filtered image is free of aberrations for the object being located at a distance greater than 10 cm from the optical assembly.
24. An optical imaging assembly having cylindrical symmetry, comprising a plurality of lenses having surfaces with curvatures and spacings between the surfaces, such that an optical image formed by the plurality of lenses at a focal plane of the assembly has a modulation transfer function (MTF) for spatial frequencies greater than 50 cycles per millimeter that is less than or equal to 0.35.
25. A method for forming an optical imaging assembly, comprising: providing a plurality of lenses having cylindrical symmetry and surfaces with curvatures and spacings between the surfaces; and arranging the lenses so that an optical image formed by the plurality of lenses has a defocus aberration coefficient greater than 0.1 at a focal plane of the assembly.
26. Apparatus for optical imaging, comprising: an optical imaging assembly comprising optical elements having surfaces; and a processor which is configured to: formulate an optical specification of the optical imaging assembly, the optical specification comprising one or more aberration coefficients, generate a virtual phase element having phase coefficients which are negatives of the aberration coefficients; generate a design of the optical assembly comprising the optical elements and the virtual phase element so as to determine curvatures and spacings of the surfaces that will focus an image of an object on a focal plane of the assembly in the presence of the virtual phase element; and fabricate the optical assembly in accordance with the curvatures and spacings determined by the design.
27. A method for forming an optical imaging assembly, comprising: providing a plurality of lenses having cylindrical symmetry and surfaces with curvatures and spacings between the surfaces; and arranging the lenses so that an optical image formed by the plurality of lenses at a focal plane of the assembly has a modulation transfer function (MTF) for spatial frequencies greater than 50 cycles per millimeter that is less than or equal to 0.35.
28. An optical imaging assembly, comprising: an aperture having a semi-diameter of 7.5 10" 1 mm; and a plurality of lenses having surfaces with curvatures and spacings between the surfaces, wherein respective surfaces of the plurality of lenses are defined by an equation
Figure imgf000035_0001
where z is a sag of the surface, r is a radius measured to an optic axis of the assembly, c, k are respective curvature and conic constants for the surface, and (X2, ... αg are parameters for the surface, the lenses comprising: a first aspheric lens following the aperture, having a refractive index 1.53, a first surface wherein c = 9.57- 10"2, k = 9.74-101, α2 = -1.38-10"2, α3 = -9.0-10"1, 0C4 = 2.83,
015 = -2.26, αg = -2.96, aη = 2.69, αg = 1.61, and a second surface wherein c = -4.46-10"!, k = 0, (X2 = -9.83-10-2, α3 = 1.24-10"1, α4 = -4.25-10"1, α5 = 7.62-10"1, αg = -8.04-10"1, α7 = 4.73-10"1, α8 = -1.28-10"1; a second aspheric lens following the first aspheric lens, having a refractive index 1.61, a third surface wherein c = -7.60-10"!, ^ = .3.37,, and absolute values of α2, α3, α_j., α5, αg, aη, αg less than I, and a fourth surface wherein c = -2.42-10"!, k = 2.38, and absolute values of U2, 03, αφ α5, αg, aη, αg less than 1; and a third aspheric lens following the second aspheric lens, having a refractive index 1.53, a fifth surface wherein c = 8.1-10"!, ^- -4 Q3, and absolute values of α2, α3, αφ, 05, αg, aη, αg less than 1, and a sixth surface wherein c = 5.17-10"l, k = -4.07, and absolute values of α2, α3, 04, α5, αg, aη, αg less than 1.
PCT/IL2006/001294 2003-01-16 2006-11-09 Optics for an extended depth of field WO2007054938A2 (en)

Priority Applications (5)

Application Number Priority Date Filing Date Title
US12/093,394 US7773316B2 (en) 2003-01-16 2006-11-09 Optics for an extended depth of field
TW096108966A TWI427554B (en) 2006-03-31 2007-03-15 Optical imaging assembly and method for forming the same, and apparatus and method for optical imaging
US12/847,651 US8014084B2 (en) 2003-01-16 2010-07-30 Optics for an extended depth of field
US13/196,051 US8294999B2 (en) 2003-01-16 2011-08-02 Optics for an extended depth of field
US13/625,390 US8611030B2 (en) 2003-01-16 2012-09-24 Optics for an extended depth of field

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
US73544105P 2005-11-10 2005-11-10
US60/735,441 2005-11-10
US11/278,255 US20070236573A1 (en) 2006-03-31 2006-03-31 Combined design of optical and image processing elements
US11/278,255 2006-03-31

Related Parent Applications (4)

Application Number Title Priority Date Filing Date
PCT/IL2004/000040 Continuation-In-Part WO2004063989A2 (en) 2003-01-16 2004-01-15 Camera with image enhancement functions
US10/541,967 Continuation-In-Part US7627193B2 (en) 2003-01-16 2004-01-15 Camera with image enhancement functions
US11/278,255 Continuation-In-Part US20070236573A1 (en) 2003-01-16 2006-03-31 Combined design of optical and image processing elements
US54196706A Continuation-In-Part 2005-10-14 2006-06-23

Related Child Applications (3)

Application Number Title Priority Date Filing Date
US12/093,394 A-371-Of-International US7773316B2 (en) 2003-01-16 2006-11-09 Optics for an extended depth of field
US12/847,651 Continuation US8014084B2 (en) 2003-01-16 2010-07-30 Optics for an extended depth of field
US13/196,051 Continuation US8294999B2 (en) 2003-01-16 2011-08-02 Optics for an extended depth of field

Publications (2)

Publication Number Publication Date
WO2007054938A2 true WO2007054938A2 (en) 2007-05-18
WO2007054938A3 WO2007054938A3 (en) 2009-02-12

Family

ID=38023666

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/IL2006/001294 WO2007054938A2 (en) 2003-01-16 2006-11-09 Optics for an extended depth of field

Country Status (1)

Country Link
WO (1) WO2007054938A2 (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2931255A1 (en) * 2008-05-15 2009-11-20 Dxo Labs Sa OPTICAL SYSTEM AND METHOD FOR DESIGNING THE SAME
US20110025877A1 (en) * 2008-04-03 2011-02-03 Gallagher Dennis J Imaging System Including Distributed Phase Modification And Associated Methods
WO2011058236A1 (en) * 2009-11-16 2011-05-19 Dxo Labs Optical system, and related design method
US7999851B2 (en) 2007-05-24 2011-08-16 Tessera Technologies Ltd. Optical alignment of cameras with extended depth of field
US9118826B2 (en) 2009-03-19 2015-08-25 Digitaloptics Corporation Dual sensor camera
US9282252B2 (en) 2009-05-04 2016-03-08 Digitaloptics Corporation Dual lens digital zoom
CN111650747A (en) * 2020-05-20 2020-09-11 中国科学院西安光学精密机械研究所 Method and device for real-time correction of defocusing aberration of collimator

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040103024A1 (en) * 2000-05-24 2004-05-27 Matchcraft, Inc. Online media exchange

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040103024A1 (en) * 2000-05-24 2004-05-27 Matchcraft, Inc. Online media exchange

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7999851B2 (en) 2007-05-24 2011-08-16 Tessera Technologies Ltd. Optical alignment of cameras with extended depth of field
US20110025877A1 (en) * 2008-04-03 2011-02-03 Gallagher Dennis J Imaging System Including Distributed Phase Modification And Associated Methods
US8922700B2 (en) * 2008-04-03 2014-12-30 Omnivision Technologies, Inc. Imaging system including distributed phase modification and associated methods
FR2931255A1 (en) * 2008-05-15 2009-11-20 Dxo Labs Sa OPTICAL SYSTEM AND METHOD FOR DESIGNING THE SAME
US9118826B2 (en) 2009-03-19 2015-08-25 Digitaloptics Corporation Dual sensor camera
US9282252B2 (en) 2009-05-04 2016-03-08 Digitaloptics Corporation Dual lens digital zoom
WO2011058236A1 (en) * 2009-11-16 2011-05-19 Dxo Labs Optical system, and related design method
CN111650747A (en) * 2020-05-20 2020-09-11 中国科学院西安光学精密机械研究所 Method and device for real-time correction of defocusing aberration of collimator

Also Published As

Publication number Publication date
WO2007054938A3 (en) 2009-02-12

Similar Documents

Publication Publication Date Title
TWI427554B (en) Optical imaging assembly and method for forming the same, and apparatus and method for optical imaging
US8014084B2 (en) Optics for an extended depth of field
US8611030B2 (en) Optics for an extended depth of field
US11721002B2 (en) Imaging system and method for imaging objects with reduced image blur
US20070236574A1 (en) Digital filtering with noise gain limit
US8126287B2 (en) Camera with image enhancement functions
US20070239417A1 (en) Camera performance simulation
WO2007054938A2 (en) Optics for an extended depth of field
CA2577735A1 (en) Extended depth of field using a multi-focal length lens with a controlled range of spherical aberration and centrally obscured aperture
EP1672912B1 (en) Method for producing an optical system including an electronic image enhancement processor
Adelsberger Design guidelines for wavefront coding in broadband optical systems
Scrymgeour et al. Advanced Imaging Optics Utilizing Wavefront Coding.
Ma Composite Likelihood Inference for Multivariate Finite Mixture Models and Application to Flow Cytometry

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application
NENP Non-entry into the national phase in:

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 12093394

Country of ref document: US

122 Ep: pct application non-entry in european phase

Ref document number: 06809854

Country of ref document: EP

Kind code of ref document: A2