WO2008084585A1 - 画像処理装置および撮像装置 - Google Patents

画像処理装置および撮像装置 Download PDF

Info

Publication number
WO2008084585A1
WO2008084585A1 PCT/JP2007/070770 JP2007070770W WO2008084585A1 WO 2008084585 A1 WO2008084585 A1 WO 2008084585A1 JP 2007070770 W JP2007070770 W JP 2007070770W WO 2008084585 A1 WO2008084585 A1 WO 2008084585A1
Authority
WO
WIPO (PCT)
Prior art keywords
component
noise
signal
signal level
image
Prior art date
Application number
PCT/JP2007/070770
Other languages
English (en)
French (fr)
Inventor
Takahiro Saito
Yuki Ishii
Nobuyuki Watanabe
Original Assignee
Kanagawa University
Olympus Corporation
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 Kanagawa University, Olympus Corporation filed Critical Kanagawa University
Publication of WO2008084585A1 publication Critical patent/WO2008084585A1/ja
Priority to US12/498,845 priority Critical patent/US8189940B2/en
Priority to US13/477,225 priority patent/US8538186B2/en

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/10Cameras or camera modules comprising electronic image sensors; Control thereof for generating image signals from different wavelengths
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/80Camera processing pipelines; Components thereof
    • H04N23/81Camera processing pipelines; Components thereof for suppressing or minimising disturbance in the image signal generation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/80Camera processing pipelines; Components thereof
    • H04N23/84Camera processing pipelines; Components thereof for processing colour signals
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N9/00Details of colour television systems
    • H04N9/64Circuits for processing colour signals
    • H04N9/646Circuits for processing colour signals for image enhancement, e.g. vertical detail restoration, cross-colour elimination, contour correction, chrominance trapping filters
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N5/00Details of television systems
    • H04N5/14Picture signal circuitry for video frequency region
    • H04N5/21Circuitry for suppressing or minimising disturbance, e.g. moiré or halo

