CN104644203B - 一种剂量调制扫描方法和装置 - Google Patents
一种剂量调制扫描方法和装置 Download PDFInfo
- Publication number
- CN104644203B CN104644203B CN201410445744.4A CN201410445744A CN104644203B CN 104644203 B CN104644203 B CN 104644203B CN 201410445744 A CN201410445744 A CN 201410445744A CN 104644203 B CN104644203 B CN 104644203B
- Authority
- CN
- China
- Prior art keywords
- scanned
- data layer
- noise reduction
- value
- current position
- 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
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000013442 quality metrics Methods 0.000 claims abstract description 51
- 238000011946 reduction process Methods 0.000 claims abstract description 30
- 239000003795 chemical substances by application Substances 0.000 claims description 50
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 25
- 241000208340 Araliaceae Species 0.000 claims 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 2
- 235000008434 ginseng Nutrition 0.000 claims 2
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002146 bilateral effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 230000032696 parturition Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 102000037983 regulatory factors Human genes 0.000 description 1
- 108091008025 regulatory factors Proteins 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/54—Control of apparatus or devices for radiation diagnosis
- A61B6/542—Control of apparatus or devices for radiation diagnosis involving control of exposure
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- 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/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- 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/10081—Computed x-ray tomography [CT]
-
- 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/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
本发明实施例公开了一种剂量调制扫描方法和装置,所述方法包括:预先设定待扫描部位的目标图像质量标准以及确定所述待扫描部位的任一数据层为参考数据层,根据所述目标图像质量标准和所述参考数据层,计算当前待扫描位置的理论扫描剂量值;在所述理论扫描剂量值下,计算所述当前待扫描位置的各个数据层的衰减值与所述参考数据层的衰减值的差异;当所述当前待扫描位置中存在差异大于预设阈值的数据层时,将所述数据层确定为待降噪对象;对所述待降噪对象进行降噪处理,以得到符合所述目标图像质量标准的扫描图像。本发明通过对待降噪对象的确定以及降噪处理,解决现有剂量调制扫描中出现的单次球管曝光扫描图像不清楚的问题,改善扫描图像的质量。
Description
技术领域
本发明涉及医疗技术领域,具体涉及一种剂量调制扫描方法和装置。
背景技术
剂量调制扫描,是指在扫描过程中,随扫描位置的变化,扫描剂量不断发生变化的一种扫描方式。采用剂量调制扫描方式进行CT临床扫描,目的是实现待扫描的各个部位噪声的一致性。
但是,由于在同一扫描位置,一次球管曝光时,曝光剂量是固定值。如果剂量调制扫描中单次球管曝光时,对应的当前扫描位置的Z向范围的衰减值变化较大,则会导致当前扫描数据的噪声有明显差异,这就违背了剂量调制扫描的设计初衷。例如,当单次球管在肩颈交接处曝光时,由于肩部明显比颈部面积大,所以,对于肩颈交接处扫描结果的Z向范围的衰减值可能变化较大。也就是说,可能导致面积较大的肩部扫描结果不清楚。
所以,目前我们需要想办法来解决剂量调制扫描时,存在的这个问题。
发明内容
本发明提供了一种剂量调制扫描方法和装置,用以解决现有剂量调制扫描中出现的单次球管曝光扫描图像不清楚的问题,改善扫描图像的质量。
本发明提供了一种剂量调制扫描方法,预先设定待扫描部位的目标图像质量标准以及确定所述待扫描部位的任一数据层为参考数据层,所述方法包括:
根据所述目标图像质量标准和所述参考数据层,计算所述待扫描部位的当前待扫描位置的理论扫描剂量值;
在所述理论扫描剂量值下,计算所述当前待扫描位置的各个数据层的衰减值与所述参考数据层的衰减值的差异;
当所述当前待扫描位置中存在差异大于预设阈值的数据层时,将所述数据层确定为待降噪对象;
对所述待降噪对象进行降噪处理,以得到符合所述目标图像质量标准的扫描图像。
优选地,所述根据所述目标图像质量标准和所述参考数据层,计算当前待扫描位置的理论扫描剂量值,包括:
基于所述目标图像质量标准,计算所述参考数据层的参考扫描剂量值;
根据所述参考扫描剂量值和所述当前待扫描位置的数据层的个数,计算所述当前待扫描位置的理论扫描剂量值。
优选地,所述在所述理论扫描剂量值下,计算所述当前待扫描位置的各个数据层的衰减值与所述参考数据层的衰减值的差异,包括:
在所述理论扫描剂量值下,计算所述当前待扫描位置的各个数据层的等效水模衰减直径;
计算所述当前待扫描位置的各个数据层的等效水模衰减直径与所述参考数据层的等效水模衰减直径的差异。
优选地,所述对所述待降噪对象进行降噪处理,以得到符合所述目标图像质量标准的扫描图像,包括:
获取所述待降噪对象的生数据;
基于所述生数据,对所述待降噪对象进行降噪,以得到符合所述目标图像质量标准的扫描图像。
优选地,所述对所述待降噪对象进行降噪处理,以得到符合所述目标图像质量标准的扫描图像,包括:
利用所述理论扫描剂量值,对所述当前待扫描位置进行扫描得到初步扫描图像;
根据所述待降噪对象,对初步扫描图像进行降噪处理,以得到符合所述目标图像质量标准的扫描图像。
本发明还提供了一种剂量调制扫描装置,所述装置包括:
设定模块,用于预先设定待扫描部位的目标图像质量标准;
第一确定模块,用于确定所述待扫描部位的任一数据层为参考数据层;
第一计算模块,用于根据所述目标图像质量标准和所述参考数据层,计算所述待扫描部位的当前待扫描位置的理论扫描剂量值;
第二计算模块,用于在所述理论扫描剂量值下,计算所述当前待扫描位置的各个数据层的衰减值与所述参考数据层的衰减值的差异;
第二确定模块,用于当所述当前待扫描位置中存在差异大于预设阈值的数据层时,将所述数据层确定为待降噪对象;
降噪模块,用于对所述待降噪对象进行降噪处理,以得到符合所述目标图像质量标准的扫描图像。
优选地,所述第一计算模块,包括:
第一计算子模块,用于基于所述目标图像质量标准,计算所述参考数据层的参考扫描剂量值;
第二计算子模块,用于根据所述参考扫描剂量值和所述当前待扫描位置的数据层的个数,计算所述当前待扫描位置的理论扫描剂量值。
优选地,所述第二计算模块,包括:
第三计算子模块,用于在所述理论扫描剂量值下,计算所述当前待扫描位置的各个数据层的等效水模衰减直径;
第四计算子模块,用于计算所述当前待扫描位置的各个数据层的等效水模衰减直径与所述参考数据层的等效水模衰减直径的差异。
优选地,所述降噪模块,包括:
获取子模块,用于获取所述待降噪对象的生数据;
第一降噪子模块,用于基于所述生数据,对所述待降噪对象进行降噪,以得到符合所述目标图像质量标准的扫描图像。
优选地,所述降噪模块,包括:
扫描子模块,用于利用所述理论扫描剂量值,对所述当前待扫描位置进行扫描得到初步扫描图像;
第二降噪子模块,用于根据所述待降噪对象,对初步扫描图像进行降噪处理,以得到符合所述目标图像质量标准的扫描图像。
本发明预先设定待扫描部位的目标图像质量标准以及确定所述待扫描部位的任一数据层为参考数据层,并根据所述目标图像质量标准和所述参考数据层,计算当前待扫描位置的理论扫描剂量值。其次,在所述理论扫描剂量值下,计算所述当前待扫描位置的各个数据层的衰减值与所述参考数据层的衰减值的差异。当所述当前待扫描位置中存在差异大于预设阈值的数据层时,将所述数据层确定为待降噪对象;最终,对所述待降噪对象进行降噪处理,以得到符合所述目标图像质量标准的扫描图像。与现有技术的剂量调制扫描方式相比,本发明通过对待降噪对像的确定以及降噪处理,可以解决现有剂量调制扫描中出现的单次球管曝光扫描图像不清楚的问题,改善扫描图像的质量。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例一提供的剂量调制扫描方法流程图;
图2为本发明实施例一提供的当前待扫描位置的任一数据层的生数据坐标示意图;
图3为本发明实施例二提供的剂量调制扫描装置结构示意图。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
目前,利用传统的剂量调制扫描方法进行扫描时,期望各个扫描位置的扫描图像噪声一致。但是实际扫描时,由于在同一扫描位置,一次球管曝光时曝光剂量是固定值,若当前扫描位置的Z向扫描范围的衰减值与参考数据层比较变化较大时,就会导致扫描图像的Z向的数据噪声差异较大。为了避免上述情况的发生,本发明实施例提供的剂量调制扫描方法,实时的将当前待扫描位置的各个数据层与参考数据层进行比较,当差异较大时,对所述数据层进行降噪处理,最终得到符合目标图像质量标准的扫描图像。
实施例一
参考图1,图1为本发明实施例提供的剂量调制扫描方法流程图,所述方法包括:
S101:预先设定待扫描部位的目标图像质量标准。
本实施例中,启用剂量调制扫描方法之前,首先设定期望的图像质量标准,即目标图像质量标准。具体的设定方式或设定参数不限,目标图像质量标准可以是期望的图像噪声水准,也可以是期望的图像信噪比,也可以是期望的剂量水准等。
以目标图像质量标准为图像信噪比为例,其中,图像信噪比的设置范围可以为0.3~3.0。参数1.0代表当前目标图像的图像信噪比等于当前扫描协议的标准图像信噪比。所谓的当前扫描协议的标准图像信噪比,是指CT扫描机出厂的默认图像信噪比设定值。同理,参数0.3和3.0分别代表目标图像的图像信噪比为当前扫描协议的标准图像信噪比的0.3或3.0倍。
S102:确定所述待扫描部位的任一数据层为参考数据层。
本实施例中,待扫描部位可以包括若干数据层,可以将其中任意一个数据层作为参考数据层,以便后续作为各个数据层的比较标准。
实际应用中,尽量将所述待扫描部位中可以代表各个数据层的平均情况的数据层作为参考数据层。具体的,可以将待扫描部位的中间数据层确定为参考数据层。
值得注意的是,S101和S102的执行顺序不做限制。
S103:根据所述目标图像质量标准和所述参考数据层,计算所述待扫描部位的当前待扫描位置的理论扫描剂量值。
本实施例中,以预先设定的目标图像质量标准作为标准,为了达到所述标准,计算所述待扫描部位的当前待扫描位置理论上需要的扫描剂量值,即理论扫描剂量值。
实际操作中,首先,基于预先设定的所述目标图像质量标准,计算所述参考数据层的参考扫描剂量值。也就是说,首先计算参考数据层为了达到所述目标图像质量标准需要的扫描剂量值,即参考扫描剂量值。其次,根据所述参考扫描剂量值和所述当前待扫描位置的数据层的个数,计算所述当前待扫描位置的理论扫描剂量值。也就是说,以所述参考数据层需要的参考扫描剂量值为单位,计算当前待扫描位置理论上需要的扫描剂量值,即理论扫描剂量值。可以将当前待扫描位置的理论扫描剂量值理解为参考扫描剂量值与当前待扫描位置的数据层的个数的乘积。
另外,本实施例的当前待扫描位置的理论扫描剂量值的计算思想为,假设CT扫描机对当前扫描位置扫描一圈的总剂量不变,根据当前圈的每个扫描投影位置的衰减变化,对剂量进行分配。此处提到的衰减变化是针对参考数据层的当前投影位置下的最大衰减值。CT扫描一般会覆盖多层检测器,本实施例计算当前待扫描位置的理论扫描剂量值时,一般按照预设的原则,选择某一层检测器的数据作为参考数据层。
具体的,计算当前待扫描位置的理论扫描剂量值可以利用如下公式:
公式(1);
N0,All=Np*N0 公式(2);
其中,N0为恒定剂量下,CT扫描仪的中心通道对应的X线发射强度。恒定剂量下,每个投影位置的N0都相等;
Np为CT扫描仪对当前位置扫描一圈的总投影位置个数;
Amax,i为第i个投影位置所有通道的最大衰减值;
Nmod,i为第i个投影位置所需的扫描剂量值。
S104:在所述理论扫描剂量值下,计算所述当前待扫描位置的各个数据层的衰减值与所述参考数据层的衰减值的差异。
本实施例中,在通过计算得到所述当前待扫描位置的理论扫描剂量值后,以所述理论扫描剂量值为前提,计算所述当前待扫描位置的各个数据层的衰减值,并确定各个数据层的衰减值与所述参考数据层的衰减值的差异。
实际操作中,可以通过将各个数据层的等效水模衰减直径与所述参考数据层的等效水模衰减直径进行比较。具体的,首先,在所述理论扫描剂量值下,计算所述当前待扫描位置的各个数据层的等效水模衰减直径。其次,计算所述当前待扫描位置的各个数据层的等效水模衰减直径与所述参考数据层的等效水模衰减直径的差异。
实际应用中,针对所述当前待扫描位置的每层数据,计算等效水模衰减直径的具体方法如下:
首先,获取所述当前待扫描位置的各个数据层的生数据,其中,所述生数据由N个检测器通道采集获得,即利用公式(3)获得;
rawdata={μ0l0,μ1l1,......μili,},i=0,1,2,3,......,N-1 公式(3);
其次,基于获取的各个数据层的生数据,计算各个数据层的总衰减面积。如图2,图2为当前待扫描位置的任一数据层的生数据坐标示意图。
具体的,计算任一数据层的总衰减面积可以利用如下公式:
公式(4);
公式(5);
其中,N为通道数;μi为物体在第i个通道上的平均衰减系数;li为物体在第i个通道上的路径;Δ为相邻检测器之间的距离;R为旋转半径;α为检测器的扇角。
在计算任一数据层的总衰减面积S时,需要考虑随床高的变化,也就是说,需要对Δ值进行修正。如果待扫描部位位于扫描仪中心,则无需修正Δ值。如果待扫描部位向上偏离中心,则Δ值需要变小。如果待扫描部位向下偏离中心,则Δ值需要变大。Δ值可以基于CT扫描仪的几何关系来进行修正。
最后,基于所述数据层的总衰减面积S,计算该数据层的等效水模衰减直径Dscan,具体如下:
Dscan=2*sqrt(mean(S)/(PI*μwater)) 公式(6)
其中,mean(s)为所述当前待扫描位置的Z方向所有平片切面衰减面积的均值;μwater为水的衰减系数。
S105:当所述当前待扫描位置中存在差异大于预设阈值的数据层时,将所述数据层确定为待降噪对象。
本实施例中,预先设置一个阈值作为差异标准,具体的,阈值的大小可以根据用户对扫描图像质量的要求设定。
实际操作中,在确定所述当前待扫描位置中各个数据层与参考数据层的差异后,判断所述差异是否大于预设的阈值,如果大于,则将大于所述阈值的数据层确定为待降噪对象。例如,当前待扫描位置为肩颈交接处时,由于肩部明显比颈部面积大,所以,对于肩颈交接处扫描结果的Z向范围的衰减值可能变化较大。此时,肩颈交接处的各个数据层与参考数据层的差异可能大于预设的阈值,可以将这些数据层确定为待降噪对象。
S106:对所述待降噪对象进行降噪处理,以得到符合所述目标图像质量标准的扫描图像。
本实施例中,在确定待降噪对象之后,需要对所述待降噪对象进行降噪处理,最终得到符合所述目标图像质量标准的扫描图像。
实际应用中,利用剂量调制扫描方法进行扫描过程中,可以针对所述待降噪对象的生数据进行降噪,并保存降噪后的生数据。最终利用降噪后的生数据得到符合所述目标图像质量标准的扫描图像。
具体的,首先获取所述待降噪对象的生数据。其次,基于所述生数据,对所述待降噪对象进行降噪,以得到符合所述目标图像质量标准的扫描图像。
另外,本实施例还可以基于扫描图像进行降噪处理。具体的,在确定待降噪对象之后,直接在所述带降噪对象的生数据的头通道中进行标记,并保存未经过降噪处理的原始的生数据。在后续的图像重建过程中,根据标记再对图像进行降噪处理。
实际应用中,首先利用所述理论扫描剂量值,对所述当前待扫描位置进行扫描得到初步扫描图像。其次,根据所述待降噪对象,对初步扫描图像进行降噪处理,以得到符合所述目标图像质量标准的扫描图像。
值得注意的是,本实施例还可以基于生数据和扫描图像的组合进行降噪处理,具体的降噪方法不做限制。
实际应用中,基于生数据的降噪方法可以为双边滤波的降噪方法。如下公式所示:
公式(7)
公式(8)
公式(9)
其中,CT扫描仪采集的生数据是由投影位置和通道组成的二维数据集,其中i、j分别表示二维的数据点位置,Qi,Qj表示数据点的值,为卷积结果数据,d,σ是调节因子。从公式中可以知道W1是以到中心点距离为自变量的高斯卷积,W2是以与中心点差值为自变量的高斯卷积,可以认为是一种置信度因子,差值越大,该值越小,即对降噪的贡献越小。
另外,基于扫描图像的降噪方法可以为自适应的维纳滤波器降噪方法。其中,自适应的维纳滤波器能根据扫描图象的局部方差来调整滤波器的输出,局部方差越大,滤波器的平滑作用越强。该方法的最终目标是使恢复图像f^(x,y)与原始图像f(x,y)的均方误差e2=E[(f(x,y)-f^(x,y)2]最小。值得注意的是,该方法对保留扫描图像的边缘和其他高频部分作用较大。
本实施例中,预先设定待扫描部位的目标图像质量标准以及确定所述待扫描部位的任一数据层为参考数据层,并根据所述目标图像质量标准和所述参考数据层确定当前待扫描位置的待降噪对象。最终,对所述待降噪对象进行降噪处理,得到符合所述目标图像质量标准的扫描图像。本实施例能够通过对待降噪对像的确定以及降噪处理,解决现有剂量调制扫描中出现的单次球管曝光扫描图像不清楚的问题,改善扫描图像的质量。
实施例二
参考图3,图3为本实施例提供的一种剂量调制扫描装置结构示意图,所述装置包括:
设定模块301,用于预先设定待扫描部位的目标图像质量标准;
第一确定模块302,用于确定所述待扫描部位的任一数据层为参考数据层;
第一计算模块303,用于根据所述目标图像质量标准和所述参考数据层,计算所述待扫描部位的当前待扫描位置的理论扫描剂量值;
第二计算模块304,用于在所述理论扫描剂量值下,计算所述当前待扫描位置的各个数据层的衰减值与所述参考数据层的衰减值的差异;
第二确定模块305,用于当所述当前待扫描位置中存在差异大于预设阈值的数据层时,将所述数据层确定为待降噪对象;
降噪模块306,用于对所述待降噪对象进行降噪处理,以得到符合所述目标图像质量标准的扫描图像。
其中,所述第一计算模块可以包括:
第一计算子模块,用于基于所述目标图像质量标准,计算所述参考数据层的参考扫描剂量值;
第二计算子模块,用于根据所述参考扫描剂量值和所述当前待扫描位置的数据层的个数,计算所述当前待扫描位置的理论扫描剂量值。
其中,所述第二计算模块可以包括:
第三计算子模块,用于在所述理论扫描剂量值下,计算所述当前待扫描位置的各个数据层的等效水模衰减直径;
第四计算子模块,用于计算所述当前待扫描位置的各个数据层的等效水模衰减直径与所述参考数据层的等效水模衰减直径的差异。
在一种优选实施例中,所述降噪模块可以包括:
获取子模块,用于获取所述待降噪对象的生数据;
第一降噪子模块,用于基于所述生数据,对所述待降噪对象进行降噪,以得到符合所述目标图像质量标准的扫描图像。
在另一种优选实施例中,所述降噪模块可以包括:
扫描子模块,用于利用所述理论扫描剂量值,对所述当前待扫描位置进行扫描得到初步扫描图像;
第二降噪子模块,用于根据所述待降噪对象,对初步扫描图像进行降噪处理,以得到符合所述目标图像质量标准的扫描图像。
本实施例中,预先设定待扫描部位的目标图像质量标准以及确定所述待扫描部位的任一数据层为参考数据层,并根据所述目标图像质量标准和所述参考数据层确定当前待扫描位置的待降噪对象。最终,对所述待降噪对象进行降噪处理,得到符合所述目标图像质量标准的扫描图像。本实施例能够通过对待降噪对像的确定以及降噪处理,解决现有剂量调制扫描中出现的单次球管曝光扫描图像不清楚的问题,改善扫描图像的质量。
对于装置实施例而言,由于其基本对应于方法实施例,所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上对本发明实施例所提供的剂量调制扫描方法和装置进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种剂量调制扫描方法,其特征在于,预先设定待扫描部位的目标图像质量标准以及确定所述待扫描部位的任一数据层为参考数据层,所述参考数据层为代表所述待扫描部位的各个数据层的平均情况的数据层,所述方法包括:
根据所述目标图像质量标准和所述参考数据层,计算所述待扫描部位的当前待扫描位置的理论扫描剂量值;
在所述理论扫描剂量值下,计算所述当前待扫描位置的各个数据层的衰减值与所述参考数据层的衰减值的差异;
当所述当前待扫描位置中存在差异大于预设阈值的数据层时,将所述数据层确定为待降噪对象;
对所述待降噪对象进行降噪处理,以得到符合所述目标图像质量标准的扫描图像。
2.根据权利要求1所述的方法,其特征在于,所述根据所述目标图像质量标准和所述参考数据层,计算当前待扫描位置的理论扫描剂量值,包括:
基于所述目标图像质量标准,计算所述参考数据层的参考扫描剂量值;
根据所述参考扫描剂量值和所述当前待扫描位置的数据层的个数,计算所述当前待扫描位置的理论扫描剂量值。
3.根据权利要求1所述的方法,其特征在于,所述在所述理论扫描剂量值下,计算所述当前待扫描位置的各个数据层的衰减值与所述参考数据层的衰减值的差异,包括:
在所述理论扫描剂量值下,计算所述当前待扫描位置的各个数据层的等效水模衰减直径;
计算所述当前待扫描位置的各个数据层的等效水模衰减直径与所述参考数据层的等效水模衰减直径的差异。
4.根据权利要求1所述的方法,其特征在于,所述对所述待降噪对象进行降噪处理,以得到符合所述目标图像质量标准的扫描图像,包括:
获取所述待降噪对象的生数据;
基于所述生数据,对所述待降噪对象进行降噪,以得到符合所述目标图像质量标准的扫描图像。
5.根据权利要求1所述的方法,其特征在于,所述对所述待降噪对象进行降噪处理,以得到符合所述目标图像质量标准的扫描图像,包括:
利用所述理论扫描剂量值,对所述当前待扫描位置进行扫描得到初步扫描图像;
根据所述待降噪对象,对初步扫描图像进行降噪处理,以得到符合所述目标图像质量标准的扫描图像。
6.一种剂量调制扫描装置,其特征在于,所述装置包括:
设定模块,用于预先设定待扫描部位的目标图像质量标准;
第一确定模块,用于确定所述待扫描部位的任一数据层为参考数据层,所述参考数据层为代表所述待扫描部位的各个数据层的平均情况的数据层;
第一计算模块,用于根据所述目标图像质量标准和所述参考数据层,计算所述待扫描部位的当前待扫描位置的理论扫描剂量值;
第二计算模块,用于在所述理论扫描剂量值下,计算所述当前待扫描位置的各个数据层的衰减值与所述参考数据层的衰减值的差异;
第二确定模块,用于当所述当前待扫描位置中存在差异大于预设阈值的数据层时,将所述数据层确定为待降噪对象;
降噪模块,用于对所述待降噪对象进行降噪处理,以得到符合所述目标图像质量标准的扫描图像。
7.根据权利要求6所述的装置,其特征在于,所述第一计算模块,包括:
第一计算子模块,用于基于所述目标图像质量标准,计算所述参考数据层的参考扫描剂量值;
第二计算子模块,用于根据所述参考扫描剂量值和所述当前待扫描位置的数据层的个数,计算所述当前待扫描位置的理论扫描剂量值。
8.根据权利要求6所述的装置,其特征在于,所述第二计算模块,包括:
第三计算子模块,用于在所述理论扫描剂量值下,计算所述当前待扫描位置的各个数据层的等效水模衰减直径;
第四计算子模块,用于计算所述当前待扫描位置的各个数据层的等效水模衰减直径与所述参考数据层的等效水模衰减直径的差异。
9.根据权利要求6所述的装置,其特征在于,所述降噪模块,包括:
获取子模块,用于获取所述待降噪对象的生数据;
第一降噪子模块,用于基于所述生数据,对所述待降噪对象进行降噪,以得到符合所述目标图像质量标准的扫描图像。
10.根据权利要求6所述的装置,其特征在于,所述降噪模块,包括:
扫描子模块,用于利用所述理论扫描剂量值,对所述当前待扫描位置进行扫描得到初步扫描图像;
第二降噪子模块,用于根据所述待降噪对象,对初步扫描图像进行降噪处理,以得到符合所述目标图像质量标准的扫描图像。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410445744.4A CN104644203B (zh) | 2014-09-02 | 2014-09-02 | 一种剂量调制扫描方法和装置 |
US14/843,754 US9836825B2 (en) | 2014-09-02 | 2015-09-02 | Dose modulated X-ray scanning |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410445744.4A CN104644203B (zh) | 2014-09-02 | 2014-09-02 | 一种剂量调制扫描方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104644203A CN104644203A (zh) | 2015-05-27 |
CN104644203B true CN104644203B (zh) | 2018-01-23 |
Family
ID=53236206
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410445744.4A Active CN104644203B (zh) | 2014-09-02 | 2014-09-02 | 一种剂量调制扫描方法和装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US9836825B2 (zh) |
CN (1) | CN104644203B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104644203B (zh) * | 2014-09-02 | 2018-01-23 | 沈阳东软医疗系统有限公司 | 一种剂量调制扫描方法和装置 |
CN104287768A (zh) * | 2014-09-30 | 2015-01-21 | 沈阳东软医疗系统有限公司 | 一种ct扫描剂量控制方法及系统 |
CN104605881A (zh) * | 2014-12-31 | 2015-05-13 | 沈阳东软医疗系统有限公司 | 一种参数优化方法及医疗设备 |
WO2016203282A1 (en) | 2015-06-18 | 2016-12-22 | The Nielsen Company (Us), Llc | Methods and apparatus to capture photographs using mobile devices |
CN106994021B (zh) * | 2016-01-22 | 2022-10-14 | 通用电气公司 | 一种计算ct影像上的噪声的方法及装置 |
EP3384847B1 (de) | 2017-04-05 | 2019-06-12 | Siemens Healthcare GmbH | Bestimmen eines referenzdosisparameters einer computertomographiebildgebung |
CN107978003B (zh) * | 2017-10-30 | 2021-03-30 | 东软医疗系统股份有限公司 | 一种ct图像金属伪影的处理方法及装置 |
WO2019090533A1 (en) | 2017-11-08 | 2019-05-16 | Shenzhen United Imaging Healthcare Co., Ltd. | System and method for diagnositic and treatment |
CN108065952B (zh) * | 2018-01-03 | 2021-08-13 | 东软医疗系统股份有限公司 | 一种剂量调制方法及装置 |
CN112089431B (zh) * | 2020-09-15 | 2023-05-30 | 上海联影医疗科技股份有限公司 | 剂量调制的优化方法、装置、电子设备和介质 |
Family Cites Families (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5379333A (en) * | 1993-11-19 | 1995-01-03 | General Electric Company | Variable dose application by modulation of x-ray tube current during CT scanning |
US5400378A (en) * | 1993-11-19 | 1995-03-21 | General Electric Company | Dynamic dose control in multi-slice CT scan |
US5450462A (en) * | 1993-11-19 | 1995-09-12 | General Electric Company | Modulation of x-ray tube current during CT scanning with modulation limit |
US5485494A (en) * | 1994-08-03 | 1996-01-16 | General Electric Company | Modulation of X-ray tube current during CT scanning |
DE19957083B4 (de) * | 1999-11-28 | 2004-11-18 | Siemens Ag | Verfahren zur Untersuchung eines eine periodische Bewegung ausführenden Körperbereichs |
AU2001288024A1 (en) * | 2001-09-05 | 2003-03-18 | Koninklijke Philips Electronics N.V. | Dose control in ct-images |
JP2004173924A (ja) * | 2002-11-27 | 2004-06-24 | Ge Medical Systems Global Technology Co Llc | X線制御方法およびx線画像撮影装置 |
CN1615797A (zh) * | 2004-12-08 | 2005-05-18 | 东软飞利浦医疗设备系统有限责任公司 | 在ct系统中保证z轴方向图像噪声一致性的方法 |
WO2007138979A1 (ja) * | 2006-05-25 | 2007-12-06 | Hitachi Medical Corporation | X線ct装置 |
US8031831B2 (en) * | 2009-05-28 | 2011-10-04 | Kabushiki Kaisha Toshiba | Voltage and or current modulation in dual energy computed tomography |
US8699658B2 (en) * | 2009-07-15 | 2014-04-15 | Mayo Foundation For Medical Education And Research | System and method for automatic tube potential selection for dose reduction in medical imaging |
CN107252318A (zh) * | 2009-10-22 | 2017-10-17 | 皇家飞利浦电子股份有限公司 | 采集协议评价设备 |
CN103180875B (zh) * | 2010-10-27 | 2018-04-13 | 皇家飞利浦电子股份有限公司 | 低剂量ct降噪 |
AU2012315530A1 (en) * | 2011-09-30 | 2014-04-17 | Cincinnati Children's Hospital Medical Center | Method for consistent and verifiable optimization of computed tomography (CT) radiation dose |
JP5774447B2 (ja) | 2011-10-28 | 2015-09-09 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置および被曝線量計算方法並びにプログラム |
US9492136B2 (en) | 2012-01-06 | 2016-11-15 | Indiana University Research And Technology Corporation | Method and apparatus that automates tube current and voltage selection for CT scans |
US20130202079A1 (en) * | 2012-02-07 | 2013-08-08 | Lifeng Yu | System and Method for Controlling Radiation Dose for Radiological Applications |
JP6265992B2 (ja) * | 2012-08-06 | 2018-01-24 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 画像ノイズの低減及び/又は画像解像度の向上 |
US9111337B2 (en) * | 2012-10-12 | 2015-08-18 | Mayo Foundation For Medical Education And Research | Low dose cardiac CT imaging with time-adaptive filtration |
CN103565460B (zh) * | 2013-09-26 | 2016-06-29 | 沈阳东软医疗系统有限公司 | 一种降低扫描剂量的扫描方法和装置 |
CN104644203B (zh) * | 2014-09-02 | 2018-01-23 | 沈阳东软医疗系统有限公司 | 一种剂量调制扫描方法和装置 |
-
2014
- 2014-09-02 CN CN201410445744.4A patent/CN104644203B/zh active Active
-
2015
- 2015-09-02 US US14/843,754 patent/US9836825B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
US20160063686A1 (en) | 2016-03-03 |
US9836825B2 (en) | 2017-12-05 |
CN104644203A (zh) | 2015-05-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104644203B (zh) | 一种剂量调制扫描方法和装置 | |
CN105534543B (zh) | 一种ct扫描参数优化方法及装置 | |
CN106355570B (zh) | 一种结合深度特征的双目立体视觉匹配方法 | |
CN101984916B (zh) | 基于数字图像处理技术的血管管径测量方法 | |
CN104574327B (zh) | 一种乳腺图像处理方法和装置 | |
CN102100562B (zh) | 一种x射线ct机的在线剂量调制方法 | |
CN101727661B (zh) | 医用图像处理装置 | |
CN104318536B (zh) | Ct图像的校正方法及装置 | |
CN104346799B (zh) | 一种ct图像中脊髓的提取方法 | |
US20060198499A1 (en) | Device and method for adapting the recording parameters of a radiograph | |
CN106880906A (zh) | 一种用于放射治疗中患者摆位的实时验证装置及方法 | |
CN105389815A (zh) | 一种乳腺图像配准方法及装置 | |
CN107702663A (zh) | 一种基于具有标志点的旋转平台的点云配准方法 | |
US9980682B2 (en) | Curved movable beam stop array and CBCT comprising thereof | |
WO2009118021A1 (en) | Radiotherapeutic treatment planning apparatus and methods | |
CN106875354A (zh) | 图像去噪方法、装置及设备 | |
CN103247046B (zh) | 一种放射治疗计划中靶区自动勾画的方法和装置 | |
JP2004512070A (ja) | マスクされた円錐ビーム投影データの簡易化された処理によって発生したイメージの誤差を識別して補正するための方法および装置 | |
CN106530277B (zh) | 一种基于小波方向相关系数的图像融合方法 | |
CN102170521A (zh) | 一种数码相机基于非均匀采样窗口的自动对焦方法 | |
CN104586421B (zh) | 一种消除摆位偏心的方法和装置 | |
Chen et al. | Diaphragm motion quantification in megavoltage cone‐beam CT projection images | |
Malian et al. | Development of a robust photogrammetric metrology system for monitoring the healing of bedsores | |
Preuhs et al. | Over-exposure correction in CT using optimization-based multiple cylinder fitting | |
CN109087357A (zh) | 扫描定位方法、装置、计算机设备及计算机可读存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address | ||
CP03 | Change of name, title or address |
Address after: 110179 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Patentee after: DongSoft Medical System Co., Ltd. Address before: Hunnan New Century Road 110179 Shenyang city of Liaoning Province, No. 16 Patentee before: Dongruan Medical Systems Co., Ltd., Shenyang |