CN114719760A - 一种动物成像分布式运动监测系统及位移监测方法 - Google Patents
一种动物成像分布式运动监测系统及位移监测方法 Download PDFInfo
- Publication number
- CN114719760A CN114719760A CN202210352475.1A CN202210352475A CN114719760A CN 114719760 A CN114719760 A CN 114719760A CN 202210352475 A CN202210352475 A CN 202210352475A CN 114719760 A CN114719760 A CN 114719760A
- Authority
- CN
- China
- Prior art keywords
- sensor
- data
- control module
- animal imaging
- motion
- 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.)
- Pending
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 92
- 238000003384 imaging method Methods 0.000 title claims abstract description 81
- 238000000034 method Methods 0.000 title claims abstract description 81
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 43
- 238000001514 detection method Methods 0.000 claims abstract description 9
- 230000001133 acceleration Effects 0.000 claims description 37
- 238000013139 quantization Methods 0.000 claims description 20
- 238000004891 communication Methods 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 16
- 230000005484 gravity Effects 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000003702 image correction Methods 0.000 claims description 7
- 230000001360 synchronised effect Effects 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 2
- 230000002093 peripheral effect Effects 0.000 claims description 2
- 238000002474 experimental method Methods 0.000 description 13
- 238000002591 computed tomography Methods 0.000 description 7
- 238000004590 computer program Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 238000012360 testing method Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000007667 floating Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 238000002603 single-photon emission computed tomography Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 238000013016 damping Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005188 flotation Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000029058 respiratory gaseous exchange Effects 0.000 description 2
- 238000011076 safety test Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000003748 differential diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000009509 drug development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 230000002503 metabolic effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003068 molecular probe Substances 0.000 description 1
- 239000002547 new drug Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 229940121896 radiopharmaceutical Drugs 0.000 description 1
- 239000012217 radiopharmaceutical Substances 0.000 description 1
- 230000002799 radiopharmaceutical effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
-
- 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P15/00—Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration
- G01P15/18—Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration in two or more dimensions
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P21/00—Testing or calibrating of apparatus or devices covered by the preceding groups
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2503/00—Evaluating a particular growth phase or type of persons or animals
- A61B2503/40—Animals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2503/00—Evaluating a particular growth phase or type of persons or animals
- A61B2503/42—Evaluating a particular growth phase or type of persons or animals for laboratory research
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- High Energy & Nuclear Physics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本申请涉及活体动物成像领域，本申请提供一种动物成像分布式运动监测系统及位移监测方法，所述系统包括：运动传感器，所述运动传感器分布于X射线源、探测器以及载物台处，用于将检测的数据传输到后端的控制模块；控制模块，所述控制模块用于接收并记录所述运动传感器传输的检测数据，并将记录的数据发送至上位机；上位机，所述上位机用于接收所述控制模块记录的数据，根据接收的数据计算出位移值，并进行图像校正，得到校正图像。实现高精度的动物成像位移监测，实现了系统的位置偏差估算，为高精度的动物成像提供了硬件基础。
Description
技术领域
本申请涉及活体动物成像领域，尤其涉及一种动物成像分布式运动监测系统及位移监测方法。
背景技术
动物SPECT/CT是一种分辨率极高的活体分子代谢成像技术，已广泛用于多种疾病的鉴别诊断、疗效评价、及新药开发等方面。其检测原理通过单光子发射计算机断层成像术(Single-Photon Emission Computed Tomography，SPECT)检测生物体内的标记分子探针的信息，从而得到放射性药物在生物体内的三维分布，电子计算机断层扫描(ComputedTomography，CT)提供高精度的物理结构图像。近些年来随着SPECT/CT设备空间性能的提升(尤其是动物CT空间分辨率的提升，已经达到了um量级)，以及应用环境的日益复杂，对系统的位置稳定性提出了进一步的要求。
与此同时，在活体动物成像领域，随着设备的空间分辨率的提升，由于被测目标的呼吸，心跳等带来的微小位移所带来的图像伪影越来越严重，为保证采集数据的质量，需要对平台进行减震处理，以减小运动伪影的产生。因实验室场地和条件限制，往往多个实验系统需要共用设备和场地，无法做到实验人员、设备，在空间上的隔离，可能对成像结果造成干扰。因此，需要对实验平台的振动情况进行检测，使用振动信息对采集的数据进行修正与取舍以进一步提高图像质量。
发明内容
为了解决上述现有技术中存在的技术问题,本申请的主要目的在于提供一种动物成像分布式运动监测系统及位移监测方法，实现了一种快速的多维度位移监测及校准算法，大大提高了系统抗震及微小运动检测的能力，解决了目前设计中由于生物体心跳，呼吸以及实验平台振动所带来的严重运动伪影的问题，提高了系统的空间分辨能力。
一种动物成像分布式运动监测系统及位移监测方法。
为实现上述目的，本申请提供一种动物成像分布式运动监测系统，包括：
运动传感器，所述运动传感器分布于X射线源、探测器以及载物台处，用于将检测的数据传输到后端的控制模块；
控制模块，所述控制模块用于接收并记录所述运动传感器传输的检测数据，并将记录的数据发送至上位机；
上位机，所述上位机用于接收所述控制模块记录的数据，根据接收的数据计算出位移值，并进行图像校正，得到校正图像。
此外，为实现上述目的，本申请还提供一种动物成像位移监测方法，包括以下步骤：
S1、系统初始化
进行定时器初始化、通信接口初始化、SPI初始化和运动传感器初始化；
S2、利用重力加速度对三轴的量化系数进行传感器校准
S201、将所述传感器水平放置，至X方向与重力方向平行，读取三轴的数值；
S202、若Y方向和Z方向的数值满足限制条件，则记录此时Xi的值；
S203、多次重复步骤S201和步骤S202后获得
轴的平均数值，则X轴的量化系数为：
S204、采用上述步骤S201至步骤S203同样的方法，获取ky的值与kz的值,得到Y轴的量化系数和Z轴的量化系数；
S3、记录初始状态的信息，形成初始成像空间；
S4、进行采集数据，并监测各传感器加速度，实验过程中的加速度监测值为实测值与初始值之差；
S5、主控模块通过同步电路将各传感器加速度与探测器数据信号进行同步，将记录的数据发送至上位机，上位机计算出位移值，并进行图像校正，得到校正图像。
此外，为实现上述目的，本申请还提供一种动物成像位置监控的校准设备，所述动物成像位置监控的校准设备包括处理器、存储器、以及存储在所述存储器上并可被所述处理器执行的动物成像位置监控的校准程序，其中所述动物成像位置监控的校准程序被所述处理器执行时，实现如上述的动物成像位移监测方法的步骤。
此外，为实现上述目的，本申请还提供一种计算机可读存储介质，所述计算机可读存储介质上存储有动物成像位置监控的校准程序，其中所述动物成像位置监控的校准程序被处理器执行时，实现如上述的动物成像位移监测方法的步骤。
本申请提供一种动物成像分布式运动监测系统及位移监测方法，利用传感器具有体积小、价格低的优势，可方便集成在监测系统中，实现对多目标或多部位的振动情况监测分析，解决了在复杂应用环境下由于微小震动造成的系统空间分辨率降低的问题。通过分布式运动检测实现高精度的动物成像位移监测特别适合高精度CT系统实现；在硬件端实现了系统的位置偏差估算，从而大幅度减少了由于物理震动带来的运动偏差，为高精度的动物成像提供了硬件基础；本申请为在精密试验、安全试验、破坏分析领域具有极大的应用价值，为研发高精度成像系统提供有效的技术保证。
本申请的这些方面或其他方面在以下实施例的描述中会更加简明易懂。应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本申请。
附图说明
为了更清楚地说明本申请实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例。在附图中：
图1为本申请第一实施例中动物成像位移监测方法的系统框图；
图2为图1中本申请动物成像分布式运动监测系统中系统连接方式示意图；
图3为本申请动物成像分布式运动监测系统中数据采集的流程示意图；
图4为本申请第二实施例中动物成像位移监测方法的流程示意图；
图5为图4中本申请动物成像位移监测方法中传感器校准的流程示意图；
图6为本申请第二实施例中动物成像位移监测方法中静态下两个实验台振动情况监测对比图；
图7为本申请第二实施例中动物成像位移监测方法中扰动下两个实验台振动情况监测对比图；
图8为本申请第三实施例中动物成像位置监控的校准设备的设备示意图。
本申请目的实现、功能特点及优点将结合实施例，参照附图做进一步说明。
具体实施方式
本部分将详细描述本申请的具体实施例，本申请之较佳实施例在附图中示出，附图的作用在于用图形补充说明书文字部分的描述，使人能够直观地、形象地理解本申请的每个技术特征和整体技术方案，但其不能理解为对本申请保护范围的限制。
在本申请的描述中，若干的含义是一个或者多个，多个的含义是两个以上，大于、小于、超过等理解为不包括本数，以上、以下、以内等理解为包括本数。
在本申请的描述中，对方法步骤的连续标号是为了方便审查和理解，结合本申请的整体技术方案以及各个步骤之间的逻辑关系，调整步骤之间的实施顺序并不会影响本申请技术方案所达到的技术效果。
本申请的描述中，除非另有明确的限定，设置等词语应做广义理解，所属技术领域技术人员可以结合技术方案的具体内容合理确定上述词语在本申请中的具体含义。
在本申请的实施方式中，提供了一种动物成像分布式运动监测系统及位移监测方法，其中，动物成像分布式运动监测系统，可以应用在动物成像位置监控的校准设备上，该动物成像位置监控的校准设备可以是PC、便携计算机、移动终端等具有显示和处理功能的设备，当然也不限于此。
其中，参见图1所示，一种动物成像分布式运动监测系统，包括运动传感器300、控制模块200以及上位机100。在本实施例中，所述运动传感器300分布于X射线源、探测器以及载物台处，并通过等长的数据线连接到后端的控制模块200，用于将检测的数据传输到后端的控制模块。在一些实施例中参见图2所示，所述运动传感器300采用六轴运动传感器，例如：BMI160，所述运动传感器300与所述控制模块200的数字接口之间采用串行外设接口(Serial Peripheral Interface，SPI)以串行方式进行通信以交换信息。
参见图1所示，所述控制模块200用于接收并记录所述运动传感器传输的检测数据，并将记录的数据发送至上位机。所述控制模块200包括FPGA芯片、时钟同步模块、数字接口及所述数字接口的通信协议，所述数字接口用于接收时钟同步信号、数据传输信号、触发信号、状态反馈信号以及地信号，还用于连接电源。在一些实施例中，参见图2所示，所述控制模块200内还包括型号为赛灵思ZYNQ7000系列的FPGA芯片，所述FPGA芯片与所述上位机100之间采用采用RS-232标准串行通信接口或网线连接。
在一些实施例中，所述控制模块200的数字接口还包括串行通信与时钟同步接口，支持多机通信模式。
在一些实施例中，所述控制模块200内设有接口电路、同步电路以及主控制电路，所述主控制电路通过FPGA芯片实现，电路进行同步采集之前，还包括分布式运动监测系统的系统初始化和传感器校准。
在一些实施例中，所述系统初始化包括定时器初始化、通信接口初始化、SPI初始化和运动传感器初始化四部分。
参见图1所示，所述上位机100用于接收所述控制模块200记录的数据，根据接收的数据计算出位移值，并进行图像校正，得到校正图像。
在一些实施例中，参见图3所示，同步电路其采集流程如下：
1、系统初始化过程包括定时器初始化，通信接口初始化、SPI初始化和运动传感器初始化四部分。
2、初始化后，考虑到传感器的差异，利用重力加速度对三轴的量化系数进行校准。具体方法是，先将传感器水平放置，保证X方向与重力方向平行，读取三轴的数值，若Y，Z的数值满足限制条件，则记录此时Xi的值。为减小误差，多次重复后获得
则X轴的量化系数为：
然后，再用同样的方法，获取ky与kz。
3、考虑到实验台实验台会略有倾斜，因此需要记录初始状态的信息，并形成初始成像空间。
4、然后系统开始采集相关数据，并监测各传感器加速度，实验过程中的加速度监测值为实测值与初始值之差。
5、由FPGA构成的主控模块通过同步电路将各传感器加速度与探测器数据信号进行同步。
6、将记录的数据发送至上位机100，计算出位移值，并进行图像校正，得到校正图像。
在一些实施例中，所述系统初始化后，进行传感器校准，所述传感器校准为利用重力加速度对三轴的量化系数进行校准，所述传感器校准的方法包括：
1)将所述传感器水平放置，至X方向与重力方向平行，读取三轴的数值；
2)若Y方向和Z方向的数值满足限制条件，则记录此时Xi的值；
3)多次重复步骤1)和步骤2)后获得
轴的平均数值，则X轴的量化系数为：
4)采用上述步骤1)至步骤3)同样的方法，获取ky的值与kz的值,得到Y轴的量化系数和Z轴的量化系数。
在一些实施例中，所述运动传感器300分布的X射线源、探测器以及载物台均设置在实验台实验台上，传感器校准之后还包括记录初始状态的信息，形成初始成像空间，同步电路后进行采集。
在一些实施例中，在系统同步进行采集数据时，还包括监测所述传感器加速度，其中，实验过程中的加速度监测值为实测值与初始值之差；所述主控模块还用于通过同步电路将各传感器加速度与探测器数据信号进行同步至上位机100。
在一些实施例中，所述上位机100还用于接收所述主控模块记录的数据，根据接收的数据计算出位移值，并进行图像校正，得到校正图像。
在一些实施例中，所述数字接口包含有系统同步时钟接口，所述系统同步电路主时钟位于控制模块200，并产生系统全局时钟并分发到各个运动传感器和探测器中保证探测器数据和运动数据的同步性控制模块时钟同步电路各个运动传感器和探测器。
在本实施例中，同步电路包含了同步时钟接口，系统电源以及数字核心控制模块200。整个时钟系统位于控制模块，并由时钟同步电路产生系统全局时钟并分发到各个运动传感器和探测器中保证探测器数据和运动数据的同步性。
下面结合附图，对本申请的一些实施方式作详细说明。在不冲突的情况下，下述的实施例及实施例中的特征可以相互组合。
请参照图4，图4为本申请第二实施例的一种动物成像位移监测方法的流程示意图。本申请的实施例中，所述一种动物成像位移监测方法包括以下步骤S1～S5：
S1、系统初始化
进行定时器初始化、通信接口初始化、SPI初始化和运动传感器初始化；
S2、利用重力加速度对三轴的量化系数进行传感器校准；
S3、记录初始状态的信息，形成初始成像空间；
S4、进行采集数据，并监测各传感器加速度，实验过程中的加速度监测值为实测值与初始值之差；
S5、主控模块通过同步电路将各传感器加速度与探测器数据信号进行同步，将记录的数据发送至上位机100，上位机100计算出位移值，并进行图像校正，得到校正图像。
在本申请的实施例中，参见图5所示，在利用重力加速度对三轴的量化系数进行传感器校准阶段，还包括步骤S201～S204：
S201、将所述传感器水平放置，至X方向与重力方向平行，读取三轴的数值；
S202、若Y方向和Z方向的数值满足限制条件，则记录此时Xi的值；
S203、多次重复步骤S201和步骤S202后获得
轴的平均数值，则X轴的量化系数为：
S204、采用上述步骤S201至步骤S203同样的方法，获取ky的值与kz的值,得到Y轴的量化系数和Z轴的量化系数。
在本申请的实施例中，在所述传感器校准阶段，校准方法，包括：
将X射线源，探测器以及载物台处的振动探测器分别定义为Detector1、Detector2和Detector3，将Detector1、Detector2和Detector3采集到的数据矩阵进行低通滤波，并分别定义为D1、D2、D3；
根据T0时刻的D1、D2的数值建立成像坐标系P(X,Y,Z)；
根据T时刻的D1、D2数值计算探测器的相对位置，建立修正坐标系为P(X’,Y’,Z’)；其中相对位置公式计算方式为：
X为D1探测器T时刻的位置，x(t)为D1探测器加速度计x方向的分量的函数，以此类推Y,Z位置的计算；
根据D3的数值计算转台的实际偏差，其中，D3的方向偏差计算为：
其中，X,Y,Z为R为探测器到旋转中心的距离，ω为角速度。
示例性的，以下为一个应用实例用所设计的传感器进行了实验，实验台实验台并对普通实验台和气浮实验台进行了振动误差判断，验证设计的有效性。
参见图6和图7所示，实线为气浮平台实验结果，虚线为普通实验台结果。Y轴为Z方向加速度，单位是mm/s2,X轴为时间，单位是s。上半部分为静止状态下，两个实验台振动情况监测对比图，下半部分是在30秒后，启动偏心振动电机后观察外界扰动对实验台的影响。
按静止状态，施加扰动前和施加扰动分类后，计算了加速度数据的标准差，极值和波动范围，通过统计分析，可以清晰的看出，在无扰动情况下，气浮实验台和普通实验台加速度的标准差、极值和波动范围都比较接近，气浮实验台的性能略优于普通实验台。但在施加扰动后，由于气浮平台具有较大的柔性减震阻尼，其加速的标准差、极值和波动范围均远小于普通实验台。通过该实验可以看到，检测系统可以实时检测到系统各部分的位移状态，从而进行图像运动校准，该方法也可以为系统的运动部件选型提供参考。
其中，上述动物成像分布式运动监测系统中各个模块与上述动物成像分布式运动监测方法实施例中各步骤相对应，其功能和实现过程在此处不再一一赘述。
参照图8，为根据本申请方式的设备图。本申请第三实施例提供了一种动物成像位置监控的校准设备包括存储器600及处理器700，其中处理器700存储有计算机程序，计算机程序用于执行上述动物成像位移监测方法实施例中各步骤：
S1、系统初始化
进行定时器初始化、通信接口初始化、SPI初始化和运动传感器初始化；
S2、利用重力加速度对三轴的量化系数进行传感器校准
S201、将所述传感器水平放置，至X方向与重力方向平行，读取三轴的数值；
S202、若Y方向和Z方向的数值满足限制条件，则记录此时Xi的值；
S203、多次重复步骤S201和步骤S202后获得
轴的平均数值，则X轴的量化系数为：
S204、采用上述步骤S201至步骤S203同样的方法，获取ky的值与kz的值,得到Y轴的量化系数和Z轴的量化系数；
S3、记录初始状态的信息，形成初始成像空间；
S4、进行采集数据，并监测各传感器加速度，实验过程中的加速度监测值为实测值与初始值之差；
S5、主控模块通过同步电路将各传感器加速度与探测器数据信号进行同步，将记录的数据发送至上位机100，上位机100计算出位移值，并进行图像校正，得到校正图像。
其中，所述传感器校准，校准方法，包括：
将X射线源，探测器以及载物台处的振动探测器分别定义为Detector1、Detector2和Detector3，将Detector1、Detector2和Detector3采集到的数据矩阵进行低通滤波，并分别定义为D1、D2、D3；
根据T0时刻的D1、D2的数值建立成像坐标系P(X,Y,Z)；
根据T时刻的D1、D2数值计算探测器的相对位置，建立修正坐标系为P(X’,Y’,Z’)；其中相对位置公式计算方式为：
X为D1探测器T时刻的位置，x(t)为D1探测器加速度计x方向的分量的函数，以此类推Y,Z位置的计算；
根据D3的数值计算转台的实际偏差，其中，D3的方向偏差计算为：
其中，X,Y,Z为R为探测器到旋转中心的距离，ω为角速度。
其中，动物成像位置监控的校准程序被执行时所实现的方法可参照本申请动物成像分布式运动监测系统的各个实施例，此处不再赘述。
应当理解的是，处理器可以是中央处理单元(Central Processing Unit，CPU)，该处理器还可以是其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现场可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。其中，通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
其中，所述处理器用于运行存储在存储器中的计算机程序，以实现本申请动物成像位移监测方法的各个实施例。
应当认识到，本申请实施例中的方法步骤可以由计算机硬件、硬件和软件的组合、或者通过存储在非暂时性计算机可读存储器中的计算机指令来实现或实施。所述方法可以使用标准编程技术。每个程序可以以高级过程或面向对象的编程语言来实现以与计算机系统通信。然而，若需要，该程序可以以汇编或机器语言实现。在任何情况下，该语言可以是编译或解释的语言。此外，为此目的该程序能够在编程的专用集成电路上运行。
此外，可按任何合适的顺序来执行本文描述的过程的操作，除非本文另外指示或以其他方式明显地与上下文矛盾。本文描述的过程(或变型和/或其组合)可在配置有可执行指令的一个或多个计算机系统的控制下执行，并且可作为共同地在一个或多个处理器上执行的代码(例如，可执行指令、一个或多个计算机程序或一个或多个应用)、由硬件或其组合来实现。所述计算机程序包括可由一个或多个处理器执行的多个指令。
此外，本申请第四实施例还提供一种计算机可读存储介质。
本申请计算机可读存储介质上存储有动物成像位置监控的校准程序，其中所述动物成像位置监控的校准程序被处理器执行时，实现如上述的动物成像位移监测方法的步骤：
S1、系统初始化
进行定时器初始化、通信接口初始化、SPI初始化和运动传感器初始化；
S2、利用重力加速度对三轴的量化系数进行传感器校准；
S201、将所述传感器水平放置，至X方向与重力方向平行，读取三轴的数值；
S202、若Y方向和Z方向的数值满足限制条件，则记录此时Xi的值；
S203、多次重复步骤S201和步骤S202后获得
轴的平均数值，则X轴的量化系数为：
S204、采用上述步骤S201至步骤S203同样的方法，获取ky的值与kz的值,得到Y轴的量化系数和Z轴的量化系数；
S3、记录初始状态的信息，形成初始成像空间；
S4、进行采集数据，并监测各传感器加速度，实验过程中的加速度监测值为实测值与初始值之差；
S5、主控模块通过同步电路将各传感器加速度与探测器数据信号进行同步，将记录的数据发送至上位机100，上位机100计算出位移值，并进行图像校正，得到校正图像。
其中，所述传感器校准，校准方法，包括：
将X射线源，探测器以及载物台处的振动探测器分别定义为Detector1、Detector2和Detector3，将Detector1、Detector2和Detector3采集到的数据矩阵进行低通滤波，并分别定义为D1、D2、D3；
根据T0时刻的D1、D2的数值建立成像坐标系P(X,Y,Z)；
根据T时刻的D1、D2数值计算探测器的相对位置，建立修正坐标系为P(X’,Y’,Z’)；其中相对位置公式计算方式为：
X为D1探测器T时刻的位置，x(t)为D1探测器加速度计x方向的分量的函数，以此类推Y,Z位置的计算；
根据D3的数值计算转台的实际偏差，其中，D3的方向偏差计算为：
其中，X,Y,Z为R为探测器到旋转中心的距离，ω为角速度。
其中，动物成像位置监控的校准程序被执行时所实现的方法可参照本申请动物成像分布式运动监测系统的各个实施例，此处不再赘述。
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。
上述本申请实施例序号仅仅为了描述，不代表实施例的优劣。
本申请可用于众多通用或专用的计算机装置环境或配置中。进一步，所述方法可以在可操作地连接至合适的任何类型的计算平台中实现，包括但不限于个人电脑、迷你计算机、主框架、工作站、网络或分布式计算环境、单独的或集成的计算机平台、或者与带电粒子工具或其它成像装置通信等等。本申请的各方面可以以存储在非暂时性存储介质或设备上的机器可读代码来实现，无论是可移动的还是集成至计算平台，如硬盘、光学读取和/或写入存储介质、RAM、ROM等，使得其可由可编程计算机读取，当存储介质或设备由计算机读取时可用于配置和操作计算机以执行在此所描述的过程。此外，机器可读代码，或其部分可以通过有线或无线网络传输。当此类媒体包括结合微处理器或其他数据处理器实现上文所述步骤的指令或程序时，本文所述的申请包括这些和其他不同类型的非暂时性计算机可读存储介质。当根据本申请所述的方法和技术编程时，本申请还包括计算机本身。
计算机程序能够应用于输入数据以执行本文所述的功能，从而转换输入数据以生成存储至非易失性存储器的输出数据。输出信息还可以应用于一个或多个输出设备如显示器。在本申请优选的实施例中，转换的数据表示物理和有形的对象，包括显示器上产生的物理和有形对象的特定视觉描绘。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在如上所述的一个存储介质(如ROM/RAM、磁碟、光盘)中，包括若干指令用以使得一台终端设备(可以是手机，计算机，服务器，空调器，或者网络设备等)执行本申请各个实施例所述的方法。
本申请提供一种动物成像分布式运动监测系统及位移监测方法，利用传感器具有体积小、价格低的优势，可方便集成在监测系统中，实现对多目标或多部位的振动情况监测分析，解决了在复杂应用环境下由于微小震动造成的系统空间分辨率降低的问题。通过分布式运动检测实现高精度的动物成像位移监测特别适合高精度CT系统实现；在硬件端实现了系统的位置偏差估算，从而大幅度减少了由于物理震动带来的运动偏差，为高精度的动物成像提供了硬件基础；本申请为在精密试验、安全试验、破坏分析领域具有极大的应用价值，为研发高精度成像系统提供有效的技术保证。
以上仅为本申请的优选实施例，并非因此限制本申请的专利范围，凡是利用本申请说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本申请的专利保护范围内。
Claims (10)
1.一种动物成像分布式运动监测系统，其特征在于，所述动物成像分布式运动监测系统包括：
运动传感器，所述运动传感器分布于X射线源、探测器以及载物台处，用于将检测的数据传输到后端的控制模块；
控制模块，所述控制模块用于接收并记录所述运动传感器传输的检测数据，并将记录的数据发送至上位机；以及
上位机，所述上位机用于接收所述控制模块记录的数据，根据接收的数据计算出位移值，并进行图像校正，得到校正图像。
2.根据权利要求1所述的动物成像分布式运动监测系统，其特征在于，所述述运动传感器采用六轴运动传感器，所述运动传感器与所述控制模块的数字接口之间采用SPI以串行方式进行通信以交换信息；
所述控制模块内还包括FPGA芯片，所述FPGA芯片与所述上位机之间采用串行通信接口或网线连接。
3.根据权利要求2所述的动物成像分布式运动监测系统，其特征在于，所述控制模块的数字接口还包括串行通信与时钟同步接口，支持多机通信模式。
4.根据权利要求3所述的动物成像分布式运动监测系统，其特征在于，所述控制模块内设有接口电路、同步电路以及主控制电路，所述主控制电路通过FPGA芯片实现，电路进行同步采集之前，还包括分布式运动监测系统的系统初始化和传感器校准。
5.根据权利要求4所述的动物成像分布式运动监测系统，其特征在于，所述系统初始化包括定时器初始化、通信接口初始化、SPI初始化和运动传感器初始化四部分；系统初始化后，进行传感器校准，所述传感器校准为利用重力加速度对三轴的量化系数进行校准，所述传感器校准的方法包括：
1)将所述传感器水平放置，至X方向与重力方向平行，读取三轴的数值；
2)若Y方向和Z方向的数值满足限制条件，则记录此时Xi的值；
3)多次重复步骤1)和步骤2)后获得X一轴的平均数值，则X轴的量化系数为：
4)采用上述步骤1)至步骤3)同样的方法，获取ky的值与kz的值,得到Y轴的量化系数和Z轴的量化系数。
6.根据权利要求5所述的动物成像分布式运动监测系统，其特征在于，所述运动传感器分布的X射线源、探测器以及载物台均设置在实验台实验台上，传感器校准之后还包括记录初始状态的信息，形成初始成像空间，在系统同步进行采集；
进行采集数据时，还包括监测所述传感器加速度，其中，实验过程中的加速度监测值为实测值与初始值之差；所述主控模块还用于通过同步电路将各传感器加速度与探测器数据信号进行同步至上位机。
7.根据权利要求6所述的动物成像分布式运动监测系统，其特征在于，所述上位机还用于接收所述主控模块记录的数据，根据接收的数据计算出位移值，并进行图像校正，得到校正图像。
8.根据权利要求3所述的动物成像分布式运动监测系统，其特征在于，所述数字接口内包含有同步时钟接口，系统同步时钟主时钟位于控制模块，并产生系统全局时钟并分发到各个运动传感器和探测器中保证探测器数据和运动数据的同步性控制模块时钟同步电路各个运动传感器和探测器。
9.一种动物成像位移监测方法，其特征在于，所述动物成像位移监测方法包括如权利要求1-8中任一项所述的动物成像分布式运动监测系统，其特征在于，该方法包括：
S1、系统初始化
进行定时器初始化、通信接口初始化、SPI通信接口初始化、SPI初始化和运动传感器初始化；
S2、利用重力加速度对三轴的量化系数进行传感器校准；
S3、记录初始状态的信息，形成初始成像空间；
S4、进行采集数据，并监测各传感器加速度，实验过程中的加速度监测值为实测值与初始值之差；
S5、主控模块通过同步电路将各传感器加速度与探测器数据信号进行同步，将记录的数据发送至上位机，上位机计算出位移值，并进行图像校正，得到校正图像。
10.根据权利要求9所述的动物成像位移监测方法，其特征在于，所述传感器校准，校准方法，包括：
将X射线源，探测器以及载物台处的振动探测器分别定义为Detector1、Detector2和Detector3，将Detector1、Detector2和Detector3采集到的数据矩阵进行低通滤波，并分别定义为D1、D2、D3；
根据T0时刻的D1、D2的数值建立成像坐标系P(X,Y,Z)；
根据T时刻的D1、D2数值计算探测器的相对位置，建立修正坐标系为P(X’,Y’,Z’)；其中相对位置公式计算方式为：
X为D1探测器T时刻的位置，x(t)为D1探测器加速度计x方向的分量的函数，以此类推Y,Z位置的计算；
根据D3的数值计算转台的实际偏差，其中，D3的方向偏差计算为：
其中，X,Y,Z为R为探测器到旋转中心的距离，ω为角速度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210352475.1A CN114719760A (zh) | 2022-04-05 | 2022-04-05 | 一种动物成像分布式运动监测系统及位移监测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210352475.1A CN114719760A (zh) | 2022-04-05 | 2022-04-05 | 一种动物成像分布式运动监测系统及位移监测方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114719760A true CN114719760A (zh) | 2022-07-08 |
Family
ID=82242755
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210352475.1A Pending CN114719760A (zh) | 2022-04-05 | 2022-04-05 | 一种动物成像分布式运动监测系统及位移监测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114719760A (zh) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108037145A (zh) * | 2017-11-28 | 2018-05-15 | 北京纳米维景科技有限公司 | 一种用于计算机断层扫描成像系统的测试平台及测试方法 |
| CN109008972A (zh) * | 2013-02-01 | 2018-12-18 | 凯内蒂科尔股份有限公司 | 生物医学成像中的实时适应性运动补偿的运动追踪系统 |
| CN109141433A (zh) * | 2018-09-20 | 2019-01-04 | 江阴市雷奥机器人技术有限公司 | 一种机器人室内定位系统及定位方法 |
| JP2019180941A (ja) * | 2018-04-12 | 2019-10-24 | キヤノンメディカルシステムズ株式会社 | Ｘ線診断装置、画像処理装置、画像処理システム、画像処理プログラム、放射線診断装置、及び再構成用データのデータ構造 |
| CN110720916A (zh) * | 2019-08-26 | 2020-01-24 | 南京睿蜂健康医疗科技有限公司 | 一种便携式呼吸门控系统及方法 |
| CN111096759A (zh) * | 2018-10-26 | 2020-05-05 | 深圳迈瑞生物医疗电子股份有限公司 | X射线摄影系统及其平板探测器、相关方法 |
| CN111132618A (zh) * | 2017-09-26 | 2020-05-08 | 皇家飞利浦有限公司 | 移动式x射线成像系统 |
| CN111789609A (zh) * | 2020-08-10 | 2020-10-20 | 上海联影医疗科技有限公司 | 一种成像系统和方法 |
| CN113100794A (zh) * | 2021-03-26 | 2021-07-13 | 深圳市深图医学影像设备有限公司 | 一种x射线平板探测器运动伪影去除方法及其装置 |
| CN215305855U (zh) * | 2021-07-07 | 2021-12-28 | 重庆理工大学 | 一种心率监测装置 |
-
2022
- 2022-04-05 CN CN202210352475.1A patent/CN114719760A/zh active Pending
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109008972A (zh) * | 2013-02-01 | 2018-12-18 | 凯内蒂科尔股份有限公司 | 生物医学成像中的实时适应性运动补偿的运动追踪系统 |
| CN111132618A (zh) * | 2017-09-26 | 2020-05-08 | 皇家飞利浦有限公司 | 移动式x射线成像系统 |
| CN108037145A (zh) * | 2017-11-28 | 2018-05-15 | 北京纳米维景科技有限公司 | 一种用于计算机断层扫描成像系统的测试平台及测试方法 |
| JP2019180941A (ja) * | 2018-04-12 | 2019-10-24 | キヤノンメディカルシステムズ株式会社 | Ｘ線診断装置、画像処理装置、画像処理システム、画像処理プログラム、放射線診断装置、及び再構成用データのデータ構造 |
| CN109141433A (zh) * | 2018-09-20 | 2019-01-04 | 江阴市雷奥机器人技术有限公司 | 一种机器人室内定位系统及定位方法 |
| CN111096759A (zh) * | 2018-10-26 | 2020-05-05 | 深圳迈瑞生物医疗电子股份有限公司 | X射线摄影系统及其平板探测器、相关方法 |
| CN110720916A (zh) * | 2019-08-26 | 2020-01-24 | 南京睿蜂健康医疗科技有限公司 | 一种便携式呼吸门控系统及方法 |
| CN111789609A (zh) * | 2020-08-10 | 2020-10-20 | 上海联影医疗科技有限公司 | 一种成像系统和方法 |
| CN113100794A (zh) * | 2021-03-26 | 2021-07-13 | 深圳市深图医学影像设备有限公司 | 一种x射线平板探测器运动伪影去除方法及其装置 |
| CN215305855U (zh) * | 2021-07-07 | 2021-12-28 | 重庆理工大学 | 一种心率监测装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3793266B2 (ja) | ポジトロンｃｔ装置およびその画像再構成方法 | |
| Goodvin et al. | Development of a real-time three-dimensional spinal motion measurement system for clinical practice | |
| Picerno et al. | A spot check for assessing static orientation consistency of inertial and magnetic sensing units | |
| CN105496436B (zh) | 用于pet装置的时间校正方法和装置 | |
| US20150297120A1 (en) | Method For Tracking Motion of Subject in Real Time and for Correcting Medical Image | |
| US20110081058A1 (en) | Reconstructing Three-Dimensional Current Sources from Magnetic Sensor Data | |
| McGrath et al. | Body-worn IMU-based human hip and knee kinematics estimation during treadmill walking | |
| US11717184B2 (en) | Tracking head motion for medical imaging | |
| WO2019121127A1 (en) | Combining image based and inertial probe tracking | |
| Cutti et al. | A simple test to assess the static and dynamic accuracy of an inertial sensors system for human movement analysis | |
| Defrise et al. | Perturbative refinement of the geometric calibration in pinhole SPECT | |
| Paloschi et al. | Validation and assessment of a posture measurement system with magneto-inertial measurement units | |
| CN108634976B (zh) | 一种图像校正方法和装置 | |
| CN104076309A (zh) | 用于磁共振图像数据组校正的方法及相应磁共振设备 | |
| CN114719760A (zh) | 一种动物成像分布式运动监测系统及位移监测方法 | |
| Lehtonen et al. | A respiratory motion estimation method based on inertial measurement units for gated positron emission tomography | |
| Bellitti et al. | Smart brace for static and dynamic knee laxity measurement | |
| Ilewicz et al. | Direct method of IMU calibration | |
| Barraza-Madrigal et al. | Instantaneous position and orientation of the body segments as an arbitrary object in 3D space by merging gyroscope and accelerometer information | |
| RU2679296C1 (ru) | Сейсмокардиоблок и способ измерения сейсмокардиоцикла | |
| JP2018520815A (ja) | 方位トラッキング・デバイス | |
| van Niekerk et al. | A method for measuring orientation within a magnetic resonance imaging scanner using gravity and the static magnetic field (VectOrient) | |
| Bhowmik et al. | A head motion measurement system suitable for 3D cone-beam tomography using markers | |
| Sgrosso et al. | Torso Motion Monitoring with an IMU Set-Up | |
| CN105575239B (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 |