旋转式乳房影像的病变检测方法及病变检测装置与流程

本发明涉及一种医疗图像处理技术，尤其涉及一种旋转式乳房影像的病变(lesion)检测方法及病变检测装置。

背景技术：

乳腺癌(mammarycarcinoma)是女性常见的恶性肿瘤之一，其主要症状包括乳房肿瘤(tumor)、异常分泌物或形状变异等。提早筛检出乳房的异常症状，将有助于尽早针对肿瘤进行治疗，以降低癌细胞恶化或扩散等问题。诸如临床或自我乳房检测、活体组织检查、乳房摄影术(mammography)、超音波(ultrasound)或磁共振(magneticresonance)显像等筛检方式已广泛在临床上使用或成为学术研究的重要议题。

据研究指出，与低密度乳房相比，拥有高密度乳房的女性具有相对高的风险罹患乳癌。因此，乳房及乳腺组织的密度分析亦是乳癌评估的重要因素之一。另一方面，虽然现今临床上已使用计算机辅助检测(computeraideddetection；cade)系统来自动辨识肿瘤、肿块或钙化点，但是都存在高伪阳性的风险。

技术实现要素：

本发明提供一种旋转式乳房影像的病变检测方法及病变检测装置，其可有效降低计算机辅助检测系统的伪阳性。

本发明提出一种旋转式乳房影像的病变检测方法，且至少包括(但不仅限于)下列步骤。取得一组旋转式乳房影像。此组旋转式乳房影像包括多张子乳房影像。对这些子乳房影像进行重组，以产生重组乳房影像。比对重组乳房影像以及未重组的组旋转式乳房影像以确认至少一个病变位置。

另一观点而言，本发明提出一种病变检测装置，其至少包括(但不仅限于)储存单元及处理单元。而处理单元耦接储存单元，且取得一组旋转式乳 房影像。而此组旋转式乳房影像包括多张子乳房影像。处理单元对这些子乳房影像进行重组，以产生重组乳房影像。处理单元比对重组乳房影像以及未重组的旋转式乳房影像以确认至少一个病变位置。

基于上述，本发明实施例所提出的旋转式乳房影像的病变检测方法及病变检测装置，其将旋转式乳房影像重组，并将重组乳房影像及旋转式乳房影像进行比对，以确认病变(例如，肿瘤、肿块或钙化点等)位置。据此，本发明实施例便能助于减少计算机辅助检测系统的伪阳性。

为让本发明的上述特征和优点能更明显易懂，下文特举实施例，并配合附图作详细说明如下。

附图说明

图1是依据本发明一实施例说明病变检测装置的方框图；

图2是依据本发明一实施例说明一种旋转式乳房影像的病变检测方法流程图；

图3是旋转式扫描的示意图范例；

图4为旋转式乳房影像于立体空间的示意图；

图5a～图5c是依据本发明一实施例说明填补作业的示意图；

图6是依据本发明一实施例说明自动化病变检测方法的流程图；

图7a是未重组的旋转式乳房影像的局部影像；

图7b是区域分割影像；

图7c所示为经图7b的区域分割影像所判断出的可疑病变区域；

图7d所示为图7c经伪阳性减少作业筛选的病变位置；

图8是一范例说明决定厚度区域的示意图；

图9a是一范例说明投影影像；

图9b是一范例说明分割后的投影影像。

附图标记：

100：病变检测装置

110：储存单元

111：影像输入模块

113：影像重组模块

115：病变判断模块

117：影像质量模块

150：处理单元

s210～s250、s610～s670：步骤

301：乳房

320：乳房取像容器

340：扫描仪

410、810：圆柱体

431、432：子乳房影像

501：区隔位置

503：缺失位置

520：矩形

530：影像集合

710：旋转式乳房影像

711：肿瘤

730：区域分割影像

731：可疑病变区域

733：病变位置

811：厚度区域

910、920：投影影像

911：目标区块

921：皮肤组织类别区块

923：拍摄失误类别区块

r：半径

θ：旋转角度

th：厚度值

h：子乳房影像的高

w：子乳房影像的宽

rd1：旋转方向

具体实施方式

图1是依据本发明一实施例说明病变检测装置的方框图。请参照图1，病变检测装置100至少包括(但不仅限于)储存单元110及处理单元150。病变检测装置100可以是服务器、客户端、桌面计算机、笔记本电脑、网络计算机、工作站、个人数字助理(personaldigitalassistant；pda)、平板个人计算机(personalcomputer；pc)、计算机辅助检测(cade)系统等电子装置，且不以此为限。

