dayjournal | Try #107 – TiTilerでCloud Optimized GeoTIFF(COG)を利用した動的タイルサーバーを構築してみた

この記事は、「MIERUNE Advent Calendar 2023」の13日目の記事です。

TiTilerでCloud Optimized GeoTIFF(COG)を利用した動的タイルサーバーを構築してみました！

近年、海外のFOSS4Gコミュニティで注目を集めているTiTilerを試してみました。

TiTilerはFastAPIとRasterio/GDALを基盤として構築された動的なタイルサーバーです。その主な特徴には、Cloud Optimized GeoTIFF(COG)のサポート・複数の投影法のサポート・様々な出力フォーマット(JPEG、JP2、PNG、WEBP、GTIFF、NumpyTile)・WMTSのサポート・バーチャルモザイクのサポートなどが含まれます。また、AWS CDKを利用したLambdaやECSのデプロイ環境も提供されています。

今後、日本でも利用が広がるツールになってくると思いますので、ぜひチェックしてみてください！

事前準備

今回利用するサンプルデータは、AWSが提供しているRegistry of Open Data on AWSから取得できるCloud Optimized GeoTIFF(COG)形式です。

「Sentinel-2」のCOGを選択しました。
https://sentinel-cogs.s3.us-west-2.amazonaws.com/sentinel-s2-l2a-cogs/54/T/WN/2023/8/S2A_54TWN_20230815_0_L2A/TCI.tif

TiTilerのインストール & 起動

TiTilerをインストールします。今回はDockerを利用しローカルで起動します。

TiTilerのリポジトリをクローンします。

git clone https://github.com/developmentseed/titiler.git
cd titiler

TiTilerを起動します。

docker-compose up --build titiler

起動後に「 http://localhost:8000 」にアクセスします。TiTilerのトップページが表示されるとOKです！

APIドキュメントを見ると実際の利用方法が確認できます。

データの基本情報取得

はじめに、データの基本情報を取得します。

「/cog/info」の項目選択 → Dataset URLにCOGのURL入力 → 「Execute」をクリック。

対象データの基本情報がJSONで取得できます。

{
  "bounds": [
    140.99975320947394,
    42.356322389620374,
    142.3544473299751,
    43.35285539188815
  ],
  "minzoom": 7,
  "maxzoom": 13,
  "band_metadata": [
    [
      "b1",
      {}
    ],
    [
      "b2",
      {}
    ],
    [
      "b3",
      {}
    ]
  ],
  "band_descriptions": [
    [
      "b1",
      ""
    ],
    [
      "b2",
      ""
    ],
    [
      "b3",
      ""
    ]
  ],
  "dtype": "uint8",
  "nodata_type": "Nodata",
  "colorinterp": [
    "red",
    "green",
    "blue"
  ],
  "driver": "GTiff",
  "count": 3,
  "width": 10980,
  "height": 10980,
  "overviews": [
    2,
    4,
    8,
    16
  ],
  "nodata_value": 0
}

ビューアーでデータ表示

次に、ビューアーで対象データを表示します。TiTilerには標準でデータを表示するビューアーが用意されています！

「/cog/viewer」の項目を選択します。そして「 http://localhost:8000/cog/viewer 」に直接アクセスすることでビューアー画面が表示されます。

ビューアーにアクセス → COGのURLを入力 → 「Apply」をクリック。

対象のデータがビューアー上に表示されます。ズームイン・アウトや各パラメーターの設定を試すこともできます。

GeoJSON形式の基本情報取得

次に、GeoJSON形式の基本情報を取得しQGISで表示します。

「/cog/info.geojson」のAPIを利用しGeoJSON形式のデータをQGISで読み込みます。

http://localhost:8000/cog/info.geojson?url=https://sentinel-cogs.s3.us-west-2.amazonaws.com/sentinel-s2-l2a-cogs/54/T/WN/2023/8/S2A_54TWN_20230815_0_L2A/TCI.tif

対象データの範囲ポリゴンと属性に基本情報が入っているのを確認できます。

XYZ Tilesを取得

次に、XYZ Tiles形式のラスタタイルを取得しQGISで表示します。

「/cog/tiles」のAPIを利用しXYZ Tiles形式のラスタタイルをQGISで読み込みます。

http://localhost:8000/cog/tiles/WebMercatorQuad/{z}/{x}/{y}.png?url=https://sentinel-cogs.s3.us-west-2.amazonaws.com/sentinel-s2-l2a-cogs/54/T/WN/2023/8/S2A_54TWN_20230815_0_L2A/TCI.tif

対象データがXYZ Tilesで動的に配信されているのを確認できます。

WMTSを取得

最後に、WMTS形式のラスタタイルを取得しQGISで表示します。

「/cog/WMTSCapabilities.xml」のAPIを利用しWMTS形式のラスタタイルをQGISで読み込みます。

http://localhost:8000/cog/WMTSCapabilities.xml?url=https://sentinel-cogs.s3.us-west-2.amazonaws.com/sentinel-s2-l2a-cogs/54/T/WN/2023/8/S2A_54TWN_20230815_0_L2A/TCI.tif