【GAMES101】作业1（提高）与框架理解
作业描述
本次作业的任务是填写一个旋转矩阵和一个透视投影矩阵。给定三维下三个点 v0(2.0, 0.0, −2.0), v1(0.0, 2.0, −2.0), v2(−2.0, 0.0, −2.0), 你需要将这三个点的坐标变换为屏幕坐标并在屏幕上绘制出对应的线框三角形 (在代码框架中，我们已经提供了 draw_triangle 函数，所以你只需要去构建变换矩阵即可)。简而言之，我们需要进行模型、视图、投影、视口等变换来将三角形显示在屏幕上。在提供的代码框架中，我们留下了模型变换和投影变换的部分给你去完成。
以下是你需要在 main.cpp 中修改的函数（请不要修改任何的函数名和其他已经填写好的函数，并保证提交的代码是已经完成且能运行的）：
• get_model_matrix(float rotation_angle):
逐个元素地构建模型变换矩阵并返回该矩阵。在此函数中，你只需要实现三维中绕 z 轴旋转的变换矩阵，而不用处理平移与缩放。• get_projection_matrix(float eye_fov, float aspect_ratio, floatzNear, float zFar): 使用给定的参数逐个元素地构建透视投影矩阵并返回该矩阵。
• [Optional] main(): 自行补充你所需的其他操作。
当你在上述函数中正确地构建了模型与投影矩阵，光栅化器会创建一个窗口显示出线框三角形。由于光栅化器是逐帧渲染与绘制的，所以你可以使用 A 和 D 键去将该三角形绕 z 轴旋转(此处有一项提高作业，将三角形绕任意过原点的轴旋转)。当你按下 Esc 键时，窗口会关闭且程序终止。
另外，你也可以从命令行中运行该程序。你可以使用以下命令来运行和传递旋转角给程序，在这样的运行方式下，是不会生成任何的窗口，输出的结果图像会被存储在给定的文件中 (若未指定文件名，则默认存储在 output.png 中)。图像的存储位置在可执行文件旁，所以如果你的可执行文件是在 build 文件夹中，那么图像也会在该文件夹内。
命令行的使用命令如下：
./Rasterizer //循环运行程序，创建一个窗口显示，且你可以使 用A键 和D键 旋 转 三 角 形。
./Rasterizer −r 20 //运行程序并将三角形旋转20度,然后将结果存在output.png中
./Rasterizer −r 20 image.png //运行程序并将三角形旋转20度，然后将结果存在image.png中。
解
Eigen::Matrix4f get_model_matrix(float rotation_angle)//模型变换矩阵
Eigen::Matrix4f get_model_matrix(float rotation_angle)//模型变换矩阵
{
Eigen::Matrix4f model = Eigen::Matrix4f::Identity();
// TODO: Implement this function
// Create the model matrix for rotating the triangle around the Z axis.
// Then return it.
Eigen::Matrix4f rotation;
double fangle = rotation_angle / 180 * MY_PI;//角度转弧度，便于计算
rotation << cos(fangle), -sin(fangle), 0, 0,
sin(fangle), cos(fangle), 0, 0,
0, 0, 1, 0,
0, 0, 0, 1;//模型旋转矩阵（绕z轴）
model = rotation * model;
return model;
}
Eigen::Matrix4f get_projection_matrix(float eye_fov, float aspect_ratio, float zNear, float zFar)//投影变换矩阵
Eigen::Matrix4f get_projection_matrix(float eye_fov, float aspect_ratio,
float zNear, float zFar)//投影变换矩阵
{
// Students will implement this function
Eigen::Matrix4f projection = Eigen::Matrix4f::Identity();
// TODO: Implement this function
// Create the projection matrix for the given parameters.
// Then return it.
Eigen::Matrix4f proj, ortho;
proj << zNear, 0, 0, 0,
0, zNear, 0, 0,
0, 0, zNear + zFar, -zNear * zFar,
0, 0, 1, 0;//透视投影矩阵
double w, h, z;
h = zNear * tan(eye_fov / 2) * 2;
w = h * aspect_ratio;
z = zFar - zNear;
ortho << 2 / w, 0, 0, 0,
0, 2 / h, 0, 0,
0, 0, 2 / z, -(zFar+zNear) / 2,
0, 0, 0, 1;//正交投影矩阵，因为在观测投影时x0y平面视角默认是中心，所以这里的正交投影就不用平移x和y了
projection = ortho * proj * projection;
return projection;
}
这次作业本身难度并不大，更多是用来温习MVP的流程和公式，具体要点的都在代码注释里写出了，这里稍微多讲一些我对提高部分的实现，首先看到提高部分的要求：在 main.cpp 中构造一个函数，该函数的作用是得到绕任意
过原点的轴的旋转变换矩阵。Eigen::Matrix4f get_rotation(Vector3f axis, float angle)
这就很明显是lecture 04里提到过的罗德里格斯旋转公式了：
根据给定的参数写出公式即可，这里要注意齐次坐标在这个公式应用时第四维的计算结果可能会影响最终成像，要记得改过来。
Eigen::Matrix4f get_rotation(Vector3f axis, float angle) //任意轴旋转矩阵
Eigen::Matrix4f get_rotation(Vector3f axis, float angle) {//任意轴旋转矩阵（罗德里格斯旋转公式，默认轴过原点）
double fangle = angle / 180 * MY_PI;
Eigen::Matrix4f I, N, Rod;
Eigen::Vector4f axi;
Eigen::RowVector4f taxi;
axi << axis.x(), axis.y(), axis.z(), 0;
taxi << axis.x(), axis.y(), axis.z(), 0;
I << 1, 0, 0, 0,
0, 1, 0, 0,
0, 0, 1, 0,
0, 0, 0, 1;
N << 0, -axis.z(), axis.y(), 0,
axis.z(), 0, -axis.x(), 0,
-axis.y(), axis.x(), 0, 0,
0, 0, 0, 1;
Rod = cos(fangle) * I + (1 - cos(fangle)) * axi * taxi + sin(fangle) * N;
Rod(3, 3) = 1;//这里要注意，非齐次坐标的公式应用在齐次坐标上时记得运算完成后把矩阵的右下角改为1，否则会导致图形比例错误
return Rod;
}
但是现在我有了函数输出的矩阵，应该怎么操作才能让它真正的应用到图像上呢？
当然就是加入一个类似mvp矩阵的传递接口，在main.cpp文件中我们可以看到类似set_projection这种函数，他的作用是将内部的投影矩阵设为给定矩阵 p，并传递给光栅化器，同样的我们也可以在rasterizer.cpp加入接口函数
void rst::rasterizer::set_rodrigues(const Eigen::Matrix4f& r)
{
rodrigues = r;
}
接下来在draw函数里加入mvp矩阵的计算：
Eigen::Matrix4f mvp = projection * view * model * rodrigues;//计算mvp矩阵，因为绕任意轴旋转也是模型变换的一种，所以放在model矩阵的相邻位置
注意这步之前要在rasterizer.hpp里声明接口函数和矩阵：
class rasterizer
{
...
public:
void set_rodrigues(const Eigen::Matrix4f& r);
private:
Eigen::Matrix4f rodrigues;
...
};
最后在main.cpp的main函数里面加入旋转轴和角度的参数输入：
std::cout << "Please enter the axis and angle:" << std::endl;
std::cin >> raxis.x() >> raxis.y() >> raxis.z() >> ra;//定义罗德里格斯旋转轴和角
加入按下键盘’r’键就绕指定轴旋转的判定：
if (rflag) //如果按下r了，就开始绕给定任意轴旋转
r.set_rodrigues(get_rotation(raxis, rangle));
else
r.set_rodrigues(get_rotation({ 0,0,1 }, 0));
...
else if (key == 'r') {//按下r，再次绕给定任意轴旋转
rflag = true;
rangle += ra;
}
效果
我们输入时设定每按一次’R’就绕z轴（0，0，1）旋转90度：
这是初始的图像：
这是按下三次‘A’(即绕z轴逆时针旋转30°)时的图像：
这是再按下四次‘D’(即顺时针旋转10°)时的图像：
这是再按下一次’R’（即逆时针旋转80°）时的图像：
可以自己尝试，比如绕（1，2，3）旋转10°的图像：
框架理解
这份作业框架包含了许多功能，提前了解的话有助于我们后续的学习，为方便c++基础一般的朋友理解，给大部分代码做了注释：
main.cpp
int main(int argc, const char** argv)
{
float angle = 0;//定义角度
bool command_line = false;//定义命令行开关标志，默认为关
std::string filename = "output.png";//定义文件名，默认为output.png"
Eigen::Vector3f raxis(0, 0, 1);
double rangle = 0, ra;
if (argc >= 3) {//接收到的参数大于三个，即检测到通过命令行传入参数时
command_line = true;//设命令行开关标志为开
angle = std::stof(argv[2]); //从命令行获取角度参数
if (argc == 4) {//接收到的参数为四个，那么说明命令行输入了文件名参数
filename = std::string(argv[3]);//从命令行获取文件名
}
}
rst::rasterizer r(700, 700);//设定700*700像素的光栅器视口
Eigen::Vector3f eye_pos = { 0, 0, 5 };//设定相机位置
std::vector<Eigen::Vector3f> pos{ {2, 0, -2}, {0, 2, -2}, {-2, 0, -2} };//设定三顶点位置
std::vector<Eigen::Vector3i> ind{ {0, 1, 2} };//设定三顶点序号,用于画图时确定需要处理几个顶点，这里表示的是三个顶点
auto pos_id = r.load_positions(pos);
auto ind_id = r.load_indices(ind);//保存多个图形的顶点和序号，本次作业只涉及一个图形，可以不管
int key = 0;//键盘输入
int frame_count = 0;//帧序号
if (command_line) {//如果命令行开关标志为开（这一段if代码是为了应用命令行传入的参数，比如初始角度和文件名）
r.clear(rst::Buffers::Color | rst::Buffers::Depth);//初始化帧缓存和深度缓存（本次作业本次作业只涉及一个图形，所以不涉及深度，可以不管）
r.set_model(get_model_matrix(angle));
r.set_view(get_view_matrix(eye_pos));
r.set_projection(get_projection_matrix(45, 1, 0.1, 50));//向光栅器传入MVP矩阵
r.set_rodrigues(get_rotation(raxis, rangle));
r.draw(pos_id, ind_id, rst::Primitive::Triangle);//开始画图
cv::Mat image(700, 700, CV_32FC3, r.frame_buffer().data());
image.convertTo(image, CV_8UC3, 1.0f);
cv::imwrite(filename, image);//写入文件名
return 0;
}
bool rflag = false;
std::cout << "Please enter the axis and angle:" << std::endl;
std::cin >> raxis.x() >> raxis.y() >> raxis.z() >> ra;//定义罗德里格斯旋转轴和角
while (key != 27) {//只要没有检测到按下ESC就循环(ESC的ASCII码是27)
r.clear(rst::Buffers::Color | rst::Buffers::Depth);
r.set_model(get_model_matrix(angle));
r.set_view(get_view_matrix(eye_pos));
r.set_projection(get_projection_matrix(45, 1, 0.1, 50));
if (rflag) //如果按下r了，就开始绕给定任意轴旋转
r.set_rodrigues(get_rotation(raxis, rangle));
else
r.set_rodrigues(get_rotation({ 0,0,1 }, 0));
r.draw(pos_id, ind_id, rst::Primitive::Triangle);
cv::Mat image(700, 700, CV_32FC3, r.frame_buffer().data());
image.convertTo(image, CV_8UC3, 1.0f);
cv::imshow("image", image);//显示图像
key = cv::waitKey(10);//等待10号码接收键盘输入，没有输入就为空，图像不做调整，保持原状
std::cout << "frame count: " << frame_count++ << '\n';//显示当前是第几帧画面
if (key == 'a') {//按下a，逆时针旋转10°
angle += 10;
}
else if (key == 'd') {//按下d，顺时针旋转10°
angle -= 10;
}
else if (key == 'r') {//按下r，绕给定任意轴旋转
rflag = true;
rangle += ra;
}
}
return 0;
}
rasterization.cpp
#include 附件
中点画线算法
中点算法会在后面的课程提到，想提前了解的朋友可以看下这个ppt，看看中点算法是怎么推导的：
直线扫描转换-中点算法
源代码
附上源代码，网页看不清或者想要自己运行尝试修改的朋友可以自行下载：
CSDN:【GAMES101】作业1（提高）
GITHUB:【GAMES101】作业合集