Python 使用GDAL和rasterio读取栅格数据，并进行仿射变换信息的相互转换

使用gdal可以读取栅格数据，并获取投影及仿射变换信息

import sys
import gdal

dataset = gdal.Open('./zhuhai.tif')
if dataset is None:
    print('Unable to open *.tif')
    sys.exit(1)

projection = dataset.GetProjection()   # 获取投影信息
### GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],
    # AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],
    # UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],
    # AXIS["Latitude",NORTH],AXIS["Longitude",EAST],AUTHORITY["EPSG","4326"]]

transform = dataset.GetGeoTransform()    # 获取放射变换信息
### (113.39159774577, 5.000000000000014e-06, 0.0, 22.058475004479, 0.0, -5.000000000000014e-06)

img_array = dataset.ReadAsArray()   # 读取波段数据信息

使用rasterio读取栅格数据，并获取投影及仿射变换信息

import rasterio

dataset = rasterio.open('./zhuhai.tif')

projection = dataset.crs
# EPSG:4326

transf