在Geodjango中通过几何交点关联两个模型

问题描述

在GeoDjango中有两个包含几何场的模型：

from django.contrib.gis.db import models 

class Country(models.Model):
    territory = models.MultiPolygonField()
    language = models.CharField(max_length=2)

class House(models.Model):
    location = models.PointField()

我要进行一个查询，返回位于英语国家/地区的所有房屋。
乡村和住宅之间的关系应该通过House.location和Country.territory相交来完成。

如何使用GeoDjango的ORM实现此目的？

---

解决方案

一个有用且相当优化的解决方案是将讲英语的国家的多边形合并为一个multipolygon(至少由两个定义明确的多边形组成的区域)。然后过滤哪些点与该区域相交。

为此，我们将使用GeoDjango的Union：

返回一个包含查询集中每个几何体的并集的GEOSGeometry对象。请注意，UNION的使用是处理器密集型的，在大型查询集上可能会花费大量时间

在Subquery中：

Houses.objects.filter(
    location__intersects=Subquery(
        Country.objects.filter(language='English')
                       .aggregate(area=Union('territory'))['area']
    )
)

或者我们可以避免子查询(对于Django版本<；1.11)：

engish_speaking_area = Country.objects.filter(language='English')
                                      .aggregate(area=Union('territory'))['area']
Houses.objects.filter(location__intersects=english_speaking_area)

---

另一种方法是在这里修改我的答案：GeoDjango query: all point that are contained into a multi polygon以满足您的需要。