地形渲染的动态LOD四叉树算法详细实现

最新推荐文章于 2024-08-10 07:08:38 发布
最新推荐文章于 2024-08-10 07:08:38 发布
阅读量4.6k 收藏 12
点赞数 2
分类专栏： OpenGL 文章标签： GIS LOD OpenGL 三维地形
OpenGL 同时被 2 个专栏收录
21 篇文章 0 订阅
订阅专栏
GIS
6 篇文章 0 订阅
订阅专栏

地形渲染的动态LOD四叉树算法详细实现

                                                                                    作者：野草 email:JerryLi@sina.com

声明：请将本文档和程序配合使用，旨在使读者费很少的时间和精力就能理解此算法。读者应该熟悉递归程序设计,以及基本的VC OpenGL编程.

渲染地形时，如果采用固定地形精度等级，则会随着地形面积的增大，渲染的三角形数量会随着面积的增大而增长。采用动态精度等级，越远处的地形画的三角形面积越大(精度越低)，越近的地形面积越小(精度越高)，坡度高的地方精度高，平坦的地方精度低，这样能有效地减少画三角形的数量，提高游戏引擎的速度.

采用四叉树算法，可有效把地形不在视野内的网格去掉,降低程序的时间复杂度

其基本思路是：先把地形一分为四，用递归的方法对每个网格渲染。对每个网格，如果达到最高精度，则退出；如果不在视野内，也退出。再对符合条件的网格递归下去。

渲染网格的方法是用GL_TRIANGLE_FAN方法，如下图

但是这样在相邻精度等级相差1倍或1倍以上时，渲染的地形会出现裂缝。

 

解决方法见下图：

 

这样相邻的网格可以相差任意级精度

参见下面的截图：看中间的那个大正方形。

 

 

 

                  

 

程序如下：

 

 

 

void CTerrain::Render()

{

    ………

//细分标记数组清0（这里为了节约内存，一个网格用一个bit表示）

m_cFlag.Reset();  

//在细分标记数组记录需要细分的网格

Update((m_nMapSize-1)/2,(m_nMapSize-1)/2,(m_nMapSize-1)/2,1);

        //根据细分标记数组递归渲染网格

        RenderQuad((m_nMapSize-1)/2,(m_nMapSize-1)/2,(m_nMapSize-1)/2);

        …….    

}

 

//渲染网格的递归程序

//参数:

//                 nXCenter: 网格中心点坐的X值

//                 nZCenter: 网格中心点坐的Y值

//                 nSize:网格大小

//返回:渲染的三角形数量

//使用前提:m_cFlag需要细分数组要先调用Update函数记录好.

int CTerrain::RenderQuad(int nXCenter,int nZCenter,int nSize)

{

                   int nCount = 0;

                   int nH0 = Height(nXCenter,nZCenter);

                   int nH1 = Height(nXCenter-nSize,nZCenter-nSize);

                   int nH2 = Height(nXCenter+nSize,nZCenter-nSize);

                   int nH3 = Height(nXCenter+nSize,nZCenter+nSize);

                   int nH4 = Height(nXCenter-nSize,nZCenter+nSize);

                   //求网格5个点中高度最大和最小值

                   int nMax = nH0,nMin = nH0;

                   if(nMax<nH1)nMax = nH1;

                   if(nMax<nH2)nMax = nH2;

                   if(nMax<nH3)nMax = nH3;

                   if(nMax<nH4)nMax = nH4;

                   if(nMin>nH1)nMin = nH1;

                   if(nMin>nH2)nMin = nH2;

                   if(nMin>nH3)nMin = nH3;

                   if(nMin>nH4)nMin = nH4;

                   SIZE size;

                   size.cx = nSize;size.cy=(nMax-nMin)/2;

                   //如果此网格不在视野内，返回

                   if(!m_pFrustum->CubeInFrustum(nXCenter,(nMax+nMin)/2,nZCenter,max(size.cx,size.cy)))

                     return 0;

                   if(m_cFlag.IsTrue(nXCenter,nZCenter))//需要细分,递归渲染4个子网格

                   {

                     nCount += RenderQuad(nXCenter-nSize/2,nZCenter-nSize/2,nSize/2);

                     nCount += RenderQuad(nXCenter+nSize/2,nZCenter-nSize/2,nSize/2);

                     nCount += RenderQuad(nXCenter+nSize/2,nZCenter+nSize/2,nSize/2);

                     nCount += RenderQuad(nXCenter-nSize/2,nZCenter+nSize/2,nSize/2);

                   }

                   else//渲染此网格

                   {

                     glBegin( GL_TRIANGLE_FAN );

                     SetTextureCoord(nXCenter,nZCenter);  // center

                     glVertex3i(nXCenter, nH0, nZCenter);

 

                     SetTextureCoord(nXCenter-nSize,nZCenter-nSize); //left top

                     glVertex3i(nXCenter-nSize,nH1 , nZCenter-nSize);

                     //为了避免出现裂缝而多画三角形扇（就是多加那些红线啦）

                     RemedyTop(nXCenter,nZCenter,nSize);

                    

                     SetTextureCoord(nXCenter+nSize,nZCenter-nSize); //right top

                     glVertex3i(nXCenter+nSize, nH2, nZCenter-nSize);

                     //为了避免出现裂缝而多画三角形扇（就是多加那些红线啦）

                     RemedyRight(nXCenter,nZCenter,nSize);

 

                     SetTextureCoord(nXCenter+nSize,nZCenter+nSize); //right bottom

                     glVertex3i(nXCenter+nSize, nH3, nZCenter+nSize);

                     //为了避免出现裂缝而多画三角形扇（就是多加那些红线啦）

                     RemedyBottom(nXCenter,nZCenter,nSize);

 

                     SetTextureCoord(nXCenter-nSize,nZCenter+nSize);//left bottom

                     glVertex3i(nXCenter-nSize, nH4, nZCenter+nSize);

                     //为了避免出现裂缝而多画三角形扇（就是多加那些红线啦）

                      RemedyLeft(nXCenter,nZCenter,nSize);

 

                     SetTextureCoord(nXCenter-nSize,nZCenter-nSize); //left top

                     glVertex3i(nXCenter-nSize, nH1, nZCenter-nSize);

                     glEnd();

                     nCount += 4;//这里为求简单，计算渲染的三角形数目并不是十分准确,比实际要小些

                   }

                   return nCount;

}

 

 

 

 

//更新细分记录数组的递归函数

//参数:

//                 nXCenter: 网格中心点坐的X值

//                 nZCenter: 网格中心点坐的Y值

//                 nSize:网格大小

//  nLevel:当前的细分等级(其实可由nSize计算出来)    

void CTerrain::Update(int nXCenter,int nZCenter,int nSize,int nLevel)

{

                   CVector3 vPos = m_pCamera->GetPos();

//                 int nSize = (m_nMapSize-1)/pow(2,nLevel);//由细节级度得到当前渲染矩形的大小

                   //求网格5个点中高度最大和最小值

                   int nMax = Height(nXCenter,nZCenter);

                   int nMin = nMax;

                   int nH1 = Height(nXCenter-nSize,nZCenter-nSize);

                   int nH2 = Height(nXCenter+nSize,nZCenter-nSize);

                   int nH3 = Height(nXCenter+nSize,nZCenter+nSize);

                   int nH4 = Height(nXCenter-nSize,nZCenter+nSize);

                   if(nMax<nH1)nMax = nH1;

                   if(nMax<nH2)nMax = nH2;

                   if(nMax<nH3)nMax = nH3;

                   if(nMax<nH4)nMax = nH4;

                   if(nMin>nH1)nMin = nH1;

                   if(nMin>nH2)nMin = nH2;

                   if(nMin>nH3)nMin = nH3;

                   if(nMin>nH4)nMin = nH4;

                   //网格的中心点坐标

                   CVector3 vDst(nXCenter,(nMax+nMin)/2,nZCenter);

                   SIZE size;

                   size.cx = nSize;size.cy=(nMax-nMin)/2;

                   //如果此网格不在视野内，返回

                   if(!m_pFrustum->CubeInFrustum(nXCenter,(nMax+nMin)/2,nZCenter,max(size.cx,size.cy)))

                     return ;

                   //视点到网格中心点的距离

                   int nDist = Dist(vPos,vDst);

                   if(nDist>2000)//距离大于2000时才将距离和误差(坡度)2个因素考虑

                   {

                     //与距离成正比,与误差(坡度)成反比,还与当前细分等级有关

                     if(70.0*Dist(vPos,vDst)/nSize/CalcError(nXCenter,nZCenter,nSize)>nLevel)

                            return;//此网格不需要细分

                   }

                   else//近距离(<2000)时只考虑误差(坡度)因素(因为动态LOD会造成地表呼吸现象,为减弱此现象不考虑距离因素)

                   {

                     if(CalcError(nXCenter,nZCenter,nSize)<30)

                            return;//此网格不需要细分

                   }

                   if(nSize!=MIN_GRID)

                   { 此网格需要细分

                     m_cFlag.Set(nXCenter,nZCenter,true);//把细分数组相应位设为true

                     //递归4个子网格

                     Update(nXCenter-nSize/2,nZCenter-nSize/2,nSize/2,nLevel+1);

                     Update(nXCenter+nSize/2,nZCenter-nSize/2,nSize/2,nLevel+1);

                     Update(nXCenter+nSize/2,nZCenter+nSize/2,nSize/2,nLevel+1);

                     Update(nXCenter-nSize/2,nZCenter+nSize/2,nSize/2,nLevel+1);

                   }

                   //此网格不需要细分

}

 

//计算误差(坡度)的函数

//参数:

//                 nXCenter: 网格中心点坐的X值

//                 nZCenter: 网格中心点坐的Y值

//                 nSize:网格大小

//返回: 误差大小

//说明:只是计算某些点细分与不细分的高度值差的和

inline int CTerrain::CalcError(int nXCenter,int nZCenter,int nSize)

{

                   int nError = 0;

                   int nH0 = Height(nXCenter,nZCenter);

                   int nH1 = Height(nXCenter-nSize,nZCenter-nSize);

                   int nH2 = Height(nXCenter+nSize,nZCenter-nSize);

                   int nH3 = Height(nXCenter+nSize,nZCenter+nSize);

                   int nH4 = Height(nXCenter-nSize,nZCenter+nSize);

                   nError += abs(Height(nXCenter,nZCenter-nSize)-(nH1+nH2)/2);

                   nError += abs(Height(nXCenter+nSize,nZCenter)-(nH2+nH3)/2);

                   nError += abs(Height(nXCenter,nZCenter+nSize)-(nH3+nH4)/2);

                   nError += abs(Height(nXCenter-nSize,nZCenter)-(nH4+nH1)/2);

 

                   nError +=abs(Height(nXCenter-nSize/2,nZCenter-nSize/2)-(nH0+nH1)/2);

                   nError +=abs(Height(nXCenter+nSize/2,nZCenter-nSize/2)-(nH0+nH2)/2);

                   nError +=abs(Height(nXCenter+nSize/2,nZCenter+nSize/2)-(nH0+nH3)/2);

                   nError +=abs(Height(nXCenter-nSize/2,nZCenter+nSize/2)-(nH0+nH4)/2);

                   return nError;

}

如下图:

当不细分时红点的高度值是它所在线段的两个端点的平均值,细分时就是它实际的高度值,这2个值是不一样的. 两者这间的高度差就是误差,把所有红点的误差相加就是此网格的误差

 

//修补网格顶部裂缝的函数

//参数:

//                 nXCenter: 网格中心点坐的X值

//                 nZCenter: 网格中心点坐的Y值

//                 nSize:网格大小

void CTerrain::RemedyTop(int nXCenter,int nZCenter,int nSize)

{

                   if(nZCenter-2*nSize>=0)//保证寻址不越界

                   {

                     if(!m_cFlag.IsTrue(nXCenter,nZCenter-2*nSize))//此网格相邻的顶部网格(同级的)没有细分

                            return;

                   }

                   else

                     return;

                   //此网格相邻的顶部网格(同级的)已经细分

                   //对顶部一分为二递归(跨级精度的情况)

                   RemedyTop(nXCenter-nSize/2,nZCenter-nSize/2,nSize/2);

                   SetTextureCoord(nXCenter,nZCenter-nSize);

                   glVertex3i(nXCenter,Height(nXCenter,nZCenter-nSize),nZCenter-nSize);

                   RemedyTop(nXCenter+nSize/2,nZCenter-nSize/2,nSize/2);

}

其它的修补裂缝的函数也是一样的流程.

 

 

 

近处没有考虑距离因素,只考虑误差(坡度)因素

 

距离越运,精度越低

 

感谢你的阅读,如果你有什么更好的思路或更好的算法(比如能消除地表呼吸现象),请和我联系:JerryLi@sina.com

