第3章 5 裁剪算法
2016新编Weiler-Atherton任意多边形裁剪算法

Weiler-Atherton任意多边形裁剪Sutherland-Hodgeman算法解决了裁剪窗口为凸多边形窗口的问题,但一些应用需要涉及任意多边形窗口(含凹多边形窗口)的裁剪。
Weiler-Atherton多边形裁剪算法正是满足这种要求的算法。
一、Weiler-Atherton任意多边形裁剪算法描述:在算法中,裁剪窗口、被裁剪多边形可以是任意多边形:凸的、凹的(内角大于180o)、甚至是带有内环的(子区),见下图。
裁剪窗口和被裁剪多边形处于完全对等的地位,这里我们称:1、被裁剪多边形为主多边形,记为A;2、裁剪窗口为裁剪多边形,记为B。
主多边形A和裁剪多边形B的边界将整个二维平面分成了四个区域:1、A∩B(交:属于A且属于B);2、A-B(差:属于A不属于B);3、B-A(差:属于B不属于A);4、A∪B(并:属于A或属于B,取反;即:不属于A且不属于B)。
内裁剪即通常意义上的裁剪,取图元位于窗口之内的部分,结果为A∩B。
外裁剪取图元位于窗口之外的部分,结果为A-B。
观察右图不难发现裁剪结果区域的边界由被裁剪多边形的部分边界和裁剪窗口的部分边界两部分构成,并且在交点处边界发生交替,即由被裁剪多边形的边界转至裁剪窗口的边界,或者反之。
由于多边形构成一个封闭的区域,所以,如果被裁剪多边形和裁剪窗口有交点,则交点成对出现。
这些交点分成两类:一类称“入”点,即被裁剪多边形由此点进入裁剪窗口,如图中a、c、e;一类称“出”点,即被裁剪多边形由此点离开裁剪窗口,如图中b、d、f。
二、Weiler-Atherton任意多边形裁剪算法思想:假设被裁剪多边形和裁剪窗口的顶点序列都按顺时针方向排列。
当两个多边形相交时,交点必然成对出现,其中一个是从被裁剪多边形进入裁剪窗口的交点,称为“入点”,另一个是从被裁剪多边形离开裁剪窗口的交点,称为“出点”。
算法从被裁剪多边形的一个入点开始,碰到入点,沿着被裁剪多边形按顺时针方向搜集顶点序列;而当遇到出点时,则沿着裁剪窗口按顺时针方向搜集顶点序列。
简单多边形裁剪算法

摘 要 :为 了尽量 降低任意 多边形裁剪复杂度,提 出了一种基 于多边形顶点遍历 的简单 多边形裁 剪算 法。该 算法将 多边形
交 点插 入 到裁 剪 多边 形 和 被 裁 减 多边 形 顶 点 矢 量数 组 中 ,通 过 记 录 交 点及 其 前 驱 、 后 继 信 息 ,可 快 速 生 成 结 果 多边 形 。其
2 .C h i n a C e n t r e f o r R e s o u r c e s S a t e l l i t e Da t a a n d Ap p l i c a t i o n , B e i j i n g 1 0 0 0 9 4 , C h i n a ;3 .C h i n a C o a l Te c h n o l o g i e s G r o u p or C p o r a t i o n ,B e i j i n g 1 0 0 0 1 3 , C h i n a ; 4 .B e i j i n g I n s t i t u t e o f A p p l i e d Me t e o r o l o g y , B e i j i n g 1 0 0 0 2 9 , C h i n a )
中,时间复杂度 为 ( ) ( ( m+k ) ×k ) ,m 是两 多边形 中顶点数较 大者 ,k是 两 多边形的 交点数 。该算 法简化 了交点的数据 结 构 ,节省 了存储 空间,降低 了算法的时间复杂度 ,具有 简单 、易于编程 实现 、运行效 率高的特点 。
关 键 词 : 多边 形 裁 剪 ; 交点 ;前 驱 ;后 继 ;矢 量数 组 中图 法 分 类 号 :TP 3 9 1 文 献 标 识 号 :A 文章 编 号 :1 0 0 0 — 7 0 2 4( 2 0 1 4 )0 1 - 0 1 9 2 - 0 6
计算机图形学(简单多边形裁剪算法)

简单多边形裁剪算法摘要:多边形裁剪算法与线性裁剪算法具有更广泛的实用意义,因此它是目前裁剪研究的主要课题。
本文主要介绍了一种基于多边形顶点遍历的简单多边形裁剪算法,它有效降低了任意多边形裁剪复杂度。
通过记录交点及其前驱、后继信息,生成结果多边形,该算法简化了交点的数据结构,节省了存储空间,降低了算法的时间复杂度,具有简单、易于编程实现、运行效率高的特点。
关键词:多边形裁剪;交点;前驱;后继;矢量数组一、技术主题的基本原理简单多边形裁剪算法综合考虑现有多边形裁剪算法的优缺点,它是一种基于多边形顶点遍历来实现简单多边形裁剪工作的。
其主要的原理是遍历多边形并把多边形分解为边界的线段逐段进行裁剪,输出结果多边形。
二、发展研究现状近年来,随着遥感绘图、CAD辅助设计、图象识别处理技术的发展,图形裁剪算法从最初在二维平面上线和图形的裁剪扩展到三维空间里体和场的裁剪,国内外相继提出不少行之有效的算法,但越来越复杂的图形和计算也对算法的速度和适用性提出了越来越高的要求。
因此,不断简化算法的实现过程,完善细节处理,满足大量任意多边形的裁剪也就成了当今算法研究的焦点之一。
以往多边形裁剪算法不是要求剪裁多边形是矩形,就是必须判断多边形顶点的顺时针和逆时针性,即存在不实用或者是增加了多边形裁剪算法的难度。
为了解决现在的问题,我们研究现在的新多边形算法,其中,裁剪多边形和被裁剪多边形都可以是一般多边形,且不需要规定多边形输入方向。
它采用矢量数组结构,只需遍历剪裁多边形和被裁剪多边形顶点即完成多边形的裁剪,具有算法简单、运行效率高的特点。
三、新算法设计1、算法的思想本算法是为了尽量降低任意多边形裁剪算法复杂度而提出的,其主要思想是采用矢量数组结构来遍历裁剪多边形和被裁多边形顶点,记录裁剪多边形和被裁减多边形交点及其前驱、后继信息,并通过记录相邻交点的线段,然后通过射线法选择满足条件的线段,之后进行线段连接,输出对应的裁剪结果。
裁剪算法——cohen-sutherland算法

裁剪算法——cohen-sutherland算法实验环境:VC6.0算法思想: 延长窗⼝的边,将⼆维平⾯分成9个区域,每个区域赋予4位编码C t C b C r C l,裁剪⼀条线段P1P2时,先求出所在的区号code1,code2。
若code1=0,且code2=0,则线段P1P2在窗⼝内,应取之。
若按位与运算code1&code2,则说明两个端点同时在窗⼝的上⽅、下⽅、左⽅或右⽅,则可判断线段完全在窗⼝外,可弃之;否则,按第三种情况处理,求出线段与窗⼝某边的交点,在交点处把线段⼀分为⼆,其中必有⼀段在窗⼝外,可弃之,再对另⼀段重复上述处理。
100110001010000100000010010********* 多边形裁剪编码程序实现:#include "stdafx.h"#include<stdio.h>#include<conio.h>#include<graphics.h>#define LEFT 1#define RIGHT 2#define BOTTOM 4#define TOP 8void midpointLine(int x0,int y0,int x1,int y1,int color)//中点画线算法{int a,b,d1,d2,d,x,y;a=y0-y1;b=x1-x0;d=2*a+b;d1=2*a;d2=2*(a+b);x=x0;y=y0;putpixel(x,y,color);while(x<x1){if(d<0){x++;y++;d+=d2;}else{x++;d+=d1;}putpixel(x,y,color);}}int encode(int x,int y,int XL,int XR,int YB,int YT)//编码{int c=0;if(x<XL) c|=LEFT;if(x>XR) c|=RIGHT;if(y<YB) c|=BOTTOM;if(y>YT) c|=TOP;return c;}void C_SLineClip(int x1,int y1,int x2,int y2,int XL,int XR,int YB,int YT){int code1,code2,code,x,y;code1=encode(x1,y1,XL,XR,YB,YT);code2=encode(x2,y2,XL,XR,YB,YT);while((code1 != 0) || (code2 != 0)){if((code1 & code2) != 0){midpointLine(x1,y1,x2,y2,RGB(0, 255, 0));//如果直线在裁剪窗⼝外就⽤绿⾊画出printf("线段在窗⼝外!");return;}if(code1 != 0) code=code1;else code=code2;if((LEFT & code) != 0){x=XL;y=y1+(y2-y1)*(XL-x1)/(x2-x1);}else if((RIGHT & code) != 0){x=XR;y=y1+(y2-y1)*(XR-x1)/(x2-x1);}else if((BOTTOM & code) != 0){y=YB;x=x1+(x2-x1)*(YB-y1)/(y2-y1);}else if((TOP & code) != 0){y=YT;x=x1+(x2-x1)*(YT-y1)/(y2-y1);}if(code == code1){x1=x; y1=y; code1=encode(x,y,XL,XR,YB,YT);}else{x2=x; y2=y; code2=encode(x,y,XL,XR,YB,YT);}}midpointLine(x1,y1,x2,y2,RGB(255,0,0));//将裁减的直线⽤红⾊标注return;}int main(int argc, char* argv[]){int gdriver=DETECT,gmode;int x1=20,y1=30,x2=250,y2=300,XL=10,XR=200,YT=400,YB=30;initgraph(&gdriver,&gmode,"c:\\tc");//setbkcolor(WHITE);cleardevice();midpointLine(x1,y1,x2,y2,RGB(0,255,0));//将被裁剪直线⽤绿⾊画出rectangle(10,400,200,30);//rectangle(int left,int top,int right,int bottom);//裁剪窗⼝ C_SLineClip(x1,y1,x2,y2,XL,XR,YB,YT);// cohen sutherland算法getch();closegraph();return0;}显⽰效果:将在窗⼝内的线段设为红⾊,窗⼝外的线段设为绿⾊。
计算机图形学 第三章 二维图形的裁剪概述

