Django gis.db.models.functions.Intersection 関数徹底解説

Djangoのgis.db.models.functions.Intersection関数

gis.db.models.functions.Intersectionは、Djangoのcontrib.gisモジュールで提供される関数です。この関数は、2つのジオメトリオブジェクトの共通部分（交差部分）を計算します。

使い方

Intersection関数は、以下の引数を受け取ります。

• geom1: 最初のジオメトリオブジェクト。

Intersection関数は、以下の出力フィールドを生成します。

• intersection: 2つのジオメトリオブジェクトの共通部分。

例

以下のコードは、2つのポリゴンオブジェクトの交差部分を見つける例です。

from django.contrib.gis.db.models.functions import Intersection

# モデル
class MyModel(models.Model):
    geom = models.PolygonField()

# クエリ
polygons = MyModel.objects.filter(
    geom__intersects=Intersection(
        # 最初のジオメトリ
        models.Point(x=1, y=1),
        # 2番目のジオメトリ
        models.Polygon(
            ((0, 0), (10, 0), (10, 10), (0, 10), (0, 0)),
        ),
    )
)

このクエリは、最初のジオメトリオブジェクト（点）と2番目のジオメトリオブジェクト（ポリゴン）が交差するすべてのMyModelオブジェクトを取得します。

詳細

Intersection関数に関する詳細は、以下のDjangoドキュメントを参照してください。

その他の関数

gis.db.models.functionsモジュールには、Intersection関数以外にも、さまざまなジオメトリ関数があります。

これらの関数を使用して、ジオメトリオブジェクトに関するさまざまな情報を取得することができます。

gis.db.models.functions.Intersection関数は、2つのジオメトリオブジェクトの交差部分を見つけるために使用できます。この関数を使用して、空間データ分析を行うことができます。

gis.db.models.functions.Intersection 関数のサンプルコード

点とポリゴンの交差

from django.contrib.gis.db.models.functions import Intersection

# モデル
class MyModel(models.Model):
    geom = models.PolygonField()

# クエリ
polygons = MyModel.objects.filter(
    geom__intersects=Intersection(
        # 点
        models.Point(x=1, y=1),
        # ポリゴン
        models.Polygon(
            ((0, 0), (10, 0), (10, 10), (0, 10), (0, 0)),
        ),
    )
)

2つのポリゴンの交差

from django.contrib.gis.db.models.functions import Intersection

# モデル
class MyModel(models.Model):
    geom = models.PolygonField()

# クエリ
polygons = MyModel.objects.filter(
    geom__intersects=Intersection(
        # 最初のポリゴン
        models.Polygon(
            ((0, 0), (5, 0), (5, 5), (0, 5), (0, 0)),
        ),
        # 2番目のポリゴン
        models.Polygon(
            ((3, 3), (8, 3), (8, 8), (3, 8), (3, 3)),
        ),
    )
)

このクエリは、2つのポリゴンの交差部分を持つすべての MyModel オブジェクトを取得します。

複数のジオメトリとの交差

from django.contrib.gis.db.models.functions import Intersection

# モデル
class MyModel(models.Model):
    geom = models.PolygonField()

# クエリ
polygons = MyModel.objects.filter(
    geom__intersects=Intersection(
        # 点
        models.Point(x=1, y=1),
        # ポリゴン
        models.Polygon(
            ((0, 0), (10, 0), (10, 10), (0, 10), (0, 0)),
        ),
        # ライン
        models.LineString(((-1, -1), (11, 11))),
    )
)

このクエリは、点 (1, 1)、ポリゴン、ラインのいずれかと交差するすべての MyModel オブジェクトを取得します。

交差部分の面積

from django.contrib.gis.db.models.functions import Intersection, Area

# モデル
class MyModel(models.Model):
    geom = models.PolygonField()

# クエリ
polygons = MyModel.objects.annotate(
    intersection_area=Area(
        Intersection(
            # 最初のポリゴン
            models.Polygon(
                ((0, 0), (5, 0), (5, 5), (0, 5), (0, 0)),
            ),
            # 2番目のポリゴン
            models.Polygon(
                ((3, 3), (8, 3), (8, 8), (3, 8), (3, 3)),
            ),
        )
    )
)

# 結果
for polygon in polygons:
    print(polygon.intersection_area)

このコードは、2つのポリゴンの交差部分の面積を計算します。

複雑な条件

from django.contrib.gis.db.models.functions import Intersection, Within

# モデル
class MyModel(models.Model):
    geom = models.PolygonField()

# クエリ
polygons = MyModel.objects.filter(
    geom__intersects=Intersection(
        # 最初のポリゴン
        models.Polygon(
            ((0, 0), (10, 0), (10, 10), (0, 10), (0, 0)),
        ),
        # 2番目のポリゴン
        models.Polygon(
            ((3, 3), (8, 3), (8, 8), (3, 8), (3, 3)),
        ),
    ),
    geom__within=models.Point(x=

gis.db.models.functions.Intersection 関数の代替方法

geom__intersects ルックアップを使用して、2つのジオメトリオブジェクトが交差しているかどうかを調べることができます。

# モデル
class MyModel(models.Model):
    geom = models.PolygonField()

# クエリ
polygons = MyModel.objects.filter(
    geom__intersects=models.Polygon(
        ((0, 0), (10, 0), (10, 10), (0, 10), (0, 0)),
    )
)

このクエリは、渡されたポリゴンと交差するすべての MyModel オブジェクトを取得します。

GeoQuerySet メソッドを使用して、2つのジオメトリオブジェクトの交差部分を見つけることができます。

# モデル
class MyModel(models.Model):
    geom = models.PolygonField()

# クエリ
polygons = MyModel.objects.filter(
    geom__intersects=models.Polygon(
        ((0, 0), (10, 0), (10, 10), (0, 10), (0, 0)),
    )
).intersection(
    models.Polygon(
        ((3, 3), (8, 3), (8, 8), (3, 8), (3, 3)),
    )
)

このクエリは、最初のポリゴンと2番目のポリゴンの両方に含まれるすべての MyModel オブジェクトを取得します。

カスタム SQL を使用して、2つのジオメトリオブジェクトの交差部分を見つけることもできます。

SELECT *
FROM my_model
WHERE geom && ST_Intersects(
    ST_GeomFromText('POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))'),
    geom
);

このクエリは、最初のポリゴンと交差するすべての my_model テーブルのレコードを取得します。

gis.db.models.functions.Intersection 関数は、2つのジオメトリオブジェクトの交差部分を見つけるための強力なツールです。ただし、他の方法もいくつかあるので、状況に応じて最適な方法を選択する必要があります。