MNE-Python中与Scipy 1.15.0兼容性问题分析

在MNE-Python 1.9.0版本中，当与Scipy 1.15.0版本一起使用时，出现了与球谐函数相关的测试失败问题。这个问题主要影响到了两个关键测试模块：预处理模块中的maxwell测试和转换模块中的transforms测试。

问题背景

Scipy 1.15.0版本中对sph_harm函数进行了重大变更，将其标记为已弃用(deprecated)，并引入了新的替代函数sph_harm_y。这一变更导致了MNE-Python中依赖球谐函数计算的部分功能出现了兼容性问题。

具体表现

测试失败主要表现在两个方面：

1. 预处理模块测试失败：在测试球面坐标转换时，系统抛出了SystemError异常，提示DeprecationWarning返回了一个带有异常设置的结果。

2. 转换模块测试失败：在测试薄板样条(Thin Plate Spline, TPS)变换时，同样出现了类似的SystemError异常。

技术分析

深入分析这些问题，我们可以发现：

1. 函数替代差异：当尝试将sph_harm替换为新的sph_harm_y函数时，计算结果出现了显著差异。测试数据显示，在预处理模块中，96.7%的数据点不匹配，最大绝对差异达到0.3449；在转换模块中，7%的数据点不匹配，最大绝对差异为0.00437。

2. 数值计算差异：新旧函数在计算球谐函数时采用了不同的算法或约定，导致计算结果存在系统性差异。这种差异不仅体现在数值大小上，还体现在结果的符号和相位等方面。

解决方案

MNE-Python开发团队已经在主分支(main)中修复了这个问题。修复方案可能包括：

1. 函数替换：将弃用的sph_harm函数替换为新的sph_harm_y函数。

2. 结果调整：根据新函数的输出特性，对计算结果进行适当的后处理或调整，以保持与旧版本行为的一致性。

3. 兼容层：可能实现了某种兼容层，使得在不同Scipy版本下都能获得一致的结果。

影响范围

这个问题主要影响：

1. 使用MNE-Python进行脑磁图(MEG)数据预处理的用户
2. 依赖球谐函数进行空间变换的用户
3. 使用Scipy 1.15.0或更新版本的用户

用户建议

对于遇到此问题的用户，建议：

1. 升级到包含修复的MNE-Python版本
2. 如果暂时无法升级，可以考虑降级Scipy版本
3. 关注MNE-Python的更新日志，了解相关修复的详细信息

这个问题很好地展示了科学计算生态系统中库之间相互依赖的复杂性，也提醒开发者在依赖其他库的核心功能时需要做好版本兼容性管理。