maxloc

pylayers.gis.ezone.maxloc(f, threshold=-1.4142135623730951)

determine local maximum above a threshold

f : np.array threshold : float

gnp.array

values of local maximum

indnp.array

index of local maximum

>>> import numpy as np
>>> t = np.arange(0,6*np.pi)
>>> f = np.sin(2*t)*cos(3*t)*sin(5*t)
>>> g,ind = maxloc(f,threshold=0.3)