VB中利用OpenGL绘制三维地质构造图

6) 初始化颜色表

预先构造好颜色表。当读入高程数据时，可以根据Z值来访问颜色表，从而计算出该顶点的颜色（不同的高度层用不同的颜色显示）。

Private Sub InitColorTable()
'初始化颜色表
For i = 0 To 255
With ColorTable(i)
.R = 255 - i: .G = i: .B = 0
End With
Next i
For i = 256 To 510
With ColorTable(i)
.R = 0: .G = 255 - (i - 255): .B = i - 255
End With
Next i
End Sub

7) 从数据文件中获取数据

从data.txt文件中获取顶点数据。data.txt数据文件中第一行是构造图横向采样点的个数，第二行是构造图纵向采样点的个数，第三行为采样点的间距，第四行起为采样点的高度值。GetDataFromFile()过程从data.txt文件中读取数据，由于是等间距采样数据，所以顶点的X和Y坐标值可以直接计算出来，只有Z坐标值从数据文件中读入。获得顶点数据之后，计算出顶点的颜色，及填充三角形的法线矢量（用三角形序列填充）。

8) 设置窗口识见区和裁剪区

设置窗口识见区和修剪空间,保持坐标系到屏幕坐标的映射不变

glViewport 0, 0, w, h '定义视见区
glMatrixMode GL_PROJECTION '定义矩阵为投影矩阵
glLoadIdentity '用于在进行矩阵处理之前"复位"坐标系
If w <= h Then '定义修剪空间
glOrtho -2000, 2000, -1000* h / w, 2000* h / w, -10000, 10000
Else
glOrtho -1000 * w / h , 2000* w / h, -1000, 2000, -10000, 10000
End If
'告诉OpenGL将来的所有变换都会影响模型
glMatrixMode GL_MODELVIEW '定义矩阵为模型变换矩阵
glLoadIdentity

9) 绘制三维地质构造图

先进行旋转、比例和平移变换，再根据顶点数据进行绘图，运行结果如图2所示。绘制出来的三维地质构造图可以绕X、Y、Z轴进行360度自由旋转，也可以进行自动旋转；并可以对三维图分别沿X、Y、Z轴进行比例缩放，也可以使三维图沿X、Y、Z轴三个方向同时进行相同比例缩放。

　图2 运行结果

10) 把绘图描述表从线程上脱离并删除

Public Sub ReleaseGLRC(hGLRC As Long)
'把绘图描述表从线程上脱离并删除
If hGLRC <> 0 Then
wglMakeCurrent 0, 0
wglDeleteContext hGLRC
End If
End Sub

3. 结束语

目前，计算机图形学已进入三维时代，三维图形在人们四周无所不在，三维图形在军事、航空、航天、医学、地质勘探、艺术造型和文化娱乐等方面都有着十分广泛的应用，而OpenGL是一种三维图形开发标准，是从事三维图形开发工作的必要工具，其稳定性、可靠性、可扩展性等特点，赋予了其强大的生命力和应用前景，OpenGL已经广泛应用于在CAD/CAM/CAE、医学图像处理、虚拟现实、娱乐、广告等领域，因而对于OpenGL的研究和应用有着重大意义。