【2017年整理】核辐射成像练习指导.docVIP

核辐射成像练习指导 整理调试 张天奎 审订指导 胡华四教授 2010年12月22日 图像处理与分析的关系框图 像素领域概念 四领域 4 4 P 4 4 八领域 8 8 8 8 P 8 8 8 8 Matlab使用技巧 help命令 图像处理工具箱（Image Processing） 图像的读取、写入和显示 例1 读取lena.jpg图像并显示，最后另存为lena.bmp （1）程序： % 读入图像 I=imread(lena.jpg); % 显示图像 figure,imshow(I); % 写出图像 imwrite(I,lena.bmp); （2）效果： 图像类型转换 例2 将2.jpg真彩图像转换成灰度图像 程序： % 读入图像 RGB=imread(2.jpg); % 真彩图像转化成灰度图像 GRAY=rgb2gray(RGB); % 显示图像 figure,subplot(221),imshow(RGB); subplot(222),imshow(GRAY); 效果 图像傅里叶变换 例3 将lena.jpg图像进行傅里叶变换，并显示其功率谱 程序： % 读入原始图像 I=imread(lena.jpg); imshow(I); % 求离散傅里叶频谱 J=fftshift(fft2(I));%对原始图像进行二维傅里叶变换，并将其中心移到零点 figure,imshow(log(abs(J)),[6 10]); 效果： Gamma校正 压缩低灰度，扩张高灰度（蓝色）； 压缩高灰度，扩张低灰度（红色） 例4 gamma校正实验 程序： I = imread(pout.tif);%pout.tif为matlab7自带图像 J = imadjust(I); % 显示图像 figure,subplot(121),imshow(I),title(I); subplot(122),imshow(J),title(J); % 灰度直方图 figure,subplot(121),imhist(I); subplot(122),imhist(J); 效果： 对应的灰度直方图 灰度直方图 灰度等级 具有灰度等级的像素数目 总的像素个数 灰度直方图均衡化：使得变换后的灰度直方图呈均匀分布。 例5 灰度直方图均衡化 程序： I = imread(tire.tif);%tire.tif为matlab7自带图像 J = histeq(I); % 显示图像 figure,subplot(121),imshow(I),title(I); subplot(122),imshow(J),title(J); % 灰度直方图 figure,subplot(121),imhist(I,64); axis square subplot(122),imhist(J,64); axis square 效果： 对应的灰度直方图 图像增强算法（空域） 滤波器 目的 负作用 高通滤波器 （锐化，Sharp） 增强边缘 放大噪声 低通滤波器 （平滑，Blur） 消除噪声 模糊边缘 均值平滑（线性平滑滤波器一种） 3*3模板 1/9 1/9 1/9 1/9 1/9 1/9 1/9 1/9 1/9 例6 不同尺度均值滤波去除椒盐噪声 程序： I=imread(eight.tif); J=imnoise(I,salt pepper,0.02);%添加椒盐噪声 figure, subplot(221),imshow(I) title(原图像); subplot(222),imshow(J) title(添加椒盐噪声图像); K1=filter2(fspecial(average,3),J);%应用3*3模板 K1= mat2gray(K1); subplot(223),imshow(K1); title(3*3窗的领域均值滤波); K2=filter2(fspecial(average,7),J);%应用7*7模板 K2= mat2gray(K2); subplot(224),imshow(K2); title(7*7窗的领域均值滤波); 效果： 中值滤波器（非线性平滑滤波器的一种） 先在领域里排序，然后挑中间的。可以有效的过滤椒盐噪声。 例7 中值滤波去除椒盐噪声 程序： I=imread(eigh