下面的代码是显示一个旋转立方体的正投影(侧视图,前视图和顶视图)。我采用的方法是齐次坐标。我采取显示立方体的每个边,至于侧视图,我没有考虑立方体顶点的x分量,只使用y和z分量绘制边。这也显示了其他侧面,我怎么能只检测左边的侧面呢?
#include<iostream>
#include "graphics.h"
#include<cmath>
using namespace std;
#define PI 22/(float)7
float x =0, y = 0, z= 0;
void matrixmult(int A[4][8], float R[4][4]){
float res[4][8];
for(int i = 0; i<4; i++){
for(int j = 0; j<8; j++){
res[i][j] = 0;
for(int k = 0; k<4; k++){
res[i][j]+=R[i][k]*A[k][j];
}
}
}
for(int i = 0; i<4; i++){
for(int j = 0; j<8; j++){
A[i][j] = res[i][j];
}
}
}
void rotMataboutX(float Rx[4][4]){
float xrad = PI*x/(float)360;
float cosx = cos(xrad), sinx = sin(xrad);
Rx[1][1]=cosx;
Rx[1][2]=-sinx;
Rx[2][1]=sinx;
Rx[2][2]=cosx;
}
void rotMatAboutY(float Ry[4][4]){
float yrad = PI*y/(float)360;
float cosy = cos(yrad), siny = sin(yrad);
Ry[0][0]=cosy;
Ry[0][2]=siny;
Ry[2][0]=-siny;
Ry[2][2]=cosy;
}
void rotMatAboutZ(float Rz[4][4]){
float zrad = PI*z/(float)360;
float cosz = cos(zrad), sinz = sin(zrad);
Rz[0][0]=cosz;
Rz[0][1]=-sinz;
Rz[1][0]=sinz;
Rz[1][1]=cosz;
}
void sideView(int A[4][8]){
int xdisp = 106, ydisp = 120;
line(xdisp+A[1][0], ydisp+A[2][0], xdisp+A[1][1], ydisp+A[2][1]);
line(xdisp+A[1][1], ydisp+A[2][1], xdisp+A[1][2], ydisp+A[2][2]);
line(xdisp+A[1][2], ydisp+A[2][2], xdisp+A[1][3], ydisp+A[2][3]);
line(xdisp+A[1][3], ydisp+A[2][3], xdisp+A[1][0], ydisp+A[2][0]);
line(xdisp+A[1][4], ydisp+A[2][4], xdisp+A[1][5], ydisp+A[2][5]);
line(xdisp+A[1][5], ydisp+A[2][5], xdisp+A[1][6], ydisp+A[2][6]);
line(xdisp+A[1][6], ydisp+A[2][6], xdisp+A[1][7], ydisp+A[2][7]);
line(xdisp+A[1][7], ydisp+A[2][7], xdisp+A[1][4], ydisp+A[2][4]);
line(xdisp+A[1][0], ydisp+A[2][0], xdisp+A[1][4], ydisp+A[2][4]);
line(xdisp+A[1][1], ydisp+A[2][1], xdisp+A[1][5], ydisp+A[2][5]);
line(xdisp+A[1][6], ydisp+A[2][6], xdisp+A[1][2], ydisp+A[2][2]);
line(xdisp+A[1][7], ydisp+A[2][7], xdisp+A[1][3], ydisp+A[2][3]);
line(xdisp+A[1][7], ydisp+A[2][7], xdisp+A[1][3], ydisp+A[2][3]);
}
void frontView(int A[4][8]){
int xdisp = 319, ydisp = 120;
line(xdisp+A[0][0], ydisp+A[1][0], xdisp+A[0][1], ydisp+A[1][1]);
line(xdisp+A[0][1], ydisp+A[1][1], xdisp+A[0][2], ydisp+A[1][2]);
line(xdisp+A[0][2], ydisp+A[1][2], xdisp+A[0][3], ydisp+A[1][3]);
line(xdisp+A[0][3], ydisp+A[1][3], xdisp+A[0][0], ydisp+A[1][0]);
line(xdisp+A[0][4], ydisp+A[1][4], xdisp+A[0][5], ydisp+A[1][5]);
line(xdisp+A[0][5], ydisp+A[1][5], xdisp+A[0][6], ydisp+A[1][6]);
line(xdisp+A[0][6], ydisp+A[1][6], xdisp+A[0][7], ydisp+A[1][7]);
line(xdisp+A[0][7], ydisp+A[1][7], xdisp+A[0][4], ydisp+A[1][4]);
line(xdisp+A[0][0], ydisp+A[1][0], xdisp+A[0][4], ydisp+A[1][4]);
line(xdisp+A[0][1], ydisp+A[1][1], xdisp+A[0][5], ydisp+A[1][5]);
line(xdisp+A[0][6], ydisp+A[1][6], xdisp+A[0][2], ydisp+A[1][2]);
line(xdisp+A[0][7], ydisp+A[1][7], xdisp+A[0][3], ydisp+A[1][3]);
line(xdisp+A[0][7], ydisp+A[1][7], xdisp+A[0][3], ydisp+A[1][3]);
}
void topView(int A[4][8]){
int xdisp = 534, ydisp = 120;
line(xdisp+A[0][0], ydisp+A[2][0], xdisp+A[0][1], ydisp+A[2][1]);
line(xdisp+A[0][1], ydisp+A[2][1], xdisp+A[0][2], ydisp+A[2][2]);
line(xdisp+A[0][2], ydisp+A[2][2], xdisp+A[0][3], ydisp+A[2][3]);
line(xdisp+A[0][3], ydisp+A[2][3], xdisp+A[0][0], ydisp+A[2][0]);
line(xdisp+A[0][4], ydisp+A[2][4], xdisp+A[0][5], ydisp+A[2][5]);
line(xdisp+A[0][5], ydisp+A[2][5], xdisp+A[0][6], ydisp+A[2][6]);
line(xdisp+A[0][6], ydisp+A[2][6], xdisp+A[0][7], ydisp+A[2][7]);
line(xdisp+A[0][7], ydisp+A[2][7], xdisp+A[0][4], ydisp+A[2][4]);
line(xdisp+A[0][0], ydisp+A[2][0], xdisp+A[0][4], ydisp+A[2][4]);
line(xdisp+A[0][1], ydisp+A[2][1], xdisp+A[0][5], ydisp+A[2][5]);
line(xdisp+A[0][6], ydisp+A[2][6], xdisp+A[0][2], ydisp+A[2][2]);
line(xdisp+A[0][7], ydisp+A[2][7], xdisp+A[0][3], ydisp+A[2][3]);
line(xdisp+A[0][7], ydisp+A[2][7], xdisp+A[0][3], ydisp+A[2][3]);
}
int main(){
int gd = DETECT, gm;
initgraph(&gd, &gm, "C:\\TC\\BGI");
int A[4][8]={{0, 50, 50, 0, 0, 50, 50, 0},
{0, 0, 50, 50, 0, 0, 50, 50},
{0, 0, 0, 0, 50, 50, 50, 50},
{1, 1, 1 ,1 ,1 ,1, 1, 1}};
float cosx = cos(x), sinx = sin(x);
float Rx[4][4]={{1, 0, 0, 0},
{0, cosx, -sinx, 0},
{0, sinx, cosx, 0},
{0, 0, 0, 1}};
float cosy = cos(y), siny = sin(y);
float Ry[4][4]={{cosy, 0, siny, 0},
{0, 1, 0, 0},
{-siny, 0, cosy, 0},
{0, 0, 0, 1}};
float cosz = cos(z), sinz = sin(z);
float Rz[4][4]={{cosz, -sinz, 0, 0},
{sinz, cosz, 0, 0},
{0, 0, 1, 0},
{0, 0, 0, 1}};
cout<<"Matrix A:\n";
for(int i = 0; i<4; i++){
for(int j = 0; j<8; j++){
cout<<A[i][j]<<" ";
}
cout<<endl;
}
sideView(A);
topView(A);
frontView(A);
while(1){
char axisofrot;
cin>>axisofrot;
switch(axisofrot){
case 'X':
case 'x':
cout<<"ROTATING ABOUT X:\n";
x+=5;
cout<<"X:\t"<<x<<"\n";
cout<<"Y:\t"<<y<<"\n";
cout<<"Z:\t"<<z<<"\n";
rotMataboutX(Rx);
cout<<"Matrix Rx:\n";
for(int i = 0; i<4; i++){
for(int j = 0; j<4; j++){
cout<<Rx[i][j]<<" ";
}
cout<<endl;
}
matrixmult(A, Rx);
break;
case 'Y':
case 'y':
cout<<"ROTATING ABOUT Y:\n";
y+=5;
cout<<"X:\t"<<x<<"\n";
cout<<"Y:\t"<<y<<"\n";
cout<<"Z:\t"<<z<<"\n";
rotMatAboutY(Ry);
cout<<"Matrix Ry:\n";
for(int i = 0; i<4; i++){
for(int j = 0; j<4; j++){
cout<<Ry[i][j]<<" ";
}
cout<<endl;
}
matrixmult(A, Ry);
break;
case 'Z':
case 'z':
cout<<"ROTATING ABOUT Z:\n";
z+=5;
cout<<"X:\t"<<x<<"\n";
cout<<"Y:\t"<<y<<"\n";
cout<<"Z:\t"<<z<<"\n";
rotMatAboutY(Rz);
cout<<"Matrix Rz:\n";
for(int i = 0; i<4; i++){
for(int j = 0; j<4; j++){
cout<<Rz[i][j]<<" ";
}
cout<<endl;
}
matrixmult(A, Rz);
break;
}
cout<<"Matrix A:\n";
for(int i = 0; i<4; i++){
for(int j = 0; j<8; j++){
cout<<A[i][j]<<" ";
}
cout<<endl;
}
cleardevice();
delay(100);
sideView(A);
topView(A);
frontView(A);
}
getch();
closegraph();
return 0;
}
最初,无旋转
沿着x、y和z轴旋转20度后
1条答案
按热度按时间lmyy7pcs1#
你需要的是使用法线和背面剔除。所以如果你定义法向量
n
指向你的立方体的每个面比:将决定你是渲染正面还是背面。两个矢量必须在同一个坐标系中,所以要么变换
n
,要么变换相机视图,或者两者都变换......不要忘记,变换矢量是在没有偏移的情况下完成的!!!如果你没有法线
n
,你可以在运行中为每个由点p0,p1,p2,...
定义的面计算它们,如下所示:但是它的方向是指向内还是指向外取决于你的面的缠绕规则,交叉操作数的顺序和你的坐标系约定。所以调整以满足你的需要(渲染法线以直观地查看它们是否正确是个好主意)
再看看这个:
所有的视图和变换都可以用4x 4矩阵来表示,你可以不用投影函数而只用矩阵。