WO2005071615A1 - Stochastic analysis of cardiac function - Google Patents
Stochastic analysis of cardiac function Download PDFInfo
- Publication number
- WO2005071615A1 WO2005071615A1 PCT/IB2005/050060 IB2005050060W WO2005071615A1 WO 2005071615 A1 WO2005071615 A1 WO 2005071615A1 IB 2005050060 W IB2005050060 W IB 2005050060W WO 2005071615 A1 WO2005071615 A1 WO 2005071615A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- images
- set forth
- samples
- functional parameter
- volume
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
Definitions
- a database representing temporal image slices of the heart is first generated.
- the database may be described as a temporal sequence of slices of the heart, or a stack of slices at temporal intervals, covering at least one whole cycle of the heart.
- An attempt is then made to delineate the cardiac contours in each heart phase from the stack of slices.
- Current methods of accomplishing this task are manual or semi-manual using visualization software. For example, a clinician marks boundary points of a heart chamber in each slice of the stack of slices at each temporal interval.
- the shape (volume) of the heart chamber is determined in each of the temporal intervals, at least for the temporal intervals corresponding to the cardiac phases of interest. From the estimated shapes, functional parameters such as end-diastolic and end-systolic left ventricular volumes, stroke volumes and ejection fraction are computed. With a typical sample of 200 slices, manually delineating the cardiac chambers is time-consuming and error prone. Methods of accomplishing the above-described procedure automatically are currently being tried in the research art. One automated procedure is to fit to the image slices a mathematical model of the heart. This model is based on some prior knowledge on the shape of the myocardium, and possibly a statistical description of the typical shape variations from patient to patient.
- the automatic segmentation technique generally consists in optimizing an energy functional made of two terms: the first one describes the fit of the model to the data, while the second one penalizes too strong deformations with respect to some reference configuration.
- the result of this procedure is one optimal segmentation with respect to the mathematical model used to represent the myocardium shape and the mathematical form of the energy used for the optimization.
- MRI scans are noisy, contain artifacts and only a few slices are acquired, with possible gaps between the slices.
- the image information does not permit delineating the cardiac chambers without ambiguity. Due to these problems and the like, the segmentation problem is characterized by an inherent uncertainty. This aspect is well known in clinical practice, where it is commonly accepted that the marked boundary points vary from clinician to clinician in manual techniques.
- a method of stochastic model-based segmentation of diagnostic images of a subject, from samples of the subject is provided.
- the samples are generated according to a Bayesian stochastic model describing the conditional probability distribution of the organ shape given the images and functional parameters are derived for each of the respective samples. An uncertainty value is estimated for each derived parameter .
- a diagnostic imaging apparatus is provided.
- the apparatus includes a means for generating the diagnostic images and a processor programmed to perform the above-described method.
- an apparatus for generating stochastic model-based segmentation of diagnostic images of a subject is provided.
- the apparatus includes a diagnostic imaging scanner configured to perform a scan of a volume region of interest of a subject.
- a processing system is provided to process scanned data of the volume region of interest from the diagnostic imaging scanner and a provided reconstruction module receives the scanning data and generate a stack of slice images of the region of interest stored in the processing system.
- a shape module computes shape samples representing multiple solutions of the segmentation of the slice images according to a probability distribution described by a Bayesian model, and a function module derives functional parameters for each of the respective samples and estimates a probability value for each derived parameter.
- a display module displays the functional parameters and respective probability values on a provided display device.
- a diagnostic imaging apparatus performs stochastic model-based segmentation of diagnostic images of a subject.
- a computing means computes a plurality of shape samples according to a probability distribution described by a Bayesian model.
- a function means derives at least one functional parameter for each of the respective samples and estimates a probability value for each derived parameter.
- One advantage resides in an improved robustness of the segmentation and the tracking of the heart chambers.
- Another advantage resides in the precise information provided to the clinician regarding the accuracy of the functional parameters generated by the automated process.
- Yet another advantage resides in the means for an improved diagnosis of cardiac function, with an improved level of confidence.
- FIGURE 1 diagrammatically shows a magnetic resonance imaging system employing automated estimating of cardiac functional parameters according to concepts of the present invention.
- a diagnostic imaging scanner 10 such as a magnetic resonance scanner, includes a housing 12 defining a generally cylindrical scanner bore 14 inside of which an associated imaging subject 16 is disposed. Details of the diagnostic imaging scanner are not shown because magnetic resonance, CT, SPECT, PET, and other suitable scanners are well known in the art. It suffices to say that the diagnostic imaging scanner 10 performs a cardiac scan and communicates scan data to a reconstruction module 18 resident in a processing system 20.
- the processing system 20 presented herein is not inherently related to any particular computer or other apparatus. In particular, various general-purpose machines may be used with program modules in accordance with the teachings herein, or it may prove more convenient to construct more specialized apparatus to perform the required method steps.
- the processing system 20 may be a single system or an interconnected distributed system of processors.
- the reconstruction module 18 generates a series of volume images at temporal intervals during the cardiac cycle.
- Each of the temporally offset volume images includes a stack of slice images 22. Due to possible poor contrast between the cardiac tissues and the surroundings of the myocardium, image noise, discrete pixel sizes, and the like, the boundary of the cardiac chamber is not without ambiguity.
- a shape module or means 24 computes multiple solutions of the heart shape rather than the single optimized solution of prior-art methods.
- the Bayesian probability distribution can be represented by a finite number of samples, for example 500.
- the multiple samples are stored in a myocardium shape sample memory 26.
- These multiple solutions can be displayed to the clinician for example in the form of an animated shape, possibly superposed on the image slices, so that a qualitative representation of the set of solutions is given to the clinician.
- a function module 28 then accesses the shape samples 26 and computes cardiac functional parameters and statistical parameters 30 from the multiple shapes 26, for example ejection fraction, end-diastolic volume, end-systolic volume, stroke volume, wall thickness, and the like.
- the exemplary functional parameter is the ejection fraction 32.
- the present invention computes a plurality of ejection fractions from the plurality of shapes in the shape memory 26 and statistical parameters from this plurality of the ejection fractions.
- a histogram is computed, having a relative frequency or probability 34 for each of the computed ejection fractions 32.
- the functional parameters and statistical parameters 30 are then presented via a display module 36 on a display device 38 for analysis by a clinician or other user of the diagnostic imaging system 10.
- the display can be in graphical form as a histogram, numerical form as a median value and standard deviation, or the like. In this manner, a clinician is presented along with the solutions, a measure of the degree of confidence in the solutions.
- an animation module 39 displays the plurality of solutions in the aforementioned form of an animated shape, possibly superposed on the images, so that a qualitative representation of the set of solutions is given to the clinician.
- the animation sequence may be ordered according to a predetermined criteria such as, for example, ordered from low probability to high probability, most far from a reference shape to the most close to the reference shape, etc.
- the animation may also be superposed on the most probable shape.
- One approach that the shape module 24 uses in computing the shape samples 26 involves the use of Markov chain algorithms as described by W.R. Gilks, S. Richardson and D.J. Spiegelhalter in Markov Chain Monte Carlo in Practice, Chapman and Hall, 1966. To deal with the temporal component of the problem, sequential approaches as described by M. Isard and A. Blake in Condensation — conditional density propagation for visual tracking, International Journal of Computer Vision, 1998.
- the general approach is to construct a Bayesian probability distribution describing the space of solutions for the segmentation problem, and to generate samples of this probability distribution.
- the prior model describes the prior knowledge about the shape; it is a probability distribution on the parameters used to describe the shape (for example the coordinates of the nodes of a mesh representing the organ shape, or the coefficients of the decomposition of the surface based on surface harmonics).
- the likelihood model describes the structure of the images (for example the statistical distribution of the grey values) for a fixed shape.
- Monte Carlo integration allows one to compute this integral straightforwardly by generating a sample ..., z n ) of the distribution ⁇ (z/y) and to use the approximation:
- Computing statistics on the functional parameters requires, therefore, generating a finite number of samples of the Bayesian probability distribution ⁇ (z/y). This can be done, for example, by generating a Markov chain which has the Bayesian probability distribution as a stationary distribution, see Markov Chain Monte Carlo in Practice for details.
- the samples span the range of possible solutions. If a temporal sequence of slices is to be segmented, the segmentation consists then in a temporal series of shapes. Such a sample can be sequentially generated, using sequential Monte Carlo techniques.
- a finite number of shapes representing volumes are generated.
- a motion model is used to predict the next set of shapes.
- a simple motion model assumes that the myocardium shape contracts by a constant factor at each time step.
- a weight of each sample which depends on the likelihood value of the sampled shape for the new time step, is then computed and used in the Monte Carlo integration.
- the sequential Monte Carlo method is well known in the arts, and used in other applications such as, for example, financial mathematics, tracking, etc.
- the method also called particle filtering, is useful for following in time objects in motion.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
- Image Generation (AREA)
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Description
Claims
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006548516A JP2007521862A (en) | 2004-01-15 | 2005-01-05 | Stochastic analysis of cardiac function |
EP05702590A EP1709590A1 (en) | 2004-01-15 | 2005-01-05 | Stochastic analysis of cardiac function |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US53676204P | 2004-01-15 | 2004-01-15 | |
US60/536,762 | 2004-01-15 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2005071615A1 true WO2005071615A1 (en) | 2005-08-04 |
Family
ID=34807047
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/IB2005/050060 WO2005071615A1 (en) | 2004-01-15 | 2005-01-05 | Stochastic analysis of cardiac function |
Country Status (4)
Country | Link |
---|---|
EP (1) | EP1709590A1 (en) |
JP (1) | JP2007521862A (en) |
CN (1) | CN1910618A (en) |
WO (1) | WO2005071615A1 (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009530008A (en) * | 2006-03-20 | 2009-08-27 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Ultrasound diagnosis by quantifying myocardial performance |
WO2010136584A1 (en) * | 2009-05-29 | 2010-12-02 | Institut Telecom-Telecom Paris Tech | Method for quantifying the development of diseases involving changes in the volume of bodies, in particular tumours |
US9549713B2 (en) | 2008-04-24 | 2017-01-24 | Boston Scientific Scimed, Inc. | Methods, systems, and devices for tissue characterization and quantification using intravascular ultrasound signals |
US10456105B2 (en) | 2015-05-05 | 2019-10-29 | Boston Scientific Scimed, Inc. | Systems and methods with a swellable material disposed over a transducer of an ultrasound imaging system |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4973162B2 (en) * | 2006-12-01 | 2012-07-11 | 株式会社ニコン | Image processing apparatus, image processing program, and observation system |
FR2965951B1 (en) * | 2010-10-11 | 2013-10-04 | Olea Medical | SYSTEM AND METHOD FOR ESTIMATING A QUANTITY OF INTEREST IN AN ARTERY / FABRIC / VEIN DYNAMIC SYSTEM |
WO2014064577A1 (en) * | 2012-10-23 | 2014-05-01 | Koninklijke Philips N.V. | Spatial configuration determination apparatus |
KR102211154B1 (en) * | 2013-08-06 | 2021-02-04 | 삼성전자주식회사 | Method and apparatus of diagnosis cardiac diseases based on simulation of cardiac motion |
CN104814754A (en) * | 2015-04-01 | 2015-08-05 | 王有彬 | Positron emission neurotomography device |
CN113139950B (en) * | 2021-05-08 | 2024-04-16 | 佳都科技集团股份有限公司 | Target object identification method and device |
-
2005
- 2005-01-05 EP EP05702590A patent/EP1709590A1/en not_active Withdrawn
- 2005-01-05 CN CNA200580002514XA patent/CN1910618A/en active Pending
- 2005-01-05 WO PCT/IB2005/050060 patent/WO2005071615A1/en active Application Filing
- 2005-01-05 JP JP2006548516A patent/JP2007521862A/en not_active Withdrawn
Non-Patent Citations (3)
Title |
---|
ISARD M ET AL: "CONDENSATION-conditional density propagation for visual tracking", INTERNATIONAL JOURNAL OF COMPUTER VISION KLUWER ACADEMIC PUBLISHERS NETHERLANDS, vol. 29, no. 1, 1998, pages 5 - 28, XP002329561, ISSN: 0920-5691 * |
K.M. HANSON, G.S. CUNNINGHAM, AND R.J. MCKEE: "Uncertainty Estimation in Reconstructed Deformable Models", MAXENT 96: PROC. MAXIMUM ENTROPY CONF., 1997, JOHANNESBURG, SOUTH AFRICA, pages 41 - 51, XP002329560 * |
MONTAGNAT J ET AL: "Space and time shape constrained deformable surfaces for 4D medical image segmentation", MEDICAL IMAGE COMPUTING AND COMPUTER-ASSISTED INTERVENTION - MICCAI 2000. THIRD INTERNATIONAL CONFERENCE. PROCEEDINGS (LECTURE NOTES IN COMPUTER SCIENCE VOL.1935) SPRINGER-VERLAG BERLIN, GERMANY, 2000, pages 196 - 205, XP002329562, ISBN: 3-540-41189-5 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009530008A (en) * | 2006-03-20 | 2009-08-27 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Ultrasound diagnosis by quantifying myocardial performance |
US9549713B2 (en) | 2008-04-24 | 2017-01-24 | Boston Scientific Scimed, Inc. | Methods, systems, and devices for tissue characterization and quantification using intravascular ultrasound signals |
WO2010136584A1 (en) * | 2009-05-29 | 2010-12-02 | Institut Telecom-Telecom Paris Tech | Method for quantifying the development of diseases involving changes in the volume of bodies, in particular tumours |
FR2946171A1 (en) * | 2009-05-29 | 2010-12-03 | Groupe Ecoles Telecomm | METHOD OF QUANTIFYING THE EVOLUTION OF PATHOLOGIES INVOLVING CHANGES IN BODY VOLUMES, IN PARTICULAR TUMORS |
US9026195B2 (en) | 2009-05-29 | 2015-05-05 | Institute Telecom-Telecom Paris Tech | Method for characterizing the development of pathologies involving changes in volumes of bodies, notably tumors |
US10456105B2 (en) | 2015-05-05 | 2019-10-29 | Boston Scientific Scimed, Inc. | Systems and methods with a swellable material disposed over a transducer of an ultrasound imaging system |
Also Published As
Publication number | Publication date |
---|---|
JP2007521862A (en) | 2007-08-09 |
CN1910618A (en) | 2007-02-07 |
EP1709590A1 (en) | 2006-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7822246B2 (en) | Method, a system and a computer program for integration of medical diagnostic information and a geometric model of a movable body | |
US7043063B1 (en) | Non-rigid motion image analysis | |
EP2618307B1 (en) | Method and system for determining the volume of epicardial fat from volumetric images, and corresponding computer program | |
Gee et al. | Processing and visualizing three-dimensional ultrasound data | |
US8600128B2 (en) | Image segmentation | |
US7813537B2 (en) | Motion-guided segmentation for cine DENSE images | |
US7155042B1 (en) | Method and system of measuring characteristics of an organ | |
CN107464231B (en) | System and method for determining optimal operating parameters for medical imaging | |
EP1163644B1 (en) | Method and apparatus for image processing | |
US9536318B2 (en) | Image processing device and method for detecting line structures in an image data set | |
US20030099391A1 (en) | Automated lung nodule segmentation using dynamic progamming and EM based classifcation | |
CN109661682B (en) | Image-based diagnostic system | |
CN101028187B (en) | System and method for image based physiological monitoring of cardiovascular function | |
EP1709590A1 (en) | Stochastic analysis of cardiac function | |
EP1772825A1 (en) | Method for registering images of a sequence of images, particularly ultrasound diagnostic images | |
US20060228014A1 (en) | Estimation of solitary pulmonary nodule diameters with reaction-diffusion segmentation | |
US7590273B2 (en) | Estimation of solitary pulmonary nodule diameters with a hybrid segmentation approach | |
Mazonakis et al. | Development and evaluation of a semiautomatic segmentation method for the estimation of LV parameters on cine MR images | |
EP2153408B1 (en) | Cardiac contour propagation | |
US9691159B2 (en) | Local contraction measurements | |
EP4407558A1 (en) | Obtaining a medical image at a target plane | |
Fabbri et al. | A Semi-Automated Approach for the Quantification of the Left Ventricle Chamber Volumes from Cine Magnetic Resonance Images | |
Nesser et al. | Volumetric analysis of regional left ventricular function with real-time 3D echocardiography: validation by magnetic resonance and clinical utility testing | |
Huang | Contour Tracking in Echocardiography via Sparse Modeling | |
Corsi et al. | Quantification of regional left ventricular function by real-time 3d echocardiography: validation by magnetic resonance imaging and clinical utility |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AK | Designated states |
Kind code of ref document: A1 Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW |
|
AL | Designated countries for regional patents |
Kind code of ref document: A1 Designated state(s): GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LT LU MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application | ||
WWE | Wipo information: entry into national phase |
Ref document number: 2005702590 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2557/CHENP/2006 Country of ref document: IN |
|
WWE | Wipo information: entry into national phase |
Ref document number: 200580002514.X Country of ref document: CN Ref document number: 2006548516 Country of ref document: JP |
|
WWP | Wipo information: published in national office |
Ref document number: 2005702590 Country of ref document: EP |