这篇文章给大家分享的是有关如何利用python GDAL库读写geotiff格式的遥感影像的内容。小编觉得挺实用的，因此分享给大家做个参考，一起跟随小编过来看看吧。
成都创新互联公司专注于夏津网站建设服务及定制，我们拥有丰富的企业做网站经验。 热诚为您提供夏津营销型网站建设，夏津网站制作、夏津网页设计、夏津网站官网定制、小程序定制开发服务，打造夏津网络公司原创品牌,更为您提供
夏津网站排名全网营销落地服务。python可以做什么
Python是一种编程语言，内置了许多有效的工具，Python几乎无所不能，该语言通俗易懂、容易入门、功能强大，在许多领域中都有广泛的应用，例如最热门的大数据分析，人工智能，Web开发等。
如下所示：
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格式的遥感影像”这篇文章就分享到这里了，希望以上内容可以对大家有一定的帮助，让大家可以学到更多知识，如果觉得文章不错，可以把它分享出去让更多的人看到吧！
当前名称：如何利用pythonGDAL库读写geotiff格式的遥感影像-创新互联
当前网址：
http://cdkjz.cn/article/dsdcod.html
