Markus Neteler wrote > Hi, > > for the upcoming GRASS GIS course at the FOSS4G [1] I have processed a > set of Sentinel-2 scenes + a Landsat time series. > Altogether it accumulates to 90GB of compressed GeoTIFFs. > Since that's not very practical for a course, I would like to loop > over all files and cut the extent by the current region. > > Any ideas? > > Maybe it would be sufficient to pass the xmin, xmax, ymin, ymax on to > the right GDAL function? > > thanks > Markus > > PS: sure, otherwise via GDAL ... > > [1] http://2016.foss4g.org/ws27.html
what I've used recently quite a lot on the GDAL side of life: feed gdal_translate with x,y of upper left and x,y of lower right and to get automagically the raster with the extent of my interest with the same pixel size of the original input raster and also the same pixel alignment of the original input raster; without resampling, but with some kind of an off adjustment of 1-2 pixels more or less regarding the input extent of interest. it would be very nice to have such an automagic also in r.in.gdal. ----- best regards Helmut -- View this message in context: http://osgeo-org.1560.x6.nabble.com/r-in-gdal-restrict-extent-of-imported-data-set-to-current-region-tp5281223p5281259.html Sent from the Grass - Dev mailing list archive at Nabble.com. _______________________________________________ grass-dev mailing list grass-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-dev