目錄

• • 官方關于產品質量控制的說明(機翻)
  • MCTK處理QA文件
  • QA文件篩選質量可靠像元
  • 篩選前后對比

官方關于產品質量控制的說明(機翻)

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

最關鍵的一點就是只有全部中間數據異常最后的溫度結果才會設置為NaN，所以根據QA文件篩選出質量好的LST像素是非常有必要的。

MCTK處理QA文件

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

QA文件篩選質量可靠像元

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

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進制轉16位2進制qc_goodmark = strmid(qc_2, 14) ;取出代表強制質量標志的兩位數字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值轉換為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

篩選前后對比

處理前后對比，發現不可靠像元還是不少的，所以使用產品前一定要根據QA文件篩選出質量可靠地像元。

總結

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

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。