2.845

2023影响因子

(CJCR)

  • 中文核心
  • EI
  • 中国科技核心
  • Scopus
  • CSCD
  • 英国科学文摘

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

结合MRF能量和模糊速度的乳腺癌图像分割方法

冯宝 陈业航 刘壮盛 李智 宋嵘 龙晚生

冯宝, 陈业航, 刘壮盛, 李智, 宋嵘, 龙晚生. 结合MRF能量和模糊速度的乳腺癌图像分割方法. 自动化学报, 2020, 46(6): 1188−1199 doi: 10.16383/j.aas.c180759
引用本文: 冯宝, 陈业航, 刘壮盛, 李智, 宋嵘, 龙晚生. 结合MRF能量和模糊速度的乳腺癌图像分割方法. 自动化学报, 2020, 46(6): 1188−1199 doi: 10.16383/j.aas.c180759
Feng Bao, Chen Ye-Hang, Liu Zhuang-Sheng, Li Zhi, Song Rong, Long Wan-Sheng. Segmentation of breast cancer on DCE-MRI images with MRF energy and fuzzy speed function. Acta Automatica Sinica, 2020, 46(6): 1188−1199 doi: 10.16383/j.aas.c180759
Citation: Feng Bao, Chen Ye-Hang, Liu Zhuang-Sheng, Li Zhi, Song Rong, Long Wan-Sheng. Segmentation of breast cancer on DCE-MRI images with MRF energy and fuzzy speed function. Acta Automatica Sinica, 2020, 46(6): 1188−1199 doi: 10.16383/j.aas.c180759

结合MRF能量和模糊速度的乳腺癌图像分割方法

doi: 10.16383/j.aas.c180759
基金项目: 广西高等学校千名中青年骨干教师培育计划项目基金(2018GXQGFB160)资助
详细信息
    作者简介:

    冯宝:2014年获华南理工大学博士学位. 现为中山大学博士后. 主要研究方向为机器学习, 模式识别及其在生物医学信号处理中的应用. E-mail: fengbao1986.love@163.com

    陈业航:桂林电子科技大学硕士研究生. 主要研究方向为生物医学信号处理与模式识别. E-mail: cyh93yl@163.com

    刘壮盛:广东省江门市中心医院副主任医师. 2010年获中山大学影像医学与核医学硕士学位. 主要研究方向为乳腺及骨骼肌肉系统影像诊断. E-mail: zhuangshengliu@126.com

    李智:2003年获电子科技大学博士学位. 桂林航天工业学院教授. 主要研究方向为智能仪器系统, 现代测试理论与技术. E-mail: cclizhi@guet.edu.cn

    宋嵘:2006年获中国香港理工大学生物医学工程博士学位. 现为中山大学生物医学工程学院教授. 主要研究方向为肌肉骨骼建模, 生物医学信号处理, 人体运动分析和机器人辅助中风康复. E-mail: songrong@mail.sysu.edu.cn

    龙晚生:广东省江门市中心医院主任医师. 1990年获四川大学放射诊断硕士学位. 主要研究方向为腹部影像鉴别诊断. 本文通信作者.E-mail: jmlws2@163.com

Segmentation of Breast Cancer on DCE-MRI Images With MRF Energy and Fuzzy Speed Function

