Python - GDAL 学习(一)

本科的时候跟着老师做过一些图像处理方面的工作,到别业设计的时候选择了图像水印技术在遥感图像方面的应用,之前打算用matlab来做,因为之前经常使用,而且参加了数学建模比赛,用的熟悉一些,但是读研究生的导师建议我用python来做,可以为以后的路打下基础,研究生的学习,我也想学习python来处理图像。

 
   
import gdal
fileName = 'f:/20010803.tif' #其实在这里存在一个问题,就是打开普通jpg图片的时候,用的是\而打开tif文件的时候,用的是/ 
dataset = gdal.Open(fileName)
if dataset == None: 
print(fileName+"文件无法打开") 
return im_width = dataset.RasterXSize #列数 
im_height = dataset.RasterYSize #行数
im_bands = dataset.RasterCount #波段数 
im_data = dataset.ReadAsArray(0,0,im_width,im_height)#获取数据 
im_geotrans = dataset.GetGeoTransform()#获取仿射矩阵信息 
im_proj = dataset.GetProjection()#获取投影信息 
band1= im_data[0,0:im_height,0:im_width] 
band2= im_data[1,0:im_height,0:im_width] 
band3= im_data[2,0:im_height,0:im_width] 
band4= im_data[3,0:im_height,0:im_width] 
band5= im_data[4,0:im_height,0:im_width] 
band6= im_data[5,0:im_height,0:im_width]#老师给我的tif文件是6个波段




但是遇到了第一个问题,将tif文件分波段后,不知道如何对tif文件进行显示和保存 




    


你可能感兴趣的:(Python学习)