由于下一步准备做辐射定标,因此需要从元数据文件中获取。
代码如下:
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]