技术领域
[0001] 本
发明涉及彩色
图像处理技术领域，尤其涉及一种低照度条件下图像增强方法及系统。
背景技术
[0002] 现实生活中，彩色图像增强在数字图像处理中比灰度图像增强起更加关键作用，因为彩色图像具有三种
感知属性(包括
色调，
饱和度和图像强度)的
色彩表达。由于彩色图像在许多场景中的广泛应用，图像和视频增强已经成为重要的研究领域，近年来它一直是一个研究热点。彩色图像常常会出现的问题之一是，由于照明或某些其他条件(如成像装置限制和天气的变化)，经常造成图像暗部区域模糊不清、图像
对比度低等现象，这样的图像首先给人带来视觉上的不舒适，同时也不利于后续的图像处理工作。尽管图像捕获设备获得极大的改进，但仍然存在各种自然和人工伪影，这导致所捕获图像的
质量较差。因此，针对原始捕获彩色图像的视觉效果和质量的改进是图像预处理必不可少的部分。增强低照度、模糊或不完整图像等的色彩仍然是一个悬而未决的问题。低照度彩色图像有两种来源：夜间和
不平衡的光线。本发明主要增强不平衡光线下的低照度彩色图像。
[0003] 低照度图像增强非常重要。现有增强方法分为三类：基于
视网膜理论的方法、直方图均衡化方法和使用dehaze模型的方法，这些方法的出现都是为了尽可能地改善图像质量。其中基于直方图均衡化的方法和基于Retinex理论的方法尤其受到很多学者的关注。直方图均衡化方法主要关注增强
图像对比度，但是常常由于过度增强图像产生了不自然的效果，这类方法存在诸如过度增强、图像变白、图像
亮度无法保留等问题。Retinex理论自提出以来就吸引了许多学者的兴趣，它主要关注提高图像亮度，假设图像是照明和反射的相互作用，通过消除照明效果，可以增强低照度图像。Retinex理论随后发展成为单尺度Retinex
算法(Single-scale Retinex，SSR)、多尺度Retinex算法(Multi-scale Retinex，MSR)和带色彩恢复的多尺度Retinex增强算法(MSRCR)。虽然这些方法增强了诸如边缘和
角落的高频信息，可以使细节增强，但它们分别处理RGB通道中的图像，当原始图像不符合“灰色世界假设”时，它们可能会导致
颜色失真，甚至出现“光晕伪影”的现象，而且它们无法有效避免对比度不均匀情况的出现。因此，研究一种不平衡光线下的低照度彩色图像增强方法，是急需解决的问题。
发明内容
[0004] 针对照度不均、整体亮度低的彩色图像，为了解决现有方法对图像增强时出现的暗部区域细节增强不明显以及局部亮区域的过度增强问题，本发明提出了一种基于自适应混沌粒子群优化的低照度彩色图像增强算法，它可以提高对比度并保持低照度彩色图像的自然效果。利用对比度受限的自适应直方图均衡化(contrast limited adaptive histogram equalization，CLAHE)方法处理L*a*b*色彩空间的L*通道图像、增强对比度，采用提出的自适应混沌粒子群
优化算法(Adaptive chaotic particle swarm optimization，ACPSO)寻找出最优伽
马值代入伽马校正函数对HSV色彩空间V通道进行处理、增强亮度，利用构建的自适应拉伸函数对S通道进行拉伸、增强饱和度，最终得到增强的彩色图像。有效改善了低照度彩色图像的亮度和对比度低、饱和度不高，图像细节信息不清晰的问题，明显改善图像中的高亮区域。处理后的图像清晰度高、增强结果自然。
[0005] 本发明所要解决的技术问题在于克服
现有技术的不足，提供一种低照度条件下彩色图像增强方法，解决现有图像增强技术针对低照度且不平衡光线环境下的低照度彩色图像进行增强所产生的图像局部亮区域过度增强以及暗部细节增强效果不佳等问题。
[0006] 本发明具体采用以下技术方案解决上述技术问题：
[0007] S1、读入低照度彩色图像；
[0008] S2、将RGB彩色图像转换至L*a*b*色彩空间，得到亮度分量L*、
色度分量a*和b*；
[0009] S3、采用CLAHE对L*a*b*色彩空间的L*分量处理得到对比度增强的分量L1*；
[0010] S4、将亮度分量L1*和色度分量a*、b*组合得到L*a*b*色彩空间图像，并转换至RGB色彩空间；
[0011] S5、将步骤S4得出的图像从RGB色彩空间转换成HSV色彩空间，得到色调分量H、饱和度分量S和亮度分量V；
[0012] S6、采用ACPSO寻找出最优伽马值，代入伽马校正函数，对步骤S5得出的亮度分量V进行伽马校正处理，得到亮度调整的V1分量；
[0013] S7、使用一种自适应非线性拉伸函数对HSV色彩空间的S分量拉伸，得到饱和度增强的S1分量；
[0014] S8、将亮度分量V1、饱和度分量S1和色调分量H组合得到HSV彩色图像，并进行反变换至RGB色彩空间；
[0015] S9、最后获得经过增强后的彩色图像；
[0016] 优选的，所述S2的从RGB色彩空间到L*a*b*色彩空间之间没有简单的转换公式，RGB色彩空间需要先转为CIEXYZ色彩空间，然后再由CIEXYZ色彩空间转为CIE L*a*b*色彩空间，转换过程具体包括：
[0017] S21、将RGB色彩空间转换到CIEXYZ色彩空间，使用的转换公式如下所示：
[0018] 其中，X、Y、Z分别为CIEXYZ图像中X分量、Y分量、Z分量的值。R、G、B分别为RGB图像中R分量、G分量、B分量的值。
[0019] S22、CIEXYZ色彩空间转换到L*a*b*色彩空间，使用的转换公式如下所示：
[0020] 其中，L*、a*、b*分别为L*a*b*图像中L*分量、a*分量、b*分量的值，Xn，Yn，Zn是标准D65照明白点，取值为Xn＝0.950456，Yn＝1.000000，Zn＝1.088754。
[0021] 优选的，所述S4的从L*a*b*色彩空间重新反转换为RGB色彩空间的具体步骤为：
[0022] S41、L*a*b*色彩空间转换到CIEXYZ色彩空间，使用的转换公式如下所示：
[0023] 这里， 分别为转换之后CIEXYZ图像中X分量、-1
Y分量、Z分量的值，g(t)是f(t)的反函数f (t)，表达形式为：
[0024] S42、CIEXYZ色彩空间转换到RGB色彩空间，使用的转换公式如下所示：
[0025] 其中， 分别为转换之后RGB图像中R分量、G分量、B分量的值。
[0026] 优选的，所述S5的从RGB色彩空间转换成HSV色彩空间，使用的转换公式如下所示：
[0027] 其中，H、S、V分别为HSV图像中的分量，H代表色调，S代表饱和度，V代表亮度， 为步骤S42所得值。
[0028] 优选的，所述S6采用ACPSO寻找出最优伽马值，代入伽马校正函数，对步骤S5得出的亮度分量V进行伽马校正处理，具体步骤为：
[0029] S61、提取步骤S5所得V通道图像；
[0030] S62、设置粒子群优化算法参数并初始化每个粒子的
位置和速度，利用评价函数的初始值初始化个体历史最优适应度值和群体历史最优适应度值，利用粒子初始位置初始化粒子个体最优位置和群体最优位置；
[0031] S63、自适应调整每个粒子的惯性权重，从而更新每个粒子的速度和位置；
[0032] S64、利用每个粒子的位置作为校正参数，带入伽马校正，对步骤S5所得V通道图像进行处理，得到各个粒子的中间增强图像；
[0033] S65、将中间增强图像的信息作为评价函数的输入信息，得到每个粒子的适应度值fitness；
[0034] S66、利用每个粒子的适应度值fitness更新个体历史最优适应度值，并将个体历史最优适应度值对应的粒子位置存储在个体最优位置中；
[0035] S67、利用所有粒子的适应度值的最大值更新群体历史最优适应度值，并将群体历史最优适应度值对应的粒子位置存储在群体最优位置中；
[0036] S68：采用Logistic映射产生混沌序列对粒子位置进行优化；
[0037] S69、判断是否满足
迭代寻优终止条件，若满足，执行步骤S610，否则，执行步骤S63，继续迭代；
[0038] S610、输出粒子群的群体最优位置，结束迭代过程；
[0039] S611、使用步骤S610得到的群体最优位置作为校正参数带入伽马校正，对V通道图像进行校正，得到增强后的图像。
[0040] 优选的 ，所述步骤S62的评价函数为每个粒子的适应度值 ，即：其中，V′是步骤S5所得V通道图像
通过伽马校正转换得到的增强图像，V″由Canny算法处理增强图像V′得到边缘化图像，H(V′)为增强图像V′熵值，n_edges(V″)是增强图像V′经过Canny
边缘检测算子所得的边缘
像素的数量，T是增强图像V′中的总像素数，E(V″)是由Canny算法得到的边缘化图像V″中灰度像素的总和。
[0041] 优选的，所述伽马校正的实现方法为： 其中，lmax为步骤S5所得V通道图像中最大的像素值，l∈[0，lmax]为步骤S5所得V通道图像的实际像素值，校正参数γ的取值为迭代过程中所得的粒子位置，T(l)为步骤S5所得V通道图像中像素值为l的灰度值经过伽马校正后得到的增强图像的灰度值。
[0042] 优选的，所述步骤S62的初始化过程为：
[0043] S621：初始化参数：初始化粒子群粒子个数为N，最大迭代次数为tmax，学习因子c1和c2，混沌系数μ；
[0044] S622：初始化每个粒子的位置和速度：随机产生一个种群，初始化各粒子的位置和速度，将粒子的位置和速度限制在一定区间内，将初始化粒子位置作为校正参数γ的值代入伽马校正函数得到初步增强图像，计算各个粒子的适应度值的初始值fitness0；
[0045] S623：初始化个体最优值和群体最优值：引用步骤S622中每个粒子的适应度值fitness0初始化该粒子个体历史最优适应度值fitnesspbest，并将该粒子的初始随机位置存储为个体最优位置pbest；用各粒子适应度值fitness0中的最大值初始化群体历史最优适应度值fitnessgbest，并将第一个粒子的位置初始化存储为群体最优位置gbest。
[0046] 优选的，所述步骤S63中的惯性权重、速度和位置的更新方法为：
[0047]这里，ωmax和ωmin分别是惯性权重的最大值和最小值，所有粒子ω的最大值和最小值分别设置为0.9和0.4。fitnessmin为当前整个种群适应度值的最小值，fitnessmean是当前整个种群适应度值的平均值，fitnessi为当前第i个粒子的适应度值，i＝(1，2，...，m)。
[0048] 粒子群中粒子i的速度和位置：
[0049] vi(t+1)＝ωtvi(t)+c1r1[pbest(t)-xi(t)]+c2r2[gbest(t)-xi(t)]，[0050] xi(t+1)＝xi(t)+vi(t+1)，
[0051] 其中，ωt表示更新后的惯性权重，r1和r2是[0，1]内的随机数，vi(t+1)和xi(t+1)分别表示粒子i的更新后的速度和位置，xi(t+1)将作为进入下一次迭代时步骤S64对应粒子i的伽马校正的校正参数，pbest(t)为t次迭代时的个体最优位置，gbest(t)为t次迭代时的群体最优位置。
[0052] 优选的，所述步骤S64中将更新后的粒子位置xi作为校正参数γ的值代入伽马校正处理步骤S5所得亮度分量V，得到中间增强图像为Vi′。
[0053] 优选的，所述步骤S65的每个粒子的适应度值fitness的计算步骤为：
[0054] S651：由Canny算法处理中间增强图像Vi′得到
边缘图像V″，计算边缘图像V″像素强度之和E(Vi″)；
[0055] S652：计算中间增强图像Vi′中的边缘像素的数量n_edges(Vi″)，则中间增强图像Vi′的边缘内容为
[0056] S653：中间增强图像Vi′的熵值Hi为： 其中，p(j)是中间增强图像Vi′灰度值为j的像素所占的比例， 其中，Ti是中间增强图像Vi′中的像素总数，kj为中间增强图像Vi′中灰度值为j的像素点的数量，i的范围为i＝1，2，…，N，N为粒子群粒子个数；
[0057] S654：粒子i的适应度值fitnessi为：
[0058]
[0059] 优选的，所述步骤S66的粒子群中每个粒子的个体历史最优适应度值和个体最优位置的更新方法为：若个体历史最优适应度值fitnessgbest小于步骤S654中粒子i的适应度值fitnessi，则更新个体历史最优适应度值fitnesspbest，即将粒子i的适应度值fitnessi赋值给个体历史最优适应度值fitnesspbest，并将该个体历史最优适应度值fitnesspbest对应的粒子的位置作为该粒子的个体最优位置pbest。
[0060] 优选的，所述步骤S67的粒子群的群体历史最优适应度值和群体最优位置的更新方法为：若群体历史最优适应度值fitnessgbest小于粒子群的个体适应度值fitnesspbest中的最大值，则更新群体历史最优适应度值fitnessgbest，即将粒子群所有粒子个体历史最优适应度值fitnesspbest中的最大值赋值给群体历史最优适应度值fitnessgbest，并将该群体历史最优适应度值fitnessgbest对应的粒子的位置作为群体最优位置gbest。
[0061] 优选的，所述步骤S68中利用Logistic混沌序列扰动粒子群，其步骤包括：
[0062] S681：将初始粒子根据公式 映射到[0，1]范围内，其中，gbest∈[ci，di]，ci为粒子i的位置最小值，di为粒子i的位置最大值，t为迭代次数。
[0063] S682：根据Logistic公式生成序列：
[0064]
[0065] 其中，当μ接近4时，Logistic处于混沌状态。可以使用该公式得到N个混沌的解空间。
[0066] S683：将步骤S682中生成的混沌空间中的向量根据式Xi＝ci+yi(di-ci)映射到原始解空间范围内，得到向量组X＝{X1，X2，...，Xn}。
[0067] 优选的，所述步骤S69的迭代寻优终止条件为：迭代次数已达到最大迭代次数tmax。
[0068] 优选的，所述S7的所使用的自适应非线性拉伸函数公式为：
[0069] 式中，S为非线性拉伸前图像饱和度，S1为非线性拉伸后图像饱和度， 和 分别为步骤S42
所得RGB色彩空间对应像素点的 颜色分量的最大值、最小值和平均值。
[0070] 优选的，所述S8中由HSV色彩空间重新变换回RGB色彩空间，使用的变换公式如下所示：
[0071]
[0072]
[0073]
[0074] 式中，C是色度，X是具有这种颜色的第二大成分的中间值。
[0075] 本发明所达到的有益效果：本发明是一种低照度彩色图像增强方法，通过对L*a*b*色彩空间中的亮度L*进行局部处理来提高图像对比度，对HSV色彩空间中的亮度分量V、饱和度分量S进行处理，来消除图像整体亮度偏低、亮度不均、暗处细节不清晰以及饱和度低等问题，解决了图像局部亮区域的过度增强以及暗部细节增强效果不佳的技术问题。与HE，BPDHE，MSR，MSRCR，NPEA，AGCWD和Kanmani等人提出的方法进行比较，其中观察到所提出的方法提供了更好的性能，本发明输出的图像亮度分布均匀，图像色彩更加饱满和自然，清晰度和对比度显著提高。
附图说明
[0076] 为了更清楚地说明本发明
实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0077] 图1为本发明中低照度彩色图像增强方法
流程图；
[0078] 图2为使用ACPSO的亮度增强算法的流程示意图；
[0079] 图3为本发明一较佳实施例中的彩色图像增强结果对比图，其中，(a)为原始图像，(b)为HE的输出结果，(c)为BPDHE的输出结果，(d)为MSR的输出结果，(e)为MSRCR的输出结果，(f)为NPEA的输出结果，(g)为AGCWD的输出结果，(h)为Kanmani提出方法的输出结果，(i)为本发明的输出结果。
具体实施方式
[0080] 为了使本技术领域的人员更好地理解本发明中的技术方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0081] 一种低照度彩色图像增强方法流程图如图1所示，采用CLAHE对L*通道图像进行处理以提高图像对比度，利用ACPSO优化的伽马校正函数对HSV色彩空间中V通道的亮度进行调整，伽马校正可有效提高图像的亮度，提供非常清晰、低噪声的图像。而且，对图像中低亮度区域的亮度提升作用更大，这正好可以弥补CLAHE对于彩色图像暗部细节增强效果不佳的缺点。使用一种自适应拉伸函数对HSV色彩空间的S分量进行拉伸，以提高图像饱和度，增强非均匀照明彩色图像，包括如下步骤：
[0082] S1、输入并读取低照度条件下彩色图像I，具体实现为：利用matlab中的函数imread读取输入的原始彩色图像；
[0083] S2、将RGB彩色图像I转换至L*a*b*色彩空间，得到亮度分量L*、色度分量a*和b*，从RGB色彩空间到L*a*b*色彩空间之间没有简单的转换公式，RGB色彩空间需要先转为CIEXYZ色彩空间，然后再由CIEXYZ色彩空间转为CIE L*a*b*色彩空间，具体转换步骤为：
[0084] S21、将RGB色彩空间转换到CIEXYZ色彩空间，具体是实施时使用的转换公式如下所示：
[0085] 其中，X、Y、Z分别为CIEXYZ图像中X分量、Y分量、Z分量的值。R、G、B分别为RGB图像中R分量、G分量、B分量的值。
[0086] S22、CIEXYZ色彩空间转换到L*a*b*色彩空间，具体是实施时使用的转换公式如下所示：
[0087] 其中，L*、a*、b*分别为CIE L*a*b*图像中L*分量、a*分量、b*分量的值，Xn，Yn，Zn是标准D65照明白点，取值为Xn＝0.950456，Yn＝1.000000，Zn＝1.088754。
[0088] S3、采用对比度受限直方图均衡方法(CLAHE)对L*a*b*色彩空间的L*分量处理得到对比度增强的亮度分量L1*，具体实现为：利用matlab中的函数adapthisteq将图像L*灰度化进行限制对比度直方图均衡化；
[0089] S4、将亮度分量L1*和色度分量a*、b*组合得到L*a*b*色彩空间图像，并转换至RGB彩色图像I1，从L*a*b*色彩空间重新反转换为RGB色彩空间的具体步骤为：
[0090] S41、L*a*b*色彩空间转换到CIEXYZ色彩空间，具体是实施时使用的转换公式如下所示：
[0091] 这里， 分别为转换之后CIEXYZ图像中X分量、Y分量、Z分量的值，g(t)是f(t)的反函数f-1(t)，表达形式为：
[0092]
[0093] S42、CIEXYZ色彩空间转换到RGB色彩空间，得到RGB彩色图像I1，具体是实施时使用的转换公式如下所示：
[0094] 其中， 分别为转换之后RGB图像中R分量、G分量、B分量的值；
[0095] S5、将图像I1从RGB色彩空间转换成HSV色彩空间，得到色调分量H、饱和度分量S和亮度分量V，具体是实施时将彩色图像I1转换成HSV色彩空间使用的转换公式如下所示：
[0096] 其中，H、S、V分别为HSV图像中的分量，H代表色调，S代表饱和度，V代表亮度， 为步骤S42所得值。
[0097] S6、图2为使用ACPSO的亮度增强算法的流程示意图，本发明采用ACPSO寻找出最优伽马值代入伽马校正函数对HSV色彩空间的V分量进行伽马校正处理得到亮度调整的V1分量，具体步骤如下：
[0098] S61、提取步骤S5所得V通道图像，具体实现为：利用matlab中的V＝hsv(：，：，3)语句提取HSV颜色空间中V通道图像；
[0099] S62、设置粒子群优化算法参数并初始化每个粒子的位置和速度，利用评价函数的初始值初始化个体历史最优适应度值和群体历史最优适应度值，利用粒子初始位置初始化粒子个体最优位置和群体最优位置；
[0100] S621：初始化参数：初始化粒子群粒子个数N＝20，最大迭代次数tmax＝100，惯性权重ω最大值ωmax＝0.9，惯性权重ω最小值ωmin＝0.4，学习因子c1＝2，学习因子c2＝2，控制参数μ＝4；
[0101] S622：初始化每个粒子的位置和速度：随机产生一个种群，初始化各粒子的位置和速度，将粒子的位置和速度限制在一定区间内，将所有粒子位置x随机设置在[0，1]之间，所有粒子速度v随机设置在[-1，1]之间；
[0102] S623：利用评价函数计算每个粒子的适应度值，即：
[0103] 其中，V′是原始V通道图像通过伽马校正转换得到的增强图像，V″由Canny算法处理增强图像V′得到，H(V′)为图像V′熵值，n_edges(V″)是增强图像V′经过Canny边缘检测算子所得的边缘像素的数量，T是增强图像V′中的总像素数，E(V″)是由Canny算法得到的图像V″中灰度像素的总和；利用matlab中函数nnz得出；图像熵值、图像边缘内容和图像边缘强度这三个测试函数可以较为全面的反映粒子群算法在各种场景设置下的优化性能，因此，粒子适应度值越大表示作为伽马校正参数的粒子位置越好；粒子适应度值的初始化：在区间[0，1]范围内随机产生粒子位置x0，将粒子位置x0作为伽马值(γ)代入伽马校正函数获取增强图像并计算获得粒子的适应度值的初始值为fitness0；
[0104] S624：初始化个体最优值和群体最优值：引用步骤S23中每个粒子的适应度值fitness0初始化该粒子个体历史最优适应度值fitnesspbest，并将该粒子的初始随机位置存储为个体最优位置pbest；用各粒子适应度值fitness0中的最大值初始化群体历史最优适应度值fitnessgbest，并将第一个粒子的位置初始化存储为群体最优位置gbest。
[0105] S63、自适应调整每个粒子的惯性权重和学习因子，从而更新每个粒子的速度和位置；
[0106] S631：惯性权重：
[0107] 其中，ωmax和ωmin分别是惯性权值的最大值和最小值，fitnessmin为当前整个种群适应度值的最小值，fitnessmean是当前整个种群适应度值的平均值，fitnessi为当前第i个粒子的适应度值，i＝(1，2，...，m)。
[0108] S632：粒子群中粒子i的速度和位置：
[0109] vi(t+1)＝ωtvi(t)+c1r1[pbest(t)-xi(t)]+c2r2[gbest(t)-xi(t)]，[0110] xi(t+1)＝xi(t)+vi(t+1)，
[0111] 其中，ωt表示更新后的惯性权重，r1和r2是[0，1]内的随机数，vi(t+1)和xi(t+1)分别表示粒子i的更新后的速度和位置，xi(t+1)将作为进入下一次迭代时步骤S64对应粒子i的伽马校正的校正参数，pbest(t)为t次迭代时的个体最优位置，gbest(t)为t次迭代时的群体最优位置。
[0112] S64、利用每个粒子的位置作为校正参数带入伽马校正对步骤S5所得V通道图像进行处理，得到各个粒子的中间增强图像Vi′，伽马校正的实现方法为： 其中，lmax为步骤S5所得V通道图像的灰度最大的像素值，可通过matlab中函数max来获取，l∈[0，lmax]为步骤S5所得V通道图像的实际像素值，校正参数γ的取值为迭代过程中所得的粒子位置，首次迭代时校正参数γ取粒子初始化位置，T(l)为步骤S5所得V通道图像中像素值为l的灰度值经过伽马校正后得到的增强图像的灰度值。
[0113] S65、将中间增强图像Vi′的信息作为评价函数的输入信息，得到每个粒子的适应度值fitness，粒子种群的粒子i的适应度值fitnessi的具体步骤为：
[0114] S651：由Canny算法处理中间增强图像Vi′得到边缘化图像V″，计算边缘化图像V″像素强度之和E(Vi″)；
[0115] S652：计算中间增强图像Vi′的边缘内容为 其中，n_edges(Vi″)表示中间增强图像Vi′经过Canny边缘检测算子所得的边缘像素的数量，具体为先调用Canny函数将增强图像Vi′边缘化，利用matlab函数sum针对边缘化后的图像计算出所有边缘像素的数量n_edges(Vi″)，Ti是中间增强图像Vi′中的像素总数；
[0116] S653：中间增强图像Vi′的熵值Hi为： 其中，p(j)是中间增强图像Vi′灰度值为j的像素所占的比例， 其中，Ti是中间增强图像Vi′中的像素总数，对于大小为n×m的中间增强图像Vi′，中间增强图像Vi′的像素总数为n×m，也即，Ti＝n×m，kj为中间增强图像Vi′中灰度值为j的像素点的数量，i的范围为i＝1，2，…，N，V为粒子群粒子个数；
[0117] S654：粒子i的适应度值fitnessi为：
[0118]
[0119] S66、利用每个粒子的适应度值fitness更新个体历史最优适应度值，并将个体历史最优适应度值对应的粒子位置存储在个体最优位置中；具体为：若个体历史最优适应度值fitnessgbest小于步骤S654中粒子i的适应度值fitnessi，则更新个体历史最优适应度值fitnesspbest，即将粒子i的适应度值fitnessi赋值给个体历史最优适应度值fitnesspbest，并将该个体历史最优适应度值fitnesspbest对应的粒子的位置(伽马校正参数)作为该粒子的个体最优位置pbest。
[0120] S67、利用所有粒子的适应度值的最大值更新群体历史最优适应度值，并将群体历史最优适应度值对应的粒子位置存储在群体最优位置中；具体为：若群体历史最优适应度值fitnessgbest小于粒子群的个体适应度值fitnesspbest中的最大值，则更新群体历史最优适应度值fitnessgbest，即将粒子群所有粒子个体历史最优适应度值fitnesspbest中的最大值赋值给群体历史最优适应度值fitnessgbest，并将该群体历史最优适应度值fitnessgbest对应的粒子的位置(伽马校正参数)作为群体最优位置gbest。
[0121] S68：采用Logistic映射产生混沌序列对粒子位置进行优化，其步骤包括：
[0122] S681：将初始粒子根据公式 映射到[0，1]范围内，其中，gbest∈[ci，di]，ci为粒子i的位置最小值，di为粒子i的位置最大值，t为迭代次数。
[0123] S682：根据Logistic公式生成序列：
[0124] 其中，当μ接近4时，Logistic处于混沌状态。可以使用该公式得到N个混沌的解空间。
[0125] S683：将步骤S682中生成的混沌空间中的向量根据式Xi＝ci+yi(di-ci)映射到原始解空间范围内，得到向量组X＝{X1，X2，...，Xn}。
[0126] S69、判断是否满足迭代寻优终止条件，若满足，执行步骤S610，否则，执行步骤S63；迭代寻优终止条件为：迭代次数已达到最大迭代次数tmax。
[0127] S610、输出粒子群的群体最优位置gbest，结束迭代过程；
[0128] S611、使用步骤S610得到的群体最优位置gbest作为校正参数γ带入伽马校正对步骤S5所得V通道图像进行校正，得到增强后的图像V1。本步骤中伽马校正的方法同步骤S64，只需选择不同的校正参数(群体最优位置gbest)即可，校正参数γ控制图像拉伸程度。
[0129] S7、使用一种自适应非线性拉伸函数对HSV色彩空间的S分量拉伸得到饱和度增强的S1分量，本发明所构建的自适应非线性拉伸函数定义为：
[0130] 式中，S为非线性拉伸前图像饱和度，S1为非线性拉伸后图像饱和度， 和 分别为步骤S42
所得RGB色彩空间对应像素点的 颜色分量的最大值、最小值和平均值。
[0131] S8、将亮度分量V1、饱和度分量S1和色调分量H组合得到HSV彩色图像，并进行反变换成RGB彩色图像，具体是实施时使用的变换公式如下所示：
[0132] 式中，C是色度，X是具有这种颜色的第二大成分的中间值。
[0133]
[0134]
[0135] S9、最后获得经过增强后的彩色图像Iout；
[0136] 在本发明的一种实施例中，从大量低照度彩色图像中选取一个场景来进行实验，图像大小为468×262，可验证本发明的适用性与有效性。原始图像的可观测性比较差，采用本发明方法对彩色图像进行增强处理，并与有代表性的对比度增强方法如直方图均衡化(HE)、亮度保持的动态直方图均衡(BPDHE)、多尺度Retinex算法(MSR)、带彩色恢复的多尺度Retinex算法(MSRCR)、非均匀照明图像的自然保持增强算法(NPEA)、基于加权分布的自适应伽马校正的有效对比度增强(AGCWD)和Kanmani等人提出的方法进行增强效果对比，对比结果如图3所示。
[0137] 图像为一群人在街头游行的场景，在这个图像中，由于
建筑物和树木的遮挡，不能清晰的观察所有人和建筑物等
背光面信息，需要同时提高图像对比度和亮度。从增强结果图像可以看出，HE算法对光亮处存在过度增强现象，这影响了原始光亮区域信息，如红色条幅和白色条幅上的标语经过HE处理后均变得模糊不清。BPDHE对图像的亮度有一定提升，但是提升能
力不足，图像整体亮度仍然很低，暗部细节仍然不清晰。MSR和MSRCR处理图像亮度明显提升，但是亮度提升过度而对比度较低，导致MSR处理的图像泛白，MSRCR结果则呈现蓝白色调。NPEA处理图像整体偏暗。AGCWD在图像中间光照处出现过度增强现象。Kanmani等人提出方法的增强效果非常明显，增强后图像的平均亮度也是所有算法中最高的，但是这也带来了增强过度的问题，图像中出现过曝和失真现象，如条幅和人群信息过度增强。本发明提出方法的增强清晰度最好，同时能够平衡处理图像整体亮度和图像细节。
[0138] 为了对本发明提出的方法在处理彩色图像方面进行有效的客观量化评价，分别采用均值(Mean value)、评估函数值(Evaluation Function)、熵(Entropy)和平均梯度(Mean Gradient)四种评价指标对实验结果进行评估。
[0139] 图像亮度为一副图像给人的一种直观感受，对于彩色图像，则跟灰度值有关。均值(Mean)反映了图像的平均亮度，值越高则图像越亮，如果图像亮度值较低的话，画面太黑什么都看不清，而画面过亮常常会令人感觉不适。均值适中(灰度值在128附近)，表明视觉效果良好。均值采用下面公式计算：
[0140] 式中，M×N为图像大小，g(x，y)为图像x行y列处像素灰度值。
[0141] 评估函数值为增强图像的熵值、边缘强度平均值和边缘数三个性能指标的综合，值越大图像增强效果越好，评估函数值(EF)计算由下面公式给出：
[0142]
[0143] 平均梯度(Mean Gradient)反映了图像的清晰度和纹理变化，值越大说明图像越清晰，平均梯度值MG采用下面公式计算：
[0144] 其中，M×N为图像大小， 表示
水平方向梯度，为垂直方向梯度。
[0145] 灰度均值(Gray mean value)反映了图像的平均亮度，值越大说明图像亮度越大。
[0146] 熵(Entropy)可以衡量图像所携带的信息量，图像的熵值越大，代表图像所含的信息量越大，细节越丰富，熵值采用下面公式计算：
[0147] 其中，p(i)某个灰度值(i)在该图像中出现的概率。
[0148] 在本发明的一种实施例中，分别采用指标Mean value(Mean)、Entropy、Evaluation Function(EF)和Mean Gradient(MG)评价本发明提出方法和参与对比的算法对图像增强处理的效果，结果见表1。
[0149] 表1不同增强方法对测试图像的Mean value、Entropy、Evaluation Function和Mean Gradient值定量比较
[0150]
[0151] 经过各种方法处理后，本发明方法和HE均值靠近128，亮度得到很大改善，HE方法会对高亮度区域过度增强，影响高亮度区域细节观察，因此熵值(Entropy)较低。BPDHE处理图像之后，均值较低，说明BPDHE对图像亮度提高不明显。NPEA处理图像平均亮度大于BPDHE，但是数值仍与128差别较大，图像整体亮度仍然偏低。MSR、MSRCR和Kanmani的方法处理图像均值都很大，说明图像过亮。AGCWD同样会出现平均亮度过小现象。本发明方法的Entropy指标在所有比较算法中排名第一，这说明相比其他方法，经本发明方法增强后图像的细节更丰富。本发明方法评估函数值最优，充分表明本发明方法处理图像能够呈现更多图片细节，对比度较高，增强效果最好。在清晰度的对比中，本发明方法的MG指标最优，说明本发明方法可提高图像清晰度。因此整体而言本发明所提出的图像增强方法增强效果明显优于其他方法。
[0152] 以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何
修改、等同替换、改进等，均应包含在本发明的保护范围之内。