储存单元110可以是任何型态的固定或可移动随机存取内存(randomaccessmemory，ram)、只读存储器(read-onlymemory，rom)、闪存(flashmemory)或类似元件或上述元件的组合。在本实施例中，储存单元110用以储存旋转式乳房影像、子乳房影像、重组乳房影像、扫描参数、程序代码、装置组态、缓冲的或永久的数据，并记录影像输入模块111、影像重组模块113、病变判断模块115及影像质量模块117等软件程序。前述模块的详细运作内容待稍后实施例详细说明。本实施例中所述的储存单元110并未限制是单一内存元件，上述的各软件模块亦可以分开储存在两个或两个以上相同或不同型态的内存元件中。

处理单元150的功能可藉由使用诸如中央处理单元(centralprocessingunit；cpu)、微处理器、微控制器、数字信号处理(digitalsignalprocessing；dsp)芯片、场可程序化逻辑门阵列(fieldprogrammablegatearray；fpga)等可程序化单元来实施。处理单元150的功能亦可用独立电子装置或集成电路(integratedcircuit；ic)实施，且处理单元150亦可用硬件或软件实施，存取储存单元110中的各模块、软件程序、程序代码等并执行对应操作。

为了方便理解本发明实施例的操作流程，以下将举诸多实施例详细说明本发明实施例中病变检测装置100进行乳房图像处理及病变检测的流程。图2是依据本发明一实施例说明一种旋转式乳房影像的病变检测方法流程图。请参照图2，本实施例的方法适用于图1中的病变检测装置100。下文中，将搭配病变检测装置100中的各项元件及模块说明本发明实施例所述的方法。本方法的各个流程可依照实施情形而随之调整，且并不仅限于此。

在步骤s210中，处理单元150通过影像输入模块111取得一组旋转式乳房影像。而此组旋转式乳房影像包括多张(或片段(slice))子乳房影像。 在本发明实施例中，这些子乳房影像是通过扫描仪绕着乳房取像容器下方进行旋转扫描一圈的过程中所分别取得。此扫描仪例如是具有基于自动乳房超音波(automatedbreastultrasound；abus)、数字乳房断层层析(digitalbreasttomosynthesis；dbt)、磁共振显影(magneticresonanceimaging：mri)等医疗影像扫描技术的探头(probe)。针对超音波扫描，此乳房取像容器中可装载液体或水溶性软膏等来作为超音波传导媒介。

举例而言，图3是旋转式扫描的示意图范例。请参照图3，使用者以俯卧方式，使其乳房301可完全或部分地置于圆柱型的乳房取像容器320中。具有探头的可移动式扫描仪340置于(固定式)乳房取像容器320的圆柱型底部下方，且可通过机械装置(未显示)来驱动其以旋转方向rd1(即，顺时针)或旋转方向rd1的相反方向(即，逆时针)旋转至少一圈(360°)，以使扫描范围能涵盖乳房301的全部或部分投影面积。在扫描的过程中，扫描仪340会随着每旋转一旋转角度(例如，3°、5°、8°等)对乳房301扫描以取得一张子乳房影像。例如，旋转角度为3°，则扫描仪340旋转一圈(360°)后便可取得120张子乳房影像。

需说明的是，在进行旋转扫描之前，可默认或手动调整扫描参数。此扫描参数至少包括(但不仅限于)影像扫描起始位置、旋转方向(顺时针或是逆时针)、旋转角度等，且可记录于储存单元110中以备后续使用。

影像输入模块111可自储存单元110、通过无线或有线通信单元(例如，wi-fi、以太网(ethernet))、直接通过图3所述医学影像扫描仪(例如，abus扫描仪器、mri扫描仪器等)或储存装置(例如，dvd、随身碟、硬盘等)取得旋转式乳房影像。

在步骤s230中，处理单元150通过影像重组模块113对这些子乳房影像进行重组，以产生重组乳房影像。为了让通过旋转方式取得的旋转乳房影像能够依照不同方向的视角(view)来观看，本发明实施例依据旋转式扫描的旋转特性来进行重组。

在一实施例中，影像重组模块113将依据扫描仪的扫描起始位置、旋转角度及旋转方向将这些子乳房影像转换成影像集合，判断各相邻的子乳房影像之间的间隔位于此影像集合中的缺失位置，且通过内插法补足此缺失位置。

具体而言，请先参照图4所示为旋转式乳房影像于立体空间的示意图。 为了方便说明，图4中仅以两张子乳房影像431、子乳房影像432进行说明，然不以此为限。每张子乳房影像的宽w即为旋转半径(例如，图3中扫描仪340单次可扫描的长度)，而其宽h为扫描仪可扫描的最大深度。影像重组模块113将旋转扫描一圈的每张子乳房影像(例如，子乳房影像431、子乳房影像432)依据先前记录或默认的扫描参数(例如，影像扫描起始位置、顺时针或是逆时针旋转、旋转角度等)，以各自对应的扫描位置依序排列，以形成圆柱体410。此时，圆柱体410中例如是两张子乳房影像431、子乳房影像432之间存在间隔(即，未存在像素或扫描影像)。

接着，影像重组模块113对这些间隔进行填补作业。图5a～图5c是依据本发明一实施例说明填补作业的示意图。请先参照图5a，影像重组模块113将图4中数张子乳房影像所形成的圆柱体410建立二维直角坐标系统。例如，圆柱体410的圆形顶所属平面的二维直角坐标系统。

影像重组模块113可将此二维直角坐标系统转换成以旋转角度及半径长度的坐标系统。具体而言，首先，定义相邻两子乳房影像中区隔位置501在直角坐标系统的坐标为(x1,y2)(以圆柱体410中圆形顶的中心点为原点(0,0)且影像扫描起始位置至圆形顶的中心点的沿线为x轴(或称水平轴))、圆柱体410中圆形顶的直径宽width(即，2*w)及高height(即，h)。接着，影像重组模块113将坐标(x1,y1)以下列公式(1)(勾股定理)及(2)进行转换，以取得坐标(r1,θ1)：

r1代表此位置相距圆柱体410中圆形顶边缘的距离，而θ1代表自扫描起始位置所旋转的角度。

藉由前述转换，影像重组模块113将部分或全部的子乳房影像映射成影像集合。请参照图5b，此坐标系统的水平轴为(自扫描起始位置起始所旋转的)旋转角度θ，而垂直轴为半径r。多张子乳房影像依旋转方向在此r-θ坐标系统上依序排列，并重组成影像集合530。例如，图4中旋转方向rd1为顺时针，则在影像集合中子乳房影像431将位于子乳房影像432左侧。而图 5a中的区隔位置501亦可映射至影像集合530中缺失位置503。

请接着参照图5c为影像集合530的局部放大示意图。影像重组模块113可通过诸如双线性内插法(bilinearinterpolation)、双立方内插法(bicubicinterpolation)、最邻近点(nearest)内插法等来补足此缺失位置503。而在补足所有相邻子乳房影像间的间隔之后，便可形成完整的(例如，无缺失位置)影像集合。影像重组模块113可以此影像集合再将自r-θ坐标系统转换回如图5a的呈现方式，以形成三维的重组乳房影像。而此重组乳房影像并不会存在如图4所示相邻两子乳房影像之间的间隔。

而处理单元150可进一步通过显示单元(未显示，例如，液晶显示器(lcd)、电浆显示器面板(pdp)、有机发光二极管(oled)等)显示此重组乳房影像，并可通过输入单元(未显示，例如，触控装置、键盘、鼠标等)接收用户的输入操作来以不同方向(例如，图5a中矩形520是用以显示的方向)的视角检视重组乳房影像。

在步骤s250中，处理单元150通过病变判断模块115比对重组乳房影像以及未重组的旋转式乳房影像以确认至少一个病变位置。在本实施例中，使用基于区域(regionbased)的病变(例如，肿瘤、肿块或钙化点等)检测方法，对旋转式乳房影像进行自动检测。通过所提出的区域筛选方式，可对可疑病变区域进行条件筛选以找出病变位置。

图6是依据本发明一实施例说明自动化病变检测方法的流程图。请参照图6，在步骤s610中，病变判断模块115对未重组的旋转式乳房影像进行区域分割，以产生区域分割影像。病变判断模块115可利用原始影像(即，未重组的旋转式乳房影像中的一片段子乳房影像)进行基于像素或纹路差异的区域分割(例如，分水岭(watershed)切割、马可夫随机场(markovrandomfield；mrf))。

例如，图7a是未重组的旋转式乳房影像的局部影像。假设旋转式乳房影像710中存在肿瘤711。病变判断模块115对旋转式乳房影像710进行区域分割后，形成如图7b所示的区域分割影像730。区域分割影像730中具有相同或相似纹路特征的像素受划分为相同区域(以相同颜色表示)。需说明的是，区域分割算法中的各项参数可依据需求自行调整，本发明实施例不对此限制。

在步骤s630中，处理单元150通过病变判断模块115判断区域分割影像中的至少一个可疑病变区域。具体而言，病变判断模块115将旋转式乳房影像切割后，区域分割影像仍存在许多不必要的区块，且目标病变区块亦包含于其中。病变判断模块115可利用各区块的像素值特征(例如，平均值、最大值、最小值、中位数、变异数等)，依据所欲检测的目标特性(例如，较暗区域、近似椭圆形、长短轴比例等)来进行初步筛选，以决定可疑病变区域，此步骤中所使用的影像各个区块的特征以及特性将可依照需求进行调整。例如，图7c所示为经图7b的区域分割影像730所判断出的可疑病变区域731。

在步骤s650中，病变判断模块115对可疑病变区域进行伪阳性减少作业，以判断病变位置。具体而言，经初步的筛选后，病变判断模块115利用已经筛选过后的可疑病变区块再次对所剩区块进行伪阳性减少(false-positivereduction)。伪阳性减少作业所使用特征可至少包括(但不仅限于)三个部分：形状(例如，面积、长短轴比率)、像素强度(例如，平均值、标准偏差)以及纹理(例如，灰阶共生矩阵(gray-levelco-occurrencematrix；glcm)、马可夫随机场(markovrandomfield；mrf)、贾柏滤波器(gaborfilter)。换句而言，病变判断模块115可基于默认或经由手动选择的特征，对所剩区块进行筛选，以比对出目标病变位置。例如，图7d所示为图7c经伪阳性减少作业筛选的病变位置733。

在步骤s670中，病变判断模块115依据病变位置寻找重组乳房影像中是否存在相连的病变区域，以确认此病变位置。具体而言，病变判断模块115将决定的病变位置(例如，图7d中的病变位置733)与三维的重组乳房影像比对。若在此重组乳房影像中具有相连的病变区域(即，与病变位置相对应的区域)，则代表实际存在此病变(或是存在的机率高于80％、90％等)。而若在此重组乳房影像中未具有相连的病变区域，则代表未存在病变(或是存在的机率低于10％、15％等)。

在一些实施例中，处理单元150可进一步通过显示单元呈现诸如病变位置、可疑病变区域、病变区域、发现病变的提示信息(例如，“发现肿瘤！”)等其中的一个或其组合，以协助医疗人员清楚得知检测状况。

此外，为了维持影像扫描的质量控制(qualitycontrol)，在一些实施例 中，处理单元150通过影像质量模块117可进一步判断旋转式乳房影像是否完整或拍摄失误过高(例如，失误比例大于70％、80％等)。影像质量模块117可对旋转式乳房影像中的部分厚度区域的影像进行垂直投影，以产生投影影像，并依据投影影像中的拍摄失误类别与皮肤组织类别的比值，判断旋转式乳房影像的影像质量。

具体而言，在步骤s210或步骤s230之后或是步骤s250之前，处理单元150可依据影像质量模块117所判断影像质量的结果，来决定是否继续进行后续的病变检测。影像质量模块117可自三维的重组乳房影像或例如图4形成圆柱体410的子乳房影像集合中决定厚度值(例如，2公分、5公分等，不同使用者可能有所不同)的厚度区域。

例如，图8是一范例说明决定厚度区域的示意图。影像质量模块117自形成圆柱体810的子乳房影像集合或是重组乳房影像上方决定厚度值th，从而决定厚度区域811。

影像质量模块117接着对已决定的厚度区域(例如，图8的厚度区域811)的影像进行垂直投影，其可选择各位置中以垂直方向的所有像素中具有最低像素值，来分别作为各位置上的值。

影像质量模块117可进一步去除皮肤组织外围不必要的区域。例如，图9a是一范例说明投影影像。请参照图9a，投影影像910中的目标区块911以外的影像视为不必要的区域。

影像质量模块117可依据目标区域内的像素特征，利用利用区域分割方法(例如，分水岭分割、马可夫随机场等)，将投影影像区分为皮肤组织类别区块与拍摄失误类别区块。例如，图9b是一范例说明分割后的投影影像。请参照图9b，经分割后的投影影像920包括皮肤组织类别区块921与拍摄失误类别区块923。影像质量模块117可利用两个类别的面积来产生比值(例如，拍摄失误类别区块923面积总和/皮肤组织类别区块921面积总和)用于决定此影像扫描质量的优劣程度。例如，若此比值大于质量门槛值(例如，30％、20％等)，则影像质量模块117判断质量不佳，进而重新进行扫描作业。例如，通过显示单元显示“建议重新扫描！”的提示信息。反之，若比值小于质量门槛值(例如，15％、30％等)，则影像质量模块117判断质量良好，且可进行例如步骤s250的病变检测作业。

综上所述，本发明实施例所提出的旋转式乳房影像的病变检测方法及病变检测装置，其是将旋转式乳房影像进一步重组，以方便医疗人员能以不同方向的视角来检视。利用基于区域的病变检测方法能协助筛选病变位置，且可进一步与三维的重组乳房影像进行比对，以降低伪阳性。此外，本发明实施例更能通过影像质量判断，从而维持影像扫描质量。

虽然本发明已以实施例揭示如上，然其并非用以限定本发明，任何所属技术领域中普通技术人员，在不脱离本发明的精神和范围内，当可作些许的改动与润饰，故本发明的保护范围当视所附权利要求界定范围为准。