开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > GIS编程 > 查看源码
DelaunayView.cpp
资源名称：VC++Delaunay.zip　[点击查看]
上传用户：azhong891
上传日期：2013-06-04
资源大小：197k
文件大小：36k
源码类别：

GIS编程

开发平台：

Visual C++

DelaunayView.cpp：源码内容
  1. // DelaunayView.cpp : implementation of the CDelaunayView class
  2. //
  4. #include "stdafx.h"
  5. #include "Delaunay.h"
  7. #include "mainfrm.h"
  8. #include "pointpos.h"
  9. #include "triangle.h"
  11. #include "DelaunayDoc.h"
  12. #include "DelaunayView.h"
  14. #ifdef _DEBUG
  15. #define new DEBUG_NEW
  16. #undef THIS_FILE
  17. static char THIS_FILE[] = __FILE__;
  18. #endif
  20. /////////////////////////////////////////////////////////////////////////////
  21. // CDelaunayView
  23. IMPLEMENT_DYNCREATE(CDelaunayView, CView)
  25. BEGIN_MESSAGE_MAP(CDelaunayView, CView)
  26. //{{AFX_MSG_MAP(CDelaunayView)
  27. ON_WM_SIZE()
  28. ON_WM_CREATE()
  29. ON_WM_LBUTTONDOWN()
  30. ON_WM_MOUSEMOVE()
  31. ON_WM_DESTROY()
  32. ON_WM_KEYDOWN()
  33. ON_COMMAND(IDM_CORLOR, OnCorlorMaterial)
  34. ON_UPDATE_COMMAND_UI(IDM_CORLOR, OnUpdateCorlorMaterial)
  35. ON_UPDATE_COMMAND_UI(IDM_LIGHT, OnUpdateLight)
  36. ON_COMMAND(IDM_LIGHT, OnLight)
  37. ON_COMMAND(ID_BUTTON_TWO, OnButtonTwo)
  38. ON_UPDATE_COMMAND_UI(ID_BUTTON_TWO, OnUpdateButtonTwo)
  39. ON_COMMAND(ID_BUTTON_THREE, OnButtonThree)
  40. ON_UPDATE_COMMAND_UI(ID_BUTTON_THREE, OnUpdateButtonThree)
  41. ON_COMMAND(ID_BUTTON_HCT, OnButtonHct)
  42. ON_UPDATE_COMMAND_UI(ID_BUTTON_HCT, OnUpdateButtonHct)
  43. ON_COMMAND(ID_BUTTON_FILL, OnButtonFill)
  44. ON_UPDATE_COMMAND_UI(ID_BUTTON_FILL, OnUpdateButtonFill)
  45. ON_COMMAND(ID_BUTTON_LINE, OnButtonLine)
  46. ON_UPDATE_COMMAND_UI(ID_BUTTON_LINE, OnUpdateButtonLine)
  47. ON_COMMAND(ID_BUTTON_WNAG, OnButtonWnag)
  48. //}}AFX_MSG_MAP
  49. END_MESSAGE_MAP()
  51. /////////////////////////////////////////////////////////////////////////////
  52. // CDelaunayView construction/destruction
  54. CDelaunayView::CDelaunayView()
  55. {
  56. // TODO: add construction code here
  58. }
  60. CDelaunayView::~CDelaunayView()
  61. {
  62. }
  64. BOOL CDelaunayView::PreCreateWindow(CREATESTRUCT& cs)
  65. {
  66. // TODO: Modify the Window class or styles here by modifying
  67. //  the CREATESTRUCT cs
  68.     cs.style |=WS_CLIPCHILDREN |WS_CLIPSIBLINGS;
  70. return CView::PreCreateWindow(cs);
  71. }
  73. /////////////////////////////////////////////////////////////////////////////
  74. // CDelaunayView drawing
  76. void CDelaunayView::OnDraw(CDC* pDC)
  77. {
  78. CDelaunayDoc* pDoc = GetDocument();
  79. ASSERT_VALID(pDoc);
  81. DrawScene();
  82. SwapBuffers(pDC->m_hDC);
  83. }
  85. /////////////////////////////////////////////////////////////////////////////
  86. // CDelaunayView diagnostics
  88. #ifdef _DEBUG
  89. void CDelaunayView::AssertValid() const
  90. {
  91. CView::AssertValid();
  92. }
  94. void CDelaunayView::Dump(CDumpContext& dc) const
  95. {
  96. CView::Dump(dc);
  97. }
  99. CDelaunayDoc* CDelaunayView::GetDocument() // non-debug version is inline
  100. {
  101. ASSERT(m_pDocument->IsKindOf(RUNTIME_CLASS(CDelaunayDoc)));
  102. return (CDelaunayDoc*)m_pDocument;
  103. }
  104. #endif //_DEBUG
  106. /////////////////////////////////////////////////////////////////////////////
  107. // CDelaunayView message handlers
  109. void CDelaunayView::OnSize(UINT nType, int cx, int cy) 
  110. {
  111. CView::OnSize(nType, cx, cy);
  112. if(cy>0)
  113. {
  114. m_oldRect.right=cx;
  115. m_oldRect.bottom=cy;
  116. glMatrixMode(GL_PROJECTION);
  117. glLoadIdentity();
  118. glFrustum(-2.0,2.0,-2.0,2.0,0.0,20.0);
  119. glLoadIdentity();//????不能不要,事关显示
  120. //gluPerspective(40.0f,(GLdouble)cx/cy,1.0f,50.0f);
  121. glMatrixMode(GL_MODELVIEW);
  122. glViewport(0,0,cx,cy);
  123. }
  125. }
  127. int CDelaunayView::OnCreate(LPCREATESTRUCT lpCreateStruct) 
  128. {
  129. if (CView::OnCreate(lpCreateStruct) == -1)
  130. return -1;
  131. InitOpenGl();
  132. m_pDoc=this->GetDocument();
  133. m_ox=50;m_oy=50;m_len=400;
  134. m_Dimension=DO_DRAW2;
  135. //m_HCT=DO_NO_HCT;
  136. m_Draw_what=DO_LINE;
  137. //OpenGl/////////////////////////////
  138. m_translate_x=-1.0;
  139. m_translate_y=-1.0;
  140. m_translate_z=0;
  141. m_rotate_angle=0;
  142. m_gl_x=1.0;m_gl_y=1.0;m_gl_z=0.0;
  143. eye.x=4.0;eye.y=5.0;eye.z=6.0;
  145. m_colorRed=float(0.2);m_colorGreen=float(0.5);m_colorBlue=float(0.8);
  146. m_bMaskRed=FALSE;m_bMaskGreen=FALSE;m_bMaskBlue=FALSE;
  148. m_materialAmb[0]=0.2;m_materialAmb[1]=0.2;m_materialAmb[2]=0.2;m_materialAmb[3]=1.0;
  149. m_materialDif[0]=0.8;m_materialDif[1]=0.8;m_materialDif[2]=0.8;m_materialDif[3]=1.0;
  150. m_materialSpe[0]=0.0;m_materialSpe[1]=0.0;m_materialSpe[2]=0.0;m_materialSpe[3]=1.0;
  151. m_materialEmi[0]=0.0;m_materialEmi[1]=0.0;m_materialEmi[2]=0.0;m_materialEmi[3]=1.0;
  152. m_matshininess=0.0f;
  154. m_lightAmb[0]=0.2;m_lightAmb[1]=0.2;m_lightAmb[2]=0.2;m_lightAmb[3]=1.0;
  155. m_lightDif[0]=1.0;m_lightDif[1]=1.0;m_lightDif[2]=0.3;m_lightDif[3]=1.0;
  156. m_lightSpe[0]=1.0;m_lightSpe[1]=1.0;m_lightSpe[2]=1.0;m_lightSpe[3]=1.0;
  157. m_lightPos[0]=1.0;m_lightPos[1]=1.0;m_lightPos[2]=1.0;m_lightPos[3]=0.0;
  158. //OpenGl/////////////////////////////
  159. return 0;
  160. }
  162. void CDelaunayView::OnLButtonDown(UINT nFlags, CPoint point) 
  163. {////////////////////////////////////////////
  164. ASSERT(m_pDoc!=NULL);
  165. CPoint pos=point;
  166. int m_plen;//当前节点个数
  167. double s;
  168. //ScreenToClient(&pos);
  169. if(m_pDoc->m_DoWhat==DO_ADD)
  170. {
  171. //'x' and 'y' are that to be stored in 'm_point'
  172. double x,y;
  173. x=double(pos.x)/double(m_oldRect.right);
  174. y=double((m_oldRect.bottom-pos.y))/double(m_oldRect.bottom);
  175. //form the original convexity when 'm_plen' is smaller than 3
  176. //the original convexity is a triangle
  177. m_plen=m_pDoc->m_point.GetSize();//==1????
  178. switch(m_plen) 
  179. {
  180. case 0 :
  181. m_pDoc->AddPoint(x,y);
  182. m_pDoc->m_con.SetAtGrow(0,0);
  183.                         break;
  184. case 1 :
  185. if(m_pDoc->m_point[0]->m_x==x &&
  186.    m_pDoc->m_point[0]->m_y==y)
  187. {
  188.   AfxMessageBox("The two point not suitable");
  189.                       break;
  190. //please 输入节点,跳出
  191. }
  192. m_pDoc->AddPoint(x,y);
  193. m_pDoc->m_con.SetAtGrow(1,1);
  194.                         break;
  195.                     case 2 ://确保凸包顶点逆时针存储
  196. m_pDoc->AddPoint(x,y);
  197. s=m_pDoc->S(0,1,2);
  198. if(s==0)
  199. {
  200. AfxMessageBox("please 输入节点 again,因三点共线,跳出");
  201. m_pDoc->m_point.RemoveAt(2,1);
  202.                           //  m_pDoc->m_n.RemoveAt(2,1);
  203. }
  205. m_plen=m_pDoc->m_point.GetSize();
  206. m_pDoc->m_con.SetAtGrow(m_plen-1,m_plen-1);                    
  207. if(s<0) 
  208. {
  209. m_pDoc->m_con[0]=1;
  210. m_pDoc->m_con[1]=0;                          
  211. }
  212. CTriangle* pTriangle;
  213. pTriangle=new CTriangle(m_pDoc->m_con[0],
  214.       m_pDoc->m_con[1],m_pDoc->m_con[2]);
  215. m_pDoc->Center(pTriangle);
  216. m_pDoc->BaryCenter(pTriangle);
  217. m_pDoc->m_tri.AddHead(pTriangle);                   
  218.                         break;
  219.                 default :
  220. m_pDoc->AddPoint(x,y);
  221.                         m_plen=m_pDoc->m_point.GetSize();
  222.                         int i=DelTriList(x,y,m_plen-1);
  223. if(i==POS_ERROR){
  224. int k=m_pDoc->m_point.GetSize();
  225. m_pDoc->m_point.RemoveAt(k-1,1);
  226. }
  227.                         break;
  228. }
  229. Invalidate();
  230. }
  231. ////////////////////////////////////////////
  232. CView::OnLButtonDown(nFlags, point);
  233. }
  235. void CDelaunayView::OnMouseMove(UINT nFlags, CPoint point) 
  236. {
  237. // TODO: Add your message handler code here and/or call default
  238. // CClientDC dc(this);
  239. double x,y;
  240. double a,b;
  241. a=double(point.x);
  242. b=double(m_oldRect.right);
  243. x=a/b;
  244. y=double((m_oldRect.bottom-point.y))/double(m_oldRect.bottom);
  246.     CString str1;
  247. str1.Empty();
  248. int i=m_pDoc->m_DoWhat;
  249.     str1.Format("%f;%f;%d;%d;%d;%d",x,y,i,m_Dimension,m_HCT,m_Draw_what);
  250. CMainFrame *pFram=(CMainFrame *)(AfxGetApp()->GetMainWnd());
  251. pFram->m_wndStatusBar.SetPaneText(0,str1);
  252. CView::OnMouseMove(nFlags, point);
  253. }
  256. array CDelaunayView::Wher(CPointPos *pos)
  257. {
  258. //Points and Convexity
  259. /*return Value:a.a
  260. POS_IN=1,a point  belong a triangle or convexity
  261. POS_ON=2,a point belong a circle or on the edge of convexity
  262. POS_OUT=0,a point out of a circle or convexity
  263. a.b , a.c: record the point's mark in m_con */
  264. //1. get a ponit in convexity:p
  265. array a;
  266. int max=m_pDoc->m_con.GetSize();
  267. int z=max/2;
  268. CPointPos *p1=new CPointPos(m_pDoc->m_point[m_pDoc->m_con[0]]->m_x,
  269. m_pDoc->m_point[m_pDoc->m_con[0]]->m_y);
  270.     CPointPos *p2=new CPointPos(m_pDoc->m_point[m_pDoc->m_con[z]]->m_x,
  271. m_pDoc->m_point[m_pDoc->m_con[z]]->m_y);
  272.     CPointPos *p3=new CPointPos(m_pDoc->m_point[m_pDoc->m_con[max-1]]->m_x,
  273. m_pDoc->m_point[m_pDoc->m_con[max-1]]->m_y);
  274. CPointPos *p=new CPointPos((p1->m_x+p2->m_x+p3->m_x)/3,
  275.          (p1->m_y+p2->m_y+p3->m_y)/3);
  276. p3->m_x=double(pos->m_x);
  277. p3->m_y=double(pos->m_y);
  278. //遍寻Convexity
  279. double s1,s2,s;
  280. for(int i=0;i<max-1;i++)
  281. {
  282. p1->m_x=m_pDoc->m_point[m_pDoc->m_con[i]]->m_x;
  283.     p1->m_y=m_pDoc->m_point[m_pDoc->m_con[i]]->m_y;
  284.       p2->m_x=m_pDoc->m_point[m_pDoc->m_con[i+1]]->m_x;
  285.     p2->m_y=m_pDoc->m_point[m_pDoc->m_con[i+1]]->m_y;
  286.         
  287. s1=m_pDoc->S(p3,p,p1);
  288. s2=m_pDoc->S(p3,p,p2);
  289. if(s1>0 && s2<0)//the point in the 扇形 con[i],con[i+1] and p formed
  290. {
  291.   a.b=i;
  292.   a.c=i+1;
  293.           s=m_pDoc->S(p1,p2,p3);
  294.   if(s>0)
  295.   {
  296.   a.a=1;//POS_IN  
  297.   }
  298.           if(s==0)
  299.   {
  300.   a.a=2;//POS_ON
  301.   }
  302.   if(s<0)
  303.   {
  304.   a.a=0;//POS_OUT
  305.   }
  306.   return a;
  307. }         
  308. }
  309.     p1->m_x=m_pDoc->m_point[m_pDoc->m_con[i]]->m_x;
  310. p1->m_y=m_pDoc->m_point[m_pDoc->m_con[i]]->m_y;
  311.     p2->m_x=m_pDoc->m_point[m_pDoc->m_con[0]]->m_x;
  312. p2->m_y=m_pDoc->m_point[m_pDoc->m_con[0]]->m_y;
  313. s1=m_pDoc->S(p3,p,p1);
  314. s2=m_pDoc->S(p3,p,p2);
  315. if(s1>0 && s2<0)
  316. {
  317.           a.b=i;
  318.   a.c=0;
  319.           s=m_pDoc->S(p1,p2,p3);
  320.   if(s>0)
  321.   {
  322.   a.a=1;//POS_IN  
  323.   }
  324.           if(s==0)
  325.   {
  326.   a.a=2;//POS_ON
  327.   }
  328.   if(s<0)
  329.   {
  330.   a.a=0;//POS_OUT
  331.   }
  332.   return a;
  333. }           
  334. }
  336. void CDelaunayView::DrawTri(CTriangle *tri)
  337. {
  338. int i;
  339. //POI normal;
  340. int p1,p2,p3;
  342. p1=tri->m_p1;
  343. p2=tri->m_p2;
  344. p3=tri->m_p3;
  345. GLfloat fLineWidth[2];
  346. glGetFloatv(GL_LINE_WIDTH_RANGE,fLineWidth);
  347. if(m_Dimension==DO_DRAW2 && m_HCT!=DO_HCT)
  348. {
  349. glColor4d(0.0,0.0,1.0,0.5);
  350.     glBegin(GL_TRIANGLES);
  351. glVertex3d(m_pDoc->m_point[p1]->m_x,m_pDoc->m_point[p1]->m_y,0.0); 
  352.      glVertex3d(m_pDoc->m_point[p2]->m_x,m_pDoc->m_point[p2]->m_y,0.0);
  353.      glVertex3d(m_pDoc->m_point[p3]->m_x,m_pDoc->m_point[p3]->m_y,0.0);
  354.      glEnd();
  355. }
  356. if(m_Dimension==DO_DRAW3 && m_HCT!=DO_HCT)
  357. {
  358. glPushMatrix();
  359. //glScaled(2,2,1.5);
  360. glColor3f(1.0f,0.0f,0.0f);         
  361. glBegin(GL_TRIANGLES);
  362. glVertex3d(m_pDoc->m_point[p1]->m_x,m_pDoc->m_point[p1]->m_y,m_pDoc->m_point[p1]->m_z); 
  363.      glVertex3d(m_pDoc->m_point[p2]->m_x,m_pDoc->m_point[p2]->m_y,m_pDoc->m_point[p2]->m_z);
  364.      glVertex3d(m_pDoc->m_point[p3]->m_x,m_pDoc->m_point[p3]->m_y,m_pDoc->m_point[p3]->m_z);
  365.      glEnd();
  366. glPopMatrix();
  367. }
  368. if(m_Dimension==DO_DRAW3 && m_HCT==DO_HCT)
  369. {
  370. m_pDoc->DrawTri(tri->m_p1,tri->m_p2,tri);
  371. m_pDoc->DrawTri(tri->m_p2,tri->m_p3,tri);
  372. m_pDoc->DrawTri(tri->m_p3,tri->m_p1,tri);
  373. /*for(i=1;i<=3;i++)
  374. {
  375. HCT(tri,i);
  376. glColor4d(0.0,1.0,0.0,0.5);
  377. if(i==2)
  378. {
  379. glBegin(GL_TRIANGLES);//anticlock    
  380.    normal=GetTriNormal(m_hct[0],m_hct[1],m_hct[2]);
  381.    glNormal3d(normal.x,normal.y,normal.z);
  382.    glVertex3d(m_hct[0].x,m_hct[0].y,m_hct[0].z);          
  383.    glVertex3d(m_hct[1].x,m_hct[1].y,m_hct[1].z);
  384.    glVertex3d(m_hct[2].x,m_hct[2].y,m_hct[2].z);
  385.    
  386.    normal=GetTriNormal(m_hct[1],m_hct[3],m_hct[2]);
  387.    glNormal3d(normal.x,normal.y,normal.z);
  388.    glVertex3d(m_hct[1].x,m_hct[1].y,m_hct[1].z);
  389.    glVertex3d(m_hct[3].x,m_hct[3].y,m_hct[3].z);
  390.    glVertex3d(m_hct[2].x,m_hct[2].y,m_hct[2].z);
  392.    normal=GetTriNormal(m_hct[3],m_hct[1],m_hct[4]);
  393.    glNormal3d(normal.x,normal.y,normal.z);    
  394.    glVertex3d(m_hct[3].x,m_hct[3].y,m_hct[3].z);
  395.    glVertex3d(m_hct[1].x,m_hct[1].y,m_hct[1].z);
  396.    glVertex3d(m_hct[4].x,m_hct[4].y,m_hct[4].z);
  398.    normal=GetTriNormal(m_hct[3],m_hct[4],m_hct[5]);
  399.    glNormal3d(normal.x,normal.y,normal.z);
  400.    glVertex3d(m_hct[3].x,m_hct[3].y,m_hct[3].z);
  401.    glVertex3d(m_hct[4].x,m_hct[4].y,m_hct[4].z);
  402.    glVertex3d(m_hct[5].x,m_hct[5].y,m_hct[5].z);
  403.    
  404.    normal=GetTriNormal(m_hct[4],m_hct[6],m_hct[5]);
  405.    glNormal3d(normal.x,normal.y,normal.z);
  406.    glVertex3d(m_hct[4].x,m_hct[4].y,m_hct[4].z);    
  407.    glVertex3d(m_hct[6].x,m_hct[6].y,m_hct[6].z);
  408.         glVertex3d(m_hct[5].x,m_hct[5].y,m_hct[5].z);
  410.    normal=GetTriNormal(m_hct[8],m_hct[7],m_hct[9]);
  411.    glNormal3d(normal.x,normal.y,normal.z);
  412.    glVertex3d(m_hct[8].x,m_hct[8].y,m_hct[8].z); 
  413.         glVertex3d(m_hct[7].x,m_hct[7].y,m_hct[7].z);
  414.    glVertex3d(m_hct[9].x,m_hct[9].y,m_hct[9].z);
  416.    normal=GetTriNormal(m_hct[8],m_hct[3],m_hct[7]);
  417.    glNormal3d(normal.x,normal.y,normal.z);
  418.    glVertex3d(m_hct[8].x,m_hct[8].y,m_hct[8].z);
  419.    glVertex3d(m_hct[3].x,m_hct[3].y,m_hct[3].z);
  420.    glVertex3d(m_hct[7].x,m_hct[7].y,m_hct[7].z);
  422.    normal=GetTriNormal(m_hct[3],m_hct[5],m_hct[7]);
  423.    glNormal3d(normal.x,normal.y,normal.z);
  424.    glVertex3d(m_hct[3].x,m_hct[3].y,m_hct[3].z);
  425.    glVertex3d(m_hct[5].x,m_hct[5].y,m_hct[5].z);
  426.    glVertex3d(m_hct[7].x,m_hct[7].y,m_hct[7].z);
  427.     
  428.    normal=GetTriNormal(m_hct[2],m_hct[3],m_hct[8]);
  429.    glNormal3d(normal.x,normal.y,normal.z);
  430.    glVertex3d(m_hct[2].x,m_hct[2].y,m_hct[2].z); 
  431.    glVertex3d(m_hct[3].x,m_hct[3].y,m_hct[3].z);
  432.    glVertex3d(m_hct[8].x,m_hct[8].y,m_hct[8].z);
  433.      glEnd();
  434. }
  435. else
  436. {
  437.    glBegin(GL_TRIANGLES);//anticlock    
  438.    normal=GetTriNormal(m_hct[0],m_hct[2],m_hct[1]);
  439.    glNormal3d(normal.x,normal.y,normal.z);
  440.    glVertex3d(m_hct[0].x,m_hct[0].y,m_hct[0].z);          
  441.    glVertex3d(m_hct[2].x,m_hct[2].y,m_hct[2].z);
  442.    glVertex3d(m_hct[1].x,m_hct[1].y,m_hct[1].z);
  443.    
  444.    normal=GetTriNormal(m_hct[1],m_hct[2],m_hct[3]);
  445.    glNormal3d(normal.x,normal.y,normal.z);
  446.    glVertex3d(m_hct[1].x,m_hct[1].y,m_hct[1].z);
  447.    glVertex3d(m_hct[2].x,m_hct[2].y,m_hct[2].z);
  448.    glVertex3d(m_hct[3].x,m_hct[3].y,m_hct[3].z);
  450.    normal=GetTriNormal(m_hct[1],m_hct[3],m_hct[4]);
  451.    glNormal3d(normal.x,normal.y,normal.z);    
  452.    glVertex3d(m_hct[1].x,m_hct[1].y,m_hct[1].z);
  453.    glVertex3d(m_hct[3].x,m_hct[3].y,m_hct[3].z);
  454.    glVertex3d(m_hct[4].x,m_hct[4].y,m_hct[4].z);
  456.    normal=GetTriNormal(m_hct[4],m_hct[3],m_hct[5]);
  457.    glNormal3d(normal.x,normal.y,normal.z);
  458.    glVertex3d(m_hct[4].x,m_hct[4].y,m_hct[4].z);
  459.    glVertex3d(m_hct[3].x,m_hct[3].y,m_hct[3].z);
  460.    glVertex3d(m_hct[5].x,m_hct[5].y,m_hct[5].z);
  461.    
  462.    normal=GetTriNormal(m_hct[4],m_hct[5],m_hct[6]);
  463.    glNormal3d(normal.x,normal.y,normal.z);
  464.    glVertex3d(m_hct[4].x,m_hct[4].y,m_hct[4].z);    
  465.    glVertex3d(m_hct[5].x,m_hct[5].y,m_hct[5].z);
  466.         glVertex3d(m_hct[6].x,m_hct[6].y,m_hct[6].z);
  468.    normal=GetTriNormal(m_hct[7],m_hct[8],m_hct[9]);
  469.    glNormal3d(normal.x,normal.y,normal.z);
  470.    glVertex3d(m_hct[7].x,m_hct[7].y,m_hct[7].z); 
  471.         glVertex3d(m_hct[8].x,m_hct[8].y,m_hct[8].z);
  472.    glVertex3d(m_hct[9].x,m_hct[9].y,m_hct[9].z);
  474.    normal=GetTriNormal(m_hct[3],m_hct[8],m_hct[7]);
  475.    glNormal3d(normal.x,normal.y,normal.z);
  476.    glVertex3d(m_hct[3].x,m_hct[3].y,m_hct[3].z);
  477.    glVertex3d(m_hct[8].x,m_hct[8].y,m_hct[8].z);
  478.    glVertex3d(m_hct[7].x,m_hct[7].y,m_hct[7].z);
  480.    normal=GetTriNormal(m_hct[3],m_hct[7],m_hct[5]);
  481.    glNormal3d(normal.x,normal.y,normal.z);
  482.    glVertex3d(m_hct[3].x,m_hct[3].y,m_hct[3].z);
  483.    glVertex3d(m_hct[7].x,m_hct[7].y,m_hct[7].z);
  484.    glVertex3d(m_hct[5].x,m_hct[5].y,m_hct[5].z);
  485.     
  486.    normal=GetTriNormal(m_hct[3],m_hct[2],m_hct[8]);
  487.    glNormal3d(normal.x,normal.y,normal.z);
  488.    glVertex3d(m_hct[3].x,m_hct[3].y,m_hct[3].z); 
  489.    glVertex3d(m_hct[2].x,m_hct[2].y,m_hct[2].z);
  490.    glVertex3d(m_hct[8].x,m_hct[8].y,m_hct[8].z);
  491.      glEnd();
  492. }//if else
  493. }//for*/
  494. }//if
  495. if(m_Dimension==DO_DRAW2 && m_HCT==DO_HCT)
  496. {
  497. for(i=1;i<=3;i++)
  498. {
  499. HCT(tri,i);
  500. glColor4d(0.0,0.0,1.0,0.5);
  501. glBegin(GL_TRIANGLES);//anticlock
  502.    glVertex3d(m_hct[0].x,m_hct[0].y,0);          
  503.    glVertex3d(m_hct[1].x,m_hct[1].y,0);
  504.    glVertex3d(m_hct[2].x,m_hct[2].y,0); 
  506.    glVertex3d(m_hct[1].x,m_hct[1].y,0); 
  507.    glVertex3d(m_hct[2].x,m_hct[2].y,0); 
  508.    glVertex3d(m_hct[3].x,m_hct[3].y,0); 
  510.    glVertex3d(m_hct[3].x,m_hct[3].y,0); 
  511.    glVertex3d(m_hct[1].x,m_hct[1].y,0); 
  512.    glVertex3d(m_hct[4].x,m_hct[4].y,0); 
  514.    glVertex3d(m_hct[3].x,m_hct[3].y,0); 
  515.    glVertex3d(m_hct[4].x,m_hct[4].y,0); 
  516.    glVertex3d(m_hct[5].x,m_hct[5].y,0); 
  518.    glVertex3d(m_hct[4].x,m_hct[4].y,0); 
  519.    glVertex3d(m_hct[6].x,m_hct[6].y,0); 
  520.         glVertex3d(m_hct[5].x,m_hct[5].y,0); 
  522.    glVertex3d(m_hct[7].x,m_hct[7].y,0);  
  523.         glVertex3d(m_hct[8].x,m_hct[8].y,0); 
  524.    glVertex3d(m_hct[9].x,m_hct[9].y,0); 
  526.    glVertex3d(m_hct[8].x,m_hct[8].y,0); 
  527.    glVertex3d(m_hct[3].x,m_hct[3].y,0); 
  528.    glVertex3d(m_hct[7].x,m_hct[7].y,0); 
  529.  
  530.    glVertex3d(m_hct[3].x,m_hct[3].y,0); 
  531.    glVertex3d(m_hct[5].x,m_hct[5].y,0); 
  532.    glVertex3d(m_hct[7].x,m_hct[7].y,0); 
  533.     
  534.    glVertex3d(m_hct[2].x,m_hct[2].y,0); 
  535.    glVertex3d(m_hct[3].x,m_hct[3].y,0); 
  536.    glVertex3d(m_hct[8].x,m_hct[8].y,0); 
  537.      glEnd();
  538. }
  539. }
  540. }
  542. int CDelaunayView::DelTriList(double x, double y,int p)
  543. {//record the position of triangles which are to be deleted
  544. //'x' and 'y' are the coordinates of the insert point
  545. //'p' is the mark or th insert point in 'm_point'
  547. CPointPos *point=new CPointPos(x,y);
  548. int j=0,i=0,k=0,max=0;
  549. array a;
  550. // GetInitEdges : ransack each trangle to record it's edges
  551. //when the piont belong it's circle
  552.     k=m_pDoc->GetInitEdges(x,y,p);// the returned value 'k/2' is the number of cirecle the point belonged??
  553. if(k==POS_ERROR){
  554. return POS_ERROR;
  555. }
  556. //默认:插入点属于某个外接圆,且属于凸包
  557.           //judge a point belong to the convexity or not
  558.   a=Wher(point);
  559.   /*return Value:a.a
  560. POS_IN=1,a point  belong a triangle or convexity
  561. POS_ON=2,a point belong a circle or on the edge of convexity
  562. POS_OUT=0,a point out of a circle or convexity
  563. a.b , a.c: record the point's mark in m_con */
  564.    max=m_pDoc->m_con.GetSize();
  565.   if(0==a.a && k>=2)//插入点属于某个外接圆,但不属于凸包
  566.   {  //delete an edge "belong to" the inserted polygon,the edge is also 
  567.  //belong  to the border of the convexity
  568.  if (max==3){
  569. int t=m_pDoc->m_edge.GetSize();       
  570. for(j=0;j<t;j++){
  571. if((m_pDoc->m_edge[j]->m_p1==m_pDoc->m_con[a.b]) & 
  572.    (m_pDoc->m_edge[j]->m_p2==m_pDoc->m_con[a.c]))
  573.    {
  574.    m_pDoc->m_edge.RemoveAt(j,1);
  575.    m_pDoc->m_con.InsertAt(a.c,p);
  576.    break;
  577.    }
  578. }
  579.  }
  580.  else
  581.  {
  582.     CWordArray con_index;//save the point's mark can be saw on the con
  583.     for(i=0;i<max-1;i++)
  584. {
  585.    int t=m_pDoc->m_edge.GetSize();       
  586.    for(j=0;j<t;j++){
  587.   if((m_pDoc->m_edge[j]->m_p1==m_pDoc->m_con[i]) & 
  588.    (m_pDoc->m_edge[j]->m_p2==m_pDoc->m_con[i+1]))
  589.   {
  590.   bool y=m_pDoc->DelEdgeOrNot(m_pDoc->m_edge[j]->m_p1,m_pDoc->m_edge[j]->m_p2,p);
  591.                           if(y){
  592.               m_pDoc->m_edge.RemoveAt(j,1);
  593.               con_index.Add(i);
  594.               con_index.Add(i+1);    
  595.   break;
  596.   }
  597.   }
  598.    }
  599. }//end for i
  600.                 int t=m_pDoc->m_edge.GetSize();       
  601.     for(j=0;j<t;j++){       
  602.    if((m_pDoc->m_edge[j]->m_p1==m_pDoc->m_con[max-1]) && 
  603.    (m_pDoc->m_edge[j]->m_p2==m_pDoc->m_con[0]))
  604.    {
  605.     bool y=m_pDoc->DelEdgeOrNot(m_pDoc->m_edge[j]->m_p1,m_pDoc->m_edge[j]->m_p2,p);
  606.                         if(y){
  607.               m_pDoc->m_edge.RemoveAt(j,1);
  608.               con_index.Add(max-1);
  609.               con_index.Add(0);    
  610.   break;
  611. }     
  612.    }
  613. }
  614. t=con_index.GetSize();
  615.     m_pDoc->EditCon(con_index[t-1],con_index[0],p);
  616.              }
  617.   }
  618.   if(k==0){//插入点不属于任意一个外接圆,显然也不属于凸包
  619.   CBorder * m_border;
  620.   m_border=new CBorder(m_pDoc->m_con[a.b],m_pDoc->m_con[a.c]);
  621.   m_pDoc->m_edge.Add(m_border);
  622.   m_pDoc->EditCon(a.c,a.b,p);
  623.   }
  624.    //Delete triangles that have been marked
  625.   m_pDoc->DelTriMarked();
  626.   //add new triangles
  627.       m_pDoc->AddTriangle(p);
  628.   max=m_pDoc->m_index.GetSize();
  629.   m_pDoc->m_index.RemoveAt(0,max);//RemoveAll : no use ???
  630.   m_pDoc->m_edge.RemoveAll();
  631. }
  633. void CDelaunayView::OnDestroy() 
  634. {
  635. HGLRC hrc;
  636. hrc=::wglGetCurrentContext();
  637. ::wglMakeCurrent(NULL,NULL);
  638. if(m_pDC)
  639. delete m_pDC;
  640. CView::OnDestroy();
  641. }
  643. void CDelaunayView::DrawScene()
  644. {
  645. glClearColor(1.0f,1.0f,1.0f,1.0f);//setup backcolor
  646. glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
  647.     //GL_COLOR_BUFFER_BIT 用当前设置的backcolor 第一次清除 color buffer
  648. // GL_DEPTH_BUFFER_BIT清除 depth buffer,以便进行depth test and 消除隐藏面
  649. glClearDepth(1.0f);
  650. glDepthFunc(GL_LESS);
  651. glEnable(GL_DEPTH_TEST);
  652. if(m_bMaterial){
  653. CreateMaterial();
  654. CreateLights();
  655. }
  656. glColorMask(!m_bMaskRed,!m_bMaskGreen,!m_bMaskBlue,GL_TRUE);
  657. glColor3f(m_colorRed,m_colorGreen,m_colorBlue);
  658. glMatrixMode(GL_MODELVIEW);
  659. glLoadIdentity();
  660. glTranslated(m_translate_x,m_translate_y,m_translate_z);
  661. if(m_rotate_angle!=0.0){
  662. glRotated(m_rotate_angle,m_gl_x,m_gl_y,m_gl_z);
  663. }
  664. glScaled(2,2,0);
  665. // glScaled((cx-2*m_ox)/cx,(cy-2*m_oy)/cy,0.0);
  666. if(m_Dimension==DO_DRAW3){
  667. // glPushMatrix();
  668.     gluLookAt(10.0,11.0,12.0,0.0,0.0,0.0,0.0,0.0,1.0);//eye.x,eye.y,eye.z
  670. // glPopMatrix();
  671. }
  672. glBegin(GL_LINE_STRIP);
  673.            glVertex3d(0.0,0.0,0.0);
  674.    glVertex3d(0.0,1.0,0.0);
  675.    glVertex3d(1.0,1.0,0.0);
  676.    glVertex3d(1.0,0.0,0.0);
  677.    glVertex3d(0.0,0.0,0.0);
  678.     glEnd(); 
  679. glPushMatrix();
  680. if(m_Dimension==DO_DRAW2)
  681. DrawPoints();
  682. if(m_Draw_what==DO_WANG && m_Dimension==DO_DRAW3)
  683. m_pDoc->Wang();
  684. else
  685.         DrawTris();
  686. glPopMatrix();
  687. glTranslated(0.5,0.5,0.0);
  689. // glPushMatrix();
  691. // auxSolidSphere(0.5);
  693. // glPopMatrix();
  695. glDisable(GL_LIGHTING);
  696. glDisable(GL_LIGHT0);
  697. glDisable(GL_DEPTH_TEST);
  698. // glFinish();
  699. }
  701. void CDelaunayView::InitOpenGl()
  702. {
  703. PIXELFORMATDESCRIPTOR pfd;
  704.     int n;
  705. HGLRC hrc;
  706. m_pDC=new CClientDC(this);
  707. ASSERT(m_pDC!=NULL);
  708. if(!bSetupPixelFormat()) return;
  709. n=::GetPixelFormat(m_pDC->GetSafeHdc());
  710. ::DescribePixelFormat(m_pDC->GetSafeHdc(),n,sizeof(pfd),&pfd);
  711. hrc=wglCreateContext(m_pDC->GetSafeHdc());
  712. //create a instance of RC: hrc
  713. wglMakeCurrent(m_pDC->GetSafeHdc(),hrc);
  714. //make hrc a current Rendering context of a thread
  715. GetClientRect(&m_oldRect);//get size of ClientRect
  716. }
  718. BOOL CDelaunayView::bSetupPixelFormat()
  719. {
  720. static PIXELFORMATDESCRIPTOR pfd={
  721. sizeof(PIXELFORMATDESCRIPTOR),
  722. 1,
  723. PFD_DRAW_TO_WINDOW |
  724. PFD_SUPPORT_OPENGL |
  725. PFD_DOUBLEBUFFER,
  726. PFD_TYPE_RGBA,
  727. 24,//24位真彩
  728. 0,0,0,0,0,0,
  729. 0,
  730. 0,
  731. 0,
  732. 0,0,0,0,
  733. 32,//32-bit z-buffer,可从物体中消除隐藏面
  734. 0,
  735. 0,
  736. PFD_MAIN_PLANE,//选择主层面
  737. 0,
  738. 0,0,0
  739. };
  740. int pixelformat;
  742. if((pixelformat=ChoosePixelFormat(m_pDC->GetSafeHdc(),&pfd))==0)
  743. {
  744. MessageBox("ChoosePixelFormat failed");
  745. return FALSE;
  746. }
  748. if(SetPixelFormat(m_pDC->GetSafeHdc(),pixelformat,&pfd)==FALSE)
  749. {
  750. MessageBox("SetPixelFormat failed");
  751. return FALSE;
  752. }
  754. return TRUE;
  755. //CClientDC clientdc(this);
  756. //int pf=ChoosePixelFormat(clientdc.m_hDC,&pfd);
  757. //BOOL rt=SetPixelFormat(clientdc.m_hDC,pf,&pfd);
  758. //hglrc=wglCreateContext(clientdc.m_hDC);
  759. }
  763. void CDelaunayView::DrawPoint(int i)
  764. {
  765. double x,y,z;
  766. x=m_pDoc->m_point[i]->m_x;
  767. y=m_pDoc->m_point[i]->m_y;
  768. z=m_pDoc->m_point[i]->m_z;
  769.     glBegin(GL_POINTS);
  770. glVertex3d(m_pDoc->m_point[i]->m_x,m_pDoc->m_point[i]->m_y,0.0);
  771.     glEnd();
  772. }
  774. void CDelaunayView::OnKeyDown(UINT nChar, UINT nRepCnt, UINT nFlags) 
  775. {
  776. switch(nChar)
  777. {
  778. case VK_LEFT:
  779. m_translate_x=m_translate_x-0.2;
  780. break;
  781. case VK_RIGHT:
  782.     m_translate_x=m_translate_x+0.2;
  783. break;
  784. case VK_UP:
  785. m_translate_y=m_translate_y+0.1;
  786.         break;
  787. case VK_DOWN:
  788. m_translate_y=m_translate_y-0.2;
  789. break;
  790. case 88://x
  791. if(m_gl_x==0.0){ 
  792.             m_gl_x=0.1;
  793. }
  794. else{
  795. m_gl_x=0.0;
  796. }
  797. break;
  798. case 89://y
  799. if(m_gl_y==0.0){ 
  800.             m_gl_y=0.1;
  801. }
  802. else{
  803. m_gl_y=0.0;
  804. }
  805. break;
  806. case 90://z
  807. if(m_gl_z==0.0){ 
  808.             m_gl_z=0.1;
  809. }
  810. else{
  811. m_gl_z=0.0;
  812. }
  813. break;
  814. case 100://小键盘的left
  815. m_rotate_angle=m_rotate_angle-10.0;
  816. break;
  817. case 102://小键盘的right
  818. m_rotate_angle=m_rotate_angle+10.0;
  819. break;
  820. case 104://小键盘的up
  821. m_rotate_angle=m_rotate_angle+90.0;
  822. break;
  823. case 98://小键盘的down
  824. m_rotate_angle=m_rotate_angle-90.0;
  825. break;
  826. case 65://a
  827. eye.x=eye.x+1.0;
  828. break;
  829.     case 83://s
  830. eye.y=eye.y+1.0;
  831. break;
  832.     case 68://d
  833. eye.z=eye.z+1.0;
  834. break;
  835. default:
  836. break;
  837. }
  838. CView::OnKeyDown(nChar, nRepCnt, nFlags);
  839.     Invalidate();
  840. }
  842. void CDelaunayView::DrawTris()
  843. {
  844. POSITION POS;
  845. CTriangle* pTriangle;
  846. POS = m_pDoc->m_tri.GetHeadPosition();
  847. if(m_Draw_what==DO_LINE)
  848. glPolygonMode(GL_FRONT_AND_BACK,GL_LINE);
  849.     else{ 
  850.     glPolygonMode(GL_FRONT,GL_FILL);
  851.     glPolygonMode(GL_BACK,GL_FILL);
  852. }
  853. while(POS != NULL ){
  854.   DrawTri(m_pDoc->m_tri.GetAt(POS));
  855.       pTriangle=m_pDoc->m_tri.GetNext(POS);
  856. }   
  857. }
  859. void CDelaunayView::DrawPoints()
  860. {
  861. int max=m_pDoc->m_point.GetSize();
  862. glPointSize(6.0);
  863. for(int i=0;i<max;i++)
  864. {
  865.    glColor4d(double(i)/double(max),1.0,0.0,0.0);
  866.    DrawPoint(i);
  867. }
  868. }
  871. void CDelaunayView::OnCorlorMaterial() 
  872. {
  873.    CLightMaterial dlg;
  874.    dlg.m_colorRed=m_colorRed;
  875.    dlg.m_colorGreen=m_colorGreen;
  876.    dlg.m_colorBlue=m_colorBlue;
  877.    dlg.m_bMaskRed=m_bMaskRed;
  878.    dlg.m_bMaskGreen=m_bMaskGreen;
  879.    dlg.m_bMaskBlue=m_bMaskBlue;
  881.    dlg.m_matAmbientR=m_materialAmb[0];
  882.    dlg.m_matAmbientG=m_materialAmb[1];
  883.    dlg.m_matAmbientB=m_materialAmb[2];
  884.    dlg.m_matDiffuseR=m_materialDif[0];
  885.    dlg.m_matDiffuseG=m_materialDif[1];
  886.    dlg.m_matDiffuseB=m_materialDif[2];
  887.    dlg.m_matSpecularR=m_materialSpe[0];
  888.    dlg.m_matSpecularG=m_materialSpe[1];
  889.    dlg.m_matSpecularB=m_materialSpe[2];
  890.    dlg.m_matEmissionR=m_materialEmi[0];
  891.    dlg.m_matEmissionG=m_materialEmi[1];
  892.    dlg.m_matEmissionB=m_materialEmi[2];
  894.    dlg.m_matShininess=m_matshininess;
  895.    dlg.m_bMaterialOnOff=m_bMaterial;
  897.    int responeDlg=dlg.DoModal();
  898.    if(responeDlg==IDOK){
  899.    m_colorRed=dlg.m_colorRed;
  900.        m_colorGreen=dlg.m_colorGreen;
  901.        m_colorBlue=dlg.m_colorBlue;
  902.        m_bMaskRed=dlg.m_bMaskRed;
  903.        m_bMaskGreen=dlg.m_bMaskGreen;
  904.        m_bMaskBlue=dlg.m_bMaskBlue;
  906.    m_materialAmb[0]=dlg.m_matAmbientR;
  907.        m_materialAmb[1]=dlg.m_matAmbientG;
  908.        m_materialAmb[2]=dlg.m_matAmbientB;
  909.        m_materialDif[0]=dlg.m_matDiffuseR;
  910.        m_materialDif[1]=dlg.m_matDiffuseG;
  911.        m_materialDif[2]=dlg.m_matDiffuseB;
  912.        m_materialSpe[0]=dlg.m_matSpecularR;
  913.        m_materialSpe[1]=dlg.m_matSpecularG;
  914.    m_materialSpe[2]=dlg.m_matSpecularB;
  915.        m_materialEmi[0]=dlg.m_matEmissionR;
  916.        m_materialEmi[1]=dlg.m_matEmissionG;
  917.        m_materialEmi[2]=dlg.m_matEmissionB;
  919.        m_matshininess=dlg.m_matShininess;
  920.        m_bMaterial=dlg.m_bMaterialOnOff;
  921.    Invalidate();
  922.    }
  923. }
  925. void CDelaunayView::OnUpdateCorlorMaterial(CCmdUI* pCmdUI) 
  926. {
  927. // TODO: Add your command update UI handler code here
  929. }
  931. void CDelaunayView::CreateLights()
  932. {
  933. glLightfv(GL_LIGHT0,GL_AMBIENT,m_lightAmb);
  934. glLightfv(GL_LIGHT0,GL_DIFFUSE,m_lightDif);
  935. glLightfv(GL_LIGHT0,GL_SPECULAR,m_lightSpe);
  936. glLightfv(GL_LIGHT0,GL_POSITION,m_lightPos);
  938. glEnable(GL_LIGHTING);
  939. glEnable(GL_LIGHT0);
  940. }
  942. void CDelaunayView::OnUpdateLight(CCmdUI* pCmdUI) 
  943. {
  944. // TODO: Add your command update UI handler code here
  946. }
  948. void CDelaunayView::OnLight() 
  949. {
  950. CLighting lightDlg;
  951. lightDlg.m_lightAmbientR=m_lightAmb[0];
  952. lightDlg.m_lightAmbientG=m_lightAmb[1];
  953. lightDlg.m_lightAmbientB=m_lightAmb[2];
  954. lightDlg.m_lightDiffuseR=m_lightDif[0];
  955. lightDlg.m_lightDiffuseG=m_lightDif[1];
  956. lightDlg.m_lightDiffuseB=m_lightDif[2];
  957. lightDlg.m_lightSpecularR=m_lightSpe[0];
  958. lightDlg.m_lightSpecularG=m_lightSpe[1];
  959. lightDlg.m_lightSpecularB=m_lightSpe[2];
  961. lightDlg.m_lightPositionX=m_lightPos[0];
  962. lightDlg.m_lightPositionY=m_lightPos[1];
  963. lightDlg.m_lightPositionZ=m_lightPos[2];
  965. int responeDlg=lightDlg.DoModal();
  966. if(responeDlg==IDOK){
  967. m_lightAmb[0]=lightDlg.m_lightAmbientR;
  968. m_lightAmb[1]=lightDlg.m_lightAmbientG;
  969. m_lightAmb[2]=lightDlg.m_lightAmbientB;
  970.         m_lightDif[0]=lightDlg.m_lightDiffuseR;
  971. m_lightDif[1]=lightDlg.m_lightDiffuseG;
  972. m_lightDif[2]=lightDlg.m_lightDiffuseB;
  973.         m_lightSpe[0]=lightDlg.m_lightSpecularR;
  974. m_lightSpe[1]=lightDlg.m_lightSpecularG;
  975. m_lightSpe[2]=lightDlg.m_lightSpecularB;
  977. m_lightPos[0]=lightDlg.m_lightPositionX;
  978. m_lightPos[1]=lightDlg.m_lightPositionY;
  979. m_lightPos[2]=lightDlg.m_lightPositionZ;
  981. Invalidate();
  982. }
  983. }
  985. void CDelaunayView::CreateMaterial()
  986. {
  987. glMaterialfv(GL_FRONT,GL_AMBIENT,m_materialAmb);
  988. glMaterialfv(GL_FRONT,GL_DIFFUSE,m_materialDif);
  989. glMaterialfv(GL_FRONT,GL_SPECULAR,m_materialSpe);
  990. glMaterialfv(GL_FRONT,GL_EMISSION,m_materialEmi);
  991. glMaterialf(GL_FRONT,GL_SHININESS,m_matshininess);
  992. }
  994. void CDelaunayView::OnButtonTwo() 
  995. {
  996. m_Dimension=DO_DRAW2;
  997. Invalidate();
  998. }
  1000. void CDelaunayView::OnUpdateButtonTwo(CCmdUI* pCmdUI) 
  1001. {
  1002. if(m_Dimension==DO_DRAW2)
  1003. pCmdUI->SetCheck(1);
  1004. else 
  1005. pCmdUI->SetCheck(0);
  1006. }
  1008. void CDelaunayView::OnButtonThree() 
  1009. {
  1010. m_Dimension=DO_DRAW3;
  1011. Invalidate();
  1012. }
  1014. void CDelaunayView::OnUpdateButtonThree(CCmdUI* pCmdUI) 
  1015. {
  1016.     if(m_Dimension==DO_DRAW3)
  1017. pCmdUI->SetCheck(1);
  1018. else 
  1019. pCmdUI->SetCheck(0);
  1020. }
  1022. void CDelaunayView::HCT(CTriangle *temp, int p)
  1023. {//每次处理1/3三角形
  1024. m_hct[6].x=temp->m_x;//重心
  1025. m_hct[6].y=temp->m_y;
  1026. m_hct[6].z=temp->o;
  1027. if(p==3){
  1028.     m_hct[0].x=m_pDoc->m_point[temp->m_p1]->m_x;//f1
  1029.     m_hct[0].y=m_pDoc->m_point[temp->m_p1]->m_y;
  1030.     m_hct[0].z=m_pDoc->m_point[temp->m_p1]->m_z;
  1032. m_hct[9].x=m_pDoc->m_point[temp->m_p2]->m_x;//f2
  1033.     m_hct[9].y=m_pDoc->m_point[temp->m_p2]->m_y;
  1034.     m_hct[9].z=m_pDoc->m_point[temp->m_p2]->m_z;
  1036.     m_hct[1]=PointOnEdge_3(m_hct[0].x,m_hct[0].y,m_hct[6].x,m_hct[6].y);//c1
  1037.     m_hct[1].z=temp->c1;//c1
  1038.     m_hct[2]=PointOnEdge_3(m_hct[0].x,m_hct[0].y,m_hct[9].x,m_hct[9].y);//d31
  1039.         m_hct[2].z=temp->d31;//d31
  1040.         m_hct[4]=PointOnEdge_3(m_hct[6].x,m_hct[6].y,m_hct[0].x,m_hct[0].y);//b1
  1041.         m_hct[4].z=temp->b1;//b1
  1042.         m_hct[3]=BaryCenter(m_hct[0].x,m_hct[0].y,0,m_hct[9].x,m_hct[9].y,0,
  1043.  m_hct[6].x,m_hct[6].y,0);//e3
  1044.         m_hct[3].z=temp->e3;//e3
  1045.         m_hct[5]=PointOnEdge_3(m_hct[6].x,m_hct[6].y,m_hct[9].x,m_hct[9].y);//b2
  1046.         m_hct[5].z=temp->b2;//b2
  1047. m_hct[7]=PointOnEdge_3(m_hct[9].x,m_hct[9].y,m_hct[6].x,m_hct[6].y);//c2
  1048.         m_hct[7].z=temp->c2;//c2
  1049. m_hct[8]=PointOnEdge_3(m_hct[9].x,m_hct[9].y,m_hct[0].x,m_hct[0].y);//d32
  1050.         m_hct[8].z=temp->d32;//d32
  1051. }
  1053. if(p==2){
  1054.     m_hct[0].x=m_pDoc->m_point[temp->m_p1]->m_x;//f1
  1055.     m_hct[0].y=m_pDoc->m_point[temp->m_p1]->m_y;
  1056.     m_hct[0].z=m_pDoc->m_point[temp->m_p1]->m_z;
  1057. m_hct[9].x=m_pDoc->m_point[temp->m_p3]->m_x;//f3
  1058.     m_hct[9].y=m_pDoc->m_point[temp->m_p3]->m_y;
  1059.     m_hct[9].z=m_pDoc->m_point[temp->m_p3]->m_z;
  1061.     m_hct[1]=PointOnEdge_3(m_hct[0].x,m_hct[0].y,m_hct[6].x,m_hct[6].y);//c1
  1062.     m_hct[1].z=temp->c1;//c1
  1063.     m_hct[2]=PointOnEdge_3(m_hct[0].x,m_hct[0].y,m_hct[9].x,m_hct[9].y);//d21
  1064.         m_hct[2].z=temp->d21;
  1065.        
  1066. m_hct[4]=PointOnEdge_3(m_hct[6].x,m_hct[6].y,m_hct[0].x,m_hct[0].y);//b1
  1067.         m_hct[4].z=temp->b1;//b1
  1068.         m_hct[3]=BaryCenter(m_hct[0].x,m_hct[0].y,0,m_hct[6].x,m_hct[6].y,0,
  1069.  m_hct[9].x,m_hct[9].y,0);//e2
  1070.         m_hct[3].z=temp->e2;//e2
  1071.        
  1072. m_hct[5]=PointOnEdge_3(m_hct[6].x,m_hct[6].y,m_hct[9].x,m_hct[9].y);//b3
  1073.         m_hct[5].z=temp->b3;//b3
  1074. m_hct[7]=PointOnEdge_3(m_hct[9].x,m_hct[9].y,m_hct[6].x,m_hct[6].y);//c3
  1075.         m_hct[7].z=temp->c3;//c3
  1076. m_hct[8]=PointOnEdge_3(m_hct[9].x,m_hct[9].y,m_hct[0].x,m_hct[0].y);//d23
  1077.         m_hct[8].z=temp->d23;//d23
  1078. }
  1079. if(p==1){
  1080.     m_hct[0].x=m_pDoc->m_point[temp->m_p2]->m_x;//f2
  1081.     m_hct[0].y=m_pDoc->m_point[temp->m_p2]->m_y;
  1082.     m_hct[0].z=m_pDoc->m_point[temp->m_p2]->m_z;
  1083. m_hct[9].x=m_pDoc->m_point[temp->m_p3]->m_x;//f3
  1084.     m_hct[9].y=m_pDoc->m_point[temp->m_p3]->m_y;
  1085.     m_hct[9].z=m_pDoc->m_point[temp->m_p3]->m_z;
  1087.     m_hct[1]=PointOnEdge_3(m_hct[0].x,m_hct[0].y,m_hct[6].x,m_hct[6].y);//c2
  1088.     m_hct[1].z=temp->c2;//c2
  1089.     m_hct[2]=PointOnEdge_3(m_hct[0].x,m_hct[0].y,m_hct[9].x,m_hct[9].y);//d12
  1090.         m_hct[2].z=temp->d12;//d12
  1091.        
  1092. m_hct[4]=PointOnEdge_3(m_hct[6].x,m_hct[6].y,m_hct[0].x,m_hct[0].y);//b2
  1093.         m_hct[4].z=temp->b2;//b2
  1094.         m_hct[3]=BaryCenter(m_hct[0].x,m_hct[0].y,0,m_hct[9].x,m_hct[9].y,0,
  1095.  m_hct[6].x,m_hct[6].y,0);//e1
  1096.         m_hct[3].z=temp->e1;//e1
  1097.         
  1098. m_hct[5]=PointOnEdge_3(m_hct[6].x,m_hct[6].y,m_hct[9].x,m_hct[9].y);//b3
  1099.         m_hct[5].z=temp->b3;//b3
  1100. m_hct[7]=PointOnEdge_3(m_hct[9].x,m_hct[9].y,m_hct[6].x,m_hct[6].y);//c3
  1101.         m_hct[7].z=temp->c3;//c3
  1102. m_hct[8]=PointOnEdge_3(m_hct[9].x,m_hct[9].y,m_hct[0].x,m_hct[0].y);//d13
  1103.         m_hct[8].z=temp->d13;//d13
  1104. }
  1105. b[0][0][3]=m_hct[6].z;b[0][3][0]=m_hct[9].z;b[3][0][0]=m_hct[0].z;
  1107. b[0][1][2]=m_hct[5].z;b[0][2][1]=m_hct[7].z;
  1109. b[1][0][2]=m_hct[4].z;b[2][0][1]=m_hct[1].z;
  1111. b[2][1][0]=m_hct[2].z;b[1][2][0]=m_hct[8].z;b[1][1][1]=m_hct[3].z;
  1112.     
  1113. for(int i=1;i<9;i++){
  1114. if(i!=6){
  1115. m_hct[i].z=Bezier(m_hct[i]);//得到所有细分点的z值
  1116. }
  1117. }
  1118. }
  1120. POI CDelaunayView::BaryCenter(double x1, double y1,double z1,double x2, double y2,double z2,double x3, double y3,double z3)
  1121. {//平面上三角形重心
  1122. POI temp;
  1123.     temp.x=(x1+x2+x3)/double(3);
  1124. temp.y=(y1+y2+y3)/double(3);
  1125. temp.z=(z1+z2+z3)/double(3);
  1126.     return temp;
  1127. }
  1129. POI CDelaunayView::PointOnEdge_3(double x1, double y1, double x2, double y2)
  1130. {
  1131. POI temp;
  1132. temp.x=x1+(x2-x1)/double(3);
  1133. temp.y=y1+(y2-y1)/double(3);
  1134. return temp;
  1135. }
  1137. void CDelaunayView::OnButtonHct() 
  1138. {
  1139. if(m_HCT==DO_HCT)
  1140.     m_HCT=-1;
  1141. else
  1142. m_HCT=DO_HCT;
  1143. }
  1145. void CDelaunayView::OnUpdateButtonHct(CCmdUI* pCmdUI) 
  1146. {
  1147.     
  1148. }
  1150. void CDelaunayView::OnButtonFill() 
  1151. {
  1152.     m_Draw_what=DO_FILL;
  1153. }
  1155. void CDelaunayView::OnUpdateButtonFill(CCmdUI* pCmdUI) 
  1156. {
  1157. if(m_Draw_what==DO_FILL)
  1158. pCmdUI->SetCheck(1);
  1159. else 
  1160. pCmdUI->SetCheck(0);
  1162. }
  1164. void CDelaunayView::OnButtonLine() 
  1165. {
  1166. m_Draw_what=DO_LINE;
  1167. }
  1169. void CDelaunayView::OnUpdateButtonLine(CCmdUI* pCmdUI) 
  1170. {
  1171.     if(m_Draw_what==DO_LINE)
  1172. pCmdUI->SetCheck(1);
  1173. else 
  1174. pCmdUI->SetCheck(0);
  1175. }
  1177. POI CDelaunayView::GetTriNormal(POI p1, POI p2, POI p3)
  1178. { //得到三角片的单位法向
  1179. POI vector1,vector2,Normal;
  1180.     /*POI p;
  1181. double s;
  1182. s=S(p1,p2,p3);
  1183. if(s<0){
  1184. p=p2;p2=p3;p3=p;
  1185. }*/
  1186. vector1.x=p2.x-p1.x;
  1187. vector1.y=p2.y-p1.y;
  1188. vector1.z=p2.z-p1.z;
  1189. vector2.x=p3.x-p1.x;
  1190. vector2.y=p3.y-p1.y;
  1191. vector2.z=p3.z-p1.z;
  1192.     //get the normal n=a*b(外积)
  1193. Normal=m_pDoc->VectorProduct(vector1.x,vector1.y,vector1.z,vector2.x,vector2.y,vector2.z);
  1194. //单位化 normal
  1195. Normal=m_pDoc->Unitization(Normal);
  1196.     return Normal;
  1197. }
  1200. double CDelaunayView::S(POI p1, POI p2, POI p3)
  1201. { //求三角形面积,以右下角为原点时,s>0为逆时针(左转),
  1202.    //s=0为三点重合
  1203.    double s;
  1204.    s=p1.x*p2.y+p2.x*p3.y+p1.y*p3.x-p2.y*p3.x-p1.y*p2.x-p1.x*p3.y;
  1205.    s=s/2.0;
  1206.    return s;
  1207. }
  1209. int CDelaunayView::Factorial(int n)
  1210. {//阶乘
  1211. int x=n;
  1212. if(n==0 ||n==1) 
  1213. return 1;
  1214. else{
  1215. for(int i=1;i<n;i++){
  1216. x=x*(n-i);
  1217. }
  1218. }
  1219.     return x;
  1220. }
  1222. double CDelaunayView::Power(double a, int e)
  1223. {
  1224. double x=1.0;
  1225. if(e==0){
  1226. return 1.0;
  1227. }
  1228. else{
  1229. if(a==0.0){
  1230. return 0.0;
  1231. }
  1232.         for(int i=1;i<=e;i++){
  1233. x=x*a;
  1234. }
  1235. }
  1236.     return x;
  1237. }
  1239. double CDelaunayView::Bezier(POI p)
  1240. {//s1,s2,s3 : 面积坐标
  1241. double s1,s2,s3,s;
  1242. s=S(m_hct[0],m_hct[9],m_hct[6]);
  1243. s1=S(p,m_hct[9],m_hct[6])/s;
  1244. s2=S(m_hct[0],p,m_hct[6])/s;
  1245. s3=S(m_hct[0],m_hct[9],p)/s;
  1246. if(fabs(s1)<0.000000000001) s1=0;
  1247. if(fabs(s2)<0.000000000001) s2=0;
  1248. if(fabs(s3)<0.000000000001) s3=0;
  1249. int i,j,k;
  1250. double B=0.0;
  1251. for(i=0;i<4;i++){
  1252. for(j=0;j<=(3-i);j++){
  1253. for(k=0;k<=(3-i-j);k++){
  1254. if((i+j+k)==3){
  1255. B=B+b[i][j][k]*double(6)/(Factorial(i)*Factorial(j)*Factorial(k))*
  1256. Power(s1,i)*Power(s2,j)*Power(s3,k);
  1257. }
  1258. }
  1259. }
  1260. }
  1261. return B;
  1262. }
  1264. void CDelaunayView::OnButtonWnag() 
  1265. {
  1266.     m_Draw_what=DO_WANG;
  1267. }