SQLServerのgeography型をポリゴンで検索する場合は左回りである必要があります

2022年3月19日

はじめに

SQLServerのgeography型をポリゴンで空間検索する際、座標が左回りでなければうまく検索できません
NetTopologySuiteを使用し、座標を左回り(時計回り)に修正する方法です

SQL

geography型のlocationを持つ、customersテーブルを空間検索するとします
下記例は左回りの座標でポリゴンを作成しているので正常に検索できます

SELECT 
    * 
FROM 
    customers 
WHERE 
    location.STIntersects(geography::STGeomFromText('POLYGON((139.781225 35.686380, 139.782368 35.683731, 139.777267 35.683319, 139.781225 35.686380))', 4326)) = 1

下記は右回りの座標でポリゴンを作成しています
エラーにはなりませんが、空間検索ができません

SELECT 
    *
FROM 
    customers 
WHERE 
    location.STIntersects(geography::STGeomFromText('POLYGON((139.781225 35.686380, 139.777267 35.683319, 139.782368 35.683731, 139.781225 35.686380))', 4326)) = 1

NetTopologySuite.Geometries.Polygonを使う

nugetでNetTopologySuiteをインストールします

あわせて読みたい

NetTopologySuite 2.5.0 The NTS Topology Suite is an API for modelling and manipulating 2-dimensional linear geometry. It provides numerous geometric predicates and functions. NTS conf...

NetTopologySuite.Geometries.Polygonを使えば、右回り座標を左回りに修正することができます
polygon.Shell.IsCCWがfalseなら右回りなので、座標を逆転させると左回りになります

// 右回りの座標 
var coordinates = new[] { 
        new Coordinate(139.781225, 35.686380), 
        new Coordinate(139.777267, 35.683319), 
        new Coordinate(139.782368, 235.683731), 
        new Coordinate(139.781225, 35.686380) 
}; 
var polygon = new Polygon(new LinearRing(coordinates)); 
if(!polygon.Shell.IsCCW){ 
    // 右回りなら座標を逆転させる 
    coordinates = coordinates.Reverse().ToArray(); 
}