一种新的构建序列影像与呼吸信号对应关系的方法技术

本发明专利技术公开了一种构建序列影像与呼吸信号对应关系的方法，在运动追踪法框架下，利用双目立体视觉系统，采集红外标记物的位移作为人体呼吸信号。采用新一代CT扫描机，获取受试者深呼吸模式下完整的肺部4D CT图像。结合体外标记物的距离移动信息，构造距离测度函数，提出距离对应关联法，通过标记物相同时间段内在4DCT中和呼吸信号中的位移量的比较，确定4D CT和呼吸信号的对应关系。该发明专利技术简化了呼吸信号获取的操作步骤，提高了数据测量精度，降低了实验时间以及受试者接收的辐射剂量。

【技术实现步骤摘要】
一种新的构建序列影像与呼吸信号对应关系的方法
本专利技术属于运动建模领域，具体涉及一种在影像导航手术中的构建序列影像与呼吸信号对应关系的方法，提高序列影像与呼吸信号对应的准确度。
技术介绍
影像导航手术中，呼吸运动是引起病灶定位不准确的重要因素之一。此外，当规划肺癌患者的放疗方案时，若缺乏对呼吸运动的认识，肿瘤会接收不到足够的剂量，周围的健康组织会接收不必要的剂量，或两者兼有。因此需要构建肺部组织的运动模型，模拟和预测呼吸运动对病灶的影响。为了构建肺部呼吸运动模型，需要同时获取4D肺部CT图像与呼吸信号，并确立两者之间的对应关系。如申请号为201410209850.2、名称为“肺部CT的三维配准方法”的专利中提出了一种三维配准方法，能够将同一位置的结节在不同CT数据中的成像位置进行自动对应，该方法首先进行刚性配准，通过互信息将CT数据中的肺部区域进行整体对齐，然后在刚性配准的基础上利用提取的匹配特征点对结合局部单点互信息(SMI)最大化的方式进行非刚性配准，得到最终的精确配准结果，然后利用配准结果提供的空间变换信息，进行相同结节成像位置的对应，从而方便结节的特征变化进行对比。获取呼吸信号时，一些方法通常采用运动追踪法获得，即在人体腹部贴上标记物，用摄像机在CT扫描过程中同步跟踪标记物，将得到的视频提取为单帧图像，再提取标记物质心，以质心坐标任意维度与时间的关系曲线作为呼吸曲线。如申请号为201410009192.2、名称为“下颌三维运动捕捉及可视化系统及方法”；申请号为201420011840.3、名称为“下颌三维运动捕捉及可视化系统与下颌运动追踪装置”的专利中均提出了三维运动捕捉及可视化系统，其标记物表面覆盖有荧光材料，且用红外摄像头带有红外发射装置，标记物可通过荧光材料反射出红外发射装置发出的红外线，以反映运动轨迹。获取4D肺部CT图像时，根据受试者呼吸模式的不同分为屏气模式和自由呼吸模式。屏气模式中，受试者在一个呼吸周期内多个时相点位置，采用屏气方式，同步获取CT图像与呼吸信号，由此可直接确定CT图像与呼吸信号的对应关系。自由呼吸模式中，受试者在自由呼吸状态下，进行CT图像和呼吸曲线的同步获取。该模式适用于任何受试者，在临床中具有广泛的应用前景。然而由于CT技术的限制，现有方法获取完整肺部的4DCT图像时，通常采用图像堆积扫描方式。扫描过程中，每个床位需要获取25或30甚至更多的图像。为了覆盖整个肺部区域，每个受试者需要扫描3到4个连续的床位，整个过程大概需要70秒。获得4D肺部CT图像和体外呼吸信号后，需要构建两者间的一一对应关系，即确定4D肺部CT图像中每一个3D肺部CT图像与体外信号在同一时刻获取，此处通常采用基于相位的关联法。利用高斯平滑滤波器对采集的呼吸曲线做平滑处理，滤除曲线上非对应于吸气末或呼气末的波峰和波谷；求取曲线上所有点的梯度，将梯度为零的点与该点附近最接近于受试者吸气末或呼气末状态的3D肺部CT对应；再将其他CT图像分布在完整的呼吸周期中。然而，上述的体外呼吸信号获取的操作方法均存在不足：1.获取体外呼吸信号的步骤繁琐，且在提取标记物质心过程中容易引入误差；2.获取CT图像时，在屏气模式下，不易控制获取数据时的呼吸状态，不适合呼吸功能障碍患者；而且屏气状态与自由呼吸状态下肺部组织运动存在差异性，导致对应误差；3.自由呼吸模式中，扫描时间长，可能引起受试者的不适，增加受试者接收的辐射剂量；4.构建序列图像与体外信号对应关系时，所用方法不能确保受试者吸气末或呼气末肺部组织的实时解剖结构与选取作为肺部组织吸气末或呼气末状态的3D肺部CT相一致，导致对应误差；5.同时，尽管通过该方法构建的肺部呼吸运动模型，可以估计呼吸周期中任意位置对应的当前的肺部状态，但是无法做到监测体外信号而实时的做出估计。
技术实现思路
为克服工业技术的限制，避免由于测量仪器以及构建序列图像与体外呼吸信号对应关系的方法的选取不足，本专利技术要研发的“一种新的构建序列影像与呼吸信号对应关系的方法”，利用立体视觉技术，简化体外呼吸信号获取的操作步骤，提高获取数据的精度；降低实验时间以及受试者接收的辐射剂量，避免采集4D肺部CT时受试者长时间屏气，克服受试者在屏气状态与自由呼吸状态下肺部活动存在差异性的问题；提高序列影像与呼吸信号对应的准确度。为了解决上述技术问题，实现上述目的，本专利技术通过以下技术方案实现：本专利技术获取呼吸信号的方法是立体视觉技术，即采用两个或两个以上的相机搭建的立体视觉测量系统，对放置在病人胸腹部的靶标进行实时空间(三维)定位，获取随时间变化的体表运动幅度曲线。本专利技术的技术方案的硬件主体包括五个部分，即红外标记物、立体视觉测量单元、数据分析单元、通信单元、控制单元；关键技术在于结合标记物在4DCT中与呼吸信号中的移位信息，构建距离测度函数，同时确定所有CT图像对应的标记物的空间坐标。该技术方案的工作流程是：使用CT扫描仪采集4D肺部CT图像，该CT扫描仪能够完成自主呼吸状态的肺部动态扫描。使用实验室自行研制的双目立体视觉系统采集体外呼吸信号，该系统可直接记录标记物的相对空间坐标，且精度高。数据采集时，受试者静卧在CT床上，将红外标记物黏贴于受试者胸部剑突位置的体表。首先开启双目立体视觉系统，采集红外标记物位移数据。接着，CT操作者通过麦克指导受试者作呼吸运动，接着启动CT扫描仪进行肺部动态扫描，实验场景如图1所示。通过上述方式获得了受试者自由深呼吸下的体外呼吸信号和完整肺部4DCT图像，结合标记物相同时间段内在4D肺部CT图像中和呼吸信号中的距离移动信息，构造距离测度函数，提出基于距离的对应关联法。通过计算距离测度函数可确定序列3D肺部CT图像与标记物空间坐标的一一对应关系，实验效果如图4所示。现有技术构建4DCT与体外呼吸信号对应关系的方法通常是采用高斯平滑滤波器对采集的呼吸曲线做平滑处理，滤除曲线上非对应于吸气末或呼气末的波峰和波谷；求取曲线上所有点的梯度，将梯度为零的点与该点附近最接近于受试者吸气末或呼气末状态的3D肺部CT对应，再将其他CT图像分布在完整的呼吸周期中。与现有技术相比，本专利技术的有益效果是提高了序列影像与呼吸信号对应的准确度，使得通过监测体外呼吸信号实时预测体内组织运动成为可能。上述说明仅是本专利技术技术方案的概述，为了能够更清楚了解本专利技术的技术手段，并可依照说明书的内容予以实施，以下以本专利技术的较佳实施例并配合附图详细说明如后。本专利技术的具体实施方式由以下实施例及其附图详细给出。本专利技术多处仅仅对做出改进的部分进行描述，而其他未说明部分可以借助本领域的现有技术实现，亦即未说明部分通过现有技术实现，在此不进行详细说明。附图说明此处所说明的附图用来提供对本专利技术的进一步理解，构成本申请的一部分，本专利技术的示意性实施例及其说明用于解释本专利技术，并不构成对本专利技术的不当限定。在附图中：图1为本专利技术的试验场景图。图1中，1和2为工业相机，3和4为红外滤光片，5和6为支架，7和8为数据连接线，9和10为相机电源线，11和12为采集卡，13为计算机，14为显示屏，16为红外LED，17为受试者，18为CT床，19为扫描架。图2为本专利技术的硬件主体工作流程图。图3为本专利技术的距离测度函数的序列的3DCT。图4为本本文档来自技高网...

