CN102103745A - 一种数字减影的软件实现方法 - Google Patents

Abstract

本发明涉及一种数字减影的软件实现方法，针对DICOM格式多帧检查影像，通过打开DICOM格式的多帧医学检查影像，装载指定帧的像素，并显示影像；获取用户指定的减影蒙片帧，或者减影图像文件；遍历蒙片与当前帧，根据指定减影区域，进行像素减影计算，获取当前帧与蒙片影像减影后效果；针对多帧影像进行减影正，倒序，循环播放，实现动态减影阅片。

Description

一种数字减影的软件实现方法

技术领域

本发明涉及一种对于多帧医学影像的数字减影实现方法，尤其是针对DICOM格式多帧检查影像的数字减影实现方法。

背景技术

在医学影像检查中，核磁检查，C型臂X光机，血管造影等影像检查设备会产生多帧连续检查影像。为了能够清晰地观察检查部位血管的实际情况，观察血管病变，血管狭窄的定位测量，诊断，及辅助介入治疗，医生通过在检查的过程中为病人注射(或者服用)血管造影剂来使得血管影像在放射照射下凸显出来，但是由于血管通常位于组织器官的内部，即使在造影放射检查下的影像，仍然会被各种不同密度的组织器官所遮挡，不能清晰地观察独立的血管影像，尤其是较细的血管影像会混淆在组织中难以通过检查影像看清楚。

数字减影(DSA Digital Subtration angiography)技术通常通过大型放射医学影像设备实时实现，但对于针对血管病情的检查诊断，对于血管位置的测量，对于介入治疗的手术规划，对于医学会诊教学，医生需要在专业的医学影像工作站上通过软件方式对于多帧检查影像进行数字减影操作，以观察，测量血管病症的详细情况，以辅助诊断，治疗。

发明内容

一种数字减影的软件实现方法，针对DICOM多帧影像进行数字减影运算，实现数字减影的影像效果，其实现方法步骤如下，实现流程如图3：

第1步：打开DICOM格式的多帧医学检查影像，装载指定帧的像素，并显示影像。

第2步：获取用户指定的减影蒙片帧，或者减影图像文件。

第3步：遍历蒙片与当前帧，根据指定减影区域，进行像素减影计算，减影前后的当前帧效果示例分别如图1，图2。

第4步：针对多帧影像进行减影正，倒序，循环播放，实现动态减影阅片。

附图说明

图1：减影前血管影像

图2：减影后血管影像

图3：数字减影实现流程

具体实施方法

一种数字减影的软件实现方法，其实现方法步骤如下，实现流程如图3：

第1步：打开DICOM格式的多帧医学检查影像，装载指定第i帧的像素，并显示影像：读取DICOM文件TAG(0028，0008)Number of Frames的值，获取多帧文件的总帧数，通过读取TAG(FFFE，E000)Image Fragment作为像素区分隔标记来读取指定帧的像素数据，并显示该帧影像。

第2步：获取用户指定的减影蒙片帧，或者减影图像文件。

如果用户指定当前多帧文件的另外一帧作为减影蒙片帧，则按照步骤1方法装载蒙片帧的像素数据，而指定减影区域座标则自动设置为左上(0，0，Height，Width)。果用户指定一个单独的影像文件作为减影蒙片，则装载该蒙片影像文件像素，并获取用户输入的减影区域座标(x，y)。

第3步：遍历蒙片与当前帧，根据指定减影区域，进行像素减影计算：

遍历当前帧影像像素，针对每个像素P_cur(x，y)做如下判断：

X_DSALeft＜X_cur＜X_DSARight与Y_DSATop＜Y_cur＜Y_DSABottom来判断当前像素是否在指定减影区域之内。

上式中各符号的含义如下：

X_DSALeft： DSA减影区域的左边界X座标；

X_cur：     当前帧影像的当前像素X座标；

X_DSARight：DSA减影区域的右边界X座标；

Y_DSALeft： DSA减影区域的左边界Y座标；

Y_cur：     当前帧影像的当前像素Y座标；

Y_DSARight：DSA减影区域的右边界Y座标；

然后通过：

P_{curDSA(x，y)}＝P_{cur(x+DSALeft，y+DSATop)}-Pcover(x，y)

上式中各符号含义如下：

P_{curDSA(x，y)}：            当前帧X，Y座标像素减影后的像素值；

P_{cur(x+DSALeft，y+DSATop)}：当前帧在减影区域中的相对应座标的像素值；

P_cover(x，y)：             蒙片帧影像x，y座标像素的像素值；

得到减影后的像素集合，减影前后的当前帧示例效果分别如图1，图2。第4步：针对多帧影像进行减影正，倒序，循环播放，实现动态减影阅片。

通过获取当前播放方向的下一帧影像，装载该帧像素，并针对该帧影像执行第3步减影计算，并在当前帧率控制帧间暂停时间，实现动态减影阅片播放。

Claims (4)

1.一种数字减影的软件实现方法，针对DICOM多帧影像通过蒙片影像与当前帧影像进行数字减影运算，实现数字减影的影像效果，其特征在于按照如下步骤进行：

(1)打开DICOM格式的多帧医学检查影像，装载指定帧的像素；

(2)获取用户指定的减影蒙片帧，或者减影图像文件；

(3)遍历蒙片与当前帧，根据指定减影区域，进行像素减影计算；

(4)针对多帧影像进行减影正，倒序，循环播放，实现动态减影阅片。

2.根据上述权利要求1所述的一种数字减影的软件实现方法，其特征在于通过如下方法判断当前像素是否在指定减影区域之内：

遍历当前帧影像像素，针对每个像素Pcur(x，y)做如下是否在减影区域之内的判断：

X_DSALeft＜X_cur＜X_DSARight与Y_DSATop＜Y_cur＜Y_DSABottom

上式中各符号的含义如下：

X_DSALeft：DSA减影区域的左边界X座标；

X_cur：当前帧影像的当前像素X座标；

X_DSARight：DSA减影区域的右边界X座标；

Y_DSALeft：DSA减影区域的左边界Y座标；

Y_cur：当前帧影像的当前像素Y座标；

Y_DSARight：DSA减影区域的右边界Y座标。

3.根据上述权利要求1所述的一种数字减影的软件实现方法，其特征在于通过如下方法计算当前影像与减影蒙片影像对应像素之间的减影后像素值：

P_{curDSA(x，y)}＝P_{cur(x+DSALeft，y+DSATop)}-P_cover(x，y)

上式中各符号含义如下：

P_{curDSA(x，y)}：当前帧X，Y座标像素减影后的像素值；

P_{cur(x+DSALeft，y+DSATop)}：当前帧在减影区域中的相对应座标的像素值；

P_cover(x，y)：蒙片帧影像x，y座标像素的像素值。

4.根据上述权利要求1所述的一种数字减影的软件实现方法，其特征在于通过如下方法实现针对多帧影像进行减影正，倒序，循环播放动态减影阅片：

通过获取当前播放方向的下一帧影像，装载该帧像素，并针对该帧影像执行减影计算，并在当前帧率控制帧间暂停时间，实现动态减影阅片播放。 

Images