我如何在R中读取MTL文件

我对 R 编程非常陌生，请有人告诉我如何读取用陆地卫星数据存档的 MTL 文件。

对于随 Landsat 场景提供的标准 MTL 文件，该文件来自地球探索号 http://earthexplorer.usgs.gov/ or Glovis http://glovis.usgs.gov/服务，您可以简单地执行以下操作：

mtl <- read.delim('L71181068_06820100518_MTL.txt', sep = '=', stringsAsFactors = F)

所以，对于这样开始的事情：

GROUP = L1_METADATA_FILE  GROUP = METADATA_FILE_INFO...

你可以使用这个：

> mtl[grep("LMAX",mtl$GROUP),]
                  GROUP L1_METADATA_FILE
64          LMAX_BAND1           293.700
66          LMAX_BAND2           300.900
68          LMAX_BAND3           234.400
70          LMAX_BAND4           241.100
72          LMAX_BAND5            47.570
74         LMAX_BAND61            17.040
76         LMAX_BAND62            12.650
78          LMAX_BAND7            16.540
80          LMAX_BAND8           243.100
84       QCALMAX_BAND1             255.0
86       QCALMAX_BAND2             255.0
88       QCALMAX_BAND3             255.0
90       QCALMAX_BAND4             255.0
92       QCALMAX_BAND5             255.0
94      QCALMAX_BAND61             255.0
96      QCALMAX_BAND62             255.0
98       QCALMAX_BAND7             255.0
100      QCALMAX_BAND8             255.0

每个服务都有提供的字典，找到here https://lta.cr.usgs.gov/landsat_dictionary.html and here http://glovis.usgs.gov/ImgViewer/landsat_dictionary.html.

来自 MTL 的信息对于应用大气和辐射校正可能至关重要。顺便一提，landsat http://cran.r-project.org/web/packages/landsat/index.html包允许您使用以下命令运行一些更典型的校正DOS() and radiocorr()功能。

您还需要以下提供的标准校准值钱德等人。 (2009) http://landsathandbook.gsfc.nasa.gov/pdfs/Landsat_Calibration_Summary_RSE.pdf.

对于更复杂的方法this http://www.fs.fed.us/pnw/pubs/journals/pnw_2006_schroeder001.pdf也许是一个好的开始。