基于猕猴感觉运动皮层Spike信号的运动解码分析不同运动参数对解码的影响

news2025/1/17 18:06:58
  1. 公开数据集中文版详细描述参考前文:https://editor.csdn.net/md/?not_checkout=1&spm=1011.2124.3001.6192
  2. 神经元Spike信号分析参考前文:https://blog.csdn.net/qq_43811536/article/details/134359566?spm=1001.2014.3001.5501
  3. 神经元运动调制分析参考前文:https://blog.csdn.net/qq_43811536/article/details/134401004?spm=1001.2014.3001.5501
  4. 离散卡尔曼滤波器算法详解及重要参数(Q、R、P)基于公开数据集解码的讨论参考前文:https://blog.csdn.net/qq_43811536/article/details/134403022?spm=1001.2014.3001.5501

目录

  • 0. 公开数据集
  • 1. 运动参数对解码结果的影响
    • 1.1 基于Spike的运动解码
    • 1.2 特征提取
    • 1.3 解码算法
    • 1.3 不同运动参数的影响


0. 公开数据集

  • 网址:Nonhuman Primate Reaching with Multichannel Sensorimotor Cortex Electrophysiology 1
  • Session:
    • " indy_20170124_01 "

1. 运动参数对解码结果的影响

从 公开数据集中文版详细描述 以及其 神经元Spike信号分析 中我们知道:猕猴的感觉运动皮层对运动的调制主要有六个参数——二维平面中x和y方向上的位移、速度和加速度。这六个参数也是我们在运动解码中感兴趣的部分。

1.1 基于Spike的运动解码

首先,什么是基于Spike的运动解码?运动解码到底是想做什么?

为了回答这个问题,让我们再次回顾 公开数据集中文版详细描述 一文,已知的数据是猕猴控制光标(cursor)追逐目标(target)时光标和猕猴手指(finger)的运动轨迹以及同步记录下来的颅内脑电信号

  • 运动轨迹:这部分数据是250Hz采样率的时间点和当前时刻的cursor和finger所在的位置,前者为二维数据,后者是三维数据(z方向固定),所以用任何一部分位置数据均可。
  • 颅内脑电:颅内记录的是同等采样率下不同神经元Spike发放的时间点,即已经sort好的Spike数据(很友好,不需要我们再动手sort)。

既然二者数据是在时间上一一对应的,那么最直观的运动解码旨在希望仅通过记录的颅内脑电信号即可得到精确的运动轨迹,也就类似于我们常说的“读心术”。

那么运动解码最重要且必要的实际意义就在于:恢复神经功能尚且完好的瘫痪患者的运动能力

1.2 特征提取

基于 神经元Spike信号分析 一文,我们主要将Spike特征运动特征定义如下:

  • Spike特征:spike_feat = [time_bins, n_neuron]。时间箱的大小可以根据实际解码意义和结果进行选取,我这里取100ms;神经元个数取符合发放率(比如大于1Hz)的即可。
  • 运动特征:motion_feat = [time_bins, n_motion]。时间箱同上,运动参数即二维平面中x和y方向上的位移、速度和加速度

1.3 解码算法

我们这里就以 离散卡尔曼滤波器算法详解及重要参数(Q、R、P)的讨论 中的卡尔曼滤波器作为解码算法。

1.3 不同运动参数的影响

我们分别讨论了以下几种情况的运动参数对最终解码结果的影响:

  • 只有位置/速度/加速度;
  • 同时使用位置、速度;
  • 同时使用位置、加速度;
  • 同时使用速度、加速度;
  • 全部考虑,即同时使用位置、速度、加速度。

在合适的Q、R、P初始值设置下(参考 重要参数(Q、R、P)的讨论 一文),我们使用卡尔曼滤波器分别对上述几种情况进行解码评估,使用十折交叉验证。我们对7种情况的解码结果的运动参数(特征)维度进行平均,然后对十折实验结果绘制了箱型图,如 Figure 1 所示:

在这里插入图片描述

Figure 1. 不同运动参数下卡尔曼滤波器的解码性能

Figure 1 的横轴分别为7种运动参数情况,纵轴为Pearson相关性。可以看到加速度对解码结果的贡献最大,仅加速度特征存在时就可以达到0.4以上的解码相关性,这和预期是一致的:因为加速度是通过位置数据的二阶差分得到,包含了最丰富的信息。而仅速度存在的情况下解码相关性最低,甚至低于仅有位置特征的情况,这与我们的预期相反,理论上速度特征包含了部分位置特征信息,预期应该介于位置和加速度解码性能之间,我们认为可能的原因是速度作为一种关键特征只有和其他特征同时存在时才能最大化预测效果,仅速度特征的解码性能最低无法说明速度不重要,因为这无法解释速度+加速度大于位置+加速度的结果,这反而说明速度应该是一个极为核心的参数,但可能无法单独发挥作用。


创作不易,麻烦点点赞和关注咯!

  1. Makin, J. G., O’Doherty, J. E., Cardoso, M. M. B. & Sabes, P. N. (2018). Superior arm-movement decoding from cortex with a new, unsupervised-learning algorithm. J Neural Eng. 15(2): 026010. doi:10.1088/1741-2552/aa9e95 ↩︎

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1225092.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

2D槽道流

之前看槽道流时,一直无法在二维槽道流里计算出湍流状态,后来了解到二维槽道流需要额外添加随机扰动,但是这个体积力的植入方式一直不知道。而且看稳定性分析中的OS方程的推导,也是基于2d的NS方程,至今还是很疑惑这个问…

Struts2 数据校验之四兄弟

现在是科技的时代,大多数人都在网上购物了, 我们都碰到过相同的问题,各大网站弄的那些各种各样的注册页面,相信大家都深有体会。 有了这验证就很好的保证了我们的信息的准确性和安全性。 接下来我给大家讲解一下用struts2怎么实…

H5ke11--3介绍本地,会话存储

