WO2011014170A1 - Method for tone mapping an image - Google Patents

Method for tone mapping an image Download PDF

Info

Publication number
WO2011014170A1
WO2011014170A1 PCT/US2009/052226 US2009052226W WO2011014170A1 WO 2011014170 A1 WO2011014170 A1 WO 2011014170A1 US 2009052226 W US2009052226 W US 2009052226W WO 2011014170 A1 WO2011014170 A1 WO 2011014170A1
Authority
WO
WIPO (PCT)
Prior art keywords
bit depth
linear space
high bit
value
values
Prior art date
Application number
PCT/US2009/052226
Other languages
English (en)
French (fr)
Inventor
Niranjan Damera-Venkata
Nelson Liang An Chang
Original Assignee
Hewlett-Packard Development Company, L.P.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hewlett-Packard Development Company, L.P. filed Critical Hewlett-Packard Development Company, L.P.
Priority to US13/258,563 priority Critical patent/US20120014594A1/en
Priority to PCT/US2009/052226 priority patent/WO2011014170A1/en
Priority to KR1020117027076A priority patent/KR20120046103A/ko
Priority to EP09847913A priority patent/EP2411962A4/en
Priority to CN2009801606961A priority patent/CN102473289A/zh
Priority to JP2012522787A priority patent/JP2013500677A/ja
Priority to TW099124715A priority patent/TW201106295A/zh
Publication of WO2011014170A1 publication Critical patent/WO2011014170A1/en

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/90Dynamic range modification of images or parts thereof
    • G06T5/92Dynamic range modification of images or parts thereof based on global image properties
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/40Picture signal circuits
    • H04N1/407Control or modification of tonal gradation or of extreme levels, e.g. background level
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N5/00Details of television systems
    • H04N5/14Picture signal circuitry for video frequency region
    • H04N5/20Circuitry for controlling amplitude response
    • H04N5/202Gamma control
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10016Video; Image sequence
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20004Adaptive image processing
    • G06T2207/20012Locally adaptive
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20172Image enhancement details
    • G06T2207/20208High dynamic range [HDR] image processing

Definitions

  • Many capture device for example scanners or digital cameras, capture images as a two dimensional array of pixeis. Each pixel will have associated intensity values in a predefined color space, for example red, green and blue.
  • the intensity values may be captured using a high bit depth for each color, for example 12 or 16 bits deep.
  • the captured intensity values are typically linearly spaced.
  • the intensity values of each color When saved as a final image, or displayed on a display screen, the intensity values of each color may be converted to a lower bit depth with a non-iinear spacing, for example 8 bits per color.
  • a final image with 8 bits per color (with three colors) may be represented as a 24 bit color image. Mapping the linear high bit depth image (12 or 16 bits per color) into the lower nonlinear bit depth image (8 bits per color) is typically done using a gamma correction tone map.
  • Multi-projector systems often require high-bit depth to prevent contouring in the blend areas (the blends must vary smoothly). This becomes a much more significant issue when correcting black offsets digitally since a discrete digital jump from 0 to I does not allow a representation of continuous values in that range. Also, in a display system the "blends" or s ⁇ bframe values are often computed in linear space with high precision (16-bit) and then gamma corrected to 8 non-iinear bits.
  • FIG. 1 is a two dimensional array of intensity values representing a small part of an image, in an example embodiment of the invention.
  • FIG.2 is a tabic showing the mapping of the intensity values of a linear 4 bit image into the intensity values of a non-linear 2 bit image with a gamma of 2.2.
  • FlG. 3 shows the image from figure 1 after having been mapped into a 2 bit
  • FIG. 4 is a flow chart showing a method for combining gamma correction with dithering in an example embodiment of the invention.
  • FIG. Sa is a table showing the intensity values of the high bit depth image in an example embodiment of the invention.
  • Sb is a table showing the intensity values of the lower bit depth image in non-iinear space and in linear space, in an example embodiment of the invention.
  • FIG.6 is a dither pattern in an example embodiment of the invention.
  • FIG. 7 is a small image, in an example embodiment of the invention.
  • FIG. 8 is a table that lists the results for overlaying the dither pattern in figure 6 onto the small image of figure 7, in an example embodiment of the invention.
  • FIG. 9 is a final image in an example embodiment of the invention.
  • FIG. 10 is a block diagram of a computer system 1000 in an example embodiment of the invention.
  • FIGS, i - 10 and the following description depict specific examples to teach those skilled in the art how to make and use the best mode of the invention. For the purpose of teaching inventive principles, some conventional aspects have been simplified or omitted. Those skilled in the art will appreciate variations from these examples that fall within the scope of the invention. Those skilled in the art will appreciate that the features described below can be combined in various ways to form multiple variations of the invention. As a result, the invention is not limited to the specific examples described below, but only by the claims and their equivalents.
  • mapping an image from a high bit depth linear image into a lower bit depth non-linear image can be done over many different bit depth levels. For example mappings may be done from 16 bits (65,536 levels) to 8 bits (256 levels), from 12 bit to S bite, from S bits to 4 bits, from 4 bits into 2 bits, or the like.
  • each intensity level in the high bit depth image is first normalized to between 0 and 1.
  • each color channel is processed independently. Normalization is done by dividing the original intensity value by the largest possible intensity value for the current bit depth.
  • Normalized Non-linear Value (Normalized Value) ⁇ ( I /gamma) Equation 1
  • the normalized Non-linear intensity value is given by raising die normalized intensity value to one over the gamma value. For a gamma of 2.2, the normalized intensity value would be raised to the power of 1/2.2 or 0.4545. The original intensity value of 50 would yield a normalized mapped value of 0.4812
  • Figure 1 is a two dimensional array of intensity values representing a small part of an image, in an example embodiment of the invention.
  • the image in figure 1 is a 4 bit image with intensity values ranging from 0 - 15.
  • Figure 2 is a table showing the mapping of the intensity values of a linear 4 bit image into the intensity values of a nonlinear 2 bit image with a gamma of 2.2.
  • Figure 3 shows the image from figure 1 after having been mapped into a 2 bit (4 level) space using a 2.2 gamma mapping.
  • Figure 3 may have visible banding between the 3 different levels.
  • a dithering step is combined with the mapping step to produce an image that may show less contouring.
  • Figure 4 is a flow chart showing a method for combining gamma correction with dithering in an example embodiment of the invention. Using the method shown in figure 4, a high bit depth linear image is represented using a smaller number of non-linear levels where the smaller number of non-linear levels are spatially modulated across the final image.
  • each intensity value in the high bit depth linear image is mapped to an intensity value in the non-linear space, hi one example embodiment of the invention, the mapping is done using gamma correction. In other example embodiments of the invention, other mapping algorithms may be used.
  • a left and right interval boundary is calculated for each of the intensity values in nonlinear space. Once the left and right interval boundaries are calculated, they are mapped into linear space.
  • a dither pattern is overlaid onto the pixels of the original image in linear space.
  • the intensity value at each pixel is snapped to one of the two closest left and right interval boundaries in linear space, based on the original linear intensity value, the left and right interval boundary values (in linear space), and the value of the dither screen at that pixel location.
  • the non-linear gamma corrected intensity value for the pixel location is determined.
  • FIG. 7 The first step is to map each intensity value in the high bit depth linear image to an intensity value in the nonlinear space.
  • Equation i is used for mapping from a linear image to a non-linear image when the mapping is done using a gamma correction function.
  • Figure Sa is a table showing the intensity values of the high bit depth image in an example embodiment of the invention.
  • the first column in figure 5a lists the normalized intensity values in 4 bit linear space.
  • the second column in figure Sa lists the normalized intensity values in non-linear space.
  • Each intensity value in column 2 was generated using equation 1 with a 2.2 gamma correction.
  • the gamma corrected value for intensity value 2 (in non-linear space) is generated by first normalizing the 4 bit value, and then raising that normalized value to the power of 1/2.2 resulting in a value of 0.40017 ( (0. 13333 > ⁇ ⁇ 3/2.2) ⁇ 0.40017).
  • the next step is to generate the left and right boundary intervals for each high bit depth intensity value.
  • the left and right boundary intervals represent the two closest lower bit depth non-linear intensity values to the current non-linear intensity value. Equations 2 and 3 are used to calculate the left and right boundary intervals respectively.
  • IntensityVal is the normalized high bit depth intensity value in non- liner space
  • MaxIV is the maximum low bit depth intensity value
  • intergerValue is a function that truncates any fractional value (i.e. it converts a floating point value into an integer value).
  • the first step in equation 1 [integerValue(intensityVal * MaxIV) ] takes the normalized high bit depth intensity value and multiplies it by the maximum quantized low bit depth intensity value. The result is converted from a floating point value into an integer. This converts the normalized high bit depth intensity value into a lower bit depth intensity value.
  • the second step in equation 1 normalizes the lower bit depth value to between zero and one by dividing by the maximum low bit depth intensity value. The calculation for the left boundary interval value in non-linear space for the 4 bit intensity value of 6 is shown below. Left ⁇ ((integerValue(0.65935• 3))/3) Left ⁇ ((integerValueO .97805))/3) Left ⁇ (1/3) Left ⁇ 0.33333
  • FIG. 5b is a table showing the intensity values of the lower bit depth image in non-linear space and in linear space, in an example embodiment of the invention.
  • the first column in figure 5b lists the intensity values of the lower bit depth image in non-linear space.
  • the second column in table 5b lists the intensity values of the lower bit depth image in linear space.
  • a dither pattern is overlaid onto the pixels of the original image in linear space.
  • a dither pattern may be a matrix of threshold intensity values, a single threshold intensity value with a pattern for propagating error to other pixels, a single threshold with a pattern of noise addition, or the like.
  • the dither pattern is shown in figure 6. Any type of dither pattern may be used, including error diffusion or random noise injection. The size of the dither pattern may also be varied.
  • the dither pattern shown in figure 6 is a 4x4 Bayer dither pattern. Before the dither pattern is overlaid onto the intensity values in the original image, the intensity values in the dither pattern are normalized to a value between 0 and 1.
  • the intensity value at each pixel is snapped to one of the two closest left and right interval boundaries in linear space, based on the original linear intensity value, the left and right interval boundary values in linear space, and the value of the dither screen at that pixel location.
  • the correct left or right interval boundary is selected using equations 4 and 5.
  • IntensityN is the original high bit depth linear intensity value for the current pixel normalized to between 0 and 1
  • left and right are the left and right boundary intervals in linear space for the current intensity value
  • Dither is the normalized dither value for the current pixel.
  • CompVal is set to zero when the expression is false and CompVal is set to one when the expression is true.
  • SelectedVal will equal the right value when CompVal is one, and will equal the left value when CompVal is a zero.
  • Figure 7 is a small section of an image, in an example embodiment of the invention.
  • Figure S is a table that lists the results for overlaying the dither pattern in figure 6 onto the small image of figure 7, in an example embodiment of the invention.
  • the first column in figure 8 lists the pixel location in the image.
  • the second column lists the normalized intensity value of the image for each pixel location.
  • the third and fourth columns list the left and right boundary intervals in linear space for each pixel location, respectively.
  • the fifth column lists the normalized dither pattern value for each pixel location.
  • the sixth column lists the calculated CompVal for each pixel location.
  • the last column lists the SelectedVal for each pixel location.
  • Equations 4 and 5 are used to calculate the last two columns in figure 8.
  • the last step is to map the selected value from the linear space to the nonlinear space. This can be done using a lookup table.
  • the lookup table in figure Sb is used tor this example.
  • Figure 9 is the final image from the example above.
  • the image can be saved or stored onto a computer readable medium.
  • a computer readable medium can comprise the following: random access memory, read only memory, hard drives, tapes, optical disk drives, non-volatile ram, video ram, and the like.
  • the image can be used in many ways, for example displayed on one or more displays, transferred to other storage devices, or the like.
  • Computer system has a processor 1002, a memory device 1004, a storage device 1006, a display 1008, and an I/O device 1010.
  • the processor 1002, memory device 1004, storage device 1006, display device 1008 and I/O device 1010 are coupled together with bus 1012.
  • Processor 1002 is configured to execute computer instruction that implement the method describe above.

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Processing (AREA)
  • Facsimile Image Signal Circuits (AREA)
PCT/US2009/052226 2009-07-30 2009-07-30 Method for tone mapping an image WO2011014170A1 (en)

Priority Applications (7)

Application Number Priority Date Filing Date Title
US13/258,563 US20120014594A1 (en) 2009-07-30 2009-07-30 Method for tone mapping an image
PCT/US2009/052226 WO2011014170A1 (en) 2009-07-30 2009-07-30 Method for tone mapping an image
KR1020117027076A KR20120046103A (ko) 2009-07-30 2009-07-30 이미지를 톤 매핑하기 위한 방법
EP09847913A EP2411962A4 (en) 2009-07-30 2009-07-30 METHOD FOR CONVERTING IMAGE TONES
CN2009801606961A CN102473289A (zh) 2009-07-30 2009-07-30 用于色调映射图像的方法
JP2012522787A JP2013500677A (ja) 2009-07-30 2009-07-30 画像をトーンマッピングするための方法
TW099124715A TW201106295A (en) 2009-07-30 2010-07-27 Method for tone mapping an image

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/US2009/052226 WO2011014170A1 (en) 2009-07-30 2009-07-30 Method for tone mapping an image

Publications (1)

Publication Number Publication Date
WO2011014170A1 true WO2011014170A1 (en) 2011-02-03

Family

ID=43529592

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2009/052226 WO2011014170A1 (en) 2009-07-30 2009-07-30 Method for tone mapping an image

