高光譜在冬蟲夏草含量及真偽鑒別中的應用

高光譜在冬蟲夏草含量及真偽鑒別中的應用

摘要：目的 以冬蟲夏草粉末為研究對象
，利用高光譜成像技術建立對冬蟲夏草粉末的真假鑒別及含量判斷的無損檢測模型。方法 將真偽樣品粉末分別按9：1、8：2、7：3
、6：4、5：5等不同比例進行混合
，通過光譜範圍為940 nm~2500 nm的高光譜成像儀分別獲取真偽樣品、不同比例真偽混合樣品的高光譜影像數據，經過主成分分析（PCA）初步判別真偽樣品的差異性，再結合偏最小二乘法（PLS）對真偽粉末樣品的含量進行分析判斷。結果 不同比例真偽混合物中冬蟲夏草粉末的識別度為97.0%-98.78%，偽冬蟲夏草粉末的識別度為83.10%-99.3%。 結論 基於高光成像技術可以實現對冬蟲夏草粉的真假辨別
，並可準確地判別出冬蟲夏草粉末的有效含量。

關鍵詞： 高光譜成像技術；冬蟲夏草粉末；真假鑒別；含量判斷

冬蟲夏草最初載於《本草從新》[1]，主產於我國四川、雲南、青海、西藏、甘肅等省區，是我國名貴中藥材之一。與人參、鹿茸並稱為補品“三寶”

，具ju有you良liang好hao的de醫yi療liao保bao健jian作zuo用yong。由you於yu冬dong蟲chong夏xia草cao人ren工gong撫fu育yu技ji術shu尚shang未wei突tu破po，其qi來lai源yuan僅jin靠kao野ye生sheng采cai挖wa
，產chan量liang較jiao少shao
，價jia格ge昂ang貴gui，少shao數shu不bu法fa分fen子zi為wei謀mou取qu暴bao利li而er采cai取qu以yi次ci充chong好hao

、chanjiachongzhendengshouduan
，daozhishichangshangdongchongxiacaodegezhongweipincengchubuqiong，yanzhongweihaixiaofeizhehehuanzhedeshentijiankang。yinci，wanshanchongcaodepinzhongjiandingyanjiu，zhunquedipandingqilaiyuan
，duiyaopindejianduhejianyanxiandeyouweibiyao。

本文運用近紅外高光譜成像技術對冬蟲夏草粉末的真假鑒別及含量進行無損檢測的判斷。與傳統的中藥檢測技術如色譜法、質譜法和光譜法[ 2-4] dengfangfaxiangbi
，gaoguangpuchengxiangjishubaerweichengxiangheguangpujishurongweiyiti，qigaoguangpushujubaohanguangpuhetuxiangxinxi，keyitongshibiaozhengbeiceduixiangdewaibutezhengheneibuxinxi，youshizaiyucaijidaodetuxiangxinxiliangfengfu，shibiedujiaogaoheshujumiaoshumoxingduo。youyuwutidefansheguangpujuyou“指紋”效xiao應ying，不bu同tong物wu不bu同tong譜pu，同tong物wu一yi定ding同tong譜pu的de原yuan理li來lai分fen辨bian不bu同tong的de物wu質zhi信xin息xi
，它ta不bu追zhui求qiu單dan一yi成cheng分fen的de控kong製zhi，其qi整zheng體ti性xing和he模mo糊hu性xing可ke以yi提ti供gong豐feng富fu的de中zhong藥yao信xin息xi，能neng夠gou更geng加jia有you效xiao地di體ti現xian中zhong藥yao成cheng分fen的de綜zong合he作zuo用yong，從cong而er更geng好hao地di鑒jian別bie中zhong藥yao真zhen偽wei
，評ping價jia中zhong藥yao質zhi量liang。高gao光guang譜pu技ji術shu近jin年nian來lai逐zhu漸jian受shou到dao生sheng物wu醫yi學xue
、精細農業

、食品安全等許多領域的重視。例如，外在品質檢測如水果
、蔬菜表麵損傷
、淤痕[5-6]
；內部品質檢測如水果的可溶性固體含量、水分含量、堅硬程度[7-8]，豬肉的嫩度[9-10]

，鱈魚的新鮮程度[11]
，菠菜葉片硝酸鹽含量等[12]
；食品安全檢測主要指食品中是否含有可能損害或威脅人體健康的物質，如蘋果
、哈密瓜表麵排泄物汙染檢測[13-15]，雞肉排泄物汙染檢測[16]，玉米、辣椒等黃曲黴毒素檢測等。

本研究利用高光譜成像儀“圖譜合一”的特點
，分析冬蟲夏草粉末、蛹蟲草粉末、地蠶粉末的光譜差異
，以期準確地判別出混合粉末中冬蟲夏草粉末的有效含量。

1  儀器與試劑

高光譜成像數據采集采用四川雙利合譜科技有限公司的GaiaSorter高光譜分選儀係統。該係統主要由高光譜成像儀

，CCD相機
、光源、暗箱
、平移台

、計算機組成。

蟲草粉樣品是由青海唐古拉藥業公司提供，其中包括3個標準樣品
，2個偽品（分別為地蠶
、蛹蟲草）。

2  實驗方法

2.1 樣品的采集與處理

  將青海唐古拉藥業公司提供的樣品粉末各取5g放置於培養皿上，標號
，用於高光譜相機的光譜采集。其中偽品2個標號為9號
、10號（分別為地蠶
、蛹蟲草），3個標準樣品6號

、7號、8號
，同時將8號與10號樣品粉末分別以9：1
、8：2、7：3、6：4、5：5等不同比例進行混合，得到1~5號樣品。

 2.2 高光譜數據的采集

樣品的采集高光譜成像係統
，掃描方式為推掃式成像技術，該儀器的光譜掃描範圍為940 nm ~2500 nm，每個樣品的采樣間隔為2.5 nm，每個樣品的測定速度一般小於 1分鍾。在獲取樣品的高光譜影像後，需要對采集的光譜圖像進行圖像校正才能得到樣品的反射率，圖像校正公式如下：

式中，R_ref 是校正過的圖像

，DN_raw是原始圖像，DN_white為白板校正圖像，DN_dark是黑板校正圖像。

樣品的典型高光譜如圖1所示。紅色為8號樣品、藍色和綠色分別為9號、10號樣品。

圖1  8號、9號、10號樣品典型光譜

Fig.1  Reflectance spectra of the 8號、9號
、10號 samples

3結果與分析

3.1 主成份分析

由於主成分圖像都是由原始數據中的各個波段下的圖像經過線性組合而成，根據 （其中

，為第m個主成分，為該主成分的權重係數，為單個波段的原始圖像）。比較該線性組合的權重係數，ruguoquanzhongxishujueduizhiyueda，duizhuchengfentuxianggongxianjiuyueda。zuidaquanzhongxishusuoduiyingdebochangxiadetuxiangweizuijiatezhengbochangtuxiang。weilexuanquzuijiadebochangzuhe

，bixubaozhengtamendequanzhongxishujueduizhijinkenengda
，tongshihaiyaobaozhengtamenzhijianyouyidingdebochangjiange。jiangceshiyuanshishujujinxingjiangzaochulihou，tongguozhuchengfenfenxibianhuan（PCA）進行背景扣除後，再次進行主成分分析變換（PCA）後結果如圖2所示。

  圖2 主成分分析的結果

  Fig. 2  The results of principal component analysis

通過主成分分析

，可以明顯將9號
、10號兩個偽樣品與其它樣品進行區分，由此可以判斷9號與10號樣品為偽樣品。

3.2 偏最小二乘法

偏最小二乘法（PLS）是(shi)一(yi)種(zhong)數(shu)學(xue)優(you)化(hua)的(de)技(ji)術(shu)，它(ta)主(zhu)要(yao)是(shi)通(tong)過(guo)最(zui)小(xiao)化(hua)誤(wu)差(cha)的(de)平(ping)方(fang)和(he)來(lai)找(zhao)到(dao)一(yi)組(zu)數(shu)據(ju)的(de)最(zui)佳(jia)函(han)數(shu)匹(pi)配(pei)，然(ran)後(hou)用(yong)最(zui)簡(jian)的(de)方(fang)法(fa)求(qiu)得(de)一(yi)些(xie)絕(jue)對(dui)不(bu)可(ke)知(zhi)的(de)真(zhen)值(zhi)

，而(er)令(ling)誤(wu)差(cha)平(ping)方(fang)之(zhi)和(he)為(wei)最(zui)小(xiao)。偏(pian)最(zui)小(xiao)二(er)乘(cheng)法(fa)（PLS）就相當於將多元線性回歸分析、典型的相關分析以及主成分分析方法融合在一起的數學方法。偏最小二乘法（PLS）主(zhu)要(yao)是(shi)通(tong)過(guo)主(zhu)成(cheng)分(fen)分(fen)析(xi)法(fa)將(jiang)多(duo)為(wei)空(kong)間(jian)數(shu)據(ju)的(de)曲(qu)線(xian)壓(ya)縮(suo)到(dao)較(jiao)低(di)維(wei)的(de)空(kong)間(jian)數(shu)據(ju)上(shang)，使(shi)其(qi)原(yuan)曲(qu)線(xian)分(fen)解(jie)為(wei)多(duo)種(zhong)主(zhu)成(cheng)分(fen)分(fen)析(xi)曲(qu)線(xian)，而(er)不(bu)同(tong)的(de)曲(qu)線(xian)的(de)主(zhu)成(cheng)分(fen)分(fen)別(bie)代(dai)表(biao)不(bu)同(tong)的(de)主(zhu)分(fen)和(he)因(yin)素(su)間(jian)對(dui)曲(qu)線(xian)的(de)貢(gong)獻(xian)率(lv)，選(xuan)取(qu)貢(gong)獻(xian)率(lv)較(jiao)大(da)的(de)主(zhu)成(cheng)分(fen)，去(qu)除(chu)有(you)幹(gan)擾(rao)組(zu)分(fen)和(he)幹(gan)擾(rao)因(yin)素(su)的(de)主(zhu)成(cheng)分(fen)，僅(jin)僅(jin)將(jiang)貢(gong)獻(xian)率(lv)較(jiao)高(gao)的(de)主(zhu)成(cheng)分(fen)與(yu)質(zhi)量(liang)參(can)數(shu)進(jin)行(xing)回(hui)歸(gui)。

本研究將8號樣品與10號樣品的粉末分別按9：1、8：2、7：3、6：4、5：5等不同比例進行混合，得到1~5號樣品。將1～5號樣品作為未知樣品進行偏最小二乘法變換（PLS）實現分類判別

，判斷8號與10號樣品的混合比例
，結果如下表1-表5所示。

從表1到表5可知
，1號到5號樣品中8號樣品冬蟲夏草粉末的成分識別度在97.0%到98.78%之間；蛹蟲草粉末的成分識別度在83.10%到99.3%之間；未識別成分所占比例在0.57%到0.84%之間。

從總體上看，基於PCA變換後的高光譜影像，利用偏最小二乘法可準確地識別出肉眼無法分辨的冬蟲夏草和蛹蟲草粉末。其中當冬蟲夏草和蛹蟲草粉末進行9：1混合時，91.1%識別為冬蟲夏草粉末，8.31%識別為蛹蟲草粉末，未識別成分占0.59%；當冬蟲夏草和蛹蟲草粉末按8：2混合時，78.26%識別為冬蟲夏草粉末
，20.14%識別為蛹蟲草粉末，未識別成分占0.60%

；當冬蟲夏草和蛹蟲草粉末按7：3混合時

，71.47%識別為冬蟲夏草粉末

，28.80%識別為蛹蟲草粉末，未識別成分占0.73%；當冬蟲夏草和蛹蟲草粉末按6：4混合時
，59.04%識別為冬蟲夏草粉末，40.36%識別為蛹蟲草粉末，未識別成分占0.57%；當冬蟲夏草和蛹蟲草粉末按5：5混合時，51.50%識別為冬蟲夏草粉末，47.66%識別為蛹蟲草粉末，未識別成分占0.84%。

表1 1號分類判別結果

Tab.1 The result of No.1 Classification

表2 2號分類判別結果

Tab.2 The result of No.2 Classification

表3 3號分類判別結果

Tab.3 The result of No.3 Classification

表4 4號分類判別結果

Tab.4 The result of No.4 Classification

表5 5號分類判別結果

Tab.5 The result of No.5 Classification

4 討論

tongguoyishangduidongchongxiacaofenmoyangpinzaijinhongwaiboduandefanshegaoguangputuxiangcaiji
，jingguozhuchengfenfenxi
，keyouxiaodeduidongchongxiacaofenmodezhenweijinxingchubudejianbie。jinyibutongguopianzuixiaoerchengfafenxiduiyangpinquyujinxingfenxipanduan
，keduiyangpindeyouxiaochengfenhanliangjinxingjianbie，qizhongdongchongxiacaofenmodechengfenshibieduzai97.0%到98.78%之間；蛹蟲草粉末的成分識別度在83.10%到99.3%之間。但由於目前樣品采樣數量較小，對於成分含量判斷的準確性還需進一步實驗驗證。

