Matlab图像处理代码及仿真结果.docx

Matlab图像处理代码及仿真结果.docx

ID:58149789

大小:1.88 MB

页数:23页

时间:2020-04-25

Matlab图像处理代码及仿真结果.docx_第1页
Matlab图像处理代码及仿真结果.docx_第2页
Matlab图像处理代码及仿真结果.docx_第3页
Matlab图像处理代码及仿真结果.docx_第4页
Matlab图像处理代码及仿真结果.docx_第5页
资源描述:

《Matlab图像处理代码及仿真结果.docx》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库

1、预习报告要求:A.绘出信号,当、时、、时的信号实部和虚部图;当时呢?此时信号周期为多少?代码:functionx=reim(z)n=1:100;x=exp(z.*n)subplot(2,1,1)plot(n,real(x));%xµÄʵ²¿title('realpartofx');subplot(2,1,2)plot(n,imag(x));%xµÄÐ鲿title('imaginarypartofx');enda=reim(-(1./12)+(pi./6).*i)a=reim((1./12)+(pi./6).*i)a=reim(1./12)a=reim(2+(pi./6).*i)a

2、=reim((pi./6).*i)B.绘出信号的频率是多少?周期是多少?产生一个数字频率为0.9的正弦序列,并显示该信号,说明其周期。频率是0.1;周期是10代码:functionx=s(a,b)n=1:100;x=a.*sin(b.*n)plot(n,x);a=s(1.5,2.*pi.*0.1)2、使用帮助功能学习square(方波),sawtooth(锯齿波)和sinc函数,并绘图1.实验目的:熟悉图像处理常用函数和方法;培养通过查阅文献解决问题的能力。2.实验内容编程实现对该图像的任意比例的放大及缩小编程实现对该图像的任意角度旋转解决锯齿等图像不平滑问题(邻近法,双线性法)将自己

3、实现的方法与matlab的方法进行比较针对黑白、灰度、彩色三类图像分别测试3.实验步骤与分析考虑到在平面直角坐标系进行坐标点的平移变换会比较直观易懂,所以先把图像原点(这里取为图像的中心点)转化为直角坐标系原点,转换矩阵和代码如下:change=[100;0-10;-0.5.*h0.5.*k1];进行图像处理操作后,要把图像的坐标原点转换回来,对图像进行还原,还原转换矩阵和代码如下:rechange=[100;0-10;0.5.*R0.5.*C1];3.1图像的旋转原图:对原点变换完成后的图像进行旋转操作,旋转矩阵和代码如下:tras=[cos(ang)-sin(ang)0;sin(a

4、ng)cos(ang)0;001];3.1.1正映射方法旋转处理图像rotimg.m下面采用正映射的方法对处理后的图像进行灰度值赋值:完整代码如下:functiona=rotimg(a,b)A=imread(a);[R,C,O]=size(A);ang=(b./180).*pi;h=double(abs(R.*cos(ang))+abs(C.*sin(ang)));k=double(abs(C.*cos(ang))+abs(R.*sin(ang)));o=O;r=uint8(h);c=uint8(k);change=[100;0-10;-0.5.*R0.5.*C1];tras=[cos

5、(ang)-sin(ang)0;sin(ang)cos(ang)0;001];rechange=[100;0-10;0.5.*h0.5.*k1];B=zeros(r,c,o);fori=1:R  forj=1:C      temp=[ij1];      temp=temp*change;      temp=temp*tras;      temp=temp*rechange;      x=uint8(temp(1,1));      y=uint8(temp(1,2));      if(x<=r)&&(y<=c)&&(x>=1)&&(y>=1)        B(x,y,:)

6、=A(i,j,:);      end  endend;imshow(uint8(B));运行结果:a=rotimg('lolicon.jpg',45)可以观察到图像的旋转和图像的大小比较符合要求,但是放大看会有有规律的像素点空洞:这是由于正向映射,代码x=uint8(temp(1,1));y=uint8(temp(1,2));会使部分像素点赋值重合,造成空洞点的问题3.1.2最邻近插值法旋转处理图片rotimg_1.m如果采用反向映射的方法,就可以保证每个像素点都有灰度值,不会出现像素点空洞的情况,完整代码如下:functiona=rotimg_1(a,b)A=imread(a);[

7、R,C,O]=size(A);ang=(-b./180).*pi;h=double(abs(R.*cos(ang))+abs(C.*sin(ang)));k=double(abs(C.*cos(ang))+abs(R.*sin(ang)));o=O;r=uint8(h);c=uint8(k);change=[100;0-10;-0.5.*h0.5.*k1];tras=[cos(ang)-sin(ang)0;sin(ang)cos(ang)0;0

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。