CN101422352A - Interactive coronary artery virtual angioscope implementation method - Google Patents
Interactive coronary artery virtual angioscope implementation method Download PDFInfo
- Publication number
- CN101422352A CN101422352A CNA2008100800020A CN200810080002A CN101422352A CN 101422352 A CN101422352 A CN 101422352A CN A2008100800020 A CNA2008100800020 A CN A2008100800020A CN 200810080002 A CN200810080002 A CN 200810080002A CN 101422352 A CN101422352 A CN 101422352A
- Authority
- CN
- China
- Prior art keywords
- rightarrow
- coronary artery
- path
- image
- ivus
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
The invention relates to a realization method for an interactive coronary artery virtual angioscope, which pertains to the technical field of medical imaging and aims at solving the visible problem of making a diagnosis and giving treatment of coronary artery. The technical proposal is that the realization method obtains the three-dimensional model of blood vessels by means of combining three-dimensional geometrical feature information of lumen by a nearly orthorhombic X-ray visualization picture of coronary artery with the data of the section of lumen by lumen endobronchial ultrasonography, then adopts virtual reality modeling language to interactively describe the model of blood vessels and realizes the visualization of coronary artery in the endoscope roaming mode. The realization method realizes interactive visit and display of the three-dimensional model of blood vessels and provides an ideal platform for the research on the development of coronary atherosclerosis pathological changes, the visible diagnosis and treatment of coronary heart diseases, the involvement in the evaluation of curative effect, and the like, and the training of medical personnel.
Description
Technical field
The present invention relates to a kind of implementation method of the interactive coronary artery virtual angioscope that merges based on many formation methods, belong to the medical imaging technology field.
Background technology
Virtual endoscopic techniques (Virtual Endoscope, VE) be to utilize medical image as initial data, technology such as comprehensive utilization Digital Image Processing, computer graphics, visualization in scientific computing, virtual reality, reconstruction of three-dimensional images forms the virtual human body tissue; Then viewpoint is inserted in the organ cavity that reconstructs, simulated the tract endoscopy realistically by navigation or roaming technology and pseudo-color technology.
The at present clinical coronary artery interventional imaging means that generally adopt are X ray coronarography (CAG, Coronary ArteryAngiography) and intravascular ultrasound (IVUS, Intravascular Ultrasound), the two carries out simultaneously, and angiography and intravascular ultrasound show the configuration of catheterbased US probe in intraluminal position and corresponding blood vessel wall respectively synchronously.CAG and IVUS have the projected outline of advantage and not enough complementary characteristics: CAG reflection lumen of vessels after by the contrast agent filling, diseases such as energy diagnosing ischemia heart disease and malformation of coronary artery, and significant, but can not provide the structural information and the lesion degree of blood vessel wall to thrombolytic in the coronary artery, the treatment of PTCA intervene operations such as (percutaneous transluminal coronary angioplasties); IVUS can clear demonstration blood vessel transverse section, carries out histological classification according to the speckle acoustic features, finds the vascular lesion that CAG can not show, observes the fuzzy pathological changes of crotch or blood vessel overlapping etc.But owing to adopt the high frequency ultrasound probe, influenced investigation depth, can only measure, can not enter serious narrow tube chamber, and can't determine the axial location and the direction in space in cross section a certain section lesion vessels.
Get involved the imaging detection in addition and also comprise the coronary artery mirror, it is an a kind of small endoscopic imaging technology of utilizing optical fiber technology.But this technology is not accepted extensively clinically, and reason comprises: the morphologic data of luminal surface can only be provided, can not observe the pathological changes DEEP STRUCTURE in the tube wall, can not carry out the quantitative analysis of stenosis and blood flow; The pathological changes that can not be used for pathological changes and anterior descending branch and the Zhi Jinduan that circles round at video picture aorta-coronary ostium place; The blood flow that enters from side opening can make the visual field fuzzy; Conduit lacks navigability, has limited the video picture scope; Need temporary transient plug flow in the checking process, can cause the generation of myocardial ischemia etc.
At present, Non-Invasive cardiovascular image check mainly comprises CTA (CT Angiography) and MRCA (Magnetic ResonanceCoronary Angiography).But the main limitation of heart CT examination is to be easy to generate pseudomorphism, influences picture quality.For MRCA, because arteria coronaria blood vessel itself is thin, distortion and structure are complicated, and has the influence of heartbeat and breathing, signals such as fatty tissue and cardiac muscular tissue can influence the result of its video picture around the arteria coronaria.Though there are not lonizing radiation in the MR checking process simultaneously, comparatively safe, the influence of noise is arranged, the safety of some metal implants (as artificial metal lobe, cardiac pacemaker etc.) also receives publicity.In a word, because deficiency and technical limitation that image-forming principle caused, feasible CTA up to the present and MRCA generally only can be used as a kind of system of selection to heart overall merit, or as the examination measure of coronary heart disease conduit contrast examination, reduce the unnecessary wound inspection, and, can not replace the image check method of intervention property on to the clinical diagnosis and treatment of coronary heart disease fully to the non-invasive follow-up investigation of operation on heart or interventional therapy effect.
In sum, the main image method that the CAG of intervention property and IVUS remain clinical diagnosis and treatment coronary heart disease, and also the two has advantage and not enough complementary characteristics.At present also do not have a kind of interactive coronary artery virtual angioscope system, can realize that the coronary artery of scope roam mode is visual based on CAG and IVUS image co-registration.
Summary of the invention
The objective of the invention is to overcome the deficiencies in the prior art, a kind of implementation method that can realize the visual interactive coronary artery virtual angioscope of coronary artery of scope roam mode is provided.
The alleged problem of the present invention realizes with following technical proposals:
A kind of implementation method of interactive coronary artery virtual angioscope, it is to merge mutually with the tube chamber cross-sectional data that is obtained by intravascular ultrasound by the tube chamber three-dimensional geometry shape information that will be obtained by the X ray coronarogram picture of nearly orthogonal, obtain the threedimensional model of blood vessel, use Virtual Reality Modeling Language (VRML) alternatively to describe vascular pattern then, the coronary artery of realizing the scope roam mode is visual, and concrete steps are as follows:
A, gather the intravascular ultrasound IVUS and the X ray coronarography CAG image of vessel segment simultaneously:
Mechanical type ultrasound catheter probe is placed the far-end of vessel segment, remove in the process of guide wire at the uniform velocity equidistant time, utilize the intravascular ultrasound imaging instrument to gather equidistant IVUS image sequence in the mode of electrocardio ECG gate, promptly with the R ripple of electrocardiosignal as triggering, only identical cardiac phase place images acquired in each cardiac cycle can solve the motion artifacts problem in the coronary artery IVUS image sequence.Simultaneously, utilize C type arm single face X ray angiography machine to return the CAG image that the starting point of removing the path is taken two width of cloth near normal angles of record same cardiac state at conduit;
The threedimensional model that B, the IVUS that utilizes collection and CAG image are set up blood vessel;
C, utilization Virtual Reality Modeling Language are realized coronary artery revascularization result's scope roam mode visual.
The concrete steps of threedimensional model that the implementation method of above-mentioned interactive coronary artery virtual angioscope, the IVUS of described utilization collection and CAG image are set up blood vessel are as follows:
A, according to the CAG image of two near normal angles, three-dimensional reconstruction goes out returning of ultrasound catheter and removes the path;
B, three-dimensional reconstruction goes out vessel lumen from the CAG image:
With the 3D rodding that reconstructs footpath two CAG imaging plane back projections to the left and right, obtain corresponding 2D path, for each point on the 2D path, by on direction perpendicular to the path, seek two maximum of shade of gray, obtain the vessel lumen left and right edges, the cross section at the hypothesis tube chamber is under the prerequisite of ellipse then, finish the three-dimensional reconstruction of whole vessel lumen, this result is used for the follow-up direction in space of determining each frame ultrasonoscopy;
The extraction of c, ivus image sequence medium vessels wall profile:
In first two field picture, manually select the several characteristic point on the inside and outside film profile of blood vessel wall, put formed polygon as initial position to connect these, obtain the profile of the inside and outside film of blood vessel wall by the snake distortion, the speckle that is partitioned into blood vessel wall and may exists, for subsequent frame, then, finish cutting apart to continuous multiple frames IVUS image with the extraction result of former frame initial position as snake;
D, determine the axial location of each frame ivus image:
According to the acquisition order and the spacing of IVUS image, the 3D conduit that the edge reconstructs returns removes the path with each frame IVUS image sequence arrangement, determines the axial location of each two field picture;
E, determine the dimensional orientation of each frame ivus image:
On the 3D conduit path after the reconstruction, set up the local coordinate system of each frame ultrasonoscopy, i.e. Frenet-Serret frame, three secondary method vector b of the tangent vector t of the coordinate axes unit of being respectively, the master of unit method vector n and unit, the position of conduit is positioned at the center of IVUS image;
Each frame ultrasonoscopy is rotated to its correct direction to determine the dimensional orientation of ivus image around conduit:
The centrifugal vector of representing the deviation of gravity center catheter position of blood vessel wall profile with ρ, the elliptic contour of the vessel lumen that goes out from the CAG image reconstruction is projected on the corresponding ultrasonoscopy, represent that with μ vessel lumen elliptic contour centrage departs from the centrifugal vector of catheter position, ε is the mould of vectorial ρ, θ is the angle of ρ and μ, determine the dimensional orientation of ultrasonic image sequence with the statistics optimization method, purpose is to make the θ minimum:
Set the moving window of a fixed width w, carry out statistical analysis in this window, there is n in the ultrasonic image sequence for the N frame is formed
wThe individual moving window of=N-(w-1), at each the window's position m place, accumulative total eccentric distance ∑ ε
m, the eccentric angle meansigma methods of weighting θ
mAnd the weighted standard deviation σ (θ of eccentric angle
m) can calculate by following formula respectively:
Utilize these numerical value, in each the window's position place, computed reliability weight factor: r
m=∑ ε
m/ σ (ε
m), give bigger weight factor in the bigger position of eccentric distance, limit σ (θ simultaneously
m) bigger position, calculate one by following formula and proofread and correct eccentric angle θ
Corr:
And apply it in all images of sequence, thereby obtain the dimensional orientation of each two field picture;
The drafting of vessel lumen surfaces externally and internally is finished in f, utilization based on the surface extraction method of nurbs surface match.
The implementation method of above-mentioned interactive coronary artery virtual angioscope, the visual concrete steps of described utilization Virtual Reality Modeling Language realization coronary artery revascularization result's scope roam mode are as follows:
1. roam the calculating in path:
Returning the three-dimensional reconstruction result of removing the path according to ultrasound catheter, is the current view point position with the coordinate Pi along the i frame IVUS image acquisition point of conduit, P
I+1Be next viewpoint position, with the direction vector of the negative semiaxis of z axle
Be the initial viewpoint direction, then in VRML, the position vector of current view point is
, rotating shaft for the current view point direction is
Unit vector
, rotating shaft is
, the anglec of rotation is
2. the IVUS pixel data is inserted in the virtual scene, and adopts translucent display mode to show;
3. develop interactively user graphical interface.
The implementation method of above-mentioned interactive coronary artery virtual angioscope, described ultrasound catheter returns the three-dimensional rebuilding method of removing the path: at first set up the perspective projection imaging model of CAG system two near normal angles, again according to the distance and the angle parameter of synchronous recording in angiographic procedure, obtain the geometric transformation of imaging system, adopt three-dimensional snake modelling technique then, snake directly is out of shape in the space, finishes conduit and returns the three-dimensional reconstruction of removing the path;
The implementation method of above-mentioned interactive coronary artery virtual angioscope, the span of angle is 60 ° to 120 ° between two acquisition angles of described coronarogram picture, and only returns the starting point of removing the path at ultrasound catheter and take a pair of contrastographic picture.
The blood vessel space geometry information that the present invention will be obtained by the single face contrastographic picture of two nearly orthogonal angles combines with the tube chamber cross section information that is obtained by ivus image, make full use of the complementarity of two kinds of imaging means, finished the accurate three-dimensional reconstruction of blood vessel, and the utilization Virtual Reality Modeling Language realizes that the coronary artery of scope roam mode is visual.The present invention has realized the interactive visit of three-dimensional vascular pattern and demonstration, for the visual diagnosis and treatment of the development of coronary atherosclerosis pathological changes, coronary heart disease, to the research of interventional therapy effect assessment etc., and medical worker's training provides an ideal platform.
Description of drawings
The invention will be further described below in conjunction with accompanying drawing.
Fig. 1 is the flow chart according to the three-dimensional reconstruction blood vessel of the inventive method;
Fig. 2 is according to CAG of the inventive method and IVUS image acquisition sketch map;
Fig. 3 is according to the radiography system of the inventive method imaging sketch map two angles;
Fig. 4 is the definite sketch map according to each frame ultrasonoscopy relative bearing of the inventive method;
Fig. 5 is according to ultrasonoscopy eccentric distance of the inventive method and eccentric angle sketch map;
Fig. 6 is the definite sketch map according to the roaming viewpoint position of the inventive method;
Fig. 7 is the definite sketch map according to the roaming viewpoint direction of the inventive method.
Each symbol is among the figure: Image A, Image B, imaging plane; s
1, s
2, the position of x-ray source focus in twice angiographic procedure; s
1x
1y
1z
1, with s
1Space coordinates for initial point; s
2x
2y
2z
2, with s
2Space coordinates for initial point; U
1V
1O
1, the rectangular coordinate system on the imaging plane A; U
2V
2O
2, the rectangular coordinate system on the imaging plane B; D
1, s
1Vertical dimension to imaging plane A; D
2, s
2Vertical dimension to imaging plane B; Point on P, the space blood vessel; p
1, the projection of P point on imaging plane A; p
2, the projection of P point on imaging plane B; u
1, p
1At coordinate system U
1V
1O
1Interior abscissa; v
1, p
1At coordinate system U
1V
1O
1Interior vertical coordinate; u
2, p
2At coordinate system U
2V
2O
2Interior abscissa; v
2, p
2At coordinate system U
2V
2O
2Interior vertical coordinate; The spatial parameter curve in c (s), expression 3D conduit path; The position of conduit in C, the ultrasonoscopy, it also is the center of ultrasonoscopy; O
C, the oval cross section profile center (promptly in the hypothesis vessel cross-sections when being ellipse, pairing blood vessel center line position in the three-dimensional reconstruction based on CAG); O
I, the tube chamber cross section profile that from ultrasonoscopy, extracts center of gravity; ρ, ρ=O
1-C is the centrifugal vector of the deviation of gravity center conduit of ultrasound contour; μ, μ=O
C-C is the centrifugal vector that vessel centerline departs from conduit; The angle of θ, ρ and μ; P
i, P
I+1, viewpoint;
Position vector;
Rotating shaft; φ
i, the anglec of rotation;
The direction vector of the negative semiaxis of z axle;
Viewpoint direction to be asked;
The unit vector of y axle also is the direction of giving tacit consent among the VRML that makes progress; ε
i,
The anglec of rotation in the x-y plane.
Used symbol: t, unit tangent vector in the literary composition; N, the master of unit method vector; B, the secondary method vector of unit; The mould of ε, vectorial ρ; W, moving window width;
The eccentric angle meansigma methods of weighting; σ (θ
m), the weighted standard deviation of eccentric angle; r
m, the reliability weight factor; θ
Corr, proofread and correct eccentric angle.
The specific embodiment
Describe step of the present invention in detail below in conjunction with accompanying drawing:
(1) image acquisition:
Collecting device comprises C type arm single face X ray angiography machine and intravascular ultrasound imaging instrument.
Referring to Fig. 2, IVUS and CAG imaging are carried out simultaneously.Conventional brachial artery puncture through right femoral artery or upper arm, the row selectivity pulse prefixing contrast.Under the guidance of radioscopy image, insert high frequency ultrasound probe conduit, to the blood vessel far-end.Ultrasonic probe is connected the pseudo-movie queen of removal with ultrasonic imaging device, at the uniform velocity equidistantly returns through motor control and remove conduit.When doing 360 ° of rotations with 1800 rev/mins, the probe conduit obtains the real-time blood vessel tangent plane picture of 30 frame/seconds continuously.Adopt clinical commonly used, allow patient return the method for holding the breath in the process of removing at conduit, reduce respirometric influence.Adopt the mode of ECG gate to gather ultrasonoscopy, thereby reduce the influence of heart movement.
Only return the starting point of removing the path, adopt the mode of ECG gate, take the contrastographic picture of a pair of near normal angle at corresponding cardiac phase place at conduit.Because adopt mechanical type ultrasound catheter probe, ultrasonic transducer is positioned at a flexible axle center head end, the axle center is at epitheca pipe internal rotation, and the sheath pipe is fixed, therefore can guarantee back to remove the stable of path.Record radiography angle and x-ray source focus are to the distance of receiving screen in the imaging process.
(2) extraction and the three-dimensional reconstruction at conduit path and tube chamber edge in the contrastographic picture:
The present invention at first sets up the perspective projection imaging model (accompanying drawing 3) of CAG system two near normal angles.Afterwards, according to the distance and the angle parameter of synchronous recording in angiographic procedure, obtain the geometric transformation of imaging system.Utilize three-dimensional snake modelling technique then, finish the three-dimensional reconstruction in conduit path.
The snake model claims movable contour model (active contour model) again, be a kind of distorted pattern technology (the Kass M that proposed in 1987 by Kass etc., Witkin A, Terzopoulos T.Snakes:active contour models.International Journal of Computer Vision, 1987,1 (4): 321-331), application is very extensive in image processing field in recent years, finishes image segmentation, coupling and motion tracking.
The specific implementation method is: the initial position of snake adopts manually gets an acquisition, promptly in a projection of conduit, manually choose some sampled points (generally choosing back starting point, terminal point and 3~4 intermediate points of removing the path gets final product), obtain these corresponding point in another projection according to outer utmost point constraint then.Obtain their three-dimensional coordinate respectively by these several groups of corresponding point, connect these 3D points with straightway, the gained broken line is as the initial position of 3D snake.
The energy function of snake model is:
C (s)=(x (s), y (s), z (s)) wherein, s ∈ [0,1] are the B-spline Curve of expression conduit.Internal energy E in the formula (1)
IntExpression formula be:
E
int(c(s))=(α|c′(s)|
2+β|c"(s)|
2)/2 (2)
Wherein c ' (s) and c " (s) be respectively single order and the second dervative of c (s).Internal energy guarantees the continuous and smooth of curve.
External energy function E
ExtBe to guarantee snake astringent external force, comprise two parts, correspond respectively to left and right sides projection, guarantee that the projection of three-dimensional curve on two angle imaging planes is positioned at corresponding conduit projection place just:
(3)
I wherein
L(u
1, v
1) and
Be respectively the gray scale and the shade of gray value of left subpoint; I
R(u
2, v
2) and
I
R(u
2, v
2) be respectively the gray scale and the shade of gray of right subpoint.Because in the contrastographic picture, the gray value of blood vessel is littler than background, thus weight coefficient γ get on the occasion of.According to the geometrical relationship and the outer polar curve restriction relation of perspective projection imaging, can derive u
1, v
1, u
2And v
2All be spatial point three-dimensional coordinate c=(x
1, y
1, z
1) function:
[u
1?v
1]
T=F
L(c),[u
2?v
2]
T=F
R(c) (4)
Afterwards, minimize by the energy function that makes formula (1), the three-dimensional axis of conduit has just been determined in the final position of snake curve.This method has been avoided having improved reconstruction precision and arithmetic speed based on the pointwise coupling between two angles of outer utmost point constraint.
According to the geometric transformation of imaging system,, obtain corresponding 2D path with the 3D rodding that reconstructs footpath two CAG imaging plane back projections to the left and right.For each point on the 2D path, by on direction, seek two maximum of shade of gray perpendicular to the path, finish extraction to the vessel lumen left and right edges.Afterwards, be under the prerequisite of ellipse at the cross section of supposing tube chamber, finish the three-dimensional reconstruction of whole tube chamber, this result is that the direction in space of follow-up definite each frame ultrasonoscopy is used.
(3) extraction of ivus image sequence medium vessels wall profile:
The present invention adopts in conjunction with the snake model of dynamic programming and finishes extraction to adventitia profile in the IVUS image sequence medium vessels wall.The operator only needs the several characteristic point on the manual select target profile in first frame, connects these and puts the initial position of formed polygon as snake.For subsequent frame, with the extraction result of former frame initial position, finish the cutting apart of continuous multiple frames IVUS image as next frame snake, can save computation time greatly.
(4) fusion of IVUS and CAG data:
Here mainly need solve two problems: 3D axial location and the dimensional orientation of determining each IVUS frame.
(4.1) determining of the three-dimensional axial location of ultrasonoscopy:
In gathering the process of ultrasonoscopy, the mode that adopts motor to drive at the uniform velocity equidistantly pulls straight conduit from far-end to near-end.The speed of conduit is pulled out in adjusting, can regulate the tangent plane spacing as required.Adopt the CAG image reconstruction to go out after the axis of conduit, according to axially each frame IVUS image sequence being arranged, can determine the axial location of each two field picture according to known tangent plane spacing.
(4.2) each frame ultrasonoscopy dimensional orientation determines
The present invention utilizes a kind of non-iterative statistics optimization method to calculate the dimensional orientation of each frame ultrasonoscopy.At first on the 3D conduit path after the reconstruction, set up the local coordinate system of each frame ultrasonoscopy, it is the Frenet-Serret frame, three secondary method vector b (accompanying drawing 4) of the tangent vector t of the coordinate axes unit of being respectively, the master of unit method vector n and unit, zero is the position of conduit in the IVUS image.Behind the three-dimensional reconstruction of finishing the conduit path, can obtain its 3D curvilinear equation c (s), according to differential geometric knowledge, t, n and b can be calculated as follows according to curvilinear equation:
Wherein " * " represents the multiplication cross of vector; "
The dot product of expression vector; C ' (s) and c " (s) be respectively single order and the second dervative of c (s).
The position of conduit is positioned at the center of IVUS image, and the center of gravity of the objective contour that is partitioned into does not generally overlap with catheter position, and as shown in Figure 5, wherein the C point is represented conduit, O
CBe the center of elliptic contour (promptly when the hypothesis vessel cross-sections is ellipse, pairing blood vessel center line position in the three-dimensional reconstruction based on CAG), O
1Center of gravity for the vascular cross-section profile that from ultrasonoscopy, extracts.Adopt centrifugal vectorial ρ to represent the degree of the deviation of gravity center catheter position of profile: ρ=O
I-C.
Because vessel centerline and conduit path do not overlap, the elliptic contour orientation that reconstructs at the ultrasonoscopy profile at blood vessel same position place with based on contrastographic picture is inconsistent, and elliptic contour is projected on the corresponding ultrasonoscopy.The same centrifugal vectorial μ of elliptic contour that adopts represents that vessel centerline departs from the degree of catheter position: μ=O
C-C.
The mould ε of the matching error availability vector ρ of ultrasonoscopy and the angle theta of ρ and μ are represented.The present invention utilizes the statistics optimization method to determine that the absolute orientation of ultrasonic image sequence, purpose are the angle theta minimums that makes between the centrifugal vector of elliptic contour and ultrasound contour.Set the moving window of a fixed width w, in this window, carry out statistical analysis.There is n in ultrasonic image sequence for the N frame is formed
WThe individual moving window of=N-(w-1).At each the window's position m place, accumulative total eccentric distance ∑ ε
m, the eccentric angle meansigma methods of weighting θ
mAnd the weighted standard deviation σ (θ of eccentric angle
m) can calculate by following formula respectively:
Utilize these numerical value, in each the window's position place, computed reliability weight factor: r
m=∑ ε
m/ σ (ε
m).Give bigger weight factor in the bigger position of eccentric distance, limit σ (θ simultaneously
m) bigger position.Calculate one by following formula and proofread and correct eccentric angle θ
Corr:
And apply it in all images of sequence, thereby obtain the dimensional orientation of each two field picture.
(5) match of lumen of vessels surfaces externally and internally
After the IVUS image sequence being finished edge extracting and being determined the locus of each frame, the present invention adopts NURBS (non-uniform rational B-spline) surface fitting to return the sampled point of removing on correct each cross section arranged in path along three-dimensional, obtains the continuous three-dimensional blood vessel surface.
(6) interactive coronary artery virtual endoscopic system
Utilize virtual reality modeling language to show the reconstructed results of scope roam mode, not only can show the overall appearance of rebuilding the back vessel segment, and can show major axis profilograph picture.
The demonstration of IVUS view data and the exploitation of interactively user graphical interface in the drafting of the tube chamber surfaces externally and internally that comprise the calculating of roaming the path, reconstructs, the virtual scene.
(6.1) calculating in roaming path
The roaming path is the sequence of forming in a series of viewpoints of target blood intracavity.For each viewpoint, all need to determine its position and direction, wherein direction adopts observer place local coordinate system to represent around the rotation of arbitrary axis.
Determining of viewpoint position: in VRML, 1 P on the roaming path
iThe position of (being viewpoint) and direction ternary number
Represent, wherein
Be position vector,
Be rotating shaft, φ
iIt is the anglec of rotation (as shown in Figure 6).
The roaming of virtual observer in target blood both can be removed the path along returning of IVUS conduit and carried out, and also can carry out along the tube chamber centrage.So
Can directly obtain according to both three-dimensional reconstruction result: for the former,
Be exactly coordinate along the i frame IVUS image acquisition point of conduit, the also i.e. centre coordinate of this two field picture; For the latter,
The barycentric coodinates of the tube chamber profile that from i frame IVUS image, is partitioned into exactly.Because the rigidity and the seriality of conduit are all good than the tube chamber centrage that calculates, so the present invention adopts first method to determine
Thereby, obtain more continuous slick animation effect.
Determining of viewpoint direction: the initial value of viewpoint direction is set at the negative semiaxis along the z axle, as shown in Figure 6, and wherein
The direction vector of the negative semiaxis of expression z axle also is the initial viewpoint direction.The path carries out because roaming is to remove along returning of conduit, therefore according to definite method of above-mentioned viewpoint position, by the next as can be known viewpoint P of current view point Pi
I+1Position vector
Thereby, obtain vector
, its unit vector is to be asked looking
The point direction
With
Cross product be rotating shaft:
Promptly
Perpendicular to
With
The plane that is determined, as shown in Figure 7, the inceptive direction of viewpoint
Around rotating shaft
Rotation φ
iThe angle promptly obtains the current viewpoint direction
, promptly by current view point P
iPoint to next viewpoint P
I+1Direction, the anglec of rotation is:
(7)
By the computing formula of vector cross product as can be known,
The z component be 0, the expression
In the x-y plane, so
Also can be by the unit vector of y axle
At x-y plane internal rotation ε
iObtain:
R wherein
z(ε
i) expression is around the spin matrix of z axle,
It is the direction of giving tacit consent among the VRML that makes progress.
Terminal point in the roaming path is promptly for viewpoint set { P
0, P
1..., P
N-1In 1 P
i, when i=n-1, owing to there is not P
I+1So unavailable preceding method calculates the rotating shaft and the anglec of rotation, at this moment:
(6.2) demonstration of IVUS pixel data among the VRML
The present invention adopts translucent display mode to the IVUS pixel data that inserts in the virtual scene, and promptly the transparence value of each pixel is not same constant in the frame IVUS image, but depends on the position of pixel in image and its gray value.In the ultrasonoscopy except blood vessel wall and speckle, other structure all should be sightless in virtual scope scene, utilize the two-dimentional segmentation result of front this moment to ultrasonoscopy, to be set to full impregnated bright for the pixel of echo signal beyond expression tube chamber and the adventitia, allows the roaming path to pass through these zones unobstructedly.And for tube wall and these area-of-interests of speckle, its transparence value depends on the gray value of pixel, for example, and the speckle that bright echo signal expresses possibility and exists, so its opacity value should be set to higher numerical value; The tube chamber of other blood vessel may be represented or not have echogenic other structure that its opacity value should be made as lower numerical value in dark space in the image.
(6.3) interactively user graphical interface
The present invention designs and develops out simple and clear clear, convenient, flexible user's control panel in the VRML environment, make it can finish following function: 1. user's opening and closing control panel at any time, and when opening, with as far as possible not the shelter target scene be principle.2. when advancing, at certain viewpoint place, can switch between different display modes by the user along the roaming path for virtual observer, for example: the IVUS image that obtains at this point according to correct direction and position display; The luminal surface (can open or close translucent IVUS frame simultaneously) that perhaps only shows this some place; Perhaps show the luminal surface of having finished coloud coding, wherein coloud coding is represented measures of quantization result etc.Simultaneously, the user can enter or withdraw from the observing pattern of virtual scope at any time, shows the overall appearance of rebuilding the back vessel segment, perhaps major axis longitudinal section image.3. can adjust the speed and the direction of roaming, virtual observer can stop in intraluminal optional position.
Claims (5)
1, a kind of implementation method of interactive coronary artery virtual angioscope, it is characterized in that, it is to merge mutually with the tube chamber cross-sectional data that is obtained by intravascular ultrasound by the tube chamber three-dimensional geometry shape information that will be obtained by the X ray coronarogram picture of nearly orthogonal, obtain the threedimensional model of blood vessel, use Virtual Reality Modeling Language (VRML) alternatively to describe vascular pattern then, the coronary artery of realizing the scope roam mode is visual, and concrete steps are as follows:
A, gather the intravascular ultrasound IVUS and the X ray coronarography CAG image of vessel segment simultaneously:
Mechanical type ultrasound catheter probe is placed the far-end of vessel segment, remove in the process of guide wire at the uniform velocity equidistant time, utilize the intravascular ultrasound imaging instrument to gather equidistant IVUS image sequence at identical cardiac phase place in the mode of ECG gate, simultaneously, utilize C type arm single face X ray angiography machine to return the CAG image that the starting point of removing the path is taken two width of cloth near normal angles of record same cardiac state at conduit;
The threedimensional model that B, the IVUS that utilizes collection and CAG image are set up blood vessel;
C, utilization Virtual Reality Modeling Language are realized coronary artery revascularization result's scope roam mode visual.
According to the implementation method of the described interactive coronary artery virtual angioscope of claim 1, it is characterized in that 2, the concrete steps of threedimensional model that the IVUS of described utilization collection and CAG image are set up blood vessel are as follows:
A, according to the CAG image of two near normal angles, three-dimensional reconstruction goes out returning of ultrasound catheter and removes the path;
B, three-dimensional reconstruction goes out vessel lumen from the CAG image:
With the 3D rodding that reconstructs footpath two CAG imaging plane back projections to the left and right, obtain corresponding 2D path, for each point on the 2D path, by on direction perpendicular to the path, seek two maximum of shade of gray, obtain the vessel lumen left and right edges, the cross section at the hypothesis tube chamber is under the prerequisite of ellipse then, finish the three-dimensional reconstruction of whole vessel lumen, this result is used for the follow-up direction in space of determining each frame ultrasonoscopy;
The extraction of c, ivus image sequence medium vessels wall profile:
In first two field picture, manually select the several characteristic point on the inside and outside film profile of blood vessel wall, put formed polygon as initial position to connect these, obtain the profile of the inside and outside film of blood vessel wall by the snake distortion, the speckle that is partitioned into blood vessel wall and may exists, for subsequent frame, then, finish cutting apart to continuous multiple frames IVUS image with the extraction result of former frame initial position as snake;
D, determine the axial location of each frame ivus image:
According to the acquisition order and the spacing of IVUS image, the 3D conduit that the edge reconstructs returns removes the path with each frame IVUS image sequence arrangement, determines the axial location of each two field picture;
E, determine the dimensional orientation of each frame ivus image:
On the 3D conduit path after the reconstruction, set up the local coordinate system of each frame ultrasonoscopy, i.e. Frenet-Serret frame, three secondary method vector b of the tangent vector t of the coordinate axes unit of being respectively, the master of unit method vector n and unit, the position of conduit is positioned at the center of IVUS image;
Each frame ultrasonoscopy is rotated to its correct direction to determine the dimensional orientation of ivus image around conduit:
The centrifugal vector of representing the deviation of gravity center catheter position of blood vessel wall profile with ρ, the elliptic contour of the vessel lumen that goes out from the CAG image reconstruction is projected on the corresponding ultrasonoscopy, represent that with μ vessel lumen elliptic contour centrage departs from the centrifugal vector of catheter position, ε is the mould of vectorial ρ, θ is the angle of ρ and μ, determine the dimensional orientation of ultrasonic image sequence with the statistics optimization method, purpose is to make the θ minimum:
Set the moving window of a fixed width w, carry out statistical analysis in this window, there is n in the ultrasonic image sequence for the N frame is formed
wThe individual moving window of=N-(w-1), at each the window's position m place, accumulative total eccentric distance ∑ ε
m, the eccentric angle meansigma methods of weighting θ
mAnd the weighted standard deviation σ (θ of eccentric angle
m) can calculate by following formula respectively:
Utilize these numerical value, in each the window's position place, computed reliability weight factor: r
m=∑ ε
m/ σ (ε
m), give bigger weight factor in the bigger position of eccentric distance, limit σ (θ simultaneously
m) bigger position, calculate one by following formula and proofread and correct eccentric angle θ
Corr:
And apply it in all images of sequence, thereby obtain the dimensional orientation of each two field picture;
The drafting of blood vessel surface is finished in f, utilization based on the surface extraction method of nurbs surface match.
According to the implementation method of claim 1 or 2 described interactive coronary artery virtual angioscopes, it is characterized in that 3, the visual concrete steps of described utilization Virtual Reality Modeling Language realization coronary artery revascularization result's scope roam mode are as follows:
1. roam the calculating in path:
Return the three-dimensional reconstruction result of removing the path according to ultrasound catheter, with coordinate P along the i frame IVUS image acquisition point of conduit
iBe current view point position, P
I+1Be next viewpoint position, with the direction vector of the negative semiaxis of z axle
Be the initial viewpoint direction, then in VRML, the position vector of current view point is
Rotating shaft for the current view point direction is
Unit vector
Rotating shaft is
The anglec of rotation is
2. the IVUS pixel data is inserted in the virtual scene, and adopts translucent display mode to show;
3. develop interactively user graphical interface.
4, according to the implementation method of the described interactive coronary artery virtual angioscope of claim 3, it is characterized in that, described ultrasound catheter returns the three-dimensional rebuilding method of removing the path: at first set up the perspective projection imaging model of CAG system two near normal angles, again according to the distance and the angle parameter of synchronous recording in angiographic procedure, obtain the geometric transformation of imaging system, adopt three-dimensional snake modelling technique then, snake directly is out of shape in the space, finishes conduit and returns the three-dimensional reconstruction of removing the path.
5, according to the implementation method of the described interactive coronary artery virtual angioscope of claim 4, it is characterized in that, the span of angle is 60 ° to 120 ° between two acquisition angles of described coronarogram picture, and only returns the starting point of removing the path at ultrasound catheter and take a pair of contrastographic picture.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100800020A CN101422352B (en) | 2008-12-10 | 2008-12-10 | Interactive coronary artery virtual angioscope implementation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100800020A CN101422352B (en) | 2008-12-10 | 2008-12-10 | Interactive coronary artery virtual angioscope implementation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101422352A true CN101422352A (en) | 2009-05-06 |
CN101422352B CN101422352B (en) | 2011-07-13 |
Family
ID=40613397
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100800020A Expired - Fee Related CN101422352B (en) | 2008-12-10 | 2008-12-10 | Interactive coronary artery virtual angioscope implementation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101422352B (en) |
Cited By (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102289814A (en) * | 2011-08-30 | 2011-12-21 | 北京理工大学 | Cardiac nuclear magnetic resonance image segmentation method |
CN102509267A (en) * | 2011-11-08 | 2012-06-20 | 华北电力大学(保定) | Retrospective off-line gating method for intravascular ultrasound image sequence |
CN102592274A (en) * | 2010-11-12 | 2012-07-18 | 通用电气公司 | Method for processing radiographic images for stenosis detection |
CN103247071A (en) * | 2013-03-29 | 2013-08-14 | 哈尔滨工业大学深圳研究生院 | Method and device for constructing three-dimensional blood vessel model |
CN105938513A (en) * | 2015-03-04 | 2016-09-14 | 三星电子株式会社 | Apparatus and method for providing reliability for computer aided diagnosis |
CN102592274B (en) * | 2010-11-12 | 2016-12-14 | 通用电气公司 | Processing procedure for the radiation image of stenosis detection |
CN106909770A (en) * | 2015-12-21 | 2017-06-30 | 佳能株式会社 | Medical image processing unit, its control method and storage medium |
CN107330236A (en) * | 2017-05-11 | 2017-11-07 | 青岛大学附属医院 | Virtual endoscope system with improvement roaming effect |
CN107468334A (en) * | 2017-08-01 | 2017-12-15 | 强联智创(北京)科技有限公司 | A kind of moulding aided design system of three-dimensional microtubular and design method |
CN107689072A (en) * | 2016-06-12 | 2018-02-13 | 中慧医学成像有限公司 | A kind of 3-D view imaging method and system |
CN107945169A (en) * | 2017-12-01 | 2018-04-20 | 中国人民解放军第三军医大学 | A kind of coronary artery image analysis methods and data structure |
CN108428210A (en) * | 2017-02-15 | 2018-08-21 | 浙江京新术派医疗科技有限公司 | Blood-vessel image method for reconstructing and reconstructing device |
CN108882854A (en) * | 2016-03-21 | 2018-11-23 | 华盛顿大学 | The virtual reality or augmented reality of 3D medical image visualize |
CN109087352A (en) * | 2018-08-16 | 2018-12-25 | 数坤(北京)网络科技有限公司 | A kind of heart coronary artery advantage type automatic distinguishing method |
CN109219810A (en) * | 2016-06-01 | 2019-01-15 | 微软技术许可有限责任公司 | Online visual angle for 3D component is searched for |
CN109427059A (en) * | 2017-08-18 | 2019-03-05 | 西门子医疗有限公司 | The plane visualization of anatomical structure |
CN109584148A (en) * | 2018-11-27 | 2019-04-05 | 重庆爱奇艺智能科技有限公司 | A kind of method and apparatus handling two-dimentional interface in VR equipment |
CN109598011A (en) * | 2017-08-31 | 2019-04-09 | 和承R&A有限公司 | The method and its recording medium of pipe layout are explained by using three-dimensional coordinate |
CN109800814A (en) * | 2019-01-25 | 2019-05-24 | 西南科技大学 | The Invariant feature point extracting method of curved measurement positioning |
CN110458853A (en) * | 2019-08-01 | 2019-11-15 | 北京灵医灵科技有限公司 | Ankle ligament separation method and separation system in a kind of medical image |
CN111009032A (en) * | 2019-12-04 | 2020-04-14 | 浙江理工大学 | Blood vessel three-dimensional reconstruction method based on improved epipolar line constraint matching |
CN111553979A (en) * | 2020-05-26 | 2020-08-18 | 广州狄卡视觉科技有限公司 | Operation auxiliary system and method based on medical image three-dimensional reconstruction |
CN112057167A (en) * | 2019-05-22 | 2020-12-11 | 深圳市德力凯医疗设备股份有限公司 | Ultrasonic navigation method and ultrasonic navigation equipment for vascular surgery |
CN112292076A (en) * | 2018-04-11 | 2021-01-29 | 弗劳恩霍夫应用研究促进协会 | Position determining apparatus for determining the position of an article within a tubular structure |
CN112529976A (en) * | 2020-11-26 | 2021-03-19 | 上海商汤智能科技有限公司 | Target display method and device, electronic equipment and storage medium |
CN112652052A (en) * | 2020-12-15 | 2021-04-13 | 山东大学 | Coronary artery three-dimensional reconstruction method and system based on blood vessel branch registration |
WO2021109122A1 (en) * | 2019-12-05 | 2021-06-10 | 苏州润迈德医疗科技有限公司 | Method and device for establishing mathematical model of blood vessel having stenotic lesion |
CN113409275A (en) * | 2021-06-22 | 2021-09-17 | 青岛海信医疗设备股份有限公司 | Method for determining thickness of transparent layer behind fetal neck based on ultrasonic image and related device |
CN114145719A (en) * | 2022-02-08 | 2022-03-08 | 天津恒宇医疗科技有限公司 | Method and system for three-dimensional fusion of dual-mode coronary vessel images |
-
2008
- 2008-12-10 CN CN2008100800020A patent/CN101422352B/en not_active Expired - Fee Related
Cited By (45)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102592274B (en) * | 2010-11-12 | 2016-12-14 | 通用电气公司 | Processing procedure for the radiation image of stenosis detection |
CN102592274A (en) * | 2010-11-12 | 2012-07-18 | 通用电气公司 | Method for processing radiographic images for stenosis detection |
CN102289814B (en) * | 2011-08-30 | 2012-12-26 | 北京理工大学 | Cardiac nuclear magnetic resonance image segmentation method |
CN102289814A (en) * | 2011-08-30 | 2011-12-21 | 北京理工大学 | Cardiac nuclear magnetic resonance image segmentation method |
CN102509267A (en) * | 2011-11-08 | 2012-06-20 | 华北电力大学(保定) | Retrospective off-line gating method for intravascular ultrasound image sequence |
CN102509267B (en) * | 2011-11-08 | 2013-11-06 | 华北电力大学(保定) | Retrospective off-line gating method for intravascular ultrasound image sequence |
CN103247071A (en) * | 2013-03-29 | 2013-08-14 | 哈尔滨工业大学深圳研究生院 | Method and device for constructing three-dimensional blood vessel model |
CN103247071B (en) * | 2013-03-29 | 2015-11-11 | 哈尔滨工业大学深圳研究生院 | A kind of structure three-dimensional blood vessel model method and apparatus |
CN105938513A (en) * | 2015-03-04 | 2016-09-14 | 三星电子株式会社 | Apparatus and method for providing reliability for computer aided diagnosis |
CN106909770A (en) * | 2015-12-21 | 2017-06-30 | 佳能株式会社 | Medical image processing unit, its control method and storage medium |
CN106909770B (en) * | 2015-12-21 | 2020-11-03 | 佳能株式会社 | Medical image processing apparatus, control method thereof, and storage medium |
CN108882854B (en) * | 2016-03-21 | 2022-05-24 | 华盛顿大学 | Virtual reality or augmented reality visualization of 3D medical images |
CN108882854A (en) * | 2016-03-21 | 2018-11-23 | 华盛顿大学 | The virtual reality or augmented reality of 3D medical image visualize |
US11468111B2 (en) | 2016-06-01 | 2022-10-11 | Microsoft Technology Licensing, Llc | Online perspective search for 3D components |
CN109219810A (en) * | 2016-06-01 | 2019-01-15 | 微软技术许可有限责任公司 | Online visual angle for 3D component is searched for |
CN107689072A (en) * | 2016-06-12 | 2018-02-13 | 中慧医学成像有限公司 | A kind of 3-D view imaging method and system |
CN108428210A (en) * | 2017-02-15 | 2018-08-21 | 浙江京新术派医疗科技有限公司 | Blood-vessel image method for reconstructing and reconstructing device |
CN108428210B (en) * | 2017-02-15 | 2021-11-16 | 浙江京新术派医疗科技有限公司 | Blood vessel image reconstruction method and reconstruction device |
CN107330236A (en) * | 2017-05-11 | 2017-11-07 | 青岛大学附属医院 | Virtual endoscope system with improvement roaming effect |
CN107468334B (en) * | 2017-08-01 | 2019-07-16 | 强联智创(北京)科技有限公司 | A kind of three-dimensional microtubular moulding aided design system and design method |
CN107468334A (en) * | 2017-08-01 | 2017-12-15 | 强联智创(北京)科技有限公司 | A kind of moulding aided design system of three-dimensional microtubular and design method |
CN109427059A (en) * | 2017-08-18 | 2019-03-05 | 西门子医疗有限公司 | The plane visualization of anatomical structure |
CN109427059B (en) * | 2017-08-18 | 2022-05-24 | 西门子医疗有限公司 | Planar visualization of anatomical structures |
CN109598011B (en) * | 2017-08-31 | 2023-08-22 | 和承R&A有限公司 | Method for interpreting tube layout by using three-dimensional coordinates and recording medium thereof |
CN109598011A (en) * | 2017-08-31 | 2019-04-09 | 和承R&A有限公司 | The method and its recording medium of pipe layout are explained by using three-dimensional coordinate |
CN107945169A (en) * | 2017-12-01 | 2018-04-20 | 中国人民解放军第三军医大学 | A kind of coronary artery image analysis methods and data structure |
CN107945169B (en) * | 2017-12-01 | 2022-02-15 | 中国人民解放军第三军医大学 | Coronary artery image analysis method |
CN112292076A (en) * | 2018-04-11 | 2021-01-29 | 弗劳恩霍夫应用研究促进协会 | Position determining apparatus for determining the position of an article within a tubular structure |
CN109087352A (en) * | 2018-08-16 | 2018-12-25 | 数坤(北京)网络科技有限公司 | A kind of heart coronary artery advantage type automatic distinguishing method |
CN109584148A (en) * | 2018-11-27 | 2019-04-05 | 重庆爱奇艺智能科技有限公司 | A kind of method and apparatus handling two-dimentional interface in VR equipment |
CN109800814A (en) * | 2019-01-25 | 2019-05-24 | 西南科技大学 | The Invariant feature point extracting method of curved measurement positioning |
CN109800814B (en) * | 2019-01-25 | 2022-08-09 | 西南科技大学 | Invariant characteristic quantity extraction method for blade curve measurement positioning |
CN112057167A (en) * | 2019-05-22 | 2020-12-11 | 深圳市德力凯医疗设备股份有限公司 | Ultrasonic navigation method and ultrasonic navigation equipment for vascular surgery |
CN110458853A (en) * | 2019-08-01 | 2019-11-15 | 北京灵医灵科技有限公司 | Ankle ligament separation method and separation system in a kind of medical image |
CN110458853B (en) * | 2019-08-01 | 2021-02-26 | 北京灵医灵科技有限公司 | Ankle ligament separation method and system in medical image |
CN111009032A (en) * | 2019-12-04 | 2020-04-14 | 浙江理工大学 | Blood vessel three-dimensional reconstruction method based on improved epipolar line constraint matching |
CN111009032B (en) * | 2019-12-04 | 2023-09-19 | 浙江理工大学 | Vascular three-dimensional reconstruction method based on improved epipolar line constraint matching |
WO2021109122A1 (en) * | 2019-12-05 | 2021-06-10 | 苏州润迈德医疗科技有限公司 | Method and device for establishing mathematical model of blood vessel having stenotic lesion |
CN111553979A (en) * | 2020-05-26 | 2020-08-18 | 广州狄卡视觉科技有限公司 | Operation auxiliary system and method based on medical image three-dimensional reconstruction |
CN111553979B (en) * | 2020-05-26 | 2023-12-26 | 广州雪利昂生物科技有限公司 | Operation auxiliary system and method based on three-dimensional reconstruction of medical image |
CN112529976A (en) * | 2020-11-26 | 2021-03-19 | 上海商汤智能科技有限公司 | Target display method and device, electronic equipment and storage medium |
CN112652052A (en) * | 2020-12-15 | 2021-04-13 | 山东大学 | Coronary artery three-dimensional reconstruction method and system based on blood vessel branch registration |
CN113409275A (en) * | 2021-06-22 | 2021-09-17 | 青岛海信医疗设备股份有限公司 | Method for determining thickness of transparent layer behind fetal neck based on ultrasonic image and related device |
CN113409275B (en) * | 2021-06-22 | 2022-07-01 | 青岛海信医疗设备股份有限公司 | Method for determining thickness of transparent layer behind fetal neck based on ultrasonic image and related device |
CN114145719A (en) * | 2022-02-08 | 2022-03-08 | 天津恒宇医疗科技有限公司 | Method and system for three-dimensional fusion of dual-mode coronary vessel images |
Also Published As
Publication number | Publication date |
---|---|
CN101422352B (en) | 2011-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101422352B (en) | Interactive coronary artery virtual angioscope implementation method | |
CN101283929B (en) | Rebuilding method of blood vessel three-dimensional model | |
US6148095A (en) | Apparatus and method for determining three-dimensional representations of tortuous vessels | |
Roelandt et al. | Three-dimensional reconstruction of intracoronary ultrasound images. Rationale, approaches, problems, and directions. | |
JP5474342B2 (en) | Anatomical modeling with 3-D images and surface mapping | |
EP1088515B1 (en) | Vascular reconstruction | |
CN100539947C (en) | Medical image demonstration and processing method, CT equipment, work station and program product | |
EP1030191A2 (en) | Semi-automated segmentation method for 3-dimensional ultrasound | |
US20100189337A1 (en) | Method for acquiring 3-dimensional images of coronary vessels, particularly of coronary veins | |
EP0884975A1 (en) | Catheter for obtaining three-dimensional reconstruction of a vascular lumen and wall, and method of use | |
Zheng et al. | Sequential reconstruction of vessel skeletons from X-ray coronary angiographic sequences | |
JP2006516440A (en) | 3D object reconstruction method | |
Di Mario et al. | Three dimensional reconstruction of cross sectional intracoronary ultrasound: clinical or research tool? | |
CN114145719B (en) | Method and system for three-dimensional fusion of dual-mode coronary vessel images | |
Weng et al. | Three-dimensional surface reconstruction using optical flow for medical imaging | |
CN116188677A (en) | Three-dimensional reconstruction method, system and device for vascular intervention operation area | |
CN112669449A (en) | CAG and IVUS accurate linkage analysis method and system based on 3D reconstruction technology | |
Weichert et al. | Virtual 3D IVUS vessel model for intravascular brachytherapy planning. I. 3D segmentation, reconstruction, and visualization of coronary artery architecture and orientation | |
Zheng et al. | Reconstruction of coronary vessels from intravascular ultrasound image sequences based on compensation of the in-plane motion | |
Lengyel et al. | Three-dimensional reconstruction and volume rendering of intravascular ultrasound slices imaged on a curved arterial path | |
Rotger et al. | Internal and external coronary vessel images registration | |
Garcia | Three-dimensional imaging for coronary interventions | |
Zheng | Three-dimensional reconstruction of vessels from intravascular ultrasound sequence and X-ray angiograms | |
Chandran et al. | Coronary arteries: imaging, reconstruction, and fluid dynamic analysis | |
US20240070855A1 (en) | Method and System of Calculating the Cross-section Area and Included Angle of Three-dimensional Blood Vessel Branch |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110713 Termination date: 20141210 |
|
EXPY | Termination of patent right or utility model |