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

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

Info

Publication number
WO2009001956A1
WO2009001956A1 PCT/JP2008/061806 JP2008061806W WO2009001956A1 WO 2009001956 A1 WO2009001956 A1 WO 2009001956A1 JP 2008061806 W JP2008061806 W JP 2008061806W WO 2009001956 A1 WO2009001956 A1 WO 2009001956A1
Authority
WO
WIPO (PCT)
Prior art keywords
color
image
pixel
sampling
unit
Prior art date
Application number
PCT/JP2008/061806
Other languages
English (en)
French (fr)
Inventor
Yasuhiro Sawada
Xianji Zhang
Takashi Masuda
Original Assignee
Acutelogic Corporation
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from JP2007166106A external-priority patent/JP5036421B2/ja
Priority claimed from JP2007335065A external-priority patent/JP5036524B2/ja
Application filed by Acutelogic Corporation filed Critical Acutelogic Corporation
Priority to EP08765855A priority Critical patent/EP2175656A1/en
Priority to KR1020107001700A priority patent/KR101356286B1/ko
Priority to US12/666,511 priority patent/US8494260B2/en
Publication of WO2009001956A1 publication Critical patent/WO2009001956A1/ja
Priority to US13/929,423 priority patent/US8934713B2/en

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4015Image demosaicing, e.g. colour filter arrays [CFA] or Bayer patterns
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/10Cameras or camera modules comprising electronic image sensors; Control thereof for generating image signals from different wavelengths
    • H04N23/12Cameras or camera modules comprising electronic image sensors; Control thereof for generating image signals from different wavelengths with one sensor only
    • 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/10Circuitry of solid-state image sensors [SSIS]; Control thereof for transforming different wavelengths into image signals
    • H04N25/11Arrangement of colour filter arrays [CFA]; Filter mosaics
    • H04N25/13Arrangement of colour filter arrays [CFA]; Filter mosaics characterised by the spectral characteristics of the filter elements
    • H04N25/134Arrangement of colour filter arrays [CFA]; Filter mosaics characterised by the spectral characteristics of the filter elements based on three different wavelength filter elements
    • 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/61Noise processing, e.g. detecting, correcting, reducing or removing noise the noise originating only from the lens unit, e.g. flare, shading, vignetting or "cos4"
    • H04N25/611Correction of chromatic aberration
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N2209/00Details of colour television systems
    • H04N2209/04Picture signal generators
    • H04N2209/041Picture signal generators using solid-state devices
    • H04N2209/042Picture signal generators using solid-state devices having a single pick-up sensor
    • H04N2209/045Picture signal generators using solid-state devices having a single pick-up sensor using mosaic colour filter
    • H04N2209/046Colour interpolation to calculate the missing colour values

Definitions

  • the present invention relates to an image processing device, an image processing method, a program, and an imaging device, and in particular, from a color mosaic image obtained by using a single-plate type power image pickup device, a plurality of colors for all pixels.
  • the present invention relates to an image processing apparatus, an image processing method, a program, and an imaging apparatus that execute a demosaicing process that interpolates luminance information to generate a color image. Background art
  • each pixel has only single color information.
  • the power image is represented by combining three separate monochromatic images. That is, to display a force image, red (R), green
  • demosic processing color interpolation processing
  • the force image generated from the color mosaic image as described above may be subjected to processing such as deformation processing and correction processing.
  • processing such as deformation processing and correction processing.
  • a power zoom image may be enlarged or reduced by digital zoom processing, or a power zoom image may be rotated by camera shake correction processing.
  • zoom lenses for digital cameras often have barrel distortion at the wide-angle end due to their configuration.
  • the force Lull surface image may be deformed.
  • chromatic aberration may be corrected.
  • the lens of the imaging optical system used for the digital camera has a different refractive index depending on the wavelength of the imaging light. For this reason, the magnification of the image differs depending on each color of R GB.
  • the size of the image formed on the image sensor varies from color component to color component, and as shown in FIG. 20, it is known that an image shift on the image sensor occurs for each color component. It has been. This is called lateral chromatic aberration (lateral chromatic aberration) of the lens. If there is double chromatic aberration, when the white point light source is photographed, especially the periphery of the screen appears to be rainbow colored and extend in the radial direction. In addition, with the image shift, the edge portion of the image appears as a color shift, causing a problem that the quality of the image quality is impaired.
  • the color corresponding to the distance from the reference position in the color image is based on the color image signal obtained by imaging.
  • an imaging apparatus that detects an amount of aberration and corrects the color difference of the color image signal based on the detected amount of chromatic aberration.
  • an imaging device that detects effective edges in a color image and detects a chromatic aberration amount according to a distance from a reference position based on the edge position.
  • an incident image from the optical lens system is picked up and R, G, B signals are output, and the R, G, B signals are subjected to chromatic aberration inherent in the optical lens system.
  • R, G, and B After correcting the chromatic aberration using the focal length correction value (correction to enlarge or reduce the image with respect to R and B with G as a reference), R, G, and B outputs corrected for chromatic aberration
  • a digital camera that synthesizes and outputs an image signal corrected for chromatic aberration (see, for example, Patent Document 3).
  • Patent Document 1 Special Table 2 0 0 4— 5 3 4 4 2 9
  • Patent Document 2 Japanese Patent Laid-Open No. 2 00 0-2 7 0 2 94
  • Patent Document 3 Japanese Patent Application Laid-Open No. 6-1 1 3 3 09 Disclosure of Invention
  • the digital camera described in Patent Document 3 is not a single plate type, but a so-called three-plate type digital camera.
  • the three-plate digital camera is R
  • G B are equipped with three image sensors, and R, G, B signals output from each image sensor are combined to obtain a color image.
  • the number of R, G, and B pixels is the same as the number of output images, so a color image can be obtained by relatively simple image composition.
  • demosic processing color interpolation processing
  • demosic processing is not required as is done with a single-plate digital camera that has fewer R, G, and B pixels than the number of pixels in the output image.
  • interpolation processing is performed when generating a color image from a color mosaic image, and another interpolation is performed when correcting lateral chromatic aberration. It is necessary to perform processing. In other words, the pixel interval must be doubled. For this reason, there was a problem that the processing load increases. In addition, since another interpolation is performed on the image generated by the interpolation, there is a problem that the image quality of the generated image is greatly deteriorated.
  • the present invention has been made to solve such a problem, and when generating a force-ratio image that has been subjected to deformation processing from a color mosaic image, with a smaller amount of processing load, The purpose is to generate a color image with little image quality degradation.
  • a color mosaic separation unit that separates a color mosaic image into a plurality of color planes that include only pixel values of the same color light, and a color mosaic from the pixel position of the power error image. Corresponding sample on the image
  • the coordinate conversion unit that calculates sampling coordinates
  • the sampling unit that interpolates and generates pixel values at the sampling coordinates for each color plane, and the interpolated values of each color plane.
  • a color generation unit that generates a color image, and in the coordinate conversion unit, on the color mosaic image corresponding to the pixel position of the color image when the image is deformed from the pixel position of the color image. Sampling coordinates are calculated.
  • the color mosaic corresponding to the pixel position of the color image when the image is deformed as a step of generating a color image from the color mosaic image.
  • the sampling coordinates on the image are calculated.
  • the color interpolation operation for obtaining the interpolation value of the sampling coordinate is performed by the sampling unit, so that each pixel value of the force-transformed image is sampled.
  • the pixel value of the ring coordinate is obtained from the color mosaic image.
  • color interpolation processing for generating a color image from a color mosaic image and deformation processing of the color image can be realized by a single interpolation operation.
  • a color plane separation unit that separates a color mosaic image into a plurality of color planes that include only pixel values of the same color light, and a color mosaic from the pixel position of the color image.
  • a coordinate conversion unit that calculates the corresponding sampling coordinates on the sample image, a sampling unit that interpolates and generates pixel values at the sampling coordinates for each of the color planes, and an interpolation value for each color plane is synthesized.
  • a color generation unit that generates a color image.
  • a coefficient having a different value for each color plane is used for each color plane.
  • the sampling coordinates that differ for each color plane are calculated from the pixel positions of the force raster image.
  • a power mosaic image is generated from a color mosaic image.
  • the image position shift of each color component caused by the chromatic aberration of magnification [the pixel position for each color plane corresponds to the sampling coordinates on the color mosaic image corresponding to the pixel position of the color image. And then calculate
  • the color interpolation calculation to obtain the interpolation value of the sampling coordinate is performed by the sampling V unit, so that each pixel value of the color image corrected for the chromatic aberration of magnification becomes the pixel value of the sampling coordinate. Therefore, it is obtained from the color mosaic image.
  • a color interpolation process (demosaic process) that generates a force image from a color mosaic image and a process for correcting the chromatic aberration of magnification of the imaging optical system can be realized by a single interpolation operation.
  • the processing load when generating a color image with corrected chromatic aberration of magnification from a color mosaic image can be reduced, and the interpolation processing can be performed twice as compared with the conventional method. Degradation of image quality can also be suppressed.
  • the coordinate conversion unit corrects the image deformation coefficient representing the image deformation with respect to the color image in addition to the chromatic aberration coefficient having different values depending on the color plane, and the blur of the imaging apparatus. Using at least one of the distortion aberration coefficient determined according to the focal length and the subject distance set in the imaging optical system for guiding the imaging light of the subject to the image sensor. The sampling coordinates that differ for each color plane are calculated from the pixel position of the image.
  • Each pixel value is obtained from the color mosaic image as the pixel value of the sampling coordinates.
  • the process of correcting the chromatic aberration of magnification of the imaging optical system and the process such as image deformation, blur correction, and distortion correction of the power image can be realized by a single interpolation calculation. This reduces the processing load when generating a color image that has been corrected for chromatic aberration of magnification and has undergone image deformation, etc., from a color mosaic image, and has been subjected to interpolation processing many times. Deterioration of image quality due to this can also be suppressed.
  • an aberration coefficient calculation unit that calculates a BIJ chromatic aberration coefficient by taking a predetermined image is provided, and a color plane is calculated from the pixel position of the color image using the color yield coefficient. Try to calculate different sampling coordinates for each.
  • an aberration coefficient calculation unit that calculates a chromatic aberration coefficient by capturing a predetermined image is provided. Therefore, a color mosaic is used by using the calculated chromatic aberration coefficient. Power capture from images-Color capture processing (demosaic processing) that generates images and processing for correcting magnification chromatic aberration of the imaging optical system can be realized by a single interpolation operation.
  • the image processing apparatus includes an image output unit that outputs a captured image to an external device connected to the outside, and a coefficient input unit that inputs the chromatic aberration coefficient from the external device. Using the coefficients, different sampling coordinates are calculated for each color plane from the pixel position of the color image.
  • an image can be taken on an external device connected to the outside.
  • An image output unit that outputs a shadow image and a coefficient input unit that inputs the chromatic aberration coefficient from the external device are provided, so that a color image can be obtained from the color mosaic image using the input chromatic aberration coefficient.
  • the generated color interpolation process demosaic process
  • the process of correcting the magnification chromatic aberration of the imaging optical system can be realized by a single interpolation operation.
  • FIG. 1 is a diagram illustrating a configuration example of a color imaging apparatus according to a first embodiment in which an image processing apparatus according to the present invention is implemented.
  • FIG. 2 is a diagram showing the primary color layer array of the color filter array according to the first to fifth embodiments.
  • FIG. 3 is a block diagram showing a functional configuration example of the demosaic unit according to the first embodiment.
  • FIG. 4 is a flowchart showing an example of the image processing operation executed by the demosaic unit of the first embodiment.
  • FIG. 5 is an image diagram for specifically explaining an example of image processing executed by the demosaic unit of the first embodiment.
  • FIG. 6 is an image diagram for explaining the bilinear interpolation of the first embodiment.
  • FIG. 7 is an image diagram for specifically explaining another example of the image processing executed by the demosaic unit of the first embodiment.
  • FIG. 8 is a diagram illustrating a configuration example of a color imaging apparatus according to the second embodiment and the third embodiment in which the image processing apparatus according to the present invention is implemented.
  • FIG. 9 is a block diagram illustrating a functional configuration example of the demosaicing unit according to the second embodiment and the third embodiment.
  • Figure 10 shows the images executed by the demosaicing unit of the second to fifth embodiments.
  • 3 is a flowchart showing an example of image processing operation.
  • FIG. 11 is an image diagram for specifically explaining an example of image processing executed by the demosaicing unit of the second embodiment and the third embodiment.
  • FIG. 12 is an image diagram for specifically explaining another example of the image processing executed by the demosaicing unit according to the second to fifth embodiments.
  • FIG. 13 is a diagram showing a configuration example of a force label imaging device according to the fourth embodiment in which the image processing device according to the present invention is implemented.
  • FIG. 14 is a diagram showing an example of an aberration coefficient table according to the fourth embodiment.
  • FIG. 15 is a block diagram showing a functional configuration example of the demosaicing unit according to the fourth embodiment.
  • FIG. 16 is a diagram showing a configuration example of the color imaging apparatus according to the fifth embodiment.
  • FIG. 17 is an explanatory diagram of chromatic aberration amount detection in the fifth embodiment, and FIG. FIG. 17 is a diagram showing the form of the chart, and FIG. 17B is a diagram showing the arrangement of the chart with respect to the image sensor.
  • FIG. 18 is an explanatory diagram of chromatic aberration amount detection in the fifth embodiment, in which FIGS. 18 (a) and (b) are explanatory diagrams when detecting intersection points, and FIG. 18 (c) is an edge for each intersection point.
  • FIG. 18 (d) is an explanatory diagram when setting a sampling pixel row for the intersection point.
  • FIG. 19 is an explanatory diagram of edge detection in the fifth embodiment.
  • FIG. 20 is a diagram for explaining the lateral chromatic aberration of the lens.
  • FIG. 21 is a diagram for explaining color misregistration due to lateral chromatic aberration.
  • Figure 22 shows the color information of R G B and Y U V color in consideration of false color suppression.
  • Figure 23 shows the relationship between the X and u coordinates.
  • Fig. 24 is a diagram showing the equations used when performing blur correction.
  • Fig. 25 is a diagram showing the equations used to perform coordinate conversion between color planes.
  • Figure 26 shows the relational expression of the pixel coordinates on the color mosaic image corresponding to the xy coordinates on each color plane.
  • FIG. 27 is a diagram showing an expression used when obtaining the chromatic aberration matrix.
  • FIG. 1 is a diagram showing a configuration example of a color imaging apparatus 100 according to the first embodiment in which the image processing apparatus according to the present invention is implemented.
  • the color imaging device 100 of this embodiment includes an imaging optical system 110, a single-plate color image sensor 120, an AD conversion unit 130, a demosic unit 140, and a visual correction unit 15. 0, a compression unit 1600 and a recording unit 1700.
  • the demosaicing unit 140 corresponds to the image processing apparatus of the present invention.
  • the single-plate color image sensor 1 2 0 has a color filter array 1 2 1 that separates the imaged light output from the imaging optical system 1 1 0 into predetermined color components, and a force filter array. And an image sensor 1 2 2 that generates a pixel signal by photoelectrically converting the imaged light that has passed through 1 2 1.
  • the image sensor 1 2 2 is composed of, for example, a CCD (Charge Coupled Device), and a certain layer is a CMOS (Complementary Metal Oxide Semiconductor) sensor.
  • the imaging optical system 1 1 0 plays a role of guiding imaging light of a subject to the single plate color image sensor 1 2 0.
  • the imaging optical system 110 includes an imaging lens, an infrared ray removal filter, and the like, including an optical low-pass filter.
  • the infrared filter is a single plate color image.
  • the optical aperture is used to block infrared rays incident on the sensor 120.
  • the single plate force color image sensor 1 2 0 color-in-rectifier 1 2 1 is regularly arranged in a predetermined pattern on the light-receiving surface of each pixel constituting the image sensor 1 2 2.
  • a color filter array 1 2 1 having a primary color bay arrangement using three colors RG and B is used as the color component. I'm going.
  • the primary color bay arrangement is an arrangement of G color filters in Himatsumatsu, and R color filters and B color filters arranged alternately in each row. .
  • G color filters those sandwiched between horizontal R color filters are called Gr color filters, and those sandwiched between horizontal B color finoleta are called Gb color filters.
  • the fe image element 1 2 2 plays a role in photoelectrically converting the received imaging light into electrical pixel information and storing it as a charge amount, and outputting it as an air signal to the AD conversion unit 1 3 0 Is.
  • the image pickup device 1 2 2 has a plurality of pixels (a pixel diode) arranged in a predetermined pattern, and the color filter array 1 2 1 is arranged on the light receiving surface of each pixel. Regularly arranged in an array
  • a subject that has been subjected to the above-described configuration is imaged on the image sensor 1 2 2 of the single plate force image sensor 1 2 0 via the ⁇ image optical system 1 1 0.
  • the imaged subject image is deteriorated. For example, a straight line on the subject becomes a curved image due to distortion.
  • the single-plate color image sensor 120 converts the object image formed on the image sensor 1222 into an analog electrical signal as a color mosaic image.
  • the color filter array 1 2 1 has a primary color bay arrangement as shown in FIG. 2, and R, G, and B color filters are arranged for each pixel of the image sensor 1 2 2. . Therefore, the imaging light of the subject reaches the imaging device 1 2 2 as imaging light in a color mosaic image state in which only the color component corresponding to each pixel is transmitted.
  • the image sensor 1 2 2 photoelectrically converts the reaching light and outputs it as an electric signal of the color mosaic image to the AD converter 1 30.
  • the AD converter 1 30 converts the analog signal of the color mosaic image photoelectrically converted by the image sensor 1 2 2 into a digital signal so that digital signal processing can be performed.
  • the color mosaic image immediately after AZD conversion by the AD conversion unit 130 is also called R AW data.
  • the demosaicing unit 140 converts the color mosaic image into a color image. In this embodiment, the image deformation process is performed at the same time, thereby correcting the image quality deterioration due to the distortion aberration of the imaging optical system 110 described above. The method of image processing by the demosic unit 140 will be described in detail later.
  • the visual correction unit 150 mainly performs processing for improving the appearance of the surface image on the color image generated by the demosic unit 140.
  • the visual correction unit 1 5 0 performs image correction processing such as tone curve (gamma) correction, saturation enhancement, and edge enhancement.
  • the compression unit 160 compresses the color image corrected by the visual correction unit 150 using a method such as JPEG (Joint Photographic Experts Group) to reduce the recording size.
  • the recording unit 170 records the compressed digital image signal on a recording medium (not shown) such as a flash memory.
  • Each configuration from the demosaicing unit 14 0 to the recording unit 1 70 may be configured as a separate device, or may be configured as a single micro processor. good. In the latter case, a single microprocessor executes processing related to each configuration from the demographic unit 14 0 to the recording unit 1 70
  • Fig. 3 is a block diagram showing an example of the functional configuration of demosaicing unit 140.
  • FIG. 4 is a flowchart showing an example of the image processing executed by the demosaicing unit 140.
  • FIG. 5 is an image diagram for specifically explaining the contents of the image processing executed by the demosaicing unit 140.
  • the demosaicing unit 140 has a functional configuration of a color plane resolution unit 14 1, a coordinate conversion unit 1 4 2, a sampling unit 1 4 3, And a color generation unit 1 4 4.
  • the color plane separation unit 1 4 1 decomposes the color mosaic image output from the AD conversion unit 1 3 0 into a plurality of color planes including only pixel values of the same color light (step in FIG. 4). S 1).
  • S 1 the color plane separation unit 1 4 1 decomposes the color mosaic image output from the AD conversion unit 1 3 0 into a plurality of color planes including only pixel values of the same color light.
  • S 1 decomposes the color mosaic image output from the AD conversion unit 1 3 0 into a plurality of color planes including only pixel values of the same color light.
  • the coordinate conversion unit 14 2 uses a coefficient representing the image deformation of the color image, and uses the color image generated from the color mosaic image to display the color image (when the deformation process is performed)
  • the sampling coordinates on the color mosaic image corresponding to the pixel position of the deformed color image are calculated (step S 2 in FIG. 4).
  • the distorted color generated from the color mosaic image is used.
  • the procedure for calculating the sampling coordinates is described in detail below. First, for the xy coordinate system, the origin is the image center, the maximum image height (maximum distance from the origin) is 1, the positive X coordinate is taken from the origin to the right of the screen, and the origin is positive from the origin to the bottom of the screen. The y coordinate of is taken.
  • the pixel coordinate (u d , V d ) of the uv coordinate system is set to (0) from the upper left corner of the screen as shown in FIG. , 0), (1, 0), (2, 0) ..., if the next line is assigned as (1, 0), (1, 1), (2, 1) ... 'uv coordinates
  • the pixel coordinates of the system (3 19.5, 2 39.5) are the origin of the xy coordinate system.
  • ⁇ k ,, k 2 ⁇ are coefficients indicating distortion aberration of the imaging optical system 110, k, is a third-order aberration coefficient, and k 2 is a fifth-order aberration coefficient.
  • a method for performing nonlinear coordinate transformation on a color image including distortion is already known, and the distortion coefficient ⁇ k ,, k 2 ⁇ is obtained by simulation or the like. Is possible.
  • This distortion coefficient ⁇ k ,, k 2 ⁇ corresponds to the image deformation coefficient of the present invention.
  • the pixel coordinates (u s , VJ are not necessarily integer values, but are generally non-integer.
  • the pixel coordinates (U s , V s ) on this color mosaic image are the sampling coordinates.
  • one of these sampling coordinates is indicated by the symbol 20 0.
  • the color mosaic image is broken down into three color planes.
  • Figure 5 shows the sampling coordinates 2 0 0 on each color plane, as mentioned above, since the sampling coordinates are non-integer values, the sampling coordinates are shifted to the position shifted from the pixel center. Coordinate 2 0 0 exists.
  • the sampling unit 1 4 3 is a pixel at the sampling coordinates 2 0 0 calculated by the coordinate conversion unit 1 4 2 for each of a plurality of color planes decomposed by the color plane decomposition unit 1 4 1.
  • the value is interpolated from the pixel values of the same color light contained in the color plane (step S3 in Fig. 4). That is, the sampling unit 14 3 calculates the pixel value of the sampling coordinate 2 0 0 from each of the R plane, the G plane, and the B plane, and outputs it.
  • the value (u s , V s ) of the sampling coordinate 2 0 0 is not necessarily an integer value, so there are four possible values surrounding the sampling coordinate 2 0 0.
  • Linear interpolation is performed from pixels (pixel values of the same color light that each color plane originally has). This linear interpolation is preferably performed by bilinear interpolation.
  • the R and B planes have valued pixels in the form of vertical and horizontal grid points, so the four valued pixels surrounding the sampling coordinates 200 are the corresponding sampling coordinates 2 It is located at each vertex of a square with a side length of 2 surrounding 0 0.
  • Respective pixel values of the relevant pixel are represented by R (1 0 0, 1 0 0), R (1 0 0, 1 0 2), R (1 0 2, 1 0 0), R (1 0 2, 1 0 2), the interpolated pixel value R (1 0 0. 8, 1 0 1.4) of sampling coordinates 2 0 0 on the R plane generated by bilinear interpolation as shown in Fig. 6 ) Is expressed by the following equation.
  • the four value pixels surrounding the sampling coordinate 2 00 have a side length of 2 surrounding the sampling coordinate 2 0 0. 4 Located at each vertex of a square tilted 5 degrees.
  • Each pixel value of the relevant pixel is set to G (1 0 0, 1 0 1), G (1 0 1, 1 0 0), G (1 0 1, 1 0 2), G (1 0 2, 1 0
  • the interpolated pixel value G (10.8.8, 10.4.4) of the sampling coordinate 2200 on the G plane is expressed by the following equation.
  • G (1 0 0 .8, 1 0 1 .4) 0 .7 * 0 .3 * G (1 0 0, 1 0 1) + 0. 3 * 0 .3 * G (1 0 1, 1 0 0) + 0. 7 * 0. 7 * G (1 0 1, 1 0 2) + 0. 3 * 0. 7 * G (1 0 2, 1 0 1)
  • the color generation unit 14 4 synthesizes the interpolated pixel values of each color plane interpolated by the sampling unit 14 3, so that each surface element value has luminance information of multiple colors.
  • a color image is generated (step S 4 in FIG. 4).
  • the interpolation value of each color plane obtained by the sampling unit 14 3 for one sampling coordinate (u s , V s ) is respectively represented by R (u s) V s ), G (u s , v J, B (u s , v s ), the color generator 1 4 4 combines these three interpolation values to obtain the sampling coordinates (u s , V s ). Generate color information.
  • the color generation unit 14 4 converts the RGB color information thus obtained into YUV color information (Y is luminance information, U and V are color information), and U and V color information is obtained.
  • a low frequency filter is applied to. Note that from RGB to YUV A known process can be applied to the conversion process and the low frequency filter process for U and V.
  • the color generation unit 144 performs the above processing for all the pixels (all sampling coordinates) of the deformed color image, and outputs the resulting deformed color image to the visual correction unit 150.
  • the processing after the visual correction unit 150 is as described above.
  • the color mosaic corresponding to the pixel position of the modified color image to be output is generated as the previous stage of generating the deformed color image from the color mosaic image.
  • the sampling coordinates 2 0 0 on the image are calculated.
  • the pixel values of the deformed color image are converted into the sampling coordinates 2
  • the pixel value of the color mosaic image is used as the pixel value at 0, and is generated by interpolation calculation.
  • the color interpolation processing for generating a power error image from the color mosaic image and the deformation processing of the color image can be realized by a single interpolation operation.
  • the processing load when generating a deformed color image from a color mosaic image can be reduced, and as before, the degradation of image quality due to repeated interpolation processing is also suppressed. can do.
  • the color mosaic image is decomposed into a plurality of color planes of R, G, and B, and interpolation pixel values of sampling coordinates 200 are obtained for each color plane.
  • the interpolated pixel values are combined to generate color information that includes luminance information of three colors in one pixel. By doing so, it is possible to obtain the catch pixel value of sampling coordinates 200 by simple linear interpolation from the pixel values of the same color light contained in the color plane, Processing load can be reduced.
  • the color plane separation unit 14 1 separates the color mosaic image into three color planes of the R plane, the G plane, and the B plane. Not limited to this.
  • the two kinds of green pixels G r and G b that are in contact with each other may be separated into different color planes.
  • the color plane separation unit 14 1 extracts only the R plane and G r component pixels from the color mosaic image output from the AD conversion unit 1 30, which extracts only the R component pixels.
  • the image is divided into four color planes: the extracted Gr plane, the Gb plane from which only the Gb component pixels are extracted, and the B plane from which only the B component pixels are extracted.
  • the color generation unit 14 4 uses the pixel value (sampling value) of the sampling coordinates 2 0 0 generated by the sampling unit 1 4 3 for the G r plane and the G b plane, respectively. Is added to generate color information in which luminance information of R, G, and B color components is contained in one pixel.
  • the sampling value is used as it is for the R component and B component, and the average value of G r and G b is used for the G component.
  • the demosaic unit 1 4 0 is connected to each color of the G r plane and G b plane.
  • the image processing apparatus further includes a false color determination unit that calculates a difference between the interpolation pixel values generated by the sampling units 1 to 4 on the plane and determines the presence or absence of a false color based on the difference between the interpolation pixel values. It's okay.
  • the interpolated pixel values obtained from the Gr plane and Gb plane should be the same value for both. It is. However, if a false color is generated, a difference occurs in the interpolated pixel values obtained by the Gr plane and Gb plane. Therefore, it is possible to detect the presence or absence of false color by looking at the difference between the interpolated pixel values. If it is determined that there is a false color, the color generator 1 4 4 uses the formula in Fig. 2 2 considering the suppression of false color when converting RGB power error information into YUV color information. Perform the conversion process.
  • the distortion correction is described as an example of the deformation process for the color image.
  • the present invention is not limited to this.
  • Demosaicing part 1 4 0 can be applied.
  • enlargement, reduction, and rotation can be expressed by, for example, affine transformation, and a coefficient representing the affine transformation is used as an image deformation coefficient.
  • FIG. 8 is a diagram showing a configuration example of a color imaging device 1 0 0 ′ according to the second embodiment in which the image processing device according to the present invention is implemented.
  • the color imaging device 1 0 0 ′ of the second embodiment includes an imaging optical system 1 1 0, a single plate color image sensor 1 2 0, and an AD conversion unit 1 3 0.
  • the demosaic unit 140 'corresponds to the image processing apparatus of the present invention.
  • those having the same reference numerals as those shown in FIG. 1 have the same functions, and therefore redundant description is omitted here.
  • the photographed subject forms an image on the image sensor 1 2 2 of the single plate color image sensor 1 2 0 via the imaging optical system 1 1 0.
  • Demosaicing unit 1 4 0 'converts color mosaic image to color image.
  • the color chromatic aberration due to the magnification chromatic aberration of the imaging optical system 110 described above is corrected by performing a correction process for the magnification chromatic aberration at the same time.
  • FIG. 9 is a block diagram showing an example of a functional configuration of the demosaic unit 140 '.
  • FIG. 10 is a flowchart showing an example of the image processing operation executed by the demosaic unit 140 0 ′.
  • FIG. 11 is an image diagram for specifically explaining the contents of the image processing executed by the demosaicing unit 140.
  • the demosaicing section 1 4 0 ′ is composed of a color plane separation section 1 4 1, a coordinate conversion section 1 4 2, and a sampling section 1 as its functional configuration.
  • the color plane separation unit 14 1 1 is the R plane that extracts only the R component pixels from the color mosaic image output from the AD conversion unit 1 3 0. , G plane extracted only G component pixels, B pixel extracted only ⁇ Plane separated into 3 color planes (step S1 1 in Fig. 10) ⁇
  • the coordinate conversion unit 1 4 2 ′ uses the pixel position on the color mosaic image corresponding to the pixel position on the power raster image when chromatic aberration correction is performed from the pixel position on the power raster image generated from the color mosaic image.
  • Calculate the sampling coordinates (step S 1 2 in Fig. 10).
  • the coordinate conversion unit 1 4 2 ′ has a chromatic aberration coefficient having a different value depending on the color plane for each of the plurality of color planes decomposed by the color plane decomposition unit 1 4 1.
  • Using calculate different sampling coordinates for each color plane.
  • the imaging optical system 1 1 0 is used for color mosaic images. In order to prevent this effect from appearing in the color image, the sampling coordinates should be given different values for the color planes corresponding to different color lights.
  • coordinate conversion taking into account the correction of the chromatic aberration of magnification of the imaging optical system 110 is performed as follows.
  • coordinate conversion is performed separately for each of the three color planes of the R, G, and B planes.
  • ⁇ k R , k B ⁇ are coefficients indicating the chromatic aberration of magnification of the imaging optical system 110, k R is the magnification of the R component with respect to the G component, and k B is the magnification of the B component with respect to the G component. ing.
  • chromatic aberration coefficients ⁇ k R , k B ⁇ can be obtained by optical simulation.
  • the color mosaic image is assumed to be composed of 16 0 0 X 1 2 0 0 square pixels, and the pixel coordinates (u s , V s ) of the uv coordinate system are allocated in the same manner as in the color image described above.
  • the pixel coordinates of the uv coordinate system (79.9.5, 59.9.5) are the origin of the xy coordinate system
  • the pixel coordinates (U sR , V sR ) (U sG) V sG ) (U sB , V sB ) on this color mosaic image are the sampling coordinates for each color plane.
  • the sampling coordinates on the R plane are denoted by reference numeral 30 1
  • the sampling coordinates on the G plane are denoted by reference numeral 30 2
  • the sampling coordinates on the B plane are denoted by reference numeral 30.
  • the sampling coordinates 3 0 1, 3 0 2, 3 0 3 of each color plane exist at positions shifted from the pixel center. .
  • the sampling unit 1 4 3 ′ has a sample for each color plane calculated by the coordinate conversion unit 1 4 2 ′ for each of the color planes decomposed by the color plane separation unit 1 4 1.
  • the pixel values (sampling values) at the sampling coordinates 3 0 1, 3 0 2, 3 0 3 are interpolated from the pixel values of the same color light contained in the color plane (step of Fig. 10). S 1 3).
  • the sampling unit 1 4 3 ′ calculates the pixel values of the sampling coordinates 3 0 1, 3 0 2, and 3 0 3 through interpolation calculation using the power of the R plane, G plane, and B plane. Output.
  • the values of sampling coordinates 3 0 1, 3 0 2, 3 0 3 (u sR ,
  • V sR (U sG , V sC ) (U sB) V sB ) is not necessarily an integer value, so the four value pixels (3 0 1, 3 0 2, 3 0 3) surrounding the sampling coordinates 3 0 1, 3 0 2, 3 0 3 Interpolation is performed from the pixel values of the same color light that each color plane originally has). This interval is preferably performed by bilinear interpolation.
  • the R plane and B plane have value pixels in vertical and horizontal grid points, so the four value pixels surrounding the sampling coordinates 3 0 1, 3 0 3 are It is located at each vertex of a square with a side length of 2 surrounding the coordinates 3 0 1 and 3 0 3.
  • the four pixels surrounding this in the R plane ( u d , V d ) (1 0 0, 1 0 0)
  • (1 0 0, 1 0 2), (1 0 2, 1 0 0), (1 0 2, 1 0 2) are the R pixel valued pixels.
  • the G plane has value pixels in a checkered pattern.
  • the four valued pixels surrounding 30 2 are located at the vertices of a square tilted 45 degrees with a side length of “/ 2” surrounding the sampling coordinate 30 2. This ⁇ ,
  • the four pixels (u d , V d,) (1 0 0, 1 0 1), (1 0 1, 1 0 0), (1 0 1, 1 0 2), (1 0 2, 1 0 1) are the value pixels of the G plane.
  • Each pixel value of the relevant pixel is represented by G (1 0 0, 1 0 1), G (1 0 1, 1 0 0), G (1 0 1, 1 0 2), G (1 0 2, 1 0 If it is expressed by 1), the interpolated pixel value G (1 0 1.0, 1 0 1.4) of the sampling coordinates 30 2 on the G plane is expressed by the following equation.
  • G (1 0 1. 0, 1 0 1. 4) 0. 7 * 0. 3 * G (1 0 0, 1 0 1) + 0. 3 * 0. 3 * G (1 0 1, 1 0 0) + 0. 7 * 0. 7 * G (1 0 1, 1 0 2) + 0. 3 * 0. 7 * G (1 0 2, 1 0 1)
  • the color generation unit 14 4 synthesizes the interpolated pixel values of the color planes interpolated and generated by the sampling unit 1 4 3 ′, so that each pixel value represents luminance information of multiple colors.
  • a color image is generated (step S 1 4 in FIG. 10). Further, the color generation unit 14 4 4 converts the RGB color information thus obtained into YUV color information (Y is luminance information, U and V are color information), and U and V colors. Apply a low frequency filter to the information. It should be noted that known processing can be applied to the conversion processing from R GB to Y U V and the low-frequency filter processing for U and V.
  • the color generation unit 144 performs the above processing for all pixels (all sampling coordinates) of the color image, and outputs the resulting deformed color image to the visual correction unit 150. The processing after the visual correction unit 150 is as described above.
  • Interpolated pixel values R (u sR , V sR ), G (u sG , v sG ), ⁇ (u sB ) of sampling coordinates 3 0 1, 3 0 2, 3 0 3 used in the processing of the color generator 14 4 , ⁇ sB ) takes into account the color shift of the chromatic aberration of magnification and shows the same part on the subject. For this reason, R, G, and B change at the same time where there is a black and white border on the subject. This makes it possible to obtain a sharp luminance signal for a color image obtained by combining the R, G, and B color components. That is, Capable of obtaining a force image equivalent to that captured and shaded by an imaging optical system without lateral chromatic aberration
  • the color position is determined from the pixel position of the output force level image. Sampling coordinates that are positions on the color mosaic image that correspond to the pixel positions of the image and that are adjusted in consideration of lateral chromatic aberration
  • 3 0 1, 3 0 2 and 3 0 3 are calculated for each color plane.
  • the pixel values at the sampling coordinates 3 0 1, 3 0 2, 3 0 3 are generated by interpolation using the pixel values of the color mosaic image.
  • the color interpolation processing for generating a color image from the color mosaic image and the magnification chromatic aberration correction processing for the power error image can be realized by a single interpolation calculation. For this reason, the processing load when generating a color image corrected for chromatic aberration from a color mosaic image can be reduced, and the interpolation process can be performed twice as before. Image quality deterioration due to
  • the color mosaic image is decomposed into a plurality of color planes R, G, and B, and the sampling coordinates 3 0 1, 3 0 2,
  • the interpolated pixel values of each color plane are synthesized to generate color information in which luminance information of three colors is included in one pixel.
  • the sampling coordinates 3 0 1, 3 0 2, and 3 are obtained by simple linear interpolation from the pixel values of the same color light contained in the color plane.
  • the interpolated pixel value of 3 0 3 can be obtained, and the processing load can be reduced.
  • the color plane decomposition unit 14 1 1 decomposes the color mosaic image into three color planes of R plane, G plane, and B plane.
  • R plane R plane
  • G plane G plane
  • B plane B component images
  • the color generator 1 4 4 has sampled coordinates 3 0 generated by the sampler 1 4 3 'in the G ⁇ plane and G b plane, respectively.
  • the two pixel values are added to generate color information in which luminance information for each color component of R, G, and B is included in one pixel IV.
  • R component For example, R component
  • the B component uses the sampling value as it is, and the G component uses the average value of G r and G b o
  • the amosaic part 1 4 0 ' is used for the G r plane and the G b plane. It further includes a false color determination unit that calculates the difference between the capture pixel values generated by the sampling unit 1 on each color plane and determines the presence or absence of a false color based on the difference between the interpolation pixel values. You may do it.
  • the difference between the interpolation pixel value of the sampling coordinate 30 2 obtained on the G r plane and the interpolation pixel value of the sampling coordinate 3 0 2 obtained on the G b plane is taken.
  • G r and G b are essentially the same G color finisher, so G ⁇ pre-
  • the interpolated pixel values found in the G and Gb planes should be the same for both. However, if a false color is generated, there will be a difference in the interpolated pixel values found in the Gr plane and Gb plane. Therefore, by looking at the difference between the interpolated pixel values, it is possible to detect the presence or absence of false colors.
  • the conversion process is performed using the formula in Fig. 22.
  • the coordinate conversion unit 1 4 2 ′ is a color plane decomposition unit.
  • a color image pixel is obtained by using a chromatic aberration coefficient having a different value depending on the color plane and an image deformation coefficient representing an image deformation of the power error image. From the position, the sampling coordinates on the color mosaic image corresponding to the surface element position of the color image when chromatic aberration correction and image deformation are performed, and different sampling coordinates for each color plane are calculated.
  • a straight line on the subject becomes a curved image due to distortion.
  • image quality degradation due to distortion is also corrected.
  • the distorted color image generated from the color mosaic image is corrected.
  • a method for performing non-linear coordinate transformation on a color image including distortion is already known, and the above distortion aberration coefficient ⁇ k ,, k 2 ⁇ can be obtained by simulation or the like. .
  • This distortion coefficient ⁇ k ,, k 2 ⁇ corresponds to the image deformation coefficient of the present invention.
  • the sampling unit 14 3 ′ uses the pixel value (sampling value) at the sampling coordinates for each color plane calculated by the coordinate conversion unit 14 2 ′ as described above as the color plane. Interpolation is generated from the pixel values of the same color light contained within. Also, the color generation unit 14 4 synthesizes the interpolated pixel values of each color plane interpolated by the sampling unit 1 4 3 ′ so that each pixel value has luminance information of a plurality of colors.
  • image optics is used as an example of image deformation.
  • image deformation can be achieved by enlarging or reducing a color image using the Digitals Theorem, and power It may be an image rotation.
  • enlargement, reduction, and rotation can be expressed by, for example, an Affine transformation, and a coefficient representing the Affine transformation is used as an image deformation coefficient.
  • FIG. 6 is a diagram showing a configuration example of a force-ra-fe image device 100, according to a fourth embodiment in which an image processing device according to the present invention is implemented.
  • FIG. 13 those given the same reference numerals as those shown in FIG. 8 have the same functions, so redundant description is omitted here.
  • the force imager 100 of the fourth embodiment includes an imaging optical system 1 1 0 ”, a single plate color image sensor 1 2 0, and an AD conversion unit 1.
  • demosic unit 14 40 visual correction unit 15 50, compression unit 16 60, recording unit 170, blur detection unit 180, aberration coefficient setting unit 190, carrier blurring It comprises a projection unit 20 0 and an aberration coefficient table storage unit 2 1 0.
  • the demosic unit 1400 and the aberration coefficient setting unit 1900 are images of the present invention. It corresponds to a processing device.
  • the photographed subject forms an image on the image sensor 1 2 2 of the single-plate force image sensor 1 2 0 through the image optical system 1 1 0 ".
  • / (3 ⁇ 4 image optical system 1 1 10 "has various aberrations, and the formed subject image is deteriorated.
  • a straight line on the subject due to distortion is a curved line. That is, an image in which image displacement (color misregistration) occurs on the image sensor 1 2 2 for each color component due to the chromatic aberration of magnification occurs in the image forming optical system 110 in the fourth embodiment.
  • the blur detection unit 1800 detects the blur of the color image pickup device 100 "and detects the blur.
  • the correction value of the subject image size due to blurring in the longitudinal direction of 0 is z
  • the correction value of the subject image rotation due to blurring of the mouth axis is 0
  • the correction value of the subject image rotation is
  • the correction value for the left-right position of the subject image is d X
  • the correction value for the vertical position of the subject image that accompanies vertical or horizontal blur is dy.
  • the coefficient z may include the digital zoom magnification.
  • the aberration coefficient setting unit 1 90 corresponds to the coefficient setting unit of the present invention. O The aberration coefficient setting unit 1 90 detects and detects the lens state of the imaging optical system 1 1 0 ". Aberration coefficient ⁇ kk 2 , p ,, P
  • 0 is a coefficient indicating distortion
  • ⁇ k ,, k 2 ⁇ are distortions in the radiation direction
  • P 1 and P 2 ⁇ indicate the tangential distortion.
  • ⁇ k R , k B ⁇ is the imaging optical system 1
  • k R denotes the ratio of B component to magnification, k component of the R component for the G component.
  • the aberration coefficient setting unit 190 detects the lens state after the change, and an aberration coefficient corresponding to the lens state is obtained from the aberration coefficient table storage unit 2 10. Reading, demosaicing part 1 4 0
  • the detection of the lens state is set to the image optical system 1 1 0 "from the controller (not shown) of the color imaging device 10 0" which controls the lens state, for example. It is possible to do this by receiving information on the lens states that have been registered.
  • the aberration coefficients for all lens states are recorded in a table. It is unrealistic to record. Therefore, in this embodiment, only the aberration coefficients for a limited lens state are recorded in the aberration coefficient table storage unit 2 10.
  • the focal length and subject distance are 3 patterns each, and a total of 9 lens states and corresponding aberration coefficient values are recorded.
  • An example of the part related to the aberration coefficient k, in this aberration coefficient table is as shown in FIG.
  • the aberration coefficient setting unit 190 has an aberration coefficient of 0.08 with an object distance of Mid: 1.0 m and a focal length of Wide: 28 mm from the aberration coefficient table.
  • Object distance is Far: Inf
  • focal length is Wide: 28 mm
  • aberration coefficient is 0.05
  • subject distance is Mid: 1.0 m
  • focal length is Mid: 5 O mm
  • k, ((1 / 2.0. 0-l Z lnf) (l / 3 5 _ l Z 5 0) * 0. 0 8 + (1 / 1.0. 0-1 / 2.0. 0) (1/3 5 -1/5 0) * 0. 0 5 + (1/2. 0-1/1 nf) (l / 2 8-l / 3 5) * 0. 0 2 + (1/1. 0 — 1 / 2.0) (1/2 8-1/3 5) * 0. 0 0) / (1 / 1.0. — L I nf) / (1/2 8-1/5 0)
  • the interpolated value calculated as shown in Fig. 6 is set in demosaicing section 1 4 0 ".
  • the interpolated values are calculated in the same way, Demoza Set to “1 4 0”.
  • the calibration unit 20 0 0 generates an aberration coefficient table value to be stored in the aberration coefficient table storage unit 2 1 0 based on the digital image signal output from the AD conversion unit 1 3 0. .
  • the calibration unit 20 0 assigns the yield coefficient ⁇ k ⁇ k 2 , p p p 2 , k R , k B ⁇ to each of the plurality of lens states to be recorded in the yield coefficient table. These are obtained and recorded in the aberration coefficient table storage unit 2 10.
  • the calibration unit 20 0 may be provided as a part of the color imaging device 10 0 ", but it is separate from the color imaging device 1 0 0". It may be configured as a revision device. When a separate device is used, the force imager device 100 “0” outputs the photographed image (digital image signal output from the AD conversion unit 130) to the calibration device, Communication means for inputting aberration coefficients from a calibration device.
  • the demosaicing unit 1 4 0 is a blur correction coefficient ⁇ z, 0, dx, dy ⁇ set by the shake detection unit 1 80 and an aberration coefficient set by the aberration coefficient setting unit 1 90 ⁇ k have k 2, p have p 2, k R, k B ⁇ while performing correction based on, performs conversion into a color image from the color Moza I click the image.
  • Fig. 15 is a block diagram showing an example of the functional configuration of demosaic unit 1 4 0 ". As shown in Fig. 1 5, demosaic unit 1 4 0" is shown as its functional configuration. A color plane separation unit 1 4 1 ′′, a coordinate conversion unit 1 4 2 ′′, a sampling unit 1 4 3 ′′, and a color generation unit 1 4 4 ′′.
  • the color plane decomposing unit 14 1 decomposes the color mosaic image output from the AD converting unit 13 30 into a plurality of color planes including only pixel values of the same color light.
  • the color plane separation unit 1 4 1 is an R plane that extracts only R component pixels, a G r plane that extracts only G r component pixels, and G G b plane that only extracted b component pixels, Separate only the B component pixels into four color planes, the B plane.
  • Coordinate conversion unit 1 4 2 " the above-described blur correction coefficient ⁇ z, 0, dx, dy ⁇ and aberration coefficients ⁇ kk 2, pp 2, k K, k B ⁇ with, from the color mosaic image From the pixel position of the generated color image, the sampling coordinates on the color mosaic image corresponding to the pixel position of the color image when chromatic aberration correction, distortion aberration correction and blur correction are applied are calculated.
  • the following is a detailed explanation of the sampling coordinate calculation procedure.
  • the origin is the image center
  • the maximum image height (maximum distance from the origin) is 1
  • the positive X coordinate is taken from the origin to the right of the screen
  • the origin is positive from the origin to the bottom of the screen.
  • the y coordinate of is taken. It is also assumed that the color image to be output consists of 1600 0 X 1 2 0 0 square pixels.
  • the pixel coordinates (u d , V d ) of the uv coordinate system are (0, 0), (1, 0), (2, 0) ⁇
  • the next row is assigned to (1, 0)
  • corresponding to pixel coordinates (u d , V d ) xy coordinates (x d, y d) are you express by equation shown in FIG 3.
  • Coordinate conversion unit 1 4 2 " apply the xy coordinate (x d, y d) with respect to, first, the above-described blur correction coefficient ⁇ z, 0, dx, dy ⁇ to cormorants good expression shown in FIG. 2 4 Then, the xy coordinates ( ⁇ ', y') after blur correction are obtained.
  • the coordinate transformation unit 14 2 "applies the coefficients ⁇ k,, k 2 , P l , p 2 ⁇ related to distortion as shown in the following equation, and X in the G r plane and G b plane Find the y coordinate (X G , y G ).
  • X G X '(l + k, r 2 + k 2 r 4 ) + 2 p 1 x' y '+ p 2 (r' 2 + 2 x ' 2 )
  • the coordinate conversion unit 14 2 "converts the coordinates between the color planes in consideration of the coefficients ⁇ k K , k B ⁇ relating to the magnification color yield of the imaging optical system 1 1 0".
  • X y coordinates in R and B planes (X K, y R) ( x B, y B) , respectively.
  • d R x, d R y is the parallel displacement of the R plane when the G plane is the reference
  • d B x, d B y is the parallel displacement of the B plane when the G plane is the reference.
  • the color mosaic image is also a square pixel of 16 0 0 X 1 2 0 0, and the pixel coordinates (u s , V s ) of the uv coordinate system are divided in the same way as the color image described above.
  • these pixel coordinates ( usR , VsR ) ( usG , vsG ) ( usB , vsB ) are not necessarily integer values, and are generally non-integer values.
  • the pixel coordinates ( usR , VsR ) ( usG) vssG ) ( usB , vsB ) on this color mosaic image are the sampling coordinates for each color plane.
  • the sampling unit 1 4 3 is provided for each color plane calculated by the coordinate conversion unit 1 4 2" for each of the plurality of color planes decomposed by the color plane separation unit 1 4 1 ".
  • the pixel values (sampling values) in the sampling coordinates of the same are generated from the pixel values of the same color light contained in the color plane.
  • the color generator 1 4 4 is the sampling unit 1 4 3 By synthesizing the interpolated pixel values of each color plane generated by “”, a color image in which each pixel value has brightness information of a plurality of colors is generated. This sample
  • the processing contents of the printing unit 14 3 "and the color generation unit 14 4" are the same as those in the second embodiment or the third embodiment.
  • the color mosaic image is displayed.
  • Color capture processing that generates a color image from the image, magnification chromatic aberration correction of the imaging optical system 1 1 0 ", deformation of the color image to correct distortion aberration of the imaging optical system 1 1 0
  • correction processing for correcting aberrations according to the lens state of the imaging optical system 110, and processing for correcting camera shake of the color imaging device 100 can be realized by a single interpolation operation. For this reason, it is possible to reduce the processing load when generating a color image from a color mosaic image in which the chromatic aberration of magnification is corrected, and the force, distortion, distortion, and camera shake are also corrected. In addition, it is possible to suppress the deterioration of image quality due to the interpolation process being performed many times.
  • the blur correction coefficients ⁇ z, 0, dx, dy ⁇ and the aberration coefficient ⁇ k are used as the coefficients used by the coordinate conversion unit 14 2 ′′ to calculate the sampling coordinates.
  • the remaining coefficients can be used in any combination.
  • the present invention is not limited to this.
  • the value of the aberration coefficient table may be fixed without providing the calibration unit 200.
  • FIG. 16 is a diagram showing a configuration example of a color imaging device 100 A according to the fifth embodiment.
  • the color imaging device 1 OOA according to the fifth embodiment includes a chromatic aberration amount detection device 20 OA that detects a chromatic aberration coefficient by photographing a predetermined chart CH, and is detected by the chromatic aberration amount detection device 20 OA.
  • the chromatic aberration coefficients are stored in the coordinate conversion unit 14 2 ′ in the second and third embodiments and the aberration coefficient setting unit 1 90 in the fourth embodiment, and the second to fourth embodiments.
  • the chart CH corresponds to the predetermined image in claims 23 and 24 of the present invention.
  • FIG. 17 is an explanatory diagram of the detection of the amount of chromatic aberration in the embodiment
  • Fig. 17 (a) shows the chart form
  • Fig. 17 (b) shows the chart for the imaging element 1 2 2. It is a figure showing arrangement
  • FIG. 18 is an explanatory diagram of chromatic aberration amount detection in the same embodiment
  • FIG. 18 (a) and (b) are explanatory diagrams when detecting intersection points
  • FIG. 18 (c) is an edge for each intersection point.
  • FIG. 18 (d) is an explanatory diagram when setting a sampling pixel row with respect to the intersection
  • FIG. 19 is an explanatory diagram of edge detection in the same embodiment.
  • the chromatic aberration amount detection device 2 0 0 of the present embodiment First, as shown in FIG. 16, the chromatic aberration amount detection device 2 0 0 of the present embodiment.
  • A calculates the chromatic aberration coefficient of the imaging optical system 110 using a digital image signal obtained by photographing the chart CH.
  • the chart CH is the first image pattern P 1 with respect to the pixel array of the image sensor 1 2 2 in the single plate force image sensor 1 2 0.
  • the second image pattern P 2 are tilted by the tilt angle ⁇ .
  • one image pickup pattern area is an image pickup element 1 2.
  • the chromatic aberration amount detection device 20 OA receives a digital image signal (a so-called pixel signal representing the luminance of the pixel) input from the AD conversion unit 13 30. Based on the pixel signals stored in the field memory 2 2 1 and field memory 2 2 1 stored for each RGB color, the first imaging pattern P 1 and the second imaging pattern P 2 Detects the edge positions of the first image pattern P 1 and the second image pattern P 2 for each RGB around the intersection point detected by the intersection detection processing unit 2 2 5, which detects multiple intersections RGB edge detection processing unit 2 2 6, edge position storage unit 2 2 8 that stores edge positions detected by RGB edge detection processing unit 2 2 6 in association with intersections, edge position storage unit 2 2 8 Aberration coefficient calculation unit 2 2 9 for calculating chromatic aberration coefficient from edge position stored in CPU 2, CPU (Central Processing Unit) 2 3 0 , ROM (Red Only Memory) 2 3 1, etc., and CPU 2 3 0 performs each process of the chromatic aberration amount
  • the field memory 2 2 1 is associated with the Bayer array and stores an R field memory 2 2 2 that stores a red (R) pixel signal and a G field that stores a green (G) pixel signal.
  • the intersection detection processing unit 2 2 5 calculates a luminance gradient using pixel values in a predetermined range centered on the target pixel, and the luminance gradient is The position of the largest target pixel is detected as the intersection point Int.
  • Fig. 18 (b) five vertical and horizontal pixels centered on the pixel of interest are set, and the intersection position is detected by adding a weight according to the pixel position. In other words, the vertical and horizontal pixel values centered on the target pixel are multiplied by the coefficients shown in Fig. 18 (b) and the results are totaled.
  • the absolute value is the evaluation value of the pixel of interest
  • the position of the pixel of interest when the evaluation value exceeds the predetermined threshold is the intersection point Int
  • multiple intersection points Int is detected in a matrix.
  • the first imaging pattern P 1 and the second imaging pattern P 2 are arranged so that the intersection points Int appear in a matrix shape at equal intervals.
  • the RGB edge detection processing unit 2 2 6 has a plurality of pixel rows H s, V located vertically and horizontally via the intersection I nt for each RGB color. s is scanned with a predetermined sampling line length, and pixel values are sequentially acquired, and the sampling position where the change amount of the pixel value is the largest relative to the adjacent sampling position is detected as an edge. .
  • the luminance (pixel value) of each pixel was obtained for each sampling as represented by the curve In in FIG. 19, and was obtained by sampling as represented by the curve SL.
  • the amount of change (gradient SL) of the pixel value is calculated based on the pixel value, and the position EP where the largest amount of change (gradient SL) appears is detected as an edge.
  • Sampling is performed for each pixel of the same color and is sampled along the left-right direction.
  • sampling H s as shown in Fig. 17 (b)
  • the sampling length SL (1 1) in the left-right direction and the number of samples in the vertical direction of the sampling are shown.
  • the ring number SN (4) is predetermined according to the required detection accuracy. Also, when performing sampling V s along the vertical direction, the sampling length and the number of samplings in the upper and lower directions are determined in advance.
  • the interval S can be obtained geometrically as shown in Fig. 18 (d), for example.
  • the distance S from the intersection I nt is L-SL no 2
  • S (W + L) X
  • the edge blur E, the inclination angle ⁇ , and the sampling line length SL are known quantities. It can be obtained by the formula of tana.
  • the start position of the sampling H sl does not fall within the imaging pattern P 1-2 in Fig. 18 (d), but is separated from the P 1-2 by the edge blur amount E. In order to do so, it is better to find the interval S.
  • the edge position storage unit 2 2 8 detects the color detected by the RGB edge detection processing unit 2 2 6 at each intersection I ntj (j is a serial number 1, 2... Assigned to each intersection). Based on each edge position, the edge positions of G (green), R (red), and B (blue) are associated with the left and right direction (u direction) and the up and down direction (direction), respectively. (U ci, G j, (. U j, R j), (U j, B j)
  • aberration coefficient calculation unit 2 2 9 d stored in or falling edge of di position storage unit 2 2 8 tree located U G j, V G j) ,, U R j, V R j,, U B j, j;
  • the chromatic aberration coefficients k R , k are calculated using this.
  • the original point is the center of the image
  • the maximum image height (maximum distance from the origin) is 1, and from the origin to the right of the screen.
  • the positive X coordinate is taken, and the positive y coordinate is taken from the origin downward in the screen.
  • the pixel coordinates (u d , V d ) of the uv coordinate system are displayed on the upper left corner of the screen as shown in FIG. 2 for a color image composed of 6 40 X 4 80 square pixels.
  • the next line is (1, 0), (1, 1), (2, 1) ...
  • the pixel coordinates (3 19.5, 2 39.5) in the uv coordinate system are used as the origin of the xy coordinate system, and the xy coordinates (u d , V d ) corresponding to the pixel coordinates (u d , V d ) are assigned for each color.
  • x d , y d ) is calculated.
  • the chromatic aberration matrix can be obtained by substituting the xy coordinates and aberration coefficients calculated by the aberration coefficient calculator 2 29 into the equation shown in Fig. 25. .
  • the matrix located on the rightmost side is a pseudo inverse matrix.
  • the edge position for each color can be detected in association with the matrix-like intersection position, and the detected color can be detected.
  • Chromatic aberration coefficients (k R , k B) can be calculated based on each edge position, and using the calculated chromatic aberration coefficients, color interpolation processing (demosaic processing) for generating a power image from a color mosaic image, The process to correct the chromatic aberration of magnification of the imaging optical system can be realized by a single interpolation operation.-3
  • the color filter array 1 2 1 1 is a primary color barrier array
  • the present invention is not limited to this.
  • a color filter array having a complementary color array may be used.
  • the imaging device 10 OA includes the chromatic aberration amount detection device 2 0 0, but the external device is not configured so that the chromatic aberration amount detection device 2 0 OA is integrated with the imaging device 1 0 OA.
  • the aberration coefficient setting unit 1 90 may be configured to be stored.
  • demosaic processing methods according to the first to fifth embodiments described above can be realized by any of a hardware configuration, a DSP (Digitital> 1 gna 1 Processor, and a software).
  • a DSP Digitital> 1 gna 1 Processor
  • a software For example, in the case of realizing the software, the demosaic units 1 4 0, 1 4 0 ′, 1 4 0 ′′ (image processing apparatus) of the first to fifth embodiments are Actually, computer
  • each of the first to fifth embodiments is merely an example of a specific example for carrying out the present invention, and the technical scope of the present invention is interpreted in a limited manner by these. It must not be. In other words, the present invention can be carried out in various forms without departing from the spirit or the main features thereof.
  • the present invention provides an image for executing a demosic process for generating a color image by interpolating luminance information of a plurality of colors into all pixels from a color mosaic image obtained using a single-plate color image sensor. It is suitable for use in a processing apparatus.

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Color Television Image Signal Generators (AREA)
  • Image Processing (AREA)

