Python gdal库读取tif文件

from osgeo import gdal    # GDAL库主要提供对栅格数据的处理，使用抽象数据模型来解析所支持的数据格式
import filename_cut as fc
import matplotlib.pyplot as plt

def read_tif(filename):
    dataset = gdal.Open(str(filename))
    if dataset == None:
        print(filename + "无法打开")
        return
    im_width = dataset.RasterXSize                              # 栅格矩阵的列数（宽）
    im_height = dataset.RasterYSize                             # 栅格矩阵的行数（高）
    im_bands = dataset.RasterCount                              # 波段数
    im_data = dataset.ReadAsArray(0, 0, im_width, im_height)    # 获取数据
    im_geotrans = dataset.GetGeoTransform()                     # 获取仿射矩阵信息
    im_proj = dataset.GetProjection()                           # 获取投影信息
    return im_data

data_path = "E:/FY_GPM_dataset/2020/FY/congxincut_384_280/0615_0015.tif"
plt.imshow(read_tif(data_path)[11, :, :])
plt.savefig('F:/hj/0615_0015_12.jpg')
plt.show()

可视化的数据结果