【技术保护点】
一种新的构建序列影像与呼吸信号对应关系的方法，其特征在于：该方法的主体至少包括五个部分：红外标记物、立体视觉测量单元、数据分析单元、通信单元、控制单元；该方法结合标记物在4D CT中与呼吸信号中的移位信息，构建距离测度函数，同时确定所有CT图像对应的标记物的空间坐标，其步骤为：(1)、使用CT扫描仪采集4D肺部CT图像；该CT扫描仪能够完成自主呼吸状态的肺部动态扫描；(2)、使用立体视觉系统采集体外呼吸信号，直接记录标记物的相对空间坐标。

【技术特征摘要】
1.一种新的构建序列影像与呼吸信号对应关系的方法，其特征在于：该方法的主体至少包括五个部分：红外标记物、立体视觉测量单元、数据分析单元、通信单元、控制单元；该方法结合标记物在4DCT中与呼吸信号中的移位信息，构建距离测度函数，同时确定所有CT图像对应的标记物的空间坐标，其步骤为：(1)、使用CT扫描仪采集4D肺部CT图像；该CT扫描仪能够完成自主呼吸状态的肺部动态扫描；(2)、使用立体视觉系统采集体外呼吸信号，直接记录标记物的相对空间坐标；所述的构造的距离测度函数F的计算公式为：式中，r＝1，2，3代表x，y，z三个方向，qr为CT图像的空间分辨率；m-1为获取m个序列3D肺部CT图像I的时间间隔个数；Ci是标记物在第i个时间点的3D肺部CT图像Ii中的坐标，Dj为双目视觉系统记录的红外标记物第j个时间点的相对空间坐标。2.根据权利要求1所述的新的构建序列影像与呼吸信号对应关系的方法，其特征在于：所述红外标记物为一个120°投射角、直径5㎜的红外发光二极管，其支撑结构的底部具有容易放置或吸附在病人体表的结构，在病人进行呼吸运动时，红外标记物会随着体表进行运动；所述立体视觉测量单元由两个或两个以上的数字相机、计算机或单片机处理器、机械结构组成。3.根据权利要求2所述的新的构建序列影像与呼...

【专利技术属性】
技术研发人员：崔学理，高欣，夏威，
申请(专利权)人：中国科学院苏州生物医学工程技术研究所，
类型：发明
国别省市：江苏;32