542 """Size of natural block reported by GDAL for bands of given dataset"""
545 for i
in range(band_from, band_to):
546 band = ds.GetRasterBand(i)
547 assert band
is not None,
"Cannot access raster band %d" % i
548 dims = band.GetBlockSize()
555 if block_dims != dims:
556 logit(
"MSG: Block sizes don't match: %s != %s\n" % (str(block_dims), str(dims)))
558 assert block_dims
is not None,
"Failed to calculate block size"
559 return (int(block_dims[0]), int(block_dims[1]))
def calculate_block_size(ds, band_from, band_to)