数字图像处理MATLAB作业(3)

3.1I=imread('bmp.bmp'); %把一个BMP灰度图像读入到变量I中I1=im2double(I); %图像灰度值的双精度化figure,imhist(I1); %打开一个新图形窗口来显示直方图I2=imadjust(I1,[0 0.3],[0 1]); %把图像的灰度范围从[0,0.3]拉伸到[0,1] figure,imshow(I2); %打开一个新图形窗口显示增强后的图像figure,imhist(I2); %打开一个新图形窗口显示增强后的直方图3.2I=imread('lf.bmp'); imshow(I);%读入并显示一幅BMP频谱图像I1=double(I);%双精度化I2=log(I1+1);%对数变换I3=mat2gray(I2);%把图像的灰度范围变换为[0 1]范围figure,imshow(I3); %打开一个新的图形窗口显示减少动态范围后的频谱图像3.4I=imread('fj.bmp'); I=rgb2gray(I);imshow(I);title('原始图像');%读入并显示一幅灰度图像figure,imhist(I);title('原始直方图');%打开一个新的图形窗口显示原始图像的直方图[I1,T]=histeq(I);%直方图均衡化figure,imshow(I1);title('均衡化后'); %打开一个新的图形窗口显示均衡化后的图像figure,imhist(I1);title('均衡化后直方图');%打开一个新的图形窗口显示均衡化后的直方图figure,plot((0:255)/255,T);title('变换曲线');%打开一个新的图形窗口画出变换曲线3.6I=imread('fj.bmp'); I=rgb2gray(I);imshow(I);title('原始图像');%读入并显示一幅灰度图像figure,imhist(I);title('原始直方图');%打开一个新的图形窗口显示原始图像的直方图[I1,T]=histeq(I);%直方图均衡化figure,imshow(I1);title('均衡化后'); %打开一个新的图形窗口显示均衡化后的图像figure,imhist(I1);title('均衡化后直方图');%打开一个新的图形窗口显示均衡化后的直方图figure,plot((0:255)/255,T);title('变换曲线');%打开一个新的图形窗口画出变换曲线3.7clc;clear all;I=imread('hua.bmp');imshow(I),title('hua.bmp'); %把一个图形窗口划分为1×2矩形显示区域，在图形窗口的左侧区域显示图像II1=double(I);%数据类型转换，不支持无符号整型计算[M,N]=size(I1);%计算图像的高和宽%产生高斯型高通滤波器m=round(M/2);n=round(N/2);%四舍五入取整hh=1.035;hl=0.7;c=3;for i=1:Mfor j=1:Nd(i,j)=sqrt((i-m)^2+(j-n)^2);endendd0=median(median(d));for i=1:Mfor j=1:Nh(i,j)=(hh-hl)*(1-exp(-c*(d(i,j)/d0)))+hl;%高通滤波器endend[i,j]=meshgrid(1:M,1:N);figure,mesh(h); title('highpass filter');%同态滤波h=h';I2=log(I1+1);%对数变换I3=fft2(I2);%傅立叶变换I3=fftshift(I3);%频谱中心化I3=reshape(I3,384,1152);I4=I3.*(h);%滤波I4=reshape(I4,384,384,3);I4=ifftshift(I4); %频谱反中心化I4=ifft2(I4);%傅立叶反变换I5=exp(I4-1);%指数变换I6=real(I5);%取幅值I7=mat2gray(I6);%恢复变换后的数值范围figure,imshow(I7);title('adjusted image intensity value');%在图形窗口的右侧区域显示结果图像3.8I=imread('lf.bmp');I=rgb2gray(I);imshow(I);%读入和显示一幅原始图像title('original image') %给图像加标题X=grayslice(I,16);%将灰度范围分成16层% X=grayslice (I,n)用阈值1/n, 2/n, ..., (n-1)/n 阈值化灰度图像I ，并返回索引图像X。

