WO2017017742A1 - 画像処理装置、画像処理方法およびプログラム - Google Patents

画像処理装置、画像処理方法およびプログラム Download PDF

Info

Publication number
WO2017017742A1
WO2017017742A1 PCT/JP2015/071175 JP2015071175W WO2017017742A1 WO 2017017742 A1 WO2017017742 A1 WO 2017017742A1 JP 2015071175 W JP2015071175 W JP 2015071175W WO 2017017742 A1 WO2017017742 A1 WO 2017017742A1
Authority
WO
WIPO (PCT)
Prior art keywords
noise
pixel
unit
image
image data
Prior art date
Legal status (The legal status 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 status listed.)
Ceased
Application number
PCT/JP2015/071175
Other languages
English (en)
French (fr)
Inventor
裕輝 丸山
学 市川
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Olympus Corp
Original Assignee
Olympus Corp
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 Olympus Corp filed Critical Olympus Corp
Priority to JP2017530484A priority Critical patent/JP6727205B2/ja
Priority to PCT/JP2015/071175 priority patent/WO2017017742A1/ja
Priority to DE112015006619.5T priority patent/DE112015006619T5/de
Priority to CN201580081299.0A priority patent/CN107710734B/zh
Publication of WO2017017742A1 publication Critical patent/WO2017017742A1/ja
Priority to US15/845,166 priority patent/US10368017B2/en
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/40Picture signal circuits
    • H04N1/409Edge or detail enhancement; Noise or error suppression
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; 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 OR CALCULATING; 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/60Control of cameras or camera modules
    • H04N23/66Remote control of cameras or camera parts, e.g. by remote control devices
    • H04N23/663Remote control of cameras or camera parts, e.g. by remote control devices for controlling interchangeable camera parts based on electronic image sensor signals
    • 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
    • H04N23/843Demosaicing, e.g. interpolating colour pixel values
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/60Noise processing, e.g. detecting, correcting, reducing or removing noise
    • H04N25/616Noise processing, e.g. detecting, correcting, reducing or removing noise involving a correlated sampling function, e.g. correlated double sampling [CDS] or triple sampling
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/60Noise processing, e.g. detecting, correcting, reducing or removing noise
    • H04N25/618Noise processing, e.g. detecting, correcting, reducing or removing noise for random or high-frequency noise
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/60Noise processing, e.g. detecting, correcting, reducing or removing noise
    • H04N25/68Noise processing, e.g. detecting, correcting, reducing or removing noise applied to defects
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/70SSIS architectures; Circuits associated therewith
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/70SSIS architectures; Circuits associated therewith
    • H04N25/76Addressed sensors, e.g. MOS or CMOS sensors
    • H04N25/78Readout circuits for addressed sensors, e.g. output amplifiers or A/D converters
    • 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
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20216Image averaging

