CN108537803B - A kind of CT image interactive segmentation method and device - Google Patents
A kind of CT image interactive segmentation method and device Download PDFInfo
- Publication number
- CN108537803B CN108537803B CN201810298411.1A CN201810298411A CN108537803B CN 108537803 B CN108537803 B CN 108537803B CN 201810298411 A CN201810298411 A CN 201810298411A CN 108537803 B CN108537803 B CN 108537803B
- Authority
- CN
- China
- Prior art keywords
- region
- drawing view
- view
- initial data
- primary segmentation
- 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
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/187—Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20036—Morphological image processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The present invention provides a kind of CT image interactive segmentation method and devices, wherein method includes: adjustment CT intensity threshold range, it obtains CT intensity value and belongs to the position to be shown within the scope of intensity threshold, and it treats display area and coloring is marked, it is shown using Volume Rendering Techniques simulation contour surface, shows initial data volume drawing view;The region to be split in initial data volume drawing view is determined by marking tool, cut zone is treated using fill tool and is calculated, and generates primary segmentation as a result, and adjusting primary segmentation result using corrosion expansion tool;On the basis of primary segmentation result, undivided region out in region to be split is chosen using marking tool, segmentation result is generated, shows dividing body drawing view.It is possible thereby to simplify cutting procedure with intuitive interactive means, a large amount of sliced time is saved, more conducively user uses.
Description
Technical field
The present invention relates to digital medical field more particularly to a kind of CT, (Computed Tomography, computer are disconnected
Layer scanning) image interactive segmentation method and device.
Background technique
In CT image segmentation field, popular software has Seg3D and the Germany of Univ Utah USA in the world at present
MITK workbench of Cancer Research Center etc..These softwares provide the tools such as filtering, threshold value, produce to the first of target object
Step segmentation.Next it needs artificially to smear every slice, erasing label, and checks other views confirmation segmentation effects, from
And obtain final segmentation.When practical segmentation data, a large amount of man-machine interactively operation is still needed to, such as need using paintbrush tool
Smearing image slice etc. one by one, and the hundreds of layers of slice number of plies Chang Duoda of a CT image, complete data using this mode and divide
Very time-consuming, the target object segmentation completed in entire CT usually requires 1~2 hour, and is difficult to divide the tiny knot such as blood vessel
Structure.And research work in recent years then proposes several automatic division method according to the difference of cutting object, but these methods are suitable
Relatively narrow with property, precision is relatively low, it is difficult to directly apply in practical operation.
Summary of the invention
The present invention is intended to provide a kind of a kind of CT image for overcoming the above problem or at least being partially solved the above problem
Interactive segmentation method and device can obtain accurate segmentation result by intuitive, simple interaction.
In order to achieve the above objectives, technical solution of the present invention is specifically achieved in that
One aspect of the present invention provides a kind of CT image interactive segmentation method, comprising: adjustment CT intensity threshold range,
It obtains CT intensity value and belongs to the position to be shown within the scope of intensity threshold, and treat display area and coloring is marked, utilize body
Rendering technique simulation contour surface is shown, shows initial data volume drawing view;Initial data volume drawing is determined by marking tool
Region to be split in view, treats cut zone using fill tool and is calculated, and generates primary segmentation as a result, and using
Corrode expansion tool and adjusts primary segmentation result;On the basis of primary segmentation result, area to be split is chosen using marking tool
Undivided region out, generates segmentation result in domain, shows dividing body drawing view.
In addition, while display initial data volume drawing view, further includes: the sagittal of display initial data volume drawing view
Face slice, coronal-plane slice and cross-sectional slices two dimension view, and the sagittal plane of initial data volume drawing view be sliced,
The region for being higher than current global threshold in coronal-plane slice and the two dimension view of cross-sectional slices carries out coloring and shows;And/or it is aobvious
While showing dividing body drawing view, further includes: sagittal plane slice, coronal-plane slice and the cross of display dividing body drawing view
The two dimension view of cross-sectional slice, and the two of the sagittal plane of dividing body drawing view slice, coronal-plane slice and cross-sectional slices
Region in dimensional view higher than current global threshold carries out coloring and shows.
In addition, determining that the region to be split in initial data volume drawing view includes: in initial data by marking tool
In volume drawing view, by dragging mouse, the region near label dragging track drags trail termination as region to be split
Position before depth mutation;Wherein, track is dragged are as follows: at a certain moment of dragging, the extended line that sight passes through mouse pointer
Several intersection points are generated with currently displayed contour surface, the midpoint of the two o'clock wherein intersected at first are taken, by the midpoint at all moment
The broken line connected into.
It is calculated in addition, treating cut zone using fill tool, generating primary segmentation result includes: to be existed using mouse
Seed point is set in initial data volume drawing view, chooses region associated with seed point, and using flood filling algorithm or
Person's algorithm of region growing calculates primary segmentation result.
In addition, the use of corrosion expansion tool adjustment primary segmentation result including: adjustment region associated with seed point
Edge executes three-dimensional morphological erosion/expansion algorithm to region associated with seed point, adjusts primary segmentation result.
In addition, transmission function includes:
Wherein, I is CT intensity value, and T (A), T (B) indicate the threshold value for being respectively mapped to two color labels;And CT intensity value
The part for meeting intensity threshold range is arranged to opaque, and the part for being unsatisfactory for intensity threshold range is arranged to transparent.
Another aspect of the present invention provides a kind of CT image interactive segmentation device, comprising: and display module is strong for adjusting CT
Threshold range is spent, CT intensity value is obtained and belongs to the position to be shown within the scope of intensity threshold, and treat display area and be marked
Coloring is shown using Volume Rendering Techniques simulation contour surface, shows initial data volume drawing view;Primary segmentation module, for leading to
It crosses marking tool and determines region to be split in initial data volume drawing view, treat cut zone using fill tool and carry out
It calculates, generates primary segmentation as a result, and adjusting primary segmentation result using corrosion expansion tool;Divide module, for preliminary
On the basis of segmentation result, undivided region out in region to be split is chosen using marking tool, segmentation result is generated, leads to
Know that display module shows dividing body drawing view.
In addition, display module, is also used to while showing initial data volume drawing view, initial data volume drawing is shown
Sagittal plane slice, the coronal-plane of view are sliced and the two dimension view of cross-sectional slices, and in the arrow of initial data volume drawing view
The region for being higher than current global threshold in the two dimension view of shape face slice, coronal-plane slice and cross-sectional slices carries out coloring and shows
Show;And/or display dividing body drawing view while, display dividing body drawing view sagittal plane slice, coronal-plane slice and
The two dimension view of cross-sectional slices, and in the sagittal plane of dividing body drawing view slice, coronal-plane slice and cross-sectional slices
Region in two dimension view higher than current global threshold carries out coloring and shows.
In addition, primary segmentation module is specifically used for determining that initial data volume drawing regards by marking tool in the following way
Region to be split in figure: primary segmentation module is also used in initial data volume drawing view, by dragging mouse, label
The region near track is dragged as region to be split, position of the dragging trail termination before depth mutation;Wherein, it drags
Dynamic rail mark are as follows: at a certain moment of dragging, if sight is generated by the extended line and currently displayed contour surface of mouse pointer
Dry intersection point takes the midpoint of the two o'clock wherein intersected at first, the broken line that the midpoint at all moment is connected into.
In addition, primary segmentation module is counted specifically for treating cut zone using fill tool in the following way
It calculates, generates primary segmentation result: primary segmentation module, be also used to that kind is arranged in initial data volume drawing view using mouse
It is sub-, region associated with seed point is chosen, and utilize flood filling algorithm or algorithm of region growing, calculates preliminary point
Cut result.
In addition, primary segmentation module is specifically used for adjusting primary segmentation knot using corrosion expansion tool in the following way
Fruit: primary segmentation module is also used to adjust the edge in region associated with seed point, to region associated with seed point
Three-dimensional morphological erosion/expansion algorithm is executed, primary segmentation result is adjusted.
In addition, transmission function includes:
Wherein, I is CT intensity value, and T (A), T (B) indicate the threshold value for being respectively mapped to two color labels;And CT intensity value
The part for meeting intensity threshold range is arranged to opaque, and the part for being unsatisfactory for intensity threshold range is arranged to transparent.
It can be seen that the CT image interactive segmentation method and device provided through the invention, can substitute current Seg3D,
It is directed to the dividing method of CT image in the softwares such as MITK workbench, simplifies cutting procedure with intuitive interactive means, saves a large amount of
Sliced time, more conducively user use.
Detailed description of the invention
In order to illustrate the technical solution of the embodiments of the present invention more clearly, required use in being described below to embodiment
Attached drawing be briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, for
For those skilled in the art, without creative efforts, it can also be obtained according to these attached drawings
His attached drawing.
Fig. 1 is the flow chart of CT image interactive segmentation method provided in an embodiment of the present invention;
Fig. 2 is the structural schematic diagram of CT image interactive segmentation device provided in an embodiment of the present invention;
Fig. 3 is the transmission function schematic diagram of display multiple color provided in an embodiment of the present invention.
Specific embodiment
Exemplary embodiments of the present disclosure are described in more detail below with reference to accompanying drawings.Although showing the disclosure in attached drawing
Exemplary embodiment, it being understood, however, that may be realized in various forms the disclosure without should be by embodiments set forth here
It is limited.On the contrary, these embodiments are provided to facilitate a more thoroughly understanding of the present invention, and can be by the model of the disclosure
It encloses and is fully disclosed to those skilled in the art.
For the segmentation problem of CT image data, it to include original that the invention proposes a kind of CT image interactive segmentation methods
The intuitive display of beginning data and segmentation result.
Fig. 1 shows the flow chart of CT image interactive segmentation method provided in an embodiment of the present invention, referring to Fig. 1, CT image
Interactive segmentation method includes:
S101 adjusts CT intensity threshold range, obtains CT intensity value and belongs to the position to be shown within the scope of intensity threshold, and
It treats display area and coloring is marked, shown using Volume Rendering Techniques simulation contour surface, display initial data volume drawing view
Figure.
Specifically, the invention proposes the sides for using Volume Rendering Techniques simulation contour surface to show, and colour for data markers
Method.
Following transmission function is designed, volume drawing can be made only to show the part greater than threshold value:
That is, the part that CT intensity value is greater than threshold value is arranged to opaque, it is arranged to transparent less than the part of threshold value.It adjusts
Save threshold value can contour surface under simulative display difference threshold value, volume drawing is calculated by video card, so calculating is very fast,
It can real-time display.
Simultaneously in order to which to there is flag data to show multiple color, the present invention constructs a special transmission function:
Wherein, I is CT intensity value, and T (A), T (B) indicate the threshold value for being respectively mapped to two color labels;And CT intensity value
The part for meeting intensity threshold range is arranged to opaque, and the part for being unsatisfactory for intensity threshold range is arranged to transparent.
Specifically, Fig. 3 is display multiple color transmission function schematic diagram:
It, only need to be by this region together with the transparent region near it when needing the region of segmentation with other color marks
CT value takes opposite number.The purpose that neighbouring transparent region is also taken opposite number is to utilize the linear interpolation of volume drawing, makes region
Boundary seems more smooth.
S102 determines the region to be split in initial data volume drawing view by marking tool, uses fill tool pair
Region to be split is calculated, and generates primary segmentation as a result, and adjusting primary segmentation result using corrosion expansion tool.
Specifically, it the present invention provides several convenient interactive segmentation tools, can be chosen in volume data using these tools
1. these voxels after choosing, can be labeled as segmentation, 2. exclude from segmentation or 3. hide by qualified certain voxels.
As an optional embodiment of the embodiment of the present invention, partition tools may include: marking tool, filling work
Tool, area locking tool, corrosion expansion tool etc..
Wherein, as an optional embodiment of the embodiment of the present invention, initial data body painting is determined by marking tool
Region to be split in view processed includes: in initial data volume drawing view, and by dragging mouse, label dragging track is attached
Close region drags position of the trail termination before depth mutation as region to be split;Wherein, track is dragged are as follows: dragging
The a certain moment, sight generates several intersection points by the extended line of mouse pointer and currently displayed contour surface, takes wherein most
The midpoint of the two o'clock first intersected, the broken line that the midpoint at all moment is connected into.Specifically, marking tool is for choosing body number
The region specified according to middle mouse.In initial data volume drawing view, is dragged, can be marked along object to be split using mouse
Drag the region near track.Wherein, dragging track definition is as follows: at a certain moment of dragging, sight passes through mouse pointer
Extended line and currently displayed contour surface will generate several intersection points, take the midpoint of the two o'clock wherein intersected at first, all moment
A series of broken line that connects into of such midpoints be to drag track.The region that the toolmark goes out is a tubular area,
The axis of pipe is above-mentioned dragging track, and diameter is the 1.5 of the above-mentioned median for intersecting a series of distances calculated by two o'clock at first
Times.Particularly, when the depth that mouse picks up coordinate mutates, dragging track is terminated at the position before depth mutation.Example
Such as, mouse is marked along blood vessel when beginning, at a time out of this angiosomes, then, pickup point will transfer to contour surface
Other positions, or without pickup point, pick up depth at this time and be mutated, dragging track will terminate and stop receiving mouse
Response.The tubular area radius is likely larger than the size of practical desire marked region, so that some background areas are included,
So the tool will only retain the region being connected with dragging track.Thus, it is possible to determine area to be split by marking tool
Domain.
As an optional embodiment of the embodiment of the present invention, cut zone is treated using fill tool and is calculated,
Generating primary segmentation result includes: that seed point is arranged in initial data volume drawing view using mouse, is chosen and seed point phase
Associated region, and flood filling algorithm or algorithm of region growing are utilized, calculate primary segmentation result.Specifically, work is filled
Tool can quickly choose a continuum associated with seed point.Seed point can use mouse directly in volume drawing view
It clicks and chooses, the point of selection is the point that mouse picks up on contour surface.General it is desired that seed point location be all located at contour surface
Inside enclosing region rather than surface, therefore, the depth of seed point is adjusted in the present invention, and method is roll mouse idler wheel, in body painting
The reference line of x, y, z axis can be shown to show seed point location in view processed;Two dimension view can also adjust display position to see
It examines.Such as the straight line that x, y, z reference axis is parallel at seed point can be indicated using red, green, blue line respectively.And judges and plant
The principle of sub- point " associated " is that flood fills (flood fill) algorithm or region growing (region growing) algorithm,
Different filling principles may be selected according to usage scenario difference.Wherein, the object of flood filling algorithm is above current global threshold
Voxel, all voxels being connected with seed point will be filled;Algorithm of region growing is unrelated with current global threshold, from kind
Son point start constantly by with have been filled with the voxel that region is connected and gray difference is little and be included, when not meeting item more
It grows and stops when the voxel of part.Meanwhile area locking tool is usually used in combination with fill tool, to solve fill tool appearance
The problem of easily choosing multizone.Other tools are first used to choose a region as the sphere of action currently paid close attention to, when the work
When tool is in initiate mode, fill tool will be restricted to only execute algorithm in sphere of action.Thus, it is possible to pass through filling
Tool obtains primary segmentation result.
Further, it as an optional embodiment of the embodiment of the present invention, is adjusted using corrosion expansion tool preliminary
Segmentation result includes: the edge in adjustment region associated with seed point, executes three-dimensional to region associated with seed point
Morphological erosion/expansion algorithm adjusts primary segmentation result.Specifically, corrosion expansion tool can be used for adjusting selected areas
Edge executes three-dimensional morphological erosion/expansion algorithm to selected areas.Selected areas edge can be made to receive using erosion algorithm
Contracting can be such that selected areas edge expands using expansion algorithm.Thus, it is possible to adjust primary segmentation knot by corrosion expansion tool
Fruit keeps it more smooth.
S103 chooses undivided area out in region to be split using marking tool on the basis of primary segmentation result
Domain generates segmentation result, shows dividing body drawing view.
Specifically, it if any undivided region out, can also continue to select undivided region, until all areas to be split
It all chooses and finishes in domain.
Meanwhile user from simulation contour surface only it is observed that the medium distribution in some threshold value of volume data, and be difficult to
Other not shown data are observed, grasp these information in cutting procedure for the convenience of the user, one as the embodiment of the present invention
A optional embodiment, while showing initial data volume drawing view, further includes: display initial data volume drawing view
The two dimension view of sagittal plane slice, coronal-plane slice and cross-sectional slices, and cut in the sagittal plane of initial data volume drawing view
The region for being higher than current global threshold in the two dimension view of piece, coronal-plane slice and cross-sectional slices carries out coloring and shows;With/
Or while display dividing body drawing view, further includes: sagittal plane slice, the coronal-plane slice of display dividing body drawing view
With the two dimension view of cross-sectional slices, and the sagittal plane of dividing body drawing view slice, coronal-plane slice and cross-sectional slices
Two dimension view in be higher than the region of current global threshold and carry out coloring and show.The present invention provides three in addition to volume drawing view
A two dimension view is respectively used to the image of display sagittal plane, coronal-plane, cross-sectional slices.In two dimension view, it is higher than current complete
Office threshold value region (data being shown in i.e. current volume drawing view) for example with green display, rest part for example with
Grayscale image is shown.
In addition, can also show these sectioning images in initial data volume drawing view and/or dividing body drawing view, show
Whole slice of display may be selected in mode, can also only show the region higher than threshold value.It in actual use, can be according to usage scenario spirit
Selection display pattern living.
It can be seen that the CT image interactive segmentation method provided through the invention, can substitute current Seg3D, MITK work
Make the dividing method for being directed to CT image in the softwares such as platform, simplify cutting procedure with intuitive interactive means, saves a large amount of segmentation
Time, more conducively user use.
A kind of splitting scheme for certain incidence CT image medium vessels presented below, cutting procedure are as follows:
(1) appropriate threshold is adjusted, shows that blood vessel clearly;
(2) in volume drawing interface, a mouse click picks up a thicker blood vessel, is calculated using the region growing of fill tool
Method generates primary segmentation, is sliced in volume drawing interface display, and adjusts cut zone using corrosion expansion tool according to slice;
(3) minute blood vessel being not yet partitioned into is found, chooses these blood vessels using marking tool;Volume drawing interface with cut
Piece checks that interface can confirm the accuracy of current pickup point, as inaccuracy can adjust at any time.
(4) segmentation is completed, and generates equivalent surface grids.
Fig. 2 shows the structural schematic diagram of CT image interactive segmentation device provided in an embodiment of the present invention, which is handed over
Mutual segmenting device is applied to the above method, only the structure of CT image interactive segmentation device is briefly described below, other are not
Matters to the greatest extent, please refer to the associated description in above-mentioned CT image interactive segmentation method, referring to fig. 2, CT image interactive segmentation device,
Include:
Display module 201, for adjusting CT intensity threshold range, obtain CT intensity value belong within the scope of intensity threshold to
Display area, and treat display area and coloring is marked, it is shown using Volume Rendering Techniques simulation contour surface, shows original number
According to volume drawing view;
Primary segmentation module 202, for determining the region to be split in initial data volume drawing view by marking tool,
Cut zone is treated using fill tool to be calculated, and generates primary segmentation as a result, and adjusting using corrosion expansion tool preliminary
Segmentation result;
Divide module 203, for being chosen in region to be split not on the basis of primary segmentation result using marking tool
The region being partitioned into, generates segmentation result, and notice display module 201 shows dividing body drawing view.
It can be seen that the CT image interactive segmentation device provided through the invention, can substitute current Seg3D, MITK work
Make the dividing method for being directed to CT image in the softwares such as platform, simplify cutting procedure with intuitive interactive means, saves a large amount of segmentation
Time, more conducively user use.
Due to user from simulation contour surface only it is observed that the medium distribution in some threshold value of volume data, and be difficult to see
Other not shown data are examined, grasp these information in cutting procedure for the convenience of the user, one as the embodiment of the present invention
Optional embodiment, display module 201 are also used to while showing initial data volume drawing view, show initial data body
Sagittal plane slice, the coronal-plane of drawing view are sliced and the two dimension view of cross-sectional slices, and in initial data volume drawing view
Sagittal plane slice, coronal-plane slice and cross-sectional slices two dimension view in be higher than current global threshold region carry out
Color is shown;And/or while display dividing body drawing view, sagittal plane slice, the coronal-plane of display dividing body drawing view are cut
The two dimension view of piece and cross-sectional slices, and cut in the sagittal plane of dividing body drawing view slice, coronal-plane slice and cross section
Region in the two dimension view of piece higher than current global threshold carries out coloring and shows.
As an optional embodiment of the embodiment of the present invention, primary segmentation module 202 is specifically used for by such as lower section
Formula determines the region to be split in initial data volume drawing view by marking tool: primary segmentation module 202 is also used to
In initial data volume drawing view, by dragging mouse, the region near label dragging track is as region to be split, dragging
Position of the trail termination before depth mutation;Wherein, track is dragged are as follows: at a certain moment of dragging, sight passes through mouse pointer
Extended line and currently displayed contour surface generate several intersection points, take the midpoint of the two o'clock wherein intersected at first, by sometimes
The broken line that the midpoint at quarter connects into.Thus, it is possible to determine region to be split by marking tool.
As an optional embodiment of the embodiment of the present invention, primary segmentation module 202 is specifically used for by such as lower section
Formula is treated cut zone using fill tool and is calculated, and generate primary segmentation result: primary segmentation module 202 is also used to benefit
Seed point is set in initial data volume drawing view with mouse, chooses region associated with seed point, and fill out using flood
Algorithm or algorithm of region growing are filled, primary segmentation result is calculated.Thus, it is possible to obtain primary segmentation knot by fill tool
Fruit.
As an optional embodiment of the embodiment of the present invention, primary segmentation module 202 is specifically used for by such as lower section
Formula adjusts primary segmentation result: primary segmentation module 202 using corrosion expansion tool, is also used to adjust associated with seed point
Region edge, execute three-dimensional morphological erosion/expansion algorithm to region associated with seed point, preliminary point of adjustment
Cut result.Thus, it is possible to adjust primary segmentation as a result, keeping it more smooth by corrosion expansion tool.
As an optional embodiment of the embodiment of the present invention, transmission function includes:
Wherein, I is CT intensity value, and T (A), T (B) indicate the threshold value for being respectively mapped to two color labels;And CT intensity value
The part for meeting intensity threshold range is arranged to opaque, and the part for being unsatisfactory for intensity threshold range is arranged to transparent.
It should be understood by those skilled in the art that, embodiments herein can provide as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the application
Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the application, which can be used in one or more,
The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) produces
The form of product.
The application is referring to method, the process of equipment (system) and computer program product according to the embodiment of the present application
Figure and/or block diagram describe.It should be understood that can be realized by computer program instructions each in flowchart and/or the block diagram
The combination of process and/or box in process and/or box and flowchart and/or the block diagram.It can provide these computers
Processor of the program instruction to general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices
To generate a machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute
For realizing the function of being specified in one or more flows of the flowchart and/or one or more blocks of the block diagram
Device.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that instruction stored in the computer readable memory generation includes
The manufacture of command device, the command device are realized in one box of one or more flows of the flowchart and/or block diagram
Or the function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer
Or the instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or box
The step of function of being specified in figure one box or multiple boxes.
In a typical configuration, calculating equipment includes one or more processors (CPU), input/output interface, net
Network interface and memory.
Memory may include the non-volatile memory in computer-readable medium, random access memory (RAM) and/
Or the forms such as Nonvolatile memory, such as read-only memory (ROM) or flash memory (flash RAM).Memory is computer-readable Jie
The example of matter.
Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by any method
Or technology come realize information store.Information can be computer readable instructions, data structure, the module of program or other data.
The example of the storage medium of computer include, but are not limited to phase change memory (PRAM), static random access memory (SRAM),
Dynamic random access memory (DRAM), other kinds of random access memory (RAM), read-only memory (ROM), electricity can
Erasable programmable read-only memory (EPROM) (EEPROM), flash memory or other memory techniques, read-only disc read only memory (CD-ROM) (CD-
ROM), digital versatile disc (DVD) or other optical storage, magnetic cassettes, tape magnetic disk storage or other magnetism are deposited
Equipment or any other non-transmission medium are stored up, can be used for storage can be accessed by a computing device information.According to boundary herein
Fixed, computer-readable medium does not include temporary computer readable media (transitory media), such as the data-signal of modulation and
Carrier wave.
The above is only embodiments herein, are not intended to limit this application.To those skilled in the art,
Various changes and changes are possible in this application.All any modifications made within the spirit and principles of the present application are equally replaced
It changes, improve, should be included within the scope of the claims of this application.
Claims (8)
1. a kind of CT image interactive segmentation method characterized by comprising
CT intensity threshold range is adjusted, the CT intensity value is obtained and belongs to the position to be shown within the scope of the intensity threshold, and
Coloring is marked to the position to be shown, is shown using Volume Rendering Techniques simulation contour surface, shows initial data volume drawing
View;
Determine the region to be split in the initial data volume drawing view by marking tool, using fill tool to it is described to
Cut zone is calculated, and generates primary segmentation as a result, and adjusting the primary segmentation result using corrosion expansion tool;
On the basis of the primary segmentation result, undivided area out in the region to be split is chosen using marking tool
Domain generates segmentation result, shows dividing body drawing view;
Wherein: described to determine that the region to be split in the initial data volume drawing view includes: by marking tool
In the initial data volume drawing view, by dragging mouse, region near label dragging track as it is described to
Cut zone, position of the dragging trail termination before depth mutation;Wherein, the dragging track are as follows: in a certain of dragging
Moment, sight generate several intersection points by the extended line and currently displayed contour surface of the mouse pointer, take wherein at first
The midpoint of the two o'clock of intersection, the broken line that the midpoint at all moment is connected into.
2. the method according to claim 1, wherein
While the display initial data volume drawing view, further includes:
Show the sagittal plane slice of the initial data volume drawing view, the two dimension view of coronal-plane slice and cross-sectional slices,
And it is higher than in the two dimension view of the sagittal plane of initial data volume drawing view slice, coronal-plane slice and cross-sectional slices
The region of current global threshold carries out coloring and shows;And/or
While the display dividing body drawing view, further includes:
Show that the sagittal plane of the dividing body drawing view is sliced, coronal-plane is sliced and the two dimension view of cross-sectional slices, and
It is higher than in the two dimension view of the sagittal plane slice of the dividing body drawing view, coronal-plane slice and cross-sectional slices current global
The region of threshold value carries out coloring and shows.
3. the method according to claim 1, wherein described carry out the region to be split using fill tool
It calculates, generating primary segmentation result includes:
Seed point is set in the initial data volume drawing view using mouse, chooses area associated with the seed point
Domain, and flood filling algorithm or algorithm of region growing are utilized, calculate the primary segmentation result.
4. according to the method described in claim 3, it is characterized in that, described adjust the primary segmentation using corrosion expansion tool
Result includes:
The edge for adjusting region associated with the seed point executes three-dimensional shape to region associated with the seed point
State corrosion/expansion algorithm adjusts the primary segmentation result.
5. a kind of CT image interactive segmentation device characterized by comprising
Display module obtains the CT intensity value and belongs within the scope of the intensity threshold for adjusting CT intensity threshold range
Position to be shown, and coloring is marked to the position to be shown, show that display is former using Volume Rendering Techniques simulation contour surface
Beginning data volume drawing view;
Primary segmentation module makes for determining the region to be split in the initial data volume drawing view by marking tool
The region to be split is calculated with fill tool, generates primary segmentation as a result, and adjusting institute using corrosion expansion tool
State primary segmentation result;
Divide module, for being chosen in the region to be split on the basis of the primary segmentation result using marking tool
Undivided region out, generates segmentation result, and the display module is notified to show dividing body drawing view;
Wherein: the primary segmentation module is specifically used for determining the initial data body painting by marking tool in the following way
Region to be split in view processed:
Primary segmentation module is also used in the initial data volume drawing view, and by dragging mouse, label dragging track is attached
Close region is as the region to be split, position of the dragging trail termination before depth mutation;Wherein, the dragging rail
Mark are as follows: at a certain moment of dragging, if sight is generated by the extended line and currently displayed contour surface of the mouse pointer
Dry intersection point takes the midpoint of the two o'clock wherein intersected at first, the broken line that the midpoint at all moment is connected into.
6. device according to claim 5, which is characterized in that
Display module is also used to show the initial data volume drawing view while showing initial data volume drawing view
Sagittal plane slice, coronal-plane slice and cross-sectional slices two dimension view, and in the arrow of the initial data volume drawing view
The region for being higher than current global threshold in the two dimension view of shape face slice, coronal-plane slice and cross-sectional slices carries out coloring and shows
Show;And/or while display dividing body drawing view, show that sagittal plane slice, the coronal-plane of the dividing body drawing view are cut
The two dimension view of piece and cross-sectional slices, and the sagittal plane of the dividing body drawing view slice, coronal-plane slice and it is transversal
The region for being higher than current global threshold in the two dimension view of face slice carries out coloring and shows.
7. device according to claim 5, which is characterized in that the primary segmentation module is specifically used in the following way
The region to be split is calculated using fill tool, generates primary segmentation result:
The primary segmentation module is also used to that seed point is arranged in the initial data volume drawing view using mouse, chooses
Region associated with the seed point, and flood filling algorithm or algorithm of region growing are utilized, calculate the primary segmentation
As a result.
8. device according to claim 7, which is characterized in that the primary segmentation module is specifically used in the following way
The primary segmentation result is adjusted using corrosion expansion tool:
The primary segmentation module is also used to adjust the edge in region associated with the seed point, to the seed point
Associated region executes three-dimensional morphological erosion/expansion algorithm, adjusts the primary segmentation result.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810298411.1A CN108537803B (en) | 2018-03-30 | 2018-03-30 | A kind of CT image interactive segmentation method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810298411.1A CN108537803B (en) | 2018-03-30 | 2018-03-30 | A kind of CT image interactive segmentation method and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108537803A CN108537803A (en) | 2018-09-14 |
CN108537803B true CN108537803B (en) | 2019-08-23 |
Family
ID=63483156
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810298411.1A Active CN108537803B (en) | 2018-03-30 | 2018-03-30 | A kind of CT image interactive segmentation method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108537803B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109872312B (en) * | 2019-02-15 | 2022-12-20 | 腾讯科技(深圳)有限公司 | Medical image segmentation method, device and system and image segmentation method |
US11423544B1 (en) | 2019-11-14 | 2022-08-23 | Seg AI LLC | Segmenting medical images |
US10762629B1 (en) | 2019-11-14 | 2020-09-01 | SegAI LLC | Segmenting medical images |
CN111275717B (en) * | 2020-01-17 | 2023-04-11 | 杭州迪英加科技有限公司 | Contrast display method for different stained sections, electronic device and readable medium |
CN111833427B (en) * | 2020-07-21 | 2021-01-05 | 推想医疗科技股份有限公司 | Method and device for volume rendering of three-dimensional image |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102945127A (en) * | 2012-11-05 | 2013-02-27 | 深圳市旭东数字医学影像技术有限公司 | Interactive method for volume rendering displaying and system thereof |
CN103049907A (en) * | 2012-12-11 | 2013-04-17 | 深圳市旭东数字医学影像技术有限公司 | Interactive image segmentation method |
CN103678847A (en) * | 2012-09-24 | 2014-03-26 | 邵雁 | Skin soft tissue expanding auxiliary system |
CN106570868A (en) * | 2015-10-10 | 2017-04-19 | 中国科学院深圳先进技术研究院 | Three-dimensional liver tumor semi-automatic segmentation method based on level set |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8200015B2 (en) * | 2007-06-22 | 2012-06-12 | Siemens Aktiengesellschaft | Method for interactively segmenting structures in image data records and image processing unit for carrying out the method |
-
2018
- 2018-03-30 CN CN201810298411.1A patent/CN108537803B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103678847A (en) * | 2012-09-24 | 2014-03-26 | 邵雁 | Skin soft tissue expanding auxiliary system |
CN102945127A (en) * | 2012-11-05 | 2013-02-27 | 深圳市旭东数字医学影像技术有限公司 | Interactive method for volume rendering displaying and system thereof |
CN103049907A (en) * | 2012-12-11 | 2013-04-17 | 深圳市旭东数字医学影像技术有限公司 | Interactive image segmentation method |
CN106570868A (en) * | 2015-10-10 | 2017-04-19 | 中国科学院深圳先进技术研究院 | Three-dimensional liver tumor semi-automatic segmentation method based on level set |
Also Published As
Publication number | Publication date |
---|---|
CN108537803A (en) | 2018-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108537803B (en) | A kind of CT image interactive segmentation method and device | |
Kowalski et al. | Art-based rendering of fur, grass, and trees | |
CN109523631B (en) | Contour line generation method and device based on B spline fitting | |
US7412362B2 (en) | River modeling | |
CN103049907B (en) | Interactive image segmentation method | |
CN109685914A (en) | Cutting profile based on triangle grid model mends face algorithm automatically | |
GB2555698B (en) | Three-dimensional model manipulation and rendering | |
CN108595858A (en) | Lightweight advanced treating working method based on BIM | |
US20090322748A1 (en) | Methods,systems, and computer program products for GPU-based point radiation for interactive volume sculpting and segmentation | |
CN106887000A (en) | The gridding processing method and its system of medical image | |
CN108876868B (en) | Computerized generation of decorative designs by placing instances of simple shapes according to direction guidance | |
Rocha et al. | Decal-maps: Real-time layering of decals on surfaces for multivariate visualization | |
Günther et al. | Decoupled opacity optimization for points, lines and surfaces | |
CN116822160A (en) | Cartesian grid generation method, device, equipment and medium | |
Rojo et al. | Fourier opacity optimization for scalable exploration | |
Breen et al. | 3d scan-conversion of CSG models into distance, closest-point and colour volumes | |
Gagnon et al. | Dynamic lapped texture for fluid simulations | |
CN106373192B (en) | A kind of non-topological coherence three-dimensional grid block tracing algorithm | |
Borgeat et al. | Layered surface fluid simulation for surgical training | |
Schmitz et al. | Efficient and quality contouring algorithms on the GPU | |
Foster et al. | Pen-and-ink for BlobTree implicit models | |
CN108538370A (en) | A kind of illumination volume drawing output method and device | |
Schein et al. | Adaptive extraction and visualization of silhouette curves from volumetric datasets | |
Balázs et al. | Efficient trimmed NURBS tessellation | |
Tekumalla et al. | Reverse engineering point clouds to fit tensor product B-spline surfaces by blending local fits |
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 | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 100083 Beijing Haidian District Zhongguancun East Road 8 Dongsheng building AB block nine floor 917 units. Applicant after: Beijing Ling Ling Medical Technology Co., Ltd. Address before: 100083 B 917, Dongsheng building, 8 Zhongguancun East Road, Haidian District, Beijing. Applicant before: Beijing Ling Ling Medical Technology Co., Ltd. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |