CN117495893B - 一种基于主动轮廓模型的颅骨剥离方法 - Google Patents
一种基于主动轮廓模型的颅骨剥离方法 Download PDFInfo
- Publication number
- CN117495893B CN117495893B CN202311788250.1A CN202311788250A CN117495893B CN 117495893 B CN117495893 B CN 117495893B CN 202311788250 A CN202311788250 A CN 202311788250A CN 117495893 B CN117495893 B CN 117495893B
- Authority
- CN
- China
- Prior art keywords
- image
- skull
- brain
- mask
- pyramid
- 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.)
- Active
Links
- 210000003625 skull Anatomy 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 36
- 210000004556 brain Anatomy 0.000 claims abstract description 78
- 238000003709 image segmentation Methods 0.000 claims abstract description 8
- 238000007781 pre-processing Methods 0.000 claims abstract description 5
- 230000011218 segmentation Effects 0.000 claims description 37
- 230000006870 function Effects 0.000 claims description 30
- 230000009466 transformation Effects 0.000 claims description 16
- 238000012545 processing Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 4
- 238000012952 Resampling Methods 0.000 claims description 3
- 238000012804 iterative process Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims 1
- 210000001519 tissue Anatomy 0.000 abstract description 4
- 230000004927 fusion Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 7
- 238000002595 magnetic resonance imaging Methods 0.000 description 6
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 4
- 210000005013 brain tissue Anatomy 0.000 description 4
- 238000002224 dissection Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 208000026062 Tissue disease Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/149—Segmentation; Edge detection involving deformable models, e.g. active contour models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/12—Edge-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
- G06T7/337—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods involving reference images or patches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20112—Image segmentation details
- G06T2207/20116—Active contour; Active surface; Snakes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30016—Brain
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种基于主动轮廓模型的颅骨剥离方法，包括步骤1：获取原始三维医学图像数据及脑模板数据；步骤2：对三维医学图像数据进行预处理；步骤3：将脑模板数据配准到原始医学图像数据空间；步骤4：基于主动轮廓模型进行金字塔图像分割；步骤5：实现颅骨剥离。本发明可以快速实现颅骨剥离，且保证其精确度，以帮助医生更好地分析大脑其他组织结构。且此方法适用性广，不局限于MR图像，同样也适用于CT图像。
Description
技术领域
本发明属于医学图像处理技术领域，涉及到一种基于主动轮廓模型的颅骨剥离方法。
背景技术
核磁共振成像（MRI）具有空间分辨率高、软组织对比度高、成像参数多、信息量大、非侵入、低损伤等优点，在各级医疗机构、科研院所以及企业研发机构得到了广泛应用，并成为了脑组织疾病诊断与治疗、脑神经科学研究与实验领域常用研究手段之一。
由于脑组织具有体积较小，组织间边缘模糊，所以直接对脑部的某个区域研究很困难。脑MRI通常能够详细的覆盖整个大脑的神经解剖学信息，在此基础上用核磁共振图像分割的目的是提取医生需要研究的特定区域，医生可以通过分割出的区域对组织、器官进行仔细的观察和诊断，并据此制定有效的治疗方案，从而有效的提高疾病诊断的正确率。
而在分析颅脑中的特定组织或结构时，颅骨的存在会对其识别产生一定的影响，而且颅骨的存在也会对增加后续的各种图像处理的运算量。因此为了后续的脑组织分割及各种图像处理更加准确、快速，有必要对颅骨进行剥离。
发明内容
为解决上述问题，本发明公开了一种基于主动轮廓模型的颅骨剥离方法。
具体方案如下：
一种基于主动轮廓模型的颅骨剥离方法，包括：
步骤1：获取原始三维医学图像数据及脑模板数据；
步骤2：对三维医学图像数据进行预处理；
步骤3：将脑模板数据配准到原始医学图像数据空间；
步骤4：基于主动轮廓模型进行金字塔图像分割；
步骤5：实现颅骨剥离。
进一步地，所述三维医学图像数据为颅脑MRI序列，为待分割的原始图像OriImage序列；所述脑模板数据包括脑模板图像TemplateImage序列及对应的掩膜标签LabelMask和模板的头颅掩膜区域AreaMask。脑模板图像TemplateImage序列为一标准常规颅脑MRI序列，掩膜标签LabelMask为对应的颅骨剥离后的大脑区域掩膜Mask数据，头颅掩膜区域AreaMask为常规颅脑序列对应的脑部区域掩膜Mask。
进一步地，所述步骤2中图像预处理，具体方法为：
步骤2.1：将步骤1中所涉及到的图像数据进行重采样，重采样后的数据为各向同性，且体素大小一致；
步骤2.2：使用头颅掩膜区域AreaMask数据对脑模板图像TemplateImage进行预分割。头颅掩膜区域AreaMask为脑部区域的掩膜Mask，是二值图像，脑部区域为1，背景为0，将头颅掩膜区域AreaMask和脑模板图像TemplateImage进行与运算，在脑模板图像TemplateImage中保留脑部区域，其他区域置0，以减少后续图像处理所需的运算量；
步骤2.3：将原始图像OriImage和脑模板图像TemplateImage的灰度值映射到0-255区间，统一压缩灰度区间，有利于减少图像内存占用，且能减少运算量。
进一步地，所述步骤3中所述脑模板数据配准到原始医学图像数据空间，其目的为通过模板匹配，获得原始图像大脑区域的Mask图像，具体方法为：
步骤3.1：将原始图像OriImage设为基准的参考图像(fixedImage)，将脑模板图像TemplateImage设为移动图像(movingImage)，通过图像配准，使脑模板图像TemplateImage与原始图像OriImage进行对齐；
步骤3.2：获取步骤3.1中配准的变换模型；
步骤3.3：基于步骤3.2中的变换模型，将掩膜标签LabelMask变换到原始图像OriImage空间。
进一步地，所述步骤4中主动轮廓算法对大脑轮廓进行分割。主动轮廓算法基于曲线演化的思想，通过自动寻找图像中感兴趣区域的边界来实现分割，为提升分割算法的效率和准确性，可采用金字塔图像分割思想，从粗糙到细致地处理图像，其具体方法为：
步骤4.1：构建图像金字塔；通过对原始图像进行多次降采样，生成一系列不同分辨率的图像，每一层图像都是上一层图像的1/4大小，直到达到所需的分辨率或金字塔的最底层；
步骤4.2：分割金字塔底层；构建初始轮廓，使用主动轮廓算法对最底层金字塔做初始分割；
步骤4.3：金字塔上采样；将步骤4.2中分割结果上采样到较高分辨率。这可以通过插值方法，如双线性插值或最近邻插值等来完成；
步骤4.4：分割结果融合；将上采样后的分割结果与高分辨率的图像进行融合，得到更精确的分割结果；
步骤4.5：迭代；重复进行步骤4.2至步骤4.4，将融合后的分割结果作为下一层金字塔的输入，直到达到金字塔的顶层；
步骤4.6：输出最终分割结果；在金字塔的顶层得到最终的分割结果。
进一步地，所述步骤4中主动轮廓算法对大脑轮廓进行分割。主动轮廓算法具体实现方法为：
步骤4-1：初始化轮廓：以步骤3中变换后的掩膜标签LabelMask作为初始化轮廓；
步骤4-2：定义能量函数：能量函数一般包括内部能量和外部能量，以曲率能量项作为内部能量函数；外部能量函数为基于图像梯度的能量项；
步骤4-3：优化能量函数：通过迭代的方式优化能量函数，调整轮廓的形状；迭代过程中，根据能量函数的最小值来更新轮廓的位置；
步骤4-4：定义停止条件：定义迭代停止的条件，如迭代次数达到阈值、轮廓收敛或能量函数的变化小于某个阈值等。
进一步地，所述步骤5中实现颅骨剥离，具体为基于步骤4中的轮廓分割结果提取掩膜Mask，并生成颅骨剥离图像，实现颅骨剥离。
本发明的有益效果在于：应用本发明的技术方案，可以快速实现颅骨剥离，且保证其精确度，以帮助医生更好地分析大脑其他组织结构。且此方法适用性广，不局限于MR图像，同样也适用于CT图像。
附图说明
图1为本发明的流程示意图。
图2为本发明中待分割的原始图像OriImage序列示意图。
图3为本发明中脑模板图像TemplateImage示意图。
图4 为本发明中脑模板图像TemplateImage与其对应的掩膜标签LabelMask融合效果示意图。
图5 为本发明中脑模板图像TemplateImage与其对应的头颅掩膜区域AreaMask融合效果示意图。
图6 为本发明中分割后的大脑轮廓示意图。
图7为本发明中基于轮廓生成的mask示意图。
图8 为本发明中生成的mask与原始图像融合显示示意图。
图9 为本发明中颅骨剥离后的大脑组织示意图。
图10为本发明中基于颅脑CT数据颅骨剥离的融合显示效果图。
具体实施方式
下面结合附图和具体实施方式，进一步阐明本发明，应理解下述具体实施方式仅用于说明本发明而不用于限制本发明的范围。
本发明提供了一种基于主动轮廓模型的颅骨剥离方法，具体流程如图1所示，包括：
步骤1：获取原始三维医学图像数据及脑模板数据。所述三维医学图像数据为颅脑MRI序列，为待分割的原始图像OriImage序列，如图2所示；所述脑模板数据包括脑模板图像TemplateImage序列（如图3所示）及对应的掩膜标签LabelMask（其与脑模板图像TemplateImage的融合显示图如图4所示）和模板的头颅掩膜区域AreaMask（其与脑模板图像TemplateImage的融合显示图如图5所示）。脑模板图像序列为一标准常规颅脑MRI序列，掩膜标签LabelMask为对应的颅骨剥离后的大脑区域掩膜Mask数据，头颅掩膜区域AreaMask为常规颅脑序列对应的脑部区域掩膜Mask。
步骤2：对三维医学图像数据进行预处理，具体方法为：
步骤2.1：将步骤1中所涉及到的图像数据进行重采样，重采样后的数据为各向同性，且体素大小一致，体素大小（长、宽、高）可设为（1,1,1），即重采样的3D图像数据分辨率为(1,1,1)；
步骤2.2：使用头颅掩膜区域AreaMask数据对脑模板图像TemplateImage进行预分割。头颅掩膜区域AreaMask为脑部区域的掩膜Mask，是二值图像，脑部区域为1，背景为0，将头颅掩膜区域AreaMask和脑模板图像TemplateImage进行与运算，在脑模板图像TemplateImage中保留脑部区域，其他区域置0，以减少后续图像处理所需的运算量；
步骤2.3：将原始图像OriImage和脑模板图像TemplateImage的灰度值映射到0-255区间，统一压缩灰度区间，有利于减少图像内存占用，且能减少运算量。具体计算方法为：
获取图像的最大灰度值max和最小灰度值min；
计算缩放因子dFactor，
；
；
其中，pixel_val为当前像素的灰度值，nPixelVal计算得到的新的像素灰度值，(int)为取整操作。
步骤3：将脑模板数据配准到原始医学图像数据空间，其目的为通过模板匹配，获得原始图像大脑区域的Mask图像，具体方法为：
步骤3.1：将原始图像OriImage设为基准的参考图像(fixedImage)，将脑模板图像TemplateImage设为移动图像(movingImage)，通过图像配准，使T脑模板图像emplateImage与原始图像OriImage进行对齐。可以采用金字塔图像配准方式进行图像配准。 其思想为使用一种粗糙的尺度对少量的图像像素进行处理,然后在下一层使用一种精确的尺度, 并用上一层的结果对其参数进行初始化. 迭代该过程, 直到达到最精确的尺度。这种由粗到细, 在大尺度上看整体, 在小尺度上看细节的方法能够极大程度地提高配准成功率；
在配准中，采用Mattes Mutual Information（Mattes互信息，MMI）来度量图像相似性，其求解公式为：
其中，A和B是待比较的两幅图像，N是采样点的数量，表示第i个采样点在两幅图像中的像素值，/>表示/>在两幅图像中同时出现的概率，/>表示/>在B中出现的概率；
采用VersorRigid3DTransformOptimizer作为优化器，其作用为作用是通过最小化图像之间的差异来寻找最佳的变换参数，以将源图像和目标图像对准。它使用梯度下降法来搜索最优参数，并通过反向传播来计算梯度。设置每次迭代中最大步长为0.05，设置最小步长为0.005，设置最大迭代步数为500；
步骤3.2：获取步骤1中配准的变换模型，记为transformModel；
步骤3.3：基于步骤3.2中的变换模型，将掩膜标签LabelMask变换到原始图像OriImage空间，计算方式为：
其中为/>基于变换模型得到的新的mask；
步骤3.4：此外，为提高配准的精度和效率，可通过先采用刚性配准后仿射配准或非刚性配准的策略，刚性变换（如平移、旋转和缩放）是一种比较简单的变换模型，它可以处理图像之间的整体几何差异。通过先进行刚性配准，可以更好地对齐图像的整体结构，使得后续的配准任务更加容易和准确。仿射变换是一种比刚性变换更灵活的变换模型，它除了平移、旋转和缩放外，还可以进行剪切和错切等变换操作。通过在刚性配准后进行仿射配准，可以更好地处理图像之间的局部形变和细微差异，提高配准的精度。通过先进行刚性配准，再进行仿射配准，可以逐步处理图像的整体结构和局部形变，提高配准的效果。
步骤4：使用主动轮廓算法对大脑轮廓进行分割。进一步地，所述步骤4中主动轮廓算法对大脑轮廓进行分割。主动轮廓算法基于曲线演化的思想，通过自动寻找图像中感兴趣区域的边界来实现分割，为提升分割算法的效率和准确性，可采用金字塔图像分割思想，从粗糙到细致地处理图像，其具体方法为：
步骤4.1：构建图像金字塔。通过对原始图像进行多次降采样，生成一系列不同分辨率的图像，每一层图像都是上一层图像的1/4大小，直到达到所需的分辨率或金字塔的最底层；
步骤4.2：分割金字塔底层。构建初始轮廓，使用主动轮廓算法对最底层金字塔做初始分割；
步骤4.3：金字塔上采样。将步骤4.2中分割结果上采样到较高分辨率。这可以通过插值方法，如双线性插值或最近邻插值等来完成；
步骤4.4：分割结果融合：将上采样后的分割结果与高分辨率的图像进行融合，得到更精确的分割结果；
步骤4.5：迭代：重复进行步骤4.2至步骤4.4，将融合后的分割结果作为下一层金字塔的输入，直到达到金字塔的顶层；
步骤4.6：输出：最终分割结果：在金字塔的顶层得到最终的分割结果。
在本实施例中，步骤4中主动轮廓算法对大脑轮廓进行分割。主动轮廓算法是一种基于能量最小化的图像分割方法。它的基本思想是将一个初始轮廓放置在图像中，并通过调整轮廓的形状和位置来逐步逼近目标边界，从而实现图像分割。主动轮廓算法为一种算法思想，有多种具体算法形态，此处以Chan-Vese(CV)水平集分割算法为例，水平集算法的优点在于它能够自适应地调整轮廓的形状和拓扑结构，从而在处理具有复杂形状和不规则边界的图像时具有较高的精度和鲁棒性。其具体实现步骤为：
步骤4-1：初始化轮廓：以步骤3中变换后的掩膜标签LabelMask的轮廓作为初始化轮廓，从掩膜标签LabelMask中提取轮廓可使用边缘提取算法（如Canny边缘检测算法）；
步骤4-2：定义能量函数：能量函数一般包括内部能量和外部能量，以曲率能量项作为内部能量函数；外部能量函数为基于图像梯度的能量项；
其能量函数为：
其中，E表示能量函数，c1和c2是两个不同的灰度值，φ是水平集函数，μ、ν、λ1和λ2是权重参数，I是输入图像，R1和R2是分割区域。Length(φ)表示轮廓的长度，可以通过计算水平集函数的梯度模长来得到；
表示分割区域的面积，可以通过计算水平集函数大于0和小于0的像素数量来得到；
RegionTerm(I, R)表示基于图像亮度的能量项，可以通过计算分割区域内的灰度均值和图像灰度值的差平方来得到。该项用于推动轮廓向图像中灰度变化较大的区域移动；
步骤4-3：优化能量函数：通过迭代的方式优化能量函数，调整轮廓的形状。迭代过程中，根据能量函数的最小值来更新轮廓的位置；
步骤4-4：定义停止条件：定义迭代停止的条件，如迭代次数达到阈值、轮廓收敛或能量函数的变化小于某个阈值等。一种可行的参数设置方案为：设置最大迭代次数为400；设置迭代停止的均方根误差阈值为0.02；设置曲率项的系数为1；设置传播项的系数为1；设置平流项的系数为1。
步骤5，输出分割结果，实现颅骨剥离。具体方法为：
步骤5.1：对步骤4中最终分割得到的轮廓进行上采样到分辨率与原始图像OriImage一致，上采样后的轮廓与原始图像OriImage融合显示效果如图6所示；
步骤5.2：生成与原始图像OriImage尺寸一致的mask图像数据，并初始化为0；
步骤5.3：以轮廓的重心为种子点，通过区域生长方法获得颅骨剥离后的大脑mask图像数据，轮廓内部像素灰度值置为1，生成的mask图像如图7所示，mask图像与原始图像OriImage融合显示效果如图8所示；
步骤5.4：将mask与原始图像OriImage进行与操作，只保留mask区域对应的原始图像OriImage的原始灰度值，其他区域像素值均为0，完成颅骨剥离，颅骨剥离后的图像如图9所示。
本专利所实现方法不局限于MR图像，同样也使用与颅脑CT图像，基于颅脑CT图像的颅骨剥离效果如图10所示。
Claims (6)
1.一种基于主动轮廓模型的颅骨剥离方法，其特征在于，包括：
步骤1：获取原始三维医学图像数据及脑模板数据；
步骤2：对三维医学图像数据进行预处理；
步骤3：将脑模板数据配准到原始医学图像数据空间；
步骤4：基于主动轮廓模型进行金字塔图像分割；
步骤5：实现颅骨剥离；
所述步骤4中的主动轮廓算法基于曲线演化的思想，通过自动寻找图像中感兴趣区域的边界来实现分割，为提升分割算法的效率和准确性，采用金字塔图像分割思想，从粗糙到细致地处理图像，其具体方法为：
步骤4.1：构建图像金字塔；通过对原始图像进行多次降采样，生成一系列不同分辨率的图像，每一层图像都是上一层图像的1/4大小，直到达到所需的分辨率或金字塔的最底层；
步骤4.2：分割金字塔底层；构建初始轮廓，使用主动轮廓算法对最底层金字塔做初始分割；
步骤4.3：金字塔上采样；将步骤4.2中分割结果上采样到高分辨率；通过插值方法来完成；
步骤4.4：分割结果融合；将上采样后的分割结果与高分辨率的图像进行融合，得到更精确的分割结果；
步骤4.5：迭代；重复进行步骤4.2至步骤4.4，将融合后的分割结果作为下一层金字塔的输入，直到达到金字塔的顶层；
步骤4.6：输出最终分割结果；在金字塔的顶层得到最终的分割结果。
2.根据权利要求1所述的一种基于主动轮廓模型的颅骨剥离方法，其特征在于，所述三维医学图像数据为颅脑MRI序列，为待分割的原始图像序列；所述脑模板数据包括脑模板图像序列及对应的掩膜标签和模板的头颅掩膜区域；脑模板图像序列为一标准常规颅脑MRI序列，掩膜标签为对应的颅骨剥离后的大脑区域掩膜数据，头颅掩膜区域为常规颅脑序列对应的脑部区域掩膜。
3.根据权利要求2所述的一种基于主动轮廓模型的颅骨剥离方法，其特征在于，所述步骤2的具体方法为：
步骤2.1：将步骤1中所涉及到的图像数据进行重采样，重采样后的数据为各向同性，且体素大小一致；
步骤2.2：使用头颅掩膜区域数据对脑模板图像进行预分割；头颅掩膜区域为脑部区域的掩膜，是二值图像，脑部区域为1，背景为0，将头颅掩膜区域和脑模板图像进行与运算，在脑模板图像中保留脑部区域，其他区域置0，以减少后续图像处理所需的运算量；
步骤2.3：将原始图像和脑模板图像的灰度值映射到0-255区间，统一压缩灰度区间，减少图像内存占用，减少运算量。
4.根据权利要求3所述的一种基于主动轮廓模型的颅骨剥离方法，其特征在于，所述步骤3为通过模板匹配，获得原始图像大脑区域的Mask图像，具体方法为：
步骤3.1：将原始图像设为基准的参考图像，将脑模板图像设为移动图像，通过图像配准，使脑模板图像与原始图像进行对齐；
步骤3.2：获取步骤3.1中配准的变换模型；
步骤3.3：基于步骤3.2中的变换模型，将掩膜标签变换到原始图像空间。
5.根据权利要求4所述的一种基于主动轮廓模型的颅骨剥离方法，其特征在于，所述步骤4中主动轮廓算法的具体方法为：
步骤4-1：初始化轮廓；以步骤3中变换后的掩膜标签作为初始化轮廓；
步骤4-2：定义能量函数；能量函数包括内部能量和外部能量，以曲率能量项作为内部能量函数；外部能量函数为基于图像梯度的能量项；
步骤4-3：优化能量函数；通过迭代的方式优化能量函数，调整轮廓的形状；迭代过程中，根据能量函数的最小值来更新轮廓的位置；
步骤4-4：定义停止条件；定义迭代停止的条件，包括迭代次数达到阈值、轮廓收敛或能量函数的变化小于某个阈值。
6.根据权利要求5所述的一种基于主动轮廓模型的颅骨剥离方法，其特征在于，所述步骤5中实现颅骨剥离，具体为基于步骤4中的轮廓分割结果提取掩膜，并生成颅骨剥离图像，实现颅骨剥离。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311788250.1A CN117495893B (zh) | 2023-12-25 | 2023-12-25 | 一种基于主动轮廓模型的颅骨剥离方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311788250.1A CN117495893B (zh) | 2023-12-25 | 2023-12-25 | 一种基于主动轮廓模型的颅骨剥离方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117495893A CN117495893A (zh) | 2024-02-02 |
| CN117495893B true CN117495893B (zh) | 2024-03-19 |
Family
ID=89683246
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311788250.1A Active CN117495893B (zh) | 2023-12-25 | 2023-12-25 | 一种基于主动轮廓模型的颅骨剥离方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117495893B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107909589A (zh) * | 2017-11-01 | 2018-04-13 | 浙江工业大学 | 一种结合C‑V水平集和GrabCut算法的牙齿图像分割方法 |
| CN109509203A (zh) * | 2018-10-17 | 2019-03-22 | 哈尔滨理工大学 | 一种半自动脑部图像分割方法 |
| CN112419340A (zh) * | 2020-12-09 | 2021-02-26 | 沈阳先进医疗设备技术孵化中心有限公司 | 脑脊液分割模型的生成方法、应用方法及装置 |
| CN115100123A (zh) * | 2022-06-10 | 2022-09-23 | 北京理工大学 | 一种结合UNet和主动轮廓模型的大脑提取方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070047790A1 (en) * | 2005-08-30 | 2007-03-01 | Agfa-Gevaert N.V. | Method of Segmenting Anatomic Entities in Digital Medical Images |
| WO2013082289A1 (en) * | 2011-11-30 | 2013-06-06 | Rush University Medical Center | System and methods for identification of implanted medical devices and/or detection of retained surgical foreign objects from medical images |
| CN109872328B (zh) * | 2019-01-25 | 2021-05-07 | 腾讯科技（深圳）有限公司 | 一种脑部图像分割方法、装置和存储介质 |
-
2023
- 2023-12-25 CN CN202311788250.1A patent/CN117495893B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107909589A (zh) * | 2017-11-01 | 2018-04-13 | 浙江工业大学 | 一种结合C‑V水平集和GrabCut算法的牙齿图像分割方法 |
| CN109509203A (zh) * | 2018-10-17 | 2019-03-22 | 哈尔滨理工大学 | 一种半自动脑部图像分割方法 |
| CN112419340A (zh) * | 2020-12-09 | 2021-02-26 | 沈阳先进医疗设备技术孵化中心有限公司 | 脑脊液分割模型的生成方法、应用方法及装置 |
| CN115100123A (zh) * | 2022-06-10 | 2022-09-23 | 北京理工大学 | 一种结合UNet和主动轮廓模型的大脑提取方法 |
Non-Patent Citations (3)
| Title |
|---|
| Skull Stripping Based on the Segmentation Models;Rasha Helmi Nayyef 等;《Journal of Engineering》;20231001;第29卷(第10期);74-89 * |
| 三维模型分割技术在数字牙齿医疗领域的应用研究;孟宏伟;《中国优秀硕士学位论文全文数据库 医药卫生科技辑》;20180815(第8期);E07-2 * |
| 基于多模态MRI图像的脑肿瘤分割方法;刘瑞;《中国优秀硕士学位论文全文数据库 医药卫生科技辑》;20180215(第2期);E072-1356 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117495893A (zh) | 2024-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP7493464B2 (ja) | ３ｄオブジェクトの正準ポーズの自動化判定、および深層学習を使った３ｄオブジェクトの重ね合わせ | |
| CN111798462B (zh) | 一种基于ct图像的鼻咽癌放疗靶区自动勾画方法 | |
| CN110310287B (zh) | 基于神经网络的危及器官自动勾画方法、设备和存储介质 | |
| CN106934821B (zh) | 一种基于icp算法和b样条的锥形束ct和ct图像配准方法 | |
| CN112967236B (zh) | 图像的配准方法、装置、计算机设备和存储介质 | |
| US7907768B2 (en) | Method and apparatus for probabilistic atlas based on shape modeling technique | |
| CN110599528A (zh) | 一种基于神经网络的无监督三维医学图像配准方法及系统 | |
| US8135189B2 (en) | System and method for organ segmentation using surface patch classification in 2D and 3D images | |
| CN108876769B (zh) | 一种左心耳ct图像分割方法 | |
| Maan et al. | Prostate MR image segmentation using 3D active appearance models | |
| Tavares | Analysis of biomedical images based on automated methods of image registration | |
| WO2008123969A1 (en) | Combined feature ensemble mutual information image registration | |
| Wirth et al. | Nonrigid mammogram registration using mutual information | |
| Abdel-Nasser et al. | Temporal mammogram image registration using optimized curvilinear coordinates | |
| CN113223004A (zh) | 一种基于深度学习的肝脏图像分割方法 | |
| CN111080658A (zh) | 基于可形变配准和dcnn的宫颈mri图像分割方法 | |
| US11334995B2 (en) | Hierarchical systems and methods for image segmentation | |
| CN115136189A (zh) | 基于图像处理的肿瘤的自动化检测 | |
| CN110782428A (zh) | 一种用于构建临床脑部ct图像roi模板的方法及系统 | |
| Jaffar et al. | Anisotropic diffusion based brain MRI segmentation and 3D reconstruction | |
| CN110570430B (zh) | 基于体配准的眼眶骨组织分割方法 | |
| Rexilius et al. | Multispectral brain tumor segmentation based on histogram model adaptation | |
| Rad et al. | Dental x-ray image segmentation and multiple feature extraction | |
| Papademetris et al. | Computing 3D non-rigid brain registration using extended robust point matching for composite multisubject fMRI analysis | |
| Alvén et al. | Überatlas: fast and robust registration for multi-atlas segmentation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |