3d.c

计算机图形学

         if(dp[i].d>dp[j].d)
         {
           temp.s = dp[i].s; temp.d = dp[i].d;
           dp[i].s = dp[j].s; dp[i].d = dp[j].d;
           dp[j].s = temp.s; dp[j].d = temp.d;
         }
      }
   }

   i = 0;

   while(i <= 4)      /*画家算法*/
  {
    switch(dp[i++].s)
    {


  case'A': B.x = A.x - L/2.0; B.y = A.y; B.z = A.z;  /*四条竖线*/
           Cft3D(B,15,15,h + 40,BB,qx,qy,qz);
           break;

  case'B': B.x = A.x - L/2.0 + L/4.0; B.y = A.y; B.z = A.z;
           Cft3D(B,15,15,h - 15,BB,qx,qy,qz);
           break;

  case'C': B.x = A.x - L/2.0 + L/2.0; B.y = A.y; B.z = A.z;
           Cft3D(B,15,15,h - 15,BB,qx,qy,qz);
           break;

  case'D': B.x = A.x - L/2.0 + 3*L/4.0; B.y = A.y; B.z = A.z;
           Cft3D(B,15,15,h - 15,BB,qx,qy,qz);
           break;

  case'E': B.x = A.x - L/2.0 + L; B.y = A.y; B.z = A.z;
           Cft3D(B,15,15,h + 40,BB,qx,qy,qz);
           break;


   }
  }

  B.x = A.x; B.y = A.y; B.z = A.z + h/2.0;        /*两条横线*/
           Cft3D(B,15,L - 15,15,BB,qx,qy,qz);


  B.x = A.x; B.y = A.y; B.z = A.z - h/2.0;
           Cft3D(B,15,L - 15,15,BB,qx,qy,qz);

}

/*****************************/
/**    栏杆(平行于Y轴)      **/
/*****************************/
void LanGan_Y(Point3D A,int L,int h,Point3D BB,float qx,float qy,float qz)
{
  int X,Y,Z,i,j;
  Point3D B;
  Depth dp[7],temp;

  X = A.x - L/2.0; Y = A.y; Z = A.z;    /*五条竖线的深度*/
     Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
  dp[0].s = 'A'; dp[0].d = Z;

  X = A.x - L/2.0+ L/4.0; Y = A.y; Z = A.z;
     Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
  dp[1].s = 'B'; dp[1].d = Z;

  X = A.x - L/2.0 + L/2.0; Y = A.y; Z = A.z;
    Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
  dp[2].s = 'C'; dp[2].d = Z;

  X = A.x - L/2.0 + 3*L/4.0 - 10; Y = A.y; Z = A.z;
     Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
  dp[3].s = 'D'; dp[3].d = Z;

  X = A.x - L/2.0 + L; Y = A.y; Z = A.z;
     Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
  dp[4].s = 'E'; dp[4].d = Z;


  for(i = 0; i <=3; i++)                 /*  深度排序  */
   {
      for(j = i+1;j<=4;j++)
      {
         if(dp[i].d>=dp[j].d)
         {
           temp.s = dp[i].s; temp.d = dp[i].d;
           dp[i].s = dp[j].s; dp[i].d = dp[j].d;
           dp[j].s = temp.s; dp[j].d = temp.d;
         }
      }
   }

   i = 0;

   while(i <= 4)      /*画家算法*/
  {
    switch(dp[i++].s)
    {


  case'A': B.x = A.x; B.y = A.y - L/2.0; B.z = A.z;  /*四条竖线*/
           Cft3D(B,15,15,h +40,BB,qx,qy,qz);
           break;

  case'B': B.x = A.x; B.y = A.y - L/2.0 + L/4.0; B.z = A.z;
           Cft3D(B,15,15,h - 15,BB,qx,qy,qz);
           break;

  case'C': B.x = A.x ; B.y = A.y- L/2.0 + L/2.0; B.z = A.z;
           Cft3D(B,15,15,h - 15,BB,qx,qy,qz);
           break;

  case'D': B.x = A.x ; B.y = A.y- L/2.0 + 3*L/4.0; B.z = A.z;
           Cft3D(B,15,15,h - 15,BB,qx,qy,qz);
           break;

  case'E': B.x = A.x ; B.y = A.y- L/2.0 + L; B.z = A.z;
           Cft3D(B,15,15,h + 40,BB,qx,qy,qz);
           break;


   }
  }
    B.x = A.x; B.y = A.y; B.z = A.z + h/2.0;
           Cft3D(B,L - 15,15,15,BB,qx,qy,qz);


    B.x = A.x; B.y = A.y; B.z = A.z - h/2.0;
           Cft3D(B,L - 15,15,15,BB,qx,qy,qz);


}





