陈舒迪 柴琴琴 张勋 黄泽豪 林羽 徐伟
摘要: 金線莲是中国珍稀中草药,不同品系的金线莲具有细微的形态差异和显著的药效差异。针对金线莲的单一特征贡献能力不足以及传统分类器泛化能力不佳的问题,提出使用形状、颜色和纹理特征对金线莲叶片图像进行特征提取与融合,再使用表现性能更优的LightGBM(轻量级梯度提升机)构建分类器,以提高金线莲识别正确率。LightGBM具有精确高效等优点,将提取得到的高层次特征导入LightGBM进行训练预测,可以有效提高分类准确性。对金线莲数据集中的6个品系共368幅叶片图像进行试验,结果表明,相比于传统的分类方法,基于多特征融合和LightGBM的模型识别效果最好,10次随机试验的平均识别率比传统方法KNN、SVM和GBDT高,并且在分类评价指标精确率、召回率、综合评价指标上有较优表现,该研究结果可为中药材品系识别提供参考。
关键词: 金线莲;多特征融合;LightGBM;叶片识别
中图分类号: TP391 文献标识码: A 文章编号: 1000-4440(2021)01-0155-08
Identification of Anoectochilus roxburghii strains based on multi feature fusion and LightGBM
CHEN Shu-di1,2, CHAI Qin-qin1,2, ZHANG Xun3, HUANG Ze-hao3, LIN Yu3, XU Wei3
(1.College of Electrical Engineering and Automation, Fuzhou University, Fuzhou 350108, China;2.Ministry of Education Key Laboratory of Medical Instrument and Pharmaceutical Technology, Fuzhou University, Fuzhou 350108, China;3.College of Pharmacy, Fujian University of Traditional Chinese Medicine, Fuzhou 350122, China)
Abstract: Anoectochilus roxburghii is a rare Chinese herbal medicine, different strains of A.roxburghii have slight morphological difference and significant variance in medicinal effects. To solve the problems of insufficient contribution ability of single-feature in A.roxburghii and poor generalization ability of traditional classifiers, it was proposed to use shape, color and texture features to extract and fuse features of A.roxburghii leaf images. Then LightGBM (light weight class elevator) with a better performance was used in building classifier, so as to improve the recognition accuracy of A.roxburghii. LightGBM had the advantages of accurate and efficient, and the prediction accuracy could be improved effectively by importing the extracted high-level features into LightGBM to forecast the training. A total of 368 leaf images from six strains of A.roxburghii dataset were trained and tested. The results showed that model based on multi-feature fusion and LightGBM had the best recognition effect compared with traditional classification methods. The average recognition rate of ten random experiments was higher than traditional methods such as KNN, SVM and GBDT, and it showed good performance in classification evaluation indices like precision, recall rate and comprehensive evaluation index. The result can provide reference for the identification of different strains of traditional Chinese medicine.
Key words: Anoectochilus roxburghii;multi feature fusion;LightGBM;leaf recognition
金线莲属于兰科开唇兰属植物[1-2],含糖类、多种氨基酸及无机元素等,常用于糖尿病、肺热咳嗽、急慢性肝炎等疾病的治疗[3-4]。由于金线莲药用价值高、市场需求量大,而野生金线莲的繁殖率低、生存条件受限等原因,导致目前市场上销售的金线莲基本是人工栽培的。不同品系的金线莲虽外观相似但生物产量积累和药用化合物组成有着很大的差异,市面上以低药用价值的金线莲冒充优质品系的金线莲的现象尤为常见。因此,金线莲品系识别是选择优质品系的重要依据,如何精确高效地对金线莲品系进行识别有着重要的意义。
常用的金线莲品系鉴别方法有中药鉴定法、解剖学法和近红外光谱检测法等[5-8],这些传统方法存在主观性强、效率低、正确率低和操作复杂等缺陷,而基于机器视觉[9]的植物品系识别技术具有速度快、准确率高和操作方便等优点,因此使用计算机进行植物品系识别受到了越来越多的关注。机器视觉可获得不同品系金线莲植株的颜色[10]、形状[11]和纹理[12]等遗传信息,其中叶片数据最易采集且具有较长存活期,故以金线莲叶片作为研究对象,提取其颜色、形状和纹理特征,并将3种特征加以融合以进行识别。然而,不同品系的金线莲叶片图像十分相似,识别难度大,因此如何基于获得的图像特征,构建合适的分类器对采集的叶片图像进行进一步品系识别仍是具有挑战性的问题。
在分类器的选择上,传统的机器算法,如支持向量机(SVM)[13]、最近邻分类器(KNN)[14]、决策树(DT)[15]及逻辑回归(LR)[16]等在植物图像识别上应用广泛,但上述分类器的应用仍具有局限性,由于金线莲脉纹复杂、种类繁多,这些传统的分类方法在金线莲数据集上的表现也差强人意。深度学习[17]算法具有覆蓋范围广、学习能力强和特征提取方便等特点,Dyrmann等[18]基于Keras学习框架使用卷积神经网络对10 413幅旱田杂草图像进行识别,取得了较好的识别效果;邓向武等[19]将深度置信网络应用于稻田杂草图像的检测和分类,识别正确率为91.13%。然而直接通过原图像构建深度学习模型通常需要大量的标记样本,且训练时间太长;在处理小样本数据集时深度学习算法往往需要对数据集进行扩充,易出现精度低、过拟合等问题。考虑到本研究样本数量、分类品系不多和普通计算机硬件资源有限等问题,本研究提出用轻量级梯度提升机(Light gradient boosting machine,LightGBM)[20]作为分类模型。LightGBM是一种基于梯度提升树(Gradient boosting decision tree,GBDT)的提升算法,比起传统的 GBDT算法,LightGBM通过直方图优化对连续特征值进行分段处理,可以节省内存并加快运行速度,决策树通过Leaf-wise的叶子生长策略进行生长,限制树的深度防止过拟合,能够有效提高模型预测的准确率。Wang等[21]将LightGBM用于乳腺肿瘤的分类和诊断,并使用模型RF、XGBoost和LightGBM分别对样本进行分类识别,试验结果表明,LightGBM模型的分类精度最高。Zhang等[22]将LightGBM用于化学毒性预测,在所有Tox21(21世纪毒理学测试)和致突变性数据集上,与其他分类算法相比,LightGBM提供了最佳的预测性能。目前,LightGBM分类算法鲜被用于植物图像识别中。
综上所述,针对不同品系金线莲叶片差异不明显、复杂生长环境下不同品系金线莲叶片形态多变等因素引起的品系鉴别困难,本研究提出了基于多特征融合和LightGBM的金线莲品系识别。对提取得到的特征进行融合,解决了金线莲叶片单一特征贡献能力不足的缺点,通过LightGBM分类器进行预测,解决了传统机器学习算法识别效果不佳的问题。
1 样本预处理及特征提取
1.1 图像预处理
在试验研究阶段,所用金线莲样品均采自福建葛园生物科技有限公司,为金线莲成熟期叶片,采集6个品系的金线莲叶片,共368个样本,在实验室环境下用尼康单反数码相机(型号为NIKON D7100)获取金线莲叶片图像。其中台湾金线莲图像69幅、红霞图像70幅、小圆叶图像50幅、 尖叶图像58幅、一株圆叶图像66幅、 大圆叶图像55幅,采集的金线莲各品系代表图像如图1所示。
针对获取的金线莲叶片图像,在进行特征提取和识别之前先对其进行预处理,突出叶片图像中的有效特征,削弱噪声和虫洞的影响。将金线莲叶片图像转换成灰度图像,接着使用中值滤波对灰度图像进行去噪处理,增强叶片的目标信息,之后使用形态学膨胀和腐蚀方法对叶片虫洞进行填充,并将叶柄与叶面分割,最后选取自适应阈值确定的大津法获取叶片二值图像,将二值图像和金线莲叶片原始图像通过掩膜操作,获得金线莲叶片前景图像(图2)。
1.2 金线莲叶片特征提取
在人工智能时代,特征数据决定了机器学习的上限,特征数据是否有效是识别模型作用能否发挥到极致的重要前提[23]。如果从金线莲叶片图像中提取到的特征数据含有大量品系差异信息,那么识别模型的识别性能将会大大提高。本研究从金线莲叶片图像中提取23维形状特征、18维颜色特征和75维纹理特征,共116维特征,作为表征金线莲品系差异的主要特征。
1.2.1 形状特征 形状特征是图像识别中广泛应用的重要特征,由图1可以看出,不同品系的金线莲叶片存在一定的形态差异。本研究选择16维傅里叶特征以及7维HU不变矩特征共23维特征作为金线莲叶片的形状特征。
1.2.1.1 傅里叶描述子 傅里叶描述子只需通过少量的描述子即可大致代表整个金线莲叶片的轮廓[24]。将傅里叶描述子归一化后,提取到的形状特征可以不受轮廓在图像中的位置、大小和角度的影响,具有良好的鲁棒性。傅里叶描述子用1个总长度为N的离散的坐标序列{x(n),y(n)|n=0,1,…,N-1}来表示金线莲的闭合轮廓,坐标序列使用复数形式表示,如公式(1)所示:
通过傅里叶变换可得a(u):
对a(u)进行傅立叶反变换可获得金线莲叶片的边缘系数,本研究通过选取16个边缘系数来描述金线莲叶片轮廓。
1.2.1.2 HU不变矩 HU不变矩通过统计矩计算得到,常被用来表示图像的形状特征[25]。对于灰度分布为f(x,y)的图像,该图像的HU不变矩(p+q)阶统计量的定义为:
其中p,q=0,1,2,…,(p+q)阶中心矩定义为:
其中矩心为:x0=m10m00,y0=m01m00,通过公式(5)计算可得归一化中心矩:
式中:r=p+q+22,p,q=0,1,2,…,本研究通过二阶和三阶归一化中心矩共构造7个HU不变矩。
1.2.2 颜色特征 颜色特征是彩色图像中最显著的特征,常用于描述颜色特征的方法有颜色矩、颜色聚合向量和颜色集等[26]。本研究通过颜色矩来提取颜色特征,共提取了18维特征,分别是RGB(红-绿-蓝)、HSV(色调-饱和度-明度)颜色空间的R、G、B、H、S和V 6个分量的均值、方差和斜度,均值(Mi1)、方差(Mi2)和斜度(Mi3)的计算方法如公式(6)~公式(8)所示:
式中,Pij为图像中灰度为j的像素在第i个颜色通道分量中出现的概率,N为图像中所有像素的个数。
1.2.3 纹理特征 纹理特征描述了物体局部特征的周期性变化及排列规则,本研究根据金线莲叶片纹理特征分别使用灰度共生矩阵(GLCM)[27]和圆形LBP(CLBP)方法[28]对叶片图像进行纹理特征提取,共得到75维纹理特征。
1.2.3.1 灰度共生矩阵 灰度共生矩阵通过计算像素之间不同灰度值组合出现的次数来反映纹理信息。本研究通过计算金线莲叶片图像4个方向(0°,45°,90°,135°)的 4 个GLCM关键特征[对比度(Con)、能量(Ent)、熵(Asm)和相关性(Idm)],共可提取得到 16 个特征值。4个关键特征的计算方法如公式(9)~公式(12)所示:
式中P(ij)为归一化后的灰度共生矩阵,N为灰度级数。
1.2.3.2 圆形LBP LBP是一种常用的纹理描述方法,具有计算简单、速度快和旋转不变性等特点。然而传统LBP算子提取的是固定领域内的灰度值,当图像的比例发生改变时,传统LBP特征将无法正确地反映纹理信息。故本研究使用圆形LBP算子来描述纹理特征,圆形LBP可适应图像不同尺度的变换,使用灵活的圆形邻域代替固定的3×3正方形邻域,可设置得到采样点为P、半径为R的圆形LBP算子。对于1个点(Xc,Yc),它鄰近的(Xa,Ya)a∈P的计算公式如公式(13)、公式(14)所示:
式中,R为所设置的圆形领域的半径,P为采样点的个数。中心为(Xc,Yc)的LBP值计算方法如公式(15)、公式(16)所示:
本研究共设置8个采样点,共产生二进制形式的LBP值28=256种,由于太大的特征维度会影响特征融合后其他特征的描述效果,故使用Uniform降维方法[29]对采样得到的特征进行降维处理,特征由256维降至59维。
2 金线莲品类识别模型
2.1 梯度提升树(GBDT)原理
梯度提升(Gradient boosting)是一种常用的机器学习算法,常被用于回归和分类任务。梯度提升通过目前模型的残差来对新的弱分类器进行训练,之后将训练好的弱分类器不断累加到现有模型中,保证损失函数不断减小。假设fi(x)为弱分类器,得到最终强分类器为:
决策树(Decision tree)为树形结构的分类方法。从决策树的根节点开始进行决策,其中的每个内部节点代表1个划分标准,每个分支代表1个划分结果。
GBDT将梯度提升思想和决策树思想相融合,以决策树为弱分类器,通过多轮迭代,不断减小损失。在迭代过程中,逐步增加弱分类器,最后将训练好的各个弱分类器的结果进行集成,从而得出最终的预测结果,GBDT训练过程如图3所示。
2.2 LightGBM
2.2.1 直方图优化 LightGBM是GBDT算法的一种改进,LightGBM没有使用传统的预排序思想,而是对特征值进行直方图优化。直方图优化首先对特征数据进行“装箱处理”,得到n个“桶”。之后对数据进行遍历,将离散化后的值作为索引在直方图中的每个“桶”中进行累计,遍历完成后,就可以得到每个“桶”的统计量,直方图统计过程如图4所示。比起传统的预排序思想,直方图优化只需保存特征离散化后的值,明显降低了内存的损耗,加快了训练速度。
2.2.2 叶子生长策略为了提高模型训练的效率,LightGBM选择Leaf-wise作为叶子生长策略,传统的GBDT
算法通常使用Level-wise叶子生长策略,Leaf-wise和Level-wise的对比如图5所示。由图5可见,Level-wise叶子生长策略每次分裂需要对同一层的所有叶子同时进行分裂,因此在很多叶子上进行了不必要的分裂和搜索,而Leaf-wise每次分裂只对分裂增益最大的1个叶子进行分裂,避免不必要的分裂和搜索,可以达到提升分类精度的目的。但Leaf-wise可能会产生1棵过深的决策树,使结果过拟合,因此LightGBM在Leaf-wise叶子分裂策略的基础上增加了决策树最大深度的限制,最大程度避免了过拟合。
2.3 识别模型
为训练LightGBM识别模型,本研究选取70%的金线莲叶片特征数据用作训练样本,30%的特征数据用作测试样本,识别模型如图6所示。
识别流程如下:
(1)特征提取及融合:对6个品系共368幅金线莲叶片图像分别提取形状、颜色和纹理特征,通过融合,共提取得到116维特征。
(2)数据集划分:本研究从368个数据样本中选择258个数据样本作为训练样本,将剩下的110个数据作为测试样本。
(3)训练过程:
(a)初始化决策树棵数为m,训练样例权重为1/m;初始化迭代次数iter=1,设置最大迭代次数imax;
(b)根据样例特征直方图计算各个切分点,并寻找最优分割点;再采用带有深度限制的Leaf-wise叶子生长策略训练每棵决策树fi(x);
(c)根据误差确定当前弱分类器的话语权ai;
(d)按公式(17)合并基分類器为强分类器;
(e)若训练误差或者最大迭代次数满足要求则停止训练,转步骤(4); 否则iter=iter+1, 转步骤(b);
(4)预测过程:使用训练后的LightGBM识别模型对测试集进行预测,并获得识别结果。
2.4 试验结果及分析
本试验均在pycharm环境下进行,试验所用电脑为Windows 10系统,配备4 GB内存,处理器为英特尔第七代酷睿i5-7200U@2.50 GHz。为了研究基于多特征融合和LightGBM的识别模型的识别性能,本研究基于融合特征分别构建KNN、SVM、GBDT和LightGBM识别模型,统计金线莲品系识别正确率。并且为了避免识别结果的不确定性,数据集划分的随机数种子是由系统随机给出的,重复测试10次,4种分类器的10次试验分类结果如图7所示。
考虑到寻优算法在寻优过程中所需要的时间和模型的自适应等原因,本研究所用的4种识别模型KNN、SVM、GBDT和LightGBM在试验过程中均未对参数进行优化,4种识别模型均选择合理参数。从图7可知:将3种特征进行融合后,LightGBM分类器10次试验的平均正确率高于其他分类器,由于SVM分类模型对于小样本数据的识别效果比较好,故通过网格寻优算法对SVM模型进行参数寻优,10次随机试验平均正确率为88.20%,正确率仍然低于未寻优的LightGBM模型,体现了LightGBM识别模型良好的性能。为验证特征融合的有效性,基于颜色、形状和纹理单一特征分别构建KNN、SVM、GBDT和LightGBM识别模型,统计不同分类器和不同特征的金线莲品系识别率,随机进行10次试验,取平均值,结果如表1所示。
由表1可以看出,在金线莲品系识别中,叶片颜色特征对识别率贡献最大,将所有特征进行融合后识别率有很大的提升,识别率高于任意单一特征向量,证明了对特征进行融合的有效性。除基于单一颜色特征的LightGBM识别模型的识别率略低于GBDT识别模型外,基于单一形状特征、单一纹理特征和融合特征的LightGBM识别模型的识别率均高于GBDT、KNN和SVM识别模型,充分体现了LightGBM识别模型性能的优越性。
此外,本研究使用常见识别性能度量指标精确率、召回率和综合评价指标来评价识别模型,LightGBM对各个金线莲品系的识别性能度量指标如表2所示。
从表2数据可得,除红霞由于与尖叶的叶片特征相似度太高,精确率只有80.77%以外,其他各个品系的精确率、召回率、综合评价指标均在85.00%以上,可见基于多特征融合和LightGBM的模型在分类评价指标方面有较优表现。
综上可得,本研究提出的方法对金线莲叶片特征的表述更加准确,对金线莲数据集的识别率达91.08%,取得了较好的识别效果。
3 结论
本研究从金线莲叶片图像的单一特征贡献能力不足以及传统分类器泛化能力不佳引起的金线莲品系识别问题出发,研究基于机器视觉的金线莲品系识别问题,提出了基于多特征融合和LightGBM的金线莲品系识别方法。对金线莲叶片图像进行预处理后,分别从叶片颜色、形状和纹理3个方面共提取116维特征,并对特征进行融合,构建了基于LightGBM的分类模型,对金线莲进行品系识别,并与传统分类算法GBDT、KNN和SVM进行比较。结果表明,基于多特征融合和LightGBM的金线莲品系识别方法能够有效提取金线莲叶片特征,具有良好的鲁棒性,能够快速识别不同品系的金线莲。但本研究依然具有局限性,研究所用金线莲叶片图像数据集均为实验室环境下采集得来,叶片特征提取相对简单,如何处理复杂背景下的叶片图像,在复杂环境下对植物叶片进行识别,将是下一步研究的重点。
参考文献:
[1] 王 勇,陈 硕,卢端萍, 等.金线莲化学成分的研究[J]. 中草药, 2017(13):36-41.
[2] 何春年,王春兰,郭顺星,等. 福建金线莲的化学成分研究Ⅱ[J]. 中国中药杂志, 2005, 39(10): 581-582.
[3] 马玉芳,郑小香,衣伟萌,等. 金线莲多糖对免疫抑制小鼠脾淋巴细胞体外增殖、 分泌NO及细胞因子的影响[J]. 天然产物研究与开发, 2018,30(1):21-26.
[4] 李 芹,周 文,刘 路,等. 金线莲喷雾剂治疗手足口病口腔疱疹临床观察[J]. 福建中医药, 2012, 43(3):9-10.
[5] 肖雪洋. 植物叶片图像识别特征的研究和在线识别系统实现[D]. 合肥:中国科学技术大学,2011.
[6] 蔡 超. 金线莲的鉴别与应用[J]. 中国中医药现代远程教育, 2016, 14(21):123-124, 127.
[7] 蔡文燕,肖华山,范秀珍. 金线莲研究进展(综述)[J]. 亚热带植物科学, 2003(3):69-73.
[8] 王海阁,许 文,张 勋,等. 林下栽培金线莲的生药鉴别[J]. 中药材, 2020, 43(2):303-308.
[9] TAO O, LIN Z, ZHANG X B, et al. Research on identification model of chinese herbal medicine by texture feature parameter of transverse section image[J]. World Science and Technology-Modernization of Traditional Chinese Medicine, 2014(12):2558-2562.
[10]沈宝国,陈树人,尹建军,等. 基于颜色特征的棉田绿色杂草图像识别方法[J]. 农业工程学报,2009,25(6):163-167.
[11]HERRERA P J, DORADO J, RIBEIRO . A novel approach for weed type classification based on shape descriptors and a fuzzy decision-making method[J]. Sensors, 2014(14):15304-15324.
[12]TURKOGLU M, HANBAY D. Leaf-based plant species recognition based on improved local binary pattern and extreme learning machine[J]. Physica A: Statistical Mechanics and its Applications, 2019, 527: 121297.
[13]PAL M , FOODY G M . Feature selection for classification of hyperspectral data by SVM[J]. IEEE Transactions on Geoence and Remote Sensing, 2010, 48(5): 2297-2307.
[14]ZHANG S C, LI X L, ZONG M,et al. Learning k for kNN classification[J]. ACM Transactions on Intelligent Systems and Technology, 2017, 8(8):1-19.
[15]KIM S Y , UPNEJA A . Predicting restaurant financial distress using decision tree and AdaBoosted decision tree models[J]. Economic Modelling, 2014, 36(1): 354-362.
[16]FAGERLAND M W , HOSMER D W , BOFIN A M . Multinomial goodness-of-fit tests for logistic regression models[J]. Statistics in Medicine, 2008, 27(21): 4238-4253.
[17]HAN J , ZHANG D , CHENG G , et al. Advanced deep-learning techniques for salient and category-specific object detection: a survey[J]. IEEE Signal Processing Magazine, 2018, 35(1):84-100.
[18]DYRMANN M, KARSTOFT H, MIDTIBY H S, et al. Plant species classification using deep convolutional neural network[J]. Biosystems Engineering,2016, 151:72-80 .
[19]鄧向武,齐 龙,马 旭,等. 基于多特征融合和深度置信网络的稻田苗期杂草识别[J]. 农业工程学报, 2018, 34(14):173-180.
[20]KE G, MENG Q, FINLEY T, et al. Lightgbm: a highly efficient gradient boosting decision tree[C]//Ulrike von Luxburg, Isabelle Guyon. Proceedings of the 31st International Conference on Neural Information Processing SystemsDecember. San Diego,CA: Neural Information Processing Systems Foundation, 2017: 3146-3154.
[21]WANG D H, ZHANG Y, ZHAO Y. LightGBM: an effective miRNA classification method in breast cancer patients [C]// Association for Computing Machinery.Proceedings of the 2017 International Conference on Computational Biology and Bioinformatics. New York: Association for Computing Machinery, 2017: 7-11.
[22]ZHANG J, MUCS D, NORINDER U, et al. LightGBM: an effective and scalable algorithm for prediction of chemical toxicity-application to the Tox21 and mutagenicity data sets[J]. Journal of Chemical Information & Modeling, 2019, 59(10):4150-4158.
[23]POST F H, VROLIJK B, HAUSER H, et al. The state of the art in flow visualisation: feature extraction and tracking[J]. Computer Graphics Forum, 2010, 22(4):775-792.
[24]胡维炜,张 武,刘连忠,等. 利用图像处理技术计算大豆叶片相对病斑面积[J]. 江苏农业学报,2016,32(4):774-779.
[25]WANG X F, HUANG D S, DU J X, et al. Classification of plant leaf images with complicated background[J]. Applied Mathematics and Computation, 2008, 205(2): 916-926.
[26]HWANG W, WANG H T, KIM H, et al. Face recognition system using multiple face model of hybrid fourier feature under uncontrolled illumination variation[J]. IEEE Transactions on Image Processing, 2011, 20(4):1152-1165.
[27]HARALICK R M, SHANMUGAM K, DINSTEIN I H. Textural features for image classification[J]. Studies in Media and Communication, 1973, 3(6): 610-621.
[28]LIU L, FIEGUTH P, GUO Y, et al. Local binary features for texture classification: taxonomy and experimental study[J]. Pattern Recognition, 2017(62): 135-160.
[29]OJALA T, PIETIKINEN M, MENP T. Multiresolution gray-scale and rotation invariant texture classification with local binary patterns[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2002, 24(7): 971-987.
(責任编辑:陈海霞)
收稿日期:2020-06-02
基金项目:福建省科技厅重大产学研项目(2019Y4009);晋江市福大科教园区发展中心科研项目(2019-JJFDKY-48)
作者简介:陈舒迪(1997-),女,福建南平人,硕士研究生,主要从事机器学习与图像处理研究。(E-mail)chshud@163.com
通讯作者:柴琴琴,(E-mail)qq.chai@fzu.edu.cn;林 羽,(E-mail)lyfjtcm@163.com