CN109994198A - A kind of method for reconstructing of dynamic heart threedimensional model - Google Patents
A kind of method for reconstructing of dynamic heart threedimensional model Download PDFInfo
- Publication number
- CN109994198A CN109994198A CN201910267439.3A CN201910267439A CN109994198A CN 109994198 A CN109994198 A CN 109994198A CN 201910267439 A CN201910267439 A CN 201910267439A CN 109994198 A CN109994198 A CN 109994198A
- Authority
- CN
- China
- Prior art keywords
- model
- heart
- dynamic
- cardiac
- reconstructing
- 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
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
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
Abstract
The present invention provides a kind of method for reconstructing of dynamic heart threedimensional model, comprising: the cardiac module of important phase in selected cardiac cycle；Time parameter is added in the mathematical model expression of step 1, indicates the cardiac module in space-time；The motion profile of each point on the cardiac module that step 2 indicates is fitted with interpolation；It checks that the model for not meeting convention in Interpolation Process interts part, and is modified for unreasonable part；As unit of a cardiac cycle, Dynamically Announce heart movement is to get dynamic heart threedimensional model.The present invention can establish the dynamic 3 D cardiac module with medical value, reproduce physiological movement of the heart within cardiac cycle, model includes atrium sinistrum, atrium dextrum, left ventricle, right ventricle, left ventricle outer wall and all cardia of entire chambers of the heart outer wall, it can be realized comprehensive analysis of heart movement feature based on this, and help to check for heart physiological dysfunction.
Description
Technical field
The invention belongs to field of Computer Graphics, are related to three-dimensional reconstruction, and in particular to a kind of three-dimensional cardiac moves
State reconstruction technique.
Background technique
Heart is one of most important organ and a structure and all extremely complex organ of movement in human body, wherein
The diastole and contraction of atrium and ventricle are the power resources of the circulatory system.Cardiovascular disease is world today's morbidity and mortality
One of highest disease can influence labour and the daily life of patient significantly.But existing medical technology can not be moved
State three-dimensionally shows the motion state of patient's heart, it is difficult to carry out accurate description to it.Dynamic three-dimensional reconstruction with image segmentation,
The technological means such as image registration and information fusion, three-dimensional visualization, motion analysis, reconstruct the Three-Dimensional Dynamic model of heart, energy
The heart movement information collection for enough realizing higher-dimension, has very important significance, establishes heart dynamic three-dimensional reconstruction theory gesture and exist
It must go.
So far, there are also research achievements for the Medical Image Processing of heart, are schemed using magnetic resonance image, CT mostly
Picture or ultrasound image model heart, wherein being that this individual chamber models to left ventricle mostly, for heart
The modeling of other three chambers or even entire heart is less.In addition, existing medical visualization technology is mostly static three-dimensional model,
It is only limited to the 3 D anatomical mechanism of description heart, the foundation of rare dynamic 3 D model, the motion information of tissue is not
It preserves, therefore is difficult to accurately and effectively evaluate dynamic heart function.
Therefore, the dynamic 3 D model for establishing heart, can show the motion information of heart, medically have important valence
Value.And still lack the modeling method easily realized in the art at present.
Summary of the invention
To solve the above problems, can be utilized quiet the present invention provides a kind of method for reconstructing of dynamic heart threedimensional model
The cardiac three-dimensional model of state generates the dynamic cardiac three-dimensional model with medical value, reproduces heart within cardiac cycle
Physiological movement.
In order to achieve the above object, the invention provides the following technical scheme:
A kind of method for reconstructing of dynamic heart threedimensional model, includes the following steps:
Step 1 selectes the cardiac module of important phase in cardiac cycle；
Step 2 adds time parameter in the mathematical model expression of step 1, indicates the cardiac module in space-time；
Step 3 fits the motion profile of each point on the cardiac module that step 2 indicates with interpolation；
Step 4 checks that the model for not meeting convention in Interpolation Process interts part, and is repaired for unreasonable part
Change；
Step 5, as unit of a cardiac cycle, Dynamically Announce heart movement is to get dynamic heart threedimensional model.
Further, the step 1 cardiac model is static cardiac module；Time parameter is added in the step 2
The point on model is expressed as the form of (x, y, z, t) afterwards, wherein x, and y, z are respectively the point coordinate on static cardiac module, when t is
Between parameter；Interpolation is carried out respectively to three coordinates of all the points on model with cubic spline interpolation method in the step 3
Fit the motion profile of each point on cardiac module.
Further, the step 1 cardiac model is to generate after carrying out principal component analysis by static cardiac three-dimensional model
Mathematical model；(x is expressed as plus model after time parameter in the step 21,x2,…xn, t) form, wherein x1,…,
xnIt is all coefficients of the phase heart mathematical model, t is time parameter；All coefficients of model are distinguished in the step 3
Interpolation is carried out to fit the motion profile of each point on cardiac module.
Further, the step 1 heart cycles include ventricular systole and ventricular diastole.
Further, important phase includes at least ventricular diastole latter stage and/or ventricular end systolic in the step 1.
Further, the step 4 specifically comprises the following steps: that the movement to heart entirety is tested, if there is model
The case where surface is interspersed or separates is modified, and is tested to time shaft, modifying if having unreasonable situation it makes it
More meet the actual motion track of heart.
Further, the step 5 further includes following process:
Social assessment is made to dynamic heart threedimensional model, using the criterion as reconstruction model superiority and inferiority.
Compared with prior art, the invention has the advantages that and the utility model has the advantages that
1. the present invention can establish the dynamic 3 D cardiac module with medical value, including atrium sinistrum, atrium dextrum, Zuo Xin
Room, right ventricle, left ventricle outer wall and all cardia of entire chambers of the heart outer wall, can be realized heart movement feature based on this
Analysis comprehensively.2. the heart dynamic model ratio obtained by the mathematical model interpolation after principal component analysis is direct by static models
The dynamic model that interpolation obtains is more smooth.
3. the Normal Human Heart's dynamic model and sufferer heart established compare, can help to check for heart physiological
Dysfunction.4. judge by superiority and inferiority of the social assessment to model makes it possess the specialty evaluation under Medicine standard.
Detailed description of the invention
Fig. 1 is the flow chart for the dynamic heart method for reconstructing three-dimensional model that the embodiment of the present invention one provides.
Fig. 2 is the static cardiac three-dimensional model schematic referred in the embodiment of the present invention one.
Fig. 3 is the flow chart of dynamic heart method for reconstructing three-dimensional model provided by Embodiment 2 of the present invention.
Fig. 4 is the static cardiac three-dimensional mathematical model schematic diagram referred in the embodiment of the present invention two.
Fig. 5 be the embodiment of the present invention one, in embodiment two interpolation method schematic diagram.
Specific embodiment
Technical solution provided by the invention is described in detail below with reference to specific embodiment, it should be understood that following specific
Embodiment is only illustrative of the invention and is not intended to limit the scope of the invention.
Embodiment one:
The present embodiment is that model interpolation is directly carried out to existing multidate static state cardiac module, complete it is whole to heart with
And the estimation of each chamber movement track, dynamic heart threedimensional model is finally obtained, specifically, the dynamic heart that this example provides
The method for reconstructing process of threedimensional model is as shown in Figure 1, comprising the following steps:
Step S110, select cardiac cycle in important phase static cardiac module, static cardiac module as shown in Fig. 2,
Including ventricular systole (isovolumic contraction period 0.05s, phase of maximum ejection 0.09s, reduced ejection period 0.13s) and ventricular diastole
(protodiastole 0.04s, isovolumic relaxation period 0.08s, phase of rapid filling 0.11s, reduced filling period 0.19s and auricular systole period
0.1s), important phase therein is selected, other phases such as diastasis and end-systole and during cardiac cycle.Wherein
The point of these models is mutual corresponding, it is known that the three-dimensional coordinate of the identical point of each phase.
Step S120, the point in static three-dimensional model add a time parameter t, indicate the locating moment, that is, use (x, y, z, t)
The point coordinate on the cardiac module in space-time is indicated, as the basis of dynamic modeling, for characterizing dynamic heart three-dimensional
Model.
Step S130 fits the motion profile of each point on cardiac module with interpolation.In interpolation method, cubic spline
Interpolation can guarantee that functional value, first derivative, second dervative are all continuous, so that model interpolation result is as smooth mellow and full hairless as possible
Thorn.Therefore cubic spline interpolation method is chosen to fit the motion profile of each point on cardiac module, guarantee that fitting track must wrap
Containing the location of known models, and meet the kinematics character and shared time scale of each phase of heart as far as possible.Interpolation
Method is as shown in Figure 5.Interpolation is carried out respectively to the three-dimensional coordinate of each point.By taking x coordinate as an example, interpolation problem can be equivalent to
Know that (wherein independent variable x is model at the moment to n+1 node of the function y=f (x) on section [a, b], and dependent variable y is three-dimensional
Coordinate), a=x0<x1<…<xn=b, yi=f (xi) (i=0,1 ..., n), interpolating function S (x) is sought, so that meeting:
Wherein, there is form
S (x)=Ci(x), xi-1≤x≤xi
Ci(x)=ai+bix+Cix2+dix3
It is intended to solve, 4n coefficient need to be calculated.The wherein available 2n equation of the interface point in each section:
And by Second Order Continuous can micro- condition 2n-2 equation can be obtained:
C″i(xi)=C "i+1(xi)→2ci+6dixi=2ci+1+6di+1xi
Again since heart is cyclically-varying, therefore assume that S (x) has periodically, meets at endpoint
S ' (a+0)=S ' (b-0), S " (a+0)=S " (b-0)
Total 4n equation, can pass through Solving Linear.Therefore interpolation model is calculated.
Step S140 tests to the movement of heart entirety, carries out if the case where having model surface interspersed or separation
Modification.It tests to time shaft, there is unreasonable part then to modify, it is made more to meet the actual motion track of heart.
Step S150, as unit of a cardiac cycle, Dynamically Announce heart movement process, and it is made with Medicine standard
Social assessment out, the criterion as dynamic model quality good or not.
Embodiment two:
As the improvement of embodiment one, the present embodiment is to carry out motion profile by the mathematical model of the PCA heart formed
Reproduction, to mathematical model carry out interpolation, to realize the foundation of four-dimensional cardiac model.Specifically, the dynamic that this example provides
The method for reconstructing process of cardiac three-dimensional model is as shown in Figure 3, comprising the following steps:
Step S210 selectes the heart mathematical model of important phase in cardiac cycle, as shown in figure 4, including ventricular contraction
Phase (isovolumic contraction period 0.05s, phase of maximum ejection 0.09s, reduced ejection period 0.13s) and ventricular diastole (protodiastole
0.04s, isovolumic relaxation period 0.08s, phase of rapid filling 0.11s, reduced filling period 0.19s and auricular systole period 0.1s).Wherein this
A little mathematical models are to carry out PCA (principal component analysis) by static cardiac three-dimensional model to generate, and the coefficient of phase and each single item is all
It is known.
Model formulation is (x by step S2201,x2,…xn, t) form, wherein x1,…,xnIt is the phase heart mathematics
All coefficients of model add time parameter t in original mathematical model expression, to characterize the cardiac module under the four-dimension.
Step S230, it is similar with step S130, only S130 be three coordinates of all the points on model are carried out respectively it is slotting
Value, as long as carrying out interpolation respectively to all coefficients of model here.Meter is greatly reduced in the case where guaranteeing model accuracy
Calculation amount, and mathematical model interpolation is more more smooth than the effect of model interpolation and the kinetic characteristic that gets close to nature.
Step S240 tests to the movement of heart entirety, if having model surface interspersed or being repaired the case where separation
Change.It tests to time shaft, it is made more to meet the actual motion track of heart.
Step S250, as unit of a cardiac cycle, Dynamically Announce heart movement process, and it is made with Medicine standard
Social assessment out, the criterion as dynamic model quality good or not.
It is demonstrated experimentally that the heart dynamic model ratio that this example is obtained based on the mathematical model interpolation after principal component analysis pass through it is quiet
The dynamic model that the direct interpolation of states model obtains is more smooth.
By the situation comparison under heart dynamic process and pathologic condition under normal circumstances, heart physiological can be intuitively found out
The unusual part of function.This method can provide more information for the diagnosis of heart disease and the evaluation of cardiac function.
Those skilled in the art should be understood that each module of the above invention or each step can use general calculating
Device realizes that they can be concentrated on a single computing device, or be distributed in network constituted by multiple computing devices
On, optionally, they can be realized with the program code that computing device can perform, it is thus possible to be stored in storage
It is performed by computing device in device, perhaps they are fabricated to each integrated circuit modules or will be more in them
A module or step are fabricated to single integrated circuit module to realize.In this way, the present invention is not limited to any specific hardware and
Software combines.
The technical means disclosed in the embodiments of the present invention is not limited only to technological means disclosed in above embodiment, further includes
Technical solution consisting of any combination of the above technical features.It should be pointed out that for those skilled in the art
For, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also considered as
Protection scope of the present invention.
Claims (7)
1. a kind of method for reconstructing of dynamic heart threedimensional model, which comprises the steps of:
Step 1 selectes the cardiac module of important phase in cardiac cycle；
Step 2 adds time parameter in the mathematical model expression of step 1, indicates the cardiac module in space-time；
Step 3 fits the motion profile of each point on the cardiac module that step 2 indicates with interpolation；
Step 4 checks that the model for not meeting convention in Interpolation Process interts part, and is modified for unreasonable part；
Step 5, as unit of a cardiac cycle, Dynamically Announce heart movement is to get dynamic heart threedimensional model.
2. the method for reconstructing of dynamic heart threedimensional model according to claim 1, it is characterised in that: the step 1 center
Dirty model is static cardiac module；The shape of (x, y, z, t) is expressed as in the step 2 plus the point on model after time parameter
Formula, wherein x, y, z are respectively the point coordinate on static cardiac module, and t is time parameter；Spline interpolation is used in the step 3
Method carries out interpolation to three coordinates of all the points on model to fit the motion profile of each point on cardiac module respectively.
3. the method for reconstructing of dynamic heart threedimensional model according to claim 1, it is characterised in that: the step 1 center
Dirty model is the mathematical model for generate after principal component analysis by static cardiac three-dimensional model；The time is added in the step 2
Model is expressed as (x after parameter1,x2,…xn, t) form, wherein x1,…,xnIt is all systems of the phase heart mathematical model
Number, t is time parameter；It is each on cardiac module to fit to carry out interpolation respectively to all coefficients of model in the step 3
The motion profile of point.
4. the method for reconstructing of dynamic heart threedimensional model described in any one of -3 according to claim 1, it is characterised in that: institute
Stating step 1 heart cycles includes ventricular systole and ventricular diastole.
5. the method for reconstructing of dynamic heart threedimensional model according to claim 4, it is characterised in that: weight in the step 1
Phase is wanted to include at least ventricular diastole latter stage and/or ventricular end systolic.
6. the method for reconstructing of dynamic heart threedimensional model described in any one of -3 according to claim 1, it is characterised in that: institute
It states step 4 and specifically comprises the following steps: that the movement to heart entirety is tested, if having model surface interspersed or isolated feelings
Condition is modified, and is tested to time shaft, and practical fortune for making it more meet heart is modified it if having unreasonable situation
Dynamic rail mark.
7. the method for reconstructing of dynamic heart threedimensional model described in any one of -3 according to claim 1, it is characterised in that: institute
Stating step 5 further includes following process:
Social assessment is made to dynamic heart threedimensional model, using the criterion as reconstruction model superiority and inferiority.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910267439.3A CN109994198A (en) | 2019-04-03 | 2019-04-03 | A kind of method for reconstructing of dynamic heart threedimensional model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910267439.3A CN109994198A (en) | 2019-04-03 | 2019-04-03 | A kind of method for reconstructing of dynamic heart threedimensional model |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109994198A true CN109994198A (en) | 2019-07-09 |
Family
ID=67132295
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910267439.3A Pending CN109994198A (en) | 2019-04-03 | 2019-04-03 | A kind of method for reconstructing of dynamic heart threedimensional model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109994198A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111816313A (en) * | 2020-07-01 | 2020-10-23 | 杭州脉流科技有限公司 | Simulation display method and device for analyzing myocardial motion based on left ventricle three-dimensional model, computer equipment and storage medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101639948A (en) * | 2009-08-20 | 2010-02-03 | 浙江工业大学 | Method for calculating characteristic value and characteristic vector of flexible body point distribution model based on interpolation algorithm |
| CN101901498A (en) * | 2009-12-31 | 2010-12-01 | 华中科技大学 | Recursion modelling method of human coronary artery |
| US20140173486A1 (en) * | 2012-05-14 | 2014-06-19 | Heartflow, Inc. | Method and system for providing information from a patient-specific model of blood flow |
| CN106600596A (en) * | 2016-12-21 | 2017-04-26 | 南昌航空大学 | Dynamic display method and system for heart cavity |
-
2019
- 2019-04-03 CN CN201910267439.3A patent/CN109994198A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101639948A (en) * | 2009-08-20 | 2010-02-03 | 浙江工业大学 | Method for calculating characteristic value and characteristic vector of flexible body point distribution model based on interpolation algorithm |
| CN101901498A (en) * | 2009-12-31 | 2010-12-01 | 华中科技大学 | Recursion modelling method of human coronary artery |
| US20140173486A1 (en) * | 2012-05-14 | 2014-06-19 | Heartflow, Inc. | Method and system for providing information from a patient-specific model of blood flow |
| CN106600596A (en) * | 2016-12-21 | 2017-04-26 | 南昌航空大学 | Dynamic display method and system for heart cavity |
Non-Patent Citations (1)
| Title |
|---|
| WOCHUBU: "插值法-三次样条插值", 《HTTPS://WWW.DOCIN.COM/P-1636017667.HTML》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111816313A (en) * | 2020-07-01 | 2020-10-23 | 杭州脉流科技有限公司 | Simulation display method and device for analyzing myocardial motion based on left ventricle three-dimensional model, computer equipment and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10299862B2 (en) | Three-dimensional quantitative heart hemodynamics in medical imaging | |
| Bermejo et al. | Intraventricular vortex properties in nonischemic dilated cardiomyopathy | |
| CN108992057A (en) | A kind of method and apparatus of determining coronary flow reserve score FFR | |
| Banerjee et al. | A completely automated pipeline for 3D reconstruction of human heart from 2D cine magnetic resonance slices | |
| US9462952B2 (en) | System and method for estimating artery compliance and resistance from 4D cardiac images and pressure measurements | |
| KR102620470B1 (en) | Patient-specific modeling of coronary hemodynamic parameters | |
| Wang et al. | A knowledge‐based cone‐beam x‐ray CT algorithm for dynamic volumetric cardiac imaging | |
| CN109961508A (en) | A kind of heart static three-dimensional model method for reconstructing | |
| CN104825133B (en) | The quasistatic ventricular heart magnetic field model being imaged based on color Doppler 3D | |
| Nguyen et al. | A semi-automated method for patient-specific computational flow modelling of left ventricles | |
| CN101116104A (en) | A method, a system and a computer program for segmenting a surface in a multidimensional dataset | |
| CN109994198A (en) | A kind of method for reconstructing of dynamic heart threedimensional model | |
| Wong et al. | Velocity-based cardiac contractility personalization from images using derivative-free optimization | |
| JP2013162921A (en) | Image processing apparatus, image processing method, and image processing program | |
| CN110457765B (en) | Method for calculating odd vein shunt ratio in half-ventricle model based on geometric multi-scale model | |
| Dahl | Numerical simulations of blood flow in the left side of the heart | |
| Lee et al. | Validation of three‐dimensional echocardiographic principal strain analysis for assessing left ventricular contractility: An animal study | |
| Behdadfar et al. | A centerline-based model morphing algorithm for patient-specific finite element modeling of the left ventricle | |
| Kounchev et al. | Application of PDE methods to visualization of heart data | |
| Wan et al. | Reconstructing Cardiac Shape via Constrained Voronoi Diagram and Cyclic Dynamic Time Wrapping From CMR | |
| Choi et al. | Motion visualization of human left ventricle with a time‐varying deformable model for cardiac diagnosis | |
| Patel | 3D Reconstruction of Left Ventricle with Mitral and Aortic Valves from Multiple Axis 2D Echocardiography at Standard Cross Sections and Flow Simulations | |
| Yang et al. | Modeling of individual left ventricle for cardiac resynchronization therapy | |
| Segars et al. | Enhanced 4D heart model based on high resolution dual source gated cardiac CT images | |
| Serinagaoglu et al. | A Bayesian approach to inclusion and performance analysis of using extra information in bioelectric inverse problems |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190709 |