代码顺序: 1.设置input,捕获input如果有多个用属性选择符例如 input[typefile]点击事件.向我们的本地存储设置键值对 2.在点击事件外面设置本地存储表示初始化的值.点击上面的事件才能修改我们想修改的值 会话(session)浏览a数据可以写到本地硬盘,关闭页面数据就没了 本地(…

2023年【P气瓶充装】报名考试及P气瓶充装复审考试

题库来源:安全生产模拟考试一点通公众号小程序 2023年P气瓶充装报名考试为正在备考P气瓶充装操作证的学员准备的理论考试专题,每个月更新的P气瓶充装复审考试祝您顺利通过P气瓶充装考试。 1、【多选题】充装过程中出现充气头漏气的主要原因为&#xff1…

Java Web——JavaScript运算符与流程语句

1. 运算符 1.1. 算数运算符 数字是用来计算的,比如:乘法 * 、除法 / 、加法 、减法 - 等等,所以经常和算术运算符一起。 算术运算符:也叫数学运算符,主要包括加、减、乘、除、取余(求模)等 …

git拉取普通idea Java项目module没有build的问题

在不断完成一个项目的时候,会有不断新加的module,我们用git拉取时会发生没有识别新module的情况。 解决方法是右键项目名称,然后点击Open Module Settings 接下来,点击Module,加号,新建Module的名字就是在g…

高效文件管理:一键批量修改文件名,并统一转换为大写扩展名

在日常生活和工作中,文件处理成为了一项必不可少的任务。无论是个人还是企业,都需要管理大量的文件,包括图片、文档、音频和视频等。这些文件的名字可能千奇百怪,格式各不相同,而且往往需要按照一定的规则进行修改或整…

MFC 对话框

目录 一、对话款基本认识 二、对话框项目创建 三、控件操作 四、对话框创建和显示 模态对话框 非模态对话框 五、动态创建按钮 六、访问控件 控件添加控制变量 访问对话框 操作对话框 SendMessage() 七、对话框伸缩功能实现 八、对话框小项目-逃跑按钮 九、小项…

十一、统一网关GateWay(搭建网关、过滤器、跨越解决)

目录 一、网关技术的实现 在SpringCloud中网关的实现包括两种: 作用: 二、搭建网关服务 1、新建模块,并添加依赖 2、新建Gateway包,并编写启动类 3、编写yml文件 4、启动服务,并在网页内测试 5、步骤 三、路由断言工厂 …

Vue3 shallowRef 和 shallowReactive

一、shallowRef 使用shallowRef之前需要进行引入: import { shallowRef } from vue; 使用方法和ref 的使用方法一致,以下是二者的区别: 1. 如果ref 和 shallowRef 都传入的是普通数据类型的数据,那么他们的效果是一样的&#x…

OpenGL 的学习之路-4(变换)

三大变换:平移、缩放、旋转(通过这三种变换,可以将图像移动到任意位置) 其实,这背后对应的数学在 闫令琪 图形学课程 中有过一些了解,所以,理解起来也不觉得很困难。看程序吧。 1.画三角形&am…

【C++】类和对象(6)--运算符重载

目录 一 概念 二 运算符重载的实现 三 关于时间的所有运算符重载 四 默认赋值运算符 五 const取地址操作符重载 一 概念 C为了增强代码的可读性引入了运算符重载,运算符重载是具有特殊函数名的函数,也具有其返回值类型,函数名字以及参数…

【使用vscode在线web搭建开发环境--code-server搭建】

官方版本下载 https://github.com/coder/code-server/releases?q4.0.0&expandedtrue使用大于版本3.8.0,因为旧版本有插件市场不能访问的情况版本太高需要更新环境依赖 拉取安装包 []# wget "https://github.com/coder/code-server/releases/download/v4.0.0/code-…

[游戏开发][Untiy]跨平台可视化Log系统

工具介绍 今天介绍的主角是LogViewer 工具运行时长这个样子,Unity的Log日志都会在这里显示 如何安装 在Unity商店搜索Log,排名第一的就是它 也可以去Github官网下载源码: Unity-Logs-Viewerhttps://github.com/aliessmael/Unity-Logs-Vie…

机器学习的医疗乳腺癌数据的乳腺癌疾病预测

项目视频讲解:基于机器学习的医疗乳腺癌数据的乳腺癌疾病预测 完整代码数据分享_哔哩哔哩_bilibili 效果演示: 代码: #第一步!导入我们需要的工具 import numpy as np import pandas as pd import matplotlib.pyplot as plt import seaborn as sns %matplotlib inlin…

Python的数据分析包Pandas?示例文章完成版来啦~

文章目录 前言一、Pandas简介二、Python Pandas的使用 总结 前言 pandas 是基于NumPy 的一种工具,该工具是为了解决数据分析任务而创建的。 Pandas 是python的一个数据分析包,最初由AQR Capital Management于2008年4月开发,并于2009年底开源…

Javaweb之Vue指令案例的详细解析

2.3.5 案例 需求: 如上图所示,我们提供好了数据模型,users是数组集合,提供了多个用户信息。然后我们需要将数据以表格的形式,展示到页面上,其中,性别需要转换成中文男女,等级需要…

Consistency Models 阅读笔记

Diffusion models需要多步迭代采样才能生成一张图片,这导致生成速度很慢。Consistency models的提出是为了加速生成过程。 Consistency models可以直接一步采样就生成图片,但是也允许进行多步采样来提高生成的质量。 Consistency models可以从预训练的扩…

制作Go程序的Docker容器

今天突然遇到需要将 Go 程序制作成 Docker 的需求,所以进行了一些研究。方法很简单,但是官方文档和教程有些需要注意的地方,所以写本文进行记录。 源程序 首先介绍一下示例程序,示例程序是一个 HTTP 服务器,会显示si…

Docker安装Zookeeper

天行健,君子以自强不息;地势坤,君子以厚德载物。 每个人都有惰性,但不断学习是好好生活的根本,共勉! 文章均为学习整理笔记,分享记录为主,如有错误请指正,共同学习进步。…