用VTK绘制DEM之一

DEM的显示是GIS的基本功能，利用GDAL和VTK强大的功能可以很容易做到这一点，而Python则使实现更为简单。

关于GDAL可以参考lilin的学习笔记http://wiki.woodpecker.org.cn/moin/lilin/gdal-index，VTK可以参考我博客里的其他文章。

先看看我的效果图，目前只显示地形，没有根据高程设置颜色，也没有叠加遥感影像。高程数据采用的是SRTM3，可以去http://www2.jpl.nasa.gov/srtm/cbanddataproducts.html下载。

先看代码再解释几个关键函数。

# !C:/Python24/python.exe
# coding=utf-8

# app.py
# author:liujunzhi
# date:20071210

import  vtk
from  vtk.wx.wxVTKRenderWindow  import  wxVTKRenderWindow

import  sys, time
from  gdal  import  gdal
from  wxPython.wx  import   *

VTK_POLYGON  =   7

def  display_dem(renWin, demName):
     # open dem file with gdal
    dataset  =  gdal.Open(demName)
    
     if  dataset  is  None:
         print   ' open DEM file failed! '
        sys.exit(0)
    
    geoTrans  =  dataset.GetGeoTransform()
    
    band  =  dataset.GetRasterBand( 1 )
     # ReadAsArray(self, xoff=0, yoff=0, win_xsize=None, win_ysize=None, bu