一霎微雨洒庭轩

2007-11-30

三维地形绘制－－四叉树递归算法 - [Algorithms]

版权声明：转载时请以超链接形式标明文章原始出处和作者信息及本声明
http://conoon.blogbus.com/logs/11308647.html

 #define EDGE_POINT 0
 #define NODE_POINT 1
 #define UNKNOWN 2

 class CDEM 
 {
 public:
 CString strFileName; //文件名
 CString strPathName; //文件路径名
 double **dDEMGrid; //原始DEM格网点高程数据
 int iCols; //研究区格网列数
 int iRows; //研究区格网行数
 double dXllcorner; //格网左下角X坐标（起始坐标）
 double dYllcorner; //格网左下角Y坐标（起始坐标）
 double dCellsize; //格网大小
 double dNodata; //无数据区域的标识数字
 double dMinZ,dMaxZ; //最大、最小高程值
 double dCenterX,dCenterY; //DEM中心的坐标
 float fScale; //范围比例
 float fZRatio; //DEM高程夸张系数
 bool IsShow; //是否显示
 CDEM();
 virtual ~CDEM();
 void GetDEMScope(); //获得DEM最大、最小高程值
 void Draw(); //绘制DEM
 bool GetPt(int i,int j,CPt3d &pt); //获得DEM中的i行j列点

 //四叉树算法实现
 int iLodRows,iLodCols; //经过LOD扩展后的行数和列数
 int** quadtree; //四叉树边点记录
 void init_quadtree(int i, int j); //初始化四叉树
 void setup_quadtree(int i, int j, int width); //建立四叉树
 void Draw(int i, int j, int width); //绘制四叉树
 void triangle(CPt3d &pt1, CPt3d &pt2, CPt3d &pt3); //绘制三角形
 };

 void CDEM::init_quadtree(int iRow, int iCol)
 {
 int i,j;
 quadtree = new int *[iRow];
 for(i=0; i<iRow; i++)
 {
 quadtree[i] = new int[iCol];
 for(j=0; j<iCol; j++)
 {
 quadtree[i][j] = UNKNOWN;
 }
 }
 }

 void CDEM::setup_quadtree(int i, int j, int width)
 {
 int width2;
 width2 = width / 2;
 if(width > 1)
 {
 quadtree[i][j] = NODE_POINT;
 setup_quadtree(i-width2,j-width2,width2);
 setup_quadtree(i-width2,j+width2,width2);
 setup_quadtree(i+width2,j-width2,width2);
 setup_quadtree(i+width2,j+width2,width2);
 }
 else
 quadtree[i][j] = EDGE_POINT;
 }

 void CDEM::Draw(int i, int j, int width)
 {
 int width2;
 width2 = width / 2;
 if((width > 1) && (quadtree[i][j] == NODE_POINT))
 {
 Draw(i-width2,j-width2,width2);
 Draw(i-width2,j+width2,width2);
 Draw(i+width2,j-width2,width2);
 Draw(i+width2,j+width2,width2);
 }
 else if(quadtree[i][j] == EDGE_POINT)
 {
 //绘制8个三角行
 CPt3d pt[9];
 GetPt(i, j, pt[0]);

 GetPt(i-width, j-width, pt[1]);
 GetPt(i, j-width, pt[2]);
 GetPt(i+width, j-width, pt[3]);
 GetPt(i+width, j, pt[4]);

 GetPt(i+width, j+width, pt[5]);
 GetPt(i, j+width, pt[6]);
 GetPt(i-width, j+width, pt[7]);
 GetPt(i-width, j, pt[8]);

 int m;
 for(m=1; m<8; m++)
 triangle(pt[0], pt[m], pt[m+1]);
 triangle(pt[0], pt[8], pt[1]);

 }
 else 
 {
 return;
 }
 }

 void CDEM::triangle(CPt3d &pt1, CPt3d &pt2, CPt3d &pt3)
 { 
 if (pt1.z != -9999.0f * fZRatio && pt2.z != -9999.0f*fZRatio &&
 pt3.z != -9999.0f * fZRatio)
 { 
 glBegin(GL_LINE_LOOP);
 glBegin(GL_TRIANGLES); 

 glColor4f(1.0f, 0.0f, 0.0f, 1.0f);
 glVertex3f(pt1.x, pt1.y, pt1.z);

 glColor4f(1.0f, 0.0f, 0.0f, 1.0f);
 glVertex3f(pt2.x, pt2.y, pt2.z);

 glColor4f(1.0f, 0.0f, 0.0f, 1.0f);
 glVertex3f(pt3.x, pt3.y, pt3.z);

 glEnd();
 glEnd();
 }
 }