Working With TileMatrixSets (other than WebMercator)¶

TiTiler has builtin support for serving tiles in multiple Projections by using rio-tiler and morecantile.

This Notebook shows how to use and display tiles with non-webmercator TileMatrixSet

Requirements¶

• ipyleaflet
• httpx

In [ ]:

# Uncomment if you need to install those module within the notebook
# !pip install ipyleaflet httpx

# Uncomment if you need to install those module within the notebook # !pip install ipyleaflet httpx

In [ ]:

import json

import httpx

from ipyleaflet import (
    Map,
    basemaps,
    basemap_to_tiles,
    TileLayer,
    WMSLayer,
    GeoJSON,
    projections
)

import json import httpx from ipyleaflet import ( Map, basemaps, basemap_to_tiles, TileLayer, WMSLayer, GeoJSON, projections )

In [ ]:

titiler_endpoint = "https://titiler.xyz"  # Developmentseed Demo endpoint. Please be kind.
url = "https://s3.amazonaws.com/opendata.remotepixel.ca/cogs/natural_earth/world.tif" # Natural Earth WORLD tif

titiler_endpoint = "https://titiler.xyz" # Developmentseed Demo endpoint. Please be kind. url = "https://s3.amazonaws.com/opendata.remotepixel.ca/cogs/natural_earth/world.tif" # Natural Earth WORLD tif

List Supported TileMatrixSets¶

In [ ]:

r = httpx.get(f"{titiler_endpoint}/tileMatrixSets").json()

print("Supported TMS:")
for tms in r["tileMatrixSets"]:
    print("-", tms["id"])

r = httpx.get(f"{titiler_endpoint}/tileMatrixSets").json() print("Supported TMS:") for tms in r["tileMatrixSets"]: print("-", tms["id"])

WebMercator - EPSG:3857¶

https://epsg.io/3857

In [ ]:

r = httpx.get(
    f"{titiler_endpoint}/cog/WebMercatorQuad/tilejson.json", params = {"url": url}
).json()

m = Map(center=(0, 0), zoom=2, basemap={}, crs=projections.EPSG3857)

layer = TileLayer(url=r["tiles"][0], opacity=1)
m.add_layer(layer)
m

r = httpx.get( f"{titiler_endpoint}/cog/WebMercatorQuad/tilejson.json", params = {"url": url} ).json() m = Map(center=(0, 0), zoom=2, basemap={}, crs=projections.EPSG3857) layer = TileLayer(url=r["tiles"][0], opacity=1) m.add_layer(layer) m

WGS 84 -- WGS84 - World Geodetic System 1984 - EPSG:4326¶

https://epsg.io/4326

In [ ]:

r = httpx.get(
    f"{titiler_endpoint}/cog/WorldCRS84Quad/tilejson.json", params = {"url": url}
).json()

m = Map(center=(0, 0), zoom=1, basemap={}, crs=projections.EPSG4326)

layer = TileLayer(url=r["tiles"][0], opacity=1)
m.add_layer(layer)
m

r = httpx.get( f"{titiler_endpoint}/cog/WorldCRS84Quad/tilejson.json", params = {"url": url} ).json() m = Map(center=(0, 0), zoom=1, basemap={}, crs=projections.EPSG4326) layer = TileLayer(url=r["tiles"][0], opacity=1) m.add_layer(layer) m

ETRS89-extended / LAEA Europe - EPSG:3035¶

https://epsg.io/3035

In [ ]:

r = httpx.get(
    f"{titiler_endpoint}/cog/EuropeanETRS89_LAEAQuad/tilejson.json", params = {"url": url}
).json()

my_projection = {
    'name': 'EPSG:3035',
    'custom': True, #This is important, it tells ipyleaflet that this projection is not on the predefined ones.
    'proj4def': '+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs',
    'origin': [6500000.0, 5500000.0],
    'resolutions': [
        8192.0,
        4096.0,
        2048.0,
        1024.0,
        512.0,
        256.0
    ]
}

m = Map(center=(50, 65), zoom=0, basemap={}, crs=my_projection)

layer = TileLayer(url=r["tiles"][0], opacity=1)
m.add_layer(layer)
m

r = httpx.get( f"{titiler_endpoint}/cog/EuropeanETRS89_LAEAQuad/tilejson.json", params = {"url": url} ).json() my_projection = { 'name': 'EPSG:3035', 'custom': True, #This is important, it tells ipyleaflet that this projection is not on the predefined ones. 'proj4def': '+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs', 'origin': [6500000.0, 5500000.0], 'resolutions': [ 8192.0, 4096.0, 2048.0, 1024.0, 512.0, 256.0 ] } m = Map(center=(50, 65), zoom=0, basemap={}, crs=my_projection) layer = TileLayer(url=r["tiles"][0], opacity=1) m.add_layer(layer) m

In [ ]: