r - 我如何在R中读取MTL文件

Question

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

Answer 1

对于从 EarthExplorer 获得的 Landsat 场景提供的标准 MTL 文件或 Glovis服务，你可以简单地做:

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和 here .

来自 MTL 的信息对于应用大气和辐射校正可能至关重要。顺便说一下，landsat软件包允许您使用 DOS() 和 radiocorr() 函数运行一些更典型的校正。

您还需要 Chander et al. (2009) 提供的标准校准值.

对于更复杂的方法 this可能是一个好的开始。