CARTO Analytics Toolbox for Snowflakeを使ってみよう

この記事はSnowflake Advent Calendar 2023シリーズ2の9日目です。
みんなでアドベントカレンダーを埋め尽くそうぜ！

はじめに

ここ1年くらいのSnowflakeは地理空間データの扱いがものすごいスピードで改良されているのですが、20年以上の歴史のあるPostGISと比べるとまだまだ足りない機能があります。
…という話をSnowVillageでしていたら「足りない関数はCARTO Analytics Toolboxにあったりするので、そっちも見るとハッピーになれるかと」というコメントが。
なにそれ知らない、ハッピーになりたい！試してみましょう。

使い方

CARTO Analytics Toolboxには2つのバージョンがあります。CARTOのアカウントなしで利用できる基本的なものがCARTO Analytics Toolbox Coreです。
これはSnowflake上のMarketplaceでNative Appとして公開されていて、無料で使用することができます。

取得ボタンを押すことでNative Appを自分のアカウントにインストールできます。
Application名はワークシートなどからCARTOの関数にアクセスするときに使うアプリケーション名、
インストールに使用したウェアハウスはインストール中に使うウェアハウスです。

1分ほどでインストールできました。

以降はワークシートからCARTOの関数を呼び出すことができます。
アプリケーション名がデフォルトのままであれば、例えば以下のように呼び出せます。

SELECT CARTO_ANALYTICS_TOOLBOX_CORE.CARTO.VERSION_CORE();

試してみよう

ではワークシートからCARTOの関数を呼び出してみましょう。
例えば東京からサンフランシスコまでの大圏コースを調べてみましょう。サンフランシスコに行きたいかー！

APIリファレンスはこちらです。

-- 東京とサンフランシスコの緯度・経度を定義してGEOGRAPHY型に変換
WITH tokyo AS (SELECT TO_GEOGRAPHY('POINT(139.6917337 35.6895014)') as geog),
     sanfrancisco AS (SELECT TO_GEOGRAPHY('POINT(-122.4192524 37.7794658)') as geog)
-- 大圏コースを計算。50個の点で近似
SELECT CARTO_ANALYTICS_TOOLBOX_CORE.CARTO.ST_GREATCIRCLE(tokyo.geog, sanfrancisco.geog, 50) AS geog
FROM tokyo, sanfrancisco;

JSON形式で出力されました。ただ、数字を見ても正しいか分かりませんね…

OpenStreetMap WKT Playgroundというサイトで地理空間データの可視化ができるのですが、入力フォーマットがWKT形式に限定されています。
WKT形式で出力されるように設定して再実行してみましょう。

-- 出力フォーマットをWKT形式に設定する
ALTER SESSION SET GEOGRAPHY_OUTPUT_FORMAT='WKT';
-- 東京とサンフランシスコの緯度・経度を定義してGEOGRAPHY型に変換
WITH tokyo AS (SELECT TO_GEOGRAPHY('POINT(139.6917337 35.6895014)') as geog),
     sanfrancisco AS (SELECT TO_GEOGRAPHY('POINT(-122.4192524 37.7794658)') as geog)
-- 大圏コースを計算。50個の点で近似
SELECT CARTO_ANALYTICS_TOOLBOX_CORE.CARTO.ST_GREATCIRCLE(tokyo.geog, sanfrancisco.geog, 50) AS geog
FROM tokyo, sanfrancisco;

出力された結果をコピペすれば、大圏コースを確認することができました！

おわりに

CARTO Analytics Toolboxには色んな関数があるので地理空間データで遊ぶのが捗りますね！
そしてNative AppによるSnowflakeの機能拡張、簡単なのにポテンシャルがすごい。