指针式仪表自动读数识别系统设计

1、精选优质文档-倾情为你奉上精选优质文档-倾情为你奉上专心-专注-专业专心-专注-专业精选优质文档-倾情为你奉上专心-专注-专业指针式仪表自动读数系统设计摘要随着模式识别技术、计算机技术等多种技术的不断完善和发展，机器视觉获得了巨大的进步与发展。目前在许多企业中，存在着大量的仪表，仪表的读数都要靠人来完成，工作量很大而且误差率相对来说比较高，基于这个原因，设计了一个工业生产线在线检测数据数字化处理系统。首先通过摄像头采集仪表图像，通过传输装置以无线的方式把图像传输到电脑上，然后在电脑上通过matlab程序设计，处理图像，读取仪表的数据，然后用labview调用matlab程序，并设计一个显示界面

2、，通过界面可以看到实时的数据与仪表图像以及在一段时间内数据的变化情况。最后通过labview访问access数据库，将读取的数据传入数据库中，便于被调用，最终实现检测数据的数字化处理。关键词：机器视觉 图像处理 Matlab LabviewOnline testing data of industrial production line digital processing system designAbstract With the development of computer technology and pattern recognition technology,machine vi

3、sion technology makes a great progress and develop-ment.At present,there are a lot of instrument in many enterprises.in-strument reading work need people to complete.So there are a lot of workto do and efficiency is very low,as the same time,error rate is quitehigh.For this reason,there design a onl

4、ine testing data of industrial production line digital processing system.First of all,there need to takePicture by camera.next,through transmission device in wireless way to transmit the image to a computer. Then by matlab programming on the compu-ter,data processing images, read the meter.And then

5、call matlab by labview and design a display interface.Through the interface can see the real-time data.At last,through labview access access database,and take the data into database.Finally,realize the testing data of the digital processing.Keywords： Machine vision image processing Matlab Labview目 录

6、TOC o 1-4 h u 1.前言1.1课题背景近年来，随着科学技术的进步，许多新兴的技术不断的产生与发展，机器视觉技术也随着图像处理等技术的发展而逐渐产生，而且发展很快。近些年来，世界范围内涌现了许多的专门从事机器视觉技术的公司，有许多的公司得到了很好的利润。现阶段，国内的机器视觉技术虽然刚起步不长的时间，但是，也有几家公司在这方面取得了不错的成绩，并且发展的速度也很快。所谓的机器视觉是以通过处理图像达到类似于人眼的视觉的效果，提高工作效率减轻人的工作负担，为最终目的。相比于人类的视觉，首先机器视觉不会出现疲劳现象，而且机器视觉识别的精度可以比人眼的识别的精度还要高一些。目前，很多领域中都

7、出现了机器视觉技术的影子，这也体现出了机器视觉技术的研究价值。主要应用的领域有:医学领域、工业检测领域、监控管理、航天领域等。当机器视觉技术应用于工业检测领域时，可以实现工业生产线在线检测数据的自动读取，提高读数的准确性，有效减轻人的负担，及时准确的发现问题、解决问题。进入21世纪后，电子类的仪表应用的越来越多，准确度更高的数字式仪表大有取代传统的指针式仪表的趋势。但是，数字式的仪表在某些情况下不能使用，对环境的适应性比较差，比如，在读数快速变化的时候就不适合使用了，而且，相比于指针式仪表来说，数字仪表价格相对较高，大批量的购买成本会上升很多。指针式仪表则可以适用于很多的场合，环境因素对其约束

8、相对来说不大。除此之外，指针式仪表还有许多其他的优点，比如，抗干扰性强，防尘，防水等，因此我们国家的许多行业中，尤其是那些严禁易燃易爆的场合，比如电力、石油、化工等，指针式仪表仍然是首选，在环境恶劣的场合，根本无法使用数字式的仪表。在一般的工厂中，指针式的仪表数量大、种类多，这些仪表的读数十分的麻烦，而且工作量很大，工作效率很低。传统上人们是检定指针式仪表的方法是人工判读，但是这种判别方法会受到许多人为因素的干扰，比如人的观测仪表角度、疲劳强度以及观测的距离，具有工作强度大，误差率大等缺点，不能实现仪表读数与检测的自动化。每个仪表都需要依靠人来瞄准，每个仪表的操作也需要依靠人来完成，这种方式的

9、工作量大，效率低。 为了缩短检测的时间，有效减轻工作人员工作强度 ，提高检测工作的自动化水平，设计工业生产线在线检测数据数字化处理系统是很有必要的。利用机器视觉技术，可以实现仪表图像的采集、识别以及仪表数据的读取，而整个设计的关键就在于如何用图像识别的方式识别仪表，读取数据。1.2课题研究现状 工业生产线在线监测数据的数字化处理系统的设计的核心是使用图像识别的方法识别仪表,指针式仪表的读数识别研究起步相对来说较晚，1994年，sablatnig等人提出了利用hough变换来检测指针式仪表的读数的方法1 Sablatnig, Robert, Kropatsch, Walter G. Automa

10、tic Reading of Analog DisplayInstruments.Conference on Pattern Recognition, 1994, 1：794-797.。1995 年，韩国的 Kyong-Ho Kim 等人成功的实现了对核装置数字表的图像分割，采用的方法是区域划分和阈值法，最终采取一系列的措施，建立一个完全自动化仪表识别系统2 Kyong-Ho Kim, Sung-Li Chen, Yong-Bum Lee, Jong Min Lee. A Study on Analog and Digital Meter Recognition Based on Image

11、Processing Technique. Journal of the Korean Institute of elematics and Eleatronics, 1995, 9(32)：79-94. 。Correa Alegria等人完成了利用机器视觉技术对指针式仪表的检定，整个过程为，首先使用摄像头采集图像，然后采用减影法处理图像，保存图像，然后利用霍夫变换检测直线，识别指针，并且读出指针的角度然，最后，根据指针角度通过计算，得到指针的读数3 F.C0rreaAlegria,A.Cruzz.Serra.Automatie Calibration of Analog and D1 Sab

12、latnig, Robert, Kropatsch, Walter G. Automatic Reading of Analog DisplayInstruments.Conference on Pattern Recognition, 1994, 1：794-797.2 Kyong-Ho Kim, Sung-Li Chen, Yong-Bum Lee, Jong Min Lee. A Study on Analog and Digital Meter Recognition Based on Image Processing Technique. Journal of the Korean

13、Institute of elematics and Eleatronics, 1995, 9(32)：79-94. 3 F.C0rreaAlegria,A.Cruzz.Serra.Automatie Calibration of Analog and Digital Measuring Instruments Using ComPuter VisionJ.IEEE Transaction on Instrumentation and Measurement,2000,49(l):94一99.4 孙凤杰，安田江，范杰清，杨春萍，徐征. 电力变压器温度表指针位置识别研究J. 中 国电机工程学报,

14、2007,27(7): 70-75.5 黎明，王厚枢.飞机零位检测系统中的座舱读数识别法团.南京航空航天大学学报，1996， l(28)：106-112.6 周洪，钟明慧.工业过程显示的数字处理与视频传输.仪表仪表学报(增刊)，2003， 4(24):359-361.577 杨晓敏，何小海，吴炜等.基于高斯混合模型的车辆字符识别算法.光电子激光，2007， 4:487-490.8 何智杰 , 张彬 , 金连文 . 高精度指针仪表自动读数识别方法 J. 计算机辅助工 程 ,2006,15(3): 9-18.9 宁志刚, 汪仁煌, 唐冬, 李貌. 基于图像处理的高精度仪表判读方法J. 传感器与微

