RGBW sensor array

Information

  • Patent Grant
  • 8878967
  • Patent Number
    8,878,967
  • Date Filed
    Monday, October 11, 2010
    13 years ago
  • Date Issued
    Tuesday, November 4, 2014
    9 years ago
Abstract
A color filter enhancement method for a portable digital image acquisition device includes digitally exposing color pixels of a color sensor array for a first digital exposure duration and digitally exposing white pixels of a color sensor array for a second digital exposure time shorter than the first digital exposure duration. A color filter enhanced digital image is generated using data from both the color pixels exposed for the first digital exposure duration and the white pixels exposed for the second digital exposure duration.
Description
BACKGROUND OF THE INVENTION

1. Field of the Invention


The present invention relates to an image processing method and apparatus.


2. Description of the Related Art


Sensor arrays for digital cameras exist for capturing color photos. Sensors known as RGBW sensors are used for capturing red, green, blue colors, and for capturing luminance information for multiple pixels of an array of pixels. The red, green and blue pixels include filters such that only certain narrow ranges of wavelengths of incident light are counted. The white pixels capture light of red, green and blue wavelengths, i.e., of a broader range of wavelengths than any of the blue, green and red pixels. Thus, the white pixels are typically brighter than any of the blue, red and green pixels if they are exposed for the same duration.


Noise removal algorithms tend to blur face regions in an undesirable manner. Noise removal algorithms are described at U.S. patent application Ser. Nos. 11/856,721 and 11/861,257 and are hereby incorporated by reference, as are Ser. Nos. 10/985,650, 11/573,713, 11/421,027, 11/673,560, 11/319,766, 11/744,020, 11/753,098, 11/752,925, and 12/137,113, which are assigned to the same assignee as the present application and are hereby incorporated by reference.


Kodak has developed a RGBW color filter pattern differing from the previously known Bayer Color Filter. The RGBW pattern of Kodak is referred to as a Color Filter Array (CFA) 2.0. One half of cells in a RGBW pattern are panchromatic, i.e. sensing all color spectrum (Y component)—usually called white cells. This way more light energy is accumulated in the same amount of time than for color pixels. A Bayer filter uses only ⅓ (˜0.33) of color spectrum energy. An RGBW filter uses 4/6 (˜0.67) of energy, where ½ comes from white cells and ⅙ from RGB cells.


CFA Array looks something like the following:

    • WBWG . . .
    • BWGW . . .
    • WGWR . . .
    • RWRW . . .


In this context, the following are incorporated by reference: U.S. Pat. Nos. 7,195,848, 7,180,238, 7,160,573, 7,019,331, 6,863,368, 6,607,873, 6,602,656, 6,599,668, 6,555,278, 6,387,577, 6,365,304, 6,330,029, 6,326,108, 6,297,071, 6,114,075, 5,981,112, 5,889,554, 5,889,277, 5,756,240, 5,756,239, 5,747,199, 5,686,383, 5,599,766, 5,510,215, 5,374,956, and 5,251,019.


Two source images nominally of the same scene may be used to produce a single target image of better quality or higher resolution than either of the source images.


In super-resolution, multiple differently exposed lower resolution images can be combined to produce a single high resolution image of a scene, for example, as disclosed in “High-Resolution Image Reconstruction from Multiple Differently Exposed Images”, Gunturk et al., IEEE Signal Processing Letters, Vol. 13, No. 4, April 2006; or “Optimizing and Learning for Super-resolution”, Lyndsey Pickup et al, BMVC 2006, 4-7 Sep. 2006, Edinburgh, UK, hereby incorporated by reference. However, in super-resolution, blurring of the individual source images either because of camera or subject motion are usually not of concern before the combination of the source images.


U.S. Pat. No. 7,072,525, incorporated by reference, discloses adaptive filtering of a target version of an image that has been produced by processing an original version of the image to mitigate the effects of processing including adaptive gain noise, up-sampling artifacts or compression artifacts.


US published applications 2006/0098890, 2007/0058073, 2006/0098237, 2006/0098891, European patent EP1779322B1, and PCT Application No. PCT/EP2005/011011, each hereby incorporated by reference, describe uses of information from one or more presumed-sharp short exposure time (SET) preview images to calculate a motion function for a fully exposed higher resolution main image to assist in the de-blurring of the main image.


Indeed many other documents, including US 2006/0187308, Suk Hwan Lim et al.; and “Image Deblurring with Blurred/Noisy Image Pairs”, Lu Yuan et al, SIGGRAPH07, Aug. 5-9, 2007, San Diego, Calif. are directed towards attempting to calculate a blur function in the main image using a second reference image before de-blurring the main image.


Other approaches, such as disclosed in US2006/0017837 have involved selecting information from two or more images, having varying exposure times, to reconstruct a target image where image information is selected from zones with high image details in SET images and from zones with low image details in longer exposure time images.


SUMMARY OF THE INVENTION

A color filter enhancement method is provided for a portable digital image acquisition device. The method includes digitally exposing color pixels of a color sensor array for a first digital exposure duration. White pixels of a color sensor array are digitally exposed for a second digital exposure time shorter than the first digital exposure duration. A color filter enhanced digital image is generated using data from both the color pixels exposed for the first digital exposure duration and the white pixels exposed for the second digital exposure duration shorter than the first digital exposure duration. The color filter enhanced digital image and/or a further processed version is stored, transmitted, communicated, displayed, and/or projected.


The second digital exposure time may be less than half of the first digital exposure time, for example, it may be approximately a third of the first digital exposure time.


The digitally exposing of the color pixels and the white pixels for different exposure times may include clocking the color pixels and the white pixels independently.


The digitally-exposing of the color pixels and the white pixels for different exposure times may involve including sensor data over different temporal ranges. The different temporal ranges may be overlapping. A first temporal range corresponding to the digitally-exposing of the color pixels may include an entire second temporal range corresponding to the digitally-exposing of the white pixels.


The color pixels may include greater motion blurring effect than the white pixels due to the color pixels being digitally-exposed for a longer duration than the white pixels. The method may further include compensating blurring in the color pixels using less-blurred data from the white pixels.


The color sensor array may include a CMOS-based sensor.


One or more processor-readable media are also provided that have code embedded therein for programming the processor to perform a color filter enhancement method in accordance with any of the methods described herein.


A portable digital image acquisition device is also provided including optics and a color sensor array for acquiring a color digital image, a processor, and one or more processor-readable media having code embedded therein for programming the processor to perform a color filter enhancement method that comprises any of the methods described herein.





BRIEF DESCRIPTIONS OF THE DRAWINGS

Several embodiments will now be described, by way of example, with reference to the accompanying drawings, in which:



FIG. 1 is a block diagram illustrating the processing of images prior to adaptive filtering according to a first embodiment of the present invention;



FIG. 2 illustrates corresponding grid points from a preview and a full resolution image used in the processing of FIG. 1;



FIG. 3 illustrates the adaptive filtering of images in R/G/B color space according to one implementation of the present invention;



FIG. 4 illustrates the adaptive filtering of images in YCbCr color space according to another implementation of the present invention;



FIGS. 5(
a) and (b) illustrate in more detail the adaptive filtering of images according to two variants of the first embodiment of the invention;



FIG. 6 illustrates a sliding vector employed in the filtering of FIG. 5 at successive iterations for L=3;



FIG. 7 is a block diagram illustrating the processing of images prior to adaptive filtering according to a second embodiment of the present invention;



FIG. 8 shows the timing involved in acquiring two images for use in a further embodiment of the present invention;



FIGS. 9(
a)-9(e) shows some image data produced during the image acquisition of FIG. 8;



FIG. 10 is a block diagram illustrating a method of color filter enhancement for a portable digital image acquisition device in accordance with certain embodiments; and



FIG. 11 is a block diagram illustrating a portable digital image acquisition device with color filter enhancement in accordance with certain embodiments.





DETAILED DESCRIPTIONS OF THE EMBODIMENT

Improved methods are described of combining a sharp image and a blurred image of differing resolution and exposure to produce a relatively high resolution, fully exposed and relatively sharp image.


Referring now to FIG. 1, in accordance with certain embodiments, a well-exposed blurred relatively low resolution image 12 and a sharp but under-exposed full resolution image 10 are available for processing with a view to combining the images to produce an improved quality full resolution image.


The size of the lower resolution image 12 is O×P and the size of the under-exposed full resolution image 10 is Q×R, with O<Q and P<R.


Where the images are acquired in a digital image acquisition device such as a digital stills camera, camera phone or digital video camera, the lower resolution image 12 may be a preview image of a scene acquired soon before or after the acquisition of a main image comprising the full resolution image 10, with the dimensions of the preview and full resolution images depending on the camera type and settings. For example, the preview size can be 320×240 (O=320; P=240) and the full resolution image can be much bigger (e.g. Q=3648; R=2736).


In accordance with certain embodiments, adaptive filtering (described in more detail later) is applied to the (possibly pre-processed) source images 10, 12 to produce an improved filtered image. Adaptive filtering requires an input image (referred to in the present specification as x(k)) and a desired image (referred to in the present specification as d(k)) of the same size, with the resultant filtered image (referred to in the present specification as y(k)) having the same size as both input and desired images.


As such, in one embodiment, the preview image is interpolated to the size Q×R of the full resolution image.


It will be seen that in interpolating the preview image, a misalignment between the interpolated image 14 and the full resolution image might exist. As such, in this embodiment, the images are aligned 16 to produce an aligned interpolated preview image 18 and an aligned full resolution image 20. Any known image alignment procedure can be used, for example, as described in Kuglin C D., Hines D C. “The phase correlation image alignment method”, Proc. Int. Conf. Cybernetics and Society, IEEE, Bucharest, Romania, September 1975, pp. 163-165, hereby incorporated by reference.


Other possible image registration methods are surveyed in “Image registration methods: a survey”, Image and Vision Computing 21 (2003), 977-1000, Barbara Zitova and Jan Flusser, hereby incorporated by reference.


Alternatively, the displacements between the images 10 and 12/14 can be measured if camera sensors producing such a measure are available.


In any case, either before or during alignment, the full resolution image can be down-sampled to an intermediate size S×T with the preview image being interpolated accordingly to produce the input and desired images of the required resolution, so that after alignment 16, the size of the aligned interpolated image and the aligned full resolution image will be S×T (S≦Q, T≦R).


These images are now subjected to further processing 22 to compute the input and desired images (IMAGE 1 and IMAGE 2) to be used in adaptive filtering after a decision is made based on the displacement value(s) provided from image alignment 16 as indicated by the line 24.


In real situations, there may be relatively large differences between the images 10, 14, with one image being severely blurred and the other one being under-exposed. As such, alignment may fail to give the right displacement between images.


If the displacement values are lower than a specified number of pixels (e.g. 20), then the full resolution aligned image 20 is used as IMAGE 1 and the aligned interpolated preview image 18 is used as IMAGE 2.


Otherwise, if the displacement values are higher than the specified number of pixels, several alternatives are possible for IMAGE 2, although in general these involve obtaining IMAGE 2 by combining the interpolated preview image 14 and the full resolution image 10 in one of a number of manners.


