一种核磁共振肿瘤区域提取方法与流程

1.本发明涉及图像数据处理技术领域，具体涉及一种核磁共振肿瘤区域提取方法。


背景技术：

2.在核磁共振波谱（mrs）中对肿瘤区域数据进行提取的过程中，需要从大量的mrs信号中选择出感兴趣的mrs信号进行分析，在提取过程中首先需要通过mri图像确定多体素波谱采集过程中的感兴趣区域（voi），之后使用核磁共振仪采用化学位移法（csi）进行感兴趣区域的多体素波谱采集。在对感兴趣区域的体素进行磁共振波谱（mrs）提取过程中，需要在感兴趣区域中将肿瘤体素的波谱数据进行提取，从而通过肿瘤区域的波谱数据进行肿瘤提取。在实际的肿瘤区域数据提取过程中，因为体素是三维的空间结构，在提取的过程中需要对感兴趣区域周围的体素波谱进行相似性判断，使得对肿瘤区域的提取更加精确。在判断感兴趣区域体素的相邻区域体素波谱与肿瘤区域的相似性的过程中，首先需要获取voi区域体素的核磁共振波谱，之后根据voi区域中的中心体素向外检测voi区域边缘以外的相邻体素波谱，通过判断波谱的谱线相似性，根据所获取的相似程度将相似的体素加入到voi区域，从而提取到完整的肿瘤区域mrs数据。
3.为了实现上述目的，本领域技术人员通过动态时间规整（dtw）对两个体素的磁共振波谱谱线进行相似性判断。但是在体素波谱的相似性判断过程中因为传统的序列相似性判断是通过固定长度的子区间进行判断，在判断过程中会因为谱线的代谢物峰偏移造成相似性判断错误的情况。


技术实现要素：