/*****************************/
/**          树             **/
/*****************************/
void Shu(Point3D A,int r,Point3D BB,float qx,float qy,float qz)
{
   int k;
   Point3D B;
   setbkcolor(2);
   B.x = A.x + r; B.y = A.y; B.z = A.z;
   k=0;
   while(k<360)
   {
     Rect_3DXY_XuanZhuan(B,40,20,7,k,A,3,BB,qx,qy,qz);
     k = k+30;
   }

   B.x = A.x; B.y = A.y; B.z = A.z + 50; /* 树干 */
   YuanZhu3D(B,15,100,BB,qx,qy,qz);

   B.x = A.x; B.y = A.y; B.z = A.z + 15;
   k = 15;
   while(B.z>=A.z)
   {
       if((int)B.z%2==0)
          Yuan_3DXY(B,k,2,BB,qx,qy,qz);
       else
          Yuan_3DXY(B,k,4,BB,qx,qy,qz);
       B.z--;
       k = k + 1;
   }



   B.x = A.x; B.y = A.y; B.z = A.z + 100; /*树叶*/
   k=80;
   while(k>=0)
   {
       Yuan_3DXY(B,k--,8,BB,qx,qy,qz);
   }

   k=100;            /*树叶1*/
   while(k>=0)
   {
      if(k%2==0)
        Yuan_3DXY(B,k,3,BB,qx,qy,qz);
      else
          Yuan_3DXY(B,k,4,BB,qx,qy,qz);

      B.z = B.z + 2;
      k--;

   }

   B.x = A.x; B.y = A.y; B.z = A.z + 200;
   k=80;
   while(k>=0)
   {
       Yuan_3DXY(B,k--,8,BB,qx,qy,qz);
   }

   k=80;            /*树叶2*/
   while(k>=0)
   {
      if(k%2==0)
        Yuan_3DXY(B,k,3,BB,qx,qy,qz);
      else
          Yuan_3DXY(B,k,4,BB,qx,qy,qz);

      B.z = B.z + 2;
      k--;

   }


}



