问题标签 [neuroscience]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
756 浏览

python - 如何使用 NiBabel (Python) 预处理 NIfTI 数据格式

使用 NiBabel 将 NIfTI 文件转换为数组后,数组具有三个维度,数字如下所示:

如何预处理此数组以进行机器学习?当只选择指数时,大部分信息在绘制图像时都会丢失,所以也许基数也很重要?

任何帮助表示赞赏。

0 投票
1 回答
234 浏览

r - R中EEG数据的巴特沃斯过滤

我对 R 和 EEG 信号非常陌生,所以如果问题的答案很明显,请原谅。

我正在尝试对 EEG 信号执行巴特沃斯滤波器以提取 Alpha 波段。当我执行过滤器时,生成的信号看起来很奇怪,完全不是我所期望的,在时间框架的开始有一个异常大的峰值。我尝试使用 eegfilter 和 bwfilter 来查看代码是否存在问题,但是当我绘制结果时两者之间几乎没有区别。我无法解释最终结果,如果有人可以向我解释特殊的最终结果,我将不胜感激。

这是我正在查看的数据中的一个示例: https ://ufile.io/1ji48wg6

采样率为 512。

我想提取 alpha 波段,所以频率在 8 到 12 Hz 之间

这是绘制时的数据图片。最左边的图像是原始数据。中心图是使用 eegfilter 应用 Butterworth 滤波器的结果。右图是使用 bwfilter 应用巴特沃斯滤波器的结果。

应用过滤器时的数据图

数据集的标题:

脑电图

-8438.876837

-8442.718979

-8441.877183

-8439.974768

-8443.436883

-8448.900711

-8452.433874

-8441.616546

0 投票
1 回答
15 浏览

neuroscience - 如何获取 NeuroML 数据库中可用模型的完整列表?

我想获得 NeuroML-DB.org 中所有细胞模型的列表。如何以编程方式执行此操作?

0 投票
1 回答
66 浏览

python - 在 Shady 中呈现预先准备好的图像数组的一部分

我有兴趣从 psychtoolbox 迁移到shady以进行刺激演示。我查看了在线文档,但我不太清楚如何复制我目前在 matlab 中所做的事情。

我所做的其实很简单。对于每次试验,

  1. 我从磁盘加载单个图像(我离线进行亮度线性化),其中包含我计划在该试验中显示的所有帧(刺激为 1000x1000 像素,我呈现 25 帧,因此图像为 5000x5000 像素。我只使用 BW 图像,所以我每个像素有一个 int8 值)。

  2. 我将整个图像从 CPU 传输到 GPU

  3. 在某些时候(外部控制)我将第一帧复制到视频缓冲区并呈现它

  4. 在其他某个点(外部控制),我触发剩余 24 帧的显示(将图像的相关部分复制到每个视频帧的视频缓冲区,然后调用翻转())。

  5. 外部控制是通过让另一台机器通过 TCP/IP 与刺激呈现代码进行通信来实现的。在控制PC向演示PC发送命令并执行后,演示PC需要向控制PC发送回确认消息。我需要发送三个 ACK​​ 消息,一个当屏幕上出现第一帧时,一个当屏幕上出现第二帧时,一个当屏幕上出现第 25 帧时(这样控制 PC 可以轻松验证是否丢帧)。在 matlab 中,我通过调用阻塞方法 flip() 来呈现一个帧,当它返回时,我将 ACK 发送到控制 PC。

而已。我怎么会在阴暗的地方做到这一点?有我应该看的例子吗?

0 投票
1 回答
16 浏览

neuroscience - 如何使用 NeuroML Database API 收集与给定通道模型关联的单元模型?

我想将 NeuroML 数据库中的通道模型子集与它们嵌入的单元模型相关联。我想知道是否有一种快速简便的方法来获取给定通道模型的所有关联单元。

