def block_reduce(self, block_size, how=np.ma.mean):
raster2 = block_reduce(self.raster, block_size, func=how)<EOL>geot = self.geot<EOL>geot = (geot[<NUM_LIT:0>], block_size[<NUM_LIT:0>] * geot[<NUM_LIT:1>], geot[<NUM_LIT:2>], geot[<NUM_LIT:3>], geot[<NUM_LIT:4>],<EOL>block_size[<NUM_LIT:1>] * geot[-<NUM_LIT:1>])<EOL>return GeoRaster(raster2, geot, nodata_value=self.nodata_value,projection=self.projection, datatype=self.datatype)<EOL>
geo.block_reduce(block_size, how=func)

Returns copy of raster aggregated to smaller resolution, by adding cells.
Default: func=np.ma.mean