利用python读取Landsat8的定标参数

最新推荐文章于 2023-06-27 16:42:22 发布

Remote Sensing

最新推荐文章于 2023-06-27 16:42:22 发布

阅读量1.1k

点赞数

分类专栏：遥感图像处理

本文链接：https://blog.csdn.net/RSstudent/article/details/107510515

版权

遥感图像处理专栏收录该内容

16 篇文章 31 订阅

订阅专栏

由于下一步准备做辐射定标，因此需要从元数据文件中获取。

代码如下：

filename =\
"LC81220352018123LGN00/LC08_L1TP_122035_20180503_20180516_01_T1_MTL.txt"
f = open(filename, 'r')
metadata = f.readlines()
f.close()

multi_parameters = []
add_parameters = []

parameter_start_line = 0
#检测定标参数的位置
for lines in metadata:
    test_line = lines.split('=')
    if test_line[0] == '    RADIANCE_MULT_BAND_1 ':
        break
    else:
        parameter_start_line = parameter_start_line + 1

for lines in range(parameter_start_line, parameter_start_line+11):
    parameter = float(metadata[lines].split('=')[1])
    multi_parameters.append(parameter)

for lines in range(parameter_start_line+11, parameter_start_line+22):
    parameter = float(metadata[lines].split('=')[1])
    add_parameters.append(parameter)
    
print(add_parameters)
print(multi_parameters)

结果如下：

[-61.79147, -63.27522, -58.30757, -49.16823, -30.0885, -7.48274, -2.52208, -55.6449, -11.75927, 0.1, 0.1]
[0.012358, 0.012655, 0.011662, 0.0098336, 0.0060177, 0.0014965, 0.00050442, 0.011129, 0.0023519, 0.0003342, 0.0003342]