#?-*-?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數(shù)據(jù)

dataset?=?gdal.Open（filePath）
adfGeoTransform?=?dataset.GetGeoTransform（）
band?=?dataset.GetRasterBand（1）???#用gdal讀寫數(shù)據(jù)