2. 程式人生 > >計算機圖形學：虛擬服裝設計

計算機圖形學：虛擬服裝設計

阿新 • • 發佈：2018-11-25

要求：

（1）讀取給定的模型檔案並能夠分別使用線模式、點模式、填充模式進行渲染。

（2）對模型進行平移、縮放和旋轉。

（3）改變視點，對模型進行三維漫遊。

（4）為模型中的衣裙新增物理模型，要求能夠拖動某個節點或面片展示區域性形變效果。

（5）為模型中的衣裙新增碰撞檢測，使得衣裙能在重力作用下懸掛在衣架上，以及模擬風力作用下的衣裙飄動效果。


一、實驗環境

1. 作業系統：windows 8.1

2. 程式設計軟體：visual studio 2013

 

二、實驗內容



l 關鍵程式碼分析

（1）讀取obj模型

void ReadPIC()

 {

int dianshu;

dianshu = 0;

 

ifstream ifs("dress.obj"); //用檔案輸入流讀入檔案

string s;

Mian *f;

POINT3 *v;

FaXiangLiang *vn;

WenLi *vt;

while (getline(ifs, s))  //從標準輸入流中讀取一個字串，儲存到字串string(物件)中

{

if (s.length()<2)continue

;

if (s[0] == 'v')

{

if (s[1] == 't') //vt 紋理 0.001992 0.001992

{

istringstream in(s); //定義一個字串輸入流的物件in,將s中所包含的字串放入in 物件中

    vt = new WenLi();

string head;

in >> head >> vt->TU >> vt->TV;

m_pic.VT.push_back(*vt);

}

else if (s[1] == 'n') //vn 法向量 0.000000 -1.000000 0.000000

{

istringstream in(s);

vn = new FaXiangLiang();

string head;

in >> head >> vn->NX >> vn->NY >> vn->NZ;

m_pic.VN.push_back(*vn);

}

else //v 點 - 0.500000 - 0.500000 0.500000

{

istringstream in(s);

v = new POINT3();

string head;

in >> head >> v->X >> v->Y >> v->Z;

m_pic.V.push_back(*v);

dianshu++;

}

//printf("dianshu=%d",dianshu);

}

else if (s[0] == 'f')//f 面 1/1/1 2/2/2 3/3/3  這個面的頂點、紋理座標、法向量的索引  //f 2443//2656 2442//2656 2444//2656 面

{

for (int k = s.size() - 1; k >= 0; k--)

{

if (s[k] == '/')s[k] = ' ';

}

istringstream in(s);

f = new Mian();

string head;

in >> head;

int i = 0;

while (i<3)

{

if (m_pic.V.size() != 0)

{

in >> f->V[i];

f->V[i] -= 1;

}

if (m_pic.VT.size() != 0)

{

in >> f->T[i];

f->T[i] -= 1;

}

if (m_pic.VN.size() != 0)

{

in >> f->N[i];

f->N[i] -= 1;

}

i++;

}

m_pic.F.push_back(*f);

}

}

}

讀取檔案dress.obj，依據檔案中的索引將資料分別存入法向量m_pic.VN，紋理座標m_pic.VT,頂點m_pic.V，依據面的索引，法向量紋理座標頂點按所在三角形面分別存入m_pic.F中，方便後面繪製以及變換操作。由於資料量大，讀取檔案函式ReadPIC()設定了標誌位，只在檔案啟動時執行1次。

（2）繪製模型

void InitScene()

