使用Rasterio读取栅格数据的实例讲解

affine.Affine(a, b, c,
    d, e, f)

import rasterio

with rasterio.open('example.tif') as ds:
 print('该栅格数据的基本数据集信息（这些信息都是以数据集属性的形式表示的）：')
 print(f'数据格式：{ds.driver}')
 print(f'波段数目：{ds.count}')
 print(f'影像宽度：{ds.width}')
 print(f'影像高度：{ds.height}')
 print(f'地理范围：{ds.bounds}')
 print(f'反射变换参数（六参数模型）：\n {ds.transform}')
 print(f'投影定义：{ds.crs}')
 # 获取第一个波段数据，跟GDAL一样索引从1开始
 # 直接获得numpy.ndarray类型的二维数组表示，如果read()函数不加参数，则得到所有波段（第一个维度是波段）
 band1 = ds.read(1)
 print(f'第一波段的最大值：{band1.max()}')
 print(f'第一波段的最小值：{band1.min()}')
 print(f'第一波段的平均值：{band1.mean()}')
 # 根据地理坐标得到行列号
 x, y = (ds.bounds.left + 300, ds.bounds.top - 300) # 距离左上角东300米，南300米的投影坐标
 row, col = ds.index(x, y) # 对应的行列号
 print(f'(投影坐标{x}, {y})对应的行列号是({row}, {col})')
 # 根据行列号得到地理坐标
 x, y = ds.xy(row, col) # 中心点的坐标
 print(f'行列号({row}, {col})对应的中心投影坐标是({x}, {y})')
 # 那么如何得到对应点左上角的信息
 x, y = (row, col) * ds.transform
 print(f'行列号({row}, {col})对应的左上角投影坐标是({x}, {y})')

该栅格数据的基本数据集信息（这些信息都是以数据集属性的形式表示的）：
数据格式：GTiff
波段数目：3
影像宽度：4800
影像高度：4800
地理范围：BoundingBox(left=725385.0, bottom=2648415.0, right=869385.0, top=2792415.0)
反射变换参数（六参数模型）：
 | 30.00, 0.00, 725385.00|
| 0.00,-30.00, 2792415.00|
| 0.00, 0.00, 1.00|
投影定义：CRS({'init': 'epsg:32649'})
第一波段的最大值：5459
第一波段的最小值：-313
第一波段的平均值：489.80300625
(投影坐标725685.0, 2792115.0)对应的行列号是(10, 10)
行列号(10, 10)对应的中心投影坐标是(725700.0, 2792100.0)
行列号(10, 10)对应的左上角投影坐标是(725685.0, 2792115.0)