系统仿真学报论文模板下载

1.文章标题：一般不超过20个汉字，必要时加副标题，并译成英文。2.作者姓名、工作单位：题目下面均应写作者姓名，姓名下面写单位名称（一、二级单位）、所在城市（不是省会的城市前必须加省名）、邮编，不同单位的多位作者应以序号分别列出上述信息。3.提要：用第三人称写法，不以“本文”、“作者”等作主语，应是一篇能客观反映文章核心观点和创新观点的表意明确、实在的小短文，切忌写成背景交代或“中心思想”，100－200字为宜。4.关键词：3－5个，以分号相隔，选择与文章核心内容相关的具有独立性的实在词。5.正文标题：内容应简洁、明了，层次不宜过多，层次序号为一、（一）、1、（1），层次少时可依次选序号。6.正文文字：一般不超过1万字，正文用小4号宋体，通栏排版。7.数字用法：执行GB/T15835－1995《出版物上数字用法的规定》，凡公元纪年、年代、年、月、日、时刻、各种记数与计量等均采用阿拉伯数字；夏历、清代及其以前纪年、星期几、数字作为语素构成的定型词、词组、惯用语、缩略语、临近两数字并列连用的概略语等用汉字数字。8.图表：文中尽量少用图表，必须使用时，应简洁、明了，少占篇幅，图表均采用黑色线条，分别用阿拉伯数字顺序编号，应有简明表题（表上）、图题（图下），表中数据应注明资料来源。9.注释：注释主要包括释义性注释和引文注释，集中列于文末参考文献之前。释义性注释是对论着正文中某一特定内容的进一步解释或补充说明；引文注释包括各种不宜列入文后参考文献的引文和个别文后参考文献的节略形式，其序号为①②③……10.参考文献：参考文献是作者撰写论着时所引用的已公开发表的文献书目，是对引文作者、作品、出处、版本等情况的说明，文中用序号标出，详细引文情况按顺序排列文尾。以单字母方式标识以下各种参考文献类型：普通图书[M]，会议论文[C]，报纸文章[N]，期刊文章[J]，学位论文[D]，报告[R]，标准[S]，专利〔P〕，汇编[G]，档案[B]，古籍[O]，参考工具[K]，其他未说明文献〔Z〕。格式与示例如下：（1）图书类格式：[序号]主要责任者。文献题名：其他题名信息（任选）[文献类型标识].其他责任者（任选）。出版地：出版者，出版年：起止页码。（2）期刊文章格式：[序号]主要责任者。文献题名[J].刊名（建议外文刊名后加ISSN号），年，卷（期）：起止页码。（3）报纸文章格式：[序号]主要责任者。文献题名[N].报纸名，出版日期（版次）。（4）古籍格式：[序号]主要责任者。文献题名[O].其他责任者（包括校、勘、注、批等）。刊行年代（古历纪年）及刊物机构（版本）。收藏机构。（5）析出文献格式：[序号]析出文献主要责任者。析出文献题名[文献类型标识]//原文献主要责任者（任选）。原文献题名。出版地：出版者，出版年：析出文献起止页码。（6）电子文献格式：主要责任者。文献题名[文献类型标识/载体类型标识].出处或可获得地址，发表或更新日期/引用日期（任选）。（7）文献重复引用标记：同一作者的同一文献被多次引用时，在文后参考文献中只出现一次，其中不注页码；而在正文中标注首次引用的文献序号，并在序号的角标外着录引文页码。11.基金项目：获得国家基金资助和省部级科研项目的文章请注明基金项目名称及编号，按项目证明文字材料标示清楚。12.作者简介：姓名（出生年-），性别，民族（汉族可省略），籍贯，现供职单位全称及职称、学位，研究方向。13.来稿请注明作者电话、E-mail，收刊人及详细地址、邮编。14.其他：请勿一稿两发，并请自留原稿，本刊概不退稿。

[1] Xiangyun Liao, Zhiyong Yuan*, Qi Zheng, Qian Yin, Dong Zhang, Jianhui Zhao. Multi-Scale and Shape Constrained Localized Region-Based Active Contour Segmentation of Uterine Fibroid Ultrasound Images in HIFU Therapy, PLoS One, Vol. 9, No. 7, e103334, 2014.(SCI收录)[2] Xiangyun Liao, Zhiyong Yuan*, Pengfei Hu and Qianfeng Lai. GPU-assisted energy asynchronous diffusion parallel computing model for soft tissue deformation. Simulation: Transactions of the Society for Modeling and Simulation International, Vol. 90, , 1199-1208, 2014. (SCI收录)[3] Weixin Si, Zhiyong Yuan*, Xiangyun Liao, Zhaoliang Duan, Jianhui Zhao. An Energy Based Free Boundary Asynchronous Diffusion Model for 3D Warping of Tissue Dynamics. Journal of Statistical Computation and Simulation, Vol. 84, No. 6, Pages: 1280-1296, 2014. (SCI收录)[4]袁志勇，童倩倩，喻思娇, 廖祥云. 基于压缩感知与LS_SVM的三维组织表面重建. 华中科技大学学报(自然科学版), 第42卷, 第8期, Pages: 7-11, 2014. (EI收录)[5]Zhiyong Yuan, Weixin Si, Xiangyun Liao, Zhaoliang Duan, Yihua Ding, Jianhui Zhao. Parallel computing of 3D smoking simulation based on OpenCL heterogeneous platform. Journal of Supercomputing, 2012, Volume 61, Number 1, Pages: 84-102. (SCI&EI收录)[6]Jianhui Zhao, Yihua Ding, Ravindra S. Goonetilleke, Shuping Xiong, Yuanyuan Zhang, Chengjiang Long, Zhiyong Yuan. Interactive Deformation Simulation of Manual Girth Measurement for Limbs, Information-An International Interdisciplinary Journal, 2012, , , Pages: 339-346. (SCI收录[7]Weixin Si, Zhiyong Yuan*, Xiangyun Liao, Zhaoliang Duan, Yihua Ding and Jianhui Zhao. 3D Soft Tissue Warping Dynamics Simulation Based on Force Asynchronous Diffusion Model. CASA2011 Special Issue of Journal of Computer Animation and Virtual Worlds, 2011, Volume 22, Issue 2-3: 251-259. (SCI收录)[8]Jianhui Zhao, Zhong Zhang, Shizhong Han, Chengzhang Qu, Zhiyong Yuan, Dengyi Zhang. SVM Based Forest Fire Detection Using Static and Dynamic Features. Computer Science and Information Systems, Vol. 8, No. 3, 2011, Pages: 821-841. (SCI收录)[9]Zhiyong Yuan, Shikun Feng, Qian Yin, Xiali Wang, Dengyi Zhang, Jianhui Zhao, Mianyun Chen. Endoscopic Image Cutting Simulation Based on Mass-Spring Model and Computational Geometry. International Journal of Circuits, Systems, and Computers,(2009), Pages: 1453-1465. (SCI收录)[10]Qian Yin, Zhi-Yong Yuan, Ping Guo. Studies on the distribution of the shortest linear recurring sequences. Information Sciences, 2009, , Pages: 2379-2389. (SCI收录)[11]Mengyi Liu, Zhiyong Yuan*, Yufeng Ma, Xingwei Chen, Qian Yin. Heterogeneous Face Recognition and Synthesis Using Canonical Correlation Analysis. Journal of Convergence Information Technology, 2012, , , Pages: 398-407. (EI)[12]Zhaoliang Duan, Zhiyong Yuan*, Xiangyun Liao, Weixin Si, Jianhui Zhao. 3D Tracking and Positioning of Surgical Instruments in Virtual Surgery Simulation. Journal of Multimedia, 2011, Volume 6, Number 6, Pages: 502-509.(EI收录)[13]袁志勇, 司伟鑫, 廖祥云, 段昭亮, 赵俭辉. 虚拟针灸视觉建模与仿真方法研究. 系统仿真学报, 2011年第23卷第8期, Pages: 1708-1713.[14]Zhiyong Yuan, Yuanyuan Zhang, Jianhui Zhao, Yihua Ding, Chengjiang Long,Lu Xiong, Dengyi Zhang, Guozhong Liang. Real-time Simulation for 3D Tissue Deformation with CUDA Based GPU Computing. Journal of Convergence Information Technology, Volume 5, Number 4, June 2010, Pages:109-119. (EI收录)[15]Zhiyong Yuan, Zhaoliang Duan, Weixin Si, Xiangyun Liao, Ruixue Mao, Jianhui Zhao. Real-time smoking simulation based on OpenCL and GPU acceleration technologies. Journal of Information and Computational Science, , , Pages: 2503-2512, December, 2010. (EI收录)[16]Yihua Ding, Jianhui Zhao, Zhiyong Yuan, Yuanyuan Zhang, Chengjiang Long, Lu Xiong. Constrained Surface Recovery Using RBF and Its Efficiency Improvements. Journal of Multimedia, 2010, Vol. 5, , Page(s): 55-62. (EI收录)[17]Zhiyong Yuan, Yuanyuan Zhang, Jianhui Zhao, Yihua Ding, Chengjiang Long,Lu Xiong, Dengyi Zhang, Guozhong Liang. Real-time Simulation for 3D Tissue Deformation with CUDA Based GPU Computing. Journal of Convergence Information Technology, Volume 5, Number 4, June 2010, Pages:109-119. (EI收录)[18]Weixin Si, Zhiyong Yuan*, Ruixue Mao, Jianhui Zhao, Xiangyun Liao, Zhaoliang Duan. Modeling and Realization of Virtual Acupuncture Training Simulation System. Journal of Digital Content Technology and its Applications, , , November 2010, Pages: 126-136. (EI收录)[19]Xiangyun Liao, Zhiyong Yuan*, Weixin Si, Zhaoliang Duan, Ruixue Mao, Jianhui Zhao. Research and Application of Parallel Computing Technologies based on CUDA and OpenCL. Journal of Convergence Information Technology, Volume 6, Number 6, June 2011, Pages: 25-37.(EI收录)[20]Yuanyuan Zhang, Zhiyong Yuan*, Yihua Ding, Jianhui Zhao. Real-time Simulation of Smoking Based on GPU Computing for Surgical Simulation. Journal of Computational Information Systems, Volume 6, Number 13, December 2010, Pages: 4263-4273. (EI收录)[21]Zhiyong Yuan, Yihua Ding, Yuanyuan Zhang, Jianhui Zhao. Real-time Simulation of Tissue Cutting with CUDA Based on GPGPU. Advanced Materials Research, (2010), Pages: 154-161. (EI收录)[22]陈二虎, 袁志勇*, 余泽江, 廖祥云, 朱炜煦. 虚拟手术系统中电磁力反馈建模仿真与实现. 系统仿真学报, 第26卷, 第9期, Pages: 2003-2008, 2014.[23]赖虔葑, 袁志勇*, 郭甲翔, 廖祥云, 郑绵仑. 基于三维均值坐标和GPU加速的弹性体形变仿真. 系统仿真学报, 第26卷, 第9期, Pages: 1907-1926, 2014.[24]胡鹏飞, 袁志勇*, 廖祥云, 郑奇, 陈二虎. 基于CPU-GPU混合加速的SPH流体仿真方法. 计算机工程与科学, 第36卷, 第7期, Pages: 1231-1237, 2014.[25]郑奇, 袁志勇*, 肖玲, 王慧玲, 王高华. 基于Cortex-M3的高速安全U盘设计. 计算机工程与设计, 第35卷, 第4期, Pages: 1214-1220, 2014.[26]袁志勇*, 郭甲翔, 杨成, 韩婕, 廖祥云. 基于粒子的流体和可形变固体双向耦合. 系统仿真学报, 第25卷，第10期, Pages: 2327-2331, 2013.[27]袁志勇, 尹乾, 胡君, 冯仕坤. 基于改进质点-弹簧模型的图像变形仿真方法.华中科技大学学报自然科学版, 2009年第37卷, 第1期, Pages: 35-37. (EI收录)[28]袁志勇, 尹乾, 冯仕坤, 胡君. 内窥镜图像实时动态流血合成与仿真. 武汉大学学报信息科学版, 2009年第34卷, 第3期, Pages: 317-320. (EI收录)[29]袁志勇, 丁乙华, 张圆圆, 赵俭辉. 基于GPU加速的虚拟内窥镜手术实时冒烟仿真. 武汉理工大学学报, 2009年第31卷第18期, Pages: 55-58.[30]赵俭辉, 龙成江, 丁乙华, 袁志勇. 一种基于立方体小栅格的K邻域快速搜索新算法. 武汉大学学报(信息科学版), 2009年第34卷,第5期, pages: 615-618.(EI收录)[31]袁志勇, 熊惠林, 陈绵云. IKE协议的研究与改进. 计算机工程, 2007年第9期.[32]袁志勇, 查桂峰, 陈绵云. 基于反对称小波的车型识别研究. 武汉大学学报(信息科学版), 2005年第30卷, 第6期, pages: 560-563. (EI收录)[33]袁志勇, 查桂峰, 陈绵云. 基于聚类的二级模糊综合评判的车型识别研究. 计算机工程与应用, 2005年第12期, Pages: 202-205.[34]袁志勇, 肖畅, 陈绵云. 一种NVF量化的小波域图像水印算法及实现. 计算机工程与科学, 2005年第8期, Pages: 44-46.[35]袁志勇, 杨土安, 夏维, 陈绵云. 一种自适应目标DCT的数字水印算法. 计算机工程与科学, 2005年第10期, Pages: 40-43.[36]袁志勇, 李敏, 陈绵云. 基于整型小波变换和嵌入式零树编码的医学图像压缩. 计算机应用研究, 2004年第6期, Pages: 110-113.[37]谭晓华, 袁志勇. XML数据存储模型及优化的研究. 微型机与应用, 2004年第1期[38]Zhiyong Yuan, Chang Xiao, Mianyun Chen. Image Watermarking Based on Wavelet Quantization. Journal of Advances in Systems Science and Applications, Pennsylvania, USA, Dec. 2003, Vol. 3, , Pages: 641-644.[39]李敏, 苏光奎, 袁志勇. 生物组织脱水机智能系统设计与实现. 微型机与应用, 2003年第12期.[40]苏光奎, 袁志勇, 杨维均. 远程智能数据采集与处理系统设计. 武汉大学学报(信息科学版), 2002年第2期.[41]袁志勇, 苏光奎, 杨维均. 利用Win32 API通信函数设计串行通信程序. 测绘信息工程, 2000年第4期.[42]袁志勇. 实用微机红外通讯技术, 实用无线电, 1998年第5期.

1. 重庆工商大学实验实习中心,重庆400033 ; 2. 重庆工商大学图书馆,重庆400033)摘 要:介绍了MATLAB 图像处理工具箱中的函数,给出了图像处理与分析的技术实现,如用直方图均衡进行图像增强,通过形态学方法进行图像特征抽取与分析,借助于分水岭图像分割实现目标检测等。关键词:灰度图像;形态学变换;标记;分割;特征抽取中图分类号:TP 317. 4 文献标识码:AMATLAB6. 1(R12. 1) 是一套功能十分强大的工程计算及数据分析软件,其应用范围涵盖了数学、工业技术、电子科学、医疗卫生、建筑、金融、数字图像处理等各个领域。许多工程师和研究人员发现,MATLAB能迅速测试其构思,综合评测系统性能,并能借此快速设计出更多的解决方案,达到更高的技术要求。MATLAB 的图像处理工具箱,功能十分强大,支持的图像文件格式丰富,如3 . BMP , 3 . JPG, 3 . JPEG,3 . GIF , 3 . TIF , 3 . TIFF , 3 . PNG, 3 . PCX , 3 . XWD , 3 . HDF , 3 . ICO , 3 . CUR 等。利用MATLAB 所提理函数,并给出用MATLAB 实现图像处理与分析的应用技术实例。1 MATLAB 的图像处理工具概述MATLAB6. 1(R12. 1) 提供了20 类图像处理函数,涵盖了图像处理的包括近期研究成果在内的几乎所有的技术方法,是学习和研究图像处理的人员难得的宝贵资料和加工工具箱。这些函数按其功能可分为:图像显示;图像文件I/ O ;图像算术运算;几何变换;图像登记;像素值与统计;图像分析;图像增强;线性滤波;线性二元滤波设计;图像去模糊;图像变换;邻域与块处理; 灰度与二值图像的形态学运算;结构元素创建与处理;基于边缘的处理; 色彩映射表操作;色彩空间变换;图像类型与类型转换。2 应用MATLAB 工具箱进行图像分析处理2. 1 用直方图均衡实现图像增强当图像对比度较低,即灰度直方图分布区间较窄时,可用直方图均衡实现灰度分布区间展宽而达到图像增强的效果。下面是实现的源程序及相关功能的注解:%源程序:test1. mX 收稿日期:2003 - 02 - 27 ;修回日期:2003 - 03 - 30作者简介:何希平(1968 - ) ,男,四川人,博士生,重庆工商大学副教授,从事多媒体数据压缩、网络信息系统研究。. 1995-2005 Tsinghua Tongfang Optical Disc Co., Ltd. All rights , close all %清除所有内存变量、图形窗口I = imread(’pout. tif’); %将图像文件pout. tif 的图像像素数据读入矩阵Iimshow( I) %显示图像I ,图像对比度低,如图1afigure , imhist ( I) %在新图形窗口中显示图像I 的直方图,如图1c。可以注意到图像%亮度范围相当狭窄,并未完全覆盖可能的范围[0 ,255 ]I2 = histeq( I) ; %对图像I 做直方图均衡补偿在整个范围内展宽亮度值并输出到矩阵I2 ,因而改进了图像I 的对比度figure , imshow( I2) %在新图形窗口中显示新图像I2 , 如图1bfigure , imhist ( I2) %在新图形窗口中显示图像I2 的直方图, 如图1dimwrite ( I2 , ’pout2. png’); %将对比度调节的结果图像写入PNG格式的文件a 原图 b 直方图均衡结果图 c 原图像的直方图 d 结果图像的直方图图1 直方图均衡补偿消去图像噪声程序运行后,可得如图1 的对比图像。2. 2 用形态学方法进行图像处理与分析以rice. tif 为图像实例,介绍用形态学方法对灰度图像进行处理与分析的技术要点,即对灰度图像进行如下操作:去除图像的不均匀背景;用设置阈值的方法(thresholding) 将结果图像转换成二值图像;通过成分标记(components labeling) 返回图像中的目标对象属性,并计算目标对象的统计数字特征。其算法步骤描述如下:(1) 用工具箱函数imread 和imshow 读取和显示8 位灰度图,如图2a 。(2) 用形态学开运算(Morphological Opening) 估计背景。通过调用imopen 并对输入图像I 执行形态学开运算, 取半径为15 的圆盘结构元素,且结构元素通过函数strel 建立。形态学开运算有消除不能完全包含在半径为15 的圆盘内的目标对象的作用。注意到图像(如图2b) 中央的背景照度(background illumina2tion) 比底部要亮。(3) 用surf 指令察看背景图像。用Surf 指令创建近似背景的彩色表面图(如图2c) ,使人可以看到在一个矩形区域上的数学函数特征。在表面图中,[0 , 0 ] 表示原点, 或图像左上角,曲面图最高部分表示背景的最亮像素(从而rice. tif 的背景的最亮像素出现在图像中央行的附近,而最暗像素出现在图像的底部) 。(4) 从原图像中减去背景图像。须用图像处理工具箱的图像算术函数imsubtract 产生均匀的背景(如图2d) 。(5) 调节图像对比度。用imadjust 指令增大图像对比度(如图2e) 。imadjust 函数需要一个输入图像且也可带两个矢量: [ low high ] 和[ bottom top ] . 输出图像通过将输入图像中low 值映射到输出图像中的bot2tom 值、high 值映射到输出图像中的top 值,并将low 与high 间的值进行线性缩放而产生。(6) 对图像进行阈值处理。先调用graythresh ,自动计算一个适当的阈值;然后使用graythresh 返回的阈值,调用im2bw 执行阈值处理,将灰度图像转换成二值图像(如图2f) 。(7) 确定图像中的目标对象并予以标记。调用bwlabel 寻找连通成份而且用惟一的数字将他们分类标记。bwlabel 接受一个二值图像和指定各目标对象的连通性的值(4 或8 ,表示4 或8 连通) 作为输入。注意: 结果的准确性依赖于许多因素,包括: 目标对象的大小; 近似背景的准确程度; 是否设定连接3 2 第2 期 何希平等: 基于MATLAB 的图像处理与分析. 1995-2005 Tsinghua Tongfang Optical Disc Co., Ltd. All rights 原图b 背景性参数为4 或8 ; 是否任何目标对象均相接(在这种情况下他们可能被标记为同一目标对象) ;在该实例中, rice 的一些谷粒正好相接,因此bwlabel 把它们视为了同一目标对象。(8) 查看标记矩阵。看一看bwlabel 产生的标记的近似形状是有用的。调用imcrop 并用鼠标选择包含某一目标对象的一部分及其一些背景的部分图,则所选部分图的像素值会在MAT2LAB 窗口中返回。若查看上面的结果,你会看到c 背景的表面图d 原图与背景的差一个对象的一角标以某数字标记k ,这意味着它是第k 个被bwlabel 分类的目标对象。imcrop 函数也可带矢量指定剪裁矩形的坐标。在这种情况下,它不执行交互式操作。举例来说,rect = [15 25 10 10 ] ; roi = imcrop (labeled ,rect)这个调用指定一个剪裁矩形的左上角坐标始于(15 ,25) ,而且高度和宽度均为10 。一种查看标记矩阵的好方法是将它显示成e 图像对比度调节结果 f 阈值处理后的二值图一种假彩色索引图像(如图2g) 。在假彩色索引图像中,将标记矩阵中区分每一对象的数字映射成了相关色彩映射矩阵中的一种不同的颜色。当把一个标记矩阵看成一个RGB 图像时,图像中的对象是比较容易区别的。为此, 使用la2bel2rgb 函数。使用该函数时,可以指定色彩映射表,背景颜色,以及标记矩阵中的对象如何映射为色彩映射表中的颜色。(9) 测量图像中的对象属性。regionprops 指令可测量图像中的对象或区域的属性,并返回一g 假彩色标记图h 谷粒大小分布图图2 形态学图像处理的对比分析结果个结构数组。当将其作用于一个图像成分的标记矩阵时,它为每个成分建立一个结构元素,而每一结构元素包含一个标记成分的一些基本属性。regionprops 函数支持对许多不同的属性予以测量, 但是设定属性参数为’basic’旨在返回最常用的三个量: 面积(Area) , 质心或块中心(Centroid) 和边框(BoundingBox) 。边框Bounding2Box 表示能容纳一个区域(所举实例中的谷粒)的最小长方形, 为四元素矢量: [ left top widthheight ] 。(10) 在图像中计算目标对象的统计特性。使用MATLAB 函数max , mean , 和hist 可计算被阈值处理的目标对象的一些统计属性(如图2h) 。图像处理工具箱也有一些统计函数,如mean2 和std2 ,适用于图像数据,因为他们对二维空间的数据返回单一值。下面是算法实现的程序代码:%程序代码:test2. mclear , close all ,I = imread(’rice. tif’); imshow( I) %读取和显示8 位灰度图rice. tif4 2 重庆工商大学学报 (自然科学版) 第20 卷. 1995-2005 Tsinghua Tongfang Optical Disc Co., Ltd. All rights = imopen( I ,strel (’disk’,15) ) ; %取半径为15 的圆盘结构元素对图像用开运算估计背景figure ,imshow(background) %显示背景图figure , surf (double (background (1 :8 : end ,1 :8 : end) ) ) ,zlim( [ 0 255 ]) ; %显示背景的彩色表面图,对8 ×8 格点取样set(gca ,’ydir’,’reverse’);I2 = imsubtract ( I ,background) ; figure , imshow( I2) %从原图像中减去背景图像,并显示结果图I3 = imadjust ( I2 , stretchlim( I2) ,[0 1 ]) ;figure , imshow( I3) ; %调节图像对比度,并显示结果图level = graythresh( I3) ;bw = im2bw( I3 ,level) ; figure , imshow(bw) %将灰度图像转换成二值图像[ labeled ,numObjects ] = bwlabel(bw ,4) ; % 成分标记,4 具体指定4 - 连通成分.grain = imcrop (labeled) % 用鼠标选取实现交互式剪裁标记成分的一部分RGB-label = label2rgb(labeled , @spring , ’c’, ’shuffle’); %把一个标记矩阵转换成一个RGB 图像figure ,imshow(RGB-label) ;graindata = regionprops(labeled ,’basic’) %调用regionprops ,为rice 的每一经阈值处理%的谷粒返回一个基本属性的结构。由BoundingBox 的域返回四元素矢量: [ left top width height ]。graindata (51) . Area , graindata(51) .BoundingBox , graindata(51) . Centroidallgrains = [graindata. Area ] %用点号存取graindata 的所有元素的面积域并将该数据存入%新的矢量allgrains。这个步骤简化了对面积量的分析,因为不必使用域名存取面积。max(allgrains) %找最大谷粒的大小。allgrains 中的数据是一维的, 故函数mean 和std 是适用的。biggrain = find(allgrains = = ans) %使用find 指令返回该最大谷粒的成分标记mean(allgrains) %求平均粒径hist (allgrains ,20) %作包含20 个方柱的显示谷粒大小分布的直方图。直方图表明,在rice 图像中谷粒最通常的%大小在300 到400 个像素的范围内(如图2h) 。2. 3 用分水岭分割法检测连通目标在一个图像中检测目标是图像分割的一个例子。为分割连通目标,时常用Watershed 变换。如果把一幅图像看做一个具有山(高亮度) 和低谷(低亮度) 的表面,那么这个变换在一幅图像中找亮度低谷。实现包括下列步骤:(1) 读图像。读入图像afmsurf . tif , 它是一幅原子能显微镜下的衣料表面图像(如图3a) 。(2) 对比度最大化。注意到图像中有许多彼此连通的不同大小的对象。为使通过watershed 变换找到的低谷数目最小,我们使感兴趣的对象的对比度达到最大。对比度增强的一个常用的技术是综合应用top- hat 和bottom - hat 变换。top - hat 变换定义为原图像和它的开之差。图像的开是一与特定结构元素匹配的图像前景部分的集合(如图3b) 。bottom - hat 变换定义为在原图像和它的闭之间的差。图像的闭是一与特定结构元素匹配的图像背景的集合(如图3c) 。通用的结构元素是正方形,长方形,圆盘,菱形,球和线。既然图像中我们感兴趣的目标对象看起来像圆盘,我们用strel 函数建立一个半径为15 个像素的圆盘形结构元素。这个圆盘尺度是图像中的目标对象的平均半径的一个估计。(3) 图像相加减。看到top - hat 图像含有与结构元素匹配的对象的”巅峰”。相反,bottom - hat 图像显示出感兴趣的目标对象之间的间隙。为使目标对象与分隔它们的间隙之间的对比达到最大,用“原图+ top - hat 图像- bottom - hat 图像”得到增强的结果图(如图3d) 。(4) 转换感兴趣的对象。调用watershed 变换找出图像的亮度”低谷”,把imcomplement 函数作用于增强过的图像上,将感兴趣的目标对象转换为亮度低谷,得到增强图的补图(如图3e) 。(5) 检测亮度低谷。对所得补图运用imextendedmin 函数检测低于某特别阈值的所有亮度低谷。imextendedmin 函数的输出是一个二值(逻辑值) 图像(如图3f) 。二值图像中重要的是区域的位置而非区域的大小。用imimposemin 函数把补图改为只含有那些由imextendedmin 函数找到的低谷,并将低谷的像素值变为0 (8 位图像可能的深谷) (如图3g) 。(6) Watershed 分割。通过watershed 变换,可找出来所有含有强加给最小值的区域。用watershed 函数实现Watershed 分割。watershed 函数返回一个标记矩阵,它含有对应于watershed 区域的非负数。凡未落入5 2 第2 期 何希平等: 基于MATLAB 的图像处理与分析. 1995-2005 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved.图3 用Watershed 分割法检测连通目标的图像渐近过程任何watershed 区域的像素均被赋予像素值0。用label2rgb把一个标记矩阵变为一幅图像(如图3h) 。(7) 从标记矩阵中抽取目标对象的特征。可用region2props 函数从标记矩阵中抽取特征。比如说,可以计算两个量(面积和方向) 并把他们看成彼此的一个函数。下面给出算法的实现代码:%程序代码:test3. mafm = imread (’afmsurf . tif’); figure , imshow(afm) , title (’surface im2age’);se = strel (’disk’, 15) ;Itop = imtophat (afm , se) ; figure , imshow( Itop , [ ]) , title (’top - hat im2age’);Ibot = imbothat (afm , se) ; figure , imshow( Ibot , [ ]) , title (’bottom - hatimage’);Ienhance = imsubtract (imadd ( Itop , afm) , Ibot) ; figure , imshow( Ien2hance) , title (’original + top - hat - bottom - hat’);Iec = imcomplement ( Ienhance) ; figure , imshow( Iec) , title (’complementof enhanced image’);Iemin = imextendedmin( Iec , 22) ; figure , imshow( Iemin) , title (’extend2ed minima image’);Iimpose = imimposemin ( Iec , Iemin) ; figure , imshow( Iimpose) , title ( ’imposed minima image’);wat = watershed( Iimpose) ;rgb = label2rgb(wat) ; figure , imshow(rgb) ;title (’watershed segmented image’);stats = regionprops (wat , ’Area’, ’Orientation’); area = [ stats ( :) .Area ] ; orient = [ stats( :) . Orientation] ;figure , plot (area , orient , ’b 3 ’); title (’Relationship of Particle Orienta2tion to Area’);xlabel (’particle area (pixels) ’); ylabel (’particle orientation (degrees) ’);参考文献:[1 ] 孙兆林.MATLAB 6. x 图像处理[M] . 北京:清华大学出版社,2002[2 ] 崔屹. 图像处理与分析———数学形态学方法及其应用[M] . 北京:科学出版社,2000[3 ] 张远鹏,董海,周文灵. 计算机图像处理技术基础[M] . 北京:北京大学出版社,1996Image processing and analysis based on MATLABHE Xi - ping1 , ZHANG Qiong - hua2(1. Center of Experiment and Practice ,ChongQing Technology and Business University ,ChongQing 400033 ,China ;2. Library , ChongQing Technology and Business University , ChongQing 400033 ,China)Abstract :This paper first introduces the functions of MATLAB image processing toolbox , then presents sometechniques in image processing and analysis , such as image enhancement by using histogram equalization , image fea2ture extracting and analysis with morphological methods , and objects detection through watershed image words : grayscale intensity image ; morphological transform; labeling ; segmentation ; feature extraction责任编辑:杨祖彬6 2 重庆工商大学学报 (自然科学版) 第20 卷. 1995-2005 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved.收稿日期:2002208224MATLAB 及其在图像处理中的应用许志影,李晋平(中国矿业大学资源学院,江苏徐州 221008)摘要: 介绍了MATLAB 的特点和功能,分析了MATLAB 在图像处理中的应用,并结合实例说明了MATLAB 在图像处理中关键词: MATLAB ; 图像处理; 边缘提取中图分类号: TN911. 73 文献标识码: AMATLAB and Its Application to Digital Image ProcessingXU Zhi2ying ,LI Jin2ping(School of Resource & Geoscience ,China University of Mining and Technology ,Xuzhou 221008 ,China)Abstract :Introduces properties and functions of MATLAB ,and analyses its applications to digital image processing ,finally ,displays the a2bility of MATLAB in image processing with an :MATLAB ;image processing ;edge detectionMATLAB 软件由美国Math Works 公司于1984 年推出,历经十几年的发展和竞争,现已成为( IEEE) 国际公认的最优秀的科技应用软件之一。作为一个跨平台的软件,MATLAB 已推出Unix、Windows 9x/ NT、Linux 和Mac 等十多种操作系统平台下的版本,大大方便了在不同操作系统平台下的研究工作。目前基于Windows 系统的最新版本已上升到MATLAB6. 5 ,它继承了以往版本的优点,非常容易使用。现在,MATLAB 已经发展成为一个系列产品:MATLAB 主包和各种工具箱(TOOLBOX) 。目前已经推出了30 多个工具箱,这些工具箱可分为两大类:功能性工具箱和学科性工具箱。功能性工具箱主要用来扩充其符号计算功能、图示建模仿真功能、文字处理功能以及硬件实时交互功能,能用于多种学科。而学科性工具箱是专业比较强的,如控制工具箱、信号处理工具箱、图像处理工具箱和小波工具箱等多个学科的专业工具箱。借助于这些工具箱,各个层次的研究人员就可方便地进行研究工作,提高工作效率。本文将简要介绍MATLAB6. 5 及其在图像处理中的应用,希望对从事图像处理工作的研究人员有所帮助。1 MATLAB 概述MATLAB 最初是作为矩阵实验室(Matrix Labora2tory) 用来提供通往LINPACK和EISPACK矩阵软件包接口的。后来,它逐渐发展成为通用科技计算和图视交互系统的程序语言,其数据的基本单元是矩阵。它的指令表达与数学、工程中常用的习惯形式十分相似,从而使许多用C 或Fortran 实现起来十分复杂和费时的问题用MATLAB 就可以轻松地解决。MAT2LAB 的典型应用包括:数学计算、算法研究、数据分析和计算结果可视化、建模与仿真等。1. 1 MATLAB的特点MATLAB 有三大特点:一是功能强大。主要包括数值计算和符号计算、计算结果和编程可视化、数学和文字统一处理、离线和在线计算。二是界面友好,编程效率高。MATLAB 是一种以矩阵为基本单元的可视化程序设计语言,语法结构简单,数据类型单一,指令表达与标准教科书的数学表达式相近。三是开放性强。MATLAB 有很好的可扩充性,可以把它当成一种更高级的语言去使用。使用它很容易编写各种通用或专用应用程序。1. 2 MATLAB的主要功能MATLAB 之所以成为世界顶级的科学计算与数学应用软件,是因为它随着版本的升级与不断完善而具有愈来愈强大的功能。