计算机图形学代码

DDA生成直线算法：void DDALine(int x₀,int y₀,int x₁,int y₁,int color)

{int i;float dx,dy,length,x,y;

if(fabs(x₁ -x₀)>=fabs(y₁ -y₀))

length=fabs(x₁-x₀);

else

length=fabs(y₁-y₀);

dx = (x₁-x₀)/length;    dy=(y₁-y₀)/length;

i = 1 ;x = x ₀;y= Yo;

while(i<=length)

{SetPixel ( int(x+0.5),int(y+0.5),color);

x=x+ dx;y=y+dy;i + + ;}〕

Bresenham画线算法：(0<k<1)

void Bresenham_Line( int x₀, int y₀, int x₁, int y₁, int color)

{int dx,dy,e,i,x,y;

dx = x₁ -x₀,dy = y₁ - y₀,e=2*dy-dx;

x = x₀,y= y₀;

for ( i= 0;i < =dx;i++)

{SetPixel(x,y,color);     x++ ;

if ( e>= 0)

{y++ ;e=e+2*dy-2*dx;}else

e=e+2*dy;}}

中点画圆法：MidPointCircle(int r,int color)

{int x,y;int e;x = 0;y = r;e = 1 - r;

CirclePoints(x,y,color);

while( x < = y)

{if( e< 0)

e+ =2 * x + 3;

else

{e+= 2 * ( x- y) + 5; y -- ;

}  x++ ;

CirclePoints(x,y,color);}}

四连通：（内点表示）

void FloodFill4(int x,int y,int oldcolor,int newcolor)

{if(GetPixel(x,y)==oldcolor)

{SetPixel(x,y,newcolor);

FloodFill4(x,y+1,oldcolor,newcolor);

FloodFill4(x,y-1,oldcolor,newcolor);

FloodFill4(x-1,y,oldcolor,newcolor);

FloodFill4(x+1,y,oldcolor,newcolor);}}

（边界表示）

void BoundaryFill4(int x,int y,int boundarycolor,int newcolor)

{int color=GetPixel(x,y);

if(color!=newcolor &&color!=boundarycolor)

{SetPixel(x,y,newcolor);

BoundaryFill4(x,y+1,boundarycolor,newcolor);

BoundaryFill4(x,y-1,boundarycolor,newcolor);

 BoundaryFill4(x-1,y,boundarycolor,newcolor);

 BoundaryFill4(x+1,y,boundarycolor,newcolor);}}

边表：void polyfill(polygon,color)

 int color;

多边形polygon；

{for(各条扫描线，标识为i)

{初始化新边表头指针 NET[i]；

把ymin=1的边放进边表NET[i]；}

y= 最低扫描线号；

初始化活性边表AET为空；

for(各条扫描线i) {把新边表NET[i]中的边结点用插入排序法插入AET表，使之按x坐标递增顺序排列；遍历AET表，把配对交点区间(左闭右开)上的像素(x，y)，用drawpixel(x，y，color)改写像素颜色值；遍历AET表，把ymax=1的结点从AET表中删除，并把ymax＞i结点的x值递增△x；若允许多边形的边自相交，则用冒泡排序法对AET表重新排序；}}

边界标志算法：void edgemark_fill(polydef，color)

多边形定义 polydef; int color;

{对多边形polydef每条边进行直线扫描转换；inside = FALSE;

for(每条与多边形polydef相交的扫描线y)

for(扫描线上每个像素x)

{if(像素x被打上边标志)

inside =! (inside);

if(inside!= FALSE)

drawpixel(x, y，color);

else  drawpixel(x, y,background);}}​

扫描线种子填充算法：typedef struct{

int x;  int y;

}Seed;

void ScanLineFil14(int x,int y,COLORREF oldcolor,COLORREF newcolor)

{int x1,xr,i;

bool spanNeedFill;

Seed pt;

setstackempty();

pt.x=x; pt.y=y;

stackpush(pt);

while(! isstackempty())

{pt = stackpop()；y=pt.y；x= pt.x；

while(getpixel(x,y)==oldcolor)

{SetPixel(x,y,newcolor);x++;}

xr=x-1；   x=pt.x-1；

while(getpixel(x,y)==oldcolor)

{ SetPixel(x,y,newcolor);x--；}

x1=x+1；x=x1；y=y+1;

while(x＜xr)

{spanNeedFill = FALSE;

while(getpixel(x,y)==oldcolor)

{spanNeedF i11 = TRUE;x++;}

if(spanNeedFill)

{ pt.x=x-1;pt.y=y;stackpush(pt);

spanNeedFi11 = FALSE;

}while(getpixel(x,y)！=oldcolor && x＜xr) x++；

}x=x1;y=y-2;while(x＜xr){spanNeedFi11 = FALSE;

while(getpixel(x，y)==oldcolor)

{spanNeedFill = TRUE;

x++；}if(spanNeedFil1)

{pt.x=x-1;pt.y=y;

stackpush(pt);

spanNeedFil1=FALSE;}while(getpixel(x，y)!=oldcolor && x<xr) x++;}}}

Z缓冲器消隐算法：

帧缓冲区置成背景色；

Z缓冲区置成某初始值，该值比场景在观察坐标系下的最小z值还小；

for(各个多边形)

{扫描转换该多边形；

for(多边形所覆盖的每个像素(x，y))

{计算该像素所对应多边形上的点在观察坐标系下的z坐标值Z(x，y)；

if(Z(x,y)大于Z缓冲区在(x,y)处的值)

{Z缓冲区中(x,y)处深度值替换为Z(x,y);

帧缓冲区中(x，y)处亮度值替换为多边形在(x，y)处的亮度值；}}}

光线跟踪：

DrawImage()

{for(帧缓冲区中的每一个像素e)

{计算像素e所对应的投影平面上的点 e′,从视点V向 e′发出一条光线,方向记为direction。

Ray-Tracing(e',direction,1.0,1,&I);

将像素e的颜色设置为 (Ir,Ig,Ib);}}

Ray-Tracing(startPoint,direction,weight,depth,I)

if(( depth>TREE_MAX_DEPTH)||( weight<CONTRIBUTE_COEFFICIENT_THRESHOLD))

{*I= 0;  return;}

将起点为startPoint、方向为direction的光线与场景中的所有物体表面求交；如有交点，记距离startPoint最近的交点P₁，记光线在相交物体上P₁处的反射方向和透射方向分别为direction_r、direction_t。相交物体表面的镜面反射系数和透射系数分别为Ks，Kt。

if(无交点)

{*I=0；return;}

else

{根据简单光照模型计算P₁处因光源直接照射引起的(-1*direction)方向的反射光强度Ic；

Ray-Tracing(P₁, direction_r, weight*Ks，depth+ 1,&Is);

Ray-Tracing(P₁, direction_t, weight*Kt,depth+1,&It);

*I=Ic+Ks·Is+Kt·It

return;}}

曲线的De Casteljau算法：

float Casteljau( int k, float controlsa[], float t)

{int r,i;  float t1, contolsa[10];    t1=1.0-t;

for(i=0；i<=k,；i++)

controlsa[]=controls[i];

for(r=1:r<=k ;i++)     for(i=0; i<= k-r;i++)

controlsa[i]=t1*controlsa[i]+t*controlsa[i];

return controlsa[0];)

void bezier_to_points(int k，int npoints,float controls[], float points[])

{float t, delt;  int i;

delt=1.0/(float)npoints;

t=0.0;

for(i=0；i<=points；i++)

points=Casteljau( k, controls,t)

 t=t + delt;}

曲面的De Casteljau算法：

DeCasteljauSurf(P, n, m, u, v, S)

{  if(n<=m)

{ for(j=0;j<=m;j++)

DeCasteljau(P[j][], n, u, Q[j])

DeCasteljau(Q, m, v, S) ; }

else {for(i=0;i<=n;i++)

DeCasteljau(P[][i], m, v, Q[i]) ;

DeCasteljau(Q, n, u, S) ; }

光线跟踪算法：

设置视点，投影平面以及窗口的参数；

For （窗口内的每一条扫描线）

    for （扫描线上的每一个像素）

    {    确定从视点指向像素中心的光线ray；

          像素的颜色=RayTracing(ray,1)；}

Color RayTracing(Ray ray, int depth)

{    求ray与物体表面最近的交点P；

      if （有交点）

      {    用局部光照明模型计算P点的Ic；

            color = Ic;

            if  （depth<给定的最大跟踪层次）

            {    计算ray的反射光线；

                   Is=RayTracing(反射光线，depth+1）；

                  if  （物体是透明的）

                  {    计算ray的透射光线；

                        It=RayTracing(透射光线，depth+1);  }

                 color = Ic + Is + It ;    }    else color = black；

       return color ;}