2017-06-02 4 views
0

私はhttps://cran.r-project.org/web/packages/water/vignettes/Landsat8.htmlの例でLandsat 8の例を実行しようとしています。私はエラーが発生したread.WSdataの例につきましたError in data.frame(date = unique(WSdata$date), radiation_sum = tapply(WSdata$radiation, : arguments imply differing number of rows: 1, 0自分のデータを使用しています - この例で提供されているデータではありません。read.WSdata data.frameのエラー

私のcsvファイルは、サンプルデータセット( "INTA.csv")と同じように編成されています。私がデータセット間で気づいた唯一の違いは、私のdatetimeは15分ごとで、サンプルデータセットはdatetimeが1時間ごとです。

ここに私のコードです。

`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")` 

私はまた、あなたがESPAからダウンロードした表面反射率MTLファイルは、元のようにMTLファイルの同じ名前を与えます?私たちは、元のMTLファイルを使用すると想定しておりませんか

MTLfile < - system.file( "rossfrk072616"、 "LC08_L1TP_039030_20160726_20170221_01_T1_MTL.txt"、パッケージ= "水")私はread.WSdataを実行した後、私はエラーを取得 WeatherStation <- read.WSdata(WSdata = csvfile,datetime.format = "%Y/%m/%d %H:%M",columns = c("datetime", "temp","RH", "pp", "radiation", "wind"),lat=43.07138, long= -112.4311, elev=1354.5, height= 2.5, MTL = MTLfile)

Error in data.frame(date = unique(WSdata$date), radiation_sum = 
tapply(WSdata$radiation,: arguments imply differing number of rows: 1, 0 
+0

は、MTLの詳細に話すことはできませんが、tapply:1、0:引数は、行の異なる数を意味するものではあり - 1行が呼び出されているが、そこには存在しないことを示唆しています。すなわち、列のデータ= c( "datatime"、&etc. – Chris

+0

データセットでエラーを検出できませんでした。 –

答えて

0

このエラーは日付の形式が異なるため発生します。あなたの最初の試みでは、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)

1

何らかの理由で、Webサイトからコードを取得してデータセットで機能させることができませんでした。しかし、私は次のコードで私の気象観測所のデータを読むことができました。 WeatherStation <- read.WSdata(WSdata = csvfile, date.format = "%d/%m/%Y", lat=43.07138, long= -112.4311, elev=1354.5, height= 2.5, MTL = MTLfile)