OpenGL实现3D自由变形

标签: opengl变形
158人阅读 评论(0) 收藏 举报
分类:

笔者介绍:姜雪伟,IT公司技术合伙人,IT高级讲师,CSDN社区专家,特邀编辑,畅销书作者,已出版书籍:《手把手教你架构3D游戏引擎》电子工业出版社和《Unity3D实战核心技术详解》电子工业出版社等。

CSDN视频网址:http://edu.csdn.net/lecturer/144

以前有朋友问我,关于变形动画是如何实现的,实现方式主要有两种,一种是通过美术人员利用Max工具自己调出动画,这种调出的动画太僵硬而且不能根据用户的要求随意变形,只能按照预先调好的动画变形,这种做法可以排除了。第二种做法是通过程序实现自由的动画变换,在这里给读者介绍利用OpenGL实现的变形动画,目的是把动画的变换原理给读者介绍清楚。

在介绍3D变形的之前,读者首先要掌握的是3D的固定流水线, 现在做游戏的开发者对于什么是3D固定流水线一无所知,市面上使用的大部分引擎都封装的非常好,开发者也沦为了只会使用工具,对于一些常识一概不知。在游戏公司做的主要事情就是写写逻辑,调用调用接口,如果需求超出接口的范畴基本上就是一筹莫展,或者说引擎的功能满足不了。面对这种情况,作为开发者必须把3D的基本知识掌握了,这样遇到3D的任何问题都不会感觉难下手,至少知道解决问题的思路。如果读者对于3D固定流水线不清楚可以自己去查阅,我也写过关于3D固定流水线的博客,在这里就不做介绍了,简单的一句话表示3D固定流水线就是将3D模型在2D屏幕上显示的过程,这个过程就是对于矩阵的换算,对应3D固定流水线是3D可编程流水线,不清楚的读者自行查阅一下,这个都必须要掌握的。

接下来介绍实现变形的基本思路:

       第一步,准备使用的模型,模型格式很多种,这个可以使用网上的,也可以自己编写插件实现自定义的模型。本项目对应的模型是后缀名为m的自定义的模型格式,模型格式内容如下所示:


上图显示的只是一部分,为了能让读者看清楚,并没有对模型进行加密,为了计算方便,这里只给出了Vertex顶点和Face面的数据。

第二步,定义模型的数据结构体

#ifndef halfedge_h
#define halfedge_h

#include <cstdlib>

using namespace std;

// Halfedge structures
struct HE_edge;
struct HE_vert;
struct HE_face;
struct HE_line;

struct HE_edge
{
	int edgeID;
	HE_vert* vertS;	// Vertex at the start of halfedge
	HE_vert* vertE;	// Vertex at the end of halfedge
	HE_edge* pair;	// Oppositely oriented halfedge
	HE_face* face;	// Incident face
	HE_edge* next;	// Next halfedge around the face
};

struct HE_vert
{
	int vertID;			// Vertex ID
	float x, y, z;		// Vertex coordinates
	float nx, ny, nz;	// Vertex normals
	bool selected;		// Selected for FFD
	HE_edge* edge;		// Halfedge emanating from vertex
};

struct HE_face
{
	int faceID;			// Face ID
	int v1, v2, v3;		// Vertex IDs
	float nx, ny, nz;	// Face normals
	HE_edge* edge;		// Halfedge bordering the face
};

struct HE_line
{
	unsigned int startVert, endVert;
	bool operator < (const HE_line &other) const
	{
		if(startVert < other.startVert)
			return true;
		else if(startVert == other.startVert)
			return endVert < other.endVert;
		else
			return false;
	}
	HE_line(){};
	HE_line(const HE_line &vert) {
		startVert = vert.startVert;
		endVert = vert.endVert;
	}
	HE_line(unsigned int v1, unsigned int v2) {
		startVert = v1;
		endVert = v2;
	}
};

#endif

       第三步,需要加载模型,并计算顶点和面的法线以及设置包围盒,头文件的内容给读者展示如下:

#ifndef mesh_h
#define mesh_h

#include <cstdlib>
#include <GL\glui.h>
//#include <GL\freeglut.h>
#include <map>
#include <Windows.h>
#include <WTypes.h>

using namespace std;
// Mesh functions
void loadMesh(std::string path);
void findMinMax(float tempX, float tempY, float tempZ);
void calcFaceNorm(void);
void calcVertNorm(void);
void setToOrigin(void);
void scaleToOrigin(void);

// Draw functions
double editLineWidth(int width);
void drawMeshWire(void);
void drawMeshSolid(void);
void drawMeshPoint(void);
void drawBbox(void);
void drawGrid(void);
void drawAxes(void);

// Dialog functions
void createFileDialog(void);
int createMsgbox(int msgboxID);
enum
{
	ERR_no_mesh = 0,
	ERR_no_pt
};

#endif
具体实现会在博客的末尾给读者展示。

第四步,开始对模型的变形处理,在这里并没有使用什么高深的算法,其实就是对它们的顶点进行矩阵变换。首先进行的是世界变换,投影变换,最后是视口变换。

变形的函数如下所示:

void deformMesh()
{
	copyVertices = new Vertex[numOfVert]();
	for(int count = 0; count < numOfVert; count++) {
		/* Create copy of vertex position */
		copyVertices[count].x = v0[count].x;
		copyVertices[count].y = v0[count].y;
		copyVertices[count].z = v0[count].z;

		/* Recalculate new minimum/maximum vertices */
		if(vertices[count].selected) {
			findSelectedMinMax(vertices[count].x, vertices[count].y, vertices[count].z);
		}
	}
	for(int count = 0; count < numOfVert; count++) {
		/* Only calculate new position of vertex if selected */
		if(vertices[count].selected) {
			/* Highlights selected vertex */
			glColor3f(0.5, 0.5, 1);
			glPushMatrix();
				glTranslatef(copyVertices[count].x, copyVertices[count].y, copyVertices[count].z);
				glutSolidSphere(0.1, 20, 20);
			glPopMatrix();
			glColor3f(1, 1, 1);

			vertices[count].x = 0;
			vertices[count].y = 0;
			//vertices[count].z = 0;

			/* Calculate the FFD position and update vertex position */
			Vertex tempResult = trivariateBernstein(copyVertices[count]);
			vertices[count].x = tempResult.x;
			vertices[count].y = tempResult.y;
			vertices[count].z = tempResult.z;
			//cout << "Vertex (old): " << vertices[count].x << ", " << vertices[count].y << ", " << vertices[count].z << "\n";
			//cout << "Vertex (new): " << vertices[count].x << ", " << vertices[count].y << ", " << vertices[count].z << "\n";
		}
	}
}
当你看到代码后,第一印象就是太简单了,确实是这样,变形其实就是将需要改变的顶点进行一系列矩阵变换,在这里使用了OpenGL的库函数接口
			glPushMatrix();
				glTranslatef(copyVertices[count].x, copyVertices[count].y, copyVertices[count].z);
				glutSolidSphere(0.1, 20, 20);
			glPopMatrix();
进行顶点的移动,同时调用了函数
trivariateBernstein

进行顶点的变换以及权值计算,代码如下所示:

Vertex trivariateBernstein(Vertex vert)
{
	Vertex stuVert;
	stuVert.x = 0;
	stuVert.y = 0;
	stuVert.z = 0;
	stuVert = convertToSTU(vert);
	/*cout << "Vertex (XYZ): " << vert.x << ", " << vert.y << ", " << vert.z << "\n";
	cout << "Vertex (STU): " << stuVert.x << ", " << stuVert.y << ", " << stuVert.z << "\n";*/
	
	double weight = 0.0;
	Vertex convertedVert;
	convertedVert.x = 0;
	convertedVert.y = 0;
	convertedVert.z = 0;

	/* Performing summations using for loops */
	for(int i = 0; i <= l; i++) {
		for(int j = 0; j <= m; j++) {
			for(int k = 0; k <= n; k++) {
				weight = 0;
				weight = bernsteinPoly(n, k, stuVert.z) * bernsteinPoly(m, j, stuVert.y) * bernsteinPoly(l, i, stuVert.x);
				convertedVert.x += weight * lattice[i][j][k].x;
				convertedVert.y += weight * lattice[i][j][k].y;
				convertedVert.z += weight * lattice[i][j][k].z;
			}
		}
	}
	return convertedVert;
}
没有使用任何算法,实现的就是模型的矩阵变换。

最后一步就是实现,首先要做的事情是选择需要变形的顶点或者区域,使用了函数如下所示:

void moveSelectedCP(int endX, int endY)
{
	
	/* Initialise movement variables */
	float moveX = 0;
	float moveY = 0;
	moveX = (endX-startX);
	moveY = (endY-startY);

	
	/* Set end point to start point to find next round of moved distance */
	startX = endX;
	startY = endY;

	if(!(moveX == 0) || !(moveY == 0)) {
		//lattice2d[sel_i][sel_j][sel_k].x = endX;
		//lattice2d[sel_i][sel_j][sel_k].y = endY;
		lattice2d[sel_i][sel_j][sel_k].x += moveX;
		lattice2d[sel_i][sel_j][sel_k].y += moveY;

		/* Convert to 3D to change position of 3D CP */
		gluUnProject(lattice2d[sel_i][sel_j][sel_k].x, lattice2d[sel_i][sel_j][sel_k].y, lattice2d[sel_i][sel_j][sel_k].z,
						modelView, projection, viewport,
						&lattice[sel_i][sel_j][sel_k].x, &lattice[sel_i][sel_j][sel_k].y, &lattice[sel_i][sel_j][sel_k].z);
		cout << "Move CP: [" << sel_i << "][" << sel_j << "][" << sel_k << "]\n";
		//drawLattice();
		//deformMesh();
	}
}
在再main函数中调用封装的函数即可。其实做起来还是比较简单的,就是对变形的部分做了矩阵的运算。大家关注上面加粗的部分。

代码下载地址:链接:http://pan.baidu.com/s/1i5klXXV 密码:m1f5


0
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:214206次
    • 积分:4152
    • 等级:
    • 排名:第7417名
    • 原创:189篇
    • 转载:1篇
    • 译文:0篇
    • 评论:142条
    著作书籍
    《Unity3D实战核心技术详解》电子工业出版社 《手把手教你架构3D游戏引擎》电子工业出版社 《Cocos2dx 3.x 图形学渲染技术讲解》电子工业出版社
    博客专栏
    最新评论