Definitions

  • the present invention relates to an image processing apparatus, an image processing method, and a program for correcting blinking defect noise in which a pixel value such as RTS noise generated in an image sensor having a plurality of pixels arranged two-dimensionally varies within a certain range. .
  • CMOS Complementary Metal Oxide Semiconductor
  • CMOS Complementary Metal Oxide Semiconductor
  • a shared pixel structure that reads a signal by sharing a plurality of pixels with one readout circuit reduces the area required for the circuit in the image sensor, and the aperture ratio ( The sensitivity is improved by improving the ratio of the light receiving part).
  • the noise generated in the image sensor includes a defective pixel whose pixel value always shows an abnormal value, and a pixel value other than a dark current shot noise caused by a dark current and a random noise caused by a thermal noise in a readout circuit.
  • RTS Random Telegraph Signal
  • As a technique for correcting this RTS noise a pixel value of a target pixel in a captured image, a pixel value of a peripheral pixel of the target pixel, and a noise level of RTS noise detected in advance for each pixel of an image sensor (hereinafter referred to as “RTS noise”).
  • RTS noise a technique for correcting this RTS noise, a pixel value of a target pixel in a captured image, a pixel value of a peripheral pixel of the target pixel, and a noise level of RTS noise detected in advance for each pixel of an image sensor (hereinafter referred to as “RTS noise”).
  • Patent Document 1 when the correlation between the pixel value of the target pixel and the pixel values of the peripheral pixels of the target image is high, RTS noise can be reduced.
  • the image quality may be deteriorated because the correlation between the pixel value of the pixel of interest and the pixel value of surrounding pixels is low.
  • the present invention has been made in view of the above, and an image processing apparatus and image processing capable of improving image quality even when the pixel value of a pixel of interest and the pixel value of peripheral pixels are low
  • An object is to provide a method and a program.
  • an image processing apparatus is arranged in a two-dimensional manner, and receives a plurality of pixels that receive light from the outside and generate a signal according to the amount of received light.
  • An image processing apparatus that corrects blinking defect noise included in image data generated by an image sensor having a plurality of readout circuits that read out the signals as pixel values, and the blinking defect noise caused by the readout circuit is generated
  • An acquisition unit that acquires position information of the readout circuit or position information of each of the plurality of pixels and the image data, and the image data acquired by the acquisition unit is image data to be corrected
  • the correction target image data based on the reference image data based on the image data acquired at a different time from the correction target image data and the noise information.
  • a correcting unit for correcting the pixel value of kicking the target pixel characterized by comprising a.
  • the image processing apparatus further includes a movement amount calculation unit that calculates a movement amount of a subject based on the correction target image data and the reference image data. Based on the movement amount calculated by the movement amount calculation unit, a reference pixel of the reference image data corresponding to the pixel of interest or a pixel near the pixel of interest is obtained, and based on a pixel value of the reference pixel, The pixel value of the target pixel is corrected.
  • the image processing apparatus is the image processing apparatus according to the above invention, wherein the correction unit reduces noise in the target pixel when the flashing defect noise occurs in the target pixel as compared with a case where the flashing defect noise does not occur. It is characterized by increasing the strength of processing.
  • the correction unit increases the strength of noise reduction processing in the time direction using the reference image data when the blinking defect noise occurs in the target pixel. It is characterized by doing.
  • the image processing apparatus is the image processing apparatus according to the above invention, wherein, when the blinking defect noise occurs in the target pixel, the correction unit performs a spatial noise reduction process using pixels around the target pixel. It is characterized by increasing.
  • the correction unit calculates a representative value corresponding to a pixel value when the blinking defect noise does not occur, based on a pixel value of the reference pixel, The pixel value of the target pixel is corrected based on the representative value.
  • the correction unit calculates the representative value based on a pixel value of the pixel of interest.
  • the correction unit is based on a pixel value of the reference pixel in which the blinking defect noise does not occur or a pixel value of the reference pixel in which the blinking defect noise is corrected.
  • the representative value is calculated.
  • the correction unit calculates the representative value based on pixel values of surrounding pixels in the target pixel in which the blinking defect noise does not occur.
  • the image processing apparatus further includes a noise amount estimation unit that estimates a random noise amount around the target pixel in the above invention, and the correction unit is configured to perform the above-described correction based on the random noise amount. A representative value is calculated.
  • the noise information corresponds to position information of the reading circuit or position information of each of the plurality of pixels in which blinking defect noise caused by the reading circuit is generated. And further includes a blinking defect noise level that is a noise level of the blinking defect noise, and the correction unit corrects the pixel value of the pixel of interest using the blinking defect noise level.
  • the noise information corresponds to position information of the reading circuit or position information of each of the plurality of pixels in which blinking defect noise caused by the reading circuit is generated. And further includes a blinking defect noise level that is a noise level of the blinking defect noise, and the correction unit increases the intensity of noise reduction processing in the target pixel as the blinking defect noise level increases.
  • the noise information corresponds to position information of the reading circuit or position information of each of the plurality of pixels in which blinking defect noise caused by the reading circuit is generated.
  • a blinking defect noise level that is a noise level of the blinking defect noise is further included, and the correction unit calculates the representative value based on the blinking defect noise level.
  • the image processing apparatus further includes a reliability calculation unit that calculates a reliability indicating the certainty of the movement amount calculated by the movement amount calculation unit in the above invention, and the correction unit includes the reliability The representative value is calculated based on the degree.
  • the reference image data is the image data generated by the image sensor immediately before the image data corresponding to the correction target image data. To do.
  • the reference image data is image data corrected by the correction unit immediately before the correction unit corrects the correction target image data. It is characterized by.
  • the reference image data is a plurality of the image data generated by the imaging element before the image data corresponding to the correction target image data. It is characterized by.
  • the blinking defect noise is random telegraph signal noise.
  • the image processing method is arranged in a two-dimensional manner, receives a light from the outside, generates a signal corresponding to the amount of received light, and a plurality of readout circuits that read the signal as a pixel value And an image processing method executed by an image processing apparatus that corrects blinking defect noise included in image data generated by an image pickup device having a position of the readout circuit where the blinking defect noise caused by the readout circuit is generated Information or noise information including position information of each of the plurality of pixels and the image data, and the correction target image data using the image data acquired in the acquisition step as the correction target image data.
  • the correction target image data based on the reference image data based on the image data acquired at a different time and the noise information Characterized in that it comprises a correction step of correcting the pixel value of the definitive target pixel.
  • a program according to the present invention is arranged in a two-dimensional manner, receives a light from the outside, generates a signal according to the amount of received light, a plurality of readout circuits that read the signal as a pixel value,
  • the image processing apparatus that corrects the blinking defect noise included in the image data generated by the image pickup device having the position information of the readout circuit or the position of each of the plurality of pixels in which the blinking defect noise due to the readout circuit is generated
  • a correction step that, characterized in that for the execution.
  • the image quality can be improved even when the pixel value of the target pixel and the pixel values of the surrounding pixels are low in correlation.
  • FIG. 1 is a block diagram schematically showing a configuration of an imaging system according to Embodiment 1 of the present invention.
  • FIG. 2 is a schematic diagram schematically showing the configuration of the main part of the image sensor included in the imaging apparatus according to Embodiment 1 of the present invention.
  • FIG. 3 is a diagram illustrating an example of fluctuations in the amplifier output output from the amplifier unit when RTS noise is generated in the case where the image sensor according to Embodiment 1 of the present invention is shielded from light so that light does not strike. It is.
  • FIG. 4 is a diagram illustrating an example of a distribution of pixel values read using the amplifier unit that generates RTS noise according to the first embodiment of the present invention.
  • FIG. 1 is a block diagram schematically showing a configuration of an imaging system according to Embodiment 1 of the present invention.
  • FIG. 2 is a schematic diagram schematically showing the configuration of the main part of the image sensor included in the imaging apparatus according to Embodiment 1 of the present invention.
  • FIG. 3 is
  • FIG. 5 is a flowchart showing an outline of processing executed by the image processing apparatus according to Embodiment 1 of the present invention.
  • FIG. 6 is a diagram illustrating an example of a noise model stored in a lookup table provided in the image processing apparatus according to Embodiment 1 of the present invention.
  • FIG. 7 is a flowchart showing an overview of the time direction NR process of FIG.
  • FIG. 8 is a diagram showing a relationship between RTS_Value and coefficient Co_we_pre set by the noise reduction processing unit included in the image processing apparatus according to Embodiment 1 of the present invention.
  • FIG. 9 is a flowchart showing an outline of the spatial direction NR process of FIG. FIG.
  • FIG. 10 is a diagram showing a relationship between RTS_Value and coefficient Co_Cb set by the noise reduction processing unit included in the image processing apparatus according to Embodiment 1 of the present invention.
  • FIG. 11 is a block diagram schematically showing the configuration of the imaging system according to Embodiment 2 of the present invention.
  • FIG. 12 is a block diagram showing a detailed configuration of the RTS noise correction unit according to Embodiment 2 of the present invention.
  • FIG. 13 is a block diagram showing a detailed configuration of the representative value calculation unit according to Embodiment 2 of the present invention.
  • FIG. 14 is a block diagram showing a detailed configuration of a correction value calculation unit according to Embodiment 2 of the present invention.
  • FIG. 11 is a block diagram schematically showing the configuration of the imaging system according to Embodiment 2 of the present invention.
  • FIG. 12 is a block diagram showing a detailed configuration of the RTS noise correction unit according to Embodiment 2 of the present invention.
  • FIG. 13 is a
  • FIG. 15 is a flowchart showing an outline of processing executed by the image processing apparatus according to Embodiment 2 of the present invention.
  • FIG. 16 is a flowchart showing an overview of the representative value calculation processing of FIG.
  • FIG. 17 is a flowchart showing an outline of the correction value calculation processing of FIG.
  • FIG. 18 is a block diagram showing a detailed configuration of the noise reduction unit according to Embodiment 3 of the present invention.
  • FIG. 19 is a block diagram showing a detailed configuration of the RTS noise correction unit according to Embodiment 3 of the present invention.
  • FIG. 20 is a flowchart showing an outline of representative value calculation processing executed by the image processing apparatus according to Embodiment 3 of the present invention.
  • FIG. 21 is a block diagram showing a detailed configuration of an RTS noise correction unit according to a modification of the third embodiment of the present invention.
  • FIG. 22 is a flowchart showing an overview of representative value calculation processing according to a modification of the third embodiment of the present invention.
  • FIG. 23 is a diagram schematically illustrating a blend rate calculation method calculated by the representative value determination unit according to the modification of the third embodiment of the present invention.
  • FIG. 24 is a block diagram showing a detailed configuration of the noise reduction unit according to Embodiment 4 of the present invention.
  • FIG. 25 is a block diagram showing a detailed configuration of the RTS noise correction unit according to Embodiment 4 of the present invention.
  • FIG. 22 is a flowchart showing an overview of representative value calculation processing according to a modification of the third embodiment of the present invention.
  • FIG. 23 is a diagram schematically illustrating a blend rate calculation method calculated by the representative value determination unit according to the modification of the third embodiment of the present invention.
  • FIG. 24 is a block diagram
  • FIG. 26 is a flowchart showing an outline of representative value calculation processing executed by the image processing apparatus according to Embodiment 4 of the present invention.
  • FIG. 27 is a block diagram showing a detailed configuration of the RTS noise correction unit according to the first modification of the fourth embodiment of the present invention.
  • FIG. 28 is a flowchart showing an outline of representative value calculation processing executed by the image processing apparatus according to the first modification of the fourth embodiment of the present invention.
  • FIG. 29 is a flowchart showing an outline of representative value calculation processing executed by the image processing apparatus according to the second modification of the fourth embodiment of the present invention.
  • FIG. 30 is a flowchart showing an outline of representative value calculation processing executed by the image processing apparatus according to the third modification of the fourth embodiment of the present invention.
  • FIG. 31 is a block diagram schematically showing the configuration of the imaging system according to Embodiment 5 of the present invention.
  • FIG. 32 is a flowchart showing an outline of processing executed by the imaging system according to Embodiment 5 of the present invention.
  • FIG. 33 is a flowchart showing an outline of the image processing of FIG.
  • FIG. 1 is a block diagram schematically showing a configuration of an imaging system according to Embodiment 1 of the present invention.
  • the imaging system 1 illustrated in FIG. 1 includes an imaging device 10, an image processing device 30, and a display device 40.
  • the imaging apparatus 10 includes an optical system 101, a diaphragm 102, a shutter 103, a driver 104, an imaging element 105, an analog processing unit 106, an A / D conversion unit 107, and an operation unit. 108, a memory I / F unit 109, a recording medium 110, a volatile memory 111, a nonvolatile memory 112, a bus 113, an imaging control unit 114, and a first external I / F unit 115.
  • the optical system 101 is configured using one or a plurality of lenses.
  • the optical system 101 is configured using, for example, a focus lens and a zoom lens.
  • the diaphragm 102 adjusts exposure by limiting the amount of incident light collected by the optical system 101.
  • the diaphragm 102 limits the amount of incident light collected by the optical system 101 under the control of an imaging control unit 114 described later. Note that the amount of incident light may be limited using an electronic shutter in the shutter 103 or the image sensor 105 without using the diaphragm 102. Note that the optical system 101 and the diaphragm 102 may be detachable from the imaging apparatus 10.
  • the shutter 103 sets the state of the image sensor 105 to an exposure state or a light shielding state.
  • the shutter 103 is configured using, for example, a focal plane shutter. Note that an electronic shutter in the image sensor 105 may be used without using the shutter 103.
  • the driver 104 drives the optical system 101, the diaphragm 102, and the shutter 103 under the control of an imaging control unit 114 described later.
  • the driver 104 moves the optical system 101 along the optical axis O1 to change the zoom magnification of the imaging apparatus 10 or adjust the focus position.
  • the imaging element 105 receives light collected by the optical system 101 under the control of an imaging control unit 114 described later, converts the light into image data (electric signal), and outputs the image data.
  • the image sensor 105 is configured using a CMOS (Complementary Metal Oxide Semiconductor) or the like in which a plurality of pixels are two-dimensionally arranged.
  • a Bayer array RGB filter is arranged in front of each pixel. Note that the image sensor 105 is not limited to the Bayer array, and may of course be a stacked type such as Fovion.
  • the filter to be used is not limited to RGB, and any filter such as a complementary color filter can be applied.
  • a light source capable of irradiating different color lights in a time-sharing manner is arranged, and no color filter is arranged in the image sensor 105, and a color image is configured using images that are sequentially captured while changing the irradiating color. You may be able to do it.
  • the image sensor 105 has an electronic shutter function that can electronically control the amount of received light.
  • FIG. 2 is a schematic diagram schematically showing the configuration of the main part of the image sensor 105.
  • 2 shows an example in which a readout circuit is shared by a plurality of pixels in order to improve sensitivity by improving the aperture ratio of the pixels.
  • 2 includes one readout circuit for 8 pixels of 2 pixels in the horizontal direction (horizontal direction) ⁇ 4 pixels in the vertical direction (vertical direction).
  • FIG. 2 an example in which one readout circuit is grouped with respect to 8 pixels of 2 pixels in the horizontal direction (horizontal direction) ⁇ 4 pixels in the vertical direction (vertical direction) has been described. It is assumed that the above-described pixels and readout circuits are arranged side by side in the horizontal direction and the vertical direction on the image sensor 105 of the first embodiment.
  • the image sensor 105 receives light by exposure and performs photoelectric conversion to generate a plurality of pixels 105a (photodiodes) that generate charges corresponding to the exposure amount, and a plurality of pixels 105a.
  • a first switch 105b that is provided in each and opens and closes according to the control of the imaging control unit 114, a vertical transfer line 105c that transfers a signal (charge) output from each of the plurality of pixels 105a in the vertical direction, and a plurality of An FD unit 105d (Floating Diffusion) that accumulates signals output from each of the pixels 105a, an amplifier unit 105e that amplifies signals output from the FD unit 105d, and a second that opens and closes according to the control of the imaging control unit 114
  • a switch 105f, a control line 105g for controlling the second switch 105f, and a transfer line 10 for transferring the electric signal amplified by the amplifier unit 105e. Includes a h, the.
  • the image sensor 105 configured as described above reads out a signal corresponding to the exposure amount in the pixels 105a (1) to 105a (8) as a pixel value
  • the FD unit 105d is set in a reset state, and the imaging control unit 114 is reset.
  • the imaging control unit 114 turns on the second switch 105f, the imaging element 105 amplifies the charge accumulated in the FD unit 105d by the amplifier unit 105e and reads (outputs) the pixel value.
  • the imaging element 105 resets the FD unit 105d, and the imaging control unit 114 turns on only the first switch 105b (2), whereby the charge generated in the pixel 105a (2) is transferred to the FD unit 105d. Forward to. Thereafter, when the imaging control unit 114 turns on the second switch 105f, the imaging element 105 amplifies the charge accumulated in the FD unit 105d by the amplifier unit 105e and reads it as a pixel value.
  • the imaging element 105 can sequentially output signals corresponding to the exposure amounts in the pixels 105a (1) to 105a (8) as pixel values by sequentially performing such readout operations.
  • the amplifier unit 105e functions as a readout circuit that reads out charges from each of the plurality of pixels 105a.
  • the analog processing unit 106 performs predetermined analog processing on the analog signal input from the image sensor 105 and outputs the analog signal to the A / D conversion unit 107. Specifically, the analog processing unit 106 performs noise reduction processing, gain increase processing, and the like on the analog signal input from the image sensor 105. For example, the analog processing unit 106 performs waveform shaping on the analog signal after reducing reset noise and the like, and further increases the gain so that the target brightness is obtained.
  • the A / D conversion unit 107 generates digital image data (hereinafter referred to as “RAW image data”) by performing A / D conversion on the analog signal input from the analog processing unit 106, and the bus 113. Is output to the volatile memory 111.
  • RAW image data digital image data
  • the A / D conversion unit 107 may directly output RAW image data to each unit of the imaging apparatus 10 to be described later.
  • the above-described analog processing unit 106 and A / D conversion unit 107 may be provided in the image sensor 105 so that the image sensor 105 directly outputs digital RAW image data.
  • the operation unit 108 gives various instructions of the imaging apparatus 10.
  • the operation unit 108 includes a power switch that switches the power state of the imaging device 10 to an on state or an off state, a release switch that gives a still image shooting instruction, an operation switch that switches various settings of the imaging device 10, and a moving image shooting. And a moving image switch for giving instructions.
  • the recording medium 110 is configured using a memory card that is mounted from the outside of the imaging device 10, and is detachably mounted on the imaging device 10 via the memory I / F unit 109. Further, the recording medium 110 may output a program and various types of information to the nonvolatile memory 112 via the memory I / F unit 109 under the control of an imaging control unit 114 described later.
  • the volatile memory 111 temporarily stores image data input from the A / D conversion unit 107 via the bus 113.
  • the volatile memory 111 temporarily stores image data that the image sensor 105 sequentially outputs for each frame via the analog processing unit 106, the A / D conversion unit 107, and the bus 113.
  • the volatile memory 111 is configured using SDRAM (Synchronous Dynamic Random Access Memory) or the like.
  • the nonvolatile memory 112 is configured using a flash memory or the like, and records various programs for operating the imaging apparatus 10 and various data used during execution of the programs.
  • the nonvolatile memory 112 includes RTS noise caused by the program recording unit 112a and the position information of the readout circuit (amplifier unit 105e) that reads out the pixel value or the position information of each of the plurality of pixels 105a and the readout circuit (amplifier unit 105e).
  • An RTS noise information recording unit 112b that records RTS noise position information associated with feature quantities related to the random noise model, and a random noise model information recording unit 112c that records one or more random noise models.
  • the feature amount is one of the amplitude of RTS noise (RTS_Value), the occurrence frequency of blinking defect noise, and the occurrence frequency of RTS noise less than the amplitude of RTS noise.
  • the bus 113 is configured by using a transmission path or the like that connects each component of the imaging device 10, and transfers various data generated inside the imaging device 10 to each component of the imaging device 10.
  • the imaging control unit 114 is configured by using a CPU (Central Processing Unit) or the like, and performs an instruction, data transfer, or the like to each unit constituting the imaging apparatus 10 in accordance with an instruction signal or a release signal from the operation unit 108.
  • the operation of the apparatus 10 is comprehensively controlled.
  • the imaging control unit 114 starts controlling the shooting operation in the imaging device 10.
  • the photographing operation in the image pickup apparatus 10 means that the analog processing unit 106 and the A / D conversion unit 107 are predetermined with respect to the exposure timing of the image pickup element 105, the output timing of the analog signal, and the analog signal output from the image pickup element 105. This is the operation to perform the process.
  • the image data thus processed is recorded on the recording medium 110 via the bus 113 and the memory I / F unit 109 under the control of the imaging control unit 114.
  • the first external I / F unit 115 outputs information input from an external device via the bus 113 to the nonvolatile memory 112 or the volatile memory 111, while the volatile memory 111 stores the information to an external device via the bus 113.
  • Information, information stored in the nonvolatile memory 112, and image data generated by the image sensor 105 are output.
  • the first external I / F unit 115 outputs RAW image data, RTS noise information, and random noise model information generated by the image sensor 105 to the image processing apparatus 30 via the bus 113.
  • the image processing device 30 includes a third external I / F unit 31, an interpolation processing unit 32, a frame memory 33, a motion detection unit 34, a noise reduction unit 35, a storage unit 36, an operation unit 37, an image A processing unit 38 and an image processing control unit 39 are provided.
  • the third external I / F unit 31 records the image data generated by the image sensor 105 via the first external I / F unit 115 of the imaging apparatus 10 and the RTS recorded by the RTS noise information recording unit 112b in the nonvolatile memory 112.
  • RTS noise information related to noise and random noise model information recorded by the random noise model information recording unit 112c are acquired, and the acquired raw image data (correction target image data), RTS noise information, and random noise model information are interpolated by the interpolation processing unit 32 and The data is output to the storage unit 36.
  • the third external I / F unit 31 and the first external I / F unit 115 are connected via, for example, a control cable capable of exchanging information bidirectionally or wireless communication.
  • the third external I / F unit 31 functions as an acquisition unit.
  • the interpolation processing unit 32 performs an interpolation process on the RAW image data generated by the image sensor 105 input from the third external I / F unit 31. Specifically, when the image sensor 105 is a Bayer array, the interpolation processing unit 32 performs an interpolation process on each pixel of a RAW image (hereinafter, “current image”) corresponding to the RAW image data. The missing pixel value (signal value) is interpolated to generate a current image having all the pixel values of the R, G, and B signals in each pixel, and this current image is output to the noise reduction unit 35.
  • current image a RAW image
  • the frame memory 33 is bi-directionally connected to the noise reduction unit 35 and temporarily stores frame data of an image to be displayed on the display device 40.
  • the motion detection unit 34 stores the current image input from the interpolation processing unit 32 and the frame memory 33, the current image whose noise is corrected by the noise reduction unit 35, and the corrected image that is temporally continuous with the current image. Based on the corrected image corresponding to the data (hereinafter referred to as “past image”), the motion vector (movement amount) of the subject is detected, and the detection result is output to the noise reduction unit 35.
  • the motion detector 34 detects a motion vector using, for example, a known block matching process.
  • the motion detector 34 detects the feature amount and the feature amount of the target pixel of the current image.
  • the pixel of the past image having the minimum distance to is detected, and the motion vector is detected based on the position of the detected pixel and the pixel position of the target pixel of the current image.
  • the feature amount of the target pixel of the current image and the feature amount of the past image corresponding to the feature amount position of the target pixel of the current image are a pixel value, a luminance value, and the like.
  • the motion detector 34 may detect a motion vector for each pixel of the current image. In the first embodiment, the motion detector 34 functions as a movement amount calculator.
  • the noise reduction unit 35 performs noise reduction processing (hereinafter referred to as “NR processing”) on the current image input from the interpolation processing unit 32. Specifically, the noise reduction unit 35 determines the still state and the operation state for each pixel of the current image, and switches the NR process according to the determination result. For example, when the current image is determined to be still, the noise reduction unit 35 selects a first NR process (hereinafter referred to as “time direction NR process”) that reduces noise in the time direction in which a high-frequency component can be held. On the other hand, when it is determined that the current image is in the operating state, a second NR process (hereinafter referred to as “spatial direction NR process”) that reduces noise in the spatial direction is selected.
  • NR processing noise reduction processing
  • the time-direction NR processing is to reduce noise by weighted average processing using a current image to be subjected to NR processing and a past image (reference image data) acquired at a time different from the current image.
  • a pixel of interest (processing target pixel) of the current image and a pixel of the past image (processing target pixel) corresponding to the position of the pixel of interest of the current image are used.
  • noise is reduced by a weighted average process using a target pixel (processing target pixel) to be processed by the NR process and peripheral pixels around the target pixel.
  • the stationary state indicates a state in which the relative positional relationship between the imaging element 105 (imaging device 10) and the subject does not change with time. Furthermore, the operating state indicates a state in which the relative positional relationship between the imaging element 105 (imaging device 10) and the subject changes with time.
  • the noise reduction unit 35 functions as a correction unit. In addition, when the subject is not moving, the noise reduction unit 35 acquires the pixel value of the pixel of the past image corresponding to the pixel near the target pixel of the current image from the storage unit 36, and the pixel of the pixel of the acquired past image The pixel value of the target pixel may be corrected based on the value. *
  • the noise reduction unit 35 includes an evaluation value calculation unit 351, an estimated noise amount acquisition unit 352, a lookup table 353, a determination unit 354, and a noise reduction processing unit 355.
  • the evaluation value calculation unit 351 determines whether or not the subject included in each of the temporally continuous current image input from the interpolation processing unit 32 and the past image input from the frame memory 33 is in a stationary state between frames. An evaluation value for determination is calculated.
  • the estimated noise amount acquisition unit 352 acquires an estimated noise amount estimated from an image corresponding to the image data input from the interpolation processing unit 32.
  • the lookup table 353 stores a noise model used when the estimated noise amount acquisition unit 352 estimates the noise amount.
  • the determination unit 354 determines that the subject in the current image input from the interpolation processing unit 32 is stationary. It is determined whether or not there is.
  • the noise reduction processing unit 355 corrects the pixel value of the target pixel by increasing the noise reduction intensity and outputs the corrected value to the image processing unit 38.
  • the storage unit 36 is configured using a volatile memory or a non-volatile memory, and stores RTS noise information and a plurality of past images output from the third external I / F unit 31.
  • the operation unit 37 accepts input of various operation signals related to the image processing device 30.
  • the operation unit 37 is configured using, for example, a cross button, a push button, a touch panel, and the like.
  • the image processing unit 38 performs predetermined image processing on the noise-reduced RAW image data and outputs it to the display device 40.
  • the predetermined image processing performs basic image processing including at least white balance adjustment processing, ⁇ correction processing, color reproduction processing, edge enhancement processing, and the like.
  • the image processing unit 38 performs image processing for reproducing a natural image based on preset image processing parameters.
  • the parameters of each image processing are values of contrast, sharpness, saturation, white balance, and gradation.
  • the image processing control unit 39 comprehensively controls each unit constituting the image processing apparatus 30.
  • the image processing control unit 39 is configured using a CPU (Central Processing Unit) or the like.
  • the image processing control unit 39 controls the transfer of instructions, data, and the like of each unit constituting the image processing apparatus 30.
  • the display device 40 displays an image corresponding to the image data input from the image processing device 30.
  • the display device 40 is configured using a display panel such as liquid crystal or organic EL (Electro Luminescence).
  • the image processing device 30 corrects RTS noise generated in the imaging element 105, and the display device 40 displays an image corresponding to the image data subjected to image processing by the image processing device 30. To do.
  • FIG. 3 is a diagram illustrating an example of fluctuations in the amplifier output output from the amplifier unit 105e when RTS noise is generated when light is blocked so that the image sensor 105 does not receive light.
  • FIG. 4 is a diagram illustrating an example of a distribution of pixel values read using the amplifier unit 105e that generates RTS noise.
  • the RTS noise is generated when charges are trapped or released at a trap timing at random timing. For this reason, as shown in FIG. 3, in the amplifier unit 105e where RTS noise occurs, the amplifier output fluctuates randomly within a range of about Vrts. In addition, the potential fluctuation does not occur in an instant and requires a short time ⁇ .
  • a correlated double sampling process (hereinafter referred to as “CDS process”) is performed in order to reduce noise from the pixel value read from the pixel 105a.
  • the imaging control unit 114 turns on the reset switch (not shown) of the imaging element 105 to reset the charge of the FD unit 105d, and the imaging control unit 114 turns on the second switch 105f, A reset state is created, and a signal (reference signal) in the reset state is read (output).
  • the imaging control unit 114 turns on only the first switch 105b (or any one of the first switches 105b (1) to 105b (8)), and charges generated in the pixel 105a are transferred to the FD unit 105d. Then, a read state (output state) is generated in which the second switch 105f is turned on, and a signal in the read state is read (output). Subsequently, in the CDS process, a signal obtained by subtracting the reset state signal (reference signal) from the readout state signal is converted as a pixel value.
  • the image sensor 105 when the image sensor 105 reads signals at time tr1 (reset state) and time ts1 (readout state) by CDS processing, the amplifier outputs V at time tr1 and time ts1 are substantially the same. For this reason, the pixel values read out are mainly influenced by random noise, and the read pixel value has a distribution centered on 0 like the distribution A shown in FIG. Similarly, the image sensor 105 has substantially the same amplifier output V at time tr2 and time ts2 at time tr2 (reset state) and time ts2 (readout state). It becomes like the distribution A shown.
  • the amplifier output at time ts3 is about Vrts lower than the amplifier output at time tr3.
  • the pixel value is shifted in the negative direction by RTS_Value corresponding to the amount of change in the amplifier output Vrts, and the read pixel value becomes a distribution B centered on -RTS_Value.
  • the amplifier output at time ts4 is about Vrts compared to the amplifier output at time tr4. Since the difference between the two signals is high, the pixel value is shifted in the positive direction by RTS_Value corresponding to the amount of change in the amplifier output Vrts, and the read pixel value becomes a distribution C centered on RTS_Value. .
  • the signal may be read out while the potential is fluctuating.
  • the amplifier output difference is larger than ⁇ Vrts and smaller than Vrts between the reset state readout time and the readout state readout time.
  • the pixel value read from the image sensor 105 is also larger than ⁇ RTS_Value and smaller than RTS_Value.
  • the time ⁇ is considered to be substantially constant if the conditions (for example, temperature, drive voltage, etc.) of the image sensor 105 are constant. Therefore, pixel values larger than ⁇ RTS_Value and smaller than RTS_Value are generated with the same probability.
  • the frequency of occurrence of these pixel values is defined as ⁇ noise.
  • each of the distribution B and the distribution C is different only in the median value, but other than that is the same distribution. Therefore, hereinafter, the ratio of the distribution B or the distribution C to the distribution A is defined as ⁇ rts. This ⁇ rts increases as the fluctuation cycle of the amplifier output of the amplifier unit 105e is shorter.
  • the pixel values read using the amplifier unit 105e in which RTS noise is generated by the CDS process have a distribution as shown in FIG.
  • the potential in the readout state changes according to the exposure amount.
  • the potential change due to the RTS noise is constant regardless of the exposure amount. That is, the RTS noise does not depend on the exposure amount, and has a characteristic that randomly varies with respect to a normal pixel value in a range of ⁇ RTS_Value or more and RTS_Value or less.
  • the distribution A, the distribution B, and the distribution C are schematically shown.
  • the RTS noise is noise caused by the readout circuit (amplifier unit 105e), as shown in FIG. 2, when each of the plurality of pixels 105a shares one readout circuit, all the shared pixels RTS noise having the same characteristics occurs in (pixels 105a (1) to 105a (8)).
  • RTS noise may also occur in a column amplifier or a source follower shared in the column direction of the image sensor 105.
  • RTS noise having the same characteristics is generated in all the pixels in the column direction sharing the same column amplifier and source follower.
  • the present invention can also be applied to RTS noise generated in a circuit other than the readout circuit (amplifier unit 105e).
  • the RTS noise is a blinking defect in which when the subject is fixed and the image is captured under the same conditions, the pixel value of the image obtained by the image capture is amplitude (fluctuates) within a certain range ( ⁇ RTS_Value or more, RTS_Value or less). It becomes a kind of noise.
  • FIG. 5 is a flowchart showing an outline of processing executed by the image processing apparatus 30, and is a flowchart of a main routine executed by the image processing apparatus 30.
  • the evaluation value calculation unit 351 uses the current image (current frame) output from the interpolation processing unit 32 and the past image (past frame) stored in the frame memory 33.
  • the difference average value mSAD (inter-frame difference value) is calculated (step S1).
  • the past image is an NR-processed image output from the noise reduction unit 35 at a time point one frame before the current image input from the interpolation processing unit 32.
  • the evaluation value calculation unit 351 calculates the difference average value mSAD using the following equation (1), where (x, y) is the coordinate of the pixel of interest that is the pixel to be subjected to NR processing.
  • min ⁇ represents a process of acquiring the minimum value in the parentheses.
  • F G — cur (x, y) is the G signal value at the coordinates (x, y) of the current image (RGB image)
  • F G — pre (x, y) is the coordinates (x of the past image)
  • Y) is the G signal value.
  • k is a natural number, and (2k + 1) corresponds to the kernel size when calculating the difference average value mSAD. Note that k may be set to a predetermined value in advance, or a user may set an arbitrary value via the operation unit 37.
  • represents a process of acquiring the absolute value of the real number A.
  • the smallest SAD in the search range is selected as mSAD.
  • the estimated noise amount acquisition unit 352 acquires an estimated noise amount N corresponding to the pixel value (signal value) of the target pixel in the current image with reference to the lookup table 353 (step S2).
  • FIG. 6 is a diagram illustrating an example of a noise model stored in the lookup table 353.
  • the vertical axis indicates the amount of noise
  • the horizontal axis indicates the pixel value.
  • the standard deviation of the pixel value is used as the amount of noise on the vertical axis, and a noise model corresponding to the characteristics of the image sensor 105 is shown.
  • the estimated noise amount acquisition unit 352 acquires the estimated noise amount N corresponding to the pixel value of the target pixel of the current image with reference to the noise model of the curve L x 1 in FIG. To do.
  • the noise model may have a characteristic approximated by an approximate expression or a broken line.
  • the determination unit 354 determines whether or not the difference average value mSAD calculated by the evaluation value calculation unit 351 is equal to or less than the noise amount N acquired by the estimated noise amount acquisition unit 352 (mSAD ⁇ N). Then, it is determined whether or not it is in a stationary state or an operating state (step S3).
  • the determination unit 354 determines that the difference average value mSAD calculated by the evaluation value calculation unit 351 is equal to or less than the noise amount N acquired by the estimated noise amount acquisition unit 352 (step S3: Yes)
  • the image processing apparatus 30 proceeds to step S4 described later.
  • step S3 when the determination unit 354 determines that the difference average value mSAD calculated by the evaluation value calculation unit 351 is not less than or equal to the noise amount acquired by the estimated noise amount acquisition unit 352 (step S3: No), image processing is performed.
  • the apparatus 30 proceeds to step S5 described later.
  • step S4 the noise reduction processing unit 355 performs time direction NR processing on the target pixel in the current image. Details of the time direction NR process will be described later.
  • step S4 the image processing apparatus 30 proceeds to step S6 described later.
  • step S5 the noise reduction processing unit 355 performs spatial direction NR processing on the target pixel in the current image. Details of the spatial direction NR processing will be described later. After step S5, the image processing apparatus 30 proceeds to step S6 described later.
  • step S6: Yes when the image processing control unit 39 determines that NR processing has been performed on all the pixels of the current image (step S6: Yes), the image processing apparatus 30 proceeds to step S7 described later.
  • step S6: No when the image processing control unit 39 determines that NR processing has not been performed on all the pixels of the current image (step S6: No), the image processing device 30 returns to step S1 described above.
  • step S7 the noise reduction processing unit 355 outputs the NR image to the frame memory 33 and the image processing unit 38.
  • the image processing apparatus 30 ends this process.
  • the determination unit 354 determines whether the difference average value mSAD calculated by the evaluation value calculation unit 351 is equal to or less than the noise amount N, so that the noise reduction processing unit 355 performs time direction NR processing or space direction NR processing.
  • the noise reduction processing unit 355 performs both the time direction NR process and the spatial direction NR process, and sets the difference average value mSAD (or the reliability indicating the certainty of the motion vector).
  • the result of each NR process may be blended based on the blend value set accordingly.
  • FIG. 7 is a flowchart showing an overview of the time direction NR process.
  • the noise reduction processing unit 355 acquires RTS noise information from the storage unit 36, and based on the RTS noise information, whether or not there is a possibility that RTS noise is generated in the target pixel of the current image. Is determined (step S10).
  • the noise reduction processing unit 355 determines that there is a possibility that RTS noise is generated in the target pixel of the current image (step S10: Yes)
  • the image processing apparatus 30 proceeds to step S11 described later.
  • the noise reduction processing unit 355 determines that there is no possibility that RTS noise occurs in the target pixel of the current image (step S10: No)
  • the image processing apparatus 30 proceeds to step S12 described later. .
  • step S11 the noise reduction processing unit 355 sets the coefficient Co_we_pre used in the time direction NR process according to the RTS_Value included in the RTS noise information acquired from the storage unit 36.
  • step S11 the image processing apparatus 30 proceeds to step S13 described later.
  • FIG. 8 is a diagram illustrating the relationship between the RTS_Value and the coefficient Co_we_pre set by the noise reduction processing unit 355.
  • the horizontal axis indicates the magnitude of RTS_Value
  • the vertical axis indicates the coefficient Co_we_pre.
  • the noise reduction processing unit 355 sets the coefficient Co_we_pre according to the magnitude of RTS_Value, as indicated by the straight line L1, the curved line L2, and the broken line L3.
  • the noise reduction processing unit 355 may set the coefficient Co_we_pre so as to increase linearly as indicated by the straight line L1 or exponentially increase as indicated by the curve L2 according to the magnitude of the RTS_Value.
  • the coefficient Co_we_pre may be set so as to be, or the coefficient Co_we_pre may be set so as to increase stepwise as indicated by a broken line L3.
  • the noise reduction processing unit 355 increases the coefficient Co_we_pre used in the time direction NR process when there is a possibility that the RTS noise is generated in the target pixel as compared with the case where the RTS noise is not generated in the target pixel.
  • step S13 the image processing apparatus 30 proceeds to step S13 described later.
  • the noise reduction processing unit 355 performs time direction NR processing using the weighting coefficient we_pre calculated in step S13 described above (step S14), and the image processing apparatus 30 returns to the main routine of FIG. Specifically, the noise reduction processing unit 355 performs time direction NR processing using the following equation (3).
  • F G_NR (x, y) is the pixel value of the current image at the coordinates (x, y) of the NR image.
  • we_cur and we_pre are weighting coefficients at the time of the weighted average process.
  • the amount of noise reduction may be increased by increasing we_pre compared to we_cur.
  • we_pre may be set to a fixed value in advance, or a user may set an arbitrary value via the operation unit 37.
  • FIG. 9 is a flowchart showing an outline of the spatial direction NR process.
  • the noise reduction processing unit 355 acquires RTS noise information from the storage unit 36, and based on the RTS noise information, whether or not there is a possibility that RTS noise is generated in the target pixel of the current image. Is determined (step S21).
  • the noise reduction processing unit 355 determines that there is a possibility that RTS noise is generated in the target pixel of the current image (step S21: Yes)
  • the image processing apparatus 30 proceeds to step S22 described later.
  • the noise reduction processing unit 355 determines that there is no possibility that RTS noise is generated in the target pixel of the current image (step S21: No)
  • the image processing apparatus 30 proceeds to step S22 described later. .
  • step S21 the noise reduction processing unit 355 sets the coefficient Co_Cb used in the spatial direction NR processing according to the RTS_Value included in the RTS noise information acquired from the storage unit 36.
  • step S21 the image processing apparatus 30 proceeds to step S23 described later.
  • FIG. 10 is a diagram illustrating a relationship between the RTS_Value and the coefficient Co_Cb set by the noise reduction processing unit 355.
  • the horizontal axis indicates the magnitude of RTS_Value
  • the vertical axis indicates the coefficient Co_Cb.
  • the noise reduction processing unit 355 sets the coefficient Co_Cb according to the magnitude of RTS_Value, as indicated by the straight line L11, the curved line L12, and the broken line L13.
  • the noise reduction processing unit 355 may set the coefficient Co_Cb so as to increase linearly as indicated by a straight line L11 according to the magnitude of RTS_Value, or may increase exponentially as indicated by a curve L12.
  • the coefficient Co_Cb may be set so as to be, or the coefficient Co_Cb may be set so as to increase stepwise as indicated by a broken line L13. In this way, the noise reduction processing unit 355 increases the coefficient Co_Cb used in the spatial direction NR process when there is a possibility that RTS noise is generated in the target pixel, compared to the case where RTS noise is not generated in the target pixel.
  • the image processing apparatus 30 proceeds to step 23 described later.
  • the noise reduction processing unit 355 calculates the weighting coefficient Cb using the coefficient Co_Cb set in Step S21 or Step S22 described above (Step S24). Specifically, the noise reduction processing unit 355 calculates the weighting coefficient Cb by the following equation (4).
  • Cb Co_Cb ⁇ Cb (4)
  • the noise reduction processing unit 355 performs the spatial direction NR process using the weighting coefficient Cb calculated in step S23 described above (step S24), and the image processing apparatus 30 returns to the main routine of FIG. Specifically, the noise reduction processing unit 355 performs the spatial direction NR processing using the following equation (5).
  • we_diff_cur (x + i, y + j) and we_diff_pre (x + i, y + j) correspond to the weighting coefficients in the weighted average process. This coefficient is given by a Gaussian distribution as shown in Equation (6) below. I and j are natural numbers. Further, m and n are SAD (m, n) selected as the difference average value mSAD in the above-described equation (1).
  • the difference between the pixel value of the target pixel of the current image and the pixel values of the peripheral pixels of the target pixel is adaptively set according to the above. Specifically, when the difference between the pixel value of the target pixel of the current image and the pixel values of the peripheral pixels of the target pixel is large, the weight during the weighted average process is small. Therefore, the pixels in the region where the pixel value changes abruptly, such as the edge portion, do not contribute to the weighted average process, and therefore, it is possible to reduce the noise component while retaining the edge portion.
  • the degree of noise reduction (smoothing) Of the image) depends on the amount of noise included in the image. Specifically, since the difference increases as the noise increases, the weighting coefficient decreases, and the contribution in the weighted average of the above equation (3) decreases. Therefore, the greater the noise, the weaker the degree of noise reduction, and the noise cannot be reduced.
  • the noise reduction processing unit 355 calculates the standard deviation ⁇ of the Gaussian distribution of the above equation (6) based on the noise amount N output in step S2 of FIG. Specifically, the noise reduction processing unit 355 calculates the standard deviation ⁇ using the following equation (7).
  • Cb ⁇ N (7)
  • the weight coefficient Cb is a positive real number.
  • the weighting factor Cb may be set to a predetermined value in advance, or may be configured to be set to an arbitrary value by the user via the operation unit 37.
  • the noise reduction processing unit 355 can perform adaptive noise reduction processing on the noise amount by calculating the standard deviation ⁇ of the Gaussian distribution based on the noise amount N. That is, the noise reduction processing unit 355 increases the standard deviation ⁇ as the amount of noise increases, and therefore the difference (for example, F G_cur (x + i, y + j) ⁇ F G_cur (x +, y)) in the above equation (6) becomes noise. Even if it becomes larger, the weight coefficient can be made larger than in the case where the standard deviation ⁇ does not depend on the estimated noise amount N. For this reason, even when noise is large, the strength of smoothing can be maintained.
  • the noise reduction unit 35 can generate RTS noise at the target pixel of the current image based on the RTS noise information input from the third external I / F unit 31. If the pixel value of the target pixel is higher than the coefficient of the NR process performed on the pixel in which RTS noise is not generated, the pixel value of the target pixel Even when the correlation with the pixel value is low, RTS noise that may occur in a moving image can be reduced, so that the image quality can be improved.
  • Embodiment 1 of the present invention when there is a possibility that RTS noise occurs in the target pixel of the current image, the NR function provided in the conventional imaging device or the like for the pixel value of the target pixel Therefore, RTS noise can be reduced with a simple configuration.
  • the SAD is used by the motion detection unit 34 to calculate the motion vector (m, n), but other SSDs (Sum of Squared Difference) or NCC (Normalized Cross-Correlation). It may be calculated by an existing method.
  • the imaging system according to the second embodiment has a configuration different from that of the image processing device 30 according to the first embodiment described above, and also a process executed by the image processing device according to the second embodiment. Specifically, the image processing apparatus according to the second embodiment calculates a representative value for correcting the target pixel based on the representative value calculated from the current image and the representative value calculated from the past image. The pixel value of the target pixel is corrected based on the representative value.
  • processing executed by the image processing apparatus according to the second embodiment will be described.
  • symbol is attached
  • FIG. 11 is a block diagram schematically showing the configuration of the imaging system according to Embodiment 2 of the present invention.
  • An imaging system 1a illustrated in FIG. 11 includes an image processing device 30a instead of the image processing device 30 of the imaging system 1 according to Embodiment 1 described above.
  • the image processing device 30a omits the configurations of the interpolation processing unit 32, the frame memory 33, and the motion detection unit 34 from the image processing device 30 according to the first embodiment described above, and replaces the noise reduction unit 35 with a noise reduction unit 35a. Is provided.
  • the noise reduction unit 35 a corrects the RTS noise for the current image input from the third external I / F unit 31 and outputs the corrected image to the image processing unit 38.
  • the noise reduction unit 35a includes an RTS noise correction unit 320.
  • the RTS noise correction unit 320 performs RTS noise correction processing for correcting RTS noise on the current image input from the third external I / F unit 31, and the current image on which the RTS noise correction has been performed is performed on the image processing unit. 38.
  • FIG. 12 is a block diagram illustrating a detailed configuration of the RTS noise correction unit 320.
  • the RTS noise correction unit 320 includes an RTS noise pixel determination unit 321, a candidate value calculation unit 322, a representative value calculation unit 323, a random noise amount estimation unit 324, and a correction value calculation unit 325. Have.
  • the RTS noise pixel determination unit 321 Based on the RTS noise information input from the third external I / F unit 31, the RTS noise pixel determination unit 321 generates RTS noise in the target pixel of the current image input from the third external I / F unit 31. It is determined whether or not there is a possibility, and the determination result is output to the candidate value calculation unit 322 and the representative value calculation unit 323. Specifically, when a position of a pixel is input to the RTS noise pixel determination unit 321, it is determined whether there is RTS noise information corresponding to the pixel (target pixel), and there is RTS noise information. While this RTS noise information (information indicating the presence of RTS noise) is output, if there is no RTS noise information, it is regarded as a pixel in which RTS noise does not occur, and RTS noise information is not output.
  • the candidate value calculation unit 322 is determined by the RTS noise pixel determination unit 321 that RTS noise may occur in the target pixel.
  • a plurality of candidate values for the correction amount for the pixel value of the target pixel, and the representative value calculation unit 323, the random noise amount estimation unit 324, and the correction value calculation for the pixel value of the target pixel and the plurality of calculated candidate values are determined by the RTS noise pixel determination unit 321 that RTS noise may occur in the target pixel.
  • a plurality of candidate values for the correction amount for the pixel value of the target pixel, and the representative value calculation unit 323, the random noise amount estimation unit 324, and the correction value calculation for the pixel value of the target pixel and the plurality of calculated candidate values To each of the units 325.
  • the representative value calculation unit 323 determines that the RTS noise pixel determination unit 321 around the target pixel has at least the RTS noise. Based on a pixel that is determined not to occur and a random noise amount corresponding to the target pixel calculated by a random noise amount estimation unit 324 described later, a representative value corresponding to a pixel value when no RTS noise occurs is obtained. calculate.
  • the representative value calculation unit 323 calculates a representative value of the current image based on the current image input from the third external I / F unit 31 and the past image stored in the storage unit 36, and calculates a correction value. Output to the unit 325.
  • FIG. 13 is a block diagram illustrating a detailed configuration of the representative value calculation unit 323.
  • the representative value calculation unit 323 includes a first reference value calculation unit 323a, a first random noise amount estimation unit 323b, an allowable range calculation unit 323c, and a representative value determination unit 323d.
  • the first reference value calculation unit 323a calculates a reference value used when estimating the random noise amount. Specifically, the first reference value calculation unit 323a uses, as a reference value, a value obtained by adding RTS_Value to the pixel value of the target pixel, or a value obtained by adding the maximum value among the correction value candidate values to the pixel value of the target pixel. calculate.
  • the first random noise amount estimation unit 323b calculates a random noise amount based on the random noise model input from the third external I / F unit 31 and the reference value calculated by the first reference value calculation unit 323a.
  • the allowable range calculation unit 323c can be used for the representative value calculation process based on the pixel value of the pixel in the calculation range set with the target pixel as a reference and the random noise amount calculated by the first random noise amount estimation unit 323b.
  • An allowable range (effective range) which is a range of pixel values is calculated.
  • the representative value determination unit 323d determines the pixel value included in the allowable range (within the effective range) for the pixel values other than the target pixel in the calculation range (in the case of the image sensor 105 using the color filter, the target pixel). If the selected pixel value is equal to or greater than a predetermined value, the median value of the selected pixel values is determined as the representative value. Further, the representative value determining unit 323d determines the representative value of the current image based on the current image and the past image stored in the storage unit 36.
  • the random noise amount estimation unit 324 estimates the random noise amount corresponding to the pixel value based on the random noise model input from the third external I / F unit 31, and the estimation result is used as the candidate value calculation unit 322, the representative value.
  • the data are output to the calculation unit 323 and the correction value calculation unit 325, respectively. That is, when a pixel value is input to the random noise amount estimation unit 324, a random noise amount corresponding to the pixel value is output.
  • the correction value calculation unit 325 is based on the plurality of candidate values calculated by the candidate value calculation unit 322.
  • the pixel value of the target pixel is corrected.
  • the correction value calculation unit 325 is based on the pixel value of the target pixel, the plurality of candidate values calculated by the candidate value calculation unit 322, and the representative value calculated by the representative value calculation unit 323.
  • the pixel value corrected for the RTS noise is calculated and output to the image processing unit 38. More specifically, the correction value calculation unit 325 selects a candidate value whose correction result is closest to the representative value calculated by the representative value calculation unit 323 from among a plurality of candidate values calculated by the candidate value calculation unit 322.
  • the correction value calculation unit 325 outputs the pixel value of the target pixel as it is.
  • the correction value calculation unit 325 functions as a correction unit.
  • FIG. 14 is a block diagram illustrating a detailed configuration of the correction value calculation unit 325.
  • the correction value calculation unit 325 includes a second reference value calculation unit 325a, a second random noise amount estimation unit 325b, a correction amount determination unit 325c, and a pixel value correction unit 325d.
  • the second reference value calculation unit 325a calculates, as a reference value, a value obtained by subtracting RTS_Value from the pixel value of the target pixel or a value obtained by subtracting the maximum value among the correction value candidate values from the pixel value of the target pixel.
  • the second random noise amount estimation unit 325b calculates a random noise amount based on the random noise model input from the third external I / F unit 31 and the reference value calculated by the second reference value calculation unit 325a.
  • the correction amount determination unit 325c calculates a value obtained by multiplying the random noise amount calculated by the second random noise amount estimation unit 325b by a certain coefficient as a threshold value, and the correction amount candidate value calculated by the candidate value calculation unit 322 After determining whether or not the maximum candidate value is equal to or larger than the threshold value, the representative value calculated by the representative value calculation unit 323 and the candidate value of the correction amount closest to the absolute value of the target pixel are selected as the correction amount. .
  • the pixel value correction unit 325d calculates the corrected pixel value of the target pixel by adding or subtracting the correction amount to the pixel value of the target pixel so as to approach the representative value.
  • FIG. 15 is a flowchart showing an outline of processing executed by the image processing apparatus 30a, and is a flowchart of a main routine executed by the image processing apparatus 30a.
  • the RTS noise correction unit 320 sets a target pixel for sequentially performing the processes of steps S102 to S105 described later (step S101). Note that the RTS noise correction unit 320 assigns integers larger than 0 as indexes in order of rasters from the upper left to the lower right for each pixel in raster order. Next, each time step S101 is executed, the RTS noise correction unit 320 increments the counter by 1 (the counter is reset to 0 when the processing of FIG. 15 is started). The RTS noise correction unit 320 sets the pixel to which the index indicated by the counter is assigned as the target pixel.
  • step S101 when step S101 is first executed by the RTS noise correction unit 320, the counter is incremented by 1 because the RTS noise correction unit 320 increments the counter by 1, so that the upper left pixel becomes the target pixel. .
  • the RTS noise correction unit 320 executes the process of step S101 twice (second time), the counter indicates 2, so the pixel on the right side of the upper left pixel becomes the target pixel.
  • the RTS noise pixel determination unit 321 includes the RTS noise information recording unit 112b of the nonvolatile memory 112 of the imaging device 10 via the third external I / F unit 31, the first external I / F unit 115, and the bus 113.
  • the RTS noise information to be recorded is acquired, and based on the acquired RTS noise information, it is determined whether or not there is a possibility that RTS noise is generated in the pixel of interest (step S102). That is, the RTS noise pixel determination unit 321 determines whether or not the position information of the shared pixel block including the target pixel is included in the RTS noise information.
  • the RTS noise pixel determination unit 321 determines whether or not the position information of the shared pixel block including the target pixel is included in the RTS noise information as a shared pixel block in which RTS noise may occur. judge. When the RTS noise pixel determination unit 321 determines that RTS noise may occur in the target pixel (determined that the position information of the shared pixel block including the target pixel is included in the RTS noise information) ( In step S102: Yes), the image processing apparatus 30a proceeds to step S103 described later. On the other hand, the RTS noise pixel determination unit 321 determines that there is no possibility of occurrence of RTS noise in the target pixel (determined that the position information of the shared pixel block including the target pixel is not included in the RTS noise information.
  • Step S102 No
  • the image processing device 30a proceeds to step S106 described later.
  • the pixel value of the target pixel is directly output to the representative value calculation unit 323 as a corrected pixel value. To do.
  • the candidate value calculation unit 322 calculates a plurality of correction amount candidate values for correcting the RTS noise. Specifically, the candidate value calculation unit 322 calculates a pixel value of 0 or more and RTS_Value or less based on RTS_Value (included in the RTS noise information output from the RTS noise pixel determination unit 321) corresponding to the target pixel. All possible values (all integers from 0 to RTS_Value when only an integer can be taken as a RAW image) are set as candidate values.
  • the representative value calculation unit 323 determines a representative value (in the case where no RTS noise is generated in the target pixel) based on the pixel values of the peripheral pixels of the target pixel of the past image and the current image stored in the storage unit 36.
  • a representative value calculation process for calculating (predicted pixel value) is executed (step S104). Note that the representative value calculation unit 323 may include pixels that generate RTS noise.
  • FIG. 16 is a flowchart showing an overview of the representative value calculation processing in step S104 of FIG.
  • the representative value calculation unit 323 sets a minimum calculation range that is a target for calculating a representative value with reference to the target pixel (step S201). Specifically, the representative value calculation unit 323 minimizes 3 ⁇ 3, which is the minimum range of 7 ⁇ 7 or less, for example, when the target range is a maximum of 7 ⁇ 7 with the target pixel as the center. Set as the calculation range.
  • the first reference value calculation unit 323a calculates, as a reference value, a value obtained by adding RTS_Value to the pixel value of the pixel of interest or a value obtained by adding the maximum correction amount candidate value to the pixel value of the pixel of interest (step) S202).
  • the first random noise amount estimation unit 323b performs random noise based on the random noise model input from the third external I / F unit 31 and the reference value calculated by the first reference value calculation unit 323a in step S202.
  • the amount is calculated (step S203).
  • the first random noise amount estimation unit 323b according to the second embodiment is based on the noise model (random noise model) of the curve L X1 in FIG. 6 and the reference value calculated by the first reference value calculation unit 323a.
  • the amount of random noise is calculated (standard deviation is calculated).
  • the random noise model may be approximated by an approximate expression or a broken line.
  • the allowable range calculation unit 323c calculates an allowable range (effective range) that is a range of pixel values that can be used for the representative value calculation processing based on the pixel values of the pixels within the calculation range (step S203). S204). Specifically, the allowable range calculation unit 323c calculates the upper limit of the allowable range (effective range) by the following equation (8). Reference value + random noise amount (standard deviation) x R + RTS_Value ... (8)
  • R is a predetermined coefficient, and is set according to how much the RTS noise can be visually grasped with respect to the random noise. For example, a value of around 2 is preferable as the coefficient of R. Further, the representative value calculation unit 323 calculates the lower limit of the allowable range by the following equation (9). Reference value-Random noise amount (standard deviation) x R-RTS_Value ... (9)
  • the maximum value of a plurality of candidate values may be used instead of RTS_Value.
  • the reference values in the equations (8) and (9) are obtained by a reference value method different from the reference value used for estimating the random noise amount by the first random noise amount estimation unit 323b in step S203 described above. It is also possible to use the reference value.
  • the allowable range calculation unit 323c can calculate the allowable range in consideration of the RTS noise of the target pixel and the random noise around the target pixel.
  • the allowable range calculation unit 323c performs the above-described steps for each of the pixel values of pixels other than the target pixel in the calculation range (in the case of the image sensor 105 using a color filter, the pixel value of the same color as the target pixel). It is determined whether it is within the allowable range calculated in S204, and the number of pixel values within this allowable range is counted (step S205). The count value obtained in step S205 tends to be larger for a flat subject and smaller for a subject including an edge. Note that pixels that may cause RTS noise within the calculation range may not be counted.
  • step S206: Yes when the count value (number of pixels) counted in step S205 described above is larger than the predetermined value ThRef (step S206: Yes), the image processing device 30a proceeds to step S209 described later.
  • the predetermined value ThRef is preferably 1 or more because the representative value calculation unit 323 calculates a representative value from the peripheral pixels of the target pixel.
  • step S206: No when the count value counted in step S205 described above is not larger than the predetermined value ThRef (step S206: No), the image processing device 30a proceeds to step S207 described later.
  • step S207 when the calculation range that is the target of the representative value calculation is the maximum (step S207: Yes), the image processing apparatus 30a proceeds to step S209 described later. On the other hand, when the calculation range that is the target of the representative value calculation is not the maximum (step S207: No), the image processing device 30a proceeds to step S208 described later.
  • step S208 the allowable range calculation unit 323c expands the calculation range for calculating the representative value (step S208). Specifically, the allowable range calculation unit 323c enlarges the calculation range that is the target of calculating the representative value within the maximum range by one or more pixels in the horizontal or vertical direction. For example, when the 3 ⁇ 3 range centered on the target pixel is set as the calculation range, the allowable range calculation unit 323c resets the 5 ⁇ 5 range centered on the target pixel as the calculation range.
  • the image processing apparatus 30a returns to step S202.
  • the allowable range calculation unit 323c has set the 3 ⁇ 3 or 5 ⁇ 5 range as the calculation range. For example, the allowable range calculation unit 323c enlarges only the horizontal or vertical range to 5 ⁇ 3 or 3 ⁇ 5. The range may be set as the calculation range.
  • the representative value determining unit 323d calculates a representative value of the current image. Specifically, the representative value determination unit 323d first uses a pixel value (using a color filter) included in the allowable range (effective range) for the pixel value of the pixel other than the target pixel in the calculation range. In the case of the image sensor 105, the pixel value of the same color as the target pixel) is selected. Thereafter, when the number of selected pixels is equal to or greater than the predetermined value ThRef, the representative value determining unit 323d calculates the median value of the selected pixel values as the representative value of the current image.
  • a pixel value using a color filter included in the allowable range (effective range) for the pixel value of the pixel other than the target pixel in the calculation range. In the case of the image sensor 105, the pixel value of the same color as the target pixel) is selected. Thereafter, when the number of selected pixels is equal to or greater than the predetermined value ThRef, the representative value determining unit 323
  • the representative value determining unit 323d calculates the median value closer to the pixel value of the target pixel as the representative value of the current image. In this case, overcorrection can be prevented.
  • the representative value determining unit 323d calculates the pixel values of pixels other than the target pixel in the calculation range having the pixel value closest to the pixel value of the target pixel of the current image. Calculate as a representative value.
  • the representative value determining unit 323d calculates the representative value using the median value, but may be calculated by another method such as an average or an intermediate value of distribution.
  • the representative value determining unit 323d may determine the edge direction within the calculation range, and calculate the peripheral pixel value in the direction with the highest correlation as the representative value of the current image based on the result of the edge direction determination. . Further, the representative value determining unit 323d may exclude pixels that may generate RTS noise in pixels within the calculation range other than the target pixel. At this time, the representative value determining unit 323d determines the pixel value in the current image of the pixel of interest when there is no pixel that is unlikely to cause RTS noise within the calculation range at the time when step S209 is executed. The representative value of the current image.
  • step S210 if the current image is the first image data for which a representative value is calculated (step S210: Yes), the representative value determining unit 323d determines the representative value of the current image calculated in step S209 described above as a representative value. (Step S211). After step S211, the image processing apparatus 30a returns to the main routine of FIG. 15 described above.
  • step S210 when the current image is not the first image data for calculating the representative value (step S210: No), the representative value determining unit 323d determines the representative value of the current image and the representative value of the past image stored in the storage unit 36. Is determined as a representative value (step S212). In this case, the representative value determination unit 323d may calculate an average value using the representative values of the plurality of past images in the storage unit 36. The representative value determining unit 323d determines the average value of the representative value of the current image and the representative value of the past image as the representative value, but may determine the weighted average or the median value as the representative value.
  • the representative value determining unit 323d is based on the representative values of the plurality of past images (for example, three or more) and the representative value of the current image. A representative value may be determined. Furthermore, the representative value determining unit 323d sets a predetermined range based on the representative value of the past image stored in the storage unit 36 in order to prevent fluctuation within a certain range of the previous result. The representative value of the current image may be determined so that the representative value of the current image is included therein. After step S212, the image processing device 30a returns to the main routine of FIG. 15 described above.
  • the representative value calculation unit 323 calculates the representative value of the current image in the above-described representative value calculation process, and then the current image is the first image data (the past image is stored in the storage unit 36). If the current image is not the first image data (when the past image is stored in the storage unit 36), the representative value of the current image and the past image are calculated. The average value of the representative values is calculated as a representative value.
  • the representative value calculation unit 323 stores the current image and an uncorrected image (hereinafter referred to as “previous image before correction”) temporally continuous with the current image in the storage unit 36, and the corrected past image.
  • the representative value may be calculated from the representative value of the current image and the representative value of the current image.
  • step S105 the correction value calculation unit 325 is based on the plurality of candidate values calculated by the candidate value calculation unit 322 in step S103 described above and the representative value calculated by the representative value calculation unit 323 in step S104 described above. Then, a correction value calculation process for calculating a pixel value in which the RTS noise in the target pixel is corrected is executed.
  • step S105 the image processing apparatus 30a proceeds to step S106 described later.
  • FIG. 17 is a flowchart showing an outline of the correction value calculation processing in step S105 of FIG.
  • the second reference value calculation unit 325a calculates, as a reference value, a value obtained by subtracting RTS_Value (maximum correction amount candidate value) from the pixel value of the target pixel (step S301).
  • the second random noise amount estimation unit 325b receives the random noise recorded in the random noise model information recording unit 112c via the third external I / F unit 31, the first external I / F unit 115, and the bus 113.
  • a model is acquired, and a random noise amount is calculated based on the acquired random noise model and the reference value calculated by the second reference value calculation unit 325a in step S301 (step S302).
  • the correction amount determination unit 325c calculates, as a threshold value, a value obtained by multiplying the random noise amount calculated by the second random noise amount estimation unit 325b in step S302 described above by a certain coefficient (step S303).
  • the threshold value is calculated by the following equation (10). Random noise amount x Rm (10) Rm is determined according to how much RTS noise is visually seen with respect to random noise. For example, the value of Rm is preferably around 2.
  • step S304 the maximum candidate value among the correction amount candidate values calculated by the candidate value calculation unit 322 in step S103 of FIG. 15 described above is equal to or greater than the threshold calculated in step S303 described above. Whether or not (step S304).
  • step S304 determines that the maximum candidate value among the correction amount candidate values calculated by the candidate value calculation unit 322 in step S103 of FIG. 15 described above is greater than or equal to the threshold value (step S304: Yes).
  • step S304 Yes
  • the image processing device 30a proceeds to step S305 described later.
  • step S304 determines that the maximum candidate value among the correction amount candidate values calculated by the candidate value calculation unit 322 in step S103 of FIG. 15 described above is not greater than or equal to the threshold (step S304). : No)
  • the image processing apparatus 30a proceeds to Step S306 to be described later.
  • step S305 the correction amount determination unit 325c selects a correction amount candidate value closest to the absolute value of the difference between the representative value calculated by the representative value calculation unit 323 in step S104 of FIG. 15 and the pixel value of the target pixel. Select as the correction amount.
  • the image processing apparatus 30a proceeds to step S307 described later.
  • step S306 the correction amount determination unit 325c sets 0 as the correction amount. After step S306, the image processing apparatus 30a proceeds to step S307.
  • the pixel value correction unit 325d corrects the corrected pixel value of the target pixel by adding or subtracting a correction amount to the pixel value of the target pixel so as to approach the representative value (step S307).
  • the image processing apparatus 30a returns to the main routine of FIG.
  • step S106 the RTS noise correction unit 320 determines whether or not the processing in steps S101 to S105 described above has been completed for all pixels.
  • step S106: Yes the image processing device 30a ends this processing.
  • step S106: No the image processing device 30a performs the above steps. Return to S101.
  • the representative value calculation unit 323 sets the representative value of the current image input from the third external I / F unit 31 and the representative value of the past image stored in the storage unit 36. Based on the representative value calculated by the representative value calculation unit 323 by the correction value calculation unit 325, a representative value corresponding to the pixel value when no RTS noise is generated is calculated. Therefore, even if the pixel value of the target pixel is constant and the pixel values of the peripheral pixels of the target pixel change in magnitude for each current image, the representative value including the representative value of the past image is calculated. Thus, fluctuations in the pixel value of the target pixel after correction can be reduced. As a result, even when the correlation between the pixel value of the target pixel and the pixel values of the peripheral pixels is low, the image quality can be improved.
  • the third embodiment corresponds to a pixel value when no RTS noise is generated after the pixel of interest of the current image and the pixel of the past image corresponding to the pixel of interest of the current image are aligned. A representative value is calculated.
  • the representative value calculation process executed by the image processing apparatus according to the third embodiment will be described.
  • symbol is attached
  • FIG. 18 is a block diagram showing a detailed configuration of the noise reduction unit according to the third embodiment.
  • the noise reduction unit 35b illustrated in FIG. 18 includes the motion detection unit 34 according to the first embodiment and the RTS noise correction unit 320b.
  • the RTS noise correction unit 320b performs an RTS noise correction process for correcting RTS noise on the current image (RAW) input from the third external I / F unit 31, and the corrected RAW image is converted into an image processing unit. 38.
  • FIG. 19 is a block diagram illustrating a detailed configuration of the RTS noise correction unit 320b.
  • the RTS noise correction unit 320 b includes a representative value calculation unit 326 instead of the representative value calculation unit 323 according to the second embodiment described above.
  • the representative value calculation unit 326 includes the motion vector detected by the motion detection unit 34, the pixel value of the target pixel of the current image (RAW image) input from the third external I / F unit 31, and the current image stored in the storage unit 36. Based on the pixel value of the target pixel of the past image corresponding to the target pixel of the (RAW image), a representative value corresponding to the pixel value when no RTS noise occurs is calculated. Specifically, the representative value calculation unit 326 uses the motion vector detected by the motion detection unit 34 to align the target pixel of the current image with the pixel of the past image corresponding to the target pixel of the current image. Based on the pixel value of the target pixel of the current image and the pixel value of the target pixel of the past image, a representative value corresponding to the pixel value when no RTS noise is generated is calculated.
  • FIG. 20 is a flowchart showing an outline of representative value calculation processing executed by the image processing apparatus 30a according to the third embodiment.
  • step S311: Yes when the current image input from the third external I / F unit 31 is the first image data for calculating the representative value (step S311: Yes), the image processing device 30a proceeds to step S312 described later. To do.
  • step S311: No when the current image input from the third external I / F unit 31 is not the first image data for calculating the representative value (step S311: No), the image processing device 30a proceeds to step S321 described later. To do.
  • Step S312 to Step S320 correspond to Step S201 to Step S209 of FIG. 16 described above, respectively.
  • the image processing device 30a returns to the main routine of FIG. 15 described above.
  • step S321 the representative value calculation unit 326 determines a representative value corresponding to a pixel value when no RTS noise occurs based on the past image and the current image stored in the storage unit 36. Specifically, the representative value calculation unit 326 aligns the target pixel of the current image and the pixel of the past image corresponding to the position of the target pixel of the current image based on the motion vector detected by the motion detection unit 34. Thereafter, the pixel value of the past image corresponding to the position of the target image of the current image is determined as the representative value.
  • the representative value calculation unit 326 desirably uses a local vector as a motion vector.
  • the representative value calculation unit 326 may use a global vector. In this case, it is possible to perform alignment with respect to movement of the entire image such as camera shake.
  • the image processing apparatus 30a returns to the main routine of FIG.
  • the representative value calculation unit 326 corresponds to the position of the target pixel of the current image and the target pixel of the current image based on the motion vector detected by the motion detection unit 34. Since the pixel value of the past image corresponding to the position of the target image of the current image is calculated as the representative value after the pixel alignment of the past image is performed, the peripheral pixels of the target pixel of the current image may change in size. However, since this change in size can be reduced, RTS noise can be reduced without degrading the image quality.
  • Embodiment 3 of the present invention even when only the pixel value of the target pixel changes in size, RTS noise can be reduced without degrading the image quality.
  • the noise reduction unit according to the modification of the third embodiment includes a reliability calculation unit that calculates the reliability of the motion vector.
  • the pixel value of the pixel of interest of the current image and the pixel value of the pixel of the past image corresponding to the pixel of interest of the current image are used as pixel values when no RTS noise occurs. The corresponding representative value is calculated.
  • FIG. 21 is a block diagram showing a detailed configuration of the noise reduction unit according to the modification of the third embodiment of the present invention.
  • the noise reduction unit 35b1 illustrated in FIG. 21 further includes a reliability calculation unit 34a in addition to the configuration of the noise reduction unit 35b according to Embodiment 3 described above.
  • the reliability calculation unit 34a calculates a reliability representing the certainty of the amount of movement of the subject calculated by the motion detection unit 34.
  • FIG. 22 is a flowchart showing an overview of representative value calculation processing according to a modification of the third embodiment of the present invention.
  • steps S401 to S410 correspond to the above-described steps S311 to S320 of FIG.
  • steps S411 to S418 correspond to steps S312 to S319 of FIG. 20 described above, respectively.
  • step S419 the representative value calculation unit 326 determines a representative value Repi. Specifically, the representative value calculation unit 326 determines the representative value Repi by performing the same processing as in step S410. After step S419, the image processing apparatus 30a proceeds to step S421 described later.
  • step S420 the representative value calculation unit 326 calculates the pixel value of the past image as the representative value Repr. After step S420, the image processing apparatus 30a proceeds to step S421 described later.
  • the representative value calculation unit 326 calculates the similarity between the representative value Repi of the current image and the representative value Repr of the past image (step S421). Specifically, the representative value calculation unit 326 calculates a difference (
  • the representative value calculation unit 326 calculates the magnitude of the motion vector input from the motion detection unit 34, the reliability calculated by the reliability calculation unit 34a, the representative value Repi of the current image calculated in step S421, and the representative value of the past image.
  • the similarity is calculated based on the difference of Repr (
  • FIG. 23 is a diagram schematically illustrating a method of calculating the blend ratio b calculated by the representative value calculation unit 326.
  • the horizontal axis indicates the similarity and the vertical axis indicates the blend ratio b.
  • each of the straight line LB1, the curved line LB2, and the broken line LB3 shows an example of the blend ratio.
  • the representative value calculation unit 326 indicates that the state of the current image is a case where the pixel value of the same pixel changes suddenly when the degree of similarity increases, and thus the representative value calculation unit 326 The pixel value is weighted.
  • the representative value calculation unit 326 multiplies the similarity between the representative value Repi of the current image and the representative value Repr of the past image by the reliability Rel of the motion vector calculated by the reliability calculation unit 34a (
  • the reliability of the motion vector calculated by the reliability calculation unit 34a is high, and the past image is likely to be the pixel value of the normal image. (Representative value) is weighted. Note that it may be assumed that the target pixel of the past image is the pixel value of the normal image.
  • the reliability Rel is a value representing the probability of each motion vector based on images that have been aligned between images. The reliability Rel is calculated by the difference average value mSAD of the first embodiment described above, the known SSD, MCC, ZMCC, or the like.
  • step S423 the representative value calculation unit 326 determines a representative value Rep.
  • the image processing device 30a returns to the main routine of FIG. Specifically, the representative value determining unit 323d determines the representative value Rep by the following equation (11).
  • Rep Repi ⁇ b + Repr ⁇ (1-b) (11)
  • the representative value calculation unit 326 may determine the representative value Rep from the following equation (12) using the pixel value Vali of the current image and the pixel value Valr of the past image.
  • Rep Vali ⁇ b + Valr ⁇ (1 ⁇ b) (12)
  • the representative value calculation unit 326 may store the current image and the past image before correction in the storage unit 36, and may use the past image before correction instead of the past image.
  • FIG. 24 is a block diagram illustrating a detailed configuration of the noise reduction unit according to the fourth embodiment.
  • the noise reduction unit 35c illustrated in FIG. 24 includes the motion detection unit 34 according to the first embodiment described above and the RTS noise correction unit 320c.
  • the RTS noise correction unit 320 c performs RTS noise correction processing for correcting RTS noise on the current image input from the third external I / F unit 31, and outputs the corrected RAW image to the image processing unit 38. To do.
  • FIG. 24 is a block diagram illustrating a detailed configuration of the RTS noise correction unit 320c.
  • the RTS noise correction unit 320c includes a representative value calculation unit 327 instead of the representative value calculation unit 323 according to the second embodiment described above.
  • the representative value calculation unit 327 focuses attention when RTS noise does not occur based on the motion vector detected by the motion detection unit 34, the current image, and a plurality of past images temporally continuous with the current image stored in the storage unit 36. A representative value corresponding to the pixel value of the pixel is calculated.
  • FIG. 26 is a flowchart showing an outline of representative value calculation processing executed by the image processing device 30a according to the fourth embodiment.
  • the representative value calculation unit 327 calculates the average value Ave of the pixel value of the target pixel in the current image and the pixel value of each target pixel of the plurality of past images (step S500). In addition to the average value, the representative value calculation unit 327 may calculate using a statistical value such as the pixel value of the target pixel in the current image and the median value of the pixel values of each target pixel in the plurality of past images. .
  • the representative value calculation unit 327 determines the average value Ave calculated in step S500 described above as a representative value (step S501). After step S501, the image processing apparatus 30a returns to the main routine of FIG. 15 described above.
  • the representative value calculation unit 323 uses the current image input from the third external I / F unit 31 and a plurality of past images stored in the storage unit 36, Since the representative value corresponding to the pixel value when no RTS noise occurs is calculated, even if the peripheral pixel of the target pixel of the current image changes in size, this change in size can be reduced. RTS noise can be reduced without degrading the image quality.
  • the representative value calculation unit 323 calculates a representative value corresponding to a pixel value when no RTS noise occurs using the current image and a plurality of past images.
  • the representative value may be calculated by using the image (hereinafter “future image”) corresponding to the shooting time (generation timing) of the current image and the past or future (new) image data.
  • future image the image corresponding to the shooting time (generation timing) of the current image
  • n of these m images may be used as past images
  • mn images may be used as future images.
  • this future image there is a possibility that a frame delay may occur. Therefore, it is preferable to use the future image within a range that allows the frame delay.
  • the representative value calculation unit 323 stores an uncorrected image that is temporally continuous with the current image in the storage unit 36, and uses the uncorrected past image instead of the past image. Also good.
  • Modification 1 of Embodiment 4 Modification 1 of Embodiment 4 of the present invention will be described.
  • the configuration of the noise reduction unit is different, and the representative value calculation process executed by the image processing apparatus is different.
  • the image processing apparatus according to the first modification of the fourth embodiment calculates the representative value using the current pixel and the past image before correction.
  • a representative value calculation process executed by the image processing apparatus according to the first modification of the fourth embodiment will be described.
  • symbol is attached
  • FIG. 27 is a block diagram showing a detailed configuration of the noise reduction unit according to the first modification of the fourth embodiment of the present invention.
  • the noise reduction unit 35c1 illustrated in FIG. 27 further includes a reliability calculation unit 34a in addition to the configuration of the noise reduction unit 35c according to Embodiment 4 described above.
  • FIG. 28 is a flowchart illustrating an overview of representative value calculation processing executed by the image processing device 30a according to the first modification of the fourth embodiment.
  • Co is a predetermined coefficient, for example, a value of 1 or more is preferable. Note that Co may be set in advance to a constant value, or the user may set an arbitrary value.
  • the representative value calculation unit 327 determines that the difference between the pixel value Val of the target pixel of the current image and the pixel value Val (i) of the target pixel of the past image before the i-th correction is equal to or smaller than the allowable range Th (
  • the representative value calculation unit 327 determines that the difference between the pixel value Val of the target pixel of the current image and the pixel value Val (i) of the target pixel of the past image before the i-th correction is less than or equal to the allowable range Th (step) (S514: Yes)
  • the image processing apparatus 30a proceeds to Step S515 described later.
  • the representative value calculation unit 327 determines that the difference between the pixel value Val of the target pixel of the current image and the pixel value Val (i) of the target pixel of the past image before the i-th correction is not less than the allowable range Th. If it has been performed (step S514: No), the image processing apparatus 30a proceeds to step S517 described later.
  • step S515 the representative value calculation unit 327 determines that the reliability Rel (i) of the motion vector of the i-th previous image calculated by the reliability calculation unit 34a is equal to or greater than a predetermined threshold Rel_Th (Rel (i ⁇ ( Rel_Th))
  • the threshold value Rel_Th may be set to a fixed value in advance, or may be set to an arbitrary value by the user.
  • the image processing device 30a determines that the reliability Rel (i) of the motion vector of the i-th previous image calculated by the reliability calculation unit 34a is greater than or equal to a predetermined threshold Rel_Th (step S515: Yes)
  • the image processing device 30a proceeds to step S516, which will be described later. If the reliability Rel of the motion vector of the image (i) is determined not greater than a predetermined threshold value Rel_Th (step S515: No), the image processing apparatus 30a proceeds to step S517 to be described later.
  • step S516 the representative value calculation unit 327 stores the pixel value Val (i) of the target pixel of the i-th past image before correction in mVal [] of the storage unit 36.
  • the representative value calculation unit 327 determines whether the counter i is smaller than the number m of past images before correction input from the storage unit 36 (i ⁇ m) (step S517).
  • the image processing device 30a performs the step The process proceeds to S519.
  • step S519 the representative value calculation unit 327 stores the pixel value Val of the target pixel of the current image in mVal [] of the storage unit 36.
  • the representative value calculation unit 327 calculates an average value Ave of all the pixel values stored in mVal [] of the storage unit 36 (Step S520).
  • the representative value calculation unit 327 may calculate an average value by weighting the pixel value of the target pixel of each image.
  • the weighting coefficient may be increased when the reliability of the motion vector is large, or the weighting coefficient may be decreased when
  • the weight may be calculated from a combination of a plurality of parameters (for example, reliability and difference of motion vectors).
  • the representative value calculation unit 327 may use other statistical values such as a median value and a mode value in addition to the average value and the average value obtained by weighting.
  • the representative value calculation unit 327 determines the average value Ave calculated in step S520 as a representative value (step S521). After step S521, the image processing apparatus 30a returns to the main routine of FIG. *
  • the representative value calculation unit 327 performs a plurality of pre-correction corresponding to the position of the target pixel in the current image based on the pixel value of the target pixel in the current image. A value extracted and calculated from the pixel value of each pixel of the past image is determined as a representative value. Thereby, even if the luminance value of the subject changes, the da representative value is calculated only from the pixel value similar to the pixel value of the current image, so that the reliability of the representative value can be improved.
  • FIG. 29 is a flowchart showing an outline of representative value calculation processing executed by the image processing device 30a according to the second modification of the fourth embodiment.
  • step S601 to step S607 respectively correspond to step S511 to step S518 of FIG. 26 described above.
  • step S608 the representative value calculation unit 327 determines whether or not the number of pixels stored in mVal [] of the storage unit 36 is larger than the threshold value ThRef (mVal []> ThRef).
  • the image processing device 30a proceeds to step S611 described later.
  • the image processing device 30a performs steps described later. The process proceeds to S609.
  • step S609 the representative value calculation unit 327 determines that the difference between the pixel value Val of the target pixel of the current image and the pixel value Val (i) of the target pixel of the past image is equal to or smaller than the allowable range Th (
  • step S609: Yes the image The processing device 30a proceeds to step S610 described later.
  • step S609: No when the representative value calculation unit 327 determines that the difference between the pixel value Val of the target pixel of the current image and the pixel value Val (i) of the target pixel of the past image is not less than the allowable range Th (step S609: No), the image processing apparatus 30a proceeds to Step S611 described later.
  • step S610 the representative value calculation unit 327 stores the pixel value Val (i) of the target pixel of the past image in mVal [] of the storage unit 36.
  • step S610 the image processing apparatus 30a proceeds to step S611.
  • Steps S611 to S613 correspond to Steps S519 to S521 of FIG. 26 described above, respectively.
  • the image processing apparatus 30a returns to the main routine of FIG.
  • the representative value calculation unit 327 when the representative value calculation unit 327 has few similar pixel values from the target pixel in the past image before correction, the pixel value of the past pixel is referred to. Since the representative value is calculated, the reliability of the representative value can be improved.
  • Modification 3 of Embodiment 4 Modification 3 of Embodiment 4 of the present invention will be described.
  • the third modification of the fourth embodiment is different in the representative value calculation process executed by the image processing device 30a according to the fourth embodiment described above. Specifically, the image processing device 30a according to the third modification of the fourth embodiment calculates the representative value using the pixel values of the peripheral pixels of the past image.
  • a representative value calculation process executed by the image processing device 30a according to the third modification of the fourth embodiment will be described.
  • FIG. 30 is a flowchart showing an outline of representative value calculation processing executed by the image processing device 30a according to the third modification of the fourth embodiment.
  • steps S701 to S704 respectively correspond to steps S201 to S204 of FIG.
  • step S707 when the number of pixels of mVal [] stored in the storage unit 36 is greater than a preset threshold ThRef (step S707: Yes), or when the counter i is greater than the number m of past images input from the storage unit 36.
  • step S709 the image processing apparatus 30a proceeds to step S709 described later.
  • Step S707 If the number of pixels of mVal [] stored in the storage unit 36 is not greater than the preset threshold ThRef in step S707, and the counter i is not greater than the number m of past images input from the storage unit 36 (step S707). : No), the image processing apparatus 30a proceeds to Step S708 to be described later.
  • Step S709 and step S710 correspond to step S612 and step S613 of FIG. 27 described above, respectively.
  • the representative value calculation unit 327 may store an uncorrected image that is temporally continuous with the current image in the storage unit 36 and use the uncorrected past image instead of the past image.
  • the representative value calculation unit 327 calculates the representative value using the pixel values of the peripheral pixels of the past image, thereby improving the reliability of the representative value. be able to.
  • FIG. 31 is a block diagram schematically showing the configuration of the imaging system 2 according to Embodiment 5 of the present invention.
  • the imaging system 2 illustrated in FIG. 31 includes a main body unit 3 and a lens unit 4 that can be detachably connected to the main body unit 3.
  • the main unit 3 includes a shutter 103, an image sensor 105, an analog processing unit 106, an A / D conversion unit 107, an operation unit 108, a memory I / F unit 109, a recording medium 110, and a volatile memory 111. , Nonvolatile memory 112, bus 113, imaging control unit 114, AE processing unit 116, AF processing unit 117, external I / F unit 118, display unit 119, driver 120, and RTS noise correction unit 320. And comprising.
  • the driver 120 drives the shutter 103 under the control of the imaging control unit 114.
  • the AE processing unit 116 acquires image data stored in the volatile memory 111 via the bus 113, and sets an exposure condition for performing still image shooting or moving image shooting based on the acquired image data. Specifically, the AE processing unit 116 calculates the luminance from the image data, and determines the automatic exposure (Auto Exposure) of the imaging system 2 by determining, for example, the aperture value, the exposure time, the ISO sensitivity, and the like based on the calculated luminance. )I do.
  • the AF processing unit 117 acquires the image data stored in the volatile memory 111 via the bus 113, and adjusts the automatic focus of the imaging system 2 based on the acquired image data. For example, the AF processing unit 117 extracts a high-frequency component signal from the image data, performs AF (Auto Focus) calculation processing on the high-frequency component signal, and determines the focus evaluation of the imaging system 2 to perform imaging. The system 2 auto focus adjustment is performed. Note that the method of adjusting the automatic focus of the imaging system 2 may acquire the phase difference signal by the imaging element 105.
  • the external I / F unit 118 can perform reading and writing of data in various blocks in the main unit 3 and control by a dedicated command or the like.
  • the external I / F unit 118 is an interface capable of controlling various blocks in the main unit 3 by connecting an external device such as a dedicated circuit or a personal computer (PC) equipped with an FPGA, DSP, GPU, or the like. .
  • the display unit 119 is configured using a display panel made of liquid crystal, organic EL (Electro Luminescence), or the like.
  • the display unit 119 displays an image corresponding to the image data generated by the image sensor 105.
  • the lens unit 4 forms an object image condensed from a predetermined field of view on the image sensor 105.
  • An optical system 101, a diaphragm 102, and a driver 104 are provided.
  • FIG. 32 is a flowchart illustrating an outline of processing executed by the imaging system 2.
  • the imaging control unit 114 initializes the imaging system 2. This is performed (step S801). Specifically, the imaging control unit 114 performs initialization to turn off a recording flag indicating that a moving image is being recorded. This recording flag is a flag that is turned on during moving image shooting and is turned off when no moving image is being shot, and is stored in the volatile memory 111.
  • step S802 when the moving image button of the operation unit 108 is pressed (step S802: Yes), the imaging control unit 114 inverts the recording flag indicating that the moving image is being recorded in the ON state (step S803).
  • the imaging control unit 114 determines whether or not the imaging system 2 is recording a moving image (step S804). Specifically, the imaging control unit 114 determines whether or not the recording flag stored in the volatile memory 111 is on.
  • step S804: Yes the imaging system 2 proceeds to step S805 described later.
  • step S804: No the imaging system 2 proceeds to step S806 described later.
  • step S805 the imaging control unit 114 generates a moving image file for recording the image data on the recording medium 110 in time series.
  • step S806 the imaging system 2 proceeds to step S806 described below.
  • step S802 when the moving image button of the operation unit 108 is not pressed (step S802: No), the imaging system 2 proceeds to step S806.
  • the imaging control unit 114 determines whether or not the imaging system 2 is recording a moving image (step S806).
  • the imaging control unit 114 determines that the imaging system 2 is recording a moving image (step S806: Yes)
  • the imaging system 2 proceeds to step S817 described later.
  • the imaging control unit 114 determines that the moving image is not being recorded (step S806: No)
  • the imaging system 2 proceeds to step S807 described later.
  • step S807 If the playback button of the operation unit 108 is pressed in step S807 (step S807: Yes), the imaging system 2 plays back and displays an image corresponding to the image data recorded on the recording medium 110 on the display unit 119. (Step S808). After step S808, the imaging system 2 proceeds to step S809 described later.
  • step S807 when the playback button of the operation unit 108 is not pressed (step S807: No), the imaging system 2 proceeds to step S809.
  • step S809 Yes
  • the imaging system 2 executes a setting process for performing various settings (step S810). After step S810, the imaging system 2 proceeds to step S811 described below.
  • step S809 when the menu button of the operation unit 108 is not pressed (step S809: No), the imaging system 2 proceeds to step S811.
  • step S811 when the release button of the operation unit 108 transitions from the off state to the 1st state (step S811: Yes), the imaging control unit 114 causes the AE processing unit 116 to adjust the exposure to the AE processing and AF processing unit 117. Each AF process for adjusting the focus is executed (step S812). Thereafter, the imaging system 2 proceeds to Step S824 described later.
  • step S811 when the release button of the operation unit 108 has not transitioned from the off state to the 1st state (step S811: No), the imaging system 2 proceeds to step S813.
  • the imaging control unit 114 performs imaging using the mechanical shutter (step S814). Specifically, the imaging control unit 114 controls the shutter 103 to cause the imaging element 105 to perform imaging.
  • the imaging system 2 performs image processing for performing predetermined processing after correcting the RTS noise on the image data generated by the imaging element 105 (step S815). Details of the image processing will be described later.
  • the imaging control unit 114 records the image data subjected to the image processing by the image processing unit 303 on the recording medium 110 (step S816). After step S816, the imaging system 2 proceeds to step S824 described later.
  • step S813 when the release button of the operation unit 108 has not transitioned to the 2nd state (step S813: No), the imaging system 2 proceeds to step S817.
  • the imaging control unit 114 causes the AE processing unit 116 to execute AE processing for adjusting exposure (step S817), and causes the AF processing unit 117 to perform AF processing for adjusting focus (step S818).
  • the imaging control unit 114 causes the imaging element 105 to electronically control the exposure time via the driver 120, so as to execute imaging with a so-called electronic shutter (step S819).
  • Image data generated by the image sensor 105 by photographing with the electronic shutter is output to the volatile memory 111 via the analog processing unit 106, the A / D conversion unit 107, and the bus 113.
  • the imaging system 2 executes image processing similar to that in step S815 (step S820). Details of the image processing will be described later.
  • the imaging system 2 causes the display unit 119 to display a live view image corresponding to the image data generated by the imaging element 105 by photographing with the electronic shutter (step S821).
  • step S822 when the imaging system 2 is recording a moving image (step S822: Yes), the imaging control unit 114 compresses the image data in an image compression / expansion unit (not shown) in the recording format set by the setting process in step S810, The compressed image data is recorded as a moving image in the moving image file created on the recording medium 110 (step S823). After step S823, the imaging system 2 proceeds to step S824.
  • step S822 when the imaging system 2 is not recording a moving image (step S822: No), the imaging system 2 proceeds to step S824.
  • step S824: Yes when the power button of the operation unit 108 is pressed and the power of the imaging system 2 is turned off (step S824: Yes), the imaging system 2 ends this process. On the other hand, when the power supply of the imaging system 2 is not turned off (step S824: No), the imaging system 2 returns to step S802.
  • FIG. 33 is a flowchart showing an overview of image processing.
  • the RTS noise correction unit 320 executes RTS noise correction processing for correcting RTS noise on the image data generated by the image sensor 105 (step S901).
  • the RTS noise correction process corresponds to the process executed by the image processing apparatus 30 according to the first embodiment or the process executed by the image processing apparatus 30a according to the second to fourth embodiments. Is omitted.
  • the image processing unit 303 performs basic image processing on the image data whose RTS noise has been corrected by the RTS noise correction unit 320 (step S902). After step S902, the imaging system 2 returns to the main routine of FIG.
  • the present invention is not limited to the above-described embodiments, and various modifications and applications are naturally possible within the scope of the gist of the present invention.
  • imaging that captures a subject through an optical device such as a portable device, video camera, endoscope, surveillance camera, or microscope equipped with an imaging device in a mobile phone or smartphone
  • the present invention can be applied to any device that can image a subject such as a device.
  • the reference value of the random noise model or the noise model to be referred to is changed according to RTS_Value, which is a feature quantity of RTS noise, and the maximum value of candidate values, which are values based on RTS_Value.
  • RTS_Value is a feature quantity of RTS noise
  • candidate values which are values based on RTS_Value.
  • a random noise model may be selected according to the combination and the RTS noise feature amount, and the reference value may be changed to estimate the random noise amount. Further, since the amount of random noise varies depending on the temperature, a random noise model may be selected according to the temperature of the image sensor when the image data is captured, and the above method may be applied.
  • the RTS noise information recording unit is provided in the imaging device.
  • the RTS noise information recording unit may be provided in the above-described image processing device or in a server capable of bidirectional communication via a network.
  • An RTS noise information recording unit may be provided to acquire RTS noise information via a network.
  • the present invention is also applicable to image data other than image data used for display or recording, such as image data of an OB area or image data of an area outside an image circle that is not optically guaranteed. Is possible.
  • each processing method by the image processing apparatus in the above-described embodiment can be stored as a program that can be executed by a control unit such as a CPU.
  • a control unit such as a CPU.
  • memory cards ROM cards, RAM cards, etc.
  • magnetic disks floppy disks (registered trademark), hard disks, etc.
  • optical disks CD-ROMs, DVDs, etc.
  • storage media such as semiconductor memories, etc.
  • a control unit such as a CPU reads the program stored in the storage medium of the external storage device, and the operation described above can be executed by the operation being controlled by the read program.
  • the present invention is not limited to the above-described embodiments and modifications as they are, and in the implementation stage, the constituent elements can be modified and embodied without departing from the spirit of the invention.
  • Various inventions can be formed by appropriately combining a plurality of constituent elements disclosed in the above-described embodiments. For example, some constituent elements may be deleted from all the constituent elements described in the above-described embodiments and modifications. Furthermore, you may combine suitably the component demonstrated by each embodiment and the modification.

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Studio Devices (AREA)
  • Image Processing (AREA)
  • Transforming Light Signals Into Electric Signals (AREA)
  • Facsimile Image Signal Circuits (AREA)
  • Picture Signal Circuits (AREA)

Abstract

 注目画素の画素値と周辺画素の画素値との相関が低い場合であっても、画質を向上させることができる画像処理装置、画像処理方法およびプログラムを提供する。画像処理装置30は、読み出し回路の位置情報または複数の画素の各々の位置情報を含むノイズ情報を取得するとともに、撮像素子105が生成した複数の画像データを取得する第3外部I/F部31と、第3外部I/F部31が取得した画像データを補正対象画像データとし、この補正対象画像データと異なる時間に取得された画像データに基づく参照画像データと、ノイズ情報とに基づいて、補正対象画像データにおける注目画素の画素値を補正するノイズ低減部35と、を備える。

Description

画像処理装置、画像処理方法およびプログラム
 本発明は、2次元状に配置された複数の画素を有する撮像素子に生じるRTSノイズのような画素値が一定範囲内で変動する点滅欠陥ノイズを補正する画像処理装置、画像処理方法およびプログラムに関する。
 近年、CMOS(Complementary Metal Oxide Semiconductor)等の撮像素子は、画素および該画素から信号を読み出す読み出し回路の微細化が進んでいる。このような微細化された読み出し回路では、感度の低下および様々なノイズの増加が問題となっている。感度の低下に対しては、複数の画素を1つの読み出し回路で共有させて信号を読み出す共有画素構造を採用することで、撮像素子における回路に必要な面積を削減し、各画素の開口率(受光部の割合)を向上させることによって、感度を向上させている。
 一方、撮像素子で発生するノイズには、暗電流による暗電流ショットノイズおよび読み出し回路での熱雑音等に起因するランダムノイズ以外にも、画素値が常に異常値を示す欠陥画素、および画素値がランダムに変動する点滅欠陥ノイズ等がある。このような点滅欠陥ノイズの中には、読み出し回路に起因するRTS(Random Telegraph Signal)ノイズがある。このRTSノイズを補正する技術として、撮影された画像における注目画素の画素値、該注目画素の周辺画素の画素値および予め撮像素子の画素毎に検出したRTSノイズのノイズレベル(以下、「RTSノイズレベル」という)に基づいて、注目画素に対してRTSノイズの影響があるか否かを判定し、RTSノイズの影響があると判定された場合、この注目画素の画素値に対してRTSノイズレベルだけ加算または減算する技術が知られている(特許文献1参照)。
特開2012-105063号公報
 ところで、上述した特許文献1では、注目画素の画素値と該注目画像の周辺画素の画素値との相関が高い場合にはRTSノイズを低減することができる。しかしながら、上述した特許文献1では、輝点または微細構造物等の被写体によっては、注目画素の画素値と周辺画素の画素値との相関が低くなることで、画質を劣化させる場合があった。
 本発明は、上記に鑑みてなされたものであって、注目画素の画素値と周辺画素の画素値との相関が低い場合であっても、画質を向上させることができる画像処理装置、画像処理方法およびプログラムを提供することを目的とする。
 上述した課題を解決し、目的を達成するために、本発明に係る画像処理装置は、2次元状に配置され、外部から光を受光し、受光量に応じた信号を生成する複数の画素と、前記信号を画素値として読み出す複数の読み出し回路と、を有する撮像素子が生成した画像データに含まれる点滅欠陥ノイズを補正する画像処理装置であって、前記読み出し回路に起因する点滅欠陥ノイズが発生する前記読み出し回路の位置情報または前記複数の画素の各々の位置情報を含むノイズ情報と、前記画像データと、を取得する取得部と、前記取得部が取得した前記画像データを補正対象画像データとし、該補正対象画像データと異なる時間に取得された前記画像データに基づく参照画像データと、前記ノイズ情報とに基づいて、前記補正対象画像データにおける注目画素の画素値を補正する補正部と、を備えたことを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記補正対象画像データと前記参照画像データとに基づいて、被写体の移動量を算出する移動量算出部と、を備え、前記補正部は、前記移動量算出部が算出した前記移動量に基づいて、前記注目画素または前記注目画素近傍の画素と対応する前記参照画像データの参照画素を取得し、該参照画素の画素値に基づいて、前記注目画素の画素値を補正することを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記補正部は、前記注目画素において前記点滅欠陥ノイズが発生する場合、前記点滅欠陥ノイズが発生しない場合に比べ、前記注目画素におけるノイズ低減処理の強度を大きくすることを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記補正部は、前記注目画素において前記点滅欠陥ノイズが発生する場合、前記参照画像データを用いた時間方向のノイズ低減処理の強度を大きくすることを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記補正部は、前記注目画素において前記点滅欠陥ノイズが発生する場合、前記注目画素周辺の画素を用いた空間方向のノイズ低減処理の強度を大きくすることを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記補正部は、前記参照画素の画素値に基づいて、前記点滅欠陥ノイズが発生しない場合の画素値に相当する代表値を算出し、該代表値に基づいて、前記注目画素の画素値を補正することを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記補正部は、前記注目画素の画素値に基づいて、前記代表値を算出することを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記補正部は、前記点滅欠陥ノイズが発生しない前記参照画素の画素値もしくは前記点滅欠陥ノイズを補正した前記参照画素の画素値に基づいて、前記代表値を算出することを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記補正部は、前記点滅欠陥ノイズが発生しない前記注目画素における周辺の画素の画素値に基づいて、前記代表値を算出することを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記注目画素周辺のランダムノイズ量を推定するノイズ量推定部と、をさらに備え、前記補正部は、前記ランダムノイズ量に基づいて、前記代表値を算出することを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記ノイズ情報は、前記読み出し回路に起因する点滅欠陥ノイズが発生する前記読み出し回路の位置情報または前記複数の画素の各々の位置情報と対応づけられた、前記点滅欠陥ノイズのノイズレベルである点滅欠陥ノイズレベルをさらに含み、前記補正部は、前記点滅欠陥ノイズレベルを用いて前記注目画素の画素値を補正することを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記ノイズ情報は、前記読み出し回路に起因する点滅欠陥ノイズが発生する前記読み出し回路の位置情報または前記複数の画素の各々の位置情報と対応づけられた、前記点滅欠陥ノイズのノイズレベルである点滅欠陥ノイズレベルをさらに含み、前記補正部は、前記点滅欠陥ノイズレベルが大きいほど、前記注目画素におけるノイズ低減処理の強度を大きくすることを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記ノイズ情報は、前記読み出し回路に起因する点滅欠陥ノイズが発生する前記読み出し回路の位置情報または前記複数の画素の各々の位置情報と対応づけられた、前記点滅欠陥ノイズのノイズレベルである点滅欠陥ノイズレベルをさらに含み、前記補正部は、前記点滅欠陥ノイズレベルに基づいて、前記代表値を算出することを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記移動量算出部が算出した前記移動量の確からしさを表す信頼度を算出する信頼度算出部を備え、前記補正部は、前記信頼度に基づいて、前記代表値を算出することを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記参照画像データは、前記補正対象画像データに対応する前記画像データの直前に前記撮像素子が生成した前記画像データであることを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記参照画像データは、前記補正部が前記補正対象画像データに対する補正を行う直前に、前記補正部が補正を行った画像データであることを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記参照画像データは、前記補正対象画像データに対応する前記画像データの前に、前記撮像素子が生成した複数の前記画像データであることを特徴とする。
 また、本発明に係る画像処理装置は、上記発明において、前記点滅欠陥ノイズは、ランダムテレグラフシグナルノイズであることを特徴とする。
 また、本発明に係る画像処理方法は、2次元状に配置され、外部から光を受光し、受光量に応じた信号を生成する複数の画素と、前記信号を画素値として読み出す複数の読み出し回路と、を有する撮像素子が生成した画像データに含まれる点滅欠陥ノイズを補正する画像処理装置が実行する画像処理方法であって、前記読み出し回路に起因する点滅欠陥ノイズが発生する前記読み出し回路の位置情報または前記複数の画素の各々の位置情報を含むノイズ情報と、前記画像データと、を取得する取得ステップと、前記取得ステップにおいて取得した前記画像データを補正対象画像データとし、該補正対象画像データと異なる時間に取得された前記画像データに基づく参照画像データと、前記ノイズ情報とに基づいて、前記補正対象画像データにおける注目画素の画素値を補正する補正ステップと、を含むことを特徴とする。
 また、本発明に係るプログラムは、2次元状に配置され、外部から光を受光し、受光量に応じた信号を生成する複数の画素と、前記信号を画素値として読み出す複数の読み出し回路と、を有する撮像素子が生成した画像データに含まれる点滅欠陥ノイズを補正する画像処理装置に、前記読み出し回路に起因する点滅欠陥ノイズが発生する前記読み出し回路の位置情報または前記複数の画素の各々の位置情報を含むノイズ情報と、前記画像データと、を取得する取得ステップと、前記取得ステップにおいて取得した前記画像データを補正対象画像データとし、該補正対象画像データと異なる時間に取得された前記画像データに基づく参照画像データと、前記ノイズ情報とに基づいて、前記補正対象画像データにおける注目画素の画素値を補正する補正ステップと、を実行させることを特徴とする。
 本発明によれば、注目画素の画素値と周辺画素の画素値との相関が低い場合であっても、画質を向上させることができるという効果を奏する。
図1は、本発明の実施の形態1に係る撮像システムの構成を模式的に示すブロック図である。 図2は、本発明の実施の形態1に係る撮像装置が備える撮像素子の要部の構成を模式的に示す概略図である。 図3は、本発明の実施の形態1に係る撮像素子に光が当たらないように遮光した場合において、RTSノイズが発生するときに、アンプ部から出力されるアンプ出力の変動の例を示す図である。 図4は、本発明の実施の形態1に係るRTSノイズが発生するアンプ部を用いて読み出された画素値の分布の例を示す図である。 図5は、本発明の実施の形態1に係る画像処理装置が実行する処理の概要を示すフローチャートである。 図6は、本発明の実施の形態1に係る画像処理装置が備えるルックアップテーブルが記憶するノイズモデルの一例を示す図である。 図7は、図5の時間方向NR処理の概要を示すフローチャートである。 図8は、本発明の実施の形態1に係る画像処理装置が備えるノイズ低減処理部が設定するRTS_Valueと係数Co_we_preとの関係を示す図である。 図9は、図5の空間方向NR処理の概要を示すフローチャートである。 図10は、本発明の実施の形態1に係る画像処理装置が備えるノイズ低減処理部が設定するRTS_Valueと係数Co_Cbとの関係を示す図である。 図11は、本発明の実施の形態2に係る撮像システムの構成を模式的に示すブロック図である。 図12は、本発明の実施の形態2に係るRTSノイズ補正部の詳細な構成を示すブロック図である。 図13は、本発明の実施の形態2に係る代表値算出部の詳細な構成を示すブロック図である。 図14は、本発明の実施の形態2に係る補正値算出部の詳細な構成を示すブロック図である。 図15は、本発明の実施の形態2に係る画像処理装置が実行する処理の概要を示すフローチャートである。 図16は、図15の代表値算出処理の概要を示すフローチャートである。 図17は、図15の補正値算出処理の概要を示すフローチャートである。 図18は、本発明の実施の形態3に係るノイズ低減部の詳細な構成を示すブロック図である。 図19は、本発明の実施の形態3に係るRTSノイズ補正部の詳細な構成を示すブロック図である。 図20は、本発明の実施の形態3に係る画像処理装置が実行する代表値算出処理の概要を示すフローチャートである。 図21は、本発明の実施の形態3の変形例に係るRTSノイズ補正部の詳細な構成を示すブロック図である。 図22は、本発明の実施の形態3の変形例に係る代表値算出処理の概要を示すフローチャートである。 図23は、本発明の実施の形態3の変形例に係る代表値決定部が算出するブレンド率の算出方法を模式的に示す図である。 図24は、本発明の実施の形態4に係るノイズ低減部の詳細な構成を示すブロック図である。 図25は、本発明の実施の形態4に係るRTSノイズ補正部の詳細な構成を示すブロック図である。 図26は、本発明の実施の形態4に係る画像処理装置が実行する代表値算出処理の概要を示すフローチャートである。 図27は、本発明の実施の形態4の変形例1に係るRTSノイズ補正部の詳細な構成を示すブロック図である。 図28は、本発明の実施の形態4の変形例1に係る画像処理装置が実行する代表値算出処理の概要を示すフローチャートである。 図29は、本発明の実施の形態4の変形例2に係る画像処理装置が実行する代表値算出処理の概要を示すフローチャートである。 図30は、本発明の実施の形態4の変形例3に係る画像処理装置が実行する代表値算出処理の概要を示すフローチャートである。 図31は、本発明の実施の形態5に係る撮像システムの構成を模式的に示すブロック図である。 図32は、本発明の実施の形態5に係る撮像システムが実行する処理の概要を示すフローチャートである。 図33は、図30の画像処理の概要を示すフローチャートである。
 以下、図面を参照して、本発明を実施するための形態(以下、「実施の形態」という)について説明する。なお、以下に説明する実施の形態によって本発明が限定されるものではない。さらに、図面の記載において、同一の部分には同一の符号を付して説明する。
(実施の形態1)
 〔撮像システムの構成〕
 図1は、本発明の実施の形態1に係る撮像システムの構成を模式的に示すブロック図である。図1に示す撮像システム1は、撮像装置10と、画像処理装置30と、表示装置40と、を備える。
 〔撮像装置の構成〕
 まず、撮像装置10の構成について説明する。撮像装置10は、図1に示すように、光学系101と、絞り102と、シャッタ103と、ドライバ104と、撮像素子105と、アナログ処理部106と、A/D変換部107と、操作部108と、メモリI/F部109と、記録媒体110と、揮発メモリ111と、不揮発メモリ112と、バス113と、撮像制御部114と、第1外部I/F部115と、を備える。
 光学系101は、単数または複数のレンズを用いて構成される。光学系101は、例えばフォーカスレンズとズームレンズとを用いて構成される。
 絞り102は、光学系101が集光した光の入射量を制限することで露出の調整を行う。絞り102は、後述する撮像制御部114の制御のもと、光学系101が集光した光の入射量を制限する。なお、絞り102を用いずに、シャッタ103や撮像素子105における電子シャッタを用いて光の入射量を制限するようにしてもよい。なお、光学系101、絞り102は、撮像装置10に対して着脱可能であってもよい。
 シャッタ103は、撮像素子105の状態を露光状態または遮光状態に設定する。シャッタ103は、例えばフォーカルプレーンシャッタ等を用いて構成される。なお、シャッタ103を用いずに、撮像素子105における電子シャッタを用いてもよい。
 ドライバ104は、後述する撮像制御部114の制御のもと、光学系101、絞り102およびシャッタ103を駆動する。例えば、ドライバ104は、光学系101を光軸O1に沿って移動させることによって、撮像装置10のズーム倍率の変更またはピント位置の調整を行う。
 撮像素子105は、後述する撮像制御部114の制御のもと、光学系101が集光した光を受光して画像データ(電気信号)に変換して出力する。撮像素子105は、複数の画素が2次元状に配置されたCMOS(Complementary Metal Oxide Semiconductor)等を用いて構成される。この各画素の前面には、ベイヤー配列のRGBフィルタが配置されている。なお、撮像素子105は、ベイヤー配列に限定されず、例えばFovionのような積層型の形式でも勿論かまわない。また、用いるフィルタはRGBに限定されず、補色フィルタ等任意のフィルタを適用できる。また、別途、異なるカラー光を時分割で照射可能な光源を配置し、撮像素子105には、フィルタを配置せず、照射する色を変更しながら順次取り込んだ画像を使用してカラー画像を構成できるようにしてもよい。また、撮像素子105は、受光量を電子的に制御可能な電子シャッタ機能を有する。
 ここで、撮像素子105の構成について詳細に説明する。図2は、撮像素子105の要部の構成を模式的に示す概略図である。なお、図2に示す撮像素子105は、画素の開口率向上により感度を向上させるため、複数の画素で読み出し回路を共有している例を示している。なお、図2に示す撮像素子105は、水平方向(横方向)に2画素×垂直方向(縦方向)に4画素の8画素に対して、1つの読み出し回路が配置されている。なお、図2においては、水平方向(横方向)に2画素×垂直方向(縦方向)に4画素の8画素に対して、1つの読み出し回路を1グループとする例を説明したが、本実施の形態1の撮像素子105上には、上述した画素および読み出し回路が、水平方向および垂直方向に並んで配置されているものとする。
 図2に示すように、撮像素子105は、露光により光を受光し、光電変換を行うことによって、露光量に対応した電荷を発生する複数の画素105a(フォトダイオード)と、複数の画素105aの各々に設けられ、撮像制御部114の制御に応じて開閉する第1スイッチ105bと、複数の画素105aの各々から出力された信号(電荷)を垂直方向に転送する垂直転送線105cと、複数の画素105aの各々から出力された信号を蓄積するFD部105d(Floating Diffusion)と、FD部105dから出力された信号を増幅するアンプ部105eと、撮像制御部114の制御に応じて開閉する第2スイッチ105fと、第2スイッチ105fを制御する制御線105gと、アンプ部105eで増幅された電気信号を転送する転送線105hと、を備える。
 このように構成された撮像素子105は、画素105a(1)~105a(8)における露光量に対応する信号を画素値として読み出す場合、まず、FD部105dをリセット状態にして、撮像制御部114が第1スイッチ105b(1)のみをオンとすることで、画素105a(1)に発生した電荷をFD部105dに転送する。その後、撮像素子105は、撮像制御部114が第2スイッチ105fをオンとすることで、FD部105dに蓄積された電荷をアンプ部105eによって増幅させて画素値として読み出す(出力する)。次に、撮像素子105は、FD部105dをリセット状態にして、撮像制御部114が第1スイッチ105b(2)のみをオンとすることで、画素105a(2)に発生した電荷をFD部105dに転送する。その後、撮像素子105は、撮像制御部114が第2スイッチ105fをオンとすることで、FD部105dに蓄積された電荷をアンプ部105eによって増幅させて画素値として読み出す。撮像素子105は、このような読み出し動作を順次行うことによって、画素105a(1)~105a(8)における露光量に対応する信号を順次画素値として出力することができる。なお、本実施の形態1では、アンプ部105eが複数の画素105aの各々から電荷を読み出す読み出し回路として機能する。
 図1に戻り、撮像装置10の構成の説明を続ける。
 アナログ処理部106は、撮像素子105から入力されるアナログ信号に対して、所定のアナログ処理を施してA/D変換部107へ出力する。具体的には、アナログ処理部106は、撮像素子105から入力されるアナログ信号に対して、ノイズ低減処理およびゲインアップ処理等を行う。例えば、アナログ処理部106は、アナログ信号に対して、リセットノイズ等を低減した上で波形整形を行い、さらに目的の明るさとなるようにゲインアップを行う。
 A/D変換部107は、アナログ処理部106から入力されるアナログ信号に対して、A/D変換を行うことによってデジタルの画像データ(以下、「RAW画像データ」という)を生成し、バス113を介して揮発メモリ111に出力する。なお、A/D変換部107は、後述する撮像装置10の各部に対してRAW画像データを直接出力するようにしてもよい。なお、上述したアナログ処理部106とA/D変換部107を撮像素子105に設け、撮像素子105がデジタルのRAW画像データを直接出力するようにしてもよい。
 操作部108は、撮像装置10の各種の指示を与える。具体的には、操作部108は、撮像装置10の電源状態をオン状態またはオフ状態に切り替える電源スイッチ、静止画撮影の指示を与えるレリーズスイッチ、撮像装置10の各種設定を切り替える操作スイッチおよび動画撮影の指示を与える動画スイッチ等を有する。
 記録媒体110は、撮像装置10の外部から装着されるメモリカードを用いて構成され、メモリI/F部109を介して撮像装置10に着脱自在に装着される。また、記録媒体110は、後述する撮像制御部114の制御のもと、メモリI/F部109を介してプログラムおよび各種情報それぞれを不揮発メモリ112に出力してもよい。
 揮発メモリ111は、バス113を介してA/D変換部107から入力される画像データを一時的に記憶する。例えば、揮発メモリ111は、アナログ処理部106、A/D変換部107およびバス113を介して、撮像素子105が1フレーム毎に順次出力する画像データを一時的に記憶する。揮発メモリ111は、SDRAM(Synchronous Dynamic Random Access Memory)等を用いて構成される。
 不揮発メモリ112は、Flashメモリ等を用いて構成され、撮像装置10を動作させるための各種プログラム、プログラムの実行中に使用される各種データを記録する。また、不揮発メモリ112は、プログラム記録部112aと、画素値を読み出す読み出し回路(アンプ部105e)の位置情報または複数の画素105aの各々の位置情報と読み出し回路(アンプ部105e)に起因するRTSノイズに関する特徴量とを対応付けたRTSノイズ位置情報を記録するRTSノイズ情報記録部112bと、一または複数のランダムノイズモデルを記録するランダムノイズモデル情報記録部112cと、を有する。ここで、特徴量とは、RTSノイズの振幅(RTS_Value)、点滅欠陥ノイズの発生頻度およびRTSノイズの振幅未満のRTSノイズにおける発生頻度のいずれかである。
 バス113は、撮像装置10の各構成部位を接続する伝送路等を用いて構成され、撮像装置10の内部で発生した各種データを撮像装置10の各構成部位に転送する。
 撮像制御部114は、CPU(Central Processing Unit)等を用いて構成され、操作部108からの指示信号やレリーズ信号に応じて撮像装置10を構成する各部に対する指示やデータの転送等を行って撮像装置10の動作を統括的に制御する。例えば、撮像制御部114は、操作部108からセカンドレリーズ信号が入力された場合、撮像装置10における撮影動作の制御を開始する。ここで、撮像装置10における撮影動作とは、撮像素子105の露光タイミング、アナログ信号の出力タイミング、および撮像素子105が出力したアナログ信号に対し、アナログ処理部106かつA/D変換部107が所定の処理を施す動作をいう。このように処理が施された画像データは、撮像制御部114の制御のもと、バス113およびメモリI/F部109を介して記録媒体110に記録される。
 第1外部I/F部115は、バス113を介して外部の機器から入力される情報を不揮発メモリ112または揮発メモリ111へ出力する一方、バス113を介して外部の機器へ揮発メモリ111が記憶する情報、不揮発メモリ112が記憶する情報および撮像素子105が生成した画像データを出力する。具体的には、第1外部I/F部115は、バス113を介して画像処理装置30に撮像素子105が生成したRAW画像データ、RTSノイズ情報およびランダムノイズモデル情報を出力する。
 〔画像処理装置の構成〕
 次に、画像処理装置30の構成について説明する。画像処理装置30は、第3外部I/F部31と、補間処理部32と、フレームメモリ33と、動き検出部34と、ノイズ低減部35と、記憶部36と、操作部37と、画像処理部38と、画像処理制御部39と、を備える。
 第3外部I/F部31は、撮像装置10の第1外部I/F部115を介して撮像素子105によって生成された画像データ、不揮発メモリ112内のRTSノイズ情報記録部112bが記録するRTSノイズに関するRTSノイズ情報およびランダムノイズモデル情報記録部112cが記録するランダムノイズモデル情報を取得し、取得したRAW画像データ(補正対象画像データ)、RTSノイズ情報およびランダムノイズモデル情報を補間処理部32および記憶部36へ出力する。第3外部I/F部31および第1外部I/F部115は、例えば双方向に情報をやり取り可能な制御ケーブルや無線通信等を介して接続されている。なお、本実施の形態1では、第3外部I/F部31が取得部として機能する。
 補間処理部32は、第3外部I/F部31から入力された撮像素子105によって生成されたRAW画像データに対して補間処理を行う。具体的には、補間処理部32は、撮像素子105がベイヤー配列の場合、RAW画像データに対応するRAW画像(以下、「現画像」という)の各画素に対して補間処理を行うことによって、欠落している画素値(信号値)を補間し、各画素においてR、GおよびB信号の全ての画素値を有する現画像を生成し、この現画像をノイズ低減部35へ出力する。
 フレームメモリ33は、ノイズ低減部35と双方向に接続され、表示装置40に表示するための画像のフレームデータを一時的に記憶する。
 動き検出部34は、補間処理部32から入力される現画像およびフレームメモリ33に記憶され、ノイズ低減部35によってノイズが補正された現画像と、この現画像と時間的に連続した補正後画像データに対応する補正後画像(以下、「過去画像」という)に基づいて、被写体の動きベクトル(移動量)を検出し、この検出結果をノイズ低減部35へ出力する。動き検出部34は、例えば公知のブロックマッチング処理等を用いて動きベクトルを検出する。例えば、動き検出部34は、現画像の注目画素の特徴量と現画像の注目画素の特徴量位置に対応する過去画像の特徴量を検出後、この特徴量と現画像の注目画素の特徴量との距離が最小となる過去画像の画素を検出し、この検出した画素の位置と現画像の注目画素の画素位置とに基づいて、動きベクトルを検出する。ここで、現画像の注目画素の特徴量および現画像の注目画素の特徴量位置に対応する過去画像の特徴量とは、画素値や輝度値等である。なお、動き検出部34は、現画像の各画素について動きベクトルを検出してもよい。また、本実施の形態1では、動き検出部34が移動量算出部として機能する。
 ノイズ低減部35は、補間処理部32から入力された現画像に対して、ノイズ低減処理(以下、「NR処理」という)を行う。具体的には、ノイズ低減部35は、現画像の各画素において静止状態と動作状態とを判別し、その判別結果に応じてNR処理を切り替える。例えば、ノイズ低減部35は、現画像が静止状態と判別した場合、高周波成分を保持可能な時間方向におけるノイズを低減する第1のNR処理(以下、「時間方向NR処理」という)を選択する一方、現画像が動作状態と判定した場合、空間方向におけるノイズを低減する第2のNR処理(以下「空間方向NR処理」という)を選択する。ここで、時間方向NR処理とは、NR処理の対象となる現画像と、その現画像とは異なる時間に取得された過去画像(参照画像データ)と、を用いた加重平均処理によりノイズを低減する。加重平均処理には、現画像の注目画素(処理対象画素)と現画像の注目画素の位置に対応する過去画像の画素(処理対象画素)が用いられる。また、空間方向NR処理とは、NR処理の処理対象となる注目画素(処理対象画素)と、この注目画素の周辺の周辺画素とを用いた加重平均処理によってノイズを低減する。また、静止状態とは、撮像素子105(撮像装置10)と被写体との相対的な位置関係が時間的に変化しない状態を示す。さらにまた、動作状態とは、撮像素子105(撮像装置10)と被写体との相対的な位置関係が時間的に変化する状態を示す。なお、本実施の形態1では、ノイズ低減部35が補正部として機能する。また、ノイズ低減部35は、被写体が移動していない場合、現画像の注目画素近傍の画素と対応する過去画像の画素の画素値を記憶部36から取得し、取得した過去画像の画素の画素値に基づいて、注目画素の画素値を補正してもよい。 
 ここで、ノイズ低減部35の詳細な構成について説明する。ノイズ低減部35は、評価値算出部351と、推定ノイズ量取得部352と、ルックアップテーブル353と、判定部354と、ノイズ低減処理部355と、を有する。
 評価値算出部351は、補間処理部32から入力された時間的に連続する現画像およびフレームメモリ33から入力された過去画像の各々に含まれる被写体がフレーム間において静止状態であるか否かを判定するための評価値を算出する。
 推定ノイズ量取得部352は、補間処理部32から入力された画像データに対応する画像で推定される推定ノイズ量を取得する。
 ルックアップテーブル353は、推定ノイズ量取得部352がノイズ量を推定する際に用いられるノイズモデルを記憶する。
 判定部354は、評価値算出部351が算出した評価値と推定ノイズ量取得部352が取得した推定ノイズ量とに基づいて、補間処理部32から入力された現画像内における被写体が静止状態であるか否かを判定する。
 ノイズ低減処理部355は、注目画素がRTSノイズを発生しえるとき、ノイズ低減の強度を強くすることによって注目画素の画素値を補正して画像処理部38へ出力する。
 記憶部36は、揮発メモリや不揮発性メモリを用いて構成され、第3外部I/F部31から出力されたRTSノイズ情報および複数の過去画像を記憶する。
 操作部37は、画像処理装置30に関する各種操作信号の入力を受け付ける。操作部37は、例えば十字ボタン、プッシュボタンおよびタッチパネル等を用いて構成される。
 画像処理部38は、ノイズ低減されたRAW画像データに対して、所定の画像処理を行って表示装置40へ出力する。ここで、所定の画像処理とは、少なくともホワイトバランス調整処理、γ補正処理、色再現処理およびエッジ強調処理等を含む基本の画像処理を行う。また、画像処理部38は、予め設定された各画像処理のパラメータに基づいて、自然な画像を再現する画像処理を行う。ここで、各画像処理のパラメータとは、コントラスト、シャープネス、彩度、ホワイトバランスおよび階調の値である。
 画像処理制御部39は、画像処理装置30を構成する各部を統括的に制御する。画像処理制御部39は、CPU(Central Processing Unit)等を用いて構成される。画像処理制御部39は、画像処理装置30を構成する各部の指示やデータ等の転送を制御する。
 〔表示装置の構成〕
 次に、表示装置40の構成について説明する。表示装置40は、画像処理装置30から入力される画像データに対応する画像を表示する。表示装置40は、液晶または有機EL(Electro Luminescence)等の表示パネル等を用いて構成される。
 以上の構成を有する撮像システム1は、画像処理装置30が撮像素子105において発生するRTSノイズを補正し、表示装置40が画像処理装置30によって画像処理が施された画像データに対応する画像を表示する。
 〔RTSノイズの発生原因と特性〕
 次に、RTSノイズの発生原因とRTSノイズの特性について説明する。
 図3は、撮像素子105に光が当たらないように遮光した場合において、RTSノイズが発生するときに、アンプ部105eから出力されるアンプ出力の変動の例を示す図である。図4は、RTSノイズが発生するアンプ部105eを用いて読み出された画素値の分布の例を示す図である。
 RTSノイズは、アンプ部105eにおけるゲート酸化膜にトラップ準位が存在した場合、ランダムなタイミングで、このトラップ準位に電荷が捕獲されたり、放出されたりすることで発生する。このため、図3に示すように、RTSノイズが発生するアンプ部105eでは、アンプ出力が約Vrtsの範囲でランダムに変動する。また、電位の変動は、一瞬で起こらず、わずかな時間τを要する。
 一般に、撮像素子105では、画素105aから読み出した画素値からノイズを低減するため、相関二重サンプリング処理(以下、「CDS処理」という)が行われる。CDS処理では、撮像制御部114が撮像素子105のリセットスイッチ(図示せず)をオンにして、FD部105dの電荷をリセットさせ、さらに、撮像制御部114が第2スイッチ105fをオンにして、リセット状態を作り、リセット状態の信号(基準信号)を読み出す(出力する)。次に、CDS処理では、撮像制御部114が第1スイッチ105b(または第1スイッチ105b(1)~105b(8)のいずれか)のみをオンにして、画素105aで発生した電荷をFD部105dに転送し、さらに第2スイッチ105fをオンにした読み出し状態(出力状態)を作り、読み出し状態の信号を読み出す(出力する)。続いて、CDS処理では、読み出し状態の信号からリセット状態の信号(基準信号)を減算することで得られる信号を画素値として変換する。
 図3に示すように、撮像素子105は、CDS処理により、時間tr1(リセット状態)および時間ts1(読み出し状態)それぞれの信号を読み出すと、時間tr1および時間ts1それぞれのアンプ出力Vがほぼ同様であるため、主にランダムノイズによる影響を受け、読み出された画素値が図4に示す分布Aのような0を中心とした分布となる。同様に、撮像素子105は、時間tr2(リセット状態)と時間ts2(読み出し状態)でも、時間tr2および時間ts2それぞれのアンプ出力Vがほぼ同様であるため、読み出された画素値が図4に示す分布Aのようになる。
 一方、撮像素子105は、CDS処理により、時間tr3(リセット状態)および時間ts3(読み出し状態)それぞれの信号を読み出すと、時間tr3のアンプ出力と比べ時間ts3のアンプ出力が約Vrts低いため、2つの信号の差をとると、アンプ出力の変化量であるVrtsに対応する画素値であるRTS_Value分マイナス方向にシフトし、読み出された画素値が-RTS_Valueを中心とした分布Bとなる。
 これに対して、撮像素子105は、CDS処理により、時間tr4(リセット状態)および時間ts4(読み出し状態)それぞれの信号を読み出すと、時間tr4のアンプ出力に比べて時間ts4のアンプ出力が約Vrts高いため、2つの信号の差をとるとアンプ出力の変化量であるVrtsに対応する画素値であるRTS_Value分プラス方向にシフトし、読み出された画素値がRTS_Valueを中心とした分布Cとなる。
 ここで、図3のアンプ出力の変動は、時間τを要して生じるため、電位が変動している途中で信号を読み出す場合もある。この場合、リセット状態の読み出し時間および読み出し状態の読み出し時間の間で、アンプ出力差が-Vrtsより大きく、Vrtsより小さい。この結果、撮像素子105から読み出された画素値も、-RTS_Valueより大きく、RTS_Valueより小さな値となる。時間τは、撮像素子105の条件(例えば温度や駆動電圧等)が一定であれば、ほぼ一定になると考えられるため、-RTS_Valueより大きくRTS_Valueより小さな画素値が同様の確率で発生する。ここでは、これらの画素値の発生頻度をαnoiseと定義する。また、分布Bおよび分布Cの各々は、中央値のみ異なるが、それ以外は同様の分布となる。このため、以下においては、分布Aに対する分布Bまたは分布Cの割合をαrtsと定義する。このαrtsは、アンプ部105eのアンプ出力の変動周期が短いほど、大きくなる。
 このように、CDS処理によりRTSノイズが発生するアンプ部105eを用いて読み出された画素値は、図4のような分布となる。なお、撮像素子105に光が当たっている条件では、読み出し状態の電位が露光量に応じて変化する。しかしながら、RTSノイズによる電位の変化は、露光量によらず一定である。即ち、RTSノイズは、露光量に依存せず、-RTS_Value以上、RTS_Value以下の範囲で正常な画素値に対してランダムに変動する特性を有する。なお、図4において、分布A、分布B、分布Cを模式的に示したが、一般には正規分布となる。
 また、RTSノイズは、読み出し回路(アンプ部105e)に起因するノイズであるため、図2に示すように、複数の画素105aの各々が1つの読み出し回路を共有している場合、全ての共有画素(画素105a(1)~105a(8))において同様の特性のRTSノイズが発生する。
 また、図2に示した読み出し回路(アンプ部105e)以外にも、撮像素子105の列方向で共有しているカラムアンプやソースフォロア等においても、RTSノイズが発生する場合がある。この場合、同じカラムアンプおよびソースフォロアを共有する列方向の全ての画素においても同様の特性のRTSノイズが発生する。本実施の形態では、読み出し回路(アンプ部105e)以外の回路で発生したRTSノイズにも適用することができる。
 このようにRTSノイズは、被写体を固定して同じ条件で撮影した場合、撮影により得られた画像の画素値が一定範囲内(-RTS_Value以上、RTS_Value以下)で振幅(変動)するような点滅欠陥ノイズの一種となる。
 〔画像処理装置の処理〕
 次に、画像処理装置30が実行する処理について説明する。図5は、画像処理装置30が実行する処理の概要を示すフローチャートであり、画像処理装置30が実行するメインルーチンのフローチャートである。
 図5に示すように、まず、評価値算出部351は、補間処理部32から出力された現画像(現フレーム)と、フレームメモリ33に記憶された過去画像(過去フレーム)と、を用いて差分平均値mSAD(フレーム間差分値)を算出する(ステップS1)。ここで、過去画像とは、補間処理部32から入力される現画像より1フレーム前の時点において、ノイズ低減部35から出力されたNR処理済みの画像である。具体的には、評価値算出部351は、NR処理の対象となる画素である注目画素の座標を(x,y)とした場合、差分平均値mSADを以下の式(1)を用いて算出し、式(1)のSAD(m,n)が最小となる(m,n)の値をノイズ低減処理部355および判定部354へ出力する。なお、以下においては、G信号についての処理を説明するが、R信号およびB信号についても同一の処理を行う。
Figure JPOXMLDOC01-appb-M000001
 ここで、上式(1)において、min{}は、括弧内の値の最小値を取得する処理を表す。また、m=-1、0、1、n=-1、0、1の場合を示すが、これに限定されず、m=-10、-9、・・・、9、10、n=-10、-9、・・・、9、10などより広い範囲を設定してもよい。また、FG_cur(x,y)は、現画像(RGB画像)の座標(x,y)におけるG信号値であり、FG_pre(x,y)は、過去画像の座標(x,y)におけるG信号値である。また、kは、自然数であり、(2k+1)は、差分平均値mSADを算出する際のカーネルサイズに相当する。なお、kは、予め一定の値を設定しておくこともできるし、操作部37を介してユーザが任意の値を設定する構成としてもよい。また、|A|は、実数Aの絶対値を取得する処理を表す。
 また、上式(1)で説明した(m,n)は、フレーム間の動きベクトルに相当し、この動きベクトルの探索範囲は、m=-1、0、1、n=-1、0、1の場合±1画素である。上式(1)では、探索範囲の中で最も小さいSADをmSADとして選択する。
 続いて、推定ノイズ量取得部352は、ルックアップテーブル353を参照して、現画像における注目画素の画素値(信号値)に対応する推定ノイズ量Nを取得する(ステップS2)。
 図6は、ルックアップテーブル353が記憶するノイズモデルの一例を示す図である。図6において、縦軸がノイズ量を示し、横軸が画素値を示す。なお、図6においては、縦軸のノイズ量として画素値の標準偏差を用い、撮像素子105の特性に応じたノイズモデルを示す。
 図6の曲線Lx1に示すように、撮像素子105におけるノイズ量は、画素値が大きくなるに従って増加する。このため、本実施の形態1では、推定ノイズ量取得部352は、図6の曲線Lx1のノイズモデルを参照して、現画像の注目画素の画素値に対応する推定ノイズ量Nを取得する。なお、図6に示した曲線以外にも、ノイズモデルを近似式や折れ線で近似した特性であってもよい。
 その後、判定部354は、評価値算出部351によって算出された差分平均値mSADが推定ノイズ量取得部352によって取得されたノイズ量N以下(mSAD≦N)であるか否かを判定することによって、静止状態であるか動作状態であるか否かの判別処理を行う(ステップS3)。判定部354が評価値算出部351によって算出された差分平均値mSADが推定ノイズ量取得部352によって取得されたノイズ量N以下であると判定した場合(ステップS3:Yes)、画像処理装置30は、後述するステップS4へ移行する。これに対して、判定部354が評価値算出部351によって算出された差分平均値mSADが推定ノイズ量取得部352によって取得されたノイズ量以下でないと判定した場合(ステップS3:No)、画像処理装置30は、後述するステップS5へ移行する。
 ステップS4において、ノイズ低減処理部355は、現画像における注目画素に対して時間方向NR処理を実行する。なお、時間方向NR処理の詳細は、後述する。ステップS4の後、画像処理装置30は、後述するステップS6へ移行する。
 ステップS5において、ノイズ低減処理部355は、現画像における注目画素に対して空間方向NR処理を実行する。なお、空間方向NR処理の詳細は、後述する。ステップS5の後、画像処理装置30は、後述するステップS6へ移行する。
 続いて、画像処理制御部39が現画像の全画素に対してNR処理が実行されたと判断した場合(ステップS6:Yes)、画像処理装置30は、後述するステップS7へ移行する。これに対して、画像処理制御部39が現画像の全画素に対してNR処理が実行されていないと判断した場合(ステップS6:No)、画像処理装置30は、上述したステップS1へ戻る。
 ステップS7において、ノイズ低減処理部355は、NR画像をフレームメモリ33および画像処理部38へ出力する。ステップS7の後、画像処理装置30は、本処理を終了する。
 判定部354は、評価値算出部351によって算出された差分平均値mSADがノイズ量N以下であるかを判定することによって、ノイズ低減処理部355が時間方向NR処理を行うか空間方向NR処理を行うか切替えていたが、ノイズ低減処理部355が時間方向NR処理と空間方向NR処理の両方を行い、差分平均値mSADの大きさ(または動きベクトルの確からしさを示す信頼度の大きさ)に応じて設定されるブレンド値に基づいて各NR処理の結果をブレンドするようにしてもよい。
 〔時間方向NR処理の概要〕
 次に、上述した図5のステップS4の時間方向NR処理について詳細に説明する。図7は、時間方向NR処理の概要を示すフローチャートである。
 図7に示すように、ノイズ低減処理部355は、記憶部36からRTSノイズ情報を取得し、このRTSノイズ情報に基づいて、現画像の注目画素においてRTSノイズが発生する可能性があるか否かを判断する(ステップS10)。ノイズ低減処理部355が現画像の注目画素においてRTSノイズが発生する可能性があると判断した場合(ステップS10:Yes)、画像処理装置30は、後述するステップS11へ移行する。これに対して、ノイズ低減処理部355が現画像の注目画素においてRTSノイズが発生する可能性がないと判断した場合(ステップS10:No)、画像処理装置30は、後述するステップS12へ移行する。
 ステップS11において、ノイズ低減処理部355は、記憶部36から取得したRTSノイズ情報に含まれるRTS_Valueに応じて、時間方向NR処理で用いる係数Co_we_preを設定する。ステップS11の後、画像処理装置30は、後述するステップS13へ移行する。
 図8は、ノイズ低減処理部355が設定するRTS_Valueと係数Co_we_preとの関係を示す図である。図8において、横軸がRTS_Valueの大きさを示し、縦軸が係数Co_we_preを示す。
 図8に示すように、ノイズ低減処理部355は、直線L1、曲線L2および折れ線L3に示すように、RTS_Valueの大きさに応じて、係数Co_we_preを設定する。例えば、ノイズ低減処理部355は、RTS_Valueの大きさに応じて、直線L1に示すように直線的に大きくなるように係数Co_we_preを設定してもよいし、曲線L2に示すように指数的に大きくなるように係数Co_we_preを設定してもよいし、折れ線L3に示すように段階的に大きくなるように係数Co_we_preを設定してもよい。このように、ノイズ低減処理部355は、注目画素にRTSノイズが発生する可能性がある場合、注目画素にRTSノイズが発生しない場合と比べて、時間方向NR処理で用いる係数Co_we_preを大きくする。
 図7に戻り、ステップS12以降の説明を続ける。
 ステップS12において、ノイズ低減処理部355は、係数Co_we_preを1(係数Co_we_pre=1)に設定する。ステップS12の後、画像処理装置30は、後述するステップS13へ移行する。
 続いて、ノイズ低減処理部355は、上述したステップS11またはステップS12で設定した係数Co_we_preを用いて、重み係数we_preを算出する(ステップS13)。具体的には、ノイズ低減処理部355は、以下の式(2)によって、重み係数we_preを算出する。
 we_per=Co_we_pre×we_pre   ・・・(2)
 その後、ノイズ低減処理部355は、上述したステップS13で算出した重み係数we_preを用いて、時間方向NR処理を実行し(ステップS14)、画像処理装置30は、図5のメインルーチンへ戻る。具体的には、ノイズ低減処理部355は、以下の式(3)を用いて時間方向NR処理を実行する。
Figure JPOXMLDOC01-appb-M000002
 ここで、上述(3)において、FG_NR(x,y)は、NR画像の座標(x,y)における現画像の画素値である。また、we_curおよびwe_preは、加重平均処理時の重み係数である。なお、本実施の形態1では、we_preをwe_curと比較して大きくすることによって、ノイズ低減量を大きくしてもよい。また、本実施の形態1では、we_preを予め一定の値に設定しておいてもよいし、操作部37を介してユーザが任意の値を設定可能な構成としてもよい。
 〔空間方向NR処理の概要〕
 次に、図5のステップS5において説明した空間方向NR処理について詳細に説明する。図9は、空間方向NR処理の概要を示すフローチャートである。
 図9に示すように、ノイズ低減処理部355は、記憶部36からRTSノイズ情報を取得し、このRTSノイズ情報に基づいて、現画像の注目画素においてRTSノイズが発生する可能性があるか否かを判断する(ステップS21)。ノイズ低減処理部355が現画像の注目画素においてRTSノイズが発生する可能性があると判断した場合(ステップS21:Yes)、画像処理装置30は、後述するステップS22へ移行する。これに対して、ノイズ低減処理部355が現画像の注目画素においてRTSノイズが発生する可能性がないと判断した場合(ステップS21:No)、画像処理装置30は、後述するステップS22へ移行する。
 ステップS21において、ノイズ低減処理部355は、記憶部36から取得したRTSノイズ情報に含まれるRTS_Valueに応じて、空間方向NR処理で用いる係数Co_Cbを設定する。ステップS21の後、画像処理装置30は、後述するステップS23へ移行する。
 図10は、ノイズ低減処理部355が設定するRTS_Valueと係数Co_Cbとの関係を示す図である。図10において、横軸がRTS_Valueの大きさを示し、縦軸が係数Co_Cbを示す。
 図10に示すように、ノイズ低減処理部355は、直線L11、曲線L12および折れ線L13に示すように、RTS_Valueの大きさに応じて、係数Co_Cbを設定する。例えば、ノイズ低減処理部355は、RTS_Valueの大きさに応じて、直線L11に示すように直線的に大きくなるように係数Co_Cbを設定してもよいし、曲線L12に示すように指数的に大きくなるように係数Co_Cbを設定してもよいし、折れ線L13に示すように段階的に大きくなるように係数Co_Cbを設定してもよい。このように、ノイズ低減処理部355は、注目画素にRTSノイズが発生する可能性がある場合、注目画素にRTSノイズが発生しない場合と比べて、空間方向NR処理で用いる係数Co_Cbを大きくする。
 図9に戻り、ステップS22以降の説明を続ける。
 ステップS22において、ノイズ低減処理部355は、係数Co_Cbを1(係数Co_Cb=1)に設定する。ステップS22の後、画像処理装置30は、後述するステップ23へ移行する。
 続いて、ノイズ低減処理部355は、上述したステップS21またはステップS22で設定した係数Co_Cbを用いて、重み係数Cbを算出する(ステップS24)。具体的には、ノイズ低減処理部355は、以下の式(4)によって、重み係数Cbを算出する。
 Cb=Co_Cb×Cb   ・・・(4)
 その後、ノイズ低減処理部355は、上述したステップS23で算出した重み係数Cbを用いて、空間方向NR処理を実行し(ステップS24)、画像処理装置30は、図5のメインルーチンへ戻る。具体的には、ノイズ低減処理部355は、以下の式(5)を用いて空間方向NR処理を実行する。
Figure JPOXMLDOC01-appb-M000003
 ここで、上式(5)において、we_diff_cur(x+i,y+j)およびwe_diff_pre(x+i,y+j)は、加重平均処理時の重み係数に相当する。この係数は、以下の式(6)に示すように、ガウス分布によって与えられる。また、iおよびjは、自然数である。さらに、mおよびnは、上述した式(1)において差分平均値mSADとして選択されたSAD(m,n)である。
Figure JPOXMLDOC01-appb-M000004
 上述した式(6)に示すように、本実施の形態1で用いられる空間方向NR処理では、現画像の注目画素の画素値と、この注目画素の周辺画素の画素値との差分の大きさに応じて、適応的に重み係数が設定される。具体的には、現画像の注目画素の画素値と、この注目画素の周辺画素の画素値との差分が大きい場合、加重平均処理時の重みが小さくなる。従って、エッジ部等の画素値が急に変化する領域の画素は、加重平均処理に寄与しなくなるため、エッジ部を保持してノイズ成分のみを低減することができる。しかしながら、この空間方向NR処理では、現画像の注目画素の画素値と、この注目画素の周辺画素の画素値との差分に応じて重み係数が制御されているため、ノイズ低減の程度(平滑化の強度)が画像に含まれるノイズ量に依存する。具体的には、ノイズが大きいほど差分が大きくなるため、重み係数が小さくなり、上式(3)の加重平均における寄与が小さくなる。そのため、ノイズが大きいほどノイズ低減の程度は、弱くなり、ノイズが低減されなくなってしまう。
 そこで、本実施の形態1では、ノイズ低減処理部355は、上式(6)のガウス分布の標準偏差σを、図5のステップS2で出力されるノイズ量Nに基づいて算出する。具体的には、ノイズ低減処理部355は、以下の式(7)を用いて標準偏差σを算出する。
 σ=Cb×N   ・・・(7)
 ここで、上式(7)において、重み係数Cbは、正の実数である。なお、重み係数Cbは、予め一定の値を設定してもよいし、操作部37を介してユーザが任意の値に設定可能な構成としてもよい。
 このように、ノイズ低減処理部355は、ノイズ量Nに基づいて、ガウス分布の標準偏差σを算出することによって、ノイズ量に対して適応的なノイズ低減処理を実行することができる。即ち、ノイズ低減処理部355は、ノイズ量が大きいほど、標準偏差σが大きくなるため、上式(6)において差分(例えばFG_cur(x+i,y+j)-FG_cur(x+,y))がノイズにより大きくなった場合でも、標準偏差σが推定ノイズ量Nに依存しない場合に比べて、重み係数を大きくできる。このため、ノイズが大きい場合でも、平滑化の強度を維持することができる。
 以上説明した本発明の実施の形態1によれば、ノイズ低減部35が第3外部I/F部31から入力されたRTSノイズ情報に基づいて、現画像の注目画素においてRTSノイズが発生する可能性がある場合、注目画素の画素値に対して、NR処理の係数をRTSノイズが発生しない画素に対して行うNR処理の係数と比べて大きくすることによって、注目画素の画素値と周辺画素の画素値との相関が低い場合であっても、動画に発生する可能性があるRTSノイズを低減することができるので、画質を向上させることができる。
 さらに、本発明の実施の形態1によれば、現画像の注目画素においてRTSノイズが発生する可能性がある場合、注目画素の画素値に対して、従来の撮像装置等に備えられたNR機能の係数を大きくするだけでよいため、簡易な構成でRTSノイズを低減することができる。
 なお、本発明の実施の形態1では、動き検出部34が動きベクトル(m、n)を算出するためにSADを用いたが、その他SSD(Sum of Squared Difference)やNCC(Normalized Cross-Correlation)など既存の方法によって算出してもよい。
(実施の形態2)
 次に、本発明の実施の形態2について説明する。本実施の形態2に係る撮像システムは、上述した実施の形態1に係る画像処理装置30と構成が異なるうえ、本実施の形態2に係る画像処理装置が実行する処理も異なる。具体的には、本実施の形態2に係る画像処理装置は、現画像から算出した代表値および過去画像で算出した代表値に基づいて、注目画素を補正する際の代表値を算出し、この代表値に基づいて注目画素の画素値を補正する。以下においては、本実施の形態2に係る撮像システムの構成を説明後、本実施の形態2に係る画像処理装置が実行する処理について説明する。なお、上述した実施の形態1に係る撮像システム1と同一の構成には同一の符号を付して説明を省略する。
 〔撮像システムの構成〕
 図11は、本発明の実施の形態2に係る撮像システムの構成を模式的に示すブロック図である。図11に示す撮像システム1aは、上述した実施の形態1に係る撮像システム1の画像処理装置30に換えて、画像処理装置30aを備える。
 〔画像処理装置の構成〕
 画像処理装置30aは、上述した実施の形態1に係る画像処理装置30から補間処理部32、フレームメモリ33および動き検出部34の構成を省略するとともに、ノイズ低減部35に換えてノイズ低減部35aを備える。
 ノイズ低減部35aは、第3外部I/F部31から入力された現画像に対して、RTSノイズを補正して画像処理部38へ出力する。ノイズ低減部35aは、RTSノイズ補正部320を有する。
 RTSノイズ補正部320は、第3外部I/F部31から入力された現画像に対して、RTSノイズを補正するRTSノイズ補正処理を行い、このRTSノイズ補正を行った現画像を画像処理部38へ出力する。
 ここで、RTSノイズ補正部320の詳細な構成について説明する。図12は、RTSノイズ補正部320の詳細な構成を示すブロック図である。
 図12に示すように、RTSノイズ補正部320は、RTSノイズ画素判定部321と、候補値算出部322と、代表値算出部323と、ランダムノイズ量推定部324と、補正値算出部325と、を有する。
 RTSノイズ画素判定部321は、第3外部I/F部31から入力されたRTSノイズ情報に基づいて、第3外部I/F部31から入力された現画像の注目画素においてRTSノイズが発生する可能性があるか否かを判定し、判定結果を候補値算出部322および代表値算出部323へ出力する。具体的には、RTSノイズ画素判定部321に対して画素が位置を入力されると、その画素(注目画素)に対応するRTSノイズ情報があるか否かを判定し、RTSノイズ情報がある場合、このRTSノイズ情報(RTSノイズが有りを示す情報)を出力する一方、RTSノイズ情報がない場合、RTSノイズが発生しない画素と見なし、RTSノイズ情報を出力しない。
 候補値算出部322は、注目画素の画素値と、RTSノイズ画素判定部321の判定結果とに基づいて、RTSノイズ画素判定部321によって注目画素においてRTSノイズが発生する可能性があると判定されている場合、注目画素の画素値に対する補正量の候補値を複数算出し、注目画素の画素値と、算出した複数の候補値を代表値算出部323、ランダムノイズ量推定部324および補正値算出部325それぞれへ出力する。
 代表値算出部323は、RTSノイズ画素判定部321によって注目画素においてRTSノイズが発生する可能性があると判定されている場合には、注目画素における周囲の少なくともRTSノイズ画素判定部321によってRTSノイズが発生しないと判定されている画素と、後述するランダムノイズ量推定部324が算出した注目画素に対応するランダムノイズ量とに基づいて、RTSノイズが発生しない場合の画素値に相当する代表値を算出する。また、代表値算出部323は、第3外部I/F部31から入力された現画像と記憶部36に記憶された過去画像とに基づいて、現画像の代表値を算出して補正値算出部325へ出力する。
 ここで、代表値算出部323の詳細な構成について説明する。図13は、代表値算出部323の詳細な構成を示すブロック図である。
 図13に示すように、代表値算出部323は、第1参照値算出部323aと、第1ランダムノイズ量推定部323bと、許容範囲算出部323cと、代表値決定部323dと、を有する。
 第1参照値算出部323aは、ランダムノイズ量を推定する際に用いる参照値を算出する。具体的には、第1参照値算出部323aは、注目画素の画素値にRTS_Valueを加算した値、または注目画素の画素値に補正量の候補値のうち最大値を加算した値を参照値として算出する。
 第1ランダムノイズ量推定部323bは、第3外部I/F部31から入力されたランダムノイズモデルと第1参照値算出部323aが算出した参照値とに基づいて、ランダムノイズ量を算出する。
 許容範囲算出部323cは、注目画素を基準に設定した算出範囲内にある画素の画素値と第1ランダムノイズ量推定部323bで算出したランダムノイズ量とに基づいて、代表値算出処理に使用可能な画素値の範囲である許容範囲(有効範囲)を算出する。
 代表値決定部323dは、算出範囲内における注目画素以外の画素値に対して、許容範囲内(有効範囲内)に含まれている画素値(カラーフィルタを用いた撮像素子105の場合、注目画素と同色の画素値)を選択し、選択した画素値が所定の値以上の場合、この選択した画素値の中央値を代表値として決定する。また、代表値決定部323dは、現画像と記憶部36に記憶された過去画像とに基づいて、現画像の代表値を決定する。
 図12に戻り、RTSノイズ補正部320の構成の説明を続ける。
 ランダムノイズ量推定部324は、第3外部I/F部31から入力されたランダムノイズモデルに基づいて、画素値に対応するランダムノイズ量を推定し、推定結果を候補値算出部322、代表値算出部323および補正値算出部325それぞれへ出力する。即ち、ランダムノイズ量推定部324に対して画素値を入力すると、その画素値に対応するランダムノイズ量が出力される。
 補正値算出部325は、RTSノイズ画素判定部321によって注目画素においてRTSノイズが発生する可能性がある画素と判定されている場合、候補値算出部322が算出した複数の候補値に基づいて、注目画素の画素値を補正する。具体的には、補正値算出部325は、注目画素の画素値と、候補値算出部322によって算出された複数の候補値と、代表値算出部323によって算出された代表値と、に基づいて、RTSノイズを補正した画素値を算出して画像処理部38へ出力する。より具体的には、補正値算出部325は、候補値算出部322が算出した複数の候補値の中から、代表値算出部323が算出した代表値に補正結果が最も近くなるような候補値に基づいて、注目画素の画素値を補正する。これに対して、補正値算出部325は、RTSノイズ画素判定部321によって注目画素においてRTSノイズが発生しない画素と判定されている場合、注目画素の画素値をそのまま出力する。なお、本実施の形態2では、補正値算出部325が補正部として機能する。
 ここで、補正値算出部325の詳細な構成について説明する。図14は、補正値算出部325の詳細な構成を示すブロック図である。
 図14に示すように、補正値算出部325は、第2参照値算出部325aと、第2ランダムノイズ量推定部325bと、補正量決定部325cと、画素値補正部325dと、を有する。
 第2参照値算出部325aは、注目画素の画素値からRTS_Valueを減算した値または注目画素の画素値から補正量の候補値のうち最大値を減算した値を参照値として算出する。
 第2ランダムノイズ量推定部325bは、第3外部I/F部31から入力されたランダムノイズモデルと第2参照値算出部325aが算出した参照値とに基づいて、ランダムノイズ量を算出する。
 補正量決定部325cは、第2ランダムノイズ量推定部325bが算出したランダムノイズ量に、一定の係数を乗じた値を閾値として算出し、候補値算出部322によって算出された補正量の候補値の中における最大候補値が閾値以上であるか否かを判定後、代表値算出部323が算出した代表値と、注目画素の絶対値に最も近い補正量の候補値を、補正量として選択する。
 画素値補正部325dは、代表値に近づくように、注目画素の画素値に対して補正量分加算または減算することによって、補正した注目画素の画素値を算出する。
 〔画像処理装置の処理〕
 次に、画像処理装置30aが実行する処理について説明する。図15は、画像処理装置30aが実行する処理の概要を示すフローチャートであり、画像処理装置30aが実行するメインルーチンのフローチャートである。
 図15に示すように、まず、RTSノイズ補正部320は、後述するステップS102~ステップS105の処理を順次行うための注目画素を設定する(ステップS101)。なお、RTSノイズ補正部320は、画素毎に左上から右下に向かってラスタ順に、0より大きい整数を1、2、3、・・・と順にインデックスとして割り当てる。次に、RTSノイズ補正部320は、ステップS101が実行される毎に、カウンタを1ずつ増加させる(図15の処理が開始される時点でカウンタは0にリセットする)。RTSノイズ補正部320は、カウンタが示しているインデックスが割り当てられている画素を注目画素として設定する。即ち、RTSノイズ補正部320によってステップS101が最初に実行されると、RTSノイズ補正部320が0にリセットされたカウンタを1増加させるため、カウンタが1を示し、左上の画素が注目画素となる。RTSノイズ補正部320がステップS101の処理を2回(2回目)実行すると、カウンタが2を示すため、左上の画素の右側の画素が注目画素となる。
 続いて、RTSノイズ画素判定部321は、第3外部I/F部31、第1外部I/F部115およびバス113を介して、撮像装置10の不揮発メモリ112のRTSノイズ情報記録部112bが記録するRTSノイズ情報を取得し、取得したRTSノイズ情報に基づいて、注目画素においてRTSノイズが発生する可能性があるか否かを判定する(ステップS102)。即ち、RTSノイズ画素判定部321は、注目画素が含まれる共有画素ブロックの位置情報が、RTSノイズ情報に含まれているか否かを判定する。具体的には、RTSノイズ画素判定部321は、注目画素が含まれる共有画素ブロックの位置情報が、RTSノイズが発生する可能性がある共有画素ブロックとしてRTSノイズ情報に含まれているか否かを判定する。RTSノイズ画素判定部321によって注目画素においてRTSノイズが発生する可能性があると判定(注目画素が含まれる共有画素ブロックの位置情報が、RTSノイズ情報に含まれていると判定)された場合(ステップS102:Yes)、画像処理装置30aは、後述するステップS103へ移行する。これに対して、RTSノイズ画素判定部321によって注目画素においてRTSノイズが発生する可能性がないと判定(注目画素が含まれる共有画素ブロックの位置情報が、RTSノイズ情報に含まれていないと判定)された場合(ステップS102:No)、画像処理装置30aは、後述するステップS106へ移行する。この場合において、RTSノイズ画素判定部321は、注目画素においてRTSノイズが発生する可能性がないと判定した場合、この注目画素の画素値をそのまま補正後の画素値として代表値算出部323へ出力する。
 ステップS103において、候補値算出部322は、RTSノイズを補正するための補正量の候補値を複数算出する。具体的には、候補値算出部322は、注目画素に対応するRTS_Value(RTSノイズ画素判定部321から出力されるRTSノイズ情報に含まれている)に基づいて、0以上RTS_Value以下の画素値として取り得る値全て(RAW画像として整数のみを取り得る場合、0以上RTS_Value以下の全ての整数)を候補値とする。なお、候補値算出部322は、撮像制御部114によって撮像素子105のカラムアンプ等に設定されたアンプゲイン値が、RTSノイズ検出時(アンプゲイン値=G0とする)と、RTSノイズ補正時(アンプゲイン値=G1とする)とで異なる場合、RTS_Valueを、RTSノイズ補正時のアンプゲイン値とRTSノイズ検出時のアンプゲイン値との比(G=G1/G0)に対してRTS_Valueを乗算した値に置き換えてもよい。また、候補値算出部322は、予め設定しうるアンプゲイン値毎のRTS_ValueをRTSノイズ情報に持たせ、この設定しているアンプゲイン値に応じたRTS_Valueを用いてもよい。
 続いて、代表値算出部323は、記憶部36が記憶する過去画像と現画像の注目画素の周辺画素の画素値に基づいて、代表値(注目画素において、RTSノイズが発生していない場合における予測される画素値)を算出する代表値算出処理を実行する(ステップS104)。なお、代表値算出部323は、RTSノイズが発生する画素を含めるようにしてもよい。
 〔代表値算出処理の概要〕
 図16は、図15のステップS104の代表値算出処理の概要を示すフローチャートである。
 図16に示すように、まず、代表値算出部323は、注目画素を基準に、代表値算出の対象とする最小の算出範囲を設定する(ステップS201)。具体的には、代表値算出部323は、例えば注目画素を中心として対象の範囲で最大7×7の範囲を算出範囲とする場合、7×7以下の最小の範囲である3×3を最小の算出範囲として設定する。
 続いて、第1参照値算出部323aは、注目画素の画素値に、RTS_Valueを加算した値または注目画素の画素値に補正量候補値の最大値を加算した値を参照値として算出する(ステップS202)。
 その後、第1ランダムノイズ量推定部323bは、第3外部I/F部31から入力されたランダムノイズモデルとステップS202で第1参照値算出部323aが算出した参照値とに基づいて、ランダムノイズ量を算出する(ステップS203)。本実施の形態2における第1ランダムノイズ量推定部323bは、上述した図6の曲線LX1のノイズモデル(ランダムノイズモデル)と第1参照値算出部323aが算出した参照値とに基づいて、ランダムノイズ量を算出する(標準偏差を算出する)。なお、図6に示した曲線以外にも、ランダムノイズモデルを近似式や折れ線で近似した特性であってもよい。
 ステップS203の後、許容範囲算出部323cは、算出範囲内にある画素の画素値に基づいて、代表値算出処理に使用可能な画素値の範囲である許容範囲(有効範囲)を算出する(ステップS204)。具体的には、許容範囲算出部323cは、許容範囲(有効範囲)の上限を以下の式(8)によって算出する。
 参照値+ランダムノイズ量(標準偏差)×R+RTS_Value
                           ・・・(8)
 ここで、Rは、所定の係数であり、ランダムノイズに対してRTSノイズが視覚的にどの程度把握できるかに応じて設定する。例えばRの係数としては、2前後の値が好ましい。また、代表値算出部323は、許容範囲の下限を以下の式(9)によって算出する。
 参照値-ランダムノイズ量(標準偏差)×R-RTS_Value
                           ・・・(9)
 なお、RTS_Valueに代えて、複数の候補値の最大値を用いてもよい。また、式(8)および式(9)における参照値は、上述したステップS203において第1ランダムノイズ量推定部323bによるランダムノイズ量を推定するために使用された参照値と異なる参照値方法により得られた参照値としてもよい。このように、許容範囲算出部323cは、注目画素のRTSノイズと、この注目画素周辺のランダムノイズを考慮した許容範囲を算出することができる。
 その後、許容範囲算出部323cは、算出範囲内において注目画素以外の画素の画素値(カラーフィルタを用いた撮像素子105の場合、注目画素と同色の画素値)の各々に対して、上述したステップS204で算出した許容範囲内であるか否かを判定し、この許容範囲内の画素値の個数をカウントする(ステップS205)。このステップS205で得られるカウント値は、平坦な被写体の場合ほど大きく、エッジを含む被写体の場合ほど小さくなる傾向がある。なお、算出範囲内においてRTSノイズが発生する可能性がある画素は、カウントしないようにしても良い。
 続いて、上述したステップS205でカウントしたカウント値(画素数)が所定の値ThRefより大きい場合(ステップS206:Yes)、画像処理装置30aは、後述するステップS209へ移行する。ここで、所定の値ThRefは、代表値算出部323が注目画素の周辺画素から代表値を算出するため、1以上とすることが好ましい。これに対して、上述したステップS205でカウントしたカウント値が所定の値ThRefより大きくない場合(ステップS206:No)、画像処理装置30aは、後述するステップS207へ移行する。
 ステップS207において、代表値算出の対象とする算出範囲が最大である場合(ステップS207:Yes)、画像処理装置30aは、後述するステップS209へ移行する。これに対して、代表値算出の対象とする算出範囲が最大でない場合(ステップS207:No)、画像処理装置30aは、後述するステップS208へ移行する。
 ステップS208において、許容範囲算出部323cは、代表値を算出する算出範囲を拡大する(ステップS208)。具体的には、許容範囲算出部323cは、代表値算出の対象とする算出範囲を最大の範囲内に収まる範囲で、水平または垂直方向に1画素以上拡大する。例えば、許容範囲算出部323cは、算出範囲として注目画素を中心とした3×3の範囲を設定している場合、注目画素を中心とした5×5の範囲を算出範囲に設定し直す。ステップS208の後、画像処理装置30aは、ステップS202へ戻る。なお、ステップS208においては、許容範囲算出部323cは、3×3または5×5の範囲を算出範囲に設定していたが、例えば、水平または垂直だけ拡大し、5×3や3×5の範囲を算出範囲に設定するようにしても良い。
 ステップS209において、代表値決定部323dは、現画像の代表値を算出する。具体的には、代表値決定部323dは、まず、算出範囲内における注目画素以外の画素の画素値に対して、許容範囲内(有効範囲内)に含まれている画素値(カラーフィルタを用いた撮像素子105の場合、注目画素と同色の画素値)を選択する。その後、代表値決定部323dは、選択した画素数が所定の値ThRef以上の場合、この選択した画素値の中央値を現画像の代表値として算出する。なお、代表値決定部323dは、選択した画素値の数が偶数の場合、注目画素の画素値に近い側の中央値を現画像の代表値として算出する。この場合、過補正を防止することができる。また、代表値決定部323dは、選択した画素数が所定の値ThRef未満の場合、注目画素の画素値に最も近い画素値を有する算出範囲内における注目画素以外の画素の画素値を現画像の代表値として算出する。なお、代表値決定部323dは、中央値を用いて代表値を算出していたが、例えば平均や分布の中間値等の他の方法で算出するようにしてもよい。また、代表値決定部323dは、算出範囲内のエッジ方向判別を行い、このエッジ方向判別の結果に基づいて、最も相関の高い方向の周辺画素値を現画像の代表値として算出してもよい。さらに、代表値決定部323dは、注目画素以外の算出範囲内の画素においてRTSノイズが発生する可能性がある画素は除外するようにしても良い。このとき、代表値決定部323dは、ステップS209が実行される時点での算出範囲内において、RTSノイズが発生する可能性がない画素が全くない場合には、注目画素の現画像における画素値を現画像の代表値とする。
 続いて、現画像が代表値を算出する最初の画像データである場合(ステップS210:Yes)、代表値決定部323dは、上述したステップS209で算出した現画像の代表値を代表値に決定する(ステップS211)。ステップS211の後、画像処理装置30aは、上述した図15のメインルーチンへ戻る。
 ステップS210において、現画像が代表値を算出する最初の画像データでない場合(ステップS210:No)、代表値決定部323dは、現画像の代表値と、記憶部36が記憶する過去画像の代表値との平均値を代表値として決定する(ステップS212)。この場合、代表値決定部323dは、記憶部36に複数の過去画像それぞれの代表値を用いて平均値を算出してもよい。また、代表値決定部323dは、現画像の代表値と過去画像の代表値との平均値を代表値として決定していたが、重み付け平均や中央値を代表値として決定してもよい。さらに、代表値決定部323dは、記憶部36に複数の過去画像が記憶されている場合、この複数の過去画像(例えば3枚以上)の各々の代表値と現画像の代表値とに基づいて代表値を決定してもよい。さらにまた、代表値決定部323dは、前の結果の一定範囲内で変動しないように防止するため、記憶部36が記憶する過去画像の代表値に基づいて、所定の範囲を設定し、この範囲内に現画像の代表値が入るように現画像の代表値を決定してもよい。ステップS212の後、画像処理装置30aは、上述した図15のメインルーチンへ戻る。
 このように、代表値算出部323は、上述した代表値算出処理において、現画像の代表値を算出した後、現画像が最初の画像データである場合(記憶部36に過去画像が記憶されていない場合)、現画像の代表値を代表値として算出する一方、現画像が最初の画像データでない場合(記憶部36に過去画像が記憶されている場合)、現画像の代表値と過去画像の代表値の平均値を代表値として算出する。
 なお、代表値算出部323は、記憶部36に現画像と、この現画像と時間的に連続した補正前の画像(以下「補正前の過去画像」)とを記憶し、補正前の過去画像の代表値と現画像の代表値から代表値を算出してもよい。
 図15に戻り、ステップS105以降の説明を続ける。
 ステップS105において、補正値算出部325は、上述したステップS103で候補値算出部322によって算出された複数の候補値と、上述したステップS104で代表値算出部323によって算出された代表値とに基づいて、注目画素におけるRTSノイズが補正された画素値を算出する補正値算出処理を実行する。ステップS105の後、画像処理装置30aは、後述するステップS106へ移行する。
 〔補正値算出処理の概要〕
 図17は、図15のステップS105の補正値算出処理の概要を示すフローチャートである。
 図17に示すように、まず、第2参照値算出部325aは、注目画素の画素値からRTS_Value(補正量候補値の最大値)を減算した値を参照値として算出する(ステップS301)。
 続いて、第2ランダムノイズ量推定部325bは、第3外部I/F部31、第1外部I/F部115およびバス113を介してランダムノイズモデル情報記録部112cに記録されているランダムノイズモデルを取得し、取得したランダムノイズモデルとステップS301で第2参照値算出部325aが算出した参照値とに基づいて、ランダムノイズ量を算出する(ステップS302)。
 その後、補正量決定部325cは、上述したステップS302で第2ランダムノイズ量推定部325bが算出したランダムノイズ量に、一定の係数を乗じた値を閾値として算出する(ステップS303)。ここで、閾値は、以下の式(10)によって算出する。
 ランダムノイズ量×Rm   ・・・(10)
Rmは、ランダムノイズに対してRTSノイズが視覚的にどの程度見えるかに応じて決定される。例えば、Rmの値は、2前後が好ましい。
 続いて、補正量決定部325cは、上述した図15のステップS103で候補値算出部322によって算出された補正量の候補値の中における最大候補値が上述したステップS303で算出した閾値以上であるか否かを判定する(ステップS304)。補正量決定部325cが上述した図15のステップS103で候補値算出部322によって算出された補正量の候補値の中における最大候補値が閾値以上であると判定した場合(ステップS304:Yes)、画像処理装置30aは、後述するステップS305へ移行する。これに対して、補正量決定部325cが上述した図15のステップS103で候補値算出部322によって算出された補正量の候補値の中における最大候補値が閾値以上でないと判定した場合(ステップS304:No)、画像処理装置30aは、後述するステップS306へ移行する。
 ステップS305において、補正量決定部325cは、上述した図15のステップS104で代表値算出部323が算出した代表値と注目画素の画素値との差の絶対値に最も近い補正量の候補値を、補正量として選択する。ステップS305の後、画像処理装置30aは、後述するステップS307へ移行する。
 ステップS306において、補正量決定部325cは、補正量として0を設定する。ステップS306の後、画像処理装置30aは、ステップS307へ移行する。
 続いて、画素値補正部325dは、代表値に近づくように、注目画素の画素値に対して補正量分加算または減算することによって、補正した注目画素の画素値を補正する(ステップS307)。ステップS307の後、画像処理装置30aは、図15のメインルーチンへ戻る。
 図15に戻り、ステップS106以降の説明を続ける。
 ステップS106において、RTSノイズ補正部320は、全画素に対して上述したステップS101~ステップS105の処理が終了したか否かを判定する。RTSノイズ補正部320が全画素に対して上述したステップS101~ステップS105の処理が終了したと判定した場合(ステップS106:Yes)、画像処理装置30aは、本処理を終了する。これに対して、RTSノイズ補正部320が全画素に対して上述したステップS101~ステップS105の処理が終了していないと判定した場合(ステップS106:No)、画像処理装置30aは、上述したステップS101へ戻る。
 以上説明した本発明の実施の形態2によれば、代表値算出部323が第3外部I/F部31から入力された現画像の代表値と記憶部36が記憶する過去画像の代表値に基づいて、RTSノイズが発生しない場合の画素値に相当する代表値を算出し、補正値算出部325が代表値算出部323によって算出された代表値に基づいて、現画像の注目画素の画素値を補正するので、注目画素の画素値が一定で、注目画素の周辺画素の画素値が現画像毎に大小変化する場合であっても、過去画像の代表値を含めて代表値を算出することで、補正後の注目画素の画素値の変動を軽減することができる。この結果、注目画素の画素値と周辺画素の画素値との相関が低い場合であっても、画質を向上させることができる。
(実施の形態3)
 次に、本発明の実施の形態3について説明する。本実施の形態3に係る撮像システムは、上述した実施の形態2に係る画像処理装置30aのノイズ低減部35aの構成が異なるうえ、本実施の形態3に係る画像処理装置が実行する代表値算出処理が異なる。具体的には、本実施の形態3は、現画像の注目画素と現画像の注目画素に対応する過去画像の画素の位置合わせを行った後に、RTSノイズが発生しない場合の画素値に相当する代表値を算出する。以下においては、本実施の形態3に係るノイズ低減部の構成を説明後、本実施の形態3に係る画像処理装置が実行する代表値算出処理について説明する。なお、上述した実施の形態1に係る撮像システム1と同一の構成には同一の符号を付して説明を省略する。
 図18は、本実施の形態3に係るノイズ低減部の詳細な構成を示すブロック図である。図18に示すノイズ低減部35bは、上述した実施の形態1に係る動き検出部34と、RTSノイズ補正部320bと、を備える。
 RTSノイズ補正部320bは、第3外部I/F部31から入力された現画像(RAW)に対してRTSノイズを補正するRTSノイズ補正処理を行い、この補正を行ったRAW画像を画像処理部38へ出力する。
 ここで、RTSノイズ補正部320bの詳細な構成について説明する。図19は、RTSノイズ補正部320bの詳細な構成を示すブロック図である。
 図19に示すように、RTSノイズ補正部320bは、上述した実施の形態2に係る代表値算出部323に換えて、代表値算出部326を有する。
 代表値算出部326は、動き検出部34が検出した動きベクトル、第3外部I/F部31から入力された現画像(RAW画像)の注目画素の画素値および記憶部36が記憶する現画像(RAW画像)の注目画素と対応する過去画像の注目画素の画素値に基づいて、RTSノイズが発生しない場合の画素値に相当する代表値を算出する。具体的には、代表値算出部326は、動き検出部34が検出した動きベクトルを用いて、現画像の注目画素と現画像の注目画素に対応する過去画像の画素の位置合わせを行った後に、現画像の注目画素の画素値と過去画像の注目画素の画素値とに基づいて、RTSノイズが発生しない場合の画素値に相当する代表値を算出する。
 〔代表値算出処理の概要〕
 次に、本実施の形態3に係る画像処理装置30aが実行する代表値算出処理について説明する。図20は、本実施の形態3に係る画像処理装置30aが実行する代表値算出処理の概要を示すフローチャートである。
 図20において、第3外部I/F部31から入力された現画像が代表値を算出する最初の画像データである場合(ステップS311:Yes)、画像処理装置30aは、後述するステップS312へ移行する。これに対して、第3外部I/F部31から入力された現画像が代表値を算出する最初の画像データでない場合(ステップS311:No)、画像処理装置30aは、後述するステップS321へ移行する。
 ステップS312~ステップS320は、上述した図16のステップS201~ステップS209にそれぞれ対応する。ステップS320の後、画像処理装置30aは、上述した図15のメインルーチンへ戻る。
 ステップS321において、代表値算出部326は、記憶部36に記憶された過去画像と現画像とに基づいて、RTSノイズが発生しない場合の画素値に相当する代表値を決定する。具体的には、代表値算出部326は、動き検出部34が検出した動きベクトルに基づいて、現画像の注目画素と現画像の注目画素の位置に対応する過去画像の画素の位置合わせを行った後に、現画像の注目画像の位置に相当する過去画像の画素値を代表値として決定する。なお、代表値算出部326は、動きベクトルとしてローカルベクトルを用いることが望ましい。また、代表値算出部326は、グローバルベクトルを用いてもよい。この場合、手ブレ等の画像全体の動きに対する位置合わせを行うことができる。ステップS321の後、画像処理装置30aは、上述した図15のメインルーチンへ戻る。
 以上説明した本発明の実施の形態3によれば、代表値算出部326は、動き検出部34が検出した動きベクトルに基づいて、現画像の注目画素と現画像の注目画素の位置に対応する過去画像の画素の位置合わせを行った後に、現画像の注目画像の位置に相当する過去画像の画素値を代表値として算出するので、現画像の注目画素の周辺画素が大小変化する場合であっても、この大小変化を軽減することができるので、画質を劣化させることなく、RTSノイズを低減することができる。
 また、本発明の実施の形態3によれば、注目画素の画素値のみが大小変化する場合であっても、画質を劣化させることなく、RTSノイズを低減することができる。
(実施の形態3の変形例)
 次に、本発明の実施の形態3の変形例について説明する。本実施の形態3の変形例は、ノイズ低減部の構成が異なるうえ、画像処理装置が実行する代表値算出処理が異なる。具体的には、本実施の形態3の変形例に係るノイズ低減部は、動きベクトルの信頼度を算出する信頼度算出部を有する。さらに、本実施の形態3の変形例では、現画像の注目画素の画素値および現画像の注目画素に対応する過去画像の画素の画素値を用いて、RTSノイズが発生しない場合の画素値に相当する代表値を算出する。以下においては、本実施の形態3の変形例に係るノイズ低減部の構成を説明後、本実施の形態3の変形例に係る画像処理装置が実行する代表値算出処理について説明する。なお、上述した実施の形態と同一の構成には同一の符号を付して説明を省略する。
 図21は、本発明の実施の形態3の変形例に係るノイズ低減部の詳細な構成を示すブロック図である。図21に示すノイズ低減部35b1は、上述した実施の形態3に係るノイズ低減部35bの構成に加えて、信頼度算出部34aをさらに備える。
 信頼度算出部34aは、動き検出部34が算出した被写体の移動量の確からしさを表す信頼度を算出する。
 〔代表値算出処理の概要〕
 図22は、本発明の実施の形態3の変形例に係る代表値算出処理の概要を示すフローチャートである。
 図22において、ステップS401~ステップS410は、上述した図20のステップS311~ステップS320にそれぞれ対応する。また、ステップS411~ステップS418は、上述した図20のステップS312~ステップS319それぞれに対応する。
 ステップS419において、代表値算出部326は、代表値Repiを決定する。具体的には、代表値算出部326は、ステップS410と同様の処理を行うことによって、代表値Repiを決定する。ステップS419の後、画像処理装置30aは、後述するステップS421へ移行する。
 ステップS420において、代表値算出部326は、過去画像の画素値を代表値Reprとして算出する。ステップS420の後、画像処理装置30aは、後述するステップS421へ移行する。
 続いて、代表値算出部326は、現画像の代表値Repiと過去画像の代表値Reprとの類似度を算出する(ステップS421)。具体的には、代表値算出部326は、現画像の代表値Repiと過去画像の代表値Reprの差(|Repi-Repr|)を算出する。
 その後、代表値算出部326は、動き検出部34から入力された動きベクトルの大きさや信頼度算出部34aが算出した信頼度およびステップS421で算出した現画像の代表値Repiと過去画像の代表値Reprの差(|Repi-Repr|)に基づいて類似度を算出し、その類似度に基づいて現画像の代表値Repiおよび過去画像の代表値Reprの各々のブレンド率bを算出する(ステップS422)。具体的には、現画像の代表値Repiと過去画像の代表値Reprの類似度を動きベクトルの大きさ|Vec|で除算した値(|Repi-Repr|/|Vec|)を類似度として算出し、類似度が大きい場合に、ブレンド率bが大きくなるように算出する。
 図23は、代表値算出部326が算出するブレンド率bの算出方法を模式的に示す図である。図23において、横軸が類似度を示し、縦軸がブレンド率bを示す。また、図23において、直線LB1、曲線LB2および折れ線LB3の各々は、ブレンド率の一例を示す。
 図23に示す場合、代表値算出部326は、類似度が大きくなるとき現画像の状態は、同画素で画素値が急激に変化した場合であるため、過去画像の画素値よりも現画像の画素値に重み付けを行う。
 また、代表値算出部326は、現画像の代表値Repiと過去画像の代表値Reprの類似度を、信頼度算出部34aによって算出された動きベクトルの信頼度Relで乗算した値(|Repi-Repr|*Rel)を類似度として算出してもよい。この場合、類似度が大きくなるとき、信頼度算出部34aによって算出された動きベクトルの信頼度が高く、かつ過去画像が正常画像の画素値である可能性が高いため、過去画像の画素値(代表値)に重み付けを行う。なお、過去画像の注目画素が正常画像の画素値であることを前提としてもよい。また、信頼度Relは、画像間に位置合わせを行った画像に基づいて、各動きベクトルの確からしさを表す値である。この信頼度Relは、上述した実施の形態1の差分平均値mSADや公知のSSD、MCCおよびZMCC等によって算出される。
 図22に戻り、ステップS423以降の説明を続ける。
 ステップS423において、代表値算出部326は、代表値Repを決定する。ステップS423の後、画像処理装置30aは、図15のメインルーチンへ戻る。具体的には、代表値決定部323dは、以下の式(11)によって、代表値Repを決定する。
 Rep=Repi×b+Repr×(1-b)   ・・・(11)
なお、代表値算出部326は、の現画像の画素値Valiと過去画像の画素値Valrを用いて、以下の式(12)から代表値Repを決定してもよい。
 Rep=Vali×b+Valr×(1-b)   ・・・(12)
 なお、代表値算出部326は、記憶部36に現画像と、補正前の過去画像とを記憶し、過去画像の代わりに補正前の過去画像を用いてもよい。
 以上説明した本発明の実施の形態3の変形例によれば、上述した実施の形態3と同様の効果を奏する。
(実施の形態4)
 次に、本発明の実施の形態4について説明する。本実施の形態4に係る撮像システムは、上述した実施の形態2に係る画像処理装置30aのノイズ低減部35aの構成が異なるうえ、本実施の形態4に係る画像処理装置が実行する代表値算出処理が異なる。具体的には、本実施の形態4では、現画像と複数の過去画像とを用いて、RTSノイズが発生しない場合の注目画素の画素値に相当する代表値を算出する。以下においては、本実施の形態4に係るノイズ低減部の構成を説明後、本実施の形態4に係る画像処理装置が実行する代表値算出処理について説明する。なお、上述した実施の形態1に係る撮像システム1と同一の構成には同一の符号を付して説明を省略する。
 図24は、本実施の形態4に係るノイズ低減部の詳細な構成を示すブロック図である。図24に示すノイズ低減部35cは、上述した実施の形態1に係る動き検出部34と、RTSノイズ補正部320cと、を備える。
 RTSノイズ補正部320cは、第3外部I/F部31から入力される現画像に対してRTSノイズを補正するRTSノイズ補正処理を行い、この補正を行ったRAW画像を画像処理部38へ出力する。
 ここで、RTSノイズ補正部320cの詳細な構成について説明する。図24は、RTSノイズ補正部320cの詳細な構成を示すブロック図である。
 図25に示すように、RTSノイズ補正部320cは、上述した実施の形態2に係る代表値算出部323に代えて、代表値算出部327を有する。
 代表値算出部327は、動き検出部34が検出する動きベクトル、現画像および記憶部36が記憶する現画像と時間的に連続する複数の過去画像に基づいて、RTSノイズが発生しない場合の注目画素の画素値に相当する代表値を算出する。
 〔代表値算出処理の処理〕
 次に、画像処理装置30aが実行する代表値算出処理について説明する。図26は、本実施の形態4に係る画像処理装置30aが実行する代表値算出処理の概要を示すフローチャートである。
 図26に示すように、代表値算出部327は、現画像における注目画素の画素値と複数の過去画像の各々の注目画素の画素値の平均値Aveを算出する(ステップS500)。なお、代表値算出部327は、平均値以外に、現画像における注目画素の画素値と複数の過去画像の各々の注目画素の画素値の中央値等の統計値を用いて算出してもよい。
 続いて、代表値算出部327は、上述したステップS500で算出した平均値Aveを代表値に決定する(ステップS501)。ステップS501の後、画像処理装置30aは、上述した図15のメインルーチンへ戻る。
 以上説明した本発明の実施の形態4によれば、代表値算出部323が第3外部I/F部31から入力された現画像と記憶部36が記憶する複数の過去画像とを用いて、RTSノイズが発生しない場合の画素値に相当する代表値を算出するので、現画像の注目画素の周辺画素が大小変化する場合する場合であっても、この大小変化を軽減することができるので、画質を劣化させることなく、RTSノイズを低減することができる。
 なお、本発明の実施の形態4では、代表値算出部323が現画像と複数の過去画像とを用いてRTSノイズが発生しない場合の画素値に相当する代表値を算出していたが、例えば現画像の撮影時刻(生成タイミング)と前後する過去画像または未来(新しい)の画像データに対応する画像(以下、「未来画像」)を用いて代表値を算出してもよい。この場合、過去画像および未来画像をm枚用いる場合、このm枚のうちn枚を過去画像とし、m-n枚を未来画像とすればよい。この未来画像を用いる場合、フレーム遅延が発生する可能性があるので、フレーム遅延を許容できる範囲で未来画像を用いることが好ましい。
 また、本発明の実施の形態4では、代表値算出部323が記憶部36に現画像と時間的に連続した補正前の画像を記憶し、過去画像の代わりに補正前の過去画像を用いてもよい。
(実施の形態4の変形例1)
 次に、本発明の実施の形態4の変形例1について説明する。本実施の形態4の変形例1は、ノイズ低減部の構成が異なるうえ、画像処理装置が実行する代表値算出処理が異なる。具体的には、本実施の形態4の変形例1に係る画像処理装置は、現画素と、補正前の過去画像を用いて代表値を算出する。以下においては、本実施の形態4の変形例1に係るノイズ低減部の構成を説明後、本実施の形態4の変形例1に係る画像処理装置が実行する代表値算出処理について説明する。なお、上述した実施の形態と同一の構成には同一の符号を付して説明を省略する。
 図27は、本発明の実施の形態4の変形例1に係るノイズ低減部の詳細な構成を示すブロック図である。図27に示すノイズ低減部35c1は、上述した実施の形態4に係るノイズ低減部35cの構成に加えて、信頼度算出部34aをさらに備える。
 〔代表値算出処理の概要〕
 図28は、本実施の形態4の変形例1に係る画像処理装置30aが実行する代表値算出処理の概要を示すフローチャートである。
 図28に示すように、許容範囲算出部323cは、注目画素のRTSノイズ情報に基づいて、代表値算出処理に使用可能な画素値の範囲である許容範囲Th(有効範囲)を算出する(ステップS511)。具体的には、許容範囲算出部323cは、許容範囲Th(有効範囲)を以下の式(13)によって算出する。
 Th=注目画素のRTS_Value×Co   ・・・(13)
 ここで、Coは、所定の係数であり、例えば1以上の値が好ましい。なお、Coは、予め一定の値を設定してもよいし、ユーザが任意の値を設定するようにしてもよい。
 続いて、代表値算出部327は、カウンタを初期化(i=0)し(ステップS512)、記憶部36から入力された複数の補正前の過去画像におけるi番目画像の注目画素の画素値Val(i)を抽出する(ステップS513)。具体的には、代表値算出部327は、記憶部36から入力されたm枚の補正前の過去画像を時系列順に、0、1、・・・m-1番目(値が大きい方が過去)の補正前の過去画像とする。なお、ステップS513で記憶部36から入力された複数の補正前の過去画像の各々は、注目画素の位置合わせの処理が施されているものとする。また、0番目の過去画像は現画像と同じである。
 その後、代表値算出部327は、現画像の注目画素の画素値Valとi番目の補正前の過去画像の注目画素の画素値Val(i)との差が許容範囲Th以下(|Val(i)-Val|≦Th)であるか否かを判断する(ステップS514)。代表値算出部327が現画像の注目画素の画素値Valとi番目の補正前の過去画像の注目画素の画素値Val(i)との差が許容範囲Th以下であると判断した場合(ステップS514:Yes)、画像処理装置30aは、後述するステップS515へ移行する。これに対して、代表値算出部327が現画像の注目画素の画素値Valとi番目の補正前の過去画像の注目画素の画素値Val(i)との差が許容範囲Th以下でないと判断した場合(ステップS514:No)、画像処理装置30aは、後述するステップS517へ移行する。
 ステップS515において、代表値算出部327は、信頼度算出部34aによって算出されたi番目の補正前の過去画像の動きベクトルの信頼度Rel(i)が所定の閾値Rel_Th以上(Rel(i≧(Rel_Th))であるか否かを判断する。ここで、閾値Rel_Thは、予め一定の値を設定してもよいし、ユーザが任意の値を設定する構成としてもよい。代表値算出部327が信頼度算出部34aによって算出されたi番目の補正前の過去画像の動きベクトルの信頼度Rel(i)が所定の閾値Rel_Th以上であると判断した場合(ステップS515:Yes)、画像処理装置30aは、後述するステップS516へ移行する。これに対して、代表値算出部327が信頼度算出部34aによって算出されたi番目の補正前の過去画像の動きベクトルの信頼度Rel(i)が所定の閾値Rel_Th以上でないと判断した場合(ステップS515:No)、画像処理装置30aは、後述するステップS517へ移行する。
 ステップS516において、代表値算出部327は、記憶部36のmVal[]にi番目の補正前の過去画像の注目画素の画素値Val(i)を記憶する。
 続いて、代表値算出部327は、カウンタiが記憶部36から入力された補正前の過去画像の枚数mより小さい(i<m)か否かを判断する(ステップS517)。代表値算出部327によってカウンタiが記憶部36から入力された補正前の過去画像の枚数mより小さいと判断した場合(ステップS517:Yes)、代表値算出部327は、カウンタをインクリメント(i=i+1)し(ステップS518)、ステップS513へ戻る。これに対して、代表値算出部327によってカウンタiが記憶部36から入力された補正前の過去画像の枚数mより小さくないと判断した場合(ステップS517:No)、画像処理装置30aは、ステップS519へ移行する。
 ステップS519において、代表値算出部327は、記憶部36のmVal[]に現画像の注目画素の画素値Valを記憶する。
 続いて、代表値算出部327は、記憶部36のmVal[]に記憶されている全ての画素値の平均値Aveを算出する(ステップS520)。なお、代表値算出部327は、各画像の注目画素の画素値に対して重み付けを行って平均値を算出してもよい。この重み付け係数の算出方法として、動きベクトルの信頼度が大きい場合に重み付け係数も大きくなるようにしてもよいし、|Val(i)-Val|が小さい場合に重み付け係数を小さくするようにしてもよいし、複数のパラメータ(例えば動きベクトルの信頼度や差等)の組み合わせから重み付けを算出するようにしてもよい。また、代表値算出部327は、平均値や重み付けを行った平均値以外に、中央値や最頻値等の他の統計値を用いてもよい。
 その後、代表値算出部327は、ステップS520で算出した平均値Aveを代表値に決定する(ステップS521)。ステップS521の後、画像処理装置30aは、図115のメインルーチンへ戻る。 
 以上説明した本発明の実施の形態4の変形例1によれば、代表値算出部327が現画像における注目画素の画素値に基づいて、現画像における注目画素の位置に対応する補正前の複数の過去画像の各々の画素の画素値から抽出して算出した値を代表値として決定する。これにより、被写体の輝度値が変化した場合であっても、現画像の画素値に類似する画素値のみからダ代表値を算出するので、代表値の信頼性を向上させることができる。
(実施の形態4の変形例2)
 次に、本発明の実施の形態4の変形例2について説明する。本実施の形態4の変形例2は、上述した実施の形態4に係る画像処理装置30aが実行する代表値算出処理のみが異なる。具体的には、本実施の形態4の変形例2に係る画像処理装置30aは、補正前の過去画像と、過去画像と、を用いて代表値を算出する。以下においては、本実施の形態4の変形例2に係る画像処理装置30aが実行する代表値算出処理について説明する。
 〔代表値算出処理の概要〕
 図29は、本実施の形態4の変形例2に係る画像処理装置30aが実行する代表値算出処理の概要を示すフローチャートである。
 図29において、ステップS601~ステップS607は、上述した図26のステップS511~ステップS518にそれぞれ対応する。なお、ステップS602においては、i=0の必要がなく、i=1とし、補正前の過去画像のみを使って代表値を算出するようにしてもよい。
 ステップS608において、代表値算出部327は、記憶部36のmVal[]に記憶された画素数が閾値ThRefより大きいか(mVal[]>ThRef)否かを判断する。代表値算出部327が記憶部36のmVal[]に記憶された画素数が閾値ThRefより大きいと判断した場合(ステップS608:Yes)、画像処理装置30aは、後述するステップS611へ移行する。これに対して、代表値算出部327が記憶部36にmVal[]に記憶された画素数が閾値ThRefより大きくないと判断した場合(ステップS608:No)、画像処理装置30aは、後述するステップS609へ移行する。
 ステップS609において、代表値算出部327は、現画像の注目画素の画素値Valと過去画像の注目画素の画素値Val(i)との差が許容範囲Th以下(|Val(i)-Val|≦Th)であるか否かを判断する(ステップS609)。代表値算出部327が現画像の注目画素の画素値Valと過去画像の注目画素の画素値Val(i)との差が許容範囲Th以下であると判断した場合(ステップS609:Yes)、画像処理装置30aは、後述するステップS610へ移行する。これに対して、代表値算出部327が現画像の注目画素の画素値Valと過去画像の注目画素の画素値Val(i)との差が許容範囲Th以下でないと判断した場合(ステップS609:No)、画像処理装置30aは、後述するステップS611へ移行する。
 ステップS610において、代表値算出部327は、記憶部36のmVal[]に過去画像の注目画素の画素値Val(i)を記憶する。ステップS610の後、画像処理装置30aは、ステップS611へ移行する。
 ステップS611~ステップS613は、上述した図26のステップS519~ステップS521それぞれに対応する。ステップS613の後、画像処理装置30aは、図15のメインルーチンへ戻る。
 以上説明した本発明の実施の形態4の変形例2によれば、代表値算出部327が補正前の過去画像における注目画素から類似する画素値が少ない場合、過去画素の画素値を参照して代表値を算出するので、代表値の信頼性を向上させることができる。
(実施の形態4の変形例3)
 次に、本発明の実施の形態4の変形例3について説明する。本実施の形態4の変形例3は、上述した実施の形態4に係る画像処理装置30aが実行する代表値算出処理が異なる。具体的には、本実施の形態4の変形例3に係る画像処理装置30aは、過去画像の周辺画素の画素値を用いて代表値を算出する。以下においては、本実施の形態4の変形例3に係る画像処理装置30aが実行する代表値算出処理について説明する。
 〔代表値算出処理の処理〕
 図30は、本実施の形態4の変形例3に係る画像処理装置30aが実行する代表値算出処理の概要を示すフローチャートである。
 図30において、ステップS701~ステップS704は、上述した図16のステップS201~ステップS204にそれぞれ対応する。
 ステップS705において、カウンタを初期化(i=0)し(ステップS705)、許容範囲Th内の周辺画素の画素値を記憶部36のmVal[]に記憶する(ステップS706)。
 続いて、記憶部36に記憶したmVal[]の画素数が予め設定した閾値ThRefより多い場合(ステップS707:Yes)または、カウンタiが記憶部36から入力された過去画像の枚数mより多い場合において、画像処理装置30aは、後述するステップS709へ移行する。
 ステップS707において、記憶部36に記憶したmVal[]の画素数が予め設定した閾値ThRefより多くない場合かつ、カウンタiが記憶部36から入力された過去画像の枚数mより多くない場合(ステップS707:No)、画像処理装置30aは、後述するステップS708へ移行する。
 ステップS708において、代表値算出部327は、カウンタをインクリメント(i=i+1)する。ステップS708の後、画像処理装置30aは、ステップS706へ戻る。
 ステップS709およびステップS710は、上述した図27のステップS612およびステップS613それぞれに対応する。
 また、代表値算出部327が記憶部36に現画像と時間的に連続した補正前の画像を記憶し、過去画像の代わりに補正前の過去画像を用いてもよい。
 以上説明した本発明の実施の形態4の変形例3によれば、代表値算出部327が過去画像の周辺画素の画素値を用いて代表値を算出するので、代表値の信頼性を向上させることができる。
(実施の形態5)
 次に、本発明の実施の形態3について説明する。上述した実施の形態1は、画像処理装置30が個別に設けられていたが、本実施の形態5では、撮像装置本体に画像処理装置を設ける。このため、以下においては、上述した実施の形態1に係る撮像システム1と同様の構成には同一の符号を付して説明を省略する。
 〔撮像システムの構成〕
 図31は、本発明の実施の形態5に係る撮像システム2の構成を模式的に示すブロック図である。図31に示す撮像システム2は、本体部3と、本体部3に着脱自在に接続可能なレンズ部4と、を備える。
 〔本体部の構成〕
 本体部3は、シャッタ103と、撮像素子105と、アナログ処理部106と、A/D変換部107と、操作部108と、メモリI/F部109と、記録媒体110と、揮発メモリ111と、不揮発メモリ112と、バス113と、撮像制御部114と、AE処理部116と、AF処理部117と、外部I/F部118と、表示部119と、ドライバ120と、RTSノイズ補正部320と、を備える。ドライバ120は、撮像制御部114の制御のもと、シャッタ103を駆動する。
 AE処理部116は、バス113を介して揮発メモリ111に記憶された画像データを取得し、この取得した画像データに基づいて、静止画撮影または動画撮影を行う際の露出条件を設定する。具体的には、AE処理部116は、画像データから輝度を算出し、算出した輝度に基づいて、例えば絞り値、露光時間、ISO感度等を決定することで撮像システム2の自動露出(Auto Exposure)を行う。
 AF処理部117は、バス113を介して揮発メモリ111に記憶された画像データを取得し、取得した画像データに基づいて、撮像システム2の自動焦点の調整を行う。例えば、AF処理部117は、画像データから高周波成分の信号を取り出し、高周波成分の信号に対してAF(Auto Focus)演算処理を行うことによって、撮像システム2の合焦評価を決定することで撮像システム2の自動焦点の調整を行う。なお、撮像システム2の自動焦点の調整方法は、撮像素子105で位相差信号を取得するものであってもよい。
 外部I/F部118は、本体部3における各種ブロックにおけるデータの読み書きや、専用のコマンド等による制御などを行える。外部I/F部118は、FPGA、DSPまたはGPU等を搭載した専用の回路やパーソナルコンピュータ(PC)等の外部機器を接続することで、本体部3における各種ブロックを制御が可能なインターフェースである。
 表示部119は、液晶または有機EL(Electro Luminescence)等からなる表示パネルを用いて構成される。表示部119は、撮像素子105が生成した画像データに対応する画像を表示する。
 〔レンズ部の構成〕
 図31に示すように、レンズ部4は、所定の視野領域から集光した被写体像を撮像素子105に結像する。光学系101と、絞り102と、ドライバ104と、を備える。
 〔撮像システムの処理〕
 次に、撮像システム2が実行する処理について説明する。図32は、撮像システム2が実行する処理の概要を示すフローチャートである。
 図32に示すように、まず、ユーザによって操作部108の電源ボタン(図示せず)が操作されて、本体部3の電源がオンになると、撮像制御部114は、撮像システム2の初期化を行う(ステップS801)。具体的には、撮像制御部114は、動画の記録中を示す記録中フラグをオフ状態にする初期化を行う。この記録中フラグは、動画の撮影中にオン状態になり、動画を撮影していないときにオフ状態となるフラグであり、揮発メモリ111に記憶されている。
 続いて、操作部108の動画ボタンが押された場合(ステップS802:Yes)、撮像制御部114は、オン状態で動画の記録中であることを示す記録中フラグを反転し(ステップS803)、撮像制御部114は、撮像システム2が動画記録中であるか否かを判断する(ステップS804)。具体的には、撮像制御部114は、揮発メモリ111に記憶された記録中フラグがオン状態であるか否かを判定する。撮像制御部114によって撮像システム2が動画記録中であると判断された場合(ステップS804:Yes)、撮像システム2は、後述するステップS805へ移行する。これに対して、撮像制御部114によって撮像システム2が動画記録中でないと判断された場合(ステップS804:No)、撮像システム2は、後述するステップS806へ移行する。
 ステップS805において、撮像制御部114は、記録媒体110に画像データを時系列に沿って記録するための動画ファイルを生成する。ステップS805の後、撮像システム2は、後述するステップS806へ移行する。
 ステップS802において、操作部108の動画ボタンが押されていない場合(ステップS802:No)、撮像システム2は、ステップS806へ移行する。
 続いて、撮像制御部114は、撮像システム2が動画の記録中であるか否かを判断する(ステップS806)。撮像制御部114によって撮像システム2が動画の記録中であると判断された場合(ステップS806:Yes)、撮像システム2は、後述するステップS817へ移行する。これに対して、撮像制御部114によって動画の記録中でないと判断された場合(ステップS806:No)、撮像システム2は、後述するステップS807へ移行する。
 ステップS807において、操作部108の再生ボタンが押された場合(ステップS807:Yes)、撮像システム2は、記録媒体110に記録された画像データに対応する画像を表示部119に再生させて表示させる(ステップS808)。ステップS808の後、撮像システム2は、後述するステップS809へ移行する
 ステップS807において、操作部108の再生ボタンが押されていない場合(ステップS807:No)、撮像システム2は、ステップS809へ移行する。
 続いて、操作部108のメニューボタンが押された場合(ステップS809:Yes)、撮像システム2は、各種設定を行う設定処理を実行する(ステップS810)。ステップS810の後、撮像システム2は、後述するステップS811へ移行する。
 ステップS809において、操作部108のメニューボタンが押されていない場合(ステップS809:No)、撮像システム2は、ステップS811へ移行する。
 ステップS811において、操作部108のレリーズボタンがオフ状態から1st状態に遷移した場合(ステップS811:Yes)、撮像制御部114は、AE処理部116に露出を調整するAE処理およびAF処理部117にピントを調整するAF処理のそれぞれを実行させる(ステップS812)。その後、撮像システム2は、後述するステップS824へ移行する。
 ステップS811において、操作部108のレリーズボタンがオフ状態から1st状態に遷移していない場合(ステップS811:No)、撮像システム2は、ステップS813へ移行する。
 続いて、操作部108のレリーズボタンが2nd状態に遷移した場合(ステップS813:Yes)、撮像制御部114は、メカシャッタによる撮影を実行する(ステップS814)。具体的には、撮像制御部114は、シャッタ103を制御することによって、撮像素子105に撮影を実行させる。
 続いて、撮像システム2は、撮像素子105が生成した画像データに対して、RTSノイズの補正を行った後に、所定の処理を行う画像処理を実行する(ステップS815)。なお、画像処理の詳細は後述する。
 その後、撮像制御部114は、画像処理部303が画像処理を施した画像データを記録媒体110に記録する(ステップS816)。ステップS816の後、撮像システム2は、後述するステップS824へ移行する。
 ステップS813において、操作部108のレリーズボタンが2nd状態に遷移していない場合(ステップS813:No)、撮像システム2は、ステップS817へ移行する。
 続いて、撮像制御部114は、AE処理部116に露出を調整するAE処理を実行させ(ステップS817)、AF処理部117にピントを調整するAF処理を実行させる(ステップS818)。
 その後、撮像制御部114は、ドライバ120を介して撮像素子105に露光時間を電子的に制御する、所謂、電子シャッタによる撮影を実行させる(ステップS819)。電子シャッタによる撮影によって撮像素子105が生成した画像データは、アナログ処理部106およびA/D変換部107およびバス113を介して揮発メモリ111に出力される。
 続いて、撮像システム2は、ステップS815と同様の画像処理を実行する(ステップS820)。なお、画像処理の詳細は後述する。
 その後、撮像システム2は、電子シャッタによる撮影によって撮像素子105が生成した画像データに対応するライブビュー画像を表示部119に表示させる(ステップS821)。
 続いて、撮像システム2が動画記録中である場合(ステップS822:Yes)、撮像制御部114は、画像データをステップS810の設定処理により設定した記録形式で図示しない画像圧縮展開部に圧縮させ、この圧縮した画像データを記録媒体110に作成された動画ファイルに動画として記録させる(ステップS823)。ステップS823の後、撮像システム2は、ステップS824へ移行する。
 ステップS822において、撮像システム2が動画記録中でない場合(ステップS822:No)、撮像システム2は、ステップS824へ移行する。
 続いて、操作部108の電源ボタンが押されて撮像システム2の電源がオフ状態になった場合(ステップS824:Yes)、撮像システム2は、本処理を終了する。これに対して、撮像システム2の電源がオフ状態になっていない場合(ステップS824:No)、撮像システム2は、ステップS802へ戻る。
 〔画像処理の概要〕
 次に、図32のステップS815およびステップS820で説明した画像処理について説明する。図33は、画像処理の概要を示すフローチャートである。
 図33に示すように、RTSノイズ補正部320は、撮像素子105が生成した画像データに対してRTSノイズを補正するRTSノイズ補正処理を実行する(ステップS901)。ここで、RTSノイズ補正処理は、上述した実施の形態1に係る画像処理装置30が実行する処理または上述した実施の形態2~4に係る画像処理装置30aが実行する処理に対応するため、説明を省略する。
 続いて、画像処理部303は、RTSノイズ補正部320がRTSノイズを補正した画像データに対して基本画像処理を実行する(ステップS902)。ステップS902の後、撮像システム2は、図32のメインルーチンへ戻る。
 以上説明した本発明の実施の形態5によれば、上述した実施の形態1と同様の効果を有する。
(その他の実施の形態)
 本発明は上述した実施の形態に限定されるものではなく、本発明の要旨の範囲内で種々の変形や応用が可能なことは勿論である。例えば、本発明の説明に用いた撮像装置以外にも、携帯電話やスマートフォンにおける撮像素子を備えた携帯機器、ビデオカメラ、内視鏡、監視カメラ、顕微鏡のような光学機器を通して被写体を撮影する撮像装置等、被写体を撮像可能ないずれの機器にも適用できる。
 本発明では、ランダムノイズモデルの参照値、もしくは参照するノイズモデルを、RTSノイズの特徴量であるRTS_Valueや、RTS_Valueに基づく値である候補値の最大値に応じて変更していたが、これらを組み合わせ、RTSノイズの特徴量に応じてランダムノイズモデルを選択し、さらに参照値を変更してランダムノイズ量を推定するようにしてもよい。また、温度に応じてランダムノイズ量が異なるため、画像データを撮影した際の撮像素子の温度に応じてランダムノイズモデルを選択し、さらに上記の方法を適用するようにしてもよい。
 また、本発明では、RTSノイズ情報記録部が撮像装置内に設けられていたが、上述した画像処理装置内に設けられていてもよいし、ネットワークを介して双方向に通信可能なサーバ内にRTSノイズ情報記録部を設け、ネットワークを介してRTSノイズ情報を取得するようにしてもよい。
 また、本発明は、表示または記録に用いる画像データ以外の画像データ、例えばOB領域の画像データまたは光学的に設計保証されていないイメージサークル外の領域の画像データ等の画像データであっても適用可能である。
 また、本明細書において、前述の各動作フローチャートの説明において、便宜上「まず」、「次に」、「続いて」、「その後」等を用いて動作を説明しているが、この順で動作を実施することが必須であることを意味するものではない。
 また、上述した実施の形態における画像処理装置による各処理の手法、即ち、各フローチャートに示す処理は、いずれもCPU等の制御部に実行させることができるプログラムとして記憶させておくこともできる。この他、メモリカード(ROMカード、RAMカード等)、磁気ディスク(フロッピーディスク(登録商標)、ハードディスク等)、光ディスク(CD-ROM、DVD等)、半導体メモリ等の外部記憶装置の記憶媒体に格納して配布することができる。そして、CPU等の制御部は、この外部記憶装置の記憶媒体に記憶されたプログラムを読み込み、この読み込んだプログラムによって動作が制御されることにより、上述した処理を実行することができる。
 また、本発明は、上述した実施の形態および変形例そのままに限定されるものではなく、実施段階では、発明の要旨を逸脱しない範囲内で構成要素を変形して具体化することができる。また、上述した実施の形態に開示されている複数の構成要素を適宜組み合わせることによって、種々の発明を形成することができる。例えば、上述した実施の形態および変形例に記載した全構成要素からいくつかの構成要素を削除してもよい。さらに、各実施の形態および変形例で説明した構成要素を適宜組み合わせてもよい。
 また、明細書または図面において、少なくとも一度、より広義または同義な異なる用語と共に記載された用語は、明細書または図面のいかなる箇所においても、その異なる用語に置き換えることができる。このように、発明の主旨を逸脱しない範囲内において種々の変形や応用が可能である。
 1,1a,2 撮像システム
 3 本体部
 4 レンズ部
 10 撮像装置
 30,30a 画像処理装置
 31 第3外部I/F部
 32 補間処理部
 33 フレームメモリ
 34 動き検出部
 34a 信頼度算出部
 35,35a,35b,35c,35b1,35c1 ノイズ低減部
 36 記憶部
 37 操作部
 38 画像処理部
 39 画像処理制御部
 40 表示装置
 101 光学系
 102 絞り
 103 シャッタ
 104 ドライバ
 105 撮像素子
 105a 画素
 105b 第1スイッチ
 105c 垂直転送線
 105d FD部
 105e アンプ部
 105f 第2スイッチ
 105h 転送線
 106 アナログ処理部
 107 A/D変換部
 108 操作部
 109 メモリI/F部
 110 記録媒体
 111 揮発メモリ
 112 不揮発メモリ
 112a プログラム記録部
 112b RTSノイズ情報記録部
 112c ランダムノイズモデル情報記録部
 113 バス
 114 撮像制御部
 115 第1外部I/F部
 116 AE処理部
 117 AF処理部
 118 外部I/F部
 119 表示部
 120 ドライバ
 303 画像処理部
 310 補間処理部
 320,320b,320c RTSノイズ補正部
 321 RTSノイズ画素判定部
 322 候補値算出部
 323,326,327 代表値算出部
 323a 第1参照値算出部
 323b 第1ランダムノイズ量推定部
 323c 許容範囲算出部
 324 ランダムノイズ量推定部
 325 補正値算出部
 325a 第2参照値算出部
 325b 第2ランダムノイズ量推定部
 325c 補正量決定部
 325d 画素値補正部
 351 評価値算出部
 352 推定ノイズ量取得部
 353 ルックアップテーブル
 354 判定部
 355 ノイズ低減処理部

Claims (20)

  1.  2次元状に配置され、外部から光を受光し、受光量に応じた信号を生成する複数の画素と、前記信号を画素値として読み出す複数の読み出し回路と、を有する撮像素子が生成した画像データに含まれる点滅欠陥ノイズを補正する画像処理装置であって、
     前記読み出し回路に起因する点滅欠陥ノイズが発生する前記読み出し回路の位置情報または前記複数の画素の各々の位置情報を含むノイズ情報と、前記画像データと、を取得する取得部と、
     前記取得部が取得した前記画像データを補正対象画像データとし、該補正対象画像データと異なる時間に取得された前記画像データに基づく参照画像データと、前記ノイズ情報とに基づいて、前記補正対象画像データにおける注目画素の画素値を補正する補正部と、
     を備えたことを特徴とする画像処理装置。
  2.  前記補正対象画像データと前記参照画像データとに基づいて、被写体の移動量を算出する移動量算出部と、
     を備え、
     前記補正部は、前記移動量算出部が算出した前記移動量に基づいて、前記注目画素または前記注目画素近傍の画素と対応する前記参照画像データの参照画素を取得し、該参照画素の画素値に基づいて、前記注目画素の画素値を補正することを特徴とする請求項1に記載の画像処理装置。
  3.  前記補正部は、前記注目画素において前記点滅欠陥ノイズが発生する場合、前記点滅欠陥ノイズが発生しない場合に比べ、前記注目画素におけるノイズ低減処理の強度を大きくすることを特徴とする請求項1または2に記載の画像処理装置。
  4.  前記補正部は、前記注目画素において前記点滅欠陥ノイズが発生する場合、前記参照画像データを用いた時間方向のノイズ低減処理の強度を大きくすることを特徴とする請求項3に記載の画像処理装置。
  5.  前記補正部は、前記注目画素において前記点滅欠陥ノイズが発生する場合、前記注目画素周辺の画素を用いた空間方向のノイズ低減処理の強度を大きくすることを特徴とする請求項3に記載の画像処理装置。
  6.  前記補正部は、前記参照画素の画素値に基づいて、前記点滅欠陥ノイズが発生しない場合の画素値に相当する代表値を算出し、該代表値に基づいて、前記注目画素の画素値を補正することを特徴とする請求項2に記載の画像処理装置。
  7.  前記補正部は、前記注目画素の画素値に基づいて、前記代表値を算出することを特徴とする請求項6に記載の画像処理装置。
  8.  前記補正部は、前記点滅欠陥ノイズが発生しない前記参照画素の画素値もしくは前記点滅欠陥ノイズを補正した前記参照画素の画素値に基づいて、前記代表値を算出することを特徴とする請求項6に記載の画像処理装置。
  9.  前記補正部は、前記点滅欠陥ノイズが発生しない前記注目画素における周辺の画素の画素値に基づいて、前記代表値を算出することを特徴とする請求項6~8のいずれか一つに記載の画像処理装置。
  10.  前記注目画素周辺のランダムノイズ量を推定するノイズ量推定部と、
     をさらに備え、
     前記補正部は、前記ランダムノイズ量に基づいて、前記代表値を算出することを特徴とする請求項6~9のいずれか一つに記載の画像処理装置。
  11.  前記ノイズ情報は、前記読み出し回路に起因する点滅欠陥ノイズが発生する前記読み出し回路の位置情報または前記複数の画素の各々の位置情報と対応づけられた、前記点滅欠陥ノイズのノイズレベルである点滅欠陥ノイズレベルをさらに含み、
     前記補正部は、前記点滅欠陥ノイズレベルを用いて前記注目画素の画素値を補正することを特徴とする請求項1~10のいずれか一つに記載の画像処理装置。
  12.  前記ノイズ情報は、前記読み出し回路に起因する点滅欠陥ノイズが発生する前記読み出し回路の位置情報または前記複数の画素の各々の位置情報と対応づけられた、前記点滅欠陥ノイズのノイズレベルである点滅欠陥ノイズレベルをさらに含み、
     前記補正部は、前記点滅欠陥ノイズレベルが大きいほど、前記注目画素におけるノイズ低減処理の強度を大きくすることを特徴とする請求項3~5のいずれか一つに記載の画像処理装置。
  13.  前記ノイズ情報は、前記読み出し回路に起因する点滅欠陥ノイズが発生する前記読み出し回路の位置情報または前記複数の画素の各々の位置情報と対応づけられた、前記点滅欠陥ノイズのノイズレベルである点滅欠陥ノイズレベルをさらに含み、
     前記補正部は、前記点滅欠陥ノイズレベルに基づいて、前記代表値を算出することを特徴とする請求項6~10のいずれか一つに記載の画像処理装置。
  14.  前記移動量算出部が算出した前記移動量の確からしさを表す信頼度を算出する信頼度算出部を備え、
     前記補正部は、前記信頼度に基づいて、前記代表値を算出することを特徴とする、請求項6~10のいずれか一つに記載の画像処理装置。
  15.  前記参照画像データは、前記補正対象画像データに対応する前記画像データの直前に前記撮像素子が生成した前記画像データであることを特徴とする請求項1~14のいずれか一つに記載の画像処理装置。
  16.  前記参照画像データは、前記補正部が前記補正対象画像データに対する補正を行う直前に、前記補正部が補正を行った画像データであることを特徴とする請求項1~14のいずれか一つに記載の画像処理装置。
  17.  前記参照画像データは、前記補正対象画像データに対応する前記画像データの前に、前記撮像素子が生成した複数の前記画像データであることを特徴とする請求項1~15のいずれか一つに記載の画像処理装置。
  18.  前記点滅欠陥ノイズは、ランダムテレグラフシグナルノイズであることを特徴とする請求項1~17のいずれか一つに記載の画像処理装置。
  19.  2次元状に配置され、外部から光を受光し、受光量に応じた信号を生成する複数の画素と、前記信号を画素値として読み出す複数の読み出し回路と、を有する撮像素子が生成した画像データに含まれる点滅欠陥ノイズを補正する画像処理装置が実行する画像処理方法であって、
     前記読み出し回路に起因する点滅欠陥ノイズが発生する前記読み出し回路の位置情報または前記複数の画素の各々の位置情報を含むノイズ情報と、前記画像データと、を取得する取得ステップと、
     前記取得ステップにおいて取得した前記画像データを補正対象画像データとし、該補正対象画像データと異なる時間に取得された前記画像データに基づく参照画像データと、前記ノイズ情報とに基づいて、前記補正対象画像データにおける注目画素の画素値を補正する補正ステップと、
     を含むことを特徴とする画像処理方法。
  20.  2次元状に配置され、外部から光を受光し、受光量に応じた信号を生成する複数の画素と、前記信号を画素値として読み出す複数の読み出し回路と、を有する撮像素子が生成した画像データに含まれる点滅欠陥ノイズを補正する画像処理装置に、
     前記読み出し回路に起因する点滅欠陥ノイズが発生する前記読み出し回路の位置情報または前記複数の画素の各々の位置情報を含むノイズ情報と、前記画像データと、を取得する取得ステップと、
     前記取得ステップにおいて取得した前記画像データを補正対象画像データとし、該補正対象画像データと異なる時間に取得された前記画像データに基づく参照画像データと、前記ノイズ情報とに基づいて、前記補正対象画像データにおける注目画素の画素値を補正する補正ステップと、
     を実行させることを特徴とするプログラム。
PCT/JP2015/071175 2015-07-24 2015-07-24 画像処理装置、画像処理方法およびプログラム Ceased WO2017017742A1 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP2017530484A JP6727205B2 (ja) 2015-07-24 2015-07-24 画像処理装置、画像処理方法およびプログラム
PCT/JP2015/071175 WO2017017742A1 (ja) 2015-07-24 2015-07-24 画像処理装置、画像処理方法およびプログラム
DE112015006619.5T DE112015006619T5 (de) 2015-07-24 2015-07-24 Bildverarbeitungsgerät, Bildverarbeitungsverfahren und Programm
CN201580081299.0A CN107710734B (zh) 2015-07-24 2015-07-24 图像处理装置、图像处理方法和记录介质
US15/845,166 US10368017B2 (en) 2015-07-24 2017-12-18 Image processing apparatus, image processing method, and computer-readable recording medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2015/071175 WO2017017742A1 (ja) 2015-07-24 2015-07-24 画像処理装置、画像処理方法およびプログラム

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US15/845,166 Continuation US10368017B2 (en) 2015-07-24 2017-12-18 Image processing apparatus, image processing method, and computer-readable recording medium

Publications (1)

Publication Number Publication Date
WO2017017742A1 true WO2017017742A1 (ja) 2017-02-02

Family

ID=57884484

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2015/071175 Ceased WO2017017742A1 (ja) 2015-07-24 2015-07-24 画像処理装置、画像処理方法およびプログラム

Country Status (5)

Country Link
US (1) US10368017B2 (ja)
JP (1) JP6727205B2 (ja)
CN (1) CN107710734B (ja)
DE (1) DE112015006619T5 (ja)
WO (1) WO2017017742A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10977803B2 (en) * 2017-04-10 2021-04-13 Olympus Corporation Correlation value calculation device

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019194106A1 (ja) * 2018-04-06 2019-10-10 オリンパス株式会社 画像処理装置、画像処理方法および画像処理プログラム
EP3629240B1 (en) * 2018-09-07 2023-08-23 Panasonic Intellectual Property Corporation of America Generative adversarial networks for local noise removal from an image
WO2020240853A1 (ja) * 2019-05-31 2020-12-03 日本電信電話株式会社 画像処理装置、画像処理方法及びプログラム
SG11202111977XA (en) * 2019-06-03 2021-11-29 Hamamatsu Photonics Kk Method for inspecting semiconductor and semiconductor inspecting device
WO2021260871A1 (ja) * 2020-06-25 2021-12-30 三菱電機株式会社 赤外線画像処理装置および赤外線画像処理方法
US11750944B2 (en) * 2021-05-28 2023-09-05 Varex Imaging Corporation Pixel noise cancellation system
CN114862804A (zh) * 2022-05-16 2022-08-05 Oppo广东移动通信有限公司 检测方法及装置、电子设备及存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010124321A (ja) * 2008-11-20 2010-06-03 Olympus Corp 動画像ノイズ低減処理装置、動画像ノイズ低減処理プログラム及び動画像ノイズ低減処理方法
JP2014107611A (ja) * 2012-11-26 2014-06-09 Toshiba Digital Media Engineering Corp 画像処理装置及び画像処理方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4586052B2 (ja) * 2007-08-08 2010-11-24 キヤノン株式会社 画像処理装置及びその制御方法
JP5315158B2 (ja) * 2008-09-12 2013-10-16 キヤノン株式会社 画像処理装置及び画像処理方法
JP5561112B2 (ja) 2010-11-10 2014-07-30 株式会社ニコン 画像処理装置、撮像装置およびプログラム
JP2013121099A (ja) * 2011-12-08 2013-06-17 Sony Corp 画像処理装置、および画像処理方法、並びにプログラム
JP6160292B2 (ja) * 2013-06-24 2017-07-12 富士通株式会社 画像補正装置、撮像装置及び画像補正用コンピュータプログラム
JP2015115922A (ja) * 2013-12-16 2015-06-22 オリンパス株式会社 撮像装置および撮像方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010124321A (ja) * 2008-11-20 2010-06-03 Olympus Corp 動画像ノイズ低減処理装置、動画像ノイズ低減処理プログラム及び動画像ノイズ低減処理方法
JP2014107611A (ja) * 2012-11-26 2014-06-09 Toshiba Digital Media Engineering Corp 画像処理装置及び画像処理方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10977803B2 (en) * 2017-04-10 2021-04-13 Olympus Corporation Correlation value calculation device

Also Published As

Publication number Publication date
US20180109745A1 (en) 2018-04-19
DE112015006619T5 (de) 2018-05-09
JP6727205B2 (ja) 2020-07-22
CN107710734B (zh) 2020-05-12
US10368017B2 (en) 2019-07-30
JPWO2017017742A1 (ja) 2018-05-17
CN107710734A (zh) 2018-02-16

Similar Documents

Publication Publication Date Title
JP6727205B2 (ja) 画像処理装置、画像処理方法およびプログラム
JP6320571B2 (ja) 画像処理装置、画像処理方法およびプログラム
JP5643563B2 (ja) 画像処理装置及びその制御方法
JP6223059B2 (ja) 撮像装置、その制御方法およびプログラム
JP6366739B2 (ja) 画像処理装置、画像処理方法およびプログラム
KR20130061083A (ko) 촬상장치, 촬상장치의 제어 방법, 및 기억매체
JP6471222B2 (ja) 画像処理装置、画像処理方法およびプログラム
US20140176754A1 (en) Image capturing apparatus and control method thereof
JP2016142999A (ja) 撮像装置及びその制御方法
JP6305290B2 (ja) 画像処理装置、撮像装置、画像処理方法
JP6516830B2 (ja) 画像処理装置、画像処理方法およびプログラム
US10425602B2 (en) Image processing apparatus, image processing method, and computer-readable recording medium
JP6080579B2 (ja) 撮像装置およびその制御方法
JP2020046489A (ja) 撮像制御装置、撮像装置および撮像制御プログラム
JP2018054815A (ja) 領域判別装置、その制御方法、および制御プログラム、並びに撮像装置
JP3206826U (ja) 画像処理装置
JP5610967B2 (ja) 画像処理装置、撮像装置、画像処理方法および画像処理プログラム
JP2013115758A (ja) 撮像装置、その制御方法及びプログラム

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: 15899576

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 112015006619

Country of ref document: DE

ENP Entry into the national phase

Ref document number: 2017530484

Country of ref document: JP

Kind code of ref document: A

122 Ep: pct application non-entry in european phase

Ref document number: 15899576

Country of ref document: EP

Kind code of ref document: A1