基于半边数据结构（翼边数据结构）的Eular操作回实现扫成 通过OpenGL进行CAD实体建模 经测试可直接运行

基于半边数据结构（翼边数据结构）的Eular操作来实现扫成 通过OpenGL进行CAD实体建模 经测试可直接运行

话说,今天重装电脑程序,把VS2010重装以后,运行了几个之前写的小程序测试了一下.不经意间,我发现2年前自己参考着资料写的一个Sweep扫成程序,当时程序能够跑,但是和翼边数据结构的含义有所不同,当时费了好大的劲也没弄明白其中一个指针的处理怎么会是那个样子的,现在时间过去两年了,自己竟然还是没有抽出时间来将其弄明白,真是惭愧.也许我当时坚持查阅资料仔细研究就不会是这个样子了,而且有些东西一旦放下，再抽出整块的时间去理解它是不现实的。现在2年过去了,如果不是重装电脑，自己早忘了有这么一回事情。这些年自己在技术方面没有多少成长不怪其他，一切源于自己没有坚持钻研和探索的精神.所以遇到挫折了不要一味的怪罪环境,最主要的还是自己不够努力.下面把程序的结果贴出来,后面还要仔细再研究一下.

下面是我实现程序的例子，第一幅图是一个线框模型，第二幅图是另一个实体模型。

这个程序是在一门选修课（忘了叫什么名字了，大概叫CAD实体建模之类的）上老师布置的大程作业，当时为了做出来费了好大的力气，现在看来这只不过是个要求通过欧拉操作来实现扫成的小程序，而且算法的思想很简单，只是需要用到实体建模中的名为翼边数据结构（又叫，半边数据结构）的东西，和一个欧拉公式。我看过网络上的一些同样的介绍半边数据结构的东西，他们写的也算比较详细，但是有个很致命的问题，那就是他们只写出了核心代码，没有具体实现的环境。对于一个编程初学者而言，你光给出核心代码，他们是没有能力把这个程序给搭建起来的，因为程序不仅涉及到核心数据结构和算法，通过MFC实现鼠标点击和画图也够一般初学者头疼脑大的了。

为了方便初学者学习，也为了分享自己的一些心得，我下面不仅给出了核心代码，而且给出了整个实现程序的源代码。想学习实体建模的同学能够通过这个程序学到很多东西，对于那些想要交图形学的代码课程作业的同学来说，你也可以直接把这个交上去了吆，因为这个代码我已经测试过了，没有任何问题，可以直接运行。当然，最好你是为了学习需要，而不是为了直接down下来应付老师，而且前提是你要部署好OpenGl的库。具体的下载链接在下面。

如果运行过程中有什么问题，欢迎留言讨论。如果大家觉得此文对您有帮助，麻烦你“顶”一下，予人玫瑰，手有余香~~

运行的环境和平台是：
Windows 7， visual studio 2010

刚装了win7和VS2010。配置OpenGL的方法如下。

glut下载地址:

http://download.****.net/detail/daringpig/4375605

glut.h ---> C:\Program Files (x86)\Microsoft SDKs\Windows\v7.0A\Include\gl

glut.dll,glut32.dll ---> C:\Windows\SysWOW64 (windows7 64位操作系统)

---> C:\Windows\System32 (windows7 32位操作系统)

glut.lib,glut32.lib ---> C:\Program Files (x86)\Microsoft Visual Studio 10.0\VC\lib

核心数据结构采用半边数据结构实现，具体的实现功能是：

通过绘制一个基本的面，然后通过该面通过指定的向量和距离就可以扫成为一个体。

具体的核心代码如下：