figure,imshow(X,hot(16));%黑红黄白浓淡色的伪彩色title('gray slice image'); %给图像加标题3.9I = imread('lf.bmp');I=rgb2gray(I);figure;subplot(1,2,1);imshow(I);title('原始图像'); [M,N] = size(I);J = double(I);L = 256;for i = 1:Mfor j = 1:Nif J(i,j) < L/4R(i,j) = 0;G(i,j) = 4*J(i,j);B(i,j) = L;elseif J(i,j) <=L/2R(i,j) = 0;G(i,j) = L;B(i,j) = -4*(J(i,j)-L/2);elseif J(i,j) <=3*L/4R(i,j) = 4*(J(i,j)-L/2);G(i,j) = L;B(i,j) = 0;elseR(i,j) = L;G(i,j) = -4*(J(i,j)-L);B(i,j) = 0;endendendOUT = repmat(I,[1 1 3]);for i = 1:Mfor j = 1:NOUT(i,j,1) = R(i,j);OUT(i,j,2) = G(i,j);OUT(i,j,3) = B(i,j);endendsubplot(1,2,2);imshow(OUT);title('灰度变换法伪彩色增强后图像')3.10I=imread('lf.bmp');%读入一幅原始图像I=rgb2gray(I);I1=imadjust(I,stretchlim(I),[]); %对原始图像进行灰度变换的增强figure,imshow(I1);% 打开一个新的图形窗口并显示增强后的图像B=im2bw(I1,graythresh(I1)); %对增强后的图像进行阈值处理figure,imshow(B); %打开一个新的图形窗口并显示阈值处理后的图像3.11[X,Y,Z]=peaks(80);n=8;%获得peaks图形数据，给定等高线分级数n subplot(1,2,1),pcolor(X,Y,Z)%伪彩色图colormap('jet'),shading interp%色图矩阵采用蓝头红尾饱和色，浓淡插值hold on,C=contour(X,Y,Z,n, 'k');%图上用黑虚线画等位线clabel(C) %以+引导水平放置标识Ccolorbar; %显示垂直色标尺hold off,subplot(1,2,2) %画图结束在右侧显示下图[C,h]=contourf('v6',X,Y,Z,n, 'k'); %黑虚线画填色等位线clabel(C,h) %沿线布置标识C3.12I=imread('hua.bmp');I=rgb2gray(I);I1 =imadjust(I,stretchlim(I),[]); %把图像灰度从原始范围拉伸到全部范围imshow(I1); %显示灰度增强后的图像colorbar('vert'); %显示颜色条colormap('jet');3.13subplot(1,3,1); sphere(30);%打开一个图形窗口并将其在水平方向分成三个子窗口，在左侧子窗口显示一球体( SPHERE(N) Generate sphere by three (N+1)-by-(N+1) matrices )axis square; %横轴和纵轴的比例是1:1shading flat %采用分层明暗处理grid;title('Flat Shading'); %除去网格，显示标题subplot(1,3,2);sphere(30);%在中间的子窗口显示一个球体axis square; %横轴和纵轴的比例是1:1shading faceted %采用小方块明暗处理grid;title('Faceted Shading');% 除去网格，显示标题subplot(1,3,3);sphere(30); %在最右边的子窗口显示一个球体axis square; %横轴和纵轴的比例是1:1shading interp %采用插补明暗处理grid;title('Interpolated Shading'); %除去网格，显示标题3.14rgb=imread('hua.bmp'); %读入一幅真彩色图像subplot(1,2,1);imshow(rgb);% 显示该真彩色图像f_rgb(:,:,1)=rgb(:,:,3);%把蓝色分量变为红色分量f_rgb(:,:,2)=rgb(:,:,1);%把红色分量变为绿色分量f_rgb(:,:,3)=rgb(:,:,2);%把绿色分量变为蓝色分量subplot(1,2,2);imshow(f_rgb);%显示假彩色处理后图像。