地理空間データ型¶

GEOGRAPHY データ型¶

GEOGRAPHY データ型は WGS 84標準に従います（空間参照 ID 4326。詳細については、https://epsg.io/4326を参照）。

地球上のポイントは、経度（-180度から+180度）および緯度（-90から+90）で表されます。Snowflakeは、小数点14桁を使用して GEOGRAPHY の座標を格納します。データにこの制限を超える小数点以下の桁数が含まれている場合、座標は指定された長さ制約に準拠するように丸められます。

現在、高度はサポートされていません。

線分は、地球表面の測地弧として解釈されます。

Snowflakeは、 GEOGRAPHY データ型を操作する 地理空間関数 も提供します。

地理空間データ（例: 緯度経度データ、 WKT、 WKB、 GeoJSON など）がある場合は、このデータを元の形式のまま VARCHAR、 VARIANT、 NUMBER 列に保存するのではなく、 GEOGRAPHY 列に変換して保存することをお勧めします。データを GEOGRAPHY 列に保存すると、地理空間機能を使用するクエリのパフォーマンスが大幅に向上します。

GEOMETRY データ型¶

GEOMETRY データ型は、平面（ユークリッド、デカルト）座標系の機能を表します。

座標は、実数（x、y）のペアとして表されます。現在、2D座標のみがサポートされています。

XとYの単位は、 GEOMETRY オブジェクトに関連付けられた 空間参照系（SRS） によって決定されます。空間参照系は、 空間参照系識別子（SRID） 番号によって識別されます。GEOMETRY オブジェクトの作成時に、または ST_SETSRID の呼び出しによって、 SRID が提供される場合を除き、 SRID は0です。

Snowflakeは、小数点14桁を使用して GEOMETRY の座標を格納します。データにこの制限を超える小数点以下の桁数が含まれている場合、座標は指定された長さ制約に準拠するように丸められます。

Snowflakeは、 GEOMETRY データ型で動作する地理空間関数 のセットを提供します。これらの関数に対して、

• ジオメトリが非平面 SRS を使用している場合でも、すべての関数は平面座標を想定しています。

• 測定関数（例: ST_LENGTH）では、座標系と同じ単位が用いられます。

• 引数として複数の GEOMETRY 式を受け付ける関数（例: ST_DISTANCE）の場合、入力式は同じ SRS で定義する必要があります。

サポートされている標準の入出力形式¶

GEOGRAPHY および GEOMETRY データ型は、次の業界標準の入出力形式をサポートしています。

• 周知のテキスト (WKT)

• 周知のバイナリ (WKB)

• 拡張された WKT と WKB （EWKT と EWKB） （ EWKT と EWKB の取り扱いに関する注意 を参照）

• IETF GeoJSON （GeoJSON の取り扱いに関する注意 を参照）

以下のOpen Geospatial ConsortiumのSimple Feature Accessのリファレンスも参考になる可能性があります。

• 共通アーキテクチャ

• SQL オプション

これらの標準から逸脱する場合については、Snowflakeのドキュメントに明記されています。

サポートされている地理空間オブジェクト型¶

GEOGRAPHY および GEOMETRY のデータ型は、次の型の地理空間オブジェクトを保存できます。

• WKT / WKB / EWKT / EWKB / GeoJSON 地理空間オブジェクト:

  • ポイント

  • MultiPoint

  • LineString

  • MultiLineString

  • ポリゴン

  • MultiPolygon

  • GeometryCollection

• これらの GeoJSON 固有の地理空間オブジェクト:

  • 機能

  • FeatureCollection

結果セットの出力形式の指定¶

セッションパラメーター GEOGRAPHY_OUTPUT_FORMAT と GEOMETRY_OUTPUT_FORMAT は、結果セットの GEOGRAPHY と GEOMETRY 列のレンダリングを制御します（それぞれ）。

これらのパラメーターは、以下の値のいずれかとなります。

EWKT および EWKB の場合、出力では SRID は常に4326です。 GEOGRAPHY 値用の EWKT および EWKB 処理 をご参照ください。

このパラメーターは、Snowflake UI や SnowSQL コマンドラインクライアント、 JDBC、 ODBC、Node.js、Pythonなどのドライバーやコネクタを含むすべてのクライアントに影響します。

たとえば、 JDBC ドライバーは、 GEOGRAPHY 型の結果列（この例では列 i）に対して次のメタデータを返します。

• GEOGRAPHY_OUTPUT_FORMAT='GeoJSON' または GEOMETRY_OUTPUT_FORMAT='GeoJSON' の場合、

  • ResultSetMetaData.getColumnType(i) は、 java.sql.Types.VARCHAR を返します。

  • ResultSetMetaData.getColumnClassName(i) は、 "java.lang.String" を返します。

• GEOGRAPHY_OUTPUT_FORMAT='WKT' または 'EWKT' の場合、または GEOMETRY_OUTPUT_FORMAT='WKT' または 'EWKT' の場合、

  • ResultSetMetaData.getColumnType(i) は、 java.sql.Types.VARCHAR を返します。

  • ResultSetMetaData.getColumnClassName(i) は、 "java.lang.String" を返します。

• GEOGRAPHY_OUTPUT_FORMAT='WKB' または 'EWKB' の場合、または GEOMETRY_OUTPUT_FORMAT='WKB' または 'EWKB' の場合、

  • ResultSetMetaData.getColumnType(i) は、 java.sql.Types.BINARY を返します。

  • ResultSetMetaData.getColumnClassName(i) は、 "[B" （バイトの配列）を返します。

GEOGRAPHY データの挿入とクエリの例¶

以下のコードは、 GEOGRAPHY データ型のサンプルの入力と出力を示しています。次の点に注意してください。

• WKT、 EWKT、 GeoJSON の座標では、緯度の前に経度が表示されます（例: POINT(lon lat)）。

• WKB および EWKB 出力では、 BINARY_OUTPUT_FORMAT パラメーターが HEX （パラメーターのデフォルト値）に設定されていると想定しています。

次の例では、 GEOGRAPHY 列のテーブルを作成し、 WKT 形式でデータを挿入し、さまざまな出力形式でデータを返します。

CREATE OR REPLACE TABLE geospatial_table (id INTEGER, g GEOGRAPHY);
INSERT INTO geospatial_table VALUES
  (1, 'POINT(-122.35 37.55)'),
  (2, 'LINESTRING(-124.20 42.00, -120.01 41.99)');

ALTER SESSION SET GEOGRAPHY_OUTPUT_FORMAT='GeoJSON';

SELECt g
  FROM geospatial_table
  ORDER BY id;

+------------------------+
| G                      |
|------------------------|
| {                      |
|   "coordinates": [     |
|     -122.35,           |
|     37.55              |
|   ],                   |
|   "type": "Point"      |
| }                      |
| {                      |
|   "coordinates": [     |
|     [                  |
|       -124.2,          |
|       42               |
|     ],                 |
|     [                  |
|       -120.01,         |
|       41.99            |
|     ]                  |
|   ],                   |
|   "type": "LineString" |
| }                      |
+------------------------+

ALTER SESSION SET GEOGRAPHY_OUTPUT_FORMAT='WKT';

SELECt g
  FROM geospatial_table
  ORDER BY id;

+-------------------------------------+
| G                                   |
|-------------------------------------|
| POINT(-122.35 37.55)                |
| LINESTRING(-124.2 42,-120.01 41.99) |
+-------------------------------------+

ALTER SESSION SET GEOGRAPHY_OUTPUT_FORMAT='WKB';

SELECt g
  FROM geospatial_table
  ORDER BY id;

+------------------------------------------------------------------------------------+
| G                                                                                  |
|------------------------------------------------------------------------------------|
| 01010000006666666666965EC06666666666C64240                                         |
| 010200000002000000CDCCCCCCCC0C5FC00000000000004540713D0AD7A3005EC01F85EB51B8FE4440 |
+------------------------------------------------------------------------------------+

ALTER SESSION SET GEOGRAPHY_OUTPUT_FORMAT='EWKT';

SELECt g
  FROM geospatial_table
  ORDER BY id;

+-----------------------------------------------+
| G                                             |
|-----------------------------------------------|
| SRID=4326;POINT(-122.35 37.55)                |
| SRID=4326;LINESTRING(-124.2 42,-120.01 41.99) |
+-----------------------------------------------+

ALTER SESSION SET GEOGRAPHY_OUTPUT_FORMAT='EWKB';

SELECt g
  FROM geospatial_table
  ORDER BY id;

