def averagePeg(pegList, planet):
'''Computes the average of a given list of pegpoints.'''
nPeg = len(pegList)
elp = planet.get_elp()
avgPeg = Peg()
for attribute in ['latitude', 'longitude', 'heading']:
setattr(avgPeg, attribute, sum([getattr(pegPt, attribute) for pegPt in pegList])/(1.0*nPeg))
avgPeg.updateRadiusOfCurvature(elp)
return avgPeg
def medianPeg(pegList, planet):
'''Computes the median of a given list of pegpoints.'''
import numpy
elp = planet.get_elp()
medPeg = Peg()
nPeg = len(peglist)
for attribute in ['latitude', 'longitude', 'heading']:
setattr(medPeg, attribute,numpy.median([getattr(pegPt, attribute) for pegPt in pegList]))
medPeg.updateRadiusOfCurvature(elp)
return medPeg