图像处理与分析-报告.docVIP

中 国 地 质 大 学 研究生课程论文封面 课程名称 图像处理与分析 教师姓名 研究生姓名 研究生学号 研究生专业 所在院 类 别: 日 期:  目录 1 Matlab图像处理操作 1 1.1 图像文件读取 1 1.2 图像变换 2 1.2.1 图像旋转 2 1.2.2 图像平移 3 1.2.3 图像缩放 4 1.3 图像增强 6 1.3.1 均值滤波和中值滤波 6 1.4 边缘提取 9 1.4.1 各种边缘提取方法 9 2 VC图像处理操作 12 2.1 直方图规定化的单映射与组映射 12 2.2 Sobel算子检测边缘 17 2.3 Susan算子检测边缘 20 2.4 Snake算子边缘检测 23 Matlab图像处理操作 图像文件读取 %图像读取 a=imread(cameraman.tif); subplot(1,3,1),imshow(a); title(tif 图像); a=imread(ImgTrans.bmp); subplot(1,3,2),imshow(a); title(bmp 图像); a=imread(greens.jpg); subplot(1,3,3),imshow(a); title(jpg 图像); 读取三种不同格式的图像 图像变换 图像旋转 主程序代码： %图像旋转 clear; %读取图像 Image=imread(ImgTrans.bmp);% angle=30; im1=rgb2gray(Image); subplot(1,2,1),imshow(im1); title(原始图像); %根据角度旋转图像 im2=FunRotate(im1,angle); %绘制新图 subplot(1,2,2),imshow(im2); title(旋转图像); 图像旋转函数代码： %旋转函数 function im2=FunRotate(im1,angle) pai=3.14; roteang=angle/180*pai; [im1X,im1Y]=size(im1); %计算旋转矩阵 rotmatrix=[cos(roteang),-sin(roteang);sin(roteang),cos(roteang)]; %计算反旋转矩阵 rrotmatrix=[cos(-roteang),-sin(-roteang);sin(-roteang),cos(-roteang)]; %计算旋转后的图像大小 im2XY1=[im1X,im1Y]*rotmatrix; im2XY2=[im1X,1]*rotmatrix; im2XY3=[1,1]*rotmatrix; im2XY4=[1,im1Y]*rotmatrix; im2p11=[min([im2XY1(1),im2XY2(1),im2XY3(1),im2XY4(1)]),min([im2XY1(2),im2XY2(2),im2XY3(2),im2XY4(2)])]; %图像框架 im2X=round(max([im2XY1(1),im2XY2(1),im2XY3(1),im2XY4(1)])-im2p11(1)); im2Y=round(max([im2XY1(2),im2XY2(2),im2XY3(2),im2XY4(2)])-im2p11(2)); %建立新图像 im2=uint8(zeros(im2X,im2Y))+255; %对于每个像素点，进行上文所述的计算。 for i=1:im2X for j=1:im2Y pij=round([i+im2p11(1)-1,j+im2p11(2)-1]*rrotmatrix); if (pij(1)=1)(pij(2)=1)(pij(1)=im1X)(pij(2)=im1Y) im2(i,j)=im1(pij(1),pij(2)); end end end; 旋转图像 图像平移 %图像平移 ImgOrg = imread(ImgTrans.bmp); J = double(ImgOrg); H = size(ImgOrg); I_moveresult = zeros(H); x_move = 50; y_move = 50; %移动参数设置 I_moveresult(x_move+1:H(1),y_move+1:H(2),1:H(3)) = J(1:H(1)-x_move,1:H(2)-x_move,1:H(3)); subplot(1,2,1),imshow(