Google earth engine对Landsat8辐射定标

• 辐射定标的本质是波段运算，y = ax+b
• 以第10波段为例，首先从元数据文件找到RADIANCE_MULT_BAND_10和RADIANCE_ADD_BAND_10，第一个是乘项 3.3420E-04，第二个是加项0.1。在gee中进行波段运算

var rad10 = image.select('B10').multiply(3.3420*0.0001).add(0.1)

• 以下是一个Landsat8影像自带的元数据文件。

GROUP = L1_METADATA_FILE
  GROUP = METADATA_FILE_INFO
    ORIGIN = "Image courtesy of the U.S. Geological Survey"
    REQUEST_ID = "0702003135965_00029"
    LANDSAT_SCENE_ID = "LC81220352020065LGN00"
    LANDSAT_PRODUCT_ID = "LC08_L1TP_122035_20200305_20200314_01_T1"
    COLLECTION_NUMBER = 01
    FILE_DATE = 2020-03-14T05:18:46Z
    STATION_ID = "LGN"
    PROCESSING_SOFTWARE_VERSION = "LPGS_13.1.0"
  END_GROUP = METADATA_FILE_INFO
  GROUP = PRODUCT_METADATA
    DATA_TYPE = "L1TP"
    COLLECTION_CATEGORY = "T1"
    ELEVATION_SOURCE = "GLS2000"
    OUTPUT_FORMAT = "GEOTIFF"
    SPACECRAFT_ID = "LANDSAT_8"
    SENSOR_ID = "OLI_TIRS"
    WRS_PATH = 122
    WRS_ROW = 35
    NADIR_OFFNADIR = "NADIR"
    TARGET_WRS_PATH = 122
    TARGET_WRS_ROW = 35
    DATE_ACQUIRED = 2020-03-05
    SCENE_CENTER_TIME = "02:48:28.1796780Z"
    CORNER_UL_LAT_PRODUCT = 37.09079
    CORNER_UL_LON_PRODUCT = 115.64750
    CORNER_UR_LAT_PRODUCT = 37.09195
    CORNER_UR_LON_PRODUCT = 118.24675
    CORNER_LL_LAT_PRODUCT = 34.97083
    CORNER_LL_LON_PRODUCT = 115.68327
    CORNER_LR_LAT_PRODUCT = 34.97190
    CORNER_LR_LON_PRODUCT = 118.21378
    CORNER_UL_PROJECTION_X_PRODUCT = 379800.000
    CORNER_UL_PROJECTION_Y_PRODUCT = 4105800.000
    CORNER_UR_PROJECTION_X_PRODUCT = 610800.000
    CORNER_UR_PROJECTION_Y_PRODUCT = 4105800.000
    CORNER_LL_PROJECTION_X_PRODUCT &#