def split_metatiles(output_dir, metatiles):
# Determine metatile size and bit depth.
mt_ds = gdal.Open(metatiles[0])
assert mt_ds.RasterXSize % UF_TILE_SIZE == 0
assert mt_ds.RasterYSize % UF_TILE_SIZE == 0
x_uf_count = mt_ds.RasterXSize / UF_TILE_SIZE
y_uf_count = mt_ds.RasterYSize / UF_TILE_SIZE
for ytile in range(y_uf_count):
print 'Process Row: ', ytile
split_row_of_unit_fields(output_dir, metatiles, ytile)
logging.info('Wrote %d piles to directory %s.',
x_uf_count * y_uf_count, output_dir)
评论列表
文章目录