Hi, I think I found the solution with the extract function of the raster package, which seems pretty fast from what I just experienced... But any other suggestions are welcome ! Best regards