微岩
关注
专栏目录
基于四叉树（QuadTree）的LOD地形实现
Gypsyy的专栏
08-20 1万+
实现基于四叉树的LOD地形时，我遇到的主要问题是如何修补地形裂缝。        本段我将描述使用LOD地形的优势，我实现LOD地形的思路，实现LOD地形核心模块的详细过程，以及修补地形裂缝的思路。        首先，LOD地形与一般地形不同：一般的地形是这样实现的：整个地形是一个三角形网格，一个513*513的地形包括513*513个顶点，512*512*2个三角形；在开始渲染地形之前，写
地形渲染的动态LOD四叉树算法
12-10
地形渲染的动态LOD四叉树算法，读者应该熟悉递归程序设计,以及基本的VC OpenGL编程.
Quadtree-js 项目使用教程
最新发布
gitblog_00795的博客
08-10 255
Quadtree-js 项目使用教程 quadtree-jsA lightweight quadtree implementation for javascript项目地址:https://gitcode.com/gh_mirrors/qu/quadtree-js 项目介绍 Quadtree-js 是一个轻量级的 JavaScript 四叉树实现，适用于高效的二维空间数据管理，特别是在处理大量对...
回顾四叉树LOD地形（上）
lgcacddn的专栏
05-27 7992
唉！~其实这是在差不多一年前实现的东西，但当时没作好记录。放了那么久了，如果不做点总结的话，好像有点对不起自己，于是·········还是做点什么吧。        本人不才，脑洞比较小，只能大量参考“潘李亮”大神的四叉树LOD地形论文和一个叫做“曾涛”大神的代码实现，小弟在此先向诸位致敬了。下面分享本人实践之后的一点心得。 为什么要用LOD地形？        也许看完了传说中的“
四叉树的LOD
qq960885333的专栏
11-15 1180
创建输出多边形的索引.这里首先从根节点开始判断照相机在哪个节点里面，通过距离来判断当前节点是否可以直接输出三角形，否则的话判断当前节点的四个孩子节点是否可以直接输出三角形； 如果当前节点可以直接输入三角形了，这个时候如果当前节点是最小节点，则直接输出三角形，然后退出，如果当前节点不是最小节点看看上下左右四个节点哪个是可以直接输出三角形的，如果四个邻居节点都可以直接输出三角形，没的说了当前节点和四
三维地形渲染：四叉树递归LOD算法详解
"这篇文档介绍了如何使用四叉树算法来实现三维地形的LOD（Level of Detail）渲染，以优化大规模地形数据的显示性能。在地形渲染中，四叉树是一种常用的分治策略，用于根据观察距离动态调整地形的细节级别，减少计算...
基于四叉树的动态LOD虚拟地形优化 (2009年)
06-12
提出了一种基于四叉树的动态LOD优化方法。...并在 LOD简化过程中根据节点分割和渲染的规律,提出合适的裂缝消除算法,去掉四叉树中不合理的分割,最后使得分割后的节点达到最优以完成四叉树网格优化的目的。
采用VC及Opengl基于四叉树及LOD算法的室外三维地形实时渲染源代码
02-17
2. LOD算法的实现：根据摄像机位置和四叉树节点的距离动态选择合适的细节层次。 3. 地形的加载和渲染：读取地形数据，使用OpenGL进行高效的渲染。 4. 摄像机控制：允许用户通过键盘或鼠标控制摄像机移动，从而改变...
qslim算法，LOD的实现
12-04
在网上找了好久了，是该算法创始人上传的资料，如果文件里有缺少的东西，麻烦大家说一声，我会找到需要的文件继续上传~~
OpenGL实现的四叉树地形LOD源码（含论文说明）
12-27
用OpenGL实现的四叉树地形LOD源码和论文说明，具体的介绍可以看https://blog.csdn.net/qq_31709249/article/details/93659972，效果视频可以看链接：https://pan.baidu.com/s/1tuwCTWVJUSqxEYJFx363vQ 提取码：jy3p。下载后环境配置问题，可直接发邮箱1126318649@qq.com咨询
回想四叉树LOD地形（上）
weixin_33896726的博客
05-31 187
       唉。~事实上这是在差点儿相同一年前实现的东西，但当时没作好记录。放了那么久了，假设不做点总结的话，好像有点对不起自己，于是·········还是做点什么吧。        我脑洞比較小，仅仅能大量參考“潘李亮”大神的四叉树LOD地形论文和一个叫做“曾涛”大神的代码实现，小弟在此先向诸位致敬了。以下分享本人实践之后的一点心得。 为什么要用LOD地形？       ...
四叉树lod结合灯塔AOI
213的博客
05-25 2348
灯塔相关结构体 1：灯塔区域结构 struct towerSpace_s { void (*callback)(void*pUserData,bool bAddTo,uint64_t watcher, uint64_t marker); -- 回调函数 void* pUserData; //用户信息 float fMin[2]; //最小位置 float fGridLength[3];//网格x y,z方向长度 float fMovefRange;//移动范围 int32_t
最新-基于四叉树的LOD地形设计
weixin_34254823的博客
03-17 222
  上次说过要在这里把最新的地形设计写出来，不过一直没时间写，也想画些UML图，但没工具画，呵呵，只能通过手写了。希望大家看得明白吧。        基本设计思想：        主要有以下几个类，        CTerrainBlock，地形一个block，每个对象拥有自己的顶点缓冲和相对应的LOD值。        HeightMap类，读入高度图，并根据策略创建一定数量的CTerrainB...
Lod技术的实现
m0_51364829的博客
11-09 2047
Lod技术的实现 Lod，多层次细节，其实就是同一个物体，制作多个不同精度的模型，根据物体的重要性和位置来显示不同精度的模型。思想是空间换时间。 实现过程，在3Dmax 中制作三个不同分段数的球体。 导入Unity，给上不同颜色的材质，红色精度最高，蓝色次之，紫色最低 新建空物体，绑上Lod Group 组件 点击Add，把不同的模型挂载到不同百分比的Lod上，再把所有模型Reset复位 使用鼠标滚轮改变视口离物体的距离，近处物体是红色高精度模型，远处是紫色低精度模型 打开Stats面板可以看
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值