当前位置： 移动技术网 > IT编程>脚本编程>Python > 利用python GDAL库读写geotiff格式的遥感影像方法

利用python GDAL库读写geotiff格式的遥感影像方法

2018年12月04日  | 移动技术网IT编程  | 我要评论

半边生鳞不生角,63岁女教师被误认模特,星空战记sf

如下所示：

from osgeo import gdal
import numpy as np
def read_tiff(inpath):
  ds=gdal.open(inpath)
  row=ds.rasterxsize
  col=ds.rasterysize
  band=ds.rastercount
  geotransform=ds.gettransform()
  proj=ds.gettransform()
  data=np.zeros([row,col,band])
  for i in range(band):
   dt=ds.getrasterband(1)
   data[:,:,i]=dt.readasarray(0,0,col,row)
  return data
 
def array2raster(outpath,array,geotransform,proj):
 cols=array.shape[1]
 rows=array.shape[0]
 driver=gdal.getdriverbyname('gtiff')
 outraster=driver.create(newrasterfn,cols,rows,1,gdal.gdt_byte)
 outraster.setgeotransform(geotransform)#参数2,6为水平垂直分辨率，参数3,5表示图片是指北的
 outband=outraster.getrasterband(1)
 outband.writearray(array)
 outraster.setprojection(proj)#将几何对象的数据导出为wkt格式
 outraster.flushcache()
 
if _name=="_main_":
 
 data,geotransform,proj=read_tiff('d:/a.tif')
 
 array2raster("d:/b.tif",np.zeros[2400,2400],geotransform,proj)

利用python gdal库读写geotiff格式的遥感影像，并生成与原影像具有相同地理坐标和投影坐标的geotiff格式图片。

以上这篇利用python gdal库读写geotiff格式的遥感影像方法就是小编分享给大家的全部内容了，希望能给大家一个参考，也希望大家多多支持移动技术网。

您可能感兴趣的文章:

如对本文有疑问，请在下面进行留言讨论，广大热心网友会与你互动！！ 点击进行留言回复

相关文章:

验证码：
最近更新的文章
大家感兴趣的文章
移动技术网