读取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]

更多推荐

python 读取TIFF信息