{

static GLint flag = 1;//設定標誌 位，第一次讀取obj檔案到陣列，後面直接顯示

if (flag == 1)

{

flag = 0;

ReadPIC();

}

 

glClearColor(1.000f, 1.000f, 1.000f, 1.0f); //Background color

// TODO: Replace the following sample code with your initialization code.

// Activate lighting and a light source

//用於啟用各種功能。具體功能由引數決定。與glDisable相對應。glDisable用以關閉各項功能。

glEnable(GL_LIGHT0);//啟用0號燈到7號燈(光源)　　光源要求由函式glLight函式來完成

glEnable(GL_LIGHTING);//啟用燈源

glEnable(GL_DEPTH_TEST);//啟用深度測試。　　根據座標的遠近自動隱藏被遮住的圖形（材料）

glEnable(GL_TEXTURE_2D);   // 啟用二維紋理

// Define material parameters

static GLfloat glfMatAmbient[] = { 0.000f, 0.450f, 1.000f, 1.0f };

static GLfloat glfMatDiffuse[] = { 0.000f, 0.000f, 0.580f, 1.0f };

static GLfloat glfMatSpecular[] = { 1.000f, 1.000f, 1.000f, 1.0f };

static GLfloat glfMatEmission[] = { 0.000f, 0.000f, 0.000f, 1.0f };

static GLfloat fShininess = 128.000f;

// Set material parameters

//指定用於光照計算的當前材質屬性。引數face的取值可以是GL_FRONT、GL_BACK或GL_FRONT_AND_BACK，指出材質屬性將應用於物體的哪面。

//void glMaterial{if}(GLenum face, GLenum pname, TYPE param);

glMaterialfv(GL_FRONT, GL_AMBIENT, glfMatAmbient);

glMaterialfv(GL_FRONT, GL_DIFFUSE, glfMatDiffuse);

glMaterialfv(GL_FRONT, GL_SPECULAR, glfMatSpecular);

glMaterialfv(GL_FRONT, GL_EMISSION, glfMatEmission);

glMaterialf(GL_FRONT, GL_SHININESS, fShininess);

}

 

void DrawScene()

 {

// TODO: Replace the following sample code with your code to draw the scene.

glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);  // 清除螢幕及深度快取

glLoadIdentity(); // 重置模型觀察矩陣

 

gluLookAt(viewer[0], viewer[1], viewer[2], 0.0f, -4.0f, -9.0f, 0.0, 1.0, 0.0);

 

glTranslatef(0.0f, -13.0f, -8.0f); // 移入螢幕 5.0

 

if (trackballMove&&translating&&!scaling&&!rotating)

{

glTranslatef(tran[0], tran[1], tran[2]);

//printf("tran[0]=%f, tran[1]=%f, tran[2]=%f\n", tran[0], tran[1], tran[2]);

}

if (trackballMove&&!translating&&scaling&&!rotating)

{

glScalef(scale[0], scale[1], scale[2]);

}

if (trackballMove&&!translating&&!scaling&&rotating)

{

glRotatef(angle0, axis[0], axis[1], axis[2]);//旋轉（角度，軸）

}

GLCube();// Draw a cube

glFlush();

}

 

void GLCube()

{

if (dian_model&&!xian_model&&!mian_model)

//if (0)

{

for (int i = 0; i < (m_pic.F.size()); i++)

//for (int i = 0; i<7959; i++)

{

glBegin(GL_POINTS);//點模式

if (qizhi0 == 1)

{

printf("dianmoshi");

qizhi0 = 0;

}

/*glVertex3f(m_pic.V[m_pic.F[3 * i].V[0]].X / YU, m_pic.V[m_pic.F[3 * i].V[0]].Y / YU, m_pic.V[m_pic.F[3 * i].V[0]].Z / YU);

glVertex3f(m_pic.V[m_pic.F[3 * i].V[1]].X / YU, m_pic.V[m_pic.F[3 * i].V[1]].Y / YU, m_pic.V[m_pic.F[3 * i].V[1]].Z / YU);

glVertex3f(m_pic.V[m_pic.F[3 * i].V[2]].X / YU, m_pic.V[m_pic.F[3 * i].V[2]].Y / YU, m_pic.V[m_pic.F[3 * i].V[2]].Z / YU);

*/

 

glVertex3f(m_pic.V[m_pic.F[i].V[0]].X / YU, m_pic.V[m_pic.F[i].V[0]].Y / YU, m_pic.V[m_pic.F[i].V[0]].Z / YU);

glVertex3f(m_pic.V[m_pic.F[i].V[1]].X / YU, m_pic.V[m_pic.F[i].V[1]].Y / YU, m_pic.V[m_pic.F[i].V[1]].Z / YU);

glVertex3f(m_pic.V[m_pic.F[i].V[2]].X / YU, m_pic.V[m_pic.F[i].V[2]].Y / YU, m_pic.V[m_pic.F[i].V[2]].Z / YU);

 

//glVertex3f(10*(model.vertices[3 * i]), 10*(model.vertices[3 * i + 1]), 10*(model.vertices[3 * i + 2]));

glEnd();

}

}

if (!dian_model&&xian_model&&!mian_model)

//if (0)

{

for (int i = 0; i<(m_pic.F.size()); i++)

{

glBegin(GL_LINES);//線模式

if (qizhi1 == 1)

{

printf("xianmoshi");

qizhi1 = 0;

}

glVertex3f(m_pic.V[m_pic.F[i].V[0]].X / YU, m_pic.V[m_pic.F[i].V[0]].Y / YU, m_pic.V[m_pic.F[i].V[0]].Z / YU);

glVertex3f(m_pic.V[m_pic.F[i].V[1]].X / YU, m_pic.V[m_pic.F[i].V[1]].Y / YU, m_pic.V[m_pic.F[i].V[1]].Z / YU);

//glVertex3f(m_pic.V[m_pic.F[i].V[2]].X / YU, m_pic.V[m_pic.F[i].V[2]].Y / YU, m_pic.V[m_pic.F[i].V[2]].Z / YU);

glEnd();

glBegin(GL_LINES);

//glVertex3f(m_pic.V[m_pic.F[i].V[0]].X / YU, m_pic.V[m_pic.F[i].V[0]].Y / YU, m_pic.V[m_pic.F[i].V[0]].Z / YU);

glVertex3f(m_pic.V[m_pic.F[i].V[1]].X / YU, m_pic.V[m_pic.F[i].V[1]].Y / YU, m_pic.V[m_pic.F[i].V[1]].Z / YU);

glVertex3f(m_pic.V[m_pic.F[i].V[2]].X / YU, m_pic.V[m_pic.F[i].V[2]].Y / YU, m_pic.V[m_pic.F[i].V[2]].Z / YU);

glEnd();

glBegin(GL_LINES);

glVertex3f(m_pic.V[m_pic.F[i].V[0]].X / YU, m_pic.V[m_pic.F[i].V[0]].Y / YU, m_pic.V[m_pic.F[i].V[0]].Z / YU);

//glVertex3f(m_pic.V[m_pic.F[i].V[1]].X / YU, m_pic.V[m_pic.F[i].V[1]].Y / YU, m_pic.V[m_pic.F[i].V[1]].Z / YU);

glVertex3f(m_pic.V[m_pic.F[i].V[2]].X / YU, m_pic.V[m_pic.F[i].V[2]].Y / YU, m_pic.V[m_pic.F[i].V[2]].Z / YU);

glEnd();

}

}

if (!dian_model&&!xian_model&&mian_model)

//if (1)

{

for (int i = 0; i<(m_pic.F.size()); i++)

{

//printf("m_pic.F.size()=%d", m_pic.F.size());//16213

glBegin(GL_TRIANGLES);//面模式// 繪製三角形

if (qizhi2 == 1)

{

printf("mianmoshi");

qizhi2 = 0;

}

if (m_pic.VT.size() != 0)glTexCoord2f(m_pic.VT[m_pic.F[i].T[0]].TU, m_pic.VT[m_pic.F[i].T[0]].TV);  //紋理

if (m_pic.VN.size() != 0)glNormal3f(m_pic.VN[m_pic.F[i].N[0]].NX, m_pic.VN[m_pic.F[i].N[0]].NY, m_pic.VN[m_pic.F[i].N[0]].NZ);//法向量

glVertex3f(m_pic.V[m_pic.F[i].V[0]].X / YU, m_pic.V[m_pic.F[i].V[0]].Y / YU, m_pic.V[m_pic.F[i].V[0]].Z / YU); // 上頂點

//printf("dian:%f,%f,%f\n",m_pic.V[m_pic.F[i].V[0]].X / YU, m_pic.V[m_pic.F[i].V[0]].Y / YU, m_pic.V[m_pic.F[i].V[0]].Z / YU);

if (m_pic.VT.size() != 0)glTexCoord2f(m_pic.VT[m_pic.F[i].T[1]].TU, m_pic.VT[m_pic.F[i].T[1]].TV);  //紋理

if (m_pic.VN.size() != 0)glNormal3f(m_pic.VN[m_pic.F[i].N[1]].NX, m_pic.VN[m_pic.F[i].N[1]].NY, m_pic.VN[m_pic.F[i].N[1]].NZ);//法向量

glVertex3f(m_pic.V[m_pic.F[i].V[1]].X / YU, m_pic.V[m_pic.F[i].V[1]].Y / YU, m_pic.V[m_pic.F[i].V[1]].Z / YU); // 左下

if (m_pic.VT.size() != 0)glTexCoord2f(m_pic.VT[m_pic.F[i].T[2]].TU, m_pic.VT[m_pic.F[i].T[2]].TV);  //紋理

if (m_pic.VN.size() != 0)glNormal3f(m_pic.VN[m_pic.F[i].N[2]].NX, m_pic.VN[m_pic.F[i].N[2]].NY, m_pic.VN[m_pic.F[i].N[2]].NZ);//法向量

glVertex3f(m_pic.V[m_pic.F[i].V[2]].X / YU, m_pic.V[m_pic.F[i].V[2]].Y / YU, m_pic.V[m_pic.F[i].V[2]].Z / YU); // 右下

glEnd(); // 三角形繪製結束

}

}

}

