计算机图形学专业课程设计实验报告.doc

计算机图形学试验汇报班级 计算机科学和技术 姓名 学号 年 6 月 2 日试验一 基础图形生成算法一、试验目标： 1、掌握中点Bresenham绘制直线原理； 2、设计中点Bresenham算法； 3、掌握八分法中点Bresenham算法绘制圆原理； 4、设计八分法绘制圆中点Bresenham算法； 5、掌握绘制1/4椭圆弧上半部分和下半部分中点Bresenham算法原理； 6、掌握下半部分椭圆偏差判别式初始值计算方法； 7、设计顺时针四分法绘制椭圆中点Bresenham算法。二、试验过程：1、 试验描述 试验1：使用中点Bresenham算法绘制斜率为0=k=1直线。 试验2：使用中点Bresenham算法绘制圆心在屏幕用户区中心圆。 试验3：使用中点Bresenham算法绘制圆心在屏幕用户区中心椭圆。2、试验过程 1）用MFC(exe)建立一个单文档工程； 2）编写对话框，生成对应对象，设置对应变量； 3）在类CLineView中申明对应函数，并在相关cpp文件中实现； 4）在OnDraw（）函数里调用函数实现绘制直线、圆、椭圆； 5）运行程序，输入对应值，绘制出图形。三、 源代码试验1：直线中点Bresenham算法1./ cline.cpp : implementation file/ cline dialogcline:cline(CWnd* pParent /*=NULL*/): CDialog(cline:IDD, pParent)/AFX_DATA_INIT(cline)m_x0 = 0;m_y0 = 0;m_x1 = 0;m_y1 = 0;/AFX_DATA_INITvoid cline:DoDataExchange(CDataExchange* pDX)CDialog:DoDataExchange(pDX);/AFX_DATA_MAP(cline)DDX_Text(pDX, IDC_x0, m_x0);DDX_Text(pDX, IDC_y0, m_y0);DDX_Text(pDX, IDC_x1, m_x1);DDX_Text(pDX, IDC_y1, m_y1);/AFX_DATA_MAPBEGIN_MESSAGE_MAP(cline, CDialog)/AFX_MSG_MAP(cline)/AFX_MSG_MAPEND_MESSAGE_MAP()2、 / LineView.hclass CLineView : public CViewpublic:CLineDoc* GetDocument();.void Mbline(double,double,double,double); /直线中点Bresenham函数.3、/ Line.cpp /*直线中点Bresenham函数*/void CLineView:Mbline(double x0, double y0, double x1, double y1) CClientDC dc(this); COLORREF rgb=RGB(0,0,255); /定义直线颜色为蓝色double x,y,d,k;x=x0; y=y0; k=(y1-y0)/(x1-x0); d=0.5-k;for(x=x0;x=x1;x+)dc.SetPixel(int)x,(int)y,rgb);if(d0)y+;d+=1-k;elsed-=k;4、/LineView.cppvoid CLineView:OnDraw(CDC* pDC)CLineDoc* pDoc = GetDocument();ASSERT_VALID(pDoc);/ TODO: add draw code for native data herecline a;a.DoModal();/初始化CLineView:Mbline(a.m_x0,a.m_y0,a.m_x1,a.m_y1);试验2：圆中点Bresenham算法 1、/cricle.cpp/ Ccricle dialogCcricle:Ccricle(CWnd* pParent /*=NULL*/): CDialog(Ccricle:IDD, pParent)/AFX_DATA_INIT(Ccricle)m_r = 0;/AFX_DATA_INITvoid Ccricle:DoDataExchange(CDataExchange* pDX)CDialog:DoDataExchange(pDX);/AFX_DATA_MAP(Ccricle)DDX_Text(pDX, r_EDIT, m_r);/AFX_DATA_MAP 2、/CcircleView.hclass CCcircleView : public CView.public:CCcircleDoc* GetDocument();void CirclePoint(double,double); /八分法画圆函数void Mbcircle(double); /圆中点Bresenham函数. 3、/CcircleView.cppvoid CCcircleView:OnDraw(CDC* pDC)CCcircleDoc* pDoc = GetDocument();ASSERT_VALID(pDoc);/ TODO: add draw code for native data hereCcricle r;r.DoModal();CCcircleView:Mbcircle(r.m_r);/画圆4、/CcircleView.cpp/*八分法画圆*/void CCcircleView:CirclePoint(double x,double y)CClientDC dc(this);COLORREF rgb=RGB(0,0,255);dc.SetPixel(int)(300+x),(int)(300+y),rgb);dc.SetPixel(int)(300-x),(int)(300+y),rgb);dc.SetPixel(int)(300+x),(int)(300-y),rgb);dc.SetPixel(int)(300-x),(int)(300-y),rgb);dc.SetPixel(int)(300+y),(int)(300+x),rgb);dc.SetPixel(int)(300-y),(int)(300+x),rgb);dc.SetPixel(int)(300+y),(int)(300-x),rgb);dc.SetPixel(int)(300-y),(int)(300-x),rgb);/*圆中点Bresenham函数*/void CCcircleView:Mbcircle(double r)double x,y,d;COLORREF rgb=RGB(0,0,255);d=1.25-r;x=0; y=r;for(x=0;xy;x+)CirclePoint(x,y); /调用八分法画圆子函数if(d0)d+=2*x+3;elsed+=2*(x-y)+5;y-;试验3：椭圆中点Bresenham算法1、/ellipse1.cpp/ Cellipse dialogCellipse:Cellipse(CWnd* pParent /*=NULL*/): CDialog(Cellipse:IDD, pParent)/AFX_DATA_INIT(Cellipse)m_a = 0;m_b = 0;/AFX_DATA_INITvoid Cellipse:DoDataExchange(CDataExchange* pDX)CDialog:DoDataExchange(pDX);/AFX_DATA_MAP(Cellipse)DDX_Text(pDX, IDC_EDIT1, m_a);DDX_Text(pDX, IDC_EDIT2, m_b);/AFX_DATA_MAP2、/EllipseView.hclass CEllipseView : public CView.public:CEllipseDoc* GetDocument();void EllipsePoint(double,double); /四分法画椭圆void Mbellipse(double a, double b); /椭圆中点Bresenham函数.3、/Ellipse.cpp/*四分法画椭圆*/void CEllipseView:EllipsePoint(double x,double y)CClientDC dc(this);COLORREF rgb=RGB(0,0,255); dc.SetPixel(int)(300+x),(int)(300+y),rgb);dc.SetPixel(int)(300-x),(int)(300+y),rgb);dc.SetPixel(int)(300+x),(int)(300-y),rgb);dc.SetPixel(int)(300-x),(int)(300-y),rgb);/*椭圆中点Bresenham函数*/void CEllipseView:Mbellipse(double a, double b)double x,y,d1,d2;x=0;y=b;d1=b*b+a*a*(-b+0.25);EllipsePoint(x,y);while(b*b*(x+1)a*a*(y-0.5)/椭圆AC弧段if(d10)if(d20)d2+=b*b*(2*x+2)+a*a*(-2*y+3);x+;elsed2+=a*a*(-2*y+3);y-;EllipsePoint(x,y);4、/EllipseView.cppvoid CEllipseView:OnDraw(CDC* pDC)CEllipseDoc* pDoc = GetDocument();ASSERT_VALID(pDoc);/ TODO: add draw code for native data hereCellipse el;el.DoModal();/初始化CEllipseView:Mbellipse(el.m_a, el.m_b);/画椭圆四、 实结果验试验1：直线中点Bresenham算法试验2：圆中点Bresenham算法试验3：椭圆中点Bresenham算法试验二 有效边表填充算法一、试验目标：1、设计有效边表结点和边表结点数据结构；2、设计有效边表填充算法；3、编程实现有效边表填充算法。二、 试验过程：1、试验描述 下图1 所表示多边形覆盖了 12 条扫描线，共有 7 个顶点和 7 条边。7 个顶点分别为：P0（7，8） ，P1（3，12） ，P2（1，7） ，P3（3，1）, P4（6，5）, P5（8，1）, P6（12，9）。在 1024768 显示分辩率下，将多边形顶点放大为 P0（500,400） ，P1（350，600） ，P2（250，350），P3（350，50）, P4（500，250）, P5（600，50）, P6（800，450）。请使用有效边表算法填充该多边形。2、 试验过程 （1）建立AET类和Bucket类； （2）初始化桶，并在建立桶结点时为其表示扫描线初始化为带头结点链表； （3）对每个桶结点进行循环，将桶内每个结点边表合并为有效边表，并进行有效边表循环； （4）根据扫描线从小到大移动次序，计算目前扫描线和多边形各边交点，然后把这些交点按X值递增次序进行排序，配对，以确定填充区间； （5）用指定颜色点亮填充区间内全部像素，即完成填充工作。三、 源代码1、/AET.hclass AET public:AET();virtual AET();double x;int yMax;double k;/替换1/kAET *next;/AET.cppAET:AET()AET:AET()2、 /Bucket.h#include AET.hclass Bucket public:Bucket();virtual Bucket();int ScanLine;AET *p;/桶上边表指针Bucket *next;/ Bucket.cppBucket:Bucket()Bucket:Bucket()3、/test2.cppvoid CTest2View:Et()/结构边表函数for(int i=0;iPointi.y)/终点比起点高while(CurrentB-Scanline!=Pointi.y)/在桶内寻求该边yMinCurrentB=CurrentB-next;/移到下一个桶结点Ei.x=Pointi.x;/计算AET表值Ei.yMax=Pointj.y;Ei.k=double(Pointj.x-Pointi.x)/(Pointj.y-Pointi.y);/代表1/kEi.next=NULL;CurrentE=CurrentB-p;/取得桶上连接边表地址if(CurrentB-p=NULL)CurrentE=&Ei;/边表起始地址CurrentB-p=CurrentE;/第一个边表直接连接到对应桶中elsewhile(CurrentE-next!=NULL)/假如目前边已连有边结点CurrentE=CurrentE-next;/移动指针到目前边最终一个结点CurrentE-next=&Ei;/把目前边接上去if(Pointj.yScanline!=Pointj.y)/在桶内寻求该边yMinCurrentB=CurrentB-next;/移到下一个桶结点Ei.x=Pointj.x;/计算AET表值Ei.yMax=Pointi.y;Ei.k=double(Pointi.x-Pointj.x)/(Pointi.y-Pointj.y);/代表1/kEi.next=NULL;CurrentE=CurrentB-p;/取得桶上连接边表地址if(CurrentE=NULL)CurrentE=&Ei;/边表起始地址CurrentB-p=CurrentE;/第一个边表直接连接到对应桶中elsewhile(CurrentE-next!=NULL)/假如目前边已连有边结点CurrentE=CurrentE-next;/移动指针到目前边最终一个结点CurrentE-next=&Ei;/把目前边接上去CurrentB=NULL;CurrentE=NULL;void CTest2View:AddEdge(AET *NewEdge)/插入临时边表函数T1=HeadE;if(T1=NULL)/边表为空，将边表置为TempEdgeT1=NewEdge;HeadE=T1;elsewhile(T1-next!=NULL)/边表不为空，将TempEdge连在该边以后T1=T1-next;T1-next=NewEdge;void CTest2View:EdegeOrder()/对边表进行排序函数T1=HeadE;if(T1=NULL)return;if(T1-next=NULL)/假如该边表没有再连边表return;/桶结点只有一条边，不需要排序elseif(T1-next-xx)/边表按x值排序T2=T1-next;T1-next=T2-next;T2-next=T1;HeadE=T2; T2=HeadE; T1=HeadE-next; while(T1-next!=NULL)/继续两两比较相连边表x值，进行排序 if(T1-next-xx) T2-next=T1-next; T1-next=T1-next-next; T2-next-next=T1;T2=T2-next; else T2=T1; T1=T1-next; void CTest2View:PolygonFill()/多边形填充函数CClientDC dc(this);COLORREF rgb=RGB(0,0,255); HeadE=NULL;for(CurrentB=HeadB;CurrentB!=NULL;CurrentB=CurrentB-next)/访问全部桶结点for(CurrentE=CurrentB-p;CurrentE!=NULL;CurrentE=CurrentE-next)/桶中全部结点AET *TempEdge= new AET;TempEdge-x=CurrentE-x;TempEdge-yMax=CurrentE-yMax; TempEdge-k=CurrentE-k; TempEdge-next=NULL; AddEdge(TempEdge);/将该边插入临时Aet表EdegeOrder();/边表根据x递增次序存放 T1=HeadE;/依据ymax抛弃扫描完边结点if(T1=NULL)return;while(CurrentB-Scanline=T1-yMax)/放弃该结点，Aet表指针后移，下闭上开T1=T1-next;HeadE=T1;if(HeadE=NULL)return;if(T1-next!=NULL)T2=T1;T1=T2-next;while(T1!=NULL)if(CurrentB-Scanline=T1-yMax)/跳过一个结点 T2-next=T1-next; T1-next=NULL; T1=T2-next;elseT2=T1;T1=T2-next;BOOL In=false;/设置一个BOOL变量In，初始值为假double xb,xe;/扫描线起点和终点for(T1=HeadE;T1!=NULL;T1=T1-next)/填充扫描线和多边形相交区间if(In=false)xb=T1-x;In=true;/，每访问一个结点，把In值取反一次else/假如In值为真，则填充从目前结点x值开始到下一个结点x值结束区间xe=T1-x-1;/左闭右开CClientDC dc(this);for(double x=xb;xScanline,rgb);/填充 Sleep(1);/延时1ms，提升填充过程可视性 In=FALSE;for(T1=HeadE;T1!=NULL;T1=T1-next)/边连贯性 T1-x=T1-x+T1-k;/x=x+1/kdelete HeadB;delete CurrentB;delete CurrentE;delete HeadE;/*初始化桶*void CTest2View:CreatBucket() int ScanMax,ScanMin;ScanMax = ScanMin = Point0.y;for(int i=1;iNumber;i+)if(Pointi.yScanMax)ScanMax = Pointi.y; /扫描线最大值for( i=ScanMin;iScanline = ScanMin;CurrentB-p = NULL; /没有连接边链表CurrentB-next = NULL;else /建立桶其它节点CurrentB-next = new Bucket; /新建一个桶节点CurrentB = CurrentB-next; /使CurrentB指向新建桶节点CurrentB-Scanline = i;CurrentB-p = NULL; /没有连接边链表CurrentB-next = NULL;/*菜单函数*/*void CTest2View:OnMenuAET() AfxGetMainWnd()-SetWindowText(多边形有效边表填充算法);/显示标题CColorDialog ccd(GetColor);if(ccd.DoModal()=IDOK)/调用调色板选择前景色GetColor=ccd.GetColor();getchar();RedrawWindow();/刷新屏幕CreatBucket();/初始化桶Et();/建立边表PolygonFill();/多边形填充*/4、/test2View.cppvoid CTest2View:OnDraw(CDC* pDC)CTest2Doc* pDoc = GetDocument();ASSERT_VALID(pDoc);/ TODO: add draw code for native data here/-pDC-Polygon(Point,7);/绘制多边形/输出多边形顶点编号pDC-TextOut(550,410,P0);pDC-TextOut(350,600,P1); pDC-TextOut(230,340,P2);pDC-TextOut(350,30,P3);pDC-TextOut(490,220,P4);pDC-TextOut(600,30,P5);pDC-TextOut(805,450,P6);getchar();RedrawWindow();/刷新屏幕CreatBucket();/初始化桶Et();/建立边表PolygonFill();/多边形填充/OnMenuAET();四、 试验结果