mnorm2d = np.ma.masked_less(norm2d, threshold)
com = ndimage.measurements.center_of_mass(mnorm2d)
plt.imshow(norm2d, origin='lower')
plt.scatter(*com[::-1])
plt.show()
center_of_mass ndimage для расчета положения пика Гаусса