Die Rolle der Vorverarbeitung in der Chemometrie
Insbesondere Nahinfrarotspektren (NIR).chemische und pharmazeutische Spektren sind hochkomplex und oft durch überlappende Absorptionsbanden, Grundlinienversätze und hochfrequentes Rauschen gekennzeichnet. Vor der Anpassung quantitativer Kalibrierungsmodelle (wie der PLS-Regression) müssen Chemometriker die rohen Spektraldaten vorverarbeiten. Das Ziel besteht darin, physikalische Artefakte – wie Streuung durch unterschiedliche Partikelgrößen oder Basislinienverschiebungen aufgrund von Änderungen der Lichtweglänge – zu beseitigen und gleichzeitig die interessierenden chemischen Signale zu erhalten.
Unter den Werkzeugen im Vorverarbeitungsarsenal ist der Savitzky-Golay (SG)-Algorithmus der Industriestandard für Glättung und Differenzierung. Es wurde 1964 von Abraham Savitzky und Marcel J. E. Golay entwickelt und ist bis heute die Grundlage moderner chemometrischer Datenpipelines.
„Im Gegensatz zu Standardfiltern mit gleitendem Durchschnitt, die spektrale Merkmale verwischen, bewahrt die Savitzky-Golay-Filterung Peakhöhen, -breiten und lokale Formen, indem lokale Polynome niedrigen Grades an die Daten angepasst werden.“
So funktioniert der Savitzky-Golay-Algorithmus
Das Kernkonzept besteht darin, ein Fenster mit einer festen Anzahl von Datenpunkten ($w$, das eine ungerade ganze Zahl sein muss) über das Spektrum zu verschieben. Bei jedem Schritt wird ein Polynom vom Grad $k$ (normalerweise $k = 2$ oder $3$, was quadratische oder kubische Modelle darstellt) mithilfe der linearen Methode der kleinsten Quadrate an die Datenpunkte innerhalb des Fensters angepasst. Der Wert des angepassten Polynoms am Mittelpunkt des Fensters wird zum neuen, geglätteten Wert für diese Wellenlänge.
Da die Anpassungskoeffizienten der kleinsten Quadrate nur von der Fenstergröße $w$ und dem Polynomgrad $k$ abhängen, können sie als feste Gewichte (sogenannte Faltungen) vorberechnet werden. Dadurch lässt sich der SG-Filter auch für Tausende langer Spektren extrem schnell berechnen.
Derivate und Basislinienkorrektur
Eine der leistungsstärksten Funktionen der Savitzky-Golay-Methode ist ihre Fähigkeit, während des Glättungsprozesses numerische Ableitungen zu berechnen. Derivate sind äußerst effektiv bei der Korrektur von Basislinienversätzen und der Auflösung überlappender Peaks:
- Erste Ableitung (SG 1.):Entfernt konstante Grundlinienversätze (Verschiebungen entlang der Y-Achse) und hebt Neigungsänderungen hervor.
- Zweite Ableitung (SG 2nd):Entfernt die lineare Basisliniendrift (Steigungen entlang des Spektrums) und isoliert überlappende Features. In einem Spektrum der zweiten Ableitung erscheinen Absorptionspeaks als nach unten gerichtete Täler, wodurch sie leichter zu identifizieren sind.
Richtlinien zur Parameterauswahl
Das Einstellen der richtigen Parameter ist entscheidend, um eine Unterglättung (Beibehaltung des Rauschens) oder eine Überglättung (Zerstörung der tatsächlichen Spitzenauflösung) zu verhindern. Die folgende Tabelle enthält Standardempfehlungen basierend auf der Spektrometerauflösung:
| Parametersatz | Fenstergröße ($w$) | Polynomordnung ($k$) | Empfohlene Anwendung |
|---|---|---|---|
| Milde Glättung | 5 oder 7 Punkte | 2 (Quadratisch) | Hochauflösende Gasphasen-NIR, rauscharme Systeme |
| Standard-Grundlinienkorrektur | 15 Punkte | 2 (Quadratisch, 1. Ableitung) | Typische flüssige chemische Prozesse, Feuchtigkeit in Pulvern |
| Tiefe Trennung | 25 Punkte | 3 (kubisch, 2. Ableitung) | Stark heterogene Feststoffe (Futtermehle, Futterpflanzen) |
Implementierung von SG im caliXSuite
Im caliX1,50 % – 6,00 %In der Vorverarbeitungspipeline wird der SG-Algorithmus automatisch während der AutoML-Modellvalidierung angewendet. Benutzer können die Fenstergröße und die Ableitungsreihenfolge im Vorverarbeitungseditor konfigurieren und beobachten, wie der Kreuzvalidierungsfehler (RMSECV) in Echtzeit aktualisiert wird, während verschiedene Einstellungen getestet werden. Diese interaktive Rückkopplungsschleife garantiert, dass Sie die mathematischen Vorverarbeitungsparameter auswählen, die die Modellgenauigkeit und die Übertragbarkeit von Einheit zu Einheit maximieren.
Referenzen
- Savitzky, A. & Golay, M. J. E. (1964). „Glättung und Differenzierung von Daten durch vereinfachte Verfahren der kleinsten Quadrate“,Analytische Chemie, 36(8), 1627-1639.
- „Wellenlängenauswahl und Rauschfilterung in NIR-Spektren unter Verwendung von Polynomableitungen“,Chemometrie und intelligente Laborsysteme, 2019.
- „Grundlinienkorrekturprotokolle für industrielle Prozessspektrometer“,Rezensionen zur angewandten Spektroskopie, 2021.
