《最新pet图像肿瘤的检测-结题报告剖析.doc》由会员分享，可在线阅读，更多相关《最新pet图像肿瘤的检测-结题报告剖析.doc（65页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、Four short words sum up what has lifted most successful individuals above the crowd: a little bit more.-author-datepet图像肿瘤的检测-结题报告剖析pet图像肿瘤的检测-结题报告剖析南京理工大学本科科研训练报告作者：XXXXXXXXXXXXXXXXXX学院：计算机科学与工程学院题目：PET图象肿瘤的检测指导老师：XX XXX年 月 日第一章 绪论1.1研究背景肿瘤的早期预防、早期发现和治疗对于人的生命是至关重要的，PET检测技术的使用为肿瘤患者诊断方案的确定提供了便利条件，它独特
2、的成像技术和PET图像的分析与研究，也越来越受到学者和研究人员的广泛关注。准确的将PET图像中的肿瘤区域分割出来具有要的临床意义。PET检测技术的引进，给癌症患者带来了一种准确的诊断方案。它能收到人们如此热切的关注，因为PETCT相比于传统的诊断手段有了许多明显的优势。在PET临床应用中，主要用到数字图像技术中的图像分割技术。运用图象处理中的分割方法能够提取影像数据中特殊组织，也是可视化实现的前提，分割后的区域能够进行组织容积的定量分析，诊断，治疗规划等，从而实现计算机辅助诊断分析。这就要求参加此项目的同学对数字图像学有一定的了解，尤其是要掌握其中的图像分割的多种方法，例如边缘检测阈值分割区域
3、分裂与合并等等，当然这离不开计算机技术的配合，该项目要求同学有较强的编程能力，最好能将图像的检测与分割及计算等用代码实现，利用计算机最终实现肿瘤的较为精确的检测和定量分析（主要是算其体积）。1.2研究目的和内容1.2.1 研究目的本项目的主要实现了利用分割算法来实现肿瘤区域的检测分割。首先学习图像处理软件库Matlab，熟悉其各项功能以及相关函数，然后针对PET序列影像数据，研究肿瘤的属性，运用多种方法对肿瘤进行检测并分割，最后实现分割肿瘤的各项特征计算，与医学专家分割的结果进行对比分析，给出定量的分析评价。通过该项目的研究，项目参与人员能够锻炼实际的编程能力，熟悉医学影像分析的流程，掌握图象
4、处理的相关知识，提高解决实际问题的能力。1.2.2 研究内容肿瘤的早期预防、早期发现和治疗对于人的生命是至关重要的。近年来，PET在诊断和指导治疗肿瘤、冠心病和脑部疾病等方面均已显示出独特的优越性，越来越受到学者和研究人员的广泛关注，准确的将PET图像中的肿瘤区域检测分割出来具有重要的临床意义。运用图象处理中的分割方法能够提取影像数据中特殊组织，也是可视化实现的前提，分割后的区域能够进行组织容积的定量分析，诊断，治疗规划等，从而实现计算机辅助诊断分析。第二章 MATLAB简单介绍2.1简介MATLAB是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算
5、的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。2.2优势特点1) 高效的数值计算及符号计算功能，能使用户从繁杂的数学运算分析中解脱出来；2) 具有完备的图形处理功能，实现计算结果和编程的可视化；3) 友好的用户界面及接近数学表达式的自然化语言，使学者易于学习和掌握；4) 功能丰富的应用工具箱(如信号处理工具箱、通信工具箱等) ，为用户提供了大量方便实用的处理工具。2.3编程环境MATLAB由一系列工具组成。这些工具方便用户使用MATLAB的函数和文件，其中许多工具采用的是图形用户界面。包括MATLAB桌面和命令窗口、历史命令窗口、编辑器和调试器、路径搜索和用于
6、用户浏览帮助、工作空间、文件的浏览器。2.4图形处理MATLAB自产生之日起就具有方便的数据可视化功能，以将向量和矩阵用图形表现出来，并且可以对图形进行标注和打印。高层次的作图包括二维和三维的可视化、图象处理、动画和表达式作图。可用于科学计算和工程绘图。新版本的MATLAB对整个图形处理功能作了很大的改进和完善，使它不仅在一般数据可视化软件都具有的功能（例如二维曲线和三维曲面的绘制和处理等）方面更加完善，而且对于一些其他软件所没有的功能（例如图形的光照处理、色度处理以及四维数据的表现等），MATLAB同样表现了出色的处理能力。同时对一些特殊的可视化要求，例如图形对话等，MATLAB也有相应的功
7、能函数，保证了用户不同层次的要求。另外新版本的MATLAB还着重在图形用户界面（GUI）的制作上作了很大的改善，对这方面有特殊要求的用户也可以得到满足。2.5 MATLAB在图像处理中的应用图像处理工具包是由一系列支持图像处理操作的函数组成的。所支持的图像处理操作有:图像的几何操作、邻域和区域操作、图像变换、图像恢复与增强、线性滤波和滤波器设计、变换(DCT变换等) 、图像分析和统计、二值图像操作等。下面就MATLAB 在图像处理中各方面的应用分别进行介绍。 （1） 图像文件格式的读写和显示。MATLAB 提供了图像文件读入函数 imread()，用来读取如:bmp、tif、tiffpcx 、
8、jpg 、gpeg 、hdf、xwd等格式图像文；图像写出函数 imwrite() ，还有图像显示函数 image()、imshow()等等。 （2） 图像处理的基本运算。MATLAB 提供了图像的和、差等线性运算 ，以及卷积、相关、滤波等非线性算。例如，conv2(I，J)实现了I，J两幅图像的卷积。（3） 图像变换。MATLAB 提供了一维和二维离散傅立叶变换(DFT) 、快速傅立叶变换(FFT) 、离散余弦变换 (DCT) 及其反变换函数，以及连续小波变换(CWT)、离散小波变换(DWT)及其反变换。（4） 图像的分析和增强。针对图像的统计计算MATLAB 提供了校正、直方图均衡、中值滤
9、波、对比度调整、自适应滤波等对图像进行的处理。（5） 图像的数学形态学处理。针对二值图像，MATLAB 提供了数学形态学运算函数；腐蚀(Erode)、膨胀(Dilate)算子，以及在此基础上的开 (Open)、闭(Close)算子、厚化 (Thicken) 、薄化 (Thin) 算子等丰富的数学形态学运算。以上所提到的 MATLAB在图像中的应用都是由相应的MATLAB函数来实现的，使用时，只需按照函数的调用语法正确输入参数即可。具体的用法可参考MATLAB丰富的帮助文档。图像边缘对图像识别和计算机分析十分有用，在MATLAB中，函数 edge()用于灰度图像边缘的提取，它支持六种不同的边缘提
10、取方法，即Sobel方法、Prewitt 方法、Robert 方法，Laplacian2Gaussian方法、过零点方法和Canny方法。第三章 数字图像处理3.1概念数字图像处理（Digital Image Processing）又称为计算机图像处理，它是指将图像信号转换成数字信号并利用。数字图像处理(Digital Image Processing)是通过计算机对图像进行去除噪声、增强、复原、分割、提取特征等处理的方法和技术。3.2主要目的（1）提高图像的视感质量，对图像进行几何变换等，以改善图像的质量。（2）提取图像中所包含的某些特征或特殊信息，这些被提取的特征或信息往往为计算机分析图像
11、提供便利。（3）图像数据的变换、编码和压缩，以便于图像的存储和传输。3.3常用方法1）图像变换：由于图像阵列很大，直接在空间域中进行处理，涉及计算量很大。因此，往往采用各种图像变换的方法，如傅立叶变换、沃尔什变换、离散余弦变换等间接处理技术，将空间域的处理转换为变换域处理，不仅可减少计算量，而且可获得更有效的处理（如傅立叶变换可在频域中进行数字滤波处理）。2 ）图像编码压缩：图像编码压缩技术可减少描述图像的数据量（即比特数），以便节省图像传输、处理时间和减少所占用的存储器容量。编码是压缩技术中最重要的方法，它在图像处理技术中是发展最早且比较成熟的技术。3 ）图像增强和复原：图像增强和复原的目的
12、是为了提高图像的质量，如去除噪声，提高图像的清晰度等。图像复原要求对图像降质的原因有一定的了解，一般讲应根据降质过程建立“降质模型”，再采用某种滤波方法，恢复或重建原来的图像。4 ）图像分割：图像分割是数字图像处理中的关键技术之一。图像分割是将图像中有意义的特征部分提取出来，其有意义的特征有图像中的边缘、区域等，这是进一步进行图像识别、分析和理解的基础。虽然已研究出不少边缘提取、区域分割的方法，但还没有一种普遍适用于各种图像的有效方法。因此，对图像分割的研究还在不断深入之中，是图像处理中研究的热点之一。5 ）图像描述：图像描述是图像识别和理解的必要前提。作为最简单的二值图像可采用其几何特性描述
13、物体的特性，一般图像的描述方法采用二维形状描述，它有边界描述和区域描述两类方法。6 ）图像分类（识别）：图像分类（识别）属于模式识别的范畴，其主要内容是图像经过某些预处理（增强、复原、压缩）后，进行图像分割和特征提取，从而进行判决分类。3.4主要优点1. 再现性好。数字图像处理与模拟图像处理的根本不同在于：它不会因图像的存储、传输或复制等一系列变换操作而导致图像质量的退化。只要图像在数字化时准确地表现了原稿，则数字图像处理过程始终能保持图像的再现。2. 数组，这主要取决于图像数字化设备的能力。现代扫描仪可以把每个像素的灰度等级量化为16位甚至更高，这意味着图像的数字化精度可以达到满足任一应用需
14、求。3适用面宽。图像可以来自多种信息源。这些来自不同信息源的图像只要被变换为数字编码形式后，均是用二维数组表示的灰度图像（彩色图像也是由灰度图像组合成的，因而均可用计算机来处理。即只要针对不同的图像信息源，采取相应的图像信息采集措施，图像的数字处理方法适用于任何一种图像。4灵活性高图像处理大体上可分为图像的像质改善、图像分析和图像重建三大部分，每一部分均包含丰富的内容。由于图像的光学处理从原理上讲只能进行线性运算，这极大地限制了光学图像处理能实现的目标。而数字图像处理不仅能完成线性运算，而且能实现非线性处理，即凡是可以用数学公式或逻辑关系来表达的一切运算均可用数字图像处理实现。第四章 用MAT
15、LAB进行PET图象肿瘤的检测4.1图像的读取与显示1、读取图像 函数imread可以从任何Matlab支持的图像文件格式中，以任意位深度读取一幅图像。格式为：MAP=imread(FILENAME,FMT)，其中：FILENAME为需要读入的图像文件名称，FMT为图像格式。2、读取图像信息 可以通过调用imfinfo函数获得与图像文件有关的信息，格式如下：INFO=imfinfo(FILENAME,FMT)其中：返回的INFO是Matlab的一个结构体。3、显示图像Matlab的图像显示函数主要使用imshow。 imshow(I,n) : 显示灰度级为n的图像，n缺省为256。imshow
16、(I,low,high) :以灰度范围low,high显示图像，如果不知道灰度范围，可以用imshow(I,)显示。imshow(BW) :显示二值图像。4.2阈值分割(整个科研实验的重点)（1）调入并显示图像rice.png。（2）用迭代式阈值选择方法分割图像：（3）用Otsu法阈值选择法分割图像；（4）用分水岭算法分割图像。（5）将上述几种处理方法的结果作比较。在一幅图像中，人们常常只对其中的部分目标感兴趣，这些目标通常占据一定的区域，并且在某些特性（如灰度、轮廓、颜色和纹理等）上和临近的图像有差别。这些特性差别可能非常明显，也可能很细微，以至肉眼察觉不出来。随着计算机图像处理技术的发展，
17、使得人们可以通过计算机来获取和处理图像信息。图像识别的基础是图像分割，其作用是把反映物体真实情况的、占据不同区域的、具有不同特性的目标区分开来，并形成数字特征。图像分割是图像理解的基础，而在理论上图像分割又依赖图像理解，它们彼此是紧密关联的。图像分割在一般意义下是十分困难的问题，目前的图像分割一般作为图像的前期处理阶段，是针对分割对象的技术，最常用到的是利用阈值化处理进行的图像分割。所谓阈值的方法实质是利用图像的灰度直方图信息得到用于分割的阈值。基于阈值的分割方法可以分为全局阈值的方法好局部阈值的方法，全局阈值的方法是利用整幅图像的灰度信息，从整个图像中得到用于分割的阈值，并且根据该阈值对图像
18、进行分割；而局部阈值的方法是根据图像中不同区域获得对应不同区域的几个阈值，利用这些得到的阈值对图像进行分割，也就是一个阈值对应图像中的一个子区域。该方法的关键是如何合理地选择阈值，人工选择法是通过人眼的观察，应用人们对图像的知识，在分析图像直方图的基础上人工选出合适的阈值，也可以在人工选出阈值后，根据分割效果不断地交换操作，从而选择出较佳的阈值；自动阈值法通常使用灰度直方图来分析图像中灰度值的分布，结合特定的应用区域知识来选取合适的阈值。（1）迭代式阈值选择迭代式阈值选择方法的基本思想是：开始时候选择一个阈值作为初始计值，然后按某种策略不断地改进这一估计值，直到满意给定的准则为止。在迭代过程中
19、，关键之处在于选择什么样的阈值改进策略。好的与之改进侧率应该具备两个特征：一是能够快速收敛；二是在每一个迭代过程中，新产生阈值优于上一次的阈值。如下是具体操作步骤：a 选择图像灰度的中值作为初始阈值；b 利用此阈值把图像分割成两个区域和，用下式计算区域和的灰度均值和 c 计算出和后，计算新的阈值d 重复步骤2和3，直到和的差小于某个给定值。下面为实现迭代式阈值选择的MATLAB代码：f=imread(rice.png);subplot(121);imshow(f);title(原始图像)f=double(f);T=(min(f(:)+max(f(:)/2;done=false;i=0;whil
20、e done r1=find(fT); Tnew=(mean(f(r1)+mean(f(r2)/2; done=abs(Tnew-T)1; T=Tnew; i=i+1;endf(r1)=0;f(r2)=1;subplot(122);imshow(f);title(迭代阈值二值化图像) 用迭代式阈值选择法二值化图像（2）Otsu法阈值选择Otsu法是一种使类间方差最大的自动确定阈值的方法，该方法具有简单、处理速度快的特点，是一种常用的阈值选取方法。MATLAB工具箱提供的graythresh函数求取阈值采用的就是这种方法。其基本思想是：设图像像素数为N，灰度范围为0，L-1，对应灰度级i的像素数
21、为，概率为 i=0,1,2,L-1 把图像中的像素按灰度值用阈值T分成两类和，由灰度值在0,T之间的像素组成，由灰度值在T+1,L-1之间的像素组成，对于灰度分布概率，整幅图像的均值为则和的均值为，其中，由上面三式可得，类间方差定义为让T在0,L-1范围依次取值，使得方差最大的T值即为Otsu法的最佳阈值。下面为用Otsu方法进行阈值选择的MATLAB代码：f=imread(rice.png);subplot(121);imshow(f);title(原始图像)T=graythresh(f);g=im2bw(f,T);subplot(122);imshow(g);title(Otsu方法二值化
22、图像)用Otsu方法二值化图像（3）分水岭算法分水岭算法（watershed）是一种借鉴了形态学理论的分割方法，在该方法中，将一幅图像看成一个拓扑地形图，其中灰度值对应地形高度值，高灰度值对应着山峰，低灰度值对应着山谷。水总是朝着地势低的地方流动，直到某一局部低洼处才停下来，这个低洼处被称为吸水盆地。最终所有的水会分聚在不同的吸水盆地，吸水盆地之间的山脊被称为分水岭。水从分水岭流下时，它朝着不同的吸水盆地流去的可能性是相等的。将这种想法应用于图像分割，就是要在灰度图像中找到不同的吸水盆地和分水岭，由这些不同的吸水盆地和分水岭组成的区域即为我们要分割的目标。 阈值T0 阈值T1 阈值T2 分水岭
23、形成示意图I=imread(eight.tif);subplot(221);imshow(I);title(原始图像)subplot(222);I=double(I);hv=fspecial(prewitt);hh=hv.;gv=abs(imfilter(I,hv,replicate);gh=abs(imfilter(I,hh,replicate);g=sqrt(gv.2+gh.2);subplot(222);L=watershed(g);wr=L=0;imshow(wr);title(分水岭)I(wr)=255;subplot(223);imshow(uint8(I);title(分割结果)
24、rm=imregionalmin(g);subplot(224);imshow(rm);title(局部极小值)不准确标记分水岭算法导致的过分割从上图可以看到对应于目标的边缘，出现了比较严重的过分割现象，原因在于分水岭算法是以梯度图的局部极小点作为吸水盆地的标记点。从第四幅图中可以看到有过多的局部极小点。下面的改进程序克服了这个缺点。I=imread(eight.tif);subplot(231);imshow(I);I=double(I);hv=fspecial(prewitt);hh=hv.;gv=abs(imfilter(I,hv,replicate);gh=abs(imfilter(I
25、,hh,replicate);g=sqrt(gv.2+gh.2);subplot(232);df=bwdist(I);imshow(uint8(df*8);L=watershed(df);em=L=0;subplot(233);imshow(em);im=imextendedmax(I,20);subplot(234);imshow(im);g2=imimposemin(g,im|em);subplot(235);imshow(g2);L2=watershed(g2);wr2=L2=0;subplot(236);I(wr2)=255;imshow(uint8(I);准确标记的分水岭算法分割过程
26、4.3图像分割1、图像分割概述 图像分割一般采用的方法有边缘检测（edge detection）、边界跟踪（edge tracing）、区域生长（region growing）、区域分离和聚合等。 图像分割算法一般基于图像灰度值的不连续性或其相似性。不连续性是基于图像灰度的不连续变化分割图像，如针对图像的边缘有边缘检测、边界跟踪等算法。相似性是依据事先制定的准则将图像分割为相似的区域，如阈值分割、区域生长等。2、边缘检测图像的边缘点是指图像中周围像素灰度有阶跃变化或屋顶变化的那些像素点，即灰度值导数较大或极大的地方。边缘检测可以大幅度的减少数据量，并且剔除不相关信息，保留图像重要的结构属性。边
27、缘检测基本步骤：平滑滤波、锐化滤波、边缘判定、边缘连接。边缘检测算法： 基于一阶导数：Roberts算子、Sobel算子、Prewitt算子 基于二阶导数：高斯-拉普拉斯边缘检测算子 Canny边缘检测算法Matlab实现1）、基于梯度算子的边缘检测BW=edge(I,type,thresh,direction,nothinning)type合法取值梯度算子sobelsobel算子prewittprewitt算子rebertsrobert算子thresh是敏感度阈值参数，任何灰度值低于此阈值的边缘将不会被检测到。默认值为空矩阵，此时算法自动计算阈值。direction指定了我们感兴趣的边缘方向
28、，edge函数将只检测direction中指定方向的边缘，其合法值如下：direction合法值边缘方向horizontal水平方向vertical竖直方向both所有方向可选参数nothinning，指定时可以通过跳过边缘细化算法来加快算法运行的速度。默认是thinning，即进行边缘细化。2）、基于高斯-拉普拉斯算子的边缘检测BW=edge(I,log,thresh,sigma)sigma指定生成高斯滤波器所使用的标准差。默认时，标准差为2。滤镜大小n*n，n的计算方法为：n=ceil(sigma*3)*2+1。3）、基于Canny算子的边缘检测BW=edge(I,canny,thresh
29、,sigma)thresh是敏感度阈值参数，默认值为空矩阵。此处为一列向量，为算法指定阈值的上下限。第一个元素为阈值下限，第二个元素为阈值上限。如果只指定一个阈值元素，则默认此元素为阈值上限，其0.4倍的值作为阈值下限。如阈值参数没有指定，则算法自行确定敏感度阈值上下限。b1=imread(nir.bmp);h58=fspecial(gaussian,5,0.8);b=imfilter(b1,h58);bw1=edge(b,sobel);%sobel算子bw2=edge(b,prewitt);%prewitt算子bw3=edge(b,roberts);%roberts算子bw4=edge(b,
30、log); %log算子bw5=edge(b,canny);%canny算子figure;imshow(bw1);imwrite(bw1,nirbwsobel.bmp);figure;imshow(bw2);imwrite(bw2,nirbwprewitt.bmp);figure;imshow(bw3);imwrite(bw3,nirbwroberts.bmp);figure;imshow(bw4);imwrite(bw4,nirbwlog.bmp);figure;imshow(bw5);imwrite(bw5,nirbwcanny.bmp);roberts算子log算子分析：1、边缘定位精度
31、方面： Roberts算子和Log算子定位精度较高。Roberts算子简单直观，Log算子利用二阶导数零交叉特性检测边缘。但Log算子只能获得边缘位置信息，不能得到边缘方向信息。2、边缘方向的敏感性： Sobel算子、Prewitt算子检测斜向阶跃边缘效果较好，Roberts算子检测水平和垂直边缘效果较好。Log算子不具有边缘方向检测功能。Sobel算子能提供最精确的边缘方向估计。3、去噪能力： Roberts算子和Log算子虽然定位精度高，但受噪声影响大。 Sobel算子和Prewitt算子模板相对较大因而去噪能力较强，具有平滑作用，能滤除一些噪声，去掉一部分伪边缘，但同时也平滑了真正的边缘
32、，降低了其边缘定位精度。总体来讲，Canny算子边缘定位精确性和抗噪声能力效果较好，是一个折中方案。4.4图像特征提取1.特征提取的一般原则：选择在同类图像之间差异较小（较小的类内距），在不同类别的图像之间差异较大（较大的类间距）的图像特征。2简单区域描绘子 周长：区域边界上的像素数目。 面积：区域中像素数目。 致密性：（周长）2/面积。 区域的质心。 灰度均值：区域中所有像素的平均值。 灰度中值：区域中所有像素的排序中值。 包含区域的最小矩形。 最小或最大灰度级。 大于或小于均值的像素数。 欧拉数：区域中的对象数减去这些对象的孔洞数。D=regionprops(L,properties);L
33、为一个标记矩阵，通过连通1区域标注函数bwlabel得到。Properties合法值含义Area区域内像素总数BoundingBox包含区域的最小矩形Centroid区域的质心ConvexHull包含区域的最小凸多边形EquivDiameter和区域有着相同面积的圆的直径EulerNumber区域中的对象数减去这些对象的孔洞数3、主成分分析（Principal Component Analysis，PCA）COEFF,SCORE,latent=pincomp(X); X为原始样本组成n*d的矩阵，其每一行是一个样本特征向量，每一列表示样本特征向量的一维。 COEFF：主成分分量，也是样本协方差
34、矩阵的本征向量。 SCORE：主成分，X的低维表示。 latent：一个包含着样本协方差矩阵本征值的向量。第五章 科研项目总结和展望5.1总结与心得通过这次科研项目，我们团队收获了很多。首先，我们掌握了课本上没有的知识，了解了有关于数字图像方面的知识，在指导老师曹国的带领和帮助下，我们利用课余时间上网查资料，用工具，大大提高了我们的动手能力。其次，我们学习到了一些有关报告与论文的书写格式与要求，学到了很多，为我们以后自己完成毕业设计的报告打下了基础，最后，关于我们进行的基于Matlab的PET图象肿瘤的检测研究，是我们在进入大学之后，研究时间最长，投入时间最多的一个项目，也深深地体会到研究的不
35、易，由于自身能力的原因，我们进行的研究层次有限。在看起来比较简单程序上我们也会耗费比较多的时间，此次科研训练对我们最重要的影响就是提高了我们的学习能力，刚开始在遇到一个完全不明白的工具的时候，我们的反应非常一致就是懵，在曹老师的指导以及督促下，我们从对这个项目一窍不通到可以用Matlab工具进行肿瘤图像处理，达到了科研训练对于我们的锻炼，也提升了我们团队每个人的个人能力，这次科研项目对我们可以说是意义重大，影响深远。在论文完成之际，我心情无法平静。大学期间，尤其是毕业论文完成期间，有可敬的老师和同学给予了我无私帮助。对此我深怀感激。首先感谢导师曹老师从始至终的关心、指导和教诲。曹老师始终如一关
36、心科研项目的进展，对我们小组的不足给与了很好的提示与指导，其兢兢业业的敬业精神一直鞭策我们前进，。其次在设计的过程中，感谢给予我帮助的小组成员们，在此对其表示感谢。 最后感谢学校无私的给我们提供了这样一个锻炼机会。5.2今后的工作这次科研训练历时一年，但由于我们利用的是课余时间而且是本科生，所以项目的成果还有很大的提升空间，如果能力可行，我们希望可以将加深对数字图像处理基础知识的了解，加强对Matlab工具的使用，力求将项目的要求完成好。今后的话，我们希望工作可以着力于完成对于图像阈值分割阈值这一方面的研究，并试着把研究方向从PET图像转向于多种来源的图像。在未来的学术生涯中有更多有关图像处理
37、方面的学习机会和收获。5.3展望随着计算机处理能力的不断增强，数字图像处理学科在飞速发展的同时，也越来越广泛地向许多其他学科快速交叉渗透，使得图像作为信息获取以及信息的利用等方面也变得越来越重要。目前数字图像处理的应用越来越广泛，已经渗透到工业、医疗保健、航空航天、军事等各个领域，在国民经济中发挥越来越大的作用。MathWorks 公司开发的MATLAB最突出的特点就是简洁。MATLAB用更直观的，符合人们思维习惯的代码，代替了C和FORTRAN语言的冗长代码。MATLAB给用户带来的是最直观，最简洁的程序开发环境。应用MATLAB友好的界面和丰富、实用、高效的指令及模块，可以使人较快地认识、
38、理解图像处理的相关概念，逐步掌握图像信号处理的基本方法，进而能够解决相关的工程和科研中的问题。利用MATLAB来实现数字图像处理使数字图像处理更加方便简单。图像处理技术在许多应用领域受到广泛重视并取得了重大的开拓性成就，属于这些领域的有航空航天、生物医学工程、工业检测、机器人视觉、公安司法、军事制导、文化艺术等，使数字图像处理成为一门引人注目、前景远大的新型学科。参考文献于万波，基于MATLAB的图像处理，清华大学出版社，2008.32飞思科技产品研发中心，MATLAB6.5辅助图像处理，电子工业出版社，2003.13章毓晋，图像处理与分析，北京清华大学出版社 2004.74张志涌，精通MATLAB6.5， 北京北航电子版，2002.125闫敬文，数字图像处理MATLAB版，国防工业出版社，2007.26周新伦，柳建，数字图像处理，国防工业出版社，1986.3-