地理空間データ型¶

GEOGRAPHY データ型¶

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

地球上のポイントは、経度（-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」）

• Well-Known Binary （「WKB」）

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

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

次の参照資料には、関連項目が記載されています。

• Open Geospatial Consortium's Simple Feature Access / Common Architecture and SQL Option:

  • https://www.opengeospatial.org/standards/sfa

  • https://www.opengeospatial.org/standards/sfs

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

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

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

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

  • ポイント

  • MultiPoint

  • LineString

  • MultiLineString

  • ポリゴン

  • MultiPolygon

  • GeometryCollection

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

  • 機能

  • FeatureCollection

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

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

GEOGRAPHY_OUTPUT_FORMAT および GEOMETRY パラメーターは、次のいずれかの値を指定できます。

パラメーター値	説明
GeoJSON （デフォルト）	GEOGRAPHY / GEOMETRY の結果は、 GeoJSON 形式の OBJECT としてレンダリングされます。
WKT	GEOGRAPHY / GEOMETRY の結果は、 WKT 形式の VARCHAR としてレンダリングされます。
WKB	GEOGRAPHY / GEOMETRY の結果は、 WKB 形式の BINARY としてレンダリングされます。
EWKT	GEOGRAPHY / GEOMETRY の結果は、 EWKT 形式の VARCHAR としてレンダリングされます。
EWKB	GEOGRAPHY / GEOMETRY の結果は、 EWKB 形式の BINARY としてレンダリングされます。

EWKT および EWKB の場合、出力では SRID は常に4326です。 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 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)');

Copy

alter session set GEOGRAPHY_OUTPUT_FORMAT='GeoJSON';

Copy

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" |
| }                      |
+------------------------+

Copy

alter session set GEOGRAPHY_OUTPUT_FORMAT='WKT';

Copy

select g
    from geospatial_table
    order by id;
+-------------------------------------+
| G                                   |
|-------------------------------------|
| POINT(-122.35 37.55)                |
| LINESTRING(-124.2 42,-120.01 41.99) |
+-------------------------------------+

Copy

alter session set GEOGRAPHY_OUTPUT_FORMAT='WKB';

Copy

select g
    from geospatial_table
    order by id;
+------------------------------------------------------------------------------------+
| G                                                                                  |
|------------------------------------------------------------------------------------|
| 01010000006666666666965EC06666666666C64240                                         |
| 010200000002000000CDCCCCCCCC0C5FC00000000000004540713D0AD7A3005EC01F85EB51B8FE4440 |
+------------------------------------------------------------------------------------+

Copy

alter session set GEOGRAPHY_OUTPUT_FORMAT='EWKT';

Copy

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) |
+-----------------------------------------------+

Copy

alter session set GEOGRAPHY_OUTPUT_FORMAT='EWKB';

Copy

select g
    from geospatial_table
    order by id;
+--------------------------------------------------------------------------------------------+
| G                                                                                          |
|--------------------------------------------------------------------------------------------|
| 0101000020E61000006666666666965EC06666666666C64240                                         |
| 0102000020E610000002000000CDCCCCCCCC0C5FC00000000000004540713D0AD7A3005EC01F85EB51B8FE4440 |
+--------------------------------------------------------------------------------------------+

Copy

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);

Copy

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

SELECT ST_TRANSFORM(geometry_expression, 4326, 28992);

Copy

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

SELECT ST_SETSRID(geometry_expression, 4326);

Copy

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.8
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進数表現を取得します（逆も同様）。

H3の SQL 関数を以下に示します。

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_GEOM_POINT(13.4814, 52.5015),
         ST_GEOM_POINT(-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 の出力を比較しています。

ST_INTERSECTS . GEOGRAPHY 入力の使用	ST_INTERSECTS . GEOMETRY 入力の使用
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 | +--------------------------+ Copy	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 | +--------------------------+ Copy

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

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

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

地理空間関数 が異なれば、無効な形状に対して 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);