4.为了解决现有技术中在体素波谱的相似性判断过程中因为传统的序列相似性判断是通过固定长度的子区间进行判断，在判断过程中会因为谱线的代谢物峰偏移造成相似性判断错误的问题，本发明提供一种核磁共振肿瘤区域提取方法，该方法通过对感兴趣区域边界的相邻体素进行代谢物相似性判断，在感兴趣区域基础上根据波谱相似性进行扩展，使得对感兴趣的肿瘤区域的波谱数据提取更加精确。
5.本发明的目的是提供一种核磁共振肿瘤区域提取方法，包括以下步骤：获取磁共振图像中的感兴趣区域；获取感兴趣区域中的中心体素；将感兴趣区域以外的区域中的体素作为待判断体素；并获取每个体素的相对浓度及波谱数据；利用中心体素和每个待判断体素的波谱数据中所有数据点获取中心体素与每个待判断体素的初始dtw距离；根据中心体素与每个待判断体素的波谱数据中代谢物数据点所处位置对中心体素与每个待判断体素的初始dtw距离进行优化获取中心体素与每个待判断体素的优化后的dtw距离；获取中心体素到每个待判断体素之间所包含的体素的相对浓度形成相对浓度序列，获取相对浓度序列的差分序列；
根据中心体素到每个待判断体素之间所包含的体素数量，中心体素到每个待判断体素之间所包含的体素相对浓度连续梯度变化的体素数量，以及相对浓度序列的差分序列，获取每个待判断体素的相对浓度渐变符合程度；根据中心体素到每个待判断体素之间所包含的体素形成的体素序列、每个待判断体素的相对浓度渐变符合程度以及优化dtw距离基于连通性离群因子算法获取每个待判断体素的平均链接距离；根据每个待判断体素的平均链接距离获取每个待判断体素的连通性离群因子；根据每个待判断体素的连通性离群因子获取每个待判断体素的离群程度；根据每个待判断体素的离群程度判断是否加入感兴趣区域获取核磁共振肿瘤区域。
6.在一实施例中，所述优化后的dtw距离是按照以下步骤获取：将每个待判断体素的波谱数据中代谢物数据点所处位置调整至与中心体素的波谱数据中相对应的代谢物数据点处于同一位置，获取每个待判断体素相对于中心体素相同代谢物数据点的化学位移偏离因子，和每个待判断体素相对于中心体素相同代谢物数据点的距离差异因子；并获取每个待判断体素相对于中心体素不同代谢物数据点的化学位移偏离因子，以及每个待判断体素相对于中心体素不同代谢物数据点的距离差异因子；根据每个待判断体素相对于中心体素相同代谢物数据点的化学位移偏离因子，及每个待判断体素相对于中心体素相同代谢物数据点的距离差异因子，对相同代谢物的数据点距离进行优化；获取相同代谢物数据点优化后的dtw距离；根据每个待判断体素相对于中心体素不同代谢物数据点的化学位移偏离因子，及每个待判断体素相对于中心体素不同代谢物数据点的距离差异因子，对不同代谢物数据点距离进行优化，获取不同代谢物数据点优化后的dtw距离；将相同代谢物数据点优化后的dtw距离和不同代谢物数据点优化后的dtw距离作为优化后的dtw距离。
7.在一实施例中，所述每个待判断体素相对于中心体素相同代谢物数据点的化学位移偏离因子是按照以下步骤获取：在每个待判断体素的波谱数据中代谢物数据点所处位置调整时，获取每个待判断体素的波谱数据中代谢物数据点所处位置相对于中心体素的波谱数据中相对应的代谢物数据点的化学位移差异；将归一化的化学位移差异作为化学位移偏离因子。
8.在一实施例中，所述每个待判断体素相对于中心体素相同代谢物数据点的距离差异因子是按照以下步骤获取：在每个待判断体素的波谱数据中代谢物数据点所处位置调整时，获取每个待判断体素的波谱数据中代谢物数据点在调整前与调整后相对于中心体素的波谱数据中相对应的代谢物数据点的动态时间规整距离差异；将归一化的动态时间规整距离差异作为距离差异因子。
9.在一实施例中，所述每个体素的相对浓度是按照以下步骤获取：将每个体素中代谢物浓度指标的作为每个体素的相对浓度；其中，cho表示代谢物为四甲基胺的浓度；cr表示代谢物为磷酸肌酸和肌酸的浓度。
10.在一实施例中，所述中心体素到每个待判断体素之间所包含的体素的相对浓度形成相对浓度序列，包括：获取中心体素到每个待判断体素之间所包含的体素形成的体素序列；将体素序列中每个体素对应的相对浓度形成的序列作为相对浓度序列；其中，体素序列是将每个待判断体素到中心体素的连线上的所有体素形成的体素序列。
11.在一实施例中，所述核磁共振肿瘤区域是按照以下步骤提取：根据每个待判断体素的离群程度获取每个待判断体素的兴趣程度；设置阈值；当待判断体素的兴趣程度大于阈值，则将该待判断体素加入感兴趣区域；依次将每个待判断体素进行判断，将大于阈值的待判断体素均加入感兴趣区域，获取核磁共振肿瘤区域。
12.在一实施例中，所述感兴趣区域的中心体素是将感兴趣区域的所有波谱数据中代谢物平均浓度最高的体素作为中心体素。
13.进一步地，所述阈值设置为0.7。
14.进一步地，所述感兴趣区域通过人为标注获得。
15.本发明的有益效果是：本发明提供的一种核磁共振肿瘤区域提取方法，该方法通过中心体素与每个待判断体素的波谱数据中代谢物数据点所处位置对中心体素与每个待判断体素的初始dtw距离进行优化获取中心体素与每个待判断体素的优化后的dtw距离，避免了对传统dtw动态规整过程中数据点处于不同代谢物子区间的对应。通过距离的优化使得dtw距离可以保证数据点之间的对应关系，从而获取待检测谱线与感兴趣区域中最近的谱线之间的dtw距离；其次，通过获取每个待判断体素的相对浓度渐变符合程度，相较于传统的通过距离归一化获取的相似程度判断是否将待判断体素放入感兴趣区域，可以通过肿瘤区域磁共振数据中的代谢物浓度信息对相似程度判断过程中的距离进行校正，从而获取到准确的待判断体素磁共振波谱数据与感兴趣区域体素的相似程度。最后，结合每个待判断体素的相对浓度渐变符合程度以及优化的dtw距离，通过基于连通性离群因子算法获取的每个待判断体素的离群程度，使得体素波谱曲线之间的相似性判断更加准确，并通过波谱相似性对感兴趣区域进行扩张判断，在核磁共振检测中提取出更加准确的肿瘤区域。
附图说明
16.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
17.图1为本发明的一种核磁共振肿瘤区域提取方法的实施例总体步骤的流程示意图。
具体实施方式
18.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完
整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
19.本发明的目的是通过感兴趣的肿瘤区域体素的波谱中代谢物峰的空间渐变情况对相似性判断的过程中子区间的确定进行优化，使得体素波谱曲线之间的相似性判断更加准确，并通过波谱相似性将相似的体素集合在一起作为感兴趣区域，直到没有待判断体素可以加入感兴趣区域中，判断结束，将最终获取的感兴趣区域作为提取的肿瘤区域，使得核磁共振检测中提取出更加准确的肿瘤区域。
20.本发明所针对的场景是在对磁共振成像（mri）中感兴趣的肿瘤区域进行波谱数据提取的过程时，因为感兴趣区域的人工标注存在着图像不准确的情况，所以需要进一步对感兴趣区域以外的待判断体素进行相似性判断，通过体素在波谱曲线中的表征进行肿瘤区域的准确提取。在体素波谱的相似性判断过程中，因为传统的序列相似性判断是通过固定长度的子区间进行判断，在判断过程中会因为谱线的代谢物峰偏移造成相似性判断错误的情况。
21.本发明提供的一种核磁共振肿瘤区域提取方法，参见图1所示，包括以下步骤：s1、获取磁共振图像中的感兴趣区域；获取感兴趣区域中的中心体素；将感兴趣区域以外的区域中的体素作为待判断体素；并获取每个体素的相对浓度及波谱数据；在本实施例中，感兴趣区域通过mri图像确定多体素波谱采集过程中的感兴趣区域，通常是通过人工标注的，在此基础上对感兴趣区域以外的体素进行判断。也就是本实施例中所提到的感兴趣区域视为肿瘤区域。
22.在本实施例中，获取波谱数据的过程中，采用的设备为philips 3.0t磁共振扫描仪，采用标准8通道头线圈。常规采用mri检测和t2wi肿瘤横断位薄层进行扫描，结合图像初步定位肿瘤区域，确定多体素磁共振波谱扫描的视野（fov）和感兴趣（voi）区域；多体素磁共振波谱数据采用3d-press序列，参数为tr:2000ms，te 32ms，反转角90度。对于感兴趣区域采集的每个体素的磁共振波谱，通过lcmodel进行量化和基线校准。
23.其中，在本实施例中，每个体素的相对浓度是按照以下步骤获取：将每个体素中代谢物浓度指标的作为每个体素的相对浓度；其中，cho表示代谢物为四甲基胺的浓度；cr表示代谢物为磷酸肌酸和肌酸的浓度。
24.需要说明的是，在感兴趣区域的波谱数据中不同代谢物的峰值会以浓度最高的体素为中心向外逐渐降低，所以将所有波谱数据中代谢物平均浓度最高的体素作为中心体素；将感兴趣区域以外的区域中的体素作为待判断体素；在进行核磁共振波谱的相似性判断时，目的主要是提取出完整的肿瘤区域波谱数据，所以相似性主要取决于对于肿瘤判断最重要的几种代谢物峰的相似性判断，需要通过不同体素的核磁共振波谱中相同代谢物峰的相似性进行判断。其中，在本实施例中，每个待判断体素的波谱数据通常以波谱曲线或者谱线来表示。
25.s2、利用中心体素和每个待判断体素的波谱数据中所有数据点获取中心体素与每个待判断体素的初始dtw距离；根据中心体素与每个待判断体素的波谱数据中代谢物数据点所处位置对中心体
素与每个待判断体素的初始dtw距离进行优化获取中心体素与每个待判断体素的优化后的dtw距离；在本实施例中，中心体素的波谱数据对应的是一个时间序列，而每个待判断体素的波谱数据也对应一个时间序列，为此，基于动态时间规整算法（dtw）获取中心体素与每个待判断体素的波谱数据中所有数据点之间的初始dtw距离。
26.需要说明的是，在磁共振波谱谱线中，不同化学位移量（ppm）位置的代谢物峰代表着相应的代谢物，在dtw距离的计算过程中，因为存在着不同代谢物峰对应的情况，所以需要在动态规整过程将同种代谢物峰进行校正，从而使得在相似程度的衡量过程中dtw距离可以正确衡量两个谱线的相似程度。对于代谢物峰区间的提取因为物理位移等原因会出现波峰的偏移，所以无法通过固定的区间进行提取。对于波谱曲线中的波峰，通过自适应窗口进行波峰整体范围的提取。
27.在本实施例中，通过自适应窗口对每个待判断体素的波谱曲线中的波峰在整体范围提取的具体过程为：首先，利用窗口右边线作为起始点从谱线最左侧（化学位移量（ppm）最高的位置）开始向右开始进行判断；其次，经过基线校正谱线会呈现为基线稳定的谱线，对于每个ppm值对应的谱信号强度，检测其相较于前一个数据点是否上升。然后，当检测到一个波动的上升，下降直到回归基线则窗口截止，并将此窗口所框选的谱线对应的数据记为第个待判断体素的第个谱线波动窗口。也就是利用一个自适应窗口将谱线中的代谢物波峰进行框选，作为第个待判断体素的第个波动窗口；同时将每个框选波峰的窗口对应的数据作为谱线的子区间；最后，检测谱线中的全部波动（只考虑从基线上升的波动）。
28.为此，获取到的每个待判断体素的核磁共振波谱中多个波动窗口，对于每个波动需要通过其所在的化学位移（ppm）位置确定其所代表的代谢物。在脑部肿瘤中，主要确定如下几种：cho：四甲基胺，位于3.2ppm。 cr：磷酸肌酸和肌酸，位于3.0ppm。 naa：n-乙酰基基团，位于2.0ppm附近。 ala：丙氨酸，一种中间代谢的代谢产物，位于1.4ppm。 la：乳酸盐，一种间接反映由于异常酶促反应产生的厌氧糖酵解或者糖酵解的指标，位于1.3ppm。通过每个窗口中波动峰值所在位置的最近邻特征ppm确定波动所对应的代谢物。获取到第个体素的核磁共振波谱中第个波动窗口所对应的代谢物。对于待对比相似程度的两个核磁共振波谱谱线，通过代谢物的对应关系确定一一对应的关系。
29.在本实施例中，根据每个待判断体素的波谱数据中代谢物数据点所处位置与中心体素的波谱数据中代谢物数据点的相对位置关系，对初始dtw距离进行优化获取优化后的dtw距离。
30.具体的，优化后的dtw距离是按照以下步骤获取：将每个待判断体素的波谱数据中代谢物数据点所处位置调整至与中心体素的波谱数据中相对应的代谢物数据点处于同一位置，获取每个待判断体素相对于中心体素相同代谢物数据点的化学位移偏离因子，和每个待判断体素相对于中心体素相同代谢物数据点的距离差异因子；并获取每个待判断体素相对于中心体素不同代谢物数据点的化学位移偏离因子，以及每个待判断体素相对于中心体素不同代谢物数据点的距离差异因子；根据每个待判断体素相对于中心体素相同代谢物数据点的化学位移偏离因子，及每个待判断体素相对于中心体素相同代谢物数据点的距离差异因子，对相同代谢物的数据
点距离进行优化；获取相同代谢物数据点优化后的dtw距离；根据每个待判断体素相对于中心体素不同代谢物数据点的化学位移偏离因子，及每个待判断体素相对于中心体素不同代谢物数据点的距离差异因子，对不同代谢物数据点距离进行优化，获取不同代谢物数据点优化后的dtw距离；将相同代谢物数据点优化后的dtw距离和不同代谢物数据点优化后的dtw距离作为优化后的dtw距离。
31.在本实施例中，通过获取的中心体素谱线中代谢物峰与每个待判断体素谱线的对应关系后，因为代谢物峰之间存在着化学位移的差异，所以在对代谢物峰之间进行相似性对比时需要将其进行配准。对于中心体素与待判断体素波谱数据中对应的波动窗口，也就是波谱数据中对应的相同代谢物波峰。需要说明的是，两个对应的波动窗口对应于波谱数据中的是两个子序列，同时两个对应的波动窗口对应两个子区间；通过dtw动态规整获取两个子序列中数据点的最近对应关系。将待判断体素的代谢物峰的波峰数据点调整为与中心体素的代谢物峰的波峰处于同一位置，并在调整之后在其两侧通过基线数值填充将两个子区间填充为等长。获取配准后的子区间中两个波动的dtw距离。波峰位置调整的ppm距离即为曲线化学位移差异，配准前和后的子区间序列dtw距离差异，对于谱线中每个对应子区间之间的ppm差异与配准前后的dtw距离差异获取待检测体素磁共振波谱的化学位移偏离因子与子区间距离差异因子。
32.具体的，每个待判断体素相对于中心体素相同代谢物数据点的化学位移偏离因子是按照以下步骤获取：在每个待判断体素的波谱数据中代谢物数据点所处位置调整时，获取每个待判断体素的波谱数据中代谢物数据点所处位置相对于中心体素的波谱数据中相对应的代谢物数据点的化学位移差异；将归一化的化学位移差异作为化学位移偏离因子。
33.每个待判断体素相对于中心体素相同代谢物数据点的距离差异因子是按照以下步骤获取：在每个待判断体素的波谱数据中代谢物数据点所处位置调整时，获取每个待判断体素的波谱数据中代谢物数据点在调整前与调整后相对于中心体素的波谱数据中相对应的代谢物数据点的动态时间规整距离差异；将归一化的动态时间规整距离差异作为距离差异因子。
34.在本实施例中，化学位移偏离因子与子区间距离差异因子计算公式如下：计算公式如下：式中，表示每个待判断体素中的第个波动窗口相对于中心体素的化学位移偏离因子；表示每个待判断体素中的第个波动窗口相对于中心体素的距离差异因子；表示每个待判断体素的波谱数据中代谢物数据点（第个波动窗口）所处位置相对于中心体素的波谱数据中相对应的代谢物数据点（第个波动窗口）的化学位移差异；表示每个待判断体素的波谱数据中代谢物数据点（第个波动窗口）在调整前与调整后相对于中心体素的波谱数据中相对应的代谢物数据点（第个波动窗口）的动态时间规整距离差异；表示归一化函数；需要说明的是，在获取到中心体素与待判断体素的磁共振波
谱曲线对应的第个波动窗口之间的两种差异因子之后，由于通过距离矩阵进行动态规整的过程中，因为不同代谢物峰之间的距离在现有的计算中还是同样的相对距离计算，所以需要对其通过化学位移偏离因子与子区间距离差异因子进行校正，使得当数据点处于不同的子区间时，数据点之间的距离会更大一些，当数据点处于相同的子区间时，数据点之间的距离会根据波峰的偏移进行校正。
35.具体的，根据每个待判断体素相对于中心体素相同代谢物数据点的化学位移偏离因子，及每个待判断体素相对于中心体素相同代谢物数据点的距离差异因子，对相同代谢物的数据点距离进行优化；获取相同代谢物数据点优化后的dtw距离；根据每个待判断体素相对于中心体素不同代谢物数据点的化学位移偏离因子，及每个待判断体素相对于中心体素不同代谢物数据点的距离差异因子，对不同代谢物数据点距离进行优化，获取不同代谢物数据点优化后的dtw距离。
36.在本实施例中，根据每个待判断体素相对于中心体素相同代谢物数据点的化学位移偏离因子，及每个待判断体素相对于中心体素相同代谢物数据点的距离差异因子，对相同代谢物的数据点距离进行优化；具体对于中心体素与待判断体素中相同代谢物的子区间之间的数据点的优化距离计算公式如下：对于中心体素与待判断体素中不同代谢物的子区间之间的数据点的优化距离计算公式如下：式中，表示第个体素中的第个数据点与中心体素中相同代谢物中数据点的初始dtw距离；表示每个待判断体素中第个波动窗口与中心体素的谱线中第个波动窗口对应代谢物与当前代谢物之间所有的波动窗口；表示每个待判断体素中的第个波动窗口相对于中心体素的化学位移偏离因子；表示每个待判断体素中的第个波动窗口相对于中心体素的距离差异因子；表示中心体素与待判断体素中相同代谢物的子区间之间的数据点的优化的dtw距离。为此，通过将待判断体素的代谢物峰的波峰数据点调整为与中心体素代谢物峰的波峰处于同一位置，对待判断体素中的数据点与中心体素中的数据点之间的dtw距离进行优化获取优化后的dtw距离，也就是通过对dtw动态规整过程中两个谱线之间数据点的对应距离的计算通过化学位移因子与距离差异因子进行校正，避免了对传统dtw动态规整过程中数据点处于不同代谢物子区间的对应。通过距离的校正使得dtw距离可以确保数据点之间的对应关系保留于子区间的对应关系，从而获取待判断体素与中心体素中谱线之间的dtw距离，并可以据此进行接下来的谱线相似度衡量。
37.s3、获取中心体素到每个待判断体素之间所包含的体素的相对浓度形成相对浓度序列，获取相对浓度序列的差分序列；根据中心体素到每个待判断体素之间所包含的体素数量，中心体素到每个待判断
体素之间所包含的体素相对浓度连续梯度变化的体素数量，以及相对浓度序列的差分序列，获取每个待判断体素的相对浓度渐变符合程度；具体的，所述中心体素到每个待判断体素之间所包含的体素的相对浓度形成相对浓度序列，包括：获取中心体素到每个待判断体素之间所包含的体素形成的体素序列；将体素序列中每个体素对应的相对浓度形成的序列作为相对浓度序列；其中，体素序列是将每个待判断体素到中心体素的连线上的所有体素形成的体素序列。
38.需要说明的是，通过对中心体素与每个待判断体素的数据点之间的距离计算进行优化，在获取到体素对应的磁共振波谱谱线之间的距离之后，对于视野域以内，感兴趣区域以外的体素，也就是所有待判断的体素，则可以进一步通过该场景下肿瘤区域的代谢物浓度特征对相似性的度量进行优化。在传统的序列相似性度量中，需要通过全部序列之间的距离进行归一化，从而获取序列之间的相似程度。但是在这个相似程度的判断过程中只考虑到谱线之间的序列距离。在磁共振波谱谱线中，对于肿瘤区域的提取还需要根据体素对应的磁共振波谱中代谢物浓度。通过浓度渐变过程的符合程度对通过谱线之间距离获取到的体素相似程度进行优化，从而获取更准确的肿瘤区域。在本实施例中，首先，通过代谢物浓度判断方法获取感兴趣区域中代谢物平均浓度最高的体素，获取代谢物浓度最高的体素与待判断体素连线上的所有体素，从浓度最高的体素开始获取每个体素与前一个体素在磁共振波谱谱线上的距离，也就是通过优化的dtw距离。其次，通过每个待判断体素对应谱线中作为代谢物浓度指标的作为每个待判断体素的相对浓度，确定浓度渐变的符合程度。最后，通过浓度渐变符合程度对距离在序列中的连通性离群检测过程中的序列距离进行调整，从而获取准确的待判断体素的离群程度，通过1减去离群程度获取该体素的兴趣程度。
39.在本实施例中，每个待判断体素的相对浓度渐变符合程度计算公式如下：式中，表示第个待判断体素的相对浓度渐变符合程度；表示中心体素到第个待判断体素之间所包含的体素的数量；表示中心体素到第个待判断体素之间所包含的体素相对浓度连续梯度变化的体素数量，其中，连续梯度变化的体素数量是通过对体素序列中从第个待判断体素向中心体素判断相对浓度差分的正负号，将连续为负号的数量作为相对浓度连续梯度变化的体素数量；表示差分序列中第t个元素值；表示差分序列中的元素数量；表示第个待判断体素与前一个体素的相对浓度差异值；表示归一化函数。为此，通过连续为负号的体素数量衡量渐变过程的渐变程度，对于序列中的浓度数值，通过将浓度序列差分获取到差分序列，对于差分序列其表示每两个体素之间的浓度差异，通过待判断体素与前一个体素的差异值与浓度差分序列的平均差异衡量体素对应的磁共振波谱中的代谢物浓度渐变状态的符合程度。
40.需要说明的是，对于浓度序列中的体素对应磁共振波谱中的代谢物浓度信息，通过浓度值序列的差分序列的符号来判断体素序列是否停止渐变，并根据待判断体素的差分数值与序列差分均值的差异判断待判断体素位置是否已经渐变平稳，从而确定待判断体素
的渐变符合程度。为此，通过待判断体素在序列中磁共振波谱数据中的代谢物浓度信息获取到渐变符合程度。相较于传统的通过距离归一化获取到的相似程度判断是否将待判断体素放入感兴趣区域可以通过肿瘤区域磁共振数据中的代谢物浓度信息对相似程度判断过程中的距离进行校正，从而获取到准确的待判断体素磁共振波谱数据与感兴趣区域体素的相似程度。
41.s4、根据中心体素到每个待判断体素之间所包含的体素形成的体素序列、每个待判断体素的相对浓度渐变符合程度以及优化dtw距离基于连通性离群因子算法获取每个待判断体素的平均链接距离；根据每个待判断体素的平均链接距离获取每个待判断体素的连通性离群因子；根据每个待判断体素的连通性离群因子获取每个待判断体素的离群程度；需要说明的是，在获取到体素磁共振波谱序列中的代谢物浓度的渐变符合程度之后，就可以在获取待判断体素在体素序列中的基于连通性的离群程度的过程中对每个数据点的路径距离进行校正。对于待判断体素的兴趣程度，可以通过其在体素序列中的离群程度进行判断，当该体素对应的磁共振波谱谱线在体素序列中基于连通性的离群程度越高，则说明该波谱谱线越不能放入感兴趣区域，也就是兴趣程度越低。
42.对于连通性的离群因子（cof）的计算，将体素序列中的体素数量作为cof计算过程中的k值，即每个体素都通过全部体素进行离群程度的计算。而对于体素序列中每两个体素之间通过优化后的dtw距离作为其路径距离。而在通过路径距离计算离群因子的过程中，因为体素对应的浓度渐变符合程度越高，则说明该体素的波谱谱线与前一个体素的谱线之间的距离因为浓度对于渐变过程的符合应该被校正的更低。通过渐变符合程度对到达该谱线的路径距离进行校正，从而使得对于每个体素在谱线序列中的离群因子更加准确。
43.具体的，每个待判断体素的平均链接距离，对于第个待判断体素优化后的平均链接距离计算方法如下：式中，表示序列中的第个体素；表示局部sbn路径中的数据点索引；表示第个路径；表示第个路径的路径长度，即谱线之间优化后的dtw距离；表示第个体素相对浓度渐变符合程度；表示对于序列中的数据点将其作为sbn路径中的第一个数据点根据对应的，即邻域范围中的平均链接距离；此处k邻域范围即为序列中的全部数据。需要说明的是，通过cof离群因子中平均链接距离的计算方式，在原有的基础上通过序列中体素对应的磁共振波谱谱线中的代谢物浓度信息对于感兴趣区域的渐变特征的符合程度对路径距离，即优化后的dtw距离进行校正。从而使得在衡量离群因子的过程中可以避免单一的谱线之间的距离差异大但是浓度符合肿瘤区域渐变而导致的满足条件的体素被判断为无法加入兴趣区域。为此，通过每个体素对应的磁共振波谱谱线中的代谢物浓度符合程度对序列离群衡量过程中到达该体素的路径长度进行校正，相较于传统的序列中平均链接距离，可以根据浓度符合程度来拉伸路径距离，使得浓度符合渐变的体素
在谱线之间距离较大时也可以被认为可以加入兴趣区域。在本实施例中，通过获取平均链接距离之后，利用该平均链接距离进一步计算每个待判断体素在序列中对应的每个待判断体素的连通性离群因子，通过每个待判断体素的每个待判断体素的连通性离群因子进行归一化获取每个待判断体素的离群程度。
44.s5、根据每个待判断体素的离群程度判断是否加入感兴趣区域获取核磁共振肿瘤区域。
45.具体的，核磁共振肿瘤区域是按照以下步骤提取：根据每个待判断体素的离群程度获取每个待判断体素的兴趣程度；设置阈值；当待判断体素的兴趣程度大于阈值，则将该待判断体素加入感兴趣区域；依次将每个待判断体素进行判断，将大于阈值的待判断体素均加入感兴趣区域，获取以中心体素为中心的核磁共振肿瘤区域。
46.在本实施例中，获取每个待判断体素的离群程度后，对于待判断的体素，通过减去离群程度获取其兴趣程度，也就是其是否可以加入兴趣区域的指标。对于判断标准，通过待判断体素之外所有体素的兴趣程度的均值作为阈值，当待判断体素的兴趣程度高于该阈值则将这个体素加入感兴趣区域，将最终获取的感兴趣区域作为提取的肿瘤区域，使得核磁共振检测中提取出更加准确的肿瘤区域。在本发明一个实施例中，阈值设置为0.7。
47.以上仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

技术特征：
1.一种核磁共振肿瘤区域提取方法，其特征在于，包括以下步骤：获取磁共振图像中的感兴趣区域；获取感兴趣区域中的中心体素；将感兴趣区域以外的区域中的体素作为待判断体素；并获取每个体素中四甲基胺相对于磷酸肌酸和肌酸的相对浓度，及波谱数据；利用中心体素和每个待判断体素的波谱数据中所有数据点获取中心体素与每个待判断体素的初始dtw距离；根据中心体素与每个待判断体素的波谱数据中代谢物数据点所处位置对中心体素与每个待判断体素的初始dtw距离进行优化获取中心体素与每个待判断体素的优化后的dtw距离；获取中心体素到每个待判断体素之间所包含的体素的相对浓度形成相对浓度序列，获取相对浓度序列的差分序列；根据中心体素到每个待判断体素之间所包含的体素数量，中心体素到每个待判断体素之间所包含的体素相对浓度连续梯度变化的体素数量，以及相对浓度序列的差分序列，获取每个待判断体素的相对浓度渐变符合程度；根据中心体素到每个待判断体素之间所包含的体素形成的体素序列、每个待判断体素的相对浓度渐变符合程度以及优化dtw距离基于连通性离群因子算法获取每个待判断体素的平均链接距离；根据每个待判断体素的平均链接距离获取每个待判断体素的连通性离群因子；根据每个待判断体素的连通性离群因子获取每个待判断体素的离群程度；根据每个待判断体素的离群程度判断是否加入感兴趣区域获取核磁共振肿瘤区域。2.根据权利要求1所述的核磁共振肿瘤区域提取方法，其特征在于，所述优化后的dtw距离是按照以下步骤获取：将每个待判断体素的波谱数据中代谢物数据点所处位置调整至与中心体素的波谱数据中相对应的代谢物数据点处于同一位置，获取每个待判断体素相对于中心体素相同代谢物数据点的化学位移偏离因子，和每个待判断体素相对于中心体素相同代谢物数据点的距离差异因子；并获取每个待判断体素相对于中心体素不同代谢物数据点的化学位移偏离因子，以及每个待判断体素相对于中心体素不同代谢物数据点的距离差异因子；根据每个待判断体素相对于中心体素相同代谢物数据点的化学位移偏离因子，及每个待判断体素相对于中心体素相同代谢物数据点的距离差异因子，对相同代谢物的数据点距离进行优化；获取相同代谢物数据点优化后的dtw距离；根据每个待判断体素相对于中心体素不同代谢物数据点的化学位移偏离因子，及每个待判断体素相对于中心体素不同代谢物数据点的距离差异因子，对不同代谢物数据点距离进行优化，获取不同代谢物数据点优化后的dtw距离；将相同代谢物数据点优化后的dtw距离和不同代谢物数据点优化后的dtw距离作为优化后的dtw距离。3.根据权利要求2所述的核磁共振肿瘤区域提取方法，其特征在于，所述每个待判断体素相对于中心体素相同代谢物数据点的化学位移偏离因子是按照以下步骤获取：在每个待判断体素的波谱数据中代谢物数据点所处位置调整时，获取每个待判断体素的波谱数据中代谢物数据点所处位置相对于中心体素的波谱数据中相对应的代谢物数据
点的化学位移差异；将归一化的化学位移差异作为化学位移偏离因子。4.根据权利要求2所述的核磁共振肿瘤区域提取方法，其特征在于，所述每个待判断体素相对于中心体素相同代谢物数据点的距离差异因子是按照以下步骤获取：在每个待判断体素的波谱数据中代谢物数据点所处位置调整时，获取每个待判断体素的波谱数据中代谢物数据点在调整前与调整后相对于中心体素的波谱数据中相对应的代谢物数据点的动态时间规整距离差异；将归一化的动态时间规整距离差异作为距离差异因子。5.根据权利要求1所述的核磁共振肿瘤区域提取方法，其特征在于，所述每个体素的相对浓度是按照以下步骤获取：将每个体素中代谢物浓度指标的作为每个体素的相对浓度；其中，cho表示代谢物为四甲基胺的浓度；cr表示代谢物为磷酸肌酸和肌酸的浓度。6.根据权利要求1所述的核磁共振肿瘤区域提取方法，其特征在于，所述中心体素到每个待判断体素之间所包含的体素的相对浓度形成相对浓度序列，包括：获取中心体素到每个待判断体素之间所包含的体素形成的体素序列；将体素序列中每个体素对应的相对浓度形成的序列作为相对浓度序列；其中，体素序列是将每个待判断体素到中心体素的连线上的所有体素形成的体素序列。7.根据权利要求1所述的核磁共振肿瘤区域提取方法，其特征在于，所述核磁共振肿瘤区域是按照以下步骤提取：根据每个待判断体素的离群程度获取每个待判断体素的兴趣程度；设置阈值；当待判断体素的兴趣程度大于阈值，则将该待判断体素加入感兴趣区域；依次将每个待判断体素进行判断，将大于阈值的待判断体素均加入感兴趣区域，获取核磁共振肿瘤区域。8.根据权利要求1所述的核磁共振肿瘤区域提取方法，其特征在于，所述感兴趣区域的中心体素是将感兴趣区域的所有波谱数据中代谢物平均浓度最高的体素作为中心体素。9.根据权利要求7所述的核磁共振肿瘤区域提取方法，其特征在于，所述阈值设置为0.7。10.根据权利要求1所述的核磁共振肿瘤区域提取方法，其特征在于，所述感兴趣区域通过人为标注获得。

技术总结
本发明公开了一种核磁共振肿瘤区域提取方法，属于图像数据处理技术领域；包括以下步骤：获取磁共振图像中的感兴趣区域；获取每个待判断体素的相对浓度渐变符合程度；根据每个待判断体素的平均链接距离获取每个待判断体素的连通性离群因子；根据每个待判断体素的连通性离群因子获取每个待判断体素的离群程度；根据每个待判断体素的离群程度判断是否加入感兴趣区域获取核磁共振肿瘤区域。本发明通过对感兴趣区域边界的相邻体素进行代谢物相似性判断，在感兴趣区域基础上根据波谱相似性进行扩展，使得对感兴趣的肿瘤区域的波谱数据提取更加精确。取更加精确。取更加精确。


技术研发人员：陈博 高燕 王晓勇 赵谦 赵广明 祁长银 白易民
受保护的技术使用者：山东奥洛瑞医疗科技有限公司
技术研发日：2023.08.14
技术公布日：2023/9/13