首页/新闻资讯/正文详情

医学影像预处理

发布时间:2026/11/12 18:28:08 来源:云帆数科 栏目:资讯中心
医学影像预处理
本操作都是在ubuntu20.04上进行的Linux系统一、转换格式核磁影像一般拍出来都是一系列的.dcm格式需要先把其压缩成NIfTI格式.nii.gz我们使用dcm2niix这款工具它可以从不同序列的DICOM文件中准确地重建出3D体积图像。可以使用dcm2niix的GUI版本或命令行工具。dcm2niix安装方法:sudo apt update sudo apt-get install dcm2niix基础命令与参数‌最小化命令‌默认压缩dcm2niix DICOM 文件夹路径输出文件位于当前目录自动压缩为.nii.gz。‌指定输出与压缩‌dcm2niix -z y -o /path/to/output /path/to/dicom_input。‌关键参数说明‌-z y启用 GZIP 压缩生成.nii.gz-z n为不压缩-z i可选 Zstandard 压缩需编译支持。-o 路径指定输出目录。-f 格式自定义文件名如%p_%s代表协议名_序列号。-b y同时生成 BIDS 兼容的.json元数据文件二、提取脑组织掩膜Mask2.1 安装工具FSL使用的工具是FSL根据FSL官网的步骤指令curl -Ls https://fsl.fmrib.ox.ac.uk/fsldownloads/fslconda/releases/getfsl.sh | sh -s这一步要等很久~很久~很久~要有耐心安装好后修改环境变量打开~/.bashrc添加下面几句指令export FSLDIRpath/to/fsl export PATH${FSLDIR}/bin:${PATH} source ${FSLDIR}/etc/fslconf/fsl.sh保存好后运行source ~/.bashrc2.2 安装工具3D Slicer下载3D Slicer Linux版, 以5.10.0版本为例tar xzvf Slicer-5.10.0-linux-amd64.tar.gz #解压缩到当前文件夹 cd Slicer-5.10.0-linux-amd64 #cd到解压缩好的文件夹 ./Slicer # 运行3D Slicer2.3 去颅骨操作针对DWI弥散加权成像的.nii.gz格式影像要进行去除颅骨和软组织的“颅骨剥离”Skull Stripping最专业且主流的工具是FSL的BET (Brain Extraction Tool)。由于DWI图像对比度与结构像T1不同直接处理时需要一些技巧。我们不直接对信噪比较低、有扭曲的DWI原始图像进行剥离而是先计算一张对比度更好的平均b0图像无弥散加权图像在其上生成mask再将mask应用到整个4D DWI数据上。PS理论上用对应的T2 flair生成mask效果最好因为对比度极佳T2 flair 抑制了脑脊液信号同时脑组织灰质/白质与颅骨、软组织之间有非常清晰的边界这让bet算法能非常准确地勾勒出脑轮廓图像质量高T2 flair 是结构像信噪比高、没有 DWI 那样的扭曲和信号衰减生成的mask边缘平滑、精确。但由于T2 flair很难和dwi完美对齐而配准时会有误差配准指令如下# 1. 对T2 flair做颅骨剥离生成蒙版 bet T2_flair.nii.gz T2_brain.nii.gz -m # 第 2 步在 T2 flair 上生成脑蒙版 bet T2_flair.nii.gz T2_brain -m # 第 3 步将 T2 flair 线性配准到 DWI 空间以 DWI 的第一个 volume 为参考 flirt -in T2_flair.nii.gz -ref DWI.nii.gz -out T2_flair_reg.nii.gz -omat T2_to_DWI.mat # 第 4 步将蒙版用同样的变换矩阵映射到 DWI 空间 # 注意这里一定要用 nearestneighbour 插值保证蒙版是二值的0和1 flirt -in T2_flair_mask.nii.gz -ref DWI.nii.gz -out mask_reg.nii.gz -init T2_to_DWI.mat -applyxfm -interp nearestneighbour # 第 5 步将配准后的蒙版应用到原始 DWI 数据 fslmaths DWI.nii.gz -mul mask_reg.nii.gz DWI_brain.nii.gz原始dwi影像如下包含颅骨及头皮组织等经T2 flair获取mask去颅骨结果如下仍有少许头皮组织未去除干净博主猜测是因为配准有误差有大神知道原因麻烦评论区告知感谢感谢博主研究了下采集到的DWI影像中有40张影像其中包含1个b值因此使用b0值计算掩膜来提取干净的脑组织。查看了下DWI的nii.gz图像维度提取b0值和三维DWI值fslroi DWI.nii.gz DWI_0.nii.gz 0 1 #从第一个volume开始提取一个 fslroi DWI.nii.gz DWI_1.nii.gz 1 1 #从第二个volume开始提取一个b0是最后一个volume,对提取的b0做颅骨剥离bet DWI_1.nii.gz b0_brain.nii.gz -m -f 0.5将蒙版应用到DWI文件fslmaths DWI_1.nii.gz -mul b0_brain_mask.nii.gz b0_dwi_brain.nii.gz结果如下干净多了但影像还是很模糊后续再研究下如何提高清晰度

