- def compute_bounding_sphere(self, obj):
- p1 = max(((v.co, v.co.length) for v in obj.data.vertices), key=lambda x:x[1])[0]
- p2 = max(((v.co, (v.co-p1).length) for v in obj.data.vertices), key=lambda x:x[1])[0]
- center = (p1+p2)/2
- radius = (p1-p2).length/2
- for v in obj.data.vertices:
- d = v.co-center
- if d.length>radius:
- center += d*(1-radius/d.length)/2
- radius = (radius+d.length)/2
-
- return center, radius
-