In a first implementation, we compute two coefficients c1 and c2 and the pixel values of IMAGE 2 are obtained by multiplying the pixel values of the full resolution image 10 with c1 and adding c2. These coefficients are computed using a linear regression and a common form of linear regression is least square fitting (G. H. Golub and C. F. Van Loan, Matrix Computations. John Hopkins University Press, Baltimore, Md., 3rd edition, 1996), hereby incorporated by reference. Referring to FIG. 2, a grid comprising for example 25 points is chosen from the preview image 12 and the corresponding 25 grid points from the full resolution image 10. If one pixel of the preview image has the coordinates (k,l), the corresponding chosen pixel from the full resolution image has the coordinates







(

(


k
·

Q
O


,

l
·

R
P



)

)

.




Therefore we obtain two 5×5 matrices, M1 that corresponds to the pixel values chosen from the preview image and M2 that corresponds to the pixel values chosen from the full resolution image. Two vectors are obtained from the pixel values of these matrices by column-wise ordering of M1 (a=(ai) and M2 b=(bi)). We therefore have pairs of data (ai,bi) for i=1, 2, . . . ,n, where n=25 is the total number of grid points from each image. We define the matrix






V
=


(





a
1


1







a
2


1







a
n


1




)

.






The coefficient vector c=[c1 c2] is obtained by solving the linear system VTVc=VTb. The linear system can be solved with any known method.


Another alternative is to amplify the pixels of the under-exposed image 10 with the ratio of average values of the 25 grid points of both images 10, 12 and rescale within the [0-255] interval for use as IMAGE 2.


In a still further alternative, IMAGE 2 is obtained by combining the amplitude spectrum of the interpolated blurred preview image 14 and the phase of the under-exposed full resolution image 10. As such, IMAGE 2 will be slightly deblurred, with some color artifacts, although it will be aligned with the under-exposed image 10. This should produce relatively fewer artifacts in the final image produced by adaptive filtering.


Alternatively, instead of computing FFTs on full resolution images to determine phase values, an intermediate image at preview resolution can be computed by combining the amplitude spectrum of the blurred image 12 and the phase of a reduced sized version of the under-exposed image 10. This can then be interpolated to produce IMAGE 2.


Another possibility is to use as IMAGE 2, a weighted combination of image 20 and image 18, e.g. 0.1*(Image 18)+0.9*(Image 20). This can be used if the preview image 12 has large saturated areas.


In any case, once the processing 22 is complete, two images of similar size are available for adaptive filtering 30. See FIGS. 3-4 in this context.


In a first implementation, the input and desired images are in RGB color space, FIG. 3, whereas in another implementation the input and desired images are in YCC space, FIG. 4. For the RGB case, one color plane (e.g. G plane) is selected from both images and the computed filter coefficients from adaptive filtering are used to update the pixel values for all color planes. The filter coefficients w(k) are obtained at each iteration of the filter 36. The updated pixel value for all color planes will be yG(k)=w(k)·xG(k), yR(k)=w(k)·xR(k), yB(k)=w(k)·xB(k), where xR(k), xG(k), xB(k) are the sliding vectors 32 for the R,G,B planes respectively. This provides a solution of reduced numerical complexity vis-à-vis filtering all three color planes.


In the YCC case, the Y plane is selected with the Cb and Cr planes being left unchanged.


Referring now to FIG. 5(a), where the adaptive filtering of FIGS. 3 and 4 is shown in more detail. Two sliding one-dimensional vectors 32, 34 with the dimension L are created, L being the length of the adaptive filter. Within the adaptive filter, the input signal x(k) is the first vector signal 32, while the desired signal d(k) is second vector 34.


In the simplest implementation, L=1 and this can be used if the original image acquisition device can provide good quality under-exposed pictures with a low exposure time. Where the acquisition device produces low quality and noisy under-exposed images, a longer filter length L should be chosen (e.g. 2 or 3 coefficients).


The sliding vectors 32, 34 are obtained from the columns of the image matrices, as illustrated at FIG. 6. The vectors scan both matrices, column by column and with each iteration of the adaptive filter, the following pixel value is added to the vector and the trailing pixel value is discarded.


When the vectors 32, 34 are combined in the adaptive filter 36, the most recent pixel value added to the first sliding vector 32 is updated. In the preferred embodiment, the updated pixel is the dot product of the filter coefficients and the L pixel values of the first vector. Any adaptive algorithm (Least Mean Square based, Recursive Least Square based) can be applied and many such algorithms can be found in S. Haykin, “Adaptive filter theory”, Prentice Hall, 1996. Preferably, the sign-data LMS described in Hayes, M, Statistical Digital Signal Processing and Modeling, New York, Wiley, 1996, incorporated by reference, is employed.


The formulae are:

