林业科学  2019, Vol. 55 Issue (1): 38-46   PDF    
DOI: 10.11707/j.1001-7488.20190105
0

文章信息

蔡小娜, 苏筱雨, 黄大庄, 沈佐锐.
Cai Xiaona, Su Xiaoyu, Huang Dazhuang, Shen Zuorui.
基于几何形态测量学的天蛾科成虫数字化分类
Digital Classification of Sphingid Moths Adults (Lepidoptera: Sphingidae) Based on Geometric Morphometry
林业科学, 2019, 55(1): 38-46.
Scientia Silvae Sinicae, 2019, 55(1): 38-46.
DOI: 10.11707/j.1001-7488.20190105

文章历史

收稿日期:2017-06-19
修回日期:2018-04-07

作者相关文章

蔡小娜
苏筱雨
黄大庄
沈佐锐

基于几何形态测量学的天蛾科成虫数字化分类
蔡小娜1, 苏筱雨2, 黄大庄2, 沈佐锐3     
1. 河北金融学院基础课教学部 保定 071051;
2. 河北农业大学林学院 保定 071000;
3. 中国农业大学 IPMist 实验室 北京 100193
摘要:【目的】对天蛾科10种成虫的前翅进行几何形态测量学分析,探讨利用几何形态测量学实现天蛾科成虫数字化分类的可行性,为逐步实现蛾类昆虫的数字化分类提供新的形态学依据。【方法】首先,进行几何形态测量学分析:以天蛾右前翅为研究对象,按特定次序选取17个翅脉交点作为标记点,并获得坐标数据;对原始标记点坐标数据进行普氏叠加分析消除标本摆放位置、方向和大小等非形状因素等信息;对普氏叠加后的标记点数据进行相对扭曲分析,得到17个标记点对10种天蛾分类作用的大小。其次,对几何形态测量学分析所得数据进行多元统计分析:利用单因素方差分析对17个标记点的差异显著性进行检验,再对普氏叠加后的标记点数据进行主成分分析,然后利用前3个主成分进行判别分析。【结果】相对扭曲分析表明标记点2、4和5对于10种天蛾成虫的分类作用相对较大;单因素方差分析显示17个标记点均具有显著的差异,即标记点在种间是有显著差异的,可以用于本文10种天蛾成虫的分类鉴定;主成分分析的前3个主成分的累计贡献率为97.7%,可作为对10种天蛾成虫进行分类鉴定的变量;判别分析结果显示回归判别和交叉判别的准确率均为100%,实现对天蛾科10种成虫的分类鉴定。【结论】研究表明几何形态测量学可应用于天蛾成虫的数字化分类鉴定,可为未来进一步实现蛾类成虫的自动识别奠定基础。
关键词:天蛾    前翅翅脉    几何形态测量学    分类鉴定    森林昆虫    
Digital Classification of Sphingid Moths Adults (Lepidoptera: Sphingidae) Based on Geometric Morphometry
Cai Xiaona1, Su Xiaoyu2, Huang Dazhuang2, Shen Zuorui3     
1. Department of Basic Courses, Hebei Finance University Baoding 071051;
2. College of Forestry, Hebei Agricultural University Baoding 071000;
3. IPMist Laboratory, China Agricultural University Beijing 100193
Abstract: 【Objective】The present paper analyzed geometric morphology of 10 species of sphingid adults using geometric morphometric method to explore the feasibility of geometric morphometry on digital classification of sphingid moths for providing new morphological foundation on realizing digital classification of moths progressively.【Method】In the first step, geometric morphometric analysis was carried out. Right forewings of 10 sphingid species were used as research objects and 17 intersections of wing veins were selected as landmarks according to specific order, and coordinate data were obtained. Procrustes superimposition was applied to data of the 17 landmarks to remove non-shape variation from the landmark coordinates such as placing location, direction and sizes; and further, relative warp analysis was used to the new landmark coordinates to obtain effects of 17 landmarks to the 10 sphingid species. In the second step, multivariate statistical analysis was applied to the data collected by geometric morphometric analysis. One-way variance analysis was used for significant test of difference of landmark coordinates firstly; then, principal component analysis of the Procrustes transformed data-set was implemented; at last, discriminant analysis was carried out on the first three components.【Result】The result showed that Landmarks 2, 4 and 5 had relative large contributions on the classification of the 10 sphingid adult species according to relative warp analysis. All of the 17 landmarks had significant difference and could be used for the classification of those 10 Sphingidae adult species by using one-way variance analysis. The accumulated rate of contribution of three principal components of the principal component analysis reached 97.7%, which could be used as variables of the classification and identification of the 10 sphingid adult species. The accuracy of original and cross validation tests reached to 100% and 99.7%, respectively. Classification and identification of the 10 sphingid moths were realized.【Conclusion】It was indicated that geometric morphometric analysis could be used to identify the sphingid adults accurately and would be useful in gradually realizing the automatic recognition of moths in the future.
Key words: Sphingidae species    forewing veins    geometric morphometry    identification    forest insect    

据文献记载,世界上已知隶属于鳞翅目(Lepidoptera)天蛾科(Sphingidae)的天蛾种类约有1 000种,我国天蛾种类约占世界种类的1/5(朱弘复等,1997)。天蛾科昆虫多为农林业重要害虫,对农林业造成了重要的危害和经济损失,准确地进行昆虫物种鉴定是害虫测报与防治的基础,同时动植物检疫和资源昆虫的保护及利用等方面也需要昆虫鉴定的工作。当前对蛾类昆虫鉴定还主要依靠传统的人工方法,这不仅要求鉴定人员具备一定的专业知识与长期的经验积累,而且鉴定工作耗时、主观性强。如何在生物分类学家人数减少的同时,能够对大量分类学资源进行快捷有效的信息提取和准确鉴定,是当今人们面临的一个重要挑战(Gaston et al., 2004张蕾等,2011)。

几何形态测量学(Geometric morphometric)是关于生物形态结构定量比较的科学计算工具(白明等,2014a),该方法通过在图形中设置具有生物学意义的同源标记点(Landmarks)或轮廓线而将图形引入二维或三维坐标系,并利用标记点的坐标值来捕获研究对象的形状信息(Slice, 2001Adams et al., 2004Zelditch et al., 2004Mitteroecker et al., 2009)。同时,通过标记点间的相对位置使样本间形状差异变得更加直观(Adams et al., 2013)。一直以来,关于几何形态测量学的报道主要集中在医学(Slice,2007)、大型哺乳动物(Fureix et al., 2011Baab et al., 2012)、古人类学(邢松等,2012张亚盟等,2016周蜜等,2016)和少数节肢动物(Gutierrez et al., 2011Bai et al., 2012Hübler et al., 2013白明等,2014b)等领域。近十年,随着几何形态测量学日臻完善和成熟,其越来越多地应用在昆虫学的相关领域,如:昆虫的形态进化、系统发育和亲缘关系(潘鹏亮等,2008闫宝荣,2011Ge et al., 2011Outomuro et al., 2013);昆虫的雌雄二态性(白义等,2015);种间或种下阶元昆虫的形态分化规律、地理环境因素对昆虫形态分化的影响(Hernández et al., 2010Kitthawee et al., 2011张秀秀等,2017)等。在国内,几何形态测量学方法在各个领域发挥着越来越重要的作用(詹庆斌,2014佟一杰等,2016李荣荣等,2016),但是应用于昆虫数字化分类的研究还不多见。

针对天蛾成虫翅型极其相似,利用数学形态特征较难实现天蛾成虫数字化分类识别的问题,本研究以天蛾科10种天蛾成虫作为研究对象,对天蛾成虫前翅进行几何形态测量学分析,探讨利用几何形态测量方法实现天蛾成虫的数字化分类。若利用几何形态测量学对10种天蛾可取得较好的分类效果,便可进一步将其应用在更多种类的昆虫分类鉴定中。本文利用天蛾成虫右前翅进行研究,而不是利用昆虫的整个标本特征,这是考虑到:蛾类翅的翅脉特征不仅是蛾类等多数昆虫的主要鉴定特征之一,且其翅的翅脉图像便于获取,同时,具有同源性的标记点容易选取,并且保证了标记点的可重复选择,便于日后添加其他种类的蛾类昆虫。如果研究证实了本文提出的方法是可行的,则可将其扩展应用到其他蛾类昆虫成虫的分类鉴定中;同时,如果研究表明仅以右前翅翅脉特征鉴定蛾类昆虫的成虫仍有不足,还可考虑利用蛾类成虫其他部位的特征,以提高鉴定的准确性。

1 材料与方法 1.1 试验材料

本文选择天蛾科10种天蛾成虫作为试验材料,以它们的右前翅为研究对象,对天蛾成虫前翅脉图像进行标记点的选取,并获得相应的坐标数据,探讨标记点数据应用于天蛾成虫分类鉴定的可行性。试验所选10种天蛾隶属于4个亚科,其中面形天蛾亚科(Acherontiinae)1种、蜂形天蛾亚科(Philampelinae)1种、斜纹天蛾亚科(Choerocampinae)3种、云纹天蛾亚科(Ambulicinae)5种。这10种天蛾均为较常见的物种,并且物种间前翅翅脉的拓扑结构很相似,而几何形态测量学侧重于对生物形态的拓扑结构进行研究,这可帮助笔者从技术角度探明几何形态测量学方法在蛾类昆虫分类鉴定中应用的可行性。试验所用10种天蛾,每个种类取30个样本。10种蛾具体信息见表 1

表 1 试验所用10种天蛾名录 Tab.1 Directory of 10 sphingid species for experiment
1.2 翅脉图像获取

由于鳞翅目成虫翅面上覆盖有大量鳞片,因此获取其翅脉图像首先要考虑去除蛾翅上的鳞片。本文利用化学浸泡法去除蛾翅面的鳞片和斑纹(蔡小娜,2013),所采用药品、仪器和具体步骤如下。

药品和仪器:盐酸(质量分数36%~38%),国药集团化学试剂有限公司;乙醇(质量分数95%),国药集团化学试剂有限公司;次氯酸钠(有效氯≥10%),国药集团化学试剂有限公司;大号培养皿(5套)、定性滤纸(中速)、水盆(2个)、镊子。

具体操作步骤:第1步,从蛾标本上取下完整的前翅,放入70%左右的乙醇溶液中约3 min或更长时间,使翅面充分浸润并具有一定的弹性,以免去鳞片不彻底或在以后的操作中折损蛾翅;第2步,用镊子把已浸润的蛾翅移入次氯酸钠溶液中1~3 min;第3步,再将翅移入盐酸溶液中约1~3 min,翅面的鳞片和斑纹会逐渐褪去。如果鳞片或斑纹没有去除干净,可重复第2和第3步,直到翅脉清晰可辨;第4步,将处理好的蛾翅移入水盆中,稀释掉残留的盐酸或次氯酸钠溶液;第5步,用纸条把翅从水中黏附出来,等翅面不再有流动的水时把其压在平整的纸层中,让其自然阴干,制成蛾翅脉标本。

蛾翅翅脉数字图像获取通过带有透扫功能的中晶Phantom v900 Plus型平板式扫描仪,在600 dpi分辨率下对翅脉标本进行负片扫描。图像保存为8位(8bit)的灰度图,在分析时使用原始图像,图像未压缩。

2 几何形态测量学 2.1 标记点选取

标记点的选择对于分析的结果有直接而且重要的影响,一般选择Bookstein分类(Bookstein,1991)中的第Ⅰ型或第Ⅱ型作为标记点的标准。这2种标记点的共同点是可以从解剖学特征上准确地在每一个标本上找到,而无不明确或复现困难的可能性。为了保证标记点的可重复选择,方便日后添加其他种类的蛾类昆虫,本研究选择第Ⅰ型标记点,即研究对象间其同源性有强烈证据支持的数学点(闫宝荣等,2010),譬如结构上的交点、昆虫翅脉的交点和端点。

利用tpsDig软件对天蛾右前翅翅脉交点按特定次序进行标记(Rohlf,1998a蔡小娜,2013)。翅面交点从径脉(R1)基部为起始点,顺时针方向依次取翅室周围所有翅脉交点,之后选取径脉R4和R5的交点和臀脉(A)基部分岔的交点,最后从径脉(R4)与翅外缘的交点开始顺时针方向选取所有翅外缘翅脉交点,共选择17个翅脉交点作为标记点,以深色白眉天蛾(Celerio gallii)为例给出所选取标记点的位置分布简图(图 1)。

图 1 深色白眉天蛾右前翅标记点分布 Fig. 1 Landmarks distribution from right forewing of Celerio gallii
2.2 普氏叠加

标记点的原始坐标数据包含了许多非形状因素,如获取图像时标本摆放的位置、方向、坐标系的选取(原点、刻度)及大小等(Rohlf,1999)。因此,利用软件tpsSuper(Rohlf,1998b)对样本数据进行普氏叠加分析,通过置中、旋转和缩放,以几何学上相应的线性代数作矩阵运算,所得运算结果消除了上述非形状因素等多余的信息,留下真正想要的形状。以深色白眉天蛾为例,进行普氏叠加分析可得到所有样本标记点的散点图和向量图(图 2),从图中可以看出通过普氏分析相应的标记点都排列在质心的周围,这也从可视化角度验证了叠加后的标记点分布情况。