/*****************************/
/**          凳子           **/
/*****************************/
void DengZi(Point3D A,int a,int b,int h,Point3D B,float qx,float qy,float qz)
{
   int X,Y,Z,i,j;
   Point3D D;
   Depth dp[12],temp;
   setbkcolor(2);
 /*  Rect_3DXY(A,a + 100,b + 100,8,B,qx,qy,qz);*//*凳子下的地板*/



                                              /* 深度初始化*/
   X =  A.x - b/2.0; Y = A.y - a/2.0; Z = A.z + h/2.0;
      Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
   dp[0].s = 'A'; dp[0].d = Z;


    X =  A.x - b/2.0; Y = A.y + a/2.0; Z = A.z + h/2.0;
      Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
   dp[1].s = 'B'; dp[1].d = Z;

   X =  A.x + b/2.0; Y = A.y - a/2.0; Z = A.z + h/2.0;
      Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
   dp[2].s = 'C'; dp[2].d = Z;


   X =  A.x + b/2.0; Y = A.y + a/2.0; Z = A.z + h/2.0;
      Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
   dp[3].s = 'D'; dp[3].d = Z;

   X =  A.x - b/2.0; Y = A.y; Z = A.z + h/2.0;
      Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
   dp[4].s = 'E'; dp[4].d = Z;

   X =  A.x + b/2.0; Y = A.y; Z = A.z + h/2.0;
      Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
   dp[5].s = 'F'; dp[5].d = Z;

   X = A.x; Y = A.y - a/2.0; Z = A.z + h;
      Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
   dp[6].s = 'H'; dp[6].d = Z;

   X = A.x; Y = A.y - a/2.0 + a/5.0; Z = A.z + h;
       Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
   dp[7].s = 'I'; dp[7].d = Z;

   X = A.x; Y = A.y - a/2.0 + 2*a/5.0; Z = A.z + h;
        Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
   dp[8].s = 'J'; dp[8].d = Z;

   X = A.x; Y = A.y - a/2.0 + 3*a/5.0; Z = A.z + h;
        Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
   dp[9].s = 'K'; dp[9].d = Z;

   X = A.x; Y = A.y - a/2.0 + 4*a/5.0; Z = A.z + h;
        Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
   dp[10].s = 'L'; dp[10].d = Z;

   X = A.x; Y = A.y - a/2.0 + a; Z = A.z + h;
        Point3D_2D(&X,&Y,&Z,B,qx,qy,qz);
   dp[11].s = 'M'; dp[11].d = Z;




    for(i = 0; i <=11; i++)                 /*  深度排序  */
   {
      for(j = i+1;j<=12;j++)
      {
         if(dp[i].d>dp[j].d)
         {
           temp.s = dp[i].s; temp.d = dp[i].d;
           dp[i].s = dp[j].s; dp[i].d = dp[j].d;
           dp[j].s = temp.s; dp[j].d = temp.d;
         }
      }
   }

   i = 0;

   while(i <= 11)      /*画家算法*/
  {
    switch(dp[i++].s)
    {
    case 'A': D.x = A.x - b/2.0; D.y = A.y - a/2.0; D.z = A.z + h/2.0;  /* 四个脚 */
              Cft3D(D,20,20,h+20,B,qx,qy,qz);
              break;


    case 'B': D.x = A.x - b/2.0; D.y = A.y + a/2.0; D.z = A.z + h/2.0;
              Cft3D(D,20,20,h+20,B,qx,qy,qz);
              break;

    case 'C': D.x = A.x + b/2.0; D.y = A.y - a/2.0; D.z = A.z + h/2.0;
              Cft3D(D,20,20,h+20,B,qx,qy,qz);
              break;


    case 'D': D.x = A.x + b/2.0; D.y = A.y + a/2.0; D.z = A.z + h/2.0;
              Cft3D(D,20,20,h+20,B,qx,qy,qz);
              break;


    case 'E': D.x = A.x - b/2.0; D.y = A.y; D.z = A.z + h;  /*上面的条木*/
              Cft3D(D,a+20,10,10,B,qx,qy,qz);
              break;


    case 'F': D.x = A.x + b/2.0; D.y = A.y; D.z = A.z + h;
              Cft3D(D,a+20,10,10,B,qx,qy,qz);
              break;

    case 'H': D.x = A.x; D.y = A.y - a/2.0; D.z = A.z + h;
              Cft3D(D,10,b+20,10,B,qx,qy,qz);
              break;

    case 'I': D.x = A.x; D.y = A.y - a/2.0 + a/5.0; D.z = A.z + h;
              Cft3D(D,10,b+20,10,B,qx,qy,qz);
              break;


    case 'J': D.x = A.x; D.y = A.y - a/2.0 + 2*a/5.0; D.z = A.z + h;
              Cft3D(D,10,b+20,10,B,qx,qy,qz);
              break;


    case 'K': D.x = A.x; D.y = A.y - a/2.0 + 3*a/5.0; D.z = A.z + h;
              Cft3D(D,10,b+20,10,B,qx,qy,qz);
              break;

    case 'L': D.x = A.x; D.y = A.y - a/2.0 + 4*a/5.0; D.z = A.z + h;
              Cft3D(D,10,b+20,10,B,qx,qy,qz);
              break;

    case 'M': D.x = A.x; D.y = A.y - a/2.0 + a; D.z = A.z + h;
              Cft3D(D,10,b+20,10,B,qx,qy,qz);
              break;

    }
  }
}

