做农产品检测、光谱分析的同行常会遇到一个头疼难题:作物早期病害藏在叶片内部,表面看不出异样,等肉眼能看清病斑时,病菌早已大面积扩散,果园损失直接拉满。
《Chemometrics and Intelligent Laboratory Systems》刊发一篇有意思的研究,把时序分析里的矩阵轮廓(Matrix Profile)算法嫁接到高光谱成像,专门针对苹果黑星病做早期无损筛查,解决传统光谱检测容易误判、病灶对比度弱的痛点。
高光谱成像能同时拿到图像空间信息与数百波段光谱数据,不用破坏样品就能分析作物内部水分、有机质变化,是农业无损检测的常用工具。但放到田间叶片检测场景,短板会被无限放大。
论文里提到,叶片本身结构复杂,叶脉、厚薄、表面粗糙度、光照散射都会产生大面积光谱干扰,形成很强的全局基线偏移。拿苹果黑星病叶片举例,感染初期病斑极小,水分细微变化带来的光谱差异,会被叶片本身的物理噪声完全盖住。
常规处理手段依赖PCA降维、光谱指数或者监督学习,问题很明显:
1. 监督模型需要提前标注大量病叶、健康叶片样本,新病害场景很难快速复用;
2. PCA直接处理原始反射光谱,叶脉、阴影都会被识别成异常,病斑和正常组织区分不开;
3. 全局乘性散射效应干扰严重,轻微光照变化就会打乱整套数据。
简单说:传统算法分不清“叶片天生纹理”和“病菌造成的病变”,误判多,早期微小病灶很难凸显。

这篇论文最巧妙的创新,是把原本用于时间序列挖掘的矩阵轮廓算法,改造成高光谱数据预处理工具,全程无监督、不用提前标注样本,核心逻辑很好理解。
时序分析里,矩阵轮廓会逐段对比序列片段,算出每一小段和全局最相似片段的距离,距离越高,代表这段数据越特殊、越反常。研究团队把二维叶片光谱像素矩阵,转换成一维长序列,把连续波段当成“时序片段”,设置 50 个波段为一组滑动窗口。
每个光谱窗口都会计算和其他窗口的最小欧式距离,生成一套全新的“距离光谱”,也就是论文里的MP距离谱。相当于给原始光谱加了一层筛选器:只保留局部光谱差异,弱化叶片整体散射带来的全局干扰。
1. 采集苹果叶片SWIR短波红外高光谱数据,圈选包含早期病斑的感兴趣区域;
2. 将三维光谱立方体拉伸为一维序列,滑动截取波段子序列,计算矩阵轮廓距离;
3. 生成像素-波段距离矩阵,截取有效波段做二次PCA分析;
4. 计算得分距离SD生成可视化热力图,高亮标记异常病斑像素。
整套流程不需要人工标注病害样本,拿到原始高光谱图像就能直接运算,对田间快速筛查场景十分友好。

研究选取接种后第4天、第5天两组叶片样本做对照,分别用原始反射光谱、矩阵轮廓变换后的距离光谱做PCA分析,两组数据差距一目了然。

直接对原始光谱做PCA得分距离成像,整张叶片布满高亮噪点,叶脉、叶片边缘、轻微阴影全部被判定为异常,细小病斑淹没在杂讯里,人工很难定位病灶。
经过矩阵轮廓变换后再成像,叶脉、叶片纹理带来的干扰信号被大幅压制,整张图只有黑星病病斑区域呈现高亮,无关杂点几乎消失,病灶边界清晰可辨。论文统计,两组样本威尔克斯λ 数值均出现明显下降,类内方差占比缩减近一半,健康叶片与病斑像素区分度显著提升。
论文通过载荷曲线分析得出关键结论:病菌侵染会改变叶片内部水分状态,1400–1450nm、1900–2000nm 两处水吸收峰两侧的光谱斜率变化,是区分健康与染病叶片的核心特征波段。
原始光谱会被整体反射强度掩盖这些细微斜率变化;矩阵轮廓变换放大了波段间局部差异,水吸收带的细微形变会转化为高距离数值,哪怕发病仅2–5天的早期叶片,也能捕捉到微弱病变信号。

同时论文验证,这套算法对不同光谱特征的异常区域兼容度较好,图像配准误差带来的伪异常、真实病斑都能同步识别,不会漏掉不同类型的光谱异常。
抛开论文的实验室场景,这套矩阵轮廓预处理思路,对做果蔬、作物高光谱无损检测的从业者有不少参考价值:
1. 降低样本标注成本:无监督算法,面对新型作物病害、未知缺陷时,不用提前采集大量标注样本,缩短设备调试周期;
2. 适配复杂田间样本:叶片、果皮天然存在纹理、厚度差异,算法可弱化物理结构干扰,提升野外、大棚现场检测稳定性;
3. 轻量化适配成像设备:矩阵轮廓计算逻辑简洁,可集成到小型高光谱成像设备、在线分选产线,实现实时预处理成像;
4. 拓展多品类检测场景:除果树叶片病害,同理可用于果实内部霉变、农产品隐性损伤、工业材料微小缺陷筛查。
目前该研究仅以苹果黑星病作为概念验证,后续还有参数灵敏度、多作物数据集拓展的研究空间,但足以证明:矩阵轮廓是一套适配高光谱数据、值得落地尝试的预处理方案。
免责声明
本内容为公开学术文献链接:https://www.sciencedirect.com/science/article/pii/S0169743926001292 整理编撰,仅用于行业技术探讨与科普学习,不作任何商业相关承诺,亦不能作为投资参考依据。文中所列各项实验数据与结论会受试验环境、样本个体区别、模型构建方案等多重变量干扰,若落地实际场景使用,相关效果需结合自身场景另行实测核验。




销售直拨:13858065387
固定电话:13858065387