Country Status (7)

Country Link
US (1) US20120014594A1 (ja)
EP (1) EP2411962A4 (ja)
JP (1) JP2013500677A (ja)
KR (1) KR20120046103A (ja)
CN (1) CN102473289A (ja)
TW (1) TW201106295A (ja)
WO (1) WO2011014170A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014178989A1 (en) * 2013-04-29 2014-11-06 Dolby Laboratories Licensing Corporation Dithering for chromatically subsampled image formats
US9767542B2 (en) 2013-02-27 2017-09-19 Thomson Licensing Method and device for selecting an image dynamic range conversion operator

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013055472A (ja) * 2011-09-02 2013-03-21 Sony Corp 画像処理装置および方法、並びにプログラム
US9955084B1 (en) 2013-05-23 2018-04-24 Oliver Markus Haynold HDR video camera
GB2520406B (en) * 2013-10-17 2015-11-04 Imagination Tech Ltd Tone mapping
US10277771B1 (en) 2014-08-21 2019-04-30 Oliver Markus Haynold Floating-point camera
US10225485B1 (en) 2014-10-12 2019-03-05 Oliver Markus Haynold Method and apparatus for accelerated tonemapping
CN108241868B (zh) * 2016-12-26 2021-02-02 浙江宇视科技有限公司 图像客观相似度到主观相似度的映射方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5377041A (en) 1993-10-27 1994-12-27 Eastman Kodak Company Method and apparatus employing mean preserving spatial modulation for transforming a digital color image signal
US20040075744A1 (en) * 2002-10-17 2004-04-22 Canon Kabushiki Kaisha Automatic tone mapping for images
US20040105106A1 (en) * 1999-08-24 2004-06-03 Miller Steven O. Reducing quantization errors in imaging systems
WO2009002121A1 (en) * 2007-06-27 2008-12-31 Core Logic, Inc. Non-linear tone mapping apparatus and method and computer readable medium stored thereon computer executable instructions for performing the method

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5963714A (en) * 1996-11-15 1999-10-05 Seiko Epson Corporation Multicolor and mixed-mode halftoning
US7054038B1 (en) * 2000-01-04 2006-05-30 Ecole polytechnique fédérale de Lausanne (EPFL) Method and apparatus for generating digital halftone images by multi color dithering
EP1264473A2 (en) * 2000-02-01 2002-12-11 Pictologic, Inc. Method and apparatus for quantizing a color image through a single dither matrix
US6637851B2 (en) * 2001-03-09 2003-10-28 Agfa-Gevaert Color halftoning for printing with multiple inks
US7079684B2 (en) * 2001-12-05 2006-07-18 Oridus, Inc. Method and apparatus for color quantization of images employing a dynamic color map
JP4208911B2 (ja) * 2006-08-31 2009-01-14 キヤノン株式会社 画像処理装置及びその方法、並びに、コンピュータプログラムおよび記録媒体

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5377041A (en) 1993-10-27 1994-12-27 Eastman Kodak Company Method and apparatus employing mean preserving spatial modulation for transforming a digital color image signal
US20040105106A1 (en) * 1999-08-24 2004-06-03 Miller Steven O. Reducing quantization errors in imaging systems
US20040075744A1 (en) * 2002-10-17 2004-04-22 Canon Kabushiki Kaisha Automatic tone mapping for images
WO2009002121A1 (en) * 2007-06-27 2008-12-31 Core Logic, Inc. Non-linear tone mapping apparatus and method and computer readable medium stored thereon computer executable instructions for performing the method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
DALE A. SCHUMACHER: "Graphics gems II", 1991, ACADEMIC PRESS, article "11.2 A comparison of digital halftoning techniques", pages: 57 - 71
ORCHARD M.: "COLOR QUANTIZATION OF IMAGES", IEEE TRANSACTIONS ON SIGNAL PROCESSING, vol. 39, no. 12, 1 December 1991 (1991-12-01), pages 2677 - 2690, XP000275119, DOI: doi:10.1109/78.107417
See also references of EP2411962A4

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9767542B2 (en) 2013-02-27 2017-09-19 Thomson Licensing Method and device for selecting an image dynamic range conversion operator
WO2014178989A1 (en) * 2013-04-29 2014-11-06 Dolby Laboratories Licensing Corporation Dithering for chromatically subsampled image formats
KR101746491B1 (ko) * 2013-04-29 2017-06-13 돌비 레버러토리즈 라이쎈싱 코오포레이션 채색적으로 서브샘플링된 이미지 포맷들을 위한 디더링
US9762876B2 (en) 2013-04-29 2017-09-12 Dolby Laboratories Licensing Corporation Dithering for chromatically subsampled image formats

Also Published As

Publication number Publication date
KR20120046103A (ko) 2012-05-09
EP2411962A4 (en) 2012-09-19
EP2411962A1 (en) 2012-02-01
US20120014594A1 (en) 2012-01-19
TW201106295A (en) 2011-02-16
JP2013500677A (ja) 2013-01-07
CN102473289A (zh) 2012-05-23

Similar Documents

Publication Publication Date Title
EP2411962A1 (en) Method for tone mapping an image
JP6614859B2 (ja) 表示装置、表示装置の制御方法、画像処理装置、プログラム、及び、記録媒体
US8897559B2 (en) Method, system and apparatus modify pixel color saturation level
KR20200050453A (ko) 번-인 통계 및 번-인 보상
JP6548517B2 (ja) 画像処理装置および画像処理方法
US20160352975A1 (en) Method for conversion of a saturated image into a non-saturated image
JP2004159344A (ja) コントラスト補正装置及びその方法
JP2013520934A5 (ja)
KR20080038414A (ko) 테이블 구성을 위한 시스템, 방법, 및 장치, 그리고 이미지프로세싱 시의 사용
WO2016031006A1 (ja) 表示装置、階調補正マップ生成装置、階調補正マップ生成方法及びプログラム
KR20170040865A (ko) 표시장치와 그 영상 렌더링 방법
US20110090243A1 (en) Apparatus and method for inter-view crosstalk reduction
US9832395B2 (en) Information processing method applied to an electronic device and electronic device having at least two image capturing units that have the same image capturing direction
US10152945B2 (en) Image processing apparatus capable of performing conversion on input image data for wide dynamic range
US7289666B2 (en) Image processing utilizing local color correction and cumulative histograms
CN109448644B (zh) 一种用于校正显示设备的灰阶显示曲线的方法、电子设备及计算机可读存储介质
US20100328504A1 (en) Adaptive method for processing digital images, and an image processing device
US10346711B2 (en) Image correction device, image correction method, and image correction program
US20120321178A1 (en) Method for stitching image in digital image processing apparatus
JP4397623B2 (ja) 階調補正装置
US20120188390A1 (en) Methods And Apparatuses For Out-Of-Gamut Pixel Color Correction
US8774553B1 (en) Advanced adaptive contrast enhancement
KR20070012017A (ko) 디스플레이 기기의 특정 색보정 방법 및 이의 장치
US20080158429A1 (en) Circuit and method of dynamic contrast enhancement
JP2004260835A (ja) 画像処理装置、画像処理方法、画像処理制御プログラムを記録した媒体

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 200980160696.1

Country of ref document: CN

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

Ref document number: 09847913

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 13258563

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 2009847913

Country of ref document: EP

ENP Entry into the national phase

Ref document number: 20117027076

Country of ref document: KR

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 2012522787

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE