# NOT RUN {
# Read a point in a Raster Brick
# define the file that has the raster brick
files <- c(system.file("extdata/raster/mod13q1/sinop-crop-ndvi.tif",
package = "sits"
))
# define the timeline
data(timeline_modis_392)
# create a data cube based on the information about the files
raster_cube <- sits_cube(
type = "BRICK", satellite = "TERRA",
sensor = "MODIS", name = "Sinop-crop",
timeline = timeline_modis_392,
bands = c("NDVI"), files = files
)
# read the time series of the point from the raster
point_ts <- sits_get_data(raster_cube,
longitude = -55.554,
latitude = -11.525
)
plot(point_ts)
# }
Run the code above in your browser using DataLab