實(shi)驗(yan)初(chu)步(bu)驗(yan)證(zheng)了(le)高(gao)光(guang)譜(pu)成(cheng)像(xiang)技(ji)術(shu)在(zai)蟲(chong)草(cao)粉(fen)末(mo)鑒(jian)別(bie)的(de)可(ke)行(xing)性(xing)。進(jin)一(yi)步(bu)還(hai)需(xu)通(tong)過(guo)實(shi)驗(yan)和(he)分(fen)析(xi)判(pan)斷(duan)出(chu)理(li)想(xiang)的(de)特(te)征(zheng)波(bo)段(duan)
，以(yi)降(jiang)低(di)數(shu)據(ju)采(cai)集(ji)量(liang)，並(bing)進(jin)一(yi)步(bu)優(you)化(hua)數(shu)據(ju)分(fen)析(xi)模(mo)型(xing)與(yu)數(shu)據(ju)處(chu)理(li)速(su)度(du)

，從(cong)而(er)達(da)到(dao)在(zai)線(xian)檢(jian)測(ce)的(de)速(su)度(du)與(yu)準(zhun)確(que)性(xing)要(yao)求(qiu)。

REFERENCES

1. 張白翟,愛華.中藥冬蟲夏草的真偽鑒別[J]. 學員天地,2005,12(3):58.
2. 聶黎行,王鋼力,李誌猛等. 近紅外光譜法在中藥生產過程分析中的應用[J]. 光學學報, 2009, 29( 2): 541-547.
3. 毛佩芝,楊凱,金葉,劉雪鬆,王龍虎.近紅外光譜法快速測定野菊花藥材中水分及蒙花苷含量[J]. 中國現代應用藥學. 2015，32(12)：1428-1432
4. 白雁,劉建營,雷敬衛,謝彩俠,陳誌紅,龔海燕. 近紅外光譜技術快速測定杞菊地黃丸的水分含量[J]. 中國現代應用藥學. 2013，30(01):52-54
5. NICOLA B,LTZE E,PEIRS A. Non-destructive measurement of bitter pit in apple fruit using NIR hyperspectral imaging[J]. Postharvest Biology and Technology, 2006, 40: 1-6.
6. XING J,BRAVO C,JANCSÓK P T,et al. Detecting bruises on “golden delicious”apples using hyperspectral imaging with multiple wavebands[J].Biosystems Engineering, 2005, 90: 27-36.
7. RAJKUMAR P, WANG N, EIMASRY G. Studies on banana fruit quality and maturity stages using hyperspectral imaging[J]. Journal of Food Engineering, 2012,108: 194-200.
8. 吳建虎,彭彥昆,陳菁菁,等. 基於高光譜散射特征的牛肉品質參數的預測研究[J]. 光譜學與光譜分析, 2010, 30( 7) :1815-1 819．
9. WU J H,PENG Y K,LI Y Y,et al. Prediction of beef quality attributes using vis/nir hyperspectral scattering imaging technique[J]．Journal of Food Engineering,2012,109( 2):267- 273．
10. SIVERTSEN A H,KIMIYA T,HEIA K. Automatic freshness assessment of cod(gadus morhua) fillets byvis/nir spectroscopy [J]. Journal of Food Engineering, 2011, 103: 317-323.
11. 薛利紅,楊林章. 基於可見近紅外高光譜的菠菜硝酸鹽快速無損測定研究[J]. 光譜學與光譜分析, 2009, 29(4):926-930.
12. VARGAS A M,KIM M S,TAO Y,et al. Detection of fecal contamination on cantaloupes using hyperspectral fluorescence imagery[J]. Journal of Food Science, 2005,70: E471-E476.
13. Lefcout A M, Kim M S. Technique for normalizing intensity histograms of images when approximate size of the target is known: detection of feces on apples using fluorescence imaging[J]. Computers and Electronics in Agriculture, 2006, 50:135-147.
14. LEFCOUT A M, KIM M S. Technique for normalizing intensity histograms of images when approximate size of the target is known: detection of feces on apples using fluorescence imaging [J]. Computers and Electronics in Agriculture,2006,50:135-147.
15. PARK B,WINDHAM W R,LAWRENCE K C, et al. Contaminant classification of poultry hyperspectral imagery using a spectral angle mapper algorithm[J]. Biosystems Engineering, 2007, 96:323-333.
16. YAO H, HRUSKA Z,BROWN R L,et al. Hyperspectral bright greenish - yellow fluorescence(bgyf) imaging of aflatoxin contaminated corn kernels[J].  Proceedings of the SPIE, 2006, 6381: 63810B1-63810B8.