Python GDAL库读写GeoTIFF遥感影像完整指南

小林子

发布时间：2026-10-13| 6人浏览过

原创

本文介绍如何利用Python GDAL库读取GeoTIFF格式遥感影像，并生成具有相同地理坐标和投影信息的新影像。通过示例代码展示如何获取影像行列数、波段数、地理变换参数及投影信息，并将数组数据写回新的GeoTIFF文件，确保输出影像与原图地理参考一致。

具体实现代码如下所示：

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格式遥感影像的方法就是小编分享给大家的全部内容了，希望能给大家一个参考，也希望大家多多支持本站。