library(reticulate) catalogclient <- import("catalogclient") cat=catalogclient$catalog$Catalog() cat$get_producttypes() date = "2017-04-20" products = cat$get_products('PROBAV_L3_S1_TOC_333M', fileformat='GEOTIFF', startdate=date, enddate=date, min_lon=0, max_lon=10, min_lat=36, max_lat=53) str(products[[1]]$files) str(products[[1]]$producttype) ndvi_files = lapply(products,function(p) p$file("NDVI")) str(ndvi_files) library(raster) str(substring(ndvi_files[1],6)) r <- raster(substring(ndvi_files[1],6)) plot(r, main='PROBA-V Tile X18 Y02')