Summary of the invention
The technical problem to be solved in the present invention is, manually counts in radioscopic image in prior art
Calculate, be difficult to measure rachiocamposis angle, it is provided that a kind of method measuring rachiocamposis angle.
The present invention solves the technical scheme of the problems referred to above and there is provided a kind of method measuring rachiocamposis angle,
It is characterized in that, the method comprises the following steps:
S1, obtain spine image by 3-D supersonic imaging;
S2, in described spine image, obtain be positioned at supraspinal curve point;
S3, generate final spinal column curve according to described curve point;
S4, according to the flex point on spinal column curve, calculate rachiocamposis angle.
The present invention measure rachiocamposis angle method in, in described step S2, manually and/
Or computer-aid method obtains and is positioned at supraspinal curve point.
The present invention measure rachiocamposis angle method in, in described step S3, use 5 rank or with
On fitting of a polynomial generate final spinal column curve.
In the method measuring rachiocamposis angle of the present invention, manually obtain and be positioned at supraspinal curve
During point, the number of described curve point at least includes 10.
In the method measuring rachiocamposis angle of the present invention, obtained by computer-aid method and be positioned at ridge
Curve point in vertebra post, specifically includes following steps:
Described volume rendering is carried out histogram equalization, obtains the enhancing image of spinal column;
Use OTSU algorithm that described enhancing image is split, and obtain threshold value;
From the bottom of described enhancing image to top, detection is gone to be positioned on vertebra according to the direction by-line of vertical coordinate
Curve point;
Calculate meansigma methods and the standard deviation thereof of the abscissa of all curve point detected, if described curve point
When being unsatisfactory for certain condition, delete this curve point.
In the method measuring rachiocamposis angle of the present invention, remove detecting position according to the direction by-line of vertical coordinate
In supraspinal curve point, wherein: use 5 rank or curve that above polynomial curve goes matching to detect
Point.
In the method measuring rachiocamposis angle of the present invention, matching is gone to detect using polynomial curve
Curve point before, define the spinal column scope that this curve is corresponding, and calculate the distance value in the range of this, and
Relatively this distance value and described threshold value.
In the method measuring rachiocamposis angle of the present invention, in described step S4, spinal column angle of bend is
It is positioned on knee of curve maximum angle between tangent line.
Implement the method measuring rachiocamposis angle of the present invention so that utilize 3-D supersonic imaging to measure ridge
Post angle of bend becomes simple possible.Manually obtain curve point simultaneously, or computer-aid method obtains
Curve point, or two ways combination, and then improve the precision measuring rachiocamposis angle.
Detailed description of the invention
The method measuring rachiocamposis angle of the present invention, by the spine image obtaining 3-D supersonic imaging
Carry out relevant treatment, obtain supraspinal curve point, so form spinal column curve, and then can be according to spinal curve
Flex point on line, calculates rachiocamposis angle so that utilize three-dimensional ultrasound pattern to measure rachiocamposis angle
Degree becomes simple possible.Manually obtain curve point simultaneously, or computer-aid method obtain curve point,
Or two ways combines, meet the demand that difference uses the user of 3-D imaging system, and then ridge is measured in raising
The precision of post angle of bend.
In order to make the purpose of the present invention, technical scheme and advantage clearer, below in conjunction with accompanying drawing and reality
Execute example, the present invention is further elaborated.Only should be appreciated that specific embodiment described herein
In order to explain the present invention, it is not intended to limit the present invention.
As it is shown in figure 1, be the flow chart of the method measuring rachiocamposis angle of the present invention, the method includes
Following steps:
S110, obtain spine image by 3-D supersonic imaging;
In this step, the spine image obtained by three-dimension ultrasonic imaging system, and 3-D supersonic imaging system
System discloses, and does not repeats them here.
S120, in spine image, obtain be positioned at supraspinal curve point;
In this step, manually and/or computer-aid method obtain be positioned at supraspinal curve point.
Concrete step will be described in detail below, and does not repeats them here.
S130, generate final spinal column curve according to the curve point obtained;
In this step, generate spinal column curve and can use 5 rank or above fitting of a polynomial.
S140, according to the flex point on spinal column curve, calculate rachiocamposis angle.
In this step, by the coordinate of flex point, and combine arctan function, calculate tangent line in this flex point
Angle, rachiocamposis angle is for being positioned on knee of curve maximum angle between tangent line.
It is below several particular preferred embodiments of the method for the rachiocamposis angle of the present invention:
As in figure 2 it is shown, be the flow chart of the method measuring rachiocamposis angle of first embodiment of the invention,
In the present embodiment, the method comprises the following steps:
S210, obtain spine image by 3-D supersonic imaging;
In this step, the spine image obtained by three-dimension ultrasonic imaging system, and 3-D supersonic imaging system
System discloses, and does not repeats them here.
S220, in spine image, manually obtain and be positioned at supraspinal curve point;
In this step, medical personnel or the operator of three-dimension ultrasonic imaging system, according to the experience of self
The number choosing curve point at least includes 10, and makes Curve Point on vertebra.
S230, according to the curve point obtained, and use 5 rank or above fitting of a polynomial to generate final ridge
Post curve;
S240, according to the flex point on spinal column curve, calculate rachiocamposis angle.
In this step, flex point is also that the operator of medical personnel or three-dimension ultrasonic imaging system is according to self
Experience choose, by the coordinate of flex point, and combine arctan function, calculate tangent line in this flex point
Angle, rachiocamposis angle is for being positioned on knee of curve maximum angle between tangent line.
As shown in Figure 3 a-3b, it is the result figure measuring rachiocamposis angle of first embodiment of the invention, its
In, Fig. 3 a is the spine image of the 3-D supersonic imaging acquisition of first embodiment of the invention, the song manually obtained
The number of line point has 15, as shown in Figure 3 b, according to flex point A1 on spinal column curve and A2, and flex point A1
It is respectively as follows: 10.4 ° and 16.6 ° with the angle of tangent line on A2, then calculates rachiocamposis angle A 1-2.
As shown in Figure 4, it is the flow chart of the method measuring rachiocamposis angle of second embodiment of the invention,
With difference in first embodiment it is: obtained by computer-aid method and be positioned at supraspinal curve point,
Making to measure rachiocamposis angle and become automatic, in the present embodiment, the method comprises the following steps:
S410, obtain spine image by 3-D supersonic imaging;
In this step, by three-dimension ultrasonic imaging system acquisition spine image, and three-dimension ultrasonic imaging system
Disclose, do not repeated them here.
S420, in spine image, by computer-aid method obtain be positioned at supraspinal curve point;Tool
Body comprises the following steps:
S421, spine image is carried out histogram equalization, obtain spinal column enhancing image f (x, y);
In this step, spine image is also carried out medium filtering.
(x, y) splits, and obtains threshold value T to strengthening image f for S422, employing OTSU algorithmh;
S423, from strengthen image f (x, bottom y) to top, according to the direction by-line of vertical coordinate go detection ridge
Curve point S of postcurve(x, y), wherein: use 5 rank or song that above polynomial curve goes matching to detect
Line point;
In this step, specifically include:
A, the spinal column scope of definition straight line y, calculate the distance value dy in the range of this, and compare this distance value
Dy and threshold value ThIf this distance value dy is less than threshold value Th, then step c is performed;
B, the intensity distributions of use polynomial curve g (x) fitting a straight line y, meeting condition g ' (x)=0 He
G " the minimum curve point that (x) > 0 obtainsAs curve point, wherein, polynomial curve g (x) is 5 rank
More than or, " (x) represents the first derivative of g (x), second dervative respectively, and straight line y corresponds to this song for g ' (x), g
Line;
C, judge that whether straight line y is to strengthen image f (x, the last item straight line y) the most then perform step
S424, otherwise detects straight line y+1 according to above-mentioned steps.
S424, calculate all curve point detectedThe meansigma methods of abscissaAnd standard
Deviationsx, when being unsatisfactory for certain condition, delete curve pointWherein, condition is
S430, according to the curve point obtained, and use 5 rank or above fitting of a polynomial to generate final ridge
Post curve;
S440, according to the multiple flex points on spinal column curve, calculate rachiocamposis angle.
In this step, first determine the multiple flex points on spinal column curve, more respectively by the coordinate of each flex point,
And combine arctan function, and calculate the angle of tangent line in this flex point, wherein, rachiocamposis angle is for being positioned at
Maximum angle between tangent line on knee of curve.
As illustrated in figs. 5 a-5 c, it is the result figure measuring rachiocamposis angle of second embodiment of the invention, its
In, Fig. 5 a is the spine image obtained by 3-D supersonic imaging of second embodiment of the invention, with the present invention
The spine image obtained by 3-D supersonic imaging of first embodiment is identical, the song obtained by above-mentioned steps
Line point, as shown in Figure 5 b, according to the angle of tangent line in the flex point on spinal column curve, then calculates rachiocamposis
Angle A 1-2, as shown in Figure 5 c, wherein, in figure labelling flex point A1 and on A2 the angle of tangent line be respectively as follows:
17.0 ° and 16.5 °.
As shown in Figure 6, it is the flow chart of the method measuring rachiocamposis angle of third embodiment of the invention,
In the present embodiment, it is also have hands not only by computer-aid method with the second embodiment difference
What dynamic method obtained is positioned at supraspinal curve point so that measures rachiocamposis angle and becomes automanual, the party
Method comprises the following steps:
S610, obtain spine image by 3-D supersonic imaging;
In this step, the spine image obtained by three-dimension ultrasonic imaging system, and 3-D supersonic imaging system
System discloses, and does not repeats them here.
S620, in spine image, manually obtain and be positioned at supraspinal curve point, and pass through area of computer aided
Method obtains and is positioned at supraspinal curve point;Wherein: the number of the curve point manually obtained is any;Passing through
Computer-aid method obtains in the concrete step being positioned at supraspinal curve point, deletes and is unsatisfactory for condition
During curve point, the curve point manually obtained being included, the curve point even manually obtained is unsatisfactory for accordingly
During condition, deleting corresponding curve point, concrete step is described in detail the most, at this
Repeat no more.
S630, according to the curve point obtained, and use medium filtering and 5 rank or above fitting of a polynomial to generate
Final spinal column curve;
S640, according to the multiple flex points on spinal column curve, calculate rachiocamposis angle.
In this step, first determine the multiple flex points on spinal column curve, more respectively by the coordinate of each flex point,
And combine arctan function, and calculate the angle of tangent line in this flex point, wherein, rachiocamposis angle is for being positioned at
Maximum angle between tangent line on knee of curve.
As shown in Fig. 7 a-7c, it is the result figure measuring rachiocamposis angle of second embodiment of the invention, its
In, Fig. 7 a is the spine image of the acquisition by 3-D supersonic imaging of second embodiment of the invention, with this
The spine image obtained by 3-D supersonic imaging of bright first embodiment and the second embodiment is identical, by upper
State the curve point that step obtains, as shown in Figure 7b, according to the angle of tangent line in the flex point on spinal column curve, then
Calculate rachiocamposis angle A 1-2, as shown in Figure 7 c, wherein, figure is cut on flex point A1 of labelling and A2
The angle of line is respectively as follows: 16.0 ° and 16.5 °.
Summary, the method measuring rachiocamposis angle of the present invention so that utilize 3-D supersonic imaging to measure
Rachiocamposis angle becomes simple possible.Manually obtain curve point simultaneously, or computer-aid method obtains
Take curve point, or two ways combines, and then improve the precision measuring rachiocamposis angle.
The above, the only present invention preferably detailed description of the invention, but protection scope of the present invention not office
Being limited to this, any those familiar with the art, can be easily in technical scope disclosed by the invention
The change expected or replacement, all should contain within protection scope of the present invention.Therefore, the protection of the present invention
Scope should be as the criterion with scope of the claims.