bfastts Output

[smithve87]

Hello,

I am trying to extract the real time-series data for an analysis separate from the model produced by bfast/bfastmonitor (basically I want to use only the real data to eliminate false positive breaks). However, I am having difficulty interpreting the time series object produced by bfastts. I subset my raster brick NDVI data by sensor and ran bfastts with an output that looks like this:

I wrote the output to a .csv file and opened it in excel which looked like this:

I understand the headings (Landsat scene ID) and values within the table (NDVI), but am unsure of the rows. What does the 2000.16, 2000.19, etc. indicate? Also, in the excel, are the rows supposed to be the pixels? There are about 6,000 rows, but the ncell of the raster brick is over 99,000. Also in the excel, as you scroll down (excel table attached), only certain rows will have values, which doesn't seem right since the cloud mask produces NA values, but clouds wouldn't be in the same pixels for each date. If someone could shed some light on this issue I would much appreciate it.

Best Regards,
-Vaughn
ndvi_TS.xlsx

[loicdtx]

With the code below I'm able to produce a time-series object for each pixel of the first row of the tura dataset. Rows correspond to dates in decimal form, and each column correspond to an individual pixel.
There are a lot more NAs than in the original raster object because bfastts produces a time-series with daily observations, while the original data are irregular with at least 8 days between each observation. Days filled in between observations to produce the daily time-series are automatically set to NA.

library(bfastSpatial)
data("tura")
time <- getZ(tura)
v <- t(getValues(tura, row = 1, nrows = 1))
ts <- bfastts(v, dates = time, type = 'irregular')

[sumanjnu]

@loicdtx , @bendv How do I smooth the bfmpixel generated graph, as it shows lots of ups and down? How do I save corresponding data of bfmpixel for doing the interactive graph in other software?