Funds: Supported by Natural Science Foundation of Guangxi Zhuang Autonomous Region (2016GXNSFBA380160)
  • 摘要: 乳腺癌灶的精确分割是乳腺癌计算机辅助诊断的重要前提. 在动态对比增强核磁共振成像(Dynamic contrast-enhanced magnetic resonance imaging, DCE-MRI)的图像中, 乳腺癌灶具有对比度低、边界模糊及亮度不均匀等特点, 传统的活动轮廓模型方法很难取得准确的分割结果. 本文提出一种结合马尔科夫随机场(Markov random field, MRF)能量和模糊速度函数的活动轮廓模型的半自动分割方法来完成乳腺癌灶的分割, 相对于专业医生的手动分割, 本文方法具有速度快、可重复性高和分割结果相对客观等优点. 首先, 计算乳腺DCE-MRI图像的MRF能量, 以增强目标区域与周围背景的差异. 其次, 在能量图中计算每个像素点的后验概率, 建立基于后验概率驱动的活动轮廓模型区域项. 最后, 结合Gabor纹理特征、DCE-MRI时域特征和灰度特征构建模糊速度函数, 将其引入到活动轮廓模型中作为边缘检测项. 在乳腺癌灶边界处, 该速度函数趋向于零, 活动轮廓曲线停止演变, 完成对乳腺癌灶的分割. 实验结果表明, 所提出的方法有助于乳腺癌灶在DCE-MRI图像中的准确分割.
  • 乳腺癌是现代女性最常见的恶性肿瘤之一, 严重威胁着女性患者的身心健康. 如今, 乳腺癌检测技术主要有乳腺X线摄影术和乳腺磁共振成像技术. 乳腺X线摄影术的原理是利用X射线的物理特性进行投影成像, 得到全乳的二维图像, 其对非致密腺体组织的乳腺具有高敏感性, 但在具有高密度腺体组织的乳腺诊断价值不大. 乳腺动态增强核磁共振成像(Dynamic contrast-enhanced magnetic resonance imaging, DCE-MRI)是一种以病变区和正常组织中的微血管系统为生理基础的无辐射磁共振功能成像技术. 患者在注射增强对比剂后, 由于乳腺癌灶比正常组织具有更高的血管密度, 在注射增强对比剂后的前中期乳腺癌灶的强化程度更高. DCE-MRI与X线摄影相比, DCE-MRI对乳腺癌病灶的检出不受腺体密度的影响, 能够反映病变内部血流特征, 且对早期乳腺癌的敏感性、乳腺癌的分期和病变范围的一致性都优于X线摄影, 是进行乳腺癌研究的主要影像手段[1-2]. 在DCE-MRI图像上, 乳腺癌灶表现为肿块型或非肿块型. 由于乳腺癌灶大多沿腺体分布, 导致非肿块型乳腺癌灶的形态极其不规则, 肿瘤区域边界模糊且对比度低, 对乳腺癌灶的准确分割带来了极大困难[3].

    目前, 国内外研究学者已经对乳腺癌灶分割算法进行了大量的研究, 主要的方法包括边界活动轮廓模型(Active contour model, ACM)和区域活动轮廓模型. ACM主要思想是建立一个能量泛函, 在内力和外力的共同作用下使能量泛函趋近于最小值, 从而轮廓曲线演变到目标边界[4-6]. 边界ACM依据图像的梯度信息构建边缘检测函数, 在目标边界处速度函数为零,轮廓曲线停止演变; 区域ACM依据目标边界内外区域信息构建能量泛函并驱动轮廓曲线演变[7-9]. 张旭梅等[10]提出一种整合全局和局部信息的区域ACM, 可以在乳腺核磁共振图像上分割多种灰度不均匀场景下的乳腺癌灶. 虞红伟等[11]利用高斯金字塔模型得到粗尺度图像再使用局部区域ACM对粗尺度图像进行局部优化, 提高乳腺癌灶的分割精度. 但是, 乳腺DCE-MRI的灰度图像可能无法在前景和背景图像统计之间提供足够大的差异(乳腺癌灶和正常乳腺具有低对比度), 无法使区域ACM的轮廓趋近于乳腺癌灶的真实边界. 范虹等[12]利用小波多尺度分解对图像进行多尺度空间分析, 提出一种多分辨率水平集的活动轮廓模型算法应用在乳腺核磁共振图像的癌灶分割上, 避免了核磁共振图像中灰度偏移场对分割精度的影响. 然而小波变换只使用灰度变化信息, 忽略了局部区域信息, 不能准确分割具有灰度不均匀性的目标. Kuo等[13]提出基于3 维径向梯度指数和水平集演化的边界ACM并成功应用于乳腺CT图像的病灶分割. Rodtook等[14]提出一个广义梯度矢量流的边界ACM应用于乳腺超声图像分割. 但由于乳腺DCE-MRI图像具有边界模糊的特性, 基于梯度信息的ACM很难搜索到准确边界, 从而极易发生边界泄露. Xie等[15]提出一种基于学生t分布的马尔科夫随机场(Markov random field, MRF)的多相水平集聚类方法对各种噪声或无噪声图像和fMRI时间序列影像的分割. Yang等[16]利用MRF能量函数, 代数多重网格方法和稀疏场方法, 提出一种具有快速性和鲁棒性的水平集方法, 应用在合成噪声图像, 自然图像和医学等图像的分割. 上述两篇文献是将MRF能量嵌入到区域ACM的分割方法研究, 也取得了较好的分割效果, 但在乳腺DCE-MRI图像的分割中, 同时将乳腺DCE-MRI图像的时域特征和MRF能量一起融入ACM中, 分割效果可能会更好. 此外, 基于深度学习的分割方法(如卷积神经网络和全卷积神经网络)也被应用在乳腺DCE-MRI图像的病灶分割上. 基于深度学习的分割方法旨在训练图像体素的分类模型, 在病灶区域, 体素被视为阳性样本, 其余的被视为阴性样本. Kallenberg等[17]把无监督深度学习应用到乳腺癌灶分割. Rouhi等[18]提出两种基于神经网络的分割方法, 一是通过训练人工神经网络确定自动区域的增长阈值达到分割乳腺癌灶的目的; 二是通过细胞式神经网络来执行分割, 而网络的参数由遗传算法得出. Ashraf等[19]使用多通道MRF对动态观测模型的乳腺DCE-MRI数据进行自动分割.

    为了解决DCE-MRI图像中乳腺癌灶边界模糊、对比度低和亮度不均匀等问题, 本文提出一种结合MRF能量和模糊速度函数的活动轮廓模型方法来完成乳腺癌DCE-MRI图像分割. 主要工作包括以下两个方面: 1)通过对乳腺DCE-MRI图像构建二维MRF, 将图像的灰度信息转化为MRF能量, 加强乳腺癌灶和正常腺体的区分. 然后在MRF能量图上计算每个像素点的后验概率, 构建活动轮廓模型的区域项. 2)由于乳腺癌灶具有边界模糊的特点, 边界处的梯度与非边界处梯度变化并不明显, 采用梯度信息构造边缘检测函数易造成边界泄露, 本文避免采用梯度信息构建边缘检测函数, 而是利用乳腺癌灶边界处的模糊特性计算模糊速度函数, 提出结合Gabor纹理特征、DCE-MRI时域特征和灰度特征的模糊聚类算法计算模型边缘检测项的模糊速度函数. 当轮廓曲线位于乳腺癌灶的边界时, 该模糊速度函数为零, 活动轮廓曲线停止演变.

    为了精确分割DCE-MRI图像中的乳腺癌灶, 本文提出一种结合MRF能量和模糊速度函数的活动轮廓模型的乳腺癌灶分割方法. 该算法流程如图1所示.

    图 1  本文分割算法流程图
    Fig. 1  The flowchart of segmentation algorithm in this paper

    在活动轮廓模型中, 通常假设$ C $是一条封闭的轮廓曲线, 曲线$ C $将图像区域$ \Omega $划分为曲线内部区域$ \Omega_1 $和曲线外部区域$ \Omega_2 $. 活动轮廓模型的能量函数可定义为

    $$ E(C, \phi) = \mu E_1(\phi) + \lambda E_2(C) + \gamma E_3(\phi) $$ (1)

    其中, $ E_1(\phi) $是控制轮廓曲线演变方向的区域项, $ E_2(C) $是控制轮廓曲线边界演化的边缘检测项, $ E_3(\phi) $是控制轮廓曲线平滑程度的惩罚项, $ \mu $, $ \lambda $$ \gamma $分别是控制区域项、边缘检测项和惩罚项的权重参数, 在本文中$ \mu $, $ \lambda $$ \gamma $的取值都为1[20]. $ \phi $是水平集函数.

    由于DCE-MRI图像主要反映乳腺腺体供血变化情况, 癌灶组织通常以肿块形式散布于正常腺体中, 在图像上表现出空间聚类特性. MRF能量在增强目标区域与背景对比度的同时, 也考虑像素间的空间相关性[21-22]. 这一特点符合乳腺癌灶在DCE-MRI图像中的空间分布情况.

    在计算MRF能量时, 通常把图像看作一个二维的MRF, 构造MRF需要建立两个场: 标记场和特征场. 标记场是用来对待测图像进行分类标记; 特征场是对已分类的区域进行特征分析. 给定大小为$ M \times N $的图像, 内部的像素点集合用$ S $表示且$ S = \{(x, y) | 1 \le x \le M, 1 \le y \le N\}; $特征场用 $A = $$ \{a_s, s \in S\} $表示, 任意像素点所对应的特征场的值为$ a_s; $标记场用$ B = \{ b_s, $$ s \in S\} $表示, 任意像素点所对应的标记场的值为$ b_s ;$相空间用$ \Lambda = \{1, 2, \cdots, L\} $表示, 即将图像分成$ L $个区域, $ \Lambda $是标记场随机变量$ b_s $的集合, 所以$ b_s \in \Lambda .$由贝叶斯理论可得最大后验概率, 可表示为

    $$ \hat{B} = \arg \max\limits_b {P(B|A)} = \arg \max\limits_b \dfrac{P(B)P(A|B)}{{P(A)}} $$ (2)

    其中, 在一幅图像中$ P(A) $是固定不变的, 所以式$ (2) $可简化为

    $$ \hat{B} = \arg \max\limits_b {P(B|A)} = \arg \max\limits_b {P(B)P(A|B) } $$ (3)

    $ \rm Harmmersley\text{-}Clifford $定理知, 式$ (3) $中的$ P(B) $可定义为

    $$ P(B) = \prod\limits_{s \in S} P(b_s) = \prod\limits_{s \in S} \frac{\exp \left(-\displaystyle\sum\limits_{s \in S} W_c\left(b_s\right)\right)}{\displaystyle\sum\limits_{b_s = 1}^L \exp \left(-\displaystyle\sum\limits_{s \in S} W_c\left(b_s\right)\right) } $$ (4)

    其中,$ W_c(b_s) $是包含有$ b_s $的势团$ c $的势函数, $ G $表示$ s $中所有势团的集合. $ W_c(b_s) $的建模(场建模)使用Potts模型. Potts模型只考虑二元势函数, 其定义为

    $$ {W_c}({b_s})={W_2}({b_x},{b_y})=\left\{{\begin{array}{*{20}{l}} {0,}&{{b_x}={b_y}}\\ {\beta,}&{{b_x}\ne{b_y}} \end{array}}\right.$$ (5)

    对于条件概率$ P(A|B) $的求解, 假设乳腺DCE-MRI图像像素的灰度服从混合高斯分布[15, 23-24], 则可得到在标记场$ B $条件下, 特征场$ A $的条件概率$ P(A|B) $, 定义为

    $$ \begin{split} P(A|B) =& \prod\limits_{s \in S} P(a_s | b_s) = \\ &\prod\limits_{s \in S} \dfrac{1}{\sqrt {2\pi } \sigma _m}\exp \left( { - \dfrac{{\left( {a_s - \mu _m } \right)^2 }}{{2\sigma _m^2 }}} \right) \end{split}$$ (6)

    其中, 参数$ \mu_m $$ \sigma_m $分别是第$ m $目标区域的均值和方差, 可通过期望最大算法进行求解得到最优参数, 如下:

    $$ \begin{split} \mu_m^{\left(t+1\right)} = \dfrac{\displaystyle\sum\limits_{s = 1}^S P\left(z_s = i | a_s, \theta_1^{\left(t\right)} \right)a_s }{\displaystyle\sum\limits_{s = 1}^S P\left(z_s = i | a_s, \theta_1^{\left(t\right)} \right)}\quad\quad\quad\quad\quad\quad\quad\;\\ \sigma_m^{\left(t+1\right)} = \dfrac{\displaystyle\sum\limits_{s = 1}^S P\left(z_s = i | a_s, \theta_1^{\left(t\right)} \right)\left(a_s - \theta_1^{\left(t\right)}\right)^2 }{\displaystyle\sum\limits_{s = 1}^S P\left(z_s = i | a_s, \theta_1^{\left(t\right)} \right)}\quad\quad\quad\\ p(z_s \!=\! i |a_s, \theta_1^{\left(t\right)}) \!=\! \dfrac{P(a_s | z_s \!=\! i, \theta_1^{\left(t\right)})P(z_s \!=\! i | \theta_1^{\left(t\right)})}{\displaystyle\sum\limits_{k = 1}^L P(a_s | z_s \!=\! k, \theta_1^{\left(t\right)})P(z_s \!=\! k | \theta_1^{\left(t\right)})} \end{split}$$ (7)

    其中, $ \theta_1\!=\! (\mu_1 \cdots \mu_L, \sigma_1 \cdots \sigma_L, P(z \!=\! 1) \cdots P(z \!=\! L)) $, $ a_s $所从属的类用$ z_s $表示, $ P(z = i) $是先验概率, 有$\sum\nolimits_{i = 1}^L P(z = i) = 1$. 结合式$ (4) $和式$ (6) $, 对式$ (3) $两边取对数, 可得

    $$ \begin{aligned} \hat{B} =& \arg \max\limits_b \left[\ln P(B)+\ln P(A|B)\right] =\\ &\arg \max\limits_b \left\{- \sum\limits_{c \in C} W(c)(b_s) \!-\! \left[ \sum\limits_{s \in S} \left(\ln \left(\sqrt{2\pi} \sigma_m \right)\right.\right.\right.\!+\\ &\left.\left.\left.\dfrac{\left( a_s - \mu _m \right)^2 }{2 \sigma _m^2 } \right)\right]\right\} = \arg \min\limits_b \left\{\sum\limits_{c \in C} W(c)(b_s)+\right.\\ &\left.\left[ \sum\limits_{s \in S} \left(\ln \left(\sqrt{2\pi} \sigma_m \right)+\dfrac{\left( a_s - \mu _m \right)^2 }{2 \sigma _m^2 } \right)\right] \right\} =\\ &\arg \min\limits_b \left(U_1 \left(B, A \right)+U_2\left(B\right)\right) \\[-15pt]\end{aligned} $$ (8)

    其中,$ U_1(B, A) $为条件能量函数(特征场能量), $ U_2(B) $为先验能量函数(标记场能量), 分别为

    $$ \begin{split} U_1(B, A) =& \ln P(A|B) = \\ & \sum\limits_{s \in S} \left( \ln \left(\sqrt{2\pi} \sigma_m \right)+\dfrac{\left( a_s - \mu _m \right)^2 }{2 \sigma _m^2 }\right) \\ U_2(B)=& \sum\limits_{c \in C} W_c(b_c) \end{split} $$ (9)

    则MRF能量可定义为

    $$ U(B, A) = U_1(B, A)+U_2(B) $$ (10)

    进一步计算MRF能量图中每一个像素点的后验概率. 假设乳腺DCE-MRI图像像素的MRF能量服从混合高斯模型, 该模型可定义为K个高斯分布的线性组合, 则

    $$ p(x) = \sum\limits_{k = 1}^K \pi _k{\rm N}(x|\mu_k, \Sigma_k) $$ (11)

    其中, $ K $取2, ${\rm N}(x|\mu_k, \Sigma_k)$表示均值为$ \mu_k $、协方差为$ \Sigma_k $的高斯分布, $ \pi_k $表示系数, 且有

    $$ \sum\limits_{k = 1}^K \pi_k = 1,\; 0 \le \pi_k \le 1 $$ (12)

    假设$ W = \{w_1, w_2, \cdots, w_n\} $是由n个MRF能量样本组成的数据集合, 这n个样本的对数似然函数可定义为

    $$ L(\theta_2) = \ln p(W|\theta_2) $$ (13)

    其中, $ \theta_2 $是参数$ \{\pi_k, \mu_k, \Sigma_k|k\in K\} $的集合, 通过期望最大法求解对数似然函数的最大值, 得到最优参数$ \theta_2 ,$表示如下:

    $$ \begin{split}& \pi_k^{\left(t+1\right)} = \dfrac {1}{n} \sum\limits_{i = 1}^n p_k\left(w_{xy}|\theta_{2k}^{\left(t\right)}\right) \\ &\mu_k^{\left(t+1\right)} = \dfrac{\displaystyle\sum\limits_{i = 1}^n p_k\left(w_{xy}|\theta_{2k}^{\left(t\right)}\right)} {\displaystyle\sum\limits_{i = 1}^n p_k\left(k|w_{xy}, \theta_{2k}^{\left(t\right)}\right)} \\ &\Sigma_k^{\left(t+1\right)} = \dfrac{\displaystyle\sum\limits_{i = 1}^n \left( w_{xy}-\mu_k^{\left(t+1\right)}\right)^2 p_k\left(k|w_{xy}, \theta_{2k}^{\left(t\right)}\right)}{\displaystyle\sum\limits_{i = 1}^n p_k\left(k|w_{xy}, \theta_{2k}^{\left(t\right)}\right)} \end{split} $$ (14)

    其中, $ w_{xy} $表示像素点$ (x, y) $处的MRF能量. 则$ w_{xy} $服从第k个高斯分量的后验概率$ p(k|w_{xy}, \theta_{2k}) $定义为

    $$ p\left(k|w_{xy}, \theta_{2k}\right) = \dfrac{\pi_kp_k\left(w_{xy}|\theta_{2k}\right)}{\displaystyle\sum\limits_{k = 1}^K \pi_kp_k\left(w_{xy}|\theta_{2k}\right)} $$ (15)

    则活动轮廓模型的区域项可表示为

    $$ \begin{split} E_1(\phi) =\ & \displaystyle\iint_\Omega p\left(1|w_{xy}, \theta\right)H\left(\phi\left(x, y\right)\right) +\\ &p\left(2|w_{xy}, \theta\right)\left[1-H\left(\phi\left(x, y\right)\right)\right]{\rm d}x{\rm d}y \end{split} $$ (16)

    其中, $ H(\cdot) $是Heaviside函数, 定义为$H(\phi) = 0.5 \times [1+$$ 0.5 \pi \arctan(\phi / \varepsilon )] $, 参数$ \varepsilon $控制$ H(\phi) $从0上升到1的快慢[25].

    图2是乳腺图像的MRF能量分布图和对应的后验概率分布图, 从图2(b)可以看出, 由于乳腺正常腺体和乳腺癌灶具有相似的灰度值, 很难从灰度图区分正常腺体和癌灶区域. 图2(c)的MRF能量分布图可区分出正常腺体和癌灶, 但癌灶区域还存在亮度不均匀的问题.图2(d)是后验概率图, 其中癌灶区域明显增强.

    图 2  乳腺DCE-MRI图像的MRF能量和后验概率
    Fig. 2  MRF energy and posterior probability of breast DCE-MRI image

    结合Gabor纹理特征、DCE-MRI时域特征和灰度特征, 计算图像像素的模糊隶属度并构建模糊速度函数, 将其作为模型的边缘检测函数. 利用Gabor函数对图像进行纹理特征的提取[26]. 图3是一个乳腺DCE-MRI图像的Gabor纹理特征的分析. 图3(a)图3(b)分别是原始的DCE-MRI图像和经Gabor滤波器后的纹理图像. 图3(c)图3(d)分别是乳腺癌灶和正常腺体像素的纹理值. 正如图3(c)图3(d)所示, 乳腺癌灶和正常腺体之间的纹理值存在着明显的差异, 因此可用Gabor特征区分乳腺癌灶和正常腺体.

    图 3  一个乳腺图像的纹理分析
    Fig. 3  Texture analysis of a breast image

    对于DCE-MRI图像, 由于乳腺癌病灶和正常组织对所注射的对比剂吸收情况有差异, 乳腺癌病灶的吸收速度大于正常组织的吸收速度. 本文选取注射增强对比剂后的前中期(T1期和T3期)计算信号增强比率作为DCE-MRI时域特征, 可表示为

    $$ Rate_{T1{\text{-}}T3} = \dfrac{f_{T3}\left(x, y\right)-f_{T1}\left(x, y\right)}{t} $$ (17)

    其中, $ f_{T1}(x, y) $$ f_{T3}(x, y) $分别是$ T_1 $期和$ T_3 $期在像素点$ (x, y) $处的灰度值, $ t $$ T_1 $期至$ T_3 $期的时间. 图4展示一个乳腺DCE-MRI图像的时域特征.图4(a)图4(b)分别是$ T_1 $期和$ T_3 $期图像, 图4(c)是DCE-MRI时域特征, 从图4(c)可以看出乳腺癌病灶的强化程度明显强于正常组织, 该特征有助于加强乳腺癌灶与正常组织的区别.

    图 4  DCE-MRI时域特征分析
    Fig. 4  Time domain characteristic analysis of DCE-MRI

    将Gabor特征、DCE-MRI时域特征和灰度值构成的特征向量, 利用模糊聚类算法计算每一个像素的隶属度, 从而计算模糊速度函数.

    模糊聚类算法属于迭代聚类方法, 通过迭代寻找最优的目标函数, 目标函数可定义为[27]

    $$ J\left(U, { G}\right) = \sum\limits_{l = 1}^n \sum\limits_{h = 1}^2 \left(u_{lh}\right)^od^2\left({ X}_l, { G}_h\right) $$ (18)

    其中, $ { X}_l $由Gabor特征、DCE-MRI时域特征和灰度特征构成的特征向量; $ U = \{u_{lh}\} $是隶属度矩阵, 其大小为2 × n, $ u_{lh} $是向量$ { X}_l $隶属$ h $类的程度; G =$\{v_1, v_2\} $是聚类中心向量; o是一个模糊控制常数; $ d({ X}_l, { G}_h) $$ { X}_l $和距离中心$ { G}_h $之间的距离, 该距离函数定义为

    $$ d\left({ X}_l, { G}_h\right) = ||{ X}_l-{ G}_h|| $$ (19)

    重复计算$ u_{lh} $$ { G}_h $, 使所有的特征被正确分离, 此时$ J(U, { G}) $目标函数达到最小值. $ u_{lh} $$ { G}_h $可定义为

    $$ \begin{split} &u_{lh} = \dfrac{1}{\sum\limits_{k = 1}^2\left[\dfrac{d\left({ X}_l, { G}_h\right)}{d\left({ X}_l, { G}_k\right)} \right]^{\frac{2}{o-1}}}, \quad h = 1, 2;l = 1, 2, \cdots, n \\ &{ G}_h = \dfrac{\sum\limits_{l = 1}^n\left(u_{lh}\right)^o{ X}_l}{\sum\limits_{l = 1}^n\left(u_{lh}\right)^o}, \quad\quad\quad\quad\;\; h = 1, 2\\[-30pt] \end{split} $$ (20)

    然后, 根据$ u_{lh} $构造隶属乳腺癌灶的模糊隶属度矩阵$ E $. 最后, 为了消除模糊隶属度矩阵中存在的噪声和保留边界信息, 采用结构元素$ R $$ E $进行模糊开运算, 得到新的模糊隶属度矩阵$ Z $. 模糊开运算定义为

    $$ Z = D[F(E, R), R] $$ (21)

    其中, $ F $$ D $分别为模糊腐蚀运算算子和模糊膨胀运算算子, 其分别定义为

    $$ \begin{split} & F(E, R)(x, y) = \\ &\qquad{\rm min}\left\{1, \min\limits_{\left(j, t\right)\in {\bf R}}\left[1+u_E\left(x-j, y-t\right)-u_R\left(j, t\right)\right] \right\} \\ & D(E, R)(x, y)= \\ &\qquad{\rm min}\left\{0, \min\limits_{\left(j, t\right)\in {\bf R}}\left[u_E\left(x-j, y-t\right)+u_R\left(j, t\right)-1\right] \right\} \end{split} $$ (22)

    其中, $ u_E(x-j, y-t) $是模糊矩阵$ E $中的元素, $ u_R(j, t) $是结构元素$ R $中的元素. 则像素点$ (x, y) $的模糊隶属度可表示为$ Z(x, y)\in(0, 1) $.

    在理想的情况下, 乳腺癌灶的模糊隶属度大于0.5, 在乳腺癌灶边界处模糊隶属度近似等于0.5, 而背景的模糊隶属度小于0.5. 模型的模糊速度函数需要有以下特点: 在乳腺癌灶边界处, 模糊速度函数等于零; 在远离乳腺癌灶边界处, 模糊速度函数越来越大; 在趋近乳腺癌灶边界处, 模糊速度函数越来越小. 则像素$ (x, y) $的模糊速度函数可定义为

    $$ V(x, y) = {\rm e}^{\alpha\left[Z\left(x, y\right)-0.5\right]^2-1} $$ (23)

    其中, $ \alpha $为常数. 当轮廓曲线处于对象的边界时, 模糊速度函数$ V(x, y) $接近于零. 当轮廓曲线远离边界时, 模糊速度函数$ V(x, y) $越来越大. 图5显示某个图像AB段模糊速度函数的变化趋势, 可明显看到在边界AB处模糊速度函数约等于零.

    图 5  模糊速度函数
    Fig. 5  Fuzzy speed function

    在所提出的活动轮廓模型中, 边缘检测项$ E_2(C) $可定义为

    $$ E_2(C) = \iint_\Omega V\left(x, y\right)\delta \left(\phi \right)\left| \nabla \phi \right| {\rm d}x{\rm d}y $$ (24)

    其中, $ \delta(\phi) = \dfrac{\varepsilon}{\pi}\dfrac{1}{\varepsilon^2+\phi^2} $$ H(\phi) $的导数.

    1.4.1   构建能量泛函

    将式(16)和式(24)代入式(1)可得到最终的活动轮廓模型的能量泛函

    $$ \begin{split} E(\phi) =\ & \mu \iint_\Omega p\left(1|w_{xy}, \theta\right)H\left(\phi\left(x, y\right)\right)+ \\& p\left(2|w_{xy}, \theta\right)\left[1-H\left(\phi\left(x, y\right)\right)\right] {\rm d}x{\rm d}y+ \\& \lambda \iint_\Omega V\left(x, y\right)\delta \left(\phi \right)\left| \nabla \phi \right|{\rm d}x{\rm d}y+ \\ & \gamma \iint_\Omega \frac{1}{2}\left(\left| \nabla \phi \right|-1\right)^2 {\rm d}x{\rm d}y \end{split} $$ (25)

    其中, $ \iint_\Omega \frac{1}{2}\left(\left| \nabla \phi \right|-1\right)^2 {\rm d}x{\rm d}y $是对水平集函数$ \phi $的惩罚项, 可避免曲线演变过程中对$ \phi $进行重新初始化, $ \phi $需要满足$ |\nabla \phi| = 1 ,$其中$ |\nabla \phi| = \sqrt{\phi_x^2+\phi_y^2} .$

    1.4.2   求解能量泛函

    乳腺癌灶的分割可归结为最小化泛函$ E(\phi) $. 首先对能量泛函$ E(\phi) $关于$ \phi $求变分, 为了简化表示, $ p(1|w_{xy}, \theta) $表示为$ p_1 $, 相应的$ p(2|w_{xy}, \theta) $表示为$ p_2 $, 忽略变量$ (x, y) $, $ E(\phi) $可简化为

    $$ \begin{array}{l}F(\phi, \phi_x, \phi_y) = \mu p_1 H\left(\phi\right)+\mu p_2 \left(1-H\left(\phi\right)\right)+ \\ \quad\quad\quad \lambda V\delta \left(\phi \right)\sqrt{\phi_x^2+\phi_y^2}+\dfrac{1}{2}\gamma\left(\sqrt{\phi_x^2+\phi_y^2}-1\right) \end{array} $$ (26)

    其中, 最小化式(25)对应于求解如下偏微分方程

    $$ \frac{\partial F}{\partial \phi}-\frac{\rm d}{{\rm d}x}\left(\frac{\partial F}{\partial \phi_x}\right)-\frac{\rm d}{{\rm d}y}\left(\frac{\partial F}{\partial \phi_y}\right) = 0 $$ (27)

    其中, $ F $$ \phi $, $ \phi_x $, $ \phi_y $的偏导数可以表示为

    $$ \begin{split}&\dfrac{\partial F}{\partial \phi} = \mu\delta \left(\phi \right)\left(p_1-p_2\right)+\lambda V\delta^\prime \left(\phi\right)\sqrt{\phi_x^2+\phi_y^2} \\ & \frac{\partial F}{\partial \phi_x} = \lambda V\delta\left(\phi\right)\dfrac{\phi_x}{\sqrt{\phi_x^2+\phi_y^2}}+\gamma\left(\phi_x-\dfrac{\phi_x}{\sqrt{\phi_x^2+\phi_y^2}}\right) \\ &\dfrac{\partial F}{\partial \phi_y} = \lambda V\delta\left(\phi\right)\dfrac{\phi_y}{\sqrt{\phi_x^2+\phi_y^2}}+\gamma\left(\phi_y-\frac{\phi_y}{\sqrt{\phi_x^2+\phi_y^2}}\right) \end{split}$$ (28)

    进一步对式(28)求导

    $$ \begin{array}{l} \dfrac{\rm d}{{\rm d}x}\!\!\left(\!\dfrac{\partial F}{\partial \phi}\!\right) \! =\! \lambda\left[\!V\delta^\prime\!\left(\phi\right)\dfrac{\phi_x^2}{\sqrt{\phi_x^2\!+\!\phi_y^2}}\!+\!V\delta\!\left(\phi\right)\!\dfrac{\phi_x}{\sqrt{\phi_x^2\!+\!\phi_y^2}}\!\right]\!+\!\\ \;\;\;\;\;\;\;\;\;\;\gamma\dfrac{\rm d}{{\rm d}x}\left(\phi_x\left(1-\dfrac{1}{\sqrt{\phi_x^2+\phi_y^2}}\right)\right)\\ \dfrac{\rm d}{{\rm d}x}\!\!\left(\!\dfrac{\partial F}{\partial\phi}\!\right)\!=\! \lambda \left[\!V\delta^\prime \!\left(\phi\right)\dfrac{\phi_y^2}{\sqrt{\phi_x^2\!+\!\phi_y^2}}\!+\!V\delta \!\left(\phi\right)\!\dfrac{\phi_y}{\sqrt{\phi_x^2\!+\!\phi_y^2}}\!\right]\!+\!\\ \;\;\;\;\;\;\;\;\;\;\gamma\dfrac{\rm d}{{\rm d}y}\left(\phi_y\left(1-\dfrac{1}{\sqrt{\phi_x^2+\phi_y^2}}\right)\right) \end{array} $$ (29)

    同时结合

    $$ \begin{split}&\dfrac{\rm d}{{\rm d}x}\left(\dfrac{\phi_x}{|\nabla \phi|}\right)+\dfrac{\rm d}{{\rm d}y}\left(\dfrac{\phi_y}{|\nabla \phi|}\right) = {\rm div}\dfrac{\nabla \phi}{|\nabla \phi|} \\ &\dfrac{\phi_x^2}{|\nabla \phi|}+\dfrac{\phi_y^2}{|\nabla \phi|} = |\nabla \phi| \end{split} $$ (30)

    将式(28)~(30)代入式(27), 可得到$ F $关于水平集函数$ \phi $的欧拉–拉格朗日方程

    $$ \begin{split}&\mu\delta\left(\phi\right)\left(p_1-p_2\right)-\lambda\delta\left(\phi\right)\left({\rm div}\left(V\frac{\nabla \phi}{|\nabla \phi|}\right)\right)-\\&\quad\;\;\; \gamma\left(\Delta\phi-{\rm div}\frac{\nabla \phi}{|\nabla \phi|}\right) = 0 \end{split} $$ (31)

    进一步由变分原理, 可得$ \phi $的梯度下降流

    $$ \begin{split} \frac{\partial \phi}{\partial t} =& -\mu\delta\left(\phi\right)\left(p_1-p_2\right)+\lambda\delta\left(\phi\right)\left({\rm div}\left(V\frac{\nabla \phi}{|\nabla \phi|}\right)\right)+\\ &\gamma\left(\Delta\phi-{\rm div}\frac{\nabla \phi}{|\nabla \phi|}\right)\\[-20pt] \end{split} $$ (32)

    在对比算法中, 金标准由影像科医生提供. 医生首先在DCE-MRI图像中执行自动阈值的预分割, 然后通过观察DCE-MRI剪影图像, 再对肿瘤区域进行手动调整, 最终获得的肿瘤区域作为金标准. 分割精度通过对算法分割结果与医生手动分割结果进行比较得到, 分别采用阳性率(True positive ratio)、假阳性率(False positive ratio)和相似度(Similarity degree)三个指标来进行评价

    $$ \begin{split}& R_{\rm{TP}} = \dfrac{|A_m \cap A_a|}{A_m} \\ & R_{\rm{FP}} = \dfrac{|A_m \cup A_a-A_m|}{A_m} \\& D_{\rm{S}} = \dfrac{|A_m \cap A_a|}{|A_m \cup A_a} \end{split} $$ (33)

    其中, $ A_m $是医生手动分割区域, $ A_a $是分割方法获取的区域. $ R_{\rm{TP}} $反映的是算法分割区域对医生手动分割区域的覆盖率, 值越大, 说明分割效果越好. $ R_{\rm{FP}} $反映的是算法分割区域的背景区域对医生手动分割区域的占比, 值越小, 说明算法分割的区域包含更少的背景区域. 当$D_{\rm{S} }$值比较大时, 说明分割方法获取的区域与医生手动分割区域非常形似.

    为了验证本文所提方法的有效性, 采用中山大学附属江门市中心医院的临床数据进行验证分析. 本次研究所使用的仪器设备主要为西门子MAGNETO MESSENZA 1.5T MRA; 同时还有4通道乳腺专用的相控阵表面线圈. 所有患者均处于俯卧位. 使用T1加权脂肪抑制快速激发梯度回波采集双乳横轴位的DCE图像(TR/TE: 4.6/1.7; 翻转角: 7度; FOV: 280 × 340 mm; 矩阵大小: 280 × 340; 截面厚度/间隙: 1.0/0 mm; 扫描时间: 75 s).

    将所提出方法与结合二值水平集和形态学的分割方法[28]、基于区域的ACM、基于局部高斯能量的ACM、基于MRF的区域ACM[29]、基于边界的ACM、深度学习框架中的U-net和V-net网络进行对比分析. V-net和U-net使用的都是基础模型, V-net的参数设置如下: batch_size$ = 10 $, 学习率$ = 0.001 $, 迭代次数$ = 400\;000 $, 损失函数为dice coefficient, 迭代终止条件是迭代次数; U-net的参数设置如下: batch_size = 10, 学习率 = 0.00001, 迭代次数 =10 000, 损失函数为binary crossentropy, 迭代终止条件是迭代次数. 其中共使用2 255幅乳腺DCE-MRI图像进行V-net和U-net网络的训练, 测试集共50 幅乳腺DCE-MRI图像. V-net和U-net在训练集中的表现如表1所示, 可以看出U-net出现了过分割的现象, V-net的表现较好.

    表 1  V-net和U-net训练集结果
    Table 1  Results of V-net and U-net in training sets
    分割
    算法
    $R_{\rm{TP}}$
    (均值±标准差)
    $R_{\rm{FP}}$
    (均值±标准差)
    $D_{\rm{S}}$
    (均值±标准差)
    U-net0.9706±0.11890.6673±0.66620.6252±0.1456
    V-net0.9408±0.06370.0732±0.12870.8836±0.0917
    下载: 导出CSV 
    | 显示表格

    首先验证算法对肿块型乳腺癌的分割性能. 由影像科医生选择5 个典型的肿块型乳腺癌DCE-MRI图像进行分析. 如图6所示, 所提出方法对乳腺癌灶的分割结果更接近于医生手动分割的结果, 优于其余7种方法. 图6(b)展示结合二值水平集和形态学分割方法的结果, 可以发现受周围正常腺体的影响, 得到的乳腺癌灶边界朝周围腺体发生泄漏. 图6(c)图6(d)展示基于边界的ACM和基于局部高斯能量的ACM的结果, 可以看出受到乳房边缘和周围腺体的影响, 分割结果无法正常收敛甚至把乳房边缘错认为乳腺癌灶边界. 图6(e)展示基于区域的ACM的结果, 由于癌灶内部存在亮度不均匀, 也影响了最终分割的准确性. 图6(f)展示基于MRF的区域ACM的结果, 由于癌灶边界具有模糊特性, 分割结果存在欠分割问题. 图6(g)图6(h)分别展示U-net和V-net的分割结果, 也存在误分割的情况. 对比图6(a)图6(i), 可以明显看出所提出方法得到的分割结果更接近于医生手动分割结果.

    图 6  肿块形乳腺癌灶分割结果((a)医生手动分割; (b)结合二值水平集和形态学的分割方法; (c)基于边界的ACM; (d)基于局部高斯能量的ACM; (e)基于区域的ACM; (f)基于MRF的区域ACM; (g) U-net; (h) V-net; (i)本文方法)
    Fig. 6  Tumor-shaped breast lesion segmentation results ((a) Expert manual segmentation; (b) Segmentation method based on binary level set and morphological operation; (c) Edge-based ACM; (d) ACM based on local Gaussian energy; (e) Region-based ACM; (f) Regional ACM based on MRF; (g) U-net; (h) V-net; (i) The method in this paper)

    部分乳腺癌在DCE-MRI图像上表现为非肿块型, 且肿瘤区域与正常腺体组织往往交错存在. 为了进一步验证算法性能, 由影像科医生从临床数据中选择了10 个典型的非肿块型乳腺癌进行分割. 从图7(b)图7(c)可以看出, 由于乳腺癌灶和周围腺体组织具有极相近的对比度, 结合二值水平集和形态学的分割方法和基于边界的ACM的结果包含了较多的正常腺体, 算法不能稳定收敛, 出现了边界泄漏问题. 图7(d)展示基于局部高斯能量的ACM, 也无法得到正确的分割结果. 图7(e)展示基于区域的ACM和图7(f)展示基于MRF的区域ACM, 也无法分割低对比度下的乳腺癌灶, 发生了边界泄漏. 图7(g)图7(h)分别展示U-net和V-net的分割结果.图7(a)图7(i)对比, 只有本文所提出方法提取出的乳腺癌灶更接近于医生手动分割结果.

    图 7  非肿块形乳腺癌灶分割结果((a)医生手动分割; (b)结合二值水平集和形态学的分割方法; (c)基于边界的ACM; (d)基于局部高斯能量的ACM; (e)基于区域的ACM; (f)基于MRF的区域ACM; (g) U-net; (h) V-net; (i)本文方法)
    Fig. 7  Non tumor-shaped breast lesion segmentation results ((a) Expert manual segmentation; (b) Segmentation method based on binary level set and morphological operation; (c) Edge-based ACM; (d) ACM based on local gaussian energy; (e) Region-based ACM; (f) Regional ACM based on MRF; (g) U-net; (h) V-net; (i) The method in this paper)

    对于基于局部高斯能量的ACM方法不能精确分割乳腺癌灶的原因主要是: 将非负的核函数$ K(·) $引入到ACM中. 当局部窗口的大小取值不合适时, 非负的核函数$ K(·) $不能显著的减小为零, 因此基于局部高斯能量的ACM会引起轮廓边界的泄露; 对于基于边界的ACM, 当图像中的癌灶有较深的凹陷边界时, 曲线$ C $在这一局部的曲率必将是负的, 模型可能使轮廓曲线$ C $停止在使能量泛函局部最小的状态.

    对于基于区域ACM和结合二值水平集和形态学的分割方法不能精确分割乳腺癌灶的原因主要是: 1)在乳腺癌灶和腺体之间具有非常相近的灰度值, 在轮廓曲线演变时仅使用了图像的灰度信息而并没有考虑乳腺癌灶的局部形状信息, 因此发生边界泄露. 2)乳腺癌灶成分和周围腺体之间具有较低的对比度, 使用梯度计算边界检测函数, 无法准确检测低对比度下的癌灶边界, 从而导致边界泄露. 而基于MRF的区域ACM并没有用到边界信息, 所以对于模糊边界的癌灶存在欠分割. 另外, 深度学习框架U-net和V-net方法也取得了不错的分割结果, 但是在有些图像样本上也出现了边界泄露问题, 特别是在非肿块形的乳腺癌灶中分割结果不理想. 原因可能是: 1)基于深度学习的分割算法是基于数据驱动的分割方法, 无法加入有效的临床先验信息(如DCE-MRI图像的动态信息), 导致分割效果有效; 2)乳腺DCE-MRI图像中存在严重的类平衡问题(病灶区域中的体素数(正类) 远小于背景中的体素数(负类)), 并且图像中混杂着多种背景(如血管结构, 正常腺体等), 导致算法的鲁棒性较差.

    为了进一步验证算法性能, 继续选择了50例乳腺癌临床样本. 分别用$ R_{\rm{TP}} $$ R_{\rm{FP}} $$ D_{\rm{S}} $的平均值作为评价指标, 对比了7种不同分割方法在50个临床样本中的分割结果. 结果如表2所示, 所提出方法能够较好地将乳腺癌灶分割出来, 结合二值水平集和形态学的分割方法, 基于边界的ACM, 基于局部高斯能量的ACM, 基于区域ACM方法和U-net方法虽然得到了更好的$ R_{\rm{TP}} $值, 但$ R_{\rm{FP}} $值太大且$ D_{\rm{S}} $值太小, 说明分割得到的结果含有大量背景. 虽然所提出方法得到的 $ R_{\rm{TP}} $值只达到了92.97%, 但在$ R_{\rm{FP}} $指标中, 本文所提出方法分割得到的结果包含的背景最少$( R_{\rm{FP}} $= 3.37%). 在$ D_{\rm{S}} $指标中, 所提出方法与医生手动分割结果外形最相似$( D_{\rm{S}} =$89.96%).

    表 2  50个临床样本分割结果
    Table 2  Segmentation results of 50 clinical samples
    分割算法$R_{\rm{TP}}$(均值±标准差)$R_{\rm{FP}}$(均值±标准差)$D_{\rm{S}}$(均值±标准差)
    结合二值水平集和形态学的分割方法0.9759±0.02760.9100±1.61390.6910±0.2537
    基于边界的ACM0.9676±0.05510.6270±0.42380.5146±0.2272
    基于局部高斯能量的ACM0.9548±0.08180.8368±1.02170.6147±0.2040
    基于区域的ACM0.9892±0.01351.0604±1.01930.5857±0.2357
    基于MRF的区域ACM0.8709±0.11780.2018±0.32710.7524±0.1438
    U-net0.9389±0.11560.5988±0.37720.6143±0.1456
    V-net0.8762±0.14080.1131±0.09610.7920±0.1402
    本文方法0.9297±0.04130.0337±0.02220.8996±0.0410
    下载: 导出CSV 
    | 显示表格

    为了分析算法在医学应用的可行性, 对算法进行时间复杂度的分析. 表3展示了各个对比算法的时间复杂度和分割50 个临床样本的平均执行时间(U-net和V-net只考虑分割50 个临床样本的时间). 从表3中看出, 本文方法的时间复杂度为$ {\rm O}(n^2) $, 平均执行时间为4.6575 s, 耗时略高于其余方法, 但在可接受范围之内.

    表 3  算法的时间复杂度
    Table 3  Time complexity of the algorithm
    分割算法时间复杂度平均执行时间 (s)
    结合二值水平集和形态学的分割方法${\rm O}(n)$3.3500
    基于边界的ACM${\rm O}(n)$3.2550
    基于局部高斯能量的ACM${\rm O}(n)$3.6790
    基于区域的ACM${\rm O}(n)$1.1188
    基于MRF的区域ACM${\rm O}(n^2)$2.7925
    U-net${\rm O}(n)$1.6102
    V-net${\rm O}(n)$3.3669
    本文方法${\rm O}(n^2)$4.6575
    下载: 导出CSV 
    | 显示表格

    为了在DCE-MRI图像上精确分割乳腺癌灶, 提出一种结合MRF能量和模糊速度函数的活动轮廓模型方法. 通过对乳腺DCE-MRI图像构建二维MRF, 将图像的灰度信息转化为MRF能量, 并在MRF能量的基础上计算每个像素点的后验概率, 构建活动轮廓模型的区域项, 在增强乳腺癌灶和正常组织差异时, 也考虑了乳腺癌灶在DCE-MRI上的空间聚类特性. 对于乳腺癌灶边缘模糊的问题, 放弃了基于梯度的边界检测函数构建方式, 转而结合乳腺DCE-MRI时域特征、灰度特征和Gabor纹理特征来计算模糊速度函数, 并将其作为边缘检测函数. 当轮廓曲线位于乳腺癌灶的边界时, 该模糊速度函数为零, 活动轮廓模型停止演变. 通过对临床数据的验证分析, 发现本文所提出方法从乳腺DCE-MRI图像中分割得到的边界更接近于医生给出的结果. 下一步, 我们将探索将深度学习方法和活动轮廓模型方法进行结合, 拟进一步提高乳腺癌DCE-MRI图像分割的准确性.

  • 图  1  本文分割算法流程图

    Fig.  1  The flowchart of segmentation algorithm in this paper

    图  2  乳腺DCE-MRI图像的MRF能量和后验概率

    Fig.  2  MRF energy and posterior probability of breast DCE-MRI image

    图  3  一个乳腺图像的纹理分析

    Fig.  3  Texture analysis of a breast image

    图  4  DCE-MRI时域特征分析

    Fig.  4  Time domain characteristic analysis of DCE-MRI

    图  5  模糊速度函数

    Fig.  5  Fuzzy speed function

    图  6  肿块形乳腺癌灶分割结果((a)医生手动分割; (b)结合二值水平集和形态学的分割方法; (c)基于边界的ACM; (d)基于局部高斯能量的ACM; (e)基于区域的ACM; (f)基于MRF的区域ACM; (g) U-net; (h) V-net; (i)本文方法)

    Fig.  6  Tumor-shaped breast lesion segmentation results ((a) Expert manual segmentation; (b) Segmentation method based on binary level set and morphological operation; (c) Edge-based ACM; (d) ACM based on local Gaussian energy; (e) Region-based ACM; (f) Regional ACM based on MRF; (g) U-net; (h) V-net; (i) The method in this paper)

    图  7  非肿块形乳腺癌灶分割结果((a)医生手动分割; (b)结合二值水平集和形态学的分割方法; (c)基于边界的ACM; (d)基于局部高斯能量的ACM; (e)基于区域的ACM; (f)基于MRF的区域ACM; (g) U-net; (h) V-net; (i)本文方法)

    Fig.  7  Non tumor-shaped breast lesion segmentation results ((a) Expert manual segmentation; (b) Segmentation method based on binary level set and morphological operation; (c) Edge-based ACM; (d) ACM based on local gaussian energy; (e) Region-based ACM; (f) Regional ACM based on MRF; (g) U-net; (h) V-net; (i) The method in this paper)

    表  1  V-net和U-net训练集结果

    Table  1  Results of V-net and U-net in training sets

    分割
    算法
    $R_{\rm{TP}}$
    (均值±标准差)
    $R_{\rm{FP}}$
    (均值±标准差)
    $D_{\rm{S}}$
    (均值±标准差)
    U-net0.9706±0.11890.6673±0.66620.6252±0.1456
    V-net0.9408±0.06370.0732±0.12870.8836±0.0917
    下载: 导出CSV

    表  2  50个临床样本分割结果

    Table  2  Segmentation results of 50 clinical samples

    分割算法$R_{\rm{TP}}$(均值±标准差)$R_{\rm{FP}}$(均值±标准差)$D_{\rm{S}}$(均值±标准差)
    结合二值水平集和形态学的分割方法0.9759±0.02760.9100±1.61390.6910±0.2537
    基于边界的ACM0.9676±0.05510.6270±0.42380.5146±0.2272
    基于局部高斯能量的ACM0.9548±0.08180.8368±1.02170.6147±0.2040
    基于区域的ACM0.9892±0.01351.0604±1.01930.5857±0.2357
    基于MRF的区域ACM0.8709±0.11780.2018±0.32710.7524±0.1438
    U-net0.9389±0.11560.5988±0.37720.6143±0.1456
    V-net0.8762±0.14080.1131±0.09610.7920±0.1402
    本文方法0.9297±0.04130.0337±0.02220.8996±0.0410
    下载: 导出CSV

    表  3  算法的时间复杂度

    Table  3  Time complexity of the algorithm

    分割算法时间复杂度平均执行时间 (s)
    结合二值水平集和形态学的分割方法${\rm O}(n)$3.3500
    基于边界的ACM${\rm O}(n)$3.2550
    基于局部高斯能量的ACM${\rm O}(n)$3.6790
    基于区域的ACM${\rm O}(n)$1.1188
    基于MRF的区域ACM${\rm O}(n^2)$2.7925
    U-net${\rm O}(n)$1.6102
    V-net${\rm O}(n)$3.3669
    本文方法${\rm O}(n^2)$4.6575
    下载: 导出CSV
  • [1] Fitzmaurice C, Akinyemiju T F, Al Lami F H, Alam T, Alizadeh-Navaei R, Allen C, et al. Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 29 cancer groups, 1990 to 2016: A systematic analysis for the global burden of disease study. JAMA Oncology, 2018, 4(11): 1553−1568 doi: 10.1001/jamaoncol.2018.2706
    [2] Uematsu T, Kasami M, Watanabe J. Is evaluation of the presence of prepectoral edema on T2-weighted with fat-suppression 3 T breast MRI a simple and readily available noninvasive technique for estimation of prognosis in patients with breast cancer? Breast Cancer, 2014, 21(6): 684−692
    [3] Cheon H, Kim H J, Lee S M, Cho S H, Shin K M, Kim G C, et al. Preoperative MRI features associated with lymphovascular invasion in node-negative invasive breast cancer: A propensity-matched analysis. Journal of Magnetic Resonance Imaging, 2017, 46(4): 1037−1044 doi: 10.1002/jmri.25710
    [4] Fasihi M S, Mikhael W B. Overview of current biomedical image segmentation methods. In: Proceedings of the 2016 International Conference on Computational Science and Computational Intelligence. Las Vegas, USA: IEEE, 2016. 803−808
    [5] van Dijk N P, Maute K, Langelaar M, van Keulen F. Level-set methods for structural topology optimization: A review. Structural and Multidisciplinary Optimization, 2013, 48(3): 437−472 doi: 10.1007/s00158-013-0912-y
    [6] 马超, 刘亚淑, 骆功宁, 王宽全. 基于级联随机森林与活动轮廓的3D MR图像分割. 自动化学报, 2019, 45(5): 1004−1014

    Ma Chao, Liu Ya-Shu, Luo Gong-Ning, Wang Kuan-Quan. Combining concatenated random forests and active contour for the 3D MR images segmentation. Acta Automatica Sinica, 2019, 45(5): 1004−1014
    [7] 孙文燕, 董恩清, 曹祝楼, 郑强. 一种基于模糊主动轮廓的鲁棒局部分割方法. 自动化学报, 2017, 43(4): 611−621

    Sun Wen-Yan, Dong En-Qing, Cao Zhu-Lou, Zheng Qiang. A robust local segmentation method based on fuzzy-energy based active contour. Acta Automatica Sinica, 2017, 43(4): 611−621
    [8] Zhang K H, Zhang L, Lam K M, Zhang D. A level set approach to image segmentation with intensity inhomogeneity. IEEE Transactions on Cybernetics, 2016, 46(2): 546−557 doi: 10.1109/TCYB.2015.2409119
    [9] Rouhi R, Jafari M. Classification of benign and malignant breast tumors based on hybrid level set segmentation. Expert Systems with Applications, 2016, 46: 45−59 doi: 10.1016/j.eswa.2015.10.011
    [10] 张旭梅, 范虹, 乔柱. 融合全局和局部信息的水平集乳腺MR图像分割. 计算机应用研究, 2015, 32(1): 307−311 doi: 10.3969/j.issn.1001-3695.2015.01.072

    Zhang Xu-Mei, Fan Hong, Qiao Zhu. Level set method combined global and local information and its application to breast MRI. Application Research of Computers, 2015, 32(1): 307−311 doi: 10.3969/j.issn.1001-3695.2015.01.072
    [11] 虞红伟, 厉力华, 徐伟栋, 刘伟, 张娟, 邵国良. 一种基于水平集的多尺度乳腺肿块分割方法. 仪器仪表学报, 2010, 31(6): 1418−1423

    Yu Hong-Wei, Li Li-Hua, Xu Wei-Dong, Liu Wei, Zhang Juan, Shao Guo-Liang. Multi-scale approach for breast mass segmentation using level set method. Chinese Journal of Scientific Instrument, 2010, 31(6): 1418−1423
    [12] 范虹, 朱艳春, 王芳梅, 张旭梅. 多分辨率水平集算法的乳腺MR图像分割. 物理学报, 2014, 63(11): Article No.118701

    Fan Hong, Zhu Yan-Chun, Wang Fang-Mei, Zhang Xu-Mei. Segmentation of breast MR images based on multiresolution level set algorithm. Acta Physica Sinica, 2014, 63(11): Article No.118701
    [13] Kuo H C, Giger M L, Reiser I, Boone J M, Lindfors K K, Yang K, et al. Level set segmentation of breast masses in contrast-enhanced dedicated breast CT and evaluation of stopping criteria. Journal of Digital Imaging, 2014, 27(2): 237−247 doi: 10.1007/s10278-013-9652-1
    [14] Rodtook A, Makhanov S S. Multi-feature gradient vector flow snakes for adaptive segmentation of the ultrasound images of breast cancer. Journal of Visual Communication and Image Representation, 2013, 24(8): 1414−1430 doi: 10.1016/j.jvcir.2013.09.009
    [15] Xie Q Y, Pan X, Zhu H Q. A multiphase level set clustering approach using MRF-based Student's-t mixture model. In: Proceedings of the 8th International Conference on Wireless Communications and Signal Processing. Yangzhou, China: IEEE, 2016. 1−5
    [16] Yang X, Gao X B, Tao D C, Li X L, Li J. An efficient MRF embedded level set method for image segmentation. IEEE Transactions on Image Processing, 2015, 24(1): 9−21 doi: 10.1109/TIP.2014.2372615
    [17] Kallenberg M, Petersen K, Nielsen M, Ng A Y, Diao P F, Igel C, et al. Unsupervised deep learning applied to breast density segmentation and mammographic risk scoring. IEEE Transactions on Medical Imaging, 2016, 35(5): 1322−1331 doi: 10.1109/TMI.2016.2532122
    [18] Rouhi R, Jafari M, Kasaei S, Keshavarzian P. Benign and malignant breast tumors classification based on region growing and CNN segmentation. Expert Systems with Applications, 2015, 42(3): 990−1002 doi: 10.1016/j.eswa.2014.09.020
    [19] Ashraf A B, Gavenonis S, Daye D, Mies C, Feldman M, Rosen M, et al. A multichannel Markov random field approach for automated segmentation of breast cancer tumor in DCE-MRI data using kinetic observation model. In: Proceedings of the 14th International Conference on Medical Image Computing and Computer-Assisted Intervention. Toronto, Canada: Springer, 2011. 546−553
    [20] Achuthan A, Rajeswari M, Ramachandram D, Aziz M E, Shuaib I L. Wavelet energy-guided level set-based active contour: A segmentation method to segment highly similar regions. Computers in Biology and Medicine, 2010, 40(7): 608−620 doi: 10.1016/j.compbiomed.2010.04.005
    [21] Shahvaran Z, Kazemi K, Helfroush M S, Jafarian N, Noorizadeh N. Variational level set combined with Markov random field modeling for simultaneous intensity non-uniformity correction and segmentation of MR images. Journal of Neuroscience Methods, 2012, 209(2): 280−289 doi: 10.1016/j.jneumeth.2012.06.012
    [22] 宋艳涛, 纪则轩, 孙权森. 基于图像片马尔科夫随机场的脑MR图像分割算法. 自动化学报, 2014, 40(8): 1754−1763

    Song Yan-Tao, Ji Ze-Xuan, Sun Quan-Sen. Brain MR image segmentation algorithm based on Markov random field with image patch. Acta Automatica Sinica, 2014, 40(8): 1754−1763
    [23] Li B, Chen K, Peng G M, Guo Y X, Tian L F, Ou S X, et al. Segmentation of ground glass opacity pulmonary nodules using an integrated active contour model with wavelet energy-based adaptive local energy and posterior probability-based speed function. Materials Express, 2016, 6(4): 317−327 doi: 10.1166/mex.2016.1311
    [24] 张喆, 韩德强, 杨艺. 基于证据马尔可夫随机场模型的图像分割. 控制与决策, 2017, 32(9): 1607−1613

    Zhang Zhe, Han De-Qiang, Yang Yi. Image segmentation based on evidential Markov random field model. Control and Decision, 2017, 32(9): 1607−1613
    [25] Li C M, Kao C Y, Gore J C, Ding Z H. Minimization of region-scalable fitting energy for image segmentation. IEEE Transactions on Image Processing, 2008, 17(10): 1940−1949 doi: 10.1109/TIP.2008.2002304
    [26] Khan S, Hussain M, Aboalsamh H, Mathkour H, Bebis G, Zakariah M. Optimized Gabor features for mass classification in mammography. Applied Soft Computing, 2016, 44: 267−280 doi: 10.1016/j.asoc.2016.04.012
    [27] Kwedlo W. A new random approach for initialization of the multiple restart EM algorithm for Gaussian model-based clustering. Pattern Analysis and Applications, 2015, 18(4): 757−770 doi: 10.1007/s10044-014-0441-3
    [28] 郑强, 董恩清. 一种新的基于二值水平集和形态学的局部分割方法. 电子与信息学报, 2012, 34(2): 375−381

    Zheng Qiang, Dong En-Qing. A new local segmentation method based on binary level set and morphological operation. Journal of Electronics and Information Technology, 2012, 34(2): 375−381
    [29] Xu J, Monaco J P, Madabhushi A. Markov random field driven region-based active contour model (MaRACel): Application to medical image segmentation. In: Proceedings of the 13th International Conference on Medical Image Computing and Computer-Assisted Intervention. Beijing, China: Springer, 2010. 197−204
  • 期刊类型引用(14)

    1. 党红宇,唐利明,徐雅雅. 基于Bayesian最小误判的活动轮廓模型. 内江师范学院学报. 2024(04): 53-64 . 百度学术
    2. 孙淑婷,刘铖枨,周广茵,韩锐,陈立超,羊月褀,许玥. 图像分割算法在医学图像中的应用综述. 现代仪器与医疗. 2024(02): 59-68 . 百度学术
    3. 蒲秋梅,殷帅,李正茂,赵丽娜. U型卷积网络在乳腺医学图像分割中的研究综述. 计算机科学与探索. 2024(06): 1383-1403 . 百度学术
    4. 刘侠,吕志伟,李博,王波,王狄. 基于多尺度残差双域注意力网络的乳腺动态对比度增强磁共振成像肿瘤分割方法. 电子与信息学报. 2023(05): 1774-1785 . 百度学术
    5. 阮旺,郝国生,王霞,胡晓婷,杨子豪. 面向目标识别的特征融合模糊模型及其应用. 计算机科学. 2023(S1): 505-511 . 百度学术
    6. 聂生东,魏传令,张小兵. 基于磁共振影像的乳腺癌计算机辅助检测方法综述. 上海理工大学学报. 2022(01): 1-10+41 . 百度学术
    7. 罗圣钦,陈金怡,李洪均. 基于注意力机制的多尺度残差UNet实现乳腺癌灶分割. 计算机应用. 2022(03): 818-824 . 百度学术
    8. 徐胜军,周盈希,孟月波,刘光辉,史亚. 基于多节点拓扑重叠测度高阶MRF模型的图像分割. 自动化学报. 2022(05): 1353-1369 . 本站查看
    9. 王小鹏,魏统艺,房超,朱生阳. 自适应非局部空间约束与K-L信息的模糊C-均值噪声图像分割算法. 控制理论与应用. 2022(07): 1261-1271 . 百度学术
    10. 王超,王永顺,狄凡. 快速自动模糊C-均值聚类彩色图像分割算法. 激光与光电子学进展. 2022(22): 82-88 . 百度学术
    11. 邱成健,刘青山,宋余庆,刘哲. 基于循环显著性校准网络的胰腺分割方法. 自动化学报. 2022(11): 2703-2717 . 本站查看
    12. 秦传波,宋子玉,曾军英,田联房,李芳. 联合多尺度和注意力-残差的深度监督乳腺癌分割. 光学精密工程. 2021(04): 877-895 . 百度学术
    13. 胡高珍,徐胜军,孟月波,刘光辉,冯峰,段中兴. 基于边缘约束局部区域MRF的图像分割方法. 计算机工程. 2021(06): 253-261+270 . 百度学术
    14. 刘壮盛,李昌林,衣利磊,李智,陈业航,李荣岗,罗学毛,龙晚生,冯宝. MRI纹理特征分析预测浸润性乳腺癌脉管浸润. 中国医学影像技术. 2020(11): 1637-1642 . 百度学术

    其他类型引用(7)

  • 加载中
图(7) / 表(3)
计量
  • 文章访问数:  1553
  • HTML全文浏览量:  315
  • PDF下载量:  186
  • 被引次数: 21
出版历程
  • 收稿日期:  2018-11-14
  • 录用日期:  2019-06-21
  • 网络出版日期:  2020-07-10
  • 刊出日期:  2020-07-10

目录

/

返回文章
返回