[收稿日期]20050228 [作者简介]黄书先(1962),女,1983年大学毕业,硕士,副教授,现主要从事石油勘探开发研究与科研管理工作。
基于MATLAB 平台的可视化图像处理系统设计方法黄书先 (长江大学科学技术处,湖北荆州434023)[摘要]以MATLAB 为工作语言和开发环境,开发了一个在M AT LAB 平台下的可视化图像处理系统,可实现包括对测井图像在内的一般图像的精细处理,并能和用户开发的程序接口。
为M AT LAB 的再开发和可视化系统的设计作了有益的探索。
[关键词]MATLAB;图像处理;可视化GUI 界面[中图分类号]TP 311111[文献标识码]A [文章编号]16731409(2005)04015803MAT LAB 的图像处理工具为自然科学各学科领域的学者、研究人员和工程师提供了一个直观的灵活的环境,用以解决复杂的图像处理问题。
用MAT LAB 语言开发的图像处理算法可以在所有支持MAT LAB 的平台上共享。
也可以将m 语言算法和现存的C 程序集成在一起或者将MATLAB 开发的m 语言算法和GU Is 编译为C/C++代码,供其他程序调用,或者发布为一个独立的应用程序。
下面笔者提出一个基于MATLAB 平台的可视化图像处理系统设计方法,可用于包括测井图像资料在内的一般图像的精细处理。
1 系统总体设计在MAT LAB 中有个重要的图像处理工具包[1],该工具包是由一系列支持图像处理操作的函数组成的,按功能可以分为以下几类:图像显示;图像文件输入与输出;几何操作;像素值统计;图像增强;图像识别;图像滤波;图像变换;邻域和块操作;二值图像操作;颜色映射和颜色空间转换;图像格式转换等。
和其他工具包一样,用户还可以根据需要书写自己的函数,以满足特定的需要;也可以将这个工具包和信号处理工具包或小波工具包等其他工具包联合起来使用。
MAT LAB 提供了交互式的GU I 开发环境[2~4],用户只需要设置各个对象相应的属性,系统自动生成与之对应的界面,大大减少了开发的难度。
本设计的MAT LAB 图像处理系统由封面、主界面和各个子功能界面组成,其框图如图1所示。
程序总流程图如图2所示。
2 封面界面设计首先用MAT LAB 编辑封面的脚本文件,生成系统的封面界面,然后再进入处理的主界面。
在设计封面时,要做好封面的总体布局,力求完美。
要插入背景,首先要找到所用的函数(这在一般的书中很难找到);其次要注意设置axes 的属性。
在显示背景图像时,不能用imshow ()函数,如果用此函数,背景图像只能显示在封面的一部分,不能整屏覆盖;而要用imagesc ()函数。
MATLAB 提供修改文本中文字的字体,它支持华文中宋、华文彩云、华文仿宋、楷书、黑体等,功能较强。
在设计封面的过程中,最主要的技术是要解决时间的显示,如果只用MATLAB 中的clock ()函数,则显示的是静态时间,而不会显示和电脑同步的动态时间。
为了解决这个问题,需做一个循环判断语句:while find (get (0,'c hildren'))==h0。
成立的条件是,只要是当前窗口循环,否则停止,这样可以减少CPU 的负担。
同时,要实时提取clock ()函数,可用fix (clock),使提取的时间更美观。
#158#长江大学学报(自科版)2005年4月第2卷第4期/理工卷第2卷第2期Jour nal of Yangtze University (Nat Sci Edit)Apr 12005Vol12No 14/Sci &Eng V,Vol 12No 12图1 MATLAB 图像处理系统原理框图 图2 MATLAB 图像处理系统程序总流程图因为用clock ()函数提取的时间是用科学计数法表示的。
提取的时间还要转换成字符显示,用set 设置属性/string 0。
3 主界面设计当按/进入0键后,封面界面会消失,接着进入图像处理的主界面。
其功能包括图像增强和图像识别两大块。
这个主界面中还有一个打开文件的问题,因为以1DA T 为后缀的图像数据文件MAT LAB 不能直接处理,需要转换成其他格式,如1BMP,1GIF,1T IF,1JPG 和灰度图像。
按/打开0键后,出现标准Windows 的打开对话框,选中要处理的文件,需要对它进行处理。
如果是1DAT 文件,要先把它转换成灰度图像,可用mat2gr ay ()函数生成128@128的一维矩阵,也可用reshape (I ,128,128)把图像矩阵转换成128@128的矩阵;然后用imread ()函数直接读取图像文件显示。
转换后的图像数据用save 命令保存,然后在各个文件中传递。
按/返回0键,本界面消失返回到封面,按/清屏0键将清除显示在图像区的所有图像。
设置其他键,具有相应的子功能。
4 图像处理子功能界面设计以图像增强处理为例。
在设计图像增强功能时,主要设计了3个方面的处理:灰度变换、直方图处理、空域滤波。
在灰度变换中,可进行增强对比度、图像压缩、图像求反等处理。
在空域滤波中,可进行线性滤波、中值滤波、自适应滤波。
对被噪声污染的图像进行线性滤波,可以去除图像中某些类型的噪声。
如采用邻域平均法的均值滤波,就非常适应于去除通过扫描得到的测井图像中的颗粒噪声。
中值滤波是抑制噪声的非线性处理方法,它运算简单,易于实现,而且能较好地保护边界;但有时会失掉图像中的细线和小块的目标区域。
自适应滤波是用wiener2函数[5]根据图像的局部方差来调整滤波器的输出,当局部方差大时,滤波器的平滑效果较小;当局部方差小时,滤波器平滑效果强。
wiener2函数提供的自适应滤波通常比线性滤波的效果好,它比相应的线性滤波器具有更好的选择性,可以更好地保存图像的边缘和高频细节信息;另外,使用起来非常方便,wiener2函数同时计算出滤波器的参数,并对图像进行滤波计算,而且wiener 2函数并不比线性滤波器需要更多的计算时间。
#159#理工卷第2卷第2期黄书先:基于MATLAB 平台的可视化图像处理系统设计方法5 运动目标提取与图像边缘检测对图像序列中运动目标的提取最常用的方法是差值法。
这种方法的思想是:一个场景中两幅或多幅运动图像,其静止背景物的灰度值是恒定的,而运动目标的灰度值是变化的;将相邻的两幅图像相/减0或者/差分0,图像的静止背景物被消除了,而运动物体保留下一个粗略轮廓;在此基础上,采用模板对准相乘法,能准确地提取出图像序列中的运动目标。
图像边缘是所要提取目标和背景的分界线,提取出边缘才能将目标和背景区分开来。
用sobel 算子检测边缘效果很差,因为实际图像中,对应景物边缘的图像灰度变化有时并不十分陡峭;另外,图像中也存在噪声。
因此直接使用微分算子提取边界后,还需要做细化处理,才能形成一条有意义的边界。
改用canny 算子检测图像边缘效果较好。
[参考文献][1]张志涌1精通MAT LAB615版[M]1北京:北京航空航天大学出版社,20031[2]周明,李长虹1MAT LAB 图形技术[M]1西安:西北工业大学出版社,19991[3]何强,何英1MAT LAB 扩展编程[M]1北京:清华大学出版社,20021[4]沈精虎,于伟,郭照宇1微机图形用户界面设计方法与技巧[M]1北京:电子工业出版社19981[5]孙兆林1MAT LAB 61x 图像处理[M ]1北京:清华大学出版社,20021[编辑] 易国华(上接第145页)3 结 语笔者运用BCGS 2FFT 方法分析了两媒质半空间和三分层媒质中埋入目标的电磁脉冲散射问题,其快速算法的引用大大降低了计算时间和内存需求,使得普通的PC 机就可以用于解决平面分层背景下大尺度电磁散射问题。
仿真结果验证了模型和算法的可靠性,为埋地目标(石油、煤矿、管线等)的探测提供一定的理论指导。
[参考文献][1]M ichalsk i K A,zhen g D 1E lectrom agnetic scattering an radiation by s urface of arbitrary shape in layered media )))Part Ò:Implemen 2tation and resu lts for contiguous half 2space [J]1IEEE T rans 1Antennas Propagat,1990,38(3):344~3521[2]Michals ki K A,Mosig J R 1Multilayered media Green p s function in integral equation formulation s [J]1IE EE T rans 1An tennas Pr opa 2gat,1997,45(3):508~5191[3]H e J Q,Yu T J,Geng N,et al 1Method of moments analys is of electromagnetic s cattering from a gen etal th ree 2dim ens ion al dielectrictarget em bedded in a mu ltilayered m edium [J ]1Radio Sci 1,2000,35(2):305~3131[4]Su C C 1T he three dimen sional algorithm of solving th e electric field integal equation using face 2centered node points ,conjugate gradi 2ent m ethod ,and FFT 0IEEE T ran s [J]1Microwave Th eory Tech 1,1993,41(6):510~5151[5]Gan H,Chew W C 1A discr ete BCGS 2FFT algorithm for solving 3D in homogeneous scatter problem s [J]1J 1Electromagn 1Waves Ap 2pl,1995,9(10):1339~13571[6]Xu X M,Liu Q H 1T he BCGS 2FFT method for electromagnetic scattering form inh omogeneous objects in a planarly layer ed medium[J ]1IEEE Antennas and Wir eless Propagat 1Lett,2002,1(5):430~4391[7]Xu X M,Liu Q H,Zh ang Z Q 1T he stabilized biconjugate gradient fast Fourier transform method for electromagnetic s cattering [J]1J 1Appl 1Computat 1Electromagn 1Soc,2002,17(1):97~1031[8]M ichalsk i K A,Zheng D 1Electromagnetic s cattering an radiation by su rface of arbitrary shape in lay ered media )))Part Ñ:Th eory[J ]1IEEE T ran s 1Antennas Propagat,1990,38(3):335~3441[编辑] 易国华#160# 长江大学学报(自科版)2005年4月band can be improved by using albefaction filter,the interference resist ing property of the dual t ap filter is much better t han that of predict ed error filt er1Key words:spread spect rum communication;trace filter;white noise;albefaction treat2 ment155Method for Predicting the High Temperature High Pressure Performance of Oil2base Drill2 ing FluidZHA O H a i2ya n(Beijing Union University,Beij ing100009)ZHA N G Jin2bo(Great W ell Drilling Co1CN PC,100724)YA N Jie2nian(Univer s ity of Petroleum,Beijing102249)Abstract:Based on the new t heoret ical predicting model both at home and abroad,a software of prediction algorithm for high t emperat ure high pressure(HT HP)perform2 ance of oil base drilling fluid is implemented by using MATLAB1The software has the function of predict ing H TH P apparent viscosity and density of the oil base fluid,by which the apparent viscosit y and density changing with t he depth can be predict ed accu2 rately,instruct ion can be provided for it s field application1Key words:oil base drilling fluid;apparent viscosity;densit y;algorithm158Method for Designing A Visualizational Image Processing System Based on MATLABH UA N G Shu2xia n(Y ang tze University,J ing zhou434023)Abstract:By using MAT LAB as a working language and environment for development,a visualizational image processing system is developed under MATLAB environment,by which precision image process including well logging images is implemented and in2 terfaced with user2developed programs1A useful approach is made on the redevelop2 ment of MAT LAB and design of visualizational systems1Key words:MAT LAB;image process;visualized GUI interface161Design and Implementation of Point to Point Softphone Based on Session Initiation Proto2 cal(SIP)TU Ji2hui,ZHON G Ba o2rong(Yan gtze University,J ing zhou434023)Abstract:How to implement a software telephone model based on SIP is int roduced in det ail,and point to point communicat ion command in SIP is analyzed1T echniques are provided for network voice acquisition,transmit ting and playing1Key words:SIP;RTP;VoIP;multimedium164Analysis on Clearance Flow and Rigidity for A Balance Disc of Multistage Pump WAN G Jian2hua,L IU M ing2ya o,WA NG Ben2de(Yang tze Univer sity,J ing zh ou434023)Abstract:Based on the clearance flow of a balance disc,t he parameters of leakage rate, pressure,balance force and sensitivit y are derived,the rigidity concept of the balance disc is proposed,the relationships between the balance disc rigidity and it s sensitivity are approached,and the range of sensit ivity value is obtained,by which a theoretical basis and an effective method are provided for the design of balance disc1Key words:mult istage pump;balance disc;clearance;rigidity;analysis#×#。