【文章內(nèi)容簡介】 1 。 //獲取第二個點編號 j3 * pPlanes+ip pW+2 。 //獲取第三個點編號 x1 pPoints+j21 x pPoints+j11 x。 y1 pPoints+j21 y pPoints+j11 y。 z1 pPoints+j21 z pPoints+j11 z。 x2 pPoints+j31 x pPoints+j21 x。 y2 pPoints+j31 y pPoints+j21 y。 z2 pPoints+j31 z pPoints+j21 z。 x0 y1*z2y2*z1。 //計算法向量 y0 x2*z1x1*z2。 z0 x1*y2x2*y1。 r float sqrt x0*x0+y0*y0+z0*z0 。 pPlanes+ip S .x x0/r。 //法向量 單位化 pPlanes+ip S .y y0/r。 pPlanes+ip S .z z0/r。 ///下面給邊唯一編號，并與點、面對應(yīng) ////////// int i,kk。 for i 0。i pPlanes+ip wN。i++ //外環(huán) 邊編號 j1 * pPlanes+ip pW+i 。 if i pPlanes+ip wN1 //獲取兩點 j1,j2 j2 * pPlanes+ip pW 。 else j2 * pPlanes+ip pW+i+1 。 for kk 0。kk k。kk++ //檢查前 k1 個邊里是否有 j1,j2 組成的邊 if pBorders+kk p1 j1amp。amp。 pBorders+kk p2 j2|| pBorders+kk p1 j2amp。amp。 pBorders+kk p2 j1 kk kk1。 //表明 j1,j2 組成的邊已存在 break。 if kk 0 * pPlanes+ip pBw+i kk。 //設(shè)置此面中第 i+1 個外環(huán)邊的編號為 kk+1 continue。 pBorders+k p1 j1。 //置 j1,j2 組成的邊編號 為k+1,即再 pBorders[]里下標為 k pBorders+k p2 j2。 * pPlanes+ip pBw+i k+1。 //設(shè)置此面中第 i+1 個外環(huán)邊的編號為 k+1 k++。 for i 0。i pPlanes+ip nN。i++ //內(nèi)環(huán)邊編號 j1 * pPlanes+ip pN+i 。 if i pPlanes+ip nN1 //獲取兩點 j1,j2 j2 * pPlanes+ip pN 。 else j2 * pPlanes+ip pN+i+1 。 for kk 0。kk k。kk++ //檢查前 k1 個邊里是否有j1,j2 組成的邊 if pBorders+kk p1 j1amp。amp。 pBorders+kk p2 j2|| pBorders+kk p1 j2amp。amp。 pBorders+kk p2 j1 kk kk1。 //表明 j1,j2 組成的邊已存在 break。 if kk 0 * pPlanes+ip pBw+i kk。 //設(shè)置此面中第 i+1 個內(nèi)環(huán)邊的編號為 kk+1 continue。 pBorders+k p1 j1。 //置 j1,j2 組成的邊編號為 k+1,即再 pBorders[]里下標為 k pBorders+k p2 j2。 * pPlanes+ip pBn+i k+1。 //設(shè)置此面中第 i+1 個內(nèi)環(huán)邊的編號為 k+1 k++。 //////////////////////////////////////////////////// void ThreeD::MarkPlane TDPoint amp。viewAngle for int i 0。i bNpN+2。i++ if * pPlanes+i + * pPlanes+i + * pPlanes+i 0 pPlanes+i mark true。 else pPlanes+i mark false。 //////////////////////////////////////////////////// void ThreeD::paintPlane int N,TDPoint amp。viewPoint,TDPoint amp。viewAngle,CClientDC amp。dc,float L if pPlanes+N1 mark false return。 int i,j1,j2,sign 0。 float x,y,z,x0,y0,z0,xv,yv,zv,xa,ya,za,lamd。 float yyx,yyy,yyz,xx0,yy0,zz0,xp,yp。 float cosa,sina,cosb,sinb,opx,opy,opz,opr。 xv 。 xa 。 yv 。 ya 。 zv 。 za 。 xx0 xv+L*xa。 yy0 yv+L*ya。 //獲取視野坐標原點 zz0 zv+L*za。 yyx xa*za。 //視野 x 軸向量 : ya,xa,0 ,不是單位向量 yyy ya*za。 //獲取視野 y 軸向量 不是單位向量 yyz ya*ya+xa*xa。 int bkk。 for i 1。i pPlanes+N1 wN。i++ //繪制外邊 j1 * pPlanes+N1 pW+i1 。 if i pPlanes+N1 wN j2 * pPlanes+N1 pW 。 else j2 * pPlanes+N1 pW+i 。 if pPoints+j11 mark false //如果此點還沒有做過投影 x pPoints+j11 x。 y pPoints+j11 y。 //獲取第一點的三維坐標 z pPoints+j11 z。 lamd L/ xa* xxv +ya* yyv +za* zzv 。 if lamd 0 x0 lamd*x+ 1lamd *xv。 y0 lamd*y+ 1lamd *yv。 //第一投影點在三維世界的坐標 z0 lamd*z+ 1lamd *zv。 opx x0xx0。 opy y0yy0。 opz z0zz0。 opr float sqrt opx*opx+opy*opy+opz*opz 。 cosa ya*opxxa*opy / float sqrt ya*ya+xa*xa /opr。 cosb yyx*opx+yyy*opy+yyz*opz / float sqrt yyx*yyx+yyy*yyy+yyz*yyz /opr。 sina float sqrt 1cosa*cosa 。 sinb float sqrt 1cosb*cosb 。 if cosa 0 sign 1。else sign 1。 xp opr*sinb*sign。 //第一點在視野中的 x 坐標 if cosb 0 sign 1。else sign 1。 yp opr*sina*sign。 //第一點在視野中的 y 坐標 * pShodowX+j11 360+int xp 。 * pShodowY+j11 260int yp 。 pPoints+j11 mark true。 //標明此點已作了投影 if pPoints+j21 mark false //如果此點還沒有做過投影 x pPoints+j21 x。 y pPoints+j21 y。 //獲取第二點的三維坐標 z pPoints+j21 z。 lamd L/ xa* xxv +ya* yyv +za* zzv 。 if lamd 0 x0 lamd*x+ 1lamd *xv。 y0 lamd*y+ 1lamd *yv。 //第二投影點在三維世界的坐標 z0 lamd*z+ 1lamd *zv。 opx x0xx0。 opy y0yy0。 opz z0zz0。 opr float sqrt opx*opx+opy*opy+opz*opz 。 cosa ya*opxxa*opy / float sqrt ya*ya+xa*xa /opr。 cosb yyx*opx+yyy*opy+yyz*opz / float sqrt yyx*yyx+yyy*yyy+yyz*yyz /opr。 sina float sqrt 1cosa*cosa 。 sinb float sqrt 1cosb*cosb 。 if cosa 0 sign 1。else sign 1。 xp opr*sinb*sign。 //第一點在視野中的 x 坐標 if cosb 0 sign 1。else sign 1。 yp opr*sina*sign。 //第一點在視野中的 y 坐標 * pShodowX+j21 360+int xp 。 * pShodowY+j21 260int yp 。 pPoints+j21 mark true。 //標明此點已作了投影 bkk * pPlanes+N1 pBw+i1 。 //獲取此邊編號 if pBorders+bkk1 mark false //如果此邊沒有繪制過 * pShodowX+j11 ,* pShodowY+j11 。 * pShodowX+j21 ,* pShodowY+j21 。 //繪制此面外環(huán)中的第 i 條邊 pBorders+bkk1 mark true。 //標明此邊已繪制 /////////////////////////////////////////外邊繪制完畢 for i 1。i pPlanes+N1 nN。i++ //繪制內(nèi)邊 j1 * pPlanes+N1 pN+i1 。 if i pPlanes+N1 nN j2 * pPlanes+N1 pN 。 else j2 * pPlanes+N1 pN+i 。