/*****************************/
/**          房子           **/
/*****************************/
void FangZi(Point3D A,int a,int b,int h,Point3D BB,float qx,float qy,float qz)
{
   int X,Y,Z,i,j,d;
   Point3D B,C;
   Depth dp[9],temp;

   X = A.x; Y = A.y; Z = A.z;         /*深度信息初始化*/
    Point3D_2D(&X,&Y,&Z,BB,qx,qy,qz);
   dp[0].s = 'A'; dp[0].d = Z;        /* 长方体的深度*/

   X = A.x + b/2.0 + 30; Y = A.y; Z = A.z - h/4.0;
     Point3D_2D(&X,&Y,&Z,BB,qx,qy,qz);
   dp[1].s = 'B'; dp[1].d = Z;        /*门框的深度*/

   X = A.x + b/2.0 + a/12.0 + 60; Y = A.y - a/6.0; Z = A.z - h/4.0;
     Point3D_2D(&X,&Y,&Z,BB,qx,qy,qz);
   dp[2].s = 'C'; dp[2].d = Z;        /*门1的深度*/

   X = A.x + b/2.0 + a/12.0 + 60; Y = A.y + a/6.0; Z = A.z - h/4.0;
      Point3D_2D(&X,&Y,&Z,BB,qx,qy,qz);
   dp[3].s = 'D'; dp[3].d = Z;        /*门2的深度*/

   X = A.x; Y = A.y + a/2.0; Z = A.z;
      Point3D_2D(&X,&Y,&Z,BB,qx,qy,qz);
   dp[4].s = 'E'; dp[4].d = Z;       /* 窗框1的深度 */

   X = A.x; Y = A.y + a/2.0 + 25; Z = A.z - h/8.0;
       Point3D_2D(&X,&Y,&Z,BB,qx,qy,qz);
   dp[5].s = 'F'; dp[5].d = Z;       /* 窗台1的深度 */

   X = A.x; Y = A.y - a/2.0; Z = A.z;
       Point3D_2D(&X,&Y,&Z,BB,qx,qy,qz);
   dp[6].s = 'H'; dp[6].d = Z;       /* 窗框2的深度 */

   X = A.x; Y = A.y - a/2.0 - 25; Z = A.z - h/8.0;
       Point3D_2D(&X,&Y,&Z,BB,qx,qy,qz);
   dp[7].s = 'I'; dp[7].d = Z;         /* 窗台2的深度 */

   X = A.x; Y = A.y; Z = A.z + h/2.0;
       Point3D_2D(&X,&Y,&Z,BB,qx,qy,qz);
   dp[8].s = 'J'; dp[8].d = Z;         /*房顶的深度*/

    for(i = 0; i <=8; i++)                 /*  深度排序  */
   {
      for(j = i+1;j<=9;j++)
      {
         if(dp[i].d>dp[j].d)
         {
           temp.s = dp[i].s; temp.d = dp[i].d;
           dp[i].s = dp[j].s; dp[i].d = dp[j].d;
           dp[j].s = temp.s; dp[j].d = temp.d;
         }
      }
   }

  i = 0;
  while(i <= 9)      /*画家算法*/
  {
    switch(dp[i++].s)
    {

       case'A': Cft3D(A,a,b,h,BB,qx,qy,qz); /* 长方体 */
                break;

       case'B': B.x = A.x + b/2.0; B.y = A.y; B.z = A.z - h/4.0; /* 画框 */
                Rect_3DYZ(B,a/3.0,h/2.0,15,BB,qx,qy,qz);
                break;

       case'C': B.x = A.x + b/2.0 + a/12.0; B.y = A.y - a/6.0; B.z = A.z - h/4.0; /* 画门1*/
                C.x = A.x + b/2.0; C.y = B.y; C.z = B.z;
                Cft3D_XuanZhuan(B,5,a/6.0,h/2.0,0,C,3,BB,qx,qy,qz);
                break;

       case'D': B.x = A.x + b/2.0 + a/12.0; B.y = A.y + a/6.0; B.z = A.z - h/4.0;  /* 画门2*/
                C.x = A.x + b/2.0; C.y = B.y; C.z = B.z;
                Cft3D_XuanZhuan(B,5,a/6.0,h/2.0,0,C,3,BB,qx,qy,qz);
                break;

       case'E': B.x = A.x; B.y = A.y + a/2.0; B.z = A.z; /* 画窗框1 */
                Rect_3DXZ(B,b/2.0,h/4.0,15,BB,qx,qy,qz);
                Rect_3DXZ(B,5,h/4.0,12,BB,qx,qy,qz); /* 竖线 */
                break;

       case'F': B.x = A.x; B.y = A.y + a/2.0 + 25; B.z = A.z - h/8.0; /* 窗台1*/
                Cft3D(B,50,b/2.0,5,BB,qx,qy,qz);
                break;

       case'H': B.x = A.x; B.y = A.y - a/2.0; B.z = A.z; /* 画窗框2 */
                Rect_3DXZ(B,b/2.0,h/4.0,15,BB,qx,qy,qz);
                Rect_3DXZ(B,5,h/4.0,12,BB,qx,qy,qz); /* 竖线 */
                break;

       case'I': B.x = A.x; B.y = A.y - a/2.0 - 25; B.z = A.z - h/8.0; /* 窗台2*/
                Cft3D(B,50,b/2.0,5,BB,qx,qy,qz);
                break;

       case'J':  B.x = A.x; B.y = A.y; B.z = A.z + h/2.0;/* 房顶 */
                 Rect_3DXY(B,a + 100,b + 100,8,BB,qx,qy,qz);

                 B.x = A.x - b/2.0 - 50; B.y = A.y; B.z = A.z + h/2.0 + (b+100)/4.0;
                 C.x = B.x; C.y = B.y; C.z = A.z + h/2.0;
                 d = (b/2.0 + 50) * 2.0/sqrt(3);
                 Rect_3DYZ_XuanZhuan(B,a + 100,d,9,70,C,2,BB,qx,qy,qz);

                 B.x = A.x + b/2.0 + 50; B.y = A.y; B.z = A.z + h/2.0 + (b+100)/4.0;
                 C.x = B.x; C.y = B.y; C.z = A.z + h/2.0;
                 d = (b/2.0 + 50) * 2.0/sqrt(3);
                 Rect_3DYZ_XuanZhuan(B,a + 100,d,9,-70,C,2,BB,qx,qy,qz);

                                                           /* 瓦*/
                 B.x = A.x - b/2.0 - 50; B.y = A.y - a/2.0 - 50; B.z = A.z + h/2.0 + (b+100)/4.0;
                 C.x = B.x; C.y = B.y; C.z = A.z + h/2.0;
                 while(B.y <= A.y + a/2.0 + 50)
                {
                 Rect_3DYZ_XuanZhuan(B,5,d,12,70,C,2,BB,qx,qy,qz);
                 B.y = B.y + 40;
                 C.y = C.y + 40;

                 }

                  B.x = A.x + b/2.0 + 50; B.y = A.y - a/2.0 - 50; B.z = A.z + h/2.0 + (b+100)/4.0;
                  C.x = B.x; C.y = B.y; C.z = A.z + h/2.0;
                  while(B.y <= A.y + a/2.0 + 50)
                    {
                    Rect_3DYZ_XuanZhuan(B,5,d,12,-70,C,2,BB,qx,qy,qz);
                    B.y = B.y + 40;
                    C.y = C.y + 40;

                    }
                                                       /* 屋脊 */
                  B.x = A.x; B.y = A.y; B.z = A.z + h/2.0 + 45;
                  Cft3D(B,a + 100,10,10,BB,qx,qy,qz);
                  break;
        }
    }

}


