当前位置: 首页 > news >正文

IEEE-TMI:张孝勇团队开发小鼠精细脑结构自动分割的深度学习算法

近日,复旦大学类脑智能科学与技术研究院青年研究员张孝勇课题组联合德国亥姆霍兹慕尼黑研究中心,在医学图像处理领域顶尖期刊《IEEE医学影像汇刊》(IEEE Transactions on Medical Imaging,TMI) 发表了题为《MouseGAN++:用于小鼠大脑MRI多模态合成和结构分割的无监督解耦及对比学习表征算法》(“MouseGAN++: Unsupervised Disentanglement and Contrastive Representation for Multiple MRI Modalities Synthesis and Structural Segmentation of Mouse Brain”)的研究成果,提出无监督解耦和对比学习表征算法解决小鼠大脑精细脑区分割的难题。

作为最重要的模式生物之一,小鼠在神经科学、肿瘤学和药学等领域发挥着不可替代的作用,是实现基础研究到临床应用转化的重要桥梁。由于小鼠和人类大脑在进化上具有保守性,小鼠大脑已被证明是理解人类大脑的重要基石之一。基于磁共振成像(MRI)数据,对小鼠大脑精细结构进行准确分割是理解小鼠大脑脑区功能的重要分析步骤。然而,迄今为止这仍是一项具有挑战性的任务。首先,脑结构的精确分割通常需要高分辨率的多模态MRI数据,但采集多模态数据的时间成本往往过高。再者,从模型设计上如何高效融合多模态MRI信息也面临着方法学上的挑战。其次,目前尚缺乏用于小鼠脑结构自动分割和标注的工具箱。

为解决这一问题,张孝勇课题组提出了一种新型深度学习框架,MouseGAN++(图1),完成“合成+分割”的任务。模型利用属性空间和个体空间的对比损失来度量样本间的相似度(图2),同时使对抗学习和对比学习在训练过程中同步交替进行,以更好地解耦出脑结构语义特征,使MRI在模态生成前后能保持结构一致性,进而促成高效的多模态图像融合来实现小鼠大脑结构的精细分割。

具体来说,针对多模态MRI自身的性质,我们将所有模态的图像分别解耦到属性空间和内容空间。现有的解耦方法通常使用高斯先验来约束属性隐空间的分布,但若用于多模态数据,可能会面临来自不同模态的属性特征之间解耦不充分的情况。为了克服这一问题,本文提出一种新的归纳偏置,即利用对比学习先验来指导网络学习属性空间的分布。

同时,考虑到脑结构的解剖特点(对称的左右脑半球的解剖结构应该具有相似的特征),若使用patch-level的训练方式,就会使网络在内容空间里区分它们,从而导致假阴性结果(错误的负样本对)。因此,本文采用轴位(slice-level)的训练方式以纳入全局语义信息。此外,对比学习通常需要设计额外的pretext-task进行预训练,pretext-task的设计质量以及它与下游任务之间的差距也会影响到模型的最终效果。

图1:MouseGAN++模型图。给定任意模态作为网络的输入,(a)模态生成模块:训练基于特征解耦和对比学习表征来合成所有模态。(b)重用内容编码器并将模态生成模块作为辅助网络插补缺失模态。

针对以上问题,MouseGAN++包含一个模态生成模块,我们在该模块引入两个新颖的对比损失函数,可将多模态图像特征投射到共享内容隐空间(用来编码模态无关的脑结构特征),以及模态相关的属性空间。随后,将内容隐空间中的向量与其他模态的属性特征向量相结合,以插补生成其他模态的图像。内容对比损失可以迫使网络在图像生成过程中避免混淆结构信息。同时,在对抗性训练期间重用属性和内容编码器,以同步优化对抗损失和对比学习损失。共享的内容空间还有助于分割模块中的解码器训练。此外,使用该模型插补模态也可以扩增训练数据集,使网络能够在多模态数据里学习模态无关的结构语义特征,从而增强多模态融合的效果。

图2:通过对比学习进行属性判别和内容判别。(左)对于属性空间,来自同一模态的样本被定义为正样本,在隐空间中相互拉近。不同模态的样本为负样本,相互推远。(右)来自同一个体的样本即使被变换了模态,都被视为正样本,它们的结构信息在内容隐空间内应保持高度一致。而相同模态下的不同个体之间则被视为负样本。

表1:MouseGAN++模型与State-of-the-art方法的性能对比。

如表1所示,与当前最先进的9种相关方法相比,以T1w和T2w为测试模态,平均DICE系数分别达到87.9%和90.0%,性能提高了约 +10%。上述图表及文中实验结果表明,MouseGAN++在生成和分割任务上都实现了显著的性能提升。本深度学习框架的源代码已开源:https://github.com/yu02019 。

在未来的工作中,一个有前景的研究方向是将我们开发的大脑提取工具(BEN)与MouseGAN++集成,形成一个鲁棒的、高通量的、端到端的神经影像处理流程。

复旦大学类脑智能科学与技术研究院张孝勇青年研究员为文章的通讯作者,德国亥姆霍兹慕尼黑研究中心彭廷莹研究员为共同通讯作者,复旦大学类脑智能科学与技术研究院冯建峰教授对本文做出重要指导。复旦大学类脑智能科学与技术研究院博士研究生余子奇为第一作者。

该研究得到了国家自然科学基金委面上项目、重大项目及上海市“脑与类脑智能基础转化应用研究”市级重大专项等的经费支持。

原文链接:

https://ieeexplore.ieee.org/document/9966838 

参考文献:

Yu, Ziqi, et al. "BEN: a generalizable Brain Extraction Net for multimodal MRI data from rodents, nonhuman primates, and humans." bioRxiv (2022).(elife 接收)

张孝勇,复旦大学类脑智能科学与技术研究院青年研究员,中国图象图形学会脑图谱专委会委员,中华放射学会国际交流工作组委员,医学图像计算青年研讨会委员。

研究方向主要聚焦在大脑的磁共振成像、人工智能分析及分子影像研究。在磁共振代谢成像技术领域做出了若干原创性贡献。目前(2022.12)已发表学术论文50余篇,其中以通讯作者/第一作者发表学术论文30余篇,代表性成果发表在Advanced Science,IEEE Transactions on Medical Imaging,Cancer Research,Cell Reports等权威学术期刊。

研究成果曾被国际医学磁共振学会(ISMRM)以研究亮点报道。作为课题负责人主持多项国家自然科学基金项目及上海市项目,并以研究骨干参与国家自然科学基金重大项目、集成项目等。

来源:复旦类脑智能科学与技术研究院

相关文章:

IEEE-TMI:张孝勇团队开发小鼠精细脑结构自动分割的深度学习算法

近日,复旦大学类脑智能科学与技术研究院青年研究员张孝勇课题组联合德国亥姆霍兹慕尼黑研究中心,在医学图像处理领域顶尖期刊《IEEE医学影像汇刊》(IEEE Transactions on Medical Imaging,TMI) 发表了题为《MouseGAN:用于小鼠大脑…...

八股文之面向对象和面向过程的区别

面向对象(Object-Oriented)和面向过程(Procedural)是两种不同的编程思想。 面向过程是以任务为中心,将程序分解成一系列步骤,在每个步骤中定义一个函数来完成特定的任务。它主要关注程序执行的过程和如何组…...

SpringBoot使用Redis实现分布式缓存

✅作者简介:2022年博客新星 第八。热爱国学的Java后端开发者,修心和技术同步精进。 🍎个人主页:Java Fans的博客 🍊个人信条:不迁怒,不贰过。小知识,大智慧。 💞当前专栏…...

Three——二、加强对三维空间的认识

Three——二、加强对三维空间的认识 接上个例子我们接着往下看 辅助观察坐标系 THREE.AxesHelper()的参数表示坐标系坐标轴线段尺寸大小,你可以根据需要改变尺寸。 使用方法: // AxesHelper:辅助观察的坐标系 const axesHelper new THRE…...

【Java】Java8接口中方法区别和使用

Java接口说明 jdk1.8之前接口只能是抽象方法。实现接口必须重写所有方法,比较麻烦。在java8中,支持default和static方法,这样,实现接口时,可以选择是否对default修饰的方法重写。 抽象方法 接口当中的抽象方法&#x…...

WPF 控件库Live Charts 折线图多折线比较问题处理

使用Live Charts功能对比多条折线时当Label不是一一对应时会发现折线无法对比如 Labels List<double> list2 new List<double>(); list2.Add(2.1); //x为0.5时 list2.Add(2.2); //x为0.6时 …...

接口优化方案

前言 最近随着国产化热潮&#xff0c;公司的用于营业的电脑全部从windows更换成了某国产化电脑&#xff0c;换成国产化之后&#xff0c;我们系统的前台web界面也由之前的jsp页面重构成vue.所以之前的一体式架构也变成了前后端分离的架构。但是在更换过程后&#xff0c;发现一些…...

《商用密码应用与安全性评估》第二章政策法规2.1网络空间安全形式与商业密码工作

一、国际国内网络空间安全形势 网络空间已成为与陆地、海洋、天空、太空同等重要的人类第五空间。 1.国际形势 网络空间安全纳入国家战略 网络攻击在国家对抗中深度应用 网络空间已逐步深入网络底层固件 2.国内形势 核心技术仍受制于人 信息产品存在巨大安全隐患 关…...

C#实现将文件、文件夹压缩为压缩包

C#实现将文件、文件夹压缩为压缩包 一、C#实现将文件、文件夹压缩为压缩包核心 1、介绍 Title&#xff1a;“基础工具” 项目&#xff08;压缩包帮助类&#xff09; Description步骤描述&#xff1a; 1、创建 zip 存档&#xff0c;该文档包含指定目录的文件和子目录&#xf…...

程序员跳槽,要求涨薪50%过分吗?

如果问在TI行业涨工资最快的方式是什么&#xff1f; 回答最多的一定是&#xff1a;跳槽&#xff01; 前段时间&#xff0c;知乎上这样一条帖子引发了不少IT圈子的朋友的讨论 &#xff0c;有网友提问 “程序员跳槽要求涨薪50%过分吗&#xff1f;” 截图来源于知乎&#xff0c;…...

Java核心技术 卷1-总结-10

Java核心技术 卷1-总结-10 通配符类型通配符概念通配符的超类型限定无限定通配符通配符捕获 通配符类型 通配符概念 通配符类型中&#xff0c;允许类型参数变化。 例如&#xff0c;通配符类型Pair<? extends Employee>表示任何泛型Pair类型&#xff0c;它的类型参数是…...

React Props

state 和 props 主要的区别在于 props 是不可变的&#xff0c;而 state 可以根据与用户交互来改变。 所以&#xff0c;有些容器组件需要定义 state 来更新和修改数据。 而子组件只能通过 props 来传递数据。 props 使用 Demo.js &#xff1a; import React from reactfunct…...

【Hello Network】协议

作者&#xff1a;小萌新 专栏&#xff1a;网络 作者简介&#xff1a;大二学生 希望能和大家一起进步 本篇博客简介&#xff1a;简单介绍下协议并且设计一个简单的网络服务器 协议 协议的概念结构化数据传输序列化和反序列化网络版计算机服务端代码协议定制客户端代码服务线程执…...

零项目零科研,本科排名倒数,一战上岸上海交大电子与通信工程

笔者来自通信考研小马哥23上交819全程班学员 本科就读于哈工大&#xff08;威海&#xff09;&#xff0c;本科成绩很差&#xff0c;专业排名62/99&#xff0c;没有科研&#xff0c;没有实验室&#xff0c;没有项目&#xff0c;连最基本大家都会参加的科技立项我四年也没有参与…...

NOIP模拟赛 T3区间

题目大意 有 n n n个数字&#xff0c;第 i i i个数字为 a i a_i ai​。有 m m m次询问&#xff0c;每次给出 k i k_i ki​个区间&#xff0c;每个区间表示第 l i , j l_{i,j} li,j​到第 r i , j r_{i,j} ri,j​个数字&#xff0c;求这些区间中一共出现了多少种不同的数字。部…...

【Python】如何用pyth做游戏脚本(太简单了吧)

文章目录 前言一、开发前景二、开发流程3.1、获取窗口句柄&#xff0c;把窗口置顶3. 2、截取游戏界面&#xff0c;分割图标&#xff0c;图片比较 二、程序核心-图标连接算法&#xff08;路径寻找&#xff09;四、开发总结五、源码总结 前言 简述&#xff1a;本文将以4399小游戏…...

【Linux】磁盘与文件系统

目录 一、磁盘的物理结构 二、磁盘逻辑抽象 三、文件系统 1、Super Block 2、Group Descriptor Table 3、inode Table 4、Data Blocks 5、inode Bitmap 6、Block Bitmap 四、Linux下文件系统 1、inode与文件名 2、文件的增删查改 2.1、查看文件内容 2.2、删除文件…...

Transformer中的注意力机制及代码

文章目录 1、简介2、原理2.1 什么是注意力机制2.2 注意力机制在NLP中解决了什么问题2.3 注意力机制公式解读2.4 注意力机制计算过程 3、单头注意力机制与多头注意力机制4、代码4.1 代码14.2 代码2 1、简介 最近在学习transformer&#xff0c;首先学习了多头注意力机制&#xf…...

ChatGPT在连续追问下对多线程和双重检查锁模式的理解--已经超越中级程序员

一、问&#xff1a; private static final Map<Method, GZHttpClientResultModel> CACHE_RESULT_MODEL new ConcurrentHashMap<>();public void abc(Method method){cacheResultMode(method);GZHttpClientResultModel model CACHE_RESULT_MODEL.get(method);}pr…...

每天一道大厂SQL题【Day22】华泰证券真题实战(四)

每天一道大厂SQL题【Day22】华泰证券真题实战(四) 大家好&#xff0c;我是Maynor。相信大家和我一样&#xff0c;都有一个大厂梦&#xff0c;作为一名资深大数据选手&#xff0c;深知SQL重要性&#xff0c;接下来我准备用100天时间&#xff0c;基于大数据岗面试中的经典SQL题&…...

Qt Http Server模块功能及架构

Qt Http Server 是 Qt 6.0 中引入的一个新模块&#xff0c;它提供了一个轻量级的 HTTP 服务器实现&#xff0c;主要用于构建基于 HTTP 的应用程序和服务。 功能介绍&#xff1a; 主要功能 HTTP服务器功能&#xff1a; 支持 HTTP/1.1 协议 简单的请求/响应处理模型 支持 GET…...

【开发技术】.Net使用FFmpeg视频特定帧上绘制内容

目录 一、目的 二、解决方案 2.1 什么是FFmpeg 2.2 FFmpeg主要功能 2.3 使用Xabe.FFmpeg调用FFmpeg功能 2.4 使用 FFmpeg 的 drawbox 滤镜来绘制 ROI 三、总结 一、目的 当前市场上有很多目标检测智能识别的相关算法&#xff0c;当前调用一个医疗行业的AI识别算法后返回…...

在web-view 加载的本地及远程HTML中调用uniapp的API及网页和vue页面是如何通讯的?

uni-app 中 Web-view 与 Vue 页面的通讯机制详解 一、Web-view 简介 Web-view 是 uni-app 提供的一个重要组件&#xff0c;用于在原生应用中加载 HTML 页面&#xff1a; 支持加载本地 HTML 文件支持加载远程 HTML 页面实现 Web 与原生的双向通讯可用于嵌入第三方网页或 H5 应…...

html css js网页制作成品——HTML+CSS榴莲商城网页设计(4页)附源码

目录 一、&#x1f468;‍&#x1f393;网站题目 二、✍️网站描述 三、&#x1f4da;网站介绍 四、&#x1f310;网站效果 五、&#x1fa93; 代码实现 &#x1f9f1;HTML 六、&#x1f947; 如何让学习不再盲目 七、&#x1f381;更多干货 一、&#x1f468;‍&#x1f…...

嵌入式学习笔记DAY33(网络编程——TCP)

一、网络架构 C/S &#xff08;client/server 客户端/服务器&#xff09;&#xff1a;由客户端和服务器端两个部分组成。客户端通常是用户使用的应用程序&#xff0c;负责提供用户界面和交互逻辑 &#xff0c;接收用户输入&#xff0c;向服务器发送请求&#xff0c;并展示服务…...

从面试角度回答Android中ContentProvider启动原理

Android中ContentProvider原理的面试角度解析&#xff0c;分为​​已启动​​和​​未启动​​两种场景&#xff1a; 一、ContentProvider已启动的情况 1. ​​核心流程​​ ​​触发条件​​&#xff1a;当其他组件&#xff08;如Activity、Service&#xff09;通过ContentR…...

在树莓派上添加音频输入设备的几种方法

在树莓派上添加音频输入设备可以通过以下步骤完成&#xff0c;具体方法取决于设备类型&#xff08;如USB麦克风、3.5mm接口麦克风或HDMI音频输入&#xff09;。以下是详细指南&#xff1a; 1. 连接音频输入设备 USB麦克风/声卡&#xff1a;直接插入树莓派的USB接口。3.5mm麦克…...

Python 训练营打卡 Day 47

注意力热力图可视化 在day 46代码的基础上&#xff0c;对比不同卷积层热力图可视化的结果 import torch import torch.nn as nn import torch.optim as optim from torchvision import datasets, transforms from torch.utils.data import DataLoader import matplotlib.pypl…...

02.运算符

目录 什么是运算符 算术运算符 1.基本四则运算符 2.增量运算符 3.自增/自减运算符 关系运算符 逻辑运算符 &&&#xff1a;逻辑与 ||&#xff1a;逻辑或 &#xff01;&#xff1a;逻辑非 短路求值 位运算符 按位与&&#xff1a; 按位或 | 按位取反~ …...

Tauri2学习笔记

教程地址&#xff1a;https://www.bilibili.com/video/BV1Ca411N7mF?spm_id_from333.788.player.switch&vd_source707ec8983cc32e6e065d5496a7f79ee6 官方指引&#xff1a;https://tauri.app/zh-cn/start/ 目前Tauri2的教程视频不多&#xff0c;我按照Tauri1的教程来学习&…...