Matlab实现Landsat8热红外辐射定标

本篇博客主要记录一下用matlab实现landsat8热红外波段辐射定标的过程。因为参加比赛需要用到这个功能，但是折腾了好久才搞出来，所以想写出来给有需要的共享一下

首先是定标公式。
公式中：Lλ为输出结果，ML为增益参数即为MTL头文件中的(RADIANCE_MULT_BAND_x，x为需要进行定标的波段序号)，Qcal即为原始图像波段，AL为偏置参数即为MTL头文件中的（RADIANCE_ADD_BAND_x)

接下来进入代码部分。

%打开遥感影像
[Data,R] = geotiffread(inputfilepath);
info = geotiffinfo(inputfilepath);

figure(1);
imshow(Data);%输入图片

Data=single(Data)；

%辐射定标函数主体
RACAImage = Data*0.0003342+0.10000;

figure(2);
imagesc(RACAImage);%输出图片
    
 %结果保存
geotiffwrite(outputfilepath,RACAImage,R,'GeoKeyDirectoryTag', info.GeoTIFFTags.GeoKeyDirectoryTag);
%geotiffwrite('文件名',需要保存的矩阵,头文件信息, 'GeoKeyDirectoryTag', info.GeoTIFFTags.GeoKeyDirectoryTag);

代码当中需要注意的地方包括两个
1、matlab当中包含了处理tiff格式文件的函数，即为代码中读取和保存所调用的函数
2、需要将输入数据转为single类型才可以，直接使用做出来的结果只有整数。原始数据用matlab打开时uint16类型，但是看着用envi做的结果打开时single，所以我就试了一下定标前后分别将数据转为single类型，结果是需要定标前转为single类型才能得到正确结果，这部分原理我也不太清楚。
结果图如下，

图里面分别是用envi做出来的结果和用这段程序跑出来的结果，数值是一样的。