I have a COG, e.g. https://sentinel-cogs.s3.us-west-2.amazonaws.com/sentinel-s2-l2a-cogs/52/J/FS/2023/5/S2A_52JFS_20230501_0_L2A/TCI.tif and I want to inspect the internal tiling structure.

How can I view the extent of each of the individual tiles in the COG?

I am comfortable with Python, QGIS, Typescript etc.


2 Answers 2


If you just want to find out the internal tiling structure i.e the blocks, you can use GDAL or rasterio commandline tools:

gdalinfo /vsicurl/https://sentinel-cogs.s3.us-west-2.amazonaws.com/sentinel-s2-l2a-cogs/52/J/FS/2023/5/S2A_52JFS_20230501_0_L2A/TCI.tif

Band 1 Block=1024x1024

rio info https://sentinel-cogs.s3.us-west-2.amazonaws.com/sentinel-s2-l2a-cogs/52/J/FS/2023/5/S2A_52JFS_20230501_0_L2A/TCI.tif

{"blockxsize": 1024, "blockysize": 1024, <snip...>

If you want to visualise the blocks, you can use rasterio.DatasetReader.block_windows to access block window offsets and rasterio.windows.bounds to convert the windows to coordinates. Then make geometry from them if you wish, e.g.

import geopandas as gpd
import rasterio as rio
from rasterio.windows import bounds
from shapely.geometry import box

tif = r"https://sentinel-cogs.s3.us-west-2.amazonaws.com/sentinel-s2-l2a-cogs/52/J/FS/2023/5/S2A_52JFS_20230501_0_L2A/TCI.tif"

with rio.open(tif) as src:
    features = {"id": [], "geometry": []}
    for id, (block_index, window) in enumerate(src.block_windows(1)):
        features["geometry"].append(box(*bounds(window, src.transform)))

    gdf = gpd.GeoDataFrame(features, crs=src.crs)

enter image description here


This is a single tile. so all you need to do is use gdal to get extent:

import gdal
from gdalconst import GA_ReadOnly

data = gdal.Open(

geoTransform = data.GetGeoTransform()
minx = geoTransform[0]
maxy = geoTransform[3]
maxx = minx + geoTransform[1] * data.RasterXSize
miny = maxy + geoTransform[5] * data.RasterYSize
print(minx, miny, maxx, maxy)
data = None
  • Hi @B-C B. what do you mean its a single tile? Do you mean my example dataset is not tiled internally?
    – staf
    Commented Jul 28, 2023 at 0:04
  • Hi, the tile name is in the file name: 52JFS. these tiles are 110km by 110km, just like the one you have there. I've never heard of internal tiles from these tiles.
    – B-C B.
    Commented Jul 29, 2023 at 22:20
  • I think that refers to the sentinel tile. TIFFs are often tiled internally, and COGs make extra use of this.
    – staf
    Commented Jul 30, 2023 at 3:19
  • Awesome, I did not know that. Thanks
    – B-C B.
    Commented Jul 31, 2023 at 12:38
  • Thats OK! Thanks for your help too
    – staf
    Commented Aug 1, 2023 at 0:26

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.