CN101667112A - 基于cordic算法的图像色调调整方法 - Google Patents

基于cordic算法的图像色调调整方法 Download PDF

Info

Publication number
CN101667112A
CN101667112A CN200910192168A CN200910192168A CN101667112A CN 101667112 A CN101667112 A CN 101667112A CN 200910192168 A CN200910192168 A CN 200910192168A CN 200910192168 A CN200910192168 A CN 200910192168A CN 101667112 A CN101667112 A CN 101667112A
Authority
CN
China
Prior art keywords
tone
delta
tan
gamma
angle
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN200910192168A
Other languages
English (en)
Inventor
谭洪舟
徐永键
张爱华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN200910192168A priority Critical patent/CN101667112A/zh
Publication of CN101667112A publication Critical patent/CN101667112A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Processing Of Color Television Signals (AREA)

Abstract

本发明公开了一种利用cordic算法对图像色调进行调整的方法,它包括1)对图像进行初始化,获取原始图像的红色色差Cr和蓝色色差Cb,以及色调角度θ,
Figure 200910192168.6_AB_0
;2)利用cordic算法计算调整角度θ与旋转因子δi之间的关系,从而确定图像数据各位的δi的取值,即
Figure 200910192168.6_AB_1
,其中γ=-θ;3)根据图像数据各位的δi的取值进行移位操作:首先根据图像调整的色调角度计算得到δi在每一位上的符号;然后从最低位开始对图像色调信息数据各位进行移位操作,对应于各位上的符号选择,若为正号,则保持不变,若为负号,则将本位上的数据取反;最后最高位上移位后操作输出的结果为最终的色调调整结果数据。本发明简化了乘法计算,减小计算复杂度,同时提高了计算精度。

Description

基于cordic算法的图像色调调整方法
技术领域
本发明涉及一种图像色调调整的方法,具体来说,这是一种利用cordic算法对图像色调进行调整的技术。
背景技术
色调控制一般按照函数
Figure A20091019216800041
来进行处理,此处[Cb Cr]与[Cb′Cr′]为处理前后的色差信号,θ表示所要的色调角度。θ的取值范围为(-180°,180°),典型的8位色调调节范围为(-30°,30°)。色调调整的计算过程中需要用到sin和cos值的计算,无法通过硬件来实现精确的计算,业界常用的方法是使用查表法来由色度调整角度θ直接查表得到对应的sin和cos值。但这种算法精度有限,而且会消耗另外的存储空间,影响最终芯片面积。在对饱和度进行调节的过程中,需要比较准确的调节效果,而这里由于精度有限也会引入误差,对调节结果产生影响。
发明内容
针对以上的不足,本发明提供了一种计算精度更好,色调调整更加精确的方法。
本发明提出的基于cordic算法的图像色调调整方法,包括:
1)对图像进行初始化,获取原始图像的红色色差Cr和蓝色色差Cb,以及色调角度θ;
2)利用cordic算法计算调整角度θ与旋转因子δi之间的关系,从而确定图像数据各位的δi的取值。
2.1)令色调角度调整变量γ=-θ,则
Cb′=Cbcosγ-Crsinγ    (1)
Cr′=Crcosγ+Cbsinγ
将式(1)转换为矩阵形式,得到
Cb ′ Cr ′ = cos γ 1 - tan γ tan γ 1 Cb Cr
2.2)对式(2)进行cordic算法变换得到
Cb ( i + 1 ) Cr ( i + 1 ) = K i 1 - δ i 2 - i δ i 2 - i 1 Cb ( i ) Cr ( i ) - - - ( 3 )
其中,迭代系数
Figure A20091019216800053
δi为旋转因子,δi的值为1或-1,若旋转的角度为正,则δi为1,旋转的角度为负,则δi为-1,i为1到n-1之间的自然数,n为色调值的位数。
在迭代计算的过程中,我们可以忽略Ki,然后在迭代计算完成后,将结果乘以如下值:
K ( n ) = Π i = 0 n - 1 K i = Π i = 0 n - 1 1 1 + 2 - 2 i - - - ( 4 )
此系数可以预先计算出来,作为常数定义在模块中,来进一步减小计算复杂度:
K = lim n → ∞ K ( n ) ≈ 0.6072529350088812561694
2.3)由以上式子可得到
γi=tan-1(2-i)    (5)
对应的,可以得到要进行控制的色调角度和各个移位系数δi之间的关系为:
γ = Σ i = 0 n - 1 δ i γ i = Σ i = 0 n - 1 δ i tan - 1 ( 2 - i ) - - - ( 6 )
由于δi的值为1或-1,则上式变成:
γ=±γ0±γ1±γ2±γ3±γ4±γ5±γ6±γ7±......
                                                                  。
  =±tan-1(2-0)±tan-1(2-1)±tan-1(2-2)±tan-1(2-3)±......
