本篇博客主要记录一下用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做出来的结果和用这段程序跑出来的结果,数值是一样的。