DEM数据三维可视化--python实现

# -*- coding: utf-8 -*-
“““
Created on Mon Jan  7 12:53:05 2019

@author: Administrator
“““
from mpl_toolkits.mplot3d import Axes3D
from matplotlib import cbook
from matplotlib import cm
from matplotlib.colors import LightSource
import matplotlib.pyplot as plt
import numpy as np
import tkinter
from tkinter import filedialog
from osgeo import gdal
gdal.AllRegister（）



filePath=““
def openfiles（）:
    fname = filedialog.askopenfilename（title=‘打开DEM文件‘ filetypes=[（‘DEM‘ ‘*.tif‘） （‘All Files‘ ‘*‘）]）
    global filePath
    filePath=fname
 
root = tkinter.Tk（）
#root.geometry（‘500x300+500+200‘）
btn1 = tkinter.Button（root text=‘打开DEM文件‘font =（“宋体“12‘bold‘）width=12height=2 command=openfiles）

btn1.pack（side=‘top‘）

root.mainloop（）
 #输入你的dem数据

dataset = gdal.Open（filePath）
adfGeoTransform = dataset.GetGeoTransform（）
band = dataset.GetRasterBand（1）   #用gdal读写数据