3)根据图像数据各位的δi的取值进行移位操作,得到调整后的色差分量Cr′和Cb′。所述图像数据各位进行移位操作包括:
3.1)根据具体图像要的色调角度计算得到δi在每一位上的符号;
3.2)从最低位开始对图像色调信息数据各位进行移位操作:对应于各位上的符号选择,若为正号,则保持不变,若为负号,则将本位上的数据取反;
3.3)最高位上移位后操作输出的结果为最终的色调调整结果数据。
本发明的有益效果:一是简化了乘法计算,减小了计算复杂度,从而减少了最终电路资源;二是用移位操作代替传统的查找表来计算三角函数,节省了存储空间,从而减小了最终电路的面积;三是提高了计算精度,使用本方法比传统的查找表计算正弦值精度更高;四是具有较高灵活性,可以通过增加色度信号的数据位数来提高色调控制的精度。
附图说明
图1为对图像色调进行cordic算法迭代运算原理示意图。
具体实施方式
下面结合附图对本发明进行进一步阐述。
本发明的引入在图像旋转操作中常用的cordic(Coordinate RotationDigital Computer)算法来简化色调控制的计算。
首先对图像进行初始化,获取原始图像的红色色差Cr和蓝色色差Cb,以及色调角度θ,得到
Cb′=Cbcosθ+Crsinθ
Cr′=Crcosθ-Cbsinθ
对于上述公式,令γ=-θ,则原式变为:
Cb′=Cbcosγ-Crsinγ
Cr′=Crcosγ+Cbsinγ
将上式转换为矩阵的形式,且将cos值提取出来,得:
Cb ′ C r ′ = cos γ 1 - tan γ tan γ 1 Cb Cr
然后利用cordic算法计算调整角度θ与旋转因子δi之间的关系,从而确定图像数据各位的δi的取值:
使用cordic算法,上式变为:
Cb ( i + 1 ) Cr ( i + 1 ) = K i 1 - δ i 2 - i δ i 2 - i 1 Cb ( i ) Cr ( i )
其中
Figure A20091019216800073
δi的值为1或-1,用以决定旋转的方向,若旋转的角度为正,则δi为1,旋转的角度为负则δi为-1。在迭代计算的过程中,我们可以忽略Ki,然后在迭代计算完成后,将结果乘以如下值:
K ( n ) = Π i = 0 n - 1 K i = Π i = 0 n - 1 1 1 + 2 - 2 i
此系数可以预先计算出来,作为常数定义在模块中,来进一步减小计算复杂度:
K = lim n → ∞ K ( n ) ≈ 0.6072529350088812561694
由上面的计算过程中可得:
γi=tan-1(2-i)
对应的,可以得到要进行控制的色调角度和各个移位系数δi之间的关系为:
γ = Σ i = 0 n - 1 δ i γ i = Σ i = 0 n - 1 δ i tan - 1 ( 2 - i )
由于δi的值为1或-1,则上式即:
γ=±γ0±γ1±γ2±γ3±γ4±γ5±γ6±γ7±......
  =±tan-1(2-0)±tan-1(2-1)±tan-1(2-2)±tan-1(2-3)±......