这是门户网站上的“第一个”频道模型 ( https://neuroml-db.org/model_info?model_id=NMLCH000001 )。它说它属于底部的细胞模型。

我已经使用https://neuroml-db.org/api/models?id=NMLCHxxxxxx收集了所有通道模型的详细信息(“类型”:“CH”),所以我有我想要查看的每个通道模型的 JSON。我找不到“关联单元格”字段。

是否有一个包含所有 NMLCLxxxxxx 相关单元格列表的地方?我可能会错过它。

0 投票
2 回答
131 浏览

python - 将直方图函数扩展到具有任意间隙的重叠箱和箱?

是否有一个直方图函数可以处理任何一般的 bin?

在我的神经科学应用程序中,我有两个一维数组:spikestime_centers和一个参数time_window。我的目标是输出名为 的数组firing_rate,其大小与time_centers定义为相同firing_rate[i]=len(abs(spikes-time_centers[i])<time_window),以便计算time_center宽度间隔中每个箱附近的尖峰数time_window

我很快起草了一个函数,但我有点担心性能,因为我需要处理大量的脉冲序列。

是否可以利用spikes数组已排序的事实?即使我使用相同大小 ( time_window) 的垃圾箱,我的垃圾箱也可能重叠,或者垃圾箱之间可能存在可变间隙。

通常spikes包含约 1000 个元素,time_centers约 50 个箱,我需要为每个会话计算 30000 个神经元。

0 投票
0 回答
30 浏览

python - 绘制 3D 点的刺激类

我是使用 PsychoPy 进行神经科学实验的新手,我想在 3D 平面上画点。我找到了 DotStim 和 ElementArrayStim 函数,但它们只采用 X 和 Y 坐标(第 3 维没有额外的 Z 坐标)。

是否有允许在 PsychoPy 中绘制 3D 点的功能,如果没有,我对其他库开放。请帮忙!谢谢!

0 投票
1 回答
52 浏览

python - 从 AllenSDK ephy 数据库中提取时间序列

我正在寻找 AllenSDK 库。可以提取他们为估计各种参数所做的记录的时间序列(电池电压)吗?我可以找到如何访问参数,但我对它们的踪迹真的很感兴趣。任何人都可以帮忙吗?

0 投票
2 回答
876 浏览

neuroscience - 如何修复未创建 template_4D_GM 文件的 fslvbm_2_template 中的错误

我正在 FSL 中运行 VBM 分析。我在 JISC 邮件档案中搜索了此错误的解决方案,但没有运气。我的笔记本电脑上有很多可用空间,所以这不是空间问题。我在多个 Mac 系统上也遇到过这个问题。任何帮助将不胜感激,我已经为此苦苦挣扎了好几天。

当我运行代码“fslvbm_2_template -n”时,我希望创建一个名为 template_4D_GM 的文件,但在输出中会出现以下错误代码:

运行“fslvbm_2_template -n”时,我发现 fslvbm2a 和 fslvbm2b 将正常运行,但我在 fslvbm2c 遇到了问题。当我从 fslvbm_2_template 运行这段特定的代码时 - FSLVBM 模板创建(http://ftp.nmr.mgh.harvard.edu/pub/dist/freesurfer/tutorial_packages/OSX/fsl_501/bin/fslvbm_2_template):

0 投票
2 回答
189 浏览

python - 文件加载错误:可用 1.7TB 存储空间不足

我正在使用以下代码在 Python 中以 NiFTI 格式加载我的文件。

少量图像工作正常,但如果我想加载更多图像,我会收到以下错误:

我认为 img.uncache() 会从内存中删除图像,因此它不会占用太多存储空间,但仍然可以使用图像数组。将这一位添加到代码中并没有改变任何东西。

有谁知道我能帮上忙吗?我正在使用的计算机有 24 核 2.6 GHz CPU,超过 52 GB 内存,工作目录有超过 1.7 TB 的可用存储空间。我正在尝试从 ADNI 数据库加载大约 1500 张 MRI 图像。

任何建议都非常感谢。