CN111184523B - 基于dr设备的三维图像重建方法及*** - Google Patents
基于dr设备的三维图像重建方法及*** Download PDFInfo
- Publication number
- CN111184523B CN111184523B CN202010052376.2A CN202010052376A CN111184523B CN 111184523 B CN111184523 B CN 111184523B CN 202010052376 A CN202010052376 A CN 202010052376A CN 111184523 B CN111184523 B CN 111184523B
- Authority
- CN
- China
- Prior art keywords
- scanning data
- correction
- ray
- pixel
- 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
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000012545 processing Methods 0.000 claims abstract description 28
- 238000001914 filtration Methods 0.000 claims abstract description 27
- 230000006870 function Effects 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 10
- 238000009499 grossing Methods 0.000 claims description 5
- 238000001514 detection method Methods 0.000 abstract description 7
- 210000000056 organ Anatomy 0.000 abstract description 5
- 238000002591 computed tomography Methods 0.000 abstract description 3
- 230000002159 abnormal effect Effects 0.000 description 8
- 230000003068 static effect Effects 0.000 description 8
- 238000003745 diagnosis Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000002601 radiography Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000002496 gastric effect Effects 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 210000001835 viscera Anatomy 0.000 description 2
- 238000001739 density measurement Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000035515 penetration Effects 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/40—Arrangements for generating radiation specially adapted for radiation diagnosis
-
- 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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
-
- 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/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- 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
-
- 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/10004—Still image; Photographic image
- G06T2207/10012—Stereo images
-
- 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/10116—X-ray image
- G06T2207/10124—Digitally reconstructed radiograph [DRR]
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Public Health (AREA)
- Radiology & Medical Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pulmonology (AREA)
- Computer Graphics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种基于DR设备的三维图像重建方法及***,方法包括:获取预设电压和预设电流条件下,不同角度的待检测体在DR设备上的X射线扫描数据;对获取到的X射线扫描数据依次进行暗校正、坏点坏线校正和空气校正,得到校正后的扫描数据;对校正后的扫描数据依次进行锥角加权、滤波处理和反投影,得到三维重建图像。本发明实现了将DR设备图像从二维扩展到三维,得到的三维图像无组织和脏器混叠,可以更加清晰、准确地反应待检测体的情况。通过一次DR扫描既可以获得二维的DR图像,还可以获得三维的CT图像,无需在进行DR检查后再进行CT扫描,大幅降低了患者接受X射线照射剂量,且可降低检查成本。
Description
技术领域
本发明涉及图像处理技术领域,尤其涉及一种基于DR设备的三维图像重建方法及***。
背景技术
DR(Digital Radiography)指在计算机控制下直接进行数字化X线摄影的一种技术,即X射线探测器把穿透人体的X射线信息转化为数字信号,并由计算机对图像进行后处理及显示。DR***主要包括X射线发生装置、X射线探测器、***控制器、影像监示器和影像处理工作站等组成。DR由于采用数字技术,因此可以根据临床需要进行各种图像后处理,实现如图像自动处理、边缘增强、放大漫游、图像拼接、兴趣区窗宽窗位调节以及距离、面积、密度测量等丰富的功能。DR技术动态范围广,X射线光量子检出效能(DQE)高,具有很宽的曝光宽容度,即使曝光条件稍差,也能获得很好的图像。此外,DR具有低剂量、高空间分辨率、短扫描时间及低成本等特点被广泛应用于体检及医学影像诊断领域,是医学影像诊断的主要设备之一。
一般来说DR设备分为两大类:静态DR和动态DR。静态DR,顾名思义,只能拍摄静态图像。医生指导病人完成摆位和设备设定后采用单次曝光的形式获取患者静态二维图像。其拍摄方式类似于日常生活中的拍摄照片,因此静态DR图像仅能反馈出患者曝光时刻的二维图像信息。如果摆位不佳或者患者在曝光过程中出现运动可能导致拍摄失败。动态DR将传统的“拍照”改为“摄像”,将连续曝光概念引入DR拍摄过程中,不但能够拍摄患者DR图像,同时能够通过观察患者情况,实现在诸多复杂***的拍摄中,一次性成功拍摄出清晰准确的图像,可以有效避免静态DR中患者需要多次反复拍摄的情况。动态DR由于加入了***,医生可以在多个不同角度下观察病变部位,进一步避免了漏诊,动态DR可以代替传统的胃肠机进行胃肠造影。
动态DR将传统的“拍照”式静态采集转变为“摄影”式透视采集,为医生提供了更多的信息用以辅助诊断。然而无论是静态DR还是动态DR其采集获得的图像都是二维图像,即在X射线穿透路径上所有器官存在着混叠现象,虽然DR设备具有较高的空间分辨率,但是无法精确的显示各个脏器及组织的实际位置及相关信息,无法提供精确的三维空间图像,限制了DR设备在临床的领用。
发明内容
本发明所要解决的技术问题是:提供一种基于DR设备的三维图像重建方法及***,可以得到无组织和脏器混叠的三维图像。
为了解决上述技术问题,本发明采用的技术方案为:
一种基于DR设备的三维图像重建方法,包括:
获取预设电压和预设电流条件下,不同角度的待检测体在DR设备上的X射线扫描数据;
对获取到的X射线扫描数据依次进行暗校正、坏点坏线校正和空气校正,得到校正后的扫描数据;
对校正后的扫描数据依次进行锥角加权、滤波处理和反投影,得到三维重建图像。
本发明采用的另一技术方案为:
一种基于DR设备的三维图像重建***,包括高压发生器、X射线发射器、带有角度测量及回传功能的旋转体、X射线探测器、机架和三维图像重建终端,所述高压发生器与所述X射线发射器电连接,所述带有角度测量及回传功能的旋转体用于放置待检测体,且所述带有角度测量及回传功能的旋转体设置于所述X射线发射器与X射线探测器之间,所述X射线发射器和X射线探测器分别固定设置于所述机架上,所述三维图像重建终端与所述X射线探测器电连接;
所述三维图像重建终端包括存储器、处理器和存储在所述存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现以下步骤:
获取预设电压和预设电流条件下,不同角度的待检测体在DR设备上的X射线扫描数据;
对获取到的X射线扫描数据依次进行暗校正、坏点坏线校正和空气校正,得到校正后的扫描数据;
对校正后的扫描数据依次进行锥角加权、滤波处理和反投影,得到三维重建图像。
本发明的有益效果在于:采集不同角度下的待检测体的X射线扫描数据,然后对扫描数据进行校正和三维重建,实现了将DR设备图像从二维扩展到三维,得到的三维图像无组织和脏器混叠,可以更加清晰、准确地反应待检测体的情况。通过一次DR扫描既可以获得二维的DR图像,还可以获得三维的CT图像,无需在进行DR检查后再进行CT扫描,大幅降低了患者接受X射线照射剂量,且可降低检查成本。
附图说明
图1为本发明实施例一的基于DR设备的三维图像重建方法的流程图;
图2为本发明实施例二的三维图像重建***的结构示意图;
图3为本发明实施例二的三维图像重建***的另一结构示意图;
图4为本发明实施例二的三维图像重建***的另一结构示意图;
图5为本发明实施例二的三维图像重建终端的示意图。
标号说明:
100、三维图像重建终端;1、存储器;2、处理器。
200、高压发生器;300、X射线发射器;400、带有角度测量及回传功能的旋转体;500、X射线探测器;600、机架;700、待检测体。
具体实施方式
为详细说明本发明的技术内容、所实现目的及效果,以下结合实施方式并配合附图予以说明。
本发明最关键的构思在于:采用DR设备采集不同角度下的待检测体的X射线扫描数据,然后对扫描数据进行校正和三维重建,实现了将DR设备图像从二维扩展到三维。
请参照图1,一种基于DR设备的三维图像重建方法,包括:
获取预设电压和预设电流条件下,不同角度的待检测体在DR设备上的X射线扫描数据;
对获取到的X射线扫描数据依次进行暗校正、坏点坏线校正和空气校正,得到校正后的扫描数据;
对校正后的扫描数据依次进行锥角加权、滤波处理和反投影,得到三维重建图像。
从上述描述可知,本发明的有益效果在于:采集不同角度下的待检测体的X射线扫描数据,然后对扫描数据进行校正和三维重建,实现了将DR设备图像从二维扩展到三维,得到的三维图像无组织和脏器混叠,可以更加清晰、准确地反应待检测体的情况。通过一次DR扫描既可以获得二维的DR图像,还可以获得三维图像,无需在进行DR检查后再进行CT扫描,大幅降低了患者接受X射线照射剂量,且可降低检查成本。
进一步的,对获取到的X射线扫描数据进行暗校正具体为:
在没有发射X射线的条件下,采集X射线探测器的数据,得到暗校正图数据;
将获取到的X射线扫描数据减去暗校正图数据,得到暗校正后的扫描数据。
由上述描述可知,进行暗校正的目的是为了去除DR设备X射线探测器自身采集数据不均性的问题。
进一步的,进行坏点坏线校正具体为:
预设像素范围;
判断暗校正后的扫描数据的一像素点是否在所述预设像素范围内;
若否,则将邻近的像素值在所述预设像素范围内的像素点的像素值作为所述一像素点的像素值。
由上述描述可知,进行坏点坏线校正的目的是为了去除X射线探测器存在的坏点和坏线。一般的DR设备X射线探测器均会存在不同数量的坏点及坏线,这些问题像素点在图像中表面为像素为0或者其他异常值,该异常值为一定值,不随扫描条件改变,如不进行校正会引起滤波出现异常,无图像或图像中出现严重的环状伪影问题。
进一步的,进行空气校正具体为:
采集不同电压、电流条件下的空气扫描数据,得到空气校正表;
根据所述空气校正表计算得到所述预设电压和预设电流条件下的入射光线强度;
分别对所述入射光线强度、坏点坏线校正后的扫描数据和空气扫描数据作ln对数运算,得到对数入射光线强度、对数扫描数据和对数空气扫描数据;
将所述对数入射光线强度减去对数扫描数据与对数空气扫描数据的差值,得到校正后的扫描数据。
由上述描述可知,进行空气校正的目的是校正X射线探测器对于X射线接收的不均性及残影的问题,通过计算X射线的入射光线强度且对扫描数据进行ln操作以获得最终可以用来三维重建的扫描数据。
进一步的,对校正后的扫描数据进行锥角加权具体为:
对校正后的扫描数据进行锥角余弦加权处理。
由上述描述可知,DR设备的X射线探测器具有较大的X射线锥角,这个锥角会导致扫描数据与实际不匹配而引起伪影。因此,需要降低锥角对扫描数据的影响,锥角越大权重越低,锥角越小权重越高。
进一步的,进行滤波处理具体为:
采用滤波器对锥角加权处理后的扫描数据进行滤波处理。
由上述描述可知,由于反投影操作需要将不同投影角度下的数据回模到三维重建图像中,这个过程相对于增强了重建图像中的低频分量,会导致图像变得模糊及不平坦。为了抑制这一现象,需要滤波器来抑制图像的低频区域同时增强高频信息。采用滤波器对扫描数据进行滤波以增强扫描数据的边界信息,获得最终的重建图像。
进一步的,所述对校正后的扫描数据依次进行锥角加权、滤波处理和反投影,得到三维重建图像之后还包括:
去除三维重建图像中的噪声点和对三维重建图像的边界进行平滑处理。
由上述描述可知,去除噪声点和平滑边界可以得到质量更好的三维重建图像,去除噪声点和平滑边界可以采用现有的已知手段进行。
请参照图2,本发明涉及的另一技术方案为:
一种基于DR设备的三维图像重建***,包括高压发生器200、X射线发射器300、带有角度测量及回传功能的旋转体400、X射线探测器500、机架600和三维图像重建终端100,所述高压发生器200与所述X射线发射器300电连接,所述带有角度测量及回传功能的旋转体400用于放置待检测体700,且所述带有角度测量及回传功能的旋转体400设置于所述X射线发射器300与X射线探测器500之间,所述X射线发射器300和X射线探测器500分别固定设置于所述机架600上,所述三维图像重建终端100与所述X射线探测器500电连接;
所述三维图像重建终端100包括存储器1、处理器2和存储在所述存储器1上并可在处理器2上运行的计算机程序,所述处理器2执行所述计算机程序时实现以下步骤:
获取预设电压和预设电流条件下,不同角度的待检测体在DR设备上的X射线扫描数据;
对获取到的X射线扫描数据依次进行暗校正、坏点坏线校正和空气校正,得到校正后的扫描数据;
对校正后的扫描数据依次进行锥角加权、滤波处理和反投影,得到三维重建图像。
进一步的,对获取到的X射线扫描数据进行暗校正具体为:
在没有发射X射线的条件下,采集X射线探测器的数据,得到暗校正图数据;
将获取到的X射线扫描数据减去暗校正图数据,得到暗校正后的扫描数据。
进一步的,进行坏点坏线校正具体为:
预设像素范围;
判断暗校正后的扫描数据的一像素点是否在所述预设像素范围内;
若否,则将邻近的像素值在所述预设像素范围内的像素点的像素值作为所述一像素点的像素值。
进一步的,进行空气校正具体为:
采集不同电压、电流条件下的空气扫描数据,得到空气校正表;
根据所述空气校正表计算得到所述预设电压和预设电流条件下的入射光线强度;
分别对所述入射光线强度、坏点坏线校正后的扫描数据和空气扫描数据作ln对数运算,得到对数入射光线强度、对数扫描数据和对数空气扫描数据;
将所述对数入射光线强度减去对数扫描数据与对数空气扫描数据的差值,得到校正后的扫描数据。
进一步的,对校正后的扫描数据进行锥角加权具体为:
对校正后的扫描数据进行锥角余弦加权处理。
进一步的,进行滤波处理具体为:
采用滤波器对锥角加权处理后的扫描数据进行滤波处理。
进一步的,所述处理器2执行所述计算机程序时还实现以下步骤:
所述对校正后的扫描数据依次进行锥角加权、滤波处理和反投影,得到三维重建图像之后还包括:
去除三维重建图像中的噪声点和对三维重建图像的边界进行平滑处理。
实施例一
请参照图1,本发明的实施例一为:
一种基于DR设备的三维图像重建方法,包括如下步骤:
S1、获取预设电压和预设电流条件下,不同角度的待检测体在DR设备上的X射线扫描数据。
本实施例中,通过DR设备采集X射线扫描数据。待检测体可以是人或动物等,待检测体可以站立或固定在一可旋转的带有角度测量及回传功能的旋转体上,通过带有角度测量及回传功能的旋转体带动待检测体转动即可获得不同角度的X射线扫描数据,旋转的速度可以根据需要进行调整,在旋转过程中待检测体需要保持不动。回传功能指的是旋转体可以将对应的角度信息记录并进行传递。预设电压和预设电流也可以根据实际情况进行调整。X射线发射器可以通过机架进行固定,机架可以是固定机架,也可以是具有滑轨的可调整机架。
S2、对获取到的X射线扫描数据依次进行暗校正、坏点坏线校正和空气校正,得到校正后的扫描数据。
本实施例中,对获取到的X射线扫描数据进行暗校正具体为:在没有发射X射线的条件下,采集X射线探测器的数据,得到暗校正图数据;将获取到的X射线扫描数据减去暗校正图数据,得到暗校正后的扫描数据。
进行暗校正的目的是为了去除DR设备X射线探测器自身采集数据不均性的问题。进行暗校正即为将得到的扫描数据中每一像素点的像素值分别减去对应的暗校正图中对应的像素点的像素值。
本实施例中,进行坏点坏线校正具体为:预设像素范围;判断暗校正后的扫描数据的一像素点是否在所述预设像素范围内;若否,则将邻近的像素值在所述预设像素范围内的像素点的像素值作为所述一像素点的像素值;若是,则所述一像素点的像素值不变。
进行坏点坏线校正的目的是为了去除X射线探测器存在的坏点和坏线。一般的DR设备X射线探测器均会存在不同数量的坏点及坏线,这些问题像素点在图像中表面为像素为0或者其他异常值,该异常值为一定值,不随扫描条件改变,如不进行校正会引起滤波出现异常,无图像或图像中出现严重的环状伪影问题。
本实施例中,进行空气校正具体为:采集不同电压、电流条件下的空气扫描数据,得到空气校正表;根据所述空气校正表计算得到所述预设电压和预设电流条件下的入射光线强度;分别对所述入射光线强度、坏点坏线校正后的扫描数据和空气扫描数据作ln对数运算,得到对数入射光线强度、对数扫描数据和对数空气扫描数据;将所述对数入射光线强度减去对数扫描数据与对数空气扫描数据的差值,得到校正后的扫描数据。本实施例中,可以通过计算空气校正表中所有像素点的均值作为入射光线强度,也可以选择任意一个或多个像素的均值作为入射光线强度。
进行空气校正的目的是校正X射线探测器对于X射线接收的不均性及残影的问题,通过计算X射线的入射光线强度且对扫描数据进行ln操作以获得最终可以用来三维重建的扫描数据。所述空气扫描数据即X射线仅穿透空气被探测器接收获得的数据。
S3、对校正后的扫描数据依次进行锥角加权、滤波处理和反投影,得到三维重建图像。
本实施例中,对校正后的扫描数据进行锥角加权具体为:对校正后的扫描数据进行锥角余弦加权处理。DR设备的X射线探测器具有较大的X射线锥角,这个锥角会导致扫描数据与实际不匹配而引起伪影。因此,需要降低锥角对扫描数据的影响,锥角越大权重越低,锥角越小权重越高。此外,还可以通过设计线性、高斯等权重来对扫描数据进行加权。
本实施例中,进行滤波处理具体为:采用滤波器对锥角加权处理后的扫描数据进行滤波处理。由于反投影操作需要将不同投影角度下的数据回模到三维重建图像中,这个过程相对于增强了重建图像中的低频分量,会导致图像变得模糊及不平坦。为了抑制这一现象,需要滤波器来抑制图像的低频区域同时增强高频信息。采用滤波器对扫描数据进行滤波以增强扫描数据的边界信息,获得最终的重建图像。
反投影过程是针对任意一个三维重建像素点计算在任意投影角度下该点对应的扫描数据,并将这个数据累加到重建点位置以获得最终重建数据的过程。
本实施例中,步骤S3之后还包括:
去除三维重建图像中的噪声点和对三维重建图像的边界进行平滑处理。最后,将处理后的三维重建图像进行显示。
实施例二
请参照图2至图5,本发明的实施例二为:
一种基于DR设备的三维图像重建***,如图2至图4所示,包括高压发生器200、X射线发射器300、带有角度测量及回传功能的旋转体400、X射线探测器500、机架600和三维图像重建终端100,所述高压发生器200与所述X射线发射器300电连接,所述带有角度测量及回传功能的旋转体400用于放置待检测体700,且所述带有角度测量及回传功能的旋转体400设置于所述X射线发射器300与X射线探测器500之间,所述X射线发射器300和X射线探测器500分别固定设置于所述机架600上,所述三维图像重建终端100与所述X射线探测器500电连接。本实施例中,回传功能指的是旋转体可以将对应的角度信息记录并传递给三维图像重建终端100。
图2中,机架600为一弧形结构,X射线发射器300和X射线探测器500分别固定设置于所述机架600的两端,带有角度测量及回传功能的旋转体400为一带有角度控制器的装置,可以根据需要调整待检测体700的角度。
图3中,机架600为滑轨式,机架600包括两条滑轨,X射线发射器300固定在其中一条滑轨上,X射线探测器500固定在另一条滑轨上,可以根据需要调整X射线发射器300和X射线探测器500在机架600上的位置。
图4中,机架600为立柱式,包括两个立柱,X射线发射器300固定在其中一个立柱上,X射线探测器500固定在另一个立柱上,可以根据需要调整两个立柱的位置。
如图5所示,所述三维图像重建终端100包括存储器1、处理器2和存储在所述存储器1上并可在处理器2上运行的计算机程序,所述处理器2执行所述计算机程序时实现以下步骤:
获取预设电压和预设电流条件下,不同角度的待检测体在DR设备上的X射线扫描数据;
对获取到的X射线扫描数据依次进行暗校正、坏点坏线校正和空气校正,得到校正后的扫描数据;
对校正后的扫描数据依次进行锥角加权、滤波处理和反投影,得到三维重建图像。
进一步的,对获取到的X射线扫描数据进行暗校正具体为:
在没有发射X射线的条件下,采集X射线探测器的数据,得到暗校正图数据;
将获取到的X射线扫描数据减去暗校正图数据,得到暗校正后的扫描数据。
进一步的,进行坏点坏线校正具体为:
预设像素范围;
判断暗校正后的扫描数据的一像素点是否在所述预设像素范围内;
若否,则将邻近的像素值在所述预设像素范围内的像素点的像素值作为所述一像素点的像素值。
进一步的,进行空气校正具体为:
采集不同电压、电流条件下的空气扫描数据,得到空气校正表;
根据所述空气校正表计算得到所述预设电压和预设电流条件下的入射光线强度;
分别对所述入射光线强度、坏点坏线校正后的扫描数据和空气扫描数据作ln对数运算,得到对数入射光线强度、对数扫描数据和对数空气扫描数据;
将所述对数入射光线强度减去对数扫描数据与对数空气扫描数据的差值,得到校正后的扫描数据。
进一步的,对校正后的扫描数据进行锥角加权具体为:
对校正后的扫描数据进行锥角余弦加权处理。
进一步的,进行滤波处理具体为:
采用滤波器对锥角加权处理后的扫描数据进行滤波处理。
进一步的,所述处理器2执行所述计算机程序时还实现以下步骤:
所述对校正后的扫描数据依次进行锥角加权、滤波处理和反投影,得到三维重建图像之后还包括:
去除三维重建图像中的噪声点和对三维重建图像的边界进行平滑处理。
综上所述,本发明提供的一种基于DR设备的三维图像重建方法及***,实现了将DR设备图像从二维扩展到三维,且得到的三维图像无组织和脏器混叠,可以更加清晰、准确地反应待检测体的情况。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等同变换,或直接或间接运用在相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (5)
1.一种基于DR设备的三维图像重建方法,其特征在于,包括:
获取预设电压和预设电流条件下,不同角度的待检测体在DR设备上的X射线扫描数据;
对获取到的X射线扫描数据依次进行暗校正、坏点坏线校正和空气校正,得到校正后的扫描数据;
对校正后的扫描数据依次进行锥角加权、滤波处理和反投影,得到三维重建图像;
对获取到的X射线扫描数据进行暗校正具体为:
在没有发射X射线的条件下,采集X射线探测器的数据,得到暗校正图数据;
将获取到的X射线扫描数据减去暗校正图数据,得到暗校正后的扫描数据;
进行坏点坏线校正具体为:
预设像素范围;
判断暗校正后的扫描数据的一像素点是否在所述预设像素范围内;
若否,则将邻近的像素值在所述预设像素范围内的像素点的像素值作为所述一像素点的像素值;
进行空气校正具体为:
采集不同电压、电流条件下的空气扫描数据,得到空气校正表;
根据所述空气校正表计算得到所述预设电压和预设电流条件下的入射光线强度;
分别对所述入射光线强度、坏点坏线校正后的扫描数据和空气扫描数据作ln对数运算,得到对数入射光线强度、对数扫描数据和对数空气扫描数据;
将所述对数入射光线强度减去对数扫描数据与对数空气扫描数据的差值,得到校正后的扫描数据。
2.根据权利要求1所述的基于DR设备的三维图像重建方法,其特征在于,对校正后的扫描数据进行锥角加权具体为:
对校正后的扫描数据进行锥角余弦加权处理。
3.根据权利要求2所述的基于DR设备的三维图像重建方法,其特征在于,进行滤波处理具体为:
采用滤波器对锥角加权处理后的扫描数据进行滤波处理。
4.根据权利要求1所述的基于DR设备的三维图像重建方法,其特征在于,所述对校正后的扫描数据依次进行锥角加权、滤波处理和反投影,得到三维重建图像之后还包括:
去除三维重建图像中的噪声点和对三维重建图像的边界进行平滑处理。
5.一种基于DR设备的三维图像重建***,其特征在于,包括高压发生器、X射线发射器、带有角度测量及回传功能的旋转体、X射线探测器、机架和三维图像重建终端,所述高压发生器与所述X射线发射器电连接,所述带有角度测量及回传功能的旋转体用于放置待检测体,且所述带有角度测量及回传功能的旋转体设置于所述X射线发射器与X射线探测器之间,所述X射线发射器和X射线探测器分别固定设置于所述机架上,所述三维图像重建终端与所述X射线探测器电连接;
所述三维图像重建终端包括存储器、处理器和存储在所述存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现以下步骤:
获取预设电压和预设电流条件下,不同角度的待检测体在DR设备上的X射线扫描数据;
对获取到的X射线扫描数据依次进行暗校正、坏点坏线校正和空气校正,得到校正后的扫描数据;
对校正后的扫描数据依次进行锥角加权、滤波处理和反投影,得到三维重建图像;
对获取到的X射线扫描数据进行暗校正具体为:
在没有发射X射线的条件下,采集X射线探测器的数据,得到暗校正图数据;
将获取到的X射线扫描数据减去暗校正图数据,得到暗校正后的扫描数据;
进行坏点坏线校正具体为:
预设像素范围;
判断暗校正后的扫描数据的一像素点是否在所述预设像素范围内;
若否,则将邻近的像素值在所述预设像素范围内的像素点的像素值作为所述一像素点的像素值;
进行空气校正具体为:
采集不同电压、电流条件下的空气扫描数据,得到空气校正表;
根据所述空气校正表计算得到所述预设电压和预设电流条件下的入射光线强度;
分别对所述入射光线强度、坏点坏线校正后的扫描数据和空气扫描数据作ln对数运算,得到对数入射光线强度、对数扫描数据和对数空气扫描数据;
将所述对数入射光线强度减去对数扫描数据与对数空气扫描数据的差值,得到校正后的扫描数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010052376.2A CN111184523B (zh) | 2020-01-17 | 2020-01-17 | 基于dr设备的三维图像重建方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010052376.2A CN111184523B (zh) | 2020-01-17 | 2020-01-17 | 基于dr设备的三维图像重建方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111184523A CN111184523A (zh) | 2020-05-22 |
CN111184523B true CN111184523B (zh) | 2023-03-10 |
Family
ID=70685030
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010052376.2A Active CN111184523B (zh) | 2020-01-17 | 2020-01-17 | 基于dr设备的三维图像重建方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111184523B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111419255A (zh) * | 2020-03-24 | 2020-07-17 | 深圳锐探科技有限公司 | 一种ct成像方法及装置 |
CN113223109B (zh) * | 2021-05-27 | 2022-08-12 | 明峰医疗***股份有限公司 | 一种基于ct轴扫描的坏点校正方法 |
CN113570705B (zh) * | 2021-07-28 | 2024-04-30 | 广州瑞多思医疗科技有限公司 | 一种三维剂量重建方法、装置、计算机设备及存储介质 |
CN113570708A (zh) * | 2021-07-30 | 2021-10-29 | 重庆市特种设备检测研究院 | 缺陷立体建模方法、装置及计算机可读存储介质 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005077627A (ja) * | 2003-08-29 | 2005-03-24 | Fuji Photo Film Co Ltd | 放射線画像情報読取装置における暗補正方法および暗補正装置 |
JP2005296340A (ja) * | 2004-04-12 | 2005-10-27 | Canon Inc | コーンビームx線ct撮影装置とそれを用いた画像取得方法 |
CN1881075A (zh) * | 2005-06-15 | 2006-12-20 | 上海医疗器械厂有限公司 | 一种x射线数字成像的校正方法 |
CN101371786A (zh) * | 2007-08-24 | 2009-02-25 | 北京师范大学珠海分校 | 一种x射线图像三维重构的方法及*** |
CN101795349A (zh) * | 2010-03-24 | 2010-08-04 | 深圳迈瑞生物医疗电子股份有限公司 | 一种处理x线平板探测器图像坏线的方法及其装置 |
CN102762151A (zh) * | 2010-01-21 | 2012-10-31 | 卡尔斯特里姆保健公司 | 四维体积成像*** |
CN103366389A (zh) * | 2013-04-27 | 2013-10-23 | 中国人民解放军北京军区总医院 | Ct图像重建方法 |
CN103767723A (zh) * | 2012-10-25 | 2014-05-07 | 南京普爱射线影像设备有限公司 | 一种基于c形臂的锥束ct三维数字成像方法 |
WO2019000601A1 (zh) * | 2017-06-27 | 2019-01-03 | 西安立人医学科技有限公司 | 锥形束ct多方向扫描仪 |
CN109730713A (zh) * | 2019-01-28 | 2019-05-10 | 深圳市纳诺艾医疗科技有限公司 | 一种二次荧光广角多功能成像设备 |
WO2019149809A1 (de) * | 2018-02-02 | 2019-08-08 | Gom Gmbh | Verfahren zur kalibrierung einer messeinrichtung sowie computerprogramm hierzu |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4172201B2 (ja) * | 2002-04-30 | 2008-10-29 | コニカミノルタホールディングス株式会社 | 放射線撮影装置及び放射線画像形成装置 |
KR102104534B1 (ko) * | 2013-06-12 | 2020-04-27 | 삼성전자주식회사 | 엑스선 촬영 장치 및 그 제어 방법 |
-
2020
- 2020-01-17 CN CN202010052376.2A patent/CN111184523B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005077627A (ja) * | 2003-08-29 | 2005-03-24 | Fuji Photo Film Co Ltd | 放射線画像情報読取装置における暗補正方法および暗補正装置 |
JP2005296340A (ja) * | 2004-04-12 | 2005-10-27 | Canon Inc | コーンビームx線ct撮影装置とそれを用いた画像取得方法 |
CN1881075A (zh) * | 2005-06-15 | 2006-12-20 | 上海医疗器械厂有限公司 | 一种x射线数字成像的校正方法 |
CN101371786A (zh) * | 2007-08-24 | 2009-02-25 | 北京师范大学珠海分校 | 一种x射线图像三维重构的方法及*** |
CN102762151A (zh) * | 2010-01-21 | 2012-10-31 | 卡尔斯特里姆保健公司 | 四维体积成像*** |
CN101795349A (zh) * | 2010-03-24 | 2010-08-04 | 深圳迈瑞生物医疗电子股份有限公司 | 一种处理x线平板探测器图像坏线的方法及其装置 |
CN103767723A (zh) * | 2012-10-25 | 2014-05-07 | 南京普爱射线影像设备有限公司 | 一种基于c形臂的锥束ct三维数字成像方法 |
CN103366389A (zh) * | 2013-04-27 | 2013-10-23 | 中国人民解放军北京军区总医院 | Ct图像重建方法 |
WO2019000601A1 (zh) * | 2017-06-27 | 2019-01-03 | 西安立人医学科技有限公司 | 锥形束ct多方向扫描仪 |
WO2019149809A1 (de) * | 2018-02-02 | 2019-08-08 | Gom Gmbh | Verfahren zur kalibrierung einer messeinrichtung sowie computerprogramm hierzu |
CN109730713A (zh) * | 2019-01-28 | 2019-05-10 | 深圳市纳诺艾医疗科技有限公司 | 一种二次荧光广角多功能成像设备 |
Also Published As
Publication number | Publication date |
---|---|
CN111184523A (zh) | 2020-05-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111184523B (zh) | 基于dr设备的三维图像重建方法及*** | |
JP3257630B2 (ja) | 冠状動脈を形状特定する表示装置 | |
US9044190B2 (en) | C-arm computerized tomography system | |
US20090060121A1 (en) | Computed tomography data acquisition apparatus and method | |
CN105030266A (zh) | X射线计算机断层摄影装置以及扫描计划设定支援装置 | |
US8249216B2 (en) | X-ray moving image radiographing apparatus | |
JP4537129B2 (ja) | トモシンセシス用途における対象物を走査するためのシステム | |
JP6640527B2 (ja) | X線ct装置 | |
JP2007089674A (ja) | 外観形状計測装置およびx線ct装置 | |
US9194965B2 (en) | System and method for X-ray image acquisition and processing | |
CN103860185B (zh) | 多能量图像的确定 | |
JPWO2005011502A1 (ja) | 放射線断層撮影装置 | |
JP3897925B2 (ja) | コーンビームct装置 | |
EP3349655B1 (en) | Tomography apparatus and controlling method for the same | |
CN101721222B (zh) | 一种修正床板和摆位辅助装置对图像质量影响的方法 | |
KR101768520B1 (ko) | 흉부의 디지털 x선 일반촬영 및 디지털 단층영상합성의 영상을 통합적 및 연속적으로 획득하기 위한 디지털 x선 촬영 시스템의 제어방법 | |
US20220071578A1 (en) | Improved method of acquiring a radiographic scan of a region-of-interest in a metal containing object | |
US20230145920A1 (en) | Systems and methods for motion detection in medical images | |
US20220414832A1 (en) | X-ray imaging restoration using deep learning algorithms | |
CN110084753B (zh) | 基于多传感器融合的动态dr图像拼接方法及终端 | |
CN109243585B (zh) | 医学图像的生成方法和医学图像处理***及其交互方法 | |
CN110730977B (zh) | 低剂量成像方法及装置 | |
US11622741B2 (en) | System and method for imaging a subject | |
JP2014195492A (ja) | X線ct装置 | |
KR20190002960A (ko) | 신체인식을 이용한 작업 기반 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 |