InitScene()函式中設定了場景的光源和模型的材質。DrawScene()中繪製模型，呼叫函式GLCube()實現具體繪製，並由標誌位設定繪製模式：點模式，線模式和麵模式；同時使用庫函式實現模型旋轉縮放平移變換操作，具體引數由虛擬球設定；同時gluLookAt(viewer[0], viewer[1], viewer[2], 0.0f, -4.0f, -9.0f, 0.0, 1.0, 0.0);函式實現三維漫遊，viewer[0], viewer[1], viewer[2]數值由鍵盤設定，後面的視點方向等用於使模型位於世界座標系正中間。

（3）物理模型

Cloth::Cloth(GLMmodel *model)

{

//

SpringConstant = 1000;//彈力

DampingFactor = 10;//阻尼

density = 1.2;//密度

drag = 0.4;//阻力

 

//set vert

Myvertex vert;

verts_.push_back(vert);

for (int i = 1; i <= 7959; i++) // m_pic.V[m_pic.F[i].V[2]].X

{

//printf("%d\n", m_pic.V.size());//7959

vert.id_ = i;

vert.position_ = point(model->vertices[3 * i], model->vertices[3 * i + 1], model->vertices[3 * i + 2]);//point: vec

//printf("[1]=%lf\n", model->vertices[3 * i]);

//printf("[2]=%lf\n", model->vertices[3 * i+1]);

//printf("[3]=%lf\n", model->vertices[3 * i+2]);

//vert.texCoord_ = Vec2f(model->)

verts_.push_back(vert);

}

 

 

//calculate hanger_aabb 衣架aabb包圍盒

for (int i = 3381; i <= model->numvertices; i++)

{

if (model->vertices[3 * i] > hanger_aabb.max_x)

{

hanger_aabb.max_x = model->vertices[3 * i];

}

if (model->vertices[3 * i] < hanger_aabb.min_x)

{

hanger_aabb.min_x = model->vertices[3 * i];

}

if (model->vertices[3 * i + 1] > hanger_aabb.max_y)

{

hanger_aabb.max_y = model->vertices[3 * i + 1];

}

if (model->vertices[3 * i + 1] < hanger_aabb.min_y)

{

hanger_aabb.min_y = model->vertices[3 * i + 1];

}

if (model->vertices[3 * i + 2] > hanger_aabb.max_z)

{

hanger_aabb.max_z = model->vertices[3 * i + 2];

}

if (model->vertices[3 * i + 2] < hanger_aabb.min_z)

{

hanger_aabb.min_z = model->vertices[3 * i + 2];

}

}

 

update_aabb(); //calculate_cloth_aabb

 

for (int i = 1; i <= model->numvertices; i++) //if points in hanger_aabb,point is static

{

if (p_in_aabb(verts_[i].position_))

{

verts_[i].is_fixed = true;

}

else

verts_[i].is_fixed = false;

}

 

//set face

Myface face;

for (int i = 0; i < model->numtriangles; i++)

{

face.id_ = i;

int nid = model->triangles[i].findex;

face.verId = ivec3(model->triangles[i].vindices[0], model->triangles[i].vindices[1], model->triangles[i].vindices[2]);

face.normal_ = Vec3f(model->facetnorms[nid * 3], model->facetnorms[nid * 3 + 1], model->facetnorms[nid * 3 + 2]);

for (int j = 0; j < 3; j++)

{

int tid = model->triangles[i].tindices[j];

face.texCoord_[j] = Vec2f(model->texcoords[tid * 2], model->texcoords[2 * tid + 1]);

}

int id[3];

for (int j = 0; j < 3; j++)

{

id[j] = face.verId[j];

face.vertex3[j] = &verts_[id[j]];

}

face.s = fabs(((verts_[id[1]].position_ - verts_[id[0]].position_)CROSS(verts_[id[2]].position_ - verts_[id[0]].position_)).length() / 2.0);

faces_.push_back(face);

 

}

 

//delete some bad faces

for (int i = 0; i < 6546; i++)

{

if ((faces_[i].normal_[0] * faces_[i].vertex3[0]->position_[0] + faces_[i].normal_[2] * faces_[i].vertex3[0]->position_[2]) < -0.05)

{

faces_[i].badface = true;

}

}

 

 

 

//delete some bad verts,and find bound

for (int i = 0; i < 6546; i++)

{

if (faces_[i].badface)

{

for (int j = 0; j < 3; j++)

{

faces_[i].vertex3[j]->kengdie = true;

faces_[i].vertex3[j]->is_bound = true;

}

 

}

}

for (int i = 0; i < 6546; i++)

{

if (!faces_[i].badface)

{

for (int j = 0; j < 3; j++)

{

faces_[i].vertex3[j]->kengdie = false;

}

}

}

 

 

 

set_neighbor();

set_weight();

//Update(float deltatime);

 

}

}

服裝模型類，設定了服裝的物理模型，為服裝添加了SpringConstant彈力，DampingFactor阻尼，density密度，drag 阻力屬性；分別找到衣架和服裝的頂點取值範圍，確定衣架和服裝的包圍盒；形變是針對服裝模型的三角形面片作變化，將模型的基本資料儲存到face中；由於原模型中裙襬雙層布料，碰撞檢測的效果不好，故刪去裙襬內襯，即badface。

 

 

（4）碰撞檢測

static void idle() {

if (phy_){ if (loop < 16){ Update(); loop++; } else loop = 0; }

}

 

void keys(unsigned char key, int x, int y)

{

/* Use x, X, y, Y, z, and Z keys to move viewer */

if (key == 'x') viewer[0] -= 1.0;

if (key == 'X') viewer[0] += 1.0;

if (key == 'y') viewer[1] -= 1.0;

if (key == 'Y') viewer[1] += 1.0;

if (key == 'z') viewer[2] -= 1.0;

if (key == 'Z') viewer[2] += 1.0;

//display();

if (!phy_){

DrawScene();

}

 

if (key == 'p'){

phy_ = true; printf("phy_=true\n");

}

 

if (key == 'o'){

wind.x += 0.5;

wind.Print("New Wind:");

}

if (key == 'i'){

wind.y += 0.5;

wind.Print("New Wind:");

}

if (key == 'u'){

wind.z += 0.5;

wind.Print("New Wind:");

}

}

 

void Cloth::Update(float deltatime)

{

//set gravity

for (int i = 1; i < verts_.size(); i++)

{

verts_[i].applyForce(Vec3f(0, -0.1, 0)*verts_[i].weight);

//verts_[i].applyForce(Vec3f(0, 0, 0.01));

//verts_[i].Update(deltatime);

}

 

 

change_wind();

 

//set wind

for (int i = 0; i < faces_.size(); i++)

{

faces_[i].updateface();

}

 

 

//set  spring

for (int i = 0; i < springDampers.size(); i++)

{

springDampers[i].CalculateForces();

}

 

collision_detection();//碰撞檢測

 

//update verts

for (int i = 1; i < verts_.size(); i++)

{

verts_[i].Update(deltatime);

}

 

//update normals

update_normal();

}

 

void Cloth::change_wind()

{

for (int i = 0; i < 6546; i++)

{

faces_[i].Changewind(wind);

}

}

 

void Cloth::update_normal()

{

for (int i = 0; i < faces_.size(); i++)

{

Vec3f p = verts_[faces_[i].verId[0]].position_;

Vec3f q = verts_[faces_[i].verId[1]].position_;

Vec3f r = verts_[faces_[i].verId[2]].position_;

Vec3f normal_update = ((p - r) CROSS(q - r));

normal_update.normalize();

 

//normal is  continuous variable

 

if (normal_update DOT faces_[i].normal_>0)

{

faces_[i].normal_ = normal_update;

 

}

else

{

faces_[i].normal_ = -normal_update;

}

}

}

 

