WO2017153410A1 - Method for generating a noise-reduced image based on a noise model of multiple images, as well as camera system and motor vehicle - Google Patents

Method for generating a noise-reduced image based on a noise model of multiple images, as well as camera system and motor vehicle Download PDF

Info

Publication number
WO2017153410A1
WO2017153410A1 PCT/EP2017/055328 EP2017055328W WO2017153410A1 WO 2017153410 A1 WO2017153410 A1 WO 2017153410A1 EP 2017055328 W EP2017055328 W EP 2017055328W WO 2017153410 A1 WO2017153410 A1 WO 2017153410A1
Authority
WO
WIPO (PCT)
Prior art keywords
image
noise
noisy
images
value
Prior art date
Application number
PCT/EP2017/055328
Other languages
French (fr)
Inventor
David Hurych
Pavel Krizek
Jiri Kula
Original Assignee
Connaught Electronics Ltd.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Connaught Electronics Ltd. filed Critical Connaught Electronics Ltd.
Publication of WO2017153410A1 publication Critical patent/WO2017153410A1/en

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10016Video; Image sequence
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30248Vehicle exterior or interior
    • G06T2207/30252Vehicle exterior; Vicinity of vehicle

Definitions

  • the invention relates to a method for generating a noise-reduced image.
  • a noisy image is captured by a camera of a motor vehicle and an image sequence of images is captured by the camera.
  • the invention also relates to a camera system for a motor vehicle as well as to a motor vehicle with a corresponding camera system.
  • US 2014/0126808 A1 describes an approach to denoise a digital camera image. The approach is based on direct measurement of the local statistic structure of natural images in a large training set, which are impaired by faked digital camera noise.
  • the disturbing pixels deviate from those of the actual image in color and brightness.
  • a so-called signal to noise ratio is a measure of the noise component.
  • the appearance of the image noise is not directly comparable to the so-called "grain" in the photography on conventional film material, but has similar effects on the technical image quality, in particular the detail resolution.
  • the image noise is dark noise in a large part. Thus, it occurs without light being incident on the sensor.
  • the reason for this image noise is the dark current of the individual light sensitive image sensor elements on the one hand, also noise of a read-out amplifier of the image sensor on the other hand (read-out noise).
  • this object is solved by a method, by a camera system as well as by a motor vehicle having the features according to the respective independent claims.
  • a noise-reduced image is generated.
  • a noisy image is captured by a camera of a motor vehicle and the following steps are performed: a) capturing an image sequence of images by the camera;
  • the invention is based on the realization that the weighting value can be more precisely determined by determining the noise model based on the images of the image sequence.
  • the weighting value in turn provides the basis for combining the noisy image with the blurred image. Thereby, more informative details are conserved in the noise-reduced image than it would be the case in combining the noisy image with the blurred image without the weighting value.
  • the noisy image has image noise.
  • the noisy image is for example generated such that an image sensor of the camera is impaired by undesired thermal influence during the capture of the noisy image or the image sensor is exposed to undesired electromagnetic radiation during the capture of the noisy image.
  • a global smoothing filter as for example a Gaussian filter
  • a global smoothing filter can for example be applied to the entire noisy image in unchanged manner - as in known methods.
  • This is disadvantageous in that informative details and structures of the noisy image, such as for example edges, thus for example contours of objects, are herein also impaired.
  • a locally adapted noise reduction is performed by the method according to the invention.
  • the weighting value can be determined in particular for various areas of the noisy image, whereby the various areas can then be smoothed in differently intense manner.
  • the area of the image is respectively described as an individual pixel, in particular as an individual intensity value of a color channel of the pixel.
  • denoising of the noisy image can thereby be effected with the accuracy of one pixel.
  • the noisy image is thereby preferably individually adapted or denoised for each pixel of the image, in particular each intensity value of an image channel of the noisy image.
  • the noise-reduced image is generated by a linear
  • the weighting value can for example be between 0 and 1 , wherein the combination of the weighted blurred image and the weighted noisy image then together preferably again arrives at a weighting value of 1 .
  • the linear combination is performed on the level of the intensity values of the image channels of the weighted blurred image and the weighted noisy image, respectively.
  • the intensity values are then proportionally added to each other depending on the weighting value to generate a new intensity value of the noise-reduced image.
  • the noise- reduced image is generated adapted to situation respectively for different image areas in simple manner. For example, this can then be formulated as follows:
  • x column coordinate of an intensity value of the respective image
  • y line coordinate of the intensity value of the respective image.
  • step b) at least each one image area of the images of the image sequence is selected based on at least one predetermined selection criterion and the noise model is accordingly only determined based on the respective image areas in step b).
  • the steps c) to g) are preferably also only performed for the respective image areas.
  • the respective image area preferably includes a plurality of pixels, in particular 1 00 pixels.
  • the image area can also include only one pixel.
  • the noise model is then preferably determined for each of the pixels.
  • the noise model is determined for each of the pixels with a different intensity value. Thereby, the noise model can be provided for each of the intensity values.
  • the noisy image can be pixel-precisely weighted and thus be pixel-precisely combined with the blurred image.
  • image areas of the images of the image sequence are determined by the selection criterion, which substantially show the same location of the environmental region unchanged over the capture period of time of the image sequence.
  • those image areas are therefore selected, which have a static background and a static foreground. This is advantageous since thereby variations of the intensity values of the image, which relate to a movement of the camera and/or a movement of an object in the images, can be effectively excluded.
  • the image areas only those pixels remain, which should not have any intensity variation apart from the noise.
  • the noise model can then be determined in step b) only from intensity values, which are associated with static pixels and thus are substantially varied by the noise to be modeled by the noise model.
  • the noise model can be more precisely and more faultlessly determined by the predetermined selection criterion.
  • a variation rate of an average intensity value, in particular an arithmetic average, of locally corresponding intensity values of the images of the image sequence over the temporal progress of the image sequence is determined as the selection criterion of the image areas, and the image areas are selected if the variation rate is less than a predetermined variation limit value.
  • the average intensity value of the locally corresponding intensity values has a variation rate, which is less than the predetermined variation limit value.
  • the locally corresponding intensity values are pixels of consecutively captured images, which depict the same location of a scene, for example of an environmental region of the motor vehicle.
  • a distribution curve of differences of respective locally corresponding intensity values of the images of the image sequence with an average intensity value, in particular an arithmetic average, of the locally corresponding intensity values of the images is determined as the selection criterion of the image areas, and the image areas are selected if the distribution curve of the differences has an average value of 0 and/or the distribution curve of the differences is symmetrical to the ordinate axis extending to the maximum of the distribution curve.
  • the differences of the respective locally corresponding intensity values are considered as distribution, wherein the distribution is represented by the distribution curve.
  • the distribution curve in turn can then be assessed based on the average value and/or the symmetry to select the image areas.
  • the steps b) to g) are respectively performed for at least one intensity value of a color channel of the noisy image, in particular an RGB color channel of the noisy image or a YUV color channel of the noisy image.
  • the steps b) to g) are performed for the intensity value of the color channel of the noisy image
  • the noise- reduced image can be locally generated in locally precise manner.
  • an intensity value-precise generation of the noise-reduced image can thereby be effected, whereby the noise in the noisy image can be considerably reduced, but at the same time contours of objects can further be obtained in sharp manner or intense in contrast in the noisy image.
  • the noisy image can also be present in other color channels than the RGB color channel or the YUV color channel.
  • the noise model is determined in step b) based on intensity values of locally corresponding pixels of the images of the image sequence.
  • the noise model can be locally precisely applied for weighting the noisy image and/or for weighting the blurred image.
  • the noise-reduced image is generated precisely and sharply for an observer and free of noise at the same time.
  • the locally corresponding pixels of the images of the image sequence are combined with respect to respective average intensity values in determining the noise model in step b) and a noise characteristic, in particular a standard deviation, for at least one value from a range of values of the intensity values of the images of the image sequence is respectively provided by the noise model.
  • the noise model is first determined preferably for each intensity value of the images of the image sequence.
  • the noise model can provide the noise characteristic preferably for each value from the range of values of the intensity values of the images, thus for example for an 8 bit image from 0 to 255.
  • the noise characteristic is in particular described by the standard deviation of the respective value.
  • an in particular arithmetic average is respectively formed from the locally corresponding pixels of the images of the image sequence and the noise model is determined depending on all of the corresponding pixels with the same average value.
  • the noise characteristic for a single intensity value is no longer determined only based on a pixel of the images of the image sequence, but preferably based on all of the pixels of the images of the image sequence with the same average value.
  • a further image is captured by the camera depending on an update criterion and the noise model is again determined additionally considering the further image.
  • the update criterion in particular a so-called online update of the noise model is performed.
  • the noise model can for example then be adapted from time to time. For example, a predetermined period of time or else a predetermined number of captured images can be preset as the update criterion.
  • the update criterion can for example also be determined by how a brightness state in an environmental region of the camera varies.
  • the update criterion can for example depend on a value of automatic exposure setting of the camera.
  • the update criterion can be preset as a variation of a signal gain value of the camera or an exposure time variation of the camera.
  • the update criterion is advantageous in that the noise model can thereby always be provided currently adapted to the noise causing circumstances and thus the noise-reduced image can always be generated in a detail conserving and informative manner.
  • the weighting value is determined depending on a magnitude of an edge of the noisy image, and the blurred image is the less severely taken into account in weighting in step e) the higher the magnitude of the edge is and/or the noisy image is the more severely weighted in weighting in step f) the higher the magnitude of the edge is.
  • the blurred image and the noisy image depending on the magnitude of the edge of the noisy image, it is allowed that areas of the noisy image are maintained, which for example contain a contour of an object, but homogeneous areas of the image, which then in particular do not contain contours of objects, are again considerably smoothed.
  • the noise reduction of the noisy image can thereby be particularly accurately effected in locally adapted manner.
  • the magnitude of the edge it is indicated how severely the edge differs from the adjacent intensity values.
  • step a) or step b) is performed at a point of time, at which an automatic exposure setting of the camera is performed and/or at which the motor vehicle is statically operated.
  • the image sequence of the images is preferably captured when at least the camera of the motor vehicle is unmoved and thus the respective selection criterion can be more easily satisfied.
  • the noise model is then respectively newly determined when the automatic exposure setting of the camera is newly performed. For example, this can be the case in a variation of a brightness state in the environmental region of the camera. Because after the variation of the brightness state in the environmental region of the camera, another noise characteristic can also be assumed, the determination of the noise model at this point of time is particularly effective. At this point of time, the new noise model then considerably differs from the old noise model.
  • step a) and/or step b) can also be performed directly after activating the camera, for example upon starting the motor vehicle.
  • the blurred image is determined in step d) by convolution of the noisy image with a Gaussian kernel, in particular with an average value of 0.
  • the blurred image in step d) is determined by the application of the Gaussian kernel or a Gaussian filter to the noisy image.
  • the Gaussian kernel is in particular present as a two- dimensional filter and can be applied in various pixel sizes, for example 3x3, 5x5 or 7x7.
  • a standard deviation of the Gaussian kernel can also be preset with various values.
  • the invention also relates to a camera system for a motor vehicle with a camera and an evaluation unit, which is formed to perform a method according to the invention.
  • the camera system can include multiple cameras.
  • the camera has a motor vehicle fixing element for fixing to the motor vehicle.
  • the camera system can be a component of a driver assistance system of the motor vehicle.
  • the invention relates to a motor vehicle with a camera system according to the invention.
  • the preferred embodiments presented with respect to the method according to the invention and the advantages thereof correspondingly apply to the camera system according to the invention as well as to the motor vehicle according to the invention.
  • Fig. 1 a schematic plan view of an embodiment of a motor vehicle according to the invention with a camera system
  • Fig. 2 a schematic representation of a noisy image captured by a camera of the camera system
  • Fig. 3 a diagrammatic representation of weighting values of a noise model for generating a noise-reduced image based on the noisy image and a blurred image generated based on the noisy image;
  • Fig. 4 a schematic representation of an image of an image sequence provided for determining the noise model, which is captured by the camera; a schematic representation of image areas of the image of the image sequence, which are selected based on two predetermined selection criteria; a schematic representation of a first image, a second image and a third image of the image sequence with locally corresponding intensity values for determining the second selection criterion; a schematic representation of a distribution curve of differences of respective locally corresponding intensity values of the images of the image sequence to determine the second selection criterion; a schematic representation of the images of the image sequence and a locally corresponding intensity value of the images of the image sequence, which is registered in a vector; a diagrammatic representation of noise characteristics of the noise model, which is determined based on the images of the image sequence; a diagrammatic illustration of exemplary standard deviations of the noise characteristic of the noise model; a schematic representation of identical intensity values of the images of the image sequence for an intensity value of 100, which are combined and collectively contribute to the noise characteristic of an intensity value of the noise model
  • a plan view of a motor vehicle 1 with a camera system 2 is schematically illustrated.
  • the camera system 2 includes a camera 3 and an evaluation unit 4.
  • the camera 3 is disposed at a front 5 of the motor vehicle 1 .
  • the arrangement of the camera 3 is variously possible, however, preferably such that an environmental region 6 of the motor vehicle 1 can at least partially be captured by the camera 3.
  • the arrangement of the evaluation unit 4 is also variously possible, however, preferably such that it can be connected to the camera 3.
  • the evaluation unit 4 can be integrated in the camera 3 or be formed separately from the camera 3.
  • the camera system 2 can for example also include multiple cameras 3.
  • the evaluation unit can 4 for example be connected to multiple cameras 3 or else also be formed by multiple partial evaluation units.
  • the camera 3 is formed as a CMOS (complementary metal-oxide- semiconductor) camera or as a CCD (charge-coupled device) camera or else as a versatile image capturing device.
  • the camera 3 is in particular formed as a video camera, which continuously provides an image sequence of frames.
  • the camera 3 has an objective 7 and an image sensor 8.
  • the image sensor 8 can add noise to a picture in generating the picture of the environmental region 6 by different impairment of heat.
  • the noise added to the picture can for example also be added to the picture by electromagnetic influences on the image sensor 8.
  • Fig. 2 shows a noisy image l N .
  • the noisy image l N is captured from the environmental region 6 of the motor vehicle 1 .
  • noise in particular caused by the image sensor 8 is contained.
  • a noise-reduced image l D is now generated from the noisy image l N .
  • a blurred image l B is generated based on the noisy image l N .
  • the noisy image l N and the blurred image l B are weighted by a weighting value w xy and linearly combined.
  • the weighting value w xy is provided by a noise model NLF.
  • i D ( > y) (i - w xy )i N ( > y) + w X y ( > y) 0 )
  • a position of an intensity value of a color channel or of a pixel of the respective image l D , l N , IB is described by the line coordinate x and the column coordinate y.
  • the blurred image l B is generated by a convolution of the noisy image l N by means of a discrete two-dimensional Gaussian kernel with an average value of 0, a filter dimension [sxs] of the Gaussian kernel and a standard deviation of ⁇ ⁇ .
  • the filter dimension [sxs] and ⁇ ⁇ of the blurred image l B as well as a gradient magnitude normalization constant g are preferably determined based on a number C of images with a respective image dimension [XxY]. This is effected by optimization of the following mathematical expression.
  • f is a denoising function
  • I is an ideal image without noise
  • I N is the j-th noisy image, which is captured by the camera 3.
  • the weighting value w xy can be mathematically described as follows:
  • the maximization of the term is provided in case that a normalized gradient is greater than NLF(l(x,y)) and thus a result would have a negative value.
  • the weighting value in particular has a range of values of w xy e (0,1) .
  • the partial derivatives are estimated by the following equations. di N (x, y)
  • the normalized gradient magnitude describes a local intensity of a texture of the noisy image l N .
  • the more intensely the texture is formed the less intensely the blurred image l B is weighted in equation (1 ) and the more intensely the noisy image l N is weighted.
  • each intensity value of the noisy image l N and of the blurred image l B obtains an own weighting value w xy , which indicates how severely it is denoised or smoothed.
  • Fig. 3 shows, which weighting values w xy are applied for the noise-reduced image
  • the noisy image l N is present in an RGB color space or in a YUV color space.
  • the method is preferably performed for each intensity value l N (x,y) of the noisy image l N .
  • the method is applied to provide the noise-reduced image l D in real time, thus in the time, in which a new noisy image is not yet captured by the camera.
  • the position of an upper left corner coordinate is described by x and y.
  • a width of the image area p a is described by W and a height of the image area p a is described by H.
  • the index a describes an ordinal number for the linkage to a certain image 9 from the image sequence 10.
  • the x and y coordinates are linked to a certain image area p a while the width W and the height H are the same for all of the image areas p a .
  • the weight W can for example be preset with a parameter value of ten pixels and the height H can for example also be preset with a parameter value of ten pixels.
  • the image areas p a are in particular separately examined for each image channel of the respective color space. In the further, the description is only presented for a respective intensity value l(x,y) of one of the images 9 of the image sequence 1 0.
  • the first selection criterion 1 1 implies that an average intensity value ⁇ ⁇ ' over a time t has a variation rate below a predetermined variation limit value ⁇ . This can be mathematically described as follows.
  • the average intensity value ⁇ ⁇ ' can be described as follows. W-l H-l
  • the first selection criterion 1 1 can then be defined as follows.
  • predetermined variation limit value ⁇ which is for example one to two percent of the value range of intensity values of the images 9.
  • the number of images 9 are specified.
  • an overall average intensity value ⁇ ⁇ ⁇ is subtracted from the average intensity value ⁇ ⁇ ' .
  • the image areas p a which satisfy the condition from equation (8) of the first selection criterion 1 1 , are forwarded for examining the second selection criterion 12.
  • the noise in the noisy image l N substantially follows a Gaussian distribution.
  • the image area p a is not suitable with respect to the second selection criterion 12.
  • the second selection criterion 12 implies that differences A j , of locally corresponding intensity values Vji should have a distribution curve over the time t, which shows an average value of 0 or else is formed symmetrically to an ordinate axis extending to the maximum of the distribution curve.
  • a vector is preset, in which the intensity values of a pixel of an image area are preset over the changing time and within an interval j and an interval I.
  • the interval j corresponds to j e ⁇ 0,...,W - 1 ⁇ and the interval I therein corresponds to l e ⁇ 0,...,H - l ⁇ .
  • the differences A j of the locally corresponding intensity values Vj, are obtained by subtracting an average vector b.
  • a jl (t) v jl (t) ⁇ j v jl (b) te ⁇ l,...,A ⁇ (10)
  • a distribution curve ⁇ which represents the distribution of the differences in particular by a discrete histogram.
  • the maximum of the sum of the distributions is searched and it is examined if it is close to 0. If this is the case, then, the distribution curve ⁇ is examined with respect to its symmetry to the ordinate axis 13 extending to the maximum of the distribution curve ⁇ .
  • the symmetry can simply be examined by subtracting an integral from the part of the distribution curve ⁇ disposed to the left of the ordinate axis 13 from an integral of a part of the distribution curve ⁇ disposed to the right of the ordinate axis 13. If the difference is less than a predetermined difference limit value, it can be assumed that the distribution curve ⁇ is symmetrical and thus corresponds a sum of Gaussian distribution.
  • Fig. 4 shows one of the images 9 of the image sequence 10.
  • the images 9 in particular occur with the number A in the image sequence 10.
  • the image areas p a are determined.
  • Fig. 5 shows the examined image areas p a of the images 9, which have been determined based on the image shown in Fig. 4.
  • the hatched image areas 14 satisfy both the first selection criterion 1 1 and the second selection criterion 12, while the non-hatched image areas 15 do not satisfy at least one of the selection criteria 1 1 , 12.
  • Fig. 6 shows a first image 16 of the images 9, a second image 17 of the images 9 and a third image 18 of the images 9.
  • the images 16, 17, 18 are captured consecutively in time
  • the locally corresponding image area p a is drawn.
  • first locally corresponding intensity values v 1 ; second locally corresponding intensity values v 2 and third locally corresponding intensity values v 3 are drawn.
  • the locally corresponding intensity values v 1 ; v 2 , v 3 are in particular preset as a vector.
  • a fourth distribution curve 19 is described by equation (15).
  • Fig. 8 now shows the images 9 of the image sequence 10, based on which the noise model NLF is generated.
  • the respective noise model NLF is determined for each color channel of the images 9 of the image sequence 10.
  • the color channels can for example be present in an RGB color space or else in a YUV color space.
  • the creation of the noise model NLF is explained in the following based on a color channel of the RGB color space, wherein the method is also applicable to the YUV color space.
  • the noise model NLF is mathematically described as follows.
  • a noise characteristic ⁇ or a standard deviation of the noise of a determined intensity value l(x,y) of an image of the images 9 of the image sequence 10 describes as an output.
  • the noise model NLF is determined for the various color channels, thus for example a red color channel, a blue color channel and a green color channel. It is shown, how the intensity values of the pixels from the images 9 over the time can be combined to a vector v.
  • Fig. 9 shows a diagram, in which an intensity value 21 of a red color channel of the images 9 of the image sequence 10 is plotted on an abscissa 20.
  • the noise characteristic ⁇ or a standard deviation of the noise is registered on an ordinate 22.
  • the noise model NLF is described for values i with the range of values of i e ⁇ 0,1,...,255 ⁇ . Thereby, thus, for each intensity value l(x,y) contained in the images 9, a noise characteristic ⁇ in the form of the standard deviation is provided by the noise model NLF.
  • a first noise characteristic ⁇ corresponds to a high standard deviation
  • the second noise characteristic ⁇ 2 corresponds to a medium standard deviation
  • the third noise characteristic ⁇ 3 corresponds to a low standard deviation.
  • the noise model NLF is calculated over all of the pixels, which are within the previously selected image areas p a and have the same average intensity value over the time t. A standard deviation over all of these pixels provides the noise model NLF in the position i. This is schematically explained based on Fig. 8.
  • the vector v also has an average value m v . It is mathematically described as follows, wherein n is a number of the images 9 of the image sequence 10.
  • the noise model NLF can for example be approximately estimated for these values or else simply be left empty.
  • a vector of the noise differences ⁇ can then be calculated as follows.
  • the noise model NLF is then determined as follows.
  • the noise model NLF is a normalized standard deviation of the noise differences ⁇ . Furthermore, a normalization constant m is indicated, which can be described as follows.
  • the noise model NLF is again determined or updated depending on an update criterion.
  • the noise model NLF is then determined as follows.
  • NLF(i) _ ⁇ ⁇ i (22)
  • T Q three scalars T 0 , T ⁇ and T 2 are introduced here, which together are described as T Q . If new noise differences r, are then to be added, it is only required to form the sums of all of the scalars described by T Q and to newly determine the noise model NLF.
  • the scalar Tc can be mathematically described as follows.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)
  • Studio Devices (AREA)

