使用DeepRF/data_synthetic.py的Noise类触发AttributeError: __exit__错误

我尝试使用DeepRF项目中的data_synthetic.py脚本，调用其中的Noise类时出现了异常，以下是最小可复现代码：

最小可复现代码（MWE） #

import numpy as np
import sys
import os

# 导入DeepRF模块
deeprf_module_path = os.path.abspath(os.path.join(os.path.dirname(__file__), '../DeepRF-main'))
sys.path.append(deeprf_module_path)
from data_synthetic import *
import data_synthetic as deeprf_data_synthetic

# 测试DeepRF
TR = 1.2  # 秒
low_frequency_noise = deeprf_data_synthetic.LowFrequency(0.9, TR )
physiological_noise = deeprf_data_synthetic.Physiological(0.9, TR)
# 此处传入RandomState类而非实例，是错误根源
system_noise = deeprf_data_synthetic.System(0.9, np.random.RandomState)
task_noise = deeprf_data_synthetic.Task(0.9, np.random.RandomState)
temporal_noise = deeprf_data_synthetic.Temporal(0.9, np.random.RandomState)

test_t = np.arange(0, 300, 1)
test_signal = np.ones_like(test_t)

# 创建Noise实例，指定噪声组件和幅度
random_generator_t = np.random.RandomState(64556) # 用于生成随机参数
random_generator_x = np.random.RandomState(23355) # 用于生成随机信号/噪声
random_generator_y = np.random.RandomState(1258566) # 用于生成预测值
noise = deeprf_data_synthetic.Noise(random_generator_y.rand(5), low_frequency_noise, physiological_noise, system_noise, task_noise, temporal_noise)

# 生成模拟fMRI时间序列（替换为你的实际数据）
fmri_data = np.random.normal(0, 1, size=300)

# 调用__call__方法给fMRI数据添加噪声
noisy_fmri_data = noise(fmri_data)

问题现象 #

执行最后一行noisy_fmri_data = noise(fmri_data)时，触发如下异常：

Exception has occurred: AttributeError: exit

问题原因 #

查看data_synthetic.py的第144行（Noise类的__call__方法），代码使用了with语句，但你传入的System、Task、Temporal这几个噪声组件的构造参数错误：这些类需要接收np.random.RandomState的实例，但你传入的是np.random.RandomState类本身，导致后续代码尝试将类作为上下文管理器使用，而类不具备__exit__方法，因此抛出异常。

修复方案 #

将创建System、Task、Temporal实例时的第二个参数改为RandomState的实例，可使用你已定义好的random_generator_x，或者新创建实例：

修复后的代码片段：

# 使用已有的RandomState实例
system_noise = deeprf_data_synthetic.System(0.9, random_generator_x)
task_noise = deeprf_data_synthetic.Task(0.9, random_generator_x)
temporal_noise = deeprf_data_synthetic.Temporal(0.9, random_generator_x)

# 或者创建新的实例
# system_noise = deeprf_data_synthetic.System(0.9, np.random.RandomState(seed=123))
# task_noise = deeprf_data_synthetic.Task(0.9, np.random.RandomState(seed=456))
# temporal_noise = deeprf_data_synthetic.Temporal(0.9, np.random.RandomState(seed=789))

替换原代码中对应的三行，再运行即可解决该异常。

内容的提问来源于stack exchange，提问作者skm