CN109636872B - 基于滑动窗口差值和条带噪声检测的cbct环形伪影消除方法 - Google Patents
基于滑动窗口差值和条带噪声检测的cbct环形伪影消除方法 Download PDFInfo
- Publication number
- CN109636872B CN109636872B CN201811504142.6A CN201811504142A CN109636872B CN 109636872 B CN109636872 B CN 109636872B CN 201811504142 A CN201811504142 A CN 201811504142A CN 109636872 B CN109636872 B CN 109636872B
- Authority
- CN
- China
- Prior art keywords
- image
- noise
- cbct
- sliding window
- difference
- 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
Images
Classifications
-
- 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
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了基于滑动窗口差值和条带噪声检测的CBCT环形伪影消除方法,该方法包括对CBCT图像中存在的特征微小的环形伪影,通过滑动窗口差值提取边缘特征,通过水平方向的梯度图像检测条带噪声,再利用边缘特征来优化条带噪声,在CBCT图像的基础上通过对噪声区域的插值方法达到消除环形伪影的目的。本发明针对CBCT图像中存在的特征微小的环形伪影,通过检测边缘特征来优化条带噪声,在不对整体图像处理的前提下,最大限度地消除环形伪影。而且对于CBCT重建图像的每一层而言,其环形伪影具有一致性,即均处在相同的位置,所以检测一层数据的环形伪影即可对每一层数据进行环形伪影消除。
Description
技术领域
本发明属于医学图像数据处理领域,涉及CBCT环形伪影消除技术,具体是基于滑动窗口差值和条带噪声检测的CBCT环形伪影消除方法。
背景技术
CBCT(Cone-Beam Computed Tomography)以其低剂量、快捷等优势应用于辅助引导质子肿瘤治疗,因此重建CBCT图像的质量尤为重要。通常,由于采集设备的响应器件老化、三维重建算法(FDK是一种近似算法)不足、以及患者扫描过程中出现的运动等等都会造成重建后的CBCT图像存在不同程度的伪影。
针对CBCT图像中的环形伪影,现有的解决办法主要是将CBCT图像从直角坐标系转到极坐标下,这样环形伪影变为垂直直线噪声,使得环形伪影的检测变的更加简单。有的方法是直接在极坐标下进行L0范型滤波,平滑掉噪声的同时也使图像的细节部分缺失严重,而且,坐标转换需要重新进行插值运算,进一步的改变了原有图像。也有方法直接把CBCT的环形伪影当成噪声图像,利用去噪方法,如DDID(Dual-Domain Image Denoising)、利用双边滤波保留边缘细节,在去除噪声的基础上保留图像细节。最终CBCT图像环形伪影消除算法转化为对整体图像进行去噪保质的过程。而当环形伪影的特征不是特别明显,转换至极坐标系下将无法准确地检测出垂直直线;单纯地使用滤波方法直接去除所有噪声也是一个整体降质的过程。
发明内容
本发明的目的在于提供一种基于滑动窗口差值和条带噪声检测的CBCT环形伪影消除方法,其提供了一种条带检测和优化的方法。
本发明解决其技术问题采取的技术方案是:对CBCT图像中存在的特征微小的环形伪影,通过滑动窗口差值提取边缘特征,通过水平方向的梯度图像检测条带噪声,再利用边缘特征来优化条带噪声,在CBCT图像的基础上通过对噪声区域的插值方法达到消除环形伪影的目的。
本发明的目的可以通过以下技术方案实现:
基于滑动窗口差值和条带噪声检测的CBCT环形伪影消除方法,该方法包括下述步骤:
步骤一:对CBCT图像中存在的特征微小的环形伪影;
步骤二:通过滑动窗口差值提取边缘特征;
步骤三:通过水平方向的梯度图像检测条带噪声;
步骤四:再利用边缘特征来优化条带噪声;
步骤五:在CBCT图像的基础上通过对噪声区域的插值方法达到消除环形伪影的目的。
进一步地,所述步骤二中通过滑动窗口差值提取边缘特征的过程包括以下步骤:
1)对CBCT图像进行灰度化并转换至极坐标下,此时环形伪影的特征表现为垂直直线;
2)在极坐标系下,通过滑动窗口,搜索窗口内的最大值和最小值,将最大值与最小值的差值赋给滑动窗口中心点,最终生成差值图像;
3)在差值图像的基础上,通过阈值法进行二值化操作,生成边缘特征图像。
进一步地,在步骤2)中,通过窗口内最大值和最小值之间的差值消除条带部分的影响,滑动窗口大小可以设置为3至9区间内的奇数。
进一步地,在步骤3)中通过阈值法进行二值化操作,阈值的大小可以设置为差值图像的平均数。
进一步地,所述步骤三中的通过水平方向的梯度图像检测条带噪声的过程包括以下步骤:
1)对CBCT图像进行灰度化并转换至极坐标下;
2)在极坐标下,利用Sobel算子求取水平方向的梯度图像,再通过区域生长法,在八邻域内去除小区域目标;
3)计算每一列中,连续不为零的直线;设置阈值,将长度小于等于阈值的列全部赋值为0,大于阈值的赋值为1,生成条带噪声图像。
进一步地,所述步骤四种利用边缘特征来优化条带噪声的过程包含为以下步骤:
1)将边缘特征图像和条带噪声图像进行叠加,在叠加的情况下,边缘特征图像与条带噪声图像在任意点具有相同的位置;
2)遍历条带噪声图像,在任意位置下,当条带噪声图像的像素值不为零,且边缘特征图像的像素值为零的情况下,视该位置点为噪声点,赋值为1,否则赋值为0;
3)对条带噪声图像进行形态学操作即进行膨胀操作,扩大边缘。
进一步地,所述步骤五中通过对噪声区域的插值方法来消除噪声的过程包含为以下步骤:
1)将条带噪声二值图像转换到直角坐标系下,此时与原始CBCT图像具有相同的坐标;
2)遍历原始CBCT图像,在条带噪声二值图像中的条带噪声区域即像素值不为零区域,利用原始CBCT图像噪声周边区域的像素点对其进行内插,最终获取去噪图像。
本发明的有益效果:
本发明针对CBCT图像中存在的特征微小的环形伪影,通过检测边缘特征来优化条带噪声,在不对整体图像处理的前提下,最大限度地消除环形伪影。而且对于CBCT重建图像的每一层而言,其环形伪影具有一致性,即均处在相同的位置,所以检测一层数据的环形伪影即可对每一层数据进行环形伪影消除。本发明提出一种快速、保真的CBCT环形伪影消除方法。
附图说明
为了便于本领域技术人员理解,下面结合附图对本发明作进一步的说明。
图1是本发明消除CBCT环形伪影的一种方法流程图;
图2-a是本发明实时案例所采用的CBCT一层切片数据,为原始数据;
图2-b是本发明CBCT原始数据转换至极坐标系下的图像;
图3-a是本发明在极坐标系下,采用Sobel算子提取水平方向的梯度图像;
图3-b是本发明在极坐标系下,通过滑动窗口差值法获取的边缘特征图像;
图4-a是本发明在极坐标系下,优化提取后的条带噪声直线;
图4-b是本发明检测的条带噪声直线转换至直角坐标系下图像;
图5是本发明在检测条带噪声的基础上进行环形伪影消除的结果图像。
具体实施方式
为能清楚说明本方案的技术特点,下面通过具体实施方式,并结合其附图,对本发明进行详细阐述。下文的公开提供了实施例用来实现本发明的结构。为了简化本发明的公开,下文中对特定例子的部件和设置进行描述。此外,本发明可以在不同例子中重复参考数字和/或字母。这种重复是为了简化和清楚的目的,其本身不指示所讨论各种实施例和/或设置之间的关系。应当注意,在附图中所图示的部件不一定按比例绘制。本发明省略了对公知组件和处理技术及工艺的描述以避免不必要地限制本发明。
如图1所示,本发明的一种基于滑动窗口差值和条带噪声检测的CBCT环形伪影消除方法,它对特征微小的环形伪影,通过滑动窗口差值提取边缘特征,通过水平方向的梯度图像检测条带噪声,再利用边缘特征来优化条带噪声,在CBCT图像的基础上通过对噪声区域的插值方法达到消除环形伪影的目的。
下面以CBCT一层切片数据为实施案例,结合附图和技术方案进一步说明本发明的具体实施方式。
通过滑动窗口差值提取边缘特征的过程包括以下步骤:
1)对CBCT图像进行灰度化并转换至极坐标下,此时环形伪影的特征表现为垂直直线,如附图2-b所示,特征微弱的环形伪影在极坐标下也同样不是特别明显;
2)在极坐标系下,通过滑动窗口,搜索窗口内的最大值和最小值,将最大值与最小值的差值赋给滑动窗口中心点,最终生成差值图像。由于条带噪声比周边要暗,所以通过滑动窗口差值将会削弱条带噪声的影响;
3)在差值图像的基础上,通过阈值法进行二值化操作,生成边缘特征图像,见附图3-b;
在步骤2)中,通过窗口内最大值和最小值之间的差值消除条带部分的影响,滑动窗口大小一般设置为3至9区间内的奇数。
在步骤3)中,阈值的大小通常设置为差值图像的平均数。
通过水平方向的梯度图像检测条带噪声的过程包括以下步骤:
1)对CBCT图像进行灰度化并转换至极坐标下;
2)在极坐标下,利用Sobel算子求取水平方向的梯度图像(见附图3-a),其中,所选择的Sobel算子为3×3矩阵,按行排列依次为[-101],[-202],[-101]。在Sobel算子的基础上再通过区域生长法,在八邻域内去除小区域目标;
3)计算每一列中,连续不为零的直线;设置阈值为100,将长度小于等于阈值的列全部赋值为0,大于阈值的赋值为1,生成条带噪声图像;
利用边缘特征来优化条带噪声的过程包含为以下步骤:
1)将边缘特征图像和条带噪声图像进行叠加,在叠加的情况下,边缘特征图像与条带噪声图像在任意点具有相同的位置;
2)遍历条带噪声图像,在任意位置下,当条带噪声图像的像素值不为零,且边缘特征图像的像素值为零的情况下,视该位置点为噪声点,赋值为1,否则赋值为0,见附图4-a;
3)对条带噪声图像进行形态学操作(膨胀),扩大边缘;
通过对噪声区域的插值方法来消除噪声的过程包含为以下步骤:
1)将条带噪声二值图像转换到直角坐标系下,此时与原始CBCT图像具有相同的坐标;
2)遍历原始CBCT图像,在条带噪声二值图像中的条带噪声区域(像素值不为零),利用原始CBCT图像噪声周边区域的像素点对其进行内插,最终获取去噪图像,见附图5;
本发明针对CBCT图像中存在的特征微小的环形伪影,通过检测边缘特征来优化条带噪声,在不对整体图像处理的前提下,最大限度地消除环形伪影。而且对于CBCT重建图像的每一层而言,其环形伪影具有一致性,即均处在相同的位置,所以检测一层数据的环形伪影即可对每一层数据进行环形伪影消除。本发明提出一种快速、保真的CBCT环形伪影消除方法。
以上内容仅仅是对本发明结构所作的举例和说明,所属本技术领域的技术人员对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,只要不偏离发明的结构或者超越本权利要求书所定义的范围,均应属于本发明的保护范围。
Claims (5)
1.基于滑动窗口差值和条带噪声检测的CBCT环形伪影消除方法,其特征在于,该方法包括下述步骤:
步骤一:对CBCT图像中存在的特征微小的环形伪影;
步骤二:通过滑动窗口差值提取边缘特征;
步骤三:通过水平方向的梯度图像检测条带噪声;
步骤四:再利用边缘特征来优化条带噪声;
步骤五:在CBCT图像的基础上通过对噪声区域的插值方法达到消除环形伪影的目的;
所述步骤四种利用边缘特征来优化条带噪声的过程包含为以下步骤:
1)将边缘特征图像和条带噪声图像进行叠加,在叠加的情况下,边缘特征图像与条带噪声图像在任意点具有相同的位置;
2)遍历条带噪声图像,在任意位置下,当条带噪声图像的像素值不为零,且边缘特征图像的像素值为零的情况下,视该位置点为噪声点,赋值为1,否则赋值为0;
3)对条带噪声图像进行形态学操作即进行膨胀操作,扩大边缘;
所述步骤五中通过对噪声区域的插值方法来消除噪声的过程包含为以下步骤:
1)将条带噪声二值图像转换到直角坐标系下,此时与原始CBCT图像具有相同的坐标;
2)遍历原始CBCT图像,在条带噪声二值图像中的条带噪声区域即像素值不为零区域,利用原始CBCT图像噪声周边区域的像素点对其进行内插,最终获取去噪图像。
2.根据权利要求1所述的基于滑动窗口差值和条带噪声检测的CBCT环形伪影消除方法,其特征在于,所述步骤二中通过滑动窗口差值提取边缘特征的过程包括以下步骤:
1)对CBCT图像进行灰度化并转换至极坐标下,此时环形伪影的特征表现为垂直直线;
2)在极坐标系下,通过滑动窗口,搜索窗口内的最大值和最小值,将最大值与最小值的差值赋给滑动窗口中心点,最终生成差值图像;
3)在差值图像的基础上,通过阈值法进行二值化操作,生成边缘特征图像。
3.根据权利要求2所述的基于滑动窗口差值和条带噪声检测的CBCT环形伪影消除方法,其特征在于,在步骤2)中,通过窗口内最大值和最小值之间的差值消除条带部分的影响,滑动窗口大小可以设置为3至9区间内的奇数。
4.根据权利要求2所述的基于滑动窗口差值和条带噪声检测的CBCT环形伪影消除方法,其特征在于,在步骤3)中通过阈值法进行二值化操作,阈值的大小可以设置为差值图像的平均数。
5.根据权利要求1所述的基于滑动窗口差值和条带噪声检测的CBCT环形伪影消除方法,其特征在于,所述步骤三中的通过水平方向的梯度图像检测条带噪声的过程包括以下步骤:
1)对CBCT图像进行灰度化并转换至极坐标下;
2)在极坐标下,利用Sobel算子求取水平方向的梯度图像,再通过区域生长法,在八邻域内去除小区域目标;
3)计算每一列中,连续不为零的直线;设置阈值,将长度小于等于阈值的列全部赋值为0,大于阈值的赋值为1,生成条带噪声图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811504142.6A CN109636872B (zh) | 2018-12-10 | 2018-12-10 | 基于滑动窗口差值和条带噪声检测的cbct环形伪影消除方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811504142.6A CN109636872B (zh) | 2018-12-10 | 2018-12-10 | 基于滑动窗口差值和条带噪声检测的cbct环形伪影消除方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109636872A CN109636872A (zh) | 2019-04-16 |
CN109636872B true CN109636872B (zh) | 2022-12-09 |
Family
ID=66072367
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811504142.6A Active CN109636872B (zh) | 2018-12-10 | 2018-12-10 | 基于滑动窗口差值和条带噪声检测的cbct环形伪影消除方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109636872B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110111318B (zh) * | 2019-04-30 | 2021-06-29 | 上海联影医疗科技股份有限公司 | 一种环形伪影的检测方法和系统 |
CN110889878B (zh) * | 2019-11-08 | 2023-06-09 | 湖北科技学院 | 结合图像平滑处理的ct图像环形伪影抑制方法 |
CN117078791B (zh) * | 2023-10-13 | 2024-01-12 | 俐玛精密测量技术(苏州)有限公司 | 一种ct环形伪影校正方法、装置、电子设备及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103593826A (zh) * | 2013-10-22 | 2014-02-19 | 南方医科大学 | 图像环形伪影校正方法 |
CN105321155A (zh) * | 2015-10-29 | 2016-02-10 | 北京理工大学 | 一种cbct图像环形伪影消除方法 |
WO2018103015A1 (zh) * | 2016-12-07 | 2018-06-14 | 深圳先进技术研究院 | 一种环形伪影修正的方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7860341B2 (en) * | 2006-07-21 | 2010-12-28 | Varian Medical Systems, Inc. | System and method for correcting for ring artifacts in an image |
-
2018
- 2018-12-10 CN CN201811504142.6A patent/CN109636872B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103593826A (zh) * | 2013-10-22 | 2014-02-19 | 南方医科大学 | 图像环形伪影校正方法 |
CN105321155A (zh) * | 2015-10-29 | 2016-02-10 | 北京理工大学 | 一种cbct图像环形伪影消除方法 |
WO2018103015A1 (zh) * | 2016-12-07 | 2018-06-14 | 深圳先进技术研究院 | 一种环形伪影修正的方法及装置 |
Non-Patent Citations (1)
Title |
---|
一种基于边缘检测的环状伪影定位算法;奚卫东等;《医疗卫生装备》;20100115(第01期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109636872A (zh) | 2019-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109636872B (zh) | 基于滑动窗口差值和条带噪声检测的cbct环形伪影消除方法 | |
US11232543B2 (en) | System and method for image correction | |
CN105321155A (zh) | 一种cbct图像环形伪影消除方法 | |
US11151414B2 (en) | Training image enhancement model and enhancing image | |
CN102068281B (zh) | 一种占位性病变超声图像的处理方法 | |
CN102800073B (zh) | 一种锥束ct环形伪影的自动判别与校正方法 | |
JP5994858B2 (ja) | 放射線断層画像生成装置および放射線断層画像生成方法 | |
CN110060316B (zh) | 一种用于ct重建中多区域分割的环伪影校正方法 | |
CN102521801B (zh) | Ct图像的环形伪影及弧形伪影校正方法 | |
CN106886982A (zh) | Cbct图像环形伪影去除方法 | |
CN112102428B (zh) | Ct锥形束扫描图像重建方法、扫描系统及存储介质 | |
CN111353958A (zh) | 图像处理方法、装置及系统 | |
CN105469366A (zh) | 一种消减ct图像金属伪影的解析方法 | |
CN106530236A (zh) | 一种医学图像处理方法及系统 | |
CN105374023B (zh) | 一种分割目标区域方法,及其图像重建方法和装置 | |
CN106651981B (zh) | 一种环形伪影修正的方法及装置 | |
CN103226815A (zh) | 一种低剂量ct图像滤波方法 | |
CN106910166A (zh) | 一种图像处理方法和装置 | |
CN112070856B (zh) | 基于非下采样轮廓波变换的有限角c型臂ct图像重建方法 | |
CN110111318B (zh) | 一种环形伪影的检测方法和系统 | |
CN102073991A (zh) | 增强图像中的空间局域化征兆的对比度的方法 | |
CN112669406B (zh) | 一种基于投影估计的相控阵线圈磁共振图像非均匀性校正方法 | |
CN115588060A (zh) | 基于投影平均图像的快速ct环形伪影校正方法及系统 | |
CN115375560A (zh) | 一种3d-dsa图像的重建方法及系统 | |
CN112446840B (zh) | 一种基于深度学习的ct图像黑带伪影消除方法及系统 |
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 |