使用Python使用大气校正法计算地表温度

使用Python使用大气校正法计算地表温度

前言

也有段时间没有跟新博客了，这次博客就是用新学的python语言来进行一个地表温度的计算，也算是承接了之前的内容吧！

---

一、具体原理及方法

这里不再赘述，有关原理及方法可以查阅我之前的博客文章https://blog.youkuaiyun.com/qq_44935981/article/details/90319487?spm=1001.2014.3001.5501

二、代码及步骤如下

1.引入库

代码如下：

import os
import numpy as np
from osgeo import gdal
import glob
import math
import matplotlib.pyplot as plt
#author sansuixueshi

2.读入数据及处理

代码如下：

input_path = 'E:/study/2018.4.09/2018.4.9/LC81300372018099LGN00/'
output_path = 'E:/study/2018.4.09/2018.4.9/LC81300372018099LGN00/out_data/'
file_all = os.listdir(input_path)
file_tail='.TIF'
radiance_scale=[1.2519E-02,1.2819E-02,1.1813E-02,9.9612E-03,6.0957E-03,1.5160E-03,5.1096E-04,1.1273E-02,2.3824E-03,3.3420E-04,3.3420E-04]
radiance_add=[-62.59278,-64.09577,-59.06370,-49.80584,-30.47869,-7.57977,-2.55479,-56.36651,-11.91176,0.10000,0.10000]
i=0

if not os.path.exists(output_path):
    os.mkdir(output_path)

for file_i in file_all:
    if file_i.endswith(file_tail):
        file_name = input_path + file_i
        (prename,suffix) = os.path.split(file_i)
        if i == 5:
            band = gdal.Open(file_name)
            red = band.ReadAsArray()*radiance_scale[3]+radiance_add[3]
            red = red.astype(np.float32)
        
        if i==6:
            band = gdal.Open(file_name)
            nir = band.ReadAsArray()*radiance_scale[4]+radiance_add[4]
            nir = nir.astype(np.float32)
        
        if i==1:
            band = gdal.Open(file_name)
            thermal = band.ReadAsArray