本發(fā)明涉及一種對(duì)三維色譜數(shù)據(jù)進(jìn)行處理的色譜數(shù)據(jù)處理裝置,所述三維色譜數(shù)據(jù)是通過對(duì)包含經(jīng)液相色譜儀等色譜儀分離的成分的試樣或經(jīng)流動(dòng)注射(flowinjection)法導(dǎo)入的試樣反復(fù)進(jìn)行分光分析或質(zhì)量分析而收集,更詳細(xì)而言,涉及一種用于判定在色譜中出現(xiàn)的來自于目標(biāo)成分的波峰中有無雜質(zhì)等其它成分重合的數(shù)據(jù)處理裝置。
背景技術(shù):在利用光電二極管陣列(photodiodearray,PDA)檢測(cè)器等多通道型檢測(cè)器作為檢測(cè)器的液相色譜儀中,通過以對(duì)流動(dòng)相注入試樣的時(shí)點(diǎn)為基點(diǎn),對(duì)從管柱(column)出口溶出的試樣液反復(fù)獲取吸光光譜,可獲得具有時(shí)間、波長及吸光度這三個(gè)維度的三維色譜數(shù)據(jù)。并且,在利用質(zhì)譜儀作為檢測(cè)器的液相色譜儀(liquidchromatograph,LC)或氣相色譜儀(gaschromatograph,GC)中,通過在質(zhì)譜儀中反復(fù)進(jìn)行規(guī)定的質(zhì)量電荷比范圍的掃描測(cè)定,可獲得具有時(shí)間、質(zhì)量電荷比及信號(hào)強(qiáng)度這三個(gè)維度的三維色譜數(shù)據(jù)。在以下的說明中,舉出利用PDA檢測(cè)器的液相色譜儀為例,但是在將質(zhì)譜儀作為檢測(cè)器的色譜儀中情況也是同樣。圖8的(a)是通過所述液相色譜儀而獲得的三維色譜數(shù)據(jù)的示意圖。通過從所述三維色譜數(shù)據(jù)中提取特定的波長(例如λ0)下的吸光度數(shù)據(jù),可制作表示所述特定的波長λ0下的測(cè)定時(shí)刻(即保持時(shí)間)與吸光度的關(guān)系的如圖8的(b)所示的波長色譜。此外,通過從三維色譜數(shù)據(jù)中提取表示在特定的時(shí)點(diǎn)(測(cè)定時(shí)刻)的吸光度的數(shù)據(jù),也可以制作表示在所述時(shí)點(diǎn)的波長與吸光度的關(guān)系的吸光光譜。在這樣的液相色譜儀中,進(jìn)行已知的目標(biāo)成分的定量分析時(shí),通常,制作與所述目標(biāo)成分相對(duì)應(yīng)的吸收波長下的波長色譜。然后,通常是找到在所述波長色譜中出現(xiàn)的來自于目標(biāo)成分的波峰的起點(diǎn)Ts及終點(diǎn)Te,計(jì)算所述波峰的面積值,并對(duì)照預(yù)先求出所述面積值的標(biāo)準(zhǔn)曲線而算出定量值。以如上所述方式對(duì)目標(biāo)成分進(jìn)行定量時(shí),如果所提取的波長色譜中出現(xiàn)的波峰是只來自于目標(biāo)成分的波峰則沒有問題,但是波峰不一定僅限于來自單一成分(目標(biāo)成分),而經(jīng)常有重合著來自于分析者不想要的雜質(zhì)(廣而言之為其它成分)的信號(hào)的情況。如果分析者未留意到這種情況而進(jìn)行定量計(jì)算,則會(huì)損害定量精度。因此,在定量計(jì)算之前,會(huì)再三進(jìn)行雜質(zhì)判定處理(或波峰純度判定處理),調(diào)查在色譜中出現(xiàn)的某個(gè)波峰是否只來自于目標(biāo)成分,或者是否含有雜質(zhì)。作為針對(duì)色譜上的波峰的雜質(zhì)判定處理,迄今為止已提出有各種方法,并且已投入實(shí)際應(yīng)用。但是,實(shí)際狀況是現(xiàn)有的方法皆有短長,而沒有建設(shè)性的的方法。例如專利文獻(xiàn)1所述的雜質(zhì)判定方法中,是通過對(duì)在各測(cè)定時(shí)刻分別獲得的吸光光譜在目標(biāo)成分的極大(或極小)吸收波長下沿波長方向進(jìn)行微分而求出波長微分系數(shù),從而制作表示所述波長微分系數(shù)的隨時(shí)間變化的微分色譜。然后,判定是否在所述微分色譜上觀測(cè)到如通常的色譜中出現(xiàn)般的波峰波形,由此判斷波長色譜上的來自于目標(biāo)成分的波峰是否含有雜質(zhì)。所述方法雖然是可通過比較簡(jiǎn)單的運(yùn)算處理來以高準(zhǔn)確度判定雜質(zhì)的有無的優(yōu)異方法,但是如以下所述,原理上存在無法檢測(cè)出雜質(zhì)的情況。圖9的(a)~(c)是表示來自于目標(biāo)成分的吸光光譜(實(shí)線)與來自于雜質(zhì)的吸光光譜(虛線)的關(guān)系的一例的圖。在所述現(xiàn)有的雜質(zhì)判定方法中,如圖9的(a)所示,是將來自于目標(biāo)成分的吸光光譜的極值的波長λ0、即波長微分系數(shù)為0的波長下的雜質(zhì)的吸光光譜曲線的波長微分系數(shù)用于雜質(zhì)判定。因此,如圖9的(a)所示,當(dāng)在雜質(zhì)的吸光光譜中吸光度表示極大值的波長與波長λ0不一致,而所述光譜曲線在波長λ0下具有傾斜度時(shí),可檢測(cè)出雜質(zhì)。但是,如圖9的(b)所示,當(dāng)來自于目標(biāo)成分的吸光光譜的極值與來自于雜質(zhì)的吸光光譜的極值出現(xiàn)在同一波長中時(shí),雜質(zhì)的吸光光譜曲線的波長微分系數(shù)大致為零,因此無法檢測(cè)出雜質(zhì)。并且,如圖9的(c)所示,當(dāng)在來自于目標(biāo)成分的吸光光譜的極值附近來自于雜質(zhì)的吸光光譜的曲線的傾斜度小(在極端的情況下為平坦)時(shí),即使制作微分色譜,來自于雜質(zhì)的波峰也過小,從而被埋沒在噪聲中,結(jié)果有可能無法檢測(cè)出雜質(zhì)。再者,利用PDA檢測(cè)器等對(duì)通過未使用管柱的流動(dòng)注射分析(FlowInjectionAnalysis,F(xiàn)IA)法而導(dǎo)入的試樣進(jìn)行檢測(cè)時(shí)所獲得的數(shù)據(jù)也成為具有時(shí)間、波長及吸光度這三個(gè)維度的三維數(shù)據(jù),從而與通過液相色譜儀而收集的三維色譜數(shù)據(jù)實(shí)質(zhì)上相同。因此,本說明書中的“三維色譜數(shù)據(jù)”是設(shè)為也包括通過FIA法而收集的三維數(shù)據(jù)?,F(xiàn)有技術(shù)文獻(xiàn)專利文獻(xiàn)專利文獻(xiàn)1:國際公開第2013/035639號(hào)
技術(shù)實(shí)現(xiàn)要素:[發(fā)明所要解決的問題]本發(fā)明是為了解決所述問題而完成,其目的在于提供一種色譜數(shù)據(jù)處理裝置,即使在所述現(xiàn)有的雜質(zhì)判定方法中,難以對(duì)色譜上的目標(biāo)的波峰適當(dāng)?shù)嘏袛嘤袩o雜質(zhì)重合的情況下,也可以確切而穩(wěn)定地判定有無雜質(zhì)重合。[解決問題的技術(shù)手段]為了解決所述問題而完成的本發(fā)明是一種色譜數(shù)據(jù)處理裝置,對(duì)針對(duì)作為測(cè)定對(duì)象的試樣而收集的包含時(shí)間、信號(hào)強(qiáng)度及其以外的第三維度的三維色譜數(shù)據(jù)進(jìn)行處理,所述色譜數(shù)據(jù)處理裝置的特征在于包括:a)濾波器制作部,算出與多維的主向量正交的一個(gè)副向量,將所述副向量確定為雜質(zhì)提取用的濾波器,所述多維的主向量是利用向量來表達(dá)如下光譜,即表示關(guān)于作為觀測(cè)對(duì)象的目標(biāo)成分的第三維度與信號(hào)強(qiáng)度的關(guān)系的光譜或可當(dāng)作所述關(guān)系的光譜;以及b)含雜質(zhì)信息獲取部,計(jì)算處理對(duì)象多維向量與副向量的內(nèi)積,基于其計(jì)算結(jié)果判斷在所述處理對(duì)象的光譜中有無所述目標(biāo)成分以外的雜質(zhì),所述處理對(duì)象多維向量是利用向量來表達(dá)如下光譜,即根據(jù)針對(duì)作為測(cè)定對(duì)象的試樣而獲得的三維色譜數(shù)據(jù)所求出的光譜或來自于所述數(shù)據(jù)的處理對(duì)象的光譜,所述副向量被確定為所述濾波器。在這里所謂“第三維度”,是指例如波長或質(zhì)量電荷比等,所謂“三維色譜數(shù)據(jù)”,是指例如通過對(duì)包含利用色譜儀(LC或GC)的管柱而在時(shí)間方向上加以分離的各種成分的試樣,利用多通道型檢測(cè)器等檢測(cè)器反復(fù)獲取吸光光譜而獲得的數(shù)據(jù),或者通過利用質(zhì)譜儀反復(fù)獲取質(zhì)譜(massspectrum)而獲得的數(shù)據(jù)。并且,在這里所謂“三維色譜數(shù)據(jù)”,也可以是不對(duì)經(jīng)過色譜儀的管柱的試樣,而對(duì)通過FIA法,不使成分分離而導(dǎo)入的試樣利用多通道型檢測(cè)器或質(zhì)譜儀而獲得的數(shù)據(jù)。在本發(fā)明的色譜數(shù)據(jù)處理裝置中,將表示第三維度與信號(hào)強(qiáng)度的關(guān)系的光譜,例如將吸光光譜或質(zhì)譜作為利用向量表達(dá)的多維向量來處理。例如,如果是吸光光譜,則由于離散的各波長的吸光度的集合為吸光光譜,因此可將吸光光譜表示為(a(λ1)、a(λ2)、a(λ3)、……、a(λn)),對(duì)以a(λm)為成分的多維向量進(jìn)行定義。在這里,a(λm)為波長λm(m=1~n)下的吸光度。現(xiàn)在,將利用向量來表達(dá)在某個(gè)測(cè)定時(shí)點(diǎn)的處理對(duì)象的光譜的處理對(duì)象多維向量設(shè)為I,將利用向量來表達(dá)目標(biāo)成分的光譜的多維向量設(shè)為A,進(jìn)而用B表示利用向量來表達(dá)雜質(zhì)的光譜的多維向量。于是,處理對(duì)象多維向量I可利用以下(1)式的向量運(yùn)算來表示。I=A+B……(1)將表示雜質(zhì)的光譜的向量B分解成與表示目標(biāo)成分的光譜的向量A平行的方向上的向量Ba、以及與向量A正交的方向上的向量Bo來考慮。并且,進(jìn)而考慮與向量A正交的多維向量F。由于彼此正交的向量的內(nèi)積為0,所以向量F與向量Ba的內(nèi)積為0。因此,處理對(duì)象多維向量I與向量F的內(nèi)積等于向量Bo與向量F的內(nèi)積。即,以下的(2)式成立。I·F=Bo·F……(2)由于表示雜質(zhì)的光譜的向量B的長度與向量Bo的長度成比例,因此所述(2)式的右邊的Bo·F與向量B的長度成比例。因此,(2)式的左邊的向量?jī)?nèi)積I·F與表示雜質(zhì)的光譜的向量B的長度成比例。由此,向量?jī)?nèi)積I·F可用作表示雜質(zhì)的量的指標(biāo)值u。因此,在本發(fā)明的色譜數(shù)據(jù)處理裝置中,濾波器制作部算出與表示目標(biāo)成分的光譜的主向量A正交的副向量F,并將其確定為雜質(zhì)提取用的濾波器。然后,含雜質(zhì)信息獲取部計(jì)算表示由三維色譜數(shù)據(jù)求出的光譜或來自于所述數(shù)據(jù)的處理對(duì)象的光譜的向量I與作為濾波器的向量F的內(nèi)積,并基于其計(jì)算結(jié)果判斷是否存在雜質(zhì)。作為典型的一個(gè)方式,所述含雜質(zhì)信息獲取部針對(duì)伴隨著時(shí)間經(jīng)過的在各測(cè)定時(shí)點(diǎn)獲得的處理對(duì)象的光譜,分別計(jì)算表示所述光譜的向量I與作為濾波器的向量F的內(nèi)積。然后,可設(shè)為如下構(gòu)成:對(duì)沿時(shí)間序列的所述內(nèi)積的值的變化進(jìn)行觀測(cè),如果例如出現(xiàn)有類似色譜波峰的波形,則判斷為存在目標(biāo)成分以外的雜質(zhì)。在本發(fā)明的色譜數(shù)據(jù)處理裝置中,所述濾波器制作部是算出與多維的主向量正交的副向量作為雜質(zhì)提取用的濾波器,但是在多維向量空間內(nèi)存在多個(gè)與某個(gè)向量正交的向量。因此,優(yōu)選的是所述濾波器制作部可以來自于雜質(zhì)的光譜的向量Bo與作為所述濾波器的副向量F的余弦相似度成為最大、即盡可能接近于“1”的方式來確定所述副向量F的方向。由此,用所述(2)式表示的雜質(zhì)量的指標(biāo)值u的信噪(signaltonoise,SN)比成為最大或與其接近的狀態(tài),從而目標(biāo)成分以外的成分的有無的判斷的準(zhǔn)確性提高。再者,想要算出所述余弦相似度,需要算出向量Bo,這可如下所述以解析方式而求出。Bo=I-αA……(3)α=(I·A)/(A·A)如上所述,當(dāng)針對(duì)在各測(cè)定時(shí)點(diǎn)獲得的處理對(duì)象的光譜分別求出向量I與向量F的內(nèi)積時(shí),也可使用在各測(cè)定時(shí)點(diǎn)分別獲得的向量F,還可使用一個(gè)或多個(gè)代表性的向量F。例如作為一實(shí)施方式,可設(shè)為如下構(gòu)成:所述濾波器制作部求出在各測(cè)定時(shí)點(diǎn)求出的雜質(zhì)提取用的濾波器即多個(gè)向量的平均向量,所述含雜質(zhì)信息獲取部利用所述平均向量,計(jì)算針對(duì)表示在各測(cè)定時(shí)點(diǎn)的處理對(duì)象的光譜的各向量的內(nèi)積。由此,可針對(duì)噪聲將魯棒(robust)的向量用作雜質(zhì)提取用的濾波器,因此對(duì)有噪聲重疊的數(shù)據(jù)也可以確切地判定有無雜質(zhì)。并且,作為另一實(shí)施方式,也可設(shè)為如下構(gòu)成:所述濾波器制作部從在各測(cè)定時(shí)點(diǎn)求出的雜質(zhì)提取用的濾波器即多個(gè)向量之中選擇范數(shù)(norm)為最大的向量,所述含雜質(zhì)信息獲取部利用所述經(jīng)選擇的向量,計(jì)算針對(duì)表示在各測(cè)定時(shí)點(diǎn)的處理對(duì)象的光譜的各向量的內(nèi)積。并且,當(dāng)雜質(zhì)為多種時(shí),在各測(cè)定時(shí)點(diǎn)的作為雜質(zhì)提取用濾波器的副向量成為混合有來自于多個(gè)光譜的信號(hào)的向量。在這種情況下,即使計(jì)算單純的平均也有可能無法適當(dāng)?shù)乇硎倦s質(zhì)的含有。因此,作為進(jìn)而另一實(shí)施方式,也可設(shè)為如下構(gòu)成:所述濾波器制作部求出在各測(cè)定時(shí)點(diǎn)求出的雜質(zhì)提取用的濾波器即多個(gè)向量的聚類(clustering)平均,所述含雜質(zhì)信息獲取部利用通過所述聚類平均而獲得的向量,計(jì)算針對(duì)表示在各測(cè)定時(shí)點(diǎn)的處理對(duì)象的光譜的各向量的內(nèi)積。作為聚類平均,可使用K-均值(K-mean)法、均值偏移(meanshift)法等。并且,除此以外,還可使用移動(dòng)平均、雙邊濾波器(bilateralfilter)、卡爾曼濾波器(Karmanfilter)、粒子濾波器等考慮到時(shí)間序列性的變動(dòng)的平滑濾波器。并且,作為本發(fā)明的色譜數(shù)據(jù)處理裝置的進(jìn)而另一實(shí)施方式,還可設(shè)為如下構(gòu)成:所述濾波器制作部將如下向量確定為雜質(zhì)提取用的濾波器,所述向量是通過在表示目標(biāo)成分的光譜的向量上乘以規(guī)定的常數(shù),并從表示處理對(duì)象的光譜的向量加以減去而求出。即,根據(jù)(3)式,在所述情況下表示濾波器的向量為Bo自身。這時(shí),所述含雜質(zhì)信息獲取部只要計(jì)算在所述濾波器制作部中作為雜質(zhì)提取用的濾波器而求出的向量的二階范數(shù),將所述二階范數(shù)替代所述內(nèi)積來判斷在處理對(duì)象的光譜中有無雜質(zhì)即可。因此,可簡(jiǎn)單地即高速地進(jìn)行雜質(zhì)量的指標(biāo)值的計(jì)算。這特別有利于針對(duì)如上所述在伴隨著時(shí)間經(jīng)過的各測(cè)定時(shí)點(diǎn)獲得的處理對(duì)象的光譜分別計(jì)算雜質(zhì)量的指標(biāo)值的情況。并且,在本發(fā)明的色譜數(shù)據(jù)處理裝置中,優(yōu)選的是可設(shè)為如下構(gòu)成:當(dāng)在所述含雜質(zhì)信息獲取部中判定為存在雜質(zhì)時(shí),將在所述濾波器制作部中作為雜質(zhì)提取用的濾波器而求出的向量,例如將用I-αA表示的光譜作為殘余光譜,將所述殘余光譜作為處理對(duì)象的光譜,反復(fù)執(zhí)行由所述濾波器制作部及所述含雜質(zhì)信息獲取部進(jìn)行的處理。根據(jù)所述構(gòu)成,即使在混入至試樣中的雜質(zhì)為多種,而且存在通過由所述濾波器制作部及所述含雜質(zhì)信息獲取部進(jìn)行的一次處理無法檢測(cè)到的雜質(zhì)的情況下,也可在反復(fù)進(jìn)行多次處理的階段中檢測(cè)出所述雜質(zhì)。在本發(fā)明的色譜數(shù)據(jù)處理裝置中,本來,理想的是將利用向量表達(dá)目標(biāo)成分的光譜自身的向量設(shè)為向量A,但是一般多數(shù)情況下,目標(biāo)成分的準(zhǔn)確光譜為未知。因此,通常使用可當(dāng)作目標(biāo)成分的光譜,而不是目標(biāo)成分的光譜自身。作為一形態(tài),所述濾波器制作部可將基于如下數(shù)據(jù)的光譜當(dāng)作目標(biāo)成分的光譜,所述數(shù)據(jù)是在針對(duì)作為測(cè)定對(duì)象的試樣而獲得的三維色譜數(shù)據(jù)之中推測(cè)為存在所述目標(biāo)成分而無雜質(zhì)的期間內(nèi)獲得,將表示所述光譜的向量設(shè)為所述主向量、即向量A。再者,存在目標(biāo)成分而無雜質(zhì)的推定也可由分析者自身來進(jìn)行,但是可通過自動(dòng)判定色譜波峰的形狀來自動(dòng)進(jìn)行。并且,作為另一形態(tài),所述濾波器制作部也可將如下向量設(shè)為所述主向量、即向量A,所述向量是在基于針對(duì)作為測(cè)定對(duì)象的試樣而獲得的三維色譜數(shù)據(jù)的光譜之中,利用向量表達(dá)時(shí)的范數(shù)成為最大的向量。由此,可不預(yù)先求出目標(biāo)成分的光譜,而進(jìn)行雜質(zhì)判定處理。再者,當(dāng)向量A本來應(yīng)該是表示不含雜質(zhì)的光譜的向量,但是所述光譜含有雜質(zhì),其結(jié)果使得所制作的雜質(zhì)提取用的濾波器含有雜質(zhì)時(shí),如果按時(shí)間序列順序來繪制向量I與向量F的內(nèi)積I·F,則會(huì)在獲得已設(shè)為向量A的光譜的測(cè)定時(shí)點(diǎn)的前后出現(xiàn)波峰。這是因?yàn)轭~外減去已設(shè)為向量A的光譜中所含的雜質(zhì)的影響出現(xiàn)在所述測(cè)定時(shí)點(diǎn)的前后。因此,也可利用所述情況,來判定在某個(gè)測(cè)定時(shí)點(diǎn)或特定的時(shí)間范圍內(nèi)是否存在雜質(zhì)。即,在本發(fā)明的色譜數(shù)據(jù)處理裝置中,也可設(shè)為如下構(gòu)成:所述濾波器制作部將基于如下數(shù)據(jù)的光譜當(dāng)作目標(biāo)成分的光譜,所述數(shù)據(jù)是在針對(duì)作為測(cè)定對(duì)象的試樣而獲得的三維色譜數(shù)據(jù)之中在特定期間內(nèi)獲得,將如下向量確定為雜質(zhì)提取用的濾波器,所述向量是在表示所述目標(biāo)成分的光譜的向量上乘以規(guī)定的常數(shù),并從表示處理對(duì)象的光譜的向量加以減去而求出,所述含雜質(zhì)信息獲取部將利用如下向量來表達(dá)的光譜作為殘余光譜,所述向量是針對(duì)在包含所述特定期間的規(guī)定的時(shí)間范圍內(nèi)獲得的處理對(duì)象的各光譜,通過所述濾波器制作部作為雜質(zhì)提取用的濾波器而求出,基于所述殘余光譜判定在關(guān)于所述規(guī)定時(shí)間范圍而制作的色譜上在所述特定期間的前后是否出現(xiàn)波峰,由此判定在所述特定期間內(nèi)是否存在雜質(zhì)。[發(fā)明的效果]根據(jù)本發(fā)明的色譜數(shù)據(jù)處理裝置,可在基于三維色譜數(shù)據(jù)而制作的色譜中,確切而穩(wěn)定地判定在目標(biāo)波峰中是否含有雜質(zhì),所述三維色譜數(shù)據(jù)是通過利用PDA檢測(cè)器等多通道型檢測(cè)器或質(zhì)譜儀作為檢測(cè)器的色譜儀而收集。特別是在利用所述微分色譜的雜質(zhì)判定處理中難以適當(dāng)?shù)嘏袛嘤袩o雜質(zhì)重合的情況下,也可以確切而穩(wěn)定地判定有無雜質(zhì)重合。附圖說明圖1是包含本發(fā)明的色譜數(shù)據(jù)處理裝置的液相色譜儀的一實(shí)施例的概略構(gòu)成圖。圖2是表示本實(shí)施例的液相色譜儀中的雜質(zhì)判定處理操作的流程圖。圖3是表示在某個(gè)測(cè)定時(shí)點(diǎn)獲得的吸光光譜的一例的圖。圖4是本發(fā)明中的雜質(zhì)判定處理的原理說明圖。圖5的(a)及(b)是表示本實(shí)施例的液相色譜儀中的色譜波形及基于殘余光譜的雜質(zhì)量指標(biāo)值的時(shí)間經(jīng)過波形的一例的圖。圖6是表示在本實(shí)施例的液相色譜儀中已多次反復(fù)進(jìn)行雜質(zhì)分離處理時(shí)的雜質(zhì)量指標(biāo)值的時(shí)間經(jīng)過波形的變化的一例的圖。圖7是本實(shí)施例的液相色譜儀中的雜質(zhì)判定處理的說明圖。圖8的(a)是表示通過液相色譜儀而獲得的三維色譜數(shù)據(jù)的示意圖,(b)是表示通過液相色譜儀而獲得的波長色譜的圖。圖9的(a)~(c)是表示用于說明利用現(xiàn)有的微分色譜的雜質(zhì)判定處理中的問題的吸光光譜的一例的圖。[符號(hào)的說明]1:LC部11:流動(dòng)相容器12:送液泵13:噴射器14:管柱15:PDA檢測(cè)器16:A/D轉(zhuǎn)換器2:數(shù)據(jù)處理部21:色譜數(shù)據(jù)存儲(chǔ)部22:色譜制作部23:波峰檢測(cè)部24:雜質(zhì)判定處理部3:輸入部4:顯示部具體實(shí)施方式參照附圖對(duì)本發(fā)明的色譜數(shù)據(jù)處理裝置的一實(shí)施例進(jìn)行說明。所述色譜數(shù)據(jù)處理裝置如上所述,具有如下功能:判定基于例如利用包含PDA檢測(cè)器的液相色譜儀收集的三維色譜數(shù)據(jù)(參照?qǐng)D8的(a))而制作的色譜(參照?qǐng)D8的(b))上的波峰是否含有雜質(zhì)。首先,對(duì)本發(fā)明的色譜數(shù)據(jù)處理裝置中的雜質(zhì)判定處理的原理進(jìn)行說明。[雜質(zhì)判定處理的原理]在所述雜質(zhì)判定處理中,通過使用伴隨著時(shí)間經(jīng)過而依次獲得的作為處理對(duì)象的光譜(在以下的說明中所謂“光譜”,是指橫軸為波長,縱軸為吸光度的吸光光譜,但如上所述在質(zhì)譜等的情況下也是同樣)及目標(biāo)成分的光譜這兩者,來以高SN比制作表示目標(biāo)成分以外的雜質(zhì)的量的指標(biāo)值的時(shí)間性變化的曲線圖,并通過判定在所述曲線圖中是否存在類似色譜波峰的信號(hào),來判斷在色譜上的波峰中是否含有雜質(zhì)?,F(xiàn)在,設(shè)為利用向量I來表示在某個(gè)測(cè)定時(shí)點(diǎn)的處理對(duì)象的光譜,利用向量A來表示目標(biāo)成分的光譜(或者可當(dāng)作目標(biāo)成分的光譜的光譜)。再者,處理對(duì)象的光譜典型的是在如圖8的(a)所示的三維色譜數(shù)據(jù)之中已提取出某個(gè)時(shí)刻的吸光度的光譜,但在如下所述反復(fù)進(jìn)行雜質(zhì)分離處理的情況下,經(jīng)分離處理后的光譜成為處理對(duì)象光譜。在這里,將如圖3所示的光譜當(dāng)作規(guī)定波長范圍內(nèi)的離散的各波長下的吸光度數(shù)據(jù)的集合,將吸光光譜表示為(a(λ1)、a(λ2)、a(λ3)、……、a(λn))。a(λm)是波長λm(m=1~n)下的吸光度。以所述方式表示的光譜可表示為在n維空間內(nèi)的向量。即,所述光譜是以a(λ1)、a(λ2)、a(λ3)、……為成分的多維向量。并且同樣地,利用向量B表示雜質(zhì)的光譜。于是,表示處理對(duì)象的光譜的向量I可如上所述通過(1)式來表示。I=A+B……(1)圖4是將n維向量空間簡(jiǎn)化為二維向量空間而表示的圖,由所述(1)式表示的向量I、A、B的關(guān)系如圖4所示。將表示雜質(zhì)的光譜的向量B分解成與表示目標(biāo)成分的光譜的向量A平行的方向上的向量Ba以及與向量A正交的方向上的向量Bo來進(jìn)行考慮。在這里,進(jìn)一步考慮到與向量A正交的向量F。向量Ba與向量A平行,向量F與向量A正交,所以向量F與向量Ba正交。由于相互正交的向量的內(nèi)積成為0,所以向量F與向量Ba的內(nèi)積為0。因此,處理對(duì)象多維向量I與向量F的內(nèi)積等于向量Bo與向量F的內(nèi)積。即,已舉出的(2)式成立。I·F=Bo·F……(2)表示雜質(zhì)的光譜的向量B的長度及向量Bo的長度當(dāng)然處于比例關(guān)系,因此(2)式的右邊的Bo·F與向量B的長度,即與雜質(zhì)的量成比例。因此,(2)式的左邊的內(nèi)積I·F可用作雜質(zhì)量的指標(biāo)值u。這時(shí),向量F是用來從表示處理對(duì)象的光譜的向量I中提取雜質(zhì)成分,因此將向量F確定為雜質(zhì)提取用的濾波器。例如,當(dāng)想要調(diào)查在波長色譜等任意的色譜中所出現(xiàn)的來自于目標(biāo)成分的波峰中是否重合有雜質(zhì)時(shí),跨越從所述波峰的起點(diǎn)至終點(diǎn)的時(shí)間范圍,制作表示所述指標(biāo)值u(=內(nèi)積I·F)的時(shí)間性變化的曲線圖,如果在所述曲線圖中出現(xiàn)有類似色譜波峰的波形,則可判斷為存在雜質(zhì)。在n值非常大的n維向量空間內(nèi),與表示目標(biāo)成分的光譜的向量A正交的向量存在大致無數(shù)個(gè),但是當(dāng)將內(nèi)積I·F用作雜質(zhì)量的指標(biāo)值u時(shí),可通過如下方式來確定表示雜質(zhì)提取用濾波器的向量F的方向。即,現(xiàn)在,如果設(shè)為在表示處理對(duì)象的光譜的向量I上重疊有白色噪聲,則由所述白色噪聲引起而包含于內(nèi)積I·F中的信號(hào)成分不取決于向量F的偏角,而與其長度成比例。來自于雜質(zhì)的向量Bo與向量F所成的角越接近于直角,提取雜質(zhì)成分后由白色噪聲引起而包含于內(nèi)積I·F中的信號(hào)成分的影響越大。反而言之,為了在內(nèi)積I·F中提高來自于雜質(zhì)的信號(hào)的SN比,可使向量F與向量Bo盡可能地平行。即,向量F優(yōu)選的是以與向量Bo的余弦相似度為最大或盡可能與之接近的方式來確定其方向。當(dāng)然,由此,需要求出向量Bo,這可通過已舉出的(3)式而以解析方式算出。Bo=I-αA……(3)α=(I·A)/(A·A)再者,在雜質(zhì)的光譜中,有時(shí)會(huì)產(chǎn)生因試樣液(液相色譜儀的情況)的酸堿度(pH)變化或檢測(cè)器的非線性等所引起的變動(dòng),這有時(shí)會(huì)帶來作為內(nèi)積I·F而求出的雜質(zhì)量的指標(biāo)值u的變動(dòng),從而在所述指標(biāo)值u的曲線圖中出現(xiàn)假的類似波峰的波形。但是,由如上所述的主要因素而引起的光譜變動(dòng)決定了某種程度變化的模式(pattern),由此指標(biāo)值u中出現(xiàn)的波形變化及因雜質(zhì)的混入而產(chǎn)生的波形變化為可識(shí)別。因此,在顯示雜質(zhì)判定處理的結(jié)果時(shí),顯示由向量?jī)?nèi)積I·F表示的雜質(zhì)量的指標(biāo)值u(或表示所述指標(biāo)值u的隨時(shí)間變化的曲線圖)及以向量Bo表示的光譜這兩者,使得分析者看到它,便可判斷是否真的重合有雜質(zhì),并且可判斷已檢測(cè)出表示何種光譜的雜質(zhì)即可。這時(shí)所顯示的利用向量Bo表示的光譜并不是表示雜質(zhì)的光譜自身,而是去除了與向量A平行的向量成分Ba。因此,在一并顯示為了基于這樣的光譜鑒定雜質(zhì)或者進(jìn)行光譜比較而收錄于數(shù)據(jù)庫中的純物質(zhì)的光譜或進(jìn)行數(shù)據(jù)庫檢索時(shí),需要注意在此之前,對(duì)純物質(zhì)的光譜也需要進(jìn)行去除與向量A平行的成分的處理。當(dāng)如上所述跨越某個(gè)時(shí)間范圍制作表示由內(nèi)積I·F表示的指標(biāo)值u的時(shí)間性變化的曲線圖時(shí),由于存在在所述時(shí)間范圍內(nèi)的各測(cè)定時(shí)點(diǎn)獲得的處理對(duì)象的光譜,因此可針對(duì)每個(gè)所述光譜計(jì)算內(nèi)積I·F。在這里,由于表示向量I、F包含時(shí)間的要素,所以將考慮到時(shí)間要素的向量I、F記作I(t)、F(t)。表示處理對(duì)象的光譜的向量I(t)在每個(gè)測(cè)定時(shí)點(diǎn)存在,但不一定在每個(gè)測(cè)定時(shí)點(diǎn)都需要表示濾波器的向量F(t)。計(jì)算各測(cè)定時(shí)點(diǎn)的內(nèi)積I(t)·F(t)時(shí)所使用的F(t)可大致劃分為以下兩者中的任一者。(1)直接使用在各測(cè)定時(shí)點(diǎn)算出的F(t),乘以表示在各測(cè)定時(shí)點(diǎn)獲得的處理對(duì)象的光譜的向量I(t),由此計(jì)算內(nèi)積I(t)·F(t)。(2)不直接使用在各測(cè)定時(shí)點(diǎn)算出的F(t),根據(jù)每個(gè)測(cè)定時(shí)點(diǎn)的F(t)的值算出用于計(jì)算內(nèi)積的F(t)'。例如,計(jì)算在某個(gè)規(guī)定時(shí)間內(nèi)的多個(gè)測(cè)定時(shí)點(diǎn)獲得的向量F(t)的平均向量,將其作為向量F而乘以表示在各測(cè)定時(shí)點(diǎn)獲得的處理對(duì)象的光譜的向量I(t),由此計(jì)算內(nèi)積I(t)·F。由此,可求出表示對(duì)噪聲魯棒性高的平均的濾波器的向量F。但是,當(dāng)含有多種雜質(zhì)時(shí),所述多種雜質(zhì)出現(xiàn)的時(shí)間并不相同,因此在各測(cè)定時(shí)點(diǎn)的F(t)成為復(fù)雜地混合有來自于多種雜質(zhì)的光譜的信號(hào)。在這樣的情況下,即使計(jì)算出所述單純的平均向量也無法獲得表示適當(dāng)?shù)臑V波器的向量F,因此可使用所謂聚類平均等而不是單純平均。除了作為聚類平均而眾所周知的K-均值法、均值偏移法等以外,還可使用移動(dòng)平均、雙邊濾波器、卡爾曼濾波器、粒子濾波器(序貫蒙特卡洛法(sequentialMonteCarlomethod))等將時(shí)間序列變動(dòng)考慮在內(nèi)的各種平滑濾波器。[目標(biāo)成分的光譜的替代]在所述說明中,為了算出雜質(zhì)量的指標(biāo)值u而使用表示目標(biāo)成分的光譜的向量A,但是在多數(shù)情況下,目標(biāo)成分的準(zhǔn)確的光譜未知,所以獲取所述光譜會(huì)花費(fèi)工夫。因此,就實(shí)用性而言,例如可通過如下所述的方式,利用通過對(duì)試樣的分析而獲得的信號(hào)(即每個(gè)測(cè)定時(shí)點(diǎn)的光譜)模擬地生成目標(biāo)成分的光譜。通常,雜質(zhì)的濃度低于目標(biāo)成分的濃度,所以如圖7所示,在色譜中雜質(zhì)的波峰寬度窄于目標(biāo)成分的波峰寬度。由此,可以說在每個(gè)通過分析而獲得的測(cè)定時(shí)點(diǎn)的光譜之中,包含有在目標(biāo)成分的光譜中混合有雜質(zhì)的光譜的狀態(tài)的光譜與只有目標(biāo)成分的光譜這兩者的可能性高。因此,例如,可設(shè)為在通過分析而獲得的色譜數(shù)據(jù)之中,提取可推測(cè)為一方面含有目標(biāo)成分,另一方面未混合有雜質(zhì)的在特定的時(shí)間范圍內(nèi)所含的數(shù)據(jù),將由此獲得的光譜當(dāng)作目標(biāo)成分的光譜。并且,也可將使通過分析而獲得的數(shù)據(jù)沿時(shí)間方向平滑化之后提取特定的時(shí)間范圍的數(shù)據(jù)而獲得的光譜,或者,對(duì)通過分析而獲得的數(shù)據(jù)之中的特定的時(shí)間范圍的平均進(jìn)行計(jì)算而獲得的光譜當(dāng)作目標(biāo)成分的光譜。這時(shí)提取數(shù)據(jù)的時(shí)間范圍既可設(shè)為可由分析者指定,或者,也可設(shè)為通過后述的辨別法來辨別包含雜質(zhì)波峰的時(shí)間,基于其結(jié)果而自動(dòng)選定被推測(cè)為雜質(zhì)少的時(shí)間范圍。進(jìn)而此外,如果只將確認(rèn)有無雜質(zhì)重合作為目標(biāo),而不需要高精度地求出雜質(zhì)的含量,那么表示雜質(zhì)量的指標(biāo)值u的隨時(shí)間變化的曲線圖中出現(xiàn)的波峰被分成兩個(gè)(理由將在后文敘述)則不是大問題。因此,也可容許混合有雜質(zhì)或出現(xiàn)光譜變動(dòng),在通過分析而獲得的光譜之中選擇單純地SN比為最佳的光譜、即通常利用向量表達(dá)時(shí)的范數(shù)為最大的光譜而當(dāng)作目標(biāo)成分的光譜。在這里,利用圖5的(a)及(b)來說明如下現(xiàn)象,所述現(xiàn)象是在基于在某個(gè)測(cè)定時(shí)點(diǎn)或某個(gè)時(shí)間范圍內(nèi)所含的測(cè)定時(shí)點(diǎn)獲得的光譜而制作的雜質(zhì)提取用濾波器是基于包含雜質(zhì)的光譜而不是由單一成分構(gòu)成的光譜而制作的雜質(zhì)提取用濾波器時(shí)所產(chǎn)生。圖5的(a)及(b)是表示色譜波形及基于殘余光譜的雜質(zhì)量指標(biāo)值u的時(shí)間經(jīng)過波形的一例的圖。圖5的(a)中以P1表示的指標(biāo)值是將在測(cè)定時(shí)點(diǎn)42獲得的光譜(在目標(biāo)成分中混合有雜質(zhì)的狀態(tài)的光譜)當(dāng)作目標(biāo)成分的光譜,求出表示濾波器的向量F,而繪制內(nèi)積I(t)·F的隨時(shí)間變化的曲線。與此相對(duì),圖5的(b)中以P2表示的指標(biāo)值是將在測(cè)定時(shí)點(diǎn)60獲得的光譜(不含雜質(zhì)的只有目標(biāo)成分的光譜)當(dāng)作目標(biāo)成分的光譜,求出表示濾波器的向量F,而繪制內(nèi)積I(t)·F的隨時(shí)間變化的曲線。在圖5的(a)中,可知夾著選擇目標(biāo)成分的光譜的測(cè)定時(shí)點(diǎn)而在其前后分別出現(xiàn)有波峰。這是上述波峰的分裂。這時(shí),表示內(nèi)積I(t)·F的隨時(shí)間變化的曲線圖上的波峰的形狀沒有準(zhǔn)確地表示雜質(zhì)量,所以難以求出雜質(zhì)量。但是,利用所述現(xiàn)象,在內(nèi)積I(t)·F的曲線圖中在目標(biāo)成分的前后兩方存在波峰時(shí),可判斷為已將包含雜質(zhì)的光譜設(shè)定為目標(biāo)成分光譜。另一方面,如圖5的(b)所示,當(dāng)選擇不含雜質(zhì)的光譜作為目標(biāo)成分的光譜時(shí),在表示內(nèi)積I(t)·F的隨時(shí)間變化的曲線圖中會(huì)出現(xiàn)高斯波形狀的波峰,這可推測(cè)為準(zhǔn)確地表示了雜質(zhì)量。[雜質(zhì)為多種的情況下的雜質(zhì)分離處理]圖5的(a)及(b)所示的示例是雜質(zhì)只有一種的情況,但如上所述,混入至試樣中的雜質(zhì)并不限于一種,也可以是多種?,F(xiàn)在作為一例,考慮如下情況:除了目標(biāo)成分a以外,還存在雜質(zhì)b及雜質(zhì)c,并且雜質(zhì)c與目標(biāo)成分a及雜質(zhì)b相比為微量。當(dāng)作為表示雜質(zhì)提取用濾波器的向量F由在規(guī)定時(shí)間范圍內(nèi)的多測(cè)定點(diǎn)的F(t)的平均值表示時(shí),向量I(t)與求得的向量F相似.所述I(t)是由表示目標(biāo)成分a的光譜的向量A和表示不純物b的光譜的向量B構(gòu)成。這時(shí),向量F成為與向量Bo平行的向量。在這里,當(dāng)設(shè)為向量B與表示雜質(zhì)c的光譜的向量C正交時(shí),表示內(nèi)積I(t)·F的隨時(shí)間變化的曲線圖中完全不可能檢測(cè)到向量C。并且,即使在每個(gè)測(cè)定時(shí)間求出表示雜質(zhì)提取用濾波器的向量F(t)的情況下,來自于雜質(zhì)c的波峰也非常小,例如,當(dāng)與來自于雜質(zhì)b的波峰的下擺部等的信號(hào)變動(dòng)大的部分重合時(shí),來自于雜質(zhì)c的波峰的檢測(cè)變得非常困難。因此,在這種情況下,即在假設(shè)混入有多種雜質(zhì),需要對(duì)各種雜質(zhì)進(jìn)行檢測(cè)時(shí),可采用如下所述的順序。即,如由所述(3)式可知,I-αA表示雜質(zhì)的量,因此(2)式所示的處理,即在處理對(duì)象的光譜上乘以濾波器的處理可以說是雜質(zhì)分離處理。這時(shí),利用I-αA或I(t)-αA的向量表示的光譜可以說是已去除目標(biāo)成分或一種以上的雜質(zhì)的剩下的殘余光譜。當(dāng)在試樣中含有多種雜質(zhì)時(shí),可將在第n次處理中算出的I(t)-αA(表示殘余光譜的向量),作為表示用于進(jìn)行第n+1次處理的處理對(duì)象的光譜的向量I(t)而執(zhí)行雜質(zhì)分離處理。在這里,將所述方法稱為多段光譜殘?jiān)?。圖6是表示基于執(zhí)行如上所述的多段光譜殘?jiān)〞r(shí)的殘余光譜的信號(hào)波形的圖。圖中的O表示原來的色譜波形,Q1~Q4分別表示第n=1~4次目標(biāo)|I(t)|。在Q1中,被認(rèn)為本來存在看起來是Q3的小的波峰,但是即使觀察Q1的波形,也難以找到應(yīng)該包含于其中的在Q3中觀測(cè)到的波峰。與此相對(duì),通過使用所述多段光譜殘?jiān)?,還可檢測(cè)到最初無法找到的雜質(zhì)的波峰。并且,在多段光譜殘?jiān)ㄖ?,可通過關(guān)于在第n次處理中獲得的|I(t)|與在第n+1次處理中獲得的|I(t)|的差(光譜殘?jiān)罘?檢測(cè)波峰的有無,來在各階段中判斷雜質(zhì)的有無。例如在圖6中,關(guān)于在以Q2表示的曲線的左側(cè)作為凸部而檢測(cè)到的雜質(zhì),由于在Q2中混合有Q3的信號(hào),所以難以判斷是真的波峰還是噪聲變動(dòng)。但是,如果通過去除Q3而求出Q4中所示的波形,則可識(shí)別出存在明顯只偏在于左側(cè)的成分。再者,在進(jìn)行所述去除操作時(shí),由于表示基于Q2的濾波器的向量F與表示基于Q3的濾波器的向量F不同,所以即使是來自于同一成分的波峰,波峰高度也會(huì)成常數(shù)倍。因此,在實(shí)際的去除操作中,可設(shè)為只著眼于波峰部分,利用最小平方法而求出最適當(dāng)?shù)某?shù),并將所述常數(shù)乘以各強(qiáng)度值之后進(jìn)行去除處理。當(dāng)然,也可使用一般的應(yīng)對(duì)基線(baseline)變動(dòng)的波峰高度推定法,所述波峰高度推定法是不使用單純的最小平方法,而使用對(duì)計(jì)算F(t)的二階微分而求出的波形應(yīng)用最小平方法,或利用將使經(jīng)提取的波峰規(guī)格化而成者作為內(nèi)核(kernel)的匹配濾波器(matchedfilter)來推測(cè)高度等。反復(fù)進(jìn)行如以上所說明的處理,直至獲得如圖6中以Q4表示的波形般不具有明確的波峰的殘?jiān)盘?hào)波形為止,由此即使在存在多種雜質(zhì)的情況下,也可將目標(biāo)成分與雜質(zhì)完全分離。如果除去因?yàn)閜H變動(dòng)或檢測(cè)器的線性程度低等所引起的假雜質(zhì)波峰的產(chǎn)生,則在利用多段光譜殘?jiān)▽?duì)針對(duì)包含m種物質(zhì)的試樣而獲得的測(cè)量信號(hào)進(jìn)行處理的情況下,只要反復(fù)執(zhí)行m+1次雜質(zhì)分離處理,便可對(duì)m種物質(zhì)進(jìn)行分離。[實(shí)施所述原理的雜質(zhì)判定處理的實(shí)施例的構(gòu)成及操作]其次,參照?qǐng)D1、圖2對(duì)包含本發(fā)明的色譜數(shù)據(jù)處理裝置的液相色譜儀的一實(shí)施例進(jìn)行說明。圖1是本實(shí)施例中的液相色譜儀的概略構(gòu)成圖。在用于收集三維色譜數(shù)據(jù)的LC部1中,送液泵12從流動(dòng)相容器11抽吸流動(dòng)相,以固定的流量送入至噴射器(injector)13。噴射器13以規(guī)定的時(shí)間將試樣液注入至流動(dòng)相中。試樣液通過流動(dòng)相而送入至管柱14,在通過管柱14期間使試樣液中的各成分沿時(shí)間方向分離,并從管柱14溶出。在管柱14的出口設(shè)置有PDA檢測(cè)器15。PDA檢測(cè)器15將來自未圖示的光源的光照射至溶出液,使穿透溶出液的光波長分散而利用線性傳感器(linearsensor)大致同時(shí)檢測(cè)各波長的光的強(qiáng)度。由所述PDA檢測(cè)器15反復(fù)獲得的檢測(cè)信號(hào)通過模擬/數(shù)字(AnalogtoDigital,A/D)轉(zhuǎn)換器16而轉(zhuǎn)換成數(shù)字?jǐn)?shù)據(jù)(digitaldata)之后,作為三維色譜數(shù)據(jù)發(fā)送至數(shù)據(jù)處理部2。數(shù)據(jù)處理部2包括:色譜數(shù)據(jù)存儲(chǔ)部21,用于存儲(chǔ)三維色譜數(shù)據(jù);色譜制作部22,利用三維色譜數(shù)據(jù)制作表示規(guī)定的波長下的吸光度的隨時(shí)間變化的波長色譜;波峰檢測(cè)部23,對(duì)所述波長色譜中的波峰進(jìn)行檢測(cè);以及雜質(zhì)判定處理部24,判定在檢測(cè)到的波峰之中由分析者指定的目標(biāo)波峰中是否存在雜質(zhì)。所述雜質(zhì)判定處理部24是進(jìn)行如上所述的特征性處理的功能塊。并且,在數(shù)據(jù)處理部2上連接有輸入部3及顯示部4,輸入部3是為了分析者輸入設(shè)定目標(biāo)成分的吸收波長等數(shù)據(jù)處理所需要的信息等而被操作。并且,顯示部4是用于顯示色譜、吸光光譜及雜質(zhì)判定結(jié)果等各種信息的構(gòu)件。再者,數(shù)據(jù)處理部2或未圖示的控制部的功能的一部分或全部可通過執(zhí)行安裝在個(gè)人計(jì)算機(jī)或工作站(workstation)中的專用的控制及處理軟件來實(shí)現(xiàn)。這時(shí),輸入部3是個(gè)人計(jì)算機(jī)或工作站的標(biāo)準(zhǔn)裝備即鍵盤或鼠標(biāo)等指向構(gòu)件(pointing)等,顯示部4是一般的液晶顯示器等。其次,參照?qǐng)D2所示的流程圖,對(duì)本實(shí)施例的液相色譜儀中的特征性的數(shù)據(jù)處理操作進(jìn)行說明。在LC部1中執(zhí)行對(duì)目標(biāo)試樣的色譜分析,將表示規(guī)定的波長范圍內(nèi)的吸光光譜的隨時(shí)間變化的三維色譜數(shù)據(jù)(參照?qǐng)D8的(a))從PDA檢測(cè)器15發(fā)送至數(shù)據(jù)處理部2,并存儲(chǔ)在色譜數(shù)據(jù)存儲(chǔ)部21中。在色譜制作部22中,制作基于這樣存儲(chǔ)的三維色譜數(shù)據(jù)而指定的特定波長或波長范圍內(nèi)的波長色譜,波峰檢測(cè)部23在所述色譜上執(zhí)行波峰檢測(cè)。當(dāng)分析者通過輸入部3而指定這樣檢測(cè)到的波峰之后,命令執(zhí)行雜質(zhì)判定處理時(shí),實(shí)施如以下所述的處理。首先,雜質(zhì)判定處理部24從色譜數(shù)據(jù)存儲(chǔ)部21讀入針對(duì)所指定的波峰的起點(diǎn)ts與終點(diǎn)te之間的各測(cè)定時(shí)點(diǎn)的色譜數(shù)據(jù)(光譜數(shù)據(jù))(步驟S1)。由此,準(zhǔn)備好表示處理對(duì)象的光譜的向量I(t)(其中t為ts~te)。其次,雜質(zhì)判定處理部24設(shè)定用于求出向量A的目標(biāo)成分的光譜(步驟S2)。如上所述設(shè)定目標(biāo)成分的光譜有若干個(gè)方法,當(dāng)在數(shù)據(jù)庫等中存儲(chǔ)有目標(biāo)成分的光譜時(shí),只要將其讀入即可。在這里,由于目標(biāo)成分的光譜為未知,并且自動(dòng)反復(fù)進(jìn)行所述光譜的設(shè)定,所以采用不需要分析者的作業(yè)或判斷,而選擇表示可高速地執(zhí)行的最大范數(shù)的光譜的方法。即,將在即將開始之前所執(zhí)行的處理的結(jié)果即雜質(zhì)量的指標(biāo)值u=I(t)·F為最大的測(cè)定時(shí)點(diǎn)獲得的吸光光譜直接當(dāng)作下一次的目標(biāo)成分的光譜來進(jìn)行設(shè)定。由此,也準(zhǔn)備好表示目標(biāo)成分的光譜的向量A。再者,當(dāng)執(zhí)行第一次處理,即第一次執(zhí)行步驟S2的處理時(shí),計(jì)算在步驟S1中準(zhǔn)備的向量I(t)的二階范數(shù),選定在二階范數(shù)為最大的測(cè)定時(shí)點(diǎn)的目標(biāo)成分的光譜。當(dāng)然,也可設(shè)為通過分析者的手動(dòng)指定來設(shè)定目標(biāo)成分的光譜,還可如上所述,設(shè)為搜索不含雜質(zhì)的光譜,設(shè)定在不含雜質(zhì)的多個(gè)光譜之中雜質(zhì)量的指標(biāo)值或二階范數(shù)為最大的目標(biāo)成分的光譜。如果處理對(duì)象的光譜(向量I(t))及目標(biāo)成分的光譜(向量A)已確定,則如上所述確定雜質(zhì)提取用的濾波器,計(jì)算內(nèi)積I(t)·F,由此求出從處理對(duì)象光譜中去除目標(biāo)成分的光譜而反映出雜質(zhì)量的殘余光譜(步驟S3)。在這里,是采用注重運(yùn)算的高速性,將在各測(cè)定時(shí)點(diǎn)的I(t)-αA直接用作向量F(t)的方法。這時(shí),可單純地對(duì)算式進(jìn)行變形,只是通過算出I(t)-αA的二階范數(shù)來代用作為雜質(zhì)量的指標(biāo)值u的向量?jī)?nèi)積I(t)·F的計(jì)算。當(dāng)然,也可不求出在每個(gè)測(cè)定時(shí)點(diǎn)的F(t),而使用向量F(t)的平均值或移動(dòng)平均等所述各種修改(modification)。通過判定如上所述而求出的殘余光譜與在前一次的處理中所獲得的殘余光譜的差,即光譜殘?jiān)罘值亩A范數(shù)或在各次計(jì)算中獲得的雜質(zhì)量的指標(biāo)值的平方根中的任一者中是否存在波峰,來判斷是否存在來自于雜質(zhì)的波峰(步驟S4)。雜質(zhì)量的平方根或光譜殘?jiān)罘值亩A范數(shù)表示對(duì)于白色噪聲為固定的分布,因此可通過確認(rèn)這些值的平均值及有無超出基于標(biāo)準(zhǔn)偏差的范圍的值,來確認(rèn)波峰的有無。當(dāng)然,并不限于所述方法,也可應(yīng)用一般的色譜波峰的檢測(cè)算法(algorithm)來確認(rèn)波峰的有無。并且,當(dāng)判定為存在雜質(zhì)波峰時(shí),從步驟S5返回到S2,反復(fù)進(jìn)行目標(biāo)成分的光譜的設(shè)定及目標(biāo)成分的光譜的去除。即,實(shí)施所述多段光譜殘?jiān)?。另一方面,?dāng)在步驟S5中判定為無雜質(zhì)波峰時(shí),基于在此之前的判定結(jié)果將最終的雜質(zhì)判定處理結(jié)果顯示在顯示部4中,當(dāng)判定為存在雜質(zhì)時(shí)將各光譜的殘?jiān)罘忠诧@示在顯示部4中(步驟S6)。由此,分析者不僅可知道在目標(biāo)波峰中是否重合有雜質(zhì),而且可知道雜質(zhì)的量。再者,所述實(shí)施例顯然是本發(fā)明的一例,即使在本發(fā)明的主旨的范圍內(nèi)適當(dāng)進(jìn)行變形、追加、修正也包含在本申請(qǐng)的權(quán)利要求書中。例如,對(duì)本發(fā)明的色譜數(shù)據(jù)處理裝置的處理對(duì)象的三維色譜數(shù)據(jù)進(jìn)行獲取的色譜儀的檢測(cè)器也可不是PDA檢測(cè)器等多通道型檢測(cè)器,而還可以是可進(jìn)行高速的波長掃描的紫外可見分光光度計(jì)、紅外分光光度計(jì)、近紅外分光光度計(jì)、熒光分光光度計(jì)等。并且,也可以是將質(zhì)譜儀作為檢測(cè)器的液相色譜質(zhì)譜儀。并且,色譜儀也可以是氣相色譜儀而不是液相色譜儀。并且,如上所述,顯然也可將本發(fā)明應(yīng)用于如下裝置,所述裝置不對(duì)利用檢測(cè)器對(duì)經(jīng)色譜儀的管柱分離的試樣進(jìn)行檢測(cè)而獲得的數(shù)據(jù)進(jìn)行處理,而對(duì)利用PDA檢測(cè)器或質(zhì)譜儀等檢測(cè)器對(duì)通過FIA法不進(jìn)行成分分離而導(dǎo)入的試樣中的成分進(jìn)行檢測(cè)而獲得的數(shù)據(jù)進(jìn)行處理。