最后根据图像数据各位的δi的取值进行移位操作,得到调整后的色差分量Cr′和Cb′:
由上可见,采用cordic算法将色调控制中的正弦和余弦函数计算用移位操作来代替,且简化了乘法操作,更适合硬件实现也简化了计算。但cordic算法需要进行n次迭代以达到n位精度,而且吞吐量低,需要更多的时钟周期延迟才能得到结果。但对于硬件的高速主频,即使几十个周期的延时也不到1ms,故不影响视频观看。嵌入cordic算法进行色调控制,可以根据需要的色调角度计算各个数据位上的移位符号。采用移位方法显然不能很精确的表示(-180°,180°)范围内的每个角度,但对于色调控制操作已经完全能满足调节的要求。
使用cordic算法来进行色调控制的操作,用移位操作来代替正余弦函数的计算和矩阵乘法计算。如图1所示。如上面原理所述:
(1)先根据具体图像要的色调角度计算得到δi在每一位上的符号,写入寄存器。
(2)然后从色度信息数据的最低位开始,根据寄存器中对应于本位上的符号选择,若为正号,则保持不变,若为负号,则将本位上的数据取反。
(3)然后将本位上的结果数据向右移一位,加在临时结果数据上。
(4)重复步骤(2)和(3),最高位上移位后输出的结果即为最终的色调调整结果数据。

Claims (3)

1、一种基于cordic算法的图像色调调整方法,它包括:
1)对图像进行初始化,获取原始图像的红色色度Cr和蓝色色差Cb,以及色调角度θ;
2)利用cordic算法计算调整角度θ与旋转因子δi之间的关系,从而确定图像数据各位的δi的取值;
3)根据图像数据各位的δi的取值进行移位操作,得到调整后的色度分量Cr′和Cb′。
2、根据权利要求1所述的基于cordic算法的图像色调调整方法,其特征在于,所述步骤2)包括:
2.1)令色调角度调整变量γ=-θ,则
Cb′=Cbcosγ-Crsinγ           (1)
Cr′=Crcosγ+Cbsinγ
将式(1)转换为矩阵形式,得到
Cb ′ Cr ′ = cos γ 1 - tan γ tan γ 1 Cb Cr - - - ( 2 )
2.2)对式(2)进行cordic算法变换得到
Cb ( i + 1 ) Cr ( i + 1 ) = K i 1 - δ i 2 - i δ i 2 - i 1 Cb ( i ) Cr ( i ) - - - ( 3 )
其中系数
Figure A2009101921680002C3
δi为旋转因子,δi的值为1或-1,若旋转的角度为正,则δi为1,旋转的角度为负,则δi为-1,i为1到n-1之间的自然数,n为色调值的位数。
在迭代计算的过程中,我们可以忽略Ki,然后在迭代计算完成后,将结果乘以如下值:
K ( n ) = Π i = 0 n - 1 K i = Π i = 0 n - 1 1 1 + 2 - 2 i - - - ( 4 )
此系数可以预先计算出来,作为常数定义在模块中,来进一步减小计算复杂度:
K = lim n → ∞ K ( n ) ≈ 0.6072529350088812561694
2.3)由以上式子可得到
γi=tan-1(2-i)          (5)
对应的,可以得到要进行控制的色调角度和各个移位系数δi之间的关系为:
γ = Σ i = 0 n - 1 δ i γ i = Σ i = 0 n - 1 δ i tan - 1 ( 2 - i ) - - - ( 6 )
由于δi的值为1或-1,则上式变成:
γ=±γ0±γ1±γ2±γ3±γ4±γ5±γ6±γ7±......=±tan-1(2-0)±tan-1(2-1)±tan-1(2-2)±tan-1(2-3)±......。
3、根据权利要求3所述的基于cordic算法的图像色调调整方法,其特征在于,所述图像数据各位进行移位操作包括:
3.1)根据图像调整色调角度计算得到δi在每一位上的符号;
3.2)从最低位开始对图像色调信息数据各位进行移位操作:对应于各位上的符号选择,若为正号,则保持不变,若为负号,则将本位上的数据取反;
3.3)最高位上移位后操作输出的结果为最终的色调调整结果数据。
CN200910192168A 2009-09-09 2009-09-09 基于cordic算法的图像色调调整方法 Pending CN101667112A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200910192168A CN101667112A (zh) 2009-09-09 2009-09-09 基于cordic算法的图像色调调整方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200910192168A CN101667112A (zh) 2009-09-09 2009-09-09 基于cordic算法的图像色调调整方法

Publications (1)

Publication Number Publication Date
CN101667112A true CN101667112A (zh) 2010-03-10

Family

ID=41803739

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200910192168A Pending CN101667112A (zh) 2009-09-09 2009-09-09 基于cordic算法的图像色调调整方法

Country Status (1)

Country Link
CN (1) CN101667112A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103986480A (zh) * 2012-12-10 2014-08-13 英特尔移动通信有限责任公司 用于迭代计算数值的方法和装置
CN111369470A (zh) * 2020-03-10 2020-07-03 昇显微电子(苏州)有限公司 图像区域色调调整方法、装置、存储介质及设备
CN112862691A (zh) * 2021-03-17 2021-05-28 苏州灵目电子科技有限公司 一种双鱼眼图像拼接全景图像接缝融合方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103986480A (zh) * 2012-12-10 2014-08-13 英特尔移动通信有限责任公司 用于迭代计算数值的方法和装置
US9354844B2 (en) 2012-12-10 2016-05-31 Intel Deutschland Gmbh Method and apparatus for iteratively calculating a value
CN103986480B (zh) * 2012-12-10 2016-08-24 英特尔移动通信有限责任公司 用于迭代计算数值的方法和装置
CN111369470A (zh) * 2020-03-10 2020-07-03 昇显微电子(苏州)有限公司 图像区域色调调整方法、装置、存储介质及设备
CN111369470B (zh) * 2020-03-10 2024-05-31 昇显微电子(苏州)股份有限公司 图像区域色调调整方法、装置、存储介质及设备
CN112862691A (zh) * 2021-03-17 2021-05-28 苏州灵目电子科技有限公司 一种双鱼眼图像拼接全景图像接缝融合方法

Similar Documents

Publication Publication Date Title
CN108229648B (zh) 匹配存储器中数据位宽的卷积计算方法和装置、设备、介质
CN101667112A (zh) 基于cordic算法的图像色调调整方法
CN106600537B (zh) 一种反距离权重的异向性三维空间插值方法
CN108055041A (zh) 一种数据类型转换电路单元及装置
Rust et al. Low complexity QR-decomposition architecture using the logarithmic number system
CN111078187B (zh) 一种针对单精度浮点数的任意次方根求解方法及其求解器
CN104135457A (zh) 基于改进的cordic算法的数字鉴相方法及装置
CN110110285B (zh) 一种用于FPGA的并行Jacobi计算加速实现方法
CN111061455A (zh) 一种三角函数cordic迭代运算协处理器
CN106814990A (zh) Cordic算法的优化方法
US20130262541A1 (en) Method and circuitry for square root determination
CN105553913A (zh) 一种基于最大化最小乘积距离准则的scma码本搜寻方法
CN116719499A (zh) 一种应用于5g最小二乘定位的自适应伪逆计算方法
CN104657553B (zh) 一种基于快速归一化互相关法的相似性测度硬件加速方法
CN105634439A (zh) 一种异步成形滤波器设计方法
CN107657078B (zh) 基于fpga的超声相控阵浮点聚焦发射实现方法
CN1983103A (zh) 三角函数值确定装置、方法及应用其的通信装置
CN110798420B (zh) 一种内插倍数可变的gmsk调制实现方法
CN107431495A (zh) 数字预失真校正方法及装置
CN103731386A (zh) 一种基于gpp和simd技术的高速调制方法
CN102566963A (zh) 一种fpga进行数据处理的方法
CN107256140A (zh) 实现基于硬件加速的非标准浮点数重建算法的方法和设备
CN110046327B (zh) 一种基于帕德逼近的通信误差函数逼近方法
CN107707178A (zh) 低通滤波器的滤波方法、低通滤波器及伺服驱动器
CN113515259B (zh) 一种适用于浮点格式的复数的近似取模实现电路及方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20100310