读取TIFF的行列数、投影转换信息、地图投影信息以及栅格像元值。
1、读取TIFF
用到的是GDAL模块的gdal.Open()函数。
dataset = gdal.Open(filepath)
2、获取TIFF的行列数
主要用到RasterXSize和RasterYSize来获取栅格数据的行列数。
col = dataset.RasterXSize
row = dataset.RasterYSize
3、获取投影转换信息
用到GetGeoTransform()函数。
geotrans = dataset.GetGeoTransform()
4、获取地图投影信息
用到GetProjection()函数。
proj = dataset.GetProjection()
5、获取栅格像元值
用到ReadAsArray()函数。
data = dataset.ReadAsArray()
全部代码:
def readtif(filepath):
dataset = gdal.Open(filepath)
col = dataset.RasterXSize
row = dataset.RasterYSize
geotrans = dataset.GetGeoTransform()
proj = dataset.GetProjection()
data = dataset.ReadAsArray()
del dataset
return [col, row, geotrans, proj, data]