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

最新推荐文章于 2025-06-23 13:18:22 发布

原创最新推荐文章于 2025-06-23 13:18:22 发布 · 1.1w 阅读

37 ·

CC 4.0 BY-SA版权

文章标签：

#python #遥感 #tiff #gdal

python 同时被 2 个专栏收录

2 篇文章

订阅专栏

遥感

1 篇文章

订阅专栏

本文介绍如何使用 Python 的 GDAL 库来读取和写入 GeoTIFF 格式的遥感影像文件，并确保输出文件具有与原始文件相同的地理坐标和投影坐标。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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格式图片