计算机图形学报告太阳系.docx

1、计算机图形学报告太阳系Solar System计算机图形学课程报告文档编号：报告名称：算法设计编 写：天添0743111185编写日期：2009-12指导老师：征审核日期：1 图形学基本原理1.1 几何变换相关原理 几何变换所涉及坐标系和变换如下： Object coordinates-modelview Matrix eye coordinatesprojectionmatrix clip coordinate-perspectiveMatrix normalized device coordinates-viewportTransformation window coordinates.

2、1.2 几何图形像素化相关原理 几何图形像素化就是几何图形经过一系列的变换后储存在帧缓存中的图形，经过其他的一些操作后，将输出到屏幕坐标。由于光栅化前这些坐标是浮点连续的。而屏幕坐标系是离散的像素点。几何图形像素画就是将这些连续的图形变为离散的像素点集并且确保正确的渲染。如果没有该步骤，几何图形是不能渲染到使用像素的屏幕坐标。1.3 光照模型相关原理 点光源：向四面八方发射光线的单点。 聚光灯光源：从特定地点射向特定方向的光源。 平行光源：从无限远处射来的点光源。 面光源：从一个平面发出平行光的光源。 最终的每个像素color输出为几何点在光珊化中环境光，漫反射和镜面光的叠加。2 系统界面与操

3、作说明（阐述如何操作你的软件）2.1 系统界面与图形显示示例 左下角：camera的数据，当前帧率。2.2 操作说明（配合图示进行说明） 在该太阳系中，提供以下几个键： F:摄像机到目标物体的距离增加5。 R:摄像机到目标物体的距离减少5。 :开启或关闭反锯齿。 鼠标左键按住不放时，左右上下移动改变摄像机相对目标物体的位置， 但是不改变距离 选取一个星球，摄像机自动飞向该星球： 未操作前截图： 按F后按R:鼠标左键不放左上拖动选取一个星球，摄像机自动飞向该星球：过程图片：到达图片：目前在反据此模式下，点击“” 关闭反锯齿：这里只有通过帧数的改变显示出改变3 源程序：3.1 太阳系概述： 太阳系

4、 （Solar System）就是我们现在所在的恒星系统。它是以太阳为中心，和所有受到太阳引力约束的天体的集合体：8颗行星冥王星已被开除、至少165颗已知的卫星，和数以亿计的太阳系小天体。这些小天体包括小行星、柯伊伯带的天体、彗星和星际尘埃。广义上，太阳系的领域包括太阳、4颗像地球的行星、由许多小岩石组成的小行星带、4颗充满气体的巨大外行星、充满冰冻小岩石、被称为柯伊伯带的第二个小天体区。在柯伊伯带之外还有黄道离散盘面、太阳圈和依然属于假设的奥尔特云。 模拟太阳系不仅仅是完成作业，也能让我们更近一步的了解太阳系，拓广知识面，对行星的运行有基本的感性和理性感知。增加我们对宇宙探索的渴望，可能培养

5、出又一批天文学家。3.2 数据收集： 要模拟太阳系的运行首先需要他们的数据，这里我收集了“9”大行星，太阳，月球数据如下：行星a(AU)轨道周期(yr)i(deg)e自转周期自转轴倾斜（度）质量（地球）半径(地球)水星0.3870.2417.000.20458.8d00.0550.383金星0.7230.6143.400.007244d-20.8150.949地球1.0001.0000.000.01723.9h23.51.0001.000火星1.521.881.900.09424.6h28.20.1070.533木星5.2011.91.300.0499.9h3.131811.2土星9.5829

6、.42.500.05710.7h26.795.29.45天王星19.283.70.780.04617.2h-82.114.54.01海王星30.1163.71.780.01116.1h28.317.13.88冥王星39.224817.20.2446.39d1220.0020.19103kmd1020kgkm月球38527.318.3-28.6Va0.0557301738 a=轨道半长轴i=轨道倾角e=轨道偏心率1AU=149600000km质量和半径是对于地球的相应值的相对值地球质量=5.977*1027次方地球赤道半径=6378km以上每列数据是每个星球必有的。通过这些属性，以地球为标准，获

7、取行星的运动方程就可以模拟出行星的运动轨迹和位置等数据。现在还需要的外部数据就是：每个行星的纹理。在纹理选择上。我使用的是分辨率为256*256 RBGA 的tga格式文件。使用自己编写的一个类resourceImage 派生于 C4_Image.读取该纹理。以在渲染时使用。在这个程序中，假设所有星球的运动轨迹是圆。而非椭圆。这样就简化了程序在编写时的复杂度。并且假设程序不考虑万有引力。3.3 类设计： 在该太阳系中类图关系如下：3.3.1 Main ClassbaseGameEntity:提供的最主要功能为：为每一个实体提供一个独特的ID.以便为绑定纹理，选择，提供方便，以免冲突。Movin

