gdal+python 黑暗像元校正

“””
暗像元校正
方法：计算波段中的最小值进行暗像元校正
过程：1.影像DN值转为辐射能量值（L = DN*gain +bias）
2.计算相对反射率（ρ=π×D2×L/（ESUNI×COS(SZ)））
3.选择波段最小值进行波段值运算（dataset.data =bands-bandsMin）
4.新建栅格数据写入数据（driver = gdal.Create(self, *args, **kwargs)）

@time：2018/4/2 9:40
@author：wujd
“””
from osgeo import gdal
import numpy as np
import os

打开数据文件

os.chdir(r”F:\Wujd\0330darkPixel\3285750”)
dataset = gdal.Open(“HJ1A-CCD1-451-63-20171202-L20003285750-1.TIF”)

获取栅格数据的行列数，波段数，放射矩阵，投影信息

im_width = dataset.RasterXSize
im_height = dataset.RasterYSize
im_bands = dataset.RasterCount
im_geotrans = dataset.GetGeoTransform()
im_proj = dataset.GetProjection()
im_data = da