9 feb 2011

Multiple Scatter Correction (MSC)

Leyendo un interesante artículo de T. Davies y T. Fearn sobre el tratamiento MSC en la revista Spectroscopy Europe, me puse a desarrollar uno de los gráficos que proponen para un mejor entendimiento de este tratamiento para la corrección del "scatter".
Para ello use los espectros del Shootout 2002, los exporté a Excel y calculé el espectro medio. Posteriormente en Excel representé el X-Y plot de la 1ª muestra del producto Calibra1, respecto al espectro medio de este mismo producto, obteniendo la siguiente gráfica:
Reading an interesting article by T. Davies & T. Fearn about MSC treatment in the journal Spectroscopy Europe, I began to develop one of the graphics proposed for a better understanding of this treatment for the correction of "scatter".
To do this use the spectra of the Shootout 2002, export to Excel and calculated the average spectrum. Later in Excel I represented the X-Y plot of the 1st sample of the product Calibra1, compared to the average spectrum of the same product, obtaining the following graph:

En esta gráfica, podemos calcular por mínimos cuadrados la pendiente y el intercepto que usaremos para corregir el scatter en este espectro en particular, pero haremos lo mismo con todos los espectros, respecto al espectro medio.
Para corregir el espectro, restaremos al valor de log1/T para cada data point el intercepto y dividiremos el resultado por la pendiente: (Xi - a)/b.
Debido a que este tratamiento se basa en el espectro medio del conjunto de datos original, si este se amplia este se debe de calcular nuevamente el espectro medio, así como la pendiente (b) e intercepto (a) de nuevo para cada espectro respecto al nuevo espectro medio, lo que representa un cierto inconveniente hasta que la base de datos sea lo suficientemente representativa.
In this graph, we can calculate least-squares slope and intercept that we will use to correct the scatter in this spectrum in particular, but will do the same with all spectra, compared to the average spectrum.
To correct a spectrum, subtract the value of log1 / T for each data point intercept and divide the result by the slope (Xi - a) / b.
Because this treatment is based on the average spectrum of the original data set, if we add more samples we should recalculate the average spectrum, the slope (b) and Intercept (a) again for each spectrum with respect to new media spectrum, representing some inconvenience until the database is sufficiently representative.

Conjunto Espectral sin ningun tratamiento:
(raw spectra)

Aplicando el tratamiento MSC:
(spectra corrected with MSC)



No hay comentarios:

Publicar un comentario