Abstract

The invention relates to a method for generating a noise-reduced image (lD), in which a noisy image (lN) is captured by a camera (3) of a motor vehicle (1) and the following steps are performed: a) capturing an image sequence (10) of images (9) by the camera (3); b) determining a noise model (NLF) based on multiple images (9) of the image sequence (10); c) determining a weighting value (wxy) for weighting the noisy image (lN) as well as at least one image (lB) depending on the noisy image (lN) based on the noise model (NLF); d) generating a blurred image (lB) based on the noisy image (lN); e) generating a weighted blurred image by weighting the blurred image (lB) depending on the weighting value (wxy); f) generating a weighted noisy image by weighting the noisy image (lN) depending on the weighting value (wxy); and g) generating the noise-reduced image (lD) based on the weighted blurred image and the weighted noisy image.

Description

Method for generating a noise-reduced image based on a noise model of multiple images, as well as camera system and motor vehicle
The invention relates to a method for generating a noise-reduced image. A noisy image is captured by a camera of a motor vehicle and an image sequence of images is captured by the camera. The invention also relates to a camera system for a motor vehicle as well as to a motor vehicle with a corresponding camera system.
Methods for generating a noise-reduced image are known from the prior art. Thus, US 2014/0126808 A1 describes an approach to denoise a digital camera image. The approach is based on direct measurement of the local statistic structure of natural images in a large training set, which are impaired by faked digital camera noise.
One designates the deterioration of a digital or electronically captured image by disturbances, which do not have any relation to the actual image content, the image signal, as image noise. The disturbing pixels deviate from those of the actual image in color and brightness. Therein, a so-called signal to noise ratio is a measure of the noise component. The appearance of the image noise is not directly comparable to the so-called "grain" in the photography on conventional film material, but has similar effects on the technical image quality, in particular the detail resolution.
In electronic image sensors, such as for example CCD (charge-coupled device) sensors or CMOS (complementary metal-oxide-semiconductor) sensors, the image noise is dark noise in a large part. Thus, it occurs without light being incident on the sensor. The reason for this image noise is the dark current of the individual light sensitive image sensor elements on the one hand, also noise of a read-out amplifier of the image sensor on the other hand (read-out noise).
It is the object of the invention to provide a method, a camera system as well as a motor vehicle, by which a more noise-reduced image can be generated in more detail conserving manner.
According to the invention, this object is solved by a method, by a camera system as well as by a motor vehicle having the features according to the respective independent claims. In a method according to the invention, a noise-reduced image is generated. A noisy image is captured by a camera of a motor vehicle and the following steps are performed: a) capturing an image sequence of images by the camera;
and as an essential idea of the invention
b) determining a noise model based on multiple images of the image sequence;
c) determining a weighting value for weighting the noisy image as well as at least one image depending on the noisy image based on the noise model;
d) generating a blurred image based on the noisy image;
e) generating a weighted blurred image by weighting the blurred image depending on the weighting value;
f) generating a weighted noisy image by weighting the noisy image depending on the weighting value; and
g) generating the noise-reduced image based on the weighted blurred image and the weighted noisy image.
The invention is based on the realization that the weighting value can be more precisely determined by determining the noise model based on the images of the image sequence. The weighting value in turn provides the basis for combining the noisy image with the blurred image. Thereby, more informative details are conserved in the noise-reduced image than it would be the case in combining the noisy image with the blurred image without the weighting value.
The noisy image has image noise. Thus, the noisy image is for example generated such that an image sensor of the camera is impaired by undesired thermal influence during the capture of the noisy image or the image sensor is exposed to undesired electromagnetic radiation during the capture of the noisy image.
In order to reduce the image noise in the noisy image, now, a global smoothing filter, as for example a Gaussian filter, can for example be applied to the entire noisy image in unchanged manner - as in known methods. This is disadvantageous in that informative details and structures of the noisy image, such as for example edges, thus for example contours of objects, are herein also impaired. In order to prevent this, a locally adapted noise reduction is performed by the method according to the invention. Thus, the weighting value can be determined in particular for various areas of the noisy image, whereby the various areas can then be smoothed in differently intense manner.
Preferably, the area of the image is respectively described as an individual pixel, in particular as an individual intensity value of a color channel of the pixel. In particular, thus, denoising of the noisy image can thereby be effected with the accuracy of one pixel. Thus, the noisy image is thereby preferably individually adapted or denoised for each pixel of the image, in particular each intensity value of an image channel of the noisy image.
Preferably, it is provided that the noise-reduced image is generated by a linear
combination, in particular by summation, of the weighted blurred image and the weighted noisy image. Thus, the weighting value can for example be between 0 and 1 , wherein the combination of the weighted blurred image and the weighted noisy image then together preferably again arrives at a weighting value of 1 . In particular, the linear combination is performed on the level of the intensity values of the image channels of the weighted blurred image and the weighted noisy image, respectively. Thus, the intensity values are then proportionally added to each other depending on the weighting value to generate a new intensity value of the noise-reduced image. By the linear combination, the noise- reduced image is generated adapted to situation respectively for different image areas in simple manner. For example, this can then be formulated as follows:
ID (x, y) = (l - wxy )IN (x, y) + wxyIB (x, y) , with lD = noise-reduced image;
lN = noisy image;
lB = blurred image;
wxy = weighting value;
x = column coordinate of an intensity value of the respective image;
y = line coordinate of the intensity value of the respective image.
Furthermore, it is preferably provided that before step b) at least each one image area of the images of the image sequence is selected based on at least one predetermined selection criterion and the noise model is accordingly only determined based on the respective image areas in step b). In particular, the steps c) to g) are preferably also only performed for the respective image areas. The respective image area preferably includes a plurality of pixels, in particular 1 00 pixels. However, the image area can also include only one pixel. The noise model is then preferably determined for each of the pixels. In particular, the noise model is determined for each of the pixels with a different intensity value. Thereby, the noise model can be provided for each of the intensity values. Thereby, the noisy image can be pixel-precisely weighted and thus be pixel-precisely combined with the blurred image. In particular, image areas of the images of the image sequence are determined by the selection criterion, which substantially show the same location of the environmental region unchanged over the capture period of time of the image sequence. Thus, by the selection criterion, in particular those image areas are therefore selected, which have a static background and a static foreground. This is advantageous since thereby variations of the intensity values of the image, which relate to a movement of the camera and/or a movement of an object in the images, can be effectively excluded. As the image areas, only those pixels remain, which should not have any intensity variation apart from the noise. Thereby, the noise model can then be determined in step b) only from intensity values, which are associated with static pixels and thus are substantially varied by the noise to be modeled by the noise model. Thus, the noise model can be more precisely and more faultlessly determined by the predetermined selection criterion.
Preferably, it is provided that a variation rate of an average intensity value, in particular an arithmetic average, of locally corresponding intensity values of the images of the image sequence over the temporal progress of the image sequence is determined as the selection criterion of the image areas, and the image areas are selected if the variation rate is less than a predetermined variation limit value. Thus, if there should neither move the camera nor an object in the image, thus, it can be assumed that the average intensity value of the locally corresponding intensity values has a variation rate, which is less than the predetermined variation limit value. Thus, the locally corresponding intensity values are pixels of consecutively captured images, which depict the same location of a scene, for example of an environmental region of the motor vehicle. By the variation rate of the average intensity value, the image areas are simply and reliably selected to the effect that the respective intensity value can be incorporated in the noise model without being affected by a movement of the camera and/or a movement of an object in the image.
Furthermore, it is preferably provided that a distribution curve of differences of respective locally corresponding intensity values of the images of the image sequence with an average intensity value, in particular an arithmetic average, of the locally corresponding intensity values of the images is determined as the selection criterion of the image areas, and the image areas are selected if the distribution curve of the differences has an average value of 0 and/or the distribution curve of the differences is symmetrical to the ordinate axis extending to the maximum of the distribution curve. Thus, it can be recognized based on the distribution curve if the respective image area is captured from a static region or if objects move in the image area or the camera moves in capture of the image area. Hereto, the differences of the respective locally corresponding intensity values are considered as distribution, wherein the distribution is represented by the distribution curve. The distribution curve in turn can then be assessed based on the average value and/or the symmetry to select the image areas.
Preferably, it is provided that the steps b) to g) are respectively performed for at least one intensity value of a color channel of the noisy image, in particular an RGB color channel of the noisy image or a YUV color channel of the noisy image. In that the steps b) to g) are performed for the intensity value of the color channel of the noisy image, the noise- reduced image can be locally generated in locally precise manner. In contrast to global or coarse image improvement methods, in the method according to the invention, an intensity value-precise generation of the noise-reduced image can thereby be effected, whereby the noise in the noisy image can be considerably reduced, but at the same time contours of objects can further be obtained in sharp manner or intense in contrast in the noisy image. However, the noisy image can also be present in other color channels than the RGB color channel or the YUV color channel.
Furthermore, it is preferably provided that the noise model is determined in step b) based on intensity values of locally corresponding pixels of the images of the image sequence. Thereby, the noise model can be locally precisely applied for weighting the noisy image and/or for weighting the blurred image. Thereby, the noise-reduced image is generated precisely and sharply for an observer and free of noise at the same time.
Furthermore, it is preferably provided that the locally corresponding pixels of the images of the image sequence are combined with respect to respective average intensity values in determining the noise model in step b) and a noise characteristic, in particular a standard deviation, for at least one value from a range of values of the intensity values of the images of the image sequence is respectively provided by the noise model. Thus, the noise model is first determined preferably for each intensity value of the images of the image sequence. Then, the noise model can provide the noise characteristic preferably for each value from the range of values of the intensity values of the images, thus for example for an 8 bit image from 0 to 255. The noise characteristic is in particular described by the standard deviation of the respective value. Thereby, a noise
characteristic, which is formed over multiple pixels of the images of the image sequence, can then for example be assigned to each intensity value of the noisy image.
Furthermore, it can preferably be provided that an in particular arithmetic average is respectively formed from the locally corresponding pixels of the images of the image sequence and the noise model is determined depending on all of the corresponding pixels with the same average value. Thus, the noise characteristic for a single intensity value is no longer determined only based on a pixel of the images of the image sequence, but preferably based on all of the pixels of the images of the image sequence with the same average value. By combining the noise characteristic with respect to the average value, a particularly reliable and significant noise characteristic for the value from the range of values of the intensity values is provided by means of the noise model. Thus, influences of erroneous individual standard deviations of a single intensity value can thereby for example be recognized and compensated for by the amount of the intensity values with the same average value. Thus, the effect of outliers in the determination of the noise model can be reduced. Thereby, the noise model can be still more precisely and flawlessly applied for weighting the blurred image and/or the noisy image.
In a further embodiment, it is preferably provided that after step a) a further image is captured by the camera depending on an update criterion and the noise model is again determined additionally considering the further image. Thus, by the update criterion, in particular a so-called online update of the noise model is performed. Thereby, the noise model can for example then be adapted from time to time. For example, a predetermined period of time or else a predetermined number of captured images can be preset as the update criterion. However, the update criterion can for example also be determined by how a brightness state in an environmental region of the camera varies. Thus, the update criterion can for example depend on a value of automatic exposure setting of the camera. For example, the update criterion can be preset as a variation of a signal gain value of the camera or an exposure time variation of the camera. The update criterion is advantageous in that the noise model can thereby always be provided currently adapted to the noise causing circumstances and thus the noise-reduced image can always be generated in a detail conserving and informative manner.
Preferably, it is provided that the weighting value is determined depending on a magnitude of an edge of the noisy image, and the blurred image is the less severely taken into account in weighting in step e) the higher the magnitude of the edge is and/or the noisy image is the more severely weighted in weighting in step f) the higher the magnitude of the edge is. By weighting the blurred image and the noisy image depending on the magnitude of the edge of the noisy image, it is allowed that areas of the noisy image are maintained, which for example contain a contour of an object, but homogeneous areas of the image, which then in particular do not contain contours of objects, are again considerably smoothed. Thus, the noise reduction of the noisy image can thereby be particularly accurately effected in locally adapted manner. By the magnitude of the edge, it is indicated how severely the edge differs from the adjacent intensity values.
Furthermore, it can be provided that step a) or step b) is performed at a point of time, at which an automatic exposure setting of the camera is performed and/or at which the motor vehicle is statically operated. Thus, the image sequence of the images is preferably captured when at least the camera of the motor vehicle is unmoved and thus the respective selection criterion can be more easily satisfied. In particular, the noise model is then respectively newly determined when the automatic exposure setting of the camera is newly performed. For example, this can be the case in a variation of a brightness state in the environmental region of the camera. Because after the variation of the brightness state in the environmental region of the camera, another noise characteristic can also be assumed, the determination of the noise model at this point of time is particularly effective. At this point of time, the new noise model then considerably differs from the old noise model. However, step a) and/or step b) can also be performed directly after activating the camera, for example upon starting the motor vehicle.
Preferably, it is provided that the blurred image is determined in step d) by convolution of the noisy image with a Gaussian kernel, in particular with an average value of 0. Thus, the blurred image in step d) is determined by the application of the Gaussian kernel or a Gaussian filter to the noisy image. The Gaussian kernel is in particular present as a two- dimensional filter and can be applied in various pixel sizes, for example 3x3, 5x5 or 7x7. A standard deviation of the Gaussian kernel can also be preset with various values.
The invention also relates to a camera system for a motor vehicle with a camera and an evaluation unit, which is formed to perform a method according to the invention. The camera system can include multiple cameras. Therein, the camera has a motor vehicle fixing element for fixing to the motor vehicle.
For example, the camera system can be a component of a driver assistance system of the motor vehicle.
Furthermore, the invention relates to a motor vehicle with a camera system according to the invention. The preferred embodiments presented with respect to the method according to the invention and the advantages thereof correspondingly apply to the camera system according to the invention as well as to the motor vehicle according to the invention.
Further features of the invention are apparent from the claims, the figures and the description of figures. The features and feature combinations mentioned above in the description as well as the features and feature combinations mentioned below in the description of figures and/or shown in the figures alone are usable not only in the respectively specified combination, but also in other combinations without departing from the scope of the invention. Thus, implementations are also to be considered as encompassed and disclosed by the invention, which are not explicitly shown in the figures and explained, but arise from and can be generated by separated feature combinations from the explained implementations. Implementations and feature combinations are also to be considered as disclosed, which thus do not have all of the features of an originally formulated independent claim. Moreover, implementations and feature combinations are also to be considered as disclosed, in particular by the explanations set out above, which extend beyond or deviate from the feature combinations set out in the relations of the claims.
Below, the embodiments of the invention are explained in more detail based on schematic drawings.
There show:
Fig. 1 a schematic plan view of an embodiment of a motor vehicle according to the invention with a camera system;
Fig. 2 a schematic representation of a noisy image captured by a camera of the camera system;
Fig. 3 a diagrammatic representation of weighting values of a noise model for generating a noise-reduced image based on the noisy image and a blurred image generated based on the noisy image;
Fig. 4 a schematic representation of an image of an image sequence provided for determining the noise model, which is captured by the camera; a schematic representation of image areas of the image of the image sequence, which are selected based on two predetermined selection criteria; a schematic representation of a first image, a second image and a third image of the image sequence with locally corresponding intensity values for determining the second selection criterion; a schematic representation of a distribution curve of differences of respective locally corresponding intensity values of the images of the image sequence to determine the second selection criterion; a schematic representation of the images of the image sequence and a locally corresponding intensity value of the images of the image sequence, which is registered in a vector; a diagrammatic representation of noise characteristics of the noise model, which is determined based on the images of the image sequence; a diagrammatic illustration of exemplary standard deviations of the noise characteristic of the noise model; a schematic representation of identical intensity values of the images of the image sequence for an intensity value of 100, which are combined and collectively contribute to the noise characteristic of an intensity value of the noise model; and a schematic representation analogously to Fig. 1 1 , but for an intensity value of 200.
In the figures, identical or functionally identical elements are provided with the same reference characters. In Fig. 1 , a plan view of a motor vehicle 1 with a camera system 2 is schematically illustrated. The camera system 2 includes a camera 3 and an evaluation unit 4. According to the embodiment, the camera 3 is disposed at a front 5 of the motor vehicle 1 . However, the arrangement of the camera 3 is variously possible, however, preferably such that an environmental region 6 of the motor vehicle 1 can at least partially be captured by the camera 3. The arrangement of the evaluation unit 4 is also variously possible, however, preferably such that it can be connected to the camera 3. For example, the evaluation unit 4 can be integrated in the camera 3 or be formed separately from the camera 3. The camera system 2 can for example also include multiple cameras 3. Then, the evaluation unit can 4 for example be connected to multiple cameras 3 or else also be formed by multiple partial evaluation units.
In particular, the camera 3 is formed as a CMOS (complementary metal-oxide- semiconductor) camera or as a CCD (charge-coupled device) camera or else as a versatile image capturing device. The camera 3 is in particular formed as a video camera, which continuously provides an image sequence of frames.
The camera 3 has an objective 7 and an image sensor 8. Therein, the image sensor 8 can add noise to a picture in generating the picture of the environmental region 6 by different impairment of heat. However, the noise added to the picture can for example also be added to the picture by electromagnetic influences on the image sensor 8.
Fig. 2 shows a noisy image lN. The noisy image lN is captured from the environmental region 6 of the motor vehicle 1 . In the noisy image lN, noise in particular caused by the image sensor 8 is contained. A noise-reduced image lD is now generated from the noisy image lN. Hereto, a blurred image lB is generated based on the noisy image lN. In order to generate the noise-reduced image lD, then, the noisy image lN and the blurred image lB are weighted by a weighting value wxy and linearly combined. The weighting value wxy is provided by a noise model NLF. The linear combination and thus the generation of the noise-reduced image lD can be mathematically represented as follows. iD ( > y) = (i - wxy )iN ( > y) + w Xy ( > y) 0 )
A position of an intensity value of a color channel or of a pixel of the respective image lD, lN, IB is described by the line coordinate x and the column coordinate y. The blurred image lB is generated by a convolution of the noisy image lN by means of a discrete two-dimensional Gaussian kernel with an average value of 0, a filter dimension [sxs] of the Gaussian kernel and a standard deviation of σΒ. By applying the Gaussian kernel to the noisy image lN, intense gradients are reduced. Thus, the noisy image lN is smoothed or rendered blurred by the Gaussian kernel.
The filter dimension [sxs] and σΒ of the blurred image lB as well as a gradient magnitude normalization constant g are preferably determined based on a number C of images with a respective image dimension [XxY]. This is effected by optimization of the following mathematical expression.
C X Y 2
arg min∑∑∑ (/, (x, y) - f (VN (x, y), s, σΒ , g\\ (2)
S,<7B ,g j=l X=l y=l
Herein, f is a denoising function, I, is an ideal image without noise and IN ] is the j-th noisy image, which is captured by the camera 3. These parameters are in particular only once determined and can then for example be used for all of the noise models NLF of any camera.
The weighting value wxy can be mathematically described as follows:
Figure imgf000012_0001
The maximization of the term is provided in case that a normalized gradient is greater than NLF(l(x,y)) and thus a result would have a negative value. The weighting value in particular has a range of values of wxy e (0,1) . The partial derivatives are estimated by the following equations. diN (x, y)
+ l, y) - IN (x - l, y) (4) dx dlN (x, y)
= IN (x, y + l) - IN (x, y -l) (5) In the equation (4) and the equation (5), the normalized gradient magnitude describes a local intensity of a texture of the noisy image lN. The more intensely the texture is formed, the less intensely the blurred image lB is weighted in equation (1 ) and the more intensely the noisy image lN is weighted. In particular, each intensity value of the noisy image lN and of the blurred image lB obtains an own weighting value wxy, which indicates how severely it is denoised or smoothed.
Hereto, Fig. 3 shows, which weighting values wxy are applied for the noise-reduced image
Preferably, the noisy image lN is present in an RGB color space or in a YUV color space. Thereby, the method is preferably performed for each intensity value lN(x,y) of the noisy image lN. In particular, the method is applied to provide the noise-reduced image lD in real time, thus in the time, in which a new noisy image is not yet captured by the camera.
The noise model NLF, based on which the weighting value wxy is determined, is now determined as follows. First, locally corresponding image areas pa of images 9 of an image sequence 10 captured by the camera 3 are determined. Therein, the image area pa preferably does not overlap with adjacent image areas and can for example be described as follows. pa = (x, y,W, H)T (6)
Herein, the position of an upper left corner coordinate is described by x and y. A width of the image area pa is described by W and a height of the image area pa is described by H. The index a describes an ordinal number for the linkage to a certain image 9 from the image sequence 10. The x and y coordinates are linked to a certain image area pa while the width W and the height H are the same for all of the image areas pa. The weight W can for example be preset with a parameter value of ten pixels and the height H can for example also be preset with a parameter value of ten pixels. The image areas pa are in particular separately examined for each image channel of the respective color space. In the further, the description is only presented for a respective intensity value l(x,y) of one of the images 9 of the image sequence 1 0.
In order to determine if the respective image area pa is suitable to be incorporated in the noise model NFL, the respective image area pa is examined with respect to a first selection criterion 1 1 and a second selection criterion 12. The first selection criterion 1 1 implies that an average intensity value δα' over a time t has a variation rate below a predetermined variation limit value ξ . This can be mathematically described as follows. The average intensity value δα' can be described as follows. W-l H-l
f (7)
Whereby the first selection criterion 1 1 can then be defined as follows.
Figure imgf000014_0001
Wherein it is stated by the equation (8) that the maximum variation rate of the average intensity value δα' over a number A of all of the images 9 is to be less than the
predetermined variation limit value ξ , which is for example one to two percent of the value range of intensity values of the images 9. By the number A, for example, ten images 9 are specified. For this purpose, an overall average intensity value δα } is subtracted from the average intensity value δα' . The image areas pa, which satisfy the condition from equation (8) of the first selection criterion 1 1 , are forwarded for examining the second selection criterion 12. For the second selection criterion 12, it is assumed that the noise in the noisy image lN substantially follows a Gaussian distribution. Thus, it can be assumed that if differences of locally corresponding intensity values of the images 9 do not have a distribution curve, which substantially corresponds to a Gaussian distribution, thus, the image area pa is not suitable with respect to the second selection criterion 12.
Thus, the second selection criterion 12 implies that differences Aj, of locally corresponding intensity values Vji should have a distribution curve over the time t, which shows an average value of 0 or else is formed symmetrically to an ordinate axis extending to the maximum of the distribution curve. The locally corresponding intensity values Vji are mathematically represented as follows. vjl (t) = I'(pax + j,pay + l) t e {l,...,A} (9)
Hereby, a vector is preset, in which the intensity values of a pixel of an image area are preset over the changing time and within an interval j and an interval I. Therein, the interval j corresponds to j e {0,...,W - 1} and the interval I therein corresponds to l e {0,...,H - l} . The differences Aj, of the locally corresponding intensity values Vj, are obtained by subtracting an average vector b.
Ajl (t) = vjl (t) ~ jvjl (b) te {l,...,A} (10)
-A b=l
Thereby, a sum of the distributions results, which is described as follows.
Figure imgf000015_0001
Herein, a distribution curve Λ is provided, which represents the distribution of the differences in particular by a discrete histogram. Herein, the maximum of the sum of the distributions is searched and it is examined if it is close to 0. If this is the case, then, the distribution curve Λ is examined with respect to its symmetry to the ordinate axis 13 extending to the maximum of the distribution curve Λ. The symmetry can simply be examined by subtracting an integral from the part of the distribution curve Λ disposed to the left of the ordinate axis 13 from an integral of a part of the distribution curve Λ disposed to the right of the ordinate axis 13. If the difference is less than a predetermined difference limit value, it can be assumed that the distribution curve Λ is symmetrical and thus corresponds a sum of Gaussian distribution.
Fig. 4 shows one of the images 9 of the image sequence 10. The images 9 in particular occur with the number A in the image sequence 10. In the images 9, the image areas pa are determined.
Fig. 5 shows the examined image areas pa of the images 9, which have been determined based on the image shown in Fig. 4. Thus, the hatched image areas 14 satisfy both the first selection criterion 1 1 and the second selection criterion 12, while the non-hatched image areas 15 do not satisfy at least one of the selection criteria 1 1 , 12.
Fig. 6 shows a first image 16 of the images 9, a second image 17 of the images 9 and a third image 18 of the images 9. The images 16, 17, 18 are captured consecutively in time In the images 16, 17, 18, the locally corresponding image area pa is drawn. In the image area pa, first locally corresponding intensity values v1 ; second locally corresponding intensity values v2 and third locally corresponding intensity values v3 are drawn. Therein, the locally corresponding intensity values v1 ; v2, v3 are in particular preset as a vector.
Fig. 7 shows a first distribution curve Λ1 ; a second distribution curve Λ2 and a third distribution curve Λ3, which are formed as follows. = {vl - mean{vl )) (12) Λ2 = (v2 - mean(v2)) (13) K = (v - mean(v )) (14)
A fourth distribution curve 19 is described by equation (15).
Figure imgf000016_0001
Thus, on the ordinate of the diagram of Fig. 7, a standard deviation of the respective locally corresponding intensity values Vj, is in particular indicated.
Fig. 8 now shows the images 9 of the image sequence 10, based on which the noise model NLF is generated. In particular, the respective noise model NLF is determined for each color channel of the images 9 of the image sequence 10. Therein, the color channels can for example be present in an RGB color space or else in a YUV color space. For simplicity, the creation of the noise model NLF is explained in the following based on a color channel of the RGB color space, wherein the method is also applicable to the YUV color space. The noise model NLF is mathematically described as follows.
NLF(I(x, y))→a (16)
Wherein a noise characteristic σ or a standard deviation of the noise of a determined intensity value l(x,y) of an image of the images 9 of the image sequence 10 describes as an output. In particular, the noise model NLF is determined for the various color channels, thus for example a red color channel, a blue color channel and a green color channel. It is shown, how the intensity values of the pixels from the images 9 over the time can be combined to a vector v. Fig. 9 shows a diagram, in which an intensity value 21 of a red color channel of the images 9 of the image sequence 10 is plotted on an abscissa 20. The noise characteristic σ or a standard deviation of the noise is registered on an ordinate 22. Thus, the noise model NLF is described for values i with the range of values of i e {0,1,...,255}. Thereby, thus, for each intensity value l(x,y) contained in the images 9, a noise characteristic σ in the form of the standard deviation is provided by the noise model NLF.
This is shown based on Fig. 10. Thus, based on the noise model NLF, a first noise characteristic σι , a second noise characteristic σ2 and a third noise characteristic σ3 are here derived from the noise model NLF. The first noise characteristic σι corresponds to a high standard deviation, the second noise characteristic σ2 corresponds to a medium standard deviation and the third noise characteristic σ3 corresponds to a low standard deviation.
Now, the noise model NLF is determined as follows. Each discrete value i of NLF (i) is calculated from the selected intensity values over the time. Herein, i=l(x,y) and ie {0,1,...,255}. The noise model NLF is calculated over all of the pixels, which are within the previously selected image areas pa and have the same average intensity value over the time t. A standard deviation over all of these pixels provides the noise model NLF in the position i. This is schematically explained based on Fig. 8. The vector v also has an average value mv. It is mathematically described as follows, wherein n is a number of the images 9 of the image sequence 10.
Figure imgf000017_0001
It is assumed that after the image areas pa satisfy the selection criteria 1 1 , 12, thus do not move, and lighting variations are not to be expected, they are suitable for modeling the noise. Thus, the variations of the intensity values in the vector v are substantially only caused by the noise. For the estimation of the noise model NLF, the noise model NLF is formed over all of the pixels of the images 9, in which the average value mv of the vector v is the same over the time and thus i=mv is true. Thus, the noise model NLF has for example been modeled for an intensity value of i=100 for all of the vectors v, which have the average value of mv=100. Pixels of the image 9, which are encompassed by the noise model NLF for the intensity value of i=100, are marked in Fig. 1 1 . Thus, all of the pixels are characterized there by an average intensity value of 100. Analogously to Fig. 1 1 , Fig. 12 shows the noise model NLF for an intensity value of i=200.
Due to the noise, it can occur that some values of the range of values of the intensity values of the images 9 do not occur in the images 9 and thus the noise model NLF either cannot be generated for these values. Thus, this is for example the case in Fig. 9 for intensity values below i=35. Then, the noise model NLF can for example be approximately estimated for these values or else simply be left empty.
Assuming k is the number of pixels with the same average intensity value (they are the emphasized pixels of Fig. 1 1 or Fig. 12), then, all of the k vectors can be combined in an overall vector γ,. ri = (v1 ,v2 ,...,vk ) (18)
A vector of the noise differences η can then be calculated as follows.
Figure imgf000018_0001
The noise model NLF is then determined as follows.
Figure imgf000018_0002
Herein, the noise model NLF is a normalized standard deviation of the noise differences η. Furthermore, a normalization constant m is indicated, which can be described as follows.
Figure imgf000018_0003
Additionally, it can be provided that the noise model NLF is again determined or updated depending on an update criterion. Herein, the noise model NLF is then determined as follows.
NLF(i) = _^~ i (22) Wherein, three scalars T0, T^ and T2 are introduced here, which together are described as TQ. If new noise differences r, are then to be added, it is only required to form the sums of all of the scalars described by TQ and to newly determine the noise model NLF. The scalar Tc can be mathematically described as follows.
Γβ =∑ r( Or, Vare {0,1,2} (23)

Claims

Claims
Method for generating a noise-reduced image (lD), in which a noisy image (lN) is captured by a camera (3) of a motor vehicle (1 ) and the following steps are performed:
a) capturing an image sequence (10) of images (9) by the camera (3);
characterized by
b) determining a noise model (NLF) based on multiple images (9) of the image sequence (10);
c) determining a weighting value (wxy) for weighting the noisy image (lN) as well as at least one image (lB) depending on the noisy image (lN) based on the noise model (NLF);
d) generating a blurred image (lB) based on the noisy image (lN);
e) generating a weighted blurred image by weighting the blurred image (lB) depending on the weighting value (wxy);
f) generating a weighted noisy image by weighting the noisy image (lN) depending on the weighting value (wxy); and
g) generating the noise-reduced image (lD) based on the weighted blurred image and the weighted noisy image.
Method according to claim 1 ,
characterized in that
the noise-reduced image (lD) is generated by a linear combination, in particular by summation, of the weighted blurred image and the weighted noisy image.
Method according to claim 1 or 2,
characterized in that
before step b) at least each one image area (pa) of the images (9) of the image sequence (10) is selected based on at least one predetermined selection criterion (1 1 , 12) and the noise model (NLF) is then determined in step b) only from the respective image areas (pa), in particular the steps c) to g) are only performed for the respective image areas (pa).
4. Method according to claim 3,
characterized in that
a variation rate of an average intensity value ( δα' ) of locally corresponding intensity values of the images (9) of the image sequence (10) over the temporal progress of the image sequence (10) is determined as the selection criterion (1 1 ) of the image areas (pa) and the image areas (pa) are selected if the variation rate is less than a predetermined variation limit value ( ξ ) .
5. Method according to claim 3 or 4,
characterized in that
a distribution curve (Λ) of differences (λ ,) of respective locally corresponding intensity values (vj,) of the images (9) of the image sequence (10) with an average intensity value of the locally corresponding intensity values (vj,) of the images (9) is determined as the selection criterion (12) of the image areas (pa) and the image areas (pa) are selected if the distribution curve (Λ) of the differences ( Aj,) has an average value of zero and/or the distribution curve (Λ) of the differences ( Aj,) is symmetrical to the ordinate axis (13) extending to the maximum of the distribution curve (Λ).
6. Method according to any one of the preceding claims,
characterized in that
the steps b) to g) are each performed for at least one intensity value of a color channel of the noisy image (lN), in particular an RGB color channel of the noisy image (lN) or a YUV color channel of the noisy image (lN).
7. Method according to any one of the preceding claims,
characterized in that
the noise model (NLF) is determined in step b) based on intensity values (I) of locally corresponding pixels (v) of the images (9) of the image sequence (10).
8. Method according to claim 7,
characterized in that
the locally corresponding pixels (v) of the images (9) of the image sequence (10) are combined with respect to their respective average intensity values in determining the noise model (NLF) in step b), and by the noise model (NLF) a noise characteristic (σ), in particular a standard deviation, is respectively provided for at least one value (i) from a range of values of the intensity values ( I) of the images (9) of the image sequence (10).
9. Method according to claim 8,
characterized in that
an average value (mv) is respectively formed from the locally corresponding pixels (v) of the images (9) of the image sequence (10), and the noise model (NLF) is determined depending on all of the corresponding pixels (v) with the same average value (mv).
10. Method according to any one of the preceding claims,
characterized in that
after step a) a further image is captured by the camera (3) depending on an update criterion and the noise model (NLF) is again determined additionally considering the further image.
1 1 . Method according to any one of the preceding claims,
characterized in that
the weighting value (wxy) is determined depending on a magnitude of an edge of the noisy image (lN), and the blurred image is the less severely taken into account in weighting in step e) the higher the magnitude of the edge is and/or the noisy image (IN) is the more severely weighted in weighting in step f) the higher the magnitude of the edge is.
12. Method according to any one of the preceding claims,
characterized in that
step a) and/or step b) are performed at a time, at which an automatic exposure setting of the camera (3) is performed and/or at which the motor vehicle (3) is statically operated.
13. Method according to any one of the preceding claims,
characterized in that
the blurred image is determined in step d) by a convolution of the noisy image (lN) with a Gaussian kernel, in particular with an average value of zero.
14. Camera system (2) for a motor vehicle (1 ), with a camera (3) and an evaluation unit (4), which is formed to execute a method according to any one of the preceding claims.
15. Motor vehicle (1 ) with a camera system (2) according to claim 14.
PCT/EP2017/055328 2016-03-07 2017-03-07 Method for generating a noise-reduced image based on a noise model of multiple images, as well as camera system and motor vehicle WO2017153410A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102016104043.6 2016-03-07
DE102016104043.6A DE102016104043A1 (en) 2016-03-07 2016-03-07 Method for generating a noise-reduced image based on a noise model of a plurality of images, and camera system and motor vehicle

