CN109685763A - Three dimensional CT boniness block automatic division method, device, computer equipment and storage medium - Google Patents
Three dimensional CT boniness block automatic division method, device, computer equipment and storage medium Download PDFInfo
- Publication number
- CN109685763A CN109685763A CN201811339882.9A CN201811339882A CN109685763A CN 109685763 A CN109685763 A CN 109685763A CN 201811339882 A CN201811339882 A CN 201811339882A CN 109685763 A CN109685763 A CN 109685763A
- Authority
- CN
- China
- Prior art keywords
- images
- dimensional
- scribble
- boniness
- initial 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 230000011218 segmentation Effects 0.000 claims abstract description 97
- 210000000988 bone and bone Anatomy 0.000 claims abstract description 65
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 39
- 238000009877 rendering Methods 0.000 claims abstract description 31
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000005295 random walk Methods 0.000 claims description 10
- 238000000926 separation method Methods 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 5
- 239000003550 marker Substances 0.000 claims description 4
- 239000003814 drug Substances 0.000 abstract description 6
- 238000012800 visualization Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 230000002452 interceptive effect Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 206010017076 Fracture Diseases 0.000 description 2
- 239000003086 colorant Substances 0.000 description 2
- 238000007794 visualization technique Methods 0.000 description 2
- 208000010392 Bone Fractures Diseases 0.000 description 1
- 208000024779 Comminuted Fractures Diseases 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- 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
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30008—Bone
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The present invention is suitable for computer field, provides a kind of three dimensional CT boniness block automatic division method, device, computer equipment and storage medium, described method includes following steps: obtaining boniness block CT images and is simultaneously parsed into three-dimensional data；Initial segmentation and three-dimensional contour surface rendering are carried out to the CT images for being parsed into three-dimensional data according to preset mode；Scribble label is carried out to the CT images that rendering is completed；Calculating is split to the CT images using the initial segmentation and scribble label according to preset algorithm, obtains segmentation result.In embodiments of the present invention, by utilizing the above-mentioned three dimensional CT boniness block automatic division method using scribble label, the bone block or sclerite that are disposably accurately partitioned into muti-piece different parts automatically from the three-dimensional CT images of bone may be implemented, high degree reduces the time cost of human-edited in three-dimensional bone medicine reconstruction.
Description
Technical field
The invention belongs to computer field more particularly to a kind of three dimensional CT boniness block automatic division methods, device, computer
Equipment and storage medium.
Background technique
Bone CT images refer to the CT scan for bone, it is X-ray beam, the γ using Accurate collimation
Ray, ultrasonic wave etc. surround the section that a certain position of human body is made one by one together with the detector high with sensitivity and sweep
It retouches, obtains this position bone CT images.
Due to the energy power limit of image documentation equipment acquisition information, when there are many boundary in three-dimensional CT images between bone and bone
It waits and unintelligible, if you need to divide and extract muti-piece bone from one group of CT images, generally requires a large amount of user's interaction, and
It can only be disposably partitioned into one piece to two pieces bone, efficiency is very low；On the other hand, for the feelings of fracture especially comminuted fracture
The quantity of shape, sclerite is often very much, and the anatomical structure that do not fix, therefore fractures without priori knowledge to assist doctor to do
The reference of block segmentation；The interactive mode of traditional medical software is all to be edited on 2d, and intuitive is very poor, it is difficult to
Different sclerites is distinguished, the mistake of segmentation result can be brought.
The 3-D image of bone CT images is split therefore, it is necessary to invent a kind of method, by all small bone blocks
Fast Segmentation comes out.
Summary of the invention
The embodiment of the present invention provides a kind of three dimensional CT boniness block automatic division method, by all in bone CT images
Bone block carries out scribble label and three-dimensional image segmentation calculates, and all small bone block Fast Segmentations are come out.
The embodiment of the invention provides a kind of three dimensional CT boniness block automatic division methods, and described method includes following steps:
It obtains boniness block CT images and is parsed into three-dimensional data；
Initial segmentation is carried out to the CT images for being parsed into three-dimensional data according to preset mode and three-dimensional is equivalent
Face rendering；
Scribble label is carried out to the CT images that rendering is completed；
Calculating is split to the CT images using the initial segmentation and scribble label according to preset algorithm, is obtained
Segmentation result.
Three dimensional CT boniness block automatic cutting device, which is characterized in that described device includes:
CT images acquiring unit, for obtaining boniness block CT images and being parsed into three-dimensional data；
Initial segmentation and rendering unit, for according to preset mode to the CT images for being parsed into three-dimensional data
Carry out initial segmentation and three-dimensional contour surface rendering；
Scribble marking unit, for carrying out scribble label to the CT images that rendering is completed；
Separation calculation unit, for using the initial segmentation and applying crow-style label to the CT shadow according to preset algorithm
As being split calculating, segmentation result is obtained.
In conclusion may be implemented by using the above-mentioned three dimensional CT boniness block automatic division method using scribble label
Disposably accurately it is partitioned into the bone block or sclerite of muti-piece different parts, high degree automatically from the three-dimensional CT images of bone
Reduce the time cost of human-edited in three-dimensional bone medicine reconstruction.
Detailed description of the invention
Fig. 1 is a kind of flow chart of three dimensional CT boniness block automatic division method provided in an embodiment of the present invention；
Fig. 2 is a kind of three dimensional CT boniness block automatic division method schematic diagram provided in an embodiment of the present invention；
Fig. 3 is the flow chart of another three dimensional CT boniness block automatic division method provided in an embodiment of the present invention；
Fig. 4 is the flow chart of another three dimensional CT boniness block automatic division method provided in an embodiment of the present invention；
Fig. 5 is a kind of structural block diagram of three dimensional CT boniness block automatic cutting device provided in an embodiment of the present invention；
Fig. 6 is the structural block diagram of another three dimensional CT boniness block automatic cutting device provided in an embodiment of the present invention；
Fig. 7 is the structural block diagram of another three dimensional CT boniness block automatic cutting device provided in an embodiment of the present invention.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to the accompanying drawings and embodiments, right
The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used in the restriction present invention.
The term used in embodiments of the present invention is only to be not intended to be limiting merely for for the purpose of describing particular embodiments
The present invention.In the embodiment of the present invention and the "an" of singular used in the attached claims, " described " and "the"
It is also intended to including most forms, unless the context clearly indicates other meaning.It is also understood that term used herein
"and/or" refers to and includes that one or more associated any or all of project listed may combine.
The present invention is applied to the computer equipment with image processing function, including but not limited to industrial computer, individual
Computer (Personal Computer, PC) etc..
Fig. 1 shows a kind of implementation process of three dimensional CT boniness block automatic division method provided in an embodiment of the present invention, in detail
It states as follows:
Step S101 obtains boniness block CT images and is parsed into three-dimensional data.
In embodiments of the present invention, the CT images comprising boniness block are obtained, the CT images are three dimensional CT sequence, such as Fig. 2-
Shown in 1., the boniness block CT images are parsed into three-dimensional data using image viewing technology, and be reconstructed into 3-D graphic.
The inventive point of described image visualization technique non-present invention, can be using any computer medical image visualization existing at present
Technology, herein with no restrictions.
Step S102, the CT images for being parsed into three-dimensional data to described according to preset mode carry out initial segmentation and
Three-dimensional contour surface rendering.
In embodiments of the present invention, it needs to being parsed into the CT images of three-dimensional data according to preset initial segmentation side
Method is split pretreatment, and the bone block in the CT images is made to generate initial segmentation, forms that several are specific, have uniqueness
The image-region of property, each bone block distinguished in the CT image that can be preliminary are convenient for subsequent processing, as Fig. 2-2. shown in.
Preferably, the preset initial segmentation method, which can be Threshold segmentation, region segmentation etc., can generate initial pictures
The method of segmentation, herein with no restrictions.
In embodiments of the present invention, the CT image that initial segmentation is completed substantially can distinguish all bone portions,
In order to apparent, the bone portion divided is distinguished and more intuitively convenient for subsequent operation, it is also necessary to according to preset
Algorithm carries out three-dimensional contour surface rendering to the CT image that initial segmentation is completed, and realizes that three-dimensional visualization shows the segmentation
Effect, as Fig. 2-3. shown in.
Preferably, the preset algorithm can be biggest advantage of light track algorithm.
Step S103 carries out scribble label to the CT images that rendering is completed.
In embodiments of the present invention, the bone block that the CT images that rendering forms three-dimensional visualization segmentation effect are completed is carried out
Corresponding interactive system can be used in scribble label, the scribble label.Scribble label is needed to every piece in the CT images
Bone block is marked, as Fig. 2-4. shown in.
Preferably, for the ease of distinguishing each bone block and subsequent calculating, the marker color of every piece of bone block is different, i.e., not
Same bone block carries out scribble label using different colors, forms apparent area to all bone blocks in the CT images with this
Point.
Preferably, in embodiments of the present invention, it in order to carry out scribble label to all bone blocks in the CT images, needs
There is a corresponding interactive operation interface, the scribble label can be to carry out on the two-dimension picture of the CT images, can also be with
It is to be carried out on the three-dimensional contour surface of the CT images, it can also the two combination progress.
Step S104 is split the CT images using the initial segmentation and scribble label according to preset algorithm
It calculates, obtains segmentation result.
In embodiments of the present invention, on the basis of the CT images have obtained initial segmentation, in order to will further own
Small bone block quickly distinguishes, it is also necessary to by initial segmentation and scribble label as input parameter, according still further to preset calculation
Method is calculated, and is obtained final segmentation result and is shown, as Fig. 2-5. shown in.
Preferably, the preset algorithm can be Random Walks algorithm.
Below by taking two labels as an example, calculated using Random Walks algorithm:
Remember p1,p2,...,pnFor all pixels in initial segmentation.Remember xiBelong to the probability of prospect, y for ith pixeliFor
Ith pixel belongs to the probability of background.λi,j=exp (- β | | pi-pj| |) it is piAnd pjBelong to same class (prospect or background it is general
Rate).
According to total probability formula:
Subscript ijJ-th of neighborhood of ith pixel is expressed,
If i-th point is prospect, then x in scribble labeli=1,
If i-th point is background, then x in scribble labeli=0.
Then according to above three condition, x can be uniquely solvedi, similar to solve yiBy comparing xiWith yiSize,
Determine that ith pixel belongs to prospect or background.
It here only need to be multiple labels by two tag extensions, the final size of the calculated result of more each label, really
Determine ith pixel and belongs to which label completes algorithm.
In conclusion may be implemented by using the above-mentioned three dimensional CT boniness block automatic division method using scribble label
Disposably accurately it is partitioned into the bone block or sclerite of muti-piece different parts, high degree automatically from the three-dimensional CT images of bone
Reduce the time cost of human-edited in three-dimensional bone medicine reconstruction.
Fig. 3 is shown suitable for another three dimensional CT boniness block automatic division method process provided in an embodiment of the present invention
Figure, compared to Figure 1, step S102 specifically includes step S201 and S202, specific as follows:
Step S201 carries out initial segmentation to the CT images according to preset mode.
In embodiments of the present invention, it needs to being parsed into the CT images of three-dimensional data according to preset initial segmentation side
Method is split pretreatment, and the bone block in the CT images is made to generate initial segmentation, forms that several are specific, have uniqueness
The image-region of property, each bone block distinguished in the CT image that can be preliminary are convenient for subsequent processing, as Fig. 2-2. shown in.
Preferably, the preset initial segmentation method, which can be Threshold segmentation, region segmentation etc., can generate initial pictures
The method of segmentation, herein with no restrictions.
Step S202 carries out three-dimensional contour surface rendering to the CT images according to preset method.
In embodiments of the present invention, the CT image that initial segmentation is completed substantially can distinguish all bone portions,
In order to apparent, the bone portion divided is distinguished and more intuitively convenient for subsequent operation, it is also necessary to according to preset
Algorithm carries out three-dimensional contour surface rendering to the CT image that initial segmentation is completed, and realizes that three-dimensional visualization shows the segmentation
Effect, as Fig. 2-3. shown in.
Preferably, the preset algorithm can be biggest advantage of light track algorithm.
Fig. 4 is shown suitable for another three dimensional CT boniness block automatic division method process provided in an embodiment of the present invention
Figure, compared to Figure 1, step S104 specifically includes step S301 and S302, specific as follows:
Step S301, by the initial segmentation and scribble label as input parameter.
In embodiments of the present invention, in order to further by all small bone blocks of the CT images for having obtained initial segmentation
Quickly distinguish, it is also necessary to the data of the CT images are calculated, it will be described in initial segmentation and scribble label conduct
The input parameter of calculating, is calculated according still further to preset algorithm.
Step S302 is split meter to the CT images using Random Walks algorithm according to the input parameter
It calculates.
In embodiments of the present invention, the initial segmentation and scribble label are utilized into Random as input parameter
Walks algorithm calculates final segmentation result and is shown, as Fig. 2-5. shown in
Below by taking two labels as an example, calculated using Random Walks algorithm:
Remember p1,p2,...,pnFor all pixels in initial segmentation.Remember xiBelong to the probability of prospect, y for ith pixeliFor
Ith pixel belongs to the probability of background.λi,j=exp (- β | | pi-pj| |) it is piAnd pjBelong to same class (prospect or background it is general
Rate).
According to total probability formula:
Subscript ijJ-th of neighborhood of ith pixel is expressed,
If i-th point is prospect, then x in scribble labeli=1,
If i-th point is background, then x in scribble labeli=0.
Then according to above three condition, x can be uniquely solvedi, similar to solve yiBy comparing xiWith yiSize,
Determine that ith pixel belongs to prospect or background.
It here only need to be multiple labels by two tag extensions, the final size of the calculated result of more each label, really
Determine ith pixel and belongs to which label completes algorithm.
By above-mentioned calculation method, quickly the bone block of the scribble label can be divided, disposable accurately automatic point
Cut out the bone block or sclerite of muti-piece different parts, high degree reduce three-dimensional bone medicine rebuild in human-edited time at
This.
Fig. 5 is shown suitable for a kind of structural frames of three dimensional CT boniness block automatic cutting device provided in an embodiment of the present invention
Figure, described device include:
CT images acquiring unit 401, for obtaining boniness block CT images and being parsed into three-dimensional data.
In embodiments of the present invention, the CT images comprising boniness block are obtained, the CT images are three dimensional CT sequence, such as Fig. 2-
Shown in 1., the boniness block CT images are parsed into three-dimensional data using image viewing technology, and be reconstructed into 3-D graphic.
The inventive point of described image visualization technique non-present invention, can be using any computer medical image visualization existing at present
Technology, herein with no restrictions.
Initial segmentation and rendering unit 402, for according to preset mode to the CT for being parsed into three-dimensional data
Image carries out initial segmentation and three-dimensional contour surface rendering.
In embodiments of the present invention, it needs to being parsed into the CT images of three-dimensional data according to preset initial segmentation side
Method is split pretreatment, and the bone block in the CT images is made to generate initial segmentation, forms that several are specific, have uniqueness
The image-region of property, each bone block distinguished in the CT image that can be preliminary are convenient for subsequent processing, as Fig. 2-2. shown in.
Preferably, the preset initial segmentation method, which can be Threshold segmentation, region segmentation etc., can generate initial pictures
The method of segmentation, herein with no restrictions.
In embodiments of the present invention, the CT image that initial segmentation is completed substantially can distinguish all bone portions,
In order to apparent, the bone portion divided is distinguished and more intuitively convenient for subsequent operation, it is also necessary to according to preset
Algorithm carries out three-dimensional contour surface rendering to the CT image that initial segmentation is completed, and realizes that three-dimensional visualization shows the segmentation
Effect, as Fig. 2-3. shown in.
Preferably, the preset algorithm can be biggest advantage of light track algorithm.
Scribble marking unit 403, for carrying out scribble label to the CT images that rendering is completed.
In embodiments of the present invention, the bone block that the CT images that rendering forms three-dimensional visualization segmentation effect are completed is carried out
Corresponding interactive system can be used in scribble label, the scribble label.Scribble label is needed to every piece in the CT images
Bone block is marked, as Fig. 2-4. shown in.
Preferably, for the ease of distinguishing each bone block and subsequent calculating, the marker color of every piece of bone block is different, i.e., not
Same bone block carries out scribble label using different colors, forms apparent area to all bone blocks in the CT images with this
Point.
Preferably, in embodiments of the present invention, it in order to carry out scribble label to all bone blocks in the CT images, needs
There is a corresponding interactive operation interface, the scribble label can be to carry out on the two-dimension picture of the CT images, can also be with
It is to be carried out on the three-dimensional contour surface of the CT images, it can also the two combination progress.
Separation calculation unit 404, for being marked using the initial segmentation and scribble to the CT according to preset algorithm
Image is split calculating, obtains segmentation result.
In embodiments of the present invention, on the basis of the CT images have obtained initial segmentation, in order to will further own
Small bone block quickly distinguishes, it is also necessary to by initial segmentation and scribble label as input parameter, according still further to preset calculation
Method is calculated, and is obtained final segmentation result and is shown, as Fig. 2-5. shown in.
Preferably, the preset algorithm can be Random Walks algorithm.
Below by taking two labels as an example, calculated using Random Walks algorithm:
Remember p1,p2,...,pnFor all pixels in initial segmentation.Remember xiBelong to the probability of prospect, y for ith pixeliFor
Ith pixel belongs to the probability of background.λi,j=exp (- β | | pi-pj| |) it is piAnd pjBelong to same class (prospect or background it is general
Rate).
According to total probability formula:
Subscript ijJ-th of neighborhood of ith pixel is expressed,
If i-th point is prospect, then x in scribble labeli=1,
If i-th point is background, then x in scribble labeli=0.
Then according to above three condition, x can be uniquely solvedi, similar to solve yiBy comparing xiWith yiSize,
Determine that ith pixel belongs to prospect or background.
It here only need to be multiple labels by two tag extensions, the final size of the calculated result of more each label, really
Determine ith pixel and belongs to which label completes algorithm.
In conclusion may be implemented by using the above-mentioned three dimensional CT boniness block automatic division method using scribble label
Disposably accurately it is partitioned into the bone block or sclerite of muti-piece different parts, high degree automatically from the three-dimensional CT images of bone
Reduce the time cost of human-edited in three-dimensional bone medicine reconstruction.
Fig. 6 shows the structure suitable for another three dimensional CT boniness block automatic cutting device provided in an embodiment of the present invention
Block diagram, the initial segmentation and rendering unit 402 are specific further include:
Initial segmentation module 501, for carrying out initial segmentation to the CT images according to preset mode.
In embodiments of the present invention, it needs to being parsed into the CT images of three-dimensional data according to preset initial segmentation side
Method is split pretreatment, and the bone block in the CT images is made to generate initial segmentation, forms that several are specific, have uniqueness
The image-region of property, each bone block distinguished in the CT image that can be preliminary are convenient for subsequent processing, as Fig. 2-2. shown in.
Preferably, the preset initial segmentation method, which can be Threshold segmentation, region segmentation etc., can generate initial pictures
The method of segmentation, herein with no restrictions.
Rendering module 502, for carrying out three-dimensional contour surface rendering to the CT images according to preset method.
In embodiments of the present invention, the CT image that initial segmentation is completed substantially can distinguish all bone portions,
In order to apparent, the bone portion divided is distinguished and more intuitively convenient for subsequent operation, it is also necessary to according to preset
Algorithm carries out three-dimensional contour surface rendering to the CT image that initial segmentation is completed, and realizes that three-dimensional visualization shows the segmentation
Effect, as Fig. 2-3. shown in.
Preferably, the preset algorithm can be biggest advantage of light track algorithm.
Fig. 7 shows the structure suitable for another three dimensional CT boniness block automatic cutting device provided in an embodiment of the present invention
Block diagram, the separation calculation unit 404 are specific further include:
Parameter acquisition module 601, for regarding the initial segmentation and scribble label as input parameter.
In embodiments of the present invention, in order to further by all small bone blocks of the CT images for having obtained initial segmentation
Quickly distinguish, it is also necessary to the data of the CT images are calculated, it will be described in initial segmentation and scribble label conduct
The input parameter of calculating, is calculated according still further to preset algorithm.
Separation calculation module 602 is used for according to the input parameter, using Random Walks algorithm to the CT images
It is split calculating.
In embodiments of the present invention, the initial segmentation and scribble label are utilized into Random as input parameter
Walks algorithm calculates final segmentation result and is shown, as Fig. 2-5. shown in
Below by taking two labels as an example, calculated using Random Walks algorithm:
Remember p1,p2,...,pnFor all pixels in initial segmentation.Remember xiBelong to the probability of prospect, y for ith pixeliFor
Ith pixel belongs to the probability of background.λi,j=exp (- β | | pi-pj| |) it is piAnd pjBelong to same class (prospect or background it is general
Rate).
According to total probability formula:
Subscript ijJ-th of neighborhood of ith pixel is expressed,
If i-th point is prospect, then x in scribble labeli=1,
If i-th point is background, then x in scribble labeli=0.
Then according to above three condition, x can be uniquely solvedi, similar to solve yiBy comparing xiWith yiSize,
Determine that ith pixel belongs to prospect or background.
It here only need to be multiple labels by two tag extensions, the final size of the calculated result of more each label, really
Determine ith pixel and belongs to which label completes algorithm.
By above-mentioned calculation method, quickly the bone block of the scribble label can be divided, disposable accurately automatic point
Cut out the bone block or sclerite of muti-piece different parts, high degree reduce three-dimensional bone medicine rebuild in human-edited time at
This.
In one embodiment it is proposed that a kind of computer equipment, the computer equipment include memory, processor and
It is stored in the computer program that can be run on the memory and on the processor, the processor executes the computer
The step of any one of present invention method is realized when program.
In one embodiment, a kind of computer readable storage medium is provided, is stored on computer readable storage medium
Computer program, when computer program is executed by processor, so that the step of processor executes any one of present invention method.
Although should be understood that various embodiments of the present invention flow chart in each step according to arrow instruction successively
It has been shown that, but these steps are not that the inevitable sequence according to arrow instruction successively executes.Unless expressly state otherwise herein,
There is no stringent sequences to limit for the execution of these steps, these steps can execute in other order.Moreover, each embodiment
In at least part step may include that perhaps these sub-steps of multiple stages or stage are not necessarily multiple sub-steps
Completion is executed in synchronization, but can be executed at different times, the execution in these sub-steps or stage sequence is not yet
Necessarily successively carry out, but can be at least part of the sub-step or stage of other steps or other steps in turn
Or it alternately executes.
Those of ordinary skill in the art will appreciate that realizing all or part of the process in above-described embodiment method, being can be with
Relevant hardware is instructed to complete by computer program, the program can be stored in a non-volatile computer and can be read
In storage medium, the program is when being executed, it may include such as the process of the embodiment of above-mentioned each method.Wherein, provided herein
Each embodiment used in any reference to memory, storage, database or other media, may each comprise non-volatile
And/or volatile memory.Nonvolatile memory may include that read-only memory (ROM), programming ROM (PROM), electricity can be compiled
Journey ROM (EPROM), electrically erasable ROM (EEPROM) or flash memory.Volatile memory may include random access memory
(RAM) or external cache.By way of illustration and not limitation, RAM is available in many forms, such as static state RAM
(SRAM), dynamic ram (DRAM), synchronous dram (SDRAM), double data rate sdram (DDRSDRAM), enhanced SDRAM
(ESDRAM), synchronization link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) directly RAM (RDRAM), straight
Connect memory bus dynamic ram (DRDRAM) and memory bus dynamic ram (RDRAM) etc..
Each technical characteristic of embodiment described above can be combined arbitrarily, for simplicity of description, not to above-mentioned reality
It applies all possible combination of each technical characteristic in example to be all described, as long as however, the combination of these technical characteristics is not deposited
In contradiction, all should be considered as described in this specification.
The embodiments described above only express several embodiments of the present invention, and the description thereof is more specific and detailed, but simultaneously
Limitations on the scope of the patent of the present invention therefore cannot be interpreted as.It should be pointed out that for those of ordinary skill in the art
For, without departing from the inventive concept of the premise, various modifications and improvements can be made, these belong to guarantor of the invention
Protect range.Therefore, the scope of protection of the patent of the invention shall be subject to the appended claims.
Claims (10)
1. a kind of three dimensional CT boniness block automatic division method, which is characterized in that described method includes following steps:
It obtains boniness block CT images and is parsed into three-dimensional data；
Initial segmentation and three-dimensional contour surface wash with watercolours are carried out to the CT images for being parsed into three-dimensional data according to preset mode
Dye；
Scribble label is carried out to the CT images that rendering is completed；
Calculating is split to the CT images using the initial segmentation and scribble label according to preset algorithm, is divided
As a result.
2. three dimensional CT boniness block automatic division method according to claim 1, which is characterized in that described according to preset side
The step of formula carries out initial segmentation and three-dimensional contour surface rendering to the CT images for being parsed into three-dimensional data, it is specific to wrap
It includes:
Initial segmentation is carried out to the CT images according to preset mode；
Three-dimensional contour surface rendering is carried out to the CT images according to preset method.
3. three dimensional CT boniness block automatic division method according to claim 1, which is characterized in that described pair is completed rendering
CT images carry out scribble label, specifically:
Scribble label is carried out on the two-dimension picture of the CT images；And/or
Scribble label is carried out on the three-dimensional contour surface of the CT images；
Wherein, the scribble marker color of each piece of bone is all different in the CT images.
4. three dimensional CT boniness block automatic division method according to claim 1, which is characterized in that described according to preset calculation
The step of method is split calculating to the CT images using the initial segmentation and scribble label, specifically includes:
By the initial segmentation and scribble label as input parameter；
According to the input parameter, calculating is split to the CT images using Random Walks algorithm.
5. a kind of three dimensional CT boniness block automatic cutting device, which is characterized in that described device includes:
CT images acquiring unit, for obtaining boniness block CT images and being parsed into three-dimensional data；
Initial segmentation and rendering unit, for being carried out according to preset mode to the CT images for being parsed into three-dimensional data
Initial segmentation and three-dimensional contour surface rendering；
Scribble marking unit, for carrying out scribble label to the CT images that rendering is completed；
Separation calculation unit, for according to preset algorithm using the initial segmentation and apply crow-style label to the CT images into
Row separation calculation, obtains segmentation result.
6. three dimensional CT boniness block automatic cutting device according to claim 5, which is characterized in that the initial segmentation and wash with watercolours
Dye unit specifically includes:
Initial segmentation module, for carrying out initial segmentation to the CT images according to preset mode；
Rendering module, for carrying out three-dimensional contour surface rendering to the CT images according to preset method.
7. three dimensional CT boniness block automatic cutting device according to claim 5, which is characterized in that the scribble marking unit
It specifically includes:
Two-dimension picture scribble module, for carrying out scribble label on the two-dimension picture of the CT images；And/or
Three-dimensional contour surface scribble module, for carrying out scribble label on the three-dimensional contour surface of the CT images；
Wherein, the scribble marker color of each piece of bone is all different in the CT images.
8. three dimensional CT boniness block automatic cutting device according to claim 5, which is characterized in that the separation calculation unit
Further include:
Parameter acquisition module, for regarding the initial segmentation and scribble label as input parameter；
Separation calculation module, for being divided the CT images using Random Walks algorithm according to the input parameter
Cut calculating.
9. a kind of computer equipment, which is characterized in that the computer equipment includes processor, and the processor is deposited for executing
The calculating formula program stored in reservoir realizes the step of method according to any of claims 1-4.
10. a kind of computer readable storage medium is stored thereon with computer program (instruction), which is characterized in that the calculating
When machine program (instruction) is executed by processor, realize method according to any of claims 1-4 the step of.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811339882.9A CN109685763A (en) | 2018-11-12 | 2018-11-12 | Three dimensional CT boniness block automatic division method, device, computer equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811339882.9A CN109685763A (en) | 2018-11-12 | 2018-11-12 | Three dimensional CT boniness block automatic division method, device, computer equipment and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109685763A true CN109685763A (en) | 2019-04-26 |
Family
ID=66185741
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811339882.9A Pending CN109685763A (en) | 2018-11-12 | 2018-11-12 | Three dimensional CT boniness block automatic division method, device, computer equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109685763A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111325745A (en) * | 2020-03-09 | 2020-06-23 | 北京深睿博联科技有限责任公司 | Fracture region analysis method and device, electronic device and readable storage medium |
| CN113506241A (en) * | 2021-05-25 | 2021-10-15 | 首都医科大学附属北京友谊医院 | Method for processing images of the ossicular chain and related product |
| CN115131388A (en) * | 2022-03-03 | 2022-09-30 | 中国人民解放军总医院第四医学中心 | Extraction method, device and equipment for bone mass directional superposition calculation |
| WO2023066364A1 (en) * | 2021-10-21 | 2023-04-27 | 苏州微创畅行机器人有限公司 | Three-dimensional image processing method and apparatus, computer device, and storage medium |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070025616A1 (en) * | 2005-08-01 | 2007-02-01 | Leo Grady | Editing of presegemented images/volumes with the multilabel random walker or graph cut segmentations |
| EP1750226A2 (en) * | 2005-08-01 | 2007-02-07 | Siemens Medical Solutions USA, Inc. | Editing of presegmented images/volumes with the multilabel random walker or graph cut segmentations |
| TW201121509A (en) * | 2009-12-30 | 2011-07-01 | Pou Yuen Technology Co Ltd | Segmentation method of multiple bone regions. |
| US20140093153A1 (en) * | 2012-09-28 | 2014-04-03 | Siemens Corporation | Method and System for Bone Segmentation and Landmark Detection for Joint Replacement Surgery |
| CN104969260A (en) * | 2012-09-27 | 2015-10-07 | 西门子制品生命周期管理软件股份有限公司 | Multi-bone segmentation for 3D computed tomography |
| CN107818566A (en) * | 2017-10-27 | 2018-03-20 | 中山大学 | A kind of image partition method based on partial structurtes information around random walk combination pixel |
-
2018
- 2018-11-12 CN CN201811339882.9A patent/CN109685763A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070025616A1 (en) * | 2005-08-01 | 2007-02-01 | Leo Grady | Editing of presegemented images/volumes with the multilabel random walker or graph cut segmentations |
| EP1750226A2 (en) * | 2005-08-01 | 2007-02-07 | Siemens Medical Solutions USA, Inc. | Editing of presegmented images/volumes with the multilabel random walker or graph cut segmentations |
| TW201121509A (en) * | 2009-12-30 | 2011-07-01 | Pou Yuen Technology Co Ltd | Segmentation method of multiple bone regions. |
| CN104969260A (en) * | 2012-09-27 | 2015-10-07 | 西门子制品生命周期管理软件股份有限公司 | Multi-bone segmentation for 3D computed tomography |
| US20140093153A1 (en) * | 2012-09-28 | 2014-04-03 | Siemens Corporation | Method and System for Bone Segmentation and Landmark Detection for Joint Replacement Surgery |
| CN107818566A (en) * | 2017-10-27 | 2018-03-20 | 中山大学 | A kind of image partition method based on partial structurtes information around random walk combination pixel |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111325745A (en) * | 2020-03-09 | 2020-06-23 | 北京深睿博联科技有限责任公司 | Fracture region analysis method and device, electronic device and readable storage medium |
| CN111325745B (en) * | 2020-03-09 | 2023-08-25 | 北京深睿博联科技有限责任公司 | Fracture region analysis method and device, electronic equipment and readable storage medium |
| CN113506241A (en) * | 2021-05-25 | 2021-10-15 | 首都医科大学附属北京友谊医院 | Method for processing images of the ossicular chain and related product |
| WO2023066364A1 (en) * | 2021-10-21 | 2023-04-27 | 苏州微创畅行机器人有限公司 | Three-dimensional image processing method and apparatus, computer device, and storage medium |
| CN115131388A (en) * | 2022-03-03 | 2022-09-30 | 中国人民解放军总医院第四医学中心 | Extraction method, device and equipment for bone mass directional superposition calculation |
| CN115131388B (en) * | 2022-03-03 | 2023-09-01 | 中国人民解放军总医院第四医学中心 | Extraction method, device and equipment for bone quantity directional superposition calculation |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kasten et al. | End-to-end convolutional neural network for 3D reconstruction of knee bones from bi-planar X-ray images | |
| US8929635B2 (en) | Method and system for tooth segmentation in dental images | |
| Folmsbee et al. | Active deep learning: Improved training efficiency of convolutional neural networks for tissue classification in oral cavity cancer | |
| US11132792B2 (en) | Cross domain medical image segmentation | |
| US8761493B2 (en) | Method and system for tooth segmentation in dental images | |
| CN102246204B (en) | Devices and methods for processing images using scale space | |
| CN109685763A (en) | Three dimensional CT boniness block automatic division method, device, computer equipment and storage medium | |
| CN108010021A (en) | A kind of magic magiscan and method | |
| US20130101197A1 (en) | Method and apparatus for generating an enhanced image from medical imaging data | |
| EP2263210B1 (en) | Radiotherapeutic treatment planning apparatus and methods | |
| EP2810217B1 (en) | Graph cuts-based interactive segmentation of teeth in 3-d ct volumetric data | |
| US11551362B2 (en) | Automatic segmentation process of a 3D medical image by several neural networks through structured convolution according to the geometry of the 3D medical image | |
| Sestini et al. | Fun-sis: A fully unsupervised approach for surgical instrument segmentation | |
| CN111317498A (en) | Image reconstruction method and device, equipment and storage medium | |
| Xia et al. | Adversarial pseudo healthy synthesis needs pathology factorization | |
| Li et al. | Synthetic-to-real domain adaptation joint spatial feature transform for stereo matching | |
| Liu et al. | Co-Net: A collaborative region-contour-driven network for fine-to-finer medical image segmentation | |
| CN110992439A (en) | Fiber bundle tracking method, computer device and storage medium | |
| Jin et al. | CWAN: Self-supervised learning for deep grape disease image composition | |
| CN111862001B (en) | Semi-automatic labeling method and device for CT images, electronic equipment and storage medium | |
| CN110738639B (en) | Medical image detection result display method, device, equipment and storage medium | |
| KR102241312B1 (en) | Apparatus and method for displaying consecutive nodule images automatically based on machine learning | |
| Lee et al. | Multi‐layer structural wound synthesis on 3D face | |
| EP2734147B1 (en) | Method for segmentation of dental images | |
| CN112767314B (en) | Medical image processing method, device, equipment and storage medium |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190426 |