在Geodjango中按几何交集关联两个模型

时间:2018-07-10 04:10:47

标签: python django postgis 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.locationCountry.territory相交来完成。

如何使用GeoDjango的ORM实现此目标?

1 个答案:

答案 0 :(得分:2)

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

为此,我们将使用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以适应您的需求。