Publications (1)

Publication Number Publication Date
WO2017153410A1 true WO2017153410A1 (en) 2017-09-14

Family

ID=58265949

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/EP2017/055328 WO2017153410A1 (en) 2016-03-07 2017-03-07 Method for generating a noise-reduced image based on a noise model of multiple images, as well as camera system and motor vehicle

Country Status (2)

Country Link
DE (1) DE102016104043A1 (en)
WO (1) WO2017153410A1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102018204881A1 (en) * 2018-03-29 2019-10-02 Siemens Aktiengesellschaft A method of object recognition for a vehicle with a thermographic camera and a modified noise filter
CN111028171A (en) * 2019-12-06 2020-04-17 北京金山云网络技术有限公司 Method, device and server for determining noise level of image

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5081692A (en) * 1991-04-04 1992-01-14 Eastman Kodak Company Unsharp masking using center weighted local variance for image sharpening and noise suppression
WO2002027656A2 (en) * 2000-09-29 2002-04-04 Hewlett-Packard Company Selective smoothing and sharpening of images by generalized unsharp masking
US20140126808A1 (en) 2012-11-02 2014-05-08 Board Of Regents, The University Of Texas System Recursive conditional means image denoising

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102005033851A1 (en) * 2005-07-20 2007-02-01 Sci-Worx Gmbh Digital image noise artifacts reducing method, involves reducing noise artifacts from actual portion of images to be processed depending on adaptive adjusted actual global variance value
US8149336B2 (en) * 2008-05-07 2012-04-03 Honeywell International Inc. Method for digital noise reduction in low light video

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5081692A (en) * 1991-04-04 1992-01-14 Eastman Kodak Company Unsharp masking using center weighted local variance for image sharpening and noise suppression
WO2002027656A2 (en) * 2000-09-29 2002-04-04 Hewlett-Packard Company Selective smoothing and sharpening of images by generalized unsharp masking
US20140126808A1 (en) 2012-11-02 2014-05-08 Board Of Regents, The University Of Texas System Recursive conditional means image denoising

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
CAO YANG ET AL: "A novel segmentation based video-denoising method with noise level estimation", INFORMATION SCIENCES, vol. 281, 10 October 2014 (2014-10-10), pages 507 - 520, XP029035207, ISSN: 0020-0255, DOI: 10.1016/J.INS.2014.05.031 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102018204881A1 (en) * 2018-03-29 2019-10-02 Siemens Aktiengesellschaft A method of object recognition for a vehicle with a thermographic camera and a modified noise filter
CN111028171A (en) * 2019-12-06 2020-04-17 北京金山云网络技术有限公司 Method, device and server for determining noise level of image
CN111028171B (en) * 2019-12-06 2023-04-18 北京金山云网络技术有限公司 Method, device and server for determining noise level of image

