直线扫描转换代码-MATLAB_第1页
直线扫描转换代码-MATLAB_第2页
直线扫描转换代码-MATLAB_第3页
直线扫描转换代码-MATLAB_第4页
全文预览已结束

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、精选优质文档-倾情为你奉上一、直线扫描转换%file1.m文件,这是主文件%x=0,y=0,r=1;%Bresenhamcircle(x,y,r);clear all;clc;p1=19 820;p2=-1609 2370;%Line_Scanner_DDA(p1(1),p1(2),p2(1),p2(2);%Line_Scanner_Bresenham(p1,p2);line_Bresenham(p1,p2);line_DDA(p1,p2);%x=3.12 -3.6;%ceil(x)%line_Bresenham.m文件,实现直线扫描转换的bresenham算法function line_Br

2、esenham(p1,p2)XX=p2(1)-p1(1);YY=p2(2)-p1(2);XYSwitch=0;%表示x和y是否交换,1表示需要交换YSign=1;%表示YY的符号,-1表示经过了y->-y的转换if XX*YY=0 if YY=0 t=sort(p1(1),p2(1); Point(:,1)=(t(1):t(2)' Point(:,2)=zeros(t(2)-t(1)+1,1)+p1(2); else t=sort(p1(2),p2(2); Point(:,2)=(t(1):t(2)' Point(:,1)=zeros(t(2)-t(1)+1,1)+p1(

3、1); endelse %如果斜率的绝对值大于1,则交换x和y if abs(YY)>abs(XX) XYSwitch=1; t(1)=p1(1); p1(1)=p1(2); p1(2)=t(1); t(1)=p2(1); p2(1)=p2(2); p2(2)=t(1); end %如果x1>x2则交换p1和p2 if p1(1)>p2(1) t=p1(1),p1(2); p1=p2(1),p2(2); p2=t; end XX=p2(1)-p1(1); YY=p2(2)-p1(2); %如果y1>y2,则做y->-y的转换 if YY<0 p1(2)=-p

4、1(2); p2(2)=-p2(2); YSign=-1; end YY=p2(2)-p1(2); %开始计算点 d=XX-2*YY; P=p1(1),p1(2);k=1; while P(1)<=p2(1) Point(k,:)=P; if d<0 P=P+ones(1,2); d=d+2*XX-2*YY; else P(1)=P(1)+1; d=d-2*YY; end k=k+1; end %恢复点 if YSign=-1 Point(:,2)=-Point(:,2); end if XYSwitch=1 t=Point(:,1); Point(:,1)=Point(:,2);

5、 Point(:,2)=t; endend%绘制点figure();scatter(Point(:,1),Point(:,2),2,'.');legend('Bresenham算法');% line_DDA.m文件,实现直线扫描转换的数值微分算法function line_DDA(P1,P2)p1=P1(1),P1(2);p2=P2(1),P2(2);XY=p2-p1;%x轴和y轴增量XYSwitch=0;%标记P1和P2的x和y是否交换过,1表示需要交换YSign=1;%标记y轴增量符号,-1表示需要进行y->-y的转换if XY(1)*XY(2)=0

6、if XY(1)=0 t=sort(p1(2),p2(2); Point(:,2)=(t(1):t(2)' Point(:,1)=zeros(t(2)-t(1)+1,1)+p1(1); else t=sort(p1(1),p2(1); Point(:,1)=(t(1):t(2)' Point(:,2)=zeros(t(2)-t(1)+1,1)+p1(2); endelse %如果斜率的绝对值大于1,则交换x和y if abs(XY(2)>abs(XY(1) XYSwitch=1; t(1)=p1(1); p1(1)=p1(2); p1(2)=t(1); t(1)=p2(1

7、); p2(1)=p2(2); p2(2)=t(1); end %如果x1>x2,则交换p1和p2 if p1(1)>p2(1) t=p1; p1=p2; p2=t; end XY=p2-p1; %如果y1>y2,则做y->-y的转换 if XY(2)<0 p1(2)=-p1(2); p2(2)=-p2(2); YSign=-1; end XY=p2-p1; %开始计算 k=XY(2)/XY(1); Point(:,1)=(p1(1):p2(1)' Point(:,2)=round(0:abs(XY(1)'*k+p1(2); %恢复变换 if YSign=-1 Point(:,2)=-Point(:,2); end if XYSwitch=1 t=Point(:,1); Point(:,1)=Point(:,2); Po

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论