CV | Medical-SAM-Adapter论文详解及项目实现

夏天｜여름이다

已于 2024-02-08 14:04:43 修改

阅读量1.3k

点赞数 23

分类专栏： - CV - 文章标签：人工智能计算机视觉医学影像分割算法图像分割

于 2024-02-07 16:00:12 首次发布

本文链接：https://blog.csdn.net/weixin_44649780/article/details/136051574

版权

- CV - 专栏收录该内容

22 篇文章 3 订阅

订阅专栏

******************************* 👩‍⚕️ 医学影像相关直达👨‍⚕️*******************************

CV | SAM在医学影像上的模型调研【20240207更新版】-CSDN博客

CV | Segment Anything论文详解及代码实现

本文主要讲解Medical-SAM-Adapter论文及项目实现~

2023.12.29第七版_Medical SAM Adapter:Adapting Segment Anything Model for Medical Image Segmentation

论文地址：2304.12620.pdf (arxiv.org)

论文代码：KidsWithTokens/Medical-SAM-Adapter: Adapting Segment Anything Model for Medical Image Segmentation (github.com)

1.论文详解

摘要

Segment Anything Model （SAM）在图像分割领域广受欢迎，因为它在各种分割任务中具有令人印象深刻的功能和基于提示的界面。然而，最近的研究和个别实验表明，由于缺乏医学特定知识，SAM在医学图像分割方面表现不佳。这就提出了一个问题，即如何增强SAM对医学图像的分割能力。在本文中，没有对SAM模型进行微调，而是提出了医疗SAM适配器（Med-SA），它使用一种轻巧而有效的适应技术将特定领域的医学知识整合到分割模型中。在Med-SA中，提出了空间深度转置（SD-Trans）来实现2D SAM与3D医学图像的适应，以及超提示适配器（HyP-Adpt）来实现快速条件适应。对不同图像模态的 17 个医学图像分割任务进行了综合评估实验。Med-SA 的性能优于几种最先进的（SOTA）医学图像分割方法，而仅更新了 2% 的参数。

方法

主要讲解方法部分

论文基于SAM,依旧使用图像编码器、提示编码器和mask解码器。

图像编码器：基于标准ViT被MAE训练，这里使用ViT_h/16变种，有14×14窗口的注意力和4个等间距的全局注意力。图像编码器的输出是对输入图像的16倍下采样嵌入。

提示编码器：可以是稀疏的点、框、文本或密集的mask；本文中只关注稀疏提示编码器，它将点和框表示为位置编码，并对每个提示类型进行学习嵌入。

mask解码器：是一个Transformer解码器，包括动态掩模预测头。

SAM 使用双向交叉注意，一个用于提示到图像的嵌入，另一个用于图像到提示的嵌入，在每个块中学习提示和图像嵌入之间的相互作用。在运行了两个块之后，SAM 上采样嵌入图像，MLP 将输出标记映射到一个动态线性分类器，从而预测给定图像的目标掩码。

Med-SA架构(Med-SA architecture)

为了对 SAM 架构进行微调以适应医学图像分割，本文没有完全调整所有参数，而是冻结了预先训练好的 SAM 参数，并在架构的特定位置插入了 Adapter 模块。

Adapter 是一个 bottleneck 结构，它依次包括：下投影、ReLU 激活和上投影。下投影使用简单的 MLP 层将给定的嵌入压缩到较小的维度；上投影使用另一个 MLP 层将压缩的嵌入扩展回其原始维度.

(在第一版里叫MSA,第七版里叫Med-SA)

2D Medical Image Adaption
在 SAM 编码器中，本文为每个 ViT 块部署了两个 Adapter。

修改标准 ViT block (a)，得到 2D Medical Image Adaption (b)

将第一个 Adapter 放在多头注意力之后、残差连接之前
将第二个 Adapter 放在多头注意力之后 MLP 层的残差路径上

紧接着第二个 Adapter 之后，按照一定的比例系数对嵌入进行了缩放
引入缩放因子 s 是为了平衡与任务无关的特征和与任务有关的特征
Decoder Adaption
在 SAM 解码器中，本文为每个 ViT 块部署了Adapter，把这个叫做Hyper-Prompting Adapter

第一个 Adapter 部署在 prompt-to-image 嵌入的多头交叉注意之后，并添加了提示嵌入的残差
本文使用了另一种向下投影来压缩提示嵌入，并在 ReLU 激活之前将其添加到 Adapter 的嵌入上
有助于 Adapter 根据提示信息调整参数，使其更加灵活和通用于不同的模式和下游任务
第二个 Adapter 的部署方式与编码器完全相同，用于调整 MLP 增强嵌入
第三个 Adapter 部署在图像嵌入的残差连接之后，以提示交叉注意
另一个残差连接和层归一化在自适应后连接，以输出最终结果
SD-Trans
尽管 SAM 可以应用于病灶的每个切片以获得最终的分割，但是它没有考虑深度维中的相关性

本文提出了一种新的适配方法，其灵感来源于 image-to-video adaptation，具体架构如 (c)

在每个 block 中，本文将注意力操作分成两个分支：空间分支和深度分支

2.项目实现

2.0.环境设置

Python3.8+docker容器（Ubuntu）

git clone https://github.com/KidsWithTokens/Medical-SAM-Adapter

cd Medical-SAM-Adapter

wget https://dl.fbaipublicfiles.com/segment_anything/sam_vit_b_01ec64.pth

# 自定义req.txt,如本文A1
pip install -r req.txt

因为我没有使用conda新建虚拟环境，自己新建了一个requestion.txt,如【A1】,按照作者要求是torch版本是1.12，我是1.14，并删掉了安装torch的命令。

2.1.数据集设置

官网：ISIC Challenge (isic-archive.com)

新建文件夹data/isic

#下载数据集
wget https://isic-challenge-data.s3.amazonaws.com/2016/ISBI2016_ISIC_Part1_Training_Data.zip

#wget https://isic-challenge-data.s3.amazonaws.com/2016/ISBI2016_ISIC_Part1_Training_GroundTruth.zip

wget https://isic-challenge-data.s3.amazonaws.com/2016/ISBI2016_ISIC_Part1_Test_Data.zip

#wget https://isic-challenge-data.s3.amazonaws.com/2016/ISBI2016_ISIC_Part1_Test_GroundTruth.zip


#下载csv文件
https://github.com/KidsWithTokens/MedSegDiff/blob/master/data/isic_csv/ISBI2016_ISIC_Part3B_Test_GroundTruth.csv

https://github.com/KidsWithTokens/MedSegDiff/blob/master/data/isic_csv/ISBI2016_ISIC_Part3B_Training_GroundTruth.csv

进入文件夹解压缩： unzip '*.zip'

格式如图

2.2.训练

训练数据集1：ISIC2016

python train.py -net sam -mod sam_adpt -exp_name msa_isic -sam_ckpt ./checkpoint/sam/sam_vit_b_01ec64.pth -image_size 1024 -b 32 -dataset isic -data_path ./data/isic

一张单卡24GGPU的情况，batch size为2,17929MiB,如果现存较小，改小batch size或者image size.

100个epoch时，IOU:0.85,DICE:0.91

采样结果

训练数据集2

python train.py -net sam -mod sam_adpt -exp_name msa_kv
asir -sam_ckpt ./sam_vit_b_01ec64.pth -image_size 1024 -b 32 -dataset Kvasir-SEG -data_path /workspace/SAM
/datasets/Kvasir-SEG

3.代码详解

3.1.论文的核心块adapter_block.py

位置：Medical-SAM-Adapter/models/ImageEncoder/vit/adapter_block.py

基于原本的transformer块添加窗口注意力和残差传播块，流程如下：

对输入张量 x 进行归一化操作。
进行注意力计算，包括多头注意力机制和相对位置嵌入。
2D图像的情况使用 Space Adapter 处理注意力计算的输出。
如果数据是 3D ，则对输入进行额外的处理，并将其与 Space Adapter 处理后的结果相加。
对处理后的张量再次进行归一化操作。
将张量传递给 MLP 模块，并加上 MLP Adapter 处理后的结果的 scaled version。
将上述两项相加，并返回最终的输出张量。

3.2.

PS

【PS1】ValueError: num_samples should be a positive integer value, but got num_samples=0

csv文件问题

重新下载csv文件就可以啦，不能使用wget方式下载

【PS2】 TypeError: unsupported operand type(s) for %: 'int' and 'NoneType'

或者出现 ZeroDivisionError: integer division or modulo by zero错误时：

训练命令需要填写或修改vis自定义参数，默认是None,不能填写0，修改为1

扩展

【A1】pip 安装第三方库

aiosignal==1.2.0
alembic==1.10.4
appdirs==1.4.4
astor==0.8.1
asttokens==2.2.1
backcall==0.2.0
beautifulsoup4==4.12.2
blinker==1.6.2
cachetools==4.2.2
certifi==2022.12.7
charset-normalizer==2.0.4
click==8.1.3
cmaes==0.9.1
colorama==0.4.6
colorlog==6.7.0
contextlib2==21.6.0
coverage==6.5.0
coveralls==3.3.1
cucim==23.4.1
cycler==0.11.0
databricks-cli==0.17.7
docker==6.1.1
docopt==0.6.2
einops==0.6.1
entrypoints==0.4
exceptiongroup==1.1.1
executing==1.2.0
filelock==3.12.0
fire==0.5.0
flask==2.3.2
fonttools==4.25.0
future==0.18.3
gdown==4.7.1
gitdb==4.0.10
gitpython==3.1.31
google-auth==2.6.0
google-auth-oauthlib==0.4.4
greenlet==2.0.2
gunicorn==20.1.0
h5py==3.8.0
huggingface-hub==0.14.1
iniconfig==2.0.0
ipython
itk==5.3.0
itk-core==5.3.0
itk-filtering==5.3.0
itk-io==5.3.0
itk-numerics==5.3.0
itk-registration==5.3.0
itk-segmentation==5.3.0
itsdangerous==2.1.2
jedi==0.18.2
jinja2==3.1.2
json-tricks==3.16.1
jsonschema==4.17.3
kornia==0.4.1
lmdb==1.4.1
lucent==0.1.0
mako==1.2.4
mlflow==2.3.1
nibabel==5.1.0
ninja==1.11.1
nni==2.10
nptyping==2.5.0
opencv-python==4.7.0.72
openslide-python==1.1.2
optuna==3.1.1
partd==1.2.0
pluggy==1.0.0
pooch==1.4.0
prettytable==3.7.0
prompt-toolkit==3.0.38
psutil==5.9.5
pyarrow==11.0.0
pyasn1==0.4.8
pyasn1-modules==0.2.8
pydicom==2.3.1
pygments==2.15.1
pynrrd==1.0.0
pyqt5-sip==12.11.0
pyrsistent==0.19.3
pytest==7.3.1
pytest-mock==3.10.0
pythonwebhdfs==0.2.3
pytorch-ignite==0.4.10
querystring-parser==1.2.4
regex==2023.5.5
requests-oauthlib==1.3.0
responses==0.23.1
rsa==4.7.2
safetensors==0.4.1
schema==0.7.5
simplejson==3.19.1
smmap==5.0.0
soupsieve==2.4.1
scikit-image
sqlalchemy==2.0.12
sqlparse==0.4.4
tabulate==0.9.0
tensorboardx==2.2
termcolor==2.3.0
threadpoolctl==2.2.0
tifffile==2021.7.2
timm==0.9.12
tokenizers==0.12.1
tomli==2.0.1

torch-lucent==0.1.8

traitlets==5.9.0
transformers==4.21.3
typeguard==3.0.2
types-pyyaml==6.0.12.9
wcwidth==0.2.6
websocket-client==1.5.1
websockets==11.0.3
werkzeug==2.3.4