I have a 0.05 deg global maps (~8000 files) in GeoTiff format. It's a kind of a mask containig "1" and "0". I need to create a point shapefile for each raster with points located over "1" values of raster.
I tried to do this in QGIS via this tutorial http://www.qgistutorials.com/en/docs/sampling_raster_data.html, but creating a shapefile with ~26M of points and sampling it over each raster is quite slow.
Is there a better way to do this, preferably with gdal/ogr?