Question

我正试图在https://cran.r-project.org/web/packages/water/vignettes/Landsat8.html示例中执行Landsat 8示例。我陷入read.WSdata示例，我得到错误DISTINCT我正在使用自己的数据 - 而不是示例中提供的数据。

我的csv文件的组织方式与示例数据集完全相同（＆＃34; INTA.csv＆＃34;）。我在数据集之间注意到的唯一区别是我的每15分钟有一个日期时间，示例数据集每小时都有一个日期时间。

这是我的代码。

DELETE

我还假设我们使用的是原始MTL文件，而不是表面反射MTL文件，当您从ESPA下载时，它会提供与原始文件相同的mtl文件名称？

MTLfile＆lt; - system.file（＆＃34; rossfrk072616＆＃34;，＆＃34; LC08_L1TP_039030_20160726_20170221_01_T1_MTL.txt＆＃34;，package =＆＃34; water＆＃34;） Error in data.frame(date = unique(WSdata$date), radiation_sum = tapply(WSdata$radiation, : arguments imply differing number of rows: 1, 0

运行read.WSdata后，我收到错误

`rm(list=ls()) 
library(water)
aoi<-createAoi(topleft=c(385387,4776577),
bottomright=c(414825,4749526), EPSG = 32612)
raw_data_folder <- system.file("rossfrk072616", package="water")
image <- loadImage(path=raw_data_folder, aoi=aoi, sat="L8")
image.SR <- loadImageSR(path=raw_data_folder, aoi=aoi)
plot(image)
plot(image.SR)
csvfile<-system.file("rossfrk072616","FTHI_L8_1.csv",package="water")`

Answer 1

出于某种原因，我无法从网站获取代码以使用我的数据集。但是，我能够使用以下代码读取我的气象站数据。 WeatherStation <- read.WSdata(WSdata = csvfile, date.format = "%d/%m/%Y", lat=43.07138, long= -112.4311, elev=1354.5, height= 2.5, MTL = MTLfile)

Answer 2

这是错误与日期的不同格式有关。在您第一次尝试时，date.format设置为'%Y/%m/%d'。

此外，您可以直接在read.WSdata()函数中指定文件，例如：

WeatherStation <- read.WSdata(WSdata = 'FTHI_L8_1.csv', date.format = "%d/%m/%Y", lat=43.07138, long= -112.4311, elev=1354.5, height= 2.5, MTL = MTLfile)

Answer 3

我使用了数据，并且效果很好

library(water) 

aoi <- createAoi(topleft = c(810927, 2134059), bottomright = c( 272751,1985845), 
                     EPSG = 32616)
plot(aoi)

csvfile <- system.file("extdata", "datos.csv", package="water")

MTLfile <- system.file("extdata", "L8.MTL.txt", package="water")


ws<- read.WSdata(WSdata = csvfile, date.format = "%d/%m/%Y", time.format="%H:%M:%S", cf= 
                  c(1,1,1),lat=18.094, long= -89.462, elev=279, height= 2,  MTL = 
                  MTLfile, columns=c("date" = 1, "time" = 2, "radiation" = 3,"wind" = 4, 
                  "RH" = 5, "temp" = 6, "rain" = 7))

read.WSdata data.frame中的错误

3 个答案: