CN106920282A - The edit methods and system of a kind of blood vessel mathematical model - Google Patents
The edit methods and system of a kind of blood vessel mathematical model Download PDFInfo
- Publication number
- CN106920282A CN106920282A CN201510991272.7A CN201510991272A CN106920282A CN 106920282 A CN106920282 A CN 106920282A CN 201510991272 A CN201510991272 A CN 201510991272A CN 106920282 A CN106920282 A CN 106920282A
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- point
- triangle
- mathematical model
- array
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/04—Indexing scheme for image data processing or generation, in general involving 3D image data
Landscapes
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Processing Or Creating Images (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Instructional Devices (AREA)
Abstract
The present invention provides the edit methods and system of a kind of blood vessel mathematical model, and the edit methods include:Obtain multiple triangle surfaces of blood vessel mathematical model;Analyze and obtain the adjacent surface information of each triangle surface;Receive user and the second point that select in screen coordinates at first point, and determine the index information of the triangle surface of first end point, the second end points and place in and the corresponding blood vessel mathematical model of second point respectively at first point;The index information and adjacent surface information of the triangle surface according to where first end point, the second end points determine the target blood section of multiple blood vessel annulus compositions;Midpoint according to each blood vessel annulus intermediate cam shape dough sheet determines the central point of each blood vessel annulus;The relative position editor target blood section of summit according to each blood vessel annulus intermediate cam shape dough sheet and respective central point.By the present invention, vascular lesion situation can be added in any position to existing blood vessel mathematical model, can greatly enrich the adoptable case mode storehouse in virtual operation field.
Description
Technical field
The present invention is on virtual operation modelling technique, specifically on a kind of edit methods of blood vessel mathematical model
And system.
Background technology
Virtual operation is to utilize medical image data, and virtual environment is set up in a computer using virtual reality technology, is made
Doctor can carry out surgical planning, training by the information in virtual environment, and operation is guided in actual operation.Pass
There is the shortcomings of risk is high, and patient's pain is big, postoperative effect is undesirable in the surgical procedure of system, surgical effect is by doctor
The influence of people's professional skill is very big.Medical personnel can be made to be immersed in virtual scene using augmented reality, led to
Cross vision, the sense of hearing, tactilely-perceptible and learn various operation practical operations, in experiencing and learning how to deal with clinical operation
Actual conditions.Expense and the time of training medical worker are which saved, makes the risk that unskilled persons are performed the operation
Property substantially reduce, it is uneven to the efficiency and quality that improve medical education and training and the development for improving medical surgery level
The present situation of weighing apparatus has special meaning.
Actual operation process is simulated using virtual reality technology, helps medical personnel to carry out the process of surgical skills training
In, the abundant degree for simulating used case mode is extremely important.First, from for the angle of trainer, rich
Be trained on rich case mode, can help trainee adapt to the operation technique flow under different pathological feature and
Skill, is that it carries out sufficient preparation when similar cases are run into during actual operation in the future.Second, take up a job as a doctor
For raw angle, possess the autonomous control power to vascular lesion position, prophylactically future can be likely to occur
Rare case is simulated in advance, while the training resource to new hand doctor can be enriched.3rd, from operating theater instruments factory
For business's angle, new operating theater instruments research and development can be carried out to any vascular lesion position and is tested, it is to avoid single
The deficiency of imitative body Model experiment.
Blood vessel mathematical model source employed in current virtual operation process mainly has following two:
1st, three-dimensional reconstruction
The method is based on the scan images such as the CT/MRI/DSA of actual case, and then is obtained using three-dimensional reconstruction
Three-dimensional blood vessel mathematical model.But by the method obtain vascular pattern be often it is coarse, be by manual intervention
Later stage treatment could be used in during virtual operation.The premise that the method obtains case mode is the figure for possessing the case
Picture, and grasp the post-processing method of the rough model obtained to three-dimensional reconstruction.
Case mode is obtained based on scan-image, three-dimensional reconstruction, it is necessary to obtain the scan image of the actual case first.
And the special case rarer for some pathologic structures, its patient data is often to be difficult to what is obtained.Especially examining
On the premise of considering patient privacy, the case data of some patients are typically protected and not disclosed.Therefore, this
The method for planting abundant lesion vesselses model library is just greatly limited in terms of data source.Additionally, three-dimensional reconstruction
The blood vessel mathematical model of acquisition generally requires to put into that substantial amounts of manpower goes to carry out manual treatment and adjustment can be only achieved void
Intend the use standard of operation, the manpower and materials cost of this work is equally also very high, and the cycle is more long.
2nd, three-dimensional software for editing.
The method is set up and is obtaining the base by the three-dimensional blood vessel mathematical model after three-dimensional reconstruction and later stage treatment
On plinth.By the three-dimensional such as 3DMax or Maya software for editing, edlin is entered to three-dimensional blood vessel mathematical model,
And then obtain different types of neuropathy model.The method can edit acquisition pathological characters on healthy blood vessel mathematical model
The abundant case mode of species.
Lesion vesselses model is obtained by three-dimensional software editor health vascular pattern, it is desirable to which operating personnel skillfully use
The three-dimensional such as 3DMax or Maya software for editing, the software person of needing to use of this type has to the form of blood vessel three-dimensional model
It is certain to understand, and operating method is complicated, left-hand seat is slower.
The content of the invention
The main purpose of the embodiment of the present invention is the edit methods and system for providing a kind of blood vessel mathematical model, to provide
A kind of operating procedure is simple, easy left-hand seat blood vessel mathematical model edit process.
To achieve these goals, the embodiment of the present invention provides a kind of edit methods of blood vessel mathematical model, described volume
The method of collecting includes:Obtain multiple triangle surfaces of blood vessel mathematical model;Analyze and obtain each triangle surface
Adjacent surface information;Receive and second point that user selects in screen coordinates at first point, and determine respectively at described first point
And the index letter of the triangle surface of the first end point, the second end points and place in the corresponding blood vessel mathematical model of second point
Breath;The index information and adjacent surface information of the triangle surface according to where the first end point, the second end points determine many
The target blood section of individual blood vessel annulus composition;Midpoint according to each blood vessel annulus intermediate cam shape dough sheet determine it is each described in
The central point of blood vessel annulus;Summit according to each blood vessel annulus intermediate cam shape dough sheet is relative with respective central point
Edit the target blood section in position.
In one embodiment, above-mentioned analysis and the adjacent surface information of each triangle surface is obtained, including:For
Each summit V on the blood vessel mathematical modeli, obtain all comprising the summit V on the blood vessel mathematical modeli's
Triangle surface, and its face index information is stored in and each summit ViCorresponding array AiIn, wherein, i is for just
Integer;To each triangle surface T on the blood vessel mathematical modeliEach summit Vj, in the summit VjIt is right
The array A for answeringjMiddle acquisition includes the triangle surface TiNext summit A(j+1) %3Triangle surface T 'iFace
Index information;By the triangle surface T 'iFace index information be stored in and the triangle surface TiCorresponding array Bi
In, the array BiIt is the triangle surface TiThe index information of the triangle surface for being abutted, wherein, j=0,1,2.
In one embodiment, it is above-mentioned determine respectively in described first point and the corresponding blood vessel mathematical model of second point
The index information of the triangle surface of end point, the second end points and place, including:Described first point and second point are existed
In coordinate transformation to three dimensional space coordinate in screen coordinate system, first three-dimensional coordinate and second point three-dimensional coordinate are generated;
Determine first three-dimensional coordinate and second point three-dimensional coordinate respectively with virtual reference coordinate constitute the first ray and
Second ray;Using the common factor of first ray and the blood vessel mathematical model intermediate cam shape dough sheet as the first end
Point, using the common factor of second ray and the blood vessel mathematical model intermediate cam shape dough sheet as second end points, and
The index information of the triangle surface where obtaining the first end point and the second end points.
In one embodiment, the above-mentioned common factor by first ray and the blood vessel mathematical model intermediate cam shape dough sheet
As the first end point, using the common factor of second ray and the blood vessel mathematical model intermediate cam shape dough sheet as institute
The second end points is stated, including:The blood vessel mathematical model is split by octree partition method, by the blood vessel number
Word model is divided into many sub-spaces, and the quantity of each subspace intermediate cam shape dough sheet is less than a threshold value;Described first
The common factor of the triangle surface in ray and subspace is the first end point, in second ray and subspace three
The common factor of angular dough sheet is second end points.
In one embodiment, the index letter of the above-mentioned triangle surface according to where the first end point, the second end points
Breath and adjacent surface information determine the target blood section of multiple blood vessel annulus compositions, including:Step a:First end point A and
Second terminal B constitutes vectorIt is determined that excessively described first end point A and perpendicular to the vectorPlane C, institute
Plane C is stated to intersect with the first triangle surface comprising the first end point A;Step b:By first triangle
The index information of shape dough sheet is stored in array 0, and marks first triangle surface;Step c:Obtain described
In the proximal surface of one triangle surface both sides, intersect with the plane C and the first labeled proximal surface, by described the
The index information of one proximal surface is stored in the array 0, and marks first proximal surface;Step d:Obtain described first
In the proximal surface of proximal surface, intersect with the plane C and the second not labeled proximal surface, by the index of second proximal surface
Information is stored in the array 0, and marks second proximal surface;Step e:Using second proximal surface as new
One proximal surface, circulation performs the step d, until the second proximal surface described in two overlaps;Step f:By in the array 0
The triangle surface composition blood vessel annulus A of preservation0;Step g:By in the proximal surface of the triangle surface in array n along institute
State vectorThe index information of the proximal surface of extension is stored in array n+1, and the triangle that will be preserved in the array n+1
Shape dough sheet composition blood vessel annulus An+1, n is nonnegative integer, and the initial value of n is 0;Step h:If array n+1
In comprising the triangle surface where second terminal B, then perform step i;If not comprising described in array n+1
Triangle surface where second terminal B, then add 1 by n, and returns to the execution step g;Step i:By institute
State the blood vessel annulus A of 0~array of array n+1 correspondence compositions0~blood vessel annulus An+1Constitute the target blood section.
In one embodiment, if the coordinate of first end point A is (xa, ya, za), the coordinate of the second terminal B is
(xb, yb, zb), then the equation of the plane C is expressed as:(x-xa)nx+(y-ya)ny+(z-za)nz=0,
Wherein, (nx, ny, nz) it is vectorDirection vector:
In one embodiment, the above-mentioned summit according to each blood vessel annulus intermediate cam shape dough sheet and respective central point
Relative position editor described in target blood section, including:According to summit P and central point C on each blood vessel annulusn
Position determine from summit P to central point CnVectorAccording to formula P=P+tPCn, t ∈ (0,1) make point P
AlongDirection move t*l sizes distance, make the blood vessel annulus AnTo central point CnDirection movement.
The embodiment of the present invention also provides a kind of editing system of blood vessel mathematical model, and the editing system includes:Gore
Piece acquiring unit, the multiple triangle surfaces for obtaining blood vessel mathematical model;Adjacent surface information acquisition unit, is used for
Analyze and obtain the adjacent surface information of each triangle surface;Index information determining unit, for receiving user in screen
First point and the second point selected in curtain coordinate, and determine and the corresponding blood vessel numeral mould of second point respectively at described first point
The index information of the triangle surface of first end point, the second end points and place in type;Target blood section determining unit,
Index information and adjacent surface information for the triangle surface according to where the first end point, the second end points determine many
The target blood section of individual blood vessel annulus composition;Central point determining unit, for according to each blood vessel annulus intermediate cam shape
The midpoint of dough sheet determines the central point of each blood vessel annulus;Target blood section edit cell, for according to each blood
The summit of pipe annulus intermediate cam shape dough sheet and target blood section described in the relative position editor of respective central point.
In one embodiment, above-mentioned adjacent surface information acquisition unit specifically for:For on the blood vessel mathematical model
Each summit Vi, obtain all comprising the summit V on the blood vessel mathematical modeliTriangle surface, and by its
Face index information is stored in and each summit ViCorresponding array AiIn, wherein, i is positive integer;To the blood vessel number
Each triangle surface T on word modeliEach summit Vj, in the summit VjCorresponding array AjMiddle acquisition
Comprising the triangle surface TiNext summit A(j+1) %3Triangle surface T 'iFace index information;Will be described
Triangle surface T 'iFace index information be stored in and the triangle surface TiCorresponding array BiIn, the array BiFor
The triangle surface TiThe index information of the triangle surface for being abutted, wherein, j=0,1,2.
In one embodiment, above-mentioned index information determining unit specifically for:User is received to select in screen coordinates
First point and second point, coordinate transformation to the three dimensions by described first point and second point in screen coordinate system sits
In mark, first three-dimensional coordinate and second point three-dimensional coordinate are generated;Determine first three-dimensional coordinate and second point three
The first ray and the second ray that dimension coordinate is constituted with virtual reference coordinate respectively;By first ray and the blood vessel
The common factor of mathematical model intermediate cam shape dough sheet as the first end point, by second ray and blood vessel numeral mould
The common factor of type intermediate cam shape dough sheet obtains the first end point and where the second end points three as second end points
The index information of angular dough sheet.
In one embodiment, the above-mentioned common factor by first ray and the blood vessel mathematical model intermediate cam shape dough sheet
As the first end point, using the common factor of second ray and the blood vessel mathematical model intermediate cam shape dough sheet as institute
The second end points is stated, including:The blood vessel mathematical model is split by octree partition method, by the blood vessel number
Word model is divided into many sub-spaces, and the quantity of each subspace intermediate cam shape dough sheet is less than a threshold value;Described first
The common factor of the triangle surface in ray and subspace is the first end point, in second ray and subspace three
The common factor of angular dough sheet is second end points.
In one embodiment, above-mentioned target blood section determining unit is specifically for performing following steps:Step a:The
End point A and the second terminal B constitute vectorIt is determined that excessively described first end point A and perpendicular to the vector's
Plane C, the plane C intersect with the first triangle surface comprising the first end point A;Step b:Will be described
The index information of the first triangle surface is stored in array 0, and marks first triangle surface;Step c:Obtain
Take in the proximal surface of the first triangle surface both sides, intersect with the plane C and the first not labeled proximal surface,
The index information of first proximal surface is stored in the array 0, and marks first proximal surface;Step d:Obtain
In the proximal surface of first proximal surface, intersect with the plane C and the second not labeled proximal surface, it is adjacent by described second
The index information in face is stored in the array 0, and marks second proximal surface;Step e:Second proximal surface is made
It is the first new proximal surface, circulation performs the step d, until the second proximal surface described in two overlaps;Step f:Will be described
The triangle surface composition blood vessel annulus A preserved in array 00;Step g:By the neighbour of the triangle surface in array n
Along the vector in faceThe index information of the proximal surface of extension is stored in array n+1, and will be protected in the array n+1
The triangle surface composition blood vessel annulus A for depositingn+1, n is nonnegative integer, and the initial value of n is 0;Step h:If
Comprising the triangle surface where second terminal B in array n+1, then step i is performed;If in array n+1 not
Comprising the triangle surface where second terminal B, then n is added 1, and return to the execution step g;Step
i:By the blood vessel annulus A of 0~array of array n+1 correspondence compositions0~blood vessel annulus AnConstitute the target blood section.
In one embodiment, if the coordinate of first end point A is (xa, ya, za), the coordinate of the second terminal B is
(xb, yb, zb), then the equation of the plane C is expressed as:(x-xa)nx+(y-ya)ny+(z-za)nz=0,
Wherein, (nx, ny, nz) it is vectorDirection vector:
In one embodiment, above-mentioned target blood section edit cell specifically for:According on each blood vessel annulus
Summit P and central point CnPosition determine from summit P to central point CnVectorAccording to formula P=P+
tPCn, t ∈ (0,1) make point P alongDirection move t*l sizes distance, make the blood vessel annulus AnTo center
Point CnDirection movement.
The beneficial effect of the embodiment of the present invention is, by the embodiment of the present invention, existing blood vessel mathematical model can be existed
Vascular lesion situation is added in any position, and can control the extent of disease and lesion degree of diseased region.This hair
The edit methods and system of the blood vessel mathematical model of bright embodiment enter edlin to blood vessel mathematical model, can be greatly rich
The rich adoptable case mode storehouse in virtual operation field.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, embodiment will be described below
Needed for the accompanying drawing to be used be briefly described, it should be apparent that, drawings in the following description are only of the invention
Some embodiments, for those of ordinary skill in the art, without having to pay creative labor, may be used also
Other accompanying drawings are obtained with according to these accompanying drawings.
Fig. 1 is the flow chart of the edit methods of the blood vessel mathematical model according to the embodiment of the present invention;
Fig. 2 is the blood vessel grid model schematic diagram according to the embodiment of the present invention;
Fig. 3 is the triangle surface difference proximal surface quantity schematic diagram according to the embodiment of the present invention;
Fig. 4 A and Fig. 4 B are to ask friendship to choose triangle surface with blood vessel mathematical model according to the ray of the embodiment of the present invention
Process schematic;
Fig. 5 is the triangle surface process schematic that target blood section is found according to the iteration of the embodiment of the present invention;
Fig. 6 is to ask for blood vessel circle ring center point process schematic diagram according to the embodiment of the present invention;
Fig. 7 is the blood vessel annulus editing process schematic diagram of the target blood section according to the embodiment of the present invention;
Fig. 8 is the target blood narrow edited result schematic diagram of section according to the embodiment of the present invention;
Fig. 9 is the hemangioma editor's schematic diagram according to another embodiment of the present invention;
Figure 10 is the structural representation of the editor of the blood vessel mathematical model according to the embodiment of the present invention.
Specific embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Site preparation is described, it is clear that described embodiment is only a part of embodiment of the invention, rather than whole embodiments.
Based on the embodiment in the present invention, what those of ordinary skill in the art were obtained under the premise of creative work is not made
Every other embodiment, belongs to the scope of protection of the invention.
The embodiment of the present invention provides the edit methods and system of a kind of blood vessel mathematical model.Below in conjunction with accompanying drawing to the present invention
It is described in detail.
The embodiment of the present invention provides a kind of edit methods of blood vessel mathematical model, as shown in figure 1, the edit methods are main
Including following steps:
Step S101:Obtain multiple triangle surfaces of blood vessel mathematical model;
Step S102:Analyze and obtain the adjacent surface information of each triangle surface;
Step S103:Receive and second point that user selects in screen coordinates at first point, and determine respectively at first point
And the index letter of the triangle surface of the first end point, the second end points and place in the corresponding blood vessel mathematical model of second point
Breath;
Step S104:Index information and the adjacent surface letter of the triangle surface according to where first end point, the second end points
Breath determines the target blood section of multiple blood vessel annulus compositions;
Step S105:Midpoint according to each blood vessel annulus intermediate cam shape dough sheet determines the central point of each blood vessel annulus;
Step S106:Compiled with the relative position of respective central point on summit according to each blood vessel annulus intermediate cam shape dough sheet
Collect target blood section.
By above-mentioned steps S101 to step S106, the edit methods of the blood vessel mathematical model of the embodiment of the present invention can
Vascular lesion situation is added in any position to existing blood vessel mathematical model, and the lesion of diseased region can be controlled
Scope and lesion degree.The operating procedure of the edit methods of the blood vessel mathematical model of the embodiment of the present invention is simple, easily
Left-hand seat.The edit methods of the blood vessel mathematical model of the embodiment of the present invention enter edlin to blood vessel mathematical model, can be very big
The adoptable case mode storehouse in abundant virtual operation field.
The edit methods of the blood vessel mathematical model of the embodiment of the present invention are described in detail below in conjunction with specific example.
Above-mentioned steps S101, obtains multiple triangle surfaces of blood vessel mathematical model.
Blood vessel mathematical model (grid model) is made up of small triangle surface, as shown in Fig. 2 each gore
Piece has 1 to 3 proximal surface.As shown in figure 3, when model non-closed and triangle surface are located at edge of opening, there is 1
Individual (white triangles shape in figure) or 2 proximal surfaces (gray triangles in figure), otherwise have 3 proximal surfaces (black in figure
Color triangle).Model file employed in virtual operation is generally OBJ or STL forms, and this document form is only carried
Necessary vertex index information when vertex information and triangle surface for constituting vascular pattern are assembled, does not provide every
The proximal surface information of individual triangle surface.Therefore, it is the proximal surface information of the triangle surface of i to need computation index first.
Step S102, analyzes and obtains the adjacent surface information of each triangle surface.
The proximal surface information of triangle is calculated, following two steps are broadly divided into:
1st, to each summit V of composition blood vessel mathematical modeli, find out the rope of all triangle surfaces comprising the summit
Draw, place it in the array A with each vertex correspondenceiIn.
2nd, the arbitrary triangle dough sheet T on blood vessel Array Modeli, to each summit on the triangle surface
Vj(j=0,1,2), in its corresponding array AjMiddle searching includes the next summit A of the triangle surface(j+1) %3Triangle
Index, and the index is put into TiCorresponding array BiIn.Last array BiIn result be exactly triangle surface TiInstitute
The index of adjacent triangle surface.
Step S103, receives and second point that user selects in screen coordinates at first point, and determine respectively at first point
And the index of the triangle surface of the first end point, the second end points and place in the corresponding blood vessel mathematical model of second point
Information.
User can choose two positions by way of mouse is clicked on pending blood vessel Array Model:End points
A and terminal B.Vessel segment between terminal A and terminal B is the target for entering edlin in subsequent processes
Vessel segment.
The visualization of three-dimensional blood vessel mathematical model is the display screen that the geometric data of three dimensions is transformed to two-dimensional space
On curtain.Object coordinates system, world coordinate system, camera coordinates system, cutting coordinate are experienced from the mapping of 2 d-to-3 d
A series of matrixing of coordinate systems such as system and screen coordinate system.User's vascular pattern for seeing on the screen of two dimension is
The result of the coordinate transform of screen coordinate system is tied to from object coordinates.
When user carries out mouse on two-dimensional screen to be clicked, selection is point that screen coordinate is fastened, it is therefore desirable to really
The point for determining user's selection is located on which triangle surface of vascular pattern actually.Assuming that the point coordinates that user clicks
It is P (x, y), by reverse coordinate transform, can be by the point P in the point transformation to world coordinate systemw(xw, yw, zw) on.
However, only determining that the triangle surface where it is very difficult by the position of the point.According to computer graphical
Learn basis and understand that each pixel on screen actually correspond to the virtual reference coordinate from " virtual camera "
The ray to blood vessel mathematical model for sending.Therefore, as shown in fig. 4 a and fig. 4b, according to point Pw(xw, yw, zw)
The virtual reference coordinate P of " virtual camera " in world coordinate systemwc(xwc, ywc, zwc), one can be obtained and penetrated
Line, calculates this ray and intersects with which triangle surface, so as to obtain the index of intersecting triangles dough sheet.
The triangle surface quantity for constituting blood vessel mathematical model is often very huge, using ray and each triangle
It is quite time-consuming that shape dough sheet iteration seeks the method for friendship, and with the increase of triangle surface quantity, its time loss
Linear growth.In order to improve the speed of this intersecting detection, the triangle surface for constituting blood vessel mathematical model can be entered
Row Octree space is split, i.e., the three dimensions for enclosing blood vessel mathematical model is divided into eight parallel with coordinate plane
Every sub-spaces, are then recursively continued segmentation by sub-spaces, until the depth of Octree reaches designated value or every
Number of triangles in sub-spaces reaches certain designated value.Wherein, the designated value of the depth of Octree is general 5~10
Between, depending on concrete numerical value is according to model complexity;The designated value of the number of triangles in per sub-spaces is general
Between 10~20, equally floated according to model complexity difference, the present invention is not limited thereto.Set up after Octree
Carrying out ray again can quickly exclude the substantial amounts of triangle surface that can not possibly intersect with the friendship detection of asking of triangle surface
Subset, so as to reach the purpose of quick calculating.
Step S104, index information and the adjacent surface letter of the triangle surface according to where first end point, the second end points
Breath determines the target blood section of multiple blood vessel annulus compositions.
1st, all triangles in target blood section are found.
As shown in figure 5, the vessel segment comprising terminal A and terminal B that above-mentioned steps S103 determines is intercepted out,
Cut open, tiling is in one plane.The terminal A and terminal B of the whole story position specified in user are where it
Triangle surface in mark out come.More accurately, terminal A and terminal B can be respectively positioned in correspondence triangle
The center of dough sheet.One vector can be obtained by terminal A to terminal BThere will necessarily be a plane C and cross end points
A and perpendicular to vectorAssuming that A point coordinates is (xa, ya, za), B point coordinates is (xb, yb, zb), then this is put down
The equation of face C is represented by:
(x-xa)nx+(y-ya)ny+(z-za)nz=0
Wherein, (nx, ny, nz) it is vectorDirection vector:
Plane C must intersect with multiple triangle surfaces of blood vessel, and the triangle comprising A points must be wherein it
One.The index of the triangle surface where A is put into array 0 first, and is made marks to the triangle surface,
Detection is intersected in expression excessively to it.Then the proximal surface of A points is obtained, and it is asked into friendship with plane C, if
The intersecting and face is not labeled, then the index in the face should be equally put into array 0, and is marked.As schemed
Shown in 5, the process is launched to the left and right sides both direction of A simultaneously, when both direction is met (i.e. without new
When face index adds array 0), circulation terminates.Now, the triangle surface for being preserved in array 0 constitutes blood vessel
One blood vessel annulus, the blood vessel annulus includes all triangle surfaces and face C phases on point A, and the blood vessel annulus
Hand over.
Then, based on the triangle surface in array 0, ring by ring ground is along vectorDirection finding three
Angular dough sheet, untill in the blood vessel annulus for finding comprising the triangle surface where B points.
Specifically, to the triangle surface in array 0, the blood vessel annulus of the proximal surface composition that can be searched out has two, one
Individual top (the vector in plane C) direction, another is in the lower section of plane C.The blood above C will be located at
Whether the index of all triangle surfaces in pipe annulus is put into array 1, and detected in array 1 comprising B points institute
Triangle surface, if being target blood section comprising, then all tri patch in array 0 and array 1
Triangle surface;If do not included, then just continue on vector based on the triangle surface in array 1
Direction finding, untill in array i comprising the triangle surface where B points.Noteworthy point is that, often
It is secondary when finding a triangle and adding in corresponding array, it is required for being marked the triangle surface.By this
The method of kind, when beginning look for proximal surface from array 1, can just reject on a direction as long as being detected to the mark
Ring, remaining ring must be along vectorExtension.That is, the walking direction of blood vessel annulus and plane C
Only carried out once in array 0.
When all triangle surfaces are all found, array 0 is obtained, array 1 ..., array i, altogether i index
Array.Triangle surface in each array separately constitutes a blood vessel annulus, and these blood vessel annulus have connected into target
Vessel segment.
Step S105, the midpoint according to each blood vessel annulus intermediate cam shape dough sheet determines the central point of each blood vessel annulus.
As shown in fig. 6, to each blood vessel annulus, each gore for constituting the blood vessel annulus can be calculated
The midpoint (shown in figure orbicular spot) of piece.These midpoints are averaging can obtain the central point C of the blood vessel annulusi.Often
Individual blood vessel annulus has such a central point, and i central point is coupled together, you can obtain target blood section
Center line.
Compiled with the relative position of respective central point on step S106, the summit according to each blood vessel annulus intermediate cam shape dough sheet
Collect target blood section.
The i center position C of blood vessel annulus is obtained in above-mentioned steps S105i.To triangle on each blood vessel annulus
The summit P of dough sheet, is present from summit P to central point CiVectorIfLength be l, as shown in fig. 7,
According to formula:P=P+tPCi, t ∈ (0,1) allow point P alongDirection move t*l sizes distance, you can allow
The direction movement at blood vessel annulus Shang Xiang centers, so as to the effect of hemadostewnosis can be formed, (a) in such as Fig. 8, (b),
Shown in (c) three groups of comparison diagrams.
Due to sharing some summits between two adjacent blood vessel annulus, same summit is moved twice can cause
The result of mistake, accordingly, it would be desirable to opposite vertexes P is marked.So before being moved to each summit, detection
The mark, if marked, need not move again, otherwise move.
And it is possible to the contraction of blood vessel is seemed smoother by parameter t, it is not lofty.Near marginal position (mesh
Mark vessel segment whole story position) blood vessel annulus can use less t values, with away from marginal position, t values
Increase.
From above content, the edit methods of the blood vessel mathematical model of the embodiment of the present invention can be in any position of blood vessel
Put the generation different lesion of stenosis, and narrow positions and healthy part transition relative smooth, it is not lofty.
In another embodiment, the present invention can be equally used for selection and the editor of blood vessel other diseased regions.With blood vessel
As a example by knurl, as shown in figure 9, starting point A is selected in the joint of hemangioma and blood vessel, along vertical with vessel directions
And point to the direction inside hemangioma any point B is taken in hemangioma, then can determine and a plane, the plane
Hemangioma is entirely cut along vascular wall, all triangle surfaces that acquisition is intersected with the plane, and along sensing blood
The direction finding proximal surface of tuberculation, you can obtain all angiomatous triangle surfaces of composition.According to all three for choosing
Angular dough sheet is calculated after angiomatous central point, can be by controlling the relative position on each summit and central point in hemangioma
Put, reach the purpose zoomed in and out to angiomatous radius.It can be seen that, by blood vessel provided in an embodiment of the present invention numeral
The edit methods of model, can change the lesion degree (angiomatous size) of the hemangioma lesion for having existed.
The embodiment of the present invention also provides a kind of editing system of blood vessel mathematical model, as shown in Figure 10, the editing system
Mainly include:Triangle surface acquiring unit 1, adjacent surface information acquisition unit 2, index information determining unit 3,
Target blood section determining unit 4, central point determining unit 5 and target blood section edit cell 6 etc..
Wherein, triangle surface acquiring unit 1, the multiple triangle surfaces for obtaining blood vessel mathematical model;It is adjacent
Face information acquisition unit 2, for analyzing and obtains the adjacent surface information of each triangle surface;Index information determining unit
3, for receiving and second point that user selects in screen coordinates at first point, and determine and second point respectively at first point
The index information of the triangle surface of first end point, the second end points and place in corresponding blood vessel mathematical model;Target
Vessel segment determining unit 4, index information and neighbour for the triangle surface according to where first end point, the second end points
Junction information determines the target blood section of multiple blood vessel annulus compositions;Central point determining unit 5, for according to each blood vessel
The midpoint of annulus intermediate cam shape dough sheet determines the central point of each blood vessel annulus;Target blood section edit cell 6, for root
According to the relative position editor target blood section of summit and the respective central point of each blood vessel annulus intermediate cam shape dough sheet.
By above-mentioned each part, the editing system of the blood vessel mathematical model of the embodiment of the present invention can be to existing blood
Pipe mathematical model adds vascular lesion situation in any position, and can control the extent of disease and disease of diseased region
Change degree.The editing system of the blood vessel mathematical model of the embodiment of the present invention enters edlin, Ke Yiji to blood vessel mathematical model
The big adoptable case mode storehouse in abundant virtual operation field.
Each part in editing system below in conjunction with specific example to the blood vessel mathematical model of the embodiment of the present invention
It is described in detail.
Above-mentioned triangle surface acquiring unit 1, the multiple triangle surfaces for obtaining blood vessel mathematical model.
Blood vessel mathematical model (grid model) is made up of small triangle surface, as shown in Fig. 2 each gore
Piece has 1 to 3 proximal surface.As shown in figure 3, when model non-closed and triangle surface are located at edge of opening, there is 1
Individual (white triangles shape in figure) or 2 proximal surfaces (gray triangles in figure), otherwise have 3 proximal surfaces (black in figure
Color triangle).Model file employed in virtual operation is generally OBJ or STL forms, and this document form is only carried
Necessary vertex index information when vertex information and triangle surface for constituting vascular pattern are assembled, does not provide every
The proximal surface information of individual triangle surface.Therefore, it is the proximal surface information of the triangle surface of i to need computation index first.
Adjacent surface information acquisition unit 2, for analyzing and obtains the adjacent surface information of each triangle surface.
The proximal surface information of triangle is calculated, following two steps are broadly divided into:
1st, to each summit V of composition blood vessel mathematical modeli, find out the rope of all triangle surfaces comprising the summit
Draw, place it in the array A with each vertex correspondenceiIn.
2nd, the arbitrary triangle dough sheet T on blood vessel Array Modeli, to each summit on the triangle surface
Vj(j=0,1,2), in its corresponding array AjMiddle searching includes the next summit A of the triangle surface(j+1) %3Triangle
Index, and the index is put into TiCorresponding array BiIn.Last array BiIn result be exactly triangle surface TiInstitute
The index of adjacent triangle surface.
Index information determining unit 3, for receiving and second point that user selects in screen coordinates at first point, and divides
Do not determine the triangle of the first end point, the second end points and place in and the corresponding blood vessel mathematical model of second point at first point
The index information of dough sheet.
User can choose two positions by way of mouse is clicked on pending blood vessel Array Model:End points
A and terminal B.Vessel segment between terminal A and terminal B is the target for entering edlin in subsequent processes
Vessel segment.
The visualization of three-dimensional blood vessel mathematical model is the display screen that the geometric data of three dimensions is transformed to two-dimensional space
On curtain.Object coordinates system, world coordinate system, camera coordinates system, cutting coordinate are experienced from the mapping of 2 d-to-3 d
A series of matrixing of coordinate systems such as system and screen coordinate system.User's vascular pattern for seeing on the screen of two dimension is
The result of the coordinate transform of screen coordinate system is tied to from object coordinates.
When user carries out mouse on two-dimensional screen to be clicked, selection is point that screen coordinate is fastened, it is therefore desirable to really
The point for determining user's selection is located on which triangle surface of vascular pattern actually.Assuming that the point coordinates that user clicks
It is P (x, y), by reverse coordinate transform, can be by the point P in the point transformation to world coordinate systemw(xw, yw, zw) on.
However, only determining that the triangle surface where it is very difficult by the position of the point.According to computer graphical
Learn basis and understand that each pixel on screen actually correspond to the virtual reference coordinate from " virtual camera "
The ray to blood vessel mathematical model for sending.Therefore, as shown in fig. 4 a and fig. 4b, according to point Pw(xw, yw, zw)
The virtual reference coordinate P of " virtual camera " in world coordinate systemwc(xwc, ywc, zwc), one can be obtained and penetrated
Line, calculates this ray and intersects with which triangle surface, so as to obtain the index of intersecting triangles dough sheet.
The triangle surface quantity for constituting blood vessel mathematical model is often very huge, using ray and each triangle
It is quite time-consuming that shape dough sheet iteration seeks the method for friendship, and with the increase of triangle surface quantity, its time loss
Linear growth.In order to improve the speed of this intersecting detection, the triangle surface for constituting blood vessel mathematical model can be entered
Row Octree space is split, i.e., the three dimensions for enclosing blood vessel mathematical model is divided into eight parallel with coordinate plane
Every sub-spaces, are then recursively continued segmentation by sub-spaces, until the depth of Octree reaches designated value or every
Number of triangles in sub-spaces reaches certain designated value.Wherein, the designated value of the depth of Octree is general 5~10
Between, depending on concrete numerical value is according to model complexity;The designated value of the number of triangles in per sub-spaces is general
Between 10~20, equally floated according to model complexity difference, the present invention is not limited thereto.Set up after Octree
Carrying out ray again can quickly exclude the substantial amounts of triangle surface that can not possibly intersect with the friendship detection of asking of triangle surface
Subset, so as to reach the purpose of quick calculating.
Target blood section determining unit 4, for the index of the triangle surface according to where first end point, the second end points
Information and adjacent surface information determine the target blood section of multiple blood vessel annulus compositions.
1st, all triangles in target blood section are found.
As shown in figure 5, the vessel segment comprising terminal A and terminal B that above-mentioned index information determining unit 3 is determined
Interception out, is cut open, and tiling is in one plane.The terminal A and terminal B of the whole story position specified in user
Marked out in triangle surface where it and.More accurately, it is right to be respectively positioned in terminal A and terminal B
Answer the center of triangle surface.One vector can be obtained by terminal A to terminal BThere will necessarily be a plane
C crosses terminal A and perpendicular to vectorAssuming that A point coordinates is (xa, ya, za), B point coordinates is (Xb, yb, zb),
The equation of so plane C is represented by:
(x-xa)nx+(y-ya)ny+(z-za)nz=0
Wherein, (nx, ny, nz) it is vectorDirection vector:
Plane C must intersect with multiple triangle surfaces of blood vessel, and the triangle comprising A points must be wherein it
One.The index of the triangle surface where A is put into array 0 first, and is made marks to the triangle surface,
Detection is intersected in expression excessively to it.Then the proximal surface of A points is obtained, and it is asked into friendship with plane C, if
The intersecting and face is not labeled, then the index in the face should be equally put into array 0, and is marked.If
Shown in Fig. 5, the process is launched to the left and right sides both direction of A simultaneously, (does not have when both direction is met
When new face index adds array 0), circulation terminates.Now, the triangle surface for being preserved in array 0 constitutes blood
One blood vessel annulus of pipe, the blood vessel annulus includes all triangle surfaces on point A, and the blood vessel annulus and face C
It is intersecting.
Then, based on the triangle surface in array 0, ring by ring ground is along vectorDirection finding three
Angular dough sheet, untill in the blood vessel annulus for finding comprising the triangle surface where B points.
Specifically, to the triangle surface in array 0, the blood vessel annulus of the proximal surface composition that can be searched out has two, one
Individual top (the vector in plane C) direction, another is in the lower section of plane C.The blood above C will be located at
Whether the index of all triangle surfaces in pipe annulus is put into array 1, and detected in array 1 comprising B points institute
Triangle surface, if being target blood section comprising, then all tri patch in array 0 and array 1
Triangle surface;If do not included, then just continue on vector based on the triangle surface in array 1
Direction finding, untill in array i comprising the triangle surface where B points.Noteworthy point is that, often
It is secondary when finding a triangle and adding in corresponding array, it is required for being marked the triangle surface.By this
The method of kind, when beginning look for proximal surface from array 1, can just reject on a direction as long as being detected to the mark
Ring, remaining ring must be along vectorExtension.That is, the walking direction of blood vessel annulus and plane C
Only carried out once in array 0.
When all triangle surfaces are all found, array 0 is obtained, array 1 ..., array i, altogether i index
Array.Triangle surface in each array separately constitutes a blood vessel annulus, and these blood vessel annulus have connected into target
Vessel segment.
Central point determining unit 5, for determining each blood vessel annulus according to the midpoint of each blood vessel annulus intermediate cam shape dough sheet
Central point.
As shown in fig. 6, to each blood vessel annulus, each gore for constituting the blood vessel annulus can be calculated
The midpoint (shown in figure orbicular spot) of piece.These midpoints are averaging can obtain the central point C of the blood vessel annulusi.Often
Individual blood vessel annulus has such a central point, and i central point is coupled together, you can obtain target blood section
Center line.
Target blood section edit cell 6, for summit and respective center according to each blood vessel annulus intermediate cam shape dough sheet
The relative position editor target blood section of point.
The i center position C of blood vessel annulus is obtained in above-mentioned central point determining unit 5i.To on each blood vessel annulus
The summit P of triangle surface, is present from summit P to central point CiVectorIfLength be l, such as Fig. 7
It is shown, according to formula:P=P+tPCi, t ∈ (0,1) allow point P alongDirection move t*l sizes distance,
The direction at blood vessel annulus Shang Xiang centers can be allowed to move, so as to the effect of hemadostewnosis can be formed, (a) in such as Fig. 8,
Shown in three groups of comparison diagrams of (b), (c).
Due to sharing some summits between two adjacent blood vessel annulus, same summit is moved twice can cause
The result of mistake, accordingly, it would be desirable to opposite vertexes P is marked.So before being moved to each summit, detection
The mark, if marked, need not move again, otherwise move.
And it is possible to the contraction of blood vessel is seemed smoother by parameter t, it is not lofty.Near marginal position (mesh
Mark vessel segment whole story position) blood vessel annulus can use less t values, with away from marginal position, t values
Increase.
From above content, the editing system of the blood vessel mathematical model of the embodiment of the present invention can be in any position of blood vessel
Put the generation different lesion of stenosis, and narrow positions and healthy part transition relative smooth, it is not lofty.
In another embodiment, the present invention can be equally used for selection and the editor of blood vessel other diseased regions.With blood vessel
As a example by knurl, as shown in figure 9, starting point A is selected in the joint of hemangioma and blood vessel, along vertical with vessel directions
And point to the direction inside hemangioma any point B is taken in hemangioma, then can determine and a plane, the plane
Hemangioma is entirely cut along vascular wall, all triangle surfaces that acquisition is intersected with the plane, and along sensing blood
The direction finding proximal surface of tuberculation, you can obtain all angiomatous triangle surfaces of composition.According to all three for choosing
Angular dough sheet is calculated after angiomatous central point, can be by controlling the relative position on each summit and central point in hemangioma
Put, reach the purpose zoomed in and out to angiomatous radius.It can be seen that, by blood vessel provided in an embodiment of the present invention numeral
The editing system of model, can change the lesion degree (angiomatous size) of the hemangioma lesion for having existed.
One of ordinary skill in the art will appreciate that all or part of step in realizing above-described embodiment method can lead to
Program is crossed to instruct the hardware of correlation to complete, the program can be stored in a computer read/write memory medium, than
Such as ROM/RAM, magnetic disc, CD.
Particular embodiments described above, has been carried out further in detail to the purpose of the present invention, technical scheme and beneficial effect
Describe in detail bright, should be understood that and the foregoing is only specific embodiment of the invention, be not used to limit this hair
Bright protection domain, all any modification, equivalent substitution and improvements within the spirit and principles in the present invention, done etc.,
Should be included within the scope of the present invention.
Claims (14)
1. a kind of edit methods of blood vessel mathematical model, it is characterised in that described edit methods include:
Obtain multiple triangle surfaces of blood vessel mathematical model;
Analyze and obtain the adjacent surface information of each triangle surface;
Receive and second point that user selects in screen coordinates at first point, and determine and second point respectively at described first point
The index information of the triangle surface of first end point, the second end points and place in corresponding blood vessel mathematical model;
The index information and adjacent surface information of the triangle surface according to where the first end point, the second end points determine many
The target blood section of individual blood vessel annulus composition;
Midpoint according to each blood vessel annulus intermediate cam shape dough sheet determines the central point of each blood vessel annulus;
Described in the relative position editor of summit according to each blood vessel annulus intermediate cam shape dough sheet and respective central point
Target blood section.
2. edit methods of blood vessel mathematical model according to claim 1, it is characterised in that analyze and obtain
The adjacent surface information of each triangle surface, including:
For each summit V on the blood vessel mathematical modeli, obtain all comprising described on the blood vessel mathematical model
Summit ViTriangle surface, and its face index information is stored in and each summit ViCorresponding array AiIn, wherein,
I is positive integer;
To each triangle surface T on the blood vessel mathematical modeliEach summit Vj, in the summit VjIt is right
The array A for answeringjMiddle acquisition includes the triangle surface TiNext summit A(j+1) %3Triangle surface T 'iFace
Index information;
By the triangle surface T 'iFace index information be stored in and the triangle surface TiCorresponding array BiIn, institute
State array BiIt is the triangle surface TiThe index information of the triangle surface for being abutted, wherein, j=0,1,2.
3. edit methods of blood vessel mathematical model according to claim 2, it is characterised in that determine institute respectively
State the triangle surface of the first end point, the second end points and place in and the corresponding blood vessel mathematical model of second point at first point
Index information, including:
In coordinate transformation to three dimensional space coordinate by described first point and second point in screen coordinate system, generation first
Point three-dimensional coordinate and second point three-dimensional coordinate;
Determine that first three-dimensional coordinate and second point three-dimensional coordinate first are penetrated with what virtual reference coordinate was constituted respectively
Line and the second ray;
Using the common factor of first ray and the blood vessel mathematical model intermediate cam shape dough sheet as the first end point, will
The common factor of second ray and the blood vessel mathematical model intermediate cam shape dough sheet obtains institute as second end points
The index information of the triangle surface where stating first end point and the second end points.
4. edit methods of blood vessel mathematical model according to claim 3, it is characterised in that by described first
The common factor of ray and the blood vessel mathematical model intermediate cam shape dough sheet as the first end point, by second ray with
The common factor of the blood vessel mathematical model intermediate cam shape dough sheet as second end points, including:
The blood vessel mathematical model is split by octree partition method, the blood vessel mathematical model is divided into many
Sub-spaces, the quantity of each subspace intermediate cam shape dough sheet is less than a threshold value;
The common factor of the triangle surface in first ray and subspace is the first end point, second ray with
The common factor of the triangle surface in subspace is second end points.
5. edit methods of blood vessel mathematical model according to claim 3, it is characterised in that according to described
The index information and adjacent surface information of the triangle surface where end point, the second end points determine multiple blood vessel annulus compositions
Target blood section, including:
Step a:First end point A and the second terminal B constitute vectorIt is determined that excessively described first end point A and vertical
In the vectorPlane C, the plane C intersects with the first triangle surface comprising the first end point A;
Step b:The index information of first triangle surface is stored in array 0, and marks first triangle
Shape dough sheet;
Step c:Obtain in the proximal surface of the first triangle surface both sides, intersect and be not labeled with the plane C
The first proximal surface, the index information of first proximal surface is stored in the array 0, and mark first proximal surface;
Step d:Obtain in the proximal surface of first proximal surface, intersect with the plane C and the second not labeled proximal surface,
The index information of second proximal surface is stored in the array 0, and marks second proximal surface;
Step e:Using second proximal surface as the first new proximal surface, circulation performs the step d, until described in two
Second proximal surface overlaps;
Step f:The triangle surface composition blood vessel annulus A that will be preserved in the array 00;
Step g:By in the proximal surface of the triangle surface in array n along the vectorThe index information of the proximal surface of extension
It is stored in array n+1, and the triangle surface composition blood vessel annulus A that will be preserved in the array n+1n+1, n is non-
Negative integer, and the initial value of n is 0;
Step h:If comprising the triangle surface where second terminal B in array n+1, performing step i;
If adding 1 by n not comprising the triangle surface where second terminal B in array n+1, and return to execution institute
State step g;
Step i:By the blood vessel annulus A of 0~array of array n+1 correspondence compositions0~blood vessel annulus An+1Composition is described
Target blood section.
6. edit methods of blood vessel mathematical model according to claim 5, it is characterised in that if first end point
The coordinate of A is (xa, ya, za), the coordinate of the second terminal B is (xb, yb, zb), then the equation table of the plane C
It is shown as:
(x-xa)nx+(y-ya)ny+(z-za)nz=0,
Wherein, (nx, ny, nz) it is vectorDirection vector:
7. edit methods of blood vessel mathematical model according to claim 5, it is characterised in that according to each described
The summit of blood vessel annulus intermediate cam shape dough sheet and target blood section described in the relative position editor of respective central point, including:
According to summit P and central point C on each blood vessel annulusnPosition determine from summit P to central point Cn's
Vector
According to formula P=P+tPCn, t ∈ (0,1) make point P alongDirection move t*l sizes distance, make institute
State blood vessel annulus AnTo central point CnDirection movement.
8. a kind of editing system of blood vessel mathematical model, it is characterised in that described editing system includes:
Triangle surface acquiring unit, the multiple triangle surfaces for obtaining blood vessel mathematical model;
Adjacent surface information acquisition unit, for analyzing and obtains the adjacent surface information of each triangle surface;
Index information determining unit, for receiving and second point that user selects in screen coordinates at first point, and respectively
Determine the triangle of the first end point, the second end points and place in and the corresponding blood vessel mathematical model of second point at described first point
The index information of shape dough sheet;
Target blood section determining unit, for the rope of the triangle surface according to where the first end point, the second end points
Fuse ceases and adjacent surface information determines the target blood section that multiple blood vessel annulus are constituted;
Central point determining unit, for determining each blood vessel according to the midpoint of each blood vessel annulus intermediate cam shape dough sheet
The central point of annulus;
Target blood section edit cell, for the summit according to each blood vessel annulus intermediate cam shape dough sheet with it is respective in
Target blood section described in the relative position editor of heart point.
9. the editing system of blood vessel mathematical model according to claim 8, it is characterised in that the adjacent surface
Information acquisition unit specifically for:
For each summit V on the blood vessel mathematical modeli, obtain all comprising described on the blood vessel mathematical model
Summit ViTriangle surface, and its face index information is stored in and each summit ViCorresponding array AiIn, wherein,
I is positive integer;
To each triangle surface T on the blood vessel mathematical modeli, each summit Vj, in the summit VjIt is right
The array A for answeringjMiddle acquisition includes the triangle surface TiNext summit A(j+1) %3Triangle surface T 'iFace
Index information;
By the triangle surface T 'iFace index information be stored in and the triangle surface TiCorresponding array BiIn, institute
State array BiIt is the triangle surface TiThe index information of the triangle surface for being abutted, wherein, j=0,1,2.
10. the editing system of blood vessel mathematical model according to claim 9, it is characterised in that the index letter
Breath determining unit specifically for:
Receive and second point that user selects in screen coordinates at first point, described first point and second point are sat in screen
In coordinate transformation to three dimensional space coordinate in mark system, first three-dimensional coordinate and second point three-dimensional coordinate are generated;
Determine that first three-dimensional coordinate and second point three-dimensional coordinate first are penetrated with what virtual reference coordinate was constituted respectively
Line and the second ray;
Using the common factor of first ray and the blood vessel mathematical model intermediate cam shape dough sheet as the first end point, will
The common factor of second ray and the blood vessel mathematical model intermediate cam shape dough sheet obtains institute as second end points
The index information of the triangle surface where stating first end point and the second end points.
The editing system of 11. blood vessel mathematical models according to claim 10, it is characterised in that by described
The common factor of one ray and the blood vessel mathematical model intermediate cam shape dough sheet as the first end point, by second ray
With the common factor of the blood vessel mathematical model intermediate cam shape dough sheet as second end points, including:
The blood vessel mathematical model is split by octree partition method, the blood vessel mathematical model is divided into many
Sub-spaces, the quantity of each subspace intermediate cam shape dough sheet is less than a threshold value;
The common factor of the triangle surface in first ray and subspace is the first end point, second ray with
The common factor of the triangle surface in subspace is second end points.
The editing system of 12. blood vessel mathematical models according to claim 10, it is characterised in that the target
Vessel segment determining unit is specifically for performing following steps:
Step a:First end point A and the second terminal B constitute vectorIt is determined that excessively described first end point A and vertical
In the vectorPlane C, the plane C intersects with the first triangle surface comprising the first end point A;
Step b:The index information of first triangle surface is stored in array 0, and marks first triangle
Shape dough sheet;
Step c:Obtain in the proximal surface of the first triangle surface both sides, intersect and be not labeled with the plane C
The first proximal surface, the index information of first proximal surface is stored in the array 0, and mark first proximal surface;
Step d:Obtain in the proximal surface of first proximal surface, intersect with the plane C and the second not labeled proximal surface,
The index information of second proximal surface is stored in the array 0, and marks second proximal surface;
Step e:Using second proximal surface as the first new proximal surface, circulation performs the step d, until described in two
Second proximal surface overlaps;
Step f:The triangle surface composition blood vessel annulus A that will be preserved in the array 00;
Step g:By in the proximal surface of the triangle surface in array n along the vectorThe index information of the proximal surface of extension
It is stored in array n+1, and the triangle surface composition blood vessel annulus A that will be preserved in the array n+1n+1, n is non-
Negative integer, and the initial value of n is 0;
Step h:If comprising the triangle surface where second terminal B in array n+1, performing step i;
If adding 1 by n not comprising the triangle surface where second terminal B in array n+1, and return to execution institute
State step g;
Step i:By the blood vessel annulus A of 0~array of array n+1 correspondence compositions0~blood vessel annulus AnComposition is described
Target blood section.
The editing system of 13. blood vessel mathematical models according to claim 12, it is characterised in that if first end
The coordinate of point A is (xa, ya, za), the coordinate of the second terminal B is (xb, yb, zb), then the equation table of the plane C
It is shown as:
(x-xa)nx+(y-ya)ny+(z-za)nz=0,
Wherein, (nx, ny, nz) it is vectorDirection vector:
The editing system of 14. blood vessel mathematical models according to claim 12, it is characterised in that the target
Vessel segment edit cell specifically for:
According to summit P and central point C on each blood vessel annulusnPosition determine from summit P to central point Cn's
Vector
According to formula P=P+tPCn, t ∈ (0,1) make point P alongDirection move t*l sizes distance, make institute
State blood vessel annulus AnTo central point CnDirection movement.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510991272.7A CN106920282B (en) | 2015-12-25 | 2015-12-25 | Editing method and system of blood vessel digital model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510991272.7A CN106920282B (en) | 2015-12-25 | 2015-12-25 | Editing method and system of blood vessel digital model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106920282A true CN106920282A (en) | 2017-07-04 |
CN106920282B CN106920282B (en) | 2020-11-20 |
Family
ID=59459899
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510991272.7A Active CN106920282B (en) | 2015-12-25 | 2015-12-25 | Editing method and system of blood vessel digital model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106920282B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108447552A (en) * | 2018-03-05 | 2018-08-24 | 四川和生视界医药技术开发有限公司 | The edit methods of retinal vessel and the editing device of retinal vessel |
CN110956698A (en) * | 2019-11-21 | 2020-04-03 | 中国科学院深圳先进技术研究院 | Method, system and electronic device for generating tumor in blood vessel digital model |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101331517A (en) * | 2005-12-19 | 2008-12-24 | 皇家飞利浦电子股份有限公司 | Method for facilitating post-processing of images using deformable meshes |
WO2009011919A8 (en) * | 2007-07-17 | 2009-05-14 | Searete Llc | Active blood vessel sleeve |
US20110235883A1 (en) * | 2010-03-26 | 2011-09-29 | Fujitsu Limited | Three-dimensional template transformation method and apparatus |
CN103310483A (en) * | 2013-05-27 | 2013-09-18 | 重庆大学 | Vascular bifurcating surface reconstructed method based on cavity area triangulation |
US20150216548A1 (en) * | 2014-02-04 | 2015-08-06 | Numata Optical Instruments Factory | Venous valve incising blades, manufacturing method of the same and artery revascularization treatment using the same |
-
2015
- 2015-12-25 CN CN201510991272.7A patent/CN106920282B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101331517A (en) * | 2005-12-19 | 2008-12-24 | 皇家飞利浦电子股份有限公司 | Method for facilitating post-processing of images using deformable meshes |
WO2009011919A8 (en) * | 2007-07-17 | 2009-05-14 | Searete Llc | Active blood vessel sleeve |
US20110235883A1 (en) * | 2010-03-26 | 2011-09-29 | Fujitsu Limited | Three-dimensional template transformation method and apparatus |
CN103310483A (en) * | 2013-05-27 | 2013-09-18 | 重庆大学 | Vascular bifurcating surface reconstructed method based on cavity area triangulation |
US20150216548A1 (en) * | 2014-02-04 | 2015-08-06 | Numata Optical Instruments Factory | Venous valve incising blades, manufacturing method of the same and artery revascularization treatment using the same |
Non-Patent Citations (2)
Title |
---|
JIANHUANG WU ET AL.: ""Scale-adaptive surface modeling of vascular structures"", 《BIOMEDICAL ENGINEERING ONLINE》 * |
钟南昌 等: ""类Delaunay三角剖分的多分叉血管建模方法"", 《计算机辅助设计与图形学学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108447552A (en) * | 2018-03-05 | 2018-08-24 | 四川和生视界医药技术开发有限公司 | The edit methods of retinal vessel and the editing device of retinal vessel |
CN110956698A (en) * | 2019-11-21 | 2020-04-03 | 中国科学院深圳先进技术研究院 | Method, system and electronic device for generating tumor in blood vessel digital model |
CN110956698B (en) * | 2019-11-21 | 2023-02-10 | 中国科学院深圳先进技术研究院 | Method, system and electronic device for generating tumor in blood vessel digital model |
Also Published As
Publication number | Publication date |
---|---|
CN106920282B (en) | 2020-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103049907B (en) | Interactive image segmentation method | |
JP5241357B2 (en) | Biological data model creation method and apparatus | |
US20170035517A1 (en) | Dynamic and interactive navigation in a surgical environment | |
CN110087577A (en) | For the method and system from denture grid removal tooth set | |
CN106934807A (en) | A kind of medical image analysis method, system and Medical Devices | |
CN106682424A (en) | Medical image adjusting method and medical image adjusting system | |
CN110033465A (en) | A kind of real-time three-dimensional method for reconstructing applied to binocular endoscope medical image | |
CN110189352A (en) | A kind of root of the tooth extracting method based on oral cavity CBCT image | |
CN104039229A (en) | Real-time display of vasculature views for optimal device navigation | |
CN110110445B (en) | Sunshine analysis method, device, equipment and storage medium | |
US9491443B2 (en) | Image processing method and image processing apparatus | |
CN104240287A (en) | Method and system for generating coronary panorama through CT images | |
CN104715475A (en) | Automatic whole dental crown segmentation method through harmonic-field-based three-dimensional tooth jaw model | |
CN109069122A (en) | The system and method for determining the size of the structure in medical image | |
CN106920282A (en) | The edit methods and system of a kind of blood vessel mathematical model | |
Berkley et al. | Creating fast finite element models from medical images | |
CN100588370C (en) | Three-dimensional modeling method for funnel chest of children | |
CN114723893A (en) | Organ tissue spatial relationship rendering method and system based on medical images | |
Linney et al. | Three-dimensional visualization of computerized tomography and laser scan data for the simulation of maxillo-facial surgery | |
CN108510506A (en) | A kind of tubular structure image partition method | |
JP2012245374A (en) | Biological data model creating method and device therefor, data structure of biological data model and data storage device for biological data model, and load distribution method for three dimensional data model and device therefor | |
JP2015128669A (en) | Living body data model creation method and device, living body data model data structure, living body data model data storage device, three dimensional data model load distribution method, and device | |
CN102609623A (en) | Ablation therapy image guide device with two-dimensional image processing device | |
CN102609620A (en) | Ablation therapy image guide device with image segmenting device | |
CN109561870A (en) | Excision process estimation device and excision process wizard system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |