CN116797683A - 一种数字减影血管造影控制系统及控制方法 - Google Patents
一种数字减影血管造影控制系统及控制方法 Download PDFInfo
- Publication number
- CN116797683A CN116797683A CN202310785144.1A CN202310785144A CN116797683A CN 116797683 A CN116797683 A CN 116797683A CN 202310785144 A CN202310785144 A CN 202310785144A CN 116797683 A CN116797683 A CN 116797683A
- Authority
- CN
- China
- Prior art keywords
- image
- contrast enhancement
- blood vessel
- artifact
- contrast
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
- 238000002583 angiography Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 39
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 115
- 238000010586 diagram Methods 0.000 claims abstract description 44
- 230000002792 vascular Effects 0.000 claims description 23
- 230000006870 function Effects 0.000 claims description 21
- 238000004590 computer program Methods 0.000 claims description 12
- 238000003860 storage Methods 0.000 claims description 11
- 238000001514 detection method Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 3
- 230000010363 phase shift Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 6
- 238000012545 processing Methods 0.000 description 8
- 238000004891 communication Methods 0.000 description 7
- 239000002872 contrast media Substances 0.000 description 5
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 230000004075 alteration Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 206010002329 Aneurysm Diseases 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- 208000032170 Congenital Abnormalities Diseases 0.000 description 1
- 206010061619 Deformity Diseases 0.000 description 1
- 208000007536 Thrombosis Diseases 0.000 description 1
- 206010057469 Vascular stenosis Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 231100000216 vascular lesion Toxicity 0.000 description 1
- 238000007631 vascular surgery Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Signal Processing (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Physiology (AREA)
- Public Health (AREA)
- Vascular Medicine (AREA)
- Cardiology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本申请提供了一种数字减影血管造影控制系统及控制方法，通过确定对比增强图像的频率域图，将该对比增强图像的频率域图转换为镜像频率域图，根据该对比增强图像的镜像频率域图确定血管伪影角度，根据所述血管伪影角度对该对比增强图像进行旋转得到水平对比增强图像，由对比增强曲线图确定血管伪影长度；根据所述血管伪影角度和所述血管伪影长度确定血管伪影模糊核；将所述对比增强图像进行对称填充，得到该对比增强图像的扩展图像，根据该对比增强图像的扩展图像和所述血管伪影模糊核对对比增强图像进行校正，将校正后的对比增强图像与所述基线图像进行数字减法运算生成数字减影图像，可降低数字减影血管造影过程中的血管伪影干扰。
Description
技术领域
本申请涉及数字减影血管造影技术领域，更具体的说，本申请涉及一种数字减影血管造影控制系统及控制方法。
背景技术
数字减影血管造影是用于观察和评估人体血管系统的结构和功能的技术，它在临床上广泛应用于多个领域，包括心血管病学、神经学、血管外科等，数字减影血管造影可以用于检测血管狭窄、畸形、血栓形成、动脉瘤等血管病变，辅助医生做出准确的诊断和治疗决策。
现有技术中，通过数字减影血管造影机生成数字减影图像时，因为患者的呼吸、心跳或手部远动等因素导致最终生成的数字减影图像中会产生与血管图像不相关的血管伪影，血管伪影的产生会干扰医生对数字减影图像的解读，严重的甚至可能会导致误诊，如何降低数字减影血管造影过程中的血管伪影干扰是业界面临的问题。
发明内容
本申请实施例提供一种数字减影血管造影控制系统及控制方法，以降低数字减影血管造影过程中的血管伪影干扰。
为解决上述技术问题，本申请采用如下技术方案：
第一方面，本申请提供一种数字减影血管造影控制方法，包括如下步骤：
获取数字减影血管造影的基线图像；
采集数字减影血管造影的对比增强图像；
确定所述对比增强图像的频率域图，将该对比增强图像的频率域图转换为镜像频率域图，根据该对比增强图像的镜像频率域图确定血管伪影角度；
根据所述血管伪影角度对该对比增强图像进行旋转，得到水平对比增强图像，由所述水平对比增强图像通过一阶差分得到对比增强梯度图像，根据所述对比增强梯度图像得到对比增强曲线图，由所述对比增强曲线图确定血管伪影长度；
根据所述血管伪影角度和所述血管伪影长度确定血管伪影模糊核；
将所述对比增强图像进行对称填充，得到该对比增强图像的扩展图像，根据该对比增强图像的扩展图像和所述血管伪影模糊核对对比增强图像进行校正，将校正后的对比增强图像与所述基线图像进行数字减法运算生成数字减影图像。
在一些实施例中，根据该对比增强图像的镜像频率域图确定血管伪影角度具体可包括：
确定该对比增强图像的镜像频率域图中血管伪影的边缘轮廓；
对于确定的血管伪影的边缘轮廓经过直线拟合确定血管伪影角度。
在一些实施例中，采用最小二乘法进行直线拟合。
在一些实施例中，血管伪影的边缘轮廓值根据下述公式确定：
其中PC(a，i)表示对比增强图像的镜像频率域图中血管伪影的边缘轮廓值，a表示对比增强图像的镜像频率域图的像素值，i表示不同的检测方向，Wi(a)表示在i方向的权重因子，P0表示噪声能量阈值，α为预设的常数，Δφni(a)表示相位偏离函数，表示如该符号中的值为正，则结果为其本身，否则为零，Ani(a)表示在i方向上的幅度值。
在一些实施例中，根据所述血管伪影角度和所述血管伪影长度确定血管伪影模糊核具体包括：
根据所述血管伪影角度和所述血管伪影长度确定血管伪影点扩散函数；
对确定的血管伪影点扩散函数进行逆滤波得到血管伪影模糊核。
在一些实施例中，根据该对比增强图像的扩展图像和所述血管伪影模糊核对对比增强图像进行校正具体包括：
将该对比增强图像的扩展图像与血管伪影模糊核进行卷积操作；
根据卷积结果，从卷积结果中提取校正后的对比增强图像。
在一些实施例中，将校正后的对比增强图像与所述基线图像进行数字减法运算生成数字减影图像具体包括：
将校正后的对比增强图像与基线图像进行空间上的配准，以确保两组图像对应的血管区域位置一致；
将配准后的对比增强图像和基线图像进行对齐，使得两组图像的血管结构在空间上完全重合；
将配准及对齐的基线图像从对比增强图像中减去，得到生成数字减影图像。
第二方面，本申请提供一种数字减影血管造影控制系统，其包括有：
基线图像获取模块，用于获取数字减影血管造影的基线图像；
对比增强图像采集模块，用于采集数字减影血管造影的对比增强图像；
血管伪影角度确定模块，用于确定所述对比增强图像的频率域图，将该对比增强图像的频率域图转换为镜像频率域图，根据该对比增强图像的镜像频率域图确定血管伪影角度；
血管伪影长度确定模块，用于根据所述血管伪影角度对该对比增强图像进行旋转，得到水平对比增强图像，由所述水平对比增强图像通过一阶差分得到对比增强梯度图像，根据所述对比增强梯度图像得到对比增强曲线图，由所述对比增强曲线图确定血管伪影长度；
血管伪影模糊核确定模块，用于根据所述血管伪影角度和所述血管伪影长度确定血管伪影模糊核；
数字减影图像生成模块，用于将所述对比增强图像进行对称填充，得到该对比增强图像的扩展图像，根据该对比增强图像的扩展图像和所述血管伪影模糊核对对比增强图像进行校正，将校正后的对比增强图像与所述基线图像进行数字减法运算生成数字减影图像。
第三方面，本申请提供一种计算机设备，所述计算机设备包括存储器和处理器，所述存储器存储有代码，所述处理器被配置为获取所述代码，并执行上述的数字减影血管造影控制方法。
第四方面，本申请提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现上述的数字减影血管造影控制方法。
本申请实施例提供的技术方案具有以下有益效果：
本申请公开的一种数字减影血管造影控制系统及控制方法中，首先获取数字减影血管造影的基线图像；采集数字减影血管造影的对比增强图像；确定所述对比增强图像的频率域图，将该对比增强图像的频率域图转换为镜像频率域图，根据该对比增强图像的镜像频率域图确定血管伪影角度；根据所述血管伪影角度对该对比增强图像进行旋转，得到水平对比增强图像，由所述水平对比增强图像通过一阶差分得到对比增强梯度图像，根据所述对比增强梯度图像得到对比增强曲线图，由所述对比增强曲线图确定血管伪影长度；根据所述血管伪影角度和所述血管伪影长度确定血管伪影模糊核；将所述对比增强图像进行对称填充，得到该对比增强图像的扩展图像，根据该对比增强图像的扩展图像和所述血管伪影模糊核对对比增强图像进行校正，将校正后的对比增强图像与所述基线图像进行数字减法运算生成数字减影图像，可降低数字减影血管造影过程中的血管伪影干扰。
附图说明
图1是根据本申请一些实施例所示的一种数字减影血管造影控制方法的示例性流程图；
图2是根据本申请一些实施例所示的一种数字减影血管造影控制系统的示例性硬件和/或软件的示意图；
图3是根据本申请一些实施例所示的实现数字减影血管造影控制方法的计算机设备的结构示意图。
具体实施方式
为了更好的理解上述技术方案，下面将结合说明书附图以及具体的实施方式对上述技术方案进行详细的说明。参考图1，该图是根据本申请一些实施例所示的数字减影血管造影控制方法的示例性流程图，该数字减影血管造影控制方法100主要包括如下步骤：
在步骤101，获取数字减影血管造影的基线图像。
基线图像是指在注射造影剂之前获取的无造影剂的图像，在数字减影血管成像中，基线图像用作背景图像，用于与注射造影剂后的图像进行减影操作，从而突出显示血管结构，具体实现时，可在进行数字减影血管成像之前通过数字减影血管造影设备采集得到基线图像，这里不再赘述。
在步骤102，采集数字减影血管造影的对比增强图像。
具体实现时，在对患者注射造影剂后即可开始连续的图像采集，得到对比增强图像，对比增强图像是通过注入造影剂或应用其他增强技术来突出血管结构的图像，需要说明的是，在采集对比增强图像过程中由于各种原因会产生血管伪影，导致血管图像模糊不清。
在步骤103，确定所述对比增强图像的频率域图，将该对比增强图像的频率域图转换为镜像频率域图，根据该对比增强图像的镜像频率域图确定血管伪影角度。
具体实现时，可通过傅里叶变换将对比增强图像转换为频率域图，这里不再赘述，在一些实施例中，对比增强图像的镜像频率域图是将对比增强图像的频率域图通过对数变换得到，本申请中对比增强图像的镜像频率域图可根据下述公式确定，即：
C＝B-1[log|H(u，v)+1|]
其中C表示对比增强图像的镜像频率域图，B-1[]表示博里叶逆变换，H(u，v)表示对比增强图像的频率域图，u和v分别表示频率的水平和垂直分量。
需要说明的是，本申请中根据该对比增强图像的镜像频率域图确定血管伪影角度具体可采用下述方式，即：
首先，确定该对比增强图像的镜像频率域图中血管伪影的边缘轮廓，其次，对于确定的血管伪影的边缘轮廓经过直线拟合即可确定血管伪影角度，具体实现时，对比增强图像中的血管伪影是图像产生模糊的部分，而将对比增强图像转换为对应的镜像频率域图后，在该对比增强图像的镜像频率域图中，图像中的模糊部分表现为一条直线的形式，即对比增强图像的镜像频率域图中的血管伪影是一条直线，本申请通过确定血管伪影边缘轮廓即可确定该直线，作为一个优选的实施例，可首先确定对比增强图像的镜像频率域图中血管伪影边缘轮廓值，根据确定的血管伪影边缘轮廓值即可最终确定血管边缘轮廓，其中，确定血管伪影的边缘轮廓值可根据下述公式确定：
其中PC(a，i)表示对比增强图像的镜像频率域图中血管伪影的边缘轮廓值，i表示检测方向，a表示对比增强图像的镜像频率域图的像素值，Wi(a)表示在i方向的权重因子，P0表示噪声能量阈值，α为预设的常数，Δφni(a)表示相位偏离函数，表示如该符号中的值为正，则结果为其本身，否则为零，Ani(a)表示在i方向上的幅度值。
需要说明的是，上述公式中对每个像素值和检测方向，需要计算权重因子wi(a)，而根据确定的权重因子wi(a)和相位偏离函数Δφni(a)可计算在i方向上的幅度值Ani(a)，这里不再赘述。
另外，需要说明的是，对于上述确定的血管伪影的边缘轮廓经过直线拟合确定血管伪影角度时可对提取的边缘轮廓进行直线拟合，以找到最佳拟合的直线，具体直线拟合可采用最小二乘法或其他直线拟合算法来实现，通过直线的斜率即可确定血管伪影角度，这里不再赘述。
在步骤104，根据所述血管伪影角度对该对比增强图像进行旋转，得到水平对比增强图像，由所述水平对比增强图像通过一阶差分得到对比增强梯度图像，根据所述对比增强梯度图像得到对比增强曲线图，由所述对比增强曲线图确定血管伪影长度。
具体实现时，可将该对比增强图像根据所述血管伪影角度进行顺时针旋转，使血管伪影的方向与水平对齐，从而在将该对比增强图像顺时针旋转所述血管伪影角度后即可得到所述水平对比增强图像。
在一些实施例中，通过在水平对比增强图像上应用一阶差分操作，计算图像中每个像素点的梯度，可得到对比增强梯度图像，其中高的梯度值表示边缘或变化较大的区域。
另外，在一些实施例中，根据所述对比增强梯度图像进行自相关运算得到对比增强曲线图可采用下述方式：
首先，对所述对比增强梯度图像遍历每一行i计算自相关值Ri(n)；
其次，计算出来的自相关值Ri(n)进行队列求和再平均，得到对比增强曲线图，其中自相关值Ri(n)可根据下述公式确定：
其中Ri(n)表示在对比增强梯度图像中的第i行的自相关值，m表示所述对比增强梯度图像的列索引，n表示列索引的偏移量，d′(i，m)表示所述对比增强梯度图像中第i行第m列的像素值，d′(i，m-n)表示所述对比增强梯度图像中第i行第m-n列的像素值；
具体实现时，根据所述对比增强曲线图中的左右两个对称的负相关峰值之间的距离，将所述距离除以二所得到的结果即为血管伪影长度。
需要说明的，在对比增强曲线图中，也可以通过分析对比增强曲线的特征确定血管伪影的长度，例如可以选择一定的阈值或曲线特征(如梯度突变)来确定血管伪影的起点和终点，从而计算血管伪影的长度，这里不再赘述。
在步骤105，根据所述血管伪影角度和所述血管伪影长度确定血管伪影模糊核。
具体实现时，上述确定了血管伪影角度和血管伪影长度，通过所述血管伪影角度与所述血管伪影长度即可以确定所述血管伪影模糊核，例如，可根据血管伪影的长度确定模糊核的长度，并根据血管伪影的角度确定模糊核的方向，使用一个矩阵或数组来表示模糊核，将模糊核的中心位置设置为血管伪影长度的中心，并根据血管伪影的角度确定模糊核的方向，在血管伪影模糊核中，可以将中心位置的像素值设置为最大值，然后从中心逐渐降低像素值，以模拟血管伪影的逐渐模糊效果。
作为一个优选的实施例，也可以先根据所述血管伪影角度和所述血管伪影长度确定血管伪影点扩散函数，然后对确定的血管伪影点扩散函数进行逆滤波即可得到血管伪影模糊核，即将血管伪影的点扩散函数进行傅里叶变换，得到频域上的表示，然后将血管伪影的原始图像进行傅里叶变换，得到频域上的表示，通过将血管伪影原始图像的频域表示与点扩散函数的频域表示相除，即进行逆滤波操作得到血管伪影模糊核，作为一个优选的实施例，血管伪影的点扩散函数可根据下述公式确定：
其中M(x，y)表示血管伪影点扩散函数，(x，y)表示对比增强图像中像素点的位置，X表示对比增强图像中像素点的水平坐标位置，y表示对比增强图像中像素点的垂直坐标位置，表示血管伪影角度，L表示血管伪影长度。
在步骤106，将所述对比增强图像进行对称填充，得到该对比增强图像的扩展图像，根据该对比增强图像的扩展图像和所述血管伪影模糊核对对比增强图像进行校正，将校正后的对比增强图像与所述基线图像进行数字减法运算生成数字减影图像。
具体实现时，对称填充即将对比增强图像的边缘进行镜像对称，对所述对比增强图像可使用边缘检测算法提取边界信息，所述边界信息为边缘像素的集合，使用所述边界信息在所述对比增强图像边界的对称位置上进行填充操作，具体可以将所述边界信息进行镜像到对称位置上实现填充，填充后得到所述对比增强图像的四边周围有四个镜像图像，再将上下两个图像按照相同步骤进行左右镜像操作，最后得到八个镜像图像，将所述对比增强图像与镜像图像进行合成，具体操作通过在对应的边界位置上将对比增强图像与镜像图像像素进行叠加或融合，得到对比增强图像的扩展图像，这里不再赘述。
在一些实施例，根据该对比增强图像的扩展图像和所述血管伪影模糊核对对比增强图像进行校正具体可采用下述方式，即：
将该对比增强图像的扩展图像与血管伪影模糊核进行卷积操作；
根据卷积结果，从卷积结果中提取校正后的对比增强图像，例如可截取卷积结果中心部分的图像或通过阈值提取感兴趣的区域，这里不做具体限定。
另外，在一些实施例中，将校正的对比增强图像与所述基线图像进行数字减法运算生成数字减影图像具体采用下述方式，即：
将校正后的对比增强图像与基线图像进行空间上的配准，以确保两组图像对应的血管区域位置一致；
将配准后的对比增强图像和基线图像进行对齐，使得两组图像的血管结构在空间上完全重合；
将配准及对齐的基线图像从对比增强图像中减去即可生成数字减影图像。
另外，在一些实施例中，参考图2，该图是根据本申请一些实施例所示的数字减影血管造影控制系统的示例性硬件和/或软件的示意图，本实施例基于数字减影血管造影控制系统200可包括：基线图像获取模块201、对比增强图像采集模块202、血管伪影角度确定模块203、血管伪影长度确定模块204、血管伪影模糊核确定模块205和数字减影图像生成模块206，分别说明如下：
基线图像获取模块201，本申请中基线图像获取模块201主要用于获取数字减影血管造影的基线图像；
对比增强图像采集模块202，本申请中对比增强图像采集模块202主要用于采集数字减影血管造影的对比增强图像；
血管伪影角度确定模块203，本申请中血管伪影角度确定模块203主要用于确定所述对比增强图像的频率域图，将该对比增强图像的频率域图转换为镜像频率域图，根据该对比增强图像的镜像频率域图确定血管伪影角度；
血管伪影长度确定模块204，本申请中血管伪影长度确定模块204主要用于根据所述血管伪影角度对该对比增强图像进行旋转，得到水平对比增强图像，由所述水平对比增强图像通过一阶差分得到对比增强梯度图像，根据所述对比增强梯度图像得到对比增强曲线图，由所述对比增强曲线图确定血管伪影长度；
血管伪影模糊核确定模块205，本申请中血管伪影模糊核确定模块205主要用于根据所述血管伪影角度和所述血管伪影长度确定血管伪影模糊核；
数字减影图像生成模块206，本申请中数字减影图像生成模块206主要用于将所述对比增强图像进行对称填充，得到该对比增强图像的扩展图像，根据该对比增强图像的扩展图像和所述血管伪影模糊核对对比增强图像进行校正，将校正后的对比增强图像与所述基线图像进行数字减法运算生成数字减影图像。
在一些实施例中，本申请还提供一种计算机设备，所述计算机设备包括存储器和处理器，所述存储器存储有代码，所述处理器被配置为获取所述代码，并执行上述的数字减影血管造影控制方法。
在一些实施例中，参考图3，该图是根据本申请一些实施例所示的实现数字减影血管造影控制方法的计算机设备的结构示意图。上述实施例中的数字减影血管造影控制方法可以通过图3所示的计算机设备来实现，该计算机设备300包括至少一个处理器301、通信总线302、存储器303以及至少一个通信接口304。
处理器301可以是一个通用中央处理器(central processing unit，CPU)、特定应用集成电路(application-specific integrated circuit，ASIC)或一个或多个用于控制本申请中数字减影血管造影控制方法的执行。
通信总线302可包括一通路，在上述组件之间传送信息。
存储器303可以是只读存储器(read-only memory，ROM)或可存储静态信息和指令的其它类型的静态存储设备，随机存取存储器(random access memory，RAM)或者可存储信息和指令的其它类型的动态存储设备，也可以是电可擦可编程只读存储器(electricallyerasable programmable read-only memory，EEPROM)、只读光盘(compact disc read-only Memory，CD-ROM)或其它光盘存储、光碟存储(包括压缩光碟、激光碟、光碟、数字通用光碟、蓝光光碟等)、磁盘或者其它磁存储设备、或者能够用于携带或存储具有指令或数据结构形式的期望的程序代码并能够由计算机存取的任何其它介质，但不限于此。存储器303可以是独立存在，通过通信总线302与处理器301相连接。存储器303也可以和处理器301集成在一起。
其中，存储器303用于存储执行本申请方案的程序代码，并由处理器301来控制执行。处理器301用于执行存储器303中存储的程序代码。程序代码中可以包括一个或多个软件模块。上述实施例中生成数字减影图像可以通过处理器301以及存储器303中的程序代码中的一个或多个软件模块实现。
通信接口304，使用任何收发器一类的装置，用于与其它设备或通信网络通信，如以太网，无线接入网(radio access network，RAN)，无线局域网(wireless local areanetworks，WLAN)等。
在具体实现中，作为一种实施例，计算机设备可以包括多个处理器，这些处理器中的每一个可以是一个单核(single-CPU)处理器，也可以是一个多核(multi-CPU)处理器。这里的处理器可以指一个或多个设备、电路、和/或用于处理数据(例如计算机程序指令)的处理核。
上述的计算机设备可以是一个通用计算机设备或者是一个专用计算机设备。在具体实现中，计算机设备可以是台式机、便携式电脑、网络服务器、掌上电脑(personaldigital assistant，PDA)、移动手机、平板电脑、无线终端设备、通信设备或者嵌入式设备。本申请实施例不限定计算机设备的类型。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
例如，在一些实施例中，本申请还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现上述的数字减影血管造影控制方法。
本发明是根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本发明的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。
Claims (10)
1.一种数字减影血管造影控制方法，其特征在于，包括如下步骤：
获取数字减影血管造影的基线图像；
采集数字减影血管造影的对比增强图像；
确定所述对比增强图像的频率域图，将该对比增强图像的频率域图转换为镜像频率域图，根据该对比增强图像的镜像频率域图确定血管伪影角度；
根据所述血管伪影角度对该对比增强图像进行旋转，得到水平对比增强图像，由所述水平对比增强图像通过一阶差分得到对比增强梯度图像，根据所述对比增强梯度图像得到对比增强曲线图，由所述对比增强曲线图确定血管伪影长度；
根据所述血管伪影角度和所述血管伪影长度确定血管伪影模糊核；
将所述对比增强图像进行对称填充，得到该对比增强图像的扩展图像，根据该对比增强图像的扩展图像和所述血管伪影模糊核对对比增强图像进行校正，将校正后的对比增强图像与所述基线图像进行数字减法运算生成数字减影图像。
2.如权利要求1所述的方法，其特征在于，根据该对比增强图像的镜像频率域图确定血管伪影角度具体包括：
确定该对比增强图像的镜像频率域图中血管伪影的边缘轮廓；
对于确定的血管伪影的边缘轮廓经过直线拟合确定血管伪影角度。
3.如权利要求2所述的方法，其特征在于，采用最小二乘法进行直线拟合。
4.如权利要求2所述的方法，其特征在于，血管伪影的边缘轮廓值根据下述公式确定：
其中PC(a，i)表示对比增强图像的镜像频率域图中血管伪影的边缘轮廓值，a表示对比增强图像的镜像频率域图的像素值，i表示不同的检测方向，Wi(a)表示在i方向的权重因子，P0表示噪声能量阈值，α为预设的常数，Δφni(a)表示相位偏离函数，表示如该符号中的值为正，则结果为其本身，否则为零，Ani(a)表示在i方向上的幅度值。
5.如权利要求1所述的方法，其特征在于，根据所述血管伪影角度和所述血管伪影长度确定血管伪影模糊核具体包括：
根据所述血管伪影角度和所述血管伪影长度确定血管伪影点扩散函数；
对确定的血管伪影点扩散函数进行逆滤波得到血管伪影模糊核。
6.如权利要求1所述的方法，其特征在于，根据该对比增强图像的扩展图像和所述血管伪影模糊核对对比增强图像进行校正具体包括：
将该对比增强图像的扩展图像与血管伪影模糊核进行卷积操作；
根据卷积结果，从卷积结果中提取校正后的对比增强图像。
7.如权利要求1所述的方法，其特征在于，将校正后的对比增强图像与所述基线图像进行数字减法运算生成数字减影图像具体包括：
将校正后的对比增强图像与基线图像进行空间上的配准，以确保两组图像对应的血管区域位置一致；
将配准后的对比增强图像和基线图像进行对齐，使得两组图像的血管结构在空间上完全重合；
将配准及对齐的基线图像从对比增强图像中减去，得到生成数字减影图像。
8.一种数字减影血管造影控制系统，其特征在于，包括有：
基线图像获取模块，用于获取数字减影血管造影的基线图像；
对比增强图像采集模块，用于采集数字减影血管造影的对比增强图像；
血管伪影角度确定模块，用于确定所述对比增强图像的频率域图，将该对比增强图像的频率域图转换为镜像频率域图，根据该对比增强图像的镜像频率域图确定血管伪影角度；
血管伪影长度确定模块，用于根据所述血管伪影角度对该对比增强图像进行旋转，得到水平对比增强图像，由所述水平对比增强图像通过一阶差分得到对比增强梯度图像，根据所述对比增强梯度图像得到对比增强曲线图，由所述对比增强曲线图确定血管伪影长度；
血管伪影模糊核确定模块，用于根据所述血管伪影角度和所述血管伪影长度确定血管伪影模糊核；
数字减影图像生成模块，用于将所述对比增强图像进行对称填充，得到该对比增强图像的扩展图像，根据该对比增强图像的扩展图像和所述血管伪影模糊核对对比增强图像进行校正，将校正后的对比增强图像与所述基线图像进行数字减法运算生成数字减影图像。
9.一种计算机设备，其特征在于，所述计算机设备包括存储器和处理器，所述存储器存储有代码，所述处理器被配置为获取所述代码，并执行如权利要求1至7任一项所述的数字减影血管造影控制方法。
10.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至7任一项所述的数字减影血管造影控制方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310785144.1A CN116797683A (zh) | 2023-06-29 | 2023-06-29 | 一种数字减影血管造影控制系统及控制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310785144.1A CN116797683A (zh) | 2023-06-29 | 2023-06-29 | 一种数字减影血管造影控制系统及控制方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116797683A true CN116797683A (zh) | 2023-09-22 |
Family
ID=88039935
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310785144.1A Withdrawn CN116797683A (zh) | 2023-06-29 | 2023-06-29 | 一种数字减影血管造影控制系统及控制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116797683A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117726744A (zh) * | 2023-12-21 | 2024-03-19 | 强联智创（北京）科技有限公司 | 生成三维数字减影血管造影图像的方法、设备和存储介质 |
-
2023
- 2023-06-29 CN CN202310785144.1A patent/CN116797683A/zh not_active Withdrawn
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117726744A (zh) * | 2023-12-21 | 2024-03-19 | 强联智创（北京）科技有限公司 | 生成三维数字减影血管造影图像的方法、设备和存储介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6120843B2 (ja) | 時間フィルタリングおよび空間フィルタリングによる一連の蛍光透視画像における雑音を低減させるための方法 | |
| US20190206070A1 (en) | Image registration method | |
| EP2881039A1 (en) | X-ray computer tomography image pick-up device and image reconstruction method | |
| CN102999884B (zh) | 图像处理设备和方法 | |
| KR20130038794A (ko) | 디지털 엑스레이 프레임들 시리즈에서의 잡음 감소의 방법 | |
| JP2007529071A (ja) | 画像におけるエッジの検出 | |
| CN116797683A (zh) | 一种数字减影血管造影控制系统及控制方法 | |
| US20110069907A1 (en) | Image reconstruction method | |
| CN111402150B (zh) | 一种ct图像金属伪影去除方法及装置 | |
| US20170004601A1 (en) | Method for reducing noise in sequences of fluoroscopic images | |
| CN107680057A (zh) | 超声图像增强的方法及装置 | |
| CN104463815B (zh) | Dsa图像生成方法及系统 | |
| Huang et al. | Reference-free learning-based similarity metric for motion compensation in cone-beam CT | |
| CN113793272A (zh) | 图像降噪方法及装置、存储介质、终端 | |
| CN105261051B (zh) | 一种获取图像掩码的方法及装置 | |
| CN104318595A (zh) | Ct图像的运动向量场的计算方法及装置 | |
| CN116703994B (zh) | 医学图像配准的方法、计算设备和计算机可读存储介质 | |
| CN111161330B (zh) | 非刚性图像配准方法、装置、系统、电子设备、存储介质 | |
| JP2021186430A (ja) | 画像処理装置及び画像処理方法。 | |
| CN111311531A (zh) | 图像增强方法、装置、控制台设备及医学成像系统 | |
| CN107886478B (zh) | 一种ct图像重建方法及系统、终端及可读存储介质 | |
| CN105574904B (zh) | 一种图像重建方法、装置及设备 | |
| JP2013017781A (ja) | 医用画像抽出装置および医用画像抽出プログラム | |
| CN109745073B (zh) | 弹性成像位移的二维匹配方法及设备 | |
| WO2013078451A1 (en) | Methods of determining local spectrum at a pixel using a rotationally invariant s-transform (rist) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20230922 |
|
| WW01 | Invention patent application withdrawn after publication |