GB/T 41534-2022 地表溫度遙感產(chǎn)品真實(shí)性檢驗(yàn)（正式版）

ICS35.240.70GB/T41534—2022國(guó)家市場(chǎng)監(jiān)督管理總局國(guó)家標(biāo)準(zhǔn)化管理委員會(huì)GB/T41534—2022 I Ⅱ 3術(shù)語和定義 4基本要求 2 25.1檢驗(yàn)方法概述 25.2檢驗(yàn)方法選擇 2 3 56檢驗(yàn)報(bào)告 96.1封面信息 96.2正文信息 6.3檢驗(yàn)報(bào)告信息簡(jiǎn)表 附錄A(規(guī)范性)點(diǎn)源尺度地表溫度相對(duì)真值獲取方法 附錄B(資料性)非均質(zhì)地表像元尺度地表溫度相對(duì)真值獲取方法 附錄C(資料性)輔助數(shù)據(jù)的獲取方法 附錄D(資料性)地表溫度遙感產(chǎn)品真實(shí)性檢驗(yàn)不確定度誤差來源分析 附錄E(資料性)地表溫度遙感產(chǎn)品真實(shí)性檢驗(yàn)報(bào)告信息簡(jiǎn)表樣例 IGB/T41534—2022本文件按照GB/T1.1—2020《標(biāo)準(zhǔn)化工作導(dǎo)則第1部分：標(biāo)準(zhǔn)化文件的結(jié)構(gòu)和起草規(guī)則》的規(guī)定起草。本文件由中國(guó)科學(xué)院提出。本文件由全國(guó)遙感技術(shù)標(biāo)準(zhǔn)化技術(shù)委員會(huì)(SAC/TC327)歸口。本文件起草單位：中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所、昆明理工大學(xué)、中國(guó)科學(xué)院地理科學(xué)與資源研究所、中國(guó)科學(xué)院空天信息創(chuàng)新研究院、北京師范大學(xué)、河北地質(zhì)大學(xué)、國(guó)家衛(wèi)星氣象中心、GB/T41534—2022地表溫度是地球表層系統(tǒng)過程中一個(gè)非常重要的特征物理量，它體現(xiàn)了陸地、海洋和大氣之間相互作用的綜合結(jié)果，是研究地表和大氣之間物質(zhì)和能量交換、全球海洋環(huán)流、氣候變化異常等科學(xué)問題不可或缺的一個(gè)重要參數(shù)，涉及眾多基礎(chǔ)學(xué)科和重大應(yīng)用領(lǐng)域。利用熱紅外遙感數(shù)據(jù)可反演出空間和時(shí)本文件針對(duì)地表溫度遙感產(chǎn)品的真實(shí)性檢驗(yàn)問題，結(jié)合國(guó)際上現(xiàn)有的各種地表溫度產(chǎn)品真實(shí)性檢品真實(shí)性檢驗(yàn)方法和評(píng)價(jià)標(biāo)準(zhǔn)，使其在保證先進(jìn)性的同時(shí)具有一定的可操作性和普適性，為地表溫度遙感產(chǎn)品的精度評(píng)價(jià)和不確定度分析提供統(tǒng)一的標(biāo)準(zhǔn)，作為與地表溫度遙感產(chǎn)品真實(shí)性檢驗(yàn)相關(guān)的各行業(yè)共同標(biāo)準(zhǔn)，為推動(dòng)我國(guó)對(duì)地觀測(cè)技術(shù)的創(chuàng)新與發(fā)展提供技術(shù)支撐。本文件的發(fā)布機(jī)構(gòu)提請(qǐng)注意，聲明符合本文件時(shí)，可能涉及到附錄A中A.3與“一種非接觸式地表溫度測(cè)定儀及方法(ZL201711179814.6)”相關(guān)的專利的使用。本文件的發(fā)布機(jī)構(gòu)對(duì)于該專利的真實(shí)性、有效性和范圍無任何立場(chǎng)。該專利持有人已向本文件的發(fā)布機(jī)構(gòu)承諾，他愿意同任何申請(qǐng)人在合理且無歧視的條款和條件下，就專利授權(quán)許可進(jìn)行談判。該專利持有人的聲明已在本文件的發(fā)布機(jī)構(gòu)備案。相關(guān)信息可以通過以下聯(lián)系方式獲得：專利持有人：中國(guó)科學(xué)院地理科學(xué)與資源研究所，中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所。地址：北京市朝陽區(qū)大屯路甲11號(hào)。請(qǐng)注意除上述專利外，本文件的某些內(nèi)容仍可能涉及專利。本文件的發(fā)布機(jī)構(gòu)不承擔(dān)識(shí)別專利的責(zé)任。1GB/T41534—2022地表溫度遙感產(chǎn)品真實(shí)性檢驗(yàn)本文件適用于航空和航天遙感數(shù)據(jù)地表溫度反演產(chǎn)品的真實(shí)性檢驗(yàn)。2規(guī)范性引用文件下列文件中的內(nèi)容通過文中的規(guī)范性引用而構(gòu)成本文件必不可少的條款。其中，注日期的引用文本文件。GB/T36296—2018遙感產(chǎn)品真實(shí)性檢驗(yàn)導(dǎo)則GB/T39468—2020陸地定量遙感產(chǎn)品真實(shí)性檢驗(yàn)通用方法3術(shù)語和定義3.5GB/T36296—2018、GB/T39468—2020界定的以及下列術(shù)語和定義適用于本文件。亮溫表征地[球]表[面]厚度等于穿透深度(范圍0.1～10倍波長(zhǎng))的表皮的綜合溫度。真實(shí)性檢驗(yàn)validation通過獨(dú)立方法評(píng)價(jià)待檢遙感產(chǎn)品與客觀實(shí)際的吻合程度并分析其不確定度的過程。測(cè)試結(jié)果或測(cè)量結(jié)果與真值間的一致程度。2GB/T41534—20223.64基本要求地表溫度遙感產(chǎn)品真實(shí)性檢驗(yàn)應(yīng)符合GB/T36296—2018中第7章的規(guī)定，并符合以下要求：a)參考對(duì)象的選取宜考慮不同下墊面、不同氣候條件，且應(yīng)有較長(zhǎng)的時(shí)間序列(大于或等于一年);b)參考對(duì)象地表溫度值應(yīng)具有較大的動(dòng)態(tài)取值范圍，包含高溫(310K)、中溫(290K)和低溫c)當(dāng)待檢驗(yàn)對(duì)象和參考對(duì)象之間存在空間尺度差異時(shí)，應(yīng)獲取與待檢驗(yàn)對(duì)象空間尺度一致的地表溫度相對(duì)真值后再進(jìn)行真實(shí)性檢驗(yàn)；e)當(dāng)待檢驗(yàn)對(duì)象中異常值較大(方差大于3σ)或較多(大于樣本數(shù)的1/3)時(shí)，可用中值代替均值來表示待檢地表溫度遙感產(chǎn)品和參考地表溫度遙感產(chǎn)品間的準(zhǔn)確度和不確定度；定各檢驗(yàn)方法的可靠性。5檢驗(yàn)方法a)直接檢驗(yàn)法是直接利用載荷過境時(shí)刻的地面同步或準(zhǔn)同步實(shí)測(cè)數(shù)據(jù)對(duì)地表溫度遙感產(chǎn)品進(jìn)行b)交叉檢驗(yàn)法是通過與精度已知的參考地表溫度遙感產(chǎn)品的比對(duì)分析，間接對(duì)待檢驗(yàn)的地表溫度遙感產(chǎn)品進(jìn)行檢驗(yàn)與評(píng)估的方法，并給出待檢驗(yàn)地表溫度遙感產(chǎn)品和參考地表溫度產(chǎn)品間的準(zhǔn)確度和不確定度的統(tǒng)計(jì)指標(biāo)；c)正向模型仿真檢驗(yàn)法是通過與入瞳處載荷仿真模擬值的比對(duì)分析，間接對(duì)待檢驗(yàn)的地表溫度遙感產(chǎn)品進(jìn)行檢驗(yàn)與評(píng)估的方法，并給出待檢驗(yàn)地表溫度遙感產(chǎn)品的準(zhǔn)確度和不確定度的統(tǒng)計(jì)指標(biāo)；d)時(shí)空變化趨勢(shì)分析檢驗(yàn)法是通過對(duì)待檢地表溫度遙感產(chǎn)品的時(shí)空變化特征和趨勢(shì)進(jìn)行分析，間接對(duì)地表溫度遙感產(chǎn)品進(jìn)行檢驗(yàn)與評(píng)估的方法，并給出待檢驗(yàn)地表溫度遙感產(chǎn)品的準(zhǔn)確度和不確定度的統(tǒng)計(jì)指標(biāo)。3GB/T41534—2022b)當(dāng)缺乏或者難以獲取像元尺度地表溫度觀測(cè)數(shù)據(jù)，但存在或者可獲取質(zhì)量可靠，且具有代表性的像元尺度地表發(fā)射率數(shù)據(jù)和大氣廓線數(shù)據(jù)時(shí)，采用正向模型仿真檢驗(yàn)法；c)當(dāng)缺乏地面觀測(cè)數(shù)據(jù)，但存在其他質(zhì)量可靠的地表溫度遙感產(chǎn)品時(shí)，選擇交叉檢驗(yàn)法；d)當(dāng)同時(shí)缺乏地面觀測(cè)數(shù)據(jù)和其他質(zhì)量可靠的地表溫度遙感產(chǎn)品時(shí)，選擇時(shí)空變化趨勢(shì)分析檢驗(yàn)法。5.3直接檢驗(yàn)法5.3.1實(shí)施條件和要求直接檢驗(yàn)法的實(shí)施條件和要求如下：a)應(yīng)具備地表溫度數(shù)據(jù)采集條件，且載荷過境時(shí)刻晴空無云，風(fēng)力應(yīng)小于4級(jí)；b)驗(yàn)證場(chǎng)地應(yīng)具有大于(3×3)個(gè)載荷像元對(duì)應(yīng)的地面面積，驗(yàn)證場(chǎng)地的地表類型均質(zhì)或者地表類型簡(jiǎn)單且分布已知的非均質(zhì)地表；c)驗(yàn)證場(chǎng)應(yīng)相對(duì)地勢(shì)平坦，對(duì)于具有一定三維結(jié)構(gòu)的地類，如稀疏植被和壟狀分布的農(nóng)田等，宜考慮載荷過境時(shí)刻太陽導(dǎo)致的陰影效應(yīng)；d)樣本應(yīng)選擇溫度空間變異性(標(biāo)準(zhǔn)差)不高于0.5K的均一地表，例如湖面、海面、茂密植被e)樣點(diǎn)布設(shè)的地表溫度采集儀器觀測(cè)精度高且穩(wěn)定，周邊無明顯影響數(shù)據(jù)采集的大目標(biāo)遮擋物；f)地表溫度測(cè)量宜采用全自動(dòng)方式，并應(yīng)配套采集其他輔助數(shù)據(jù)，包括大氣下行輻射數(shù)據(jù)和地表發(fā)射率；其中，大氣下行輻射數(shù)據(jù)應(yīng)同步采集，地表發(fā)射率的測(cè)量時(shí)刻宜接近真實(shí)性檢驗(yàn)時(shí)刻，且宜考慮地表發(fā)射率的時(shí)效性，樣本不應(yīng)發(fā)生降雨或者林火等改變發(fā)射率數(shù)值的事件；g)地表溫度的時(shí)間采樣間隔應(yīng)小于1min;特殊情況下，地表溫度觀測(cè)時(shí)間與載荷數(shù)據(jù)獲取時(shí)間的最大差值不應(yīng)超過30min;h)地表溫度采集儀器應(yīng)進(jìn)行定標(biāo)，且能溯源到國(guó)際單位制參考標(biāo)準(zhǔn)；定標(biāo)精度的分析結(jié)果應(yīng)包括在真實(shí)性檢驗(yàn)報(bào)告中；i)對(duì)于短期(小于6個(gè)月)的真實(shí)性檢驗(yàn)，地表溫度采集儀器在試驗(yàn)前后都應(yīng)進(jìn)行定標(biāo)；對(duì)于長(zhǎng)期的真實(shí)性檢驗(yàn)，紅外輻射計(jì)應(yīng)周期性(1次/6月)的進(jìn)行輻射定標(biāo)；j)宜同步觀測(cè)氣溫、氣壓、風(fēng)和濕度等常規(guī)氣象要素，以配合選擇真實(shí)性檢驗(yàn)的有效數(shù)據(jù)。5.3.2空間異質(zhì)性分析和空間抽樣方法空間異質(zhì)性分析和空間抽樣方法如下：a)空間異質(zhì)性的分析宜考慮高空間分辨率的地表溫度遙感產(chǎn)品或者其他輔助參量的空間分布特征，例如植被覆蓋度或者歸一化植被指數(shù)等，按照GB/T39468—2020中4.2的規(guī)定來確定樣本的異質(zhì)性；b)對(duì)于均質(zhì)地表，應(yīng)按照GB/T39468—2020中4.3.1a)的規(guī)定選擇隨機(jī)抽樣方法或系統(tǒng)抽樣方c)對(duì)于非均質(zhì)地表，宜考慮地表類型差異以及陰影的分布狀況，將樣點(diǎn)細(xì)分為若干相對(duì)均質(zhì)的子區(qū)域(端元),并按照GB/T39468—2020中4.3.1b)的規(guī)定選擇分層抽樣方法，對(duì)每個(gè)端元開展地面數(shù)據(jù)采集。5.3.3地表溫度相對(duì)真值的獲取方法地表溫度相對(duì)真值的獲取方法如下：a)點(diǎn)源尺度地表溫度相對(duì)真值的獲?。喊凑?.3.2的規(guī)定，在樣點(diǎn)上按照附錄A的規(guī)定獲取點(diǎn)源4GB/T41534—2022尺度地表溫度相對(duì)真值；設(shè)備采集的數(shù)據(jù)代表了觀測(cè)視場(chǎng)角內(nèi)的離地輻射，需結(jié)合配套的輔助設(shè)備獲取所需的地表溫度；b)像元尺度地表溫度相對(duì)真值的獲?。嚎紤]樣本的異質(zhì)性狀況，并結(jié)合點(diǎn)源尺度獲取的地表溫度相對(duì)真值來獲??；對(duì)于均質(zhì)地表，可直接采用單個(gè)采樣點(diǎn)獲取的地表溫度來表示，也可采用多個(gè)采樣點(diǎn)獲取的地表溫度的平均值來表示；對(duì)于非均質(zhì)地表，應(yīng)采用空間尺度擴(kuò)展的方法(例如：面積加權(quán)平均法或者陸面過程模型法等)來確定像元尺度地表溫度相對(duì)真值，具體方法見附錄B。5.3.4檢驗(yàn)操作流程檢驗(yàn)流程應(yīng)符合GB/T36296—2018中8.1的規(guī)定，主要操作流程見圖1,具體步驟如下：a)驗(yàn)證場(chǎng)地選擇：按照5.3.1的要求確定驗(yàn)證場(chǎng)地；b)地表空間異質(zhì)性分析：按照5.3.2a)的規(guī)定，利用地表溫度采集儀器分析驗(yàn)證場(chǎng)里待檢驗(yàn)遙感像元的地表空間異質(zhì)性；c)地表溫度相對(duì)真值的獲?。航Y(jié)合地表空間異質(zhì)性分析結(jié)果，按照5.3.3的規(guī)定，分別獲取點(diǎn)源尺度相對(duì)真值和像元尺度相對(duì)真值；d)時(shí)空一致匹配點(diǎn)數(shù)據(jù)獲?。焊鶕?jù)待檢驗(yàn)地表溫度遙感產(chǎn)品的獲取時(shí)間和位置，選擇對(duì)應(yīng)的地面實(shí)測(cè)的像元尺度地表溫度的相對(duì)真值；e)準(zhǔn)確度評(píng)價(jià)：在像元尺度上，按GB/T36296—2018中6.1的準(zhǔn)確度評(píng)價(jià)指標(biāo)，定量表達(dá)地表溫度遙感產(chǎn)品的準(zhǔn)確度，準(zhǔn)確度評(píng)價(jià)指標(biāo)至少應(yīng)包括均方根誤差、相關(guān)系數(shù)；f)不確定度分析：在像元尺度上，不確定度的來源分析見附錄D和GB/T39468—2020中第6章，按GB/T36296—2018中6.2規(guī)定的不確定度評(píng)價(jià)指標(biāo)定量表達(dá)地表溫度遙感產(chǎn)品的不確定度，不確定度評(píng)價(jià)指標(biāo)至少應(yīng)包括標(biāo)準(zhǔn)差。驗(yàn)證場(chǎng)地選擇驗(yàn)證場(chǎng)地選擇地表空問像元尺度相對(duì)真值時(shí)空一致匹配點(diǎn)數(shù)據(jù)狀取準(zhǔn)確度評(píng)價(jià)不確定度分析待檢驗(yàn)地表溫度遙感產(chǎn)品圖1直接檢驗(yàn)法操作流程圖5GB/T41534—20225.4間接檢驗(yàn)法交叉檢驗(yàn)法的實(shí)施條件和要求如下：a)應(yīng)具備精度已知的其他來源的地表溫度遙感產(chǎn)品作為參考對(duì)象，且其包括完備的產(chǎn)品描述、云和質(zhì)量控制等標(biāo)識(shí)；b)待檢驗(yàn)和參考對(duì)象獲取的地表溫度產(chǎn)品應(yīng)在時(shí)間、空間和觀測(cè)天頂角上存在一定交集；c)同一地理位置篩選出的交叉驗(yàn)證的匹配點(diǎn)對(duì)應(yīng)的觀測(cè)天頂角不應(yīng)大于45°,且應(yīng)具有相似的觀測(cè)條件，例如待檢驗(yàn)和參考載荷的觀測(cè)天頂角差異小于10°,觀測(cè)時(shí)間差異小于10min;d)參考載荷和待檢驗(yàn)載荷的空間分辨率有差異，且系統(tǒng)幾何定位精度可能不準(zhǔn)確時(shí)應(yīng)進(jìn)行像素聚合處理，將二種載荷反演的地表溫度遙感產(chǎn)品重采樣至更低一級(jí)的空間分辨率。檢驗(yàn)流程應(yīng)符合GB/T36296—2018中8.2.1的規(guī)定，主要操作流程見圖2,具體步驟如下：a)數(shù)據(jù)提?。菏紫雀鶕?jù)云掩模和質(zhì)量標(biāo)識(shí)，提取出待檢驗(yàn)地表溫度產(chǎn)品和參考地表溫度產(chǎn)品具有較高質(zhì)量的地表溫度數(shù)據(jù)、坐標(biāo)位置數(shù)據(jù)、觀測(cè)天頂角數(shù)據(jù)以及觀測(cè)時(shí)間數(shù)據(jù)；b)空間聚合判斷：根據(jù)待檢驗(yàn)地表溫度產(chǎn)品與參考地表溫度產(chǎn)品的空間分辨率以及系統(tǒng)幾何定c)空間聚合：如果待檢驗(yàn)地表溫度產(chǎn)品和參考地表溫度產(chǎn)品的空間分辨率不一致，或者待檢驗(yàn)地表溫度產(chǎn)品的系統(tǒng)幾何定位精度低于參考地表溫度產(chǎn)品的空間分辨率時(shí)，則應(yīng)對(duì)分別從待檢驗(yàn)和參考地表溫度產(chǎn)品提取出來的地表溫度數(shù)據(jù)、坐標(biāo)位置數(shù)據(jù)、觀測(cè)天頂角數(shù)據(jù)、觀測(cè)時(shí)間數(shù)據(jù)進(jìn)行空間聚合處理，如采用算數(shù)平均法；d)匹配點(diǎn)篩選：挑選出觀測(cè)天頂角不大于45°,且具有相似的觀測(cè)條件的匹配點(diǎn)，如觀測(cè)天頂角差小于10°,觀測(cè)時(shí)間差小于10min;f)不確定度分析：在像元尺度上，按GB/T36296—2018中6.2的不確定度評(píng)價(jià)指標(biāo)，同時(shí)考慮參考地表溫度產(chǎn)品的質(zhì)量控制標(biāo)識(shí)，定量表達(dá)地表溫度遙感產(chǎn)品的不確定度，不確定度評(píng)價(jià)指標(biāo)至少應(yīng)包括標(biāo)準(zhǔn)差。6GB/T41534—2022待檢驗(yàn)地表溫度待檢驗(yàn)地表溫度遙感產(chǎn)品數(shù)據(jù)提取否匹配點(diǎn)篩選準(zhǔn)確度評(píng)價(jià)不確定度分析參考地表溫度遙感產(chǎn)品空問聚合是圖2交叉檢驗(yàn)法操作流程圖5.4.2正向模型仿真檢驗(yàn)法正向模型仿真檢驗(yàn)法實(shí)施條件和要求如下：a)應(yīng)具備地表發(fā)射率和大氣廓線數(shù)據(jù)采集條件，且載荷過境時(shí)刻晴空無云；b)場(chǎng)地的地表類型應(yīng)簡(jiǎn)單且分布已知；c)應(yīng)利用相關(guān)儀器直接采集大氣廓線，其中大氣廓線與待檢驗(yàn)地表溫度遙感產(chǎn)品兩者獲取的時(shí)間間隔不應(yīng)超過10min,兩者距離的空間位置不應(yīng)超過10km;d)宜考慮地表類型差異，采集不同地類的發(fā)射率，像元尺度地表發(fā)射率為不同地類發(fā)射率的面積加權(quán)值；e)地表發(fā)射率的采集時(shí)刻宜接近載荷過境時(shí)刻，且宜考慮地表發(fā)射率的時(shí)效性，樣區(qū)不應(yīng)發(fā)生降雨或者林火等改變發(fā)射率數(shù)值的事件；f)宜同步觀測(cè)氣溫、氣壓、風(fēng)和濕度等常規(guī)氣象要素，以配g)應(yīng)對(duì)大氣廓線的適用性進(jìn)行質(zhì)量評(píng)估，例如通過比較具有不同吸收特征的相鄰兩通道的實(shí)測(cè)亮度溫度與模擬仿真亮度溫度之間的亮溫差值來判斷大氣廓線的可靠性；h)應(yīng)在真實(shí)性檢驗(yàn)報(bào)告中明確檢驗(yàn)所用的大氣輻射傳輸模型的模式，即逐線積分模式或者譜帶模式。5.4.2.2像元尺度地表溫度相對(duì)真值的獲取方法利用成熟的大氣輻射傳輸模型，將待檢驗(yàn)地表溫度遙感產(chǎn)品、采集的載荷過境時(shí)刻大氣廓線數(shù)據(jù)、地表發(fā)射率、衛(wèi)星觀測(cè)角度以及通道響應(yīng)函數(shù)作為正向仿真模擬的輸入?yún)?shù)，模擬出大氣層頂載荷入瞳處的亮度溫度。通過比對(duì)模擬仿真值與載荷實(shí)際觀測(cè)值，如有差異則不斷調(diào)整輸入的地表溫度值，循環(huán)迭代，直到兩者一致，調(diào)整后的最終地表溫度則是像元尺度地表溫度的相對(duì)真值。具體計(jì)算過程見7GB/T41534—2022待檢驗(yàn)地表溫度遙感產(chǎn)品迭代初值地表溫度否是否相等是大氣輻射傳輸模型圖3基于正向模型仿真檢驗(yàn)法的像元尺度地表溫度相對(duì)真值確定示意圖檢驗(yàn)流程應(yīng)符合GB/T36296—2018中8.2.2的規(guī)定，主要操作流程見圖4,具體步驟如下：a)像元尺度地表溫度相對(duì)真值獲取：采集載荷過境時(shí)刻大氣廓線和地表發(fā)射率，利用成熟的大氣輻射傳輸模型，按照5.4.2.2的規(guī)定，獲取像元尺度的地表溫度相對(duì)真值；b)大氣廓線適用性評(píng)估：通過比較實(shí)測(cè)亮度溫度與模擬仿真亮度溫度之間在具有不同吸收特征的兩相鄰?fù)ǖ赖牧翜夭畹牟钪祦砼袛啻髿饫€的可靠性，如果亮溫差的差值小于給定的閾值，如0.3K或者0.5K,則大氣廓線能準(zhǔn)確代表載荷過境時(shí)刻的大氣狀況；c)時(shí)空匹配數(shù)據(jù)選擇：根據(jù)大氣廓線適用性評(píng)估的結(jié)果，刪除大氣廓線無法準(zhǔn)確表達(dá)載荷過境時(shí)刻的數(shù)據(jù)，從待檢驗(yàn)地表溫度遙感產(chǎn)品和像元尺度地表溫度相對(duì)真值中選擇對(duì)應(yīng)的時(shí)空匹配數(shù)據(jù)；e)不確定度分析：按照5.3.4f)的要求分析不確定度。8GB/T41534—2022像元尺度地表溫度像元尺度地表溫度和對(duì)真值獲取待檢驗(yàn)地表溫度遙感產(chǎn)品時(shí)空匹配數(shù)據(jù)選擇準(zhǔn)確度評(píng)價(jià)圖4正向模型仿真檢驗(yàn)法操作流程圖時(shí)空變化趨勢(shì)分析檢驗(yàn)法的實(shí)施條件和要求如下：a)具備不同區(qū)域或者全球范圍的待檢驗(yàn)地表溫度遙感產(chǎn)品，或者具備長(zhǎng)時(shí)間序列的待檢驗(yàn)地表溫度遙感產(chǎn)品；b)具備待檢驗(yàn)地表溫度遙感產(chǎn)品時(shí)空變化趨勢(shì)的先驗(yàn)知識(shí)，或者具備與待檢驗(yàn)地表溫度有相同時(shí)間或者空間尺度，有一定相關(guān)性，且質(zhì)量有保障的其他地表特征參數(shù)，如空氣溫度、植被指檢驗(yàn)流程應(yīng)符合GB/T36296—2018中8.2.3的規(guī)定，主要操作流程見圖5,具體步驟如下：a)時(shí)空變化趨勢(shì)分析：分別對(duì)待檢地表溫度遙感產(chǎn)品和對(duì)應(yīng)的地表溫度影響因子(如空氣溫度、植被指數(shù)和葉面積指數(shù)等)進(jìn)行時(shí)間變化或空間變化趨勢(shì)分析，可按GB/T36296—2018中8.2.3的規(guī)定，在一定時(shí)空范圍內(nèi)獲取它們的時(shí)空變化規(guī)律；b)變化趨勢(shì)比對(duì)：將分析獲得的地表溫度遙感產(chǎn)品的時(shí)空變化規(guī)律與地表溫度遙感產(chǎn)品的時(shí)空變化趨勢(shì)的先驗(yàn)知識(shí)，或者與地表溫度影響因子的時(shí)空變化規(guī)律相比較，分析待檢驗(yàn)地表溫度遙感產(chǎn)品時(shí)空分布趨勢(shì)的合理性；d)不確定度分析：不確定度的來源分析見GB/T39468—2020中第6章，按GB/T36296—2018中6.2的不確定度評(píng)價(jià)指標(biāo)，定量表達(dá)地表溫度遙感產(chǎn)品的不確定度，不確定度評(píng)價(jià)指標(biāo)至少應(yīng)包括標(biāo)準(zhǔn)差。9GB/T41534—2022待檢驗(yàn)地表溫度待檢驗(yàn)地表溫度遙感產(chǎn)品村空變化趨勢(shì)分析地表溫度的時(shí)空變化趨勢(shì)影響因子的時(shí)空變化趨勢(shì)變化趨勢(shì)比對(duì)準(zhǔn)確度評(píng)價(jià)不確定度分析地表溫度的影響因子6檢驗(yàn)報(bào)告a)檢驗(yàn)報(bào)告編號(hào)；b)檢驗(yàn)報(bào)告名稱；d)檢驗(yàn)核對(duì)人；e)檢驗(yàn)簽發(fā)人；f)檢驗(yàn)單位；g)送檢單位；h)檢驗(yàn)時(shí)間。b)產(chǎn)品主要算法以及算法的優(yōu)缺點(diǎn)。b)驗(yàn)證數(shù)據(jù)集的處理與質(zhì)量控制；c)驗(yàn)證數(shù)據(jù)集的適用性。GB/T41534—2022a)檢驗(yàn)方法概述；b)檢驗(yàn)流程；c)檢驗(yàn)中的注意事項(xiàng)；d)檢驗(yàn)結(jié)果的評(píng)價(jià)指標(biāo)；e)檢驗(yàn)過程的記錄要求；f)檢驗(yàn)結(jié)果的存檔。d)對(duì)產(chǎn)品的評(píng)價(jià)與建議：通過對(duì)各項(xiàng)指標(biāo)的評(píng)價(jià)和分析，給出所檢驗(yàn)地表溫度產(chǎn)品中存在的誤對(duì)地表溫度遙感產(chǎn)品真實(shí)性檢驗(yàn)過程中的非常規(guī)問題進(jìn)行說明與描述。待檢驗(yàn)地表溫度遙感產(chǎn)品的檢驗(yàn)報(bào)告信息簡(jiǎn)表的編制見GB/T36296—2018中的附錄D,樣例見附錄E。GB/T41534—2022(規(guī)范性)點(diǎn)源尺度地表溫度相對(duì)真值獲取方法A.1基于紅外輻射計(jì)的獲取方法紅外輻射計(jì)通過探測(cè)某一特定方向的離開物體表面的紅外輻射，并結(jié)合地表發(fā)射率和大氣下行輻射，實(shí)現(xiàn)目標(biāo)物表面溫度的獲取。紅外輻射計(jì)的觀測(cè)波長(zhǎng)范圍一般為熱紅外譜段較窄的大氣窗口區(qū)，例等?；诩t外輻射計(jì)的點(diǎn)源尺度的地表溫度按公式(A.1)估算：STground=B;1{[Lground,;一(1—E;)·Ldown.;]/e;}………………(A.1)式中：獲取的點(diǎn)源尺度地表溫度，單位為開爾文(K);B-1——普朗克函數(shù)對(duì)應(yīng)窄波段i的逆函數(shù)；Lground.;——紅外輻射計(jì)測(cè)量獲取的窄波段i的紅外輻射亮度，單位為瓦每平方米球面度微米E;——窄波段i的地表發(fā)射率；Ldown.;——對(duì)應(yīng)的窄波段i的大氣下行和周圍環(huán)境下行的輻射亮度，單位為瓦每平方米球面度微米[W/(m2·sr·μm)]。A.2基于長(zhǎng)波輻射表的獲取方法長(zhǎng)波輻射表通過探測(cè)半球出射的離開物體表面的長(zhǎng)波輻射能，并結(jié)合地表發(fā)射率和大氣下行長(zhǎng)波輻射照度，實(shí)現(xiàn)目標(biāo)物表面溫度的獲取。長(zhǎng)波輻射表的觀測(cè)波長(zhǎng)范圍一般為較寬的長(zhǎng)波區(qū)，例如3.5μm～50μm。觀測(cè)視場(chǎng)范圍隨儀器架設(shè)的高度變化，一般為10m～100m不等?；陂L(zhǎng)波輻射表的點(diǎn)源尺度的地表溫度按公式(A.2)估算：STround={[Rground—(1—∈)·Rdown]/(E·σ)}1/4……(A.2)式中：長(zhǎng)波輻射表測(cè)量獲取的寬波段的目標(biāo)表面上行長(zhǎng)波輻射照度，單位為瓦每平方米E——寬波段的地表發(fā)射率；Rdown——長(zhǎng)波輻射表測(cè)量獲取的寬波段的大氣下行長(zhǎng)波輻射照度，單位為瓦每平方米(W/m2);σ-—是斯特藩-玻耳茲曼常數(shù)，取值為5.670374419×10-?W/(m2·K?)。A.3結(jié)合紅外輻射計(jì)與長(zhǎng)波輻射表的獲取方法利用紅外輻射計(jì)獲取某一特定方向的離開物體表面的紅外輻射，同時(shí)利用長(zhǎng)波輻射表獲取大氣下行長(zhǎng)波輻射照度，并利用經(jīng)驗(yàn)統(tǒng)計(jì)關(guān)系進(jìn)行紅外輻射計(jì)和長(zhǎng)波輻射表的光譜匹配，結(jié)合地表發(fā)射率來實(shí)現(xiàn)目標(biāo)物表面溫度的獲取。結(jié)合紅外輻射計(jì)與長(zhǎng)波輻射表的點(diǎn)源尺度的地表溫度按公式(A.3)估算：其中ST?見公式(A.4),ST?見公式(A.5),Lap,;見公式(A.6),Lap.;2見公式(A.7):GB/T41534—2022ST?=B;1{[Lground,i—(1—Emax)·(a·Rdown十b)]/emax}………………(A.4)ST?=B;1{[Lground,;一(1—εmin)·(a·Rdown十b)]/emin}………………(A.5)式中：ST?———按照窄通道內(nèi)最大地表發(fā)射率設(shè)定的地表溫度參考錨點(diǎn)值，單位為開爾文(K);ST?——按照窄通道內(nèi)最小地表發(fā)射率設(shè)定的地表溫度參考錨點(diǎn)值，單位為開爾文(K);Lap.i1——對(duì)應(yīng)地表溫度參考錨點(diǎn)值ST?的地表離地輻射亮度參考錨點(diǎn)值，單位為瓦每平方米球面度微米[W/(m2·sr·μm)];Lap.;2——對(duì)應(yīng)地表溫度參考錨點(diǎn)值ST?的地表離地輻射亮度參考錨點(diǎn)值，單位為瓦每平方米球面度微米[W/(m2·sr·μm)];Rdown——長(zhǎng)波輻射表測(cè)量獲取的寬波段的大氣下行長(zhǎng)波輻射照度，單位為瓦每平方米(W/m2);a——長(zhǎng)波凈輻射表寬波段的輻射照度轉(zhuǎn)換至對(duì)應(yīng)紅外輻射計(jì)窄波段i的輻射亮度的斜率；b——長(zhǎng)波凈輻射表寬波段的輻射照度轉(zhuǎn)換至對(duì)應(yīng)紅外輻射計(jì)窄波段i的輻射亮度的截距；Emax——窄波段內(nèi)的地表發(fā)射率最大值；Emin——窄波段內(nèi)的地表發(fā)射率最小值；λ?——紅外輻射計(jì)窄波段通道響應(yīng)函數(shù)的起始波長(zhǎng)，單位為微米(μm);λ——紅外輻射計(jì)窄波段通道響應(yīng)函數(shù)的截止波長(zhǎng)，單位為微米(μm);ex——窄波段內(nèi)在波長(zhǎng)λ處的地表發(fā)射率；Bx——波長(zhǎng)λ處的普朗克函數(shù)；f(λ)——紅外輻射計(jì)窄波段的通道響應(yīng)函數(shù)。GB/T41534—2022(資料性)非均質(zhì)地表像元尺度地表溫度相對(duì)真值獲取方法B.1面積加權(quán)平均方法對(duì)于非均質(zhì)地表，宜考慮地表類型差異以及陰影的分布狀況，結(jié)合高空間分辨率的遙感產(chǎn)品(如葉面積指數(shù)、植被覆蓋度和地面反射率等)或者熱像儀等設(shè)備獲取的高空間分辨率的地表亮溫?cái)?shù)據(jù)，將樣區(qū)細(xì)分為若干相對(duì)均質(zhì)的子區(qū)域(端元),采用分層抽樣方法在不同端元采集點(diǎn)源尺度地表溫度相對(duì)真值以及對(duì)應(yīng)的地表發(fā)射率，通過面積加權(quán)的方式獲取像元尺度的等效輻射亮度以及等效地表發(fā)射率，最后反算出像元尺度的地表溫度相對(duì)真值，按公式(B.1)估算：……(B.1)式中：Bx1——普朗克函數(shù)在波長(zhǎng)λ處的逆函數(shù)；N端元的個(gè)數(shù)；Ex.i——第i個(gè)端元在波段λ處的地表發(fā)射率；STground.i——獲取的第i個(gè)端元的點(diǎn)源尺度的地表溫度相對(duì)真值，單位為開爾文(K);S;——?dú)w一化后的端元組分比率；<ex>——在波段λ處像元尺度等效的地表發(fā)射率，按公式(B.2)估算：…………(B.2)B.2陸面過程模型法陸面過程模型具有嚴(yán)格物理機(jī)理，能夠?qū)㈥憵庵g能量、水分和碳交換有機(jī)耦合在一起，準(zhǔn)確的描述控制陸地表面發(fā)生的各種動(dòng)量、熱量及水分等交換的過程。通過實(shí)測(cè)的地面局地大氣強(qiáng)迫數(shù)據(jù)(如大氣溫度、相對(duì)濕度、地表風(fēng)速、地表短波輻射和地表長(zhǎng)波輻射等),并結(jié)合高空間分辨率的遙感產(chǎn)品(如葉面積指數(shù)、植被覆蓋度和地面反射率等),驅(qū)動(dòng)陸面過程模型即可對(duì)單個(gè)像元內(nèi)不同地類的溫度進(jìn)行表征，由此就可結(jié)合面積加權(quán)平均的方法獲取像元尺度地表溫度的相對(duì)真值。在采用陸面過程模型法對(duì)像元尺度地表溫度相對(duì)真值進(jìn)行估算的時(shí)候，首先需要對(duì)陸面過程模型進(jìn)行率定，以保證其模擬各個(gè)端元地表溫度的準(zhǔn)確性。GB/T41534—2022(資料性)輔助數(shù)據(jù)的獲取方法C.1地表發(fā)射率地表發(fā)射率的獲取方法主要包括以下三種：a)實(shí)驗(yàn)室測(cè)量：如根據(jù)能量守恒定律及基爾霍夫定律，將已知強(qiáng)度的輻射能投射到被測(cè)的不透明樣本表面上，并用積分球反射計(jì)測(cè)出表面反射能量，即可求得樣本的反射率，進(jìn)而計(jì)算發(fā)射b)外場(chǎng)單蓋或雙蓋發(fā)射率箱體法：如通過箱體蓋子的打開或者關(guān)閉的方式來改變環(huán)境輻射，利用多次觀測(cè)(二次測(cè)量或者四次測(cè)量)構(gòu)建的方程組來直接求解發(fā)射率；c)多光譜或者高光譜地表溫度和發(fā)射率分離法：如基于目標(biāo)發(fā)射率波譜形狀建立發(fā)射率關(guān)系模型以增加方程個(gè)數(shù)，利用多光譜或者高光譜觀測(cè)數(shù)據(jù)，實(shí)現(xiàn)地表溫度和發(fā)射率的分離，進(jìn)而獲取地表發(fā)射率。第一種方法對(duì)樣本結(jié)構(gòu)進(jìn)行了破壞，會(huì)導(dǎo)致一定的不確定度，宜優(yōu)先選擇第二種或者第三種方法。C.2大氣下行輻射亮度大氣下行輻射亮度的獲取方法主要包括以下三種：a)外場(chǎng)采用架設(shè)紅外輻射計(jì)的方法，其觀測(cè)天頂角為近似53°;b)外場(chǎng)采用反射率已知的漫反射板的方法，通過同步測(cè)量漫反射板的溫度，結(jié)合輻射傳輸模型反算大氣下行輻射亮度；c)采用大氣再分析數(shù)據(jù)(如：歐洲中期天氣預(yù)報(bào)中心，EuropeanCentreforMedium-RangeWeatherForecasts,ECMWF或者美國(guó)氣象環(huán)境預(yù)報(bào)中心，NationalCentersforEnvironmentalPrediction,NCEP等)或者大氣廓線數(shù)據(jù)，通過輻射傳輸模型進(jìn)行估算。對(duì)于長(zhǎng)期的真實(shí)性檢驗(yàn)，宜采用第一種和第三種方法；對(duì)于晴空無云以及半球無遮擋的情況，三種方法獲取的大氣下行輻射亮度精度相當(dāng)。其他情況下，由于第三種方法存在一定不確定度，因此優(yōu)先采用第一種或者第二種方法。C.3大氣下行長(zhǎng)波輻射照度大氣下行長(zhǎng)波輻射照度的獲取方法主要包括以下三種：a)外場(chǎng)對(duì)著天空架設(shè)長(zhǎng)波輻射表直接進(jìn)行測(cè)量；b)采用大氣再分析數(shù)據(jù)(如：ECMWF或者NCEP等)或者大氣廓線數(shù)據(jù)，通過輻射傳輸模型進(jìn)行估算；c)采用外場(chǎng)氣象臺(tái)站測(cè)量的大氣溫度和大氣濕度進(jìn)行估算。對(duì)于晴空無云以及半球無遮擋的情況，三種方法獲取的大氣下行輻射精度相當(dāng)。其他情況下，由于后兩種方法存在一定不確定度，因此優(yōu)先采用第一種方法。C.4大氣廓線大氣廓線的獲取方法主要包括以下三種：a)利用無線電探空儀直接進(jìn)行采集；c)通過提取大氣再分析數(shù)據(jù)的方法獲取，其中大氣再分析數(shù)據(jù)可以選擇ECMWF或者NCEP等。GB/T41534—2022(資料性)準(zhǔn)化測(cè)量規(guī)范予以減弱。預(yù)處理中幾何校正精確程度和空間均勻性程度而不同，可通過幾何精校正和選擇均勻驗(yàn)證場(chǎng)地