Definitions

  • the present invention relates to noise included in an image signal, and particularly to noise that depends on the signal level. Background
  • a, b and c are terms and are given to the chain.
  • the simple term is used in the future.
  • the accuracy of estimating the amount of noise was inferior.
  • the frequency characteristics of filtering are controlled from the amount of noise. This filtering is performed in the region where the amount of noise is estimated to be large from the signal level in order to make the flat part and the edge part of the image equal without ⁇ 3 ⁇ 4 ⁇ . Edge force s deteriorated and the preservation of the original signal was low. Furthermore, it cannot cope with the color noise generated between each color signal.
  • JP2001-175843 ⁇ separates the input signal into 3 ⁇ 4 signal and color 3 ⁇ 4 ⁇ and obtains the edge bow based on these signals.
  • the wording is skeletonized, and this ossification process is performed fixedly regardless of the signal level.
  • skeletal defects cannot be controlled well, and there is a possibility that color noise components remain or the original signal deteriorates.
  • JP2006-023959 A proposes to obtain an allotment of hue in a noise coring ring based on the noise model, and to divide the area into regions and apply a noise model to the region. Furthermore, JP2006-101006 A uses the average value of the pixel of interest and its surroundings to estimate the & level for noise coring processing based on the noise model, and calculates the coring threshold based on the noise model. Furthermore, it proposes to segment the image again according to the class of pixels belonging to the threshold and apply the noise model to the area IJ. Disclosure of the invention
  • the conventional average value is calculated as follows: If the noise energy is large, the noise is reduced, and the noise amount is removed by the average filter. Finale is needed.
  • this ⁇ there is a signal that cannot be applied correctly due to the influence of the pixel level in the vicinity of the edge, and the correct constant value cannot be obtained due to the influence of the neighboring pixel level. .
  • segmentation according to the noise model class does not consider the structural elements of the image, so there is a ⁇ that does not show the effect of waiting for the behavioral power near the edge.
  • the present invention has been made in view of such conventional St. technique 3 ⁇ 43 ⁇ 4 ⁇ .
  • the purpose is to reduce noise components that depend on the noise more effectively and to produce high-quality noise.
  • the image signal is converted into a first component that is a skeletal component
  • the component separation means that includes multiple components, including the signal level unit that acquires the signal level of the first component, and the noise parameter based on the signal level of the first component.
  • the noise of the second component can be made effective. High-quality noise can be avoided. Since the skeletal component is separated as the first component, the skeletal component is not inferior due to noise.
  • the S ⁇ ® image signal is compared with the pixel of interest in the vicinity of the pixel of interest.
  • FIG. 1 is a system configuration diagram of the apparatus according to the first embodiment.
  • Figure 2 shows the signal level-one noise variance model after processing of the Digitano W symbol.
  • FIG. 3 is a flowchart corresponding to the processing from the component component to the component in the first mode.
  • Figure 4 shows the relationship between the signal level immediately after A / D conversion and the noise variance.
  • Figure 5 is a special figure showing the input / output characteristics of the Digitano W code processor.
  • Fig. 6 is a diagram for explaining the soft decision threshold.
  • FIG. 7 is a diagram showing the noise that is reduced by noise.
  • FIG. 8 is a system configuration diagram of an imaging apparatus according to the second embodiment.
  • FIG. 9 is a flowchart corresponding to the processing from the component separation unit to the processing unit in the second embodiment.
  • FIG. 10 is a system configuration diagram of an imaging apparatus according to the third embodiment.
  • FIG. 11 is a flowchart corresponding to the key from the component separation unit to the component synthesis processing unit in the third type.
  • FIG. 12 is a system configuration diagram of the apparatus according to the fourth embodiment.
  • Figure 13 shows the signal level 1 noise variance model before digitalo W signal processing.
  • FIG. 14 is a system configuration diagram of an apparatus according to the fifth embodiment.
  • FIG. 15 is a diagram showing an example of a one-dimensional image signal having no correlation between color components.
  • FIG. 16 is a diagram showing an example of a one-dimensional image signal having a correlation between color components.
  • FIG. 17 is a system configuration diagram of the system according to the sixth ⁇ state.
  • FIG. 18 is a system configuration diagram of the tree according to the seventh aspect.
  • FIG. 19 is a diagram showing an example of a one-dimensional image signal having no correlation between touching minutes.
  • FIG. 20 is a diagram showing an example of a one-dimensional image signal correlated with “ ⁇ minutes”.
  • Figure 21 shows the flow chart corresponding to from the component to the component M in the seventh 3 ⁇ 4 ⁇ state.
  • FIG. 22 is a system configuration diagram of the imaging apparatus according to the eighth embodiment. BEST MODE FOR CARRYING OUT THE INVENTION
  • FIG. 1 is a system configuration diagram of an imaging apparatus according to the first embodiment of the present invention.
  • represents a data signal ⁇
  • represents a control If word ⁇ , and each processing unit is controlled by the system controller 1000.
  • the original image is converted into an electrical signal through a lens 1001 and a solid 3 ⁇ 4 element 1002, such as a CCD or CMOS, and the CDS (Correlated Double Sampling) / ⁇ ] Sampling, analog gain reduction, etc. are performed. Thereafter, the A / D conversion unit 1004 converts it to a digital M code, and the digital W code key unit 1005 converts it to a predetermined reverse image signal source image signal I).
  • the fixed element 1002 is a monochromatic mouth element
  • the image signal I is a monochromatic mouth signal.
  • image signal I is divided into first component U and second component V.
  • the first component U is the flat component (the component that changes to the crunch) and the skeletal component of the image signal I including the edge component (the image structure), and the second component V is a fine texture-like component. This is the part of the first component U that includes the noise ⁇ and the image signal I, including the noise ⁇ . Since the original image signal I is a monochrome signal, the image signal I force U and V are both 3 ⁇ 4g components.
  • the noise variance ⁇ (noise 3 ⁇ 4 is obtained corresponding to the signal level of the first component U with reference to the signal level-noise distribution model or its approximation data shown in Fig. 2. Since the first component U is a skeletal component of the image signal I, the noise variance ⁇ obtained from the signal level of the first component U is approximately equal to the noise variance ⁇ contained in the image signal I. Then, the noise parameters Tl and ⁇ 2 corresponding to the noise variance ⁇ are set, and are sent to the noise processor 3002.
  • the noise parameters Tl and ⁇ 2 are expressed by, for example, the following equation (1):
  • is a coefficient, for example, 1/2. Details of the signal level and noise model will be described later.
  • the noise processing unit 3002 performs a soft decision threshold ⁇ correlation (coring processing S) for the second component V using the noise processing parameters Tl and ⁇ 2. Since the second component V includes noise equivalent to the noise included in the image signal I, this process can effectively iS ⁇ T the noise included in the second component V. Then, the obtained value is sent to the processing unit 3005 as the corrected second component V. Details of the soft decision threshold will be discussed later.
  • the first component U and the corrected second component V are synthesized at a predetermined ratio, for example, a ratio of 1: 1, and the game component I, in which noise is degenerated, is obtained with respect to the image signal I.
  • the game part ⁇ is made into a display element 2001 such as an LCD, and through a circuit 2002 to a storage medium 2003 composed of a flash memory or the like.
  • Figure 3 shows the flow chart corresponding to the components up to 3005.
  • the image signal I is divided into three components and the age of the image signal I is represented by a broken line, but here it is explained as being divided into two components.
  • step S 11 j ® image signal I is separated rr into first component u and second component V.
  • step S 12 the noise variance ⁇ is obtained based on the signal level of the first component U with reference to the signal level-noise variance model shown in FIG. Then, the noise parameters Tl and ⁇ 2 are set according to the noise variance ⁇ .
  • step S 1 3 soft decision threshold ⁇ ⁇ 3 ⁇ 4 using noise 0 «0 « parameter Tl, ⁇ 2 is applied to the second component V, and the noise contained in the second component V is printed and corrected. Obtain the second component V.
  • step S 14 the first component U and the corrected second component V are combined to obtain a sum ⁇ in which noise is deducted from the image signal I.
  • the image signal I is expressed as the sum of the first component U and the second component V as shown in the following equation (2).
  • the property of the first component U required as early is as follows: Bounded Variation Function Space (BV) composed of multiple "smooth brightness variation shells"
  • BV Bounded Variation Function Space
  • the energy of the first component U is defined by the TV (Total Variation) norm J (U) in the following equation (3).
  • the second component V which is separated from the i ⁇ f image signal I force, is the force affected by noise.
  • the first component U is hardly affected by noise, and the skeletal component without dulling the edge.
  • Force S is extracted.
  • addition-type components include literary separation using a linear filter, adaptive smoothing using a median filter, a morphological filter, etc.
  • the following shows an example of 1 ⁇ 2fe of the addition type separation that does not depend on the bounded 3 ⁇ 41) function.
  • Example 1 The first component U is the low-order component of the bottom development in the image signal I, and the second component V is the image of the first component U in the image signal I.
  • Example 2 The first component U is the result of median filtering of the image signal I, and the second component V is the image of the first component u of the image signal I.
  • Example 3 The first component U is the result of applying the morphing filter to the image signal I, and the second component V is the wisteria of the first component U that is 3 ⁇ 4 ⁇ ⁇ to the image signal I.
  • Example 4 The first component U is the result of performing reduction processing on the image signal I and further enlargement processing, and the second component V is the image of the first component U in the original image signal I.
  • First component U is the result of applying Bilateral Filter to image signal I
  • second component V is the wisteria of first component U in image signal I ⁇ 3 ⁇ 4
  • the male of the multiplication type separation will be described.
  • the function space of the logarithmic first component u in Eq. (7) is derived from multiple “sliding force 3 ⁇ 4 small shells of change” separated by unbounds, like the first component U of the additive separation model described above.
  • Raseki G is the space of the function expressed as the equation (9) by the fibrin function gl , g 2 , and its energy is defined as G Nonorm // V // G in the equation (10) .
  • V (x, y) d x8l (x, y) + d xS2 ( x , y), g V 82 e L D
  • the second component V separated by the image signal I is affected by noise. 1
  • the first component U is hardly affected by noise and the skeletal component without dulling the edges. (Appreciation ⁇ Image structure) Force extraction.
  • the noise contained in the image signal I is Gaussian noise
  • the noise variance ⁇ (noise S) increases in a quadratic curve with respect to the signal level Lo immediately after A / D conversion.
  • this is expressed as a quadratic function as disclosed in JP 2005-175718, the following equation (12) Force S can be obtained.
  • ⁇ , ⁇ , ⁇ are terms.
  • the noise variance ⁇ varies not only with the signal level, but also with the 3 ⁇ 4 and gain of the element.
  • Fig. 4 plots the noise variance ⁇ for three ISO sensitivities (gains) of 100, 200, and 400 related to the gain at a certain temperature t. Each curve has the form shown in Eq. (12), but the number depends on the ISO sensitivity related to the gain. In is g, and the signal level-one noise dispersion model is converted into a ⁇ ;
  • Digitano W signal processing unit 1005 force 12bit input signal is converted to 8bit output signal! 3 ⁇ 4 3 ⁇ 4
  • Digitano W signal processing unit 1005 has input / output signal as shown in Fig. 5.
  • L (12) represents the signal level immediately after A / D conversion
  • L (8) represents the signal level after processing the digitano W symbol.
  • the logic parameter setting 3004 obtains the noise variance ⁇ by referring to the signal level and the noise ⁇ [model shown in Fig. 2, and the noise Meter Tl, ⁇ 2 is set. For the sake of simplicity, instead of Fig. 2, you may refer to data that approximates Fig. 2 with a broken line.
  • the noise i ⁇ m 3002 is included in the second component V by the soft decision threshold theory described below (coring processing). ⁇ ⁇ ⁇ ⁇ through the noise.
  • the following signal processing is represented by the following equation (14).
  • FIG. 6 shows how the soft decision threshold processing is performed.
  • the soft decision threshold agglomeration is performed on the 17 Aimoto signal A shown in Fig. 6 (a)
  • the signal D force S shown in Fig. 6 (b) is obtained.
  • base B uses the moving average of Shin A.
  • the above-described judgment threshold fi ⁇ is applied to the second component V separated from the original image signal I.
  • the second component V is the J ⁇ lli image signal I force, excluding the first component, which is the skeletal component.
  • Fig. 7 shows the separation of the above, since the noise «2; 1 »is shown.
  • Figure 7 (a) shows the ig® image signal I expressed in 17 fire sources and contains noise.
  • Figure 7 (b) shows the first component U separated from the image signal I by the component 3001.
  • the first component U is a tf case component including the normal job and the edge component, and contains almost no noise.
  • Figure 7 (c) shows the second component V obtained from the component 3001.
  • the second component V is a portion of the first component U in the image signal I, and includes noise equivalent to the noise contained in the S image signal I.
  • Figure 7 (d) shows the result of soft decision threshold for the second component V.
  • the corrected second component V and noise are reduced.
  • the noise variance ⁇ is obtained using the signal level of the first component U.
  • the noise ⁇ ⁇ was obtained by referring to data similar to the signal level-one noise variance model shown in Fig. 2 by a polygonal line.
  • FIG. 7 (e) shows the result of synthesizing the first component U of FIG. 7 (b) and the corrected second component V of FIG. 7 (d) in a 1: 1 ratio by the composition processing unit 3005.
  • the sum ⁇ preserves the edge component of the image signal I while noise is reduced.
  • the mixing ratio at the time of synthesis can be set arbitrarily.
  • the noise contained in the image signal I can be accurately estimated. it can.
  • the second component V which is the component of the first component U in the image signal I, separates noise equivalent to that contained in the image signal I. Therefore, based on the signal level of the first component U, set the noise parameter Tl, ⁇ 2 for noise of the second component V, and use this to set the noise of the second component V If this is done, the noise of the second component V can be effectively reduced, and high-quality noise can be reduced. Since the edge component is included in the first component, which is a skeletal component, the edge is not degraded by the noise treatment.
  • any noise processing can be used as long as the signal processing depends on the signal level.
  • the threshold value and the second component V according to the noise parameter Tl, ⁇ 2 And the soft decision threshold that corrects the second component V can be used.
  • any of the above-described addition type separation and multiplication type separation may be used. In addition to separation into two components, it may be divided into three or more components. Les.
  • the noise from the 2nd to the Nth component is not limited to the above-mentioned judgment threshold threshold collection, and if it is a process that produces noise based on the signal level, it is possible. W
  • FIG. 8 is a system configuration diagram of the imaging apparatus according to the second embodiment, and the solid-state imaging device 1012 is a full-color device having a color filter array arranged on the front surface.
  • the fixed image element 1012 may be either W3 ⁇ 4 or multiple.
  • ⁇ minutes «3008 separates the first component U into ⁇ minutes (R component, G component, B component) and obtains the signal level for each minute.
  • the noise iSMM parameter setting 3004 refers to the signal level-one noise variance model shown in Fig. 2 or its approximate data based on the ⁇ signal level, and obtains the noise variance ⁇ in J3 ⁇ 4. Then, the noise parameter Tl, ⁇ 2 for each machine is set to a value corresponding to the noise variance ⁇ of the corresponding ⁇ , for example, a value proportional to the noise variance ⁇ .
  • the second component V obtained in the component portion « 33001 is divided into ⁇ minutes.
  • Fig. 9 is a flowchart showing the process up to the component separation 3001 force generation processing unit 3005.
  • the process when separating into three components is indicated by ⁇ .
  • step S 21 the image signal I is divided into a first component U and a second component V.
  • step S 22 the first component U is separated into ⁇ minutes, and the signal level for each ⁇ is obtained.
  • step S 23 based on the signal level of each color component, the signal level—noise dispersion model shown in FIG. Then, according to the variance ⁇ of the noise of each tactile component, the noise parameter Tl, ⁇ 2 for each ⁇ component of the second component V is set.
  • step S 2 the second component V is separated into each minute tf.
  • step S 2 5 for each ⁇ component of the second component V, the soft decision threshold is independently set to ⁇ using the corresponding minute noise « ⁇ theoretical parameters Tl, ⁇ 2, and the soft decision threshold ⁇
  • the corrected value is the corrected second component V.
  • step S 26 the first component U and the corrected second component V are combined to obtain a sum ⁇ in which noise is reduced with respect to the image signal I.
  • the image signal I force S includes a plurality of ⁇ components: Noise can be performed.
  • FIG. 10 is a system configuration diagram of the device related to the third ⁇ ⁇ separation. Components that are the same as those in the previous embodiment are given the same reference numerals, and descriptions thereof are omitted.
  • a Y / C separation unit 3018 is provided instead of the color component separation unit 3008.
  • the first component U is separated into boat component Y and color components Cb and Cr. Get Y signal level.
  • the conversion from the R component, G component, and ⁇ component to the luminosity ⁇ , color integrity Cb, and Cr constituting the first component U is performed by the following equation (15).
  • step 3004 the noise level ⁇ of each color component is obtained according to the signal level of component Y (refer to the signal level / noise distribution modeler or its approximation data shown in Fig. 2). Then, the noise SML theoretical parameters Tl and ⁇ 2 for each part are set to values corresponding to the noise variance ⁇ of the corresponding color component, for example, values proportional to the noise variance ⁇ .
  • the second component V obtained with component »3001 is divided into each component.
  • the noise part 3002 performs the above EiC decision threshold independently on ⁇ 3 ⁇ 4 ⁇ , and makes noise included in each ⁇ minute.
  • Figure 11 is a flowchart showing the processing from the component separation city 3001 to the component integration unit 3005.
  • the process of dividing into three components is represented by ⁇ , but here we will explain by the age divided into two components.
  • step S 3 the image signal I is divided into a first component U and a second component V.
  • step S 3 2 the first component U is separated into Y / C components, and the signal level of brightness Y is obtained.
  • step S 3 depending on the signal level of the key component Y of the first component U,
  • the noise variance ⁇ is obtained by referring to the signal level 1 noise variance model or its approximated data, and the noise parameters Tl and ⁇ 2 are set for each touch of the second component V accordingly.
  • step S 3 the second component V is separated into each color component.
  • step S 3 5 the soft decision threshold using the noise parameter Tl, ⁇ 2 is independently applied to each component of the second component V to obtain the corrected second component V. .
  • step S 36 the first component U and the corrected second component V are combined to obtain a total ⁇ 3 ⁇ 43 ⁇ 4 of the noise ⁇ for the image signal I.
  • the image signal I includes a plurality of color components, it is possible to perform high-quality noise processing while suppressing the deterioration of the edge components as in the case of the first » ⁇ . It can be carried out.
  • any one of the component separation methods may be used.
  • the process of removing the noise of the second component V is not limited to the soft decision threshold, but any signal process can be applied as long as it is a signal level-dependent signal process.
  • the separation component ⁇ is separated from the first component, and the noise parameters Tl and ⁇ 2 are set according to the signal level of the key component ⁇ .
  • the signal level of the signal obtained by interpolating the G component that most reflects the separated component ⁇ may be used.
  • FIG. 12 is a system configuration diagram of the device according to the fourth saddle configuration.
  • the previous ⁇ state The same reference numerals are assigned to the common components and the description thereof is omitted.
  • the component 3001 to the component synthesis processing unit 3005 are arranged after the digitano M code processing unit 1005.
  • these elements are placed in the previous stage of the Digitano W signal processing unit 1005, and component separation and female noise processing are performed on the signal immediately after A / D conversion (the signal before digital signal processing).
  • the signal obtained through the color filter array is demosaiced by the interpolation processing unit 4000.
  • the noise reduction processing parameters Tl and ⁇ 2 are set to color using the signal level-one noise variance model To do.
  • the signal is separated and noise ⁇ ⁇ ⁇ is performed before the digital signal processing. Therefore, the signal level-to-noise variance model is not the model shown in Fig. 2, but Refer to the model represented by (12) and (13) and insult the noise ⁇ .
  • the noise of the second component V by the soft decision threshold ⁇ ⁇ ⁇ , the sum of the first component U and the corrected second component V) is the same as in the second state.
  • the original image signal I includes a plurality of color components, and, as in the first! ⁇ It can be performed.
  • component separation and noise processing are performed in the previous stage of the digital screen processing unit, the components of the skeleton component and noise component in the component component 3001 are improved. It should be noted that either ⁇ ⁇ H-type separation or multiplication-type separation of component separation may be used.
  • the noise that causes the noise of the second component V is not limited to the soft decision threshold, but any signal processing that depends on the signal level can be applied.
  • the noise variance is obtained based on the signal level of ⁇ and the noise parameters Tl and ⁇ 2 are set accordingly. (It is obtained by converting the entire machine immediately after A / D conversion instead of the level.)
  • the signal level of Si component ⁇ or the signal obtained by complementing the G component that most reflects component Y A signal level may be used.
  • the local intercolor correlations and numbers (hereinafter referred to as “interferences and numbers”) in the 3 ⁇ 3 near-in region show positive values close to 1.
  • the correlation coefficient that appears in the near-field region decreases because it includes pixels that have a decreased number of phases. As a result, it can be said that the smaller the value of the correlation coefficient, the more noise is contained in the image signal.
  • the noise parameters Tl and ⁇ 2 of the soft decision threshold are corrected according to the correlation between the touch minutes. This reduces the range of soft decision threshold in regions where the correlation is high and works to reduce the image signal I, while the soft decision threshold in regions where the correlation is low.
  • the width signal S is widened, and the image signal I can be flattened, enabling more accurate noise reduction.
  • FIG. 14 is a system configuration diagram of the device according to the fifth embodiment. Components that are the same as those in the previous embodiment are given the same reference numerals, and descriptions thereof are omitted.
  • the digital signal I (full color code) output from the digitano W signal processing unit 1005 is output to the color component separation unit 3022, the correlation calculation processing unit 3006, and the noise ifi processing unit 3002.
  • Minute 8 ⁇ 3022 separates the image signal I into each touch, and the representative word level calculation unit 3023 calculates a representative signal level in the vicinity of a predetermined pixel position by each separated color component.
  • nonlinear filters such as linear filters, median filters, non-linear such as bilateral filters? Filters, morphology filters, etc. can be used.
  • the noise variance ⁇ corresponding to the metabolic symbol level of each color component is obtained by referring to the signal level-one noise model shown in Fig. 2 or its approximation data. Then, the noise reduction processing parameters Tl and ⁇ 2 for each color component are set to values corresponding to the noise variance ⁇ of the corresponding color component, for example, values proportional to the noise variance ⁇ .
  • the correlation calculation processing unit 3006 calculates the number of relatives of each touch including the target pixel position and outputs the maximum number r which is the minimum value.
  • the parameter correction unit 3007 calculates a correction coefficient C based on the maximum / J4 difficulty number r, and multiplies this by the noise ⁇ » ⁇ 3 ⁇ 4 parameter Tl, ⁇ 2 for each touch, and noise management for each ⁇ Correct the parameters Tl and ⁇ 2.
  • the iS ⁇ L logic unit 3002 performs the soft decision threshold using the corrected noise ® »theoretic parameters Tl and ⁇ 2 independently for each color component, and generates noise depending on the signal level contained in the original image signal I. . Correlation calculation processing and noise correction of parameters Tl and ⁇ 2
  • Figures 15 and 16 show the relationship between the pixel position of the image signal I and the signal level of each touch. In the example shown in Fig. 15, there is no correlation between touches, and the signal level of each touch at each pixel position is
  • the maximum / J4th relation number r is -0.
  • the maximum / J4 number r is 0. 848117.
  • the noise factor Tl and ⁇ 2 are corrected by the following equation (16) using the maximum number r.
  • the noise determination parameter Tl, ⁇ 2 of the soft decision threshold ⁇ is corrected according to the correlation between the touching minutes, so that the component of the image signal I ») Among them, it is possible to discriminate between components due to noise and components due to the original image structure such as texture, and it is possible to prevent the latter component from being degraded by noise.
  • the noise ⁇ is calculated from the metabolic symbol level and the noise parameters Tl and ⁇ 2 are set accordingly.
  • the signal level of the TO component ⁇ or the signal level of the signal obtained by interpolating the G component can be used.
  • FIG. 17 is a system configuration diagram of the device according to the sixth ⁇ state.
  • the same reference numerals are assigned to the same configuration as the previous ⁇ state, and the description is omitted.
  • the signal obtained through the color filter array is demosaiced by the interpolation processing unit 4000 in order to obtain the full color ⁇ ⁇ signature image signal I).
  • the model shown in Fig. 2 is used as a model for the ⁇ of the signal level for each tactile. Instead, obtain the noise variance ⁇ for each ⁇ by referring to the model shown in Fig. 13 or Equations (12) and (13), and set the noise «processing parameters Tl and ⁇ 2 corresponding to this color. Set for each component. Noise «Processing parameters Tl and ⁇ 2 are set to values proportional to the noise variance ⁇ , for example.
  • the parameter correction unit 3007 performs the same processing as in the fifth example according to the equation (16).
  • the swipe processing unit 3002 performs soft decision threshold processing using the corrected noise processing parameters Tl and ⁇ 2. This is performed on the image signal I, and the signal level-dependent noise contained in the original image signal I is iS ⁇ T.
  • the color filter array signal extraction unit 5000 extracts a signal at a position corresponding to the arrangement of the color filter array from noise and a post-processing image signal, and outputs the image signal to the digital signal processing unit 1005.
  • the signal after the digital signal processing is sent to the display element 2001 and the storage media 2003 in the subsequent stage.
  • the noise parameters ⁇ 1 and ⁇ 2 of the soft decision threshold ffi are corrected according to the machine-to-machine correlation as in the fifth 3 ⁇ 4 »state, so that the i ⁇ image signal I Among the included 3 ⁇ 413 ⁇ 4 components, a component caused by noise and a component caused by the original structure of the image such as texture can be discriminated, and the latter component can be prevented from deteriorating due to noise.
  • the noise variance ⁇ is calculated from the level of the male power, and instead of the power level, the noise parameters ⁇ and ⁇ 2 are set accordingly. You can use the signal level of the component ⁇ or the signal level of the signal obtained by interpolating the G component that most reflects the ⁇ component Y.
  • FIG. 18 is a system configuration diagram of the apparatus according to the seventh embodiment. Constituent elements that are the same as in the previous embodiment are given the same reference numerals, and descriptions thereof are omitted.
  • the soft decision threshold ghost of the component component 3 ⁇ 4 ⁇ shown in the second ⁇ state and the noise S »theory corresponding to the minute correlation shown in the fifth and sixth ⁇ states Correct the parameters Tl and ⁇ 2.
  • the correlation calculation for touching is performed on the second component V obtained from the image signal I force.
  • Image signal I force The separated first component U is sent to the tactile divider 3008, and in the mechanical component ⁇ 133008, the signal level for each ⁇ is acquired.
  • the noise ffi ⁇ ® parameter summary 3004 as in the second 3 ⁇ 4 » based on the signal level of each machine, refer to the model of signal level-to-noise variance for each noise ⁇ . get. Then, the noise ifi ⁇ t logic parameters Tl, ⁇ 2 for each color component are set to a value corresponding to the noise variance ⁇ of the corresponding color component, for example, a value proportional to the noise variance.
  • the second component V separated from the J ⁇ image signal I force is separated by ⁇ and 3009, and then separated by the correlation calculation processing unit 3006.
  • the correlation calculation unit 3006 the second component V Correlation calculation processing of each color signal is performed on the minute V in the same manner as in the fifth and sixth 3 ⁇ 41 ⁇ states, and the minimum correlation coefficient is calculated.
  • the parameter correction unit 3007 calculates the noise (g »theoretical parameter Tl, correction coefficient C of ⁇ 2 based on the maximum / J4th M 3 ⁇ 4 r, and multiplies this by the noise sickle processing parameter Tl, ⁇ 2 of each color component. Correct the noise ifi ⁇ rational parameter Tl, ⁇ 2 for each ⁇ .
  • Noise i® »Sensor 3002 corrected noise for each touch of 2nd component V
  • the soft decision threshold ffi ⁇ using the physical parameter Tl, ⁇ 2 is applied independently to touch the noise contained in each color component.
  • the obtained signal is output as the corrected second component V.
  • the first component ⁇ and the corrected second component V are synthesized at a predetermined ratio, for example, a ratio of 1: 1, and the game component ⁇ in which noise is deceived with respect to the image signal I is obtained. obtain.
  • the second component V is one-dimensional. Since the second component V is the portion of the first component U from 3 ⁇ 4 to ⁇ in the image signal I, including the fine structure 3 ⁇ 4 ⁇ and noise, such as texture, the second component V is centered on zero It is a component that fluctuates as
  • Figures 19 and 20 show the relationship between the pixel position of the second component V and the signal level for each ⁇ .
  • the maximum / J4th base, the number r is -0.75148.
  • the maximum / J4 number r is 0.8873891.
  • the noise reduction processing parameters Tl and ⁇ 2 are respectively corrected by the equation (16) used in the fifth ⁇ state.
  • the noise correction factor C of the noise parameters Tl and T2 is (t r) 2 and the correction factor C is 3.0 ⁇ in Figure 19 and 3.0 ⁇ in Figure 20 0. 01590. Therefore, as shown in Fig. 19, the correlation between machine parts is corrected so that the range between the upper threshold and lower threshold of the soft decision threshold is increased, and as a result, the signal becomes more skeletonized. To do.
  • is corrected so that the width between the upper threshold and lower threshold of the soft decision threshold becomes smaller, and there is a slight correction that preserves the signal more. Done.
  • FIG. 21 is a flow chart corresponding to the processing from the component separation unit 3001 to the composition unit 3005.
  • the process of separating J ⁇ ® image signal I into three components is represented by stones, but here we will explain how to separate into two components.
  • step S 71 the image signal I is divided into a first component U and a second component V.
  • step S72 the first component U is separated into each color component, and the signal level of each color component is obtained.
  • step S73 based on the signal level of each color component, the noise variance ⁇ is obtained by referring to the signal level-noise variance model shown in Fig. 2 or its approximation data, and Set the noise parameter Tl, ⁇ 2.
  • step S74 the second component V is divided into ⁇ minutes.
  • step S 75 the correlation of each tactile component of the second component V is calculated to calculate the maximum / ”fourth number r.
  • step S 76 based on the maximum / J number r, the noise ⁇ ⁇ theoretical parameter Tl obtained in step S 73 and the correction coefficient C of ⁇ 2 are calculated.
  • step S 77 the noise parameters Tl and ⁇ 2 are corrected by multiplying the noise parameters Tl and ⁇ 2 by the correction coefficient C. And for each machine of second component V On the other hand, perform the soft decision threshold using the processed noise «processing parameter Tl, ⁇ 2 to obtain the corrected second component V.
  • step S 78 the first component U and the corrected second component V are combined to obtain a game component ⁇ in which noise is deceived from the J ⁇ ® image signal I.
  • the seventh 3 ⁇ 4 ⁇ ⁇ state even if the image signal I force S includes multiple touches, the deterioration of the edge component is suppressed as in the first state, and the high-quality noise ⁇ » ⁇ is reduced. It can be carried out.
  • ⁇ ⁇ H type separation or multiplication type separation of component separation may be used.
  • the processing for noise generation of the second component V is not limited to the soft decision threshold value, and any signal processing can be applied as long as it is signal level-dependent signal processing.
  • the noise variance ⁇ is obtained based on the signal level of each color signal of the first component U, and the noise parameters Tl and ⁇ 2 are set accordingly.
  • Power Each color signal of the first component ⁇ Instead of the first signal level, the signal level of the signal obtained by separating the first component from the first signal component of d3 ⁇ 43 ⁇ 4J3 ⁇ 4, or the G component that most reflects the luminance component Y is used. May be.
  • Figure 22 shows the system configuration of the »g device in the 8th ⁇ state.
  • the same components as those in the previous embodiment are denoted by the same reference numerals and description thereof is omitted.
  • the interpolation processing unit 4000 demos the signal obtained through the color filter array.
  • the noise parameter setting ⁇ 13 3004 is shown in Fig. 2 as a model of noise level-to-noise dispersion for each touch.
  • the corrected second component V output from the noise processing unit 3002 is combined with the first component U by the component combining processing unit 3005.
  • the color filter array signal extraction unit 5000 extracts a signal at a position corresponding to the arrangement of the color filter array from the sum ⁇ , and magnetizes the image signal to the Digitano W No. Mabe Fuji. According to the eighth ⁇ state, even if the original image signal I force S includes multiple mechanical elements ⁇ , the deterioration of the edge components is suppressed as in the first 3 ⁇ 4 » Noise ifi ⁇ can be done.
  • the noise reduction processing parameters Tl and ⁇ 2 in the soft decision threshold processing according to the correlation between the color components, among the fluctuation components included in the second component V, the noise-induced components and It is possible to discriminate a component caused by the original structure of an image such as a texture, and to suppress the latter component from being inferior due to noise.
  • the separation of the skeleton component and the noise component in the component 3001 is improved. It should be noted that either the ⁇ 3 ⁇ 43 ⁇ 4u arithmetic type separation or the multiplication type separation of the component ⁇ may be used. Further, the process of ifi ⁇ r removing the noise of the second component V is not limited to the soft decision threshold, and any signal can be applied as long as it is a signal level dependent signal.
  • the noise variance ⁇ is obtained based on the signal level of each color signal of the first component U, and the noise parameters Tl, ⁇ 2 are set accordingly.Each color of the first component ⁇ (It is obtained by converting the entire signal immediately after A / D conversion instead of the signal level of the signal.) 3 ⁇ 43 ⁇ 43 ⁇ 4 Minute signal level, or signal obtained by interpolating the G component that best reflects ⁇ The signal level may be used.

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)
  • Color Television Image Signal Generators (AREA)
  • Picture Signal Circuits (AREA)
  • Processing Of Color Television Signals (AREA)
  • Facsimile Image Signal Circuits (AREA)

