I want to extract raster value from point shapefile geometry. Raster size is 45924 x 61671 and shapefile has 11.000.000 points. I'm using
extract function from Raster package in R but is not very fast.
I'd like avoid non-terminal methods like ArcGIS or QGIS and prefer some terminal options like gdal_xx.py mode at prompt.
As a result I need an object with 11.000.000 length in any format (CSV, txt, shp, etc).
I need repeat this procedure over several similar rasters.
Here's a nice example but I don't know how to get or manipulate the resulting object and I'm not sure about their efficiency.