西葫芦叶片在不同施肥处理下叶绿素荧光与色素测定的像素级高光谱(pixel-level hyperspectral)指标研究

《Journal of Agriculture and Food Research》：Pixel-level hyperspectral metrics for chlorophyll fluorescence and pigment assays in zucchini leaves across fertilisation treatments

　　

摘要：植物胁迫可在肉眼可见前即已发生，而堆肥(compost)等肥料广泛用于育苗基质时常促进生长但亦可能引发生理异常。在此背景下，基于高光谱成像(hyperspectral imaging, HSI)的无损检测指标可用于捕捉叶片变化并与破坏性分析相关联。研究人员评估了一系列HSI指数及光谱描述子能否捕获堆肥剂量依赖性响应：选用21株西葫芦(Cucurbita pepo L.)植株，设三种堆肥体积比处理(0%、5%、10% v:v)，采集叶片高光谱影像并验证其与实验室破坏性检测结果的一致性。通过区域生长算法(region-growing algorithm)在像素级提取均匀叶面积，每一样本获取4个感兴趣区(Region of Interest, ROI)，进行去噪与偏差校正后，计算文献中已有光谱指数及自定义度量，并与各处理下实验室测定值做Spearman秩相关分析。结果表明光谱响应呈显著非线性——5%堆肥处理下多数光谱指标指示生理状态改善，且与光化学及色素实验室测定值具一致关联；基于可见光波段光谱积分(visible-range spectral integrals)的指标对区分各堆肥剂量下生理模式最具判别力(ρ = +0.536 于 C0，ρ = ?0.286 于 C5，ρ = +0.643 于 C10)，绿光波段(green region)指标亦然(ρ = +0.714 于 C0，ρ = ?0.214 于 C5，ρ = +0.857 于 C10)。此外，10%堆肥处理下若干描述子量值偏移，与实验室识别的胁迫症状或不良叶况相符，类胡萝卜素/叶绿素吸收窗基于连续统去除(continuum-removal, CR)的指标ρ值最高达0.714。上述发现支持利用HSI追踪叶片堆肥剂量—响应关系而无需破坏性采样，为养分管理及堆肥优化提供依据，并可推广至其他作物。

论文解读：西葫芦叶片堆肥梯度下像素级高光谱指标与叶绿素荧光及色素测定的关联研究

《Journal of Agriculture and Food Research》刊载的此研究针对当前园艺生产中堆肥部分替代泥炭育苗基质日益普遍，但过量堆肥可能引发盐分胁迫、而植物早期生理变化肉眼不可见的问题，指出现有高光谱数据分析多依赖整叶平均或理论辐射传输模型（如PROSPECT模型），忽略叶片内部异质性（如叶脉与叶缘区色素、厚度差异），且线性模型难以捕捉施肥剂量引起的非线性 hormetic（毒物兴奋）效应。研究人员以西葫芦大叶为模型（具明显中央—边缘异质性），结合像素级ROI提取、多类光谱指数（波段比、连续统去除、红边与绿光形态、微分微积分描述子）与叶绿素荧光OJIP测试及丙酮萃取色素测定，探究HSI能否无损检测堆肥0%、5%、10%(v:v)处理下的叶片生理剂量—响应关系，最终证实5%堆肥出现非线性最优（甜点）响应，可见光与绿光波段积分、红边斜率、色素窗连续统去除指标最具诊断价值，支持HSI用于苗期基质优化的无损筛查。

主要关键技术方法：研究人员设西葫芦(Cucurbita pepo L. cv. Tevere)幼苗三组堆肥替代泥炭基质处理——C0(0% compost, 80% peat + 20% perlite)、C5(5% compost, 75% peat + 20% perlite)、C10(10% compost, 70% peat + 20% perlite)，每处理7株共21株。用Specim IQ?便携式高光谱成像系统(400–1000 nm, 光谱分辨率7 nm)采集完全展开叶影像，经暗/白参考校正得反射率R(λ)；以Python编写带GUI的交互式区域生长(region-growing)算法在550 nm附近波段手动选种，于每叶取4个均质无阴影ROI(每ROI 3000–7000像素)，Madden平滑(窗口11波段)去噪后算各ROI平均反射谱。提取光谱指数包括：三波段指数(Three-Band Index, 3BI)、归一化差红边指数(Normalized Difference Red Edge Index, NDRE)、叶绿素指数(Chlorophyll Index, CI)、近红外/红边比值(NIR/RE ratio)；连续统去除(Continuum Removal, CR)于叶绿素吸收窗(640–680 nm)及类胡萝卜素吸收窗(470–500 nm)求带深(D_b)、吸收面积、偏度(chl_depth, chl_area, chl_skew; car_depth, car_area, car_skew)；红边斜率(701–722 nm, m_dyn, q_dyn)、红边位置(Red Edge Position, REP)及曲率指数(Curvature Index, CI)、绿光形态不对称(Green Shape Asymmetry)；可见光积分(VIS Integral, ∫₄₀₀⁷⁰⁰R(λ)dλ)、525 nm至红边积分及绝对导数和(Absolute Derivative Sum Index, ADS_sum)。同期用Handy PEA荧光仪测暗适应叶片快相荧光(OJIP)得F_v/F_o、F_v/F_m、Φ_Po、Φ_Eo、Φ_Do、Φ_Pav、ABS/RC、TR_o/RC、ET_o/RC、DI_o/RC等；用80%丙酮萃取—分光光度法(Lichtenthaler方程)测叶绿素a(Chl a)、叶绿素b(Chl b)及类胡萝卜素(Carotenoids, Car)。HSI各指标与实验室参数按处理分组做Spearman秩相关，置换检验(10,000次)求p值，Benjamini–Hochberg法控制错误发现率(False Discovery Rate, FDR, q < 0.05)。

研究结果：

3.1. Laboratory analyses（实验室分析）

C0→C5→C10处理下，最大PSII光化学效率(F_v/F_m、Φ_Po)微升，电子传递量子产额(Φ_Eo)略增，热耗散产额(Φ_Do)下降；单位反应中心吸收通量(ABS/RC)、捕获通量(TR_o/RC)呈递减趋势， dissipated能量(DI_o/RC)降低；色素含量(Chl a、Chl b、Car)在各重复间波动无严格单调，显示处理间存在生理差异但非简单线性变化。

3.2. Hyperspectral metric data（高光谱指标数据）

平均反射光谱显示C5整体反射率最低（吸收更强，暗示更高色素或更致密叶结构），C0与C10较高；一阶导数曲线差异较小，说明变化弥散分布于波段而非突变。可见光积分(VIS Integral)与绿光积分(Green Integral)随堆肥剂量呈C0中等→C5最小→C10回升趋势，红边起始波长(λ_start_lin)与终止波长(λ_end_lin)在C5略提前，印证C5具更锐利红边过渡。

3.3. Data correlations（数据相关性）

3.3.1. Linear correlations（线性相关）

少数指标（红边线性插值m_dyn、q_dyn）与实验室参数呈单调弱相关；VIS Integral与Green Integral对F_v/F_m、Φ_Po、F_v/F_o在C0为正相关，C5转负相关，C10再转正——体现非线性；NIR/RE ratio与Chl b相关性由C0弱负转C10强正；CR指标与Chl a相关性C0→C10递减。线性相关性总体偏弱且仅部分指标显著。

3.3.2. Non-linear correlations（非线性相关）

大多数HSI指标—实验室参数对呈显著非线性（U型或倒U型），C5常现相关性峰值或符号翻转，空白格多为已报线性相关而非无相关。积分类、CR类、红边动态类指标整体在C5表现出与荧光及色素参数最强关联。

3.3.3. Integral-based light-partitioning metrics（基于积分的光分配指标）

VIS Integral对F_v/F_o、F_v/F_m、Φ_Po的Spearman ρ为C0 = +0.536、C5 = ?0.286、C10 = +0.643；Green Integral(520–600 nm)对同组参数ρ为C0 = +0.714、C5 = ?0.214、C10 = +0.857；525 nm至红边积分(∫₅₂₅^~RER(λ)dλ)对F_v/F_o等ρ为C0 = +0.679、C5 = ?0.250、C10 = +0.643。说明C5时光合光化学较优（反射低、吸收高致积分偏小负相关），过高或零堆肥偏离该最佳态。

3.3.4. Continuum-removal metrics on pigment windows（色素窗连续统去除指标）

类胡萝卜素窗CR指标(car_area、car_depth、car_skew)与Chl a、DI_o/RC、TR_o/RC在C5呈近零或弱相关，C0与C10呈负或正中等相关；叶绿素窗CR指标(chl_area、chl_depth、chl_skew)对TR_o/RC在C0 ρ = +0.786(p = 0.036)，C5转负，C10减弱；对Chl a在C10 ρ = ?0.857(p = 0.014, FDR q = 0.240)。表明CR形状描述子可捕捉C5最佳色素—结构配置及C10潜在胁迫信号。

3.3.5. Red-edge dynamics（红边动态）

701–722 nm红边斜率(m_dyn)对TR_o/RC、ABS/RC、DI_o/RC在C5呈极值趋势（C0近0，C5轻微负，C10负向增强），佐证C5处红边陡峭度最大，对应较高叶绿素含量及有序叶肉散射结构。

讨论与结论翻译：

研究表明5%堆肥输入使西葫芦叶片产生独特的高光谱—生理特征，表现为HSI指标与实验室生理参数相关性增强及相关系数符号翻转，反映接近最优生理运作状态，可解释为hormetic效应（低剂量堆肥生物刺激，高剂量致盐分/养分失衡早期胁迫）。非线性剂量—响应中C5峰值多经FDR校正仍显著，可见形状类指标(波段积分、比值、偏度、红边斜率)比单波段指数重复性与生物学对应性更好。像素级区域生长ROI流程在均质叶区内稳健，排除叶脉与镜面反射干扰。局限性含需人工初选种子点，未来应发展全自动ROI与高通量表型适配。结论：HSI与荧光参数的关联可有效检测西葫芦叶片堆肥处理下的胁迫与最适生理状态；多数相关呈非单调剂量响应并在5%堆肥(C5)出现明确最优；可见光、绿光及红边区域光谱描述子（尤其是VIS/Green Integral、Green/VIS ratio、色素窗CR指标chl_area/chl_depth/chl_skew与car_area/car_depth、701–722 nm红边斜率）与叶绿素荧光及实验室色素测定吻合最佳，C0与C10偏离此态且C10显现胁迫征兆；基于形态学处理的像素级ROI流程提升跨处理可重复性，为养分管理与堆肥优化提供无损快速工具。