相关推荐

Qwen-RobotWorld:通过语言条件视频生成统一具身世界建模
Qwen-RobotWorld:通过语言条件视频生成统一具身世界建模

26年7月来自阿里千问的论文“Qwen-RobotWorld Technical Report: Unifying Embodied World Modeling through Language-Conditioned Video Generation”。 QWEN-ROBOTWORLD,是一个面向具身智能、由语言驱动的视频世界模型。该模型以自然语言作为统一的动作接口&… · 2026/11/12 18:28:08

大模型压缩常用方法
大模型压缩常用方法

一、为什么需要压缩Llama 3.1 405B FP16 ≈ 810 GB 显存。单张 H100 仅 80 GB。压缩目标:减少参数量 / 降低精度 / 迁移知识,在精度损失可控的前提下降低部署门槛。三种核心方法:方法原理适用阶段剪枝移除贡献小的权重或神经元训练后 / 训练中… · 2026/10/4 6:26:45

万能提示词模板:AI交互与内容创作的高效指南
万能提示词模板:AI交互与内容创作的高效指南

1. 什么是万能提示词模板?在内容创作和AI交互领域,提示词(Prompt)就像一把打开知识宝库的钥匙。我从业十年来发现,90%的创作者在使用AI工具时都卡在了"不知道该怎么问"这个环节。万能提示词模板就是经过实战… · 2026/10/2 23:36:22

随州电工证在哪里考 郑州报考避坑指南:证书过期了不知道怎么复审延期
随州电工证在哪里考 郑州报考避坑指南:证书过期了不知道怎么复审延期

随州电工证在哪里考 郑州报考避坑指南:证书过期了不知道怎么复审延期 证书过期了不知道怎么复审延期?你不是一个人。很多人拿到电工证后,以为一劳永逸,结果几年后发现证过期了,找人问又不清楚流程,甚至被坑了钱还没拿到证。这篇文章就帮你搞清楚【随州… · 2026/11/12 18:27:27

龙海哪里可以办电工证值不值得考?学历不高也能拿证指南
龙海哪里可以办电工证值不值得考?学历不高也能拿证指南

龙海哪里可以办电工证值不值得考?学历不高也能拿证指南 你是不是也遇到过这样的问题? 学历不高怕报不上名 ,想考个电工证,但又担心自己是不是不够资格?其实,别担心,很多像你一样学历不高但经验丰富的老电工,都成功考取了电工证,而且现在政策也更加… · 2026/11/12 18:27:27

上海电工证要取消吗?全国通用的实操考试怎么过
上海电工证要取消吗?全国通用的实操考试怎么过

上海电工证要取消吗?全国通用的实操考试怎么过 实操考试心里没底怕挂科?上海电工证要取消吗?很多人在报名前都在问,尤其是在全国通用的政策下,到底还值不值得花时间去考。作为从业多年的老手,我深知这种焦虑,但更清楚电工证的实际价值,下面我就带大家… · 2026/11/12 18:27:27

考了电工证要年审吗 电子证书怎么转异地
考了电工证要年审吗 电子证书怎么转异地

考了电工证要年审吗 电子证书怎么转异地 你是不是也像很多人一样, 不知道去哪报名怕被中介坑 ?电工证考下来了,但对年审流程一无所知,更别说电子证书怎么转异地了。别急,下面我会从郑州金水区项目现场管理员的实际角度出发,带你理清楚【考了电工证要… · 2026/11/12 18:27:27

电工证紧缺工种有哪些通过率揭秘:工地太忙根本没时间复习考试
电工证紧缺工种有哪些通过率揭秘:工地太忙根本没时间复习考试

电工证紧缺工种有哪些通过率揭秘:工地太忙根本没时间复习考试 工地太忙根本没时间复习考试,这是很多一线电工的真实写照。尤其是南阳卧龙区公路工程从业者,天天穿梭在钢筋水泥之间,连吃饭时间都成了奢侈,更别说抽出时间备考。但你知道吗?… · 2026/11/12 18:27:27

平舆考电工证在哪网上查询证书过期怎么办
平舆考电工证在哪网上查询证书过期怎么办

平舆考电工证在哪网上查询证书过期怎么办 证书过期了不知道怎么复审延期,你是不是也遇到过这种情况?别急,今天就带你从【平舆考电工证在哪】这个关键词出发,一步步搞清楚如何在网上查询、如何复审、怎么避免掉坑,尤其是对于南阳卧龙区的在职建筑工人来说… · 2026/11/12 18:26:27

