CN111447426A - 一种影像色彩校正方法以及装置 - Google Patents

一种影像色彩校正方法以及装置 Download PDF

Info

Publication number
CN111447426A
CN111447426A CN202010402951.7A CN202010402951A CN111447426A CN 111447426 A CN111447426 A CN 111447426A CN 202010402951 A CN202010402951 A CN 202010402951A CN 111447426 A CN111447426 A CN 111447426A
Authority
CN
China
Prior art keywords
window
pixel
processed
image
calculating
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.)
Granted
Application number
CN202010402951.7A
Other languages
English (en)
Other versions
CN111447426B (zh
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.)
China Toprs Technology Co ltd
Original Assignee
China Toprs Technology Co ltd
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 China Toprs Technology Co ltd filed Critical China Toprs Technology Co ltd
Priority to CN202010402951.7A priority Critical patent/CN111447426B/zh
Publication of CN111447426A publication Critical patent/CN111447426A/zh
Application granted granted Critical
Publication of CN111447426B publication Critical patent/CN111447426B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N9/00Details of colour television systems
    • H04N9/64Circuits for processing colour signals
    • H04N9/68Circuits for processing colour signals for controlling the amplitude of colour signals, e.g. automatic chroma control circuits
    • H04N9/69Circuits for processing colour signals for controlling the amplitude of colour signals, e.g. automatic chroma control circuits for modifying the colour signals by gamma correction

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Image Processing (AREA)

Abstract

本发明提供了影像色彩校正方法以及装置,该方法包括:根据预设的窗口尺寸将所获取的正射纠正影像的四至范围划分为多个窗口;计算每个窗口的局部均值和目标值;将所有待处理影像位于某窗口内所有像素的均值作为该窗口的局部均值;计算参考影像的像素均值、所有待处理影像的像素均值,两均值之差作为像素偏移量;根据像素偏移量,以及每个窗口对应的局部均值,计算每个窗口的目标值;最后利用每个窗口的局部均值和目标值对待处理影像进行逐像素Gamma校正。本发明实施例在不增加算法复杂度的前提下,改进目标值的获取方法,提高色彩校正的质量,改善处理结果中的偏色现象。

Description

