Essas bandas de espectro costumavam ser julgadas a olho nu, como fazer programaticamente?

Operadores usados ​​para examinar o espectro, conhecendo a localização elargura de cada pico e julgar a peça a que pertence o espectro. No novo modo, a imagem é capturada por uma câmera em uma tela. E a largura de cada banda deve ser calculada de forma programática.

Sistema antigo: espectroscópio -> olho humano Novo sistema: espectroscópio -> câmera -> programa

O que é um bom método paracalcular a largura de cada banda, dadas as suas posições aproximadas no eixo X; dado que esta tarefa costumava ser executada perfeitamente a olho, e agora deve ser realizada por programa?

Desculpe se estou com falta de detalhes, mas eles são escassos.

Listagem de programas que gerou o gráfico anterior; Espero que seja relevante:

import Image
from scipy import *
from scipy.optimize import leastsq

# Load the picture with PIL, process if needed
pic         = asarray(Image.open("spectrum.jpg"))

# Average the pixel values along vertical axis
pic_avg     = pic.mean(axis=2)
projection  = pic_avg.sum(axis=0)

# Set the min value to zero for a nice fit
projection /= projection.mean()
projection -= projection.min()

#print projection

# Fit function, two gaussians, adjust as needed
def fitfunc(p,x):
    return p[0]*exp(-(x-p[1])**2/(2.0*p[2]**2)) + \
        p[3]*exp(-(x-p[4])**2/(2.0*p[5]**2))
errfunc = lambda p, x, y: fitfunc(p,x)-y

# Use scipy to fit, p0 is inital guess
p0 = array([0,20,1,0,75,10])
X  = xrange(len(projection))
p1, success = leastsq(errfunc, p0, args=(X,projection))
Y = fitfunc(p1,X)

# Output the result
print "Mean values at: ", p1[1], p1[4]

# Plot the result
from pylab import *
#subplot(211)
#imshow(pic)
#subplot(223)
#plot(projection)
#subplot(224)
#plot(X,Y,'r',lw=5)
#show()

subplot(311)
imshow(pic)
subplot(312)
plot(projection)
subplot(313)
plot(X,Y,'r',lw=5)
show()

questionAnswers(3)

yourAnswerToTheQuestion