はじめに

PostGISってなに？、環境構築？から１歩進んで、実際にPostGISを試す流れを紹介します。

環境

以下の記事を参考にDocker環境を構築します。
dockerでPostGISを入れたPostgreSQL環境構築

多角形内にある点を探す

※以下は、先の参考記事のファイル名やデータベース構成を元に作成しています。

検索対象となる点を登録します

対象となる点は、緯度経度を示す点を 3 x 3 のマトリクスとします。

緯度経度	133	134	135
35	POINT(133, 35)	POINT(134, 35)	POINT(135, 35)
34	POINT(133, 34)	POINT(134, 34)	POINT(135, 34)
33	POINT(133, 33)	POINT(134, 33)	POINT(135, 33)

insert

INSERT INTO test_table (gid, geom) VALUES
(1, ST_GeomFromText('POINT(135 35)', 4612)),
(2, ST_GeomFromText('POINT(135 34)', 4612)),
(3, ST_GeomFromText('POINT(135 33)', 4612)),
(4, ST_GeomFromText('POINT(134 35)', 4612)),
(5, ST_GeomFromText('POINT(134 34)', 4612)),
(6, ST_GeomFromText('POINT(134 33)', 4612)),
(7, ST_GeomFromText('POINT(133 35)', 4612)),
(8, ST_GeomFromText('POINT(133 34)', 4612)),
(9, ST_GeomFromText('POINT(133 33)', 4612));

抽出条件となる多角形は？

先の 3 x 3 のマトリクスから中央のPOINT(134, 34)を抽出するとして、それを囲む多角形を定義します。

抽出イメージ

抽出対象の点

緯度経度	133	134	135
35
34		POINT(134, 34) ↑これを抽出したい
33

抽出条件となる多角形

「◇」みたいなひし形をイメージしています

緯度経度	133	134	135
35		POINT(134, 35)
34	POINT(133, 34)		POINT(135, 34)
33		POINT(134, 33)

検索クエリ

select

SELECT
  gid,
  ST_Astext(ST_Transform(geom, 4612))
FROM test_table
WHERE ST_Contains(
    ST_Transform(
      ST_GeomFromText('POLYGON((135 34, 134 35, 133 34, 134 33, 135 34))', 4612),
      4612
    ),
    geom
  );

※読み方

ST_Astext(ST_Transform(geom, 4612))
- geomフィールドを、空間参照系識別コード（SRID）を4612：JGD2000（地理座標系）に変換して、テキストに変換します。
ST_Contains(A, B)
- Aに含まれるBを探す
ST_GeomFromText('POLYGON((135 34, 134 35, 133 34, 134 33, 135 34))', 4612),
- 始点と終点が同じある輪で囲まれた多角形を示すテキストを、4612：JGD2000（地理座標系）に変換します。

検索結果

想定した値が返ってきましたね！

蛇足

ST_Astext(ST_Transform(geom, 4612))の個所を
ST_Transform(geom, 4612)とすると

何を言っているかさっぱり

PostGISで多角形内にある点を探す