一种影像色彩校正方法以及装置
技术领域
本发明涉及影像处理技术领域,具体而言,涉及一种影像色彩校正方法以及装置。
背景技术
数字正射纠正影像是以数字高程模型为辅助数据,对数字航空影像,经逐个像元进行投影差改正获取的影像数据。
航空影像是通过航空摄影获得的。由于航空摄影过程中,随着摄影时间的推移和拍摄位置的不断变化,导致外界的光照环境会不断的变化,因此造成不同的航空影像内部,或者航空影像之间存在亮度不均匀和色彩差异的问题。当对航空影像进行正射纠正时,亮度不均匀和色彩差异的问题会出现在正射纠正影像中。为了解决上述问题,需要对由航空影像生成的正射纠正影像进行色彩一致性处理。
当前对数字正射影像进行色彩一致性处理一般基于Xiuguang Zhou在《MultipleAuto-Adapting Color Balancing for Large Number of Images》中提出了一种基于常平面的Gamma色彩校正方法。该方法首先按照指定的窗口尺寸划分所有待处理影像的范围,统计窗口的像素均值进行线性内插获取每个像素对应的局部均值;然后将参考影像的均值或所有待处理影像的均值作为目标值;最后利用局部均值和目标值构成Gamma校正函数,完成输入影像的逐像素调整。
这种方法虽然在一定程度上解决了影像的色彩差异问题,但由于目标值是利用参考影像的均值来获得,忽略了待处理影像的原始信息,导致最终的处理结果存在一定的偏色现象。
发明内容
有鉴于此,本发明实施例的目的在于提供一种影像色彩校正方法以及装置,充分利用待处理影像的信息,在不增加算法复杂度的前提下,改进目标值的获取方法,提高色彩校正的质量,改善处理结果中的偏色现象。
第一方面,本发明实施例提供了一种影像色彩校正方法,包括:
根据预设的窗口尺寸将所获取的所述四至范围划分为多个窗口;每个所述窗口覆盖至少一个所述待处理影像;
计算每个所述窗口所覆盖的所述待处理影像中有效像素的局部均值;
根据参考影像中的像素值,以及所述待处理影像中像素值,计算所述参考影像以及所述待处理影像之间的像素点偏移量;
根据所述像素点偏移量,以及每个窗口对应的局部均值,计算与每个所述窗口对应的目标值;
使用每个所述窗口对应的局部均值和每个所述窗口对应的所述目标值对所述待处理影像进行逐像素Gamma校正。
第二方面,本发明实施例还提供一种影像色彩校正装置,该装置包括:
区域划分模块,用于根据预设的窗口尺寸将所获取的所述四至范围划分为多个窗口;每个所述窗口覆盖至少一个所述待处理影像;
局部均值计算模块,用于计算每个所述窗口所覆盖的所述待处理影像中像素局部均值;
偏移量计算模块,用于根据参考影像的像素值,以及所述待处理影像的像素值,计算所述参考影像以及所述待处理影像之间的像素偏移量;
目标值计算模块,用于根据所述像素偏移量,以及每个窗口对应的局部均值,计算与每个所述窗口对应的目标值;
校正模块,用于使用每个所述窗口对应的局部均值和每个所述窗口对应的所述目标值对所述待处理影像进行逐像素Gamma校正。
第三方面,本申请实施例还提供一种计算机设备,包括:处理器、存储器和总线,所述存储器存储有所述处理器可执行的机器可读指令,当计算机设备运行时,所述处理器与所述存储器之间通过总线通信,所述机器可读指令被所述处理器执行时执行如第一方面所述的方法的步骤。
第四方面,本申请实施例还提供一种计算机可读存储介质,该计算机可读存储介质上存储有计算机程序,该计算机程序被处理器运行时执行如第一方面所述的方法的步骤。
本发明实施例所提供的影像色彩校正方法以及装置,在对局部均值进行计算的时候,是对于覆盖多幅待处理影像的窗口,其局部均值为覆盖的所有影像的像素值均值;在计算目标值的时候,首先要根据所述根据参考影像的像素值,以及所述待处理影像的像素值,计算所述参考影像以及所述待处理影像之间的像素偏移量,然后再根据像素偏移量,以及每个窗口对应的局部均值,计算每个窗口对应的目标值,最后使用每个所述窗口对应的局部均值和每个所述窗口对应的所述目标值对所述待处理影像进行逐像素Gamma校正。这种方法考虑了每个窗口内重叠影像的信息以及窗口周围的信息,实现了局部色差调整的目的;利用参考影像获取偏移量,实现了整体色调的调整,较之现有技术能够提高色彩校正的质量,获得处理结果的偏色现象较少。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1示出了本发明实施例所提供的一种影像色彩校正方法的流程图;
图2示出了本发明实施例所提供的影像色彩校正方法中,窗口与待处理图像的重叠示意图;
图3示出了本发明实施例所提供的影像色彩校正方法中,计算参考影像以及待处理影像之间的像素偏移量的具体方法流程图;
图4示出了本发明实施例所提供的影像色彩校正方法中,对影像进行逐像素Gamma校正的具体方法流程图;
图5示出了本发明实施例所提供的一种影像色彩校正装置的结构示意图
图6示出了本发明实施例所提供的一种计算机设备的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
目前在对数字正射影像进行色彩一致性处理时,一般采用基于常平面的Gamma色彩校正方法,该方法在获得目标值时,直接将参考影像的像素均值或者所有待处理影像的像素均值作为目标值,造成目标值的获取方式过于简单,忽略了待处理影像的原始信息,导致对数字正射影像进行一致性的色彩校正处理所生成的影像仍然存在偏色问题。基于此,本申请提供的一种影像色彩校正方法以及装置,可以充分利用待处理影像的原始信息,在不增加算法复杂度的前提下,提高色彩校正的质量。
为便于对本实施例进行理解,首先对本发明实施例所公开的一种影像处理方法进行详细介绍,该方法除了能够对数字正射纠正影像进行色彩校正之外,还能够对其他的包含地理坐标的影像进行色彩校正。由于数字正射纠正影像可能是多波段影像,因此在进行处理时,要在不同的波段下分别进行处理。本申请实施例以在一个波段下对数字正射纠正影像进行处理为例,对本申请所提供的影像色彩校正方法加以说明。
参见图1所示,本申请实施例所提供的影像色彩校正方法,包括:
S101:获取多个待处理影像,并获取待处理的正射纠正影像的四至范围。
在具体实现的时候,待处理影像即为要进行色彩校正的正射纠正影像。航空摄影时,根据航摄设计方案,通过航摄设备多次曝光,可以获取航向旁向都存在重叠的多幅航空影像。这些航空影像镶嵌前要进行正射纠正,得到的正射纠正影像即为待处理影像。由于待处理的正射纠正影像包含地理坐标,因此四至范围是指所有待处理影像所到达的东南西北四个方向的边界。
S102:根据预设的窗口尺寸将所获取的所述四至范围划分为多个窗口;每个所述窗口覆盖至少一个所述待处理影像。
在具体实现的时候,由于太阳入射角、大气以及光照条件等因素的存在,会造成航空影像影像内部色差和待处理影像之间的局部色差,正射纠正后的影像中会遗留这些问题为了消除上述的色差问题,需要首先按照预设窗口尺寸将待处理影像的四至范围划分为多个窗口,然后以窗口为单位统计色彩校正信息,最后将待处理影像进行逐像素校正,从而改善局部色差,使整体色调趋于一致。
此处,需要注意的是,由于相邻的航空影像之间会存在重叠的部分,因此窗口对应的地理范围可能涉及一幅或多幅影像。。如图2所示,不同粗细的黑虚线为待处理影像A、待处理影像B和待处理影像C,三者相互重叠;覆盖三幅影像的红色虚线框为四至范围,根据预设窗口尺寸对范围进行划分,得到红色网格;W1、W2、W3分别代表三种类型的窗口,W1涉及一幅影像A,W2涉及两幅影像A和C,W3涉及三幅影像A、B、C。
S103:将每个所述窗口所覆盖的所述待处理影像的像素均值作为局部均值。
在具体实现的时候,首先要获取每个窗口中所覆盖的所有待处理影像的所有像素值。然后统计像素点的数量。最后根据窗口中所有像素值和像素点数量,计算该窗口的像素值均值。
其中,每一个窗口的局部均值M(m,n)满足下述公式(1):
M(m,n)=1/cm,nij vin(i,j)|(i,j)∈ωm,n (1)
其中cm,n为所有待处理影像位于所述窗口ωm,n的像素总数,vin(i,j)为所述窗口ωm,n中(i,j)处的像素值。
如果待处理影像有多个波段,则每个波段中像素值均值均满足上述公式(1)。
S104:根据参考影像的像素值,以及所述待处理影像的像素值,计算所述参考影像以及所述待处理影像之间的像素偏移量。
在具体实现的时候,参考影像是色彩均匀、与待处理影像覆盖的地物类似的影像,该参考影像可以从待处理影像中选择,也可以从已有资料中获取。
参考影像以及待处理影像之间的像素偏移量,是指参考影像每个波段的像素值统计量与待处理影像对应波段的像素值统计量之间的差异。
具体地,参见图3所示,本申请实施例提供一种计算参考影像以及待处理影像之间的像素偏移量的具体方法,包括:
S301:计算参考影像的像素点均值,并计算所有所述待处理影像的像素点均值。
在具体实现的时候,参考影像的像素点均值P(m,n)满足下述公式(2):
P(m,n)=1/sm,nij win(i,j)|(i,j)∈wm,n (2)
其中sm,n为参考影像wm,n的像素点总数,win(i,j)为所述参考影像wm,n中(i,j)处的像素值。
待处理影像的像素点均值Q(m,n)满足下述公式(3)
Q(m,n)=1/fm,nij uin(i,j)|(i,j)∈um,n (3)
其中fm,n为所有待处理影像um,n的中的有效像素总数,uin(i,j)为所述待处理影像um,n中(i,j)处的像素值。
S302:将所述参考影像的像素点均值以及所有所述待处理影像的像素点均值之间的差作为所述像素偏移量。
像素偏移量diff满足公式(4):
diff=P(m,n)-Q(m,n) (4)
S105:根据所述像素偏移量,以及每个窗口的局部均值,计算与每个所述窗口对应的目标值。
在具体实现的时候,当计算出像素偏移量后,每个窗口对应的目标值为该窗口的所述局部均值与该像素偏移量之和。
具体地,每个窗口对应的目标值T(m,n)满足下述公式(5):
T(m,n)=M(m,n)+diff (5)
S106:使用每个所述窗口的所述局部均值和每个所述窗口对应的所述目标值对所述待处理影像进行逐像素Gamma校正。
在具体实现的时候,逐像素Gamma校正的处理流程参见图4所示,包括:
S401:根据每个所述窗口的局部均值和所述目标值,计算每个窗口对应的Gamma值;
S402:对每个像素所在窗口和邻近窗口对应的Gamma值进行线性内插,获取每个像素点对应的Gamma值;
S403:根据每个像素点对应的Gamma值构建所述Gamma校正函数。
在具体实现的时候,根据每个所述窗口对应的局部均值和每个窗口对应的所述目标值构建的Gamma校正函数如下述公式(6)所示:
vout(i,j)=α×vin(i,j)γ(i,j) (6)
其中,vout(i,j)为像素点(i,j)处的输出像素值;vin(i,j)为像素点(i,j)处的输入像素值,一般被标准化为0.0~1.0;γ(i,j)为像素点(i,j)处的Gamma值;α为常数,用于控制输出像素值的范围。
其中,像素点(i,j)处的Gamma值γ(i,j),根据该像素点所在窗口以及周围窗口线性内插获得,且γ(i,j)满足下述公式(7):
Figure BDA0002490191670000081
g(m,n)=log(T(m,n))/log(M(m,n)) (8)
其中,BI为线性内插算子,g(m,n)为窗口(m,n)处的Gamma值,符号±由像素(i,j)与窗口ωm,n中心的位置确定,T(m,n)为窗口ωm,n的目标值,M(m,n)为窗口ωm,n的局部均值。
本申请实施例所提供的影像色彩校正方法,在对局部均值进行计算的时候,是对于覆盖多幅待处理影像的窗口,其局部均值为覆盖的所有影像的像素点均值;在计算目标值的时候,首先要根据所述参考影像的像素值,以及所述待处理影像的像素值,计算所述参考影像以及所述待处理影像之间的像素偏移量,然后再根据像素偏移量,以及每个窗口对应的局部均值,计算每个窗口对应的目标值,最后使用每个所述窗口对应的局部均值和每个所述窗口对应的目标值对所述待处理影像进行逐像素Gamma校正。这种方法充分利用了待处理影像的重叠信息,实现了局部色差调整的目的;利用参考影像获取偏移量,实现了整体色调的调整,较之现有技术能够提高色彩校正的质量,改善了偏色现象。
本发明又一实施例还提供一种影像色彩校正装置,参见图5所示,本发明实施例所提供的影像色彩校正装置包括:
区域划分模块501,用于根据预设的窗口尺寸将所获取的所述四至范围划分为多个窗口;每个所述窗口覆盖至少一个所述待处理影像;
局部均值计算模块502,用于计算每个所述窗口所覆盖的所述待处理影像的局部均值;
偏移量计算模块503,用于根据参考影像像素值,以及所述待处理影像的像素值,计算所述参考影像以及所述待处理影像之间的像素偏移量;
目标值计算模块504,用于根据所述像素偏移量,以及每个窗口对应的局部均值,计算与每个所述窗口对应的目标值;
校正模块505,用于使用每个所述窗口对应的局部均值和每个所述窗口对应的所述目标值对所述待处理影像中的每一个所述窗口进行Gamma校正。
本申请实施例所提供的影像色彩校正装置,在对局部均值进行计算的时候,窗口可能覆盖多幅待处理影像,其局部均值为覆盖的所有影像的像素均值;在计算目标值的时候,首先要根据所述根据参考影像的像素值,以及所述待处理影像的像素值,计算所述参考影像以及所述待处理影像之间的像素偏移量,然后再根据像素偏移量,以及每个窗口对应的局部均值,计算每个窗口对应的目标值,然后使用每个所述窗口对应的局部均值和每个所述窗口对应的所述目标值对所述待处理影像进行逐像素Gamma校正。这种方法充分利用了待处理影像的重叠信息,实现了局部色差调整的目的;利用参考影像获取偏移量,实现了整体色调的调整,较之现有技术能够提高色彩校正的质量,改善了偏色现象。
在一可选实施例中,还包括:获取模块,用于在根据预设的窗口尺寸将所获取的所述四至范围划分为多个窗口之前,获取多个待处理影像,并获取多个待处理影像的四至范围。
在一可选实施例中,局部均值计算模块502,具体用于获取每个窗口中所覆盖的所有待处理影像的像素值;
根据每个窗口中所覆盖的所有待处理影像的像素值,以及所述窗口内的像素点数量,计算所述窗口对应的局部均值。
在一可选实施例中,所述窗口的局部均值M(m,n)满足下述公式(1):
(1)M(m,n)=1/cm,nijvin(i,j)|(i,j)∈ωm,n
其中cm,n为所有待处理影像位于所述窗口ωm,n的像素总数,vin(i,j)为所述窗口ωm,n中(i,j)处的像素值。
在一可选实施例中,偏移量计算模块503具体用于:计算参考影像的像素值均值,并计算所有所述待处理影像的像素点均值;
将所述参考影像的像素值均值以及所有所述待处理影像的像素点均值之间的差作为所述像素偏移量。
在一可选实施例中,目标值计算模块504具体用于:将所述像素偏移量与每个所述窗口对应的所述局部均值的和,作为每个所述窗口对应的目标值。
在一可选实施例中,校正模块505,具体用于:使用每个所述窗口对应的局部均值和所述目标值,构建所述窗口的Gamma值;
g(m,n)满足下述公式(2):
(2)g(m,n)=log(T(m,n))/log(M(m,n));
其中T(m,n)为所述窗口对应的目标值;M(m,n)为所述窗口中像素点的像素值均值。
在一可选实施例中,校正模块505,具体用于:根据每个所述窗口和邻近窗口的Gamma值,计算每个像素对应的Gamma值;
对每个所述窗口对应的Gamma值进行线性内插,获取每个像素点对应的Gamma值,满足下述公式(3):
(3)
Figure BDA0002490191670000111
其中,BI为线性内插算子,g(m,n)为窗口(m,n)的Gamma值,符号±由像素(i,j)与窗口ωm,n中心的位置确定,T(m,n)为窗口ωm,n的目标值,M(m,n)为窗口ωm,n的局部均值。
在一可选实施例中,利用每个像素的Gamma值,构建Gamma校正函数。
根据所述Gamma校正函数计算所述窗口内每个所述像素点的输出像素值;
将所述输出像素值作为该像素点在所述待处理影像中的Gamma校正结果。
Gamma校正函数满足下述公式(6):
(6)vout(i,j)=α×vin(i,j)γ(i,j)
其中,vout(i,j)为像素点(i,j)处的输出像素值;vin(i,j)为像素点(i,j)处的输入像素值;γ(i,j)为像素点(i,j)处的Gamma值;α为常数。
本申请实施例还提供了一种计算机设备60,如图6所示,为本申请实施例提供的计算机设备60结构示意图,包括:处理器61、存储器62、和总线63。所述存储器62存储有所述处理器61可执行的机器可读指令,当计算机设备60运行时,所述处理器61与所述存储器62之间通过总线63通信,所述机器可读指令被所述处理器61执行时执行如下处理:
根据预设的窗口尺寸将所获取的所述四至范围划分为多个窗口;每个所述窗口覆盖至少一个所述待处理影像;
计算每个所述窗口所覆盖的所述待处理影像中像素的局部均值;
根据参考影像中的像素值,以及所述待处理影像中的像素值,计算所述参考影像以及所述待处理影像之间的像素偏移量;
根据所述像素偏移量,以及每个窗口对应的局部均值,计算与每个所述窗口对应的目标值;
使用每个所述窗口对应的局部均值和每个所述窗口对应的所述目标值对所述待处理影像进行逐像素伽马Gamma校正。
一种可能的实施方式中,处理器61执行的指令中,根据预设的窗口尺寸将所获取的所述四至范围划分为多个窗口之前,还包括:
获取多个待处理正射纠正影像,并获取所有待处理影像的四至范围。
一种可能的实施方式中,处理器61执行的指令中,所述计算每个所述窗口所覆盖的所述待处理影像中像素的局部均值,具体包括:
获取每个窗口中所覆盖的所有待处理影像的所有像素值;
根据每个窗口中所覆盖的所有待处理影像的所有像素值,以及像素个数,计算所述窗口对应的局部均值;
将所述窗口内像素均值作为所述窗口的局部均值。
一种可能的实施方式中,处理器61执行的指令中,所述窗口的局部均值M(m,n)满足下述公式(1):(1)M(m,n)=1/cm,nij vin(i,j)|(i,j)∈ωm,n
其中,cm,n为所有待处理影像位于所述窗口ωm,n的像素总数,vin(i,j)为所述窗口ωm,n中(i,j)处的像素值。
一种可能的实施方式中,处理器61执行的指令中,所述根据参考影像中像素值,以及所述待处理影像中像素值,计算所述参考影像以及所述待处理影像之间的像素偏移量,具体包括:
计算参考影像中在所有像素均值,并计算所有所述待处理影像中所有像素均值;
将所述参考影像中像素均值以及所有所述待处理影像中像素均值之间的差作为所述像素偏移量。
一种可能的实施方式中,处理器61执行的指令中,根据所述像素偏移量,以及每个窗口对应的局部均值,计算与每个所述窗口对应的目标值,具体包括:
将所述像素偏移量与每个所述窗口对应的所述局部均值的和,作为每个所述窗口对应的目标值。
一种可能的实施方式中,处理器61执行的指令中,所述使用每个所述窗口对应的局部均值和每个所述窗口对应的所述目标值对所述待处理影像进行逐像素Gamma校正,具体包括:
使用每个所述窗口对应的局部均值和所述目标值,计算所述窗口的Gamma值;
根据像素所在窗口和邻近窗口的Gamma值,采用线性内插方法获取每个像素的Gamma值;
根据Gamma校正函数计算每个所述像素点的输出像素值;
将所述输出像素值作为所述待处理影像对应该像素点的Gamma校正结果。
一种可能的实施方式中,处理器61执行的指令中,所述窗口的Gamma校正函数满足下述公式(6):
(6)vout(i,j)=α×vin(i,j)γ(i,j)
其中,vout(i,j)为像素点(i,j)处的输出像素值;vin(i,j)为像素点(i,j)处的输入像素值;γ(i,j)为像素点(i,j)处的Gamma值;α为常数;
且γ(i,j)满足下述公式(7):
(7)
Figure BDA0002490191670000141
其中,BI为线性内插算子,g(m,n)为窗口(m,n)的Gamma值,符号±由像素(i,j)与窗口ωm,n中心的位置确定,T(m,n)为窗口ωm,n的目标值,M(m,n)为窗口ωm,n的局部均值;
g(m,n)满足下述公式(8):
(8)g(m,n)=log(T(m,n))/log(M(m,n));
且T(m,n)为所述窗口对应的目标值;M(m,n)为所述窗口的局部均值。
本申请实施例还提供一种计算机可读存储介质,该计算机可读存储介质上存储有计算机程序,该计算机程序被处理器运行时执行上述方法实施例中所述的影像色彩校正方法的步骤。
本申请实施例所提供的影像色彩校正方法的计算机程序产品,包括存储了程序代码的计算机可读存储介质,所述程序代码包括的指令可用于执行上述方法实施例中所述的影像色彩校正方法的步骤,具体可参见上述方法实施例,在此不再赘述。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的系统和装置的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。

Claims (11)

1.一种影像色彩校正方法,其特征在于,包括:
根据预设的窗口尺寸将所获取的所述四至范围划分为多个窗口;每个所述窗口覆盖至少一个所述待处理影像;
计算每个所述窗口所覆盖的所述待处理影像中像素的局部均值;
根据参考影像中的像素值,以及所述待处理影像中的像素值,计算所述参考影像以及所述待处理影像之间的像素偏移量;
根据所述像素偏移量,以及每个窗口对应的局部均值,计算与每个所述窗口对应的目标值;
使用每个所述窗口对应的局部均值和每个所述窗口对应的所述目标值对所述待处理影像进行逐像素伽马Gamma校正。
2.根据权利要求1所述的方法,其特征在于,根据预设的窗口尺寸将所获取的所述四至范围划分为多个窗口之前,还包括:
获取多个待处理正射纠正影像,并获取所有待处理影像的四至范围。
3.根据权利要求2所述的方法,其特征在于,所述计算每个所述窗口所覆盖的所述待处理影像中像素的局部均值,具体包括:
获取每个窗口中所覆盖的所有待处理影像的所有像素值;
根据每个窗口中所覆盖的所有待处理影像的所有像素值,以及像素个数,计算所述窗口对应的局部均值;
将所述窗口内像素均值作为所述窗口的局部均值。
4.根据权利要求3所述的方法,其特征在于,所述窗口的局部均值M(m,n)满足下述公式(1):(1)M(m,n)=1/cm,nijvin(i,j)|(i,j)∈ωm,n
其中,cm,n为所有待处理影像位于所述窗口ωm,n的像素总数,vin(i,j)为所述窗口ωm,n中(i,j)处的像素值。
5.根据权利要求1所述的方法,其特征在于,所述根据参考影像中像素值,以及所述待处理影像中像素值,计算所述参考影像以及所述待处理影像之间的像素偏移量,具体包括:
计算参考影像中在所有像素均值,并计算所有所述待处理影像中所有像素均值;
将所述参考影像中像素均值以及所有所述待处理影像中像素均值之间的差作为所述像素偏移量。
6.根据权利要求5所述的方法,其特征在于,根据所述像素偏移量,以及每个窗口对应的局部均值,计算与每个所述窗口对应的目标值,具体包括:
将所述像素偏移量与每个所述窗口对应的所述局部均值的和,作为每个所述窗口对应的目标值。
7.根据权利要求1所述的方法,其特征在于,所述使用每个所述窗口对应的局部均值和每个所述窗口对应的所述目标值对所述待处理影像进行逐像素Gamma校正,具体包括:
使用每个所述窗口对应的局部均值和所述目标值,计算所述窗口的Gamma值;
根据像素所在窗口和邻近窗口的Gamma值,采用线性内插方法获取每个像素的Gamma值;
根据Gamma校正函数计算每个所述像素点的输出像素值;
将所述输出像素值作为所述待处理影像对应该像素点的Gamma校正结果。
8.根据权利要求7所述的方法,其特征在于,所述窗口的Gamma校正函数满足下述公式(6):
(6)vout(i,j)=α×vin(i,j)γ(i,j)
其中,vout(i,j)为像素点(i,j)处的输出像素值;vin(i,j)为像素点(i,j)处的输入像素值;γ(i,j)为像素点(i,j)处的Gamma值;α为常数;
且γ(i,j)满足下述公式(7):
(7)
Figure FDA0002490191660000031
其中,BI为线性内插算子,g(m,n)为窗口(m,n)的Gamma值,符号±由像素(i,j)与窗口ωm,n中心的位置确定,T(m,n)为窗口ωm,n的目标值,M(m,n)为窗口ωm,n的局部均值;
g(m,n)满足下述公式(8):
(8)g(m,n)=log(T(m,n))/log(M(m,n));
且T(m,n)为所述窗口对应的目标值;M(m,n)为所述窗口的局部均值。
9.一种影像色彩校正装置,其特征在于,该装置包括:
区域划分模块,用于根据预设的窗口尺寸将所获取的所述四至范围划分为多个窗口;每个所述窗口覆盖至少一个所述待处理影像;
局部均值计算模块,用于计算每个所述窗口所覆盖的所述待处理影像中像素的局部均值;
偏移量计算模块,用于根据参考影像中的像素值,以及所述待处理影像中像素值,计算所述参考影像以及所述待处理影像之间的像素偏移量;
目标值计算模块,用于根据所述像素偏移量,以及每个窗口对应的局部均值,计算与每个所述窗口对应的目标值;
校正模块,用于使用每个所述窗口对应的局部均值和每个所述窗口对应的所述目标值对所述待处理影像进行逐像素Gamma校正。
10.一种计算机设备,其特征在于,包括:处理器、存储器和总线,所述存储器存储有所述处理器可执行的机器可读指令,当计算机设备运行时,所述处理器与所述存储器之间通过总线通信,所述机器可读指令被所述处理器执行时执行如权利要求1至8任一所述的方法的步骤。
11.一种计算机可读存储介质,其特征在于,该计算机可读存储介质上存储有计算机程序,该计算机程序被处理器运行时执行如权利要求1至5任意一项所述的方法的步骤。
CN202010402951.7A 2020-05-13 2020-05-13 一种影像色彩校正方法以及装置 Active CN111447426B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010402951.7A CN111447426B (zh) 2020-05-13 2020-05-13 一种影像色彩校正方法以及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010402951.7A CN111447426B (zh) 2020-05-13 2020-05-13 一种影像色彩校正方法以及装置

Publications (2)

Publication Number Publication Date
CN111447426A true CN111447426A (zh) 2020-07-24
CN111447426B CN111447426B (zh) 2021-12-31

Family

ID=71653736

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010402951.7A Active CN111447426B (zh) 2020-05-13 2020-05-13 一种影像色彩校正方法以及装置

Country Status (1)

Country Link
CN (1) CN111447426B (zh)

Citations (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5373376A (en) * 1992-10-29 1994-12-13 Linotype-Hell Ag Method for processing color values for reproduction scale modification and shapness correction in reproduction of color originals
US5546084A (en) * 1992-07-17 1996-08-13 Trw Inc. Synthetic aperture radar clutter reduction system
EP1292113A2 (en) * 2001-08-23 2003-03-12 Eastman Kodak Company Tone scale adjustment
US20030053690A1 (en) * 2001-07-06 2003-03-20 Jasc Software, Inc. Automatic contrast enhancement
US20070116375A1 (en) * 2004-04-12 2007-05-24 Nikon Corporation Image processing device having color shift correcting function, image processing program and electronic camera
US20090154793A1 (en) * 2007-12-17 2009-06-18 Electronics And Telecommunications Research Institute Digital photogrammetric method and apparatus using intergrated modeling of different types of sensors
US20110157435A1 (en) * 2009-12-31 2011-06-30 Eastman Kodak Company Generating column offset corrections for image sensors
CN102129669A (zh) * 2011-02-24 2011-07-20 武汉大学 一种航空遥感影像的最小二乘区域网匀色方法
US20120002113A1 (en) * 2010-06-30 2012-01-05 Sony Corporation Image processing device, image processing method, and program
US20130002908A1 (en) * 2011-06-29 2013-01-03 Tessera Technologies Ireland Limited Axial chromatic aberration correction
CN103280188A (zh) * 2013-06-14 2013-09-04 电子科技大学 Oled器件老化补偿系统及方法
KR20140048528A (ko) * 2012-10-16 2014-04-24 현대모비스 주식회사 임계값을 이용한 이미지 센서의 픽셀 간 오프셋 보정 장치 및 그 보정 방법
US8891981B2 (en) * 2010-03-31 2014-11-18 Canon Kabushiki Kaisha Image forming apparatus capable of correcting image information
CN104574337A (zh) * 2015-01-21 2015-04-29 山东科技大学 基于双边伽马校正和多尺度图像融合的图像增强方法
US20150237360A1 (en) * 2014-02-18 2015-08-20 Electronics And Telecommunications Research Institute Apparatus and method for fast sample adaptive offset filtering based on convolution method
CN105005767A (zh) * 2015-07-01 2015-10-28 东北林业大学 一种基于微波遥感影像森林类型识别方法
US20160125660A1 (en) * 2006-05-02 2016-05-05 Digitalglobe, Inc. Automated registration of three-dimensional vectors to three-dimensional linear features in remotely-sensed data
CN106875364A (zh) * 2017-02-23 2017-06-20 深圳飞马机器人科技有限公司 一种真正射影像生成方法
US20170371046A1 (en) * 2014-12-15 2017-12-28 Koninklijke Philips N.V. Pixel based dead time correction
US20180089843A1 (en) * 2016-09-26 2018-03-29 Digitalglobe, Inc. Techniques for image co-registration
CN109035152A (zh) * 2018-05-23 2018-12-18 电子科技大学 一种合成孔径雷达图像非局部均值滤波方法
CN109063711A (zh) * 2018-07-06 2018-12-21 航天星图科技(北京)有限公司 一种基于llts框架的卫星影像正射纠正算法
CN109272465A (zh) * 2018-09-15 2019-01-25 武汉智觉空间信息技术有限公司 一种航空影像色彩一致性处理算法
CN109493332A (zh) * 2018-11-06 2019-03-19 北京京航计算通讯研究所 基于envi的优化矢量选取roi遥感影像预处理系统
CN111127318A (zh) * 2019-12-05 2020-05-08 南京莱斯电子设备有限公司 一种机场环境下的全景图像拼接方法

Patent Citations (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5546084A (en) * 1992-07-17 1996-08-13 Trw Inc. Synthetic aperture radar clutter reduction system
US5373376A (en) * 1992-10-29 1994-12-13 Linotype-Hell Ag Method for processing color values for reproduction scale modification and shapness correction in reproduction of color originals
US20030053690A1 (en) * 2001-07-06 2003-03-20 Jasc Software, Inc. Automatic contrast enhancement
US6826310B2 (en) * 2001-07-06 2004-11-30 Jasc Software, Inc. Automatic contrast enhancement
EP1292113A2 (en) * 2001-08-23 2003-03-12 Eastman Kodak Company Tone scale adjustment
US20070116375A1 (en) * 2004-04-12 2007-05-24 Nikon Corporation Image processing device having color shift correcting function, image processing program and electronic camera
US20160125660A1 (en) * 2006-05-02 2016-05-05 Digitalglobe, Inc. Automated registration of three-dimensional vectors to three-dimensional linear features in remotely-sensed data
US20090154793A1 (en) * 2007-12-17 2009-06-18 Electronics And Telecommunications Research Institute Digital photogrammetric method and apparatus using intergrated modeling of different types of sensors
US20110157435A1 (en) * 2009-12-31 2011-06-30 Eastman Kodak Company Generating column offset corrections for image sensors
US8891981B2 (en) * 2010-03-31 2014-11-18 Canon Kabushiki Kaisha Image forming apparatus capable of correcting image information
US20120002113A1 (en) * 2010-06-30 2012-01-05 Sony Corporation Image processing device, image processing method, and program
CN102129669A (zh) * 2011-02-24 2011-07-20 武汉大学 一种航空遥感影像的最小二乘区域网匀色方法
US20130002908A1 (en) * 2011-06-29 2013-01-03 Tessera Technologies Ireland Limited Axial chromatic aberration correction
KR20140048528A (ko) * 2012-10-16 2014-04-24 현대모비스 주식회사 임계값을 이용한 이미지 센서의 픽셀 간 오프셋 보정 장치 및 그 보정 방법
CN103280188A (zh) * 2013-06-14 2013-09-04 电子科技大学 Oled器件老化补偿系统及方法
US20150237360A1 (en) * 2014-02-18 2015-08-20 Electronics And Telecommunications Research Institute Apparatus and method for fast sample adaptive offset filtering based on convolution method
US20170371046A1 (en) * 2014-12-15 2017-12-28 Koninklijke Philips N.V. Pixel based dead time correction
CN104574337A (zh) * 2015-01-21 2015-04-29 山东科技大学 基于双边伽马校正和多尺度图像融合的图像增强方法
CN105005767A (zh) * 2015-07-01 2015-10-28 东北林业大学 一种基于微波遥感影像森林类型识别方法
US20180089843A1 (en) * 2016-09-26 2018-03-29 Digitalglobe, Inc. Techniques for image co-registration
CN106875364A (zh) * 2017-02-23 2017-06-20 深圳飞马机器人科技有限公司 一种真正射影像生成方法
CN109035152A (zh) * 2018-05-23 2018-12-18 电子科技大学 一种合成孔径雷达图像非局部均值滤波方法
CN109063711A (zh) * 2018-07-06 2018-12-21 航天星图科技(北京)有限公司 一种基于llts框架的卫星影像正射纠正算法
CN109272465A (zh) * 2018-09-15 2019-01-25 武汉智觉空间信息技术有限公司 一种航空影像色彩一致性处理算法
CN109493332A (zh) * 2018-11-06 2019-03-19 北京京航计算通讯研究所 基于envi的优化矢量选取roi遥感影像预处理系统
CN111127318A (zh) * 2019-12-05 2020-05-08 南京莱斯电子设备有限公司 一种机场环境下的全景图像拼接方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
H. SINGH等: "A novel gamma correction approach using optimally clipped sub-equalization for dark image enhancement", 《2016 IEEE INTERNATIONAL CONFERENCE ON DIGITAL SIGNAL PROCESSING (DSP)》 *
JIAYUAN LI等: "Optimal Illumination and Color Consistency for Optical Remote-Sensing Image Mosaicking", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 *
SHUBHI KANSAL等: "Adaptive gamma correction for contrast enhancement of remote sensing images", 《MULTIMEDIA TOOLS AND APPLICATIONS (2019)》 *
XIUGUANG ZHOU: "Multiple Auto-Adapting Color Balancing for Large Number of Images", 《THE INTERNATIONAL ARCHIVES OF THE PHOTOGRAMMETRY, REMOTE SENSING AND SPATIAL INFORMATION SCIENCES》 *
汪荣贵等: "基于Retinex理论的JPEG压缩方法研究 ", 《中国科学技术大学学报》 *
王晓宁等: "多山地区合成孔径雷达数据预处理方法 ", 《云南地理环境研究》 *
王颖洁等: "基于光学影像匹配技术的地震形变场提取研究 ", 《矿山测量》 *

Also Published As

Publication number Publication date
CN111447426B (zh) 2021-12-31

Similar Documents

Publication Publication Date Title
US8374428B2 (en) Color balancing for partially overlapping images
CN107590840B (zh) 基于网格划分的颜色阴影校正方法及其校正系统
EP3100449B1 (en) Method for conversion of a saturated image into a non-saturated image
US10223775B2 (en) Array camera image combination with feature-based ghost removal
US8284271B2 (en) Chroma noise reduction for cameras
US10521887B2 (en) Image processing device and image processing method
US20120093432A1 (en) Image processing device, image processing method and storage medium storing image processing program
JP2010278708A (ja) 画像処理装置及び方法、及びコンピュータプログラム
EP3602400A1 (en) System and methods for correcting overlapping digital images of a panorama
KR20110082903A (ko) 항공사진의 정사보정 및 최적화된 정사영상 제작 방법
CN113592739B (zh) 一种镜头阴影校正的方法、装置及存储介质
CN111598777A (zh) 天空云图的处理方法、计算机设备和可读存储介质
CN108307098A (zh) 鱼眼相机阴影校正参数确定方法、校正方法及装置、存储介质、鱼眼相机
CN102129669A (zh) 一种航空遥感影像的最小二乘区域网匀色方法
CN113496474A (zh) 一种图像处理方法、装置、环视系统、汽车和存储介质
CN114866754A (zh) 自动白平衡方法、装置及计算机可读存储介质和电子设备
CN110175967B (zh) 图像去雾处理方法、系统、计算机设备和存储介质
CN111447426B (zh) 一种影像色彩校正方法以及装置
JP2020091910A (ja) 画像処理装置、画像処理方法及びプログラム
CN116385370A (zh) 鱼眼图像处理方法、装置、电子设备和存储介质
CN116228563A (zh) 一种基于暗通道的图像去雾校正方法及装置
JP4765667B2 (ja) 画像処理プログラムおよび画像処理装置
CN117897723A (zh) 用于学习局部图像增强的色调曲线的系统和方法
CN110458768B (zh) 图像处理方法、计算机设备和可读存储介质
UA115360C2 (uk) Спосіб обробки цифрового зображення для зниження рівня шуму

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant