数字图像处理-matlab代码

subplot(1,2,2),imhist(I)%输出原图直方title('原始图像直方图')%在原图直方图上加标题I=imread('cameraman.tif');theta=K=imrotate(I,theta);%Tryvaryingtheangle,theta.figure,imshow(K)I=imread('cameraman.tif');程序实现如下： 程序实现如下：axison; axis J=imadjust(I1,[0.10.5局部拉伸，把[0.10.5]内的灰度拉伸为[01]grid axis K=imadjust(I1,[0.30.7局部拉伸，把[0.30.7]内的灰度拉伸为[01]grid axis 程序实现如下：title('灰度图像');grid axison; title('对数变换图像');grid axis 程序实现如下：用实现领域平均法抑制噪声程序I1=imnoise(I,'salt&pepper',0.02);title('添加椒盐噪声的图像 %进行3*3模板平滑滤波 进行 模板平滑滤波 %进行7*7模板平滑滤波 %进行9*9模板平滑滤波subplot(233),imshow(k1);title('3*3模板平滑滤波');subplot(234),imshow(k2);title('5*5模板平滑滤波');subplot(235),imshow(k3);title('7*7模板平滑滤波');用实现中值滤波程序如下 %进行3*3模板中值滤波 %进行5*5模板中值滤波 %进行7*7模板中值滤波 %进行9*9模板中值滤波subplot(233),imshow(k1);title('3*3模板中值滤波');subplot(234),imshow(k2);title('5*5模板中值滤波');subplot(235),imshow(k3);title('7*7模板中值滤波');subplot(236),imshow(k4);title('9*9模板中值滤波');用Sobel算子和对图像锐化：grid axison; grid axison; %选择sobel算子 %卷积运算title('sobel算子锐化图像');grid axis h=[010,1-41,01 %算 grid axis 用实现如下：grid axison; grid axison; grid axison; title('sobel算子分割结果');grid axison; title('Prewitt');grid axis LOG用程序实现如下：title('log算子分割结果Canny用程序实现如下边界（bwtraceboundary函数clearall %将彩像转化灰度图像 BW=im2bw(I1,threshold); col=round(dim(2)/2)-90;%计算起始点列坐标 holdon;plot(contour(:,2),contour(:,1g','LineWidth'2);HoughI=imread('xian.bmp');gridon;axistitle('prewittgridon;axisaxison,axisnormal,holdon;gridon;axison;holdon;fork=1:length(lines)用实现直方图阈值法：title('灰度图像')grid axis fork=0:255 %GP中相应 ylabel('出现概率')title('150的分割图像')grid axison; title('200的分割图像')grid axis 自动阈值法：Otsu 实现Otsu算法clearallgrid axison; title('Otsu法阈值分割图像')grid axis grid axison; title('膨胀后图像');grid axis 实现腐蚀操作 grid axison; grid axis 用实现开启和闭合操作 axison; axis I2=imopen(I1,se);%开启操作I3=imclose(I1,se);%闭合操作axison; axis axison; axison; I2=imopen(I1,se);%开启操作I3=imclose(I1,se);%闭合操作axison; axison; axison; axis 利用实现如下： grid axis grid axison; gridon;ax