3.2.3 梁友栋-Barsky裁剪算法
式中,Δx=x2-x1,Δy=y2-y1,参数u在0~1 之间取值,P(x,y)代表了该线段上的一个 点,其值由参数u确定,由公式可知,当u=0 时,该点为P1(x1,y1),当u=1时,该点 为P2(x2,y2)。如果点P(x,y)位于由 坐标(xwmin,ywmin)和(xwmax,ywmax)所 确定的窗口内,那么下式成立: xwmin≤x1+ u· Δx≤xwmax ywmin≤y1+ u· Δy≤ywmax(3-10)
3.2.1 Cohen-Sutherland算法
► Code(int ►{
x,int y,int *c)
*c=0; if(y>ymax) /*(xmin,ymin)和(xmax,ymax) 为窗口左下角、右上角坐标。*/ *c=*c|0x08; else if(y<ymin) *c=*c|0x04; if(x>xmax) *c=*c|0x02; else if(x<xmin) *c=*c|0x01; }
► 根据直线两点式方程:
►
2)
(3-
3.2 线段的裁剪
► 整理后得通用交点公式: ►
(3-3)
► ►
1、与上边界的求交公式: (3-4)
3.2 线段的裁剪
► ►
2、与下边界的求交公式:
(3-5)
►
► ►
3、与右边界的求交公式:
(3-6) 4、与左边界的求交公式:
►
(3-7)
3.2.1 Cohen-Sutherla2、判别 根据C1和C2的具体值,可以有三种情况: (1)C1=C2=0,表明两端点全在窗口内,因而 整个线段也在窗内,应予保留。 (2)C1&C2≠0(两端点代码按位作逻辑乘不为 0),即C1和C2至少有某一位同时为1,表明两端点 必定处于某一边界的同一外侧,因而整个线段全在 窗外,应予舍弃。 (3)不属于上面两种情况,均需要求交点。
计算机图形学的裁剪算法

计算机图形学的裁剪算法
计算机图形学的裁剪算法是图形学的一种重要算法,它的基本思想是将一个完整的几何图形(如线段、多边形、圆圈等)按照指定的裁剪窗口(矩形)进行裁剪,只保留在窗口内的部分,而把窗口外的部分抛弃掉。
由于裁剪算法的应用非常广泛,像图形显示系统、图形设备接口(GDI)和图形处理器(GPU)等都广泛使用裁剪算法。
计算机图形学的裁剪算法可以分为两种:2D裁剪算法和
3D裁剪算法。
2D裁剪算法是基于二维空间的,它将一个几何
图形投影到一个平面上,然后按照指定的窗口裁剪;而3D裁
剪算法是基于三维空间的,它将一个几何图形投影到一个三维空间,然后按照指定的窗口裁剪。
2D裁剪算法的基本步骤如下:首先,将要裁剪的几何图
形投影到平面上;其次,计算出投影后的几何图形以及裁剪窗口之间的交点;最后,将裁剪窗口内的部分保留,而把窗口外的部分抛弃掉。
3D裁剪算法的基本步骤如下:首先,将要裁剪的几何图
形投影到三维空间;其次,计算出投影后的几何图形以及裁剪窗口之间的交点;最后,将裁剪窗口内的部分保留,而把窗口外的部分抛弃掉。
计算机图形学的裁剪算法在图形处理中有着重要的作用,它不仅能够有效减少图形处理时间,而且还可以节约存储空间。
此外,它还可以有效提高图形处理效率,提高图形显示效果。
但是,它也存在着一定的局限性,比如,当几何图形的运动变得复杂时,它就会变得费时费力,这就对性能产生了一定的影响。
总之,计算机图形学的裁剪算法是图形学的重要算法,它的应用非常广泛,在图形处理中有着重要的作用。
虽然它也存在着一定的局限性,但是它仍然是一种有效的图形处理算法。
第3章 5 裁剪算法

80年代初提出了著名的Liang-Barsky裁剪算法,通过 线段的参数化表示,实现快速裁剪,至今仍是计算机 图形学中最经典的算法之一; 80年代末到90年代,梁友栋先生致力于几何连续性的 研究,提出了一系列几何连续性方面的理论和方法,成 为国际上几何连续性研究的重要力量; 1991年梁友栋先生为首完成的成果“计算机图形生 成与几何造型研究”获国家自然科学三等奖;
1
2 3 4
Sutherlan-Cohen算法 中点分割算法
梁友栋-Barsky算法 Sutherlan-Hodgman逐边裁剪算法
裁剪的意义 为了描述图形对象,我们必须存储它的全部信息,但有时 为了达到分区描述或重点描述某一部分的目的,往往将 要描述的部分置于一个窗口内,而将窗口外的部分“剪 掉”,这个处理过程叫做裁剪,裁剪在计算机图形处理中 具有十分重要的意义。 裁剪就是将指定窗口作为图形边界,将窗口内的图形保 留,而窗口外的图形则被舍弃。 裁剪处理过程 1、图元在窗口内外的判别; 2、图形元素与窗口的求交。
直线裁剪方法
◦ Cohen-Sutherland裁剪算法
◦ 中点分割算法 ◦ 梁友栋-barskey算法
多边形裁剪方法
◦ Sutherland-Hodgman逐次多边形裁剪算法 ◦ Weiler-Atherton多边形裁剪算法
对于线段的裁剪基本思想是:
线段是否全不在窗口里,若是,转5 ② 线段是否全在窗口,若是,转4 ③ 计算该线段与窗口边界的交点,以此将线段分为两部分, 丢弃不可见的部分,剩余线段转2 ④ 保留并显示该线段 ⑤ 算法结束 可以看到算法的核心有两个:
区域编码
◦ 由窗口四条边所在直线把二维平面分成9个区域,每个区域 赋予一个四位编码,CtCbCrCl,上下右左
计算机图形学实验报告实验2裁剪算法实验

一、实验目的:直线段的裁剪:编码裁剪算法,中点分割裁剪算法。
二、实验内容://BasicGraph.cpp//请将下列裁剪程序补充完整,并用注释说明是何种裁剪算法void Encode (int x,int y,int *code,int XL,int XR,int YB,int YT) {//请将此程序补充完整int c=0;if(x<XL) c=c|LEFT;else if(x>XR) c=c|RIGHT;if(y<YB) c=c|BOTTOM;else if(y>YT) c=c|TOP;(*code)=c;}//编码裁剪算法:void C_S_Line(POINT &p1,POINT &p2,int XL,int XR,int YB,int YT) {//请将此程序补充完整int x1,x2,y1,y2,x,y,code1,code2,code;x1=p1.x; x2=p2.x; y1=p1.y; y2=p2.y;Encode(x1,y1,&code1,XL,XR,YB,YT);Encode(x2,y2,&code2,XL,XR,YB,YT);while(code1!=0||code2!=0){if((code1&code2)!=0) return;code=code1;if(code1==0) code=code2;if((LEFT&code)!=0){x=XL;y=y1+(y2-y1)*(XL-x1)/(x2-x1);}else if((RIGHT&code)!=0){x=XR;y=y1+(y2-y1)*(XR-x1)/(x2-x1);}if((BOTTOM&code)!=0){y=YB;x=x1+(x2-x1)*(YB-y1)/(y2-y1);}else if((TOP&code)!=0){y=YT;x=x1+(x2-x1)*(YT-y1)/(y2-y1);}if(code==code1){x1=x;y1=y;Encode(x,y,&code1,XL,XR,YB,YT);}else{x2=x;y2=y;Encode(x,y,&code2,XL,XR,YB,YT);}}p1.x=x1;p1.y=y1;p2.x=x2;p2.y=y2;}int IsInArea(POINT point,int XL,int XR,int YB,int YT){//请将此程序补充完整if(point.x>=XL && point.x<=XR && point.y>YB && point.y<YT) return 1;else return 0;}int NotIntersect(POINT begin,POINT end,int XL,int XR,int YB,int YT) {//请将此程序补充完整int maxx,maxy,minx,miny;maxx=(begin.x>end.x)?begin.x:end.x;minx=(begin.x<end.x)?begin.x:end.x;maxy=(begin.y>end.y)?begin.y:end.y;miny=(begin.y<end.y)?begin.y:end.y;if(maxx<XL|| minx>XR||maxy<YB||miny>YT) return 1;else return 0;}//中点裁剪算法:POINT ClipMid(POINT begin,POINT end,int XL,int XR,int YB,int YT){//请将此程序补充完整POINT mid,temp;if(IsInArea(begin,XL,XR,YB,YT)) temp=begin;else if(NotIntersect(begin,end,XL,XR,YB,YT)) temp=begin;else{mid.x=(begin.x+end.x)/2;mid.y=(begin.y+end.y)/2;if(abs(mid.x-end.x)<=1&& abs(mid.y-end.y)<=1) temp=mid;else{if(NotIntersect(begin,mid,XL,XR,YB,YT))temp=ClipMid(mid,end,XL,XR,YB,YT);elsetemp=ClipMid(begin,mid,XL,XR,YB,YT);}}return temp;}//Liang-Barsky直线裁剪算法:void ClipParameter(POINT &p1,POINT &p2,int XL,int XR,int YB,int YT) {float u1=0.0,u2=1.0;float dx=p2.x-p1.x,dy=p2.y-p1.y;if(clipTest(-dx,p1.x-XL,&u1,&u2))if(clipTest(dx,XR-p1.x,&u1,&u2))if(clipTest(-dy,p1.y-YB,&u1,&u2))if(clipTest(dy,YT-p1.y,&u1,&u2)){if(u2<1.0){p2.x=p1.x+u2*dx;p2.y=p1.y+u2*dy;}if(u1>0.0){p1.x=p1.x+u1*dx;p1.y=p1.y+u1*dy;}}}int clipTest(float p,float q,float *u1,float *u2){float r;int remainFlag=1;if(p<0.0){r=q/p;if(r>*u2) remainFlag=0;else if(r>*u1) *u1=r;}else if(p>0.0){r=q/p;if(r<*u1) remainFlag=0;else if(r<*u2) *u2=r;}else //*p=0if(q<0.0) remainFlag=0;return remainFlag;}//逐边裁剪算法://typedef struct tRes { int yes,isIn; POINT pout;} Res;Res TestIntersect(int edge,int type,POINT p1,POINT p2){//判断p2是否在所裁剪的窗边edge的内侧,是否与p1点分别在窗边edge的异侧float dx,dy,m;Res res;int isIn=0,yes=0;POINT pout;dy=p2.y-p1.y;dx=p2.x-p1.x;m=dy/dx;switch(type){case 1: /*right*/if(p2.x<=edge){isIn=1;if(p1.x>edge)yes=1;}else if(p1.x<=edge)yes=1;break;case 2: /*bottom*/if(p2.y>=edge){isIn=1;if(p1.y<edge)yes=1;}else if(p1.y>=edge)yes=1;break;case 3: /*left*/if(p2.x>=edge){isIn=1;if(p1.x<edge)yes=1;}else if(p1.x>=edge)yes=1;break;case 4: /*top*/if(p2.y<=edge){isIn=1;if(p1.y>edge)yes=1;}else if(p1.y<=edge)yes=1;default: break;}if(yes){if((type==1) || (type==3)){ pout.x=edge;pout.y=p1.y+m*(pout.x-p1.x);}if((type==2) || (type==4)){ pout.y=edge;pout.x=p1.x+(pout.y-p1.y)/m;}}res.isIn=isIn;res.yes=yes;res.pout=pout;return res;}int clipSingleEdge(int edge,int type,int nin,POINT pin[50],POINT pout[50])/*对多边形pin与窗边edge进行裁剪,返回裁剪后的多边形pout及点数*/ {int i,k=0;POINT p;Res res;p.x=pin[nin-1].x;p.y=pin[nin-1].y;for(i=0;i<nin;i++){res=TestIntersect(edge,type,p,pin[i]);if(res.yes){ pout[k].x=res.pout.x;pout[k].y=res.pout.y;k++;} if(res.isIn){ pout[k].x=pin[i].x;pout[k].y=pin[i].y;k++;}p.x=pin[i].x;p.y=pin[i].y;}return k;}void ClipEdgePolygon(POINT ps[50],int &n,int XL,int XR,int YB,int YT) { /*对多边形ps进行逐边裁剪*/int n1=0,n2=0;POINT pt[50];n1=clipSingleEdge(XR,1,n,ps,pt);n2=clipSingleEdge(YB,2,n1,pt,ps);n1=clipSingleEdge(XL,3,n2,ps,pt);n2=clipSingleEdge(YT,4,n1,pt,ps);n=n2;}//多边形编码裁剪算法:void ClipEncodePolygon(POINT ps[50],int &n,int XL,int XR,int YB,int YT) {POINT tp[50];int k=0,m;int code1,code2,code;int x,y;for(int i=0;i<n-1;i++){Encode(ps[i].x,ps[i].y,&code1,XL,XR,YB,YT);Encode(ps[i+1].x,ps[i+1].y,&code2,XL,XR,YB,YT);code=code1;m=i;for(int j=0;j<2;j++){if((code1 & code2)!=0) //线段两端都在窗口外的同一侧{switch(code){case 1:x=XL;y=ps[m].y;break;case 2:x=XR;y=ps[m].y;break;case 4:x=ps[m].x;y=YB;break;case 5:x=XL;y=YB;break;case 6:x=XR;y=YB;break;case 8:x=ps[m].x;y=YT;break;case 9:x=XL;y=YT;break;case 10:x=XR;y=YT;break;}tp[k].x=x;tp[k].y=y;k++;}else if((code1 & code2)==0) //线段两端不在窗口的同一侧{if(code==0){tp[k]=ps[m];k++;}else if ((LEFT & code) !=0) //线段与左边界相交 {x=XL;y=ps[i].y+(ps[i+1].y-ps[i].y)*(XL-ps[i].x)/(ps[i+1].x-ps[i].x);if(y>YB && y<YT){tp[k].x=x;tp[k].y=y;k++;}}else if((TOP & code)!=0) //线段与上边界相交{y=YT;x=ps[i].x+(ps[i+1].x-ps[i].x)*(YT-ps[i].y)/(ps[i+1].y-ps[i].y);if(x>XL && x<XR){tp[k].x=x;tp[k].y=y;k++;}}else if((RIGHT & code)!=0) //线段与右边界相交 {x=XR;y=ps[i].y+(ps[i+1].y-ps[i].y)*(XR-ps[i].x)/(ps[i+1].x-ps[i].x);if(y>YB && y<YT){tp[k].x=x;tp[k].y=y;k++;}}else if((BOTTOM & code) != 0) //线段与下边界相交 {y=YB;x=ps[i].x+(ps[i+1].x-ps[i].x)*(YB-ps[i].y)/(ps[i+1].y-ps[i].y);if(x>XL && x<XR){tp[k].x=x;tp[k].y=y;k++;}}}code=code2;m++;}//for(j)}//for(i)for(i=0;i<k;i++)ps[i]=tp[i];n=k;}//函数的调用,裁剪窗口的调整//DrawView.cpp文件//裁剪窗口的调整CDrawView::CDrawView(){/************请在此函数中将裁剪窗口大小调整为长度100单位像素,宽度50单位像素的矩形********/// TODO: add construction code here//m_pWidth=1;m_pStyle=PEN_STYLE_SOLID;m_pColor=RGB(0,0,0);m_FFlag=0;m_FColor=RGB(0,0,0);m_HFlag=0;CurrentDraw=DRAW_VCLINE;m_Num=0;m_Drag=0;m_HCursor=AfxGetApp()->LoadStandardCursor(IDC_CROSS);//DrawType=0;ClipFlag=0;ClipType=-1;XL=200;XR=300;YB=150;YT=200;//XL=200;XR=500;YB=150;YT=400;ClipWindowColor=RGB(192,192,50);}void CDrawView::OnDraw(CDC* pDC){CDrawDoc* pDoc = GetDocument();ASSERT_VALID(pDoc);// TODO: add draw code for native data hereif(ClipFlag){CPen NewPen,*pOldPen;NewPen.CreatePen(PS_DASH,1,ClipWindowColor);pOldPen=pDC->SelectObject(&NewPen);pDC->MoveTo(XL,YB);pDC->LineTo(XR,YB);pDC->LineTo(XR,YT);pDC->LineTo(XL,YT);pDC->LineTo(XL,YB);}int index;index=pDoc->GetShapeNumber();for(int i=0;i<index;i++)pDoc->GetShape(i)->Drawing(pDC);}void CDrawView::OnInitialUpdate(){CSize sizeTotal;sizeTotal.cx = 640; sizeTotal.cy = 480;SetScrollSizes(MM_TEXT, sizeTotal);// TODO: Add your specialized code here and/or call the base class }void CDrawView::OnLButtonDown(UINT nFlags, CPoint point){// TODO: Add your message handler code here and/or call defaultCClientDC dc(this);OnPrepareDC(&dc);dc.DPtoLP(&point);m_pPrev=point;m_pOrigin=point; //点击鼠标左键作为拖动绘图的第一点m_Drag=1;SetCapture();RECT rect;GetClientRect(&rect);ClientToScreen(&rect);ClipCursor(&rect);CScrollView::OnLButtonDown(nFlags, point);}//函数调用处void CDrawView::OnLButtonUp(UINT nFlags, CPoint point){// TODO: Add your message handler code here and/or call defaultif(m_Drag){m_Drag=0;ReleaseCapture();ClipCursor(NULL);CDrawDoc *pDoc=GetDocument();CShape *pShape;POINT p1,p2;if(CurrentDraw==DRAW_VCLINE || CurrentDraw==DRAW_DDALINE ||CurrentDraw==DRAW_MIDLINE || CurrentDraw==DRAW_BSHLINE){if(ClipFlag){switch(ClipType){/****************编码裁剪函数调用处*************/case CLIP_ENCODE:C_S_Line(m_pOrigin,m_pPrev,XL,XR,YB,YT); break; /****************中点分割裁剪函数调用处************/case CLIP_MIDPOINT: ClipMid(m_pPrev,m_pOrigin,XL,XR,YB,YT);p1=ClipMid(m_pPrev,m_pOrigin,XL,XR,YB,YT);p2=ClipMid(m_pOrigin,m_pPrev,XL,XR,YB,YT);m_pOrigin=p1;m_pPrev=p2;break;case CLIP_PARAMETER:ClipParameter(m_pOrigin,m_pPrev,XL,XR,YB,YT);break;}}pShape=newCLine(m_pOrigin,m_pPrev,m_pWidth,m_pStyle,m_pColor,DrawType);pDoc->AddShape(pShape);}if(CurrentDraw==DRAW_RECTANGLE){if(ClipType==CLIP_WINDOW){XL=m_pOrigin.x;XR=m_pPrev.x;YB=m_pOrigin.y;YT=m_pPrev.y;}else{pShape=newCRectangle(m_pOrigin,m_pPrev,m_pWidth,m_pStyle,m_pColor,m_FFlag,m_FColor,m_HFlag,m_Hatch);pDoc->AddShape(pShape);}}if( CurrentDraw==DRAW_VCCIRCLE || CurrentDraw==DRAW_MIDCIRCLE || CurrentDraw==DRAW_BSHCIRCLE){pShape=newCCircle(m_pOrigin,m_pPrev,m_pWidth,m_pStyle,m_pColor,m_FFlag,m_FColor,m_HFlag,m_Hatch,DrawType);pDoc->AddShape(pShape);}if(CurrentDraw==DRAW_VCELLIPSE || CurrentDraw==DRAW_MIDELLIPSE) {pShape=newCEllipse(m_pOrigin,m_pPrev,m_pWidth,m_pStyle,m_pColor,m_FFlag,m_FColor,m_HFlag,m_Hatch,DrawType);pDoc->AddShape(pShape);}pDoc->UpdateAllViews(NULL);}CScrollView::OnLButtonUp(nFlags, point);}三实验结果:四、实验总结通过这次试验使我了解到如何运用计算机程序对窗口进行剪裁,了解到编码剪裁算法直观方便,速度较快,中点分割剪裁算法不用进行乘除运算,剪裁效率高,Liang-Barsky直线裁剪算法更快。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
Cohen-Sutherland算法的实现
算法优点: 1) 简单将不需剪裁的直线舍弃。法则是:若是 一条直线的两头在同一地区,则该直线不需裁剪, 不然,该直线为也许剪裁直线。 2) 对也许剪裁的直线缩小了与之求交的边框范 畴。法则是:若是直线的一个端点在上(下、左、 右)域,则此直线与上边框求交,然后移除上边框 以上的局部。该法则对直线的另一端点也实用。这 样,一条直线至多只需与两条边框求交。
p1
梁友栋-Barsky算法(Liang-Barsky算法)
写入图形学教科书的唯一中国人的算 法
梁有栋教授的贡献
Liang-Barsky算法 几何连续理论 从几何学与纤维缠绕理论到基因工程
梁友栋先生出生于1935年7月,福建福州人. 19561960年于复旦大学作为研究生,师从苏步青先生学 习几何理论,1960年研究生毕业后任教于浙江大学 数学系,1984-1990年任数学系主任. 一直致力于计算机辅助几何设计与计算机图形学方 面的研究.
如:CD C:0000 D:1000 Dt位等于1,与上边相交
9
8
10
1 2
AB:与窗口左边相交,求交点H, AH和BH显然不可见
EF:与窗口上边交点I,弃EI 与窗口下边交点K,弃KF 与窗口右边交点J,弃KJ 保留IJ
0
5 4 6
E
I
L M J
C
CD:与窗口上边交点L, 弃CL DO G 与窗口下边交点O, 弃OD H 与窗口右边交点M, 弃LM 与窗口左边交点N, 弃NO 保留MN 求交次数最多 可达?
区域编码
◦ 由窗口四条边所在直线把二维平面分成9个区域,每个区域 赋予一个四位编码,CtCbCrCl,上下右左
1 Ct 0 1 Cb 0 1 Cr 0 1 Cl 0 当y y max else 当y y min else 当x x max else 当x x min else
第三章
裁剪是抽取数据的一部分,或者识别一个指定区域 内部或外部的画面或图片的成分的过程。 按照被裁剪的图形可以分为:二维裁剪算法和三维 裁剪算法 按照裁剪区域的形状可以分为:规则区域和不规则 区域 本章二维裁剪算法只考虑举行和多边形,裁剪的对 象只考虑点、线段、多边形、字符等。 三维裁剪算法只考虑正方体和平截头正四棱锥两种 裁剪盒,裁剪对象只考虑线段。
y yT yB P0 xL A B
P1 C
D
设要裁剪的线段是P0P1。 P0P1和
窗口边界交于A,B,C,D四点.
算法的基本思想是从A,B和P0三点
xB
x
中找出最靠近P1的点,图中要找 的点是B; 从C,D和P1中找出最靠 近P0的点。图中要找的点是C点。 那么BC就是P0P1线段上的可见部 分。
1.
p2 示例: b
对于A线段, p1 c P1(0000).p2(0000),显示出来 对于线段B p1 p2 a p1 P1(1000),p2(1000),非零,丢弃 p2 d 对于线段C P1(1000),p2(0010),“与”为零,不可判断,对P1 开始,pa=p1,pb=p2,求得pm1(1000),pb与pm1 p2 “与”为零,不可判断, pa=pm1,求得pm2(0010),Pb与pm2 “与”非零,丢弃,pb变为pm2。直到线段长度小于一个给定的值,计算该 点的 编码。对于p2做同样处理,最后得出线段C不可见。 对于选段D: P1(0001),p2(0100), “与”为零,不可判断。 从p1开始, pa=p1,pb=p2,pm1(0100),pbpm1丢弃,pm1代替pb,求 得pm2(0000),pm2与pb“与”,为零 ,pa=pm2,计算pm3(0000), 最后汇聚成一点。计算该点的编码。对 p2进行同样处理。
区域编码
◦ 由窗口四条边所在直线把二维平面分成9个区域,每个区域 赋予一个四位编码,CtCbCrCl,上下右左
1 当y y max 例如:给出直线 AB的 Ct 端点编码。 0 else 1 当y y min A : 1001 问题 2: 给定一条线段 Cb 0 else B:0010 端点的区域码 ,该如何 1 当x x max 判断属于哪一种情况? Cr 0 1 Cl 0 else 当x x min else
A
N B
K F
最多求交次数 为4
Cohen—SutherLand直线裁剪算法 直线和窗口的关系可以分为如下三类: 1) 整条直线在窗口之内,此时,不需裁剪,显 现整条直线。 2) 整条直线在窗口之外,此时,不需裁剪,不 显现整条直线。 3) 局部直线在窗口之内,局部直线在窗口之外, 此时,必要求出直线与窗框之交点,并将窗口外的 直线局部剪裁掉,显现窗口内的局部。
直线裁剪方法
◦ Cohen-Sutherland裁剪算法
◦ 中点分割算法 ◦ 梁友栋-barskey算法
多边形裁剪方法
◦ Sutherland-Hodgman逐次多边形裁剪算法 ◦ Weiler-Atherton多边形裁剪算法
对于线段的裁剪基本思想是:
线段是否全不在窗口里,若是,转5 ② 线段是否全在窗口,若是,转4 ③ 计算该线段与窗口边界的交点,以此将线段分为两部分, 丢弃不可见的部分,剩余线段转2 ④ 保留并显示该线段 ⑤ 算法结束 可以看到算法的核心有两个:
中心思想: 1. 算法采用了裁剪算法相同的编码和检查办法 2. 先判断完全可见和显然不可见的线段 3. 不能判断的线段使用中点分割方法。 4. 算法从线段的两个端点出发,分别求出离本端点 最远的可见点。 5. 两个可见点之间的线段就是所求的可见线段
算法的思想:
①
② ③ ④
对线段p1p2进行编码,如果p1p2可以直接保留或者舍弃, 算法结束,否则进入步骤2. 从p1出发寻找离p1最远的可见点B; 从p2出发寻找离p2最远的可见点A; AB之间的线段就是可见的线段。
A
B
若P1P2完全在窗口内, code1=0,且code2=0,则“取”; 若P1P2明显在窗口外code1&code2≠0,则“弃” ; 若P1P2明显在窗口外code1&code2=0,则不确定 ; 在交点处把线段分为两段。其中一段完全在窗口外,可弃
之。然后对另一段重复上述处理。
算法的缺点?
第一,当code1&code2≠0,能判断出直 线段显然在窗口之外. 然而这个条件 并没有包含所有的在窗口之外的直线段,比如: code1=0001,code2=0100,如图 中直线段P3P4,此时, code1&code2=0,但直线段也是完全在窗口之外.这就意 味着这类型的直线必须计算与窗口的边线求交,其实是 一种无意义操作;
线段p1(x1,y1)p2(x2,y2)的参数方程可以表示为:
x=x1+t△x y=y1+t△y 其中,△x=x2-x1, △y=y2-y1,0<=t<=1
P(x,y)代表了该线段上的一个点,其值由参数t确定。 由公式可知,当t=0时,该点为P1(x1,y1),当t=1时,该 点为P2(x2,y2)。 如果点P(x,y)位于由坐标(xmin,ymin)和(xmax,ymax) 所确定的矩形窗口内,那么下式成立:
如何改进?
第二,对于p1p2,当code 码中同时有两位不等于0 时,求 交运算的次数多达四次. 如图2中直线P1P2 与窗口的四 条边的交点分别为A、B、C、D。而实际上,直线只与窗 口的两条边相交,另外两个交点发生在边的延长线上。
பைடு நூலகம்
编码裁剪算法需要求线段与边界的交点 如果使用将线段一分为二,那么求交点的过程可以 用二分法逐渐逼近来得到 中点分割算法比较适合硬件实现 如果用软件来实现,中点分割算法的运行速度反而 比编码裁剪算法慢
九十年代后期,六十多岁的他积极开展纤维缠绕几何 设计的研究,为我国几何设计与计算机图形学的研究 做出了杰出贡献. 鉴于梁友栋先生的杰出贡献和成就,中国工业与应用 数学学会几何设计与计算专委会授予梁友栋先生第 二届“中国几何设计与计算贡献奖”。 梁友栋和Barsky独立地提出了更快的参数化线段裁 剪算法
①
① 判断线段与窗口的关系。 ② 获取线段与窗口的交点。
基本思想: 对于每条线段算法分成三 种情况处理:
第一种:若线段完全在窗口内,则显示该线段称 “完全可见”,需要“取”之; 第二种:若明显在窗口外,称“完全不可见”.则 丢弃该线段。 第三种:若线段不满足“取”或 “弃”的条件,则 在交点处把线段分为两段. 其中一段完全在窗口 外,可弃之. 然后对另一段重复上述处理。 问题1:用什么方法解决该问题?
80年代初提出了著名的Liang-Barsky裁剪算法,通过 线段的参数化表示,实现快速裁剪,至今仍是计算机 图形学中最经典的算法之一; 80年代末到90年代,梁友栋先生致力于几何连续性的 研究,提出了一系列几何连续性方面的理论和方法,成 为国际上几何连续性研究的重要力量; 1991年梁友栋先生为首完成的成果“计算机图形生 成与几何造型研究”获国家自然科学三等奖;
xmin≤x0+ t· Δx≤xmax ymin≤y0+ t· Δy≤ymax
t· Dk ≤Mk , k=1,2,3,4 其中: ① D1=-Δx, M1=x0-xmin ② D2= Δx, M2=xmax-x0 ③ D3=-Δy, M3=y0-ymin ④ D4= Δy, M4=ymax-y0 进行公式验证: 公式1表示为: t*- Δx<=x0-xmin-xmin<=x0+t Δx 公式2表示为:t Δx<=xmax-x0- x0+t Δx<=xmax 公式3表示为:t*- Δy<=y0-ymin ymin<=y0+t Δy 公式4表示为: t Δy<=ymax-y0 y0+t Δy<=ymax