MODIS地表温度数据11A2的质量控制

地表温度是地球系统的重要参数，最近对地表温度数据进行了详细的了解，这里记录一下。

下载就不说了：Level-1 and Atmosphere Archive & Distribution System Distributed Active Archive Center - LAADS DAAC

MODIS官方的地表温度数据如下：

最近主要关注MODIS11A2（https://lpdaac.usgs.gov/products/myd11a2v006/）。11A2产品分Terra和Aqua两个传感器，Terra的过境时间是上午10:30和晚上22:30，Aqua的过境时间是下午1:30和凌晨1:30，因此Aqua传感器的MYD11A2数据相对来说更能代表白天和夜晚的温度，这个产品应该是当前文章中用的最多的了，各种质量的文章都有。

但是使用之前一定要搞清楚质量控制，官方文档可以直接下载，这里放一个网页版的MODIS LST Products Users' Guide。

官方的质量控制表如下：

这个表格其实代表的是一个8位的2进制数，在进行数据筛选时，需要首先把10和11值去除，也就是2、3，分别表示由于云或其他因素的影响没有反演的地区（例如海洋的地区就是3），TBD的意思应该是To Be Determined。这个质量控制不仅给了温度的误差，还给了发射率的误差，发射率就是地物的实际热辐射与同温同波长下黑体热辐射的比值，这个值对于我们反演地表温度非常重要，但是我觉得对于用户来说关注点就是地表温度的话不用太考虑二进制的第5&4bit的值

如果想要获得LST error<=1K的数据，那么QC的最大值为00111111（十进制63），因此可以直接选择QC<63的区域，

如果想要获得LST error<=2K的数据，那么QC的最大值为01111111（十进制127），因此可以直接选择QC<127的区域。

关于地表温度反演，MYD2A2是基于daily的MYD12A1产品求8天均值得到的，MYD12A1的daily数据是使用最常见的分裂窗算法得到的，关于地表温度反演，可看下面这篇遥感学报综述文章，讲的很仔细，也可以直接看MODIS算法的创立者Zhengming Wan在RSE的系列文章：

李召良,段四波,唐伯惠,吴骅,任华忠,阎广建,唐荣林,冷佩.热红外地表温度遥感反演方法研究进展[J].遥感学报,2016,20(05):899-920.

Zhengming Wan,New refinements and validation of the collection-6 MODIS land-surface temperature/emissivity product, Remote Sensing of Environment, 2014, 140(1):36-45..https://doi.org/10.1016/j.rse.2013.08.027

（最近又看了一下et的数据，如果大家要用et的话QC的最大值为00100110=38）