利用python进行Landsat8数据的辐射定标

上次用于处理Landsat8数据的类，增加了辐射定标的功能。目前是用于可见光波段的辐射定标。

这段辐射定标的代码，首先从元文件中读取定标参数，然后进行辐射定标。对NaN数据做了处理，使其一直保持为NaN。最后可以利用Landsat8类中的write函数，将辐射定标后的数据写出。

当然，这段代码的有效性，仅仅保证在我的数据上好用。

具体辐射定标的实现，参看Landsat8Reader类的radiometric_calibration函数。

import os
from osgeo import gdal
from osgeo import gdal_array
import numpy as np
from show import TwoPercentLinear
from matplotlib import pyplot as plt
import cv2 as cv

class Landsat8Reader(object):
    def __init__(self):
        self.base_path =\
        "LC81220352018123LGN00/LC08_L1TP_122035_20180503_20180516_01_T1"
        self.bands = 7
        self.band_file_name = []
        self.nan_position = []
    
    def read(self):
        for band in range(self.bands):
            band_name = self.base_path + "_B" + str(band+1) + ".tif" 
            self.band_file_name.append(band_name)
               
        ds = gdal.Open(self.band_file_name[0])
        image_dt = ds.GetRasterBand(1).DataType
        image = np.zeros((ds.RasterYSize, ds.RasterXSize, self.bands),
                        dtype = np.float)
        for band in range(self.bands):
            ds = gdal.Open(self.band_file_name[band