CSolid* CEulerDoc::mvfs(GLdouble *point) { CSolid *newSolid; CFace *newFace; CLoop *newLoop; CVertex *newVert; newSolid = new CSolid; newFace = new CFace; newLoop = new CLoop; newVert = new CVertex; // Set the coordinate of the new vertex newVert->SetCoord(point); // Add the new face to the new solid newSolid->AddFace(newFace); // Add the new vertex to the new solid newSolid->AddVertex(newVert); // Add the new loop to the new face newFace->AddLoop(newLoop); return newSolid; } CHalfEdge* CEulerDoc::mev(GLdouble *point1, GLdouble *point2, CLoop* lp) { CSolid *pSolid = lp->host_f->host_s; CVertex *v1, *newVert; CEdge *newEdge; CHalfEdge *he1, *he2; newVert = new CVertex; newEdge = new CEdge; he1 = new CHalfEdge; he2 = new CHalfEdge; // Find the start vertex of the edge v1 = pSolid->FindVertex(point1); if(!v1) { AfxMessageBox("Can't find the point1!"); return NULL; } // Set the coordinate of the new vertex newVert->SetCoord(point2); // Set v1 as the start vertex of halfedge1 he1->SetVertex(v1); // Set v2 as the start vertex of halfedge2 he2->SetVertex(newVert); // Add two new halfedges to the loop lp->AddHalfEdge(he1, he2); // Set the new edge to new halfedges he1->edge = he2->edge = newEdge; newEdge->he1 = he1; newEdge->he2 = he2; // Add the new edge to the solid pSolid->AddEdge(newEdge); // Add the new vertex to the solid pSolid->AddVertex(newVert); return he1; } CLoop* CEulerDoc::mef(GLdouble *point1, GLdouble *point2, CLoop* lp) { CSolid *pSolid = lp->host_f->host_s; CVertex *v1, *v2; // Find the two given vertex v1 = pSolid->FindVertex(point1); if(!v1) { AfxMessageBox("Can't find the point1!"); return NULL; } v2 = pSolid->FindVertex(point2); if(!v2) { AfxMessageBox("Can't find the point2!"); return NULL; } // two vertex must in the same loop if(!lp->IsVertexIn(v1) && !lp->IsVertexIn(v2)) { AfxMessageBox("该两点不位于同一个loop中"); return NULL; } CFace *newFace; CLoop *newLoop; CEdge *newEdge; CHalfEdge *he1, *he2, *temphe1, *temphe2; newFace = new CFace; newLoop = new CLoop; newEdge = new CEdge; he1 = new CHalfEdge; he2 = new CHalfEdge; // Find two halfedges start with two vertexs temphe1 = lp->FindHostHalfEdge(v1); temphe2 = lp->FindHostHalfEdge(v2); // Change halfedges' host loop to new loop temphe1->prev->next = NULL; newLoop->AddHalfEdge(temphe1, NULL); while(temphe1 != temphe2) temphe1 = temphe1->next; temphe1->prev->next = he2; he2->prev = temphe1->prev; he2->next = newLoop->GetHalfEdgeHead(); newLoop->GetHalfEdgeHead()->prev = he2; temphe1->prev = NULL; // Add halfedge start with vertex one to the old loop % close this loop lp->AddHalfEdge(he1, NULL); he1->next = temphe2; temphe2->prev = he1; // Set two halfedges' start vertex and adjacent he1->SetVertex(v1); he2->SetVertex(v2); he1->adj = he2; he2->adj = he1; // Set the new edge and add to the solid newEdge->he1 = he1; newEdge->he2 = he2; pSolid->AddEdge(newEdge); // Add new face to the solid pSolid->AddFace(newFace); // Add new loop to the new face newFace->AddLoop(newLoop); return newLoop; } CLoop* CEulerDoc::kemr(GLdouble *point1, GLdouble *point2, CLoop* lp) { CSolid *pSolid = lp->host_f->host_s; CVertex *v1, *v2; // Find the two given vertex v1 = pSolid->FindVertex(point1); if(!v1) { AfxMessageBox("Can't find the point1!"); return NULL; } v2 = pSolid->FindVertex(point2); if(!v2) { AfxMessageBox("Can't find the point2!"); return NULL; } // two vertex must in the same loop if(!lp->IsVertexIn(v1) && !lp->IsVertexIn(v2)) return NULL; CLoop *newLoop; CHalfEdge *he1, *he2; CEdge *edge; newLoop = new CLoop; // Find two halfedges start with two vertexs he1 = lp->GetHalfEdgeHead(); while(he1) { if(he1->GetVertex() == v1 && he1->next->GetVertex() == v2) break; he1 = he1->next; } he2 = he1->adj; // Make a newloop and delete two halfedges newLoop->AddHalfEdge(he1->next, NULL); he2->prev->next = he1->next; he1->next->prev = he2->prev; he1->prev->next = he2->next; he2->next->prev = he1->prev; he1->next = he1->prev = he2->next = he2->prev = NULL; delete he1; delete he2; // Find the edge and delete it edge = pSolid->GetEdgeHead(); while(edge) { if(edge->he1 == he1) break; edge = edge->next; } edge->prev->next = edge->next; edge->next->prev = edge->prev; edge->next = edge->prev = NULL; // Add new loop to the face lp->host_f->AddLoop(newLoop); return newLoop; } void CEulerDoc::kfmrh(CLoop* outlp, CLoop* lp) { if(!outlp || !lp) { AfxMessageBox("The loop is NULL!"); return; } CSolid *pSolid = lp->host_f->host_s; // Add the loop to the face outlp->host_f->AddLoop(lp); // Get the face need to be killed CFace* pFace = pSolid->GetFaceEnd(); // Delete it pFace->prev->next = NULL; delete pFace; } void CEulerDoc::sweep(CFace* pFace, GLdouble dx, GLdouble dy, GLdouble dz) { CFace *pEnd; CHalfEdge *pHe, *pHead; CLoop *pLoop, *newLoop, *pOutLoop; GLdouble *pPoint; GLdouble point1[3], point2[3], first[3], last[3]; BOOL bOut = TRUE; // Show that if is the outloop of the top face // Remember the last inside loop's face pEnd = pFace->host_s->GetFaceEnd(); // Start with the second face, because the first face is for buttom pFace = pFace->next; while(pFace) { // Get the first point of the loop newLoop = pLoop = pFace->GetLoopHead(); pHe = pHead = pLoop->GetHalfEdgeHead(); pPoint = pHe->GetVertex()->vcoord; point1[0] = pPoint[0]; point1[1] = pPoint[1]; point1[2] = pPoint[2]; // first[] is used for close the top face // last[] is used for side face first[0] = last[0] = point2[0] = point1[0] + dx; first[1] = last[1] = point2[1] = point1[1] + dy; first[2] = last[2] = point2[2] = point1[2] + dz; // Make the new edge mev(point1, point2, pLoop); // Goto next halfedge pHe = pHe->next; while(pHe->GetVertex() != pHead->GetVertex()) { // Get the point pPoint = pHe->GetVertex()->vcoord; point1[0] = pPoint[0]; point1[1] = pPoint[1]; point1[2] = pPoint[2]; point2[0] = point1[0] + dx; point2[1] = point1[1] + dy; point2[2] = point1[2] + dz; // Make the new edge mev(point1, point2, newLoop); // Make a new side face newLoop = mef(point2, last, newLoop); // Remember the lastest point last[0] = point2[0]; last[1] = point2[1]; last[2] = point2[2]; pHe = pHe->next; } // Close the top face newLoop = mef(first, last, newLoop); if(bOut) { pOutLoop = newLoop; bOut = FALSE; } else kfmrh(pOutLoop, newLoop); if(pFace == pEnd) break; pFace = pFace->next; } } CSolid* CEulerDoc::MakeSolid(GLdouble point[][3], GLint* pointnum, GLint loopnum, GLdouble* vector, GLdouble len) { CSolid *pSolid, *pTemp; CLoop *pLoop, *pHead; int i, j, pos; // Create the solid pSolid = mvfs(point[1]); pHead = pLoop = pSolid->GetFaceHead()->GetLoopHead(); // Create the out loop for(i = 1;i < pointnum[1]-1;i ++) { mev(point[i], point[i+1], pLoop); } mef(point[i], point[1], pLoop); pos = i+1; // Create the inside loops for(i = 2;i < loopnum; i ++) { mev(point[1], point[pos+1], pHead); for(j = 1;j < pointnum[i]-1;j ++) { mev(point[pos+j], point[pos+j+1], pHead); } mef(point[pos+j], point[pos+1], pHead); kemr(point[1], point[pos+1], pHead); pos = pos+j+1; } sweep(pSolid->GetFaceHead(), vector[0] * len, vector[1] * len, vector[2] * len); // Insert the new solid if(!solid) solid = pSolid; else { pTemp = solid; while(pTemp->next) pTemp = pTemp->next; pTemp->next = pSolid; } return pSolid; } void CEulerView::RenderScene(void) { CEulerDoc *pDoc = GetDocument(); CSolid *pSolid; CFace *pFace; CLoop *pLoop; CHalfEdge *pHalfEdge, *pHeHead; GLint i, j, pos = 0; GLfloat draw_ambient[] = { 1.0, 1.0, 1.0}, mat_ambient[] = { 0.5, 0.5, 0.5 }, mat_specular[] = { 1.0, 1.0, 1.0, 1.0 }, shininess[] = { 50.0 }; GLfloat position1[] = { 10.0, 10.0, 10.0, 0.0 }; glLoadIdentity(); glTranslatef(m_xPos, m_yPos, m_zPos); glRotatef(m_xAngle, 1.0f,0.0f,0.0f); glRotatef(m_yAngle, 0.0f,1.0f,0.0f); if(m_bCoord) RenderCoordinate(1); glLightfv(GL_LIGHT0, GL_POSITION, position1); ////////////////////////////////////////////////////////////////////////// // Draw glMaterialfv(GL_FRONT, GL_AMBIENT_AND_DIFFUSE, draw_ambient); for(i = 1; i <= m_iLoopNum ;i ++) { pos += m_iLoopPointNum[i-1]; glBegin(GL_LINE_STRIP); for(j = 1;j <= m_iLoopPointNum[i]; j++) glVertex3d(m_dPoint[pos + j][0], m_dPoint[pos + j][1], m_dPoint[pos + j][2]); glEnd(); } ////////////////////////////////////////////////////////////////////////// pSolid = pDoc->solid; if(m_bWire) { glMaterialfv(GL_FRONT, GL_AMBIENT_AND_DIFFUSE, draw_ambient); glLineWidth(2.0); glBegin(GL_LINES); while(pSolid) { pFace = pSolid->GetFaceHead(); while(pFace) { pLoop = pFace->GetLoopHead(); while(pLoop) { pHalfEdge = pHeHead = pLoop->GetHalfEdgeHead(); while(pHalfEdge) { glVertex3dv(pHalfEdge->GetVertex()->vcoord); glVertex3dv(pHalfEdge->next->GetVertex()->vcoord); pHalfEdge = pHalfEdge->next; if(pHalfEdge == pHeHead) break; } pLoop = pLoop->next; } pFace = pFace->next; } pSolid = pSolid->next; } glEnd(); } else { glMaterialfv(GL_FRONT, GL_AMBIENT_AND_DIFFUSE, mat_ambient); glMaterialfv(GL_FRONT, GL_SPECULAR, mat_specular); glMaterialfv(GL_FRONT, GL_SHININESS, shininess); glEnable(GL_NORMALIZE); while(pSolid) { pFace = pSolid->GetFaceHead(); while(pFace) { gluTessBeginPolygon(tessobj, NULL); glNormal3dv(pFace->feq); pLoop = pFace->GetLoopHead(); while(pLoop) { gluTessBeginContour(tessobj); pHalfEdge = pHeHead = pLoop->GetHalfEdgeHead(); while(pHalfEdge) { gluTessVertex(tessobj, pHalfEdge->GetVertex()->vcoord, pHalfEdge->GetVertex()->vcoord); pHalfEdge = pHalfEdge->next; if(pHalfEdge == pHeHead) break; } gluTessEndContour(tessobj); pLoop = pLoop->next; } gluTessEndPolygon(tessobj); pFace = pFace->next; } pSolid = pSolid->next; } } }

OpenGL库和Eular的源代码下载地址为：

http://download.****.net/detail/daringpig/4375605
如果运行过程中有什么问题，欢迎留言讨论。

http://download.****.net/detail/daringpig/4375599

如果大家觉得此文对您有帮助，麻烦你“顶”一下，谢谢

1楼daringpig3天前 15:10

mm