电工证丢失还能补吗?千万别踩坑,丢了证等于失业
电工证丢失还能补吗?千万别踩坑,丢了证等于失业

电工证丢失还能补吗?千万别踩坑,丢了证等于失业 你正在找工作,但手里的电工证突然丢了,这下怎么办?别急,先别慌, 电工证丢失还能补吗 ?答案是肯定的,但 千万别踩坑 ,补办流程如果不了解,可能会耽误你上岗、甚至影响到你的收入。… · 2026/11/12 14:09:01

电力电缆作业证直播报名入口与办理流程良心建议
电力电缆作业证直播报名入口与办理流程良心建议

电力电缆作业证直播报名入口与办理流程良心建议 你是不是也担心学历不高报不了名?别急,今天就给你讲讲 电力电缆作业证直播 的报名流程和那些你必须知道的干货, 良心建议 在这里,别错过。 电力电缆作业证直播报名需要什么条件?… · 2026/11/11 19:07:25

英德电工证报考费多少钱,怕考不过白交培训费?看完这篇别再踩坑
英德电工证报考费多少钱,怕考不过白交培训费?看完这篇别再踩坑

英德电工证报考费多少钱,怕考不过白交培训费?看完这篇别再踩坑 怕考不过白交培训费?英德电工证报考费多少钱,你真的了解吗?现在市面上培训机构五花八门,价格参差不齐,很多人一不小心就花了冤枉钱,结果还考不过。别急,今天就来给你说清楚英德电工证的… · 2026/11/11 17:14:54

电工证设年审有效期多久官方回应来了,工地太忙根本没时间复习考试
电工证设年审有效期多久官方回应来了,工地太忙根本没时间复习考试

电工证设年审有效期多久官方回应来了,工地太忙根本没时间复习考试 工地太忙根本没时间复习考试,这是很多一线电工、焊工、高处作业人员的真实写照。尤其是像你这样刚毕业的鹤壁工程类毕业生,一边要找工作,一边要考证,时间紧张,复习效率低。但你可能不知… · 2026/11/12 0:00:56

新乡电工证实操怕挂科?电工作业动漫图+报考避坑全解析
新乡电工证实操怕挂科?电工作业动漫图+报考避坑全解析

新乡电工证实操怕挂科?电工作业动漫图+报考避坑全解析 实操考试心里没底怕挂科?你不是一个人。新乡的水利工程从业者,很多在报考电工证的时候都担心自己操作不熟练、考试不通过, 千万别踩坑 ,选对资料和方法真的能提高通过率。… · 2026/11/12 0:00:56

网电工证千万别踩坑:实操考试心里没底怕挂科怎么破?
网电工证千万别踩坑:实操考试心里没底怕挂科怎么破?

网电工证千万别踩坑:实操考试心里没底怕挂科怎么破? 实操考试心里没底怕挂科,这是很多准备报考【网电工证】的学员最真实的感受。尤其是在高压电工、低压电工等工种上,操作环节的容错率极低,一旦失误就可能直接被判定不合格。很多人因为担心挂科而放弃考… · 2026/11/12 0:00:56

双极步进电机控制实战:DRV8818PWPR与PIC18F86K22驱动调试
双极步进电机控制实战:DRV8818PWPR与PIC18F86K22驱动调试

/* MD / 富文本中的 .toc(含博客园搬家等嵌套结构);.toc-box 在侧栏,不受影响 */#content_views .toc,/* 编辑器常在目录前后插入空 p(:empty 仍占 20px),一并去掉避免顶空隙 */#content_views.markdown_views > p:empty:has(+ .toc),#content_views.markdown_views … · 2026/11/12 4:35:34

STM32F207ZG + DRV8818PWPR 步进电机驱动实战解析
STM32F207ZG + DRV8818PWPR 步进电机驱动实战解析

/* MD / 富文本中的 .toc(含博客园搬家等嵌套结构);.toc-box 在侧栏,不受影响 */#content_views .toc,/* 编辑器常在目录前后插入空 p(:empty 仍占 20px),一并去掉避免顶空隙 */#content_views.markdown_views > p:empty:has(+ .toc),#content_views.markdown_views … · 2026/11/12 4:36:01

DRV8818PWPR+dsPIC30F4013工业级双极步进电机控制方案
DRV8818PWPR+dsPIC30F4013工业级双极步进电机控制方案

/* MD / 富文本中的 .toc(含博客园搬家等嵌套结构);.toc-box 在侧栏,不受影响 */#content_views .toc,/* 编辑器常在目录前后插入空 p(:empty 仍占 20px),一并去掉避免顶空隙 */#content_views.markdown_views > p:empty:has(+ .toc),#content_views.markdown_views … · 2026/11/11 17:15:20

了解更多?预约专属演示

我们的顾问将为您一对一讲解产品与方案

企业微信二维码