hist, bins = np.histogram(norm2d.ravel(), normed=True, bins=100)
threshold = bins[np.cumsum(hist) * (bins[1] - bins[0]) > 0.8][0]
mnorm2d = np.ma.masked_less(norm2d, threshold)
com = ndimage.measurements.center_of_mass(mnorm2d)
center_of_mass ndimage для расчета положения пика Гаусса