/***************************/
/**　　　　　　　　　　　 **/
/**       3D综合画面      **/
/**                       **/
/***************************/
void My3DGraphics_A(Point3D A,Point3D BB,float qx,float qy,float qz) /* 观察方向为(qx,qy,qz)*/
{
   float a,b,h,d;
   int X,Y,Z,i,j;
   Point3D B,C;
   Depth dp[20],temp;
   setbkcolor(2);

   a = 220;
   b = 150;
   h = 200;

   B.x = A.x + b/2.0 + 250; B.y = A.y + a/2.0 + 100; B.z = A.z - 100;
   C.x = A.x + b/2.0 + 600; C.y = A.y; C.z = A.z - 100;
   i = 0;
   while(B.y>=A.x - a/2.0 - 1000)   /*亭前小路*/
   {
       Rect_3DXY_XuanZhuan(B,50,100,7,i,C,3,BB,qx,qy,qz);
       B.y = B.y - 80;
       i = i + 5;
   }


                                  /*深度信息初始化*/

   X = A.x; Y = A.y; Z = A.z;            /*底台的深度*/
      Point3D_2D(&X,&Y,&Z,BB,qx,qy,qz);
   dp[0].s = 'A'; dp[0].d = Z;

   X = A.x; Y = A.y; Z = A.z + 100 + h/2.0;    /*房子的深度*/
      Point3D_2D(&X,&Y,&Z,BB,qx,qy,qz);
   dp[1].s = 'B'; dp[1].d = Z;

   X = A.x - a/2.0 - 500; Y = A.y + b/2.0 + 350; Z = A.z+250;
      Point3D_2D(&X,&Y,&Z,BB,qx,qy,qz); /* 树1的深度*/
   dp[2].s = 'C'; dp[2].d = Z;

   X = A