8、gEntity：提供实体的up,front,speed等参数。Galaxy:是星系的抽象体，他具有个天空球（CSkySphere）和他所包含的星球(star)的集合。Star:是一个星球的抽象实体。他具有表中的属性，除此还有一个CSkySphere。这里CSkySphere作为该星球的球体渲染时的纹理和顶点信息。由于每个星球可能也有他的卫星，他的卫星也可能还有卫星。所以每个star中都有一个star的集合.ChaseCamera：提供围绕一个星球旋转的功能； 提供调整到星球距离的功能；C4_Image:从文件加载32位，24位或8位tga格式的纹理。保存到一个unsigned char 数组中

9、。ResourceImage：对C4_Image进行扩展。可以加载DDS_DX1_decals_compressed,DDS_DX1_CUBMAP两类文件。但是由于DDS_DX1_CUBMAP 显示分辨率太低，所以就没有用cubmap.这使得没有使用cubmap而是使用自己写的天空球来模拟环境纹理。3.3.2 辅助类：类图如下：这6个类为辅助类，每一个需要渲染的物体都有一个C4_Material.matrix3D matrix4D,vector3D,vector4D 提供基本的几何运算，不用多说。C4_Light 提供创建8个光源实例。在程序中使用了1个光源。3.3.3 模块增强类： 所谓模块

10、增强类，就是提供更好的图形显示或更多的功能。类图如下：该太阳系提供3个增强模块，3个功能：1），C4_Selector :对派生于baseGameEntity的物体提供鼠标选择功能。2），C4_Antialiasing :提供全局的antialiasing,反锯齿（反走样），使用的是glaccum.3), C4_CgEffect :提供高级顶点，片段编程。基于NVADIA的Cg2.0.（opengl version:=1.4）以上为对程序的一些类的简介。不过大体就这15个类。（2，3运行截图请参看附录）3.4 功能设计：3.4.1 反锯齿： 在太阳系中使用了2组jitter系数，让projec

11、tion matrix产生微小位移，然后渲染物体。将两次渲染的结果累计到缓存区，然后渲染出来。使用前后对比： 前 后 使用的样本数据： double j22=0.25,0.75,0.75,0.25; （参考the red book 累积缓冲区）在使用glCalllist 调用 一个非常冗余的 圆的渲染时 在开启反锯齿后 在同样的情况下，每秒的帧数只有6帧/秒，而改进之后 帧数达到了9帧/秒左右。 在原来情况下每个顶点 要被渲染4次。改进后每个点被渲染2次。如果要进一步的该进，就要使用 索引。由于时间关系就不再优化。以下为改进前和改进后的性能对比。 改进前：改进后对比。由该表可以看出glAccu

12、m所使用的时间之巨大，也可以让我们了解到渲染时顶点越少越好。3.4.2 帧率 但是如果不开启 反锯齿，改进前和改进后的区别基本上没有。都是20帧左右。但是，不过为什么无论我怎么调，该帧数都不变。例如将太阳系中的行星删除掉80，禁用光源，但是帧数依然不变。最终发现 是glutTimerFunction中设置问题。开始设置的1/24.0。而后当我设置成1/40.0时。帧数变为32帧左右。当我继续减少该值的时候。无反锯齿下，最多可以达到66帧，反锯齿下反而有所降低到89帧。但是此时又有一个问题。当帧数载66帧左右的时候。有某个瞬间尤其不稳定可能突然渲染时间降低到0.0001以下或突然上升到0.1s。

13、此时.就需要一个采样的队列对时间进行采样使得物体和镜头的移动更加平滑。在程序中我使用了一个最简单的方法，求最近n次的平均帧长。这里应该设置一个比率。越早的数据权重应该越低。以下为部分代码。#define FRAMESAMPLESIZE 100int curframeSample=0;float frameFRAMESAMPLESIZE=0; curframeSample=(curframeSample+1)%FRAMESAMPLESIZE; framecurframeSample=time_eclipse; updateTime=0; for(int i=0;im_pAtmosphere-Re

14、nder(); glDisable(GL_ALPHA_TEST); if(!bEnable) glDisable(GL_BLEND); if(bCullEnable) glEnable(GL_CULL_FACE); glDepthMask(true); 截图如下：从左到右分别为：太阳，月球，金星，地球。由于没找到好的云图片，效果就只有这样了。不过这依然还需要改进。3.4.4 太阳特效 太阳应该看起来有一点流动效果，并且太阳的周围应该还有飞出的火焰，也就是类似于光晕的光圈。这就有2个问题。1，流动，2，光圈（glow）。3.4.4.1 流动 流动如果在常规的管道流程中来生成，对太阳的每个顶点进行

