GeoDjango：确定多边形的面积

Question

在我的模型中，我有一个定义的多边形字段

polygon = models.PolygonField(srid=4326, geography=True, null=True, blank=True)

当我想确定多边形的区域时，我打电话

area_square_degrees = object.polygon.area

但是如何使用GeoDjango将square degrees中的结果转换为m2？ This answer不起作用，因为area没有方法sq_m。有内置转换吗？

Answer 1

您需要将数据转换为正确的空间参照系。

area_square_local_units = object.polygon.transform(srid, clone=False).area

在英国，您可以使用英国国家电网SRID 27700，它使用电表。

area_square_meters = object.polygon.transform(27700, clone=False).area

您可能想要也可能不想克隆几何体，具体取决于您是否需要在未转换状态下对其进行任何其他操作。

文件在这里https://docs.djangoproject.com/en/1.8/ref/contrib/gis/geos/

Answer 2

我一直在努力解决这个问题，因为我找不到一个干净的解决方案。诀窍是你必须使用postgis功能（因此它只使用postgis ..）：

from django.contrib.gis.db.models.functions import Area
loc_obj = Location.objects.annotate(area_=Area("poly")).get(pk=??)  
# put the primary key of the object

print(loc_obj.area_)  # distance object, result should be in meters, but you can change to the unit you want, e.g  .mi for miles etc..

models.py：

class Location(models.Model):
    poly = gis_models.PolygonField(srid=4326, geography=True)

如果您必须处理地理坐标而不是投影，我认为这是最好的方法。它确实处理地球的曲线计算，即使距离/面积很大，结果也很精确

Answer 3

我需要一个应用程序来获取全球多边形区域，如果我使用了无效的国家/地区投影我得到了错误qazxsw poi

我最后使用OGRException: OGR failure使用4326投影（默认投影）来避免涉及每个国家/地区的特定投影。这是我的代码：

OpenLayers implementation

只需将几何图形传递给

import math

def getCoordsM2(coordinates):
    d2r = 0.017453292519943295  # Degrees to radiant
    area = 0.0
    for coord in range(0, len(coordinates)):
        point_1 = coordinates[coord]
        point_2 = coordinates[(coord + 1) % len(coordinates)]
        area += ((point_2[0] - point_1[0]) * d2r) *\
            (2 + math.sin(point_1[1] * d2r) + math.sin(point_2[1] * d2r))
    area = area * 6378137.0 * 6378137.0 / 2.0
    return math.fabs(area)

def getGeometryM2(geometry):
    area = 0.0
    if geometry.num_coords > 2:
        # Outer ring
        area += getCoordsM2(geometry.coords[0])
        # Inner rings
        for counter, coordinates in enumerate(geometry.coords):
            if counter > 0:
                area -= getCoordsM2(coordinates)
    return area

函数即可完成！我在我的GeoDjango模型中使用此函数作为属性。希望能帮助到你！

Answer 4

如果它所说的地球表面积，1平方度有12,365.1613平方公里。所以多个你的平方度并乘以10 ^ 6转换为米。