CN104794758B - A kind of method of cutting out of 3-D view - Google Patents
A kind of method of cutting out of 3-D view Download PDFInfo
- Publication number
- CN104794758B CN104794758B CN201510185590.4A CN201510185590A CN104794758B CN 104794758 B CN104794758 B CN 104794758B CN 201510185590 A CN201510185590 A CN 201510185590A CN 104794758 B CN104794758 B CN 104794758B
- Authority
- CN
- China
- Prior art keywords
- dimensional coordinate
- voxel
- dimensional
- coordinate system
- screen
- 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.)
- Active
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
- G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
Landscapes
- Engineering & Computer Science (AREA)
- Architecture (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)
- Image Generation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The invention discloses a kind of method of cutting out of 3-D view.Including：Image iso-surface patch is carried out according to the face of object element, three dimensional rendered images of the object under screen two-dimensional coordinate system are obtained, each voxel of the object corresponding two-dimensional coordinate under the screen two-dimensional coordinate system is obtained；Or, image volume drawing is carried out according to the voxel of object, each voxels of three dimensional rendered images and the object of the object under screen two-dimensional coordinate system corresponding two-dimensional coordinate under the screen two-dimensional coordinate system is obtained；Determine to treat clipping region on the three dimensional rendered images；Give up the corresponding voxel of two-dimensional coordinate where the three dimensional rendered images pixel treated in clipping region, and its corresponding face element；Face element according to remaining re-starts image iso-surface patch or re-starts image volume drawing according to the voxel remained.
Description
Technical field
The present invention relates to technical field of image processing, more particularly to a kind of method of cutting out of 3-D view.
Background technology
Three-dimensional (3D) image of perfection either in order to obtain object, it is desired nonetheless to simulated on the 3D rendering of object to this
The operation of object, is required for cutting 3D rendering.
By taking medical scene as an example, complicated accurately surgical operation increasingly be unable to do without 3D medical image analysis and computer aided manufacturing
Help surgery systems.It is difficult once to generate perfect 3D images in 3D medical image analysis systems, it is right from 3D images to generally require
It is cut, to form perfect 3D images.In review of computer aided surgery system, if necessary to directly enterprising in 3D images
Row surgery planning and surgical simulation, need also exist for cutting 3D images.
Therefore, the cutting to 3D rendering how is realized, is the problem of needing to solve at present.
The content of the invention
It is an object of the invention to provide a kind of method of cutting out of 3-D view, to realize the cutting to 3D rendering.
The purpose of the present invention is achieved through the following technical solutions：
A kind of method of cutting out of 3-D view, including：
Image iso-surface patch is carried out according to the face of object element, Three-dimensional Display of the object under screen two-dimensional coordinate system is obtained
Image, obtains each voxel of the object corresponding two-dimensional coordinate under the screen two-dimensional coordinate system；Or, according to object
Voxel carry out image volume drawing, obtain three dimensional rendered images and the object of the object under screen two-dimensional coordinate system
The each voxel corresponding two-dimensional coordinate under the screen two-dimensional coordinate system；
Determine to treat clipping region on the three dimensional rendered images；
Give up the corresponding body of two-dimensional coordinate where the three dimensional rendered images pixel treated in clipping region
Element, and its corresponding face element；
Face element according to remaining re-starts image iso-surface patch or re-starts image according to the voxel remained
Volume drawing.
Method provided in an embodiment of the present invention, realizes and 3-D view is cut.Also, due to obtaining object
Each voxel corresponding two-dimensional coordinate under screen two-dimensional coordinate system, therefore can be according to voxel obtained above and two-dimensional coordinate
Corresponding relation give up where three dimensional rendered images pixel in clipping region is treated the corresponding voxel of two-dimensional coordinate and its
Corresponding face element.Due to the compared pixels point position relationship under two-dimensional coordinate system, its operand is smaller, so as to improve cutting speed
Degree, realizes that real-time 3-D view is cut.
Brief description of the drawings
Fig. 1 is a kind of method flow diagram provided in an embodiment of the present invention；
Fig. 2 is the method flow diagram provided in an embodiment of the present invention for carrying out Image Rendering；
Fig. 3 is the method flow diagram provided in an embodiment of the present invention that projective transformation is carried out to voxel；
Fig. 4 be three dimensional rendered images pixel provided in an embodiment of the present invention whether the determination methods stream in closed curve
Cheng Tu；
Fig. 5 is another method flow diagram provided in an embodiment of the present invention.
Embodiment
Before the embodiment of the present invention is described in detail, the part noun first to being related in the embodiment of the present invention
It is explained.
Object, refers to the three-dimensional body shown by 3-D view.By taking medical domain as an example, refer to some histoorgan, example
Such as liver, heart etc..
World coordinate system, the coordinate system of definition space where three-dimensional body is a three-dimensional system of coordinate.Assuming that for certain
Object, is made up of in image sequence, the embodiment of the present invention N sectioning images, and the upper left corner for defining first sectioning image is generation
The origin of coordinates of boundary's coordinate system, the line direction of sectioning image is the x-axis positive direction of world coordinate system, and the column direction of sectioning image is
The y-axis positive direction of world coordinate system, the direction of image sequence is z-axis positive direction, and long measure is 1 pixel.It is of the invention real
Apply in example, the coordinate in world coordinate system and the coordinate in the coordinate system of three dimensional rendered images are corresponded.
Screen two-dimensional coordinate system, the upper left corner for defining on-screen display (osd) area is the origin of screen two-dimensional coordinate system, and screen shows
The line direction for showing region is the x-axis positive direction of screen two-dimensional coordinate system, and the column direction of on-screen display (osd) area is screen two-dimensional coordinate
The y-axis positive direction of system, long measure is 1 pixel.
A kind of voxel, 3 d image data can be obtained from the image sequence of object.
Face element, another 3 d image data can be obtained from the image sequence of object, the phase under world coordinate system
There is corresponding relation between the voxel at coordinate and face element.
Three dimensional rendered images, the image of the 3-D effect shown on screen.
Fig. 1 show the method for cutting out of 3-D view provided in an embodiment of the present invention, and its specific implementation includes as follows
Operation：
Step 100, the face element progress image iso-surface patch according to object, obtain the object under screen two-dimensional coordinate system
Three dimensional rendered images, obtain each voxel of the object corresponding two-dimensional coordinate under the screen two-dimensional coordinate system；Or,
Image volume drawing is carried out according to the voxel of object, three dimensional rendered images and institute of the object under screen two-dimensional coordinate system are obtained
State each voxel of object corresponding two-dimensional coordinate under the screen two-dimensional coordinate system.
Carry out image volume drawing or iso-surface patch includes object on screen in three dimensions.
In addition, the embodiment of the present invention does not limit the projective transformation mode that drawing process is used yet, it can both use parallel
Projection, it would however also be possible to employ perspective projection.
Step 110, determine to treat clipping region on above-mentioned three dimensional rendered images.
Step 120, give up positioned at this treat three dimensional rendered images pixel in clipping region where two-dimensional coordinate it is corresponding
Voxel, and the corresponding face element of these voxels.
Step 130, re-started according to the face element that remains image iso-surface patch or according to the voxel remained again
Carry out image volume drawing.
The step is that the object after 3-D view cutting is shown on screen.
Method provided in an embodiment of the present invention, realizes and 3-D view is cut.Also, due to obtaining object
Each voxel corresponding two-dimensional coordinate under screen two-dimensional coordinate system, therefore can be according to voxel obtained above and two-dimensional coordinate
Corresponding relation give up where three dimensional rendered images pixel in clipping region is treated the corresponding voxel of two-dimensional coordinate and its
Corresponding face element.Due to the compared pixels point position relationship under two-dimensional coordinate system, its operand is smaller, so as to improve reduction speed
Degree, realizes that real-time 3-D view is cut.Further, since under screen two-dimensional coordinate system comparison position relation, Jin Ergen
According to each voxel of object, corresponding two-dimensional coordinate gives up corresponding voxel and face element under screen two-dimensional coordinate system, it is seen that be not required to
Build and treat the corresponding 3D region in clipping region, and directly 3 d image data is entered with 3D region under three-dimensional system of coordinate
Row position judgment.Wherein, above-mentioned 3D region needs to be built along the opposite direction of volume drawing or the projecting direction of iso-surface patch,
System operations ability is limited to, parallel projection can only be used when carrying out volume drawing or iso-surface patch.Because the embodiment of the present invention is carried
The scheme of confession, which need not be built, treats the corresponding 3D region in clipping region, also avoids the need for limiting and carries out volume drawing or iso-surface patch
Projection pattern, can both use parallel projection, it would however also be possible to employ perspective projection.Therefore it can select suitable according to actual needs
Projection pattern, in addition, if using perspective projection, additionally it is possible to simulate the visual signature of human eye " near big and far smaller ", brings preferably
Consumer's Experience.
Each step to the embodiment of the present invention is described in detail separately below.
The specific implementation of step 100：
Image iso-surface patch is carried out according to the face of object element or the specific reality of image volume drawing is carried out according to the voxel of object
Existing mode has a variety of, and this is not limited by the present invention.Only it is illustrated by taking the flow shown in Fig. 2 as an example, as shown in Fig. 2
Step 100 specifically includes following operation：
Step 200, structure world coordinate system；
Step 210,3 d image data (voxel or face element) carries out to model conversion, and model conversion used in recording
Matrix.
Specifically, be rotated, Pan and Zoom, put appropriate location in the scene similar to three-dimensional body.
Step 220, by after model is converted 3 d image data carry out projective transformation, and record used in project
Transformation matrix.
In the step, parallel projection can be both carried out, perspective projection can also be carried out.
Step 230, by by model convert and projective transformation after obtained data display in the designated area of screen window,
And record viewport transform matrix.
It should be pointed out that above-mentioned processing procedure is not only suitable for volume drawing, iso-surface patch is also applied for.
Wherein, if volume drawing, then model conversion, projective transformation and the viewport transform are carried out according to voxel successively；If
It is iso-surface patch, then model conversion, projective transformation and the viewport transform is carried out according to face element successively.
200~step 230 of above-mentioned steps can with but be not limited only to existing Image Rendering platform (such as OpenGl put down
Platform) on realize.
If carrying out iso-surface patch, accordingly, the realization of voxel corresponding two-dimensional coordinate under screen two-dimensional coordinate system is obtained
Mode also has a variety of.
For example, carrying out projective transformation to the voxel of above-mentioned object, each voxel correspondence under screen two-dimensional coordinate system is obtained
Two-dimensional coordinate.The process namely converted by simulated projections, sets up the mapping of world coordinate system and screen two-dimensional coordinate system
Relation.No matter volume drawing or iso-surface patch are used, can use this implementation.During in particular by iso-surface patch, due to face
Voxel corresponding two-dimensional coordinate under screen two-dimensional coordinate system can not be obtained in drawing process, therefore the realization side can be used
Formula.The implementation is as shown in figure 3, specifically include following operation：
Step 300, three-dimensional coordinate of each voxel in world coordinate system for obtaining object.
Step 310, each above-mentioned three-dimensional coordinate is multiplied by model transformation matrix, projective transformation matrix and the viewport transform successively
Matrix, obtains each voxel of object corresponding two-dimensional coordinate under above-mentioned screen two-dimensional coordinate system.
Wherein, model transformation matrix, projective transformation matrix and viewport transform matrix are projective parameter.
It should be pointed out that regarding model transformation matrix, projective transformation matrix and viewport transform matrix as projective parameter
Be citing and it is non-limiting.
The embodiment of the present invention is not limited the acquisition modes of projective parameter.For example, it is possible to use step 100 image surface is painted
Projective parameter used in system or volume drawing process.Object, can completely, accurately be shown as by existing image rendering technologies
Three dimensional rendered images, this gives the credit to used projective parameter during it.Therefore, using used in image rendering process
Projective parameter carries out simulated projections, each voxel accurately can be mapped in screen two-dimensional coordinate system, it is ensured that voxel is with regarding
The coordinate of three dimensional object in feel is corresponded all the time, it is ensured that the follow-up accuracy for cutting process.
If carrying out image volume drawing, the projective transformation during volume drawing is utilized, you can obtain each voxel of object
The corresponding two-dimensional coordinate under above-mentioned screen two-dimensional coordinate system.Existing image rendering technologies, can be by object completely, accurately
Three dimensional rendered images are shown as, this gives the credit to used projective parameter during it.Therefore, obtained using in image rendering process
To the accurate mapping relations in screen two-dimensional coordinate system of each voxel, it is ensured that the follow-up accuracy for cutting process.
In order to further improve in processing speed, above-mentioned steps 100, if carrying out image iso-surface patch, it is possible to use CUDA
Each voxel of the accelerated mode acquisition object corresponding two-dimensional coordinate under the screen two-dimensional coordinate system.For example,
Speeded up to by GPU.
The specific implementation of step 110：
It is determined that treating the implementation of clipping region has a variety of.For example, it is aobvious in three-dimensional to obtain user according to mouse motion track
That is chosen on diagram picture treats the closed curve of clipping region, accordingly, treats that clipping region just refers to the region in closed curve.Again
For example, obtaining the closed curve that user chooses on three dimensional rendered images according to mouse motion track, anti-selection operation is then carried out,
Accordingly, it is exactly the region outside closed curve to treat clipping region.Etc..
Because the situation that mouse motion track has discontinuity point occurs, that is, non-closed curve is, it is necessary to carry out three times
Spline interpolation.Even if the region chosen manually is not closure, row interpolation formation closed curve can be also entered automatically, it is not interfered with
Judgement afterwards.For example：
It suppose there is n+1 back end (x0,y0),(x1,y1),(x2,y2)……(xn,yn)
Material calculation hi=xi+1-xi(i=0,1 ... ..., n-1)；
Back end and the first end-point condition specified are substituted into matrix equation；
Dematrix equation, tries to achieve second differential value mi, the matrix is triple diagonal matrix；
Calculate the coefficient of SPL：
ai=yi
In each subinterval xi≤x≤xi+1Middle establishment equation：
gi(x)=ai+bi(x-xi)+ci(x-xi)2+di(x-xi)3
Existing implementation is referred to by the specific implementation of cubic spline interpolation formation closed curve, the present invention
Repeat no more.
The specific implementation of above-mentioned steps 120：
Before step 120 is performed, the step of method provided in an embodiment of the present invention also includes judging.Judge on each
State whether the three dimensional rendered images pixel at two-dimensional coordinate is treated in clipping region above-mentioned.
Wherein, if step 110 is specifically to obtain the closure for treating clipping region that user chooses on three dimensional rendered images
Curve, or to user choose closed curve invert after determination treat clipping region.So, in the step of judgement, it is specifically
Judge the three dimensional rendered images pixel at each above-mentioned two-dimensional coordinate whether in the closed curve.It is right exemplified by shown in Fig. 4
The specific implementation of judgment step is illustrated.As shown in figure 4, specifically including following operation：
Step 400, under above-mentioned screen two-dimensional coordinate system build it is same with the viewing area of above-mentioned three dimensional rendered images big
Small rectangular image.
Wherein, the pixel value of whole rectangular image pixels of the rectangular image is the first pixel value.
For example, the pixel value of all matrix image pixels point of the matrix diagram elephant is set into 0.
Because the matrix diagram is as setting up in screen two-dimensional coordinate system, therefore, the matrix image pixel of the matrix diagram elephant
Corresponded with three dimensional rendered images pixel.
For example, the viewing area size of three dimensional rendered images is 100 pixel *, 100 pixels, then, on
The size for stating matrix diagram elephant is also 100 pixel *, 100 pixels.
Step 410, by above-mentioned rectangular image at the corresponding two-dimensional coordinate of above-mentioned closed curve (treating clipping region)
The pixel value of rectangular image pixel is revised as the second pixel value, obtains contour curve.
Wherein, the value of the first pixel value is different from the value of the second pixel value.
Specifically, searching the outline of above-mentioned closed curve using contour detecting, and it is mapped to the corresponding area of rectangular image
Domain, the second pixel value is changed to by the pixel value of the corresponding region, and such as 255, obtain contour curve.
Step 420, judge whether each rectangular image pixel is stated in contour curve at this in above-mentioned rectangular image.
Specifically, being compared judgement using ray method.If ray and profile that a rectangular image pixel is drawn
The intersection point number of curve is even number, then it represents that the rectangular image pixel is in the outside of contour curve, if a rectangular image
The ray and the intersection point number of contour curve that pixel is drawn are odd number, then it represents that the rectangular image pixel is in contour curve
It is internal.
It can specifically be realized to further improve the step of processing speed, above-mentioned judgement using CUDA accelerated modes.Example
Such as, speeded up to by GPU.
The specific implementation of step 130：
Projective transformation mode used in drawing process is not limited, and can both use parallel projection, it would however also be possible to employ
Perspective projection.
It should be pointed out that can be cooperated between the specific implementation of each above-mentioned step.
Based on above-mentioned any means embodiment, before the 3 d image data of object is carried out into Image Rendering, this method
Also include：
Image sequence comprising the object is split, the voxel and face element of the 3-D view of the object is obtained.
Specifically, splitting to the image sequence comprising object, object images sequence is obtained, and then according to the object diagram
As sequence obtains the voxel and face element of object.
By being split, due to not interfere with each other between the three-dimensional coordinate of multiple objects, therefore can be by image sequence
Multiple Object Segmentations of display are different object images sequences, and then are directed to different segmentation purpose and application scenarios, simultaneously
Image cropping processing is carried out to single or multiple object parallels, so as to improve processing speed.Similar to simulation surgical procedure, operation
Knife can cut off the lesion locations of a certain organ-tissue, also cut away adjacent organ-tissue lesion locations simultaneously.Cutting object
It can carry out according to actual needs, this is a very crucial step in simulation surgical procedure, is also what is tallied with the actual situation.
Below by taking concrete application scene as an example, method provided in an embodiment of the present invention is described in detail.Such as Fig. 5 institutes
Show, specifically include following operation：
Step 500, the image sequence in the thoracic cavity for obtaining use CT scan are split, and obtain liver image sequence, enter
And the voxel and face element of liver are obtained according to liver image sequence.
Wherein, the image sequence in thoracic cavity can with but be not limited only to obtain from DICOM.
It should be pointed out that in step 500, segmentation can not only obtain liver image sequence, can also obtain other devices
The image sequence of official, the image sequence of Different Organs is independently preserved.The present embodiment is only entered with the cutting of liver image sequence
Row is illustrated.
If segmentation obtains the image sequence of multiple organs, and these organs are required to carry out 3-D view cutting, then may be used
The 3-D view of these organs is cut with realizing simultaneously.
Step 510, structure world coordinate system.
Assuming that there is N liver image sequences, using the upper left corner of first liver image as the origin of coordinates, line direction is X-axis
Positive direction, column direction is Y-axis positive direction, and the direction of image sequence builds world coordinate system for Z axis positive direction.
Step 520, the voxel of liver or face element carry out to model conversion, and model transformation matrix used in recording.
The conversion of so-called model, i.e., rotated, Pan and Zoom, class to three-dimensional object model (by voxel or face element structure)
It is similar to three-dimensional body and puts appropriate location in the scene.Model transformation matrix used in model conversion is referred to as modeMatrix.
Step 530, the data after model is converted carry out to projective transformation, and projective transformation square used in recording
Battle array.
So-called projective transformation, by taking parallel projection as an example, the what comes into a driver's body of projection is the flat pattern pipeline of a rectangle, no matter object
How far is distance, article size, size constancy after projection.Projective transformation matrix is referred to as used in projective transformation
projMatrix。
Illustrated above with parallel projection, the present embodiment can also use perspective projection, its specific implementation is repeated no more.
Step 540, by by model convert and projective transformation after obtained liver three dimensional rendered images be shown in screen window
The designated area of mouth, and record viewport transform matrix (also known as viewport vector).
Specifically, after model conversion and projective transformation, in addition it is also necessary to by cutting conversion, that is, find a view.After finding a view
Liver three dimensional rendered images are shown in the designated area of screen window.Assuming that the coordinate of each voxel or face element is objcoor
[objx, objy, objz, 1.0], is multiplied by modelview matrix modelMatrix, projective transformation matrix projMatrix successively,
Finally the three dimensional rendered images pixel of homogeneous expression is limited in current viewport using viewport transform matrix.
It should be pointed out that step 510~step 540 is the process for carrying out Image Rendering, before Image Rendering is carried out,
Also need to carry out the conversion (changing into STL forms) of data format.
Step 550, three-dimensional coordinate of each voxel in world coordinate system for obtaining liver.
Step 560, the three-dimensional coordinate of each voxel of liver is multiplied by successively model transformation matrix, projective transformation matrix and
Viewport transform matrix, obtains each voxel of liver corresponding two-dimensional coordinate under above-mentioned screen two-dimensional coordinate system.
Step 570, according to mouse motion track obtain user chosen on liver three dimensional rendered images treat clipping region
Closed curve.
Step 580, under above-mentioned screen two-dimensional coordinate system build and above-mentioned liver three dimensional rendered images viewing area it is same
The rectangular image of sample size.
Step 590, by the rectangular image pixel in above-mentioned rectangular image at the corresponding two-dimensional coordinate of above-mentioned closed curve
Pixel value be revised as the second pixel value, obtain contour curve.
Step 5100, judge whether each rectangular image pixel is stated in contour curve at this in above-mentioned rectangular image.
Step 5110, the liver three dimensional rendered images pixel place two-dimensional coordinate correspondence given up in the closed curve
Voxel, and these voxels corresponding face element.
The liver voxel or face element that step 5120, basis are remained re-start Image Rendering.
It should be understood by those skilled in the art that, embodiments of the invention can be provided as method, system or computer program
Product.Therefore, the present invention can be using the reality in terms of complete hardware embodiment, complete software embodiment or combination software and hardware
Apply the form of example.Moreover, the present invention can be used in one or more computers for wherein including computer usable program code
The computer program production that usable storage medium is implemented on (including but is not limited to magnetic disk storage, CD-ROM, optical memory etc.)
The form of product.
The present invention is the flow with reference to method according to embodiments of the present invention, equipment (system) and computer program product
Figure and/or block diagram are described.It should be understood that can be by every first-class in computer program instructions implementation process figure and/or block diagram
Journey and/or the flow in square frame and flow chart and/or block diagram and/or the combination of square frame.These computer programs can be provided
The processor of all-purpose computer, special-purpose computer, Embedded Processor or other programmable data processing devices is instructed to produce
A raw machine so that produced by the instruction of computer or the computing device of other programmable data processing devices for real
The device for the function of being specified in present one flow of flow chart or one square frame of multiple flows and/or block diagram or multiple square frames.
These computer program instructions, which may be alternatively stored in, can guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works so that the instruction being stored in the computer-readable memory, which is produced, to be included referring to
Make the manufacture of device, the command device realize in one flow of flow chart or multiple flows and/or one square frame of block diagram or
The function of being specified in multiple square frames.
These computer program instructions can be also loaded into computer or other programmable data processing devices so that in meter
Series of operation steps is performed on calculation machine or other programmable devices to produce computer implemented processing, thus in computer or
The instruction performed on other programmable devices is provided for realizing in one flow of flow chart or multiple flows and/or block diagram one
The step of function of being specified in individual square frame or multiple square frames.
, but those skilled in the art once know basic creation although preferred embodiments of the present invention have been described
Property concept, then can make other change and modification to these embodiments.So, appended claims are intended to be construed to include excellent
Select embodiment and fall into having altered and changing for the scope of the invention.
Obviously, those skilled in the art can carry out the essence of various changes and modification without departing from the present invention to the present invention
God and scope.So, if these modifications and variations of the present invention belong to the scope of the claims in the present invention and its equivalent technologies
Within, then the present invention is also intended to comprising including these changes and modification.
Claims (8)
1. a kind of method of cutting out of 3-D view, it is characterised in that including：
Image iso-surface patch is carried out according to the face of object element, three-dimensional display of the object under screen two-dimensional coordinate system is obtained
Picture, obtains each voxel of the object corresponding two-dimensional coordinate under the screen two-dimensional coordinate system；Or, according to object
Voxel carries out image volume drawing, obtains the every of three dimensional rendered images and the object of the object under screen two-dimensional coordinate system
The individual voxel corresponding two-dimensional coordinate under the screen two-dimensional coordinate system；
Determine to treat clipping region on the three dimensional rendered images；
An equal amount of rectangular image in viewing area with the three dimensional rendered images is built under the screen two-dimensional coordinate system,
The pixel value of whole rectangular image pixels of the rectangular image is the first pixel value；
It will treat that the pixel value of the pixel of clipping region is revised as the second pixel value described in the rectangular image, obtain profile bent
Line；Wherein, first pixel value is different from second pixel value；
Judge the pixel of the rectangular image whether in the contour curve；
Give up two-dimensional coordinate where the three dimensional rendered images pixel treated in clipping region in the contour curve
The corresponding voxel, and its corresponding face element；
Face element according to remaining re-starts image iso-surface patch or re-starts image body painting according to the voxel remained
System.
2. according to the method described in claim 1, it is characterised in that if carrying out image iso-surface patch according to the face of object element, obtain
To three dimensional rendered images of the object under screen two-dimensional coordinate system, each voxel of the acquisition is in screen two dimension
Corresponding two-dimensional coordinate under coordinate system, including：
Projective transformation is carried out to the voxel of the object, each voxel of the object is obtained under the screen two-dimensional coordinate system
Corresponding two-dimensional coordinate.
3. method according to claim 2, it is characterised in that the voxel to the object carries out projective transformation, obtains
Each voxel of the object corresponding two-dimensional coordinate under the screen two-dimensional coordinate system is taken, including：
Obtain each three-dimensional coordinate of the voxel in world coordinate system；
Each three-dimensional coordinate is multiplied by model transformation matrix, projective transformation matrix and viewport transform matrix successively, obtains every
The individual voxel corresponding two-dimensional coordinate under the screen two-dimensional coordinate system.
4. method according to claim 3, it is characterised in that the voxel to the object carries out projective transformation, obtains
Take each voxel under the screen two-dimensional coordinate system before corresponding two-dimensional coordinate, this method also includes：
Obtain model transformation matrix, projective transformation matrix according to used in the face element of the object carries out image surface drawing process
With viewport transform matrix.
5. method according to claim 2, it is characterised in that the voxel to the object carries out projective transformation, obtains
Each voxel of the object corresponding two-dimensional coordinate under the screen two-dimensional coordinate system is taken, including：
Projective transformation is carried out to the voxel of the object using CUDA accelerated modes, each voxel of the object is obtained described
Corresponding two-dimensional coordinate under screen two-dimensional coordinate system.
6. according to the method described in claim 1, it is characterised in that judge the three dimensional rendered images at each two-dimensional coordinate
Whether pixel is treated in clipping region described, including：
Using the three dimensional rendered images pixel at each two-dimensional coordinate of CUDA accelerated modes judgement whether described to be cut
Cut in region.
7. the method according to any one of Claims 1 to 4, it is characterised in that image surface is carried out according to the face of object element and painted
Before before system or according to the voxel of object progress image volume drawing, this method also includes：
Image sequence comprising the object is split, the voxel and face element of the object is obtained.
8. the method according to any one of Claims 1 to 4, it is characterised in that on the determination three dimensional rendered images
Treat clipping region, including：
Determine the closed curve for treating clipping region that user chooses on the three dimensional rendered images.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510185590.4A CN104794758B (en) | 2015-04-17 | 2015-04-17 | A kind of method of cutting out of 3-D view |
| PCT/CN2015/081450 WO2016165209A1 (en) | 2015-04-17 | 2015-06-15 | 3d image cropping method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510185590.4A CN104794758B (en) | 2015-04-17 | 2015-04-17 | A kind of method of cutting out of 3-D view |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104794758A CN104794758A (en) | 2015-07-22 |
| CN104794758B true CN104794758B (en) | 2017-10-03 |
Family
ID=53559535
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510185590.4A Active CN104794758B (en) | 2015-04-17 | 2015-04-17 | A kind of method of cutting out of 3-D view |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN104794758B (en) |
| WO (1) | WO2016165209A1 (en) |
Families Citing this family (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107092430B (en) * | 2016-02-18 | 2020-03-24 | 纬创资通（中山）有限公司 | Space drawing scoring method, device and system for scoring space drawing |
| CN105976394B (en) * | 2016-05-31 | 2019-01-29 | 深圳市恩普电子技术有限公司 | A kind of method of cutting out of adaptive adjustment fetus three-D ultrasound data |
| CN107194925A (en) * | 2017-05-31 | 2017-09-22 | 上海联影医疗科技有限公司 | Image processing method and system |
| EP3627442A4 (en) | 2017-05-31 | 2020-05-20 | Shanghai United Imaging Healthcare Co., Ltd. | Method and system for image processing |
| CN107424211B (en) * | 2017-06-15 | 2021-01-05 | 彭志勇 | WebGL volume reconstruction method |
| CN108492299B (en) * | 2018-03-06 | 2022-09-16 | 天津天堰科技股份有限公司 | Cutting method of three-dimensional image |
| CN108830936B (en) * | 2018-05-24 | 2022-07-05 | 讯飞幻境(北京)科技有限公司 | 3D model jitter prevention method and device |
| CN108984262B (en) * | 2018-07-12 | 2021-04-13 | 宁波视睿迪光电有限公司 | Three-dimensional pointer creating method and device and electronic equipment |
| CN109147061B (en) * | 2018-07-20 | 2022-04-01 | 清华大学 | Method for carrying out human-computer interaction editing on segmentation result of volume data |
| CN111899351A (en) * | 2019-05-05 | 2020-11-06 | 中国石油化工股份有限公司 | Screening method for objects of three-dimensional visual scene |
| CN110703963B (en) * | 2019-09-09 | 2021-07-06 | 稿定（厦门）科技有限公司 | 3D character display method, medium, equipment and device |
| CN111325826A (en) * | 2020-01-10 | 2020-06-23 | 郑州航空工业管理学院 | Engineering drawing model projection system and method |
| CN111681315B (en) * | 2020-06-10 | 2022-10-04 | 中国电子科技集团公司第二十八研究所 | High altitude and profile plotting positioning method based on three-dimensional GIS map |
| CN111784823B (en) * | 2020-07-03 | 2023-10-03 | 江苏徐工工程机械研究院有限公司 | Light-weight three-dimensional model mapping-prevention display method and device and storage medium |
| CN112765706B (en) * | 2020-12-31 | 2024-02-20 | 杭州群核信息技术有限公司 | Home decoration material moving method and device, computer equipment and storage medium |
| CN113409333B (en) * | 2021-06-16 | 2022-07-22 | 青岛海信医疗设备股份有限公司 | Three-dimensional image cutting method and electronic equipment |
| CN113656954A (en) * | 2021-08-10 | 2021-11-16 | 北京首钢自动化信息技术有限公司 | Cutting map generation method and device |
| CN113706658B (en) * | 2021-08-18 | 2023-09-05 | 江苏红豆工业互联网有限公司 | Discretization parameter drawing method based on clothing circular arc curve |
| CN114459397B (en) * | 2022-03-01 | 2024-02-20 | 公安部第一研究所 | CT image size measuring method |
| CN115131371B (en) * | 2022-07-12 | 2024-02-09 | 中国建筑西南设计研究院有限公司 | Two-dimensional and three-dimensional hybrid cutting method and terminal for building structural member based on area |
| CN115147440B (en) * | 2022-07-12 | 2024-04-19 | 中国建筑西南设计研究院有限公司 | Automatic generation method of basic slope line based on surface area clipping method |
| CN116228949B (en) * | 2023-02-28 | 2024-03-08 | 如你所视(北京)科技有限公司 | Three-dimensional model processing method, device and storage medium |
| CN116431095B (en) * | 2023-03-23 | 2024-04-30 | 北京凯视达科技股份有限公司 | Panoramic display method, panoramic display device, electronic equipment and storage medium |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102074039A (en) * | 2010-09-29 | 2011-05-25 | 深圳市蓝韵实业有限公司 | Method for drawing volume rendering cutting surface |
| CN102915555A (en) * | 2011-08-03 | 2013-02-06 | 通用电气公司 | Method and system for cropping a 3-dimensional medical dataset |
| CN103049930A (en) * | 2012-12-06 | 2013-04-17 | 南方医科大学 | Method for drawing dynamic Gaussian volume model for three-dimensional medical images |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101593357B (en) * | 2008-05-28 | 2015-06-24 | 中国科学院自动化研究所 | Interactive volume cutting method based on three-dimensional plane control |
| GB201020073D0 (en) * | 2010-11-26 | 2011-01-12 | Siemens Medical Solutions | Anatomically-aware MIP shading |
| CN103810752B (en) * | 2014-02-18 | 2017-01-11 | 海信集团有限公司 | Liver segmenting method based on medical image and liver segmenting system thereof |
| CN104050689A (en) * | 2014-06-10 | 2014-09-17 | 黄淮学院 | Visual object two-dimensional cutting generation method in virtual scene |
| CN104463942B (en) * | 2014-10-31 | 2017-07-28 | 上海联影医疗科技有限公司 | 3-D view method of cutting out and device |
-
2015
- 2015-04-17 CN CN201510185590.4A patent/CN104794758B/en active Active
- 2015-06-15 WO PCT/CN2015/081450 patent/WO2016165209A1/en active Application Filing
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102074039A (en) * | 2010-09-29 | 2011-05-25 | 深圳市蓝韵实业有限公司 | Method for drawing volume rendering cutting surface |
| CN102915555A (en) * | 2011-08-03 | 2013-02-06 | 通用电气公司 | Method and system for cropping a 3-dimensional medical dataset |
| CN103049930A (en) * | 2012-12-06 | 2013-04-17 | 南方医科大学 | Method for drawing dynamic Gaussian volume model for three-dimensional medical images |
Non-Patent Citations (1)
| Title |
|---|
| 医学图像三维可视化技术研究;胡永祥;《中国优秀硕士学位论文全文数据库，信息科技辑》;20050515(第5期);第6，49-55页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104794758A (en) | 2015-07-22 |
| WO2016165209A1 (en) | 2016-10-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104794758B (en) | A kind of method of cutting out of 3-D view | |
| AU2012390266B2 (en) | A computer graphics method for rendering three dimensional scenes | |
| US9177414B2 (en) | Apparatus and method for rendering point cloud using voxel grid | |
| US10529119B2 (en) | Fast rendering of quadrics and marking of silhouettes thereof | |
| EP1775685B1 (en) | Information processing device and program | |
| CN103530907B (en) | Complicated three-dimensional model drawing method based on images | |
| WO2013069215A1 (en) | Image processing device, image processing method and image file data structure | |
| JP2012505445A (en) | Method and system for rendering a 3D model of a 3D object | |
| IL256458A (en) | Fast rendering of quadrics | |
| KR20160088902A (en) | Volume rendering polygons for 3-d printing | |
| US7893938B2 (en) | Rendering anatomical structures with their nearby surrounding area | |
| CN102968791A (en) | Interactive method for three-dimensional (3D) medical image/graphic display and system thereof | |
| Macedo et al. | A semi-automatic markerless augmented reality approach for on-patient volumetric medical data visualization | |
| Haubner et al. | Virtual reality in medicine-computer graphics and interaction techniques | |
| Jung et al. | Using X3D for medical training simulations | |
| Petkov et al. | Interactive visibility retargeting in vr using conformal visualization | |
| Wang et al. | Interactive Exploration of 4D Geometry with Volumetric Halos. | |
| CN114820980A (en) | Three-dimensional reconstruction method and device, electronic equipment and readable storage medium | |
| Stewart et al. | Rebuilding the visible man | |
| Çalışkan et al. | Overview of Computer Graphics and algorithms | |
| Lin et al. | Fusion of disconnected mesh components with branching shapes | |
| i Bartrolı et al. | Visualization techniques for virtual endoscopy | |
| EP4273811A1 (en) | Technique for optimizing rendering parameters of overlays of medical images | |
| Keeler et al. | The spherical visibility map | |
| JP2002150311A (en) | Image processing device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |