(完整版)信息熵在图像处理特别是图像分割和图像配准中的应用——信息与计算科学毕业设计

摘要

信息论是人们在长期通信实践活动中，由通信技术与概率论、随机过程、数理统计等学科相结合而逐步发展起来的一门新兴交叉学科。而熵是信息论中事件出现概率的不确定性的量度，能有效反映事件包含的信息。随着科学技术，特别是信息技术的迅猛发展，信息理论在通信领域中发挥了越来越重要的作用，由于信息理论解决问题的思路和方法独特、新颖和有效，信息论已渗透到其他科学领域。随着计算机技术和数学理论的不断发展，人工智能、神经网络、遗传算法、模糊理论的不断完善，信息理论的应用越来越广泛。在图像处理研究中，信息熵也越来越受到关注。为了寻找快速有效的图像处理方法，信息理论越来越多地渗透到图像处理技术中。本文通过进一步探讨概论率中熵的概念，分析其在图像处理中的应用，通过概念的分析理解，详细讨论其在图像处理的各个方面：如图像分割、图像配准、人脸识别，特征检测等的应用。

本文介绍了信息熵在图像处理中的应用，总结了一些基于熵的基本概念，互信息的定义。并给出了信息熵在图像处理特别是图像分割和图像配准中的应用，最后实现了信息熵在图像配准中的方法。

关键词：信息熵，互信息，图像分割，图像配准

Abstract

Information theory is a new interdisciplinary subject developed in people long-term communication practice, combining with communication technology, theory of probability, stochastic processes, and mathematical statistics. Entropy is a measure of the uncertainty the probability of the occurrence of the event in the information theory, it can effectively reflect the information event contains. With the development of science and technology, especially the rapid development of information technology, information theory has played a more and more important role in the communication field, because the ideas and methods to solve the problem of information theory is unique, novel and effective, information theory has penetrated into other areas of science. With the development of computer technology and mathematical theory, continuous improvement of artificial intelligence, neural network, genetic algorithm, fuzzy theory, there are more and more extensive applications of information theory. In the research of image processing, the information entropy has attracted more and more attention. In

order to find the fast and effective image processing method, information theory is used more and more frequently in the image processing technology. In this paper, through the further discussion on

concept of entropy, analyzes its application in image processing, such as

image segmentation, image registration, face recognition, feature detection etc.

This paper introduces the application of information entropy in

image processing, summarizes some basic concepts based on the definition of entropy, mutual information. And the information entropy

of image processing especially for image segmentation and image registration. Finally realize the information entropy in image registration.

Keywords：I nformation entropy, Mutual information, Image segmentation,

Image registration

目录

摘要.......................................................................................................................... ...... ABSTRACT .........................................................................................................................目录.............................................................................................................................

1 引言...................................................................................................................................

1.1信息熵的概念.............................................................................................................

1.2信息熵的基本性质及证明.........................................................................................

1.2.1 单峰性..................................................................................................................

1.2.2 对称性..................................................................................................................

1.2.4 展开性..................................................................................................................

1.2.5 确定性.................................................................................................................. 2基于熵的互信息理论 .......................................................................................................

2.1 互信息的概述............................................................................................................

2.2 互信息的定义............................................................................................................

2.3 熵与互信息的关系....................................................................................................

3 信息熵在图像分割中的应用...........................................................................................

3.1图像分割的基本概念 ...............................................................................................

3.1.1图像分割的研究现状 ..........................................................................................

3.1.2 图像分割的方法..................................................................................................

3.2 基于改进粒子群优化的模糊熵煤尘图像分割.........................................................

3.2.1 基本粒子群算法..................................................................................................

3.2.2 改进粒子群优化算法..........................................................................................

3.2.3 Morlet变异 ..........................................................................................................

3.2.4改建粒子群优化的图像分割方法.......................................................................

3.2.5 实验结果及分析..................................................................................................

3.3 一种新信息熵的定义及其在图像分割中的应用 ....................................................

3.3.1香农熵的概念及性质..........................................................................................

3.3.2一种信息熵的定义及证明..................................................................................

3.3.3信息熵计算复杂性分析......................................................................................

3.3.4二维信息熵阈值法..............................................................................................

3.3.5二维信息熵阈值法的复杂性分析......................................................................

3.3.6 结论及分析.........................................................................................................

4 信息熵在图像配准中的应用...........................................................................................

4.2基于互信息的图像配准.............................................................................................

4.3P OWELL算法 ..............................................................................................................

4.4变换.............................................................................................................................

4.4.1平移变换...............................................................................................................

4.4.2旋转变换...............................................................................................................

4.5基于互信息的图像配准的设计与实现.....................................................................

4.5.1总体设计思路和图像配准实现 .........................................................................

4.5.2直方图 ...................................................................................................................

4.5.3联合直方图...........................................................................................................

4.5.4灰度级差值技术 ...................................................................................................

4.4.5优化搜索办法级结论 .......................................................................................... 5结语...............................................................................................................................致谢...............................................................................................................................参考文献...........................................................................................................................

1 引言

1.1.信息熵的概念

1948年，美国科学家香农(C．E．Shannon)发表了一篇著名的论文

《通信的数学理论》。他从研究通信系统传输的实质出发，对信息做了

科学的定义，并进行了定性和定量的描述。

他指出，信息是事物运动状态或存在方式的不确定性的描述。其通

信系统的模型如下所示：

图1.1 信息的传播

信息的基本作用就是消除人们对事物的不确定性。信息熵是信

息论中用于度量信息量的一个概念。假定是随机变量的集合，表示其概率密度，计算此随机变量的信息熵的公式是：

表示一对随机变量的联合密度函数，他们的联合熵可以表示为：

),(log ),(),(Y X p y x p Y X H Y

y x ∑∑∈∈-=χ 信息熵描述的是信源的不确定性，是信源中所有目标的平均信息量。

信息量是信息论的中心概念，将熵作为一个随机事件的不确定性或信息量的量度，它奠定了现代信息论的科学理论基础，如果一条信息是由个字符连成的字符串组成，并且每个字符有种可能，那么这条信息就有种不同的排列情况，那么可以用度量信息量，但这时的信息量随着消息的长度按指数增加，为了使信息量的度量值按线性增加，Hartley 给出了取对数的信息量的定义：

(1.1)

由上式可以看出，信息量随着消息的可能性组合增多而增多，如果

消息只有一种可能性时即事件为必然事件时，那么消息中包含的信息量为零。因此可以看出，可能收到的不同消息越多，对收到哪条消息的不确定性就越大;相反，收到只有一种可能性的消息，不确定性为零，Hartley 对消息的度量实际是对不确定性的度量。

Hartley 度量方法的不足之处是他所定义信息量是假定所有符号发

生的概率相同，但实际情况各符号并不一定都等概发生，为此，Shannon 用概率加权来衡量消息出现的可能性，对Hartley 的度量方法做出改进。

设某一随机过程中有k 种可能的情况，每种情况发生的概率分别

是，，…，，Shannon 给出了熵的如下定义：

∑∑-==i i i

i p p p p H 22l o g 1l o g (1.2)

当所有可能的事件均以相等的概率发生时，上式就成了Hartley 定

义的熵，并且这时熵取得最大值，即

∑∑==-=n n n

n n m m m m m H 222log log 11log 1 (1.3)

所以，Hartley 熵是,Shannon 熵的特殊情形，而Shannon 更具有一般性。

Shannon 熵包含三种含义：第一种含义是度量信息量，事件发生概

率与获得的信息量成反比，即概率越大，信息量越少，又由式(1.3)知，概率越大，信息量越少，熵越小，所以可用熵的大小来度量信息量，熵

越大，信息量越大;第二是度量事件概率分布的分散度，概率集中分布

时熵值小，分散性越强，熵越大；三含义是度量事件发生的不确定性，概率越大，事件的不确定性越小，熵越小。利用上面第三个含义，可以

用Shannon 熵，来度量图像包含的信息量，图像灰度值的概率分布是每

灰度值出现的次数除以图像中所有灰度值出现的总次数，此时图像的信

息量可依据这个概率分布来计算，一幅图像中不同的灰度值较少，各灰

度值出现的概率较高，则对应的灰度值较低，意味着这幅图像含有的信

息量很少。反之，如果一幅图像中含有很多不同的灰度值，且各灰度值

发生的概率又基本一致，则它的熵值会很高，那么这幅图像包含的信息

量很大。

1.2信息熵的基本性质及证明

1.2.1单峰性

信息熵的单峰性可表述为：先考察由、两个事件构成的概率系统，

其产生的概率分别为和则该系统的信息

)).1(log )1(log (22P P P P H --+-=

通过求极限不难证明：

(1) 当时，.0))01(log )01(log 0(22=--+-=H 这是一种产生的概率

为0，产生的概率为1 的确定系统。

(2) 当时.0))11(log )11(1log 1(22=--+-=H 这是一种产生的概率

为1，产生的概率为0 的确定系统。

(3) 对函数)).1(log )1(log (22P P P P H --+-=可以通过求导数的方

式寻找其极值点。该函数的一阶导数为令则有，求得为该函数的驻点。因为二阶导数当时，恒小于0 ， 所以当时函数有极大值。这说明当、

两事件产生的概率相同时，具有最大值，这是一种不确定性最大的不确定系统。

(4) 若概率系统中有个事件，当每一事件产生的概率相同（均为）

时，则系统的信息熵具有最大值。该结论可以通过以下的讨论来证明：

具有个事件的概率系统其信息熵可表示为，这是在约束条件下的极

值问题。应用因子法，设：

∑∑==-+-=n i n

i i i i P P P H 1120).1(log λ

将对事件的概率求一阶偏导数，并令使用约束条件确定值，可求得

（常数）。同理有（常数），即当时，有极大值。

1.2.2 对称性

信息熵的对称性可表述为：设某一概率系统中n 个事件的概率分布

为当对事件位置的顺序进行任意置换后，得到新的概率分布为，并有以

下关系成立：

).,,(),(2121n n P P

P H P P P H '''= 它表示概率系统中事件的顺序虽不同，但概率系统的熵H 是不变

的，即概率系统的熵与事件的顺序无关。

1.2.3 渐化性

信息熵的渐化性可表述为：设概率为的事件可分解为概率分别为和

的两个事件，则有：

).,()(),,,()

,,,(),,(121121121r

q r r q q H r q r q P P P H r q P P P P H P P P P H n n n n n ++?++=+=--- 1.2.4 展开性

信息熵的展开性可表述为：设某一概率系统的概率分布为则系统的

信息熵具有展开性质：

).0,,(),(2121n n P P P H P P P H =

在此基础上，进一步展开有：

根据上述展开性不难证明。

1.2.5 确定性

信息熵的确定性可表述为：设信息系统中，任一事件产生的概率为1，

则其他事件产生的概率为0。这是一种确定的系统，对于这样的系统有： .0)1,0,,0()0,0,1,0,,0()0,,0,1(=== H H H

根据很容易证明上述性质。

2 基于熵的互信息理论

2.1 互信息的概述

互信息（Mutual Information ）来自于信息论，是信息论中的一个

基本概念，是两个随机变量统计相关性的测度。当两幅图像达到最佳配

准，它们对应像素的灰度互信息应达到最大。该测度不需要对不同成像

模式下图像灰度间的关系作任何假设，也不需要对图像进行分割或任何

预处理，具有自动化程度高的特点。因此，最近几年将互信息作为图像

配准过程的相似性测度，利用最大互信息法进行图像配准成为了图像处

理领域的研究热点。

互信息是基于概率统计论提出的，具有统计特性，它被多数研究者

公认为是一个很好的图像配准准则，许多图像配准算法的研究均是在互

信息的基础上加以改进的。互信息作为医学图像配准的一个相似性测

度，多模态医学图像的配准很实用，其配准原理是两幅基于共同人体解

剖结构的图像在配准时具有最大的互信息值。

2.2 互信息定义

定义1：随机变量和之间的互信息定义为：

或定义互信息为：

可以证明二者是相等的，即=。因此，和是随机变量和之间相互提供的

信息。

另一种定义：也可以采用直接定义与之间的互信息为：

∑∑===K k J j j k j k j k b p a p b a p b a p Y X I 11)()(),(log

),();(

=

))

(

)

(

/)

,

(

log(

)

,

(

))

(

)

(

|)

,

(

(Y

p

X

p

Y

X

p

Y

X

E

y

p

x

p

y

x

p

D=

可直接导出)

(

)

(

)

(

)

,

(XY

H

Y

H

X

H

Y

X

I-

+

=及

)).

(

),

(

min(

)

;

(

0Y

H

X

H

Y

X

I≤

≤

2.3 熵与互信息的关系

(1) 独立：)

(

)

|

(

)

(

)

|

(Y

H

X

Y

H

X

H

Y

X

H=

=，，有

(2) 确定：，则.

从而，互信息是随机变量之间相互依存度的度量信息。

互信息是信息论中的一个基本概念，通常用于描述两个系统间的相关性，或者是一个系统中所包含的另一个系统信息的多少，是两个随机变量和之间统计相关性的量度，或是一个变量包含另一个变量的信息量的量度。

它可以用熵和来描述以及联合熵，

)

|

(

)

(

)

|

(

)

(

)

,

(

)

(

)

(

)

,

(A

B

H

B

H

B

A

H

A

H

B

A

H

B

H

A

H

B

A

I-

=

-

=

-

+

=

(2.1)

其中和分别是系统和的熵，是，的联合熵，表示一直系统时的条件熵和一直系统时的条件熵。上述各种熵可分别表示为：

∑-=

a

A

A

a

p

a

p

A

H)

(

log

)

(

)

(

2

(2.2)

∑-=

b

B

B

b

p

b

p

B

H)

(

log

)

(

)

(

2

(2.3)

∑-=b

a AB AB

b a p b a p B A H ,2),(log ),(),(

(2.4)

∑-=b

a B A AB

b a p b a p B A H ,|2)|(log ),()|(

(2.5)

∑-=b

a A B AB a

b p b a p A B H ,|2)|(log ),()|(

(2.6)

其中，， 和分别是系统和系统完全独立时的边缘概率密度，是系统

和的联合概率分布，是已知系统时的条件概率分布，是已知系统时的条件概率分布，如果联合概率分布密度满足，则随机变量和相互独立;如

果和满足某映射关系使))(,())(()(a T a p a T p a p AB B A ==，则随机变量和最大相关。

在通信系统中，信源和信宿是相互联系的，因此，收到的条件下，

对信源具有一定的了解，但仍然对有不确定度，即条件熵，但总小于绝对熵。对信源的了解程度（确定度）为得到结论：差值度量了确定度。同样，在确值信源发送X 的条件下，差值度量了对的了解程度。

3 信息熵在图像分割中的应用

3.1 图像分割的基本概念

图像分割是图像处理和分析的关键步骤，也是一种基本的计算机视

觉技术。当今信息熵主要应用在图像分割技术中。为了识别和分析目标，图像分割把图像分各具特性的区域。这些特性可以是灰度、颜色、纹理等，目标可以对应单个区域，也可以对应多个区域。基于熵的图像分割方法，尽可能减少了图像信息的损失，因此可用于复杂背景，而且这种

方法有很多。随着计算机技术和数学理论的不断发展，人工智能、神经网络、遗传算法、模糊理论的不断完善，以及处理的图像越来越复杂，单一的方法已不能满足人们的需求，因此，研究多方法的结合是这一领域的趋势。

3.1.1 图像分割的研究现状

图像分割是图像处理中的一项关键技术，也是一经典难题，自20世纪70年代起一直受到人们的高度重视，至今已提出了上千种分割算法。但发展至今仍没有找出一个通用的分割理论，现提出的分割算法大都是针对具体问题的，并没有一种适合所有图像的通用分割算法。另外，也还没有制定出判断分割算法好坏和选择适用分割算法的标准，这给图像分割技术的应用带来许多实际问题。

3.1.2图像分割的方法

(1) 基于阈值的分割

这是一种最常用的区域分割技术，阈值是用于区分不同目标的灰度值。如果图像只有目标和背景两大类，那么只需选取一个阈值称为单阈值分割。这种方法是将图像中每个像素的灰度值和阈值比较，灰度值大于阈值的像素为一类，灰度值小于阈值的像素为另一类。如果图像中有多个目标，就需要选取多个阈值将各个目标分开，这种方法称为多阈值分割。为区分目标还需要对多个区域进行标记。阈值又可分为全局阈值，局部阈值和动态阈值，阈值分割的结果依赖于阈值的选取，确定阈值是阈值分割的关键，阈值分割实质上就是按照某个标准求出最佳阈值的过程。常用的全局阈值选取方法有利用图像灰度直方图的峰谷法，最小误差法，最大类间方差法，最大熵自动阈值法以及其他一些方法。

(2) 基于区域的分割

基于区域的分割技术有两种基本形式：区域生长和分裂合并。前者是从单像素出发，逐渐合并以形成所需的分割结果。后者是从整个图像出发，逐渐分裂或合并以形成所需要的分割结果。与阈值方法不同，这类方法不但考虑了像素的相似性，还考虑了空间上的邻接性，因此可以有效地消除孤立噪声的干扰，具有很强的鲁棒性。而且，无论是分裂还是合并，都能将分割深入达到像素级。

(3) 基于边缘的分割

基于边界的分割方法是利用不同区域间像素灰度不连续的特点检测出区域间的边缘，从而实现图像分割。边界的像素灰度值变化往往比较剧烈。首先检测图像中的边缘点，在按一定策略连接成轮廓，从而构成分割区域。边缘检测技术可以按照处理的顺序分为串行边缘检测及并行边缘检测。在穿性边缘检测中，当前像素点是否属于欲检测的边缘取决于先前像素的验证结果；而在并行边缘检测技术中，一个像素点是否属于欲检测的边缘，取决于当前正在检测的像素点以及该像素点的一些相邻像素点，这样该模型可以同时用于检测图像中的所有像素点。

3.2 基于改进粒子群优化的模糊熵煤尘图像分割

对煤尘图像进行有效的分割是煤尘浓度测量的重要研究内容之一，在煤尘图像分析和识别中具有重要意义。阈值法是最常用的图像分割方法，其关键是阈值的选取，用图像灰度模糊熵来确定分割阈值是一种有效的阈值确定方法。图像灰度模糊熵中模糊参数的寻优实际上是一个优化问题。解决优化问题的方法通常有穷举法、遗传算法、蚁群算法、粒子群算法等，其中Kennedy 和Eberhart提出的粒子群优化算法(PSO) 因其优越性而成为研究的热点。普通粒子群优化算法存在易陷入局部最优以及过早收敛的缺点，使得该算法难以得到理想的优化效果。近年

来出现了不少改进的PSO 算法，改进算法主要有对惯性因子的改进，

以及引入遗传算法中的交叉、变异或进化思想对部分粒子进行相应的操

作。Li 等人提出的高斯变异粒子群(GMPSO) 算法取得了不错的分割效

果，但该算法的分割精度还有待进一步提高。用一种综合Morlet 变异

和惯性因子自适应的改进粒子群优化算法，让该算法和模糊熵结合应用

于图像分割，利用改进粒子群优化(IPSO) 算法来搜索，使模糊熵最大

时的参数值得到模糊参数的最优组合，进而确定图像的分割阈值。实验

结果表明，该算法取得了令人满意的分割结果，算法运算时间较小，能满足对煤尘浓度实时精确测量的要求。

3.2.1 基本粒子群算法

粒子群优化(PSO) 算法是一种进化计算技术，最早由Kenney 和

Eberhart 于1995 年提出的。源于对鸟群捕食行为研究的PSO 算法同

遗传算法类似，是一种基于迭代的优化工具。在PSO 算法中，每个个

体称为粒子，所有的粒子都有一个由被优化的函数决定的适应值，每个

粒子还有一个速度决定他们飞翔的方向和距离，然后粒子就追随当前的

最优粒子在解空间搜索。假定粒子规模为搜索空间为维，则第个粒子的

位置表示为，第个粒子的速度表示为，每个粒子具有最好适应值的位置

称为个体最好位置，记为，整个种群中的最好适应值位置称为全局最好

位置，记为。 在找到这两个最优值时，每个粒子根据如下公式来更新

自己的速度和位置：

()())),()(())()(()()1(21t x t p rand c t x t p rand c t wv t v ij gj ij ij ij ij -+-+=+

(3.1)

(3.2)

其中：表示第个粒子；表示粒子的第维；表示第次迭代；为加速常

数，通常在间取值；为均匀分布在(0， 1) 上的随机数；为惯性因子。

3.2.2 改进粒子群优化算法

普通粒子群优化算法存在易陷入局部最优以及过早收敛的缺点，使

得该算法难以得到理想的优化效果。近年来出现了不少改进的PSO 算

法，改进算法主要有对惯性因子的改进，以及引入遗传算法中的交叉、变异或进化思想对部分粒子进行相应的操作。本文针对惯性因子的改进

提出了惯性因子自适应算法，同时引入了Morlet 变异操作，克服了普

通粒子群优化算法存在易陷入局部最优以及过早收敛的缺点。惯性因子

自适应粒子群算法，当粒子群中大多数粒子在连续的迭代中未找到最优

值前停止更新时，就会出现过早收敛的现像。当惯性因子较小或固定时

也会出现这种现像，从式(3.1) 可以看出， 当较小并且 和很小时，也

很小，即相应的粒子失去搜索能力。这种情况通常会出现在当粒子本身

是全局最优时即和等于零时的迭代早期阶段，这样在以后的迭代中粒子

就失去了多样性。为了解决该问题，一般将设为：

./)(max min max t t ax m ωωωω-?-=

(3.3)

其中：表示总迭代次数，和分别表示最大和最小惯性因子。本文对

上述方法作了改进，根据粒子距离全局最优值之间的距离对进行调节，即

))./(1(max 0dist dist i -?=ωω

(3.4)

其中：为当前粒子距离全局最优值之间的欧几里得距离，即 (4)为

最大距离。这样对进行调节能保证粒子在偏离全局最优时，粒子和全局

最优值之间的吸引力将保证粒子不会偏离最优值太远，从而避免出现过

早收敛的现像。

3.2.3 Morlet 变异

为了克服过早收敛，还有一种方法就是引入遗传算法中的变异操

作，即当用基本PSO 算法对粒子的位置和速度进行更新后，再对部分

粒子进行变异操作，使得粒子种群呈现多样性。 一般可用均匀变异或

非均匀变异来进行变异操作，Natsuki 引入了高斯变异操作，有：

)).(1()())((σgaussian t x t x mut ij ij +?=

(3.5)

其中：为高斯变异后的位置；为的高斯分布，

Natsuki 指出可以对粒子以一定概率进行高斯变异，也可以当粒子

的位置停止更新时进行高斯变异。本文中的Morlet 变异能对粒子起到微

调的作用，

每个粒子变异的概率为的大小根据粒子群的维数决定。Morlet 变

异的方程式如下：

?????≤-?+>-?+=.0),)(()(;0)),(()())((min max σσσσx t x t x t x x t x t x m ut ij ij

ij ij ij (3.6)

其中：为变异后的，和分别为的最大最小值，的计算公式如下：

(3.7)

其中：

.)ln()1)(ln(m ax g t t g wm e a +--=ξ

(3.8)

这里： 为上式单调递增方程的形状参数，为的上限值，为当前迭

代次数，为最大迭代次数。

3.2.4 改进粒子群优化的模糊熵图像分割算法

图像的最大模糊熵：

根据模糊理论，图像可看成是一个模糊事件。根据模糊熵理论，若

分割阈值将原始图像的像素分成黑和亮两个模糊集，则这两个集合的隶

属函数和模糊熵分别为：

???????????<≤

<-?--≤<-?---

≤.

0)()()()()(1,1)

(22

k c c k b b c a c c k b k a a b a c a k a k k d ，；

，）（；

，；

μ

(3.9)

???????????<≤<-?---

≤<-?--≤.

1)()(1)()()(,0)(2

2

k c c k b b c a c c k b k a a b a c a k a k k b ，；

，）（；

，；

μ

(3.10) ().)

(ln )

(255

0∑==k d

d k d d k d p k p p k p H μμ

(3.11) ().)(ln )(255

0∑==k b

b k b b k b p k p p k p H μ

μ

(3.12)

其中：参数满足)/(2/)(n m n p c a b k k ?=+=；这里为中元素的个数。模

糊事件的总模糊熵：

(3.13)

由熵理论可知，为了实现目标与背景的最佳分割，模糊事件的模

糊熵应为最大，

即得到使总模糊熵达到最大值时对应，并据此确定最优阈值

(3.14)

改进粒子群优化的模糊熵图像分割算法：根据最大模糊熵原理，基

于最大模糊熵的图像分割算法其本质是在图像的整个灰度空间上搜索

一组参数使图像的总模糊熵取最大值的优化问题。并且将改进粒子群优

化(IPSO) 算法用于搜索一组最优参数，提高了算法的分割性能。

算法的基本步骤如下：

Step 1： 初始化。初始化粒子群的位置矩阵和速度矩阵，设定粒子

群规模

和维数 (由于需寻优2 个参数，)，设定

，)()(min max min rand x x x x ij ?-+= (3.15)

(3.16)

，

)(2max max rand v v v ij ?+-= (3.17)

(3.18)

其中：为均匀分布在(0， 1) 上的随机数；为的最大值;和分别为的最大、最小值，一般取，这里和分别为图像的最大、最小灰度。

Step 2： 选择式(3.13) 作为粒子群算法的适应度函数，计算粒子

群中每个粒子的适应值， 并根据适应值选择每个粒子的当前最好位置

Pi 和粒子群的全局最好位置。

Step 3：根据式(3.3) 计算权重因子，再根据式(3.1) 和(3.2) 更新粒子的速度和位置。

Step 4：根据式(3.6) 以一定概率对部分粒子进行Morlet 变异。

Step 5：若达到最大迭代次数，则算法结束；否则，转Step 2。

Step 6：求出全局最优解对应的参数组合，计算分割阈值对图像进行分割。

3.2.5 实验结果及分析

利用本文算法对不同类型图像进行分割实验，并与其他算法的结果进行对比。实验中粒子群算法相关参数选择如下：粒子群规模，维数，最大迭代次，普通PSO 算法中的惯性因子，学习因子.Morlet变异参数为：随即抽取50% 的粒子进行变异操作，即变异概率实验中采用的图像分别为Lena，Boat，共3幅真实煤尘图像，它们代表几种不同类型的图像。

(a) Lena 图像(b) Boat 图像

(c) 煤尘图像1 (d) 煤尘图像2

(e) 煤尘图像3

图3.1 实验图像直方图

它们的灰度直方图，Lena 图像呈多峰模式；Boat 图像为明显的双峰; 3 幅煤尘图像为单峰模式。利用本文(IPSO) 算法和基本PSO 算法，