图 2 深色白眉天蛾标记点普氏叠加 Fig. 2 Procrustes superimposition on landmarks of C. gallii a. 30个样本的标记点散点图Scatter plot for landmarks of 30 specimen;b.深色白眉天蛾30个样本的标记点向量图Vector plot for landmarks of 30 specimen.
2.3 相对扭曲

由相对扭曲分析得到17个标记点的贡献率如表 2所示,由表 2可知,标记点2、4和5的贡献率相对较大,分别为34.54%、22.86%和11.85%,这也说明10种天蛾昆虫前翅基部的变异大于端部,或者可以说翅室周围的变异相对较大。

表 2 深色白眉天蛾右前翅17个标记点的贡献率 Tab.2 Contribution of 17 landmarks from right forewing of C. gallii
3 多元统计分析与结果 3.1 标记点差异显著性检验

单因素方差分析公式为$F = \frac{{{{\overline S }_{\rm{A}}}}}{{{{\overline S }_{\rm{E}}}}} = \frac{{{S_{\rm{A}}}}}{{s - 1}}/\frac{{{S_{\rm{E}}}}}{{n - s}}$(盛骤等,2001),式中:SA为效应平方和,SE为误差平方和,n为样本总和,s为因素水平。

通过相对扭曲分析,可得到17个标记点的效应平方和SA与误差平方和SE(表 3)。由本文样本数据可知n=300,s=10。由单因素方差分析公式可计算得到各个标记点的F值(表 3)。根据F分布表,F0.05(9, 290)=1.88、F0.01(9, 290)=2.41。由此可知在显著性水平为0.05和0.01时F值均小于各个标记点的F值,所以10种天蛾的17个标记点间均具有显著甚至极显著的差异,即标记点在种间是有显著差异的,可以用于本文10种天蛾的分类鉴定。

表 3 10种天蛾翅脉标记点单因素方差分析结果 Tab.3 Results of one-way variance analysis of landmarks on veins of 10 sphingid species
3.2 主成分分析

利用软件SPSS13.0对10种天蛾300个样本经普氏叠加所得标记点坐标数据进行主成分分析,得到主成分的方差解释率(表 4)和主成分的得分系数矩阵。由表 4可知,第1个主成分的贡献率已达到86.01%,而一般认为若k个主成分达到85%,表明取前k个主成分基本包含了全部测量指标所具有的信息(于秀林等,1999)。

表 4 主成分的方差解释率 Tab.4 Variance explained of the principal components

根据表 4的方差累计百分比,本文确定提取前3个主成分(累计贡献率为97.7%),由主成分分析得到的前3个主成分的得分系数矩阵如表 5所示。根据该得分系数矩阵以及变量的观察值可计算出主成分得分,因通过表 4已确定提取3个主成分,此处将相应产生3个主成分得分,它们可作为对10种天蛾成虫进行分类鉴定的新变量。

表 5 前3个主成分的得分系数矩阵 Tab.5 Score coefficient matrix of the first three principal components
3.3 判别分析

利用软件SPSS13.0对主成分分析所得的前3个主成分得分矩阵进行判别分析,得到10种天蛾300个样本的二维散点图(图 3)。从图 3中可知,10种天蛾明显地分为了10组,并且每种天蛾都很好地聚在了一起。300个样本的回归判别和交叉判别结果如表 6所示,回归判别和交叉判别的正确率均为100%,即利用几何形态测量学数据对10种天蛾成虫进行分类鉴定,达到了百分之百的准确率。

图 3 10种天蛾的判别散点分布 Fig. 3 Discriminant scatter plot of 10 sphingid species 1:白薯天蛾H. convolvuli; 2:红天蛾P. elpenor lewisi; 3:黄脉天蛾A. amurensis; 4:蓝目天蛾S. planus planus; 5:葡萄天蛾A. rubiginosa rubiginosa; 6:雀纹天蛾T. japonica; 7:深色白眉天蛾C. gallii; 8:洋槐天蛾C. deucalio; 9:榆绿天蛾 C. tatarinovi; 10:枣桃六点天蛾 M. gaschkewitschi gaschkewitschi.
表 6 10种天蛾300个样本的判别分析结果 Tab.6 Discriminant analysis results from 300 specimen of 10 sphingid species
4 讨论

天蛾昆虫中有很多种类为害虫,以后可以利用几何形态测量学逐渐扩大对其种类的研究,实现常见天蛾害虫的数字化分类。但随着需要鉴别的天蛾种类的增加,可能需要在翅脉上选取更多标记点才能满足要求。此外,在今后的研究中,如果仅用前翅特征不能满足需求,还可将蛾后翅、头部、前胸背板、腹部、体躯等部位的数字化特征应用到昆虫分类中。昆虫自动分类鉴定的一个很重要的研究步骤就是对特征参数的选取,选取一组好的参数特征,既能方便提取和应用,又可使昆虫分类的准确率和有效性得以提高。随着自动分类系统中昆虫种类的不断增加,需要不断探索新的特征参数。

本文供试样本雌雄蛾均有,且是随机抽取个体,在分析过程中采用了与翅大小尺度和方向均无关的形态特征,因此虽然同一种蛾在雌雄之间会有一定差异,但是本文所用的这些特征对10种天蛾昆虫进行研究时并未受到影响。随着今后研究蛾类昆虫种类的不断增多,蛾翅在雌雄之间的差异会逐渐体现出来,这是需要注意与克服的问题。

由于天蛾科所属的鳞翅目其自身的特点,翅面上覆盖有大量鳞片,在生产实际中对于蛾翅的处理略显繁琐。为了便于处理试验所用的300个样本,本文采用了化学去鳞的方法且用扫描仪获取图像。在实际应用中,也可考虑采用毛笔刷去鳞片或漂白粉溶液使翅脱色等简便方法去除鳞片。将去鳞后的蛾翅脉标本放置在颜色较深的背景下,用普通相机即可获取蛾翅脉图像。相信随着获取图像的技术手段的不断改进,会有更简便的方法直接获取蛾翅脉图像,以便使昆虫翅脉特征更便捷地应用到蛾类昆虫的分类鉴定中。

本文利用几何形态测量学方法对10种天蛾进行数字化分类的准确率达到了100%,这说明几何形态测量学方法具有较高的分类准确率。但是,应用几何形态测量学方法对鳞翅目昆虫进行数字化分类也有其缺点:因鳞翅目昆虫翅面上覆盖有大量鳞片,研究翅脉的拓扑结构特征需要先去除翅面上的鳞片。在今后的研究中,可将数学形态学、几何形态测量学和颜色特征等多种数字化鉴定方法进行对比研究,明确这些方法的优缺点,同时积极探索新的方法,在保证较高准确率的前提下,更快、更简便地实现昆虫的数字化分类鉴定。

将几何形态测量学应用于昆虫系统发育研究已引起了大家的广泛兴趣(闫宝荣等,2010),今后笔者也将进一步探索利用几何形态测量学对天蛾甚至更多种类的蛾类昆虫进行系统发育的研究。随着对昆虫系统发育研究的各种技术方法不断发展,可以将此方法和分子系统学以及传统分类学放到一起进行探讨,并最终将几种方法结合起来推进昆虫分类学的发展。

几何形态测量学分析主要是以二维数据为主,从二维数据到三维数据的分析也已经可行,由于三维数据具有维度高、信息量大等特点(白明等,2014b),可帮助解决某些二维数据无法涉足的问题。在人类学研究和生物医学科学领域,三维地标点方法已成为公认的标准方法(Adamset al., 2004)。本研究中的蛾翅是平面结构,较易获取其二维图像,利用蛾翅二维标记点数据能够很好地实现10种天蛾的分类鉴定。三维数据的获取需要特殊且较昂贵的设备支持,数据处理也相对繁琐,但是随着计算机处理海量数据的能力不断提高以及研究昆虫种类的不断拓展,三维数据也将会被应用到昆虫的分类鉴定中,几何形态测量学方法可能会由二维数据为主的研究向三维的方向发展,并且三维数据将在分类鉴定、系统发育以及其他各个领域得到广泛应用。

几何形态测量学方法目前在医学、考古学、人类学和脊椎动物学等领域已有广泛应用,但由于昆虫个体相对较小,该方法在昆虫学领域的应用仍处于起步阶段。如今几何形态测量学方法通过与其他研究手段相结合,开始向其他研究领域渗透(白明等,2007),如非对称性结构的研究、内部解剖对称结构的任意横向距离的研究、不定型非对称性结构的研究。随着几何形态测量学方法的不断改进和成熟,加之很多学者开发出大量免费的分析软件以及相关研讨培训活动的开展,此方法将会得到越来越广泛的应用。

5 结论

本文利用天蛾成虫右前翅翅脉的标记点数据实现了对试验所用10种天蛾成虫的数字化分类鉴定,回归判别和交叉判别的正确率均为100%,这表明利用标记点法来实现天蛾成虫的分类鉴定是可行的,为未来逐步实现蛾类成虫的自动识别鉴定提供了简便有效的方法。

本文所用10种天蛾成虫前翅翅脉拓扑结构极其相似,利用几何形态测量学方法对其进行分类研究,准确率达到了100%,这充分证明了此方法提取的特征可应用于蛾类成虫的分类鉴定。此外,几何形态测量学方法中的相对扭曲分析还可实现昆虫分类鉴定的可视化,可视化作为几何形态测量学的特点之一,能够使不同物种之间的形态差异清楚可见。今后,笔者将继续利用几何形态测量学对蛾类其他昆虫进行分类研究,将实现昆虫分类鉴定的可视化作为研究方向。

本文利用几何形态测量学对天蛾成虫进行分类鉴定,关于标记点的选取与分析须注意如下2点:1)标记点的选择对于分析的结果有直接而且重要的影响,为了方便以后添加其他种类的蛾类昆虫,一般会选择容易复现的标记点,即要保证标记点的可重复选择。2)标记点的原始坐标数据由于包含了许多非形状因素,如获取图像时标本摆放的位置、方向、坐标系的选取(原点、刻度)等,因此不能直接对其进行数据分析,要经过普世叠加分析消除这些非形状因素等多余的信息后才能进行相应的分析。

参考文献(References)
白明, 杨星科. 2007. 几何形态测量法在生物形态学研究中的应用. 昆虫知识, 44(1): 143-147.
(Bai M, Yang X K. 2007. Application of geometric morphometrics in biological researches. Chinese Bulletin of Entomology, 44(1): 143-147. DOI:10.3969/j.issn.0452-8255.2007.01.035 [in Chinese])
白明, 杨星科, 李静, 等. 2014a. 几何形态学:关于形态定量比较的科学计算工具. 科学通报, 59(10): 887-894.
(Bai M, Yang X K, Li J, et al. 2014a. Geometric morphometrics, a super scientific computing tool in morphology comparison. Chin Sci Bull, 59(10): 887-894. [in Chinese])
白明, 杨星科. 2014b. 三维几何形态学概述及其在昆虫学中的应用. 昆虫学报, 57(9): 1105-1111.
(Bai M, Yang X K. 2014b. A review of three-dimensional (3D) geometric morphometrics and its application in entomology. Acta Entomologica Sinica, 57(9): 1105-1111. [in Chinese])
白义, 戴董峰, 包克瓯, 等. 2015. 基于几何形态测量学的黄蜻翅的雌雄二态性研究. 应用昆虫学报, 52(2): 363-369.
(Bai Y, Dai D F, Bao K O, et al. 2015. Using geometric morphometrics to quantify the sexual dimorphism of Pantala flavescens. Chinese Journal of Applied Entomology, 52(2): 363-369. [in Chinese])
蔡小娜. 2013.基于数字图像的主要蛾类害虫分类识别研究.保定: 河北农业大学博士学位论文.
(Cai X N. 2013. Classification research on moths based on digital images.Baoding: PhD thesis of Agricultural University of Hebei.[in Chinese]) http://cdmd.cnki.com.cn/Article/CDMD-10086-1014106210.htm
李荣荣, 李生才, 张虎芳. 2016. 几何形态测量学及其在半翅目中的研究进展. 山西农业大学学报:自然科学版, 36(4): 235-241.
(Li R R, Li S C, Zhang H F. 2016. Geometric morphometrics and its advances in Hemiptera. J Shanxi Agric Univ:Natural Science Edition, 36(4): 235-241. [in Chinese])
潘鹏亮, 杨红珍, 沈佐锐, 等. 2008. 翅脉的数学形态特征在蝴蝶分类鉴定中的应用研究. 昆虫分类学报, 30(2): 151-160.
(Pan P L, Yang H Z, Shen Z R, et al. 2008. Research on applying vein feature for mathematical morphology in classification and identification of butterflies (Lepidoptera:Rhopalocera). Entomotaxonomia, 30(2): 151-160. DOI:10.3969/j.issn.1000-7482.2008.02.013 [in Chinese])
盛骤, 谢式千, 潘承毅. 2001. 概率论与数理统计. 北京: 高等教育出版社.
(Sheng Z, Xie S Q, Pan C Y. 2001. Probability and mathematical statistics. Beijing: Higher Education Press. [in Chinese])
佟一杰, 杨海东, 马德英, 等. 2016. 蜣螂后胸叉骨的几何形态学分析及其适应进化研究. 昆虫学报, 59(8): 871-879.
(Tong Y J, Yang H D, Ma D Y, et al. 2016. Geometric morphometrics evaluation and adaptive evolution study of scarabaeine metendosternites. Acta Entomologica Sinica, 59(8): 871-879. [in Chinese])
邢松, 张银运, 刘武. 2012. 周口店直立人3号与5号头骨形态特征对比及其演化速率所反映的群体隔离. 人类学学报, 31(3): 250-258.
(Xing S, Zhang Y Y, Liu W. 2012. Morphological comparison of the ZKD 3 and 5 skulls and the probable population isolation as reflected by evolutionary rates. Acta Anthropologica Sinica, 31(3): 250-258. [in Chinese])
闫宝荣. 2011.基于几何形态测量学的蝎蛉科昆虫系统发育研究.杨凌: 西北农林科技大学硕士学位论文.
(Yan B R. 2011. Phylogeny of Panorpidae (Mecoptera) based on geometric morphometrics.Yangling: MS thesis of Northwest A & F University.[in Chinese]) http://cdmd.cnki.com.cn/Article/CDMD-10712-1011462413.htm
闫宝荣, 花保祯. 2010. 几何形态测量学及其在昆虫分类学和系统发育中的应用. 昆虫分类学报, 32(4): 313-320.
(Yan B R, Hua B Z. 2010. Geometric morphometrics and its application in the systematics and phylogenetics of insects. Entomotaxonomia, 32(4): 313-320. [in Chinese])
于秀林, 任雪松. 1999. 多元统计分析. 北京: 中国统计出版社.
(Yu X L, Ren X S. 1999. Multivariate statistical analysis. Beijing: China Statistics Press. [in Chinese])
詹庆斌. 2014.中国蚁蛉科翅几何形态学研究及其四属的分类修订.北京: 中国农业大学博士学位论文.
(Zhan Q B. 2014. Geometric morphometric analysis of wing shape in Chinese antlion and taxonomic review of four genera (Neuroptera, Myrmeleontidea). Beijing: PhD thesis of China Agricultural University[in Chinese]) http://cdmd.cnki.com.cn/Article/CDMD-10019-1014223566.htm
张蕾, 陈小琳, 侯新文, 等. 2011. 实蝇科果实蝇属昆虫数字图像自动识别系统的构建和测试. 昆虫学报, 54(2): 184-196.
(Zhang L, Chen X L, Hou X W, et al. 2011. Construction and testing of automated fruit fly identification system-Bactrocera Macquart (Diptera:Tephritidae). Acta Entomologica Sinica, 54(2): 184-196. [in Chinese])
张秀秀, 管德龙, 覃凯华, 等. 2017. 秦岭地区中华蜜蜂5个地理种群翅形态差异的几何形态测量学研究. 陕西师范大学学报:自然科学版, 45(2): 71-74.
(Zhang X X, Guan D L, Qin K H, et al. 2017. Geometric morphometric analyses on the wing-shape variation in five Apis cerana cerana populations from Qinling Mountains. Journal of Shaanxi Normal University:Natural Science Edition, 45(2): 71-74. [in Chinese])
张亚盟, 魏偏偏, 吴秀杰. 2016. 现代人头骨断面轮廓的性别鉴定-基于几何形态测量的研究. 人类学学报, 35(2): 172-180.
(Zhang Y M, Wei P P, Wu X J. 2016. Sex determination of cross-sectional outlines in modern human skulls:A study based on geometric morphometrics. Acta Anthropologica Sinica, 35(2): 172-180. [in Chinese])
周蜜, 崔娅铭, 邢松. 2016. 中国近代人群上颌前臼齿釉质-齿质连接面形状的三维几何形态测量. 人类学学报, 35(4): 585-597.
(Zhou M, Cui Y M, Xing S. 2016. Analyzing enamel-dentine junction(EDJ) shape in recent Chinese upper premolars using 3D geometric morphometrics. Acta Anthropologica Sinica, 35(4): 585-597. [in Chinese])
朱弘复, 王琳瑶. 1997. 中国动物志:昆虫纲.第十一卷, 鳞翅目:天蛾科. 北京: 科学出版社.
(Zhu H F, Wang L Y. 1997. China Fauna:Insecta. Volume 11, Lepidoptera:Sphingidae. Beijing: Science Press. [in Chinese])
Adams D C, Rohlf F J, Slice D E. 2004. Geometric morphometrics:ten years of progress following the 'revolution'. Ital J Zool, 71(1): 5-16. DOI:10.1080/11250000409356545
Adams D C, Rohlf F J, Slice D E. 2013. A field comes of age:geometric morphometrics in the 21st century. Hystrix, 24(1): 7-14.
Baab K L, McNulty K P, Rohlf F J. 2012. The shape of human evolution:a geometric morphometrics perspective. Evol Anthropol, 21(4): 151-165. DOI:10.1002/evan.v21.4
Bai M, Beutel R G, Song K Q, et al. 2012. Evolutionary patterns of hind wing morphology in dung beetles (Coleoptera:Scarabaeinae). Arthropod Struct Dev, 41(1): 505-513.
Bookstein F L. 1991. Morphometric Tools for Landmark data:Geometry and Biology. New York: Cambridge Univ Press.
Fureix C, Hausberger M, Seneque E, et al. 2011. Geometric morphometrics as a tool for improving the comparative study of behavioural postures. Naturwissenschaften, 98(7): 583-592. DOI:10.1007/s00114-011-0803-2
Gaston K J, O'Neill M A. 2004. Automated species identification: why not? Phil Trans R Soc Lond B, 359(1444): 655-667.
Ge D Y, Chesters D, Gómez-Zurita J, et al. 2011. Anti-predator defense drives parallel morphological evolution in flea beetles. Proc Biol Sci, 278(1715): 2133-2141. DOI:10.1098/rspb.2010.1500
Gutierrez B L, MacLeod N, Edgecombe G D. 2011. Detecting taxonomic signal in an under-utilised character system:geometric morphometrics of the forcipular coxae of Scutigeromorpha (Chilopoda). ZooKeys, 156(4): 49-66.
Hernández L N, Barragán A R, Dupas S, et al. 2010. Wing shape variations in an invasive moth are related to sexual dimorphism and altitude. Bulletin of Entomological Research, 100(5): 529-541.
Hübler N, Klass K D. 2013. The morphology of the metendosternite and the anterior abdominal venter in Chrysomelinae (Insecta:Coleoptera:Chrysomelidae). Arthropod Syst Phyl, 71(1): 3-41.
Kitthawee S, Rungsri N. 2011. Differentiation in wing shape in the Bactrocera tau (Walker) complex on a single fruit species of Thailand. ScienceAsia, 37(4): 308-313. DOI:10.2306/scienceasia1513-1874.2011.37.308
Mitteroecker P, Gunz P. 2009. Advances in geometric morphometrics. Evolutionary Biology, 36(2): 235-247. DOI:10.1007/s11692-009-9055-x
Outomuro D, Adams D C, Johansson F. 2013. Wing shape allometry and aerodynamics in calopterygid damselflies: a comparative approach. BMC Evolutionary Biology, 13.[2017-05-12].http://bmcevolbiol.biomedcentral.com/articles/10.1186/1471-2148-13-118. https://link.springer.com/article/10.1186/1471-2148-13-118
Rohlf F J. 1998a. TpsDig. Ecology and evolution, SUNY. New York: Stony Brook.[2016-08-25]. http://life.bio.sunysb.edu/morph.
Rohlf F J. 1998b. TpsSuper. Ecology and evolution, SUNY. New York: Stony Brook.[2016-08-25]. http://life.bio.sunysb.edu/morph.
Rohlf F J. 1999. Shape statistics:Procrustes superimpositions and tangent spaces. Journal of Classification, 16(4): 197-223.
Slice D E. 2007. Geometric morphometrics. Annu Rev Anthropol, 36(2): 261-281.
Slice D E. 2001. Landmark coordinates aligned by procrustes analysis do not lie in Kendall's shape space. Syst Biol, 50(1): 141-149. https://www.researchgate.net/publication/11261396_Landmark_Coordinates_Aligned_by_Procrustes_Analysis_Do_Not_Lie_in_Kendall's_Shape_Space
Zelditch M L, Swiderski D, Sheets H D. 2004. Geometric Morphometrics for Biologists:A Primer. New York and London: Elsevier Academic Press, 1-437.