CN111839562B - 射线数据的处理方法、装置、存储介质和电子设备 - Google Patents
射线数据的处理方法、装置、存储介质和电子设备 Download PDFInfo
- Publication number
- CN111839562B CN111839562B CN202010568969.4A CN202010568969A CN111839562B CN 111839562 B CN111839562 B CN 111839562B CN 202010568969 A CN202010568969 A CN 202010568969A CN 111839562 B CN111839562 B CN 111839562B
- Authority
- CN
- China
- Prior art keywords
- data
- oscillation
- initial
- ray
- ray data
- 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
- 238000003672 processing method Methods 0.000 title claims description 12
- 230000010355 oscillation Effects 0.000 claims abstract description 347
- 238000012545 processing Methods 0.000 claims abstract description 42
- 238000000034 method Methods 0.000 claims abstract description 39
- 238000005070 sampling Methods 0.000 claims abstract description 23
- 238000009499 grossing Methods 0.000 claims abstract description 21
- 230000008859 change Effects 0.000 claims description 29
- 238000001228 spectrum Methods 0.000 claims description 21
- 238000004590 computer program Methods 0.000 claims description 6
- 230000000149 penetrating effect Effects 0.000 claims description 6
- 230000009514 concussion Effects 0.000 claims description 5
- 230000005855 radiation Effects 0.000 description 32
- 238000010586 diagram Methods 0.000 description 10
- 238000004891 communication Methods 0.000 description 7
- 230000005236 sound signal Effects 0.000 description 4
- 230000000737 periodic effect Effects 0.000 description 3
- 238000010894 electron beam technology Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Images
Classifications
-
- 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/06—Diaphragms
-
- 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
- A61B6/035—Mechanical aspects of CT
-
- 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- 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
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- 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]
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Pathology (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- High Energy & Nuclear Physics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Theoretical Computer Science (AREA)
- Animal Behavior & Ethology (AREA)
- Optics & Photonics (AREA)
- Public Health (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Pure & Applied Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本公开涉及一种射线数据的处理方法、装置、存储介质和电子设备,应用于CT设备,CT设备包括射线发射器、准直器和探测器,该方法包括:通过至少两个时间窗,对探测器采集到的初始射线数据进行平滑处理,以获取每个时间窗对应的平滑射线数据,每个时间窗的窗口长度不相同,初始射线数据包括指定时长内按照采样频率采集的射线经过待测物体的衰减值,根据初始射线数据和任一个时间窗对应的平滑射线数据,确定初始射线数据中包括的初始震荡数据,初始震荡数据能够反映由于射线被准直器遮挡所产生的震荡,根据每个时间窗对应的平滑射线数据和初始震荡数据,对初始射线数据进行去除震荡处理,以得到目标射线数据。能够提高射线数据的准确性。
Description
技术领域
本公开涉及电子信息技术领域,具体地,涉及一种射线数据的处理方法、装置、存储介质和电子设备。
背景技术
随着图像处理技术的不断发展,CT(英文:Computed Tomography,中文:电子计算机断层扫描)设备在医疗领域得到了广泛应用。CT设备是利用精确准直的射线(例如:X射线、γ射线、超声波等),与灵敏度极高的探测器一同围绕人体的某一部位作断面扫描,以得到该部位的CT图像。CT设备具有扫描时间快,图像清晰等特点,能够辅助医生观察该部位的状况。CT设备中的射线发射器通常为旋转阳极X线球管,旋转阳极X线球管的转轴在高速旋转过程中,由于轴承间隙和重力作用的影响,会发生周期性的位置偏转,位于转轴上的阳极靶面也会发生偏转,导致X射线偏离理想位置,而CT设备中的准直器与旋转阳极X线球管之间的位置是固定的,因此X射线会被准直器周期性地遮挡,导致CT设备中探测器采集到的X射线数据不准确。
发明内容
本公开的目的是提供一种射线数据的处理方法、装置、存储介质和电子设备,用以解决现有技术中存在的探测器采集到的X射线数据不准确的问题。
为了实现上述目的,根据本公开实施例的第一方面,提供一种射线数据的处理方法,应用于CT设备,所述CT设备包括射线发射器、准直器和探测器,所述射线发射器用于发出射线,射线通过所述准直器射向待测物体,所述探测器用于采集透过所述待测物体的射线;所述方法包括:
通过至少两个时间窗,对所述探测器采集到的初始射线数据进行平滑处理,以获取每个所述时间窗对应的平滑射线数据,每个所述时间窗的窗口长度不相同,所述初始射线数据包括指定时长内按照采样频率采集的射线经过所述待测物体的衰减值;
根据所述初始射线数据和任一个所述时间窗对应的所述平滑射线数据,确定所述初始射线数据中包括的初始震荡数据,所述初始震荡数据能够反映由于射线被所述准直器遮挡所产生的震荡;
根据每个所述时间窗对应的所述平滑射线数据和所述初始震荡数据,对所述初始射线数据进行去除震荡处理,以得到目标射线数据。
可选地,所述根据每个所述时间窗对应的所述平滑射线数据和所述初始震荡数据,对所述初始射线数据进行去除震荡处理,以得到目标射线数据,包括:
根据每个所述时间窗对应的所述平滑射线数据,确定分界区间,以按照所述分界区间将所述初始震荡数据分为平缓震荡数据和突变震荡数据,所述平缓震荡数据的峰峰值小于所述突变震荡数据的峰峰值;
根据所述平缓震荡数据,替换所述初始震荡数据中的所述突变震荡数据,以得到目标震荡数据;
将所述初始射线数据与所述目标震荡数据的差,作为所述目标射线数据。
可选地,所述根据每个所述时间窗对应的所述平滑射线数据,确定分界区间,以按照所述分界区间将所述初始震荡数据分为平缓震荡数据和突变震荡数据,包括:
根据所述至少两个时间窗中每两个所述时间窗对应的所述平滑射线数据的差值,确定所述分界区间,所述分界区间包括平稳区间和突变区间;
将所述初始震荡数据中在所述突变区间内采集的数据作为所述突变震荡数据,将所述初始震荡数据中在所述平稳区间内采集的数据作为所述平缓震荡数据。
可选地,所述根据所述平缓震荡数据,替换所述初始震荡数据中的所述突变震荡数据,以得到目标震荡数据,包括:
将所述平缓震荡数据转换为平缓震荡频谱;
根据所述平缓震荡频谱,确定由于射线被所述准直器遮挡所产生的震荡对应的震荡周期;
从所述平缓震荡数据取出时长为所述震荡周期的替换数据,并利用所述替换数据替换所述初始震荡数据中的所述突变震荡数据,以得到所述目标震荡数据。
可选地,所述根据所述平缓震荡频谱,确定由于射线被所述准直器遮挡所产生的震荡对应的震荡周期,包括:
确定所述平缓震荡频谱中峰值的频点;
根据所述峰值的频点、所述采样频率和所述指定时长,确定所述震荡周期。
可选地,所述根据所述初始射线数据和任一个所述时间窗对应的所述平滑射线数据,确定所述初始射线数据中包括的初始震荡数据,包括:
将所述初始射线数据与任一个所述时间窗对应的所述平滑射线数据的差,作为所述初始震荡数据。
可选地,在所述根据每个所述时间窗对应的所述平滑射线数据和所述初始震荡数据,对所述初始射线数据进行去除震荡处理,以得到目标射线数据之后,所述方法还包括:
根据所述目标射线数据,生成所述待测物体的CT图像。
根据本公开实施例的第二方面,提供一种射线数据的处理装置,应用于CT设备,所述CT设备包括射线发射器、准直器和探测器,所述射线发射器用于发出射线,射线通过所述准直器射向待测物体,所述探测器用于采集透过所述待测物体的射线;所述装置包括:
平滑处理模块,用于通过至少两个时间窗,对所述探测器采集到的初始射线数据进行平滑处理,以获取每个所述时间窗对应的平滑射线数据,每个所述时间窗的窗口长度不相同,所述初始射线数据包括指定时长内按照采样频率采集的射线经过所述待测物体的衰减值;
确定模块,用于根据所述初始射线数据和任一个所述时间窗对应的所述平滑射线数据,确定所述初始射线数据中包括的初始震荡数据,所述初始震荡数据能够反映由于射线被所述准直器遮挡所产生的震荡;
去除震荡模块,用于根据每个所述时间窗对应的所述平滑射线数据和所述初始震荡数据,对所述初始射线数据进行去除震荡处理,以得到目标射线数据。
可选地,所述去除震荡模块包括:
分界确定子模块,用于根据每个所述时间窗对应的所述平滑射线数据,确定分界区间,以按照所述分界区间将所述初始震荡数据分为平缓震荡数据和突变震荡数据,所述平缓震荡数据的峰峰值小于所述突变震荡数据的峰峰值;
替换子模块,用于根据所述平缓震荡数据,替换所述初始震荡数据中的所述突变震荡数据,以得到目标震荡数据;
去除子模块,用于将所述初始射线数据与所述目标震荡数据的差,作为所述目标射线数据。
可选地,所述分界确定子模块用于:
根据所述至少两个时间窗中每两个所述时间窗对应的所述平滑射线数据的差值,确定所述分界区间,所述分界区间包括平稳区间和突变区间;
将所述初始震荡数据中在所述突变区间内采集的数据作为所述突变震荡数据,将所述初始震荡数据中在所述平稳区间内采集的数据作为所述平缓震荡数据。
可选地,所述替换子模块用于:
将所述平缓震荡数据转换为平缓震荡频谱;
根据所述平缓震荡频谱,确定由于射线被所述准直器遮挡所产生的震荡对应的震荡周期;
从所述平缓震荡数据取出时长为所述震荡周期的替换数据,并利用所述替换数据替换所述初始震荡数据中的所述突变震荡数据,以得到所述目标震荡数据。
可选地,所述替换子模块用于:
确定所述平缓震荡频谱中峰值的频点;
根据所述峰值的频点、所述采样频率和所述指定时长,确定所述震荡周期。
可选地,所述确定模块用于:
将所述初始射线数据与任一个所述时间窗对应的所述平滑射线数据的差,作为所述初始震荡数据。
可选地,所述装置还包括:
生成模块,用于在所述根据每个所述时间窗对应的所述平滑射线数据和所述初始震荡数据,对所述初始射线数据进行去除震荡处理,以得到目标射线数据之后,根据所述目标射线数据,生成所述待测物体的CT图像。
根据本公开实施例的第三方面,提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现本公开实施例的第一方面所述方法的步骤。
根据本公开实施例的第四方面,提供一种电子设备,包括:
存储器,其上存储有计算机程序;
处理器,用于执行所述存储器中的所述计算机程序,以实现本公开实施例的第一方面所述方法的步骤。
通过上述技术方案,本公开中探测器采集初始射线数据,其中,初始射线数据包括指定时长内按照采样频率采集的射线经过待测物体的衰减值。首先通过至少两个窗口长度不同的时间窗对初始射线数据进行平滑处理,以得到每个时间窗对应的平滑射线数据,之后根据初始射线数据和任一个时间窗对应的平滑射线数据,确定初始射线数据中包括的,能够反映由于射线被准直器遮挡所产生的震荡的初始震荡数据,最后根据每个时间窗对应的平滑射线数据和初始震荡数据,对初始射线数据进行去除震荡处理,以得到目标射线数据。本公开通过多个时间窗对采集到的初始射线数据进行平滑处理,从而得到能够反映震荡的初始震荡数据,最后对初始射线数据进行去除震荡处理,以去除由于射线被准直器周期性地遮挡所产生的震荡,提高了射线数据的准确性。
本公开的其他特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本公开的进一步理解,并且构成说明书的一部分,与下面的具体实施方式一起用于解释本公开,但并不构成对本公开的限制。在附图中:
图1是根据一示例性实施例示出的一种CT设备的示意图;
图2是根据一示例性实施例示出的一种射线数据的处理方法的流程图;
图3是根据一示例性实施例示出的一种初始射线数据;
图4是根据一示例性实施例示出的一种初始射线数据和平滑射线数据;
图5是根据一示例性实施例示出的一种初始震荡数据;
图6是根据一示例性实施例示出的另一种射线数据的处理方法的流程图;
图7是根据一示例性实施例示出的一种两个时间窗对应的平滑射线数据的差值;
图8是根据一示例性实施例示出的一种平缓震荡频谱;
图9是根据一示例性实施例示出的一种目标震荡频谱;
图10是根据一示例性实施例示出的一种目标射线频谱;
图11是根据一示例性实施例示出的另一种射线数据的处理方法的流程图;
图12是根据一示例性实施例示出的CT图像;
图13是根据一示例性实施例示出的一种射线数据的处理装置的框图;
图14是根据一示例性实施例示出的另一种射线数据的处理装置的框图;
图15是根据一示例性实施例示出的另一种射线数据的处理装置的框图;
图16是根据一示例性实施例示出的一种电子设备的框图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本公开相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本公开的一些方面相一致的装置和方法的例子。
本公开通过多个时间窗对采集到的初始射线数据进行平滑处理,以从初始射线数据中提取出能够反映震荡的初始震荡数据,再根据初始震荡数据和经过平滑处理的数据,对初始射线数据进行去除震荡处理,以去除由于射线被准直器周期性地遮挡所产生的震荡,能够提高探测器采集到的射线数据的准确性。
在介绍本公开提供的射线数据的处理方法、装置、存储介质和电子设备之前,首先对本公开各个实施例所涉及的应用场景进行介绍。该应用场景为CT设备,如图1中的(a)所示,CT设备包括射线发射器、准直器和探测器(射线发射器、准直器和探测器共同组成扫描架),射线发射器用于发出射线,射线通过准直器射向待测物体,探测器用于采集透过待测物体的射线。需要说明的是,本公开实施例中所提及的射线可以为X射线、γ射线或者超声波等,本公开对此不作具体限定。CT设备还可以包括射线高压发生器、ADC(英文:Analog/Digital Converter,中文:模拟/数字转换器)、和控制单元,其中射线高压发生器用于使射线发射器发出射线,ADC用于将探测器采集到的模拟信号转换为数字信号(可以理解为下文所述的初始射线数据),并发送给控制单元。控制单元用于对数字信号进行处理,并根据数字信号生成CT图像。进一步的,CT设备还可以包括图像显示、存储、记录模块(例如:显示器、存储器、工作站等)。本公开所提供的射线数据的处理方法的执行主体,可以理解为上述控制单元。控制单元例如可以是MCU(英文:Microcontroller Unit,中文:微控制单元)、ECU(英文:Electronic Control Unit,中文:电子控制单元),也可以是终端设备(例如:笔记本电脑、平板电脑或者台式计算机)或者服务器。
在射线发射器为旋转阳极X线球管时,阳极端由阳极靶面、定子、转子、转轴和轴承(未在图中标出)组成,如图1中的(b)所示。阳极靶面为圆盘形状,圆盘中心固定在位于转子中心的转轴一端,定子线圈产生磁场引起转子转动会带动阳极靶面高速旋转,可以使阴极灯丝发射的电子束轰击阳极靶面产生的热量分布在连续运动的阳极靶面的圆环区域内,从而降低阳极靶面的温度。然而转子在高速旋转过程中,由于轴承间隙和重力作用的影响,会发生周期性的位置偏转,位于转轴上的阳极靶面也会发生偏转,导致X射线偏离理想位置。而CT设备中的准直器与旋转阳极X线球管之间的位置是固定的,如图1中的(c)所示,因此X射线会被准直器周期性地遮挡,即X射线在z轴方向(即电子束发射的方向)上被准直器周期性地遮挡,导致探测器采集到的X射线数据中存在震荡。本公开所提供的射线数据的处理方法就是用于去除X射线数据中存在的震荡。
图2是根据一示例性实施例示出的一种射线数据的处理方法的流程图,如图2所示,应用于上述CT设备。该方法包括:
步骤101,通过至少两个时间窗,对探测器采集到的初始射线数据进行平滑处理,以获取每个时间窗对应的平滑射线数据,每个时间窗的窗口长度不相同,初始射线数据包括指定时长内按照采样频率采集的射线经过待测物体的衰减值。
举例来说,射线发射器发出射线,射线通过准直器射向待测物体(例如可以是人体的某个部位),探测器采集透过待测物体的射线,由于待测物体会吸收部分射线,因此探测器采集到的射线,与射线发射器发出射线不同,可以根据探测器采集到的射线经过待测物体的衰减值来确定待测物体的密度,从而生成CT图像。探测器可以按照指定的采样频率在指定时长内采集射线经过待测物体的衰减值,即初始射线数据。以采样频率为1000Hz(对应的采样周期为1ms),指定时长为6s来举例,那么初始射线数据中包括了6s中每1ms探测器采集的射线经过待测物体的衰减值(即初始射线数据中包括了6000个采样时刻采集的数据),如图3所示,纵轴表示射线射出后经过待测物体的衰减值,即被待测物体吸收的光子数,可以根据探测器接收到的光子数,和射线发射器发射出的光子数求差得到,横轴表示时间,单位为ms。
可以预先设置至少两个时间窗,每个时间窗的窗口长度(也称为窗宽)均不相同。分别利用每个时间窗对探测器采集到的初始射线数据进行平滑处理,得到每个时间窗对应的平滑射线数据。每个时间窗对应的平滑射线数据,能够在该时间窗对应的时间尺度(即窗口长度)上反映初始射线数据中包括的真实数据(即去除震荡后的数据)。例如,可以设置有第一时间窗,对应的窗口长度为50ms,和第二时间窗,对应的窗口长度为100ms。然后利用第一时间窗对初始射线数据进行平滑处理,得到第一平滑射线数据,利用第二时间窗对初始射线数据进行平滑处理,得到第二平滑射线数据。
步骤102,根据初始射线数据和任一个时间窗对应的平滑射线数据,确定初始射线数据中包括的初始震荡数据,初始震荡数据能够反映由于射线被准直器遮挡所产生的震荡。
示例的,可以从多个时间窗中,选取任一个时间窗,然后根据该时间窗对应的平滑射线数据,从初始射线数据中提取能够反映由于射线被准直器遮挡所产生的震荡对应的初始震荡数据。在一种可实现的方式中:可以将初始射线数据与任一个时间窗对应的平滑射线数据的差,作为初始震荡数据。可以理解为,任一个时间窗对应的平滑射线数据,能够在该时间窗对应的时间尺度上反映初始射线数据中包括的真实数据,那么从初始射线数据中减去平滑射线数据,可以得到能够反映震荡的初始震荡数据。如图4所示,其中,与图3中相同的曲线为初始震荡数据,与初始震荡数据整体趋势相同,但更加平稳的曲线为第一时间窗对应的第一平滑射线数据,那么初始震荡数据如图5所示。在另一种可实现的方式中:可以先将初始射线数据与任一个时间窗对应的平滑射线数据求差,然后利用预设的第一衰减阈值对求差的结果进行筛选,以去除超过第一衰减阈值的数据。由于,射线会被准直器周期性地遮挡,初始射线数据中包括的震荡也应当是周期性的,并且射线经过待测物体的衰减值应当在指定的范围内,因此可以利用第一衰减阈值排除掉由于其他外界噪声产生的干扰,提高初始震荡数据的有效性。
步骤103,根据每个时间窗对应的平滑射线数据和初始震荡数据,对初始射线数据进行去除震荡处理,以得到目标射线数据。
示例的,在确定初始震荡数据之后,可以根据每个时间窗对应的平滑射线数据,进一步对初始震荡数据进行处理,从而更加准确地反映由于射线被准直器遮挡所产生的震荡。由于,射线会被准直器周期性地遮挡,初始射线数据中包括的震荡也是周期性的,如果初始震荡数据中存在震荡剧烈的部分,可能是由于射线经过待测物体的衰减值本身产生了震荡。例如,在开始使用CT设备对待测物体进行扫描时,射线发射器发射出的射线并不稳定,相应的射线经过待测物体的衰减值也会出现震荡。在一定时间之后,射线才趋于稳定,相应的,探测器采集到的射线经过待测物体的衰减值才会趋于稳定。因此,可以根据每个时间窗对应的平滑射线数据,将初始震荡数据分为震荡剧烈的数据,和震荡平缓的数据。之后,可以用震荡平缓的数据,替换震荡剧烈的数据,以得到震荡平缓的目标震荡数据。例如,可以根据震荡剧烈的数据对应的时间长度,从震荡平缓的数据中,取出该时间长度的数据,以替换震荡剧烈的数据,得到目标震荡数据。最后可以将目标震荡数据从初始射线数据中去除,以完成去除震荡处理,得到的数据即为目标射线数据。目标射线数据相比于初始射线数据,去除了由于射线被准直器遮挡所产生的震荡,能够准确反映射线经过待测物体的衰减值,提高了射线数据的准确性。
综上所述,本公开中探测器采集初始射线数据,其中,初始射线数据包括指定时长内按照采样频率采集的射线经过待测物体的衰减值。首先通过至少两个窗口长度不同的时间窗对初始射线数据进行平滑处理,以得到每个时间窗对应的平滑射线数据,之后根据初始射线数据和任一个时间窗对应的平滑射线数据,确定初始射线数据中包括的,能够反映由于射线被准直器遮挡所产生的震荡的初始震荡数据,最后根据每个时间窗对应的平滑射线数据和初始震荡数据,对初始射线数据进行去除震荡处理,以得到目标射线数据。本公开通过多个时间窗对采集到的初始射线数据进行平滑处理,从而得到能够反映震荡的初始震荡数据,最后对初始射线数据进行去除震荡处理,以去除由于射线被准直器周期性地遮挡所产生的震荡,提高了射线数据的准确性。
图6是根据一示例性实施例示出的另一种射线数据的处理方法的流程图,如图6所示,步骤103的实现方式可以包括以下步骤:
步骤1031,根据每个时间窗对应的平滑射线数据,确定分界区间,以按照分界区间将初始震荡数据分为平缓震荡数据和突变震荡数据,平缓震荡数据的峰峰值小于突变震荡数据的峰峰值。
举例来说,可以根据至少两个时间窗中每两个时间窗对应的平滑射线数据的差值,确定分界区间。可以理解为CT设备对待测物体进行扫描时,如果射线发射器刚开始发射出射线,射线经过待测物体的衰减值可能会出现震荡,或者由于射线穿过了形状不规则、密度不统一的物体(例如:人体佩戴的首饰、皮带等与人体器官密度差别较大的物体),射线经过待测物体的衰减值也可能出现震荡。而射线会被准直器周期性地遮挡,相应产生的震荡也是周期性的。因此,初始震荡数据中,存在一部分即包括了由于射线被准直器遮挡所产生的震荡,又包括了射线经过待测物体的衰减值本身产生的震荡的数据,即突变震荡数据,另一部分只包括由于射线被准直器遮挡所述产生的震荡的数据(此时射线经过待测物体的衰减值已经稳定),即平缓震荡数据。其中,突变震荡数据的峰峰值要大于平缓震荡数据的峰峰值。
分界区间用于指示初始震荡数据中哪些数据属于突变震荡数据,哪些数据属于平缓震荡数据。因此可以将分界区间分为两种时间区间:平稳区间和突变区间。其中,可能存在多个平稳区间,和多个突变区间,平稳区间和突变区间间隔排列,每个平稳区间包括了该区间的起始时刻和结束时刻,同样的,每个突变区间也包括了该区间的起始时刻和结束时刻。那么在确定分界区间之后,可以将初始震荡数据中在突变区间内采集的数据作为突变震荡数据,将初始震荡数据中在平稳区间内采集的数据作为平缓震荡数据。
例如,若步骤101中采用了三个以上的时间窗对初始射线数据进行平滑处理,那么可以先将至少两个时间窗两两分组,然后将每两个时间窗对应的平滑射线数据求差,然后通过阈值法对差值进行划分,求得每两个时间窗对应的分界区间。之后将每两个时间窗对应的分界区间求平均,作为用于划分初始震荡数据的分界区间。例如,存在第一时间窗、第二时间窗、第三时间窗,那么可以分别计算第一时间窗和第二时间窗对应的分界区间,第二时间窗和第三时间窗对应的分界区间,第一时间窗和第三时间窗对应的分界区间,然后将三个分界区间的平均值,作为用于划分初始震荡数据的分界区间。若步骤101中采用了两个时间窗对初始射线数据进行平滑处理,那么可以直接将这两个时间窗对应的平滑射线数据求差,通过阈值法对差值进行划分,求得分界区间。
下面,以对两个时间窗对应的平滑射线数据求差,通过阈值法对差值进行划分,以确定分界区间来进行具体说明。以第一时间窗对应的窗口长度为50ms,第二时间窗对应的窗口长度为100ms来举例,用第一时间窗对应的平滑数据减去第二时间窗对应的平滑数据,得到的差值如图7所示,可以看出经过第一时间窗对应的平滑数据与第二时间窗对应的平滑数据求差,可以将数据搬移到0附近,便于差值中的每个数据和预设的第二衰减阈值(例如可以是50)进行比较。之后,遍历差值中的每个数据,与第二衰减阈值进行比较,可以确定在80ms至6000ms之间采集的数据的绝对值均小于第二衰减阈值,在0至80ms之间采集的数据中存在多个数据的绝对值大于或等于第二衰减阈值,那么可以将0至80ms作为突变区间,将80ms至6000ms作为平稳区间。以图5所示初始震荡数据为例,可以将0至80ms内采集的数据作为突变震荡数据,将80ms至6000ms内采集的数据作为平缓震荡数据。
步骤1032,根据平缓震荡数据,替换初始震荡数据中的突变震荡数据,以得到目标震荡数据。
步骤1033,将初始射线数据与目标震荡数据的差,作为目标射线数据。
示例的,在将初始震荡数据划分为突变震荡数据和平缓震荡数据之后,可以先确定突变震荡数据对应的时间长度,然后从平缓震荡数据中,取出与突变震荡数据对应的时间长度相等的一部分数据,以替换突变震荡数据,再将被替换后的突变震荡数据(即平缓震荡数据中的一部分数据)和平缓震荡数据按照时间顺序拼接起来,得到目标震荡数据,这样,相比于初始震荡数据,目标震荡数据去除了由于射线不稳定,导致射线经过待测物体的衰减值产生的震荡,能够更加准确地反映射线被准直器周期性地遮挡所产生的震荡。最后,可以将初始射线数据减去目标震荡数据,从初始射线数据中去除由于射线被准直器遮挡所产生的震荡,得到目标射线数据,能够准确反映射线经过待测物体的衰减值,提高了射线数据的准确性。
需要说明的是,在存在多个平稳区间和多个突变区间的场景中,根据多个平稳区间和多个突变区间,可以将初始震荡数据分为多组突变震荡数据和多组平缓震荡数据,突变震荡数据和平缓震荡数据间隔排列。那么在执行步骤1032时,针对每一组突变震荡数据进行替换的过程,可以根据与该组突变震荡数据相邻的一组平缓震荡数据,来完成替换。也可以根据多个平稳区间中,时间长度最长的平稳区间对应的平缓震荡数据,来完成替换,本公开对此不作具体限定。
在一种实现方式中,步骤1032的实现方式可以包括:
步骤3)将平缓震荡数据转换为平缓震荡频谱。
示例的,以图5所示初始震荡数据中80ms之后采集的数据作为平缓震荡数据为例,可以先对平缓震荡数据进行傅里叶变换,得到平缓震荡频谱,如图8所示。
步骤4)根据平缓震荡频谱,确定由于射线被准直器遮挡所产生的震荡对应的震荡周期。
之后,在平缓震荡频谱中找到峰值的频点,再根据峰值的频点、采样频率和指定时长,确定震荡周期。例如,采样频率SampleFreq=1/IntegrationTime,其中IntegrationTime为采样周期。峰值的频点为MaxIndex,其中,MaxIndex的取值范围为0到N/2,其中N为采样数量(即指定时长*SampleFreq)。那么,转子发生位置偏转的频率为:WobblingFreq=SampleFreq*MaxIndex/N,相应的,由转子发生位置偏转造成射线被准直器遮挡所产生的震荡,对应的震荡周期为:Circle=SampleFreq/WobblingFreq。
步骤5)从平缓震荡数据取出时长为震荡周期的替换数据,并利用替换数据替换初始震荡数据中的突变震荡数据,以得到目标震荡数据。
最后,为了保持震荡的周期性,可以从平缓震荡数据取出时长为震荡周期的替换数据,根据突变震荡数据对应的时间长度,利用替换数据替换初始震荡数据中的突变震荡数据。例如,震荡周期Circle为40ms,突变震荡数据对应的时间长度为80ms,那么可以将替换数据重复两次(即80ms/40ms),以替换突变震荡数据,然后将重复了两次的替换数据和平缓震荡数据按照时间顺序拼接起来,得到目标震荡数据,如图9所示。可以看出,相比于图5所示的初始震荡数据,图9中的目标震荡数据的周期性更明显,更能够体现由于射线被准直器周期性地遮挡所产生的震荡。相应的,将初始射线数据减去目标震荡数据,从初始射线数据中去除了由于射线被准直器遮挡所产生的震荡,得到目标射线数据,如图10所示,能够准确反映射线经过待测物体的衰减值,提高了射线数据的准确性。
图11是根据一示例性实施例示出的另一种射线数据的处理方法的流程图,如图11所示,在步骤103之后,该方法还包括:
步骤104,根据目标射线数据,生成待测物体的CT图像。
举例来说,在得到能够准确反映射线经过待测物体的衰减值的目标射线数据后,可以根据目标射线数据,生成待测物体的CT图像。例如,可以根据预设的射线经过待测物体的衰减值与图像的明暗程度的映射关系,将目标射线数据映射为CT图像。由于目标射线数据中,去除了由于射线被准直器遮挡所产生的震荡,因此能够有效去除CT图像中的横纹伪影。图12中的(a)是直接根据初始射线数据生成的CT图像,能够看出图像中存在大量的横纹伪影。图12中的(b)是根据目标射线数据生成的CT图像,相比于图12中的(a),能够有效去除横纹伪影。
综上所述,本公开中探测器采集初始射线数据,其中,初始射线数据包括指定时长内按照采样频率采集的射线经过待测物体的衰减值。首先通过至少两个窗口长度不同的时间窗对初始射线数据进行平滑处理,以得到每个时间窗对应的平滑射线数据,之后根据初始射线数据和任一个时间窗对应的平滑射线数据,确定初始射线数据中包括的,能够反映由于射线被准直器遮挡所产生的震荡的初始震荡数据,最后根据每个时间窗对应的平滑射线数据和初始震荡数据,对初始射线数据进行去除震荡处理,以得到目标射线数据。本公开通过多个时间窗对采集到的初始射线数据进行平滑处理,从而得到能够反映震荡的初始震荡数据,最后对初始射线数据进行去除震荡处理,以去除由于射线被准直器周期性地遮挡所产生的震荡,提高了射线数据的准确性。
图13是根据一示例性实施例示出的一种射线数据的处理装置的框图,如图13所示,该装置200应用于CT设备,CT设备包括射线发射器、准直器和探测器,射线发射器用于发出射线,射线通过准直器射向待测物体,探测器用于采集透过待测物体的射线。包括以下模块:
平滑处理模块201,用于通过至少两个时间窗,对探测器采集到的初始射线数据进行平滑处理,以获取每个时间窗对应的平滑射线数据,每个时间窗的窗口长度不相同,初始射线数据包括指定时长内按照采样频率采集的射线经过待测物体的衰减值。
确定模块202,用于根据初始射线数据和任一个时间窗对应的平滑射线数据,确定初始射线数据中包括的初始震荡数据,初始震荡数据能够反映由于射线被准直器遮挡所产生的震荡。
去除震荡模块203,用于根据每个时间窗对应的平滑射线数据和初始震荡数据,对初始射线数据进行去除震荡处理,以得到目标射线数据。
图14是根据一示例性实施例示出的另一种射线数据的处理装置的框图,如图14所示,去除震荡模块203可以包括:
分界确定子模块2031,用于根据每个时间窗对应的平滑射线数据,确定分界区间,以按照分界区间将初始震荡数据分为平缓震荡数据和突变震荡数据,平缓震荡数据的峰峰值小于突变震荡数据的峰峰值。
替换子模块2032,用于根据平缓震荡数据,替换初始震荡数据中的突变震荡数据,以得到目标震荡数据。
去除子模块2033,用于将初始射线数据与目标震荡数据的差,作为目标射线数据。
在一种实现方式中,分界确定子模块2031用于执行以下步骤:
步骤1)根据至少两个时间窗中每两个时间窗对应的平滑射线数据的差值,确定分界区间,分界区间包括平稳区间和突变区间。
步骤2)将初始震荡数据中在突变区间内采集的数据作为突变震荡数据,将初始震荡数据中在平稳区间内采集的数据作为平缓震荡数据。
在另一种实现方式中,替换子模块2032用于执行以下步骤:
步骤3)将平缓震荡数据转换为平缓震荡频谱。
步骤4)根据平缓震荡频谱,确定由于射线被准直器遮挡所产生的震荡对应的震荡周期。
步骤5)从平缓震荡数据取出时长为震荡周期的替换数据,并利用替换数据替换初始震荡数据中的突变震荡数据,以得到目标震荡数据。
具体的,替换子模块2032可以用于:
首先,确定平缓震荡频谱中峰值的频点。
之后,根据峰值的频点、采样频率和指定时长,确定震荡周期。
在另一种实现方式中,确定模块202可以用于:
将初始射线数据与任一个时间窗对应的平滑射线数据的差,作为初始震荡数据。
图15是根据一示例性实施例示出的另一种射线数据的处理装置的框图,如图15所示,该装置200还包括:
生成模块204,用于在根据每个时间窗对应的平滑射线数据和初始震荡数据,对初始射线数据进行去除震荡处理,以得到目标射线数据之后,根据目标射线数据,生成待测物体的CT图像。
关于上述实施例中的装置,其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述,此处将不做详细阐述说明。
综上所述,本公开中探测器采集初始射线数据,其中,初始射线数据包括指定时长内按照采样频率采集的射线经过待测物体的衰减值。首先通过至少两个窗口长度不同的时间窗对初始射线数据进行平滑处理,以得到每个时间窗对应的平滑射线数据,之后根据初始射线数据和任一个时间窗对应的平滑射线数据,确定初始射线数据中包括的,能够反映由于射线被准直器遮挡所产生的震荡的初始震荡数据,最后根据每个时间窗对应的平滑射线数据和初始震荡数据,对初始射线数据进行去除震荡处理,以得到目标射线数据。本公开通过多个时间窗对采集到的初始射线数据进行平滑处理,从而得到能够反映震荡的初始震荡数据,最后对初始射线数据进行去除震荡处理,以去除由于射线被准直器周期性地遮挡所产生的震荡,提高了射线数据的准确性。
图16是根据一示例性实施例示出的一种电子设备300的框图。如图16所示,该电子设备300可以包括:处理器301,存储器302。该电子设备300还可以包括多媒体组件303,输入/输出(I/O)接口304,以及通信组件305中的一者或多者。
其中,处理器301用于控制该电子设备300的整体操作,以完成上述的射线数据的处理方法中的全部或部分步骤。存储器302用于存储各种类型的数据以支持在该电子设备300的操作,这些数据例如可以包括用于在该电子设备300上操作的任何应用程序或方法的指令,以及应用程序相关的数据,例如联系人数据、收发的消息、图片、音频、视频等等。该存储器302可以由任何类型的易失性或非易失性存储设备或者它们的组合实现,例如静态随机存取存储器(Static Random Access Memory,简称SRAM),电可擦除可编程只读存储器(Electrically Erasable Programmable Read-Only Memory,简称EEPROM),可擦除可编程只读存储器(Erasable Programmable Read-Only Memory,简称EPROM),可编程只读存储器(Programmable Read-Only Memory,简称PROM),只读存储器(Read-Only Memory,简称ROM),磁存储器,快闪存储器,磁盘或光盘。多媒体组件303可以包括屏幕和音频组件。其中屏幕例如可以是触摸屏,音频组件用于输出和/或输入音频信号。例如,音频组件可以包括一个麦克风,麦克风用于接收外部音频信号。所接收的音频信号可以被进一步存储在存储器302或通过通信组件305发送。音频组件还包括至少一个扬声器,用于输出音频信号。I/O接口304为处理器301和其他接口模块之间提供接口,上述其他接口模块可以是键盘,鼠标,按钮等。这些按钮可以是虚拟按钮或者实体按钮。通信组件305用于该电子设备300与其他设备之间进行有线或无线通信。无线通信,例如Wi-Fi,蓝牙,近场通信(Near FieldCommunication,简称NFC),2G、3G、4G、NB-IOT、eMTC、或其他5G等等,或它们中的一种或几种的组合,在此不做限定。因此相应的该通信组件305可以包括:Wi-Fi模块,蓝牙模块,NFC模块等等。
在一示例性实施例中,电子设备300可以被一个或多个应用专用集成电路(Application Specific Integrated Circuit,简称ASIC)、数字信号处理器(DigitalSignal Processor,简称DSP)、数字信号处理设备(Digital Signal Processing Device,简称DSPD)、可编程逻辑器件(Programmable Logic Device,简称PLD)、现场可编程门阵列(Field Programmable Gate Array,简称FPGA)、控制器、微控制器、微处理器或其他电子元件实现,用于执行上述的射线数据的处理方法。
在另一示例性实施例中,还提供了一种包括程序指令的计算机可读存储介质,该程序指令被处理器执行时实现上述的射线数据的处理方法的步骤。例如,该计算机可读存储介质可以为上述包括程序指令的存储器302,上述程序指令可由电子设备300的处理器301执行以完成上述的射线数据的处理方法。
以上结合附图详细描述了本公开的优选实施方式,但是,本公开并不限于上述实施方式中的具体细节,在本公开的技术构思范围内,可以对本公开的技术方案进行多种简单变型,这些简单变型均属于本公开的保护范围。
另外需要说明的是,在上述具体实施方式中所描述的各个具体技术特征,在不矛盾的情况下,可以通过任何合适的方式进行组合,为了避免不必要的重复,本公开对各种可能的组合方式不再另行说明。
此外,本公开的各种不同的实施方式之间也可以进行任意组合,只要其不违背本公开的思想,其同样应当视为本公开所公开的内容。
Claims (9)
1.一种射线数据的处理方法,其特征在于,应用于CT设备,所述CT设备包括射线发射器、准直器和探测器,所述射线发射器用于发出射线,射线通过所述准直器射向待测物体,所述探测器用于采集透过所述待测物体的射线;所述方法包括:
通过至少两个时间窗,对所述探测器采集到的初始射线数据进行平滑处理,以获取每个所述时间窗对应的平滑射线数据,每个所述时间窗的窗口长度不相同,所述初始射线数据包括指定时长内按照采样频率采集的射线经过所述待测物体的衰减值;
根据所述初始射线数据和任一个所述时间窗对应的所述平滑射线数据,确定所述初始射线数据中包括的初始震荡数据,所述初始震荡数据能够反映由于射线被所述准直器遮挡所产生的震荡;
根据每个所述时间窗对应的所述平滑射线数据和所述初始震荡数据,对所述初始射线数据进行去除震荡处理,以得到目标射线数据;
所述根据每个所述时间窗对应的所述平滑射线数据和所述初始震荡数据,对所述初始射线数据进行去除震荡处理,以得到目标射线数据,包括:
根据每个所述时间窗对应的所述平滑射线数据,确定分界区间,以按照所述分界区间将所述初始震荡数据分为平缓震荡数据和突变震荡数据,所述平缓震荡数据的峰峰值小于所述突变震荡数据的峰峰值;
根据所述平缓震荡数据,替换所述初始震荡数据中的所述突变震荡数据,以得到目标震荡数据;
将所述初始射线数据与所述目标震荡数据的差,作为所述目标射线数据。
2.根据权利要求1所述的方法,其特征在于,所述根据每个所述时间窗对应的所述平滑射线数据,确定分界区间,以按照所述分界区间将所述初始震荡数据分为平缓震荡数据和突变震荡数据,包括:
根据所述至少两个时间窗中每两个所述时间窗对应的所述平滑射线数据的差值,确定所述分界区间,所述分界区间包括平稳区间和突变区间;
将所述初始震荡数据中在所述突变区间内采集的数据作为所述突变震荡数据,将所述初始震荡数据中在所述平稳区间内采集的数据作为所述平缓震荡数据。
3.根据权利要求1所述的方法,其特征在于,所述根据所述平缓震荡数据,替换所述初始震荡数据中的所述突变震荡数据,以得到目标震荡数据,包括:
将所述平缓震荡数据转换为平缓震荡频谱;
根据所述平缓震荡频谱,确定由于射线被所述准直器遮挡所产生的震荡对应的震荡周期;
从所述平缓震荡数据取出时长为所述震荡周期的替换数据,并利用所述替换数据替换所述初始震荡数据中的所述突变震荡数据,以得到所述目标震荡数据。
4.根据权利要求3所述的方法,其特征在于,所述根据所述平缓震荡频谱,确定由于射线被所述准直器遮挡所产生的震荡对应的震荡周期,包括:
确定所述平缓震荡频谱中峰值的频点;
根据所述峰值的频点、所述采样频率和所述指定时长,确定所述震荡周期。
5.根据权利要求1所述的方法,其特征在于,所述根据所述初始射线数据和任一个所述时间窗对应的所述平滑射线数据,确定所述初始射线数据中包括的初始震荡数据,包括:
将所述初始射线数据与任一个所述时间窗对应的所述平滑射线数据的差,作为所述初始震荡数据。
6.根据权利要求1所述的方法,其特征在于,在所述根据每个所述时间窗对应的所述平滑射线数据和所述初始震荡数据,对所述初始射线数据进行去除震荡处理,以得到目标射线数据之后,所述方法还包括:
根据所述目标射线数据,生成所述待测物体的CT图像。
7.一种射线数据的处理装置,其特征在于,应用于CT设备,所述CT设备包括射线发射器、准直器和探测器,所述射线发射器用于发出射线,射线通过所述准直器射向待测物体,所述探测器用于采集透过所述待测物体的射线;所述装置包括:
平滑处理模块,用于通过至少两个时间窗,对所述探测器采集到的初始射线数据进行平滑处理,以获取每个所述时间窗对应的平滑射线数据,每个所述时间窗的窗口长度不相同,所述初始射线数据包括指定时长内按照采样频率采集的射线经过所述待测物体的衰减值;
确定模块,用于根据所述初始射线数据和任一个所述时间窗对应的所述平滑射线数据,确定所述初始射线数据中包括的初始震荡数据,所述初始震荡数据能够反映由于射线被所述准直器遮挡所产生的震荡;
去除震荡模块,用于根据每个所述时间窗对应的所述平滑射线数据和所述初始震荡数据,对所述初始射线数据进行去除震荡处理,以得到目标射线数据;
所述去除震荡模块包括:分界确定子模块,用于根据每个所述时间窗对应的所述平滑射线数据,确定分界区间,以按照所述分界区间将所述初始震荡数据分为平缓震荡数据和突变震荡数据,所述平缓震荡数据的峰峰值小于所述突变震荡数据的峰峰值;
替换子模块,用于根据所述平缓震荡数据,替换所述初始震荡数据中的所述突变震荡数据,以得到目标震荡数据;
去除子模块,用于将所述初始射线数据与所述目标震荡数据的差,作为所述目标射线数据。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现权利要求1-6中任一项所述方法的步骤。
9.一种电子设备,其特征在于,包括:
存储器,其上存储有计算机程序;
处理器,用于执行所述存储器中的所述计算机程序,以实现权利要求1-6中任一项所述方法的步骤。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010568969.4A CN111839562B (zh) | 2020-06-19 | 2020-06-19 | 射线数据的处理方法、装置、存储介质和电子设备 |
US17/236,977 US11471120B2 (en) | 2020-06-19 | 2021-04-21 | Method and device for processing ray data, storage medium and electronic device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010568969.4A CN111839562B (zh) | 2020-06-19 | 2020-06-19 | 射线数据的处理方法、装置、存储介质和电子设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111839562A CN111839562A (zh) | 2020-10-30 |
CN111839562B true CN111839562B (zh) | 2023-04-11 |
Family
ID=72987750
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010568969.4A Active CN111839562B (zh) | 2020-06-19 | 2020-06-19 | 射线数据的处理方法、装置、存储介质和电子设备 |
Country Status (2)
Country | Link |
---|---|
US (1) | US11471120B2 (zh) |
CN (1) | CN111839562B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4891829A (en) * | 1986-11-19 | 1990-01-02 | Exxon Research And Engineering Company | Method and apparatus for utilizing an electro-optic detector in a microtomography system |
EP1000582A2 (en) * | 1998-11-12 | 2000-05-17 | Picker International, Inc. | Apparatus and method for imaging |
WO2011030460A1 (ja) * | 2009-09-14 | 2011-03-17 | 三菱重工業株式会社 | 放射線断層撮影方法 |
CN105103005A (zh) * | 2013-04-08 | 2015-11-25 | 艾尔默斯半导体股份公司 | 尤其是作为泊车辅助的、借助超声波测量的方法,以及超声波测量系统 |
WO2016158138A1 (ja) * | 2015-04-01 | 2016-10-06 | 株式会社日立製作所 | X線ct装置、再構成演算装置およびx線ct画像生成方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11051772B2 (en) * | 2016-04-08 | 2021-07-06 | Rensselaer Polytechnic Institute | Filtration methods for dual-energy X-ray CT |
CN107981877B (zh) * | 2016-10-27 | 2021-01-05 | 北京东软医疗设备有限公司 | 扫描设备的采样方法和采样装置 |
-
2020
- 2020-06-19 CN CN202010568969.4A patent/CN111839562B/zh active Active
-
2021
- 2021-04-21 US US17/236,977 patent/US11471120B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4891829A (en) * | 1986-11-19 | 1990-01-02 | Exxon Research And Engineering Company | Method and apparatus for utilizing an electro-optic detector in a microtomography system |
EP1000582A2 (en) * | 1998-11-12 | 2000-05-17 | Picker International, Inc. | Apparatus and method for imaging |
WO2011030460A1 (ja) * | 2009-09-14 | 2011-03-17 | 三菱重工業株式会社 | 放射線断層撮影方法 |
CN105103005A (zh) * | 2013-04-08 | 2015-11-25 | 艾尔默斯半导体股份公司 | 尤其是作为泊车辅助的、借助超声波测量的方法,以及超声波测量系统 |
WO2016158138A1 (ja) * | 2015-04-01 | 2016-10-06 | 株式会社日立製作所 | X線ct装置、再構成演算装置およびx線ct画像生成方法 |
Also Published As
Publication number | Publication date |
---|---|
US11471120B2 (en) | 2022-10-18 |
CN111839562A (zh) | 2020-10-30 |
US20210393228A1 (en) | 2021-12-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11154257B2 (en) | Imaging control device, imaging control method, and imaging control program | |
JP2004033778A (ja) | X線データを取得するためのシステム及び方法 | |
US10342503B2 (en) | Medical image processing apparatus, X-ray diagnostic apparatus, and X-ray computed tomography apparatus | |
JP2017042438A (ja) | 画像処理方法、画像処理装置及び放射線断層撮影装置並びにプログラム | |
US20120099709A1 (en) | Device and method to generate x-ray radiation with two spectra | |
JP2009536434A (ja) | サドル形状陽極を備える回転陽極x線管 | |
CN111839562B (zh) | 射线数据的处理方法、装置、存储介质和电子设备 | |
US10080540B2 (en) | Systems and methods for emission tomography/computed tomography data acquisition | |
JP2010506646A (ja) | 軸心臓ctに関する高周波サドル軌跡 | |
JP7250532B2 (ja) | X線ct装置及び撮影計画装置 | |
JP7242410B2 (ja) | 医用画像処理装置、x線ct装置及び学習用データの生成方法 | |
EP3970179A1 (en) | Maintaining a given focal spot size during a kvp switched spectral (multi-energy) imaging scan | |
US10984564B2 (en) | Image noise estimation using alternating negation | |
US20050031088A1 (en) | Radiographic apparatus and radiation detection signal processing method | |
JP2020049059A (ja) | 医用画像処理装置および方法 | |
JP2018175866A (ja) | X線ct装置 | |
JP7170469B2 (ja) | X線ct装置 | |
JP7356816B2 (ja) | 医用情報処理装置 | |
JP2017093646A (ja) | X線画像診断装置及び画像処理方法 | |
JP6873831B2 (ja) | 医用画像診断装置、医用画像処理装置及び医用画像処理プログラム | |
JP2014014418A (ja) | X線ct装置 | |
US20240078723A1 (en) | Medical image processing apparatus and medical image processing method | |
JP2009240577A (ja) | X線ctスキャナ及びx線管装置 | |
JP7206163B2 (ja) | X線ct装置、医用情報処理装置、及び医用情報処理プログラム | |
JP2020089497A (ja) | X線撮影装置 |
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 |