15、系统,2006,25(3): 32-37.1.3设计的目的与意义 对于大多数的指针式仪表，尤其是准确度比较高的仪表的检验，仍然需要通过人来实现，这其中夹杂了许多的人为因素，不同的人对于同一个仪表可能会读出不同的数据，对于准确度比较高的仪表，会有比较大的误差率，而且人会出现视觉疲劳，增加了数值被读错的可能性，造成人为的误差。所以，设计一种能自动读取指针式仪表读数的系统，就可以大大减少人的工作量 ，降低工作的强度，提高读数的准确性，提高工作效率。工业设备能耗参数指标体系建立后，通过计算机就可以清楚的了解到各设备的能耗情况以及工作状态及时准确的发现问题，这样可以减少能源的消耗，优化能源结构，提高能源

16、利用率，提高生产效率，具有良好的经济效益和社会效益。1.4设计需要解决的问题对组成工业生产线在线检测数据数字化处理系统的各种软硬件进行综合设计，硬件部分包括摄像头、图像传输设备、图像接受设备的选择;软件部分包括对接收到的图像进行处理的相关程序，读数显示界面的相关程序。研究准确读取仪表数据的方法：获得较高质量表盘图像，有效的把指针检测出来，提高读数的准确性。.数据显示界面的设计：能够看到此时此刻的仪表图像以及对应的表盘读数，能够看出仪表的读数在一段时间内的变化情况。工业生产线在线检测数据数字化处理系统总体方案设计2.1设计方案工业生产线在线检测数据数字化处理系统主要由图像采集装置、图像发射与接受

17、装置、图像处理部分以及客户端等部分组成。设计方案如下所示： 目标仪表 图像采集 发射器 接收器 计算机 图像滤波去噪 图像增强 图像二值化 数据库 客户端 边缘检测 图像处理 图像分割 指针识别 读数识别 图2-1 总体方案 图像发射与接收装置采用的是无线影音发射器，图像处理部分是由matlab来完成的，而客户端则是由labview来设计的。2.2可行性分析通过摄像头可以采集到清晰的图像，无线影音传输器可实现图像的远距离传输，用matlab程序可以处理图像，读取数据，labview是一种程序开发环境,可以实现显示界面的搭建，所以，该方案可行。工业生产线在线检测数据数字化处理系统设计综述3.1图

18、像采集与传输装置设计3.1.1图像采集表盘图像的采集选择的是探头摄像头，该摄像头的参数如下图所示，该摄像头拥有六颗直径5mm的红外灯，可以实现夜拍的功能，视频线选择铜芯铜网的全铜视频线，保证传输图像不受干扰。 图3-1 摄像头参数3.1.2图像传输系统图像传输选取的是无线影音传输器，无线影音传输器是一种可以把图像和声音无线传输到电脑上的无线设备，该装置如下图所示。 图3-2 无线影音传输器无线影音传输器共有六个部分：发射机、接受机、全向天线、音视频线、12V2A电源以及12V1A电源。表3-1是发射机的技术参数,表3-2为接收机的参数。 表3-1 发射机参数 发射频率3W可用频道8CH工作频率

19、2370MHz-2510MHz频率调制方式FM/FSK视频制式PAL/NTSC音频输入阻抗1K视频输入阻抗75视频输入电平1Vp_p工作电压12V工作电流550m-600mA工作温度-20C-120C工作湿度85%RH 表3-2 接收机参数接收天线7DBi全向天线可用频道8CH工作频率2370MHz-2510MHz频率调制方式FM/FSK接受灵敏度-90dbm音频输入阻抗1K视频输入阻抗75视频输入电平1Vp_p工作电压12V工作电流150mA-200mA工作温度-20C-120C工作湿度85%RH 无线影音传输器的特点如下：具有微波抗干扰功能比较好；传输和接受的品质好，接收到的图像清晰稳定；

20、多种频率定制，适应不同的电磁环境；可以在同一地方使用八套，信号不会相互干扰可以传送的信号包括DVD、VCD、CD、电视、录像等等；具有广泛的应用范围；传输距离远；经济实惠；维护简单方便，可配合移动电源供电。该传输器的安装方法为：把发射器（2A）和电源适配器连接起来将，然后通电，如果红灯亮了，则说明发射器通电成功。将另一电源适配器接于接收器输入端（1A），接上电源，红灯亮，则证明接收机通电成功。把摄像头和发射器用音/视频线连接起来，接收器连在显示器上，把发射器和接收器的频道选择开关调到同一位置，通过显示器可以看到摄像头传来的画面，完成安装。 图3-3 传输器安装示意图视频采集卡选用的是4路USB

21、视频采集卡，如下图所示，该视频采集卡具有以下功能：该采集卡带有4路视频，4路音频，录像D1，回放D1；支持电脑操作系统，XP系统，Win7系统；支持网络远程监控，自带动态域名，设置简单；可以支持多种录像模式，如手动录像，定时录像等；每路可单独设置工作时间表，支持可编程定时录像，循环录像；支持画面捕捉、画面保存以及画面备份功能；支持远程回放功能。 图3-4 视频采集卡 3.1.3图像采集与传输系统的设计图像采集与传输系统的设计有两种方案，如下图所示。方案一：摄像头+电源+支架+BNC接头+无线接收器+USB采集卡+电脑+显示器方案二：摄像头+电源+支架+BNG接头+无线接收器+USB采集卡+硬盘

22、录像机+硬盘+显示器在电脑上操作更加灵活方便，而且图像处理也是在电脑上进行的，所以，选择方案一。 图3-5 传输系统设计示意图经过在工厂的实地测试，该图像采集与传输系统可以实现无线传输的功能，并且具有一定的抗干扰性，下图所示的图片是在工厂中距离仪表大约100米远的楼上所接收到的图片。 图3-6 工厂实地测试接收画面3.2图像处理与仪表数据的读取3.2.1图像预处理 在图像采集的过程中，会受到许多的干扰，包括自身的和外界的，例如，外界光照条件的影响，摄像头像素的高低，仪表表盘的清洁程度以及各种振动及噪声的干扰。受到诸多的干扰后得到的照片可能会是模糊的，可能无法获取我们需要的信息，要想得到有用的信

23、息就必须想办法去掉这些干扰。我们通过摄像头的到的图像都是彩色的，里面蕴含了大量的信息，处理起来相当繁琐，而且计算也非常复杂，容易出现错误，所以为了处理简单，必须要把彩色图片转化为蕴含信息比较少的灰度图片。这些都是要通过图像预处理来实现，图像预处理的目的是去除图像上的噪声，减少图像的信息，简化运算，为图像处理打好基础。图像预处理包括：图像滤波去噪，图像增强，图像二值化，边缘检测与图像检测。3.2.1.1图像滤波去噪 图像在采集、获取和传输的过程中，几乎所有的图像都会收到不同程度的噪声污染，这些噪声有的是可见的，有的是不可见的，有的噪声污染在图像上体现的很明显，有的则是难以察觉的。对于不同的地点，

24、信噪比是不同的，信噪比越大，图像的质量则越好，但当信噪比低于一定的数值时，噪声将会变成人眼可见的颗粒状，图片就会变得模糊起来，图片质量也会相应的下降。更重要的是，噪声所产生的颗粒有可能会掩盖我们需要的信息细节，由于在工厂中采集图像的过程中，噪声污染是不可避免的，而且噪声污染体现在视觉上也是比较明显的，所以必须要对图像进行滤波去噪的处理，因为摄像头所采集到的照片是彩色的，数据量很大，计算复杂，处理起来容易出错，而颜色信息对于最终的读数是没有影响的，所以，为了简化运算，使图像的处理变得先对简单，在滤波去噪之前，先要对图像进行灰度化处理，然后再进行滤波，在本设计中采用的是中值滤波，中值滤波10 于殿

25、泓. 图像检测与处理技术M. 西安: 西安电子科技大学出版社, 2006: 103-169.起初是用于一维信号的处理，后来很快被应用到二维图像的处理之中。所谓的中值滤波，10 于殿泓. 图像检测与处理技术M. 西安: 西安电子科技大学出版社, 2006: 103-169. 中值滤波是在很多地方都可以被用到，相对来说，原理也比较简单，是一种非线性的滤波，它采用的是一个包含奇数个点的滑动窗门，最终整个窗口的灰度值是窗口中的灰度值的中间值。它的原理如下，对于一个一维的序列选取的滑动窗口的长度为m，对这个一维的序列进行中值滤波，从输入的一维序列之中按照顺序取出m个元素，其中，i代表滑动窗口的中心位置，

26、=，将这m个数按数值由大到小排列起来，找出中间的数作为滤波的输出数值，数学表达式为： (3-1) 上式中，表示的是取序列的中值。例如，有这么一个一维的序列11，20，16，24，29，次序列由小到大排列的结果为11，16，20，24，29，那么，中值滤波输出的结果为20。假如灰度值是29的像素是被噪声污染的像素，那么，经过中值滤波后被噪声污染的像素将被过除掉。按照这个原理，中值滤波同样可以适用于二维的序列，在对二维序列进行中值滤波时，滤波的滑动窗口也将是二维的，把窗口内的像素按由小到大，或者由小到大进行排序，生成数据序列，二维中值滤波的结果如下（3-2）所示 (3-2) 一般情况下，二维中值滤

27、波的效果更加的明显。不同的滑动窗口可以产生不同的效果，滑动窗口的形状是多种多样的，滑动窗口的选择与所处理的图像和所要达到的目的有关，中值滤波是一种非线性的滤波，涉及到的计算也比较复杂，对于中值是0的正态分布随机噪声，中值滤波的输出噪声方差为： (3-3)上式中： 是输入的噪声功率； 是中值滤波器的窗口大小； 是输入噪声的均值； 是输入噪声的密度函数。那么，均值滤波的输出方差为： (3-4)对比上面的的几个公式可以看出，而均值滤波的输出和输入的分布是没有关系的，而中值滤波的输出的输入噪声的密度分布是相关的。中值滤波与均值滤波各有特点，均值滤波对于随机噪声的抑制能力更强，而中值滤波则更适用于脉冲干

28、扰。图3-7为加入盐噪声的图片和进行中值滤波 处理之后的图片，可以看出，中值滤波的效果还是很明显的。中值滤波的主要程序为：J= imread(1.jpg);K = imnoise(J,salt & pepper,0.02);imshow(J);figure;imshow(K);L =medfilt2(rgb2gray(K);imshow(K), figure, imshow(L); 图3-7 加入盐噪声的图片与滤波后的图片3.2.1.2图像增强 图像增强是以改善图像的视觉效果为根本目的的，在采集图像的过程中，特别是在工厂中，光照条件对图片质量的影响十分明显，而拍摄的过程中可能会出现曝光过度或者

29、曝光不足的现象，这样得到的图像将会是模糊的，表盘上的指针刻度区域与背景区域的区分不是很明显，不利于仪表图像的识别。为了使表盘上指针刻度区域与背景区域的界限变得明显，需要对图像进行空间域的图像增强。空间域的的图像增强主要有两种方式，直方图的均衡化处理和灰度变换，本设计中采用的是灰度变换。灰度变换的变换函数有三种，包括线性变换、非线性变换以及分段线性变换，三种变换函数各有特点，应根据图像的特点已经应用的场合选择相应的函数。本设计中采用的是分段线性变换法进行图像灰度的对比度增强。分段线性变换函数中最常被用到的是三段线性变换，其数学表达式如下: (3-5) 上式中是原图像的灰度值，是变换后的图像的灰度

30、值，是确定的分段区域的阈值。图3-8为增强后的图片，图像增强的主要程序为：j=imread(1.jpg);x=RGB2gray(j);subplot(1,2,1);imshow(x);title(原图像);f=double(x);m,n=size(f);h=fspecial(gaussian,25,25,80);q=imfilter(f,h,same);s=log(f+0.03)-log(q+0.03);r=exp(s);max_r=max(r(:)*0.27;min_r=min(r(:);r=(r-min_r)/(max_r-min_r);index=find(r1);r(index)=1;

31、R=mat2gray(r);subplot(1,2,2);imshow(R);title(处理后的图像);G=im2bw(R,0.7);imshow(G);I=uint8(G);bw=edge(I,sobel);imshow(bw); 图3-8 增强后的图片3.2.1.3二值化处理进行完图像增强后的图像是灰度图像，灰度图像的每个像素都可以在0-255之间取值，每个像素都可以有如此多的取值，这样的图像计算和处理起来十分的麻烦，而且容易出错，所以，就要对图像进行二值化处理，使图像的每个像素只能去0与1两个值，“0”作为目标区域，“1”作为背景区域，反之，亦可以，这样计算将变得相对简单，也不容易出错

32、。用来划分背景区域和目标区域的那个灰度值就称作为阈值，大于阈值的像素值全部设为1，小于阈值的则全部设为0.图像的二值化11 田村秀行(著), 金喜子, 乔双(译). 计算机图像处理M. 北京: 科学出版社,2004： 117-121.处理是图像处理中最常用的技术，在图像处理中具有至关重要的地位。主要有三个原因，首先，在彩色图像或者灰度图像的处理过程中，都要先将图像进行二值化处理，然后再进行更深层次的计算。其次，在理论上和方法上图像二值化比灰度图像处理更加系统化，也更加稳定。 第三，当图像处理的对象是文字或者图纸的时候，11 田村秀行(著), 金喜子, 乔双(译). 计算机图像处理M. 北京:

33、科学出版社,2004： 117-121. 1； 当时 (3-6) 0； 当时 通常情况下，在二值化图像中，数值是0的表示背景区域，数值是1的代表目标区域。而阈值t则是需要求解的，求解阈值t的方法称为阈值选择，常用的阈值选择方法有以下几种： 判别分析法。用阈值t将图像的灰度值直方图中的灰度值的分布分为两组，求出是两组分离达到最好的阈值t，这就是判别分析法。用判别分析法确定阈值t的方法如下所示，假设阈值是t，给定的图像具有L级灰度值，把大于t的灰度值的像素和小于t的像素分成两组，设为1组和2组，设i组 ，的像素数为，平均灰度值为，方差为，全体像素的平均值为，则组内方差为 (3-7)组间的方差为 (

34、3-8)如果设全部像素的灰度值的方差为，那么，则有 (3-9)因此可以得到判别比为 (3-10)由于全方差是一个常数，与阈值没有关系，由此可以得到，只需让最大就可使判别比取得最大值，所以只需使t的值不断变化并求出最大时对应的t的值。p-参数法。若要被分离的图像的面积大约等于，其与整个图像的面积S之比为p=。在选择阈值时，要使灰度直方图中大于t的像素对其他所有像素的比值亦为p为此，因此，要想求出阈值，只需要在图像的灰度直方图中，从灰度值高的一边开始计算，求出累积相对读数分布，当累积相对读数的值达到为p时的灰度值就是所求的阈值，这种方法也经常被用到。状态法。状态法只是用于灰度直方图呈现双峰型分布的

35、条件下，这时，阈值t就是两个峰间的低谷上的灰度值。在复杂的图像和瘦到很多干扰的图像中，直方图不能呈现明显的双峰型，所以，状态法有时并不适用。二值化后的图像如图3-9所示，图像二值化的主要程序为：T=0.5*(double(min(I(:)+double(max(I(:);done=false;while doneg=I=T;Tnext=0.5*(mean(I(g)+mean(I(g);done=abs(T-Tnext)=T);J(K)=255;K=find(J max_len) max_len = len; xy_long = xy; end end 图3-11 霍夫变换检测直线 图3-12

36、指针识别3.2.3读数识别 仪表表盘的读数识别是整个读数系统中最为重要的的部分，目前读数识别的方法多种多样，经过在工厂的实地考察，发现该工厂的有些材料生产车间环境相对较差，表盘难以时刻保持清洁，总是沾带有大量的粉尘，通过距离法读取表盘的读数相对来说比较困难，也不够稳定，而选取通过指针角度来读取表盘数据则可以更好的适应环境，提高读数的稳定性。读数识别的主要程序如下：plot(xy_long(:,1),xy_long(:,2),LineWidth,2,Color,cyan); k=(xy(2,2)-xy(1,2)/(xy(2,1)-xy(1,1); theta=pi/2+atan(k); if(x

37、y(1,1)+xy(2,1)/2 max_len) max_len = len; xy_long = xy; end end% highlight the longest line segment% plot(xy_long(:,1),xy_long(:,2),LineWidth,2,Color,cyan); k=(xy(2,2)-xy(1,2)/(xy(2,1)-xy(1,1); theta=pi/2+atan(k); if(xy(1,1)+xy(2,1)/21);r(index)=1;R=mat2gray(r);subplot(1,2,2);imshow(R);title(处理后的图像);

38、G=im2bw(R,0.7);imshow(G);I=uint8(G);bw=edge(I,sobel);imshow(bw);A=imread(1.jpg); I=rgb2gray(A); T=0.5*(double(min(I(:)+double(max(I(:);done=false;while done g=I=T;Tnext=0.5*(mean(I(g)+mean(I(g);done=abs(T-Tnext)=T);J(K)=255;K=find(J=L); r=0; s=0; for j=1:1:t r=r+LP(j); s=s+L(j)*LP(j); end W0(i)=r; W

39、1(i)=1-r; U0(i)=s/r; U1(i)=(m-s)/(1-r); end for i0=1:1:M BB(i0)=W0(i0)*W1(i0)*(U1(i0)-U0(i0)2); end for i=1:1:M if pbest(i,2)=gbest2 gbest2=MAX; gbest1=X(CC); end GG(k)=gbest2; for i=1:1:M V(i)=round(w(k)*V(i)+c1*rand*(pbest(i,1)-X(i)+c2*rand*(gbest1-X(i); X(i)=V(i)+X(i); end end for i=1:1:a for j=1

40、:1:b if J(i,j)gbest1 J(i,j)=250; else J(i,j)=0; end end end kk=1:1:G; gbest1; figure(1); imshow(J); figure(2); plot(kk,GG) tt=etime(clock,t1);endI=imread(1.jpg);tmin=min(I(:);tmax=max(I(:);th=(tmin+tmax)/2;ok=true;while ok g1=I=th; g2=Ith; u1=mean(I(g1); u2=mean(I(g2); tnew=(u1+u2)/2; if abs(th-tnew

41、)1 ok=0; endendth=tnew;th=floor(th);Inew=im2bw(I,th/255);subplot(1,2,1)imshow(I);title(原始图像);subplot(1,2,2)imshow(Inew);t=迭代法分割后的图像，阈值= num2str(th);title(t);A0=imread(1.jpg);seed=100,220;thresh=15;%相似性选择阈值A=rgb2gray(A0);A=imadjust(A,min(min(double(A)/255,max(max(double(A)/255,);A=double(A); B=A;r,c=

42、size(B);n=r*c;pixel_seed=A(seed(1),seed(2);q=seed(1) seed(2);top=1;M=zeros(r,c);M(seed(1),seed(2)=1;count=1; while top=0; r1=q(1,1); c1=q(1,2); p=A(r1,c1); dge=0; for i=-1:1 for j=-1:1 if r1+i0 & c1+j0 if abs(A(r1+i,c1+j)-p)=n top=1; end q=q(2:top,:); top=top-1;endsubplot(1,2,1),imshow(A,);subplot(1,2,2),imshow(B,