摘 要:提出一種相關(guān)系數(shù)模糊判別法鑒別相似的墨粉傅里葉變換紅外(FTIR)光譜圖形,改變目測圖形相似程度和圖形吸收峰波數(shù)簡單比較的主觀粗略的鑒別方法,以達到客觀準確認定激光打印機型的目的。它能廣泛用于與激光打印機有關(guān)的司法案件的鑒定。首先運用FTIR技術(shù)獲取相關(guān)激光打印機墨粉的FTIR光譜圖形;再將FTIR圖形的各吸收峰波數(shù),及其一階、二階導(dǎo)數(shù)和自卷積作為指標,計算圖形的相關(guān)系數(shù);最后,以相關(guān)系數(shù)為圖形相似的指標和建立模糊關(guān)系函數(shù)的基礎(chǔ),用模糊數(shù)學(xué)綜合判別法確定激光打印機型的鑒別結(jié)果。實驗結(jié)果達到算法設(shè)計的目的。
關(guān)鍵詞:墨粉;傅里葉變換紅外光譜圖形;相關(guān)系數(shù);模糊鑒別;激光打印機
中圖分類號:TP391 文獻標志碼:A
文章編號:1001-3695(2008)07-2200-03
Fuzzy identification of laser printer based on similarblack-powder FTIR spectrum graphics
WANG Ning1, 2, HAN Guo-qiang1
(1.School of Computer Science Engineering, South China University of Technology, Guangzhou 510641, China;2. Guangzhou Public Security Bureau, Guangzhou 510030, China)
Abstract:This paper proposed a new method of correlation coefficient fuzzy identification to evaluate the similar Fourier transform infrared (FTIR) spectrum graphics. It improved the subjective and rough identification in eyeballing the similar degree of graphics and the simple comparison of wavenumbers of graphics. So it could identify the type of laser printer objectively and accurately and could be used widely in the case identification involved with laser printer. First, the FTIR spectrum graphics of black-powder of laser printer was got by FTIR technique. Then, as the base of graphics similarity and fuzzy relationship function, the correlation coefficient of graphics was calculated by the wavenumbers, the first derivative, the second derivative and the self-deconvolution of FTIR spectrum graphics. Finally, the type of laser printer was identified by fuzzy synthetic evaluation. The result of experiment shows that the algorithm can meet the aim of design.
Key words: black-powder; FTIR spectrum graphics; correlation coefficient; fuzzy identification; laser printer
在一些司法案件鑒定中,需要快速判斷打印文件(字)的打印機類型。這類方法至今甚少報道。傅里葉變換紅外(FTIR)光譜檢測是一種快速、精確、可靠的檢測打印機墨粉成分的方法[1]。目前,在運用FTIR光譜圖形對激光打印機墨粉進行成分鑒定時,主要采用主觀目測其FTIR光譜圖形是否相似和光譜吸收峰波數(shù)簡單比較是否接近作為判斷標準。這兩個指標隨意性大,不夠精確和客觀。鑒別結(jié)果,也常隨鑒定者的知識技術(shù)水平和經(jīng)驗有所差異。特別是微量物證鑒定,差別十分微小的鑒別有時難以作出判斷,如圖1、2所示。
本研究希望利用FTIR技術(shù)獲取激光打印機墨粉的FTIR光譜圖形,通過計算圖形間的相關(guān)系數(shù)并運用模糊綜合判別法判斷圖形的相似程度,從而達到鑒別激光打印機類型的目的,其數(shù)學(xué)模型由計算機程序?qū)崿F(xiàn)。在研究中,充分利用了光譜圖形的多種信息,以待檢測物與各標準樣品FTIR圖形對應(yīng)的各吸收峰波數(shù)值,及其一階、二階導(dǎo)數(shù)和自卷積值作為指標,計算FTIR圖形間的相關(guān)系數(shù)。以此作為FTIR圖形相似的指標和建立模糊關(guān)系函數(shù)的基礎(chǔ),用模糊數(shù)學(xué)綜合判別法[2]確定鑒別結(jié)果。將整個過程用計算機程序?qū)崿F(xiàn),達到自動鑒別激光打印機類型的目的。實驗中,以方正文景A230、利盟E120、富士施樂3117、佳能LBP2900、惠普CLJ5500dn、三星ML2010、聯(lián)想LJ2000、兄弟2040、愛普生6200L和柯尼卡美能達1400W,十種常用激光打印機墨粉為樣本,進行了實際案例檢測,驗證算法的正確性和程序功能模塊的性能,結(jié)果達到設(shè)計的預(yù)期要求。
1 激光打印機鑒別方法
1.1 檢測激光打印機墨粉的FTIR光譜圖形
1.1.1 儀器和操作方法
用尼高力Impact410型傅里葉變換紅外光譜儀,對每一激光打印機墨粉樣品掃描128次,重復(fù)3次。經(jīng)自動基線和自動平滑校正,扣除水及CO2的影響后,得到墨粉的標準FTIR光譜圖形,并進行一階導(dǎo)數(shù)、二階導(dǎo)數(shù)變換處理,以及自卷積處理。
1.1.2 檢測樣品墨粉FTIR圖形
以十種常用激光打印機墨粉為標準樣品,檢測這些標準樣品的FTIR吸收峰波數(shù)及其一階、二階導(dǎo)數(shù)和自卷積值,并分別建立相應(yīng)的譜圖庫和各指標數(shù)據(jù)庫如表1所示。
1.1.3 檢測待測墨粉FTIR圖形
檢測方法同樣品。
1.2 計算待測墨粉與標準墨粉FTIR圖形間的相關(guān)系數(shù)
1.2.1 指標設(shè)定
同一種待測物即使在不同時間檢測,其FTIR光譜圖形及其一階、二階導(dǎo)數(shù)和自卷積的吸收峰的值,應(yīng)是對應(yīng)且相同的。不同待測物的FTIR吸收峰,多數(shù)是不完全對應(yīng)的[1]。在FTIR圖形上,許多吸收峰不對應(yīng),這就會造成檢測結(jié)果的差異,恰好也是檢測要達到的目的。不過,在計算FTIR圖形的相關(guān)系數(shù)時,必須要求輸入相互對應(yīng)的吸收峰波數(shù)值,波數(shù)值相差很大時,可以認為是不對應(yīng)的吸收峰波數(shù)值。這種情況下可以視對方的值為0。但是,多大的波數(shù)差才可以視對方的值為0?根據(jù)以往使用FTIR技術(shù)中的主觀判斷界限,均不超過12,所以這里暫定FTIR吸收峰波數(shù)差h對應(yīng)值區(qū)間為
區(qū)間內(nèi)視為對應(yīng)值,否則對應(yīng)值為0。FTIR圖形的一階、二階導(dǎo)數(shù)和自卷積值,也以該區(qū)間為準。以待測激光打印機墨粉與方正文景A230墨粉FTIR各吸收峰波數(shù)對應(yīng)值為例,說明取值方法如表2所示。
1.2.2 待測墨粉與樣品墨粉FTIR圖形相關(guān)系數(shù)計算[3,4]
設(shè)待測墨粉與樣品墨粉的FTIR吸收峰波數(shù)分別以a和b表示。ak′、 ak″、bk′、 b
1.2.3 待測墨粉與樣品墨粉相關(guān)系數(shù)表
用Ai表示FTIR的波數(shù)、及其一階、二階導(dǎo)數(shù)和自卷積值的四個指標,即令:FTIR波數(shù)為A1 , FTIR波數(shù)一階導(dǎo)數(shù)為A2 ,F(xiàn)TIR二階導(dǎo)數(shù)為A3 ,F(xiàn)TIR自卷積為A4 ,制作待測墨粉與N種樣品墨粉各指標相關(guān)系數(shù)的絕對值rij如表4所示[5]。
1.3 確定指標的權(quán)重系數(shù)
FTIR波數(shù)及其一階、二階導(dǎo)數(shù)和自卷積的權(quán)重,以統(tǒng)計學(xué)中的差異顯著性級別設(shè)計。完全相似的圖形應(yīng)有如下數(shù)據(jù):FTIR波數(shù)權(quán)重為0.951,一階導(dǎo)數(shù)權(quán)重為0.039,二階導(dǎo)數(shù)權(quán)重為0.005,自卷積權(quán)重為0.005 。
FTIR波數(shù)權(quán)重與一階導(dǎo)數(shù)權(quán)重之和為0.99,F(xiàn)TIR波數(shù)權(quán)重與一階導(dǎo)數(shù)權(quán)重、二階導(dǎo)數(shù)權(quán)重及自卷積權(quán)重之和為1.00。FTIR波數(shù)及其一階、二階導(dǎo)數(shù)和自卷積的權(quán)重系數(shù)這四個指標如表5所示。
a為A上的模糊集合,其隸屬函數(shù)為μA (ai ),為指標的隸屬函數(shù)[6]。在實際應(yīng)用上,構(gòu)成相應(yīng)指標的權(quán)重系數(shù),可用權(quán)重系數(shù)矩陣A表示。由表5建立權(quán)重矩陣:
A= (A1A2 A3 A4 )= ( 0.9510.0390.005 0.005) (5)
1.4.2 樣品墨粉隸屬函數(shù)及隸屬矩陣
設(shè)已知標準樣品墨粉的集合:
B ={bj}( j = 1,2,…, n ) (6)
b為B上的模糊集合,其隸屬函數(shù)μB ( bj), 為待檢測墨粉隸屬標準樣品墨粉的隸屬函數(shù)[7]。其判別隸屬矩陣B為
B =(B1 B2 … Bn )(7)
1.4.3 模糊關(guān)系函數(shù)及模糊關(guān)系矩陣[8,9]
通過模糊變換,由集合A上的元素a變換為集合B上的b 中間需要一個與集合A和集合B均有關(guān)的模糊關(guān)系集合R,它的函數(shù)形式就是模糊關(guān)系函數(shù)R(ai , bj )。由表4可建立模糊關(guān)系矩陣R
1.4.4 模糊變換
指標隸屬函數(shù)μA(ai )與模糊關(guān)系函數(shù)R(ai , bj )的乘積,進行一次模糊變換,變換為樣品墨粉隸屬函數(shù)μB ( bj)。根據(jù)該樣品墨粉隸屬函數(shù)實現(xiàn)模糊綜合評判[10~12]。
這個過程由矩陣A與矩陣R的積來實現(xiàn):
1.4.5 模糊判別[13]
由最大隸屬性原則可得到模糊判別的結(jié)果:
2 激光打印機模糊鑒別實驗
2.1 計算待測墨粉與樣品墨粉FTIR圖形的相關(guān)系數(shù)
根據(jù)式(2)計算待檢測激光打印機墨粉FTIR圖形與十種樣品激光打印機墨粉FTIR圖形的相關(guān)系數(shù), 計算出的相關(guān)系數(shù)按表1順序及表4的關(guān)系進行排列如表6所示。
2.2 建立待測墨粉與樣品墨粉的模糊關(guān)系矩陣
由模糊關(guān)系函數(shù)R(ai , bj ),式(4)以及表6的數(shù)據(jù)可以建立待測激光打印機墨粉FTIR圖形與樣品激光打印機墨粉FTIR圖形指標的模糊關(guān)系矩陣:
2.3 激光打印機模糊鑒別
根據(jù)式(9)和(10),權(quán)重矩陣A與上述模糊關(guān)系矩陣R(即式(12))進行模糊變換,可以得到模糊判別矩陣B:
由最大隸屬性原則及式(11)和(13)可得到激光打印機類型模糊判別的結(jié)果,即
B0=max(B1 B2 B3 B4 B5 B6 B7 B8 B9 B10)=B4=0.951(14)
結(jié)果,待檢測墨粉,被判斷為聯(lián)想LJ2000激光打印機墨粉,與最終結(jié)果一致。
3 激光打印機模糊鑒別軟件設(shè)計方案
激光打印機模糊鑒別軟件架構(gòu)于微軟的.NET平臺,以組件模式設(shè)計并實現(xiàn)。這樣可大大提高軟件程序的運行效率和性能,增強軟件的穩(wěn)定性和安全性,提高程序模塊的可移植性、可維護性和重用性。該軟件設(shè)計方案如圖3所示。
4 結(jié)束語
本文用FTIR光譜圖形的各吸收峰波數(shù)值及其一階、二階導(dǎo)數(shù)和自卷積值作為指標,綜合進行相似性計算。指標間相關(guān)系數(shù)作為圖形相似的指標和建立模糊關(guān)系函數(shù)的基礎(chǔ)。最后用模糊數(shù)學(xué)綜合判別法確定鑒別結(jié)果。該算法實用,可操作性好。實驗結(jié)果表明,基于激光打印機墨粉FTIR光譜圖形相似的模糊鑒別法鑒別激光打印機的類型是可行的,算法數(shù)學(xué)模型是正確的。
參考文獻:
[1] 吳瑾光.近代傅里葉變換紅外光譜技術(shù)及應(yīng)用[M].北京:科學(xué)技術(shù)文獻出版社,1994:3-131.
[2]胡寶清.模糊理論基礎(chǔ)[M].武漢:武漢大學(xué)出版社, 2004:179- 218.
[3]CHEN S M. Similarity measure between vague sets and between elements[J]. IEEE Trans on Systems, Man and Cybernetics, 1997, 27(1):153-158.
[4]CHEN S M,YEH M S,HISIAO P Y. A comparison of similarity measures of fuzzy values[J]. Fuzzy Sets and Systems, 1995, 72(1):79-80.
[5]REZAEE R M, GOEDHART B,LELIEVELDT B P F, et al. Fuzzy feature selection[J]. Pattern Recognition, 1999,32(3):2011-2019.
[6]李云,吳中福,葉春曉,等.基于擴張矩陣的模糊特征選擇算法[J].模式識別與人工智能, 2004, 17(2):141-145.
[7]BEZDEK J C,PALS K.Fuzzy models for pattern recognition methods that search for patterns in data[M].New York:IEEE Press,1992.
[8]PAPPS C P, KARACAPILIDIS N I. A comparative assessment of measures of similarity of fuzzy values[J]. Fuzzy Sets and Systems, 1993,56(2):171-174.
[9]劉華文.模糊模式識別的基礎(chǔ)—相似度量[J].模式識別與人工智能,2004,17(2):141-145.
[10]SETNES M, BABUSKA R. Fuzzy relation classifier trained by fuzzy cluster[J]. IEEE Trans on Systems, Man and Cybernetics-Part B: Cybernetics, 1999, 29(5):619-625.
[11]DUDA R O,HART P E,STORK D G. 模式分類[M].李宏東,姚天翔,等譯.2版. 北京:機械工業(yè)出版社, 2005:157-160.
[12]GYSTAFSON D E, KESSEL W C. Fuzzy clustering with a fuzzy covariance matrix[C]// Proc of IEEE Conference on Decision and Control.California:[s.n.], 1979:761-766.
[13]李士勇.工程模糊數(shù)學(xué)及應(yīng)用[M].哈爾濱:哈爾濱工業(yè)大學(xué)出版社, 2004:25-136.
注:“本文中所涉及到的圖表、注解、公式等內(nèi)容請以PDF格式閱讀原文。”