获取一定半径内的点

Question

我在 PostgreSQL 中有一个纬度和经度列。我想获得一定半径内的所有点。我知道他们是可以实现这一目标的库和扩展，但我想手动完成。 这是我从数据库中获取纬度和经度所做的：

lat1 = Places.query.with_entities(Places.latitude).all()
lon1 = Places.query.with_entities(Places.longitude).all()

接下来我编写了一个函数来将纬度和经度转换为以公里为单位的距离：

result = []
for i in range(len(lat1)):

    x = Decimal(111.12) * (Decimal(28.616700) - lat1[i]['latitude'])
    y = Decimal(111.12) * (Decimal(77.216700) - lon1[i]['longitude']) * Decimal(
    cos(lat1[i]['latitude'])) / Decimal(92.215)
    result.append(sqrt(x * x + y * y))
print result

如何在这个函数中包含我想要的点的半径。另外，lat1 和 lon1 的输出看起来像这样：

{u'latitude': Decimal('28.633300')}

我发现很难转换并不断出错：

TypeError: unsupported operand type(s) for -: 'float' and 'list'

我为 lat1lat1 和 lon1lon1 创建了一个 pastebin，以显示 lat1 和 lon1 的确切输出。

Answer 1

您的 lat1 和 lon1 是字典列表。因此，当您访问它们的实际值时，请在循环中执行以下操作，

result = []
for i in range(len(lat1)):
    x = Decimal(111.12) * (Decimal(28.616700) - lat1[i]['latitude']);
    y = Decimal(111.12) * (Decimal(77.216700) - lon1[i]['longitude']) * Decimal(math.cos(lat1[i]['latitude'])) / Decimal(92.215);
    result.append(math.sqrt(x * x + y * y));

以下是示例案例的控制台输出，

>>> lat1=[{u'latitude': Decimal('28.633300')}]
>>> lon1=[{u'longitude': Decimal('28.633300')}]
>>> x = Decimal(111.12) * (Decimal(28.616700) - lat1[0]['latitude']);
>>> y = Decimal(111.12) * (Decimal(77.216700) - lon1[0]['longitude']) * Decimal(math.cos(lat1[0]['latitude'])) / Decimal(92.215);
>>> result = math.sqrt(x * x + y * y);
>>> result
54.84298321157983