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

  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.hqwc.cn/news/192323.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

CUDA编程一、基本概念和cuda向量加法

目录 一、cuda编程的基本概念入门 1、GPU架构和存储结构 2、cuda编程模型 3、cuda编程流程 二、cuda向量加法实践 1、代码实现 2、代码运行和结果 有一段时间对模型加速比较感兴趣,其中的一块儿内容就是使用C和cuda算子优化之类一起给模型推理提速。之前一直…

HAL库STM32串口开启DMA接收数据

STM32CubeMx的配置 此博客仅仅作为记录,这个像是有bug一样,有时候好使,有时候不好,所以趁现在好使赶紧记录一下,很多地方用到串口接收数据,DMA又是一种非常好的接收方式,可以节约CPU的时间&…

Unity - Cinemachine

动态获取Cinemachine的内部组件 vCam.GetCinemachineComponent<T>() 动态修改Cinemachine的Transposer属性 var vCamComp transfrom.GetComponent<CinemachineVirtualCamera>(); var transposerComp vCamComp.GetCinemachineComponent<CinemachineTransposer&…

2023最新最全【OpenMV】 入门教程

1. 什么是OpenMV OpenMV 是一个开源&#xff0c;低成本&#xff0c;功能强大的 机器视觉模块。 OpenMV上的机器视觉算法包括 寻找色块、人脸检测、眼球跟踪、边缘检测、标志跟踪 等。 以STM32F427CPU为核心&#xff0c;集成了OV7725摄像头芯片&#xff0c;在小巧的硬件模块上&a…

我们应该如何理解Java集合框架的关键知识点?

我们应该如何理解Java集合框架的关键知识点&#xff1f; java集合是教存取数据的一个容器&#xff0c;涵盖了各种存和取的方式&#xff0c;应用在不同的工作场景中&#xff0c;要想了解java集合的相关知识&#xff0c;建议先好好学习一下数据结构这本书。最近很多小伙伴找我&am…

使用Lychee搭建个人图片存储系统并进行远程访问设置实现公网访问本地私人图床

文章目录 1.前言2. Lychee网站搭建2.1. Lychee下载和安装2.2 Lychee网页测试2.3 cpolar的安装和注册 3.本地网页发布3.1 Cpolar云端设置3.2 Cpolar本地设置 4.公网访问测试5.结语 1.前言 图床作为图片集中存放的服务网站&#xff0c;可以看做是云存储的一部分&#xff0c;既可…

qt library创建和使用

1、创建library 2、修改library中的代码 3、把library进行编译&#xff0c;编译后会生成相关文件 4、把编译后的文件拷贝到主程序目录下面。 5、并把library中的testlib头文件拷贝到主程序&#xff0c;并在pro文件加入&#xff08;这里在后面使用library中函数有关&#xff0…

【C++】类和对象(7)--友元, static成员

目录 一 友元 1 友元概念 2 友元函数 3 友元类 二 static成员 1 概念 2 用法 3 static成员特性 4 例题 一 友元 1 友元概念 友元提供了一种突破封装的方式&#xff0c;有时提供了便利。但是友元会增加耦合度&#xff0c;破坏了封装&#xff0c;所以 友元不宜多用。 …

Vue3 函数式弹窗

运行环境 vue3vitetselement-plus 开发与测试 1. 使用h、render函数创建Dialog 建议可在plugins目录下创建dialog文件夹&#xff0c;创建index.ts文件&#xff0c;代码如下 import { h, render } from "vue";/*** 函数式弹窗* param component 组件* param opti…

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

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

紫色调城市和奔跑人物剪影背景工会工作总结汇报PPT模板

这是一套紫色调城市和奔跑人物剪影背景工会工作总结汇报PPT模板&#xff0c;共33页&#xff1b; PPT模板封面&#xff0c;使用了蓝天白云、城市剪影、奔跑人物剪影背景图片。中间填写工会工作总结汇报PPT标题。界面色彩丰富充满活力。 PowerPoint模板内容页&#xff0c;由31张…

2020年06月 Scratch(二级)真题解析#中国电子学会#全国青少年软件编程等级考试

Scratch等级考试(1~4级)全部真题・点这里 一、单选题(共25题,每题2分,共50分) 第1题 如下图所示脚本运行的结果是()? A:画一条直线 B:画一个三角形 C:画一个圆形 D:画一条虚线 答案:D 第2题 运行如下图所示脚本,下面选项中说法错误的是? A:“笔的颜色”…