def getNeighbourDists(self, forceRetriang = False):
from PYME.Analysis.LMVis import visHelpers
if forceRetriang or not 'neighbourDistances' in self.GeneratedMeasures.keys():
statNeigh = statusLog.StatusLogger("Calculating mean neighbour distances ...")
self.GeneratedMeasures['neighbourDistances'] = np.array(visHelpers.calcNeighbourDists(self.getTriangles(forceRetriang)))
return self.GeneratedMeasures['neighbourDistances']
def execute(self, namespace):
from matplotlib import delaunay
from PYME.Analysis.LMVis import visHelpers
pos = namespace[self.inputPositions]
x, y = pos['x'], pos['y']
#triangulate the data
T = delaunay.Triangulation(x + .1*np.random.normal(size=len(x)), y + .1*np.random.normal(size=len(x)))
#find the average edge lengths leading away from a given point
res = np.array(visHelpers.calcNeighbourDists(T))
namespace[self.outputName] = {self.key:res}
def getNeighbourDists(self, forceRetriang=False):
from PYME.Analysis.LMVis import visHelpers
if forceRetriang or not 'neighbourDistances' in self.GeneratedMeasures.keys(
):
statNeigh = statusLog.StatusLogger(
"Calculating mean neighbour distances ...")
self.GeneratedMeasures['neighbourDistances'] = np.array(
visHelpers.calcNeighbourDists(
self.getTriangles(forceRetriang)))
return self.GeneratedMeasures['neighbourDistances']