+--------------------------------------------------------------------------------------------+
| G                                                                                          |
|--------------------------------------------------------------------------------------------|
| 0101000020E61000006666666666965EC06666666666C64240                                         |
| 0102000020E610000002000000CDCCCCCCCC0C5FC00000000000004540713D0AD7A3005EC01F85EB51B8FE4440 |
+--------------------------------------------------------------------------------------------+

GEOMETRY で異なる SRIDs を使用した場合の効果について¶

GEOMETRY 列には、異なる SRIDs を持つオブジェクトを挿入できます。列に複数の SRID が含まれている場合は、重要なパフォーマンス最適化の一部が適用されません。これにより、特に地理空間述語で結合する場合に、クエリが遅くなる可能性があります。

GEOMETRY オブジェクトの空間参照システム（SRS）および SRID の変更¶

既存の GEOMETRY オブジェクトの SRS と SRID を変更するには、 ST_TRANSFORM 関数を呼び出し、新しい SRID を渡します。この関数は、SRS を使用するために変換された新しい SRID と座標を持つ、新しい GEOMETRY オブジェクトを返します。たとえば、SRID 32633の SRS を使用する geometry_expression の GEOMETRY オブジェクトを返すには、以下のステートメントを実行します。

SELECT ST_TRANSFORM(geometry_expression, 32633);

既存の GEOMETRY オブジェクトで元の SRID が正しく設定されていない場合は、追加の引数として元の SRID を指定します。たとえば、 geometry_expression が SRID 4326を使用する GEOMETRY オブジェクトで、これを SRID 28992を使用するように変換する場合は、以下のステートメントを実行します。

SELECT ST_TRANSFORM(geometry_expression, 4326, 28992);

GEOMETRY オブジェクトが SRS に対して正しい座標を使用しているにもかかわらず、 SRID が間違っている場合は、 ST_SETSRID 関数を呼び出すと SRID を修正できます。たとえば、次のステートメントは、座標を変更せずに geometry_expression の SRID を4326に設定します。

SELECT ST_SETSRID(geometry_expression, 4326);

GEOGRAPHY 列および GEOMETRY 列での DML 操作の実行¶

GEOGRAPHY または GEOMETRY の列が DML 操作のターゲットである場合（INSERT、 COPY、 UPDATE、 MERGE、または CREATE TABLE AS ...）、列のソース式は、次の型のいずれかです。

• GEOGRAPHY または GEOMETRY: 型 GEOGRAPHY または GEOMETRY の式は通常、解析関数、コンストラクター関数、または既存の GEOGRAPHY または GEOMETRY の列の結果です。サポートされている関数と関数のカテゴリの包括的なリストについては、 地理空間関数 をご参照ください。

• VARCHAR: WKT、 WKB （16進数形式）、 EWKT、 EWKB （16進数形式）、または GeoJSON 形式の文字列として解釈されます（TO_GEOGRAPHY(VARCHAR) を参照）。

• BINARY: WKB バイナリとして解釈されます（TO_GEOGRAPHY(BINARY) および

• TO_GEOMETRY(BINARY) を参照）。

• VARIANT: GeoJSON オブジェクトとして解釈されます（TO_GEOGRAPHY(VARIANT) および TO_GEOMETRY(VARIANT) を参照）。

ステージからの地理空間データのロード¶

ステージ内の CSV または JSON/AVRO ファイルからのデータは、直接（つまり、コピー変換なしで） GEOGRAPHY 列にロードできます。

• CSV ：対応する CSV 列の文字列値は、 GeoJSON、 WKT、 EWKT、 WKB、または EWKB として解析されます（ TO_GEOGRAPHY （VARCHAR） を参照）。

• JSON/AVRO: ファイルの JSON 値は GeoJSON として解釈されます（TO_GEOGRAPHY （VARIANT） を参照）。

  GEOGRAPHY 値用の GeoJSON 処理 もご参照ください。

他のファイル形式（Parquet、 ORC など）からのデータのロードは、 COPY 変換で可能です。

Java UDFs での地理空間データの使用¶

Java UDFs は、 GEOGRAPHY 型を引数および戻り値として許可します。詳細については SQL-Javaデータ型マッピング および GEOGRAPHY 値のインラインJava UDF への引き渡し をご参照ください。

JavaScript UDFs での地理空間データの使用¶

JavaScript UDFs は、 GEOGRAPHY または GEOMETRY 型を引数および戻り値として許可します。

JavaScript UDF に型 GEOGRAPHY または GEOMETRY の引数がある場合、その引数は UDF 本文内で JSON オブジェクトとして GeoJSON 形式で表示されます。

JavaScript UDF が GEOGRAPHY または GEOMETRY を返す場合、 UDF 本体は JSON オブジェクトを GeoJSON 形式で返すことが期待されています。

たとえば、次の2つの JavaScript UDFs は、組み込み関数 ST_X および ST_MAKEPOINT とほぼ同等です。

CREATE OR REPLACE FUNCTION my_st_x(g GEOGRAPHY) RETURNS REAL
LANGUAGE JAVASCRIPT
AS
$$
  if (G["type"] != "Point")
  {
     throw "Not a point"
  }
  return G["coordinates"][0]
$$;

CREATE OR REPLACE FUNCTION my_st_makepoint(lng REAL, lat REAL) RETURNS GEOGRAPHY
LANGUAGE JAVASCRIPT
AS
$$
  g = {}
  g["type"] = "Point"
  g["coordinates"] = [ LNG, LAT ]
  return g
$$;

Python UDFs での地理空間データの使用¶

Python UDFs は、 GEOGRAPHY および GEOMETRY 型を引数および戻り値として許可します。

Python UDF が GEOGRAPHY または GEOMETRY 型の引数を持つ場合、その引数は GeoJSON オブジェクトとして表現され、 UDF 本文の中でPython dict オブジェクトに変換されます。

Python UDF が GEOGRAPHY または GEOMETRY を返す場合、 UDF 本文は GeoJSON の構造に準拠したPython dict オブジェクトを返すことが予測されます。

たとえば、このPython UDF は、複合 GEOGRAPHY 型を構成する明確なジオメトリの数を返します。

CREATE OR REPLACE FUNCTION py_numgeographys(geo GEOGRAPHY)
RETURNS INTEGER
LANGUAGE PYTHON
RUNTIME_VERSION = 3.10
PACKAGES = ('shapely')
HANDLER = 'udf'
AS $$
from shapely.geometry import shape, mapping
def udf(geo):
    if geo['type'] not in ('MultiPoint', 'MultiLineString', 'MultiPolygon', 'GeometryCollection'):
        raise ValueError('Must be a composite geometry type')
    else:
        g1 = shape(geo)
        return len(g1.geoms)
$$;

Python UDFs のさらなるサンプルについては、 Snowflakeラボ をご確認ください。中には、複雑な空間操作を可能にしたり、データのインジェスチョンを簡素化したりするものもあります。たとえば、 この UDF は、シェープファイル（.SHP）、 TAB、 KML、 GPKG など、ネイティブでサポートされていない形式を読み込むことができます。

H3での GEOGRAPHY オブジェクトの使用¶

H3 は 階層的な地理空間インデックス で、世界を 離散的なグローバルグリッドシステム の六角形のセルに分割します。

Snowflakeには SQL 関数が用意されており、 GEOGRAPHY オブジェクトでH3を使用できます。これらの関数は次のように使用できます。

• ポイントを表す GEOGRAPHY オブジェクトのH3セル ID （インデックス）を取得します（逆も同様）。

• オブジェクト GEOGRAPHY をカバーするH3セルの最小セットの IDs を取得します。

• ポリゴンを表す GEOGRAPHY オブジェクト内に重心を持つH3セルの IDs を取得します。

• H3セルの境界を表す GEOGRAPHY オブジェクトを取得します。

• 特定のH3セルの親と子を取得します。

• H3セルの重心の経度と緯度を取得します（逆も同様）。

• H3セルの 解像度 を取得します。

• H3セル ID の16進数表現を取得します（逆も同様）。

これらの関数の詳細については、 地理空間関数 をご参照ください。

GEOGRAPHY と GEOMETRY の違いについて¶

GEOGRAPHY と GEOMETRY のデータ型両方とも地理空間機能を定義しますが、これらの型は異なるモデルを使用します。次のテーブルは、相違点をまとめたものです。

GEOGRAPHY と GEOMETRY のデータ型を比較する例¶

次の例では、 GEOGRAPHY と GEOMETRY のデータ型を入力として使用した場合の地理空間関数の出力を比較しています。

SELECT ST_DISTANCE(
    ST_POINT(13.4814, 52.5015),
    ST_POINT(-121.8212, 36.8252))
  AS distance_in_meters;

+--------------------+
| DISTANCE_IN_METERS |
|--------------------|
|   9182410.99227821 |
+--------------------+

SELECT ST_DISTANCE(
    ST_GEOMPOINT(13.4814, 52.5015),
    ST_GEOMPOINT(-121.8212, 36.8252))
  AS distance_in_degrees;

+---------------------+
| DISTANCE_IN_DEGREES |
|---------------------|
|       136.207708844 |
+---------------------+

SELECT ST_AREA(border) AS area_in_sq_meters
  FROM world_countries
  WHERE name = 'Germany';

+-------------------+
| AREA_IN_SQ_METERS |
|-------------------|
|  356379183635.591 |
+-------------------+

SELECT ST_AREA(border) as area_in_sq_degrees
  FROM world_countries_geom
  WHERE name = 'Germany';

+--------------------+
| AREA_IN_SQ_DEGREES |
|--------------------|
|       45.930026848 |
+--------------------+

例3: ベルリンからサンフランシスコまでの線に重なる国名のクエリ¶

次のテーブルは、 GEOGRAPHY 型と GEOMETRY 型の ST_INTERSECTS の出力を比較しています。

. GEOGRAPHY 入力を使用した ST_INTERSECTS	. GEOMETRY 入力を使用した ST_INTERSECTS
SELECT name FROM world_countries WHERE ST_INTERSECTS(border, TO_GEOGRAPHY( 'LINESTRING(13.4814 52.5015, -121.8212 36.8252)' )); Copy +--------------------------+ | NAME                     | |--------------------------| |                  Germany | |                  Denmark | |                  Iceland | |                Greenland | |                   Canada | | United States of America | +--------------------------+	SELECT name FROM world_countries_geom WHERE ST_INTERSECTS(border, TO_GEOMETRY( 'LINESTRING(13.4814 52.5015, -121.8212 36.8252)' )); Copy +--------------------------+ | NAME                     | |--------------------------| |                  Germany | |                  Belgium | |              Netherlands | |           United Kingdom | | United States of America | +--------------------------+

入力データ検証の違いについて¶

入力形状の GEOMETRY または GEOGRAPHY オブジェクトを作成するには、 `Simple Featuresの OGC 規則`__ に従って、整形式で有効な形状を使用する必要があります。次のセクションでは、入力データの有効性が GEOMETRY と GEOGRAPHY でどのように異なるかについて説明します。

有効な GEOGRAPHY の場合でも無効な GEOMETRY の形状¶

指定された形状が、有効な GEOGRAPHY オブジェクトでも、無効な GEOMETRY オブジェクトである可能性があります（逆も同様）。

たとえば、自己交差ポリゴンは OGC 規則では許可されていません。与えられた点のセットは、直交領域では交差しても、球面上では交差しない辺を定義する場合があります。次のポリゴンを考えてみましょう。

POLYGON((0 50, 25 50, 50 50, 0 50))

Copy

デカルト領域では、この多角形は線に劣化し、結果として無効になります。

しかし、球面上では、この同じ多角形はそれ自身と交差せず、有効です。

無効な形状が地理空間関数に与える影響について¶

地理空間関数 が異なれば、無効な形状に対して GEOGRAPHY または GEOMETRY オブジェクトを渡す場合の影響も異なります。

無効な形状の操作¶

次のセクションでは、関数が無効な形状を作成できるようにする方法と、GEOGRAPHY または GEOMETRY オブジェクトが無効な形状または修復された形状を表しているかどうかを判断する方法について説明します。

TO_GEOGRAPHY( <input> [, <allowInvalid> ] )

ST_GEOGFROMWKB( <input> [, <allowInvalid> ] )

ST_GEOGFROMWKT( <input> [, <allowInvalid> ] )

TO_GEOMETRY( <input> [, <allowInvalid> ] )

ST_GEOMFROMWKB( <input> [, <allowInvalid> ] )

ST_GEOMFROMWKT( <input> [, <allowInvalid> ] )

SELECT TO_GEOMETRY('LINESTRING(100 102,100 102)', TRUE);

SELECT TO_GEOMETRY('LINESTRING(100 102,100 102)', TRUE) AS g, ST_ISVALID(g);