數(shù)據(jù)產(chǎn)品的質(zhì)量與NDVI和LST產(chǎn)品的質(zhì)量有關(guān),本研究中主要通過控制所使用的NVDI與LST數(shù)據(jù)產(chǎn)品的質(zhì)量來保證所生產(chǎn)數(shù)據(jù)產(chǎn)品的質(zhì)量。
對(duì)于LST產(chǎn)品,由于Aqua衛(wèi)星其過境掃描時(shí)間為下午,通常為云量最少的時(shí)間,最終選擇了基于Aqua衛(wèi)星的LST產(chǎn)品,考慮到在夜間地表溫度會(huì)顯著下降,繼而增加干旱檢測(cè)的不確定性,因此僅選擇使用MYD11A1.061產(chǎn)品中的日間LST數(shù)據(jù)。再依據(jù)LST產(chǎn)品中的QC_Day(日間質(zhì)量控制)波段,提取高質(zhì)量像元,該提取流程完全參照MODIS Collection(C61)Product User Guide進(jìn)行,流程如下:
(1)質(zhì)量控制波段QC_Day中Bits 0–1(Mandatory QA flags)為0的像元標(biāo)定為高質(zhì)量像元(不必進(jìn)行QA質(zhì)量確認(rèn));
(2)QC_Day中Bits 0–1為2和3的像元標(biāo)定為低質(zhì)量像元(由于云層或其他因素影響);
(3)QC_Day中Bits 0–1為1的像元標(biāo)定為待定像元(需要進(jìn)行進(jìn)一步的QA);
(4)對(duì)于待定像元,將其QC_Day中相應(yīng)像元Bits 2–7位的值進(jìn)行比較以確定該像元的質(zhì)量:當(dāng)Bits 2–3(Data quality flag)為0時(shí),為高質(zhì)量像元;當(dāng)Bits 2–3為1時(shí),查看其對(duì)應(yīng)的Bits 4–5(Emissivity error flag),和Bits 6–7(LST error flag)。對(duì)應(yīng)像元的Bits 4–5為0,Bits 6–7為0的像元為高質(zhì)量像元,其他像元為低質(zhì)量像元;
(5)將已標(biāo)定的像元進(jìn)行掩膜處理,僅提取高質(zhì)量像元,獲得新的LST數(shù)據(jù)。
考慮到云層對(duì)遙感地表溫度影響導(dǎo)致的數(shù)據(jù)缺失,以及去云算法在計(jì)算年時(shí)間尺度上造成的影響,本文采用ATC模型,確保LST數(shù)據(jù)集的準(zhǔn)確性和時(shí)空連續(xù)性。
對(duì)于NDVI產(chǎn)品,考慮到黃河流域位于北半球,從太陽(yáng)高度角考慮,Terra衛(wèi)星獲取的數(shù)據(jù)比Aqua衛(wèi)星的光譜數(shù)據(jù)效果好,同時(shí)也為保持與LST數(shù)據(jù)相同的空間分辨率,使用了基于Terra衛(wèi)星的MOD13A2.061 NDVI產(chǎn)品。按照MODIS的NDVI產(chǎn)品數(shù)據(jù)說明,依據(jù)產(chǎn)品中的SummaryQA波段對(duì)NVDI像元進(jìn)行質(zhì)量區(qū)分并提取高質(zhì)量像元。該提取流程完全參照MODIS Collection 6.1 (C61) Vegetation Index Product User Guide進(jìn)行。流程如下
[17]:
(1)將像元質(zhì)量保證波段SummaryQA中值為0的像元標(biāo)定為高質(zhì)量像元(Good Data);
(2)SummaryQA中值為1的像元標(biāo)定為邊際像元(Marginal Data);
(3)冰/雪覆蓋像元,在SummaryQA中的像元值為2,標(biāo)定為低質(zhì)量像元;
(4)云層遮蔽、不可見像元在SummaryQA中像元值為3,標(biāo)定為低質(zhì)量像元;
(5)掩膜處理,剔除邊際像元和低質(zhì)量像元,僅保留高質(zhì)量像元,獲得新的NDVI數(shù)據(jù)。