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
精彩评论