医学图像处理考试程序参考.doc

%直方图均衡化 clc,clear i=imread(cameraman.tif); %读入图像 a=histeq(i); %直方图均衡化 subplot(2,2,1); imshow(i); %显示图像 title(原图); subplot(2,2,2); imhist(i); %绘制直方图 title(原图的直方图); subplot(2,2,3); imshow(a); title(均衡化后的图像); subplot(2,2,4); imhist(a); title(直方图均衡化); %彩色图像转变成灰度图像 I=imread(lena.jpg); subplot(1,2,1); imshow(I); title(原彩色图) A=rgb2gray(I); %RGB-gray subplot(1,2,2); imshow(A); title(灰度图) subplot(1,2,1); %这里128为索引颜色种类，可以自己定义， imshow(I); 但最大为256 [X,map]=rgb2ind(I,128); %RGB图像转换索引图像，X存放灰度图像数 subplot(1,2,2); 据，map存放颜色索引表（具体可以看matlab imshow(X,map); 里面的workspace窗口的变量X和map） %imwrite(X,map,lenas,jpg); %图像旋转（绕图像中心旋转） I=imread(lena.jpg); subplot(221); imshow(I);title(原图); A=imrotate(I,30); %最邻近插值，逆时针旋转30度 subplot(222); imshow(A);title(逆时针旋转30度); B=imrotate(I,-30,bilinear,crop); %双线性插值， subplot(223); imshow(B);title(双线性插值顺时针30度,和原图像大小相等); C=imrotate(I,-30,bicubic); %双三次插值，顺时针 subplot(224); imshow(C);title(双三次插值顺时针30度); %图像的傅里叶变换 I1=imread(lena.bmp); %读入原图像文件 subplot(1,2,1); imshow(I1); %显示原图像 fftI1=fft2(I1); %二维离散傅立叶变换 sfftI1=fftshift(fftI1); %直流分量移到频谱中心 RR1=real(sfftI1); %取傅立叶变换的实部 II1=imag(sfftI1); %取傅立叶变换的虚部 A1=sqrt(RR1.^2+II1.^2); %计算频谱幅值 A1=(A1-min(min(A1)))/(max(max(A1))-min(min(A1)))*225; %归一化 subplot(1,2,2); imshow(A1); %显示原图像的频谱 %图像锐化 %Sobel算子 h1=[-1 0 1;-2 0 2;-1 0 1]; % h2=[1 2 1;0 0 0;-1 -2 -1]; M=imread(meet.jpg); N=rgb2gray(M); I=double(N); %将数据转换为双精度格式，便于后面计算 figure;imshow(N);title(原图); J1=conv2(I,h1,same); %进行卷积运算，从结果中取出与I大小相同的一部分 J2=conv2(I,h2,same); J=abs(J1)+abs(J2); %绝对值相加 T=uint8(J); %将数据格式再转变回图像格式 figure;imshow(T);title(sobel锐化); figure; imshow(255-T);title(sobel反片); % Roberts for j=2:581 for i=2:419 g(i,j)=abs(I(i+1,j+1)-I(i,j))+abs(I(i+1,j)-I(i,j+1)); % r