15、变换或纹理。是比较耗费时间的。所以这里在程序中使用了NVIDIA的cg（支持opengl1.4及其以上）。我的电脑最多支持opengl1.4,所以就只有用cg所支持的fragment shader和 vertex shader.这里只是对顶点进行放缩。 float displacement = scaleFactor * 0.5 *sin(position.y * frequency * time)+1; float4 displacementDirection = float4(normal.x, normal.y, normal.z, 0); float4 newPosition = po

16、sition + displacement * displacementDirection; oPosition = mul(modelViewProj, newPosition); 这里，对于每个顶点，都会在其法向量方向进行震动。具有相同y值的顶点会有相同的幅度，在法线方向的长度相同。这样又y值的不一样，在视觉上就有起伏感觉。而在modelviewproj构造时会有一个根据时间增加的旋转值。这样看起来每个顶点都在起伏并且移动。在贴了纹理后就有点流体的感觉。 这是在使用cg时的太阳 这是在没有使用时的图片但是太阳的反锯齿由于cg的特性（在渲染cg前的所有其他非cg的函数操作都对其没有作用）而变

17、没了。3.4.4.2 Glow: 目前我在opengl下还没弄出来不过班上有同学在dx下已经弄出来了。3.4.5 光属性 如果你查看附录的图像时，你会发现及时是很远的行星向光面都很亮，因为光的衰减函数如下： F(x)=1/(1+0.001x);Y2.0 x 0 1000AU3.5 代码实现：只将关键代码罗列如下：3.5.1 Star 类：/* Star.h*/#pragma once#include movingentity.h#include SkySphere.h#include C4_Material.h#include /* a star has a skysphere to show

18、 its outside pattern*/class Star : public MovingEntityprivate: CSkySphere *m_sSphere; C4_Material m_c4Mat; std:vector *m_pvSatellite; char m_szStarName10; CSkySphere *m_pAtmosphere;public: Star(int id, float inMass,float inRadiu, float inAU,float inYr,float inDeg,float inE, const char*texFileName, f

19、loat inRPeriod=0,float inRotInc=0, const char*starName=NONE); Star(void);public: /those method should be implemented virtual void Update(float ftime_eclipse); virtual void Render(); virtual bool HandleMessage(const Telegram& msg); CSkySphere* getSkySphere()constreturn m_sSphere; void setMaterial(con

20、st C4_Material& m)m_c4Mat=m; void addSatellite(Star*inStar)if(inStar)m_pvSatellite-push_back(inStar); Star* getSatelliteStar(int id); void setChildEnvMode(int model);protected: float m_fAU;/轨道半长轴以AU为单位1AU=149600000km float m_fYr;/轨道周期以地球的周期比。地球：.0 float m_fDeg;/轨道倾斜角deg float m_fE;/偏心率shaft eccentri

21、city float m_fRotationPeriod;/自转周期 float m_fRotInc;/rotational inclination 自转轴倾角 /mass /radius will be handle down in the base. /simple rotation float m_fTheta; float m_fSelfRotationTheta; void RenderSatellite();public: float getAU()constreturn m_fAU; void setAU(const float& inAU)m_fAU=inAU; float g

22、etYr()constreturn m_fYr; void setYr(const float& inYr)m_fYr=inYr; float getDeg()constreturn m_fDeg; void setDeg(const float& inDeg)m_fDeg=inDeg; float getE()constreturn m_fE; void setE(const float& inE)m_fE=inE; float getRotationPeriod()constreturn m_fRotationPeriod; void setRotationPeriod(const flo

23、at& inRotationPeriod)m_fRotationPeriod=inRotationPeriod; float getRotInc()constreturn m_fRotInc; void setRotInc(const float& inRotInc)m_fRotInc=inRotInc;/* Star.cpp*/*#pragma warning (disable:4996)#include Star.h#include utility.h#include Matrix4D.h#include ChaseCamera.h#include #include C4_math.h#i

24、nclude C4_Antialiasing.h#ifdef CG_ENABLEstatic C4_CgEffect * GlowEffect=NULL;#endiffloat getDegJump(const Vector3D&v) double length=v.lengthSq(); if(length=2500) return 6; else if(length200*200) return 12; else if(length500*500) return 12; else if(lengthgetSkySphere()-Render();#endifStar:Star(int id

25、, float inMass,float inRadiu, float inAU,float inYr,float inDeg,float inE, const char*texFileName, float inRPeriod,float inRotInc, const char*starName) :MovingEntity(id) m_fSelfRotationTheta=0; this-setMass(inMass); this-setRadius(inRadiu); this-setAU(inAU); this-setYr(inYr); this-setE(inE); this-setDeg(inDeg); this-setRotationPeriod(inRPeriod); this-setRotInc(inRotInc); m_sSphere=new CSkySphere(); m_sSphere-loadSkyImag