Phenotypic variations in wild populations of Paris polyphylla var. chinensis
-
摘要:目的
综合分析华重楼Paris polyphylla var. chinensis (Franch.) Hara野生种群的表型变异规律,为其种质资源保护、药材资源开发利用和繁育改良提供理论支撑。
方法以我国分布区内7个野生华重楼地理种群为研究对象,对其17个表型性状进行系统测量,利用方差分析、变异系数分析、多重比较分析和相关性分析等多种数理统计分析方法,探讨其表型变异规律及其与地理气候因子的关系。
结果除叶片数和花萼片数外,野生华重楼其他15个表型性状在种群间和种群内有显著或极显著的差异;种群内的变异(方差分量占比为52.19%)大于种群间的变异(方差分量占比为26.39%),表型性状的平均分化系数为34.27%,种群内变异是华重楼表型变异的主要来源;各表型性状的平均变异系数为17.72%,变异幅度为6.83%~29.95%,花梗长、叶面积、花萼面积和株高的变异幅度远大于平均变异系数,叶片数、花萼片数、果球横径、果球纵径、果球横纵比和种子千粒质量的变异幅度远小于平均变异系数,表明华重楼花、叶的数量指标及果球、种子的性状稳定性大于叶片形态性状的稳定性;相关分析表明,大多数表型因子与地理气候因子之间存在显著或极显著的相关性。
结论广东清远地理种群变异最为丰富,可作为种质资源收集保存的优选对象。
Abstract:ObjectiveTo comprehensively analyze the phenotypic variations in wild populations of Paris polyphylla var. chinensis (Franch.) Hara, and provide a theoretical support for the protection of germplasm resources, development and utilization of medicinal materials, and breeding improvement.
MethodThe seven wild populations of P. polyphylla var. chinensis in China were used as materials. Based on the measurement of 17 phenotypic traits, by the analyses of variance, variation coefficient, multiple comparisons and correlation, the phenotypic variation and its relationship with geographical climate factors were discussed.
ResultExcept the numbers of leaf and calyx, the other 15 phenotypic traits of wild P. polyphylla var. chinensis showed significant or extremely significant differences among and within populations. The variation within populations (Percentage of variance component was 52.19%) was greater than that among populations (Percentage of variance component was 26.39%). The mean phenotypic differentiation coefficient of phenotypic traits was 34.27% and the main source of variation was within populations. The average variation coefficient of each phenotypic trait was 17.72% and variation coefficient ranged from 6.83% to 29.95%. The variation ranges of peduncle length, leaf area, calyx area and plant height were higher than the average variation coefficients of these five traits, while those of leaf number, calyx number, fruit transverse diameter, fruit longitudinal diameter, fruit transverse diameter to longitudinal diameter ratio and seed 1000-grain weight were far less than the mean variation coefficients, indicating that the stabilities of flower and leaf quantitative trait indexes, fruit and seed traits were greater than those of leaf shape traits. Correlation analysis showed that most of the 17 phenotypic traits had significant or extremely significant correlations with the geographical and climatic factors.
ConclusionThe wild populations of P. polyphylla var. chinensis in Qingyuan City, Guangdong Province have the most abundant geographical population variations and can be used as the optimal target for germplasm collection and preservation.
-
玉米Zea mays L.作为我国三大主要粮食作物之一,具有种植成本低、周期短、收益快的优势,保证玉米作物的高质量生产是保障国家粮食安全、促进农民增收以及农业经济发展的有效措施[1]。 然而,在玉米生长过程中面临着多种挑战,包括玉米生长状况监测、病虫害预测和防治等问题。因此,高效、准确、快速、及时地获取玉米作物的生长信息对玉米生长监测和管理具有重要意义。地上部生物量和叶面积指数(Leaf area index,LAI)是反映作物生长状况的重要农艺参数,且二者间存在密切的联系。在玉米不同的生长发育阶段,由于作物冠层结构、叶片形态以及生理生态特征的变化,作物冠层在图像中表现出的光谱特征亦有不同[2, 3]。目前,测量作物LAI主要包括地面采样以及遥感与地面测量结合这2种方法,其中,遥感与地面测量结合的方法更为精确,可用于验证遥感结果的准确性[4]。广泛应用于作物长势监测研究的光谱植被指数法为玉米监测提供了有效的工具。植被指数(Vegetation index,VI)是利用2个或多个波段的反射率经过线性或非线性组合运算增强植被信息,以削弱环境背景对植被光谱特征的干扰,如归一化植被指数(Normalized difference vegetation index, NDVI)、增强植被指数(Enhanced vegetation index, EVI)和差值植被指数(Difference vegetation index, DVI)等。目前,光谱植被指数法已经作为一种快速获取地表信息的手段,被国内外学者用于作物长势监测研究[5]。杨贵军等[6]通过无人机搭载多传感器,实现了对作物产量、LAI、冠层温度等多种信息的监测。邵国敏等[7]基于无人机获取多光谱数据,在水分胁迫条件下,植被指数的随机森林回归模型可以较好地估算LAI,决定系数(R2)在0.74~0.87之间,表明使用随机森林回归算法估算多种灌溉条件下的夏玉米LAI是可行的。张瑾[8]基于无人机多光谱遥感数据和采用支持向量机构建了LAI估测模型,模型的R2可达到0.701;并通过粒子群算法进行模型优化,优化后模型的R2提高至0.822,均方根误差(RMSE)为0.267,证明粒子群算法可以有效地提高模型的预测精度。孙诗睿等[9]基于无人机多光谱影像,通过对比不同植被指数建模方法分析了冬小麦叶面积指数的反演精度,发现基于多植被指数的随机森林算法预测具有更好的拟合效果。鲸鱼优化算法(Whale optimization algorithm ,WOA)以其高效的参数优化能力和出色的全局搜索能力而闻名,能够更好地优化数据模型以适应玉米监测需求。此外,WOA具有自适应性,可以自动调整搜索策略,从而提高模型的鲁棒性。He等[10]为了研究土壤湿度对玉米产量的影响,将WOA与正弦混沌图相结合,并与支持向量机(Support vector machine, SVM)结合使用,结果显示,SVM-WOA相对于SVM在60 cm深度的平均绝对误差(MAE)、均方根误差(RMSE)、平均绝对百分比误差(MAPE)和平均偏差(MBE)方面分别提高了14.0%、13.0%、41.5%和14.0%,表明经过WOA算法优化后的模型可以更好地为智能农业和精准灌溉提供指导。将WOA算法应用于无人机遥感预测农业数据领域方面的研究鲜见报道,为进一步提高多光谱数据处理和植被指数计算的准确性,本研究对WOA算法进行优化,通过多光谱遥感无人机获取玉米作物冠层多光谱信息构建植被指数,研究植被指数与玉米LAI之间的定量关系;并对筛选出的植被指数分别使用多元线性回归、SVM、随机森林回归算法(Random Forest, RF)和基于WOA优化的随机森林算法(WOA-RF)构建玉米LAI预测模型,选择最优模型以提高玉米LAI的预测精度,为玉米监测和农田管理提供指导和数据参考。
1. 材料与方法
1.1 试验材料
本试验在广东省广州市增城区华南农业大学教学科研基地(113°40′E,23°24′N )展开研究,如图1所示,a为华南农业大学科研基地,b为本次试验的试验田,试验地海拔约为20 m,属于亚热带气候,温暖湿润,适合种植多种作物。玉米品种为‘华珍’甜玉米,播种日期为2022年8月15日,玉米的种植密度相同;播种前施入适量的基础肥料,拔节期和抽雄期追施氮肥。将试验田分为6个等份区域,每个区域标记3株玉米并实测玉米叶面积。田间玉米群体采集系统主要有计算机(配套软硬件)、大疆精灵4多光谱版一体化无人机(Phantom4-M,P4M)、电池若干、漫反射标准参照板(MFB99-16-16)以及活体叶面积测量仪(YMJ-D)。所用软件有ENVI、Pix4Dmapper、ArcMap、SPSS以及Matlab。
1.2 无人机影像获取
图2中,a为无人机航线,图中红色点为拍摄点;b为图像拍摄所在位置和角度,每次航拍作业可获取1 800张原始多光谱影像。采用P4M获取试验区域多光谱影像数据,该设备集成1个可见光传感器通道和5个多光谱传感器通道(蓝光、绿光、红光、红边和近红外),每次拍摄可以获得6张分辨率为208万像素的影像,最高飞行速度14 m/s,最大续航时间27 min,搭配有TimeSync时间同步系统,可获得厘米级定位精度。此外,P4M顶部集成光强传感器可捕捉太阳辐照度数据,用于影像的光照补偿,以排除环境光对数据的干扰,提高不同时段采集数据的准确度与一致性。北京时间上午10:30—11:30期间拍摄。飞行高度设置为30 m,航向、旁向重叠率分别设置为70%和80%,在保持光源稳定、无风的情况下,获取图像序列。
1.3 图像处理
1.3.1 多光谱影像预处理
利用漫反射标准参照板对无人机拍摄的原始多光谱影像进行辐射标定。通过漫反射标准参照板反射比数据获取各个波段的波长,计算参照板反射率,获取同一时期、相同条件下的多光谱影像,共同完成辐射校正[11]。在Envi软件中,对多光谱无人机拍摄的一组参照板影像进行勾选,通过已知参照板的反射率,对多光谱5个波段进行校正,可以获得多光谱影像的各个反射率影像,从而得到试验区域内被校正后垂直于地面的多光谱影像。在Envi软件中将影像叠加为多波段后,查看典型植被的反射率,以检验校正后数据是否正常。辐射校正可以有效地减少地形变形[12]、投影偏移[13]等问题,以便得到更真实的试验数据。
1.3.2 影像拼接与光谱信息提取
将矫正后的影像导入Pix4Dmapper软件中进行拼接,得到完整的试验区多光谱影像。将拼接后的多光谱反射率影像导入Arcgis软件中,设置正确的图像投影,确定坐标系统与影像坐标系统一致。选择WSG_1984坐标系建立Shp图层,勾画研究区内多光谱反射率影像[14],分成6个待提取小区的范围,规划小区1作为验证区,其余5个小区作为训练区(图3)。使用Envi软件中Band Math功能,通过选择多光谱影像中的5个波段,分别计算它们的反射率,获得各区域植物光谱信息。在此基础上,计算出各小区平均反射率,作为每个区域的光谱反射率,从而求得各植被指数。
1.4 地面实测数据获取
为了研究试验区内玉米植株LAI变化情况,本试验使用了活体叶面积测量仪。手持扫描仪器在试验区内扫描测量了玉米植株地上部所有叶片的叶面积,将这些测量结果相加得到叶面积总和。这一步骤旨在获取试验区内玉米植株的实际叶面积数据,为后续的分析和研究提供重要数据支持。通过使用活体叶面积测量仪,可以准确地获取玉米植株的叶面积信息,进而对其生长状态和生理特征进行详细分析。
1.5 植被指数的选取
多光谱影像是一种常用的监测和评估玉米生长状况的工具,其包含红(R)、绿(G)、蓝(B)、红边(RE)和近红外(NIR)5个光谱通道,通过计算多种植被指数来反映玉米的生长状态。为了确定最适合的植被指数,在玉米不同生长阶段选择了10种植被指数,并对其进行了比较和分析。通过对比不同植被指数的表现,最终选取了最相关的NDVI[15]、NDRE、EVI和 CIG (Chlorophyll Index Green)[16]4种光谱指数,以确保获得最准确的玉米生长状况数据。本试验使用的10个植被指数计算公式如下:
$$ {\text{NDVI = (NIR} - \text{Red) / (NIR + Red),}} $$ (1) $$ {\text{NDRE = (NIR} - \text{RE) / (NIR + RE),}} $$ (2) $$ {\text{PVI = (Red} - \text{Green) / (Red + Green),}} $$ (3) $$ {\text{GNDVI = (NIR} - \text{Green) / (NIR + Green),}} $$ (4) $$ \text{SAVI = (NIR} - \text{Red) / (NIR + Red + 0}\text{.5) × 1}\text{.5}\text{,} $$ (5) $$ \text{OSAVI = (NIR} - \text{Red) / (NIR + Red + 0}\text{.16) × 1}\text{.16}\text{,} $$ (6) $$ \begin{split} \text{EVI} = & 2\text{.5 × [(NIR} - \text{Red) / (NIR + 6 × Red }-\\ & 7\text{.5 × Blue + 1)]}\text{,} \end{split} $$ (7) $$ {\text{CIG = NIR / Green} - \text{1}}, $$ (8) $$ \text{RDVI = (NIR} - \text{Red) / (NIR + Red)}\text{,} $$ (9) $$ \text{TSAVI = 1}\text{.5 × (NIR} - \text{Red) / (NIR + Red + 0}\text{.5)}\text{,} $$ (10) 式中,NIR、Red、RE、Green、Blue分别表示近红外波段、红色波段、红边波段、绿色波段、蓝色波段的反射率。
1.6 模型构建与精度验证
采用3种不同的算法构建玉米LAI反演模型。多元线性逐步回归(Stepwise LR)方法用于优选植被指数和玉米LAI反演模型,通过筛选自变量以构建准确预测模型[17]。SVM用于评估其在LAI反演模型中的适用性,是一种统计学习算法,可用于分类和回归问题[18]。RF分析玉米LAI反演模型基于决策树的集成学习思想,通过随机特征选择和多数投票等策略构建强大的分类模型[19]。对比分析选出最优预测模型后,使用WOA优化,得到最终预测模型。
1.6.1 WOA优化
WOA作为本试验的模型优化算法,是一种自然启发式算法,灵感来自于鲸鱼捕食行为。它通过模拟鲸鱼群体的行为来寻找最优解。WOA的核心是维护一群鲸鱼,这些鲸鱼在搜索空间中游动,并根据其位置和适应度函数值来调整自己的位置。WOA通过迭代的方式搜索超参数空间,尝试不同的超参数组合。使用一群鲸鱼(空间中的候选超参数)来共同寻找最佳解[20]。
WOA在优化过程中,会尝试不同的超参数组合,不断迭代以寻找更好的超参数值。这个过程会一直进行,直到达到预定的最大迭代次数或其他停止条件。返回最佳的超参数设置,这些参数将用于构建RF模型,以使模型在给定的数据上达到最佳性能。它可以自动选择超参数,不需要手动调整,从而提高了模型的效果。
如图4所示,首先,初始化所需的参数,包括种群数量(pop)、最大迭代次数(Max_iter)、超参数维度(dim)、上下边界(lb & ub)。通过WOA函数对RF回归模型的超参数进行优化,最终得到最佳的超参数。
1.6.2 模型验证
采用决定系数(R2)、RMSE和MAE这3种指标来检验和评价模型性能。R2衡量拟合度,越接近1表示模型拟合效果越好;RMSE衡量预测值与真实值差异,较小值表示较好的预测效果;MAE评估平均绝对差异,值越小表示预测效果越好。使用R2、RMSE和MAE这3种指标来检验和评价模型的拟合效果,判断模型是否适用于未来结果的预测。
2. 结果与分析
2.1 玉米叶面积的实测结果
试验时间为2022年9月10日至2022年11月9日,试验区玉米经历了拔节期、抽雄期、喇叭口期、穗期、乳熟期和成熟期共6个生长发育阶段。采集数据频率为每10 d采集1次,在研究区域内每个小区采集3株,共18株玉米叶面积样本,7次采集共得到126株玉米叶面积总样本。这些数据将用于分析玉米叶面积随生长发育阶段的变化情况,为研究玉米生长过程中叶面积的动态变化提供数据支持。在试验区域内采集的126株玉米叶面积样本,能够为后续的数据分析和模型建立提供充分的数据支持。
使用活体叶面积测量仪作为测量工具来测量叶面积,测量试验区内已标记的18株玉米地上部。首先,测量所有叶片的叶面积总和(S1);然后,根据小区的行数(m)、行长(L)以及行间距(D)计算出小区的面积(S2);最后,根据下式计算出每个小区对应的LAI:
$$ S_2=m \times L\times (L+D), $$ (11) $$ {\mathrm{LAI}}=S_1/S_2。 $$ (12) 以上步骤将被重复执行6次,通过计算LAI来了解不同试验小区内玉米植株的生长情况,获得各个试验区内玉米植株LAI的变化趋势。6个试验区的平均LAI变化如图5所示。玉米生长周期可分为拔节期、抽雄期、喇叭口期、穗期、乳熟期和成熟期6个阶段。在拔节期,植株快速生长,叶片逐渐展开,导致LAI逐渐上升;随后进入抽雄期,生长速度依然较快,叶片覆盖面积继续增加,LAI持续上升;接着是喇叭口期,植株生长达到高峰,叶片面积达到最大值;随后进入穗期,虽然生长速度减缓,但LAI仍保持在较高水平;乳熟期中,植株生长停滞,穗开始成熟,叶片老化导致LAI略微下降;最终进入成熟期,植株生长停止,叶片逐渐凋落,LAI持续下降。
2.2 植被指数的筛选结果
通过从多光谱影像中提取的5个波段下的玉米冠层反射率值,可以计算出各植被指数,并将玉米LAI数据与相应的植被指数数据输入SPSS软件进行相关性分析。表1可知,不同生育时期玉米LAI与10个植被指数之间均存在正相关关系;通过分析对比,选择出4种高度相关的植被指数EVI、CIG、NDRE和NDVI,这4个植被指数与玉米LAI的相关系数分别为0.884、0.850、0.836和0.824,均达极显著的正相关(P < 0.01)。这4种植被指数对植被的不同生理和生态学特征具有敏感性,包括叶绿素含量、植被覆盖度和生长状态等。NDVI、NDRE和EVI都使用了红外波段的信息,植被在红外波段通常具有较高的反射率,而在可见光波段有较低的反射率,这使得这些指数能够更好地区分植被和非植被区域。所选植被指数能够更精确地估算玉米的LAI,并可作为构建玉米LAI估算模型的参数。
表 1 植被指数与叶面积指数相关性分析Table 1. Correlation analysis between vegetation index and leaf area index植被指数
Vegetation index皮尔逊相关系数
Pearson correlation coefficient显著性
Sig.EVI 0.884 0.001 CIG 0.850 0.001 NDRE 0.836 0.003 NDVI 0.824 0.002 SAVI 0.793 0.003 OSAVI 0.774 0.004 PVI 0.765 0.007 GNDVI 0.763 0.007 RDVI 0.643 0.010 TSAVI 0.571 0.012 将玉米LAI作为因变量y,4个植被指数NDVI、NDRE、EVI和CIG分别作为自变量x,通过SPSS软件进行一元线性回归分析。结果发现,不同生育时期的4个植被指数模型的R2分别为0.680、0.704、0.775、0.726,均大于0.6,表明模型的拟合度高,可用于预测玉米的生长状态和产量。其中, EVI与玉米LAI的相关性最好,R2 达到0.775,表明EVI能够更准确地反映玉米生长状态与LAI之间的关系。EVI综合了红外、近红外和蓝波段的反射率,并对大气干扰和土壤背景噪声具有一定的抵抗能力,可以更好地反映玉米生长期间的生长状态和LAI变化。
2.3 不同模型的预测结果
2.3.1 多元线性回归模型
根据逐步回归模型的分析,得到了最优的多元线性回归模型,用于预测LAI,结果见表2。由表2可见,植被指数EVI、CIG、NDVI和NDRE与LAI之间均存在显著相关(Sig. < 0.05)。另外,模型还进行了多重共线性的检验,结果显示这些植被指数之间的方差膨胀因子(VIF)均小于4,没有出现严重的多重共线性问题,可以用来构建玉米LAI的预测模型。在多元线性回归模型中,常量代表截距,表示当所有自变量(植被指数)为0时,因变量(LAI)的预期值。在表2中,常量B和Sig.分别表示截距项的回归系数和显著性水平,B和Sig.在模型中的作用是为LAI提供一个基准预测值,并检验这个基准值的统计置信度。线性回归模型在描述LAI与多光谱植被指数之间的关系方面表现出较高的拟合度(R2=0.873)和较小的均方根误差(RMSE=0.277),表明该模型能够很好地描述LAI与多光谱植被指数之间的关系,并且对未知数据的预测能力较强。
表 2 多元线性回归模型最优植被指数1)Table 2. Optimal vegetation index for multiple linear regression model模型Model B Sig. VIF R2 RMSE 常量 Intercept −1.021 0.092 0.873 0.277 EVI 3.314 0.001 3.916 CIG 0.263 0.001 3.082 NDVI 3.308 0.007 2.828 NDRE 2.871 0.020 3.262 1) B:回归系数,Sig.:显著性水平,VIF:多重共线性检验中的方差膨胀因子
1) B:Regression coefficient, Sig.: Significance level, VIF: Variance inflation factor in multicollinearity test2.3.2 支持向量机(SVM)模型
采集了126株玉米样本,为保证模型的准确性,将其中的105株样本用于训练模型,剩下21株样本用于测试模型的精度。将读取的每株玉米叶面积数据和对应的植被指数数据作为输入,用SVM算法进行训练,预测玉米LAI。结果(表3和图6)表明,训练集数据R2 为0.895 5、RMSE为0.151 1、MAE为0.222 6,测试集数据R2为0.878 0、RMSE为0.182 1、MAE为0.236 5。表明该模型在解释训练集和测试集数据的变异性方面表现出很好的能力,同时也具有较小的平均误差,RMSE较小。SVM模型表现出良好性能,对训练数据的解释能力强,同时在未见过的测试数据上也表现出良好的预测能力,表明其在预测玉米LAI方面具有高性能和泛化能力。
表 3 不同模型性能指标对比Table 3. The comparison of performance metrics of various models模型 Model 样本 Sample R2 RMSE MAE 多元线性回归 MLR 0.873 2 0.277 5 0.193 2 SVM 训练集 Training set 0.895 5 0.222 6 0.151 1 测试集 Test set 0.878 0 0.236 5 0.182 1 RF 训练集 Training set 0.963 0 0.132 0 0.151 2 测试集 Test set 0.917 7 0.209 0 0.134 9 WOA-RF 训练集 Training set 0.966 8 0.125 0 0.124 3 测试集 Test set 0.940 8 0.128 7 0.134 5 2.3.3 随机森林回归模型
原始数据集包含126个玉米样本的LAI以及与之相关的4个植被指数。为了确保随机性,将原始数据集按照8∶2的比例进行随机划分,得到训练数据集(100个样本)和预测数据集(26个样本)。训练数据集将用于训练回归模型,而预测数据集则用于验证模型的准确性。RF模型的超参数选择构成森林的决策树数量(ntree)和最小叶子节点数(Min_samples_leaf),增加ntree值可以降低RF模型的误差,如果决策树的数量过多,计算时间会增加,导致模型的运行时间被浪费。因此,需要选择合适的ntree值进行优化,通过使用Matlab编程建立RF模型,并分别设定ntree值为100、200、300和400,分析不同ntree值时的拟合系数R2 (图7a)。图7a表明,当ntree为300时, R2最高,效果最好。分别设定最小叶子节点数为1、2、3、4、5和6,并使用Matlab中的4折交叉验证方法来评估不同最小叶子节点数下的模型性能,结果(图7b)表明,不同最小叶子节点数对应的R2值范围为0.8~0.9,且最小叶子节点数为4的R2最高,效果最好。通过对ntree和Min_samples_leaf进行调优,选择了ntree为300、Min_samples_leaf为4的超参数配置,建立了一个RF模型,该模型在数据集上表现最佳。
确定模型的超参数后,在Matlab上进行回归分析。首先,导入数据,将数据划分为训练集和测试集。对数据进行归一化处理,以便更好地训练模型;然后,使用训练集数据训练模型,并使用测试集数据进行仿真测试。完成预测后,将预测结果反归一化为原始数据,通过计算RMSE和相关指标对模型进行评估。结果如表3和图8所示,RF模型训练集数据R2为0.9630、RMSE为0.1320、MAE为0.1512,测试集数据R2为0.9177、RMSE为0.2090、MAE为0.1349,该模型在训练集和测试集上都表现出较好的性能,具有较高的R2和较小的MAE和RMSE。表明该模型具有出色的数据拟合和准确的测试数据预测能力,能够很好地拟合训练数据,并在测试数据上有较准确的预测。
2.3.4 WOA-RF优化算法
为了确保随机性,采用了8∶2的比例将原始数据集随机分成2部分,即训练数据集(包含100个样本)和预测数据集(包含26个样本)。通过WOA-RF优化算法寻找随机森林算法最佳的超参数,得到了最佳随机森林的决策树数量和最佳叶子节点的最小样本数。如表3和图9所示,WOA-RF模型训练集数据R2为0.966 8、RMSE为0.125 0、MAE为 0.1243,测试集数据R2为0.940 8、RMSE 为0.1287、MAE 为0.1345;该模型在训练集和测试集上都表现出色的性能,具有较高的R2以及较小的MAE和RMSE,表明模型能够很好地拟合训练数据,获得较准确的预测。
2.4 不同模型性能比较分析
由表3可知,多元线性回归模型(MLR)在描述LAI与多光谱植被指数之间的关系方面表现出较高的拟合度和较小的均方根误差,表明该模型能够很好地描述LAI与多光谱植被指数之间的关系,并且对未知数据的预测能力较强。SVM模型表现出良好性能,对训练数据的解释能力强,同时在未见过的测试数据上也表现出良好的预测能力,表明其在预测玉米LAI方面具有高性能和泛化能力。RF模型表现出出色的性能,拥有高的R2和低的MAE和RMSE,不论是在训练集还是测试集上,证明了其出色的数据拟合和准确的测试数据预测能力。
3. 讨论与结论
多元线性回归、SVM和RF这 3种回归模型都在不同程度表现出了良好的性能,但RF模型在训练和测试集上都表现出了更高的准确性,具有更强的泛化能力,因此被认为是最适合用于预测植被指数和与LAI之间关系的模型。通过WOA算法对RF模型进行优化,可得到更高效准确的预测模型,各评价指标均高于单一模型预测精度。R2提高了约2.31%,RMSE减小了约38.49%,MAE减小了约0.30%。在优化机器学习模型中,WOA算法可以用于调整模型参数,提高模型性能。如张伟萍等[21]等提出一种针对非平稳水声信号的改进WOA-VMD小波阈值去噪方法,利用WOA选取水声信号变分模态分解(VDM)的最佳参数组合,通过仿真试验和实测数据的降噪分析,证明了使用WOA优化算法后去噪性能优于其他方法。WOA优化算法在农业领域也逐渐热门起来,赵丙秀等[22]利用WOA对BP神经网络模型进行优化来预测马铃薯产量,结果表明WOA-BP神经网络模型精度明显高于BP神经网络模型。本文通过WOA算法对RF进行优化后,预测模型的性能得到了显著改善,WOA-RF模型比RF模型有更高的预测精度,更适合基于多光谱无人机对玉米LAI的预测研究。
目前,通过无人机多光谱对作物LAI进行预测的研究很多,但大多采用的是单一模型预测,本文通过多光谱遥感无人机获取玉米作物冠层多光谱信息以构建植被指数,研究植被指数与玉米LAI之间的定量关系,对筛选后的植被指数分别使用不同的机器学习方法构建LAI预测模型,对玉米LAI的预测精度进行预测研究。结果表明, WOA算法对表现最优的RF模型进行优化后,各评价指标均高于单一模型预测精度,优化后的模型R2为0.9408、MAE为0.1345、RMSE为0.1287,预测精度可以满足玉米生产的实际需要,为农业种植管理和决策提供有益参考。然而,许多因素会对无人机多光谱在玉米LAI反演中的准确度产生影响,包括季节和生长阶段、波段选择、辐射定标和大气校正、数据处理算法等。本研究数据采集仅覆盖了拔节期至抽雄期,未考虑喇叭口期和穗期的变化,可能影响模型的泛化能力。因此,在将这一技术推广至其他地区或不同种植条件时,必须考虑外部因素对模型性能的潜在影响。
-
表 1 7个野生华重楼地理种群的地理及主要气候因子情况
Table 1 Sample collection information of Paris polyphylla var. chinensis from different distribution regions and climatic data at the sampling sites
样品采集点
Sample site种群
Population经度
Longitude纬度
Latitude海拔/m
Altitude年均温度/℃
Mean annual temperature无霜期/d
Frost-free period年日照/h
Annual sunshine hour年降水量/mm
Annual precipitation江西新余
Xinyu in JiangxiJXXY 114°37′22″E 27°31′45″N 376 17.7 270 1 655.4 1 594.8 浙江丽水
Lishui in ZhejiangZJLS 119°04′18″E 28°38′28″N 1 103 11.3 240 1 635.1 1 598.9 四川雅安
Yaan in SichuanSCYA 103°01′57″E 30°30′24″N 1 450 17.0 295 1 519.2 1 350.0 广东清远
Qingyuan in GuangdongGDQY 112°20′11″E 25°10′10″N 648 18.0 210 2 290.0 1 329.0 福建三明
Sanming in FujianFJSM 117°29′18″E 26°33′24″N 548 18.2 260 1 727.1 1 688.0 湖南娄底
Loudi in HunanHNLD 112°02′62″E 27°48′05″N 290 17.1 270 1 726.0 1 455.6 湖北宜昌
Yichang in HubeiHBYC 111°21′19″E 30°27′44″N 400 16.9 275 1 710.5 1 215.5 表 2 野生华重楼种群间和种群内全株、果、叶及种子表型性状的方差分析结果1)
Table 2 Variance analysis of phenptypic traits of wholeplant, leaf, fruit and seeds among and within Paris polyphylla var. chinensis
性状
Trait均方 Mean square 随机误差
Random
errorF 种群间
Among
populations种群内
Within
populations种群间
Among
populations种群内
Within
populations株高 Height 1 060.94 226.26 266.13 4.69** 5.97** 地径 Diameter 6.39 1.32 1.62 4.83** 7.13** 叶片数 No. of leaf 0.53 0.78 0.08 0.68 1.07 叶长 Leaf length 5 728.49 601.72 264.04 9.52** 7.18** 叶宽 Leaf width 1 043.13 43.62 51.48 23.91** 20.99** 叶面积 Leaf area 2 831.39 142.52 138.48 19.87** 3.56* 叶片长宽比 Blade aspect ratio 3.73 0.84 0.15 4.45** 5.39** 花萼片数 No. of calyx 4.52 0.55 1.22 0.82 0.79 花萼长 Calyx length 762.13 49.75 46.97 15.32** 10.32** 花萼宽 Calyx width 119.22 9.90 6.86 12.04** 11.31** 花萼面积 Calyx area 46.00 3.32 2.81 13.86** 7.15** 花萼长宽比 Calyx aspect ratio 0.55 0.13 0.03 4.12** 5.01** 花梗长 Pedicel length 152.91 70.67 25.22 3.16* 6.25** 果球横径 Fruit transverse diameter 89.74 18.79 5.95 4.78** 6.17* 果球纵径 Fruit longitudinal diameter 256.18 25.35 19.35 10.10** 8.33** 果球横纵比 Fruit transverse diameter
to longitudinal diameter ratio0.24 0.11 0.02 21.60** 17.05** 种子千粒质量 Seed 1 000 -grain weight 843.51 69.95 229.20 12.06** 10.87** 1)“*”和“**”分别表示在0.05和0.01水平差异显著(单因素方差分析)
1)“*” and “**” mean significant differences at 0.05 and 0.01 levels respectively(One-way ANOVA)表 3 野生华重楼7个地理种群20个表型性状的平均值、标准差及多重比较1)
Table 3 Mean, standard deviation and multiple comparison of 20 phenotypic traits in seven geographical populations of Paris polyphylla var. chinensis
种群2)
Population株高/cm
Height地径/mm
Diameter叶片数
No. of leaf叶长/mm
Leaf length叶宽/mm
Leaf width叶面积/cm2
Leaf areaJXXY 53.67±17.64ab 8.22±1.17ab 6.33±0.58a 134.91±21.64ab 44.89±10.01b 48.84±17.70b ZJLS 66.43±8.46a 7.93±0.70abc 6.67±1.15a 121.90±13.43bc 48.31±4.10ab 46.48±7.92b SCYA 69.20±12.51a 8.63±1.33a 7.00±1.00a 148.38±12.29a 51.85±4.90a 60.28±5.99a GDQY 49.63±24.01abc 6.47±0.84bcd 6.67±1.53a 114.58±55.20c 32.73±7.84d 30.50±21.37cd FJSM 35.33±5.81bc 5.29±0.78d 6.33±0.58a 123.02±10.82bc 35.80±3.96cd 34.73±5.96c HNLD 22.13±2.25c 6.05±0.84cd 6.67±0.58a 112.69±10.17c 38.91±5.91c 34.57±6.84c HBYC 25.83±2.25c 5.02±0.59d 5.67±0.58a 93.57±6.93d 35.70±7.53cd 26.35±6.28d 种群2)
Population叶片长宽比
Blade aspect ratio花萼片数
No. of calyx花萼长/mm
Calyx length花萼宽/mm
Calyx width花萼面积/cm2
Calyx area花萼长宽比
Calyx aspect ratioJXXY 3.06±0.37bc 5.67±0.58a 49.91±7.87a 19.13±4.13a 7.68±2.56a 2.66±0.34ab ZJLS 2.53±0.24c 6.67±1.15a 49.34±5.34a 18.46±1.74ab 7.20±1.29a 2.68±0.26a SCYA 2.89±0.41b 5.80±0.84a 52.06±8.90a 19.23±3.43a 8.03±2.56a 2.73±0.35a 续表 3 Continued table 3 种群2)
Population叶片长宽比
Blade aspect ratio花萼片数
No. of calyx花萼长/mm
Calyx length花萼宽/mm
Calyx width花萼面积/cm2
Calyx area花萼长宽比
Calyx aspect ratioGDQY 3.79±2.18a 5.33±0.58a 38.80±8.31b 13.79±1.90c 4.29±1.37b 2.81±0.42a FJSM 3.46±0.38ab 4.67±0.58a 36.71±8.41b 13.42±2.52c 3.97±1.59b 2.76±0.52a HNLD 2.95±0.45bc 5.67±0.58a 39.23±5.24b 16.56±2.48b 5.18±1.42b 2.39±0.25bc HBYC 2.72±0.51c 4.33±0.58a 34.19±2.86b 13.04±4.75c 4.26±1.16b 2.26±0.48c 种群2)
Population花梗长/cm
Pedicel length果球横径/mm
Fruit transverse diameter果球纵径/mm
Fruit longitudinal diameter果球横纵比
Fruit transverse diameter to longitudinal diameter ratio种子千粒质量/g
Seed 1 000-grain weightJXXY 18.17±5.03ab 30.45±5.25bc 35.93±7.28ab 0.86±0.10d 100.16±3.00a ZJLS 19.60±2.15ab 34.09±2.84ab 39.94±3.37a 0.85±0.04d 72.25±2.55b SCYA 18.60±3.47ab 29.50±4.41c 34.40±4.34bc 0.86±0.08d 73.25±2.55b GDQY 23.00±15.86a 29.79±3.46bc 26.33±4.40d 1.14±0.12b 77.49±10.10b FJSM 7.83±0.91ab 28.92±2.65c 30.47±3.37cd 0.96±0.88cd 95.76±4.17a HNLD 10.43±4.29ab 36.83±6.65a 36.06±5.95ab 1.02±0.11c 109.99±17.00a HBYC 5.80±1.06b 32.68±2.79abc 26.64±3.94d 1.25±0.19a 68.14±3.37b 1)同列数据后的不同小写字母表示差异显著(P<0.05,Duncan’s法);2)JXXY:江西新余;ZJLS:浙江丽水;SCYA:四川雅安;GDQY:广东清远;FJSM:福建三明;HNLD:湖南娄底;HBYC:湖北宜昌
1) Different lowercase letters in the same column indicated significant differences among different geographical populations (P<0.05, Duncan’s method);2) JXXY: Xinyu in Jiangxi; ZJLS: Lishui in Zhejiang; SCYA: Ya’an in Sichuan; GDQY: Qingyuan in Guangdong; FJSM: Sanming in Fujian; HNLD: Loudi in Hunan; HBYC: Yichang in Hubei表 4 野生华重楼表型性状的方差分量及种群间表型分化系数
Table 4 Variance component and differentiation coefficients of phenotypic traits among and within populations in Paris polyphylla var. chinensis
性状
Trait方差分量
Variance component方差分量占比/%
Percentage of variance component表型分化
系数/%
Phenotypic differentiation
coefficient种群间
Among populations种群内
Within populations随机误差
Random error种群间
Among populations种群内
Within populations随机误差
Random error株高 Height 352.96 464.74 266.13 32.57 42.88 24.55 43.16 地径 Diameter 2.13 2.77 1.62 32.67 42.48 24.85 43.47 叶片数 No. of leaf 0.18 0.72 0.08 18.37 73.47 8.16 20.00 叶长 Leaf length 301.64 829.57 264.04 21.62 59.46 18.92 26.67 叶宽 Leaf width 52.36 88.04 51.48 27.29 45.88 26.83 37.29 叶面积 Leaf area 144.13 262.02 138.48 26.46 48.11 25.43 35.49 叶片长宽比 Blade aspect ratio 0.19 0.97 0.15 14.50 74.05 11.45 16.38 花萼片数 No. of calyx 1.47 1.63 1.22 34.03 37.73 28.24 47.42 花萼长 Calyx length 53.20 90.08 46.96 27.96 47.35 24.68 37.13 花萼宽 Calyx width 7.75 15.81 6.86 25.48 51.97 22.55 32.89 花萼面积 Calyx area 3.14 5.74 2.81 26.86 49.10 24.04 35.36 花萼长宽比 Calyx aspect ratio 0.04 0.16 0.03 17.39 69.57 13.04 20.00 花梗长 Pedicel length 44.11 93.10 25.22 27.16 57.32 15.53 32.15 果球横径 Fruit transverse diameter 8.46 23.80 5.95 22.14 62.29 15.57 26.22 果球纵径 Fruit longitudinal diameter 26.53 41.64 19.35 30.31 47.58 22.11 38.92 果球横纵比 Fruit transverse diameter to longitudinal diameter ratio 0.02 0.03 0.02 28.57 42.86 28.57 40.00 种子千粒质量 Seed 1 000-grain weight 272.61 271.75 229.20 35.24 35.13 29.63 50.08 平均值 Mean 26.39 52.19 21.42 34.27 表 5 野生华重楼种群17个表型性状变异系数
Table 5 Variation coefficients of 17 phenotypic traits in Paris polyphylla var. chinensis populations
% 性状
Trait种群1) Population 平均值
MeanJXXY ZJLS SCYA GDQY FJSM HNLD HBYC 株高 Height 32.86 12.73 18.08 48.40 16.44 72.02 8.73 29.89 地径 Diameter 14.22 8.86 15.42 27.17 14.79 13.89 11.80 15.16 叶片数 No. of leaf 9.12 17.32 14.29 17.28 9.12 9.12 10.19 12.35 叶长 Leaf length 16.04 11.02 8.28 48.18 8.79 9.03 7.40 15.53 叶宽 Leaf width 22.31 8.48 9.46 23.96 11.05 15.18 21.10 15.93 叶面积 Leaf area 36.24 17.04 9.93 70.08 17.15 19.78 23.81 27.72 叶片长宽比 Blade aspect ratio 12.06 9.35 14.31 57.54 10.96 15.21 18.87 19.76 花萼片数 No. of calyx 10.19 17.32 10.83 14.43 15.75 10.19 17.32 13.72 花萼长 Calyx length 15.76 10.81 17.10 21.42 22.90 13.35 8.38 15.67 花萼宽 Calyx width 21.60 9.42 17.83 13.77 18.75 15.00 36.47 18.98 花萼面积 Calyx area 33.36 17.87 31.85 31.98 39.97 27.35 27.26 29.95 花萼长宽比 Calyx aspect ratio 12.86 9.70 12.68 15.02 19.14 10.61 21.25 14.47 花梗长 Pedicel length 27.71 10.98 18.63 68.72 11.58 41.07 18.25 28.13 续表 5 Continued table 5 性状
Trait种群1) Population 平均值
MeanJXXY ZJLS SCYA GDQY FJSM HNLD HBYC 果球横径 Fruit diameter 17.26 8.34 14.96 11.63 9.16 18.07 8.55 12.57 果球纵径 Fruit longitudinal diameter 20.25 8.45 12.62 16.72 11.05 16.49 14.79 14.34 果球横纵比 Fruit aspect ratio 11.71 4.73 9.25 10.52 9.26 10.85 15.01 10.19 种子千粒质量 Seed 1 000-grain weight 3.00 3.52 3.52 13.03 4.35 15.45 4.94 6.83 平均值 Mean 18.62 10.94 14.06 29.99 14.72 19.57 16.12 17.72 1) JXXY:江西新余;ZJLS:浙江丽水;SCYA:四川雅安;GDQY:广东清远;FJSM:福建三明;HNLD:湖南娄底;HBYC:湖北宜昌
1) JXXY: Xinyu in Jiangxi;ZJLS: Lishui in Zhejiang;SCYA: Ya’an in Sichuan;GDQY: Qingyuan in Guangdong;FJSM: Sanming in Fujian;HNLD: Loudi in Hunan;HBYC: Yichang in Hubei表 6 野生华重楼17个表型性状与生态地理因子的相关系数1)
Table 6 Correlation coefficients between 17 phenotypic traits of Paris polyphylla var. chinensis and ecological geography factors
性状
Trait经度
Longitude纬度
Latitude海拔
Altitude年均温度
Mean annual temperature无霜期
Frost-free period年日照时数
Annual sunshine
hour年降水量
Annual precipitation株高 Height −0.123 0.085 0.671** −0.353 −0.072 −0.109 0.135 地径 Diameter −0.213 0.172 0.530* −0.293 0.135 −0.267 0.157 叶片数 No. of leaf −0.134 −0.074 0.286 −0.077 −0.041 0.007 0.088 叶长 Leaf length −0.206* 0.072 0.354** 0.010 0.200* −0.288** 0.188* 叶宽 Leaf width −0.222** 0.444** 0.530** −0.365** 0.373** 0.535** 0.133 叶面积 Leaf area −0.277** 0.301** 0.512** −0.179* 0.333** −0.425** 0.154 叶片长宽比 Blade aspect ratio 0.010 −0.340** −0.114 0.284** −0.219* 0.324** 0.008 花萼片数 No. of calyx 0.078 −0.230 0.299 −0.317 −0.344 0.152 0.151 花萼长 Calyx length −0.108 0.267** 0.449** −0.291** 0.199* −0.372** 0.210* 花萼宽 Calyx width −0.075 0.199* 0.314** 0.256** 0.226* −0.385** 0.255** 花萼面积 Calyx area −0.136 0.343** 0.386** −0.268** 0.287** −0.435** 0.170 花萼长宽比 Calyx aspect ratio 0.022 −0.186 0.229* −0.011 −0.176 0.118 0.146 花梗长 Pedicel length −0.071 −0.274 0.312 −0.114 −0.377 0.311 −0.036 果球横径 Fruit transverse diameter 0.059 0.081 −0.176 −0.230* 0.037 −0.063 −0.053 果球纵径 Fruit longitudinal diameter −0.027 0.328** 0.115 −0.362** 0.390** −0.495** 0.266* 果球横纵比 Fruit transverse diameter to longitudinal diameter ratio 0.098 −0.337** −0.317** 0.210 −0.490** 0.594** −0.423** 种子千粒质量 Seed 1 000-grain weight 0.232 −0.272 −0.535** 0.306 0.178 −0.135 0.498* 1)“*”和“**”分别表示在0.05和0.01水平(双侧)上显著相关
1)“*”and“**”indicated significant correlations at 0.05 and 0.01 levels (double-tail), respectively -
[1] 汪洋, 陈文学, 明安觉, 等. 湖北红椿天然种群小叶表型性状变异研究[J]. 植物资源与环境学报, 2019, 28(2): 96-105. doi: 10.3969/j.issn.1674-7895.2019.02.13 [2] 周新华, 曾平生, 武晓玉, 等. 野生多花黄精种群表型变异研究[J]. 中南林业科技大学学报, 2019, 39(6): 99-106. [3] 李洪果, 陈达镇, 许靖诗, 等. 濒危植物格木天然种群的表型多样性及变异[J]. 林业科学, 2019, 55(4): 69-83. doi: 10.11707/j.1001-7488.20190408 [4] 郑孙元, 朱弘, 金思雨, 等. 桂花表型变化的环境依赖特征[J]. 南京林业大学学报(自然科学版), 2019, 43(2): 38-46. [5] 李因刚, 柳新红, 马俊伟, 等. 浙江楠种群表型变异研究[J]. 植物生态学报, 2014, 38(12): 1315-1324. [6] 罗建勋, 李晓清, 孙鹏, 等. 云杉天然群体的表型变异[J]东北林业大学学报, 2003, 31(1): 9-11. [7] 彭丽平, 成仿云, 钟原, 等. 凤丹栽培群体的表型变异研究[J]. 植物科学学报, 2018, 36(2): 170-180. doi: 10.11913/PSJ.2095-0837.2018.20170 [8] 赖猛, 胡冬南, 易敏, 等. 晚松人工林株间株内球果及种子形态变异研究[J]. 中南林业科技大学学报, 2018, 38(1): 22-28. [9] 冯延芝, 乔杰, 王保平, 等. 南方低山丘陵区泡桐无性系主要性状的综合选择[J]. 林业科学研究, 2017, 30(6): 969-976. [10] 尹明宇, 朱绪春, 刘慧敏, 等. 西伯利亚杏种质资源花表型变异[J]. 西北农林科技大学学报(自然科学版), 2018, 46(2): 92-103. [11] 柳江群, 尹明宇, 左丝雨, 等. 长柄扁桃天然种群表型变异[J]. 植物生态学报, 2017, 41(10): 1091-1102. [12] 钱磊, 吴忠锋, 唐昌亮, 等. 16种羊蹄甲属植物种子表型变异初探[J]. 热带农业科学, 2017, 37(9): 45-51. [13] 李珊, 甘小洪, 憨宏艳, 等. 濒危植物水青树叶的表型性状变异[J]. 林业科学研究, 2016, 29(5): 687-697. doi: 10.3969/j.issn.1001-1498.2016.05.010 [14] 杨旭, 杨志玲, 王洁, 等. 濒危植物凹叶厚朴种实特性[J]. 生态学杂志, 2012, 31(5): 1077-1081. [15] 沈涛, 申仕康, 张霁, 等. 三七表型变异及其对气候因子的响应[J]. 热带亚热带植物学报, 2017, 25(5): 445-455. doi: 10.11926/jtsb.3724 [16] 国家药典委员会. 中国药典: 一部[M]. 北京: 化学工业出版社, 2020: 271-272. [17] MAN S L, GAO W Y, ZHANG Y J, et al. Antitumor and antimetastatic activities of rhizoma paridis saponins[J]. Steroids, 2009, 74(13/14): 1051-1056.
[18] 尹显梅, 张开元, 蒋桂华, 等. 华重楼皂苷类成分的动态分布规律对药材质量的影响[J]. 中草药, 2017, 48(6): 1199-1204. [19] 刘琴, 樊小莉, 杨莉, 等. 四川盆周山地重楼属植物药用资源的筛选及评价[J]. 中草药, 2019, 50(2): 526-534. doi: 10.7501/j.issn.0253-2670.2019.02.035 [20] 中国科学院中国植物志编委会. 中国植物志: 第15卷[M]. 北京: 科学出版社, 1978. [21] 赵猛, 李燕敏, 王鹏飞, 等. 华重楼地上部分化学成分研究[J]. 中国药学杂志, 2018, 53(16): 1342-1346. [22] 陈铁柱, 文飞燕, 张涛, 等. 华重楼皂苷类化学成分与生态因子的相关性[J]. 中国实验方剂学杂志, 2018, 24(9): 46-51. [23] 周玲, 姚振, 郭永兵, 等. 重楼属植物组织培养研究进展[J]. 北方园艺, 2017(20): 184-189. [24] 王跃华, 蒋婷婷, 付伟, 等. 重楼植物的组织培养研究[J]. 时珍国医国药, 2012, 23(8): 2014-2016. doi: 10.3969/j.issn.1008-0805.2012.08.076 [25] 段宝忠, 马维思, 刘玉雨, 等. 滇重楼无公害栽培关键技术[J]. 世界中医药, 2018, 13(12): 2975-2979. doi: 10.3969/j.issn.1673-7202.2018.12.006 [26] 王林娜, 胡培, 杨光义, 等. 6种重楼属植物不同部位中9种皂苷类成分含量比较[J]. 时珍国医国药, 2018, 29(11): 2635-2638. [27] 金琳, 吴钰颖, 戴雪雯, 等. 基于多指标成分含量测定及HPLC指纹图谱的多茎重楼品质评价[J]. 中草药, 2019, 50(13): 3178-3186. doi: 10.7501/j.issn.0253-2670.2019.13.026 [28] 樊艳荣, 陈双林, 杨清平, 等. 毛竹林下多花黄精种群生长和生物量分配的立竹密度效应[J]. 浙江农林大学学报, 2013, 30(2): 199-205. doi: 10.11833/j.issn.2095-0756.2013.02.007 [29] 李翔国, 全炳武, 李虎林, 等. 人参皂苷量变异的研究进展[J]. 中草药, 2012, 43(11): 2300-2304. [30] 杨再超, 左经会, 向红, 等. 贵州特有药用植物毕节小檗的形态变异[J]. 热带亚热带植物学报, 2016, 24(1): 29-36. doi: 10.11926/j.issn.1005-3395.2016.01.004 [31] 张君诚, 宋育红, 张杭颖. 福建三明珍贵药材八角莲的群落结构及物种多样性[J]. 贵州农业科学, 2014, 42(12): 48-52. doi: 10.3969/j.issn.1001-3601.2014.12.011 -
期刊类型引用(1)
1. 刘蕾,杨易,徐金瑞. 结核分枝杆菌感染对NLRP3炎性小体活化的调控作用研究进展. 微生物学杂志. 2021(04): 91-97 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 856
- HTML全文浏览量: 1
- PDF下载量: 1332
- 被引次数: 1