Also Published As

Publication number Publication date
DE102016104043A1 (en) 2017-09-07

Similar Documents

Publication Publication Date Title
US11244432B2 (en) Image filtering based on image gradients
US9019402B2 (en) Dynamic range extension by combining differently exposed hand-held device-acquired images
EP1924966B1 (en) Adaptive exposure control
JP6159298B2 (en) Method for detecting and removing ghost artifacts in HDR image processing using multi-scale normalized cross-correlation
KR101518722B1 (en) Image enhancement processing method and apparatus for distortion correction by air particle like fog
US11017509B2 (en) Method and apparatus for generating high dynamic range image
TWI462054B (en) Estimation Method of Image Vagueness and Evaluation Method of Image Quality
Park et al. Generation of high dynamic range illumination from a single image for the enhancement of undesirably illuminated images
CN112200848B (en) Depth camera vision enhancement method and system under low-illumination weak-contrast complex environment
KR20140118031A (en) Image processing apparatus and method thereof
Kinoshita et al. Automatic exposure compensation for multi-exposure image fusion
Kinoshita et al. Pseudo multi-exposure fusion using a single image
CN115883755A (en) Multi-exposure image fusion method under multi-type scene
WO2017153410A1 (en) Method for generating a noise-reduced image based on a noise model of multiple images, as well as camera system and motor vehicle
WO2014206503A1 (en) Automatic noise modeling for ghost-free image reconstruction
Han et al. Automatic illumination and color compensation using mean shift and sigma filter
van Beek Improved image selection for stack-based hdr imaging
CN109961422B (en) Determination of contrast values for digital images
EP3547251B1 (en) Dynamic range extension of partially clipped pixels in captured images
Vavilin et al. Recursive HDR image generation from differently exposed images based on local image properties
CN112258434A (en) Detail-preserving multi-exposure image fusion algorithm in static scene
Yao et al. Half-quadratic regularization based de-noising for high dynamic range image synthesis
Kimura et al. High dynamic range imaging based on camera response estimation and multi-gradients fusion
Park et al. High dynamic range image acquisition from multiple low dynamic range images based on estimation of scene dynamic range
JP2019140674A (en) Image processing apparatus, image processing method, and program

Legal Events

Date Code Title Description
NENP Non-entry into the national phase

Ref country code: DE

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17709936

Country of ref document: EP

Kind code of ref document: A1

122 Ep: pct application non-entry in european phase

Ref document number: 17709936

Country of ref document: EP

Kind code of ref document: A1