matlab 的 QSMART 回声消除代码挺有意思的,尤其是你要搞 MRI 图像里的磁化率估计时。这套东西是两阶段的反演流程,主要就是干掉那种恼人的伪影和斑纹,在静脉边上效果还不错。它不靠传统的 SHARP 腐蚀那一套,而是直接用空间相关滤波皮质和血管的 mask,挺省事。
QSMART 的安装稍微麻烦点,依赖项不少,不过你直接把它们丢进QSMART_toolbox_v1.0
就行。再加上你得装几个常用的工具包,像 NIfTI 工具那类。只要 MATLAB 装得干净,基本没啥坑。
用法也比较直接,主要是跑Demo_QSMART.m
这个脚本。你只要指定一下 DICOM 的相位和幅度文件夹路径,剩下就靠默认参数先跑起来看看,反正参数针对 Siemens 7T 的数据已经调过了,蛮稳的。
如果你做的是脑影像、QSM、或者想试试非 SHARP 方法的 MRI 反演,这套代码还蛮值得一试的。嗯,记得检查你自己的数据格式匹不匹配,有些 DICOM 结构不标准的话得手动调下。