CN109124662A - Rib cage center line detecting device and method - Google Patents
Rib cage center line detecting device and method Download PDFInfo
- Publication number
- CN109124662A CN109124662A CN201810770823.0A CN201810770823A CN109124662A CN 109124662 A CN109124662 A CN 109124662A CN 201810770823 A CN201810770823 A CN 201810770823A CN 109124662 A CN109124662 A CN 109124662A
- Authority
- CN
- China
- Prior art keywords
- rib cage
- region
- rib
- center line
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
Abstract
The invention belongs to technical field of image processing, specifically provide a kind of rib cage center line detecting device and method, it is intended to solve prior art poor robustness and the slow technical problem of calculating speed.For this purpose, the present invention provides a kind of rib cage center line detecting methods, comprising: based on rib cage region in preset rib cage identification model recognition target image, and obtain the corresponding rib region of every rib cage according to rib cage region；The corresponding multiple specific planes of each rib region in the corresponding three-dimension space image of target image are obtained, and calculate the center of gravity of corresponding rib region according to specific plane, all centers of gravity of each rib region are connected and obtain the center line of corresponding rib cage.Based on above-mentioned steps, rib cage center line detecting method provided by the invention can accurately and efficiently identify rib cage region in target image, be adapted to a variety of abnormal conditions, robustness is good, and calculating speed is fast.The device of the invention equally has above-mentioned beneficial effect.
Description
Technical field
The invention belongs to technical field of image processing, and in particular to a kind of rib cage center line detecting device and method.
Background technique
The doctor of image department would generally detect the rib cage transfer of fracture of rib and cancer by CT image.Each patient's
The quantity of CT image is likely to be breached over one hundred, and in order to check lesion present in CT image and exception, doctor is generally required successively
CT image is checked, dynamic change of the every rib cage on piece of different layers is tracked, visually to track rib cage
The change of sectional area.But manually CT image check it is relatively time consuming, and manually check be easy to cause fault and miss
Rib cage is abnormal.
Extract rib cage center line can be used to enhance expansion rib cage visualization, can allow doctor more efficiently and effectively
Check that the correlation circumstance of rib cage, the center line of rib cage can be also used for positioning organ, registration lesion and guide interval change point
The reference of correspondence between the serial thorax CT scan of analysis.
The detection method of existing rib cage center line may include following several methods:
1, rib cage is modeled as elongated tubular construction, and intrinsic using hessian matrix (Hessian) or structure tensor
Network analysis carries out the detection of back voxel, detects tubulose knot by calculating the characteristic value of the hessian matrix of each point in 3d space
Structure.But this method spent in terms of calculating it is very high, and possibly can not all obtain consistent result to all patients.
2, carry out rib cage center lines matching by template: the section of rib cage usually has certain rule, rib cage center
Spinal cord is partially dark, and the bone of rib edge is partially bright, and rib edge is oval, can match rib cage according to above-mentioned regular design template
And then center line is found, further according to needing to be modified center line.But in many cases, rib cage marrow may be than rib cage side
Boundary is darker, so that back voxel consistently can not be detected as rib cage central point.
3, center line is identified by the method for machine learning: extracts the Haar feature of CT image, passes through the number manually marked
According to the prediction model of training centre's line, to each pixel prediction, it is the probability of center line.
4, pass through the method tracing center line of curve tracing: artificial setting seed automatically selects seed, using such as card
The method for tracing of Kalman Filtering is tracked detected rib cage central point from a slice to next slice.But
This method needs artificial initial seed point, On Local Fuzzy caused by the rib cage lesion for fracture etc or discontinuous height
It is sensitive.
The robustness of prior art is poor, and for medical image, especially the data comprising lesion often change more
Sample, the rule being manually arranged cannot handle abnormal conditions well；Calculating speed is slow, due to needing to handle more complicated exception
Situation, the method for the prior art needs to be arranged more rule and step, slow so as to cause calculating speed.
Therefore, how to propose that a kind of technical solution to solve the above problems is that those skilled in the art need to solve at present
Problem.
Summary of the invention
In order to solve the above problem of the prior art, in order to solve, prior art poor robustness and calculating speed are slow to be asked
Topic, the first aspect of the present invention provide a kind of rib cage center line detecting device, comprising:
Rib cage identification module is configured to rib cage region in preset rib cage identification model recognition target image,
And the corresponding rib region of every rib cage is obtained according to the rib cage region；
Rib cage center line obtains module, is configured to obtain each institute in the corresponding three-dimension space image of the target image
The corresponding multiple specific planes of rib region are stated, and calculate the center of gravity of corresponding rib region according to the specific plane, by institute
All centers of gravity for stating each rib region, which are connected, obtains the center line of corresponding rib cage；
Wherein, the rib cage identification model is based on the Bone images demarcated in advance and using constructed by machine learning algorithm
Neural network model；The specific plane is that Plane Angle is the plane preset polar angle and be parallel to preset reference axis,
The angle value of the polar angle depends on the center of gravity according to the rib cage region to two-dimensional projection's point in corresponding rib region
Carry out the two-dimensional projection's point pole angular region obtained after polar coordinate transform.
In the optimal technical scheme of above-mentioned apparatus, the rib cage center line acquisition module is further configured to execute as follows
Operation:
The center of gravity of the rib cage region and each rib region are projected into two-dimensional surface, respectively obtain center of gravity projection
The rib cage view field of point and each rib region；
Polar coordinate transform is carried out to two-dimensional projection's point in each rib cage view field centered on the center of gravity subpoint；
It is corresponding that corresponding rib region is obtained according to the polar coordinates of each two-dimensional projection's point in each rib cage view field
Two-dimensional projection's point pole angular region.
In the optimal technical scheme of above-mentioned apparatus, the rib cage center line obtains module and is additionally configured to execute following behaviour
Make:
According to each corresponding two-dimensional projection's point pole angular region of rib region and the extraction of preset angle step
The corresponding multiple polar angles of each rib region.
In the optimal technical scheme of above-mentioned apparatus, the rib cage center line detecting device further includes center line amendment mould
Block, the center line correction module are configured that
Center line acquired in module is obtained to the rib cage center line using preset curve smoothing method to carry out smoothly
Processing.
In the optimal technical scheme of above-mentioned apparatus, the curve smoothing method is based on preset markov random file
The curve processing method of model.
In the optimal technical scheme of above-mentioned apparatus, the center line correction module is further configured to shown according to the following formula
Method center line is smoothed:
Wherein, the E (X) indicates that the energy function of the Markov random field model, the X indicate the center line
Some dimension (x, y or the dimension where z-axis) where upper all points, xiIt indicates to tie up in X for i-th point on the center line
The coordinate of degree, NiIndicate the number of point adjacent with i-th point on the center line, xjIndicate on the center line with xiPoint phase
J-th point adjacent of coordinate, x 'iIndicate the i-th point of initial coordinate on X-dimension, γ indicates x 'iShared weight, d
Indicate preset distance.
In the optimal technical scheme of above-mentioned apparatus, the rib cage identification module is further configured to perform the following operations:
Obtain multiple connections that the rib cage region described in the corresponding three-dimension space image of the target image includes
Region, and the corresponding rib region of every rib cage is obtained according to each connected region.
In the optimal technical scheme of above-mentioned apparatus, the rib cage identification module is additionally configured to perform the following operations:
Obtain the three-dimensional space point that all Euclidean distances in the rib cage region are less than pre-determined distance；
The connected region of every rib cage is obtained according to acquired three-dimensional space point.
The second aspect of the present invention provides a kind of rib cage center line detecting method, comprising:
Based on rib cage region in preset rib cage identification model recognition target image, and according to where the rib cage
Region obtains the corresponding rib region of every rib cage；
It is corresponding multiple specific to obtain each rib region in the corresponding three-dimension space image of the target image
Plane, and the center of gravity of corresponding rib region is calculated according to the specific plane, by all centers of gravity of each rib region
It is connected and obtains the center line of corresponding rib cage；
Wherein, the rib cage identification model is based on the Bone images demarcated in advance and using constructed by machine learning algorithm
Neural network model；The specific plane is that Plane Angle is the plane preset polar angle and be parallel to preset reference axis,
The angle value of the polar angle depends on the center of gravity according to the rib cage region to two-dimensional projection's point in corresponding rib region
Carry out the two-dimensional projection's point pole angular region obtained after polar coordinate transform.
In the optimal technical scheme of the above method, " obtaining in the corresponding three-dimension space image of the target image
The corresponding multiple specific planes of each rib region, and calculate according to the specific plane weight of corresponding rib region
Before the step of heart ", the method also includes:
The center of gravity of the rib cage region and each rib region are projected into two-dimensional surface, respectively obtain center of gravity projection
The rib cage view field of point and each rib region；
Polar coordinate transform is carried out to two-dimensional projection's point in each rib cage view field centered on the center of gravity subpoint；
It is corresponding that corresponding rib region is obtained according to the polar coordinates of each two-dimensional projection's point in each rib cage view field
Two-dimensional projection's point pole angular region.
In the optimal technical scheme of the above method, " according to each two-dimensional projection in each rib cage view field
After the step of corresponding two-dimensional projection's point pole angular region of the corresponding rib region of polar coordinates acquisition of point ", the method also includes:
According to each corresponding two-dimensional projection's point pole angular region of rib region and the extraction of preset angle step
The corresponding multiple polar angles of each rib region.
In the optimal technical scheme of the above method, " all centers of gravity of each rib region is being connected and are obtaining phase
Answer the center line of rib cage " the step of after, the method also includes:
Center line acquired in module is obtained to the rib cage center line using preset curve smoothing method to carry out smoothly
Processing.
In the optimal technical scheme of the above method, the curve smoothing method is based on preset markov random file
The curve processing method of model.
In the optimal technical scheme of the above method, " the rib cage center line is obtained using preset curve smoothing method
Center line acquired in modulus block is smoothed " the step of include according to the following formula shown in method center line is carried out it is smooth
Processing:
Wherein, the E (X) indicates that the energy function of the Markov random field model, the X indicate the center line
Some dimension (x, y or the dimension where z-axis) where upper all points, xiIt indicates to tie up in X for i-th point on the center line
The coordinate of degree, NiIndicate the number of point adjacent with i-th point on the center line, xjIndicate on the center line with xiPoint phase
J-th point adjacent of coordinate, x 'iIndicate the i-th point of initial coordinate on X-dimension, γ indicates x 'iShared weight, d
Indicate preset distance.
In the optimal technical scheme of the above method, " the corresponding rib of every rib cage is obtained according to the rib cage region
The step of bone region " includes:
Obtain multiple connections that the rib cage region described in the corresponding three-dimension space image of the target image includes
Region, and the corresponding rib region of every rib cage is obtained according to each connected region.
In the optimal technical scheme of the above method, " institute in the corresponding three-dimension space image of the target image is obtained
State multiple connected regions that rib cage region includes " the step of include:
Obtain the three-dimensional space point that all Euclidean distances in the rib cage region are less than pre-determined distance；
The connected region of every rib cage is obtained according to acquired three-dimensional space point.
The third aspect of the present invention provides a kind of storage device, wherein be stored with a plurality of program, described program be suitable for by
Processor is loaded to execute rib cage center line detecting method described in any of the above embodiments.
The fourth aspect of the present invention provides a kind of control device, including processor and storage equipment, the storage equipment
Suitable for storing a plurality of program, described program is suitable for being loaded by the processor to execute rib cage center line described in any of the above embodiments
Detection method.
Compared with the immediate prior art, above-mentioned technical proposal is at least had the following beneficial effects:
1, rib cage region in target image, rib cage identification can accurately and efficiently be identified by rib cage identification model
Model is adapted to a variety of abnormal conditions in target image, and the rib cage region identified is more robust；
2, it obtains in target image passing through polar coordinate transform behind accurate rib cage region based on rib cage identification model
Rib region center line can be converted to the computational problem of two-dimensional surface by method from the computational problem of three-dimensional planar, simplify meter
Treatment process is calculated, calculating speed is improved.
Detailed description of the invention
Fig. 1 is the key step schematic diagram of an embodiment of the present invention rib cage center line detecting method；
Fig. 2 is the schematic diagram of the single rib cage connected region of an embodiment of the present invention；
Fig. 3 is the schematic diagram of the single rib cage center line of an embodiment of the present invention；
Fig. 4 is the schematic diagram of the revised center line of an embodiment of the present invention；
Fig. 5 is the primary structure schematic diagram of the rib cage center line detecting device of an embodiment of the present invention.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art
Every other embodiment obtained without making creative work, shall fall within the protection scope of the present invention
The preferred embodiment of the present invention described with reference to the accompanying drawings.It will be apparent to a skilled person that this
A little embodiments are used only for explaining technical principle of the invention, it is not intended that limit the scope of the invention.
Refering to attached drawing 1, Fig. 1 illustratively gives the key step of rib cage center line detecting method in the present embodiment.Such as
Shown in Fig. 1, rib cage center line detecting method includes the following steps: in the present embodiment
Step S101: based on rib cage region in preset rib cage identification model recognition target image.
In the embodiment of the present invention, rib cage identification model can be from the region where identifying rib cage in target image.Specifically
Ground, target image can be CT image, and rib cage identification model can be based on the Bone images demarcated in advance and utilize engineering
The neural network model of algorithm building is practised, which can use UNet-like network structure, which includes on 4
Sampling structure and 4 down-sampling structures.Target image is inputted into rib cage identification model, the convolutional layer of rib cage identification model can be adopted
Convolution operation is carried out with target image of the convolution kernel of 3*3 to input, extracts the feature of the target image of input convolutional layer；Rib cage
ReLU activation primitive can be used in identification model, removes the redundancy in image data, retains the feature of data most possibly；Rib
The upsample operation of the maxpooling and size=2 of stride=2 can be used in the pond layer of bone identification model, will be defeated
The characteristic pattern for entering pond layer reduces, and reduces the feature vector of convolutional layer output；The convolution kernel that 1*1 finally can be used carries out convolution
Operation is connected all features by full articulamentum, and the feature after connection is inputted softmax classifier, and target figure is exported
The rib cage region as in.
In order to rib cage identification model can region preferably in recognition target image where rib cage, can be by marking in advance
Fixed Bone images are trained rib cage identification model.Specifically, image translation, figure can be carried out to original target image
The operation such as image rotation turn, image scaling and Image Reversal, under the premise of not changing target image content, to the picture of target image
Element carries out geometric transformation.To image pixel carry out geometric transformation after, to after geometric transformation image carry out gamma transformation, with
The operations such as machine cutting, Gaussian Blur and normalization, so that rib cage identification model to be trained can preferably identify rib cage institute
In region.In the training process, rib cage identification model output to be trained can be judged using dice loss loss function
The degree of closeness of rib cage region and true rib cage region.Rib cage identification model after training can be accurately from mesh
The region where rib cage is identified in logo image, it is subsequent rapidly to be obtained in rib cage based on accurate rib cage region
Heart line improves the speed of operation.In addition, rib cage identification model is adapted to a variety of abnormal conditions, the rib cage institute identified
It is more robust in region.
The rib cage identification model of the embodiment of the present invention can also use VGG in addition to that can use UNet-like network structure
Network structure or ResNet network structure etc., the present invention is it is not limited here.
Step S102: multiple connected regions that rib cage region includes in the corresponding three-dimension space image of target image are obtained
Domain, and the corresponding rib region of every rib cage is obtained according to each connected region.
With continued reference to attached drawing 2, Fig. 2 illustratively gives the schematic diagram of single rib cage connected region.It is identified by rib cage
Multiple connected regions that region in model recognition target image where rib cage includes, the corresponding rib cage of each connected region.
In order to obtain the center line of every rib cage, every rib cage in the region where rib cage can be separated, to every rib cage list
Solely analyzed.Specifically, the three-dimensional space point in available rib cage region, and judge any two in three-dimensional space point
Whether the Euclidean distance of point is less than or equal to pre-determined distance (Euclidean distance of such as two points is 1), if so, can be determined that this two o'clock
It is connection, the point of several connections may be constructed the connected region of single rib cage.
Step S103: obtaining the corresponding multiple specific planes of each rib region, calculates corresponding rib cage according to specific plane
All centers of gravity of each rib region are connected and obtain the center line of corresponding rib cage by the center of gravity in region.
In embodiments of the present invention, after obtaining rib cage region, the first center of gravity of rib cage region is obtained, by rib cage
The first center of gravity and each rib region of region project to two-dimensional surface, respectively obtain the first center of gravity subpoint and each
First center of gravity subpoint is denoted as C by the rib cage view field of rib region, and two-dimensional surface can be x/y plane.By each rib cage
After region projections to two-dimensional surface, the starting point of rib cage view field and the angle of end point are obtained, α and β can be denoted as respectively,
Centered on the first center of gravity subpoint C, pole seat is carried out according to polar angle range (α, β) of the predetermined angle step-length to rib cage view field
Mark transformation, obtains the corresponding multiple polar angles of each rib region.Assuming that current angular is θ, predetermined angle step-length is ω, preset angle
Degree step-length shares 10 corresponding values in current angular θ, solves corresponding second center of gravity of this 10 values, and so on,
Available multiple second centers of gravity in polar angle range (α, β).
The corresponding multiple specific planes of each rib region in the corresponding three-dimension space image of target image are obtained, with current
For angle is θ, the plane parallel with z-axis is constructed based on corresponding second center of gravity of the angle, which is specific plane, should
The intersection point of plane and rib region is the z-axis coordinate of the second center of gravity, and multiple second centers of gravity of each rib region are mutually got continuously
To the center line of the rib cage.With continued reference to attached drawing 3, Fig. 3 illustratively gives the schematic diagram of single rib cage center line.Pass through pole
The method of coordinate transform can ask rib region center line from the calculating that the computational problem of three-dimensional planar is converted to two-dimensional surface
Topic, simplifies algorithm process process, improves calculating speed.
Step S104: center line is smoothed using preset curve smoothing method.
The above-mentioned center line being calculated is independently to calculate the separation of every rib cage, there is no in view of different rib cages it
Between correlation, certain points may not on centerline on the center line being finally calculated, therefore center line is not smooth.
In addition, the rib region of rib cage identification model identification is likely to occur deviation, it may appear that some abnormal points, these abnormal points can deviate
Normal position is farther out.In order to be modified to center line, it can use preset curve smoothing method and center line carried out smoothly
Processing.Specifically, curve smoothing method, which can be, is modified center line based on preset Markov random field model, can
To be modified according to method shown in formula (1) to center line:
Wherein, E (X) indicates that the energy function of Markov random field model, X indicate where point all on center line
Some dimension (x, y or the dimension where z-axis), xiIndicate the i-th point of coordinate in X-dimension on center line, NiIndicate center line
The number of the upper point adjacent with i-th point, xjIndicate center line on xiJ-th point adjacent of coordinate of point, x 'iIt indicates i-th
Initial coordinate of the point on X-dimension, γ indicate x 'iShared weight, d indicate preset distance.
The first part (part on the plus sige left side) of energy function requires being closer for adjacent central point, it is adjacent in
The number of heart point is to determine a factor of line smoothing degree；Second part (part on the right of the plus sige) requirement of energy function
It is too far that the distance of central point cannot deviate initial point, and d indicates preset distance, can be to avoid deviateing too far different of normal central line
Chang Dian is also only the distance of d at a distance from center line if deviateing too far point.It, can be with by the constraint of energy function
Keep center line more smooth and insensitive to abnormal point.With continued reference to attached drawing 4, Fig. 4 illustratively gives revised center
The schematic diagram of line.
Although each step is described in the way of above-mentioned precedence in above-described embodiment, this field
Technical staff is appreciated that the effect in order to realize the present embodiment, executes between different steps not necessarily in such order,
It (parallel) execution simultaneously or can be executed with reverse order, these simple variations all protection scope of the present invention it
It is interior.
Based on above-mentioned rib cage center line detecting embodiment of the method, the present invention also provides a kind of rib cage center line detecting dresses
It sets.With reference to the accompanying drawing, which is illustrated.
Refering to attached drawing 5, Fig. 5 illustratively gives the primary structure of rib cage center line detecting device in the present embodiment.Such as
Shown in Fig. 5, the present embodiment rib cage center line detecting device includes that rib cage identification module 1 and rib cage center line obtain module 2.
Rib cage identification module 1 is configured to rib cage region in preset rib cage identification model recognition target image,
And the corresponding rib region of every rib cage is obtained according to rib cage region；
Rib cage center line obtains module 2, is configured to obtain each rib cage in the corresponding three-dimension space image of target image
The corresponding multiple specific planes in region, and the center of gravity of corresponding rib region is calculated according to specific plane, by each rib region
All centers of gravity be connected obtain the center line of corresponding rib cage；
Wherein, rib cage identification model is based on the Bone images demarcated in advance and to utilize mind constructed by machine learning algorithm
Through network model；Specific plane is that Plane Angle is default polar angle and the plane for being parallel to preset reference axis, the angle of polar angle
Angle value depends on after carrying out polar coordinate transform to two-dimensional projection's point in corresponding rib region according to the center of gravity of rib cage region
Obtained two-dimensional projection's point pole angular region.
In a preferred embodiment of the embodiment of the present invention, rib cage center line obtains module 2 and is further configured to hold
The following operation of row:
The center of gravity of rib cage region and each rib region are projected into two-dimensional surface, respectively obtain center of gravity subpoint and
The rib cage view field of each rib region；
Polar coordinate transform is carried out to two-dimensional projection's point in each rib cage view field centered on center of gravity subpoint；
Corresponding rib region corresponding two is obtained according to the polar coordinates of two-dimensional projection's point each in each rib cage view field
Tie up subpoint polar angle range.
In a preferred embodiment of the embodiment of the present invention, rib cage center line obtains module 2 and is additionally configured to execute such as
Lower operation:
Each rib cage is extracted according to the corresponding two-dimensional projection's point pole angular region of each rib region and preset angle step
The corresponding multiple polar angles in region.
In a preferred embodiment of the embodiment of the present invention, rib cage center line detecting device further includes center line amendment
Module, center line correction module are configured that
Center line acquired in module 2 is obtained to rib cage center line using preset curve smoothing method smoothly to be located
Reason.
In a preferred embodiment of the embodiment of the present invention, curve smoothing method be based on preset Markov with
The curve processing method of airport model.
In a preferred embodiment of the embodiment of the present invention, center line correction module is further configured to according to the following formula
Shown in method center line is smoothed:
Wherein, E (X) indicates that the energy function of Markov random field model, X indicate where point all on center line
Some dimension (x, y or the dimension where z-axis), xiIndicate the i-th point of coordinate in X-dimension on center line, NiIndicate center line
The number of the upper point adjacent with i-th point, xjIndicate center line on xiJ-th point adjacent of coordinate of point, x 'iIt indicates i-th
Initial coordinate of the point on X-dimension, γ indicate x 'iShared weight, d indicate preset distance.
In a preferred embodiment of the embodiment of the present invention, rib cage identification module 1 is further configured to execute as follows
Operation:
Multiple connected regions that rib cage region includes in the corresponding three-dimension space image of target image are obtained, and
The corresponding rib region of every rib cage is obtained according to each connected region.
In a preferred embodiment of the embodiment of the present invention, rib cage identification module 1 is additionally configured to perform the following operations:
Obtain the three-dimensional space point that all Euclidean distances in rib cage region are less than pre-determined distance；
The connected region of every rib cage is obtained according to acquired three-dimensional space point.
Further, it is based on above method embodiment, the present invention also provides a kind of storage devices, wherein the storage device
It is stored with a plurality of program and these programs may be adapted to be loaded by processor to execute the rib cage center of above method embodiment
Line detecting method.
Still further, being based on above method embodiment, the present invention also provides a kind of control device, the control device packets
Include processor and storage equipment；Storage equipment may be adapted to store a plurality of program and these programs may be adapted to be added by processor
It carries to execute the rib cage center line detecting method of above method embodiment.
Person of ordinary skill in the field can be understood that for convenience of description and succinctly, the present invention fills
The specific work process and related description for setting embodiment, can refer to corresponding processes in the foregoing method embodiment, and with it is upper
Method beneficial effect having the same is stated, details are not described herein.
Those skilled in the art should be able to recognize that, side described in conjunction with the examples disclosed in the embodiments of the present disclosure
Method step, can be realized with electronic hardware, computer software, or a combination of the two, in order to clearly demonstrate electronic hardware and
The interchangeability of software generally describes each exemplary composition and step according to function in the above description.These
Function is executed actually with electronic hardware or software mode, specific application and design constraint depending on technical solution.
Those skilled in the art can use different methods to achieve the described function each specific application, but this reality
Now it should not be considered as beyond the scope of the present invention.
It should be noted that description and claims of this specification and term " first " in above-mentioned attached drawing, "
Two " etc. be to be used to distinguish similar objects, rather than be used to describe or indicate specific sequence or precedence.It should be understood that this
The data that sample uses can be interchanged in appropriate circumstances, so that the embodiment of the present invention described herein can be in addition at this
In illustrate or description those of other than sequence implement.
So far, it has been combined preferred embodiment shown in the drawings and describes technical solution of the present invention, still, this field
Technical staff is it is easily understood that protection scope of the present invention is expressly not limited to these specific embodiments.Without departing from this
Under the premise of the principle of invention, those skilled in the art can make equivalent change or replacement to the relevant technologies feature, these
Technical solution after change or replacement will fall within the scope of protection of the present invention.
Claims (18)
1. a kind of rib cage center line detecting device, characterized by comprising:
Rib cage identification module is configured to rib cage region in preset rib cage identification model recognition target image, and
The corresponding rib region of every rib cage is obtained according to the rib cage region；
Rib cage center line obtains module, is configured to obtain each rib in the corresponding three-dimension space image of the target image
The corresponding multiple specific planes in bone region, and the center of gravity of corresponding rib region is calculated according to the specific plane, it will be described every
All centers of gravity of a rib region, which are connected, obtains the center line of corresponding rib cage；
Wherein, the rib cage identification model is based on the Bone images demarcated in advance and to utilize mind constructed by machine learning algorithm
Through network model；The specific plane is that Plane Angle is the plane preset polar angle and be parallel to preset reference axis, described
The angle value of polar angle depends on carrying out two-dimensional projection's point in corresponding rib region according to the center of gravity of the rib cage region
The two-dimensional projection's point pole angular region obtained after polar coordinate transform.
2. rib cage center line detecting device according to claim 1, which is characterized in that the rib cage center line obtains module
It is further configured to perform the following operations:
The center of gravity of the rib cage region and each rib region are projected into two-dimensional surface, respectively obtain center of gravity subpoint and
The rib cage view field of each rib region；
Polar coordinate transform is carried out to two-dimensional projection's point in each rib cage view field centered on the center of gravity subpoint；
Corresponding rib region corresponding two is obtained according to the polar coordinates of each two-dimensional projection's point in each rib cage view field
Tie up subpoint polar angle range.
3. rib cage center line detecting device according to claim 2, which is characterized in that the rib cage center line obtains module
It is additionally configured to perform the following operations:
It is extracted according to each corresponding two-dimensional projection's point pole angular region of rib region and preset angle step described each
The corresponding multiple polar angles of rib region.
4. rib cage center line detecting device according to claim 1, which is characterized in that the rib cage center line detecting device
It further include center line correction module, the center line correction module is configured that
Center line acquired in module is obtained to the rib cage center line using preset curve smoothing method to be smoothed.
5. rib cage center line detecting device according to claim 4, which is characterized in that the curve smoothing method is to be based on
The curve processing method of preset Markov random field model.
6. rib cage center line detecting device according to claim 5, which is characterized in that
The center line correction module be further configured to according to the following formula shown in method center line is smoothed:
Wherein, the E (X) indicates that the energy function of the Markov random field model, the X indicate institute on the center line
Some dimension (x, y or the dimension where z-axis) where some points, xiIndicate on the center line i-th point in X-dimension
Coordinate, NiIndicate the number of point adjacent with i-th point on the center line, xjIndicate on the center line with xiPoint is adjacent
J-th point of coordinate, x 'iIndicate the i-th point of initial coordinate on X-dimension, γ indicates x 'iShared weight, d are indicated
Preset distance.
7. rib cage center line detecting device according to claim 1 to 6, which is characterized in that the rib cage identification
Module is further configured to perform the following operations:
Multiple connected regions that the rib cage region described in the corresponding three-dimension space image of the target image includes are obtained,
And the corresponding rib region of every rib cage is obtained according to each connected region.
8. rib cage center line detecting device according to claim 7, which is characterized in that the rib cage identification module also configures
To perform the following operations:
Obtain the three-dimensional space point that all Euclidean distances in the rib cage region are less than pre-determined distance；
The connected region of every rib cage is obtained according to acquired three-dimensional space point.
9. a kind of rib cage center line detecting method, characterized by comprising:
Based on rib cage region in preset rib cage identification model recognition target image, and according to the rib cage region
Obtain the corresponding rib region of every rib cage；
The corresponding multiple specific planes of each rib region in the corresponding three-dimension space image of the target image are obtained,
And the center of gravity that corresponding rib region is calculated according to the specific plane mutually gets all centers of gravity of each rib region continuously
To the center line of corresponding rib cage；
Wherein, the rib cage identification model is based on the Bone images demarcated in advance and to utilize mind constructed by machine learning algorithm
Through network model；The specific plane is that Plane Angle is the plane preset polar angle and be parallel to preset reference axis, described
The angle value of polar angle depends on carrying out two-dimensional projection's point in corresponding rib region according to the center of gravity of the rib cage region
The two-dimensional projection's point pole angular region obtained after polar coordinate transform.
10. rib cage center line detecting method according to claim 9, which is characterized in that " obtaining in the target image
The corresponding multiple specific planes of each rib region in corresponding three-dimension space image, and according to the specific plane meter
Calculate the center of gravity of corresponding rib region " the step of before, the method also includes:
The center of gravity of the rib cage region and each rib region are projected into two-dimensional surface, respectively obtain center of gravity subpoint and
The rib cage view field of each rib region；
Polar coordinate transform is carried out to two-dimensional projection's point in each rib cage view field centered on the center of gravity subpoint；
Corresponding rib region corresponding two is obtained according to the polar coordinates of each two-dimensional projection's point in each rib cage view field
Tie up subpoint polar angle range.
11. rib cage center line detecting method according to claim 10, which is characterized in that " according to each rib cage
The step of the corresponding two-dimensional projection's point pole angular region of the corresponding rib region of polar coordinates acquisition of each two-dimensional projection's point in view field "
After rapid, the method also includes:
It is extracted according to each corresponding two-dimensional projection's point pole angular region of rib region and preset angle step described each
The corresponding multiple polar angles of rib region.
12. rib cage center line detecting method according to claim 9, which is characterized in that " by each rib region
All centers of gravity be connected obtain the center line of corresponding rib cage " the step of after, the method also includes:
Center line acquired in module is obtained to the rib cage center line using preset curve smoothing method to be smoothed.
13. rib cage center line detecting method according to claim 12, which is characterized in that the curve smoothing method is base
In the curve processing method of preset Markov random field model.
14. rib cage center line detecting method according to claim 13, which is characterized in that " utilize preset curve smoothing
Method to the rib cage center line obtain module acquired in center line be smoothed " the step of include according to the following formula shown in
Method center line is smoothed:
Wherein, the E (X) indicates that the energy function of the Markov random field model, the X indicate institute on the center line
Some dimension (x, y or the dimension where z-axis) where some points, xiIndicate on the center line i-th point in X-dimension
Coordinate, NiIndicate the number of point adjacent with i-th point on the center line, xjIndicate on the center line with xiPoint is adjacent
J-th point of coordinate, x 'iIndicate the i-th point of initial coordinate on X-dimension, γ indicates x 'iShared weight, d are indicated
Preset distance.
15. the rib cage center line detecting method according to any one of claim 9-13, which is characterized in that " according to described
Rib cage region obtains the corresponding rib region of every rib cage " the step of include:
Multiple connected regions that the rib cage region described in the corresponding three-dimension space image of the target image includes are obtained,
And the corresponding rib region of every rib cage is obtained according to each connected region.
16. rib cage center line detecting method according to claim 15, which is characterized in that " obtain in the target image
The step of multiple connected regions that rib cage region described in corresponding three-dimension space image includes " includes:
Obtain the three-dimensional space point that all Euclidean distances in the rib cage region are less than pre-determined distance；
The connected region of every rib cage is obtained according to acquired three-dimensional space point.
17. a kind of storage device, wherein being stored with a plurality of program, which is characterized in that described program be suitable for by processor load with
Perform claim requires rib cage center line detecting method described in any one of 9-16.
18. a kind of control device, including processor and storage equipment, the storage equipment are suitable for storing a plurality of program, feature
It is, described program, which is suitable for being loaded as the processor, requires rib cage center line inspection described in any one of 9-16 with perform claim
Survey method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810770823.0A CN109124662B (en) | 2018-07-13 | 2018-07-13 | Rib center line detection device and method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810770823.0A CN109124662B (en) | 2018-07-13 | 2018-07-13 | Rib center line detection device and method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109124662A true CN109124662A (en) | 2019-01-04 |
| CN109124662B CN109124662B (en) | 2021-12-03 |
Family
ID=64800670
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810770823.0A Active CN109124662B (en) | 2018-07-13 | 2018-07-13 | Rib center line detection device and method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109124662B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109830289A (en) * | 2019-01-18 | 2019-05-31 | 上海皓桦科技股份有限公司 | Bone images display device |
| CN110807770A (en) * | 2019-10-30 | 2020-02-18 | 杭州依图医疗技术有限公司 | Medical image processing, recognizing and displaying method and storage medium |
| CN111080573A (en) * | 2019-11-19 | 2020-04-28 | 上海联影智能医疗科技有限公司 | Rib image detection method, computer device and storage medium |
| CN111985512A (en) * | 2019-05-22 | 2020-11-24 | 杭州依图医疗技术有限公司 | Method and device for establishing human body rib model |
| CN111986138A (en) * | 2019-05-22 | 2020-11-24 | 杭州依图医疗技术有限公司 | Method and device for obtaining rib positioning |
| CN112862793A (en) * | 2021-02-20 | 2021-05-28 | 北京医准智能科技有限公司 | Rib center line extraction method and device based on rib three-dimensional shape and distribution |
| CN113643176A (en) * | 2021-07-28 | 2021-11-12 | 沈阳先进医疗设备技术孵化中心有限公司 | Rib display method and device |
| CN114638828A (en) * | 2022-05-18 | 2022-06-17 | 数聚(山东)医疗科技有限公司 | Radiological image intelligent segmentation method based on computer vision |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004167109A (en) * | 2002-11-21 | 2004-06-17 | Osaka Industrial Promotion Organization | Three-dimensional measuring method, three-dimensional measuring system, image processor, and computer program |
| US20050010107A1 (en) * | 2003-06-23 | 2005-01-13 | Hong Shen | Thoracic cage coordinate system for recording pathologies in lung CT volume data |
| US20070110295A1 (en) * | 2005-10-17 | 2007-05-17 | Siemens Corporate Research Inc | System and method for enhanced viewing of rib metastasis |
| US20070249910A1 (en) * | 2004-10-28 | 2007-10-25 | Kiraly Atilla P | 2D visualization for rib analysis |
| US20080137932A1 (en) * | 2006-12-11 | 2008-06-12 | Siemens Corporation Research, Inc. | System and Method for Automatic Detection of Rib Metastasis in Computed Tomography Volume |
| CN101452577A (en) * | 2008-11-26 | 2009-06-10 | 沈阳东软医疗***有限公司 | Rib auto-demarcating method and device |
| US20100177945A1 (en) * | 2009-01-09 | 2010-07-15 | Fujifilm Corporation | Image processing method, image processing apparatus, and image processing program |
| CN102567625A (en) * | 2010-10-29 | 2012-07-11 | 西门子公司 | Automated rib ordering and pairing |
| US20130077841A1 (en) * | 2011-09-27 | 2013-03-28 | Siemens Corporation | Method and System for Automatic Rib Centerline Extraction Using Learning Base Deformable Template Matching |
| CN103054602A (en) * | 2011-09-19 | 2013-04-24 | 西门子公司 | Method and system for up-vector detection for ribs in computed tomography volume |
| US20130150704A1 (en) * | 2011-12-12 | 2013-06-13 | Shuki Vitek | Magnetic resonance imaging methods for rib identification |
| CN104252708A (en) * | 2013-06-28 | 2014-12-31 | 深圳先进技术研究院 | X-ray chest radiographic image processing method and X-ray chest radiographic image processing system |
| CN106611411A (en) * | 2015-10-19 | 2017-05-03 | 上海联影医疗科技有限公司 | Rib segmentation method in medical image and medical image processing device |
| CN106683090A (en) * | 2016-12-31 | 2017-05-17 | 上海联影医疗科技有限公司 | Rib positioning method in medical image and system thereof |
| CN107154038A (en) * | 2016-04-22 | 2017-09-12 | 孔德兴 | A kind of visual fracture of rib aided diagnosis method of rib |
| CN107798682A (en) * | 2017-08-31 | 2018-03-13 | 深圳联影医疗科技有限公司 | Image segmentation system, method, apparatus and computer-readable recording medium |
-
2018
- 2018-07-13 CN CN201810770823.0A patent/CN109124662B/en active Active
Patent Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004167109A (en) * | 2002-11-21 | 2004-06-17 | Osaka Industrial Promotion Organization | Three-dimensional measuring method, three-dimensional measuring system, image processor, and computer program |
| US20050010107A1 (en) * | 2003-06-23 | 2005-01-13 | Hong Shen | Thoracic cage coordinate system for recording pathologies in lung CT volume data |
| US20070249910A1 (en) * | 2004-10-28 | 2007-10-25 | Kiraly Atilla P | 2D visualization for rib analysis |
| US20070110295A1 (en) * | 2005-10-17 | 2007-05-17 | Siemens Corporate Research Inc | System and method for enhanced viewing of rib metastasis |
| US20080137932A1 (en) * | 2006-12-11 | 2008-06-12 | Siemens Corporation Research, Inc. | System and Method for Automatic Detection of Rib Metastasis in Computed Tomography Volume |
| CN101452577A (en) * | 2008-11-26 | 2009-06-10 | 沈阳东软医疗***有限公司 | Rib auto-demarcating method and device |
| US20100177945A1 (en) * | 2009-01-09 | 2010-07-15 | Fujifilm Corporation | Image processing method, image processing apparatus, and image processing program |
| CN102567625A (en) * | 2010-10-29 | 2012-07-11 | 西门子公司 | Automated rib ordering and pairing |
| CN103054602A (en) * | 2011-09-19 | 2013-04-24 | 西门子公司 | Method and system for up-vector detection for ribs in computed tomography volume |
| US20130077841A1 (en) * | 2011-09-27 | 2013-03-28 | Siemens Corporation | Method and System for Automatic Rib Centerline Extraction Using Learning Base Deformable Template Matching |
| CN103218800A (en) * | 2011-09-27 | 2013-07-24 | 西门子公司 | Method and system for automatic rib centerline extraction using learning base deformable template matching |
| US20130150704A1 (en) * | 2011-12-12 | 2013-06-13 | Shuki Vitek | Magnetic resonance imaging methods for rib identification |
| CN104252708A (en) * | 2013-06-28 | 2014-12-31 | 深圳先进技术研究院 | X-ray chest radiographic image processing method and X-ray chest radiographic image processing system |
| CN106611411A (en) * | 2015-10-19 | 2017-05-03 | 上海联影医疗科技有限公司 | Rib segmentation method in medical image and medical image processing device |
| CN107154038A (en) * | 2016-04-22 | 2017-09-12 | 孔德兴 | A kind of visual fracture of rib aided diagnosis method of rib |
| CN106683090A (en) * | 2016-12-31 | 2017-05-17 | 上海联影医疗科技有限公司 | Rib positioning method in medical image and system thereof |
| CN107798682A (en) * | 2017-08-31 | 2018-03-13 | 深圳联影医疗科技有限公司 | Image segmentation system, method, apparatus and computer-readable recording medium |
| CN107808379A (en) * | 2017-08-31 | 2018-03-16 | 深圳联影医疗科技有限公司 | Image processing system, method, apparatus and computer-readable medium |
Non-Patent Citations (4)
| Title |
|---|
| JIAZHENG SHI等: "Pulmonary nodule registration in serial CT scans based on rib anatomy and nodule template matching", 《MEDICAL PHYSICS》 * |
| MARCO LOOG ET AL: "Segmentation of the Posterior Ribs in Chest Radiographs Using Iterated Contextual Pixel Classification", 《IEEE TRANSACTIONS ON MEDICAL IMAGING》 * |
| 刘驰等: "一种新的肝脏CTA图像肋骨和脊椎骨分割方法", 《电脑知识与技术》 * |
| 杨杰等: "骨扫描图像预处理", 《医学影像分析和三维重建及其应用》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109830289B (en) * | 2019-01-18 | 2021-04-06 | 上海皓桦科技股份有限公司 | Rib image display device |
| CN109830289A (en) * | 2019-01-18 | 2019-05-31 | 上海皓桦科技股份有限公司 | Bone images display device |
| CN111986138B (en) * | 2019-05-22 | 2023-08-25 | 杭州依图医疗技术有限公司 | Method and device for acquiring rib positioning |
| CN111985512A (en) * | 2019-05-22 | 2020-11-24 | 杭州依图医疗技术有限公司 | Method and device for establishing human body rib model |
| CN111986138A (en) * | 2019-05-22 | 2020-11-24 | 杭州依图医疗技术有限公司 | Method and device for obtaining rib positioning |
| CN111985512B (en) * | 2019-05-22 | 2024-02-02 | 杭州依图医疗技术有限公司 | Method and device for establishing human rib model |
| CN110807770A (en) * | 2019-10-30 | 2020-02-18 | 杭州依图医疗技术有限公司 | Medical image processing, recognizing and displaying method and storage medium |
| CN111080573A (en) * | 2019-11-19 | 2020-04-28 | 上海联影智能医疗科技有限公司 | Rib image detection method, computer device and storage medium |
| CN111080573B (en) * | 2019-11-19 | 2024-02-27 | 上海联影智能医疗科技有限公司 | Rib image detection method, computer device and storage medium |
| CN112862793A (en) * | 2021-02-20 | 2021-05-28 | 北京医准智能科技有限公司 | Rib center line extraction method and device based on rib three-dimensional shape and distribution |
| CN113643176A (en) * | 2021-07-28 | 2021-11-12 | 沈阳先进医疗设备技术孵化中心有限公司 | Rib display method and device |
| CN114638828B (en) * | 2022-05-18 | 2022-09-23 | 数聚(山东)医疗科技有限公司 | Radiological image intelligent segmentation method based on computer vision |
| CN114638828A (en) * | 2022-05-18 | 2022-06-17 | 数聚(山东)医疗科技有限公司 | Radiological image intelligent segmentation method based on computer vision |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109124662B (en) | 2021-12-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109124662A (en) | Rib cage center line detecting device and method | |
| US11576645B2 (en) | Systems and methods for scanning a patient in an imaging system | |
| JP6947759B2 (en) | Systems and methods for automatically detecting, locating, and semantic segmenting anatomical objects | |
| EP3195257B1 (en) | Systems and methods for segmenting medical images based on anatomical landmark-based features | |
| CN105518744B (en) | Pedestrian recognition methods and equipment again | |
| CN106611411B (en) | Method for segmenting ribs in medical image and medical image processing device | |
| CN106340015B (en) | A kind of localization method and device of key point | |
| USRE47609E1 (en) | System for detecting bone cancer metastases | |
| US20220157047A1 (en) | Feature Point Detection | |
| US20200268251A1 (en) | System and method for patient positioning | |
| CN108986891A (en) | Medical imaging processing method and processing device, electronic equipment and storage medium | |
| JP2023507109A (en) | Automated tumor identification and segmentation from medical images | |
| EP3016027A2 (en) | Human body part detection system and human body part detection method | |
| US11386553B2 (en) | Medical image data | |
| US9092867B2 (en) | Methods for segmenting images and detecting specific structures | |
| CN103294883A (en) | Method and system for intervention planning for transcatheter aortic valve implantation | |
| CN110415792A (en) | Image detecting method, device, computer equipment and storage medium | |
| CN114022554B (en) | Massage robot acupoint detection and positioning method based on YOLO | |
| US20220207771A1 (en) | Heart Position Estimation | |
| US9675245B2 (en) | Method and device for determining the eye torsion | |
| CN110634554A (en) | Spine image registration method | |
| CN114998374A (en) | Image segmentation method, device and equipment based on position prior and storage medium | |
| CN113706514A (en) | Focus positioning method, device and equipment based on template image and storage medium | |
| CN109087357A (en) | Scan orientation method, apparatus, computer equipment and computer readable storage medium | |
| CN110570929B (en) | Structure display method and device, computer 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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |