开发者

Furthest distance between any two points with GeoDjango

开发者 https://www.devze.com 2022-12-08 15:47 出处:网络
I have a collection of Point objects. I want to find the furthest distance between any two points. Imagine a circle enclosing all of these points. I want to know the diameter of that circle. How can I

I have a collection of Point objects. I want to find the furthest distance between any two points. Imagine a circle enclosing all of these points. I want to know the diameter of that circle. How can I do this in GeoDjango开发者_JAVA技巧?

edit: this is what I have so far:

>>> r=Route.objects.get(pk=1)
>>> a=Airport.objects.filter(routebase__route=r)

>>> # this route represents a few hundred miles flight into mexico and back
>>> a
[<Airport: MMDO>, <Airport: KELP>, <Airport: KELP>, <Airport: MMCU>]

>>> # a multipoint object with all the airports
>>> mpoint = a.collect()

>>> # a polygon that represents a ring around the multipoint
>>> p = mpoint.envelope

>>> #now I just need to get the diameter of this envelope
>>> p.length
19.065994262694986

???

what units is that in? Is this even the value I'm after?

edit2: OK I'm going to try this another way:

>>> r=Route.objects.get(pk=1)
>>> a=Airport.objects.filter(routebase__route=r)

>>> # this route represents a few hundred miles flight into mexico
>>> a
[<Airport: MMDO>, <Airport: KELP>, <Airport: MMCU>]

>>> # a multipoint object with all the airports
>>> mpoint = a.collect()

>>> # get the center point of the route polygon and get the
>>> # distance between each point and the centroid
>>> # the largest should be the diameter of the ring, right?
>>> cen = mpoint.centroid

>>> dist = []
>>> for p in mp:
        dist.append(LineString(p, cen).length)

>>> dis
[0.54555421739245946,
 0.61638306853425906,
 0.53442640535933494,
 0.54555421739245946]

>>> max(dist)
0.61638306853425906

?? again, what are these units?


OK, I figured it out.

def overall_distance(route):
    a = Airport.objects.filter(routebase__route=route).distinct()
    mp = a.collect()
    ct = mp.envelope.centroid
    na = a.distance(ct)

    dist = []
    for p in na:
        dist.append(p.distance)

    diameter = max(dist) * 2

    return diameter.nm
0

精彩评论

暂无评论...
验证码 换一张
取 消