Abstract

 カラー画像の画素位置から、変形処理が施された場合のカラー画像の画素位置に対応する色モザイク画像上の対応するサンプリング座標を算出する座標変換部142と、色モザイク画像を分解した複数の色プレーン毎にサンプリング座標における画素値を補間生成するサンプリング部143と、各色プレーンの補間値を合成することによりカラー画像を生成する色生成部144とを備え、変形処理の施されたカラー画像の各画素値を、サンプリング座標の画素値として色モザイク画像から補間演算により求めることにより、色モザイク画像からカラー画像を生成する色補間処理と当該カラー画像の変形処理とを一度の補間演算によって実現することができるようにする。

Description

明 細 書 理装置、 画像処理方法、 プロ グラムおよび撮像装置 技術分野
本発明は画像処理装置 、 画像処理方法、 プロ グラムおよび撮像装置に 関 し、 特に、 単板式の力 ラ一撮像素子を用レ、て得られた色モザイ ク画像 から、 全画素に複数色の輝度情報を補間 してカ ラー画像を生成するデモ ザイ ク 処理を実行する画像処理装置、 画像処理方法 、 プロ グラムおよび 撮像装置に関する。 背景技術
近年、 消費者へのデジタルカ メ ラの普及が着実に進んでいる。 その理 由の 1 つに、 デジタノレ力 メ ラのコス トが下が り 続け 、 価格が多く の消費 者の手に十分届 く 範囲になつてきている こ と が挙げられる 。 コ ス ト を抑 えるために、 多く のデジタルカ メ ラでは、 レ、わゆる単板式の撮像素子を 採用 している。 単板式のテジタノレカ メ ラでは 、 力 ラ —画像内の各画素の 色情報を取り 込むために、 1 つの撮像素子のみを使用 している。
単板式の撮像素子では、 各画素が単色の色情報 しか持たない 。 と こ ろ が 、 力 ラ ―画像は、 3 つの別々 の単色画像を組み合わせる こ と によ り 表 される。 すなわち、 力 ラ一画像を表示するためには、 赤色 ( R ) 、 緑色
( G ) 、 青色 ( B ) の値すベてが各画素で必要である 。 こ のため単板式 のデジタル力 メ ラは 、 各画素が R , G , B成分の う ちの何れかのみを有 する色モザィ ク画像をも と に して 、 いわゆるデモザィ ク処理 (色補間処 理と も言 0 ) を行つてレゝる (例えば、 特許文献 1, 2 参照) 。 デモザィ ク処理は 、 色モザィ ク画像の各画素の単色情報に対してその周辺画素か ら集めた他の足 り ない色の輝度情報を用いて捕間演算を行う こ とによ り
、 各画素がそれぞれ R , G , B成分の全てを有する力ラー画像を生成す る処理である。
さ らに、 以上のよ う にして色モザィ ク画像から生成した力ラ一画像に 対して変形処理、 補正処理等の加工を施すこ とがある 。 例えば、 変形処 理の例と して、 デジタルズ一ム処理によ り 力ラ一画像を拡大または縮小 したり 、 手ブレ補正処理によ り 力ラ一画像を回転したりするこ とがある
。 また、 デジタルカメ ラ用のズームレンズでは、 その構成上 、 広角端で 樽型の歪曲収差を持つこ とが多い。 樽型の歪曲収差と糸卷型の歪曲収差 とが混在する、 いわゆる陣笠タイブの歪曲収差もある 。 この歪曲収差を 補正するために力ラー面像を変形処理するこ と もある 。 これらの画像変 形処理は、 変形後の補間位置に当たる各画素の情報を 、 その補間位置の 周辺にある変形前のカラ一画像が持つ複数画素の情報から補間するこ と によって行ってレゝる。
また、 補正処理の例と して、 色収差の補正を行う こ とがある。 デジタ ルカメ ラに用いられる結像光学系のレンズは 、 撮像光の波長によって屈 折率が異なる。 このため 、 R G Bの各色によつて像の倍率が異なってく る。 これによ り 、 撮像素子に結像する像の大き さが色成分毎に異なり 、 図 2 0に示すよ う に、 色成分毎に撮像素子上への結像ズレが発生するこ とが知られている。 これをレンズの倍率色収差 (横色収差) という。 倍 率色収差がある と、 白色点光源を撮影したと き、 特に画面の周辺部が虹 色に色づいて放射方向に伸びるよ う に見えてしま う。 また、 像ズレに 伴って、 画像のエッジ部分が色ズレと して現れ、 画質の品位を損なう と いう問題も生じる。
こ の よ うな色ズ レの発生を抑制するために、 撮像によ り得られたカラ —画像信号に基づいて、 カラー画像内の基準位置からの距離に応じた色 収差量を検出し、 検出した色収差量に基づいて、 カラ一画像信号に色収 差の補正を加えるよ う にした撮像装置が提供されている。 また、 カラー 画像内における有効なエッジを検出 して、 エッジ位置に基づいて、 基準 位置からの距離に応じた色収差量を検出する撮像装置もある。
さ らに、 光学レ ンズ系からの入射像を撮像して R , G, Bの信号を出 力し、 当該 R , G , Bの信号に対して光学レ ンズ系に固有の色収差に応 じた焦点距離補正値を用いてそれぞれ色収差補正 (Gを基準と して、 R , Bに対して画像の拡大または縮小を行う補正) を行った後、 色収差補 正ざれた R , G , B出力を合成して色収差補正された画像信号を出力す るよ う にしたデジタルカメ ラも提案されている (例えば、 特許文献 3参 照) 。
特許文献 1 : 特表 2 0 0 4— 5 3 4 4 2 9号公報
特許文献 2 : 特開 2 0 0 0— 2 7 0 2 9 4号公報
特許文献 3 : 特開平 6 — 1 1 3 3 0 9号公報 発明の開示
しかしながら、 上記従来の技術では、 色モザイ ク画像からカラー画像 を生成する際に補間処理を行う と と もに、 その後でカラー画像を変形する 際にも別の補間処理を行っている。 すなわち、 画素補間を二里に行つて レゝる 。 このため、 処理負荷が大き く なる とレ、う問題があつた。 また 、 補 間によ り生成した画像に対して更に別の補間を施すので 、 生成される画 像の画質劣化が大きく なるという問題もあつた。
また、 上記従来の技術では、 カラー画像內におけるェクジ位置に基づ いて色収差量を検出する際に、 R, G , Bの相関が正し < 取れず 、 ェッ ジ位置を検出するこ とができない場合がある。 すなわち 、 倍率色収差の 乗ハた色モザイ ク画像においては、 図 2 1 のよ う に 、 R G , Bの色成 分によつてサンプリ ング位置に対する輝度値が異なるため、 各色成分の ェッジ位置が一致しておらず、 高い解像度が得られない。 そのため 、 本 来のェッジ位置がエッジと して認識できず、 エッジ位置を利用 して モ ザィ クを正しく補正できない場合があるという問題があった。
なお 、 特許文献 3に記載されたデジタルカメ ラは単板式ではなく 、 い わゆる 3板式のデジタルカメ ラである。 3板式のデジタルカメ ラは 、 R
, G Bのそれぞれ用に 3つの撮像素子を備え、 それぞれの撮像素子か ら出力される R , G , Bの信号を合成してカラー画像を得るよ う に成さ れている。 3板式の場合、 R , G , Bの画素数は何れも出力画像の画 数と一致しているため、 比較的単純な画像合成でカラー画像を得るこ と がでさる。 すなわち、 出力画像の画素数に比べて R , G , Bの画素数が 何れち少ない単板式のデジタルカメ ラで行っているよ うなデモザィ ク処 理 (色補間処理) は不要である。
これに対して、 単板式のデジタルカメ ラの場合は、 色モザイ ク画像か らカラ一画像を生成する際に補間処理を行う と と もに、 倍率色収差を補 正する際にも別の補間処理を行う必要がある。 すなわち、 画素捕間を二 重に行わなければならない。 このため、 処理負荷が大き く なる という問 題があつた。 また、 補間によ り生成した画像に対して更に別の補間を施 すので 、 生成される画像の画質劣化が大き く なるという問題もあった。 本発明は、 このよ うな問題を解決するために成されたものであり 、 色 モザィ ク画像から変形処理の施された力ラ一画像を生成する際に、 よ り 少なレ、処理負荷で、 画質劣化の少ないカラ一画像を生成できるよ う にす ることを目的とする。
上記した課題を解決するために、 本発明では 、 色モザイ ク画像を同一 の色光の画素値だけを含む複数の色プレーンに分解する色プレーン分解 部と、 力ラー画像の画素位置から色モザィ ク画像上の対応するサンプリ ング座標を算出する座標変換部と 、 複数の色プレーン毎にサンプ y ング 座標における画素値を補間生成するサンプ リ ング部と 、 各色プ レ一ンの 補間値を合成する こ と によ り カ ラ一画像を生成する色生成部と を備え、 座標変換部において、 カ ラ一画像の画素位置から、 画像変形が施された 場合のカ ラー画像の画素位置に対応する色モザィ ク画像上のサンプリ ン グ座標を算出するよ う に している。
こ の よ う に構成 した本発明に れば、 色モザィ ク画像からカ ラー画像 を生成する刖段階と して 、 画像変形が施された場合の力 ラ一画像の画素 位置に対応する色モザィ ク画像上のサンプリ ング座標が算出される。 そ の上で、 当該サンプ ング座標の補間値を求める色補間演算がサンプ リ ング部によ り 行われる こ と によ り 、 画像変形の施された力 ラ一画像の各 画素値が、 サンプ リ ング座標の画素値と して色モザイ ク画像から求め ら れる こ と となる。
これによ り 、 色モザィ ク画像からカ ラ ―画像を生成する色補間処理 ( テモザィ ク処理) と 、 当該カ ラ一画像の変形処理と を一度の補間演算に つて実現する こ と がでさ る。 のため 、 色モザイ ク画像から変形処理 の施された力 ラー画像を生成する際の処理負荷を軽減する こ と ができ る と と もに、 従来のよ ラ に補間処理を二重に行 う こ と によ る画質の劣化も 抑制する とができ る <
また、 本発明の他の態様では、 色モザイ ク画像を同一の色光の画素値 だけを含む複数の色プレーンに分解する色プレ一ン分解部と 、 カ ラ一画 像の画素位置から色モザィ ク画像上の対応するサンプリ ング座標を算出 する座標変換部と 、 複数の色プレーン毎にサンプリ ング座標における画 素値を補間生成するサンプ リ ング部と 、 各色プレー ンの補間値を合成す る こ と によ り カ ラ一画像を生成する色生成部と を備え、 座標変換部にお いて、 複数の色プレ一ン毎に、 色プレーンによつて値の異なる係数を用 いて、 力ラ― 像の画素位置から、 色プレ—ン毎に異なるサンプジ ング 座標を算出するよ うにしてレ、る。
こ のよ う に構成した本発明によれば、 色モザィ ク画像から力ラ ―画像 ブ
を生成す 段階と して、 倍率色収差によ り 生じる各色成分の 像位置 ズレ【こ応じた色プレーン毎の画素位置が 、 カラー画像の画素位置に対応 する色モザィ ク画像上のサンプリ ング座標と して算出 ォしる その上で
、 当該サンプ ング座標の補間値を求める色補間演算がサンプ V ング部 によ り行われるこ とによ り 、 倍率色収差の補正されたカラー画像の各画 素値が、 サンプジ ング座標の画素値と して色モザィ ク画像から求められ ること となる。
これによ り 、 色モザィ ク画像から力ラ一画像を生成する色補間処理 ( デモザィ ク処理 ) と、 結像光学系の倍率色収差を補正する処理とを一度 の補間演算によつて実現するこ とがでさる。 このため、 色モザイ ク画像 から倍率色収差が補正されたカラー画像を生成する際の処理負荷を軽減 するこ とができる と と もに 、 従来のよ ラ に補間処理を二 に行う こ とに よる画質の劣化も抑制するこ とができる。
また、 本発明の他の態様では、 座標変換部は、 色プレ一ンによつて値 の異なる色収差係数に加えて、 カラー画像に対する画像変形を表す画像 変形係数、 撮像装置のブレを補正するためのブレ補正係数 、 撮像 子に 被写体の撮像光を導く ための結像光学系に設定されている焦点距離お び被写体距離に応じて定まる歪曲収差係数の少なく と も 1 つを用いて 、 力ラー画像の画素位置から 、 色プレーン毎に異なるサンプ ング座標を 算出するよ うにしてレ、る。
このよ う に構成した本発明によれば 、 色モザィ ク画像からカラ一画像 を生成する前段階と して、 倍率色収差によ り 生じる色ズレの補正に加え て、 画像変形、 ブ レ補正、 レンズステ ト の歪曲修正などに応じた色プ レーンの画素値が、 カラー画像の画素位置に対応する色モザィ ク画像上 のサンプリ ング座標と して算出され、 当該色ズレの補正に加えて画像変 形などの処理の施されたカラ一画像の各画素値が、 サンプリ ング座標の 画素値と して色モザイク画像から求められるこ と となる。
これによ り 、 色モザイ ク画像からカラ一画像を生成する色補間処理と
、 結像光学系の倍率色収差を補正する処理と、 さ らに力ラ一画像の画像 変形、 ブレ補正、 歪曲修正などの処理とを一度の補間演算によって実現 するこ とができる。 このため、 倍率色収差が補正され、 かつ 、 画像変形 などが施されたカラー画像を色モザィ ク画像から生成する際の処理負荷 を軽減するこ とができる と と もに、 補間処理を何回も行う こ とによる画 質の劣化も抑制することができる。
また、 本発明の他の態様では、 所定の画像を撮影するこ とによって BIJ 記色収差係数を算出する収差係数算出部を備え 、 この色収 係数を用い て、 カラー画像の画素位置から、 色プレーン毎に異なるサンプリ ング座 標を算出するよ うにしてレ、る。
このよ う に構成した本発明によれば、 所定の画像を撮影するこ とによ つて色収差係数を算出する収差係数算出部を備えているので 、 この算出 された色収差係数を用いて、 色モザイク画像から力ラ ―画像を生成する 色捕間処理 (デモザイ ク処理) と、 結像光学系の倍率色収差を補正する 処理とを一度の補間演算によって実現することができる。
また、 本発明の他の態様では、 外部に接続された外部機器に撮影画像 を出力する画像出力部と、 前記外部機器から前記色収差係数を入力する 係数入力手段とを備え、 この入力された色収差係数を用いて、 カラー画 像の画素位置から、 色プレーン毎に異なるサンプリ ング座標を算出する よ うにしている。
このよ う に構成した本発明によれば、 外部に接続された外部機器に撮 影画像を出力する画像出力部と、 前記外部機器から、 前記色収差係数を 入力する係数入力手段と を備えているのでこ の入力された色収差係数を 用いて、 色モザイ ク画像からカラ一画像を生成する色補間処理 (デモザ イ ク処理) と、 結像光学系の倍率色収差を補正する処理とを一度の補間 演算によって実現するこ とができる。 図面の簡単な説明
図 1 は、 本発明に係る画像処理装置を実施した第 1 の実施形態による カラー撮像装置の構成例を示す図である。
図 2は、 第 1 〜第 5 の実施形態によるカラーフ ィ ルタ ア レイ の原色べ ィャ配列を示す図である。
図 3は、 第 1 の実施形態によるデモザィ ク部の機能構成例を示すプロ ック図である。
図 4は、 第 1 の実施形態のデモザィ ク部によ り実行される画像処理の 動作例を示すフローチヤ一トである。
図 5 は、 第 1 の実施形態のデモザイク部によ り実行される画像処理の 一例を具体的に説明するためのイメージ図である。
図 6は、 第 1 の実施形態のバイ リ ニア補間を説明するためのイ メージ 図である。
図 7 は、 第 1 の実施形態のデモザィ ク部によ り 実行される画像処理の 他の例を具体的に説明するためのイメージ図である。
図 8は、 本発明に係る画像処理装置を実施した第 2の実施形態および 第 3の実施形態によるカラー撮像装置の構成例を示す図である。
図 9は、 第 2の実施形態および第 3の実施形態によるデモザイ ク部の 機能構成例を示すブロ ック図である。
図 1 0は、 第 2〜第 5の実施形態のデモザイ ク部によ り 実行される画 像処理の動作例を示すフロ一チヤ一 トである。
図 1 1 は、 第 2の実施形態および第 3の実施形態のデモザイ ク部によ り実行される画像処理の一例を具体的に説明するためのイ メージ図であ る。
図 1 2は、 第 2〜第 5の実施形態のデモザイ ク部によ り実行される画 像処理の他の例を具体的に説明するためのイメージ図である。 図 1 3は、 本発明に係る画像処理装置を実施した第 4 の実施形態によ る力ラ ―撮像装置の構成例を示す図である。
図 1 4は、 第 4 の実施形態による収差係数テープルの一例を示す図で ある
図 1 5は、 第 4の実施形態によるデモザイ ク部の機能構成例を示すブ 口 ック図である。
図 1 6 は、 第 5 の実施形態によるカラー撮像装置の構成例を示す図で ある 図 1 7は、 第 5 の実施形態における色収差量検出の説明図であって、 図 1 7 ( a ) がチャー トの形態を表す図、 図 1 7 ( b ) が撮像素子に対 するチャ ― トの配置を表す図である。
図 1 8 は、 第 5の実施形態における色収差量検出の説明図であって、 図 1 8 ( a ) ( b ) が交点検出の際の説明図、 図 1 8 ( c ) が交点毎に エッジを検出する際の説明図、 図 1 8 ( d ) が交点に対してサンプリ ン グ画素列を設定する際の説明図である。
図 1 9は、 第 5の実施形態におけるエッジ検出の説明図である。
図 2 0は、 レンズの倍率色収差を説明するための図である。
図 2 1 は、 倍率色収差による色ズレを説明するための図である。
図 2 2は、 偽色の抑制を考慮して R G Bのカラー情報を Y U Vのカラ
—情報に変換する際に用いる式を示す図である
図 2 3は、 X 座標と u 座標との関係式を示す図である。 図 2 4は、 ブレ補正を行う際に用いる式を示す図である。
図 2 5は、 色プ レーン間の座標変換を行う際に用いる式を示す図であ る。
図 2 6 は、 各色プ レーン上の X y座標に対応する色モザイ ク画像上の 画素座標の関係式を示す図である。
図 2 7は、 色収差行列を求める際に用いる式を示す図である。 発明を実施するための最良の形態
(第 1 の実施形態)
以下、 本発明の一実施形態を図面に基づいて説明する。 図 1 は、 本発 明に係る画像処理装置を実施した第 1 の実施形態によるカラー撮像装置 1 0 0の構成例を示す図である。 本実施形態のカラ一撮像装置 1 0 0は 、 結像光学系 1 1 0、 単板カラーイメージセ ンサ 1 2 0 、 A D変換部 1 3 0、 デモザイ ク部 1 4 0、 視覚補正部 1 5 0、 圧縮部 1 6 0および記 録部 1 7 0 を備えて構成されている。 この う ち、 デモザイ ク部 1 4 0が 本発明の画像処理装置に相当する。
単板カラーイ メージセンサ 1 2 0は、 結像光学系 1 1 0から出力され た撮像光を所定の色成分に分解するカラ一フ ィ ルタ ア レイ 1 2 1 と、 力 ラ一フ ィ ルタ ア レイ 1 2 1 を通過した撮像光を光電変換して画素信号を 生成する撮像素子 1 2 2 と を備えている。 撮像素子 1 2 2は、 例えば C C D ( Charge Coupled Device ) 、 あるレヽは C M O S ■ ( Complementary Metal Oxide Semiconductor) センサ等で構成される。
結像光学系 1 1 0は、 被写体の撮像光を単板カラーイメージセンサ 1 2 0に導く役割を果たすものである。 例えば、 結像光学系 1 1 0は、 光 学的ロ ーパスフィルタを含め、 撮影レンズや、 赤外線除去フィルタ等か ら構成されている。 なお、 赤外線除去フ ィ ルタは、 単板カラ一イメージ センサ 1 2 0 に入射する赤外線を遮断するためのものであ り 、 光学的口
―パスフィルタの前方に配置され、 1枚のガラスブロ クク と して構成さ れている
単板力 ラ一ィ メ ージセンサ 1 2 0 のカ ラーフ イ ノレタァ レィ 1 2 1 は、 撮像素子 1 2 2 を構成する各画素の受光面上に所定パタ ―ンで規則的に 配置され 像光を所定の色成分にフィ ルタ ングする役割を果たすも のである 本実施形態では、 色成分と して R G, B の 3色を用いた原 色べィャ配列のカラーフィルタアレイ 1 2 1 を使用してレ、る。
原色べィャ配列は、 図 2 に示すよ う に、 G色フィルタを巿松模 で配 置する と と もに、 R色フィルタおよび B色フ ィ ルタを各行に交互に配置 したものである。 なお 、 G色フィルタのう ち、 水平方向の R色フィルタ どう しで挟まれたものを G r 色フイ ノレタ、 水平方向の B色フィノレタ どう しで挟まれたものを G b色フィルタ と称する。
fe像素子 1 2 2は、 受光した撮像光を電気的な画素情報に光電変換し て電荷量と して貯え、 これを 気信号と して A D変換部 1 3 0に出力す る役割を果たすものである。 撮像素子 1 2 2は、 所定パターンで配列さ れた複数の画素 (フ才 トダイォ — ド) を有しており 、 この各画素の受光 面上にカラーフイ ノレタ ア レイ 1 2 1 が上述の原色ベィャ配列で規則的に 配置されている
以上の構成にぉレ、て 、 ¾ ®ίされた被写体は 、 Ψα像光学系 1 1 0 を経て 単板力ラーイ メ —ジセンサ 1 2 0 の撮像素子 1 2 2上に結像する。 この 際、 結像光学系 1 1 0が持つ各種の収差によ り 、 像した被写体像は劣 化したものとなる。 例えば、 被写体上では直線であるものが歪曲収差に よ り 曲線となつた画像となる。
単板カラ一ィメージセンサ 1 2 0は、 撮像素子 1 2 2上に結像した被 写体像を色モザィ ク画像と してアナ口 グ電気信号に変換する。 すなわち 、 カラーフィルタア レイ 1 2 1 は、 図 2 に示すよ う な原色べィャ配列で あり 、 撮像素子 1 2 2の各画素に対して R , G , Bの各色フィルタが並 ベられている。 このため、 被写体の撮像光は、 画素毎に対応した色成分 のみが透過した色モザイ ク画像状態の撮像光となって撮像素子 1 2 2 に 達する。 撮像素子 1 2 2は、 その到達光を光電変換し、 色モザイ ク画像 の電気信号と して AD変換部 1 3 0に出力する。
A D変換部 1 3 0は、 撮像素子 1 2 2 によ り光電変換された色モザィ ク画像のアナ口 グ信号を、 デジタル信号処理が行えるよ う にデジタル信 号に変換する。 なお、 A D変換部 1 3 0で AZD変換された直後の色モ ザイ ク画像は R AWデータ と も呼ばれる。 デモザイ ク部 1 4 0は、 色モ ザイ ク画像からカラ一画像への変換を行う。 本実施形態では、 このと き 同時に画像の変形処理を行う こ とによ り 、 上述した結像光学系 1 1 0の 歪曲収差による画像の画質劣化を補正する。 このデモザィ ク部 1 4 0 に よる画像処理の方法は後に詳しく説明する。
視覚補正部 1 5 0は、 デモザイ ク部 1 4 0 によ り生成されたカラー画 像に対して、 主と して面像の見栄えを良く するための処理を行う。 例え ば、 視覚補正部 1 5 0は、 トーンカーブ (ガンマ) 補正、 彩度強調、 ェ ッジ強調といった画像補正処理を行う。 圧縮部 1 6 0は、 視覚補正部 1 5 0で補正されたカラ一画像を J P E G ( Joint Photographic Experts Group) 等の方法で圧縮し、 記録時のサイズを小さ くする。 記録部 1 7 0 は、 圧縮されたデジタル画像信号をフラ ッシュメ モ リ等の記録媒体 (図 示せず) に記録する。
なお、 デモザイ ク部 1 4 0から記録部 1 7 0までの各構成は、 それぞ れを別体のデバイスと して構成しても良いし、 単一のマイ ク ロプロセッ サで構成しても良い。 後者の場合は、 単一のマイ ク ロプロセッサが、 デ モザイ ク部 1 4 0から記録部 1 7 0までの各構成に係る処理を実行する 図 3 は、 デモザイ ク部 1 4 0 の機能構成例を示すプロ ッ ク図である。 図 4 は、 デモザイ ク部 1 4 0 によ り 実行される画像処理の動作例を示す フ ロ ーチャー ト である。 図 5 は、 デモザイ ク部 1 4 0 によ り 実行される 画像処理の内容を具体的に説明するためのィ メ ージ図である。 図 3 に示 すよ う に、 デモザイ ク部 1 4 0 は、 その機能構成と して、 色プ レーン分 解部 1 4 1 と 、 座標変換部 1 4 2 と 、 サンプリ ング部 1 4 3 と 、 色生成 部 1 4 4 と を備えている。
色プ レーン分解部 1 4 1 は、 A D変換部 1 3 0 よ り 出力 される色モザ ィ ク画像を、 同一の色光の画素値だけを含む複数の色プレーンに分解す る (図 4 のステ ップ S 1 ) 。 本実施形態では、 図 5 に示すよ う に、 R成 分の画素のみを取 り 出 した Rプレーン、 G成分の画素のみを取 り 出 した Gプレーン、 B成分の画素のみを取 り 出 した Bプレーンの 3 つの色プレ —ンに分解する。 分解された各色プ レーンは、 サンプ リ ング部 1 4 3 の 処理に利用される。
座標変換部 1 4 2 は、 カ ラー画像に対する画像変形を表す係数を用い て、 色モザイ ク画像から生成されるカ ラー画像の画素位置から、 変形処 理が施された場合のカ ラー画像 (以下、 変形カ ラ一画像と称する) の画 素位置に対応する色モザィ ク画像上のサンプリ ング座標を算出する (図 4 のステ ップ S 2 ) 。
例えば、 結像光学系 1 1 0 の収差によ り 歪曲 しているカ ラー画像を補 正するために画像変形を行 う場合には、 色モザィ ク画像から生成された 歪曲 しているカ ラー画像に対して非線形の座標変換を行 う。 こ の座標変 換の方法は既知であ り 、 変形カ ラ一画像上の画素位置およびそれが色モ ザイ ク画像上でどの位置 (サンプリ ング座標) に該当するのかについて は、 計算によ り 求められる。 以下に、 サンプリ ング座標の算出手順を詳しく説明する。 まず x y座 標系と して、 原点を画像中心、 最大像高 (原点からの最大距離) を 1 と し、 原点から画面右方向に正の X座標をと り , 原点から画面下方向に正 の y座標をと るものとする。 この場合、 6 4 0 X 4 8 0の正方画素力、ら なるカラー画像に対して、 u v座標系の画素座標 ( u d, V d) を図 2のよ うに画面左上から右方向に ( 0, 0 ) , ( 1 , 0 ) , ( 2 , 0 ) · · · 、 次の行を ( 1 , 0 ) , ( 1 , 1 ) , ( 2, 1 ) · · ' と割り振れば、 u v 座標系の画素座標 ( 3 1 9 . 5 , 2 3 9 . 5 ) が x y座標系の原点とな る。 また、 画素座標系の対角長の半分の長さ 4 0 0 = ( 6 4 0 2 + 4 8 0 2 ) 1/2Z 2が X y座標系の最大像高に対応するため、 画素座標 ( u d, V d ) に対応する X y座標 ( X d, y d) は、
x d= ( u d- 3 1 9 . 5 ) / 4 0 0
y d= ( v d- 2 3 9 . 5 ) / 4 0 0
の関係式で表される。
この x y座標 ( x d, y d) に対して、 結像光学系 1 1 0の歪曲収差の補 正を考慮した座標変換を以下の式のよ うに行う。
= d V k , r 2+ k 2 r 4)
y = y d ( k , r 2+ k 2 r 4)
(ただし、 r 2= X d 2 + y d 2)
なお、 { k ,, k 2} は結像光学系 1 1 0の歪曲収差を示す係数であり 、 k , は 3次収差係数、 k 2 は 5次収差係数を示している。 上述のよ うに、 歪 曲収差を含むカラー画像に対して非線形の座標変換を行う方法は既知で あり、 歪曲収差係数 { k ,, k 2} はシ ミ ュ レーシ ョ ン等によ り求めること が可能である。 この歪曲収差係数 { k ,, k 2} が、 本発明の画像変形係数 に相当する。
一方、 色モザイ ク画像は 1 6 0 0 X I 2 0 0 の正方画素から成るもの と し、 u v座標系の画素座標 ( u s, v s) が上述のカ ラ一画像と同様に割 り 振られている と する と 、 u v座標系の画素座標 ( 7 9 9 . 5, 5 9 9 . 5 ) が x y座標系の原点と な り 、 画素座標系の対角長の半分の長さ 1 0 0 0 = ( 1 6 0 02+ 1 2 0 02) 1/2 2 が X y座標系の最大像高に対応 する。 このため、 変形カ ラー画像上の x y座標 ( X , y ) に対応する色 モザイ ク画像上の画素座標 ( u s, V s) は、
u s= 1 0 0 0 * x + 7 9 9 . 5
v s= 1 0 0 0 * y + 5 9 9 . 5
と なる。
上述の計算の結果、 画素座標 ( u s, V J は整数値になる とは限らず、 一般には非整数と なる。 この色モザイ ク画像上の画素座標 ( U s, V s) が サンプリ ング座標である。 図 5 では、 こ のサンプ リ ング座標の 1 っを符 号 2 0 0 で示 している。 上述のよ う に、 色モザイ ク画像は 3 つの色プレ ーンに分解されている。 図 5 では、 それぞれの色プレーン上にサンプリ ング座標 2 0 0 を図示 している。 上述のよ う に、 サンプリ ング座標の値 は非整数と なるので、 画素中心からずれた位置にサンプリ ング座標 2 0 0 が存在している。
サンプリ ング部 1 4 3 は、 色プレーン分解部 1 4 1 によ り 分解された 複数の色プレーン毎に、 座標変換部 1 4 2 によ り 算出 されたサンプ リ ン グ座標 2 0 0 における画素値 (サンプリ ング値) を、 色プレーン内に含 まれる同一の色光の画素値から補間生成する (図 4 のステ ップ S 3 ) 。 すなわち、 サンプリ ング部 1 4 3 は、 Rプレーン、 Gプレーンおよび B プレーンのそれぞれから、 サンプリ ング座標 2 0 0 の画素値を補間演算 によ り算出 して出力する。
上述のよ う に、 サンプリ ング座標 2 0 0 の値 ( u s, V s) は必ずしも整 数値と は限らないため、 当該サ ンプリ ング座標 2 0 0 を囲む 4 つの有値 画素 (各色プレーンが元々持っている同一色光の画素値) から線形補間 を行 う。 この線形補間は、 好ま しく はバイ リ ニア補間によって行う。 図 5 に示すよ う に、 Rプ レーンおよび Bプ レーンは縦横の格子点状に 有値画素を持っため、 サンプ リ ング座標 2 0 0を囲む 4つの有値画素は 、 当該サンプリ ング座標 2 0 0 を囲む一辺の長さが 2 の正方形の各頂点 に位置する。 例えば、 サンプ リ ング座標 2 0 0力 ( u s, V s) = ( 1 0 0 . 8, 1 0 1 . 4 ) であれば、 Rプレーンにおいてこれを囲む 4つの画 素 ( u d, V d) = ( 1 0 0 , 1 0 0 ) , ( 1 0 0, 1 0 2 ) , ( 1 0 2, 1 0 0 ) , ( 1 0 2 , 1 0 2 ) が Rプレーンの有値画素と なる。
当該有値画素の各画素値を R ( 1 0 0, 1 0 0 ) , R ( 1 0 0, 1 0 2 ) , R ( 1 0 2, 1 0 0 ) , R ( 1 0 2 , 1 0 2 ) で表すとする と 、 図 6 のよ う なバイ リ ニア補間によって生成される Rプレーン上でのサンプリ ング座標 2 0 0の補間画素値 R ( 1 0 0. 8 , 1 0 1 . 4 ) は、 次の式 で表される。
R ( 1 0 0. 8 , 1 0 1 . 4 ) = 0. 6 * 0. 3 * R ( 1 0 0 , 1
0 0 ) + 0. 6 * 0. 7 * R ( 1 0 0 , 1 0 2 ) + 0. 4 * 0. 3 * R ( 1 0 2, 1 0 0 ) + 0. 4 * 0. 7 * R ( 1 0 2 , 1 0 2 )
また、 Bプレーンにおいてサンプリ ング座標 2 0 0の位置 ( u s, V s) = ( 1 0 0. 8 , 1 0 1 . 4 ) を囲む 4つの画素 ( u d, V d) = ( 9 9 , 1 0 1 ) , ( 9 9 , 1 0 3 ) , ( 1 0 1 , 1 0 1 ) , ( 1 0 1 , 1 0 3 ) が Βプ レーンの有値画素と なる。 当該有値画素の各画素値を Β ( 9 9, 1
0 1 ) , Β ( 9 9 , 1 0 3 ) , Β ( 1 0 1 , 1 0 1 ) , Β ( 1 0 1 , 1 0 3 ) で表すとする と 、 Βプレーン上でのサンプ リ ング座標 2 0 0の捕間画 素値 Β ( 1 0 0. 8 , 1 0 1 . 4 ) は、 次の式で表される。
Β ( 1 0 0. 8, 1 0 1 . 4 ) = 0. 1 * 0. 8 * Β ( 9 9, 1 0 1 ) + 0. 1 * 0. 2 * Β ( 9 9 , 1 0 3 ) + 0. 9 * 0. 8 * Β ( 1 0 1 , 1 0 1 ) + 0 . 9 * 0. 2 * B ( 1 0 1 , 1 0 3 )
一方、 Gプ レーンは市松状に有値画素を持っため、 サンプ リ ング座標 2 0 0 を囲む 4 つの有値画素は、 当該サンプ リ ング座標 2 0 0 を囲む一 辺の長さが 2の 4 5度傾いた正方形の各頂点に位置する。 この場合、 Gプレーンにおいてサンプリ ング座標 2 0 0 の位置 ( u s, V s) = ( 1 0 0 . 8, 1 0 1 . 4 ) を囲む 4 つの画素 ( u d, V d) = ( 1 0 0, 1 0 1 ) , ( 1 0 1 , 1 0 0 ) , ( 1 0 1 , 1 0 2 ) , ( 1 0 2 , 1 0 1 ) 力 プ レーンの有値面素となる。
当該有値画素の各画素値を G ( 1 0 0 , 1 0 1 ) , G ( 1 0 1 , 1 0 0 ) , G ( 1 0 1 , 1 0 2 ) , G ( 1 0 2 , 1 0 1 ) で表すとする と、 Gプ レーン上でのサンプ リ ング座標 2 0 0 の補間画素値 G ( 1 0 0 . 8 , 1 0 1 . 4 ) は、 次の式で表される。
G ( 1 0 0 . 8 , 1 0 1 . 4 ) = 0 . 7 * 0 . 3 * G ( 1 0 0 , 1 0 1 ) + 0 . 3 * 0 . 3 * G ( 1 0 1 , 1 0 0 ) + 0 . 7 * 0 . 7 * G ( 1 0 1 , 1 0 2 ) + 0 . 3 * 0 . 7 * G ( 1 0 2 , 1 0 1 )
色生成部 1 4 4は、 サンプリ ング部 1 4 3 によ り補間生成された各色 プ レーンの補間画素値を合成するこ とによ り 、 各面素値がそれぞれ複数 色の輝度情報を有したカラ一画像を生成する (図 4のステップ S 4 ) 。 例えば、 ある 1つのサンプリ ング座標 ( u s, V s) に対してサンプリ ング 部 1 4 3によ り求められた各色プレーンの補間値をそれぞれ R ( u s) V s ) , G ( u s, v J , B ( u s, v s) と した場合、 色生成部 1 4 4は、 こ れら 3つの補間値を合成することによ り、 サンプリ ング座標 ( u s, V s) におけるカラ一情報を生成する。
さ らに、 色生成部 1 4 4は、 こ う して求めた R G Bのカラー情報を Y U Vのカラー情報に変換し (Yは輝度情報、 U , Vは色情報) 、 U , V の色情報に対して低周波フ ィ ルタを施す。 なお、 R G Bから Y U Vへの 変換処理および U , Vに対する低周波フィルタの処理は、 公知の処理を 適用するこ とが可能である。 色生成部 1 4 4 は、 以上の処理を変形カラ 一画像の全ての画素 (全てのサンプリ ング座標) について行い、 その結 果得られた変形カラー画像を視覚補正部 1 5 0に出力する。 視覚補正部 1 5 0以降の処理は上述した通りである。
以上詳し <説明したよ う に、 第 1 の実施形態では、 色モザィ ク画像か ら変形カラ一画像を生成する前段階と して、 出力すべき変形カラ一画像 の画素位置に対応する色モザィ ク画像上のサンプリ ング座標 2 0 0 を算 出する。 そして、 変形カラー画像の画素値を、 当該サンプリ ング座標 2
0 0における画素値と して色モザィ ク画像の画素値を用レ、た補間演算に よ り生成するよ うにしてレ、る。
これによ り 、 色モザィ ク画像から力ラー画像を生成する色補間処理と 当該カラー画像の変形処理とを一度の補間演算によって 現するこ とが できる。 このため、 色モザィ ク画像から変形カラ一画像を生成する際の 処理負荷を軽減するこ とができる と と もに、 従来のよ う に - 補間処理を 重に行う ことによる画質の劣化も抑制することができる。
また、 第 1 の実施形態では、 色モザィ ク画像を R , G , B の複数の色 プレーンに分解し 、 色プレーン毎にサンプリ ング座標 2 0 0の補間画素 値を求めてから、 各色プ レーンの補間画素値を合成して 1 画素に 3色の 輝度情報が含まれるカラ一情報を生成するよ うにしている。 のよ う—に するこ とによ り 、 色プレーン内に含まれる同一の色光の画素値から単純 な線形補間によ りサンプリ ング座標 2 0 0の捕間画素値を求めるこ とが でき、 処理負荷を軽減することができる。
なお、 上記第 1 の実施形態では、 色プレーン分解部 1 4 1 は 、 色モザ イ ク画像を R プレーン、 G プレーン、 B プ レーンの 3 つの色プ レ一ンに 分解する例について説明 したが、 これに限定されない 。 例えば 、 図 7 に 示すよ う に、 互いに角 を接する 2 種類の緑色画素 G r , G b を各々異な る色プレーンに分解する よ う に して も良い。 すなわち、 色プレーン分解 部 1 4 1 は、 A D変換部 1 3 0 よ り 出力 される色モザイ ク画像を、 R成 分の画素のみを取 り 出 した Rプレーン、 G r 成分の画素のみを取 り 出 し た G r プレーン、 G b 成分の画素のみを取 り 出 した G b プレー ン、 B成 分の画素のみを取り 出 した Bプレーンの 4 つの色プレーンに分解する。
この場合に色生成部 1 4 4 は、 G r プレーンおよび G b プレーンでそ れぞれサンプ リ ング部 1 4 3 によ り 生成されたサンプリ ング座標 2 0 0 の画素値 (サンプリ ング値) を加算 して、 1 画素内に R , G , Bの各色 成分の輝度情報が含まれるカ ラ一情報を生成する。 例えば、 R成分や B 成分はサンプ リ ング値をそのまま用い、 G成分は G r と G b と の平均値 を用いる。
また、 G成分を G r プレーンと G b プレーンと に分解してサンプリ ン グ座標 2 0 0 の画素値を算出する場合、 デモザイ ク部 1 4 0 は、 G r プ レーンと G b プレーンの各色プレーン上でそれぞれサンプリ ング部 1 4 3 によ り 生成された補間画素値の差分を算出 し、 当該補間画素値の差分 に基づいて偽色の有無を判定する偽色判定部を更に備える よ う に して も 良い。
べィャ配列の単板カ ラ一イ メ ージセンサ 1 2 0 では、 ナイ キス ト周波 数近傍の白黒の縞模様に対し、 赤や青の偽色が発生する と い う 問題があ る。 これに対し、 G r プレーン上で求めたサンプリ ング座標 2 0 0 の補 間画素値と G b プレーン上で求めたサンプリ ング座標 2 0 0 の補間画素 値と の差分を と る こ と で、 縞模様上での偽色の有無を検出する こ と がで き、 偽色がある場合にはこれを抑制する こ とができ る。
すなわち、 G r も G b も本来は同 じ G色フ ィルタなので、 G r プレー ンおよび G b プレーンで求めた補間画素値は双方と も同 じ値になるはず である。 しカゝし、 偽色が生じている と、 G r プレーンおよび G b プレ一 ンで求めた補間画素値に差分が生じる。 よって、 補間画素値の差分を見 るこ と によ り 、 偽色の発生の有無を検出するこ とができる。 偽色が存在 する と判定した場合には、 色生成部 1 4 4 は、 R G B の力ラー情報を Y U Vのカラー情報に変換する際に、 偽色の抑制を考慮した図 2 2 の式に よって変換処理を行う。
また、 上記第 1 の実施形態では、 カラ一画像に対する変形処理の例と して歪曲修正を挙げて説明 したが、 これに限定されない。 例えば 上記 歪曲修正のための面像変形に代えて、 またはこれに加えて、 力ラ 画像 の拡大、 縮小、 回転のう ち少なく と も 1つを含む処理を行う場合にち、 本実施形態のデモザイ ク部 1 4 0 を適用するこ とが可能である。 の場 合、 拡大、 縮小、 回転は例えばァフ ィ ン変換で表すこ とが可能であり 、 そのァフィ ン変換を表す係数を画像変形係数と して用いる。
(第 2 の実施形態)
以下、 本発明の第 2の実施形態を図面に基づいて説明する。 図 8は、 本発明に係る画像処理装置を実施した第 2 の実施形態によるカラー撮像 装置 1 0 0 ' の構成例を示す図である。 図 8 に示すよ う に、 第 2 の実施 形態のカラ一撮像装置 1 0 0 ' は、 結像光学系 1 1 0、 単板カラ一ィメ一 ジセンサ 1 2 0 、 A D変換部 1 3 0、 デモザイ ク部 1 4 0 ' 、 視覚補正部 1 5 0、 圧縮部 1 6 0および記録部 1 7 0を備えて構成されている。 この う ち、 デモザイク部 1 4 0 ' が本発明の画像処理装置に相当する。 なお、 図 8において、 図 1 に示した符号と同一の符号を付したものは同一の機能 を有するものであるので、 こ こでは重複する説明を省略する。
第 2 の実施形態において、 撮影された被写体は、 結像光学系 1 1 0 を 経て単板カラーイ メージセンサ 1 2 0 の撮像素子 1 2 2上に結像する。 こ の際、 結像光学系 1 1 0 の倍率色収差によ り 、 R, G, B の色成分毎 に撮像素子 1 2 2上での結像ズレ (色ズ レ) が生じた画像となる。 デモ ザイ ク部 1 4 0 ' は、 色モザイ ク画像からカラ一画像への変換を行う。 本実施形態では、 このと き同時に倍率色収差の補正処理を行う こ とによ り 、 上述した結像光学系 1 1 0 の倍率色収差による色ズレを補正する。 図 9 は、 デモザイク部 1 4 0 ' の機能構成例を示すブロ ック図である 。 図 1 0は、 デモザイク部 1 4 0 ' によ り実行される画像処理の動作例を 示すフ ローチャー トである。 図 1 1 は、 デモザイ ク部 1 4 0, によ り実行 される画像処理の内容を具体的に説明するためのイメージ図である。 図 9 に示すよ う に、 デモザイ ク部 1 4 0 ' は、 その機能構成と して、 色プレーン分解部 1 4 1 と、 座標変換部 1 4 2 , と、 サンプジ ング部 1
4 3 ' と 、 色生成部 1 4 4 とを備えている。 なお 、 図 9 におレ、て、 図 3 ίこ示した符号と同一の符号を付したものは同一の機能を有するものであ るので、 ここでは異なる部分を重点的に説明する。
色プ レ —ン分解部 1 4 1 は、 図 1 1 に示すよ う に、 A D変換部 1 3 0 よ り 出力される色モザイ ク画像を、 R成分の画素のみを取り 出した Rプ レーン、 G成分の画素のみを取り 出した Gプレ一ン、 B成分の画素のみ を取り 出した Β プ レーンの 3 つの色プレーンに分解する (図 1 0 のステ ップ S 1 1 ) ο
座標変換部 1 4 2 ' は、 色モザイ ク画像から生成される力ラ一画像の 画素位置から、 色収差補正が施された場合の力ラ一画像の画素位置に対 応する色モザィ ク画像上のサンプリ ング座標を算出する (図 1 0のステ ップ S 1 2 ) 。 本実施形態では特に、 座標変換部 1 4 2 ' は 、 色プレ一 ン分解部 1 4 1 によ り分解された複数の色プレ一ン毎に、 色プレ一ンに よって値の異なる色収差係数を用いて、 色プレ一ン毎に異なるサンプリ ング座標を算出する。 すなわち、 色モザイ ク画像には結像光学系 1 1 0 の倍率色収差の影響があるが、 カラ一画像にはその影響が現れないよ う にするため、 サンプリ ング座標は、 異なる色光に対応する色プ レーンに 対して異なる値を与えるよ うにする。
以下に、 第 2 の実施形態によるサンプリ ング座標の算出手順を詳しく 説明する。 なお、 こ こ で考慮する x y座標系は第 1 の実施形態と同様で あり、 画素座標 ( u d, V d) に対応する x y座標 ( x d, y d) は、
x d= ( u d- 3 1 9 . 5 ) / 4 0 0
Y d= ( v d- 2 3 9 . 5 ) / 4 0 0
の関係式で表される。
この x y座標 ( x d, y d) に対して、 結像光学系 1 1 0 の倍率色収差の 補正を考慮した座標変換を以下の式のよ う に行う。 こ こ では、 Rプレー ン、 Gプ レーン、 B プ レーンの 3 つの色プレーン毎に座標変換を別々に 行っている。
X R= X d * k R ¾ y R= y d * k R
x G= x d * 1 、 y c= y a * 1
X B= X d * k B , y B= Υ d * k B
なお、 { k R, k B} は結像光学系 1 1 0 の倍率色収差を示す係数であり 、 k R は G成分に対する R成分の倍率、 k B は G成分に対する B成分の倍 率を示している。 これらの色収差係数 { k R, k B} は光学シミ ュ レーショ ン等によ り求めるこ とが可能である。
一方、 色モザイ ク画像は 1 6 0 0 X 1 2 0 0 の正方画素から成るもの と し、 u v座標系の画素座標 ( u s, V s) が上述のカラ一画像と同様に割 り振られている とする と、 u v座標系の画素座標 ( 7 9 9 . 5 , 5 9 9 . 5 ) が x y座標系の原点となり 、 画素座標系の対角長の半分の長さ 1 0 0 0 = ( 1 6 0 02+ 1 2 0 02) >/2/ 2が X y座標系の最大像高に対応 する。 このため、 色収差補正されたカラー画像の各色プレーン上の x y に対応する色モザィ ク画像上 画素座標 ( u sK , V ( u sG sG
u SR= 1 0 0 0 * X f + 7 9 9 . 5
sR= 1 0 0 0 * y ,+ 5 9 9 . 5
U sG= 1 0 0 0 * X ;+ 7 9 9 . 5
sG= 1 0 0 0 * y ;+ 5 9 9 . 5
U sB= 1 0 0 0 * X B ,+ 7 9 9 . 5
sB= 1 0 0 0 * y B ;+ 5 9 9 • 5
となる。
上述の計算の結果、 画素座標 ( U sR, V sR) ( U sG, V sG) ( U sB ) V sB
) は何れも整数値になる と は限らず、 一般には非整数と なる。 こ の色モ ザイ ク画像上の画素座標 ( U sR, V sR) ( U sG ) V sG) ( U sB, V sB) が色 プレーン毎のサンプリ ング座標である。 図 1 1 では、 Rプレーン上のサ ンプ リ ング座標を符号 3 0 1 で示し、 Gプレーン上のサンプリ ング座標 を符号 3 0 2 で示 し、 B プレーン上のサンプ リ ング座標を符号 3 0 3 で 示している。 上述のよ う に、 サンプリ ング座標の値は非整数と なるので 、 画素中心からずれた位置に各色プ レーンのサンプリ ング座標 3 0 1 , 3 0 2 , 3 0 3 が存在してレ、る。
サンプ リ ング部 1 4 3 ' は、 色プレーン分解部 1 4 1 によ り 分解され た複数の色プレーン毎に、 座標変換部 1 4 2 ' によ り 算出された色プレ ーン毎のサンプリ ング座標 3 0 1 , 3 0 2 , 3 0 3 における画素値 (サ ンプ リ ング値) を、 色プレーン内に含まれる同一の色光の画素値から補 間生成する (図 1 0 のステ ップ S 1 3 ) 。 すなわち、 サンプ リ ング部 1 4 3 ' は、 Rプレーン、 Gプレーンおよび B プレーンのそれぞれ力ゝら、 サンプリ ング座標 3 0 1 , 3 0 2 , 3 0 3 の画素値を補間演算によ り 算 出 して出力する。 上述のよ う に、 サンプリ ング座標 3 0 1 , 3 0 2 , 3 0 3の値 ( u sR
V sR) ( U sG, V sC) ( U sB ) V sB) は必ずしも整数値と は限らないため、 当該サンプリ ング座標 3 0 1 , 3 0 2 , 3 0 3 を囲む 4つの有値画素 ( 各色プレーンが元々持っている同一色光の画素値) から補間を行 う。 こ の捕間は、 好ま しく はバイ リ ニア補間によって行う。
図 1 1 に示すよ う に、 Rプレーンおよび Bプレーンは縦横の格子点状 に有値画素を持っため、 サンプリ ング座標 3 0 1 , 3 0 3 を囲む 4つの 有値画素は、 当該サンプリ ング座標 3 0 1 , 3 0 3 を囲む一辺の長さが 2の正方形の各頂点に位置する。 例えば、 Rプレーンのサンプリ ング座 標 3 0 1 力 S ( u sR, V sR) = ( 1 0 0. 8, 1 0 1 . 4 ) であれば、 Rプ レーンにおいてこれを囲む 4つの画素 ( ud, V d) = ( 1 0 0 , 1 0 0 )
, ( 1 0 0 , 1 0 2 ) , ( 1 0 2 , 1 0 0 ) , ( 1 0 2, 1 0 2 ) が Rプレ 一ンの有値画素と なる。
当該有値画素の各画素値を R ( 1 0 0 , 1 0 0 ) , R ( 1 0 0 , 1 0 2
) , R ( 1 0 2 , 1 0 0 ) , R ( 1 0 2 , 1 0 2 ) で表すとする と 、 図 6 のよ う なバイ リ ニァ補間によって生成される Rプレーン上でのサンプリ ング座標 3 0 1 の補間画素値 R ( 1 0 0. 8 , 1 0 1 . 4 ) は 、 次の式 で表される。
R ( 1 0 0 . 8, 1 0 1 . 4 ) = 0. 6 * 0. 3 * R ( 1 0 0 , 1
0 0 ) + 0. 6 * 0. 7 * R ( 1 0 0 , 1 0 2 ) + 0. 4 * 0 3 * R
( 1 0 2 , 1 0 0 ) + 0. 4 * 0. 7 * R ( 1 0 2 , 1 0 2 )
一方、 Gプレーンは市松状に有値画素を持っため、 サンプ ング座標
3 0 2 を囲む 4つの有値画素は、 当該サンプ リ ング座標 3 0 2 を囲む一 辺の長さ が/" 2 の 4 5度傾いた正方形の各頂点に位置する。 この ½ α、
Gプレーンにおいてサンプリ ング座標 3 0 2力 S ( u sG, V sC) ( 1 0 1
. 0 , 1 0 1 . 4 ) であれば、 これを囲む 4つの画素 ( ud, V d ,) = ( 1 0 0 , 1 0 1 ) , ( 1 0 1 , 1 0 0 ) , ( 1 0 1 , 1 0 2 ) , ( 1 0 2, 1 0 1 ) が Gプレーンの有値画素となる。
当該有値画素の各画素値を G ( 1 0 0 , 1 0 1 ) , G ( 1 0 1 , 1 0 0 ) , G ( 1 0 1 , 1 0 2 ) , G ( 1 0 2 , 1 0 1 ) で表すとする と、 Gプ レーン上でのサンプリ ング座標 3 0 2の補間画素値 G ( 1 0 1 . 0 , 1 0 1 . 4 ) は、 次の式で表される。
G ( 1 0 1 . 0 , 1 0 1 . 4 ) = 0 . 7 * 0 . 3 * G ( 1 0 0 , 1 0 1 ) + 0 . 3 * 0 . 3 * G ( 1 0 1 , 1 0 0 ) + 0 . 7 * 0 . 7 * G ( 1 0 1 , 1 0 2 ) + 0. 3 * 0. 7 * G ( 1 0 2 , 1 0 1 )
色生成部 1 4 4 は、 サンプリ ング部 1 4 3 ' によ り補間生成された各 色プレーンの補間画素値を合成するこ とによ り 、 各画素値がそれぞれ複 数色の輝度情報を有したカラ一画像を生成する (図 1 0のステ ップ S 1 4 ) 。 さ らに、 色生成部 1 4 4は、 こ う して求めた R G B のカラ一情報 を Y U Vのカラー情報に変換し (Yは輝度情報、 U , Vは色情報) 、 U , Vの色情報に対して低周波フィルタを施す。 なお、 R G Bから Y U V への変換処理および U, Vに対する低周波フィルタの処理は、 公知の処 理を適用するこ とが可能である。 色生成部 1 4 4は、 以上の処理をカラ 一画像の全ての画素 (全てのサンプリ ング座標) について行い、 その結 果得られた変形カラー画像を視覚補正部 1 5 0に出力する。 視覚補正部 1 5 0以降の処理は上述した通りである。
色生成部 1 4 4 の処理で用いるサンプリ ング座標 3 0 1 , 3 0 2, 3 0 3の補間画素値 R ( u sR, V sR) , G ( u sG, v sG) , Β ( u sB, ν sB) は倍率色収差の色ズレが考慮されており 、 被写体上の同一部分を示すよ うになつている。 このため、 被写体上に白黒の境界がある部分で、 R , G , Bは同時に変化する。 これによ り 、 R, G , Bの各色成分を合成し たカラ一画像は、 シャープな輝度信号を得るこ とができる。 すなわち、 倍率色収差のない結像光学系で撮 ,影されるのと同等の力フ 画像を得る こ とができる
以上詳し <説明したよ う に、 第 2 の実施形態では 、 色モザィ ク画像か らカラー画像を生成する前段階と して 、 出力すベさ力ラ一画像の画素位 置から、 当該カラ一画像の画素位置に対応する色モザィ ク画像上の位置 であって、 倍率色収差を考慮して調 した位置であるサンプ y ング座標
3 0 1 , 3 0 2, 3 0 3 を色プレ一ン毎に算出する 。 そして 、 当該サン プリ ング座標 3 0 1, 3 0 2, 3 0 3 における画素値を 、 色モザイ ク画 像の画素値を用いた補間演算によ り生成するよ うにしている。
これによ り 、 色モザイ ク画像からカラー 像を生成する色補間処理と 当該力ラー画像の倍率色収差補正処理とを一度の補間演算によって実現 するこ とができる。 このため、 色モザイ ク画像から色収差補正された力 ラ一画像を生成する際の処理負荷を軽減するこ とができる と と もに、 従 来のよ う に補間処理を二重に行う こ とによる画質の劣化も抑制するこ と ができ
また 、 第 2の実施形態では、 色モザイ ク画像を R , G , Bの複数の色 プレーンに分解し、 色プレーン毎にサンプ y ング座標 3 0 1, 3 0 2,
3 0 3 の補間画素値を求めてから 、 各色プレ一ンの補間画素値を合成し て 1 画素に 3色の輝度情報が含まれるカラ一情報を生成するよ う にして いる。 このよ うにするこ とによ り 、 色プレ ―ン内に含まれる同一の色光 の画素値から単純な線形補間によ りサンプ ング座標 3 0 1 , 3 0 2 ,
3 0 3 の補間画素値を求めるこ とができ、 処理負荷を軽減するこ とがで さる。
なお 、 上記第 2 の実施形態では 、 色プ レ ―ン分解部 1 4 1 は、 色モザ イ ク 画像を R プ レーン、 G プ レーン 、 B プ レ ―ンの 3 つの色プ レーンに 分解する例について説明 したが、 れに限定されない。 例えば、 図 1 2 に示すよ ラ に、 互いに角を接する 2 緑色画素 G r, G b を各々異 なる色プレ一ンに分解する よ う に しても良い。 すなわち 、 色プレーン分 解部 1 4 1 は 、 A D変換部 1 3 0 よ り 出力 される色モザィ ク画像を、 R 成分の画 のみを取 り 出 した Rプレ ―ン、 G r 成分の画素のみを取 り 出 した G r プレ ―ン、 G b 成分の画素のみを取 り 出 した G b プレーン、 B 成分の画 のみを取 り 出 した B プレ ―ンの 4 つの色プレ ―ンに分解する m
o
口 に色生成部 1 4 4 は、 G Γ プレーンおよび G b プレ一ンでそ れぞれサンプ リ ング部 1 4 3 ' によ り 生成されたサンプ y ング座標 3 0
2 の画素値 (サンプリ ング値) を加算 して 、 1 画素內に R , G , B の各 色成分の輝度情報が含まれるカ ラ一情報を生成する。 例えば、 R成分や
B成分はサンプリ ング値をそのまま用い、 G成分は G r と G b と の平均 値を用いる o
また、 G成分を G r プレーンと G b プレーンと に分解してサンプリ ン グ座標 3 0 2 の画素値を算出する場 α 、 ァモザイ ク部 1 4 0 ' は、 G r プレーン と G b プレーンの各色プレ一ン上でそれぞれサンプリ ング部 1 によ り 生成された捕間画素値の差分を算出 し、 当該補間画素値の 差分に基づいて偽色の有無を判定する偽色判定部を更に備える よ う に し ても良い。
べィ ャ配列の単板カ ラ —イ メ ージセンサ 1 2 0 では、 ナイ キス ト周波 数近傍の白黒の縞模様に対し、 赤や青の偽色が発生する と い ラ 問題があ る。 これに対 し 、 G r プレーン上で求めたサンプリ ング座標 3 0 2 の補 間画素値と G b プ レーン上で求めたサンプリ ング座標 3 0 2 の補間画素 値と の差分を と る こ と で 、 縞模様上での偽色の有無を検出する こ と がで き、 偽色がある場合にはこれを抑制する こ と ができ る。
すなわち、 G r も G b も本来は同 じ G色フ イノレタなので、 G Γ プレ ― ンおよび G b プ レーンで求めた補間画素値は双方と も同じ値になるはず である。 しかし、 偽色が生じている と、 G r プレ一ンおよび G bプレ一 ンで求めた補間画素値に差分が生じる。 よつて 、 補間画素値の差分を見 るこ とによ り 、 偽色の発生の有無を検出するこ とができる た、 色生 成部 1 4 4は、 R G B のカラー情報を Y U Vの力ラ一情報に変換する際 に、 図 2 2の式によって変換処理を行う。
(第 3の実施形態)
次に、 本発明の第 3 の実施形態を図面に基づいて説明す 本発明に 係る画像処理装置を実施した第 3の実施形態に るカラ一撮像装置の構 成は、 図 8 と同様である。 また、 デモザイ ク部 1 4 0 ' の機能構成は、 図 9 と同様である。 ただし、 デモザイク部 1 4 0 , が備える座標変換部
1 4 2 ' の処理内容が第 2の実施形態と異なつてレ、る。 以下では、 第 2 の実施形態と異なる部分を中心に説明する。
第 3 の実施形態において 、 座標変換部 1 4 2 ' は 、 色プ レーン分解部
1 4 1 によ り分解された複数の色プ レーン毎に、 色プ レー ンによって値 の異なる色収差係数と、 力ラー画像に対する画像変形を表す画像変形係 数と を用いて、 カラー画像の画素位置から、 色収差補正および画像変形 が施された場合のカラー画像の面素位置に対応する色モザィ ク画像上の サンプリ ング座標であつて色プレーン毎に異なるサンプリ ング座標を算 出す
ここ では画像変形の例と して、 結像光学系 1 1 0の歪曲収差を補正す るための変形処理を例に挙げて説明する。 結像光学系 1 1 0に歪曲収差 カ ある と、 結像光学系 1 1 0 を経て単板力ラ一ィ メージセンサ 1 2 0 の tR像素子 1 2 2上に結像した被写体像は劣化したものとなる。 すなわち
、 被写体上では直線であるも のが歪曲収差によ り 曲線となった画像とな る。 第 3の実施形態では、 倍率色収差による画質劣化に加えて、 歪曲収 差による画質劣化も同時に補正する。
そのために座標変換部 1 4 2 ' は、 色モザイ ク画像から生成される力 ラー画像の画素位置から、 倍率色収差および歪曲収差の補正を考慮した 色モザイ ク画像上のサンプリ ング座標を色プレーン毎に算出する。 具体 的には、 座標変換部 1 4 2 ' は、 x y座標 ( x d, y d) に対して、 結像光 学系 1 1 0の倍率色収差と歪曲収差の補正を考慮した座標変換を以下の 式のよ うに行う (ただし、 r 2= X d 2+ y d 2) 。
x 1 — X d ( k + k 1 r 2+ k 2 r 4) + 2 p , x < j y d + p 2 ( r 2 + 2 x d 2) y 1 Υ d ( k + k . Γ 2+ k 2 r 4) + 2 p 2 x < i y d + p i ( r 2 + 2 y d 2) ( ;一 X d ( 1 + k ( r : !+ k 2 : r 4) + 2 p , x d y d + p 2 ( r 2 + 2 x d : 2) y , y d ( 1 + k , r : !+ k 2 r 4) + 2 p 2 x d y d + p i ( r 2 + 2 y d : 2) 1 ― X d ( k B + k 1 r 2+ k 4) + 2 p , X c i y a + p 2 ( r 2 + 2 x d 2) y i y d ( k B + k 1 Γ 2+ k 2 r 4) + 2 p 2 X c i y d + p i ( r 2 + 2 y d 2) こ こで、 { k ,, k 2} は結像光学系 1 1 0 の歪曲収差を示す係数であり 、 k ,は 3次収差係数、 k 2は 5次収差係数を示している。
なお、 結像光学系 1 1 0の収差によ り歪曲しているカラー画像を補正 するために画像変形を行う場合には、 色モザィ ク画像から生成された歪 曲しているカラー画像に対して非線形の座標変換を行う。 歪曲収差を含 むカラー画像に対して非線形の座標変換を行う方法は既知であり 、 上述 の歪曲収差係数 { k ,, k 2} はシミ ュ レーショ ン等によ り求めることが可 能である。 この歪曲収差係数 { k ,, k 2} が、 本発明の画像変形係数に相 当する。
これらの x y座標 ( X R, y R) ( x c, y c) ( x B. y B) から色モザィ ク画像上の画素座標 ( u sR, V sR) ( u sG, v sG) ( u sB, v sB) を求める 演算は、 第 2 の実施形態において x y座標 ( X , y ) から色モザイ ク画 像上の画素座標 ( u s , s) を求めた演算と同様であり 、 以下の式のよ う に求められる。
u s 1 0 0 0 * X 7 9 9 . 5
1 0 0 0 * 5 9 9 . 5
s 1 0 0 0 * G + 7 9 9 . 5
v s 1 0 0 0 * y , G + 5 9 9 . 5
1 0 0 0 * 7 9 9 . 5
1 0 0 0 * 5 9 9 . 5
次に、 サンプリ ング部 1 4 3 ' は、 以上のよ う にして座標変換部 1 4 2 ' によ り算出された色プレーン毎のサンプリ ング座標における画素値 (サンプリ ング値) を、 色プレーン内に含まれる同一の色光の画素値か ら補間生成する。 また、 色生成部 1 4 4は、 サンプリ ング部 1 4 3 ' に よ り補間生成された各色プ レーンの補間画素値を合成するこ とによ り 、 各画素値がそれぞれ複数色の輝度情報を有したカラー画像を生成する。 以上詳しく説明したよ う に、 第 3の実施形態によれば、 色モザイ ク画 像からカラー画像を生成する色補間処理と、 結像光学系 1 1 0 の倍率色 収差を補正する処理と、 さ らに結像光学系 1 1 0 の歪曲収差を補正する ためのカラ一画像の画像変形処理とを一度の補間演算によって実現する ことができる。 このため、 倍率色収差が補正され、 かつ、 歪曲収差も補 正された変形力ラー画像を色モザイ ク画像から生成する際の処理負 軽減するこ とができる と と もに、 補間処理を何回も行ラ こ とによる の劣化も抑制するこ とができる。
なお、 上記第 3の実施形態では 、 画像変形の例と して 、 像光学
1 0の歪曲収差を補正するための変形処理を例に挙げて説明 したが れに限定される ものではない。 例えば、 画像変形は、 テジタルズ一 理によるカラー画像の拡大または縮小、 手ブレ補正処理による力ラ 像の回転な どであって も良い。 この場合、 拡大、 縮小、 回転は例えばァ フ ィ ン変換で表すこ と が可能であ り 、 そのァフ ィ ン変換を表す係数を画 像変形係数と して用いる。
(第 4 の実施形態)
次に、 本発明の第 4 の実施形態を図面に づいて説明する。 図 1 3 は
、 本発明に係る画像処理装置を実施 した第 4 の実施形態によ る力 ラ ~ fe 像装置 1 0 0,, の構成例を示す図である。 なお、 図 1 3 において 、 図 8 に示 した符号と 同一の符号を付したものは同一の機能を有する ちのであ るので、 こ こでは重複する説明を省略する。
図 1 3 に示すよ う に、 第 4 の実施形態の力 ラー撮像装置 1 0 0 は、 結像光学系 1 1 0 " 、 単板カ ラ一イ メ ージセンサ 1 2 0、 A D変換部 1
3 0 、 デモザィ ク部 1 4 0 " 、 視覚補正部 1 5 0 、 圧縮部 1 6 0 、 記録 部 1 7 0 、 ブ レ検出部 1 8 0 、 収差係数設定部 1 9 0 、 キヤ リ ブレ一シ ョ ン部 2 0 0 および収差係数テーブル記憶部 2 1 0 を備えて構成されて レヽる 。 この う ち 、 デモザイ ク部 1 4 0 " お び収差係数設定部 1 9 0 が 本発明の画像処理装置に相当する。
図 1 3 において、 撮影された被写体は、 像光学系 1 1 0 " を経て単 板力 ラーイ メ ージセンサ 1 2 0 の撮像素子 1 2 2 上に結像する この際
、 /(¾像光学系 1 1 0 " が持つ各種の収差に り 、 結像 した被写体像は劣 化したもの と なる。 例えば、 歪曲収差によ り 被写体上では直線である も のが曲線と な り 、 倍率色収差によ り 色成分毎に撮像素子 1 2 2 上での結 像ズレ (色ズ レ ) が生じた画像と なる。 なお 、 第 4 の実施形態において 結像光学系 1 1 0 " は、 焦点距離 (ズーム ) や被写体距離 (フ ォ ―力ス
) といった レンズステ一 ト を変更可能なものとする。
ブレ検出部 1 8 0 は、 カ ラー撮像装置 1 0 0 " のブレを検出 し 、 ブレ を補正するためのブレ補正係数 { z , Θ , d x, d y } をデモザイ ク部 1 4
0 に設定する。 ブ レの検出にはジャ イ ロ セ ンサを用いる方式や、 複数 枚撮影 した面像間の特徴点の変化量な どを測定する方式があるが、 本発 明はブレの検出法を限定する ものではない。 こ こ で、 カ ラー撮像装置 1
0 0 の前後方向のブレに伴 う 被写体像の大き さの補正値を z 、 口一ル 軸のブ レに伴 う 被写体像の回転の補正値を 0 、 左右方向あるいはョ ーの ブレに伴 う被写体像の左右位置の補正値を d X 、 上下方向あるいはピク チのブレに伴 う被写体像の上下位置の補正値を d y と している。 なお 、 係数 z にはデジタルズームの倍率を含めても良レ、。
収差係数設定部 1 9 0 は、 本発明の係数設定部に相当する ものである o 収差係数設定部 1 9 0 は、 結像光学系 1 1 0 " の レンズステー ト を検 出 し 、 検出 した レンズステー 卜に応じた適切な収差係数 { k k 2, p ,, P
2 J k , R, k B} を収差係数テーブル記憶部 2 1 0 から読み出 し、 デモザィ ク 部 1 4 0 " に設定する。 こ こ で、 { k い k 2, ρ ,, ρ 2} は結像光学系 1 1
0 の歪曲収差を示す係数であ り 、 { k ,, k 2} は放射線方向の歪み、 {
P 1, P 2} は接線方向の歪みを示 している。 { k R, k B} は結像光学系 1
,,
1 0 の倍率色収差を示す係数であ り 、 k R は G成分に対する R成分の倍 率 、 k 成分に対する B成分の倍率を示している。
なお 、 レンズステー トが変化する場合、 収差係数設定部 1 9 0 は、 変 化後の レンズステー ト を検出 し、 その レンズステー ト に対応する収差係 数を収差係数テーブル記憶部 2 1 0 よ り 読み込み、 デモザイ ク部 1 4 0
,,
に Ιχ定する。 レンズステー ト の検出は、 例えば、 レンズステー ト の制 御を行つているカ ラ一撮像装置 1 0 0 " の コ ン ト ローラ (図示せず) か ら 、 像光学系 1 1 0 " に設定された レンズステー ト の情報を受ける こ と によつて行う こ と が可能である。
と ろで、 あ らゆる レンズステー ト に対する収差係数をテ一ブルに記 録しておく のは非現実的である。 そこで本実施形態では、 収差係数テー ブル記憶部 2 1 0 には限られたレンズステー トに対する収差係数だけを 記録しておく。 例えば、 焦点距離、 被写体距離を各々 3パター ン、 計 9 つの レンズステー ト とそれに対応する収差係数値とを記録しておく。 こ の収差係数テーブルの収差係数 k , に関わる部分を例示すると、 図 1 4 の よ うになる。
図 1 4 において、 例えば結像光学系 1 1 0 " に設定されている レンズ ステー ト が被写体距離 2 . O m、 焦点距離 3 5 m mである とする と、 図 1 4 の収差係数テーブルに適合する値はない。 そこで、 収差係数設定部 1 9 0は、 収差係数テーブルよ り被写体距離が M i d : 1 . 0 m、 焦点 距離が W i d e : 2 8 m mの収差係数 0 . 0 8 と、 被写体距離が F a r : I n f 、 焦点距離が W i d e : 2 8 m mの収差係数 0 . 0 5 と、 被写 体距離が M i d : 1 . 0 m、 焦点距離が M i d : 5 O m mの収差係数 0 . 0 2 と、 被写体距離が F a r : I n f 、 焦点距離が M i d : 5 0 m m の収差係数 0 . 0 0 と を読み出し、 これ らの 4 つの収差係数を補間する こ こ で、 焦点距離および被写体距離に関しては逆数に対して補間する のが望ま しく 、
k , = ( ( 1 / 2 . 0 — l Z l n f ) ( l / 3 5 _ l Z 5 0 ) * 0 . 0 8 + ( 1 / 1 . 0 - 1 / 2 . 0 ) ( 1 / 3 5 - 1 / 5 0 ) * 0 . 0 5 + ( 1 / 2 . 0 - 1 / 1 n f ) ( l / 2 8 - l / 3 5 ) * 0 . 0 2 + ( 1 / 1 . 0 — 1 / 2 . 0 ) ( 1 / 2 8 - 1 / 3 5 ) * 0 . 0 0 ) / ( 1 / 1 . 0 — l I n f )/ ( 1 / 2 8 - 1 / 5 0 )
= 0 . 0 4
のよ う に計算した補間値をデモザイ ク部 1 4 0 " に設定する。 他の収差 係数 { k 2, p p 2, k R, k B} に対しても、 同様に補間値を計算し、 デモザ イ ク部 1 4 0 " に設定する。
キャ リ ブ レーシ ョ ン部 2 0 0 は、 A D変換部 1 3 0 よ り 出力 されるデ ジタル画像信号に基づいて、 収差係数テーブル記憶部 2 1 0 に記憶する 収差係数テーブルの値を生成する。 キャ リ ブ レーシ ョ ン部 2 0 0 は、 収 差係数テーブルに記録すべき複数の レンズステー ト に対してそれぞれ収 差係数 { kい k 2, p い p 2, k R, k B} を求め、 これらを収差係数テーブル記 憶部 2 1 0 に記録する。
なお、 キャ リ ブ レーシ ョ ン部 2 0 0 はカ ラー撮像装置 1 0 0 " の一部 と して内部に備えても良いが、 カ ラ一撮像装置 1 0 0 " と は別体のキヤ リ ブレ一シ ヨ ン装置と して構成しても良い。 別体の装置とする場合、 力 ラ一撮像装置 1 0 0 " は、 撮影画像 (A D変換部 1 3 0 よ り 出力 される デジタル画像信号) をキャ リ ブ レーシ ョ ン装置に出力 した り 、 キヤ リ ブ レ一シ ヨ ン装置から収差係数を入力 した り する通信手段を備える。
デモザイ ク部 1 4 0 " は、 ブレ検出部 1 8 0 によ り 設定されるブレ補 正係数 { z , 0 , d x , d y } および収差係数設定部 1 9 0 によ り 設定され る収差係数 { kい k 2, p い p 2, k R, k B} に基づく 補正を行いつつ、 色モザ ィ ク画像からカラー画像への変換を行う。
図 1 5 は、 デモザイ ク部 1 4 0 " の機能構成例を示すブロ ッ ク図であ る。 図 1 5 に示すよ う に、 デモザイ ク部 1 4 0 " は、 その機能構成と し て、 色プレーン分解部 1 4 1 " と 、 座標変換部 1 4 2 " と 、 サンプ リ ン グ部 1 4 3 " と、 色生成部 1 4 4 " と を備えている。
色プ レーン分解部 1 4 1 " は、 A D変換部 1 3 0 よ り 出力 される色モ ザィ ク画像を、 同一の色光の画素値だけを含む複数の色プ レーンに分解 する。 例えば、 色プ レーン分解部 1 4 1 " は、 図 1 2 に示 したよ う に、 R成分の画素のみを取 り 出 した Rプレーン、 G r 成分の画素のみを取 り 出 した G r プレーン、 G b成分の画素のみを取り 出 した G b プレーン、 B成分の画素のみを取 り 出 した B プレーンの 4 つの色プレーンに分解す る。
座標変換部 1 4 2 " は、 上述のブレ補正係数 { z , 0, d x , d y } およ び収差係数 { k k 2, p p 2, k K, k B} を用いて、 色モザイ ク画像から生 成されるカ ラ一画像の画素位置から、 色収差補正、 歪曲収差補正および ブレ補正が施された場合のカ ラー画像の画素位置に対応する色モザイ ク 画像上のサンプ リ ング座標を算出する。 以下に、 サンプリ ング座標の算 出手順を詳しく 説明する。
まず x y座標系 と して、 原点を画像中心、 最大像高 (原点からの最大 距離) を 1 と し、 原点から画面右方向に正の X座標を と り , 原点から画 面下方向に正の y座標を と る ものとする。 また、 出力するカ ラー画像が 1 6 0 0 X 1 2 0 0 の正方画素からなる もの と する。 この場合、 カ ラ一 画像に対して、 u v座標系の画素座標 ( u d, V d) を図 2 のよ う に画面左 上から右方向に ( 0, 0 ) , ( 1 , 0 ) , ( 2 , 0 ) · · · 、 次の行を ( 1 , 0 ) , ( 1 , 1 ) , ( 2 , 1 ) · · · と割 り 振れば、 画素座標 ( u d, V d) に対応する x y座標 ( x d, y d) は、 図 2 3 に示す関係式で表され る。
座標変換部 1 4 2 " は、 この x y座標 ( x d, y d) に対して、 まず上述 のブレ補正係数 { z, 0 , d x , d y } を図 2 4 に示す式のよ う に適用 し、 ブレ補正後の x y座標 ( χ ' , y ' ) を得る。
さ らに、 座標変換部 1 4 2 " は、 歪曲収差に係る係数 { k ,, k 2, P l, p 2} を次式のよ う に適用 し、 G r プレーンおよび G b プレーンにおける X y座標 ( X G, y G) を求める。
X G= X ' ( l + k , r 2+ k 2 r 4) + 2 p 1 x ' y ' + p 2 ( r ' 2 + 2 x ' 2 ) y G= y ' ( 1 + k , r 2+ k 2 r 4) + 2 p 2 x ' y ' + p , ( r ' 2 + 2 y ' 2 ) (ただし、 r ' 2= x ' 2+ y ' 2) さ らに、 座標変換部 1 4 2 " は、 結像光学系 1 1 0 " の倍率色収 に 係る係数 { k K, k B} を考慮 した色プ レーン間の座標の変換を図 2 5 に示 す式のよ う に行 う こ と によ り 、 Rプレーンおよび B プレーンにおける X y座標 ( X K, y R) ( x B, y B)をそれぞれ求める。 こ こ で、 d R x , d R y は Gプレーンを基準に した際の Rプレーンの平行ズレ量、 d Bx , d By は Gプレーンを基準に した際の Bプ レーンの平行ズレ量を示す係数である 一方、 色モザイ ク画像も 1 6 0 0 X 1 2 0 0 の正方画素と し、 u v座 標系の画素座標 ( u s, V s) が上述のカラー画像と 同様に割 り振られてい る とする と、 各色プレーン上の x y座標 ( X R, y R) ( x G) y G) ( x B. y B) に対応する色モザイ ク画像上の画素座標 ( u sR, V sR) ( u sG> v sG ) ( u sB, v sB) は、 図 2 6 の式に示すよ う になる。
上述の計算の結果、 これらの画素座標 ( u sR, V sR) ( u sG, v sG) ( u sB, v sB) は何れも整数値になる と は限らず、 一般には非整数と なる。 こ の色モザイ ク画像上の画素座標 ( u sR, V sR) ( u sG) v sG) ( u sB, v sB ) が色プレーン毎のサンプリ ング座標である。
サンプ リ ング部 1 4 3 " は、 色プレーン分解部 1 4 1 " によ り 分解さ れた複数の色プ レーン毎に、 座標変換部 1 4 2 " によ り 算出された色プ レーン毎のサンプリ ング座標における画素値 (サンプ リ ング値) を、 色 プ レーン内に含まれる同一の色光の画素値から捕間生成する。 色生成部 1 4 4 " は、 サンプ リ ング部 1 4 3 " によ り 補間生成された各色プレー ンの補間画素値を合成する こ と によ り 、 各画素値がそれぞれ複数色の輝 度情報を有したカ ラ一画像を生成する。 こ のサンプ リ ング部 1 4 3 " お よび色生成部 1 4 4 " の処理内容は、 第 2 実施形態または第 3 の実施形 態と同様である。
以上詳 し く 説明 したよ う に、 第 4 の実施形態によれば、 色モザイ ク画 像からカラ一画像を生成する色捕間処理、 結像光学系 1 1 0 " の倍率色 収差を補正する処理、 結像光学系 1 1 0 の歪曲収差を補正するための カラ一画像の変形処理、 結像光学系 1 1 0 のレンズステ一トに応じた 収差を補正するための処理 、 およびカラ ― 像装置 1 0 0 の手ブレを 補正する処理を一度の補間演算によって実現するこ とがでさる このた め、 倍率色収差が補正され 、 力、つ、 歪曲収差ゃ手ブレなども補正された カラ一画像を色モザィ ク画像から生成する際の処理負荷を軽減するこ と ができる と と もに、 補間処理を何回も行ラ こ とによる画質の劣化ち抑制 することができる。
なお、 上記第 4 の実施形態では、 座標変換部 1 4 2 " がサンプリ ング 座標を算出する際に使用する係数と して、 ブレ補正係数 { z , 0 , d x, d y } および収差係数 { kい k 2, pい p 2 , k R, k B } を用いる例について説明 したが、 これらを全て用いなく ても良い。 すなわち、 色収差係数 { k R, k B } を用いるこ とは必須と して、 残りの係数は任意に組み合わせて用い ることが可能である。
また、 上記第 4の実施形態では、 キヤ リ ブレ シヨ ン部 2 0 0 を設け て収差係数テ―ブノレの値を可変とする例につレ、て説明 したが、 これに限 定されない。 例えば 、 キャ リ ブ レーシ ョ ン部 2 0 0は設けずに収差係数 テ一ブルの値を固定とするよ うにしても良い。
(第 5の実施形態)
次に、 本発明の第 5の実施形態を図面に基づレ、て説明する。 図 1 6 は
、 第 5 の実施形態によるカラ一撮像装置 1 0 0 Aの構成例を示す図であ る なお 、 図 1 6 において、 図 8 に示した符号と同一の符号を付したち のは同一の機能を有するものであるので、 ここでは重複する説明を省略 する。 第 5の実施形態におけるカラ一撮像装置 1 O O Aは、 所定のチャー ト C Hを撮影するこ とによって色収差係数を検出する色収差量検出装置 2 0 O Aを備え、 色収差量検出装置 2 0 O Aによって検出された色収差係 数が、 第 2、 第 3の実施形態における座標変換部 1 4 2 ' や、 第 4 の実 施形態における収差係数設定部 1 9 0に格納され、 第 2〜第 4の実施形 態と同様に、 デモザイ ク部 1 4 0 ' , 1 4 0 " において、 カラー画像の 画素位置から、 色収差補正が施された場合のカラー画像の画素位置に対 応する色モザイ ク画像上のサンプリ ング座標を、 色プレーン毎に算出す るよ うに構成されている。 なお、 チャー ト C Hが、 本発明の請求の範囲 第 2 3項および第 2 4項における所定の画像に相当する。
以下、 色収差量検出装置 2 0 O Aの構成と作用を、 図 1 6〜図 1 9 に 基づいて説明する。 図 1 7が同実施形態における色収差量検出の説明図 であって、 図 1 7 ( a ) がチャー トの形態を表す図、 図 1 7 ( b ) が撮 像素子 1 2 2に対するチャー トの配置を表す図である。 また、 図 1 8が同実施形態における色収差量検出の説明図であって、 図 1 8 ( a ) ( b ) が交点検出の際の説明図、 図 1 8 ( c ) が交点毎に エッジを検出する際の説明図、 図 1 8 ( d ) が交点に対してサンプリ ン グ画素列を設定する際の説明図、 図 1 9が同実施形態におけるエッジ検 出の説明図である。
まず、 図 1 6 に表したよ うに、 本実施形態の色収差量検出装置 2 0 0
Aは、 チヤ一 ト C Hを撮影して得られるデジタル画像信号を用いて結像 光学系 1 1 0の色収差係数を算出する。
チャー ト C Hは、 図 1 7 ( a ) ( b ) に表したよ う に、 単板力ラーイ メ一ジセンサ 1 2 0における撮像素子 1 2 2の画素配列に対して、 第一 撮像模様 P 1 と第二撮像模様 P 2の配列が、 傾斜角 ひ の分だけ傾斜して いる。 また、 本実施形態では、 一つの撮像模様の領域が、 撮像素子 1 2
2 で読み込まれる約 1 0 0 ピクセル (画素) に相当する。 次に、 色収差量検出装置 2 0 O Aは、 図 1 6 に表したよ う に、 A D変 換部 1 3 0から入力されたデジタル画像信号 (所謂、 画素の輝度を表す 画素信号である) を R G Bの色毎に記憶するフ ィ ール ドメ モ リ 2 2 1 、 フィ ール ドメ モ リ 2 2 1 に記憶された画素信号に基づいて第一撮像模様 P 1 及び第二撮像模様 P 2の交点を複数検出する交点検出処理部 2 2 5 、 交点検出処理部 2 2 5 で検出された交点周 り において R G B毎に第一 撮像模様 P 1 と第二撮像模様 P 2 と のエッジ位置を検出する R G Bエツ ジ検出処理部 2 2 6 、 R G Bエッジ検出処理部 2 2 6 で検出されたエツ ジ位置を交点に対応付けて記憶するエッジ位置記憶部 2 2 8 、 エッジ位 置記憶部 2 2 8 に記憶されたエッジ位置から色収差係数を算出する収差 係数算出部 2 2 9 、 C P U ( C e n t r a l P r o c e s s i n g U n i t ) 2 3 0 、 R O M ( R e d O n l y M e m o r y ) 2 3 1 等によって構成され、 C P U 2 3 0が R O M 2 3 1 に格納された制御用 プロ グラムに従って、 当該色収差量検出装置 2 0 0 Aの各処理を制御す る。
フィール ドメモリ 2 2 1 は、 B a y e r配列に対応付けて、 赤 (R ) の画素信号を記憶する R フィ ール ドメ モ リ 2 2 2 と、 緑 (G ) の画素信 号を記憶する Gフィール ドメモ リ 2 2 3 と、 青 ( B ) の画素信号を記憶 する B フィール ドメモ リ 2 2 4 と、 によって構成されている。
交点検出処理部 2 2 5は、 図 1 8 ( a ) ( b ) に表したよ う に、 注目 画素を中心と した所定の範囲における画素値を用いて輝度勾配を算出 し 、 その輝度勾配が最も大き く なる注目画素の位置を交点 I n t と して検 出する。 こ こでは、 図 1 8 ( b ) に表したよ う に、 注目画素を中心と し た縦横 5画素を設定し、 画素位置に応じた重み付けを付加して交点位置 を検出している。 つま り 、 注目画素を中心と した上下左右の画素値に対 して図 1 8 ( b ) に示す係数を乗算してその結果を集計し、 集計結果の 絶対値を注目 画素の評価値と し、 評価値が所定の閾値を超えた際の注目 画素の位置を交点 I n t と し、 図 1 8 ( a ) に表 したよ う に、 複数の交 点 I n t をマ ト リ ク ス状に検出 している。 また、 本実施形態では、 交点 I n t がマ ト リ タ ス状に等間隔で表れる よ う に、 第一撮像模様 P 1 及び 第二撮像模様 P 2を配置している。
R G Bエ ッ ジ検出処理部 2 2 6は、 図 1 8 ( c ) に表 したよ う に、 R G Bの色毎に交点 I n t を介して上下及び左右に位置する複数の画素列 H s 、 V s を所定のサンプリ ングライ ン長で走査 して、 順次、 画素値を 取得する と と もに隣接するサンプリ ング位置に対して最も画素値の変化 量が大きいサンプリ ング位置をエッジと して検出する。
詳し く は、 図 1 9の曲線 I nに表 したよ う にサンプ リ ング毎に各画素 の輝度 (画素値) を求め、 曲線 S Lに表 したよ う にサンプリ ングによ つ て求め られた画素値にも とづいて画素値の変化量 (勾配 S L ) を算出 し 、 最も変化量 (勾配 S L ) が大き く 表れる位置 E Pをエ ッ ジと して検出 する。
また、 エ ッ ジ E Pを求める際には、 図 1 8 ( c ) に表したよ う に、 交 点 I n t を介 して、 上下の画素範囲において、 夫々 、 複数列のサンプ リ ング (H s ) を行って、 列毎にエ ッ ジを検出 し、 次いで、 上部で検出 さ れたエ ッ ジ位置の平均値と 下部で検出されたエ ッ ジ位置の平均値と の平 均を算出 し、 交点 I n t における左右方向のエッジ位置とする。
また、 交点 I n t を介 して左右の画素範囲においても、 夫々 、 複数列 のサンプ リ ング (V s ) を行って、 列毎にエ ッ ジを検出 し、 次いで、 左 部で検出 されたエ ッ ジ位置の平均値と右部で検出 されたエ ッ ジ位置の平 均値と の平均を算出 し、 交点 I n t における上下方向のエ ッ ジ位置とす る。
また、 サンプ リ ングは、 同色の画素毎に行われ、 左右方向に沿ってサ ンプ リ ング H s を行 う 際には、 図 1 7 ( b ) に表 したよ う に、 左右方向 のサンプリ ング長 S L ( 1 1 ) と 、 サンプリ ングの上下方向の列の数を 表すサンプリ ング数 S N ( 4 ) が、 所要の検出精度に応 じて予め定め ら れている。 また、 上下方向に沿ってサンプ リ ング V s を行 う 際にも、 上 下方向のサンプリ ング長とサンプリ ング数が予め定められている。
また、 図 1 8 ( d ) に表 したよ う に交点 I n t を介 して上方に延出す るエ ッ ジの左右方向のエ ッ ジ位置を検出する際、 サンプリ ング H s 1 の 位置が交点 I n t に近接しすぎる と 、 エ ッ ジ E Hが傾斜しているので交 点 I n t から左方に位置する撮像模様 P I _ 2 の影響を受けて、 エ ッ ジ 検出が困難になるため、 サンプリ ングライ ン H s 1 と 交点 I n t との間 に適切な間隔 S をもたせる こ とが好ま しい。
間隔 S は、 例えば図 1 8 ( d ) に表 したよ う に幾何学的に求める こ と ができ る。 すなわち、 交点 I n t からの間隔 S は、 エ ッ ジのぼやけ量 E 、 傾斜角度 α 、 サンプリ ングライ ン長 S L を既知の量と して、 L - S L ノ 2 、 S = (W + L ) X t a n a 、 の算式によ って求める こ と ができ る 。 つま り 、 サンプリ ング H s l の開始位置が、 図 1 8 ( d ) 中の撮像模 様 P 1 — 2 内に入る こ と な く 、 P 1 — 2 からエ ッ ジのぼやけ量 Eだけ離 間するよ う に、 間隔 S を求める と よい。
次に、 エッ ジ位置記憶部 2 2 8 は、 各交点 I n t j ( j は各交点に付与 した連番 1 , 2 …である) において、 R G Bエ ッ ジ検出処理部 2 2 6 で 検出 した色毎のエ ッ ジ位置に基づいて、 G (緑) 、 R (赤) 、 B (青) のエ ッ ジ位置を、 夫々左右方向 ( u方向) 及び上下方向 ( 方向) に対 応付けて、 ( U c i , G j 、 (. U j , R j ) 、 ( U j , B j ) と し て記憶, Ό ο
次に、 収差係数算出部 2 2 9 は、 エ ッ ジ位置記憶部 2 2 8 に記憶され たエ ツ シ位置 U G j , V G j ) 、 、 U R j , V R j 、 、 U B j , j ; を 用いて、 色収差係数 k R、 k „を算出する。
詳しく は、 まず、 第 2 の実施形態と同じ様に、 x y座標系と して、 原 点を画像中心、 最大像高 (原点からの最大距離) を 1 と し、 原点から画 面右方向に正の X座標をと り, 原点から画面下方向に正の y座標をと る ものとする。 そして、 第 2 の実施形態と同じく 、 6 4 0 X 4 8 0の正方 画素からなるカラー画像に対して、 u v座標系の画素座標 ( u d, V d) を 図 2のよ う に画面左上から右方向に ( 0 , 0 ) , ( 1 , 0 ) , ( 2 , 0 ) . · . 、 次の行を ( 1 , 0 ) , ( 1 , 1 ) , ( 2 , 1 ) . . . と割り振り 、 u v座標系の画素座標 ( 3 1 9 . 5 , 2 3 9 . 5 ) を x y座標系の原 点と し、 色毎に、 画素座標 ( u d, V d) に対応する x y座標 ( x d, y d) を算出する。
こ こでは、 第 2 の実施形態で用いた x d= ( u d- 3 1 9 . 5 ) / 4 0 0 、 y d= ( v d- 2 3 9 . 5 ) 4 0 0 の演算式において、 x d及び y dを、 x G i及び y G j 、 x R j及び y R j 、 x B j及び y B jに置き替えると共に、 u d及び V dを 、 U j及び V j -> u R』及び V R j 、 U B j及び V B j に置き 替えて、 次の演算式を用い、 R G . B毎の x y座標を算出する。
X G j ( U G i - 3 1 9 . 5 ) Z 4 0 0
y G j ( G i - 2 3 9 . 5 ) / 4 0 0
X R i " ( U R i - 3 1 9 . 5 ) / 4 0 0 y j 一 ( R i - 2 3 9 . 5 ) / 0 0
X B j 一 ( U B i - 3 1 9 . 5 ) / 4 0 0 y j ' ( B i - 2 3 9 . 5 ) / 4 0 0
次いで 、 Rの色収差係数 k 、 k ∑ i ( X i 2 + y κ i 2 ) / ∑ i (
A R j Λ G j + y i G j ) の演算式で算出し、 Bの色収差係数 k Bを、 k B
2
= ∑ i ( X B j + y , ∑ i ( x , j X j + y a j y j ) の演算式で算 出する。 次に、 こ こで算出された色収差係数 k K及び k Bをデモザイク部 1 4 0
, 1 4 0 " における座標変換部 1 4 2 ' , 1 4 2 " に格納し 、 以下 、 第 2〜第 4 の実施形態と同様に、 座標変換部 1 4 2 ' , 1 4 2 におレ、 て 、 ¾像画像に対して色収差の補正を考慮したサンプリ ング座標を算出 する。
また、 図 2 7 に示す式のよ うに、 図 2 5 に示した式に対して収差係数 算出部 2 2 9で算出した X y座標及び収差係数を代入し、 色収差行列を 求めるこ とができる。 なお、 図 2 7 の式において、 最右辺に位置する行 列は、 擬似逆行列を示している。
以上のよ う に、 第 5の実施形態のカラ一撮像装置 1 0 O Aに れば 、 マ 卜 リ クス状の交点位置に対応付けて色毎のエッジ位置を検出できて 、 この検出された色毎のエッジ位置にもとづき色収差係数 ( k R、 k B ) を 算出でき 、 この算出された色収差係数を用いて、 色モザイ ク画像から力 ラ一画像を生成する色補間処理 (デモザイク処理) と、 結像光学系の倍 率色収差を補正する処理と を一度の補間演算によって実現するこ とがで き -3
以上、 本発明の一実施形態について説明 したが、 本発明は、 肓 IJ記実施 形態に限定されるものでなく 、 種々の態様をとることができる。
例えば、 上記第 1 〜第 5 の実施形態では、 カラーフ ィ ルタアレイ 1 2 1 が原色べィャ配列の場合を例にとって説明したが、 これに限定されな い。 例えば、 補色配列のカラ一フィルタアレイを用いても良い。
また、 第 5 の実施形態では、 撮像装置 1 0 O Aに色収差量検出装置 2 0 0 を備えたが、 色収差量検出装置 2 0 O Aを撮像装置 1 0 O Aと一体 に構成するこ となく外部機器と して備え、 この外部機器に撮影画像を出 力する画像出力部と、 外部機器から色収差係数を入力する係数入力手段 とを備え、 入力 した色収差係数を座標変換部 1 4 2 ' , 1 4 2 " または 収差係数設定部 1 9 0に格納するよ うに構成してもよレ、。
また、 以上に説明した第 1 〜第 5の実施形態によるデモザィ ク処理の 手法は、 ハ — ドウエア構成、 D S P (Di igital ϊ> 1 gna 1 Processorノ 、 ソ フ 卜 ゥェァの何れによっても実現するこ とが可能である。 例えばソフ 卜 ウ ェ アに つて実現する場合、 上記第 1 〜第 5の実施形態のデモザィ ク 部 1 4 0 , 1 4 0 ' , 1 4 0 " (画像処理装置) は、 実際にはコンピュ
—タの C P Uあるいは M P U、 R AM, R OMなどを備えて構成され 、
R A Mや R O Mに記憶されたプロ グラムが動作するこ とによって実現で さる。
その他 、 上記第 1 〜第 5の実施形態は 、 何れも本発明を実施するにあ たっての具体化の一例を示したものに過ぎず、 これらによって本発明の 技術的範囲が限定的に解釈されてはならないものである。 すなわち 、 本 発明はその精神、 またはその主要な特徴から'逸脱するこ となく 、 様々な 形で実施することができる。 産業上の利用可能性
本発明は、 単板式のカラ一撮像素子を用いて得られた色モザィ ク画像 から、 全画素に複数色の輝度情報を補間してカラ一画像を生成するデモ ザィ ク処理を実行する画像処理装置に用いて好適である。

Claims

請 求 の 範 囲
1 . 異なる複数の色光を光電変換する画 *を有する単板式の撮像素子に よって得られた、 各画素が単色の輝度情報を有した色モザィ ク画像を用 いて、 各画素が複数色の輝度情報を有したカラー画像を生成する画像処 理装置であって、
上記色モザィ ク画像を、 同一の色光の画素値だけを含む複数の色プレ
—ンに分解する色プレーン分解部と、
上記色モザィ ク画像から生成されるカラー画像に対する画像加ェを表 す係数を用いて 、 上記カラー画像の画素位置から、 上記画像加 Xが施さ れた場合のカラ —画像の画素位置に対応する上記色モザィ ク画像上のサ ンプ リ ング座標を算出する座標変換部と、
上記色プレ一ン分解部によ り分解された複数の色プ レーン毎に 、 上 座標変換部によ り算出された上記サンプリ ング座標における画素値を、 上記色プレーン内に含まれる同一の色光の画素値から補間生成するサン プリ ング部と、
上記サンプリ ング部によ り補間生成された各色プ レーンの補間値を合 成するこ とによ り上記カラー画像を生成する色生成部とを備えたこ と を 特徴とする画像処理装置。
2 . 上記座標変換部は、 上記画像加工を表す係数と して 、 上記色モザィ ク画像から生成されるカラ一画像に対する画像変形を表す画像変形係数 を用いて、 上記カラー画像の画素位置から、 上記画像変形が施された場 合のカラ一画像の画素位置に対応する上記色モザィ ク画像上のサンプリ ング座標を算出するこ と を特徴とする請求の範囲第 1 項に記載の画像処 理装置。
3 . 上記色モザイ ク画像は、 上記撮像素子の各画素に対応する位置にそ れぞれ赤色 ( R ) 、 緑色 (G r , G b ) 、 青色 ( B ) の何れかが配され たべィャ配列のカラーフ ィ ノレタおよび上記撮像 子によって得られる画 像であり、
上記色プ レーン分解部は、 互いに角を接する 2種類の緑色画素 G r , G b をそれぞれ異なる色プ レーンに分解するこ と を特徴とする請求の範 囲第 1項に記載の画像処理装置。
4 . 上記色生成部は、 上記 2種類の緑色画素 G Γ , G b毎に分解された 各色プレー ン上でそれぞれ上記サンプ リ ング部によ り生成された画素値 を加算して上記カラ一画像を生成するこ と を特徵とする請求の範囲第 3 項に記載の画像処理装置。
5 . 上記 2種類の緑色画素 G r , G b毎に分解された各色プ レーン上で それぞれ上記サンプリ ング部によ り生成された画素値の差分を算出し、 該画素値の差分に基づいて偽色の有無を判定する偽色判定部を更に備 えたこ とを特徴とする請求の範囲第 3項に記載の画像処理装置。
6 . 上記サンプリ ング部は、 上記サンプリ ング座標における画素値を、 その周辺に存在する上記色プレーン上の画素値から補間によって生成す ることを特徴とする請求の範囲第 1項に記載の画像処理装置。
. 上記画像変形は、 上記カラー画像の拡大、 縮小、 回転 、 歪曲修正の ラ ち少なく と も 1 つを含む処理であるこ と を特徴とする請求の範囲第 2 項に記載の画像処理装置。
8 . 上記座標変換部は、 上記画像加 Xを表す係数と して 、 上記色プレ一 ン分解部によ り分解された複数の色プレー ン毎に、 色プレ一ンによって 値の異なる係数を用いて、 上記色モザィ ク画像から生成されるカラ一画 像の画素位置から、 色収差補正が施された場合の力ラー画像の画素位置 対応する上記色モザィ ク画像上のサンプリ ング座標であつて上記色プ レ —ン毎に異なるサンプ リ ング座標をそれぞれ算出する と を特徴とす る請求の範囲第 1項に記載の画像処 3里装置。
9 . 上記座標変換部は、 上記画像加ェを表す係数と して、 上記色プレー ン分解部によ り分解された複数の色プレーン毎に 、 色プレーンによって 値の異なる色収差係数と、 上記カラ一画像に対する画像変形を表す画像 変形係数と を用いて、 上記カラ一画像の画素位置から、 上記色収差補正 および上記画像変形が施された場合の力ラー画像の画素位置に対応する 上記色モザィ ク画像上のサンプリ ング座標であって上記色プレーン毎に 異なるサンプリ ング座標を算出するこ と を特徴とする請求の範囲第 1 項 に記載の画像処理装置。
1 0 . 上記座標変換部は、 上記画像加ェを表す係数と して、 上記色プレ ーン分解部によ り分解された複数の色プレーン毎に、 色プ レーンによつ て値の異なる色収差係数と、 上記撮像素子に被写体の撮像光を導く ため の結像光学系に設定されている焦点距離および被写体距離に応じて定ま る歪曲収差係数とを用いて、 上記力ラ一画像の画素位置から、 上記色収 差補正および歪曲収差補正が施された場合のカラー画像の画素位置に対 応する上記色モザィ ク画像上のサンプ V ング座標であって上記色プレ ― ン毎に異なるサンプリ ング座標を算出するこ とを特徴とする請求の範囲 第 1項に記載の画像処理装置。
1 1 . 上記結像光学系に設定されている焦点距離および被写体距離に応 じて、 上記座標変換部で使用する上記 曲収差係数を設定する係数設定 部を更に備えたこ とを特徴とする 5H求の範囲第 8項に記載の画像処理装 置。
1 2 . 異なる複数の色光を光電変換する画素を有する単板式の撮像素子 によって得られた、 各画素が単色の輝度情報を有した色モザィ ク画像を 用いて、 各画素が複数色の輝度情報を有したカラ一画像を生成する画像 処理方法であって、 上記色モザィ ク画像を、 同一の色光の画素値だけを含む複数の色プ レ ーンに分解する色プレーン分解ステ ップと 、
上記色モザィ ク画像から生成されるカ ラー画像に対する画像加ェを表 す係数を用いて 、 上記力 ラー画像の画素位置から、 上記画像加 Xが施さ れた場合の力 ラ一画像の画素位置に対応する上記色モザィ ク画像上のサ ンプリ ング座標を算出する座標変換ステ ップと 、
上記色プレ一ン分解ステ ップで分解された複数の色プレ一ン毎に、 上 記座標変換ステ シプで算出された上記サンプリ ング座標における画素値 を、 上記色プ レー ン内に含まれる同一の色光の画素値から補間生成する サンプリ ングステ ップと、
上記サンプリ ングステ ップで補間生成された各色プレーンの補間値を 合成する こ と によ り 上記カ ラー画像を生成する色生成ステ ップと を有す るこ と を特徴とする画像処理方法。
1 3 . 上記座標変換ステ ップでは、 上記画像加ェを表す係数と して 、 上 記色モザイ ク画像から生成される力 ラー画像に対する画像変形を表す画 像変形係数を用いて、 上記カ ラー画像の画素位置から、 上記画像変形が 施された場合のカ ラ一画像の画素位置に対応する上記色モザィ ク画像上 のサンプリ ング座標を算出する こ と を特徴とする請求の範囲第 1 2項に 記載の画像処理方法。
1 4 . 上記座標変換ステ ッ プでは、 上記画像加ェを表す係数と して 、 上 記色プレーン分解ステ ップで分解された複数の色プレーン毎に、 色プレ ーンによ って値の異なる係数を用いて、 上記色モザィ ク画像から生成さ れるカ ラー画像の画素位置から、 色収差補正が施された場合のカ ラー画 像の画素位置に対応する上記色モザィ ク画像上のサンプリ ング座標であ つて上記色プレーン毎に異なるサンプリ ング座標をそれぞれ算出する こ と を特徴とする請求の範囲第 1 2項に記載の画像処理方法。
1 5 . 異なる複数の色光を光冪変換する画素を有する撮像素子によって 得られた、 各画素が単色の輝度情報を有した色モザィ ク画像を 、 同一の 色光の画素値だけを含む複数の色プレー ンに分解する色プレ ―ン分解手 段、
上記色モザィ ク画像から生成されるカラ 画像に対する 像加ェを表 す係数を用いて 、 上 力ラ一画像の画素位置から、 上記画像加工が施さ れた場合の力ラ一画像の画 置に対応する上記色モザィ ク画像上のサ ンプリ ング座標を算出する座標変換手段、
上記色プレーン分解手段によ り分解された複数の色プレーン毎に、 上 記座標変換手段によ り算出された上記サンプリ ング座標における画素値 を、 上記色プレーン內に含まれる同一の色光の画素値から補間生成する サンプリ ング手段、 および
上記サンプリ ング手段によ り補間生成された各色プレーンの補間値を 合成するこ とによ り 、 1 つの画素値がそれぞれ複数色の輝度情報を有し たカラー画像を生成する色生成手段、
と してコ ン ピュータを機能させるためのコ ン ピュータ読み取り 可能なプ ログラム。
1 6 . 上記座標変換手段は、 上記画像加工を表す係数と して、 上記色モ ザイ ク画像から生成されるカラー画像に対する画像変形を表す画像変形 係数を用いて、 上記カラー画像の画素位置から、 上記画像変形が施され た場合のカラー画像の画素位置に対応する上記色モザイ ク画像上のサン プリ ング座標を算出するこ とを特徴とする請求の範囲第 1 5項に記寧の コンピュータ読み取り可能なプログラム。
1 7 . 上記座標変換手段は、 上記画像加工を表す係数と して、 上記色プ レーン分解手段によ り分解された複数の色プ レー ン毎に、 色プ レーンに よって値の異なる係数を用いて、 上記色モザィ ク画像から生成される力 ラ―画像の画素位置力、 ら、 色収差捕正が施された場合の力ラ一画像の面 素位置に対応する上記色モザィ ク画像上のサンプリ ング座標であって上 記色プレーン母に異なるサンプリ ング座標をそれぞれ算出するこ とを特 徵とする請求の範囲第 1 5項に記載のコンピ タ ¾πみ取り 可能なプロ グラム ο
1 8 異なる複数の色光を光電変換する画素を有する単板式の撮像素子 を備えた撮像装置であって、
上記単板式の撮像素子と、
上記撮像素子によ り 光電変換された色モザイ ク画像のアナロ グ信号を デジタル信号に変換する A Z D変換器と、
上記 A Z D変換器よ り上記デジタル信号と して出力される 各画 が 単色の輝度情報を有した色モザィ ク画像を用レ、て 各画素が複数色の 情報を有したカラ一画像を生成する請求の 囲第 1 項に記載の画像処 理装置とを備えたこ とを特徴とする撮像装置。
1 9 . 上記請求の範囲第 1項に記載の画像処理装置に代えて請求の 囲 第 2項に記載の画像処理装置を備えたこ とを特徴とする g青求の範囲第 1
8項に記載の撮像装置。
2 0 . 上記請求の範囲第 1 項に記載の画像処理装置に代えて a*求の 囲 第 8に記載の画像処理装置を備えたこ と を特徵とする請求の in囲第 1 8 項に記載の撮像装置。
2 1 . 上記撮像素子に被写体の撮像光を導く ための結像光学系を更に備
X.
上記請求の範囲第 1項に記載の画像処理装置に代えて請求の範囲第 9 項または第 1 0項に記載の画像処理装置を備えたこ と を特徵とする 水 第 1 8項に記載の撮像装置。
上記撮像素子に被写体の撮像光を導く ための結像光学系を更に備 え、
上記請求の範囲第 1項に記載の画像処理装置に代えて 求の範囲第 1 1項に記載の画像処理装置を備えたこ と を特徴とする請
Figure imgf000053_0001
の範囲第 1 8 項に記载の撮像装置。
2 3 . 所定の画像を撮影するこ とによって上記色収差係数を算出する収 差係数算出部を備え、
上記収差係数算出部によって算出された上記色収差係数 上記座標変 換部に格納するこ とを特徴とする請求の範囲第 2 0項また 第 2 1 項に 記載の撮像装置。
2 4 . 所定の画像を撮影するこ とによって上記色収差係数を算出する収 差係数算出部を備え、
上記収差係数算出部によって算出された上記色収差係数を上記係数設 定部に格納するこ とを特徴とする請求の範囲第 2 2項に s i載の撮像装置
2 5 . 外部に接続された外部機器に撮影画像を出力する画像出力部と、 上記外部機器から上記色収差係数を入力する係数入力手段とを備え、 上記係数入力手段によ り入力した上記色収差係数を上記座標変換部に 格納するこ とを特徴とする請求の範囲第 2 0項または第 2 1 項に記載の 撮像装置。
2 6 . 外部に接続された外部機器に撮影画像を出力する画像出力部と、 上記外部機器から上記色収差係数を入力する係数入力手段とを備え、 上記係数入力手段によ り入力した上記色収差係数を上記係数設定部に 格納するこ とを特徵とする請求の範囲第 2 2項に記載の撮像装置。
PCT/JP2008/061806 2007-06-25 2008-06-24 画像処理装置、画像処理方法、プログラムおよび撮像装置 WO2009001956A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
EP08765855A EP2175656A1 (en) 2007-06-25 2008-06-24 Image processing device, image processing method, program, and imaging device
KR1020107001700A KR101356286B1 (ko) 2007-06-25 2008-06-24 화상 처리 장치, 화상 처리 방법, 프로그램 및 촬상 장치
US12/666,511 US8494260B2 (en) 2007-06-25 2008-06-24 Image processing device, image processing method, program, and imaging device
US13/929,423 US8934713B2 (en) 2007-06-25 2013-06-27 Image processing device, image processing method, program, and imaging device

Applications Claiming Priority (6)

Application Number Priority Date Filing Date Title
JP2007-166106 2007-06-25
JP2007-166190 2007-06-25
JP2007166106A JP5036421B2 (ja) 2007-06-25 2007-06-25 画像処理装置、画像処理方法、プログラムおよび撮像装置
JP2007166190 2007-06-25
JP2007-335065 2007-12-26
JP2007335065A JP5036524B2 (ja) 2007-06-25 2007-12-26 画像処理装置、画像処理方法、プログラムおよび撮像装置

Related Child Applications (2)

Application Number Title Priority Date Filing Date
US12/666,511 A-371-Of-International US8494260B2 (en) 2007-06-25 2008-06-24 Image processing device, image processing method, program, and imaging device
US13/929,423 Continuation US8934713B2 (en) 2007-06-25 2013-06-27 Image processing device, image processing method, program, and imaging device

Publications (1)

Publication Number Publication Date
WO2009001956A1 true WO2009001956A1 (ja) 2008-12-31

Family

ID=41819271

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2008/061806 WO2009001956A1 (ja) 2007-06-25 2008-06-24 画像処理装置、画像処理方法、プログラムおよび撮像装置

Country Status (5)

Country Link
US (1) US8494260B2 (ja)
EP (1) EP2175656A1 (ja)
KR (1) KR101356286B1 (ja)
TW (1) TW200901747A (ja)
WO (1) WO2009001956A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114051125A (zh) * 2021-11-11 2022-02-15 锐芯微电子股份有限公司 图像处理方法及系统

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5757099B2 (ja) 2010-02-15 2015-07-29 株式会社ニコン 焦点調節装置、及び焦点調節プログラム
US8692906B2 (en) * 2010-07-13 2014-04-08 Canon Kabushiki Kaisha Luminance signal generation apparatus, luminance signal generation method, image pickup apparatus, and storage medium
JP5867669B2 (ja) * 2011-02-22 2016-02-24 富士ゼロックス株式会社 色処理装置および色処理プログラム
KR20130057327A (ko) * 2011-11-23 2013-05-31 엘지이노텍 주식회사 스테레오 매칭 시스템의 전처리 장치
EP2800380B1 (en) * 2011-12-28 2018-02-28 Fujifilm Corporation Imaging device and imaging processing method
JP2013192023A (ja) * 2012-03-14 2013-09-26 Sony Corp 色距離測定装置、色距離測定方法およびプログラム
WO2013169384A1 (en) 2012-05-11 2013-11-14 Intel Corporation Systems, methods, and computer program products for compound image demosaicing and warping
KR102512521B1 (ko) 2015-10-12 2023-03-21 삼성전자주식회사 텍스쳐 처리 방법 및 장치
GB2585050A (en) 2019-06-26 2020-12-30 Light Blue Optics Ltd Image demosaicing
US11367166B1 (en) * 2019-09-09 2022-06-21 Apple Inc. Method and device for correcting chromatic aberration using chromatic characteristic values
KR102431479B1 (ko) * 2020-09-18 2022-08-11 그린스펙(주) 이미지의 밝기정보와 색깔정보의 보정기능을 갖춘 내시경기기
KR102462804B1 (ko) * 2021-06-28 2022-11-04 주식회사 아이닉스 베이어 패턴 이미지에 기반한 렌즈 왜곡 보정 장치 및 그 보정 방법
CN115829867A (zh) * 2022-11-25 2023-03-21 广州城建职业学院 基于标定的广角图像畸变和色差处理方法、装置和介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06113309A (ja) 1992-09-28 1994-04-22 Sony Corp ディジタルビデオカメラ
JP2000270294A (ja) 1998-10-01 2000-09-29 Seiko Epson Corp デジタルカメラおよび画像データ処理方法
JP2001186533A (ja) * 1999-12-22 2001-07-06 Olympus Optical Co Ltd 画像処理装置
JP2002199410A (ja) * 2000-12-25 2002-07-12 Canon Inc 撮像装置、画像劣化補正処理方法及び記憶媒体
JP2004241991A (ja) * 2003-02-05 2004-08-26 Minolta Co Ltd 撮像装置、画像処理装置及び画像処理プログラム
JP2004534429A (ja) 2001-03-20 2004-11-11 ヒューレット・パッカード・カンパニー 色不連続性均等化を使用して生データイメージを非対称にデモザイク処理するシステムおよび方法
JP2004357298A (ja) * 2003-05-27 2004-12-16 Hewlett-Packard Development Co Lp 空間的可変デモザイク処理を提供するシステムおよび方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20030059399A (ko) * 2001-12-29 2003-07-10 엘지전자 주식회사 모자이크 영상 생성장치 및 방법과 모자이크 영상 기반동영상 재생방법
US6919892B1 (en) * 2002-08-14 2005-07-19 Avaworks, Incorporated Photo realistic talking head creation system and method
US6928194B2 (en) * 2002-09-19 2005-08-09 M7 Visual Intelligence, Lp System for mosaicing digital ortho-images
JP4479457B2 (ja) * 2004-05-27 2010-06-09 ソニー株式会社 画像処理装置、および画像処理方法、並びにコンピュータ・プログラム
US20070250197A1 (en) * 2006-04-21 2007-10-25 Mosaic4U Ltd. Method and apparatus for generating mosaic from an image
US8724895B2 (en) * 2007-07-23 2014-05-13 Nvidia Corporation Techniques for reducing color artifacts in digital images

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06113309A (ja) 1992-09-28 1994-04-22 Sony Corp ディジタルビデオカメラ
JP2000270294A (ja) 1998-10-01 2000-09-29 Seiko Epson Corp デジタルカメラおよび画像データ処理方法
JP2001186533A (ja) * 1999-12-22 2001-07-06 Olympus Optical Co Ltd 画像処理装置
JP2002199410A (ja) * 2000-12-25 2002-07-12 Canon Inc 撮像装置、画像劣化補正処理方法及び記憶媒体
JP2004534429A (ja) 2001-03-20 2004-11-11 ヒューレット・パッカード・カンパニー 色不連続性均等化を使用して生データイメージを非対称にデモザイク処理するシステムおよび方法
JP2004241991A (ja) * 2003-02-05 2004-08-26 Minolta Co Ltd 撮像装置、画像処理装置及び画像処理プログラム
JP2004357298A (ja) * 2003-05-27 2004-12-16 Hewlett-Packard Development Co Lp 空間的可変デモザイク処理を提供するシステムおよび方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114051125A (zh) * 2021-11-11 2022-02-15 锐芯微电子股份有限公司 图像处理方法及系统
CN114051125B (zh) * 2021-11-11 2023-07-04 锐芯微电子股份有限公司 图像处理方法及系统

Also Published As

Publication number Publication date
EP2175656A1 (en) 2010-04-14
KR101356286B1 (ko) 2014-01-28
US8494260B2 (en) 2013-07-23
TW200901747A (en) 2009-01-01
US20100208106A1 (en) 2010-08-19
KR20100051056A (ko) 2010-05-14

Similar Documents

Publication Publication Date Title
WO2009001956A1 (ja) 画像処理装置、画像処理方法、プログラムおよび撮像装置
JP5036421B2 (ja) 画像処理装置、画像処理方法、プログラムおよび撮像装置
JP5222472B2 (ja) 画像処理装置、画像復元方法およびプログラム
JP5036482B2 (ja) 画像処理装置及び画像処理方法、画像処理プログラム
JP2931520B2 (ja) 単板式カラービデオカメラの色分離回路
JP5151075B2 (ja) 画像処理装置及び画像処理方法、撮像装置、並びにコンピュータ・プログラム
US8885067B2 (en) Multocular image pickup apparatus and multocular image pickup method
US7834907B2 (en) Image-taking apparatus and image processing method
JP5438579B2 (ja) 画像処理装置及びその制御方法
JP2004241991A (ja) 撮像装置、画像処理装置及び画像処理プログラム
JP2002084547A (ja) 画像データサイズ変換処理装置、電子スチルカメラ、および画像データサイズ変換処理用記録媒体
WO2013008517A1 (ja) 撮像装置及び画像生成方法
US8982248B2 (en) Image processing apparatus, imaging apparatus, image processing method, and program
KR20180118790A (ko) 화상 처리 장치, 화상 처리 방법, 기록 매체, 프로그램 및 촬상 장치
US8934713B2 (en) Image processing device, image processing method, program, and imaging device
US20110141321A1 (en) Method and apparatus for transforming a lens-distorted image to a perspective image in bayer space
US20080012964A1 (en) Image processing apparatus, image restoration method and program
JP5268321B2 (ja) 画像処理装置及び画像処理方法、画像処理プログラム
JP4635769B2 (ja) 画像処理装置、画像処理方法および撮像装置
JP5036524B2 (ja) 画像処理装置、画像処理方法、プログラムおよび撮像装置
JP4687454B2 (ja) 画像処理装置および撮像装置
JP4962293B2 (ja) 画像処理装置、画像処理方法、プログラム
US9635260B2 (en) Image processing apparatus, image pickup apparatus, and image processing method
JP2003348382A (ja) 輪郭強調回路
JP5454820B2 (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: 08765855

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 376/DELNP/2010

Country of ref document: IN

ENP Entry into the national phase

Ref document number: 20107001700

Country of ref document: KR

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 2008765855

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 12666511

Country of ref document: US