x(k)=[x(k),x(k−1) . . . x(k−L+1)],
w(k)=[w(k),w(k−1) . . . w(k−L+1)],
y(k)=w(kx(k),
e(k)=d(k)−y(k),
w(k+1)=w(k)+μ(ke(k)·sign(x(k))=w(k)+μ(ke(k),

where

  • w(k) are the filter coefficients calculated within the filter 36,
  • μ(k) is the step size (fixed or variable),
  • x(k) is the most recent pixel value(s) of the sliding vector 32 from Image 1 (it has always positive values),
  • d(k) is the most recent pixel value(s) of the sliding vector 34 from Image 2,
  • y(k) is the scalar product of the sliding vector 32 and the filter coefficients vector w,
  • e(k) is the error signal computed as the difference between d(k) and y(k).


Other considered variants were:

w(k+1)=w(k)+μ(ke(kx(k)(standard LMS) or
w(k+1)=w(k)+μ(ke(k)/(1+x(k))


The term 1+x(k) is used above to avoid the division by zero. Alternatively, the formula:







w


(

k
+
1

)


=


w


(
k
)


+


μ


(
k
)


·


e


(
k
)



x


(
k
)










could be used, with any zero-valued x pixel value replaced with a 1.


In a further variant, the step size μ(k) is variable as follows:







μ


(
k
)


=




1
-
α


x


(
k
)








or






μ


(
k
)



=



1
-
α


max


(

β
,

x


(
k
)



)



.






So, using the above formula:

w(k+1)=w(k)+μ(ke(k)·sign(x(k))=w(k)+μ(ke(k)

this gives:








w


(
k
)


+



1
-
α


x


(
k
)





(


d


(
k
)


-


w


(
k
)


·

x


(
k
)




)



=



w


(
k
)


+


d


(
k
)



x


(
k
)



-

w


(
k
)


-

α
·


d


(
k
)



x


(
k
)




+

α
·

w


(
k
)




=


α
·

w


(
k
)



+


(

1
-
α

)

·


d


(
k
)



x


(
k
)










If μ(k)=μ=1−α, α very close to 1 (e.g. 0.99999), for L=1, we have








w


(

k
+
1

)


=


w


(
k
)


+


μ


(
k
)


·


e


(
k
)



x


(
k
)






,





with vectors being replaced with scalars. Therefore, for this particular fixed step size, the sign-data LMS and the previous equation are equivalent.


The β parameter can be used in order to avoid division by zero and to over-amplify any black pixels. β is preferably in the interval [1 . . . 10], and preferably in the interval [5 . . . 10], particularly if the under-exposed image is too dark. If not, β=1 is enough.


Some thresholds or resetting for the filter coefficients w(k) or output values y(k) can be imposed in order to avoid artifacts in the filtered image 38. An upper threshold, δ, is imposed for the values that can be allowed for the coefficients of w(k) (i.e. wi(k)=δ for any i=1 . . . L, if its computed value at iteration k is above δ). A suitable threshold value for the mentioned LMS algorithm, can be chosen as







δ
=

1
+


b
_


4
·

a
_





,





where b and ā are the average values of above mentioned vectors b and a respectively. Also, the filter output can be forced to be within the [0 255] interval if uint8 images are used. As can be seen, the updated pixel values y(k) replace the old pixel values x(k) and can be taken into account for the next sliding vectors.


The updated color matrix 38 is completed when the last pixel from the last column has been updated. If filtering has been performed in RGB space, then a final reconstructed image 40 is obtained by concatenating the R/G/B updated matrices. Alternatively, if filtering has been performed in YCC space, the concatenated updated Y plane, i.e. matrix 38, with unchanged Cb and Cr planes of the under-exposed image 10 can be converted back to RGB color space.


The filtering can be repeated with the reconstructed image 40 replacing the under-exposed image, i.e. IMAGE 1.


In this case, adaptive filtering can be performed on the Y plane of an image converted from RGB space, if previous filtering had been performed in RGB space; or alternatively filtering can be performed on an RGB color plane of an image converted from YCC space, if previous filtering had been performed on the Y plane.


It will also be seen that filtering can be operated column wise or row wise. As such, adaptive filtering can be performed first column or row wise and subsequently in the other of column or row wise.


In each case where filtering is repeated, it has been found that the quality of the reconstructed image after two filtering operations is superior than for each individual filtering result.


Referring to FIG. 5(b), in some cases saturation problems might appear in the filtered image, especially when the coefficient c1 has a large value (e.g. when using a very dark under-exposed image and very light blurred image). This saturation can be avoided using, for example, techniques described in Jourlin, M., Pinoli, J. C.: “Logarithmic image processing, the mathematical and physical framework fro the representation and processing of transmitted images” Advances in Imaging and Electron Physics 115 (2001) 129-196; or Deng, G., Cahill, L. W., Tobin, G. R.: “The study of logarithmic image processing model and its application to image enhancement”. IEEE Trans. on Image Processing 4 (1995) 506-512, hereby incorporated by reference.


Therefore, the pixel value of the filtered image z(k) is generated by the following formula:







z


(
k
)


=

D
-


D


(

1
-


x


(
k
)


D


)



w


(
k
)









where D is the maximum permitted value (e.g. 255 for a 8 bit representation of images). The adaptive filter provides the first filter coefficient w(k) computed using the error signal e(k). Another alternative to reduce saturation problems is to reduce the value of the step size μ(k).


Referring now to FIG. 7, in a second embodiment of the invention, an under-exposed relatively-sharp low resolution image and a full resolution blurred image 72 are available. The low resolution image, for example, a preview image as before, is interpolated and aligned with the full resolution image to produce image 70.


A PSF estimation block 74 computes a PSF for the blurred image 72, from the interpolated preview 70 and the full resolution image 72, using any suitable method such as outlined in the introduction.


The blurred 72 image is then deblurred using this estimated PSF to produce a relatively deblurred image 76. Examples of deblurring using a PSF are disclosed in “Deconvolution of Images and Spectra” 2nd. Edition, Academic Press, 1997, edited by Jannson, Peter A. and “Digital Image Restoration”, Prentice Hall, 1977 authored by Andrews, H. C. and Hunt, B. R., hereby incorporated by reference.


Prior to adaptive filtering, the average luminance of the interpolated preview image 70 is equalized in processing block 78 with that of the full resolution (relatively) deblurred image 76. Preferably, this comprises a gamma (γ) amplification of the under-exposed image. The exact value of gamma is determined by obtaining a ratio of average luminance ( Y in YCC format) for the blurred full resolution and the preview image, and then using this ratio as an index for a look-up table to return γ. The deblurred full resolution image 76 is then chosen as IMAGE 2 and the interpolated/aligned/luminance equalized preview image produced by the processing block 78 is chosen as IMAGE 1.


Adaptive filtering is then applied and re-applied if necessary to IMAGE 1 and IMAGE 2 as in the first embodiment. Again when repeating adaptive filtering, the under-exposed image, i.e. IMAGE 1 is replaced with the reconstructed one.


In another embodiment, the quality of the reconstructed image 76 produced by adaptive filtering may not be good enough, especially if the PSF is relatively large. In such cases, de-blurring using the PSF may not be used, because can it introduce significant ringing.


In cases such as this, re-applying adaptive filtering as in the first embodiment can attenuate the blurring artifacts in the original image 72 and improve the quality of the image to some extent.


Again, the adaptive filtering can be performed on Y plane if RGB filtering had been performed previously and on the RGB color space if Y filtering had been performed previously.


Again, filtering can be operated on columns or rows, and sequentially on columns and rows.


It has also been found that the second embodiment is useful, if the ratio between the full resolution image 72 and the preview image sizes is less than three and the preview image is not too noisy. If this is not the case, the filtered image can have a lower quality than that obtained by deblurring the blurred image with a very good PSF estimation such as described in the introduction.


In any of the above embodiments, a single preview image is described as being interpolated to match the resolution of the full resolution image. However, it will also be appreciated that super-resolution of more than 1 preview image, nominally of the same scene, could also be used to generate the interpolated images 14, 70 of the embodiments.


In the above embodiments, the short-exposure time (presumed sharp) image is described as comprising a preview image acquired either soon before or after acquisition of a main high resolution image.


In a further refined embodiment, the two images are acquired within the longer time period of acquisition of the relatively blurred image. In an implementation of this embodiment, an image acquisition device including a CMOS sensor which allows for a non-destructive readout of an image sensor during image acquisition is employed to acquire the images.


A schematic representation of the timing involved in acquiring these images is explained in relation to FIG. 8. For a dark scene, the exposure time Tlong required to expose the image F properly can result in motion blur caused by hand jitter. Nonetheless, using a non-destructive sensor, it is possible to have an intermediate reading at Tshort providing an under-exposed (noise prone), but sharp image G.


In certain embodiments, the read-out of the under-exposed image is placed mid-way through the longer exposure period, i.e between T0 and T0+Tshort. As such, the actual exposing scheme goes as follows:

    • At t=0 start exposing
    • At t=T0 take the first readout to obtain G′
    • At t=T0+Tshort take the second readout to obtain G″
    • The short exposed image is G=G′−G″
    • At t=tlong take the third (last) readout to obtain the well-exposed frame, F.
    • Reset the image sensor.


This means that statistically, the chances of content differences between the short exposure and the long exposure images G and F are minimized. Again, statistically, it is therefore more likely that the differences are caused only by the motion existing in the period [0, Tlong]. The well exposed picture is blurred by the motion existing in its exposure period, while the other is not moved at all, i.e. the motion blur makes the content differences.


Referring now to FIGS. 9(a)-9(e), a still image of a scene is recorded. The period T0 is chosen to be long enough so that motion appears in the image G′ read at time T0, FIG. 9(c).


The values of the PSF for this image are shown in FIG. 9(a). From T0 to Tshort there is not enough time for extra motion to appear. However, the entire interval, [0; T0+Tshort], is long enough so that the resulting image G″, FIG. 9(d), will be blurred as can be seen from the corresponding PSF values of FIG. 9(b). The resulting under-exposed image, G=G″−G′, FIG. 9 (e), is not blurred as can be seen from the small difference between the PSF values for the original images G″ and G′.


The image G can now be combined with the image F through adaptive filtering as described above and in particular in relation to the second embodiment, luminance enhancement can be performed on the image G before being combined with the image F.


Subsequent to producing the filtered image 40 through one or more steps of adaptive filtering, the filtered image can be subjected to further processing to improve its quality further.


The noise correction of the filtered image can be performed using a modified version of the Lee Least mean square (LLMSE) filter. In the following example, G1 is the filtered image, G1x is the convolution of G1 with an X×X uniform averaging kernel; so G13 is the convolution of G1 with a 3×3 uniform averaging kernel; and G17 is the convolution of G1 with a 7×7 uniform averaging kernel.


The noise cleared picture is: G2=αG1x+(1−α)G1

    • where






α
=


s
n



s
n

+

s
F









    • SG1 is the filtered image standard deviation computed for a 5×5 vicinity of a pixel;

    • SF is the well-exposed image squared standard deviation computed for a 3×3 vicinity of the corresponding pixel; and

      Sn=|SF−SG1|





If SF is smaller than a predetermined threshold (meaning that the current pixel in a perfectly uniform area) then G1x=G17, otherwise (in the current pixel neighborhood there is an edge) G1x=G13. It will therefore be seen that where the variation around a pixel is high, G2 is approximately equal to G1.


As discussed, the under-exposed acquired image has intensities in the lower part of the range (darkness range). The spectral characteristics of the cameras, in this area, differ from those of normally exposed areas. Therefore, the adaptively filtered image, G1 or G2, may have deviations in color depending on whether noise filtering has been applied or not. To compensate for these deviations, a rotation or a translation in the (Cb,Cr) plane can be applied. The parameter values for these operations will depend on the camera and number of exposure stops between the well-exposed and the under-exposed images. One exemplary scheme for color correction in RBG space is as follows:

    • Compute the average luminance: ( YF, YG2)
    • Compute the color averages ( RF, RG2, GrF, GrG2BF, BG2)
    • Correct G2 to obtain G3 as follows:

      ΔR=( RG2YG2)+( YFRF);RG3(i,j)=RG2(i,j)−ΔR
      ΔGr=( GrG2YG2)+( YFGrF);GrG3(i,j)=GrG2(i,j)−ΔGr
      ΔB=( BG2YG2)+( YFBF);BG3(i,j)=BG2(i,j)−ΔB
      with ∀(i,j)ε[1,M]×[1,N]


In an alternative embodiment, an advantageous implementation may use a modified CMOS sensor array. In particular, a form of imaging sensor referred to as the Color Filter Array (CFA) incorporates an alternative arrangement of light sensitive elements to the more conventional Bayer pattern array employed by most standard sensors. Specifically, a CFA sensor incorporates an array of pixels arranged in the pattern:

    • WBWG . . .
    • BWGW . . .
    • WGWR . . .
    • RWRW . . . and so on . . .


This contrasts with the more conventional Bayer pattern, viz:

    • BGBG . . .
    • GRGR . . .
    • BGBG . . .
    • GRGR . . .


Note that the CFA arrangement includes special pixels which are sensitive to white light, rather than being specifically tuned to a Red, Green or Blue spectral response. As these white pixels are sensitive to white light, which combines all three color sensitivities, they do tend to charge more rapidly than the conventional pixels which are sensitized to Red, Green or Blue spectral frequencies.


In this alternative embodiment, the sensor array is designed to enable the white pixels which are comingled with the RGB pixels to be clocked independently of those color specific pixels.


Such an arrangement implies that these W pixels are arranged on a separate data bus and are clocked independently of the RGB pixels. Thus, it is possible to extract a “luminance-only” image from the sensor independently of the “chrominance image” provided by the RGB pixels. This involves additional real estate on the sensor IC, but fortunately, as this embodiment may use a smaller image, the density of W pixels employed by the standard CFA sensor may be advantageous reduced. In practice a 4:1 or 8:1 scaling-down of the “luminance image” is adequate and implies that only 1 in 16, or even 1 in 64 pixels should be W pixels.


Now as the luminance-only image obtained from the white pixels will charge more quickly it can be offloaded form the sensor as a fully exposed image about three times faster than the RGB pixels. Where an underexposed image is sufficient the speed of acquisition of this luminance image can be significantly more than an order of magnitude faster than the corresponding RGB image. Further, this luminance-only image derived from the W pixels can be offloaded while the RGB pixels are still charging. This ensures that sharp, underexposed images can be obtained at a time which overlaps the acquisition of the normally exposed, blurred main “chrominance image”.


The offloading of the sharp, underexposed image may be timed to coincide with the mid-point of the acquisition of the normally exposed, blurred color image.


Other advantages of this modified sensor arrangement should also be noted. For example, it is possible to implement an embodiment described in US published patent application US20050041121, hereby incorporated by reference, where a first luminance image is obtained immediately prior to the main RGB acquisition and the actuation of the flash; a second luminance image is obtained during the main RGB acquisition and these two luminance images are used to determine the location of potential red-eye defects. The final determination and correction of red-eye defects can be performed on the main RGB image. This approach has the advantage of eliminating the need to perform alignment between the pair of luminance images and the final RGB image can be enhanced (sharpened) by combining it with one of the two temporally proximate luminance images.


In another embodiment, a face tracking algorithm can be run on the luminance images acquired form such a modified sensor. Such face tracking algorithm may be based on classifier cascades trained on luminance images this approach enables the face tracker to operate directly on the luminance images acquired from the W pixels. Advantageously, where such face tracker algorithm operates using a fixed size classifier chain as described, for example, at U.S. Pat. Nos. 7,315,631, 7,403,643, 7,317,815, 7,269,292, and 7,316,630, and US published application no. 2008/0037840, 2008/0037839 and 2008/0037827, and U.S. application Ser. No. 12/063,089, each assigned to the same assignee and incorporated by reference, the main tracking algorithm for detecting new faces can operate on the luminance images obtained from the sensor, independent of a normal preview stream obtained from the main array of RGB pixels. A refined analysis can be subsequently performed on the corresponding RGB image.


A particular advantage of this configuration arises where low-lighting conditions occur. As noted above the W pixels charge at a rate at least three times faster than the conventional RGB pixels and thus they can obtain better images in low light conditions. This in turn leads to an extension of the lower lighting levels at which a face tracker algorithm can successfully operate and enables higher frame rates to be achieved in extreme conditions.


Now referring to FIGS. 10-11, further embodiments are provided of enhanced methods, hardware and software involving color sensor arrays such as RGBW sensor arrays that include separate color RGB pixels and white W pixels or luminance pixels. In particular, the W pixels tend to saturate or at least are more greatly exposed than the RGB pixels. The RGB pixels receive on average ⅓ of the light intensity of the W pixels, and thus tend to be underexposed or at least less exposed than the W pixels on average by a third and at most more than ⅓. In accordance with certain embodiments, the RGB pixels may be clocked independently of the W pixels and thus can have a longer exposure time.



FIG. 10 illustrates a method in accordance with certain embodiments. Color pixels of a color sensor array are digital exposed for a first digital exposure duration at (2), while white pixels of the array are exposed at (4) for a second digital exposure time shorter than the first digital exposure duration. The second digital exposure time may be half or less than the first digital exposure time, e.g., a third of the first digital exposure time. The exposure times for the red R, blue B, and green G pixels may be the same as or different from each other, for example, in accordance with certain determined color enhancements of the image. The exposure times for any of the red R, blue B, green G and/or white W pixels may vary over the multiple pixels of the sensor array according to spatially-dependent color enhancements, contrast enhancements, and/or intensity enhancements. The pixels may be further adjusted based on other digital image enhancement methods as may be described in other references incorporated by reference herein.


The different digital exposure durations may be provided either by gathering light for different exposure durations, and/or by using data gathered for different periods of time, i.e., while not using other data gathered at times outside a determined period particularly for the white W pixels, and/or by using computed data based on gathered data such as by multiplying or otherwise increasing color pixel data such as by adding data from surrounding pixels in a stand alone approach or in tandem with a dedicated or unrelated smoothing process. Alternatively, a color sensor may be provided that includes a ratio of more or less color pixels to white pixels such as to contribute to intensity balancing and/or certain pixels may be turned off or data gathered from certain pixels may simply not be used in pre-balancing the intensities, or may be discarded in a post-balancing filter. Certain pixels may be blocked, reduced or filtered mechanically, optically, electronically, or digitally such as in accordance with program code, and certain pixels may be likewise enhanced as to the extent of their respective contributions.


The color pixels and the white pixels make be clocked independently or differently such as to provide the different exposure durations. The color pixels and white pixels may be clocked over different temporal ranges which are proximate in time and may or may not overlap. In one example, a temporal range of exposure of the white pixels may lie entirely or partially within a longer temporal range of exposure of the color pixels or a subset thereof, or a temporal range of exposure of the white pixels may just precede or just follow a longer temporal range of exposure of the color pixels or subset thereof.


The color pixels will generally exhibit a greater motion blurring effect than the white pixels due to the color pixels being digitally-exposed for a longer duration than the white pixels. Thus the method of FIG. 10 may include at (5) compensating blurring in the color pixels, for example, by using less-blurred data from the white pixels.


A color filter enhanced digital image is generated at (6) using data from both the color pixels digitally or otherwise exposed for the first digital exposure duration and the white pixels digitally or otherwise exposed for the second digital exposure duration. The color filter enhanced digital image or a further processed version is stored, transmitted, communicated, displayed, and/or projected.



FIG. 11 illustrates in block form a portable digital imaging apparatus such as a portable camera or camera phone. The apparatus may include optional display 12 and flash 14, as well as optics 16 and digital camera sensor 18. Alternatively, methods according to certain embodiments may be applied to film images that are digitized and processed accordingly. The color sensor includes several white pixels W, as well as red pixels R, blur pixels B and red pixels R, although not necessarily in the ratios or positions indicated in the example sensor (18) of FIG. 11. The color sensor array (18) may be or may include a CMOS-based sensor or sensor elements.


The portable digital imaging apparatus include a processor 20 and memory 22 including program code for image processing and handling of image acquisition. The digital image acquisition duration for color pixels RGB is increased compared with that for white pixels W in the sensor array 18, or that for white pixels W is decreased compared to color pixels RGB. There are multiple ways that the apparatus of FIG. 11 can achieve this. In one aspect, an image is exposed to the sensor for the same time period for both color pixels RGB and white pixels W, and then data gathered by the white pixels for only a fraction of that duration is used along with the color pixel data for the entire duration or at least a longer duration than that of the white pixels W. The apparatus may be programmed such that the ratio of the durations used may be proportional to the sensor configuration, such as to depend on the ratio white pixels W to color pixels RGB. For example, the array 18 shows a 1:1 ration of white pixels W to color pixels RGB, but that ratio could be 1:2 or 1:3 or even 2:1 or otherwise. The white pixels W may be the same or different size as the color pixels RGB.


In another aspect, clocking component 24, which can include a single clock or multiple clocks, turns on the white pixels for gathering image data for a different, shorter period than it does for the color pixels. For example, the clocking component 24 may trigger the beginning of gathering data by both the color pixels RGB and the white pixels W at time t0. Then, the clocking component may trigger the stopping of gathering data by the white pixels at time t1. The color pixels RGB would continue to gather data after time t1. Then, the clocking component 24 would trigger the stopping of gathering data by the color pixels RGB at time t2. Alternatively, the clocking component can trigger the color pixels RGB and the white pixels W to begin gathering data at different times t01 and t02, and may stop the gathering of data by the white pixels W and color pixels RGB at the same time or such that the color pixels RGB are stopped before the white pixels W, and the durations may overlap entirely, partially overlap or not overlap at all. In the latter case, it is desired that the later period of the two data gathering periods closely follows the earlier period, so that the color pixels gather data of substantially the same scene as the white pixels.


The clocking component 24 could be used in the earlier aspect too. The data would be gathered by the white pixels W and color pixels RGB for the same duration between time t0 and t1 according to the clocking component. However, the clocking component would tag other times t2, t3, and perhaps times, between time t0 and t1. Then, only data gathered by the white pixels between times t2 and t3 would be used, while the data gathered by the color pixels RGB between times t0 and t1, or some other duration longer than t3 minus t2, would be used.


In a further aspect, the sensor design may be modified so that the acquisition time can be independently increased for color pixels RGB compared to white pixels W. White pixels W in this aspect would be exposed during a normal period of time (or slightly under exposed) depending on the given scene condition. As color pixels RGB use only ⅓ of light energy comparing to white pixels W, the exposure time for the color pixels RGB would be increased, e.g., by 3 times, without suffering significant over exposure assuming that white pixels W were not overexposed and that the color distribution is approximate equal for the red R, blue B and green G components of the scene.


A factor is that there may be movement of scene objects or sensor shake, jitter or panning. Then, the RGB pixels RGB would be more motion-blurred, e.g. by 3 times, compared with the white pixels W, because the exposure duration of the RGB pixels is so much longer. For this, the RGB blurring can be compensated by a software component that will use the correction from white pixels W. A process described at U.S. application Ser. No. 11/856,721, which is assigned to the same assignee as the present application and is incorporated by reference, may be used to enable the white pixels W in such a sensor to act as a reference image, or one or more dedicated reference images may be used, where the color pixels RGB provide a blurred color image. A single sharp image can be provided by such an array. In the context of adjusting images with motion blur effects, U.S. patent application Ser. Nos. 11/566,180, 10/985,657, 10/985,650, 10/986,562, 11/753,098, 12/137,113, 11/690,836, 11/856,721, 11/764,578, PCT/US2008/061328, and 61/023,774, which are assigned to the same assignee as the present application are also incorporated by reference.


This process may be performed and/or provided according to a dynamic range compression algorithm or similar. This may allow enhancement of dynamic range of an image as two sub-images are acquired with different exposures. Different acquisition durations of white pixels W and color pixels RGB is possible with CMOS based sensors.


Embodiments of the invention enable a full exposure range of a RGBW sensor array to be used. While exemplary drawings and specific embodiments of the present invention have been described and illustrated, it is to be understood that that the scope of the present invention is not to be limited to the particular embodiments discussed. Thus, the embodiments shall be regarded as illustrative rather than restrictive, and it should be understood that variations may be made in those embodiments by workers skilled in the arts without departing from the scope of the present invention as set forth in the appended claims, and structural and functional equivalents thereof.


In addition, in methods that may be performed according to preferred embodiments herein and that may have been described above, the operations have been described in selected typographical sequences. However, the sequences have been selected and so ordered for typographical convenience and are not intended to imply any particular order for performing the operations, except for those where a particular order may be expressly set forth or where those of ordinary skill in the art may deem a particular order to be necessary.


In addition, all references cited above herein, as well as the background, invention summary, abstract and brief description of the drawings, are all incorporated by reference into the Detailed Descriptions of the Embodiments as disclosing alternative embodiments.

Claims
  • 1. A color filter enhancement method for a portable digital image acquisition device, comprising: using optics, a color sensor array and a processor to acquire and process digital images;obtaining a first relatively underexposed and sharp image of a scene (“sharp image”) by exposing a first set of pixels of the color sensor array;obtaining a second relatively well exposed and blurred image of the same scene (“blurred image”) by exposing a second set of pixels of the color sensor array for a longer duration than the sharp image; andapplying color information from the blurred image to the sharp image, thereby constructing a color-enhanced version of the sharp image;wherein applying the color information from the blurred image to the sharp image includes changing one or more pixels of the sharp image based on the color information of the blurred image.
  • 2. The method of claim 1, wherein obtaining the blurred image comprises digitally exposing color pixels of the color sensor array for a first digital exposure duration; and wherein obtaining the sharp image comprises digitally exposing white pixels of the color sensor array for a second digital exposure duration shorter than the first digital exposure duration.
  • 3. The method of claim 2, wherein the constructing of the color-enhanced version of the sharp image comprises using data from both the color pixels exposed for the first digital exposure duration and the white pixels exposed for the second digital exposure duration.
  • 4. The method of claim 2, wherein the digitally exposing of the color pixels and the white pixels for different exposure times comprises clocking the color pixels and the white pixels independently.
  • 5. The method of claim 2, wherein the digitally-exposing of the color pixels and the white pixels for different exposure times comprises including sensor data over different temporal ranges.
  • 6. The method of claim 5, wherein the different temporal ranges are overlapping.
  • 7. The method of claim 6, wherein a first temporal range corresponding to the digitally-exposing of the color pixels comprises an entire second temporal range corresponding to the digitally-exposing of the white pixels.
  • 8. The method of claim 2, wherein the color pixels comprise greater motion blurring effect than the white pixels due to the color pixels being digitally-exposed for a longer duration than the white pixels, and the method further comprises compensating blurring in the color pixels using less-blurred data from the white pixels.
  • 9. The method of claim 2, wherein at least a segment of the first digital exposure duration and the second digital exposure duration temporally overlap.
  • 10. The method of claim 2, wherein the second digital exposure duration comprises not more than half of the first digital exposure duration.
  • 11. The method of claim 2, wherein the second digital exposure duration comprises approximately a third of the first digital exposure duration.
  • 12. The method of claim 1, wherein the color sensor array comprises a CMOS-based sensor.
  • 13. The method of claim 1, further comprising storing, transmitting, or displaying the color-enhanced version of the sharp image, or a further processed version, or combinations thereof.
  • 14. A portable digital image acquisition device, comprising: optics and a color sensor array for acquiring a color digital image;a processor;one or more processor-readable media having code embedded therein for programming the processor to perform a color filter enhancement method that comprises the following: obtaining a first relatively underexposed and sharp image of a scene (“sharp image”) by exposing a first set of pixels of the sensor array;obtaining a second relatively well exposed and blurred image of the same scene (“blurred image”) by exposing a second set of pixels of the color sensor array for a longer duration than the sharp image; andapplying color information from the blurred image to the sharp image, thereby constructing a color-enhanced of the sharp image;wherein applying the color information from the blurred image to the sharp image includes changing one or more pixels of the sharp image based on the color information of the blurred image.
  • 15. The device of claim 14, wherein obtaining the blurred image comprises digitally exposing color pixels of the color sensor array for a first digital exposure duration; and wherein obtaining the sharp image comprises digitally exposing white pixels of the color sensor array for a second digital exposure duration shorter than the first digital exposure duration.
  • 16. The device of claim 15, wherein the constructing of the color-enhanced version of the sharp image comprises using data from both the color pixels exposed for the first digital exposure duration and the white pixels exposed for the second digital exposure duration.
  • 17. The device of claim 15, wherein the digitally exposing of the color pixels and the white pixels for different exposure times comprises clocking the color pixels and the white pixels independently.
  • 18. The device of claim 15, wherein the digitally-exposing of the color pixels and the white pixels for different exposure times comprises including sensor data over different temporal ranges.
  • 19. The device of claim 18, wherein the different temporal ranges are overlapping.
  • 20. The device of claim 19, wherein a first temporal range corresponding to the digitally-exposing of the color pixels comprises an entire second temporal range corresponding to the digitally-exposing of the white pixels.
  • 21. The device of claim 15, wherein the color pixels comprise greater motion blurring effect than the white pixels due to the color pixels being digitally-exposed for a longer duration than the white pixels, and the method further comprises compensating blurring in the color pixels using less-blurred data from the white pixels.
  • 22. The device of claim 15, wherein at least a segment of the first digital exposure duration and the second digital exposure duration temporally overlap.
  • 23. The device of claim 15, wherein the second digital exposure duration comprises not more than half of the first digital exposure duration.
  • 24. The device of claim 15, wherein the second digital exposure duration comprises approximately a third of the first digital exposure duration.
  • 25. The device of claim 14, wherein the color sensor array comprises a CMOS-based sensor.
  • 26. The device of claim 14, wherein the method further comprises storing, transmitting, or displaying the color-enhanced version of the sharp image, or a further processed version, or combinations thereof.
  • 27. One or more non-transitory processor-readable media having code embedded therein for programming the processor to perform a color filter enhancement method that comprises: obtaining a first relatively underexposed and sharp image of a scene (“sharp image”) by exposing a first set of pixels of the color sensor array;obtaining a second relatively well exposed and blurred image of the same scene (“blurred image”) by exposing a second set of pixels of the color sensor array for a longer duration than the sharp image; andapplying color information from the blurred image to the sharp image, thereby constructing a color-enhanced version of the sharp image;wherein applying the color information from the blurred image to the sharp image includes changing one or more pixels of the sharp image based on the color information of the blurred image.
  • 28. The one or more non-transitory processor-readable media of claim 27, wherein obtaining the blurred image comprises digitally exposing color pixels of the color sensor array for a first digital exposure duration; and wherein obtaining the sharp image comprises digitally exposing white pixels of the color sensor array for a second digital exposure duration shorter than the first digital exposure duration.
  • 29. The one or more non-transitory processor-readable media of claim 28, wherein the constructing of the color-enhanced version of the sharp image comprises using data from both the color pixels exposed for the first digital exposure duration and the white pixels exposed for the second digital exposure duration.
  • 30. The one or more non-transitory processor-readable media of claim 28, wherein the digitally exposing of the color pixels and the white pixels for different exposure times comprises clocking the color pixels and the white pixels independently.
  • 31. The one or more non-transitory processor-readable media of claim 28, wherein the digitally-exposing of the color pixels and the white pixels for different exposure times comprises including sensor data over different temporal ranges.
  • 32. The one or more non-transitory processor-readable media of claim 31, wherein the different temporal ranges are overlapping.
  • 33. The one or more non-transitory processor-readable media of claim 32, wherein a first temporal range corresponding to the digitally-exposing of the color pixels comprises an entire second temporal range corresponding to the digitally-exposing of the white pixels.
  • 34. The one or more non-transitory processor-readable media of claim 28, wherein the color pixels comprise greater motion blurring effect than the white pixels due to the color pixels being digitally-exposed for a longer duration than the white pixels, and the method further comprises compensating blurring in the color pixels using less-blurred data from the white pixels.
  • 35. The one or more non-transitory processor-readable media of claim 28, wherein at least a segment of the first digital exposure duration and the second digital exposure duration temporally overlap.
  • 36. The one or more non-transitory processor-readable media of claim 28, wherein the second digital exposure duration comprises not more than half of the first digital exposure duration.
  • 37. The one or more non-transitory processor-readable media of claim 28, wherein the second digital exposure duration comprises approximately a third of the first digital exposure duration.
  • 38. The one or more non-transitory processor-readable media of claim 28, wherein the color sensor array comprises a CMOS-based sensor.
  • 39. The one or more non-transitory processor-readable media of claim 28, wherein the method further comprises storing, transmitting, or displaying the color-enhanced version of the sharp image, or a further processed version, or combinations thereof.
  • 40. The method of claim 1, further comprising: applying luminance information from the blurred image to the sharp image, thereby constructing a luminance-enhanced version of the sharp image;wherein applying the luminance information from the blurred image to the sharp image includes changing one or more pixels of the sharp image based on the luminance information of the blurred image.
  • 41. The method of claim 40, wherein: obtaining the blurred image comprises digitally exposing color pixels of the color sensor array for a first digital exposure duration;obtaining the sharp image comprises digitally exposing white pixels of the color sensor array for a second digital exposure duration shorter than the first digital exposure duration; andconstructing the luminance-enhanced version of the sharp image comprises using data from both the color pixels exposed for the first digital exposure duration and the white pixels exposed for the second digital exposure duration.
  • 42. The method of claim 40, further comprising storing, transmitting, or displaying the luminance-enhanced version of the sharp image, or a further processed version, or combinations thereof.
  • 43. The device of claim 14, wherein the method further comprises: applying luminance information from the blurred image to the sharp image, thereby constructing a luminance-enhanced version of the sharp image;wherein applying the luminance information from the blurred image to the sharp image includes changing one or more pixels of the sharp image based on the luminance information of the blurred image.
  • 44. The device of claim 43, wherein: obtaining the blurred image comprises digitally exposing color pixels of the color sensor array for a first digital exposure duration;obtaining the sharp image comprises digitally exposing white pixels of the color sensor array for a second digital exposure duration shorter than the first digital exposure duration; andconstructing of the luminance-enhanced version of the sharp image comprises using data from both the color pixels exposed for the first digital exposure duration and the white pixels exposed for the second digital exposure duration.
  • 45. The device of claim 43, wherein the method further comprises storing, transmitting, or displaying the luminance-enhanced version of the sharp image, or a further processed version, or combinations thereof.
  • 46. The one or more non-transitory processor-readable media of claim 28, wherein the method further comprises: applying luminance information from the blurred image to the sharp image, thereby constructing a luminance-enhanced version of the sharp image;wherein applying the luminance information from the blurred image to the sharp image includes changing one or more pixels of the sharp image based on the luminance information of the blurred image.
  • 47. The one or more non-transitory processor-readable media of claim 46, wherein: obtaining the blurred image comprises digitally exposing color pixels of the color sensor array for a first digital exposure duration;obtaining the sharp image comprises digitally exposing white pixels of the color sensor array for a second digital exposure duration shorter than the first digital exposure duration; andconstructing of the luminance-enhanced version of the sharp image comprises using data from both the color pixels exposed for the first digital exposure duration and the white pixels exposed for the second digital exposure duration.
  • 48. The one or more non-transitory processor-readable media of claim 46, wherein the method further comprises storing, transmitting, or displaying the luminance-enhanced version of the sharp image, or a further processed version, or combinations thereof.
PRIORITY

This application is a Continuation of U.S. patent application Ser. No. 12/330,719, filed Dec. 9, 2008; which is a Continuation-in-Part (CIP) of U.S. patent application Ser. No. 11/856,721, filed Sep. 18, 2007; which claims priority to U.S. provisional patent application No. 60/893,116, filed Mar. 5, 2007. These applications are related to PCT/EP2007/009939, filed Nov. 14, 2007, and published as WO 2009/036793. Each of these applications are incorporated by reference.

US Referenced Citations (208)
Number Name Date Kind
5251019 Moorman et al. Oct 1993 A
5374956 D'luna Dec 1994 A
5392088 Abe et al. Feb 1995 A
5428723 Ainscow et al. Jun 1995 A
5510215 Prince et al. Apr 1996 A
5599766 Boroson et al. Feb 1997 A
5686383 Long et al. Nov 1997 A
5747199 Roberts et al. May 1998 A
5751836 Wildes et al. May 1998 A
5756239 Wake May 1998 A
5756240 Roberts et al. May 1998 A
5802220 Black et al. Sep 1998 A
5889277 Hawkins et al. Mar 1999 A
5889554 Mutze Mar 1999 A
5909242 Kobayashi et al. Jun 1999 A
5981112 Roberts Nov 1999 A
6028960 Graf et al. Feb 2000 A
6035072 Read Mar 2000 A
6041078 Rao Mar 2000 A
6061462 Tostevin et al. May 2000 A
6081606 Hansen et al. Jun 2000 A
6114075 Long et al. Sep 2000 A
6122017 Taubman Sep 2000 A
6124864 Madden et al. Sep 2000 A
6134339 Luo Oct 2000 A
6198505 Turner et al. Mar 2001 B1
6269175 Hanna et al. Jul 2001 B1
6297071 Wake Oct 2001 B1
6297846 Edanami Oct 2001 B1
6326108 Simons Dec 2001 B2
6330029 Hamilton et al. Dec 2001 B1
6360003 Doi et al. Mar 2002 B1
6365304 Simons Apr 2002 B2
6381279 Taubman Apr 2002 B1
6387577 Simons May 2002 B2
6407777 DeLuca Jun 2002 B1
6535244 Lee et al. Mar 2003 B1
6555278 Loveridge et al. Apr 2003 B1
6567536 McNitt et al. May 2003 B2
6599668 Chari et al. Jul 2003 B2
6602656 Shore et al. Aug 2003 B1
6607873 Chari et al. Aug 2003 B2
6618491 Abe Sep 2003 B1
6625396 Sato Sep 2003 B2
6643387 Sethuraman et al. Nov 2003 B1
6741960 Kim et al. May 2004 B2
6863368 Sadasivan et al. Mar 2005 B2
6892029 Tsuchida et al. May 2005 B2
6940545 Ray et al. Sep 2005 B1
6947609 Seeger et al. Sep 2005 B2
6961518 Suzuki Nov 2005 B2
7019331 Winters et al. Mar 2006 B2
7072525 Covell Jul 2006 B1
7084037 Gamo et al. Aug 2006 B2
7160573 Sadasivan et al. Jan 2007 B2
7177538 Sato et al. Feb 2007 B2
7180238 Winters Feb 2007 B2
7195848 Roberts Mar 2007 B2
7269292 Steinberg Sep 2007 B2
7292270 Higurashi et al. Nov 2007 B2
7315324 Cleveland et al. Jan 2008 B2
7315630 Steinberg et al. Jan 2008 B2
7315631 Corcoran et al. Jan 2008 B1
7316630 Tsukada et al. Jan 2008 B2
7316631 Tsunekawa Jan 2008 B2
7317815 Steinberg et al. Jan 2008 B2
7336821 Ciuc et al. Feb 2008 B2
7369712 Steinberg et al. May 2008 B2
7403643 Ianculescu et al. Jul 2008 B2
7453493 Pilu Nov 2008 B2
7453510 Kolehmainen et al. Nov 2008 B2
7460695 Steinberg et al. Dec 2008 B2
7469071 Drimbarean et al. Dec 2008 B2
7489341 Yang et al. Feb 2009 B2
7548256 Pilu Jun 2009 B2
7551755 Steinberg et al. Jun 2009 B1
7565030 Steinberg et al. Jul 2009 B2
7593144 Dymetman Sep 2009 B2
7623153 Hatanaka Nov 2009 B2
7636486 Steinberg et al. Dec 2009 B2
7639888 Steinberg et al. Dec 2009 B2
7639889 Steinberg et al. Dec 2009 B2
7660478 Steinberg et al. Feb 2010 B2
7676108 Steinberg et al. Mar 2010 B2
7692696 Steinberg et al. Apr 2010 B2
7697778 Steinberg et al. Apr 2010 B2
7738015 Steinberg et al. Jun 2010 B2
7912285 Steinberg et al. Mar 2011 B2
8264576 Susanu et al. Sep 2012 B2
8351726 Albu et al. Jan 2013 B2
20010036307 Hanna et al. Nov 2001 A1
20020006163 Hibi et al. Jan 2002 A1
20030032448 Bulthuis et al. Feb 2003 A1
20030052991 Stavely et al. Mar 2003 A1
20030058361 Yang Mar 2003 A1
20030091225 Chen May 2003 A1
20030103076 Neuman Jun 2003 A1
20030151674 Lin Aug 2003 A1
20030152271 Tsujino et al. Aug 2003 A1
20030169818 Obrador Sep 2003 A1
20030193699 Tay Oct 2003 A1
20030219172 Caviedes et al. Nov 2003 A1
20040066981 Li et al. Apr 2004 A1
20040076335 Kim Apr 2004 A1
20040090532 Imada May 2004 A1
20040120598 Feng Jun 2004 A1
20040120698 Hunter Jun 2004 A1
20040130628 Stavely Jul 2004 A1
20040145659 Someya et al. Jul 2004 A1
20040169767 Norita et al. Sep 2004 A1
20040212699 Molgaard Oct 2004 A1
20040218057 Yost et al. Nov 2004 A1
20040218067 Chen et al. Nov 2004 A1
20040247179 Miwa et al. Dec 2004 A1
20050010108 Rahn et al. Jan 2005 A1
20050019000 Lim et al. Jan 2005 A1
20050031224 Prilutsky et al. Feb 2005 A1
20050041121 Steinberg et al. Feb 2005 A1
20050041123 Ansari et al. Feb 2005 A1
20050047672 Ben-Ezra et al. Mar 2005 A1
20050052553 Kido et al. Mar 2005 A1
20050057687 Irani et al. Mar 2005 A1
20050068446 Steinberg et al. Mar 2005 A1
20050068452 Steinberg et al. Mar 2005 A1
20050140801 Prilutsky et al. Jun 2005 A1
20050140829 Uchida et al. Jun 2005 A1
20050146622 Silverstein Jul 2005 A9
20050201637 Schuler et al. Sep 2005 A1
20050219391 Sun et al. Oct 2005 A1
20050231625 Parulski et al. Oct 2005 A1
20050248660 Stavely et al. Nov 2005 A1
20050259864 Dickinson et al. Nov 2005 A1
20050270381 Owens et al. Dec 2005 A1
20050281477 Shiraki et al. Dec 2005 A1
20060006309 Dimsdale et al. Jan 2006 A1
20060017837 Sorek et al. Jan 2006 A1
20060038891 Okutomi et al. Feb 2006 A1
20060039690 Steinberg et al. Feb 2006 A1
20060093212 Steinberg et al. May 2006 A1
20060098237 Steinberg et al. May 2006 A1
20060098890 Steinberg et al. May 2006 A1
20060098891 Steinberg et al. May 2006 A1
20060119710 Ben-Ezra et al. Jun 2006 A1
20060120599 Steinberg et al. Jun 2006 A1
20060125938 Ben-Ezra et al. Jun 2006 A1
20060133688 Kang et al. Jun 2006 A1
20060140455 Costache et al. Jun 2006 A1
20060170786 Won Aug 2006 A1
20060171464 Ha Aug 2006 A1
20060187308 Lim et al. Aug 2006 A1
20060204034 Steinberg et al. Sep 2006 A1
20060204054 Steinberg et al. Sep 2006 A1
20060204110 Steinberg et al. Sep 2006 A1
20060279639 Silverstein et al. Dec 2006 A1
20060285754 Steinberg et al. Dec 2006 A1
20070009169 Bhattacharjya Jan 2007 A1
20070025714 Shiraki Feb 2007 A1
20070058073 Steinberg et al. Mar 2007 A1
20070083114 Yang et al. Apr 2007 A1
20070086675 Chinen et al. Apr 2007 A1
20070097221 Stavely et al. May 2007 A1
20070110305 Corcoran et al. May 2007 A1
20070147820 Steinberg et al. Jun 2007 A1
20070171987 Trimeche Jul 2007 A1
20070189748 Drimbarean et al. Aug 2007 A1
20070201724 Steinberg et al. Aug 2007 A1
20070234779 Hsu et al. Oct 2007 A1
20070268388 Watanabe et al. Nov 2007 A1
20070269108 Steinberg et al. Nov 2007 A1
20070296833 Corcoran et al. Dec 2007 A1
20080012969 Kasai et al. Jan 2008 A1
20080037827 Corcoran et al. Feb 2008 A1
20080037839 Corcoran et al. Feb 2008 A1
20080037840 Steinberg et al. Feb 2008 A1
20080043121 Prilutsky et al. Feb 2008 A1
20080175481 Petrescu et al. Jul 2008 A1
20080211943 Egawa et al. Sep 2008 A1
20080218611 Parulski et al. Sep 2008 A1
20080219581 Albu et al. Sep 2008 A1
20080219585 Kasai et al. Sep 2008 A1
20080220750 Steinberg et al. Sep 2008 A1
20080231713 Florea et al. Sep 2008 A1
20080232711 Prilutsky et al. Sep 2008 A1
20080240555 Nanu et al. Oct 2008 A1
20080240607 Sun et al. Oct 2008 A1
20080292193 Bigioi et al. Nov 2008 A1
20080309769 Albu et al. Dec 2008 A1
20080309770 Florea et al. Dec 2008 A1
20090003652 Steinberg et al. Jan 2009 A1
20090009612 Tico et al. Jan 2009 A1
20090080713 Bigioi et al. Mar 2009 A1
20090080796 Capata et al. Mar 2009 A1
20090080797 Nanu et al. Mar 2009 A1
20090167893 Susanu et al. Jul 2009 A1
20090179999 Albu et al. Jul 2009 A1
20090185041 Kang et al. Jul 2009 A1
20090185753 Albu et al. Jul 2009 A1
20090190803 Neghina et al. Jul 2009 A1
20090196466 Capata et al. Aug 2009 A1
20090303342 Corcoran et al. Dec 2009 A1
20090303343 Drimbarean et al. Dec 2009 A1
20100053349 Watanabe et al. Mar 2010 A1
20100126831 Ceelen May 2010 A1
20100201826 Steinberg et al. Aug 2010 A1
20100201827 Steinberg et al. Aug 2010 A1
20100302393 Olsson et al. Dec 2010 A1
20110090352 Wang et al. Apr 2011 A1
20110102642 Wang et al. May 2011 A1
Foreign Referenced Citations (32)
Number Date Country
3729324 Mar 1989 DE
10154203 Jun 2002 DE
10107004 Sep 2002 DE
0 944 251 Sep 1999 EP
944251 Apr 2003 EP
1429290 Jun 2004 EP
1583033 Oct 2005 EP
1779322 Jan 2008 EP
1429290 Jul 2008 EP
10285542 Oct 1998 JP
11327024 Nov 1999 JP
2005-352721 Dec 2005 JP
2008-520117 Jun 2008 JP
2010211485 Sep 2010 JP
WO9843436 Oct 1998 WO
WO 0113171 Feb 2001 WO
WO0113171 Feb 2001 WO
WO0245003 Jun 2002 WO
WO03071484 Aug 2003 WO
WO04001667 Dec 2003 WO
WO2004036378 Apr 2004 WO
WO 2006050782 May 2006 WO
WO2006050782 May 2006 WO
WO2007093199 Aug 2007 WO
WO2007093199 Aug 2007 WO
WO2007142621 Dec 2007 WO
WO2007143415 Dec 2007 WO
WO2008017343 Feb 2008 WO
WO2007143415 May 2008 WO
WO2008131438 Oct 2008 WO
WO 2009036793 Mar 2009 WO
WO2010066381 Jun 2010 WO
Non-Patent Literature Citations (159)
Entry
Final Office Action mailed May 27, 2010, for U.S. Appl. No. 11/753,098, filed May 24, 2007.
Non-Final Office Action mailed Jan. 20, 2010, for U.S. Appl. No. 11/690,836, filed Mar. 25, 2007.
Non-Final Office Action mailed Nov. 4, 2009, for U.S. Appl. No. 11/753,098, filed May 24, 2007.
Notice of Allowance mailed Apr. 29, 2010, for U.S. Appl. No. 11/690,836, filed Mar. 25, 2007.
Notice of Allowance mailed Dec. 1, 2009, for U.S. Appl. No. 12/199,680, filed Aug. 27, 2008.
Notice of Allowance mailed Dec. 11, 2009, for U.S. Appl. No. 12/199,710, filed Aug. 27, 2008.
Notice of Allowance mailed Nov. 2, 2009, for U.S. Appl. No. 10/985,657, filed Nov. 10, 2004.
Notice of Allowance mailed Nov. 17, 2009, for U.S. Appl. No. 11/566,180, filed Dec. 1, 2006.
Notice of Allowance mailed Oct. 29, 2009, for U.S. Appl. No. 10/985,650, filed Nov. 10, 2004.
Notice of Allowance mailed Oct. 29, 2009, for U.S. Appl. No. 10/986,562, filed Nov. 10, 2004.
PCT Notification of Transmittal of the International Search Report and the Written Opinion of the International Searching Authority, or the Declaration for PCT application No. PCT/EP2009/008674, dated Mar. 29, 2010, 10 pages.
Eric P. Bennett, Leonard McMillan, Video Enhancement Using Per-Pixel Virtual Exposures, International Conference on Computer Graphics and Interactive Techniques, ACM Siggraph 2005, 2005, pp. 845-852, XP002477673, New York, NY, USA.
Chen-Kuei Y et al., Color image sharpening by moment-preserving technique, Signal Processing, Elsevier Science Publishers B.V. Amsterdam, NL, vol. 45, No. 3, Sep. 1995, pp. 397-403, XP004002600, ISSN: 0165-1684.
Oppenheim A V et al., The Importance of Phase in Signals, Proceedings of the IEEE, IEEE. New York, US, vol. 69, No. 5, May 1981 pp. 529-541, XP008060042, ISSN: 0018-9219.
Michael Elad et al., Superresolution Restoration of an Image Sequence: Adaptive Filtering Approach, IEEE Transactions on Image Processing, IEEE Service Center, Piscataway, NJ, US, vol. 8, No. 3, Mar. 1999, XP011026284, ISSN: 1057-7149.
Fujita K et al., An edge-adaptive iterative method for image restoration, Singapore ICCS/ISITA '92. Communications on the Move, Singapore 16-20, Nov. 1992, New York, NY, USA, IEEE, US, Nov. 16, 1992, pp. 361-365, XP010066997, ISBN: 0-7803-0803-4.
Notification of Transmittal of the International Search Report and the Written Opinion of the International Searching Authority, or the Declaration, for PCT/EP2007/009939, dated May 21, 2008, 13 pages.
Co-pending U.S. Appl. No. 10/764,339, filed Jan. 22, 2004.
Co-pending U.S. Appl. No. 10/985,650, filed Nov. 10, 2004.
Co-pending U.S. Appl. No. 11/319,766, filed Dec. 27, 2005.
Co-pending U.S. Appl. No. 11/421,027, filed May 30, 2006.
Co-pending U.S. Appl. No. 11/673,560, filed Feb. 10, 2007.
Co-pending U.S. Appl. No. 11/744,020, filed May 3, 2007.
Co-pending U.S. Appl. No. 11/752,925, filed May 24, 2007.
Co-pending U.S. Appl. No. 11/753,098, filed May 24, 2007.
Co-pending U.S. Appl. No. 11/856,721, filed Sep. 18, 2007.
Co-pending U.S. Appl. No. 11/861,257, filed Sep. 25, 2007.
Co-pending U.S. Appl. No. 11/861,854, filed Jun. 26, 2007.
Co-pending U.S. Appl. No. 12/042,104, filed Mar. 4, 2008.
Co-pending U.S. Appl. No. 12/042,335, filed Mar. 5, 2008.
Co-pending U.S. Appl. No. 12/137,113, filed Jun. 11, 2008.
Co-pending U.S. Appl. No. 12/336,416, filed Dec. 16, 2008.
Non-Final Office Action mailed Apr. 8, 2011, for U.S. Appl. No. 12/330,719, filed Dec. 9, 2008.
Final Office Action mailed Jul. 8, 2011, for U.S. Appl. No. 12/137,113, filed Jun. 11, 2008.
Final Office Action mailed Jul. 18, 2011, for U.S. Appl. No. 12/755,338, filed Apr. 6, 2010.
Final Office Action mailed Apr. 25, 2011, for U.S. Appl. No. 11/856,721, filed Sep. 18, 2007.
Final Office Action mailed Mar. 25, 2011, for U.S. Appl. No. 11/764,578, filed Jun. 18, 2007.
Non-Final Office Action mailed Mar. 21, 2011, for U.S. Appl. No. 12/702,092, filed Feb. 8, 2010.
Non-Final Office Action mailed Apr. 8, 2011, for U.S. Appl. No. 12/901,577, filed Oct. 11, 2010.
Non-Final Office Action mailed Dec. 7, 2010, for U.S. Appl. No. 11/764,578, filed Jun. 18, 2007.
Non-Final Office Action mailed Jan. 21, 2011, for U.S. Appl. No. 12/137,113, filed Jun. 11, 2008.
PCT International Preliminary Report on Patentability, Chapter I, for PCT Application No. PCT/EP2009/008674, dated Jun. 14, 2011, 7 pages.
PCT International Preliminary Report on Patentability, Chapter I (IB/373), for PCT Application No. PCT/EP2008/004729, dated Dec. 17, 2009, 10 pages.
PCT Notification of the Transmittal of the International Preliminary Report on Patentability, Chapter II, for PCT Application No. PCT/EP2007/009939, dated Nov. 16, 2009, 12 pages.
PCT Notification of the Transmittal of the International Search Report and the Written Opinion of the International Searching Authority, or the Declaration, for PCT Application No. PCT/EP2007/009939, dated May 21, 2008, 13 pages.
PCT Notification of the Transmittal of the International Search Report and the Written Opinion of the International Searching Authority, or the Declaration, for PCT Application No. PCT/EP2008/004729, dated Oct. 29, 2008, 14 pages.
PCT Notification of the Transmittal of the International Search Report and the Written Opinion of the International Searching Authority, or the Declaration, for PCT Application No. PCT/US2007/069638, dated Mar. 5, 2008, 7 pages.
PCT Notification of the Transmittal of the International Search Report and the Written Opinion of the International Searching Authority, or the Declaration, for PCT Application No. PCT/EP2005/011011, dated Apr. 24, 2006, 12 pages.
PCT Notification of the Transmittal of the International Search Report and the Written Opinion of the International Searching Authority, or the Declaration, for PCT Application No. PCT/EP2010/056999, dated Sep. 1, 2010, 10 pages.
PCT Written Opinion of the International Search Authority, for PCT Application No. PCT/EP2008/004729, dated Dec. 17, 2009, 10 pages.
Notice of Allowance, dated Feb. 11, 2013, for U.S. Appl. No. 11/856,721, filed Sep. 18, 2007.
Notice of Allowance, dated Jul. 3, 2012, for U.S. Appl. No. 12/330,719, filed Dec. 9, 2008.
Patent Abstracts of Japan, publication No. 2005-352721, publication date: Dec. 22, 2005, for: Imaging Device.
Wu et al., “Blind Image Blur Identification in Cepstrum Domain”, IEEE Proceedings of 16th International Conference on Computer Communications and Networks, 2007, pp. 1166-1171.
Patent Abstracts of Japan, publication No. 2001-346093, publication date: Dec. 14, 2001, for Blurred Image Correction Device, Blurred Image Correction Method, and Recording Medium for Recording Blurred Image Correction Program.
Patent Abstracts of Japan, publication No. 2004-304380, publication date: Oct. 28, 2004, for Photographing Apparatus.
Patent Abstracts of Japan, publication No. 2004-289709, publication date: Oct. 14, 2004, for Image Pickup Device, and Image Pickup Method.
Patent Abstracts of Japan, publication No. 2004-221506, publication date: Aug. 5, 2004, for Solid-State Image Sensing Element and Its Manufacturing Method.
Demir B et al: “Block motion estimation using adaptive modified two-bit transform” Jun. 5, 2007, vo1.1, No. 2, pp. 215-222.
Andrews, H.C. et al., “Digital Image Restoration”, Prentice Hall, 1977.
Bates et al., “Some Implications of Zero Sheets for Blind Deconvolution and Phase Retrieval”, J. Optical Soc. Am. A, 1990, pp. 468-479, vol. 7.
Ben Ezra, Moshe et al., “Motion Deblurring Using Hybrid Imaging”, Proceedings IEEE Computer Society Conference on Computer Vision and Pattern Recognition, 2003.
Ben-Ezra, M. el al., “Motion-Based Motion Deblurring”, IEEE Transactions on Pattern Analysis and Machine Intelligence, 2004, pp. 689-698, vol. 26—Issue 6.
Bhaskaran, V. et al., “Motion estimation using a computation-constrained criterion”, Digital Signal Processing Proceedings , 1997, pp. 229-232, vol. 1.
Bones et al., “Deconvolution and Phase Retrieval With Use of Zero Sheets”, J. Optical Soc. Am. A, 1995, pp. 1,842-1,857, vol. 12.
Cannon M., “Blind Deconvolution of Spatially Invariant Image Blurs with Phase,” IEEE Transactions on Acoustics, Speech, and Signal Processing, 1976, vol. ASSP-24, No. 1.
Caron J. N. et al., “Blind data restoration with an extracted filter function,” Optical Letters, 2001, 26(15), 1164-1166.
Co-pending U.S. Appl. No. 12/026,484, filed Feb. 5, 2008.
Co-pending U.S. Appl. No. 12/116,140, filed May 6, 2008.
Co-pending U.S. Appl. No. 12/354,707, filed Jan. 15, 2009.
Co-pending U.S. Appl. No. 12/755,338, filed Apr. 6, 2010.
Co-pending U.S. Appl. No. 12/820,002, filed Jun. 21, 2010.
Co-pending U.S. Appl. No. 12/820,034, filed Jun. 21, 2010.
Co-pending U.S. Appl. No. 12/820,086, filed Jun. 21, 2010.
Co-pending U.S. Appl. No. 12/907,921, filed Oct. 19, 2010.
Co-pending U.S. Appl. No. 11/573,713 US-2009-0273685 A1 Nov. 5, 2009 Mihai Ciuc et al.
Co-pending U.S. Appl. No. 12/063,089 US-2009-0263022 A1 Oct. 22, 2009 Petrescu et al.
Co-pending U.S. Appl. No. 12/330,719 US-2009-0167893 A1 Jul. 2, 2009 Susanu et al.
Crowley, J. et al., “Multi-modal tracking of faces for video communication, http://citeseer.ist.psu.edu/crowley97multimodal.html”, In Computer Vision and Patent Recognition, 1997.
Dalong Li et al., “Blur identification based on kurtosis minimization,” ICIP , 2005, 905-908.
Deepak Turaga, Mohamed Alkanhal, Search Algorithms for Block-Matching in Motion Estimation, Mid-Term project, 18-899, Spring, 1998, retrieved on Oct. 21, 2010 from the Internet, URL: http://www.ece.cmu.edu/˜ee899/project/deepak—mid.htm, 14 pages.
Deever, A., “In-camera all-digital video stabilization”, Proceedings of the International Conference on Decision Support Systems.Proceedings of ISDSS, 2006, pp. 190-193.
Deller J. et al., “Discrete-Time Processing of Speech Signals,” 1999, 2nd. Edition, Wiley-IEEE Press.
Deller, John R. Jr et al., “Discrete-Time Processing of Speech Signals”, 1993, 908 pages, IEEE Press Classic Reissue (Hardcover).
Demir, B. et al., “Block motion estimation using adaptive modified two-bit transform”, 2007, pp. 215-222, vol. 1—Isuue 2.
Deng G. et al., “The study of logarithmic image processing model and its application to image enhancement,” IEEE Trans. on Image Processing, 1995, pp. 506-512, vol. 4.
Dufournaud et al., “Matching Images With Different Resolutions”, IEEE Conference Proceedings on Computer Vision and Pattern Recognition, 2000.
Elad et al., “Restoration of a Single Superresolution Image from several Blurred, Noisy and Undersampled Measured Images”, IEEE Transactions on Image Processing, 1997, vol. 6—Issue 12.
Elad, Michael et al., “Superresolution Restoration of an Image Sequence: Adaptive Filtering Approach”, IEEE Transactions on Image Processing, 1999, pp. 529-541, vol. 8—Issue 3.
Favaro, Paolo, “Depth from focus/defocus, http://homepages.inf.ed.ac.uk/rbf/Cvonline/LOCAL—COPIES/FAVARO1/dfdtutorial.html.”, 2002.
Feng, J. et al., “Adaptive block matching motion estimation algorithm using bit plane matching”, ICIP, 1995, pp. 496-499.
Fergus R. et al., “Removing Camera Shake from a Single Photograph,” ACM Trans. on Graphics,Proc. SIGGRAPH , 2006, 25 (3), 787-794.
Final Office Action mailed Nov. 5, 2008, for U.S. Appl. No. 10/986,562, filed Nov. 10, 2004.
Final Office Action mailed Nov. 13, 2008, for U.S. Appl. No. 10/985,657, filed Nov. 10, 2004.
Final Office Action mailed Nov. 13, 2008, for U.S. Appl. No. 11/566,180, filed Dec. 1, 2006.
Final Office Action mailed Nov. 24, 2008, for U.S. Appl. No. 10/985,650, filed Nov. 10, 2004.
Fujita K. et al., An edge-adaptive iterative method for image restoration, Singapore ICCS/ISITA '92.“Communications on the Move” Singapore 16-20, Nov. 1992, New York, NY, USA, IEEE, US, Nov. 16, 1992,pp. 361-365, XP010066997, ISBN: 0-7803-0803-4.
Furht Borko et al., “Motion Estimation Algorithms for Video Compression Series,” The Springer International Series in Engineering and Computer Science, 1997, 379, 176.
George-Othon Glentis, “An efficient affine projection algorithm for 2-D FIR adaptive filtering and linear prediction,” Signal Processing, 2006, 86 (1), 98-116.
Golub G. H. et al., “Matrix Computations,” 1996, 3rd edition, John Hopkins University Press, Baltimore.
Gunturk et al., “High-Resolution Image Reconstruction from Multiple Differently Exposed Images,” IEEE Signal Processing Letters, 2006, vol. 13, No. 4.
Hayes M., “Statistical Digital Signal Processing and Modeling,” 1996, Wiley.
Haykin S., “Adaptive filter theory,” 1996, Prentice Hall.
James Miskin et al., “Ensemble Learning for Blind Image Separation and Deconvolution,” Adv. in Independent Component Analysis, 2000, M. Girolani,Springer-Verlag.
Jannson, Peter A., “Deconvolution of Images and Spectra,” 1997, 2nd. Edition, Academic Press.
Jansson, Peter A., “Chapter 1: Convolution and Related Concepts”, Deconvolution of Images and Spectra, 1997, 2nd. Edition, Academic Press.
Jourlin M. et al., “Logarithmic image processing the mathematical and physical framework for the representation and processing of transmitted images,” Advances in Imaging and Electron Physics, 2001, pp. 129-196, vol. 115.
Ko, S. et al., “Fast digital image stabilizer based on gray-coded bit-plane matching”, IEEE Transactions on Consumer Electronics, 1999, pp. 598-603, vol. 45—Issue 3.
Kuglin C. D. et al., “The phase correlation image alignment method,” Proc. Int. Conf. Cybernetics and Society, 1975, pp. 163-165, IEEE, Bucharest, Romania.
Kundur D. et al., “Blind Image Deconvolution Revisited,” IEEE Signal Processing Magazine, 1996, 13 (6), 61-63.
Lagendijk R. L. et al., “Iterative Identification and Restoration of Images,” 1991, Kluwer Academic.
Lane et al., “Automatic multidimensional deconvolution”, J. Opt. Soc. Am. A, 1987, pp. 180-188, vol. 4—Issue 1.
Lhuillier, M. et al., A quasi-dense approach to surface reconstruction from uncalibrated images, IEEE Transactions on Pattern Analysis and Machine Intelligence, 2005, pp. 418-433, vol. 27—Issue 3, IEEE Comput. Soc. http://rlinks2.dialog.com/NASApp/ChannelWEB/DialogProServlet?ChName=engineering.
Mahmood R. et al., “Two-Dimensional Block Diagonal LMS Adaptive Filtering,” IEEE Transaction on Signal Processing, 1994, 42 (9), 2420-2429.
Mase, Mitsuhito et al., “A Wide Dynamic Range CMOS tmage Sensor with Multiple Exposure-Time Signal Outputs and 12-bit Column-Parallel Cyclic A/D Converters”, IEEE Journal of Solid-State Circuits, 2005, vol. 40—Issue 12.
Natarajan B. et al., “Low-complexity block-based motion estimation via one-bit transforms”, IEEE Trans. Circuit Syst. Video Technol, 1997, pp. 702-706, vol. 7—Issue 5.
Nie Y. et al., “Adaptive rood pattern search for fast block-matching motion estimation,” IEEE Trans. Image Processing, 2002, 11(12).
Non-Final Office Action mailed Apr. 3, 2008, for U.S. Appl. No. 10/985,650, filed Nov. 10, 2004.
Non-Final Office Action mailed Apr. 4, 2008, for U.S. Appl. No. 10/986,562, filed Nov. 10, 2004.
Non-Final Office Action mailed Apr. 21, 2009, for U.S. Appl. No. 10/985,650, filed Nov. 10, 2004.
Non-Final Office Action mailed Apr. 22, 2009, for U.S. Appl. No. 10/986,562, filed Nov. 10, 2004.
Non-Final Office Action mailed Jun. 1, 2009, for U.S. Appl. No. 10/985,657, filed Nov. 10, 2004.
Non-Final Office Action mailed Jun. 10, 2009, for U.S. Appl. No. 11/566,180, filed Dec. 1, 2006.
Non-Final Office Action mailed Mar. 18, 2008, for U.S. Appl. No. 10/985,657, filed Nov. 10, 2004.
Non-Final Office Action mailed Mar. 21, 2008, for U.S. Appl. No. 11/566,180, filed Dec. 1, 2006.
Non-Final Office Action mailed May 11, 2009, for U.S. Appl. No. 12/199,680, filed Aug. 27, 2008.
Non-Final Office Action mailed May 29, 2009, for U.S. Appl. No. 12/199,710, filed Aug. 27, 2008.
Oppenheim, A.V. et al., “The Importance of Phase in Signals, XP008060042, ISSN: 0018-9219.”, Proceedings of the IEEE, 1981, pp. 529-541, vol. 69—Issue 5.
Park, Sung Cheol et al., “Super-resolution image reconstruction: a technical overview, ISSN: 1053-5888. DOI: 10.1109/MSP.2003.1203207.”, Signal Processing Magazine, 2003, pp. 21-36, vol. 20—Issue 3, IEEE Publication.
Patti A. et al., “Super-Resolution video reconstruction with arbitrary sampling lattices and non-zero aperture time http://citeseer.ist.psu.edu/patti97super.html”, In IEEE Transactions on Image Processing, 1997, pp. 1064-1078.
PCT International Preliminary Report on Patentability, for PCT Application No. PCT/EP2005/011011, dated Jan. 22, 2007, 8 pages.
PCT International Preliminary Report on Patentability, for PCT Application No. PCT/EP2007/009939, dated Nov. 16, 2009, 11 pages.
PCT International Preliminary Report on Patentability for PCT Application PCT/US2007/069638, dated Dec. 10, 2008, 5 pages.
PCT International Search Report and the Written Opinion of the International Searching Authority, for PCT Application No. PCT/EP2005/011011, dated Oct. 11, 2005, 12 pages.
PCT Notification of Transmittal of the International Search Report and the Written Opinion of the International Searching Authority, or the Declaration, for PCT Application No. PCT/EP2005/011011, dated Apr. 24, 2006, 12 pages.
PCT Notification of Transmittal of the International Search Report and the Written Opinion of the International Searching Authority, or the Declaration, for PCT Application No. PCT/EP2007/009939, dated May 21, 2008, 13 pages.
PCT Notification of Transmittal of the International Search Report and the Written Opinion of the International Searching Authority, or the Declaration, for PCT Application No. PCT/EP2008/004729, dated Oct. 29, 2008, 9 pages.
PCT Notification of Transmittal of the International Search Report and the Written Opinion of the International Searching Authority, or the Declaration for PCT application No. PCT/US2007/069638, dated Mar. 5, 2008, 9 pages.
PCT Written Opinion of the International Preliminary Examining Authority, PCT/EP2005/011011, dated Oct. 24, 2006, 4 pages.
PCT Written Opinion of the International Searching Authority, for PCT Application No. PCT/EP2005/011011, dated Apr. 21, 2006, 6 pages.
PCT Written Opinion of the International Searching Authority, for PCT Application No. PCT/EP2007/009939, dated May 28, 2009, 9 pages.
Peng, She et al., “Motion View Reconstruction Method with Real Object Image based on Virtual Object Movement, http://rlinks2.dialog.com/NASApp/ChannelWEB/DialogProServlet?ChName=engineering”, Eizo Joho Media Gakkai Gijutsu Hokoku, 2005, pp. 67-70, vol. 29—Issue 17.
Pickup, Lyndsey et al., “Optimizing and Learning for Super-resolution,” BMVC, Sep. 4-7 2006.
Pulli, Kari et al., “Robust Meshes from Multiple Range Maps, http://citeseer.ist.psu.edu/pulli97robust.html”, In Proc. IEEE International Conference on Recent Advances in 3-D Digital Imaging and Modeling, 1997.
Rahgozar et al., “Motion Estimation Based on Time-Sequentially Sampled Imagery”, IEEE Transactions on Image Processing, 1995.
Rav-Acha, A. et al., “Restoration of Multiple Images with Motion Blur in Different Directions, XP002375829, ISBN: 0-7695-0813-8”, Proceedings Fifth IEEE Workshop on Applications on Computer Vision IEEE Comput. Soc, 2000, pp. 22-28.
Sasaki et al., “A Wide Dynamic Range CMOS Image Sensor with Multiple Short-Time Exposures”, IEEE Proceedings on Sensors, 2004, pp. 967-972, vol. 2.
Sauer, K. et al., “Efficient Block Motion Estimation Using Integral Projections”, IEEE Trans. Circuits, Systems for video Tech, 1996, pp. 513-518, vol. 6—Issue 5.
Schultz, Richard R. et al., “Extraction of High-Resolution Frames from Video Sequences, http://citeseer.ist.psu.edu/schultz96extraction.html”, IEEE transactions on image processing, 1996, pp. 996-1011.
Seldin et al., “Iterative blind deconvolution algorithm applied to phase retrieval”, J. Opt. Soc. Am. A, 1990, pp. 428-433, vol. 7—Issue 3.
Siu, Angus et al., “Image registration for image-based rendering, http://rlinks2.dialog.com/NASApp/ChannelWEB/DialogProServlet?ChName=engineering”, IEEE Transactions on Image Processing , 2005, pp. 241-252, vol. 14—Issue 2.
Trussell, H.J. et al., “Identification and restoration of spatially variant motion blurs in sequential images, XP002375828”, IEEE Transactions on Image Processing, 1992, pp. 123-126, vol. 1—Issue 1.
Uomori, K. et al., “Automatic image stabilizing system by fulldigital signal processing” IEEE Transactions on Consumer Electronics, 1990, vol. 36, No. 3, pp. 510-519.
Wei, Jiang et al., “Dense Panoramic Imaging and 3D Reconstruction by Sensors Fusion, http://rlinks2.dialog.com/NASApp/ChannelWEB/DialogProServlet?ChName=engineering”, Japan Sci. and Technol. Agency, JPN(JST); National Inst. Industrial Safety, JPN Nippon Kikai Gakkai Robotikusu, Mekatoronikusu Koenkai Koen Ronbunshu (CD-ROM), 2006, pp. 2P1-C15.
Yang, Chen-Kui et al., “Color image sharpening by moment-preserving technique”, Signal Processing, 1995, pp. 397-403, vol. 45—Issue 3, Elsevier Science Publishers.
Yitzhaky Y. et al., “Identification of Blur Parameters from Motion Blurred Images,” Graphical Model and Image Processing, 1997, 59 (5), 310-320.
Yuan, Lu et al., “Image Deblurring with Blurred/Noisy Image Pairs,” SIGGRAPH07, Aug. 5-9, 2007.
Zhang, Junping et al., “Change detection for the urban area based on multiple sensor information fusion, http://rlinks2.dialog.com/NASApp/ChannelWEB/DialogProServlet?ChName=engineering”, IEEE International Geoscience and Remote Sensing Symposium, 2005, p. 4, IEEE.
Zitova, Barbara et al., “Image registration methods: a survey,” Image and Vision Computing, 2003, pp. 977-1000, vol. 21.
Related Publications (1)
Number Date Country
20110102638 A1 May 2011 US
Provisional Applications (1)
Number Date Country
60893116 Mar 2007 US
Continuations (1)
Number Date Country
Parent 12330719 Dec 2008 US
Child 12901577 US
Continuation in Parts (1)
Number Date Country
Parent 11856721 Sep 2007 US
Child 12330719 US