python GDAL笔记

原创已于 2022-03-25 16:52:14 修改 · 1.1k 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#python

于 2022-03-25 16:51:12 首次发布

本文介绍如何使用GDAL库读取地理影像文件的基本信息，包括图像格式、分辨率、波段数量、投影信息及坐标转换模型等。通过Python代码演示了如何获取这些关键参数。

from osgeo import gdal
from osgeo.gdalconst import *
gdal.AllRegister()#一次性注册所有数据驱动，只能读不能写
ds=gdal.Open(fr"C:\Users\Fang\Desktop\merge1_rgbz_grids_Clip11.tif",GA_ReadOnly)
print(ds.GetDriver().ShortName)#输出图像的格式信息
rows=ds.RasterYSize#栅格数据的高度（Y方向上的像元个数）
cols=ds.RasterXSize
bands=ds.RasterCount#栅格数据波段数，在GDAL中每一个波段都是一个数据集，是属性不是方法
print(f"rows:{rows}")
print(f"cols:{cols}")
print(f"bands:{bands}")
print(ds.GetProjectionRef())#输出栅格数据的投影信息
adfGeoTransform = ds.GetGeoTransform()#输出图像的坐标和分辨率信息。含6个参数的坐标转换模型，坐标、分辨率和旋转等信息
print(adfGeoTransform[0])#左上角坐标
print(adfGeoTransform[3])
print(adfGeoTransform[1])#X方向像元大小，正值。分辨率
print(adfGeoTransform[5])#Y方向像元大小，为负值
band = ds.GetRasterBand(1)#读取第一行数据(以二进制的形式)，注意波段索引的开始值是1而不是0
print(band.DataTyped)#查看band这一波段的数据类型
print(band.GetRasterColorInterpretation())#查看波段的颜色类型
data = band.ReadRaster(0,0,ds.RasterXSize,1,ds.RasterXSize,1,band.DataType)#读取图像数据（以numpy数组的形式，类型为int16）