【SpatiaLite】不正なジオメトリを正しくするMakeValid

WKTでPOLYGONを作ったときに、うっかり始点と終点を一致させず、経路を閉じ忘れることがあります。

SELECT
    IsValid(g) AS IsValid,
    IsValidReason(g) As IsValidReason
FROM (
    SELECT
        GeomFromText(
            'POLYGON((0 0, 0 1, 1 1, 1 0))',
            4326
        ) AS g
)
;

-- 0, Invalid: Unclosed Rings were detected

MakeValidはこういった不正なジオメトリを、形を崩さずに正しいジオメトリに直そうとしてくれます。

SELECT
    AsText(g2) AS geom_text,
    IsValid(g2) AS IsValid
FROM (
    SELECT
        MakeValid(g) AS g2
    FROM
    (
        SELECT
            GeomFromText(
                'POLYGON((0 0, 0 1, 1 1, 1 0))',
                4326
            ) AS g
    )
)
;

-- POLYGON((0 0, 0 1, 1 1, 1 0, 0 0)),  1

また、POLYGONで与えた頂点が描く軌跡は、交わってはいけません。交わってしまう状態を自己交差 (self-intersection) といいます。

SELECT
    IsValid(g) AS IsValid,
    IsValidReason(g) AS IsValidReason
FROM (
    SELECT
        GeomFromText(
            'POLYGON((0 0, 0 1, 2 1, 2 2, 1 2, 1 0, 0 0))',
            4326
        ) AS g
)
;

-- 0, Self-intersection[1 1]

この場合、座標 (1, 1) で自己交差がおきています。
このようなPOLYGONもMakeValidで正常化できます。

SELECT
    AsText(g2) AS geom_text,
    IsValid(g2) AS IsValid
FROM (
    SELECT
        MakeValid(g) AS g2
    FROM (
        SELECT
            GeomFromText(
                'POLYGON((0 0, 0 1, 2 1, 2 2, 1 2, 1 0, 0 0))',
                4326
            ) AS g
    )
)
;

-- MULTIPOLYGON(((0 0, 0 1, 1 1, 1 0, 0 0)), ((1 1, 1 2, 2 2, 2 1, 1 1))), 1

自己交差が起こっていた座標 (1, 1) で2つのPOLYGONに分けられ、MULTIPOLYGONになりました。