在 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()
我想进行一个查询,返回位于说英语的国家/地区的所有房屋。
Country 和 House 之间的关系应该通过交叉来完成House.location
with Country.territory
.
我如何使用 GeoDjango 的 ORM 来实现这一点?
一个有用且相当优化的解决方案是将英语国家的多边形组合成一个多多边形 https://wiki.openstreetmap.org/wiki/Relation:multipolygon(由至少 2 个明确定义的多边形生成的区域)。然后过滤哪些点与该区域相交。
为此,我们将使用 GeoDjangoUnion https://docs.djangoproject.com/en/2.1/ref/contrib/gis/geoquerysets/#django.contrib.gis.db.models.Union:
返回一个 GEOSGeometry 对象,其中包含查询集中每个几何图形的并集。请注意,Union 的使用是处理器密集型的,并且在大型查询集上可能会花费大量时间
里面一个Subquery https://docs.djangoproject.com/en/2.0/ref/models/expressions/#subquery-expressions:
Houses.objects.filter(
location__intersects=Subquery(
Country.objects.filter(language='English')
.aggregate(area=Union('territory'))['area']
)
)
或者我们可以避免子查询(对于 Django 版本
engish_speaking_area = Country.objects.filter(language='English')
.aggregate(area=Union('territory'))['area']
Houses.objects.filter(location__intersects=english_speaking_area)
另一种方法是在这里修改我的答案:GeoDjango 查询:包含在多多边形中的所有点 https://stackoverflow.com/questions/37166891/geodjango-query-all-point-that-are-contained-into-a-multi-polygon/47263632#47263632满足您的需求。
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)