CN100476876C - Method for computer-assisted rebuilding heart mitral annulus - Google Patents

Method for computer-assisted rebuilding heart mitral annulus Download PDF

Info

Publication number
CN100476876C
CN100476876C CNB2007100391293A CN200710039129A CN100476876C CN 100476876 C CN100476876 C CN 100476876C CN B2007100391293 A CNB2007100391293 A CN B2007100391293A CN 200710039129 A CN200710039129 A CN 200710039129A CN 100476876 C CN100476876 C CN 100476876C
Authority
CN
China
Prior art keywords
annulus
dimensional
mitral valve
sequence
bicuspid valve
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.)
Expired - Fee Related
Application number
CNB2007100391293A
Other languages
Chinese (zh)
Other versions
CN101034480A (en
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.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong 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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CNB2007100391293A priority Critical patent/CN100476876C/en
Publication of CN101034480A publication Critical patent/CN101034480A/en
Application granted granted Critical
Publication of CN100476876C publication Critical patent/CN100476876C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention relates to a heart bicuspid valve linkcomputer assistance reconstruction method of computer application technology area, first uses man-machine interactive way select bicuspid valve link characteristic point, next when reorganized above bicuspid valve link characteristic point, uses the triaxiality ordering method, forms a sealed ordering triaxiality lattice sequence, finally based on this, uses the non-uniform three B spline curve achieve reconstruct of the bicuspid valve link structure. This invention beneficial effect is: (1) based on the man-machine interactive way direct select spot on bicuspid valve link, compare to use the algorithm division, can enhance point accuracy of selecting cingulum characteristic; (2) uses the non-uniform three B spline curve reconstruct triaxiality bicuspid valve linkstructure, because NURBS may maintain independence of the edge characteristic, and not easily influent by other control points, therefore can accurately display the bicuspid valve link such triaxiality complex structure.

Description

The method of computer-assisted rebuilding heart mitral annulus
Technical field
What the present invention relates to is the method in a kind of Computer Applied Technology field, specifically, is a kind of method of computer-assisted rebuilding heart mitral annulus.
Background technology
At present, the three-D ultrasonic video picture can be used the Computerized 3 D Model Reconstruction that turns to three-dimensional body from computer graphics.Relatively the typical models method for reconstructing is volume elements modelling (a Voxel model), this method is divided into three-dimensional body the small cubes that is arranged in order, a small cubes is exactly a volume elements, volume elements can be thought the extension of pixel in three dimensions, the volume elements of some is arranged by corresponding locus can constitute three-dimensional image, thereby reaches the purpose of 3 d objects reconstruction.
Through the literature search of prior art is found, the typical application of volume elements modelling in medical image proposes to utilize volume elements model three-dimensional reconstruction heart inner tissue structure as (" influence of pulmonary hypertension to right ventricle volume and myocardial mass estimated in the ultrasonoscopy of volume elements model three-dimensional reconstruction " that Yang Ming etc. deliver) in the 401st page of " Chinese ultrasonic image magazine " 2000 the 9th the 7th phase of volume.Though this method can the dynamic three-dimensional reconstruction heart, blood vessel structure organize gray level information, and arbitrarily dynamically heart is cut from all angles, but because ultrasonic image noise is many, the gray level information difference is not obvious, be not easy the delicate tissues that annulus of mitral valve is such and other devices and separate, so the method accuracy in three-dimensional ultrasound pattern reconstruction lobe ring is not high.
Summary of the invention
The objective of the invention is at the deficiencies in the prior art, a kind of method of rebuilding heart mitral annulus based on man-machine conversation is provided, the accuracy rate that annulus of mitral valve is rebuild has improved 20%, and reconstruction quality is not subjected to the influence of ultrasonograph quality substantially.
The present invention is achieved by the following technical solutions, and the present invention at first adopts interactive means to choose the annulus of mitral valve unique point, improves the accuracy of gathering the annulus of mitral valve feature.Secondly when the above-mentioned annulus of mitral valve unique point of reorganization, employing be the method for three-dimensional ordering, form the orderly three-dimensional lattice sequence of a sealing.Based on this orderly three-dimensional lattice sequence, utilize non-homogeneous B-spline Curve (NURBS) to reach the purpose that reconstructs the annulus of mitral valve structure at last.
Among the present invention, obtain the section plane position in each frame in the cardiac cycle, volume data is resampled along the section plane position.The doctor can select and can clear view cut into slices to mitral sagittal plane long axis direction section of heart and sagittal plane short-axis direction by man-machine interaction mode with computing machine, in every layer of section, click several obligatory points, thereby in every frame data, obtain the feature constraint point on annulus of mitral valve.
The invention has the beneficial effects as follows: (1) directly chooses point above the annulus of mitral valve based on interactive means, want than cutting apart with algorithm, can improve the accuracy of choosing lobe ring unique point, the accuracy rate that annulus of mitral valve is rebuild has improved 20%, and reconstruction quality is not subjected to the influence of ultrasonograph quality substantially; (2) utilize non-homogeneous B-spline Curve (NURBS) reconstruction of three-dimensional annulus of mitral valve structure, because NURBS can keep the independence of edge feature, and be difficult for, so can accurately show the such 3 D complex structure of annulus of mitral valve by the influence of other reference mark.That the present invention can be applicable to is ultrasonic, the reconstruction field of complicated, meticulous three-dimensional wire institutional framework on the CT, MRI image.
Description of drawings
Fig. 1 is the three-dimensional ultrasound pattern of gathering
Fig. 2 obtains the section plane sheet from ultrasonoscopy
Fig. 3 is 20 annulus of mitral valve unique points that doctor and computing machine obtain by man-machine interaction mode
Fig. 4 is the annulus of mitral valve structure of rebuilding
Fig. 5 is the annulus of mitral valve structural drawing that reconstructs from upper left corner view
Fig. 6 is the annulus of mitral valve structural drawing that reconstructs from upper right corner view
Embodiment
Below in conjunction with accompanying drawing embodiments of the invention are elaborated: present embodiment has provided detailed embodiment and process being to implement under the prerequisite with the technical solution of the present invention, but protection scope of the present invention is not limited to following embodiment.
Following examples environment of applications is a Philips SONOS 7500 three-D ultrasonic acquisition systems and an Intel Pentium IV 2.4GHz, the desktop computer of 1G internal memory, condition be the ultrasound acquisition system acquisition to three-dimensional ultrasound pattern in, can tell the difference of annulus of mitral valve and background.
At first adopt alternant way to choose unique point on the annulus of mitral valve: to obtain the section plane position in each frame in a cardiac cycle (three-dimensional ultrasound pattern of representing a certain frame as Fig. 1), along the section plane position volume data is resampled, select respectively and can clear view arrive 5 (S1 of the mitral sagittal plane long axis direction section of heart ... S5) and 5 (C1 of sagittal plane short-axis direction section ... C5), in every layer of section, click 2 obligatory points (as Fig. 2), in every frame data, just obtain 20 feature constraint points on annulus of mitral valve like this.Fig. 3 is exactly 20 obligatory point P that extract i=[x i, y i, z i] (i=1,2 ..., 20).
Next above-mentioned annulus of mitral valve unique point of recombinating, the orderly three-dimensional lattice sequence of a sealing of formation:
(1) all obligatory points is projected on the XOY plane, obtain K i=[x i, y i] (i=1,2 ..., 20)
(2) find the some K of coordinate X maximum in 20 subpoints m=[x m, y m] and the some K of X minimum n=[x n, y n] (0<m, n<20), and on XOY plane, calculate through above-mentioned 2 straight-line equation Y = KX + b , ( k = y m - y n x m - x n ) (K is the slope of straight-line equation);
(3) with subpoint K i=[x i, y i] (i=1,2 ..., 20) substitution Equation f=KX-Y+b, when f>0, represent this subpoint on straight-line equation Y=KX+b top, otherwise this subpoint is in the straight line bottom.
(4) will be in XOY plane all S the pairing S of subpoint three-dimensional feature obligatory points on straight-line equation top sort from big to small according to directions X and obtain a sequence P s, in like manner can obtain a sequence P from little to big ordering according to directions X by T T corresponding three-dimensional feature obligatory point of point of straight line bottom T, with P sAnd P TCouple together the three-dimensional order sequence P that can obtain an end to end closure i(i=1,2 ..., 20).
Utilize non-homogeneous B-spline Curve to construct the annulus of mitral valve structure at last: the three-dimensional order sequence P that can obtain an end to end closure by above step i(i=1,2 ..., 20), utilize these orderly feature constraint points as data point, do the three-dimensional curve match based on three NURBS.Utilize the cubic B-spline of NURBS method construct fit curve only with front and back 4 three-dimensional constraining points (P I-1, P i, P I+1, P I+2) relevant, because it is a kind of zonal interpolating method, rather than all points can produce the effect of control to the edge, can keep the independence of edge feature like this, and are difficult for being produced by the influence of other reference mark the mistake on can not estimating or show.Fig. 4 is the annulus of mitral valve structural drawing that builds, Fig. 5 and Fig. 6 can see three-dimensional clearly " saddle " type lobe ring structure from upper left visual angle and upper right visual angle respectively, reach the purpose of reproducing the solid shape of annulus of mitral valve under the physiological status, also provide the foundation for next step kinematic parameter of measuring the lobe ring.
Implementation result is as shown in the table:
Select the unique point number Make up the annulus of mitral valve method Reconstructed results is influenced by picture quality The annulus of mitral valve reconstruction quality
The man-machine conversation reconstruction method 95 NURBS Little Clear
The volume elements modelling 20 Automatically develop Greatly Fuzzy, discontinuous
Selection unique point number in the table is meant needs the annulus of mitral valve unique point of gathering before rebuilding annulus of mitral valve, quantitative index is the number of unique point, other parameter constants.The reconstruction structure is influenced by picture quality and is meant the situation that structure is subjected to the good bad influence of ultrasonograph quality of rebuilding.

Claims (3)

1, a kind of method of computer-assisted rebuilding heart mitral annulus, it is characterized in that, at first adopt interactive means to choose the annulus of mitral valve unique point, secondly when the above-mentioned annulus of mitral valve unique point of reorganization, adopt the method for three-dimensional ordering, form the orderly three-dimensional lattice sequence of a sealing,, utilize non-homogeneous B-spline Curve to realize reconstructing the annulus of mitral valve structure at last based on this orderly three-dimensional lattice sequence;
The described method that adopts three-dimensional ordering, the orderly three-dimensional lattice sequence of a sealing of formation is specially:
(1) all obligatory points is projected on the XOY plane, obtain a sequence K who forms by 20 pairs of bidimensional coordinates i=[x i, y i], i=1,2 ..., 20;
(2) find the some K of coordinate X maximum in 20 subpoints m=[x m, y m] and the some K of X minimum n=[x n, y n], 0<m, n<20, and on XOY plane, calculate through above-mentioned 2 straight-line equation Y=KX+b, k = y m - y n x m - x n , K is the slope of straight-line equation;
(3) with subpoint K iSubstitution Equation f=KX-Y+b when f>0, represent this subpoint on straight-line equation Y=KX+b top, otherwise this subpoint is in the straight line bottom;
(4) will be in XOY plane all S the pairing S of subpoint three-dimensional feature obligatory points on straight-line equation top sort from big to small according to directions X and obtain a sequence P s, obtain a sequence P from little to big ordering according to directions X with T T corresponding three-dimensional feature obligatory point of point of reason straight line bottom T, with P sAnd P TCouple together the three-dimensional order sequence P that obtains a new end to end closure i, i=1,2 ..., 20.
2, the method for computer-assisted rebuilding heart mitral annulus according to claim 1, it is characterized in that, described employing alternant way is chosen the annulus of mitral valve unique point, be meant: obtain the section plane position in each frame in a cardiac cycle, along the section plane position volume data is resampled, selection can clear view be cut into slices 5 to 5 of mitral sagittal plane long axis direction sections of heart and sagittal plane short-axis direction respectively, in every layer of section, click 2 obligatory points, in every frame data, just obtain 20 feature constraint points on annulus of mitral valve like this.
3, the method for computer-assisted rebuilding heart mitral annulus according to claim 1 is characterized in that, describedly utilizes non-homogeneous B-spline Curve to realize reconstructing the annulus of mitral valve structure, is meant: by the three-dimensional order sequence P of the end to end closure that obtains i, i=1,2 ..., 20, utilize these orderly feature constraint points as data point, do the three-dimensional curve match based on non-homogeneous B-spline Curve, the cubic B-spline of structure fit curve only with 4 the three-dimensional constraining point P in front and back I-1, P i, P I+1, P I+2Relevant, kept the independence of edge feature.
CNB2007100391293A 2007-04-05 2007-04-05 Method for computer-assisted rebuilding heart mitral annulus Expired - Fee Related CN100476876C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2007100391293A CN100476876C (en) 2007-04-05 2007-04-05 Method for computer-assisted rebuilding heart mitral annulus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2007100391293A CN100476876C (en) 2007-04-05 2007-04-05 Method for computer-assisted rebuilding heart mitral annulus

Publications (2)

Publication Number Publication Date
CN101034480A CN101034480A (en) 2007-09-12
CN100476876C true CN100476876C (en) 2009-04-08

Family

ID=38731019

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2007100391293A Expired - Fee Related CN100476876C (en) 2007-04-05 2007-04-05 Method for computer-assisted rebuilding heart mitral annulus

Country Status (1)

Country Link
CN (1) CN100476876C (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100547617C (en) * 2007-12-29 2009-10-07 浙江工业大学 Heart three dimensional representation method based on NURBS
US20130016907A1 (en) * 2011-07-11 2013-01-17 Jarliden Andreas E Scribble segmentation method and apparatus
CN103284758B (en) * 2013-05-20 2015-04-01 上海交通大学 Automatic cardiac mitral valve annulus detecting method
CN109069114B (en) 2017-02-16 2021-10-22 深圳迈瑞生物医疗电子股份有限公司 Ultrasonic medical detection equipment, imaging method, imaging system and display terminal
CA3055394C (en) * 2017-03-07 2024-01-02 Cd Med S.R.L. Method for generating a mitral repair ring, and mitral repair ring
CN109087357B (en) * 2018-07-26 2021-06-29 上海联影智能医疗科技有限公司 Scanning positioning method and device, computer equipment and computer readable storage medium
CN112991522B (en) * 2021-03-30 2023-03-24 华南理工大学 Personalized automatic modeling method, system and equipment for mitral valve

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
三维超声心动图空间任意方向断层图像的重建. 周伟等.上海交通大学学报,第36卷第12期. 2002 *
体元模型法三维超声成像及其应用. 李秀平.安阳师范学院学报. 2006 *
基于切片图像的血管三维重建方法. 韩西安等.装备指挥技术学院学报,第13卷第3期. 2002 *

Also Published As

Publication number Publication date
CN101034480A (en) 2007-09-12

Similar Documents

Publication Publication Date Title
CN100476876C (en) Method for computer-assisted rebuilding heart mitral annulus
CN107392994B (en) Three-dimensional rebuilding method, device, equipment and the storage medium of coronary artery blood vessel
US5435310A (en) Determining cardiac wall thickness and motion by imaging and three-dimensional modeling
US8079957B2 (en) Synchronized three or four-dimensional medical ultrasound imaging and measurements
DE102018216296A1 (en) Determination of the measuring point in medical diagnostic imaging
US7853304B2 (en) Method and device for reconstructing two-dimensional sectional images
CN107230206A (en) A kind of 3D Lung neoplasm dividing methods of the super voxel sequence lung images based on multi-modal data
CN101625767B (en) Method for rebuilding point-cloud type three-dimensional surface of nonparallel outline medical image
US8970581B2 (en) System and method for interactive contouring for 3D medical images
CN110060337B (en) Carotid artery ultrasonic scanning three-dimensional reconstruction method and system
CN103679810A (en) Method for three-dimensional reconstruction of liver computed tomography (CT) image
CN107427279A (en) Use the Ultrasonic Diagnosis of the cardiac function of the cardiac module chamber with user's control
CN104574520A (en) Modification of a hollow organ representation
Scheipers et al. 3-D ultrasound volume reconstruction using the direct frame interpolation method
CN110807770A (en) Medical image processing, recognizing and displaying method and storage medium
CN109242947A (en) Three-dimensional ultrasound pattern display methods
CN101996415B (en) Three-dimensional modeling method for eyeball
CN105078492B (en) The bearing calibration of sawtooth artifact and device in digital galactophore cross sectional reconstruction
Eulzer et al. Temporal views of flattened mitral valve geometries
US10891776B2 (en) Imaging system and method
EP1903503A2 (en) Method, system and computer program product for providing user-customizable standardized anatomical viewing protocols for volumetric data
CN108510506A (en) A kind of tubular structure image partition method
Subramanyan et al. A procedure to average 3D anatomical structures
CN109636915B (en) Medical image three-dimensional reconstruction isosurface grid extraction method based on edge growth
CN101174331A (en) Maximum density projection generating method for medical image

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20090408

Termination date: 20120405