CN112545543B - 基于病床运动信息的扫描运动监测方法、***及存储介质 - Google Patents
基于病床运动信息的扫描运动监测方法、***及存储介质 Download PDFInfo
- Publication number
- CN112545543B CN112545543B CN202110188324.2A CN202110188324A CN112545543B CN 112545543 B CN112545543 B CN 112545543B CN 202110188324 A CN202110188324 A CN 202110188324A CN 112545543 B CN112545543 B CN 112545543B
- Authority
- CN
- China
- Prior art keywords
- motion
- scanning
- motion vector
- image
- camera
- 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 52
- 238000012544 monitoring process Methods 0.000 title claims abstract description 27
- 239000013598 vector Substances 0.000 claims abstract description 52
- 230000008569 process Effects 0.000 claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims description 15
- 230000009466 transformation Effects 0.000 claims description 12
- 238000002591 computed tomography Methods 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 6
- 238000003384 imaging method Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 5
- 238000001514 detection method Methods 0.000 claims description 5
- 230000003287 optical effect Effects 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 2
- 238000002059 diagnostic imaging Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 208000012661 Dyskinesia Diseases 0.000 description 1
- 208000015592 Involuntary movements Diseases 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000017311 musculoskeletal movement, spinal reflex action Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000000241 respiratory effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000021542 voluntary musculoskeletal movement 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/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
- A61B6/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
-
- 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]
-
- 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/547—Control of apparatus or devices for radiation diagnosis involving tracking of position of the device or parts of the device
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及医学成像技术领域,公开了一种基于病床运动信息的扫描运动监测方法、***及存储介质,方法包括步骤:在扫描孔径附近设置用于采集病人待检测部位图像的摄像机,对摄像机进行几何校准;扫描过程中获取待检测部位的彩色图像和深度图像;选择相邻帧获取的彩色图像和深度图像,进行图像配准和特征点匹配;计算特征点在世界坐标系中的位置和运动向量以及待检测部位随病床运动的运动向量,两者差值与预设的阈值相比较,作为判断待检测部位是否存在运动的依据。本发明能够根据监测结果对扫描过程进行及时调整,解决因自主或非自主运动破环投影数据的一致性和完整性的问题,提高图像质量,同时能够适用于轴扫或螺旋扫描多种扫描方式。
Description
技术领域
本发明涉及CT图像技术领域,具体涉及一种基于病床运动信息的扫描运动监测方法、***及存储介质。
背景技术
CT医学成像***自20世纪70年代发明了之后经过了长足的进步,扫描速度从开始要几分钟到现在的0.2秒。探测器排数也从开始的单排双排,到现在的64排,128排,甚至256排。这其中的变化不单单是***硬件的升级换代,***的图像重建技术也带来了革命性的变化。初期的CT***由于只有一排探测器,所以X射线束是扇形束,所用到的重建技术也都是二维的扇形束重建技术。因为每次只能扫描一层,整个扫描需要很长的时间,后来螺旋多排CT引入就是为了加快扫描的速度,比如16排,32排,64排等***。
轴扫和螺旋是CT扫描常用的两种扫描模式。理想情况下,整个CT扫描过程中,患者相对病床位置保持静止,随着病床一起运动。但是实际扫描过程中总有可能存在运动,比如头部不经意的转动,呼吸运动和心跳等等。这些自主和非自主的运动破环了投影数据的一致性和完整性,如果不进行更正将导致重建图像严重伪影,影响医生诊断。因此监测病人运动对CT成像有重要的意义。
减小运动伪影最有效的方法是提高投影数据的采集频率,也就是提高机架的旋转速度,例如0.25s每圈。但是由于机械和电子方面硬件的限制,继续提高转速面临巨大的挑战。为了抑制运动伪影减少病人不必要的辐射剂量,一些方法被提出用来解决CT扫描过程中运动的问题。如申请号为201510473092.X的专利“一种基于投影图像的CT扫描实时运动监测方法”,用当前投影的互信息和经验投影的互信息差异来判断是否有运动,经验值的选取以及病人非自主运动都会影响判断结果。申请号为201611238467.5的专利“用于计算机断层扫描的运动检测方法和装置”中,通过计算相隔180度整数倍的两幅投影的质心偏差估计是否运动,该方法同样不适用于螺旋扫描的情况。
发明内容
技术目的:针对上述技术问题,本发明公开了一种基于病床运动信息的扫描运动监测方法、***及存储介质,同时能够适用于轴扫或螺旋扫描多种扫描方式、有利于消除运动伪影以提高CT图像质量。
技术方案:为实现上述技术目的,本发明采用了如下技术方案:
一种基于病床运动信息的扫描运动监测方法,其特征在于,包括步骤:
步骤一、相机校准:在CT扫描机架的预设位置上设置用于采集扫描孔径内的病人待检测部位的图像的摄像机,对摄像机进行几何校准,确定摄像机的内部参数和外部参数,内部参数为摄像机的内部几何和光学特征,外部参数为摄像机坐标系三维位置和方向,建立成像模型,确定摄像机坐标***与世界坐标系的变换矩阵;
步骤二、图像监测:扫描过程中,通过摄像机对扫描孔径内的待检测部位进行连续监测和图像采集,得到彩色图像和深度图像;
步骤三、图像配准及特征点匹配:选择相邻帧获取的彩色图像和深度图像,首先对每帧数据的彩色图像和深度图像配准对齐,然后对配准后的图像进行特征点匹配;
步骤四、计算第一运动向量:计算特征点在世界坐标系中的位置和运动向量,记为第一运动向量;
步骤五、根据运动方向和运动速度来估计是否存在运动:计算扫描孔径内的待检测部位随CT病床运动的运动向量,记为第二运动向量,第一运动向量和第二运动向量之差作为待检测部位的绝对运动向量,将绝对运动向量与预设的阈值相比较,判断待检测部位是否存在运动;若存在运动则终止扫描或对扫描数据进行校正重建,若不存在运动则继续扫描或对扫描数据进行常规重建。
步骤一中采用的校准方法为棋盘格方法,首先输入棋盘格标定图,然后找到棋盘格中角点,再计算相机内参、外参和畸变参数,根据成像模型和变换矩阵,对相机采集到的图像进行校准,最后再对图像进行畸变校正;经过畸变校正后图像坐标到世界坐标的转换公式如下:
具体地,所述步骤四中,描述特征点在世界坐标系中的位置和运动的第一运动向量的计算公式为:
步骤五中,描述扫描孔径内的待检测部位随CT病床运动的第二运动向量为:
任意维度的大于阈值T表示该方向存在运动。
具体地,所述步骤一中采用的校准方法为直线线性变换DLT方法、RAC方法或平面标定方法中的任一种。
具体地,所述步骤三中采用SURF方法完成定位局部特征点和特征点匹配。
具体地,理想情况下,阈值T为0;若存在噪声和测量误差,阈值T根据实际环境实验测定。
一种CT扫描***,其特征在于,包括:
CT投影数据获取模块,用于获取病床上病人待检测部位的投影数据;
相机图像获取模块,用于CT检测过程中,连续获取病床上待检测部位的彩色图像和深度图像;
计算模块,用于对彩色图像和深度图像进行图像配准和特征点匹配,并计算出用于描述特征点在世界坐标系中位置和运动变化的第一运动向量;且根据病床速度信息,计算用于描述待检测部位随病床运动变化的第二运动向量;然后计算出第一运动向量和第二运动向量的差值,作为绝对运动向量,并比较绝对运动向量与阈值的差异;
控制模块,发出对应的控制指令给其它功能模块;
运动检测模块,用于根据绝对运动向量与阈值的差异判断被检测对象是否运动;并根据判断结果向控制模块发出是否进行中止扫描或继续扫描的信息。
一种存储介质,其特征在于:所述存储介质存储有至少一个可被处理器执行的指令,其中至少一个指令被处理器执行时,用于执行所述基于病床运动信息的扫描运动监测方法。
有益效果:由于采用了上述技术方案,本发明具有如下技术效果:
本发明的基于病床运动信息的运动监测方法,利用放置在扫描孔径附近的摄像采集病人图像,并结合当前病床的运动速度,实现对视野内待检测部位的运动监测,进一步实现扫描过程的及时调整,解决现有技术中自主和非自主的运动破环投影数据的一致性和完整性的问题,同时能够适用于轴扫或螺旋扫描多种扫描方式,适用范围广。
附图说明
图1为本发明的基于病床运动信息的扫描运动监测方法的示意图;
图2为本发明中摄像头/机布置在扫描机架上的示意图;
图3为本发明采用棋盘格方法对相机进行校准和图像处理的流程;
图4为对每帧数据的彩色图像和深度图像进行关键点匹配的示意图;
图5为运动向量估计的示意图。
具体实施方式
下面结合附图对本发明进行详细的描述。
本发明是利用放置在扫描孔径附近的摄像头来对病人的运动进行监测,并结合当前病床运动的速度,对视野内的物体的运动进行监测,流程如图2所示。摄像头可以安装在机架孔径上方,也可以安装在满足需求的其他位置。
本发明提供的基于病床运动信息的扫描运动监测方法,包括如下步骤:
1、对光学相机的几何进行校准。
校准的目的是确定摄像机内部几何和光学特性(内部参数),以及和相对世界坐标系对应的摄像机坐标系的三维位置和方向(外部参数),并建立成像模型,以便确定空间坐标系中物体点同它所在像平面上像点的对应关系,即确定变换矩阵。相机校准方法包含直接线性变换DLT方法、RAC方法、平面标定方法等。
本案例选用棋盘格方法,处理流程如图3所示:首先输入棋盘格标定图,然后找到棋盘格中角点,再计算相机内参、外参和畸变参数,根据成像模型和变换矩阵,对相机采集到的图像进行校准,最后再对图像进行畸变校正。
经过畸变校正后图像坐标到世界坐标的转换公式(1):
2、扫描过程中对扫描孔径内的物体进行连续监测和采集,得到彩色信息和深度信息。采集与CT投影采集同步可以定位到发生运动的CT采样角度位置。深度图像也被称为距离影像,是指将从图像采集器到场景中各点的距离(深度)作为像素值的图像,它能够直接反映被采集对像可见表面的几何形状。
3、根据相邻帧获取的RBG彩色图像和深度图像D计算运动向量。首先对每帧数据的彩色图像和深度图像配准对齐,然后对配准后的图像进行特征点匹配,如SURF方法完成定位局部特征点和特征点匹配。
4、根据公式(1),计算特征点在世界坐标系中的位置以及其运动向量:
t为相邻两帧采集图像之间的时间间隔。如图5所示,则物体绝对运动向量为:
任意维度的大于阈值T表示该方向存在运动。理想情况下T=0,由于噪声和测量精度的影响,阈值T根据实际环境实验测定。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (8)
1.一种基于病床运动信息的扫描运动监测方法,其特征在于,包括步骤:
步骤一、相机校准:在CT扫描机架的预设位置上设置用于采集扫描孔径内的病人待检测部位的图像的摄像机,对摄像机进行几何校准,确定摄像机的内部参数和外部参数,内部参数为摄像机的内部几何和光学特征,外部参数为摄像机坐标系三维位置和方向,建立成像模型,确定摄像机坐标***与世界坐标系的变换矩阵;
步骤二、图像监测:扫描过程中,通过摄像机对扫描孔径内的待检测部位进行连续监测和图像采集,得到彩色图像和深度图像;
步骤三、图像配准及特征点匹配:选择相邻帧获取的彩色图像和深度图像,首先对每帧数据的彩色图像和深度图像配准对齐,然后对配准后的图像进行特征点匹配;
步骤四、计算第一运动向量:计算特征点在世界坐标系中的位置和运动向量,记为第一运动向量;
步骤五、根据运动方向和运动速度来估计是否存在运动:计算扫描孔径内的待检测部位随CT病床运动的运动向量,记为第二运动向量,第一运动向量和第二运动向量之差作为待检测部位的绝对运动向量,将绝对运动向量与预设的阈值相比较,判断待检测部位是否存在运动;若存在运动则终止扫描或对扫描数据进行校正重建,若不存在运动则继续扫描或对扫描数据进行常规重建。
4.根据权利要求1所述的基于病床运动信息的扫描运动监测方法,其特征在于:所述步骤一中采用的校准方法为直线线性变换DLT方法、RAC方法或平面标定方法中的任一种。
5.根据权利要求1所述的基于病床运动信息的扫描运动监测方法,其特征在于:所述步骤三中采用SURF方法完成定位局部特征点和特征点匹配。
6.根据权利要求4所述的基于病床运动信息的扫描运动监测方法,其特征在于:理想情况下,预设的阈值为0;若存在噪声和测量误差,预设的阈值根据实际环境实验测定。
7.一种CT扫描***,其特征在于,包括:
CT投影数据获取模块,用于获取病床上病人待检测部位的投影数据;
相机图像获取模块,用于CT检测过程中,连续获取病床上待检测部位的彩色图像和深度图像;
计算模块,用于对彩色图像和深度图像进行图像配准和特征点匹配,并计算出用于描述特征点在世界坐标系中位置和运动变化的第一运动向量;且根据病床速度信息,计算用于描述待检测部位随病床运动变化的第二运动向量;然后计算出第一运动向量和第二运动向量的差值,作为绝对运动向量,并比较绝对运动向量与阈值的差异;
控制模块,发出对应的控制指令给其它功能模块;
运动检测模块,用于根据绝对运动向量与阈值的差异判断被检测对象是否运动;并根据判断结果向控制模块发出是否进行中止扫描或继续扫描的信息。
8.一种存储介质,其特征在于:所述存储介质存储有至少一个可被处理器执行的指令,其中至少一个指令被处理器执行时,用于执行如权利要求1至6中任一所述基于病床运动信息的扫描运动监测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110188324.2A CN112545543B (zh) | 2021-02-19 | 2021-02-19 | 基于病床运动信息的扫描运动监测方法、***及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110188324.2A CN112545543B (zh) | 2021-02-19 | 2021-02-19 | 基于病床运动信息的扫描运动监测方法、***及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112545543A CN112545543A (zh) | 2021-03-26 |
CN112545543B true CN112545543B (zh) | 2021-05-11 |
Family
ID=75035925
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110188324.2A Active CN112545543B (zh) | 2021-02-19 | 2021-02-19 | 基于病床运动信息的扫描运动监测方法、***及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112545543B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114391861B (zh) * | 2021-12-16 | 2024-06-04 | 纳米维景(上海)医疗科技有限公司 | Z型扫描床位置精度的校正方法、存储介质及终端设备 |
CN114677424B (zh) * | 2022-05-26 | 2022-09-06 | 浙江天新智能研究院有限公司 | 一种无人值守螺旋卸船机点云数据处理方法 |
CN115381471B (zh) * | 2022-10-26 | 2023-03-17 | 南方医科大学南方医院 | 一种基于运动检测的ct扫描辅助***及方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6257962B2 (ja) * | 2012-09-04 | 2018-01-10 | 東芝メディカルシステムズ株式会社 | X線ct装置 |
CN103315739B (zh) * | 2013-05-22 | 2015-08-19 | 华东师范大学 | 基于动态跟踪技术免除运动伪影的磁共振影像方法和*** |
KR101716421B1 (ko) * | 2013-06-21 | 2017-03-14 | 삼성전자주식회사 | 정보 제공 방법 및 정보 제공을 위한 의료 진단 장치 |
US20160247293A1 (en) * | 2015-02-24 | 2016-08-25 | Brain Biosciences, Inc. | Medical imaging systems and methods for performing motion-corrected image reconstruction |
CN105006004B (zh) * | 2015-08-05 | 2018-04-03 | 天津金曦医疗设备有限公司 | 一种基于投影图像的ct扫描实时运动监测方法 |
CN106821407A (zh) * | 2016-12-28 | 2017-06-13 | 上海联影医疗科技有限公司 | 用于计算机断层扫描的运动检测方法和装置 |
CN108969002A (zh) * | 2018-07-26 | 2018-12-11 | 沈阳东软医疗***有限公司 | 医疗扫描设备、医疗扫描***及扫描监测方法 |
US11510629B2 (en) * | 2018-12-26 | 2022-11-29 | General Electric Company | Systems and methods for detecting patient state in a medical imaging session |
CN110197496A (zh) * | 2019-05-31 | 2019-09-03 | 上海联影医疗科技有限公司 | 一种医疗***控制方法、装置、设备及介质 |
CN110974372A (zh) * | 2020-01-03 | 2020-04-10 | 上海睿触科技有限公司 | 一种手术过程中病人运动位置实时跟踪装置 |
CN111462168B (zh) * | 2020-04-22 | 2023-09-19 | 上海联影医疗科技股份有限公司 | 运动参数估计方法和运动伪影校正方法 |
-
2021
- 2021-02-19 CN CN202110188324.2A patent/CN112545543B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112545543A (zh) | 2021-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112545543B (zh) | 基于病床运动信息的扫描运动监测方法、***及存储介质 | |
CN112336366B (zh) | 移动ct***的扫描重建方法、存储介质和ct扫描装置 | |
CN109498051B (zh) | 一种ct病床机架自动位置校准方法及其*** | |
US5251128A (en) | Motion artifact reduction in projection imaging | |
EP2146321B1 (en) | Dynamic error correction in radiographic imaging | |
US9538974B2 (en) | Methods and systems for correcting table deflection | |
US10638996B2 (en) | System and method for increasing the accuracy of a medical imaging device | |
JP2015526708A (ja) | 医用撮像のための動き補正 | |
WO2010015957A1 (en) | Automatic pre-alignment for registration of medical images | |
CN111657986A (zh) | 一种新型的移动ct中扫描机架驱动***及驱动方法 | |
CN114848001A (zh) | 一种超声宽景成像方法 | |
US20110075798A1 (en) | Method for correcting truncated projection data | |
CN112617877B (zh) | 移动ct***的自主扫描方法、存储介质和ct扫描装置 | |
US11954768B2 (en) | Method for gating in tomographic imaging system | |
JP7166964B2 (ja) | π線最適化を使用するコンピュータ断層撮影投影の幾何学的整合、被検体動き補正および強度正規化 | |
JP7011086B2 (ja) | X線システムの誤差トラッキングおよび較正 | |
EP4018215B1 (en) | Tomographic imaging with motion detection | |
JP6821839B1 (ja) | 二軸デジタルトモシンセシスシステムに用いられる幾何学的補正方法及びそのシステム | |
CN113570710A (zh) | 一种具有c型臂抖动校正的三维成像***及方法 | |
CN113229840B (zh) | 一种口腔cbct拍摄图像运动补偿重建方法 | |
WO2023196184A1 (en) | Pose-based three-dimensional structure reconstruction systems and methods | |
CN116128817A (zh) | 一种自动矫正医学图像中头部数据角度的方法及*** | |
JP2023040212A (ja) | 断層画像撮影システムにおけるゲーティング方法 | |
Lerner et al. | Vision-based tracking system for head motion correction in fMRI images | |
CN118261949A (zh) | 一种二维医学图像形变配准方法及*** |
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 |