添加链接 注册    登录
link管理
链接快照平台
  • 输入网页链接，自动生成快照
  • 标签化管理网页链接
相关文章推荐
爱看书的盒饭  ·  nsurlerrordomain ...·  3 月前    · 
茫然的烈马  ·  语法范畴：性，数，格，时，体，级……_腾讯新闻·  3 月前    · 
傻傻的刺猬  ·  Linux C ...·  6 月前    · 
重感情的酸菜鱼  ·  家用普通投影机选购指南_投影仪该怎么挑_什么值得买·  9 月前    · 
含蓄的瀑布  ·  Uninstall Flutter ...·  9 月前    · 
link管理  ›  readMeta fails to read Landsat MSS metadata (*MTL.txt) files · Issue #7 · bleutner/RStoolbox · GitHu
landsat
https://github.com/bleutner/RStoolbox/issues/7
被表白的饺子
5 月前

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement . We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Thanks bleutner for this package, it has been very helpful so far in processing Landsat scenes!

An issue I have encountered is that readMeta fails to properly load the *MTL.txt metadata files from the various Landsat MSS scenes I have tried to process.

Here is a basic example of trying to load the metadata and the error it throws: 

meta <- readMeta("./MSS/1978/LM30520251978217PAC03_MTL.txt")
Error in `rownames<-`(`*tmp*`, value = "B6_dn") : 
attempt to set 'rownames' on an object with no dimensions

and if I use the example you provide in the documentation, readMeta throws a different error because mtlFile is only an empty character vector.


mtlFile <- system.file("./MSS/1978/LM30520251978217PAC03_MTL.txt", package = "RStoolbox")
print(mtlFile)
[1] ""
meta <- readMeta(mtlFile)
Error: Metadata file does not exist. Looking for:

I have attached one of the metadata .txt files as obtained when the MSS scene was downloaded from the USGS.  Any suggestions on how to get around this issue if I am trying to run radCor and topCor on the MSS data?


LM30520251978217PAC03_MTL.txt


Thanks,

Andras

          
Hi Andras,


thanks for reporting. I think it should work now if you install the current Github version.


Out of interest: where did you download this data from? I wasn't aware that USGS are now offering reflectance scaling coefficients for MSS data, which is pretty cool. I tried to find it on EarthExplorer but couldn't.


Let me know if it works


Thank you for your work to provide the update.  readMeta does indeed work with MSS metadata now.  However, another issue popped up as I re-ran my processing script with the MSS data.  radCor only converts the first of the four MSS bands, regardless of whether any other processing has been done, like topCor or if it's directly fed from stackMeta.


For example:


mss_sref <- radCor(mss_stack, mtl, method = "costz", bandSet = "full",
                                verbose = TRUE)
Bands to convert to reflectance: B4_dn
Excluding bands: B5_dn, B6_dn, B7_dn
hazeValues was not provided -> Estimating hazeValues automatically
hazeValues estimated as: 19
Selecting atmosphere: 'moderate'
Processing radiance / reflectance

Explicitly setting the bands for bandSet makes no difference either.  I'm unsure if this might be intentional (if there is some reason that invalidates the other bands from having radiometric corrections applied), or a bug, or user error.  On a side note, I also applied topCor to the MSS data and it worked fine.


As for your question, I downloaded the MSS data from the USGS's bulk data downloader:

Landsat Data Bulk Download *login required


And for the Data Set: "Landsat Archive:L1-5 MSS"

and I provided a Scene List with the scene ID's.

          
There was a look-up mismatch with the haze look-up tables which expected bands 1 to 4, but in Landsat 1-3 these bands are designated bands 4 to 7...

It should be fixed now.

          
I have tried installing RStoolbox as a zip file from github and it installed but when I load the MTL file I get an empty "" FILE. Iam using R version 3.3.4  in windows 10. Please help

Sethi

Kind regards
 
推荐文章
爱看书的盒饭  ·  nsurlerrordomain code=-1004 error · Issue #2863 · react-native-webview/react-native-webview · GitHub
3 月前
茫然的烈马  ·  语法范畴：性，数，格，时，体，级……_腾讯新闻
3 月前
傻傻的刺猬  ·  Linux C mktime进行时间转换的陷阱 | BewareMyPower的博客
6 月前
重感情的酸菜鱼  ·  家用普通投影机选购指南_投影仪该怎么挑_什么值得买
9 月前
含蓄的瀑布  ·  Uninstall Flutter with Ease: Step-by-Step Tutorial
9 月前
Link管理   ·   51好读   ·   Sov5搜索   ·   小百科
link管理 - 链接快照平台