MatLab图像处理基本技术初探.pdfVIP

较。 2.脚本程序处理图像：(1)在MatLab 下，使用相关函数来读取Dicom 格式文件，核心函数有 Dicomread()、imshow()，由于Dicom 图像是特殊数据类型，在进行处理之前还需要转换图像数据类型， 采用mat2gray()转化成灰度图，以备后续程序处理。(2)1.直方图分析：直方图显示Dicomread()函数用于 读入各种图像文件， imshow()函数用于显示图像，imhist()函数用于显示直方图；hgram()函数是规定化 直方图的范围，可调整；histeq()函数实现对图像的增强。(3)空域波滤处理：medfilt2(gray_data，[5 5]) 函数进行中值滤波处理，[5 5]为滤波模板，可以通过改变范围，得到满足要求的滤波效果。应用Laplacian 算子和Prewitt 模板进行锐化滤波处理，观察滤波效果。(4)边缘检测处理：应用微分算子（Robert，Sobel， Prewitt 等）、拉普拉斯算子（Log）、Canny 算子对图像进行边缘检测处理。对以上算子分别以不同阈 值输入进行对比，得出有价值的相关边缘检测数据。 结 果 一、图像处理工具箱处理 1.图像处理工具箱是可视化处理界面，直接利用菜单及工具就可操作对Dicom格式文件的读取、处 理等。并可获取获取准确、全面的图片信息，这是后期医学图像处理的基础。 2.Dicom图像的信息概要主要包括图像细（image details）和元数据（metadate）两部分。图像细节 主要包括宽度（列）、高度（行）、等级、图像类型、最小强度和最大强度。元数据主要包括文件名、格 式、序列号、日期、患者信息、医院、机器类型等等信息。 3.通过改变窗口技术，针对观察不同的组织，选择合适的窗位、窗宽对观察部位进行灰度放大，获 取合适的窗口图像，才能清晰准确地观察到不同部位的组织机构。图像的密度越高，灰度值越高，CT 值也越高，灰度值与CT值之间具有对应关系。 二、DICOM文件读取、转换与显示 1.Dicom文件的读取及显示:主要函数有：Dicomread(‘filename’)，figure，imshow(A，[]) 2.转换为灰度图像并显示以备后继处理:主要函数有：J=mat2gray(A)，figure，imshow(J) 3.核心脚本处理程序文件（read_dcm.m）：function read_dcm() READ_DCM 功能读取DICOM文件.转换并显示read_dcm() global gray_data %声明全局变量clc; filename=input(Dicom filename:，s);init_dcm=Dicomread(filename); subplot(121)，imshow(init_dcm，[])，title(原始Dicom图像) gray_data=mat2gray(init_dcm); %转换为灰度图像并显示以备后继处理 subplot(122)，imshow(gray_data)，title(转换化后的灰度图像) 三、直方图显示与规范化核心脚本文件 clf;read_dcm()；global gray_data； subplot(221)，imshow(gray_data)，title(原始图像) subplot(222)，imhist(gray_data)，title(原始图像直方图) subplot(122)，imshow(gray_data，[0.65，0.8])，title(调窗后图像) 150 % hgram是规定化直方图的范围，可调整 hgram=0.5:0.01:1;JH=histeq(gray_data，hgram);% figure; subplot(223)，imshow(JH，[0.5，1.0])，title(规范化后图像) subplot(224)，imhist(JH)，title(规范化后直方图) 四、滤波处理 1.中值滤波处理核心脚本文件：clf;read_dcm()；global gray_data med_data=medfilt2(gray_data，[5 5]);% [5 5]为滤波模板（可调整），可以通过改变范围，得到满足要 求的滤波效果。