Abstract

原画像信号Iを、骨格成分である第1成分Uと、原画像信号Iに対する第1成分Uの残差から得られる第2成分Vと、を含む複数の成分に分離し(S11)、第1成分Uの信号レベルに基づきノイズ低減処理パラメータT1、T2を設定し(S12)、ノイズ低減処理パラメータT1、T2を用いて第2成分Vのノイズを低減する(S13)。

Description

明 細 書 画舰職置およ U¾i ^置 娜扮野
本発明は、 画像信号に含まれるノイズの纖、 特に信号レベルに依存するノイズ の«に関する。 背景嫌
J P2001-157057 Aでは、 ノイズ量 Nと濃度値に変換した信号レ ベル Dの隱、を、 N=abcDにより関数化している。 a、 b、 cは 項であり、 鎖に 与えられる。 しかしながら、鍵ノイズ量は、 ¾時の? ½、 露光時間、 ゲインな どの要因により動的に変化するため、 簡な 項を用いるこの 去では、 ¾ ^時 のノィズ量に合^:た関数化に対応することができず、 ノィズ量の推定精度が劣つ ていた。 また、 ノイズ量からフィルタリングの周波数特性を制御するが、 このフィ ルタリングは画像上の平坦部、 ェッジ部を κ¾ϋすることなく同等に するため、 信号レベルからノイズ量が大と推定された領域にあるエッジ部力 s劣化してしまレヽ、 原信号の保存性が低かった。 さらに、 各色信号間に発生する色ノイズに関しては対 応できない。
また、 JP2001— 175843 Αでは、 入力信号を) ¾ 信号と色 ¾ί言号とに 分離し、 これらの信号に基づきエッジ弓娘を求め、 エッジ部!^の平坦 ¾ 域で色 言号に 骨ィ 理を行っており、 この 骨化処理は信号レベルに力かわらず固定 的に行われる。 しかしながら、 色ノイズ量は信号レベルにより異なるため、 骨化 麵を纖に制御することができず、 色ノィズ成分の残存や原信号の劣化など力発 生してしまう可能性がある。
そこで、 J P2006— 023959 Aでは、 ノイズモデ'ノレに基づくノイズコア リング鍾において、 色相の類献を求め、 領域に分割し、 領¾¾にノイズモデル を適用する 去を提案している。 さらに、 JP2006— 101006 Aでは、 ノ ィズモデルに基づくノイズコアリング処理にぉ、て、 &レベルを推定するために 注目画素とその周辺の平均値を用いてノイズモデルに基づくコアリングの閾値を求 め、 さらに、 閾値内に属する画素の類 ¾に応じて、 再度画像のセグメンテーショ ンを行レ、、 領鄉 IJにノィズモデルを適用する を提案している。 発明の開示
しかしながら、 従来の平均値を求める: であると、 特にノイズのエネノレギ一が 大きレ、画像にぉレ、ては、 そのノィズ量除レ、た値を平均値フィルタで除くために大き なカーネノレのフィノレタが必要となる。 この^^、 副ィ乍用として、 エッジ部近辺にお レ、て、 近傍の画素レベルの影響によって正しレ^ i定値が得られず、 ノイズモデルを 正しく適用することができなレヽ がある。 また、 ノイズモデルの類 ί に応じた セグメンテーションにおいても、 画像の構造要素を考慮していないので、 エッジ部 付近での挙動力溯待通りの効果を示さなレ、^^がある。
本発明は、 このような従来 St于の技 ¾¾ ^^を鑑みてなされたもので、 信号レべ ルに依存するノイズ成分をより効果的に纖し、 高画質なノイズ を行える ようにすることを目的とする。
本発明によれば、 画舰 β置において、 願像信号を、 骨格成分である第 1成 分と、
Figure imgf000005_0001
を含 む複数の成分に » る成分分離手段と、 歸己第 1成分の信号レベルを取得する信 号レベル取 段と、 編己第 1成分の信号レベルに基づきノィズ« ^理パラメ一 タを設定するパラメータ設定手段と、 ffilSノィズ 理パラメータを用レヽて嫌己 第 2成分のノイズを ί»τるノイズ 手段と、 を備える。
第 1成分の信号レベルに基づきノィズ パラメ→を設定し、 これを用レ、 て第 2成分のノィズを ifi^Tるようにすれば、 第 2成分のノィズを効果的に ることができ、 高画質なノイズ 理を撒することができる。 骨格成分は第 1 成分として分離されてレ、るので、 ノィズ により骨格成分が劣ィ匕することも ない。
また、 本発明の別の崎によれば、 複数の械分を含む顾像信号を とした 画像処«置において、 S^®像信号に対して、 注目画素とその近傍の^ §域での 信号レベルを得る信号レベル取 ^段と、 tilt己信号レベルに基づきノィズ ノ ラメータを設定するノ ラメータ設定手段と、 ΙΐίΙΕノィズ、 理ノ ラメータを用 レ、て tiHB画像信号のノィズを低減するノィズ低減手段と、 IB^画像信号に対し て、 注目画素を含 貝域での械分間の相隱数を演算する相難数演算手段 と、 嫌 目隱数に応じて tiif己ノィズ 理パラメータを補正するパラメータ補 正手段と、 を備免る。 械分間の相関に応じてノィズ 理パラメータを補正するようにしたことに より、 画像信号に含まれる麵成分のうち、 ノイズに起因する成分と、 テクスチャ のような画像本来の構造に起因する成分を判別することができ、 ノィズ に よって後者の成分が劣化するのを抑えることができる。 図面の簡単な説明
図 1は、 第 1の¾»態に係る 置のシステム構成図である。
図 2は、 デジタノ W言号処理後の信号レベル一ノイズ分散モデルである。
図 3は、 第 1の 態における、 成分分 »から成 部までの処理に 対応するフローチャートである。
図 4は、 A/D変換直後の信号レベルとノィズの分散の関係を示した特 図であ る。
図 5は、 デジタノ W言号処理部の入出力特性を示す特 図である。
図 6は、 軟判定閾 理を説明するための図である。
図 7は、 ノイズ によりノイズが赚される ίを示した図である。 図 8は、 第 2の実 »態に係る撮 «置のシステム構成図である。
図 9は、 第 2の実»態における、 成分分離部から成 理部までの処理に 対応するフローチヤ一トである。
図 1 0は、 第 3の実¾ ^態に係る撮^ ¾置のシステム構成図である。
図 1 1は、 第 3の魏形態における、成分分离鄉から成分合成処理部までの鍵 に対応するフローチヤ一トである。 図 1 2は、 第 4の実 態に係る 置のシステム構成図である。
図 1 3は、 デジタノ W言号処理前の信号レベル一ノィズ分散モデノレである。
図 1 4は、 第 5の実¾ ^態に係る 置のシステム構成図である。
図 1 5は、 色成分間に相関のない 1次元の画像信号の例を示した図である。
図 1 6は、 色成分間に相関のある 1次元の画像信号の例を示した図である。
図 1 7は、 第 6の^^態に係る |»¾置のシステム構成図である。
図 1 8は、 第 7の 態に係る樹«置のシステム構成図である。
図 1 9は、 触分間に相関のない 1次元の画像信号の例を示した図である。
図 2 0は、 "^分間に相関のある 1次元の画像信号の例を示した図である。
図 2 1は、 第 7の ¾ ^^態における、 成分分 から成^^ M部までの に対応するフローチヤ一トである。
図 2 2は、 第 8の ¾ ^態に係る撮 «置のシステム構成図である。 発明を実施するための最良の形態
以下、 樹寸図面を参照しながら本発明の雄形態について説明する。
<第 1の¾»維〉
図 1は本発明の第 1の実«態に係る撮^ ¾置のシステム構成図である。 図 1に おいて、 Λはデータ信号翁、 は制摘 If言号錄を表しており、 各々の処理部は、 システムコントローラー 1000で淛御される。
原画像は、 レンズ 1001を通して CCD、 CMOS等の固体 ¾素子 1002で電気信号に 変換され、 馳信号処理部 1003において CDS (Correlated Double Sampling)/^] サンプリング、 アナログゲインの讓等カ行われる。 その後、 A/D変換部 1004で デジタノ M言号に変換され、 デジタノ W言号鍵部 1005で所定の翻の画像信号 源 画像信号 I) に変換される。 第 1の 態においては、 固^ »素子 1002 はモ ノク口用 素子、 像信号 Iはモノク口信号である。
成分分 » 3001では、麵像信号 Iを第 1成分 Uと第 2成分 Vに分 る。 第 1成分 Uは、 平載分 ( かに変化する成分) とエッジ成分を含む顾像信号 I の骨格成分 (謝可^)画像構造) であり、 第 2成分 Vは、 テクスチャの様な細力い 構 it ^分とノイズを含む、 顾像信号 Iに ¾"Τる第 1成分 Uの 分である。 原 画像信号 Iがモノクロ信号であるので、 顾像信号 I力 分離される各成分 U、 V はいずれも) ¾g成分となる。
なお、 ここでは、説明を簡単にするために、願像信号 Iを 2成分に分 る場 合にっレ、て説明するが、 3成分以上に分离 f ることも可能である。 また、成分爐 の^去に ロ算型分離と乗算型分離があり、 その詳細にっレ、ては後财る。
ノイズ ί®^®パラメータ設¾¾ 3004では、 図 2に示す信号レベル一ノイズ分 散モデノレまたはその近似データを参照して第 1成分 Uの信号レベルに対応するノィ ズの分散 σ (ノイズ ¾ を取得する。 第 1成分 Uは顾像信号 Iの骨格成分である ので、 第 1成分 Uの信号レベルから得られるノィズの分散 σは、顾像信号 Iに含 まれていたノイズの分散 σにほぼ等しレヽ。 そして、 ノイズの分散 σに応じたノイズ パラメータ Tl、 Τ2を設定し、 ノイズ 理部 3002に{ する。 ノイズ 理パラメータ Tl、 Τ2は、 例えば、 次式 (1) :
Tl= ka T2= -ka
… ) により、 それぞれノイズの分散 σに比例する値に設定される。 kは係数で、 例えば 1/2である。 信号レベル一ノィズ獵モデルの詳細にっレ、ては後 る。
ノイズ 理部 3002 では、 第 2成分 V に対してノイズ 理パラメータ Tl、 Τ2 を用いた軟判定閾 {麵理 (コアリング処 S) を行う。 第 2成分 Vは 像 信号 Iに含まれていたノイズと同等のノイズを含むので、 この処理によれば、 第 2 成分 Vに含まれるノイズを効果的に iS^Tることができる。 そして、 得られた値を 補正第 2成分 Vとして成 理部 3005に する。 軟判定閾 の詳細に ついては後财る。
成 理部 3005では、 第 1成分 Uと補正第 2成分 Vを所定の比率、 例えば 1 : 1 の比率で合成し、 顾像信号 I に対してノイズが纖された合戯分 I,を得 る。 合戯分 Γは、 LCD等の表示素子 2001に繊されるとともに、 回路 2002 を介して、 フラッシュメモリ等で構成される記憶メディア 2003に される。 図 3は、 成分分 « 3001力ら成^ ^纏 3005までの麵に対応するフロー チヤ一トである。 図中参考のため顾像信号 Iを 3成分に分 る齢の麵を破 線で表してレ、るが、 ここでは 2成分に分 g る につレ、て説明する。
ステップ S 1 1では、 j ®像信号 Iを第 1成分 uと第 2成分 Vに分离 rrる。
ステップ S 1 2では、 第 1成分 Uの信号レベルに基づき、 図 2に示した信号レべ ル一ノイズ分散モデルまたはその近似データを参照してノイズの分散 σを取得す る。 そして、 ノイズ 理パラメータ Tl、 Τ2 をノイズの分散 σに応じて設定す る。
ステップ S 1 3では、 第 2成分 Vに対してノイズ ί«0«パラメータ Tl、 Τ2を 用いた軟判定閾 ίί^¾を行い、 第 2成分 Vに含まれているノイズを ί職し、補正第 2成分 Vを得る。
ステップ S 1 4では、 第 1成分 Uと補正第 2成分 V,を合成し、願像信号 Iに 対してノイズが赚された合 分 Γを得る。
続いて、成分分 3001において行う i i像信号 Iの成分^隹の詳細について 説明する。 麵像信号 Iから第 1成分 Uを抽出するとき、 その蘇として、 平職 分 ( ^かに変化する成分) やエッジ成分を含む顾像信号 Iの骨格成分が鹏さ れるような処理を用いる。 成分分離の 去には、 上記の通り、カロ算型分離と乗算型 分離があり、 以下各手法について順に説明する。
加算型分離
カロ算型分離では、顾像信号 Iが次式 (2)で示すように、 第 1成分 Uと第 2成分 Vの和として表される。
I = U + V
… ) ここでは、 まず、 有界変動関数とノルムを用いた分离 i^f去を説明する。
分離を行うために、 以下の文献に開示される A¾C変分モデル (Aujol-Aubert- Blanc— Feraud~Chambolle mode上) ¾ ffiレヽる。
文献 : Jean-Francois Aujol, uy Gilboa, Tony Chan & Stanley Osher, Structure— Texture Image Decomposition ~ Modeling, Algorithms, and Parameter Selection, International Journal of Computer Vision, Volume 67, Issue 1 (April 2006) Pages : 111 - 136 Year of Publication: 2006
ft ^早として求められる第 1成分 Uの性質としては、 不 によって区分さ れた複数の "滑らかな輝度変ィ匕の小部 貝 から構成された有界変動関数空間 BV (Bounded Variation Function Space)としてモデル化され、 第 1成分 Uのェネル ギ一は次式 (3)の TV (Total Variation)ノルム J (U)で定義される。
Figure imgf000011_0001
- (3) 一方、 式 (2)中の第 2成分 Vの関 間は 関輕間 Gとモデル化される。 振 動関經間 Gは、 謹母関数 gl、 g2によって次式 (4)のように表現された関数の空 間であり、 そのエネルギーは次式 (5)の Gノノレム// V // Gとして定義される。 ー(4)
\nG =mf l(s g2f\\L v = dxgl +dxg2]
gl.g2 L∞
- (5) J^iii像信号 I の分離問題は、 エネルギー汎関数を最 匕する次式 (6)の変分問題 として 化される。 この変分問題は、 Charabolleの Projection法によって角军く ことができる。
inf \J(u -^-\\i-u-vi2)
Figure imgf000011_0002
… ) これにより i^ f像信号 I力ら 離される第 2成分 Vはノィズの影響を受ける力 第 1成分 Uはノィズの影響をほとんど受けず、 ェッジを鈍らせることなく骨格成分 (謝可 画像構造) 力 S抽出される。
その他の加算型分 の例としては、線形フィルタによる ί職分離、 メジアン フィルタ、 モホロジフィルタ等による適応的平滑化を行うようにしてもよレ、。 以 下、 有界¾1)関数によらない加算型分離の ½feの例を示す。
例 1 :第 1成分 Uを顾像信号 Iに る直 底展開の低次の成分とし、 第 2成分 Vを願像信号 Iに ¾"Τる第 1成分 Uの纖とする 去
例 2:第 1成分 Uを顾像信号 Iにメジアンフィルタ麵した結果とし、 第 2 成分 Vを顾像信号 Iに ¾ ~る第 1成分 uの纖とする^ f去
例 3:第 1成分 Uを願像信号 Iに雜モホ口ジフィルタを掛けた結果とし、 第 2成分 Vを顾像信号 Iに ¾~Τる第 1成分 Uの藤とする 去
例 4: 第 1成分 Uを顾像信号 Iに縮小処理を施し、 さらに拡大処理を施し た結果とし、 第 2成分 Vを原画像信号 Iに る第 1成分 Uの歹錢とする ^去 例 5:第 1成分 Uを顾像信号 Iに Bilateral Filterを施した結果とし、 第 2成分 Vを顾像信号 Iに る第 1成分 Uの藤とする^ ¾
乗算型分離
続いて乗算型分離の雄にっレ、て説明する。 乗算型分離では、 顾像信号 Iが第 1成分 Uと第 2成分 Vの積によって表されるが、顾像信号 Iを対数変換して対数 原画像信号 f とすると、 7夂式 (7)のように加算型分離問題に変換することができ る。 I = U * V
f = u + v;
f = log/, u = logU, v = logV
…ひ) ここでは、 前述と同様に有界麵関数とノルムを用いた 去を説明する。 乗 算型 ^隹問題は対数領域において加算型分離と同様に A¾C変分モデルを用いて解 くことができる。 以下では、乗算型分離を «として、 対数領域における A¾C変 デルについて簡単に述べる。
式 (7)中の対数第 1成分 uの関数空間は、 前述の加算型分離モデルの第 1成分 U と同様に不 界によって区分された複数の "滑ら力な ¾ 変化の小部 貝 ' から構成された有界変動関数空間 BV (Bounded Variation Function Space)として モデル化され、 対 1成分 u のエネルギーは次式 (8)の TV (Total Variation)ノ ルム J(u)で定義される。 J(u) = J |Vw|^ dy
…⑧ 一方、 式ひ)中の対蘭 2成分 vの関輕間は、 mm Gとモデル化され る。 羅関 間 Gは、 繊母関数 gl、 g2によって式 (9)のように表現された関数 の空間であり、 そのエネルギーは式 (10)の Gノノレム// V //Gとして定義される。
V(x,y) = dx8l(x,y) + dxS2(x,y), gV 82 e LD
…⑨
IHIG = inf )2 + (S2 f Iし; = dxgl + dxg2 } · ' · (10) したがって、 対 像信号 fの分離問題は、 エネルギー汎関数を最 匕する次 式 (11)の変分問題として^;化される。
Figure imgf000014_0001
- (11) これにより顾像信号 Iカゝら分離される第 2成分 Vはノィズの影響を受けている 1 第 1成分 Uはノイズの影響をほとんど受けず、 エッジを鈍らせることなく骨格 成分 (謝可^画像構造) 力抽出される。
続いて、 ノイズ iJS»L パラメータ設¾3 3004で参照する信号レベル一ノイズ 分散モデルにっレ、て説明する。
顾像信号 Iに含まれるノイズをガウスノイズとすると、 ノイズの分散 σ (ノィ ズ S) は、 A/D変換直後の信号レベル Loに対して 2次曲線的に増加する。 これを特 開 2005-175718公報に開示されているように信号レベル一ノィズ モデルを 2次 関数で表すと、 次式 (12)力 S得られる。
σ = aL2 0 + pL。 + /
- - - (12) ここで、 ひ、 β、 γは 項である。 しかしながら、 ノイズの分散 σは信号レべ ノレだけではなく、 素子の ¾やゲインによっても変化する。 図 4は、 一例として、 ある温度 tにおけるゲインに関連する 3種頃の ISO感度 (ゲイン) 100、 200、 400 に るノイズの分散 σをプロットしている。 個々の曲線は式 (12)に示される形態 をしているが、 そ (^数はゲインに関連する ISO感度により異なる。 ?U¾を t、 ゲ インを gとし、 上記を考慮した形で信号レベル一ノイズ分散モデルの^;化を行う と、
〜(13) となる。 ここで、 ひ 、 y^ ^ ゲイン g に応じて決まる定数項であ る。 フル力ラ一の:^、 このノィズモデルは各舰立に適用することができる。 しかしながら、 第 1の¾»鄉では、 デジタノ M言号麵部 1005 の後段で成分分 離、 ノイズ iW理を行うので、 上記モデルをそのまま用いることはできなレヽ。 そ こで、 上記信号レベル一ノイズ分散モデノレの特性に加え、 デジタル信号処理部 1005 の特性を考慮し、 デジタル信号処理後の信号レベルについて信号レベル一ノ ィズ分散モデルを得る。
例えば、 デジタノ W言号処理部 1005力 12bitの入力信号を 8bit出力信号に! ¾ ¾ 変換する Knee処理又は Ί変換である場合、 デジタノ W言号処理部 1005は図 5の様な 入出力信 性を有している。 図中、 L (12)は A/D変換直後の信号レベル、 L(8)は デジタノ W言号処理後の信号レベルを表す。
したがって、 式 (12)あるいは (13)の A/D変換直後の信号レベル一ノィズ分散の特 性と、 図 5の様なデジタノ 言号処理部 1005 の特性を考慮すると、 デジタノ W言号処 理後の信号レベル Y (8)とノイズの分散 σは、 図 2に示した単峰の曲線の関係 (信 号レベル一ノイズ分散モデル) となる。
ノイズ iWQ;理パラメータ設 3004では、 図 2に示した信号レベル一ノイズ ^[モデルを参照してノイズの分散 σを取得し、 これに応じたノイズ 理パラ メータ Tl、 Τ2を設定している。 なお、 簡 匕のために、 図 2に代えて、 図 2を折 れ線で近似したデータを参照するようにしても.よレ、。
続いて、 ノイズ 理パラメータ Tl、 Τ2を用いた軟判定閾 理について説 明する。
願像信号 I力ら分離された第 2成分 Vは上記の通りノイズを含むので、 ノイズ i^ m 3002では、 以下に説明する軟判定閾ィ麵理 (コアリング処 により 第 2成分 Vに含まれるノィズを ί©^τる。
軟判定閾値処理は、 基準値を Β、 上限閾値を C1=B+T1 (Τ1>0)、 下限閾値を C2=B+T2 (T2く 0)とし、 鍵前の信 "^直 Aと処理後の信^ f直 Dの関係が次式 (14)で表 される信号処理である。
B + T1 ぐ (DM^- ΰ = A - Tl
B + T2 < A < B + T1の D = Β
Α ぐ B + Τ2の^^ D = Α - T2
ー(14) 図 6は軟判定閾値処理が行われる様子を示している。 図 6 (a)に示す 1 7火元の信 ^直 Aに対して軟判定閾ィ集理を行うと、 図 6 (b)に示す信 ^直 D力 S得られる。 こ のとき基 Bには信^直 Aの移動平均を用いている。
第 1の¾&¾¾態では、原画像信号 Iから分離された第 2成分 Vに対して上言 夂判 定閾 fi^理を行う。 第 2成分 Vは、 J^lli像信号 I力 骨格成分である第 1成分を除 いたものであるので、 直 Bは 0と考えて良レ、。
図 7は、 上言誠分分離、 ノイズ «2;理により麵像信号 Iのノイズが嫌され る 1»を示している。
図 7 (a)は、 17火元で表した ig®像信号 Iを示しており、 ノイズを含んでいる。 図 7 (b)は、 成分分 3001によって顾像信号 Iから分離された第 1成分 Uを 示している。 第 1成分 Uは平職分とエッジ成分を含 tf 格成分であり、 ノイズを 殆ど含んでいない。
図 7 (c) は、 成分分 3001によって得られた第 2成分 Vを示している。 第 2 成分 Vは顾像信号 Iに る第 1成分 Uの歹^)¾分であり、 S像信号 Iに含ま れてレヽたノイズと同等のノィズを含んでレ、る。
図 7 (d)は、 第 2成分 Vに対して軟判定閾 を行った結果である。 補正第 2 成分 V, は第 2成分 Vに比べてノイズが赚されている。 信号レベル一ノイズ獵 モデルを参照してノイズの分散 σを取得するためには、 通常、 ノイズの無い原信号 を推定する必要があるが、 第 1の¾»態では、 これに代えて骨格成分である第 1 成分 Uの信号レベルを用いてノイズの分散 σを取得するようにしている。 なお、 こ こでは、 簡単のために、 図 2に示した信号レベル一ノイズ分散モデルを折れ線で近 似したデータを参照してノィズの纖 σを取得した。
図 7 (e)は、 成 成処理部 3005で図 7 (b)の第 1成分 Uと図 7 (d)の補正第 2成 分 Vとを 1 : 1 で合成した結果である。 図 7の (a)と(e)を比較するとわかるよう に、合舰分 Γは顾像信号 Iに対してエッジ成分を保存しつつ、 ノイズが纖 されている。 なお、 合成時の混合比は任意に設定することができる。
続レ、て第 1の 態の作用効果にっレ、て説明する。
原画像信号 I に含まれる信号レベル依存のノイズを信号レベルから推定する場 合、 ノイズのなレ、顾像信号 Iを推定する必要があるが、 骨格成分である第 1成分 Uの信号レベルを用レ、れば 像信号 Iに含まれるノィズを精度良く推定すること ができる。 その一方で、 像信号 Iに る第 1成分 Uの 分である第 2成 分 Vには、 顾像信号 Iに含まれていたノイズと同等のノイズが分離される。 したがって、 第 1成分 Uの信号レベルに基づき、 第 2成分 Vのノイズを る ためのノイズ パラメータ Tl、 Τ2を設定し、 これを用いて第 2成分 Vのノ ィズを ί»Τるようにすれば、 第 2成分 Vのノイズを効果的に ί することがで き、 高画質なノイズ 理を徵することができる。 エッジ成分は骨格成分であ る第 1成分に含まれるので、 ノィズ 処理によってエツジが劣化することはな い。
なお、 ノイズ の; としては、 信号レベル依存の信号処理であれば何れ の;^去でもよく、 例えば、 上記の通り、 ノイズ ί©¾^理パラメータ Tl、 Τ2 に応じ た閾値と第 2成分 Vを比較して第 2成分 Vを補正する軟判定閾 を利用するこ とができる。
また、 成分分離の雜は、 上記した加算型分離、 乗算型分離の何れを用いても良 く、 また、 2成分への分離のみならず、 3成分以上に分 gt るようにしても構わな レ、。 また第 2〜第 N成分に "るノイズ «^理は、 上言嫩判定閾 ί集理に限定さ れず、 信号レべノレに基づきノイズを ί»Τる処理であれば、 f可れの; W去でもよレヽ。
<第 2の 態 >
第 2の実 »態以降の 態はフルカラーに対応した実沲形態である。 第 1の mmと共通の構成にっレヽては同じ参照符号を付して説明を省略する。 図 8は第 2の実施形態に係る撮像装置のシステム構成図であり、 固体撮像素子 1012 は色フィルタアレイを前面に配置したフルカラー用 素子である。 固^ t 像素子 1012は、 W¾、 複 のいずれでも良い。
^^分分 « 3008では第 1成分 Uを各 ^^分 (R成分、 G成分、 B成分) に分離 し、 各 分の信号レベルを取得する。
ノイズ iSMMノ ラメータ設 3004では、 各 ^^分の信号レベルに基づき、 図 2に示した信号レベル一ノィズ分散モデルあるいはその近似データを参照し、 ノ ィズの分散 σを J¾飾に取得する。 そして、 各械分のノイズ 理パラメ一 タ Tl、 Τ2 を、 対応する ^^分のノイズの分散 σに応じた値、 例えば、 ノイズの分 散 σに比例した値に設定する。
触分分 » 3009では、 成分分 «33001で得られた第 2成分 Vを各^^分に分 ¾H "る。
ノイズ «»!部 3002では、 第 2成分 Vの各械分に対して、 第 1の魏形態 で説明した軟判定閾 を ^に独立に行レ \ 各 ^分に含まれるノィズを ig«~る。
図 9は成分分离鄉3001力 成^成処理部 3005までの処理の を示したフロ 一チヤ一トである。 参考のため 3成分に分離する場合の処理を■で表している 力 ここでは 2成分に分离 f る齢にっレ、て説明する。
ステップ S 2 1では、 像信号 Iを第 1成分 Uと第 2成分 Vに分 る。 ステップ S 2 2では、 第 1成分 Uを各 ^^分に分離し、 各 ^^分の信号レベルを 得る。 ステップ S 2 3では、 各色成分の信号レべノレに基づき、 図 2に示した信号レベル —ノィズ分散モデ、ノレあるレヽはその近似データを参照してノィズの分散ひを^^ に取得する。 そして各触分のノィズの分散 σに応じて第 2成分 Vの各 ^^分に対 するノィズ« ^理パラメータ Tl、 Τ2を設定する。
ステップ S 2 4では、 第 2成分 Vを各 分に分离 tfる。
ステップ S 2 5では、 第 2成分 Vの各^^分に対して、 対応する 分のノイズ « ^理パラメータ Tl、 Τ2 を用いた軟判定閾 を^^^に独立に行い、 軟 判定閾 ί讓理後の値を補正第 2成分 Vとする。
ステップ S 2 6では、 第 1成分 U と補正第 2成分 Vとを合成し、 像信号 I に対してノイズが謹された合舰分 Γを得る。
この第 2の¾»維によれば、 顾像信号 I力 S複数の^^分を含む:^であって も、 第 1の¾»態と同様にエッジ成分の劣化を抑えつつ、 高画質なノイズ 理を行うことができる。
なお、 成分分離の雜〖 H算型分離、 乗算型分離の何れを用いても良い。 また、 第 2成分 Vのノイ
Figure imgf000020_0001
信号レ ベル依存の信号 であれば、 何れの;^去も適用できる。
<第 3の魏形態 >
図 1 0は第 3の^ β離に係る 置のシステム構成図である。 先の 態 と共通の構成にっレヽては同じ参照符号を付して説明を省略する。
第 3の雄形態では色成分分離部 3008に代えて Y/C分離部 3018を備える。 Y/C 分 « 3018では、 第 1成分 Uを艇成分 Yと色誠分 Cb、 Crに分離し、 成分 Yの信号レベルを取得する。 第 1成分 Uを構^ 1~る R成分、 G成分、 Β成分から輝 舰分 Υ、 色誠分 Cb、 Crへの変換は、 次式 (15)により行われる。
Y = 0. 29900R +0. 58700G +0. 11400B
Cb = -0. 16874R -0. 33126G +0. 50000B
Cr = 0. 50000R -0. 41869G ~0· 08131B
- (15)
Figure imgf000021_0001
3004では、 図 2に示した信号レベル一ノイズ 分散モデレレまたはその近似データを参照して) 分 Yの信号レべノレに応じた各色 成分のノイズの分散 σを取得する。 そして、 各 分のノイズSML理パラメータ Tl、 Τ2 を、 対応する色成分のノイズの分散 σに応じた値、 例えば、 ノイズの分散 σに比例した値に設定する。
^ ^^ ^ 3009では、成分分 » 3001で得られた第 2成分 Vを各械分に分 る。 ノイズ 部 3002は上言 EiC判定閾 理を^ ¾ ^に独立に行い、 各 ^^分に含まれるノィズを る。
図 1 1は成分分離都 3001から成分合 理部 3005までの処理を示したフローチ ヤートである。 参考のため 3成分に分 g る の処理を赚で表しているが、 こ こでは 2成分に分 る齢にっレ、て説明する。
ステップ S 3 1では、顾像信号 Iを第 1成分 Uと第 2成分 Vに分 g る。 ステップ S 3 2では、 第 1成分 Uを Y/C成分に分離し、輝 ^^分 Yの信号レベル を得る。
ステップ S 3 3では、 第 1成分 Uの鍵成分 Yの信号レベルに応じ、 図 2に示し た信号レベル一ノィズ分散モデノレまたはその近似データを参照してノィズの分散 σ を取得し、 これに応じた第 2成分 Vの各触分に财るノィズ 理パラメータ Tl、 Τ2をそれぞれ設定する。
ステップ S 3 4では、 第 2成分 Vを各色成分に分離する。
ステップ S 3 5では、 第 2成分 Vの各械分に対して、 ノイズィ 理パラメ一 タ Tl、 Τ2を用いた軟判定閾 を色成^に独立に行い、 補正第 2成分 V,を得 る。
ステップ S 3 6では、 第 1成分 U と補正第 2成分 V,とを合成し、 麵像信号 I に対してノイズが βされた合 β¾¾分 を得る。
この第 3の実¾ ^態によれば、 像信号 Iが複数の色成分を含む場合であって も、 第 1の »纖と同様にエッジ成分の劣化を抑えつつ、 高画質なノイズ 理を行うことができる。
なお、 成分分離の雜《¾Π算型分離、 乗算型分離の何れを用いても良い。 また、 第 2成分 Vのノイズをィ »Τる処理の 去は軟判定閾 ίίΜに限定されず、 信号レ ベル依存の信号処理であれば、 何れの; W去も適用できる。
また、 この雄形態においては、 第 1成分から離成分 Υを分離し、 鍵成分 Υ の信号レベルに応じてノイズ 理パラメータ Tl、 Τ2 を設定するようにしてい るが、 ^分 Υの信号レベルに代えて、 離成分 Υを最も反映している G成分を 補間して得られる信号の信号レベルを用いるようにしても構わなレ、。
<第 4の¾^態 >
図 1 2は第 4の纖形態に係る 置のシステム構成図である。 先の^^態 と共通の構成につレヽては同じ参照符号を付して説明を省略する。
第 1の¾»態では、 デジタノ M言号処理部 1005の後段に、 成分分 3001〜成. 分合成処理部 3005 を配置しているが、 第 4の実施形態では、 図 1 2に示すよう に、 デジタノ W言号処理部 1005の前段にこれらの要素を配置し、 A/D変換直後の信 号 (デジタル信号処理を行う前の信号) に対して成分分離、 ノイズ雌処理を行 う。 この^^、 成分分離を ft應として、 一旦フルカラ 言号 «像信号 I) を得 るために、 補間処理部 4000で、 色フィルタアレイを通して得られた信号をデモザ ィキングする。
そして、 第 4の¾»態では、 第 2の ¾5¾ ^態と同様に、 各触分の信号レベル に基づき、 信号レベル一ノイズ分散モデルを用いてノイズ低減処理パラメータ Tl、 Τ2 を色成 に設定する。 第 4の実施形態では、 デジタル信号処理の前段で 信"^離、 ノイズ β^Μを行うので、 信号レベル一ノイズ分散モデルとして、 図 2に示したモデルではなく、 図 1 3あるレヽは式 (12)、 (13)で表されるモデルを参照 してノィズの 散 σを取 ί辱する。
その後の、 軟判定閾 ίί^Μによる第 2成分 Vのノイズ «Μ、 第 1成分 Uと補 正第 2成分 V,の合)^理に関しては、 第 2の 態と同様である。
この第 4の実施形態によれば、 原画像信号 Iが複数の色成分を含む場合に適応 し、 第 1の!^態と同様にエッジ成分の劣化を抑えつつ、 高画質なノイズ®» 理を行うことができる。 また、 デジタル画脑理部の前段において成分分離、 ノィ ズ« ^理を行うので、 成分分 » 3001 における骨格成分とノイズ成分の分 力が向上する。 なお、成分分離の雜〖 H算型分離、 乗算型分離の何れを用いても良い。 また、 第 2成分 Vのノイズを is t る «の は軟判定閾ィ幽理に限定されず、信号レ ベル依存の信号処理であれば、何れの 去も適用できる。
また、 この第 4の¾»態では各^ ^分の信号レベルに基づきノイズの分散ひを 取得し、 これに応じてノイズ パラメータ Tl、 Τ2を設定するようにしてい るが、 各触分の信号レベルに代えて、 A/D変換直後の全械分を変換して得られ る) Si 成分 Υの信号レベル、 あるいは、 成分 Yを最も反映している G成分を補 間して得られる信号の信号レベルを用いるようにしてもよい。
<第 5の 態 >
ノイズを含まないカラー画像において、 3 X 3の近イ^域で言«した局所色間相 互相隱、数 (以下、 「相賺、数」 という。 ) は 1に近い正値を示す一方で、 ノイズ を含むカラー画像では、 相 Γ蕭数が低下する画素を含むため、 近 ί練域で籠した 相関係数は低下する。 このこと力 、相関係数の値が小さいほど画像信号に含まれ るノイズが多いといえる。
そこで、 第 5の魏形態では、 触分間の相関に応じて軟判定閾讓理のノイズ パラメータ Tl、 Τ2を補正する。 これにより、 相関が高レ、領域では軟判定 閾ィ維理での幅が狭まり、願像信号 Iを ί藉する様に作用する一方で、 相関の低 レ、領域では軟判定閾 ίί 理での幅力 S広がって、 ]¾1ί像信号 Iを平坦化することがで き、 より高精度のノイズ低 理が できる。
図 1 4は第 5の¾»態に係る it ^置のシステム構成図である。 先の 態 と共通の構成にっレヽては同じ参照符号を付して説明を省略する。 デジタノ W言号処理部 1005 カ ら出力された顾像信号 I (フルカラ^ ί言号) は色 成分分離部 3022、 相関演算処理部 3006及びノイズ ifi ^理部 3002 に 5¾され る。
分^ 8鄉 3022は顾像信号 Iを各触分に分離し、 代 ¾ί言号レベル算出部 3023では、 分離された各色成分にっレ、て所定の画素位置近傍の代表信号レベルを 算出する。 代 ¾ί言号レベルを算出するにあたっては、 線形フィルタ、 メジアンフィ ルタ等の非線形フィルタ、 バイラテラルフィルタ等の非線? フィルタ、 モホロ ジフィルタ等を體することができる。
ノイズ 理パラメータ設 3004では、 図 2に示した信号レベル一ノイズ モデルあるいはその近似データを参照して、 各色成分の代謝言号レベルに応じ たノイズの分散 σを取得する。 そして、 各色成分のノイズ低減処理パラメータ Tl、 Τ2 を、 対応する色成分のノイズの分散 σに応じた値、 例えば、 ノイズの分散 σに比例した値に設定する。
相関演算処理部 3006では、 後财る様に、 注目画素位置を含 應での各 触分の相赚数を演算し、 その最小値である最 湘隱数 rを出力する。
ノ ラメータ補正部 3007では、 最 /J4目難数 rに基づき補正係数 Cを算出し、 こ れを各触分のノイズ ίδ»ί¾パラメータ Tl、 Τ2 に掛けることで各^^分のノィ ズ 理パラメータ Tl、 Τ2を補正する。
ノイズ iS^L理部 3002では、 補正後のノイズ ®»理パラメータ Tl、 Τ2を用い た軟判定閾値 を色成分毎に独立に行い、 原画像信号 Iに含まれる信号レベル依 存のノイズを る。 相関演算処理とノイズ 理パラメータ Tl、 Τ2の補正
触分間の相関に応じたノイズ 理パラメータ Τ1、 Τ2の補正の内容を、 具 体例を挙げて説明する。 ここでは、説明を簡単にするために^ ®像信号 Iが 1次元 であるとして説明する。
図 1 5、 図 1 6は、願像信号 Iの画素位置と各触分の信号レベルとの関係を 示している。 図 1 5に示す例では、 触分間に相関がなく、 この 、 各画素位置 の各触分の信号レベルは、
画素位置 R G _Β_
1 150 180 130
2 220 120 140
3 160 150 100
4 200 140 160
5 140 100 200
である。 このとき、 R~G間、 G~B間、 B - R間の各相,数を求めると、
R-G GzB B-R
相関係数 -0. 15836 -0. 69668 -0. 14106
となり、 最 /J4目関係数 rは、 -0. 69668となる。
一方、 図 1 6に示す例では、 色成分間に相関があり、 各画素位置の各色成分の信 号レベルは、
画素位置 _R_ _G_ _B_
1 100 120 150 2 150 160 180
3 170 160 180
4 100 120 90
5 90 100 100
である。 このとき、 R~G間、 0"B間、 B-R間の各相 Π、数を求めると、
R-G G-B B-R
相^^数 0. 962099 0. 866459 0. 848117
となり、 最 /J4目聽数 rは、 0. 848117となる。
第 5の 態では、 最 隱数 rを用いてノイズ 理パラメータ Tl、 Τ2をそれぞれ、 次式 (16)により補正する。
Tl<-c (l-r)2Tl
T2^c (l-r)2T2
- (16) c=l とすると、 ノイズ ^理パラメータ Tl、 Τ2の補正係数 Cは (1- r)2とな り、 補正係数 Cは、 図 1 5の 、 2. 878739で、 図 1 6の 、 0. 023069 とな る。
したがって、 図 1 5の様に触分間に相関がない:^は、軟判定閾値の上限閾値 と下限閾値との幅が大きくなる様に補正され、 その結果、 原信号をより 骨化する ように作用する。 一方で、 図 1 6の様に 分間に相関がある^は、 軟判定閾値 の上限閾値と下限閾値との幅が小さくなる様に補正され、 より信号を保存するよう な僅かな補正力 S行われる。 このように、 第 5の難形態によれば、 触分間の相関に応じて軟判定閾 ίίΜ のノイズ 理パラメータ Tl、 Τ2 を補正するようにしたことで、 像信号 I に含まれる »)成分のうち、 ノイズに起因する成分と、 テクスチャのような画像本 来の構造に起因する成分を判別することができ、 ノィズ によつて後者の成 分が劣ィ匕するのを抑えることができる。
なお、 ここでは代謝言号レベルからノィズの σを求め、 これに応じたノイズ 理パラメータ Tl、 Τ2 を設定するようにしているが、 代 言号レベルに代え て、 全触分を変換して得られる) TO成分 Υの信号レベル、 あるいは、 離成分 Y を最も反映してレ、る G成分を補間して得られる信号の信号レベルを用レ、てもよレ、。
<第 6の ^態〉
図 1 7は第 6の ^態に係る 置のシステム構成図である。 先の ^態 と共通の構成にっレ、ては同じ参照符号を付して説明を省略する。
6の¾ ^態では、 第4の¾»維と同様に、 A/D変換直後の信号 (デジタル 信号処理を行う前の信号) に対して軟判定閾値処理を行う。 したがって、 第 6の実 層維では、 一旦フルカラ ^{言号 麵像信号 I) を得るために、 補間処理部 4000 で色フィルタァレイを通して得られた信号をデモザィキングする。
また、 デジタノ W言号処理を施す前に成分分離、 ノイズ 理を行うので、 ノィ ズ^^理パラメータ設 3004では、 各触分の信号レベル一ノイズの σ のモデノレとして、 図 2に示したモデルではなく、 図 1 3あるいは式 (12)、 (13)で表 されるモデルを参照して各^^分のノイズの分散 σを取得し、 これに応じたノイズ «処理パラメータ Tl、 Τ2 を色成分毎に設定する。 ノイズ «処理パラメータ Tl、 Τ2は、 例えば、 ノイズの分散 σに比例する値に設定する。
^^分間の相関による軟判定閾 tt^理のノイズ低»¾パラメータ Tl、 Τ2の補 正は、 ノ、。ラメータ補正部 3007において第 5の¾»態と同様に式 (16)により行わ れる、 また、 スィズ纖処理部 3002 では、 補正後のノイズ«処理パラメータ Tl、 Τ2を用いた軟判定閾値処理を原画像信号 Iに対して行い、 原画像信号 Iに含 まれる信号レベル依存のノィズを iS^Tる。
色フィルタアレイ信号抽出部 5000では、 ノイズ、«^理後の画像信号から色フ ィルタァレイの配列に対応する位置の信号を取り出し、 デジタル信号処理部 1005 に画像信号を する。 デジタル信号処理後の信号は後段の表示素子 2001、 記憶 メディァ 2003へと される。
第 6の ^態では、 第 5の¾»態と同様に械分間の相関に応じて軟判定閾 ffi 理のノイズ 理パラメータ Τ1、 Τ2 を補正するようにしたことで、 i^ 像 信号 Iに含まれる ¾1¾成分のうち、 ノイズに起因する成分と、 テクスチャのような 画像本来の構造に起因する成分を判別することができ、 ノィズ により後者 の成分が劣化するのを抑えることができる。
なお、 ここでは代雄号レベルからノイズの分散 σを求め、 これに応じたノイズ « ^理パラメータ Π、 Τ2 を設定するようにしている力 代 言号レベルに代え て、 A/D変換直後の全 ^^分を変換して得られる) ¾ 成分 Υの信号レベル、 あるい は、 贿成分 Yを最も反映している G成分を補間して得られる信号の信号レベルを 用いてもよレヽ。
<第 7の ¾¾^態> 図 1 8は第 7の¾ ^態に係る 置のシステム構成図である。 先の 態 と共通の構成については同じ参照符号を付して説明を省略する。
第 7の¾5^態では、 第 2の ^態に示した成分分 ¾ ^の軟判定閾ィ幽理と、 第 5、 第 6の ^態に示した 分間の相関に応じたノィズS»理パラメータ Tl、 Τ2 の補正を ¾ίϋ~る。 ただし、 第 5、 第 6の実施形態と異なり、 触分間の 相関演算は顾像信号 I力ゝら僵された第 2成分 Vに対して行う。
顾像信号 I力 分離された第 1成分 Uは触分分讓 3008に送られ、 械分 ^133008では各 ^^分の信号レベルを取得する。
ノイズ ffi^®パラメータ取纏 3004では、 第 2の¾»態と同様に、 各械 分の信号レベルに基づき、 信号レベル一ノイズ分散のモデルを参照してノイズの分 散 σを触分毎に取得する。 そして、 各色成分のノイズ ifi^t理パラメータ Tl、 Τ2 を、 対応する色成分のノイズの分散 σに応じた値、 例えば、 ノイズの分散ひに 比例した値に設定する。
J^ 像信号 I力、ら分離された第 2成分 Vは、 ^^分, 3009で各 分に分 離された後、 相関演算処理部 3006に{ され、 相関演算 部 3006では、 第 2成 分 Vに対して第 5、 第 6の ¾1 ^態と同様に各色信号の相関演算処理を行い、 最小 相関係数 を算出する。
ノ ラメータ補正部 3007では、 最 /J4目 M ¾ rに基づきノイズ (g»理パラメ一 タ Tl、 Τ2 の補正係数 C を演算し、 これを各色成分のノイズ鎌処理パラメータ Tl、 Τ2に掛けて各^ ^分のノイズ ifi^理パラメータ Tl、 Τ2を補正する。
ノイズ i®»理部 3002では、 第 2成分 Vの各触分に対して、 補正したノイズ 理パラメータ Tl、 Τ2を用いた軟判定閾 ffi^ を触^に独立に行い、各 色成分に含まれるノイズを麵する。 そして、 得られた信号を補正第 2成分 V,と して出力する。
成^^ 理部 3005では、 第 1成分 ϋと補正第 2成分 Vを所定の比率、例えば 1 : 1の比率で合成し、顾像信号 Iに対してノイズが纖された合戯分 Γを 得る。
相関演算処理とノィズfi ^理パラメータ Tl、 Τ2の補正
械分間の相関に応じたノイズ パラメータ Tl、 Τ2の補正の内容を、 具 体例を挙げて説明する。 ここでは、 説明を簡単にするために第 2成分 Vが 1次元で あるとして説明する。 第 2成分 Vがテクスチャの様な細力ゝぃ構 ¾ ^分とノイズを含 む、顾像信号 Iに ¾~Τる第 1成分 Uの 分であるので、 第 2成分 Vはゼロを 中心として変動する成分である。
図 1 9、 図 2 0は第 2成分 Vの画素位置と各 ^^分の信号レベルとの関係を示し ている。
図 1 9に示す例では、 色成分間に相関がなく、 各画素位置の各色成分の信号レべ ルは、
画素位置 _R_ _G_ _B
1 -5 7 2
2 4 -2 2
3 -2 4 4
4 1 -5 5 5 6 -1 2
である。 R~G間、 G-B間、 B- R間の各相 Γ蕭、数を求めると、
Figure imgf000032_0001
相 Μ 、数 -0. 75148 -0. 36622 -0. 19914
となり、 最 /J4目斷、数 rは、 -0. 75148となる。
一方、 図 2 0の例では、 触分間に相関があり、 各画素位置の各触分の信号レ ベルは、
画素位置 R G B
1 1 0 -1
2 4 3 3
3 5 3 4
4 -2 -1 0
5 -1 0 -1
である。 R~G間、 G~B間、 B-R間の各相関係数を求めると、
R-G G=B BzR
相関係数 0. 964023 0. 911685 0. 873891
となり、 最 /J4目隱数 rは、 0. 873891となる。
以上のようにして最小相関係数 r を算出したら、 ノイズ低減処理パラメータ Tl、 Τ2は、 第 5の^^態で用いた式 (16)によりそれぞ; ^t正される。 c=lとする と、 ノィズ« ^理パラメータ Tl、 T2の補正係数 Cは (ト r) 2となり、補正係数 C は、 図 1 9の:^で 3. 0677、 図 2 0の:^で 0. 01590となる。 したがって、 図 1 9の様に械分間に相関がなレ、 は、 軟判定閾値の上限閾値 と下限閾値との幅が大きくなる様に補正され、 その結果、 信号をより 骨化するよ うに作用する。 一方で、 図 2 0の様に触分間に相関がある:^は、 軟判定閾値の 上限閾値と下限閾値との幅が小さくなる様に補正され、 より信号を保存するような 僅かな補正が行われる。
図 2 1は、 成分分離部 3001から成 合 部 3005までの処理に対応するフ口 一チャートである。 参考のため J^®像信号 Iを 3成分に分離する の処理を石 で表してレ、るが、 ここでは 2成分に »Τる につレヽて説明する。
ステップ S 7 1では、 顾像信号 Iを第 1成分 Uと第 2成分 Vに分 ¾ る。 ステップ S 7 2では、 第 1成分 Uを各色成分に分離し、 各色成分の信号レベルを 得る。
ステップ S 7 3では、 各色成分の信号レべノレに基づき、 図 2に示した信号レベル —ノイズ分散モデルまたはその近似データを参照してノイズの分散 σを取得し、 こ れに応じて、 にノイズ 理パラメータ Tl、 Τ2を設定する。
ステップ S 7 4では、 第 2成分 Vを各 ^分に分 る。
ステップ S 7 5では、 第 2成分 Vの各触分の相関を演算し、 最 /』4目聽数 rを 算出する。
ステップ S 7 6では、 最 /J湘赚数 rに基づき、 ステップ S 7 3で得られた各色 信号のノイズ « ^理パラメータ Tl、 Τ2の補正係数 Cを算出する。
ステップ S 7 7では、 ノイズ 理パラメータ Tl、 Τ2に補正係数 Cを掛けて ノイズ 理パラメータ Tl、 Τ2を補正する。 そして、 第 2成分 Vの各械分に 対して補正後のノイズ «処理パラメータ Tl、 Τ2 を用いた軟判定閾 ίί 理を行 い、補正第 2成分 V,を得る。
ステップ S 7 8では、 第 1成分 Uと補正第 2成分 Vを合成して、 J^®像信号 I に対してノイズが赚された合戯分 Γを得る。
第 7の ¾ϋ ^態によれば、願像信号 I力 S複数の触分を含む であっても、 第 1の 態と同様にエッジ成分の劣化を抑えつつ、 高画質なノイズ ίβ»Μを 行うことができる。
また、 色成分間の相関に応じて軟判定閾値処理のノィズ低減処理パラメータ Tl、 Τ2を補正するようにしたことにより、 第 2成分 Vに含まれる β成分のうち ノイズに起因する成分と、 テクスチャのような画像本来の構造に起因する成分を判 別することができ、 ノィズ によって後者の成分が劣化するのを抑えること ができる。
なお、成分分離の雜〖 H算型分離、乗算型分離の何れを用いても良い。 また、 第 2成分 Vのノイズを る処理の は軟判定閾 ίίΜに限定されず、信号レ ベル依存の信号処理であれば、 何れの; W去も適用できる。
なお、 ここでは第 1成分 Uの各色信号の信号レベルに基づきノィズの分散 σを取 得し、 これに応じたノイズ 理パラメータ Tl、 Τ2を設定するようにしている 力 第 1成分 ϋの各色信号の信号レベルに代えて、 第 1成分から分離し d¾¾J¾分 Yの信号レべノレ、 あるいは、輝 «分 Yを最も反映している G成分を捕間して得ら れる信号の信号レベルを用いてもよい。
<第 8の難形態 > 図 2 2は第 8の^^態に係る »g置のシステム構成図である。 先の 態 と共通の構成にっレ、ては同じ参照符号を付して説明を省略する。
第 8の¾»態では、 第 4、 第 6の ^態と同様に、 A/D変換直後の信号 (デ ジタノ W言号麵を行う前の信号) に対して成分分離、 ノイズ 理を行う。 この 、 成分分離を前提として、 一旦フルカラー信号 (原画像信号 I) を得るため に、 補間処理部 4000で、 色フィルタアレイを通して得られた信号をデモザィキン グする。
また、 デジタノ W言号処理を施す前に成分分離、 ノイズ « ^理を行うので、 ノィ ズ パラメータ設^ 13 3004は、 各触分の信号レベル一ノイズ分散のモデ ノレとして、 図 2に示したモデルではなく、 図 1 3あるいは式 (12)、 (13)で表される モデルを用いてノイズの纖ひを^ ¾ ^に取得する。 そして、 各械分のノイズ
« ^理パラメータ Tl、 Τ2 を、 対応する械分のノイズの分散 σに応じた値、 例 えば、 ノィズの分散 σに比例した値に設定する。
その後の軟判定閾 «L理によるノィズ 理は他の 態と同様である。 ま た、 ^^分間の相関によるノイズ パラメータ Tl、 Τ2 の補正は第 5の 形態と同様であり、 最 /J湘関係数 rに応じて式 (16)によりノィズ« ^理パラメ一 タ Tl、 Τ2力 S捕正される。
ノイズ 理部 3002力、ら出力される補正第 2成分 V,は、 成分合成処理部 3005 で第 1成分 Uと合成される。 色フィルタアレイ信号抽出部 5000では、 合舰分 Γ から色フィルタアレイの配列に対応する位置の信号を取り出し、 デジタノ W言号摩 部藤に画像信号を磁する。 この第 8の^ ^態によれば、原画像信号 I力 S複数の械分を含む ^^であって も、 第 1の¾»態と同様にエッジ成分の劣化を抑えつつ、 高画質なノイズ ifi^ 理を行うことができる。
また、 色成分間の相関に応じて軟判定閾値処理のノィズ低減処理パラメータ Tl、 Τ2 を補正するようにしたことにより、 第 2成分 Vに含まれる変動成分のう ち、 ノイズに起因する成分と、 テクスチャのような画像本来の構造に起因する成分 を判別することができ、 ノィズ によつて後者の成分が劣、ィ匕するのを抑える ことができる。
また、 デジタノ 言号処理部 1005 の前段において成分分離、 ノイズ 理を行 うので、 成分分 3001における骨格成分とノイズ成分の分 ϋΙ^Ιが向上する。 なお、 成分彌の雜 ¾¾u算型分離、乗算型分離の何れを用いても良い。 また、 第 2成分 Vのノイズを ifi^r る処理の 去は軟判定閾 ίί ®に限定されず、信号レ ベル依存の信号麵であれば、 何れの;^去も適用できる。
また、 ここでは第 1成分 Uの各色信号の信号レべノレに基づきノイズの分散 σを取 得し、 これに応じたノイズ パラメータ Tl、 Τ2を設定するようにしている 力 第 1成分 ϋの各色信号の信号レベルに代えて、 A/D変換直後の全 分を変換 して得られる) ¾¾¾分 Υの信号レベル、 あるいは、 ^分 Υを最も反映している G成分を補間して得られる信号の信号レベルを用レ、てもよい。
以上、 本発明の実»態について説明したが、 上記 «?1 ^態は本発明の適用例の 一部を示したに過ぎず、 本発明の技«範囲を上記実»態の具体^成に限定す る趣旨ではなレ、。

Claims

請 求 の 範 囲
1. 画像処 置において、
顾像信号を、 骨格成分である第 1成分と、
Figure imgf000037_0001
1成 分の藤から得られる第 2成分と、 を含む複数の成分に分 ¾·τる成分分離手段と、 Ι ΙΞ第 1成分の信号レベルを取得する信号レベル取 段と、
編己第 1成分の信号レベルに基づきノィズ¾»¾パラメータを設定するパラメ ータ設定手段と、
廳己ノィズ« ^理パラメータを用レ、て ΐίιΐΕ第 2成分のノィズを ig^ るノイズ ί 手段と、
を備えたことを mとする画像^^ ¾
2. 請求項 1に言2¾の画像処«置において、
膽己ノィズ纖手段は、 謝己ノィズ パラメータに応じて設定される閾値 と膽己第 2成分との比^果に基づき iif己第 2成分のノィズを »Τることを if« とする画^ 0«
3. 請求項 1または 2に f¾feの画 置において、
Figure imgf000037_0002
クロであり、
嫌己第 1成分及ひ rt己第 2成分はそれぞれ) 分であることを赚とする画像 処職
4. 請求項 1または 2に記載の画像 ¾置において、
tirts^ i像信号は複数の ^^分力ゝら構成され、
膽己信号レベル取 段は、 廳己第 1成分の各触分の信号レベルを取得し、 ΙίίΙΕ^ラメータ設定手段は、 ttJlH第 1成分の各触分の信号レベルに基づき、 前 記第 2成分の各 分に対応するノィズ 理パラメータを設定し、
ΙΕノイズ ί繊手段は、 辦己第 2成分の各触分のノイズを、 対応する触分の ノィズ ί®»理パラメータを用レ、て »Τる、
ことを とする画像処 ¾
5. 請求項 1または 2に言 Ξ«の画^ «置において、
ΙΙίΙΕ^ ί像信号は複数の^^分力ゝら構成され、
嫌己信号レベル取 段は、 嫌己第 1成分の鍵成分の信号レベルあるいは G成 分を補間して得られる信号の信号レベルを取得し、
爾 Ε/、°ラメータ設定手段は、 魔己第 1成分の鍵成分の信号レベルあるいは G成 分を補間して得られる信号の信号レベルに基づき、 嫌己第 2成分の各 ^^分に対応 するノィズ パラメータを設定し、
嫌己ノイズ «手段は、 ΙΕ第 2成分の各^^分のノイズを、 対応する^ ^分の ノィズ 理パラメータを用レ、て る、
ことを とする画
6. 請求項 1または 2に言 S«の画 置にぉレヽて、
ttjiai ii像信号は複数の触分から構成され、
嫌己画腕職置は、
ΙίίΙΒ第 2成分に対して、 注目画素を含 域での触分間の相難数を演算 する相関演算 手段と、
tiri¾¾赚数に基づき膽己ノィズ パラメータを補正するパラメータ補正 手段と、
を備えたことを纖とする画 7. 請求項 6に言 2«の画 置において、
Figure imgf000039_0001
辦己注目画素を含 i錄域での触分間の相隱 数に対して負の相関を持つ補正係数を算出し、 ΙίίΙΞノィズ ifi^Mパラメータに前 言 ¾i正係数を掛けて嫌己ノィズ te^ パラメータを補正することを とする画 舰職齓
8. 請求項 7に ts¾の画 置において、
膽 、°ラメータ補正手段は嫌 Bffi正係数を次式:
C = (1-r)2
r:注目画素を含 tifii 域での^^分間の相関係数
により算出することを [とする画^
9. 請求項 6から 8のレ、ずれか一つに言 2<feの画»^置において、
嫌 21^®像信号は 3種以上の械分を含み、
1«¾関演算処理手段は、 IflfEtS隱数として複数の舗分間の相隱数のうち 最 /jィ直を職することを難とする画
1 0 · 請求項 6から 9のレヽずれか一つ言 Ξ«の画像処¾置において、
ΙΐίΙΕ信号レベル取 段は、 悲第 1成分の各触分の信号レベルを取得し、 嫌 S/ ラメータ設定手段は、 廳己第 1成分の各械分の信号レベルに基づき、 前 記第 2成分の各械分に対応するノィズ パラメータを設定し、
tiff己ノイズ 手段は、 l己第 2成分の各 分のノイズを、 対応する ^^分の ノィズWOf理パラメータを用レ、て « る、
ことを とする画^^
1 1. 請求項 6から 9のレヽずれか一つに言 の画 置において、
編己信号レベル取 ^段は、 廳己第 1成分の) ¾¾¾分の信号レベルあるいは G成 分を補間して得られる信号の信号レベルを取得し、
嫌 S/、°ラメータ設定手段は、 嫌己第 1成分の鍵成分の信号レベルあるいは G成 分を補間して得られる信号の信号レベルに基づき、 膽己第 2成分の各 ^^分に対応 するノィズ 理パラメータを設定し、
嫌己ノイズ 手段は、 ΙίϊΙΒ第 2成分の各 分のノイズを、 対応する 分の ノィズ iSl^理パラメータを用レヽて ifi^Tる、 ことを «とする画像^ Mft
1 2. 請求項 1力、ら 1 1のレ、ずれか一つに識の画難難置にぉレヽて、
嫌己第 1成分と、 ノィズ ί繊手段によってノィズが鎌された 第 2成分 との合成を行う合成手段を備えたことを赚とする画像処¾¾
1 3. 請求項 1から 1 2のレヽずれか一つに |5«の画 置において、
嫌己第 1成分は、顾像の平職分を含む成分であることを赚とする画 装
1 4. 請求項 1力、ら 1 3のレ、ずれか一つに の画 »職置にぉレ、て、
膽己第 1成分は、願像のエッジ成分を含むことを赚とする u^^am
1 5. 請求項 1力、ら 1 4のレヽずれか一つに言 B¾の画^ «置におレ、て、
ta ^分分離手段 卩算型分離により ia^®像信号を tin己第 1成分と MI己第 2 成分を含む複数の成分に分 ることを樹数とする » «¾5
1 6. 請求項 1から 1 4のレヽずれか一つに の画 置において、
Ιΐ«分分離手段は乗算型分離により嫌 画像信号を肅己第 1成分と l己第 2 成分を含む複数の成分に分 g rることを とする MfLeSo
1 7. 請求項 1から 1 6のレヽずれか一つに言 の画^; ¾置において、
,変換を行うデジタノ 言号^ M手段を備え、
tiHB^®像信号は tiff己デジタル ί言号 通手段による ΜΙΕ,変換前の信号である ことを とする画舰«
1 8. 請求項 1 2に言 Sttの画 置におレヽて、
tins合成手段によつて得られた信号のうち、 所定の色フィルタの配置に対応した 色信号を取り出す信号 «手段を備えたことを とする画^^ ¾ 1 9. 置において、
纖素子と、
ffifS搬素子から得られた顾像信号を、 骨格成分である第 1成分と、 藤 fl^® 像信号に ¾~ る謝己第 1成分の から得られる第 2成分と、 を含む複数の成分に る成 離手段と、
tins第 1成分の信号レベルを取得する信号レベル取辭段と、
嫌己第 1成分の信号レベルに基づきノィズ ig»理パラメータを設定するパラメ ータ設定手段と、
tilfSノィズ« ^理パラメータを用レヽて ilfE第 2成分のノィズを るノイズ 赚手段と、
を備えたことを»とする
2 0. 請求項 1 9に fSttの 置において、
tiifBノィズ ί職手段は、 歸己ノィズ 理パラメータに応じて設定される閾値 と tiff己第 2成分との比維果に基づき tfffS第 2成分のノィズを «1 "ることを赚 とする
2 1. 請求項 1 9または 20に! Ξ¾の 置にぉレヽて、
te^ i像信号はモノクロであり、
tilt己第 1成分及ひ曹己第 2成分はそれぞれ 分であることを とする謙 装 o
2 2. 請求項 1 9または 2 0に の 置にぉレヽて、
嫌己謹素子は色フィルタァレイを前面に配置した謙素子であり、
ttna^ i像信号は複数の^^分力ら構成され、
編己信号レベル取 段は、 編己第 1成分の各触分の信号レベルを取得し、 ラメータ設定手段は、 嫌己第 1成分の各触分の信号レベルに基づき、 前 記第 2成分の各色成分に対応するノィズ パラメータを設定し、
i己ノイズ «手段は、 I己第 2成分の各 ^^分のノイズを、 対応する械分の ノィズ 理パラメータを用レ、て » る、
ことを mとする
2 3 · 請求項 1 9または 20に言5¾の衝蝶置にお!/ヽて、 嫌己鐘素子は色フィルタァレイを前面に配置した撤素子であり、
ΙϋΙΕ^ ί像信号は複数の触分から構成され、
嫌己信号レベル取 段は、 flit己第 1成分の 分の信号レベルあるいは G成 分を補間して得られる信号の信号レベルを取得し、
嫌 5/ ラメータ設定手段は、 嫌己第 1成分の) ¾¾¾分の信号レベルあるいは G成 分を補間して得られる信号の信号レベルに基づき、 tin己第 2成分の各 ^^分に対応 するノィズ 理パラメータを設定し、
嫌己ノイズ繊手段は、 嫌己第 2成分の各 ^^分のノイズを、 対応する^^分の ノィズ ί®|^理パラメータを用レ、て ί»Τる、
ことを とする
2 4. 請求項 1 9または 2 0に |S¾の 置にぉレヽて、
嫌己 »素子は色フィルタ了レイを前面に配置した «素子であり、
嫌己 置は、
前記第 2成分に対して、 注目画素を含む近傍領域での色成分間の相関係数を演算 する相関演算処理手段と、
ttflB目難数に基づき tillSノィズ パラメータを補正するパラメータ ffi!E 手段と、
を備えたことを とする
2 5 · 請求項 2 4に言識の 置にぉレヽて、 ΙΐΠΞ/、°ラメータ補正手段は、 IflfE注目画素を含 tHfi^g域での 分間の相 Μ 、 数に対して負の相関を持つ補正係数を算出し、 tinaノィズ パラメ→に前 言 正係数を掛けて膽己ノィズ パラメータを補正することを樹敫とする撮
2 6. 請求項 2 5に f識の¾¾置にぉレヽて、
IB^ラメータ補正手段は Ι Μ正係数を次式:
C = (1-r)2
r:注目画素を含 tiifi^域での 分間の相関係数
により算出することを mとする
2 7. 請求項 2 4から 2 6のレヽずれか一つに «の 置におレ、て、
MtEi^®像信号は 3種以上の 分を含み、
tinsta関演算 ^手段は、 ΜΙΕ¾,数として複数の "^分間の相 π、数のうち 最 /Jィ直を職することを とする
2 8. 請求項 2 4から 2 7のいずれか一つ ISttの 置において、
tfiiB信号レベル取 段は、 tins第 1成分の各 ^^分の信号レベルを取得し、 嫌 ラメータ設定手段は、 tirta第 1成分の各械分の信号レベルに基づき嫌己 第 2成分の各触分に対応するノィズ 理パラメータを設定し、
嫌己ノイズ «手段は、 tiff己第 2成分の各 ^分のノイズを、 対応する ^^分の ノィズ 理パラメータを用レ、て ig» る、
ことを ί数とする
2 9. 請求項 2 4から 2 7のレ、ずれ力一つに言 Ε«の 置において、
嫌己信号レベル取 段は、 ΙίίΙΒ第 1成分の鍵成分の信号レベルあるいは G成 分を補間して得られる信号の信号レベルを取得し、
!S/、°ラメータ設定手段は、 編己第 1成分の離成分の信号レベルあるいは G成 分を補間して得られる信号の信号レベルに基づき、 fflfE第 2成分の各 ^^分に対応 するノィズ 理パラメータを設定し、
IfJlBノイズ «手段は、 Ml己第 2成分の各 ^^分のノイズを、 対応する^ ^分の ノィズ 理パラメータを用レ、て る、
ことを樹敫とする
3 0. 請求項 1 9力 ら 2 9のレ、ずれか一つに言 Bttの 置におレ、て、
l己第 1成分と、 ΙΐίΙΒノイズ 手段によってノイズが «された tiifS第 2成分 との合成を行う合成手段を備えたことを難とする .
3 1. 請求項 1 9から 3 0のレヽずれか一つに言 2¾の»¾置において、
tin己第 1成分は、 原画像の平 a ^分を含む成分であることを とする撮像装 ft
3 2. 請求項 1 9から 3 1のレヽずれか一つに言 Sttの 置において、
ΙίίΙΞ第 1成分は、顾像のエッジ成分を含むことを赚とする ¾|^¾齓
3 3. 請求項 1 9から 3 2のレヽずれか一つに の 置において、
肅诚分分離手段 ロ算型分離により嫌 a^®像信号を編己第 1成分と嫌己第 2 成分を含む複数の成分に分 81 " る
Figure imgf000047_0001
3 4. 請求項 1 9から 3 2のレヽずれか一つに言 の 置にぉレ、て、
tf«分分離手段は乗算型分離により嫌 Si ®像信号を flit己第 1成分と嫌己第 2 成分を含む複数の成分に分 i ることを «とする
3 5. 請求項 1 9力ら 3 4のレヽずれか一つに言 の 置にぉレ、て、
翻変換を行うデジタノ W言号処理手段を備え、
tina^®像信号は tin己デジタノ w言号処理手段による isra変換前の信号である ことを とする
3 6. 請求項 3 0に言識の 置にぉレヽて、
膽己合成手段によつて得られた信号のうち、所定の色フィルタの配置に対応した 色信号を取り出す信号 ^^手段を備えたことを赚とする
3 7. 複数の 分を含む 像信号を ¾ ^とした画像処«置において、 tina^®像信号に対して、 注目画素とその近傍の領域での信号レベルを得る信号 レベル取鮮段と、
tin己信号レベルに基づきノィズ 理パラメータを設定するパラメータ設定手 段と、
歸己ノィズ パラメータを用レ、て嫌己画像信号のノィズを ig るノイズ 繊手段と、
ΙβΙ^ϋ像信号に対して、 注目画素を含 ififfl 貝域での^ ^分間の相隱数を演 算する相隱欽演算手段と、
膽¾¾隱数に応じて MISノィズ パラメータを補正するパラメータ補正 手段と、 を備えたことを ^ [とする画»«
3 8. 請求項 3 7に«の 置において、
膽己ノィズ ί赚手段は、 嫌己ノィズ« ^理パラメータに応じて設定される閾値 と fS^®像信号との比 ¾ ^果に基づき |ίίΙ¾¾1ί像信号のノィズを ig» ることを | とする
3 9. 請求項 3 7または 3 8に !Ettの画 i ^Sg置において、
ffHS^ラメータ補正手段は、 嫌己注目画素を含¾©麵域での触分間の相難 数に対して負の相関を持つ補正係数を算出し、 嫌己ノィズ 理パラメータに前 S i正係数を掛けて膽己ノィズ 理パラメータを補正することを とする画 蝶¾ 4 0 · 請求項 3 9に言 2«の画飾離置にぉレ、て、
膽 E/、°ラメータ補正手段は嫌 ME係数を次式:
C = (1-r)2
r:注目画素を含む近 ί 域での 分間の相関係数
により算出することを W [とする画像処 3¾¾
4 1. 請求項 3 7から 4 0のレヽずれか一つ言 Ξ¾の画^! ¾置において、
魔 fij^®像信号は 3種以上の 分を含み、
tfit¾¾関演算処理手段は、 数として複数の"^分間の相 Π、数のうち 最 /j、値を することを とする画
4 2. 置において、
色フィルタァレイを前面に配置した 素子と、
ttHB 素子から得られた複数の触分を含む i^®像信号に対して、 注目画素と その近傍の領域での信号レベルを得る信号レベル取 段と、
ΜΙΕ信号レベルに基づきノィズ 理パラメータを設定するパラメータ設定手 段と、
膽己ノィズS^L理パラメータを用レ、て嫌 ΕΙ^Ιί像信号のノィズを ί»Τるノィ ズ ig^手段と、
嫌 a^ i像信号に対して、 tfr!E注目画素を含む近 ί赚域での触分間の相醒数 を演算する相関係数演算手段と、
tilfS目隱数に応じて ΙΞノィズ ifi^ パラメータを補正するパラメータ栖 手段と、 を備えたことを雖とする 4 3 . 請求項 4 2に言 3«の 置において、
廳己ノィズ赚手段は、 膽己ノィズ^^理パラメータに応じて設定される閾値 と tirtarai像信号との比維果に基づき te^ i像信号のノィズを ί«τることを とする 4 4 . 請求項 4 2または 4 3に記載の 置において、
Figure imgf000050_0001
ータ補正手段は、 ΙΐίΐΕ注目画素を含む近 ί魏域での触分間の相難 数に対して負の相関を持つ補正係数を算出し、 嫌己ノィズ¾»¾パラメータに前 言 Bffi正係数を掛けて tiHEノィズ iSl^ パラメータを補正することを赚とする撮 離氤
4 5 . 請求項 4 4に言 S«の 置において、
爾 S ラメータ補正手段は嫌 2«正係数を次式:
C 二 (1-r)2
r:注目画素を含む近傍領域での色成分間の相関係数
により算出することを とする衞
4 6. 請求項 4 2力、ら 4 5のレ、ずれか一つ言識の 置にぉレ、て、 ΙΗΙ^ ί像信号は 3種以上の^^分を含み、
嫌 関演算処理手段は、 ttltEtg,数として複数の 分間の相,数のうち 最 /Jィ直を することを «とする »¾go
PCT/JP2007/070770 2007-01-10 2007-10-18 画像処理装置および撮像装置 WO2008084585A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US12/498,845 US8189940B2 (en) 2007-01-10 2009-07-07 Image processing apparatus, imaging apparatus, and image processing method
US13/477,225 US8538186B2 (en) 2007-01-10 2012-05-22 Image processing apparatus, imaging apparatus, and image processing method

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2007-002655 2007-01-10
JP2007002655A JP4523008B2 (ja) 2007-01-10 2007-01-10 画像処理装置および撮像装置

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US12/498,845 Continuation US8189940B2 (en) 2007-01-10 2009-07-07 Image processing apparatus, imaging apparatus, and image processing method

Publications (1)

Publication Number Publication Date
WO2008084585A1 true WO2008084585A1 (ja) 2008-07-17

Family

ID=39608483

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2007/070770 WO2008084585A1 (ja) 2007-01-10 2007-10-18 画像処理装置および撮像装置

Country Status (3)

Country Link
US (2) US8189940B2 (ja)
JP (1) JP4523008B2 (ja)
WO (1) WO2008084585A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010211552A (ja) * 2009-03-11 2010-09-24 Rohm Co Ltd 画像処理方法及びコンピュータプログラム

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5034003B2 (ja) * 2007-06-25 2012-09-26 オリンパス株式会社 画像処理装置
JP2010147568A (ja) * 2008-12-16 2010-07-01 Olympus Corp 画像処理装置、画像処理方法、および、画像処理プログラム
KR101641543B1 (ko) 2010-02-05 2016-07-21 삼성전자주식회사 영상잡음 제거장치 및 그 방법
US8600188B2 (en) 2010-09-15 2013-12-03 Sharp Laboratories Of America, Inc. Methods and systems for noise reduction and image enhancement
US8588535B2 (en) 2010-09-15 2013-11-19 Sharp Laboratories Of America, Inc. Methods and systems for estimation of compression noise
US8532429B2 (en) * 2010-09-28 2013-09-10 Sharp Laboratories Of America, Inc. Methods and systems for noise reduction and image enhancement involving selection of noise-control parameter
US8538193B2 (en) 2010-09-28 2013-09-17 Sharp Laboratories Of America, Inc. Methods and systems for image enhancement and estimation of compression noise
JP5488621B2 (ja) 2012-01-11 2014-05-14 株式会社デンソー 画像処理装置、画像処理方法、及びプログラム
JP6287100B2 (ja) * 2013-11-20 2018-03-07 株式会社リコー 画像処理装置、画像処理方法、プログラム及び記憶媒体
WO2020127031A1 (en) * 2018-12-18 2020-06-25 Agfa Nv Method of decomposing a radiographic image into sub-images of different types

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05347723A (ja) * 1992-06-15 1993-12-27 Toshiba Corp ノイズ低減回路
JPH06178163A (ja) * 1992-12-04 1994-06-24 Casio Comput Co Ltd ノイズリダクション回路
JPH0787418A (ja) * 1993-09-13 1995-03-31 Hitachi Ltd Edtv垂直補強信号の復調方式
JPH08172552A (ja) * 1994-09-30 1996-07-02 Samsung Electron Co Ltd 映像のノイズ除去及びディテール改善装置
JPH09224172A (ja) * 1996-02-19 1997-08-26 Matsushita Electric Ind Co Ltd 映像信号処理回路

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2586687B2 (ja) * 1990-04-20 1997-03-05 松下電器産業株式会社 雑音除去回路
JP3044501B2 (ja) * 1991-11-29 2000-05-22 松下電器産業株式会社 映像信号処理方法および装置
JPH0950519A (ja) 1995-06-01 1997-02-18 Canon Inc 画像処理装置及び方法
KR0160718B1 (ko) 1995-08-29 1999-01-15 김광호 색신호에 포함된 윤곽을 보정하는 방법 및 이를 칼라 비디오기기에서 구현하기 위한 회로
JP2001157057A (ja) 1999-11-30 2001-06-08 Konica Corp 画像読取装置
JP3689607B2 (ja) 1999-12-15 2005-08-31 キヤノン株式会社 画像処理方法、装置および記憶媒体
US6904169B2 (en) * 2001-11-13 2005-06-07 Nokia Corporation Method and system for improving color images
JP4479429B2 (ja) 2003-09-09 2010-06-09 コニカミノルタホールディングス株式会社 撮像装置
JP2006023959A (ja) * 2004-07-07 2006-01-26 Olympus Corp 信号処理システム及び信号処理プログラム
JP4547223B2 (ja) 2004-09-28 2010-09-22 オリンパス株式会社 撮像システム、ノイズ低減処理装置及び撮像処理プログラム
JP5009004B2 (ja) * 2006-07-31 2012-08-22 株式会社リコー 画像処理装置、撮像装置、画像処理方法、および画像処理プログラム

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05347723A (ja) * 1992-06-15 1993-12-27 Toshiba Corp ノイズ低減回路
JPH06178163A (ja) * 1992-12-04 1994-06-24 Casio Comput Co Ltd ノイズリダクション回路
JPH0787418A (ja) * 1993-09-13 1995-03-31 Hitachi Ltd Edtv垂直補強信号の復調方式
JPH08172552A (ja) * 1994-09-30 1996-07-02 Samsung Electron Co Ltd 映像のノイズ除去及びディテール改善装置
JPH09224172A (ja) * 1996-02-19 1997-08-26 Matsushita Electric Ind Co Ltd 映像信号処理回路

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010211552A (ja) * 2009-03-11 2010-09-24 Rohm Co Ltd 画像処理方法及びコンピュータプログラム

Also Published As

Publication number Publication date
US20120293714A1 (en) 2012-11-22
JP4523008B2 (ja) 2010-08-11
US8538186B2 (en) 2013-09-17
JP2008172431A (ja) 2008-07-24
US20090324066A1 (en) 2009-12-31
US8189940B2 (en) 2012-05-29

Similar Documents

Publication Publication Date Title
WO2008084585A1 (ja) 画像処理装置および撮像装置
JP4941285B2 (ja) 撮像装置、撮像システム、撮像方法及び画像処理装置
US6108453A (en) General image enhancement framework
Bennett et al. Multispectral bilateral video fusion
CN100433785C (zh) 用于从移动场景的多个曝光中生成高动态范围图像的系统和过程
JP4455897B2 (ja) 画像処理方法および装置並びにプログラム
EP2187620B1 (en) Digital image processing and enhancing system and method with function of removing noise
JP5144202B2 (ja) 画像処理装置およびプログラム
JP2003163939A (ja) 適応モザイク減少法を実装するディジタル画像処理方法
JP2010531075A (ja) パンクロマチック画像を用いたノイズ低減されたカラーの画像
JP2010509841A (ja) パンクロマチックおよびカラーの画像のノイズ低減
CN101273623A (zh) 摄像装置
JP5012967B2 (ja) 画像処理装置及び方法、並びにプログラム
CN103873836B (zh) 图像处理设备和图像处理方法
WO2016114148A1 (ja) 画像処理装置、画像処理方法、及び記録媒体
JP2006180268A (ja) 画像処理装置、画像処理方法、プログラム、及び記録媒体
Wei et al. DA-DRN: A degradation-aware deep Retinex network for low-light image enhancement
JP2005517315A (ja) 色域圧縮
JP4456096B2 (ja) ノイズ除去方法および装置並びにプログラム
EP3547251B1 (en) Dynamic range extension of partially clipped pixels in captured images
JP5178933B1 (ja) 画像処理装置
JP5277388B2 (ja) 画像処理装置および撮像装置
WO2009081709A1 (ja) 画像処理装置、画像処理方法、および、画像処理プログラム
JP2008021219A (ja) 画像処理装置及び方法
Yang et al. The spatial correlation problem of noise in imaging deblurring and its solution

Legal Events

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

Ref document number: 07830504

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 07830504

Country of ref document: EP

Kind code of ref document: A1