CN112085833A - Analysis method for in-vivo three-dimensional motion of cervical vertebra by combining cone beam CT and image fusion - Google Patents
Analysis method for in-vivo three-dimensional motion of cervical vertebra by combining cone beam CT and image fusion Download PDFInfo
- Publication number
- CN112085833A CN112085833A CN202010853773.XA CN202010853773A CN112085833A CN 112085833 A CN112085833 A CN 112085833A CN 202010853773 A CN202010853773 A CN 202010853773A CN 112085833 A CN112085833 A CN 112085833A
- Authority
- CN
- China
- Prior art keywords
- cervical vertebra
- image fusion
- model
- cone beam
- change
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/50—Depth or shape recovery
- G06T7/55—Depth or shape recovery from multiple images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/73—Determining position or orientation of objects or cameras using feature-based 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]
-
- 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/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30008—Bone
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Multimedia (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The invention discloses a method for analyzing the in-vivo three-dimensional motion of cervical vertebra by combining cone beam CT and image fusion, which comprises the following steps: s1, shooting images of different cervical vertebra functional positions by a cone beam CT machine; s2, importing images for segmentation and 3D modeling; s3, determining a local coordinate system for the 3D cervical vertebra model in the neutral position according to the anatomical position of the human body; s4, fusing images and recording the change of space position parameters of the cervical vertebra models of the segments with different functional positions relative to the neutral position; s5, performing threshold judgment on the accuracy of image fusion to determine whether the image fusion is successful; s6, further calculating the space position parameter change between different sections to obtain the relative change of the two sections, relating to the technical field of CT influence. The invention solves the problem that the conventional X-ray imaging for perspective shooting is not clear; the dynamic motion abnormality cannot be seen by adopting nuclear magnetic resonance MRI or CT; MRI also has long scanning time and inconvenient 3D modeling; CT also has the problems of large radiation dose and the like.
Description
Technical Field
The invention relates to the technical field of CT influence, in particular to a method for analyzing cervical vertebra in-vivo three-dimensional motion by combining cone beam CT and image fusion.
Background
The cervical vertebrae are divided into 7 segments C1-C7. The rotation, flexion and extension of the human head are all achieved by the compound movement of the segments. And the spatial position abnormality (such as straightening and curvature loss) and the movement abnormality of the segments are important factors for diagnosing, treating and judging the curative effect of clinically relevant diseases. X-ray is conventionally used for fluoroscopy, a two-dimensional static position is obtained, C1 is shielded by a skull, and C7 is shielded by a clavicle and cannot be imaged clearly. And the 3D static structure can be seen by adopting the magnetic resonance MRI or CT, but the dynamic motion abnormity can not be seen, and the two scans are usually in a supine position without load, and the cervical vertebra function of the human body can not be truly reflected. MRI also has long scanning time and inconvenient 3D modeling; CT also has the problems of large radiation dose and the like.
Disclosure of Invention
In order to solve the problem that X-ray is conventionally adopted for perspective shooting, a two-dimensional static position is obtained, C1 is shielded by a skull, and C7 is shielded by a clavicle and cannot be imaged clearly; the 3D static structure can be seen by adopting the magnetic resonance MRI or the CT, but the dynamic motion abnormity cannot be seen, and the two scans are usually in a supine position without load, so that the cervical vertebra function of the human body cannot be truly reflected; MRI also has long scanning time and inconvenient 3D modeling; the CT also has the problems of large radiation dose and the like, and the invention aims to provide a method for analyzing the in-vivo three-dimensional motion of the cervical vertebra by combining cone beam CT and image fusion.
In order to achieve the purpose, the invention adopts the following technical scheme: a method for analyzing the in-vivo three-dimensional motion of cervical vertebra by combining cone beam CT and image fusion comprises the following steps:
s1, shooting images of different cervical vertebra functional positions by a cone beam CT machine;
s2, importing images for segmentation and 3D modeling;
s3, determining a local coordinate system for the 3D cervical vertebra model in the neutral position according to the anatomical position of the human body;
s4, fusing images and recording the change of space position parameters of the cervical vertebra models of the segments with different functional positions relative to the neutral position;
s5, performing threshold judgment on the accuracy of image fusion to determine whether the image fusion is successful;
and S6, further calculating the change of the space position parameter between different sections to obtain the relative change of the two sections.
Preferably, the image fusion algorithm in S4 is specifically:
I. dispersing the neutral position and the functional position in the 3D model into point sets A and B distributed in a space;
II. Each point A finds a point B which is closest to the point A, and the distances are summed;
III, adjusting the position posture of the model B, repeating the step II, and performing optimization solution, wherein the optimization objective is to minimize the distance summation;
IV, finding the position posture of the optimized B model, and enabling the sum of the position posture and the point distance corresponding to the A to be minimum, namely the fusion degree to be highest;
v, recording the angle and displacement difference values in the original position and the optimized position space of the model B, namely the change of the space position parameters;
VI, recording the average distance between the optimized B model and the A, and judging whether the average distance is larger than a required distance threshold value Xmm.
Compared with the prior art, the invention has the following beneficial effects: the analysis method for the in-vivo three-dimensional motion of the cervical vertebra by combining the cone beam CT and the image fusion is high in speed, and 3D and perspective scanning can be completed within a few seconds; the method for analyzing the cervical vertebra in-vivo three-dimensional motion by combining the cone beam CT and the image fusion is objective and automatic in the whole flow, and has no step of artificial judgment; the analysis method of the cone beam CT and image fusion combined cervical vertebra in-vivo three-dimensional motion is low in radiation dose, and low-dose image equipment which is safe in clinic is adopted.
Drawings
The invention is described in further detail below with reference to the following figures and detailed description:
FIG. 1 is a schematic diagram of different functional positions of cervical vertebrae for cone-beam CT scanning according to the present invention;
FIG. 2 is a schematic diagram of 3D modeling of an image after fast scan according to the present invention;
FIG. 3 is a schematic view of the present invention for parametric registration of spatial positional movement by determining the orientation (local coordinates) of different segments of the cervical spine based on the anatomical surface of the body;
FIG. 4 is a schematic diagram of the present invention, in which different functional positions are spatially adjusted, then fusion matching is performed with a neutral position, and fusion accuracy is calculated, and when the fusion accuracy is lower than a certain threshold (e.g., 1mm), the fusion is successful;
fig. 5 is a schematic diagram showing the motion changes of adjacent segments obtained after the fusion according to the present invention, in which the lower segment of C7 is fused, and then the motion changes are sequentially generated and analyzed in the upper segment.
Detailed Description
The following description of the embodiments of the present invention is provided for illustrative purposes, and other advantages and effects of the present invention will become apparent to those skilled in the art from the present disclosure.
Please refer to fig. 1 to 5. It should be understood that the structures, ratios, sizes, and the like shown in the drawings and described in the specification are only used for matching with the disclosure of the specification, so as to be understood and read by those skilled in the art, and are not used to limit the conditions under which the present invention can be implemented, so that the present invention has no technical significance, and any structural modification, ratio relationship change, or size adjustment should still fall within the scope of the present invention without affecting the efficacy and the achievable purpose of the present invention. In addition, the terms "upper", "lower", "left", "right", "middle" and "one" used in the present specification are for clarity of description, and are not intended to limit the scope of the present invention, and the relative relationship between the terms and the terms is not to be construed as a scope of the present invention.
The invention provides a technical scheme that: a method for analyzing the in-vivo three-dimensional motion of cervical vertebra by combining cone beam CT and image fusion comprises the following steps:
s1, shooting images of different cervical vertebra functional positions by a cone beam CT machine;
s2, importing images for segmentation and 3D modeling;
s3, determining a local coordinate system for the 3D cervical vertebra model in the neutral position according to the anatomical position of the human body;
s4, fusing images and recording the change of space position parameters of the cervical vertebra models of the segments with different functional positions relative to the neutral position;
s5, performing threshold judgment on the accuracy of image fusion to determine whether the image fusion is successful;
and S6, further calculating the change of the space position parameter between different sections to obtain the relative change of the two sections.
The image fusion algorithm in S4 specifically includes:
I. dispersing the neutral position and the functional position in the 3D model into point sets A and B distributed in a space;
II. Each point A finds a point B which is closest to the point A, and the distances are summed;
III, adjusting the position posture of the model B, repeating the step II, and performing optimization solution, wherein the optimization objective is to minimize the distance summation;
IV, finding the position posture of the optimized B model, and enabling the sum of the position posture and the point distance corresponding to the A to be minimum, namely the fusion degree to be highest;
v, recording the angle and displacement difference values in the original position and the optimized position space of the model B, namely the change of the space position parameters;
VI, recording the average distance between the optimized B model and the A, and judging whether the average distance is larger than a required distance threshold value Xmm.
The foregoing embodiments are merely illustrative of the principles and utilities of the present invention and are not intended to limit the invention. Any person skilled in the art can modify or change the above-mentioned embodiments without departing from the spirit and scope of the present invention. Accordingly, it is intended that all equivalent modifications or changes which can be made by those skilled in the art without departing from the spirit and technical spirit of the present invention be covered by the claims of the present invention.
Claims (2)
1. A method for analyzing the in-vivo three-dimensional motion of cervical vertebra by combining cone beam CT and image fusion is characterized by comprising the following steps:
s1, shooting images of different cervical vertebra functional positions by a cone beam CT machine;
s2, importing images for segmentation and 3D modeling;
s3, determining a local coordinate system for the 3D cervical vertebra model in the neutral position according to the anatomical position of the human body;
s4, fusing images and recording the change of space position parameters of the cervical vertebra models of the segments with different functional positions relative to the neutral position;
s5, performing threshold judgment on the accuracy of image fusion to determine whether the image fusion is successful;
and S6, further calculating the change of the space position parameter between different sections to obtain the relative change of the two sections.
2. The method for analyzing the in-vivo three-dimensional motion of the cervical vertebrae by combining the cone beam CT and the image fusion as claimed in claim 1, wherein: the image fusion algorithm in S4 specifically includes:
I. dispersing the neutral position and the functional position in the 3D model into point sets A and B distributed in a space;
II. Each point A finds a point B which is closest to the point A, and the distances are summed;
III, adjusting the position posture of the model B, repeating the step II, and performing optimization solution, wherein the optimization objective is to minimize the distance summation;
IV, finding the position posture of the optimized B model, and enabling the sum of the position posture and the point distance corresponding to the A to be minimum, namely the fusion degree to be highest;
v, recording the angle and displacement difference values in the original position and the optimized position space of the model B, namely the change of the space position parameters;
VI, recording the average distance between the optimized B model and the A, and judging whether the average distance is larger than a required distance threshold value Xmm.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010853773.XA CN112085833B (en) | 2020-08-24 | 2020-08-24 | Analysis method for cone beam CT and image fusion combined in-vivo three-dimensional motion of cervical vertebra |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010853773.XA CN112085833B (en) | 2020-08-24 | 2020-08-24 | Analysis method for cone beam CT and image fusion combined in-vivo three-dimensional motion of cervical vertebra |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112085833A true CN112085833A (en) | 2020-12-15 |
| CN112085833B CN112085833B (en) | 2022-09-06 |
Family
ID=73729012
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010853773.XA Active CN112085833B (en) | 2020-08-24 | 2020-08-24 | Analysis method for cone beam CT and image fusion combined in-vivo three-dimensional motion of cervical vertebra |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112085833B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112991387A (en) * | 2021-03-10 | 2021-06-18 | 上海逸动医学科技有限公司 | Noninvasive three-dimensional joint motion analysis method |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1911175A (en) * | 2005-08-10 | 2007-02-14 | 史荣 | Image fusion method used for accurate determing of tumor shape |
| TW201114440A (en) * | 2009-10-26 | 2011-05-01 | Iner Aec Executive Yuan | A kind of radiolabelling method used as hepatic receptor imaging agent with multivalent glycoside clusters |
| TW201228673A (en) * | 2011-01-07 | 2012-07-16 | Iner Aec Executive Yuan | A radiolabeled nucleoside analogue, a method for preparing the same and the use thereof |
| CN102727329A (en) * | 2012-06-28 | 2012-10-17 | 深圳清华大学研究院 | Artificial cervical intervertebral disc |
| CN103720476A (en) * | 2013-12-13 | 2014-04-16 | 天津大学 | Stability evaluation method under curved path movement mode |
| CN107530094A (en) * | 2015-02-26 | 2018-01-02 | 钩椎关节有限责任公司 | Uncovertebral joint stabilizer and related system and method |
| CN107595387A (en) * | 2017-07-28 | 2018-01-19 | 浙江大学 | A kind of vertebra image generation system and spinal operation navigation positioning system based on ultrasonic rubbing technology |
| CN107607151A (en) * | 2017-07-17 | 2018-01-19 | 金季春 | One kind measurement weights and center of gravity of moving segments of human body parametric technique |
| CN109767410A (en) * | 2018-12-22 | 2019-05-17 | 复旦大学 | A kind of lung CT and MRI image blending algorithm |
| CN109919885A (en) * | 2019-02-15 | 2019-06-21 | 中国人民解放军陆军军医大学 | CVH image and Image registration fusion method based on B-spline and mutual information |
| CN109998543A (en) * | 2019-05-07 | 2019-07-12 | 黄光� | A kind of neck motor root conduction time measuring device and method |
| CN210582487U (en) * | 2018-10-15 | 2020-05-22 | 中国人民解放军第二军医大学第二附属医院 | Cervical vertebra three-dimensional mobility intelligent measuring instrument |
| CN111340800A (en) * | 2020-03-18 | 2020-06-26 | 联影智能医疗科技(北京)有限公司 | Image detection method, computer device, and storage medium |
-
2020
- 2020-08-24 CN CN202010853773.XA patent/CN112085833B/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1911175A (en) * | 2005-08-10 | 2007-02-14 | 史荣 | Image fusion method used for accurate determing of tumor shape |
| TW201114440A (en) * | 2009-10-26 | 2011-05-01 | Iner Aec Executive Yuan | A kind of radiolabelling method used as hepatic receptor imaging agent with multivalent glycoside clusters |
| TW201228673A (en) * | 2011-01-07 | 2012-07-16 | Iner Aec Executive Yuan | A radiolabeled nucleoside analogue, a method for preparing the same and the use thereof |
| CN102727329A (en) * | 2012-06-28 | 2012-10-17 | 深圳清华大学研究院 | Artificial cervical intervertebral disc |
| CN103720476A (en) * | 2013-12-13 | 2014-04-16 | 天津大学 | Stability evaluation method under curved path movement mode |
| CN107530094A (en) * | 2015-02-26 | 2018-01-02 | 钩椎关节有限责任公司 | Uncovertebral joint stabilizer and related system and method |
| CN107607151A (en) * | 2017-07-17 | 2018-01-19 | 金季春 | One kind measurement weights and center of gravity of moving segments of human body parametric technique |
| CN107595387A (en) * | 2017-07-28 | 2018-01-19 | 浙江大学 | A kind of vertebra image generation system and spinal operation navigation positioning system based on ultrasonic rubbing technology |
| CN210582487U (en) * | 2018-10-15 | 2020-05-22 | 中国人民解放军第二军医大学第二附属医院 | Cervical vertebra three-dimensional mobility intelligent measuring instrument |
| CN109767410A (en) * | 2018-12-22 | 2019-05-17 | 复旦大学 | A kind of lung CT and MRI image blending algorithm |
| CN109919885A (en) * | 2019-02-15 | 2019-06-21 | 中国人民解放军陆军军医大学 | CVH image and Image registration fusion method based on B-spline and mutual information |
| CN109998543A (en) * | 2019-05-07 | 2019-07-12 | 黄光� | A kind of neck motor root conduction time measuring device and method |
| CN111340800A (en) * | 2020-03-18 | 2020-06-26 | 联影智能医疗科技(北京)有限公司 | Image detection method, computer device, and storage medium |
Non-Patent Citations (2)
| Title |
|---|
| 任一 等: "苗医丰氏推拿手法对颈椎病患者颈椎节段空间位置影响的临床研究", 《基层医学论坛》 * |
| 曹盛楠 等: "神经根型颈椎病C3～C7三维有限元模型的建立与意义", 《山东医药》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112991387A (en) * | 2021-03-10 | 2021-06-18 | 上海逸动医学科技有限公司 | Noninvasive three-dimensional joint motion analysis method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112085833B (en) | 2022-09-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6280676B2 (en) | Spine arrangement estimation device, spine arrangement estimation method, and spine arrangement estimation program | |
| Penney et al. | Validation of a two‐to three‐dimensional registration algorithm for aligning preoperative CT images and intraoperative fluoroscopy images | |
| US8731268B2 (en) | CT device and method based on motion compensation | |
| US8891847B2 (en) | Automatic implant detection from image artifacts | |
| US7844094B2 (en) | Systems and methods for determining geometric parameters of imaging devices | |
| WO2019020048A1 (en) | Spinal image generation system based on ultrasonic rubbing technique and navigation positioning system for spinal surgery | |
| US20080037843A1 (en) | Image segmentation for DRR generation and image registration | |
| US8121380B2 (en) | Computerized imaging method for a three-dimensional reconstruction from two-dimensional radiological images; implementation device | |
| Baka et al. | Statistical shape model-based femur kinematics from biplane fluoroscopy | |
| Munbodh et al. | Automated 2D‐3D registration of a radiograph and a cone beam CT using line‐segment enhancement a | |
| Novosad et al. | Three-dimensional (3-D) reconstruction of the spine from a single X-ray image and prior vertebra models | |
| JP6316078B2 (en) | Patient positioning system and method of calculating displacement | |
| Vo et al. | A semi-automatic 3D ultrasound reconstruction method to assess the true severity of adolescent idiopathic scoliosis | |
| Lin et al. | Intervertebral anticollision constraints improve out‐of‐plane translation accuracy of a single‐plane fluoroscopy‐to‐CT registration method for measuring spinal motion | |
| CN112085833B (en) | Analysis method for cone beam CT and image fusion combined in-vivo three-dimensional motion of cervical vertebra | |
| Haque et al. | Hierarchical model-based tracking of cervical vertebrae from dynamic biplane radiographs | |
| Eldib et al. | A motion artifact reduction method for dental CT based on subpixel-resolution image registration of projection data | |
| Whitmarsh et al. | 3D bone mineral density distribution and shape reconstruction of the proximal femur from a single simulated DXA image: an in vitro study | |
| Bögel et al. | Respiratory motion compensation using diaphragm tracking for cone-beam C-arm CT: A simulation and a phantom study | |
| Chang et al. | A Photogrammetric system for 3D reconstruction of a scoliotic torso | |
| Banks | Model based 3D kinematic estimation from 2D perspective silhouettes: application with total knee prostheses | |
| Koo et al. | In vivo measurement of the human lumbar spine using magnetic resonance imaging to ultrasound registration | |
| Fujinuma et al. | Quantitative evaluation method for tibiofibular joint alignment in ankle osteoarthritis based on 3D bone size | |
| Cukovic et al. | Development of 3D kinematic model of the spine for idiopathic scoliosis simulation | |
| Miyabe et al. | New algorithm to simulate organ movement and deformation for four‐dimensional dose calculation based on a three‐dimensional CT and fluoroscopy of the thorax |
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 |