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 PDF

Info

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
heart
dynamic
cardiac
reconstructing
threedimensional
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
Application number
CN201910267439.3A
Other languages
Chinese (zh)
Inventor
王雁刚
朱婷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201910267439.3A priority Critical patent/CN109994198A/en
Publication of CN109994198A publication Critical patent/CN109994198A/en
Pending legal-status Critical Current

Links

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

A kind of method for reconstructing of dynamic heart threedimensional model
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.
CN201910267439.3A 2019-04-03 2019-04-03 A kind of method for reconstructing of dynamic heart threedimensional model Pending CN109994198A (en)

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)

Citations (4)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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

Similar Documents

Publication Publication Date Title
US10871536B2 (en) Automated cardiac volume segmentation
Singh et al. LA strain for categorization of LV diastolic dysfunction
US10483005B2 (en) System and method for characterization of electrical properties of the heart from medical images and body surface potentials
US9805168B2 (en) Method and system for sensitivity analysis in modeling blood flow characteristics
Mittal et al. Computational modeling of cardiac hemodynamics: current status and future outlook
CN105518684B (en) For predicting device, system and the non-transitory computer-readable medium of coronary artery pathological changes relevant information
Lamata et al. An accurate, fast and robust method to generate patient-specific cubic Hermite meshes
Nelson et al. Visualization of 3D ultrasound data
CN102525443B (en) For the method and apparatus to cardiovascular circulation modeling based on medical image
Chen et al. Modeling, analysis, and visualization of left ventricle shape and motion by hierarchical decomposition
US10311978B2 (en) Method and system for patient specific planning of cardiac therapies on preoperative clinical data and medical images
Aguado-Sierra et al. Patient-specific modeling of dyssynchronous heart failure: a case study
Seed et al. Feasibility of quantification of the distribution of blood flow in the normal human fetal circulation using CMR: a cross-sectional study
Cheng et al. Fluid-structure coupled CFD simulation of the left ventricular flow during filling phase
Krishnamurthy et al. Patient-specific models of cardiac biomechanics
US10354744B2 (en) Non-invasive functional assessment of coronary artery stenosis including simulation of hyperemia by changing resting microvascular resistance
Hunter et al. Modeling total heart function
US9129053B2 (en) Method and system for advanced measurements computation and therapy planning from medical data and images using a multi-physics fluid-solid heart model
Carlsson et al. Total heart volume variation throughout the cardiac cycle in humans
Papademetris et al. Estimation of 3-D left ventricular deformation from medical images using biomechanical models
Huang et al. Spatio-temporal tracking of myocardial deformations with a 4-D B-spline model from tagged MRI
US8527251B2 (en) Method and system for multi-component heart and aorta modeling for decision support in cardiac disease
KR20140091741A (en) Method and system for multi-scale anatomical and functional modeling of coronary circulation
Thavendiranathan et al. Feasibility, accuracy, and reproducibility of real-time full-volume 3D transthoracic echocardiography to measure LV volumes and systolic function: a fully automated endocardial contouring algorithm in sinus rhythm and atrial fibrillation
KR101855369B1 (en) Systems and methods for estimating physiological heart measurements from medical images and clinical data

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