“””
暗像元校正
方法:计算波段中的最小值进行暗像元校正
过程: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