一种乳房内增强背景组织的提取方法和装置的制造方法

一种乳房内增强背景组织的提取方法和装置的制造方法
【技术领域】
[0001] 本发明涉及图像处理领域，尤其涉及一种乳房内增强背景组织的提取方法和装 置。
【背景技术】
[0002] 乳腺癌是中老年女性最常见的恶性肿瘤，发病率和死亡率正逐年上升，而对乳腺 癌的风险预测与早期诊断可W有效地提高乳腺癌的治愈率。在临床中，乳房动态对比度增 强（dynamiccontrastenhanced)磁共振技术已经被广泛地用于风险预测和乳腺癌早期诊 断中，尤其是对打入造影剂之前和之后的磁共振图像的增强情况进行比较观察。
[0003] 在临床中，乳房动态对比度增强磁共振图像通常包括若干个磁共振图像扫描序 列，W满足放射科医师的不同需求。为了便于比较观察磁共振图像中乳房内组织打入造影 剂之前和之后的增强情况，该磁共振图像扫描序列通常采用T1加权的压脂图像序列。T1加 权的压脂图像序列较其它类型图像序列具有较高的噪声水平、较强的伪影干扰和较低的图 像灰度对比度，而计算机辅助诊断方案有助于提高放射科医师对H维乳房动态对比度增强 磁共振图像阅片的准确性和效率。因此，如何自动提取乳房内增强背景组织化ackground parenchymalenhancement)是一项非常重要的工作

【发明内容】

[0004] 本发明解决的问题是提供一种乳房内增强背景组织的提取方法和装置，不仅准确 和有效地实现打入造影剂前后所有时间点的磁共振图像序列中乳房组织的全自动分割，还 能得到准确的增强背景组织。
[0005] 为了解决上述问题，本发明提供了一种乳房内增强背景组织的提取方法，包括如 下步骤：
[000引 （1)输入图像A和图像B，对所述图像A和图像B进行配准，得到配准图像序列，并 对所述配准图像序列进行融合，得到融合图像序列；
[0007] (2)基于所述融合图像序列，获取乳房的皮肤线和胸壁线，并根据所述皮肤线和胸 壁线，得到所述乳房组织的分割图像；
[0008] (3)根据所述分割图像和所述配准图像序列，得到所述乳房区域中增强的背景组 织。
[0009] 上述所述一种乳房内增强背景组织的提取方法，其中，所述图像A和图像B分别为 打入造影剂前和打入造影剂后的乳房磁共振序列图像。
[0010] 上述所述一种乳房内增强背景组织的提取方法，其中，通过对所述配准图像序列 加权的方法，得到所述融合图像序列。
[0011] 上述所述一种乳房内增强背景组织的提取方法，其中，得到所述乳房的分割图像 的过程为：基于所述融合图像序列，对所述皮肤线和胸壁线进行线性连接，形成封闭区域， 再对所述封闭区域进行填充。
[0012] 上述所述一种乳房内增强背景组织的提取方法，其中，得到所述乳房区域中增强 的背景组织的过程为：
[0013] (1)将所述图像A和图像B配准后得到的配准图像序列进行相减，得到减影图像；
[0014] (2)根据所述分割图像，在所述减影图像中定位乳房组织区域；
[0015] (3)根据预定阔值，对所述减影图像进行分割，得到所述乳房区域中增强的背景组 织。
[0016] 上述所述一种乳房内增强背景组织的提取方法，其中，通过理论模型得到所述预 定阔值。
[0017] 上述所述一种乳房内增强背景组织的提取方法，其中，所述预定阔值为3^/I<7，其 中，所述0为所述图像A或图像B的噪声水平。
[0018] 为了解决上述问题，本发明还提供了一种乳房内增强背景组织的提取装置，包 括：
[0019] 融合图像序列获取单元，适用于对输入图像A和图像B进行配准，W得到配准序列 图像，并对所述配准序列图像进行融合，W得到融合图像序列；
[0020] 分割图像获取单元，适用于基于所述融合图像序列，获取乳房的皮肤线和胸壁线， 并根据所述皮肤线和胸壁线，得到所述乳房组织的分割图像；
[0021] 增强背景组织获取单元，适用于根据所述分割图像和所述配准图像序列，得到所 述乳房区域中增强的背景组织。
[0022] 上述所述一种乳房内增强背景组织的提取装置，其中，所述融合图像序列获取单 元包括加权单元，所述加权单元适于对所述配准图像序列进行加权，得到所述融合图像序 列。
[0023] 上述所述一种乳房内增强背景组织的提取装置，其中，所述增强背景组织获取单 元包括相减单元、定位单元和获取单元，其中，所述相减单元适于将所述图像A和图像B配 准后得到的配准图像序列进行相减，得到减影图像；所述定位单元适于根据所述分割图像， 在所述减影图像中定位乳房区域；所述获取单元适于根据预定阔值，对所述减影图像进行 分割，得到所述乳房区域中增强的背景组织。
[0024] 与现有技术相比，本发明通过配准方法融合了打入造影剂前后所有时间点的磁共 振图像序列的信息，提高图像的信号噪声比，并准确和有效地实现打入造影剂前后所有时 间点的磁共振图像序列中乳房组织的全自动分割；
[00巧]进一步地，通过理论模型确定阔值，得到乳房区域中增强的背景组织，并基于所述 增强背景组织，对受试者的乳腺癌风险进行量化分析。
【附图说明】
[0026] 图1所示为本发明实施例一种乳房内增强背景组织的提取方法的流程示意图；
[0027] 图2所示为本发明实施例得到所述乳房区域中增强的背景组织的流程示意图； [002引图3所示为本发明实施例一种乳房内增强背景组织的提取装置的结构示意图。
【具体实施方式】
[0029] 在下面的描述中阐述了很多具体细节W便于充分理解本发明。但是本发明能够W 很多不同于在此描述的其它方式来实施，本领域技术人员可w在不违背本发明内涵的情况 下做类似推广，因此本发明不受下面公开的具体实施的限制。
[0030] 其次，本发明利用示意图进行详细描述，在详述本发明实施例时，为便于说明，所 述示意图只是实例，其在此不应限制本发明保护的范围。
[0031] 下面结合附图和实施例对本发明一种乳房内增强背景组织的提取方法和装置进 行详细地说明。本发明实施例乳房内增强背景组织提取方法如图1所示，首先，执行步骤 S1，输入图像A和图像B，并对所述图像A和图像B进行配准，得到配准图像序列，并对所述 配准图像序列进行融合，得到融合图像序列。其中，所述图像A和图像B分别为打入造影剂 前和打入造影剂后的乳房磁共振序列图像。临床中，放射科医师通过比较观察打入造影剂 之前和造影剂之后磁共振图像序列的增强情况进行诊断。由于病患在磁共振图像拍摄过 程中可能存在轻微运动，导致不同时间点的磁共振图像可能存在位置上的偏移，因此，首先 需要对打入造影剂前后所有时间点的磁共振图像序列进行配准。具体地，在本实施例中， 根据H维乳房组织的特点，采用一种已有的非刚体配准方法实现打入造影剂之后每个时间 点的磁共振图像序列与打入造影剂之前的磁共振图像序列的配准，所述非刚体配准方法的 具体过程详参化uang,X.H.，Arridge,S.，Hawkes,D.J.，andOurselin,S.，''Anonrigid registrationframeworkusingspatiallyencodedmutualinformationandfree-form deformations, "I邸ETransactionsonMedicalImaging30 (10) ,1819-1828 (2011)。通过 上述方法，得到打入造影剂之前和之后所有时间点的配准图像序列。
[0032] 基于上述得到的配准图像序列，将所述打入造影剂之前和之后所有时间点的磁共 振图像序列进行融合，得到融合图像序列Ifusi。。。在本实施例中，通过对所述配准图像序列 进行加权的方法，得到所述融合图像序列Ifusi。。，进行融合的公式如下：
[0033]
[0034 ] 其中，Ii为配准后打入造影剂之前任一时间点的磁共振图像序列，12. ..I。为配准 后打入造影剂之后的每个时间点的磁共振图像序列，wei曲ti为图像进行融合时所占的权 重。具体地，本实施例中，将配准后的打入造影剂之前图像序列中的每幅图像都和配准后的 打入造影剂之后每个时间点的磁共振图像序列中的每幅图像进行加权，wei曲ti= 1，即融 合图像序列Ifud。。为配准后打入造影剂前后所有时间点的磁共振图像序列的平均图像。
[0035] 接着，执行步骤S2,基于所述融合图像序列，获取乳房的皮肤线和胸壁线，并根据 所述皮肤线和胸壁线，得到所述乳房组织的分割图像。具体地，首先建立人体坐标系：左右 方向的冠状轴为X轴，前后方向的矢状轴为Y轴，上下方向的垂直轴为Z轴。在本实施例 中，对于所述步骤S1中得到的融合图像序列，根据所述乳房组织在图像中的特征点和特征 位置，确定初始参数。其中，所述初始参数包括融合图像序列的初始层和终止层，W及所述 融合图像序列在对应极坐标系转换中的射线扫描起始点P、射线扫描初始半径R、射线扫描 起始角度aW及射线扫描终止角度目和分段位置。具体地，将所述融合图像序列沿Z轴 做平均值投影，得到平均值投影图。在所述平均值投影图上，计算所述平均值投影图中每一 行所有像素点的灰度值之和，选取所述灰度值之和最大的行作为下边界。该下边界的位置 也为最大值投影图的下边界lb的位置，其中，最大值投影图为所述融合图像序列沿Z轴做 最大值投影得到。在最大值投影图上，从左到右扫描下边界lb位置上的所有像素点，第一个 灰度值不为0的像素点所在的位置为左边界ii，最后一个灰度值不为0的像素点所在的位 置为右边界1,。在最大值投影图上，从上到下进行扫描，检测到y值最小的非零像素点，经 过该像素点且平行于最大值投影图X轴的直线为上边界It。在最大值投影图上，通过canny 边缘检测得到边缘线，左边界li和右边界If的中也线将所述边缘线分为两段，分别为左边 的边缘线和右边的边缘线，再分别在两段边缘线上获取特征点Si和S2。
[0036] 在最大值投影图上，射线扫描起始点的坐标位置为（(X1+X2) /2,0)，其中，xi和X,分 别为两个特征点Si和S,在最大值投影图中X轴的坐标值；射线扫描起始点P与两个特征点 连线中的较长连线的长度为射线扫描的初始半径R;左边界li与下边界lb的交点Qi和射线 扫描起始点P的连线确定了射线扫描的起始角度a，右边界If与下边界Ib的交点Q2和射 线扫描起始点P的连线确定了射线扫描的终止角度目；两个特征点Si和S,和射线扫描起 始点P的连线确定了曲线跳跃处理的分段位置。其中，最大值投影图为融合图像序列沿Z 轴做最大值投影得到。接着，根据所述射线扫描起始点、射线扫描初始半径和射线扫描的起 始角度和终止角度，基于极坐标转换和射线扫描方法将所述中间层图像转换为二维转换图 像；通过动态规划的方法获取所述二维转换图像中的目标边界线，所述目标边界线即为所 述中间层的胸壁线经过转换后在二维转换图像中的位置；判断所述目标边界线是否具有波 动，若具有波动，则通过分段处理的方法消除所述波动；基于所述目标边界线，通过极坐标 逆转换的方法获取所述中间层的胸壁线。接着，在所述除了中间层W外的层中，所述中间层 两侧分别按离中间层由近及远的次序依次选择图像作为当前层，获取所述当前层的二维转 换图像和图像梯度；基于高斯核函数和所述当前层前一层或后一层中胸壁线的检测结果， 对所述当前层图像梯度进行改进；通过动态规划的方法，得到目标边界线，所述目标边界线 为所述当前层的胸壁线在所述二维转换图像中的位置；判断所述当前层的胸壁线是否具有 波动，若具有波动，则通过所述分段处理的方法消除所述波动；将所述目标边界通过极坐标 逆转换的方法转换为所述当前层图像中的胸壁线。
[0037] 同样地，对于所述步骤S1中得到的融合图像序列，根据所述乳房组织在图像中的 特征点和特征位置，确定左、右乳房的感兴趣区域。接着，根据初始皮肤线灰度的统计分布， 增强所述融合图像序列的乳房皮肤线。接着，根据射线扫描起始点、射线扫描半径和扫描射 线数目，基于极坐标转换和射线扫描方法分别将所述左、右乳房感兴趣区域中间层图像转 换为二维转换图像；其中，所述射线扫描起始点为所述左、右乳房的中也点在所述中间层图 像中所述左、右乳房感兴趣区域下边界上的投影点；所述射线扫描起始点距离所述左、右乳 房乳头的距离确定所述射线扫描半径；通过动态规划的方法获取所述二维转换图像中的目 标边界线，所述目标边界线即为所述中间层图像上的皮肤线经过转换后在二维转换图像中 的位置；判断所述目标边界线是否具有波动，若具有波动，则通过第二级别饱和化处理消除 所述波动；基于所述目标边界线，通过极坐标逆转换的方法获取所述中间层图像上的皮肤 线。接着，在所述除中间层W外的层中，所述中间层两侧分别按离中间层由近及远的次序 依次选择图像层作为当前层，获取所述当前层的二维转换图像和图像梯度；基于高斯核函 数和所述当前层前一层或后一层图像上皮肤线的检测结果，对所述当前层图像梯度进行改 进；通过动态规划的方法，得到目标边界线，所述目标边界线为所述当前层图像上的皮肤线 在所述二维转换图像中的位置；判断所述当前层的皮肤线是否具有波动，若具有波动，则通 过第二级别饱和化处理消除所述波动；将所述目标边界通过极坐标逆转换的方法转换为所 述当前层图像上的皮肤线。
[0038] 通过上述方法，基于融合后的图像序列，获取乳房的皮肤线和胸壁线的检测结果。 通过检测到的乳房皮肤线和胸壁线的有效结合，实现打入造影剂前后所有时间点的磁共振 图像中乳房组织的全自动分割。由于打入造影剂前后所有时间点的磁共振图像序列经过了 配准，所W基于融合的图像序列的乳房分割结果适合于打入造影剂前后所有时间点的磁共 振图像序列，而不需要再对每个时间点的磁共振图像序列上重复执行乳房分割过程。具体 地，在本实施例中，对检测到的乳房皮肤线和胸壁线的相应端点先进行线性连接，形成封闭 区域，再对封闭区域进行填充处理，即完成对乳房组织的分割，得到分割图像（二值图像）。 在融合图像序列中实现对乳房组织的分割，由于经过图像配准融合后获得的图像序列的图 像信号噪声比高于每个时间点的原始图像序列的图像信号噪声比，因此，可W有效地消除 噪声和伪影对乳房分割的影响。
[0039] 接着，执行步骤S3,根据所述分割图像和所述配准图像序列，得到所述乳房区域中 增强的背景组织。乳房内背景组织的增强程度是乳腺癌风险的一个重要量化指标，基于已 获得的乳房分割结果和已配准的打入造影剂前后两个磁共振序列图像，可W获得乳房内增 强的背景组织。具体地，得到所述乳房区域中增强的背景组织的过程如图2所示，首先，执 行步骤S201，将所述图像A和图像B配准后得到的配准图像序列进行相减，得到减影图像。 具体地，将步骤S1中得到的打入造影剂之前和之后配准图像序列进行相减，得到减影图像 序列。接着，执行步骤S202,根据所述分割图像，在所述减影图像中定位乳房区域。具体 地，在本实施例中，根据步骤S2中得到的乳房组织分割结果，根据分割图像中乳房组织的 位置，在所述减影图像中定位出与所述位置所对应的乳房区域。
[0040] 接着，执行步骤S203,根据预定阔值，对所述减影图像进行分割，得到所述乳房区 域中增强的背景组织。由于原始磁共振图像序列中存在噪声，导致减影图像序列中也存在 噪声，因此，需要设定阔值去除减影图像序列中的噪声部分，实现基于乳房组织内背景组织 增强的全自动乳腺癌风险量化分析。具体地，在本实施例中，通过理论模型得到所述预定阔 值。通过理论模型推导预定阔值的具体过程为：
[0041] 定义打入造影剂之前时间点的图像为Ib。,。,。，打入造影剂之后某个时间点的图像 为lafter，Ibetae和lafter分别可W表示为Ibetae=Ibetae'+n〇ise和lafter=lafter'+n〇ise, 其中，Ibefore'和Irfter'为图像的真实像素值，noise为图像的噪声， 假设图像Ibeftre和Lfter中的noise来源相同且均近似符合正态分布，即noise~N(y，〇2)，所述0为所述图像A 或图像B的噪声水平。则打入造影剂前后两个磁共振图像序列的减影图像序列Lub可W表 不为：
[0042] 1日血=lafter-Ib址。re=lafter'-Ib址。re' +n〇ise曰血
[004引其中，noise,ub~N(0,202)。在本实施例中，所述预定阔值为3^/1^，该阔值能 压制减影图像中99. 6%的噪声。需要说明的是，所述预定阔值也可W采用与图像A或 图像B噪声水平相关的其他阔值。对于噪声水平0,通过已有的噪声估计方法获得，详 参Coupe，P.，Yger，P.，Prima，S.，Hellier，P.，Kervrann，C.，andBarillot，C.，"An optimizedblockwisenonlocalmeansdenoisingfilterfor3-DMagneticResonance Images, "I邸EhansactionsonMedicalImaging27(4)，425-441 (2008)。需要说明的是， 所述噪声水平o也可W通过计算背景平坦区域中像素点灰度值的标准差得到。
[0044] 通过上述方法对定位出乳房区域的减影图像的噪声进行了抑制，得到所述乳房组 织区域中增强的背景组织。需要说明的是，根据所述增强背景组织的特征，可W对受试者进 行乳腺癌风险的量化分析。由于通过上述步骤阔值的处理，减影图像大部分噪声被压制，可 W计算乳房内增强背景组织的相关特征，从而进行乳腺癌风险的量化分析。比如，通过计算 增强的背景组织体积与乳房组织体积的比例对所述受试者进行乳腺癌风险的量化分析，从 而医生可W根据该结果对受试者的乳腺癌风险进行进一步地分析。
[0045] 对应于上述乳房内增强背景组织的提取方法，本发明实施例还提供了一种乳房内 增强背景组织的提取装置，如图3所示，包括融合图像序列获取单元1、分割图像获取单元2 和增强背景组织获取单元3。
[0046] 所述融合图像序列获取单元1适于对输入图像A和图像B进行配准，W得到配准 序列图像，并对所述配准序列图像进行融合，W得到融合图像序列。其中，所述图像A和图 像B分别为打入造影剂前和打入造影剂后的乳房磁共振序列图像；所述融合图像序列获取 单元包括加权单元，所述加权单元适于对所述配准图像序列进行加权，得到所述融合图像 序列。
[0047] 所述分割图像获取单元2适于基于所述融合图像序列，获取乳房的皮肤线和胸壁 线，并根据所述皮肤线和胸壁线，得到所述乳房组织的分割图像。具体地，得到所述乳房的 分割图像的过程为：对所述皮肤线和胸壁线进行线性连接，形成封闭区域，再对所述封闭区 域进行填充。
[0048] 所述增强背景组织获取单元3适于根据所述分割图像和所述配准图像序列，得到 所述乳房区域中增强的背景组织。所述获取增强背景组织单元3包括相减单元、定位单元 和获取单元。所述相减单元适于将所述图像A和图像B配准后得到的配准图像序列进行 相减，得到减影图像；所述定位单元适于根据所述分割图像，在所述减影图像中定位乳房区 域；所述获取单元适于根据预定阔值，对所述减影图像进行分割，得到所述乳房区域中增强 的背景组织。通过计算增强的背景组织体积与乳房体积的比例对所述乳房进行乳腺癌风险 的量化分析，医生可W根据该结果对受试者的乳腺癌风险进行进一步地分析。
[0049] 上述乳房内增强背景组织提取装置中各单元的配合及工作过程可W参考上述乳 腺癌风险量化分析方法的说明，在此不再费述。
[0050] 本发明虽然已W较佳实施例公开如上，但其并不是用来限定本发明，任何本领域 技术人员在不脱离本发明的精神和范围内，都可W利用上述掲示的方法和技术内容对本发 明技术方案做出可能的变动和修改，因此，凡是未脱离本发明技术方案的内容，依据本发明 的技术实质对W上实施例所作的任何简单修改、等同变化及修饰，均属于本发明技术方案 的保护范围。
【主权项】
1. 一种乳房内增强背景组织的提取方法，其特征在于，包括如下步骤： (1) 输入图像A和图像B，对所述图像A和图像B进行配准，得到配准图像序列，并对所 述配准图像序列进行融合，得到融合图像序列； (2) 基于所述融合图像序列，获取乳房的皮肤线和胸壁线，并根据所述皮肤线和胸壁 线，得到所述乳房组织的分割图像； (3) 根据所述分割图像和所述配准图像序列，得到所述乳房区域中增强的背景组织。2. 如权利要求1所述一种乳房内增强背景组织的提取方法，其特征在于，所述图像A和 图像B分别为打入造影剂前和打入造影剂后的乳房磁共振序列图像。3. 如权利要求1所述一种乳房内增强背景组织的提取方法，其特征在于，通过对所述 配准图像序列加权的方法，得到所述融合图像序列。4. 如权利要求1所述一种乳房内增强背景组织的提取方法，其特征在于，得到所述乳 房的分割图像的过程为：基于所述融合图像序列，对所述皮肤线和胸壁线进行线性连接，形 成封闭区域，再对所述封闭区域进行填充。5. 如权利要求1所述一种乳房内增强背景组织的提取方法，其特征在于，得到所述乳 房区域中增强的背景组织的过程为： (1) 将所述图像A和图像B配准后得到的配准图像序列进行相减，得到减影图像； (2) 根据所述分割图像，在所述减影图像中定位乳房组织区域； (3) 根据预定阈值，对所述减影图像进行分割，得到所述乳房区域中增强的背景组织。6. 如权利要求5所述一种乳房内增强背景组织的提取方法，其特征在于，通过理论模 型得到所述预定阈值。7. 如权利要求6所述一种乳房内增强背景组织的提取方法，其特征在于，所述预定阈 值为3V^t，其中，所述。为所述图像A或图像B的噪声水平。8. -种乳房内增强背景组织的提取装置，其特征在于，包括： 融合图像序列获取单元，适用于对输入图像A和图像B进行配准，以得到配准序列图 像，并对所述配准序列图像进行融合，以得到融合图像序列； 分割图像获取单元，适用于基于所述融合图像序列，获取乳房的皮肤线和胸壁线，并根 据所述皮肤线和胸壁线，得到所述乳房组织的分割图像； 增强背景组织获取单元，适用于根据所述分割图像和所述配准图像序列，得到所述乳 房区域中增强的背景组织。9. 如权利要求8所述一种乳房内增强背景组织的提取装置，其特征在于，所述融合图 像序列获取单元包括加权单元，所述加权单元适于对所述配准图像序列进行加权，得到所 述融合图像序列。10. 如权利要求8所述的乳房内增强背景组织的提取装置，其特征在于：所述增强背景 组织获取单元包括相减单元、定位单元和获取单元，其中，所述相减单元适于将所述图像A 和图像B配准后得到的配准图像序列进行相减，得到减影图像；所述定位单元适于根据所 述分割图像，在所述减影图像中定位乳房区域；所述获取单元适于根据预定阈值，对所述减 影图像进行分割，得到所述乳房区域中增强的背景组织。
【专利摘要】本发明提供了一种乳房内增强背景组织的提取方法和装置，所述方法包括如下步骤：(1)输入图像A和图像B，对所述图像A和图像B进行配准，得到配准图像序列，并对所述配准图像序列进行融合，得到融合图像序列；(2)基于所述融合图像序列，获取乳房的皮肤线和胸壁线，并根据所述皮肤线和胸壁线，得到所述乳房组织的分割图像；(3)根据所述分割图像和所述配准图像序列，得到所述乳房区域中增强的背景组织。本技术方案不仅准确和有效地实现打入造影剂前后所有时间点的磁共振图像序列中乳房组织的全自动分割，还能得到准确的增强背景组织。
【IPC分类】G06T5/50, A61B5/055, G06T7/00
【公开号】CN104881858
【申请号】CN201410070338
【发明人】姜娈, 李强
【申请人】上海联影医疗科技有限公司
【公开日】2015年9月2日
【申请日】2014年2月27日