目錄

• • 官方關(guān)于產(chǎn)品質(zhì)量控制的說明(機(jī)翻)
  • MCTK處理QA文件
  • QA文件篩選質(zhì)量可靠像元
  • 篩選前后對比

官方關(guān)于產(chǎn)品質(zhì)量控制的說明(機(jī)翻)

質(zhì)量指標(biāo) 在生產(chǎn)過程中生成的CoreMetadata.0全局屬性QA 中的元數(shù)據(jù)對象以及質(zhì)量控制（QC）SDS中給出，或者在數(shù)據(jù)產(chǎn)品的產(chǎn)品后科學(xué)和質(zhì)量檢查中給出。CoreMetadata.0中的 QA元數(shù)據(jù)對象全局屬性是AutomaticQualityFlag和ScienceQualityFlag及其相應(yīng)的說明。根據(jù)運(yùn)行LST算法期間遇到的數(shù)據(jù)條件，根據(jù)規(guī)則設(shè)置AutomaticQualityFlag。此質(zhì)量檢查標(biāo)志的設(shè)置是完全自動(dòng)化的。設(shè)置它的規(guī)則是自由的；幾乎所有數(shù)據(jù)或中間計(jì)算都必須異常，才能將其設(shè)置為“失敗”。通常，它將設(shè)置為“通過”或“可疑”。“可疑”是指違反了執(zhí)行約束的某些界限，應(yīng)進(jìn)行進(jìn)一步分析以確定其來源。在對數(shù)據(jù)產(chǎn)品執(zhí)行自動(dòng)質(zhì)量檢查程序之后或在合格的LST調(diào)查人員檢查數(shù)據(jù)產(chǎn)品之后，可以在生產(chǎn)后設(shè)置ScienceQualityFlag。如果存在，請務(wù)必進(jìn)行檢查。將對產(chǎn)品抽樣進(jìn)行檢查。抽樣可能是隨機(jī)的，以支持野外活動(dòng)或事件驅(qū)動(dòng)。
數(shù)據(jù)產(chǎn)品中的QC SDS提供了有關(guān)每個(gè)像素算法結(jié)果的更多信息。QC SDS無符號16位數(shù)據(jù)作為位標(biāo)志存儲(chǔ)在SDS中。可以通過讀取16位無符號整數(shù)中的位來提取此QC信息。QC SDS的目的是為用戶提供有關(guān)可以在空間上下文中查看的每個(gè)像素的算法結(jié)果的信息。QC信息表明算法結(jié)果是否正常，異常，或者是否遇到了像素的其他定義條件。QC信息應(yīng)用于幫助確定LST數(shù)據(jù)對于用戶需求的有用性。表8中列出了QC SDS中的位標(biāo)志。

最關(guān)鍵的一點(diǎn)就是只有全部中間數(shù)據(jù)異常最后的溫度結(jié)果才會(huì)設(shè)置為NaN，所以根據(jù)QA文件篩選出質(zhì)量好的LST像素是非常有必要的。

MCTK處理QA文件

由于MCTK默認(rèn)將投影的影像數(shù)值類型設(shè)置為單精度浮點(diǎn)型，所以默認(rèn)投影QA文件時(shí)會(huì)使QA文件數(shù)值失去原來的意義，但MCTK的開發(fā)者肯定考慮到了這一點(diǎn)，果然查看手冊時(shí)發(fā)現(xiàn)了關(guān)于QA的處理說明。即選擇使用雙精度處理即可使QA 文件保持原始16位整形數(shù)據(jù)(十進(jìn)制)。

QA文件篩選質(zhì)量可靠像元

MCTK處理得到的QA影像文件是十進(jìn)制數(shù)值，我們還要將其轉(zhuǎn)化為二進(jìn)制才能進(jìn)一步使用。這里我選擇了QA文件0-1位為00的像元作為最后使用結(jié)果。

pro QA_dealcompile_opt idl2e = envi()ticcd,'F:\ENVI_Tempfn'qc_raster = e.openraster('qc_Swath_2D_1_georef.dat')lst_raster = e.openraster('MYD_Swath_2D_1_georef.dat')lst_data = lst_raster.getdata(bands=0)qc_10 = qc_raster.getdata(bands=0)qc_2 = string(qc_10, format='(b016)') ; 10進(jìn)制轉(zhuǎn)16位2進(jìn)制qc_goodmark = strmid(qc_2, 14) ;取出代表強(qiáng)制質(zhì)量標(biāo)志的兩位數(shù)字good = '00'qc_bad = where(~strcmp(good, qc_goodmark))lst_data[qc_bad] = 999lst_data = lst_data * float(lst_data ne 999) / (lst_data ne 999) ; 將無效999值轉(zhuǎn)換為NaNlst_fn = 'F:\ENVI_Tempfn\QC_lst_.dat'spatialref = lst_raster.spatialrefraster = e.CreateRaster(lst_fn, lst_data, SpatialRef = spatialref)raster.saveview = e.GetView()layer = view.CreateLayer(raster)toc end

篩選前后對比

處理前后對比，發(fā)現(xiàn)不可靠像元還是不少的，所以使用產(chǎn)品前一定要根據(jù)QA文件篩選出質(zhì)量可靠地像元。

總結(jié)

以上是生活随笔為你收集整理的MODIS产品质量控制文件使用方法的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。