void Cloth::set_neighbor()

{

for (int i = 0; i < faces_.size(); i++)

{

if (!faces_[i].badface)

{

int vid1;

int vid2;

for (int j = 0; j < 3; j++)

{

vid1 = faces_[i].verId[j];

vid2 = faces_[i].verId[(j + 1) % 3];

verts_[vid1].neighborIdx.push_back(vid2);

float dis;

dis = (verts_[vid1].position_ - verts_[vid2].position_).length();

verts_[vid1].nei_dis.push_back(dis);

verts_[vid1].neighborfaceId.push_back(i);

if (verts_[vid1].is_bound&&verts_[vid2].is_bound)

{

verts_[vid2].neighborIdx.push_back(vid1);

verts_[vid2].nei_dis.push_back(dis);

}

}

 

}

}

 

for (int i = 1; i < verts_.size(); i++)

{

verts_[i].degree_ = verts_[i].neighborIdx.size();

}

 

for (int i = 0; i < faces_.size(); i++)

{

if (!faces_[i].badface)

{

int vid1;

int vid2;

for (int j = 0; j < 3; j++)

{

vid1 = faces_[i].verId[j];

vid2 = faces_[i].verId[(j + 1) % 3];

springDampers.push_back(springDamper(&verts_[vid1], &verts_[vid2], SpringConstant, DampingFactor));//彈簧阻尼

 

}

}

}

 

}

 

void Cloth::collision_detection()

{

update_aabb();

for (int i = 0; i < 40; i++)

for (int j = 0; j < 40; j++)

for (int k = 0; k < 40; k++)

{

voxels[i][j][k].clear();

}

 

 

for (int i = 0; i < 6546; i++)

{

faces_[i].voxel_pos.clear();

}

for (int i = 0; i < 6546; i++)

{

if (!faces_[i].badface)

{

aabb ab_ = faces_[i].get_aabb();

int l = floor((ab_.min_x - cloth_aabb.min_x) / (cloth_aabb.max_x - cloth_aabb.min_x)*40.f);

int m = floor((ab_.min_y - cloth_aabb.min_y) / (cloth_aabb.max_y - cloth_aabb.min_y)*40.f);

int n = floor((ab_.min_z - cloth_aabb.min_z) / (cloth_aabb.max_z - cloth_aabb.min_z)*40.f);

if (l < 0)

l = 0;

if (m < 0)

m = 0;

if (n < 0)

n = 0;

voxels[l][m][n].push_back(i);

faces_[i].voxel_pos = ivec3(l, m, n);

}

}

 

for (int i = 0; i < 6546; i++)

{

ivec3 voxel = faces_[i].voxel_pos;

for (int j = 0; j < voxels[voxel[0]][voxel[1]][voxel[2]].size(); j++)

{

int face_id = voxels[voxel[0]][voxel[1]][voxel[2]][j];

bool b = tri_collision_detection(i, face_id);

}

}

/*for (int i = 0; i < 6546; i++)

for (int j = 0; j < 6546; j++)

{

tri_collision_detection(i, j);

}*/

 

}

 

bool Cloth::tri_collision_detection(int id1, int id2)

{

if ((faces_[id1].normal_ DOT faces_[id2].normal_)>0)

return false;

if (faces_[id1].badface || faces_[id2].badface)

return false;

Vec3f dis(0, 0, 0);

for (int i = 0; i < 3; i++)

{

dis += (faces_[id1].vertex3[i]->position_ - faces_[id2].vertex3[i]->position_) / 3.0f;

}

if (dis.length() > 0.001)

return false;

 

else

{

Vec3f n = Vec3f(0, 0, 0);

 

for (int i = 0; i < 3; i++)

{

n = (faces_[id2].normal_ DOT(faces_[id1].vertex3[i]->position_ - faces_[id2].vertex3[1]->position_))*faces_[id2].normal_;

if (n.length() < 0.01)

{

n.normalize();

if (n DOT faces_[id1].vertex3[i]->velocity < 0)

{

Vec3f vv = Vec3f(0, 0, 0);

vv = (n DOT faces_[id1].vertex3[i]->velocity)*-1.1f*n;

faces_[id1].vertex3[i]->velocity += vv;

/*for (int j = 0; j < 3; j++)

{

faces_[id2].vertex3[j]->velocity += vv / -1.1f  / 3.0f;

}*/

//faces_[id1].vertex3[i]->velocity = Vec3f(0, faces_[id1].vertex3[i]->velocity[1], 0);

 

}

}

}

return true;

}

 

}

 

由鍵盤設定是否開啟碰撞檢測，按下’p’鍵即開啟碰撞檢測，風力大小也由鍵盤設定。碰撞檢測具體實現流程見上述程式碼。

 

 


 


 

 

 

 

 

相關推薦

計算機圖形學：虛擬服裝設計

要求： （1）讀取給定的模型檔案並能夠分別使用線模式、點模式、填充模式進行渲染。 （2）對模型進行平移、縮放和旋轉。 （3）改變視點，對模型進行三維漫遊。 （4）為模型中的衣裙新增物理模型，要求能夠拖動某個節點或面片展示區域性形變效果。 （5）為模型中的衣裙新增碰撞檢測，使

計算機圖形學：任意點相對於任意平面的反射矩陣

問題：已知空間中任意平面n*(x-x0，y-y0，z-z0)=0和空間中任意點V，其中n為該平面的法向量，x0，y0，z0為該平面上的一點。 求：該點相對於該平面的映象點v'和V到V'的轉換矩陣，如下圖所示： 該平面方程也可以表示為n*p+d=0，其中p=(x，

計算機圖形學：基本二維幾何變換

其中，4個元素rjk是多重旋轉選項，元素trx和try是平移項。座標位置的剛體變化有時也稱剛體運動（rigid-motion）變換。變換後坐標位置間的所有角度和距離都不變化。此外，矩陣具有其左上角的2*2矩陣是一個正交矩陣的特性。這說明，假如將子矩陣 每一行（或者每一列）作為一個向量，那麼兩個行向量，那麼兩個

計算機圖形學：基於3D遊戲開發——第二章 頂點處理機制

渲染管線：頂點處理過程——光柵化過程——片元處理過程——輸出合併操作。 光柵化階段：根據頂點對多邊形加以整合，並將各個多邊形轉換為片元集合。片元定義為一組資料，並對顏色緩衝區中的畫素進行更行。 片元處理：對各個片元進行操作（如紋理操作），進而進行各種操作確定片元的顏色值。

清華髮布：2018計算機圖形學研究報告（附下載）

導讀：報告包括概述篇、技術篇、人才篇、會議篇、應用篇、趨勢篇，本文擷取概述篇部分內容。在公眾號後

【計算機圖形學】基本圖形元素：直線的生成演算法

直線的DDA演算法【演算法介紹】設直線之起點為（x1，y1），終點為（x2，y2），則斜率m為： 直線中的每一點座標都可以由前一點座標變化一個增量（Dx, Dy）而得到，即表示為遞迴式：並有關係：Dy = m • Dx。遞迴式的初值為直線的起點（x1, y1），這樣，就可以用

計算機圖形學 學習筆記（一）：概述，直線掃描轉換演算法：DDA，中點畫線演算法，Bresenham演算法

前言 感謝中國農大 趙明老師的分享~ 現在我要為我自己走向遊戲程式設計打下基石~ 1 計算機圖形學概論 1.1 計算機圖形學課程簡介 《計算機圖形學》是計算機、地理資訊系統、應用數學、機械、建築等專業本科教學中的一門重要的專業基礎課 如影

計算機圖形學 學習筆記（七）：二維圖形變換：平移，比例，旋轉，座標變換等

在圖形學中，有兩大基本工具：向量分析，圖形變換。本文將重點講解向量和二維圖形的變換。 5.1 向量基礎知識 我們所使用的所有點和向量都是基於某一座標系定義的，比如左手座標系或者右手座標系。 從幾何的角度來看，向量是具有長度和方向的實體，但是沒有

計算機圖形學 學習筆記（二）：多邊形掃描轉換：X掃描線演算法 和 改進的X掃描線演算法

光柵圖形學演算法 2.4 多邊形掃描轉換-X掃描線演算法 多邊形的掃描轉換和區域填充這個問題是怎麼樣在離散的畫素集上表示一個連續的二維圖形？ 多邊形有兩種重要的表示方法：頂點表示和點陣表示 頂點表示 頂點表示是用多邊形的頂點序列來表

【演算法】計算機圖形學的一些經典小題：判斷點在多邊形內，隨機生成三角形內的點，判斷兩個矩形是否相交等

前幾天面試的時候被問到了，如何隨機在三角形內生成點，我按照我的想法回答了一遍，但覺得回答的不夠好。最後面試官說了一個最優的方法。覺得不錯，順帶總結一下最近看到的一些關於計算機圖形學方面的經典小題，知乎上看到的還有Leetcode上的 1.判斷一個點

計算機圖形學 學習筆記（十二）：顏色模型，簡單 / 增量 光照模型

接下來將介紹真實感圖形學的內容。 真實感圖形學包括：顏色模型，簡單光照模型，增量光照模型，區域性光照模型和投射模型，整體光照模型，紋理對映和陰影處理。 真實感圖形學 真實感圖形學研究的是什麼？簡單來說，就是希望用計算機生成像照相機拍的照片一樣逼真地

【3D計算機圖形學】變換矩陣、歐拉角、四元數

遭遇 unit 額外 star 應該 detail 導致 print uic 【3D計算機圖形學】變換矩陣、歐拉角、四元數 旋轉矩陣、歐拉角、四元數主要用於：向量的旋轉、坐標系之間的轉換、角位移計算、方位的平滑插值計算。 一、變換矩陣： 首先要區分旋轉矩陣和變換矩陣：

[計算機圖形學 with OpenGL] Chapter10 OpenGL三維觀察程序示例

chap 而不是 max argv func open position style windows 　　10.10節書中給出了一個程序示例，有一個填充正方形，從側面的角度觀察並畫到屏幕上。 　　圖0 　　這裏進一步畫出一個立方體，將相機放入立方體中心，旋轉相機，達到在

計算機圖形學（一） 視頻顯示設備_1_CRT原理

http color size 安裝 ref p s 這一 計算機圖形學 指定 第 1 章 圖形系統概述 如今。計算機圖形學的作用與應用已經得到了廣泛承認。大量的圖形硬件和軟件系統已經應用 到了差點兒全部的領域。通用計算機甚至很多手持計算器也已經

MFC計算機圖形學（2）

mct tid spc DdGzS cin html uem ubd dcs sdsdzi狗聘毫渤口毫http://huiyi.docin.com/hnbkw203d1e5gw濫良瘟侍探蝗http://weibo.com/p/10050563731520645atr4g回救

計算機圖形學-圖形幾何變換

red test position glbegin mage += logs depth window 內容：金字塔的平移以及旋轉的實現，通過鼠標控制金字塔的轉速以及運行窗口的退出 #include <GL/glut.h> #include <stdli

計算機圖形學（二）輸出圖元_3_畫線算法_2_DDA算法

通過 程序 之間 tro 取整 xen git 方程 class DDA算法? ? ? ? 數字微分分析儀(digital differential analyzer, DDA)方法是一種線段掃描轉換算法。基於使用等式(3

計算機圖形學（第2版 於萬波 於碩 編著）第45頁的Bresenham算法有錯誤

str mage mov 步長 分享圖片 圖片 方法總結 tro 計算 計算機圖形學（第2版 於萬波 於碩 編著）第45頁的Bresenham算法有錯誤： 書上本來要寫的是以x為階越步長的方法，但是他寫的是用一部分y為階越步長的方法（其實也寫的不對），最後以x為階越

計算機圖形學-mac系統下Xcode中OpenGL開發環境配置。

配置步驟 ctf 註意 posit 圖片 pen ret open 方式 mac系統下Xcode中OpenGL開發環境配置。 這學期有計算機圖形學的課程，需要用到OpenGL，最近著手開始配置開發環境了，老師上課給的安裝包都是基於windows系統的。網上也是window

CS184.1X 計算機圖形學導論 第7講 V1-3 學習筆記

線上 物體 創建 strong 公式推導 導論 幾何 ng- 解決方法 L7V1：OPENGL 著色：學習動機 1.光照的重要性 1）能夠真正顯示出形狀感知的外觀； 2）準確的著色和光照對對傳達物體的形狀非常重要； 3）著色的方式也十分重要：平面著色（GL_FLAT）、平滑