CN101991420B - Method and system for positioning suspected blood vessel narrow area - Google Patents
Method and system for positioning suspected blood vessel narrow area Download PDFInfo
- Publication number
- CN101991420B CN101991420B CN 200910169381 CN200910169381A CN101991420B CN 101991420 B CN101991420 B CN 101991420B CN 200910169381 CN200910169381 CN 200910169381 CN 200910169381 A CN200910169381 A CN 200910169381A CN 101991420 B CN101991420 B CN 101991420B
- Authority
- CN
- China
- Prior art keywords
- point
- blood vessel
- suspected
- slope
- curve
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The invention discloses a method for positioning a suspected blood vessel narrow area. The method comprises the following steps of: acquiring an image of a blood vessel to be examined; carrying out curve fitting on the diameter of the cross section of the blood vessel according to the image of the blood vessel to be examined; extracting the low point of the fitted curve as a suspected blood vessel narrow point; and determining the suspected blood vessel narrow area where the suspected blood vessel narrow point is located. Meanwhile, the invention also discloses a system for positioning the suspected blood vessel narrow area. The method and the system of the invention can be used for increasing the positioning speed and accuracy of the suspected blood vessel narrow area.
Description
Technical field
The present invention relates to Medical Imaging, particularly a kind of localization method of suspected blood vessel narrow area and system.
Background technology
Medical Imaging provides a kind of new noninvasive test means, medical imaging device, for example NMR (Nuclear Magnetic Resonance)-imaging (MRI) scanner, computer X-ray fault imaging (CT) scanner, the image in examine zone is scanned and obtains in patient's examine zone, and then the doctor can analyze image.
In clinical, sometimes need patient's suspected blood vessel narrow area is positioned, then the doctor further judges suspected blood vessel narrow area, thereby determines real angiostenosis zone.In the prior art, the localization method of suspected blood vessel narrow area is: at first, medical imaging device scans the examine blood vessel, and obtains the image of examine blood vessel; Secondly, the doctor observes the shape of blood vessel in the image of examine blood vessel, and according to experience location suspected blood vessel narrow area, then relies on manual mode to make a mark at suspected blood vessel narrow area, for follow-up judgement.
This shows, in the prior art, the location of suspected blood vessel narrow area is that the personal experience according to the doctor finishes, this might occur because doctors experience is not enough slow to the locating speed of suspected blood vessel narrow area, even also might omit suspected blood vessel narrow area, thereby locating speed and the precision of suspected blood vessel narrow area have been reduced.
Summary of the invention
The invention provides the localization method in a kind of angiostenosis zone, to improve locating speed and the precision of suspected blood vessel narrow area.
The present invention also provides the navigation system in a kind of angiostenosis zone, to improve locating speed and the precision of suspected blood vessel narrow area.
A kind of localization method of suspected blood vessel narrow area, the method comprises:
Obtain the image of examine blood vessel;
Image according to the examine blood vessel carries out curve fitting to the vessel cross-sections diameter;
The low spot of the curve of extraction institute match is as the suspected blood vessel narrow point;
Determine the suspected blood vessel narrow area at described suspected blood vessel narrow point place.
The method further comprises: the image at the examine blood vessel indicates suspected blood vessel narrow area.
Further comprise: the curve to institute's match carries out smothing filtering.
Described curve fitting comprises:
For each vessel cross-sections in the image of examine blood vessel, obtain the position of this vessel cross-sections on the blood vessel bearing of trend, and calculate the diameter of this vessel cross-sections;
With the position of vessel cross-sections as the diameter of X coordinate, vessel cross-sections as the Y coordinate, fit to curve.
The low spot of the curve of described extraction institute match comprises:
The first derivative of every bit and second dervative on the calculated curve, with first derivative be 0 and second dervative be that 0 point is as flex point;
N adjacent point of flex point left side on the trade-off curve, N adjacent point of flex point right side on the trade-off curve, wherein, N is positive integer, is preferably more than 3 and less than 10 positive integer;
Whether the Y coordinate of judging flex point less than the Y coordinate of each point in selected 2N the point, if, then with this flex point as low spot.
The suspected blood vessel narrow area at described definite suspected blood vessel narrow point place comprises:
Determine left margin and the right margin of the suspected blood vessel narrow area at suspected blood vessel narrow point place.
The left margin of determining the suspected blood vessel narrow area at suspected blood vessel narrow point place comprises: the adjacent point of choosing suspected blood vessel narrow point left side on the described curve; Calculate the tangent slope of selected point, if the absolute value of this slope is less than or equal to a preset value, the corresponding vessel position of X coordinate that then will put is as the left margin of suspected blood vessel narrow area; Otherwise, choose next point left, then repeat step and the subsequent step thereof of above-mentioned slope calculations;
The right margin of determining the suspected blood vessel narrow area at suspected blood vessel narrow point place comprises: the adjacent point of choosing suspected blood vessel narrow point right side on the described curve; Calculate the tangent slope of selected point, if this slope is less than or equal to preset value, the corresponding vessel position of X coordinate that then will put is as the right margin of suspected blood vessel narrow area; Otherwise, choose next point to the right, then repeat step and the subsequent step thereof of above-mentioned slope calculations.
The left margin of determining the suspected blood vessel narrow area at suspected blood vessel narrow point place comprises: the adjacent point of choosing suspected blood vessel narrow point left side on the described curve; Calculate the tangent slope of selected point, if the absolute value of the slope of a continuous N point is less than or equal to a preset value, the corresponding vessel position of X coordinate of the point (preferably first point in this M point) in then this M being put is as the left margin of suspected blood vessel narrow area; Otherwise, choose again next point left, then repeat step and the subsequent step thereof of above-mentioned slope calculations;
The right margin of determining the suspected blood vessel narrow area at suspected blood vessel narrow point place comprises: the adjacent point of choosing suspected blood vessel narrow point right side on the described curve; Calculate the tangent slope of selected point, if the absolute value of the slope of a continuous N point is less than or equal to a preset value, the corresponding vessel position of X coordinate of the point (preferably first point in this M point) in then this M being put is as the right margin of suspected blood vessel narrow area; Otherwise, choose again next point to the right, then repeat step and the subsequent step thereof of above-mentioned slope calculations; Wherein, M is positive integer, is preferably more than 1 and less than 5 positive integer.
The present invention also provides a kind of navigation system of suspected blood vessel narrow area, and this system comprises: an image acquisition unit, a match unit, a low spot extraction unit and a narrow zone determining unit; Wherein,
Described image acquisition unit is for the image that obtains the examine blood vessel;
Described match unit is used for according to the image of examine blood vessel the vessel cross-sections diameter being carried out curve fitting;
Described low spot extraction unit be used for to extract the low spot of curve of institute's match as the suspected blood vessel narrow point;
Described narrow zone determining unit is for the suspected blood vessel narrow area of determining described suspected blood vessel narrow point place.
This system comprises that further indicates a unit,
Be used for indicating suspected blood vessel narrow area at the image of examine blood vessel.
This system further comprises a smothing filtering unit, is used for the curve of institute's match is carried out smothing filtering.
Described match unit comprises: a position acquisition subelement, a diameter computation subunit and a curve fitting subelement; Wherein,
Described position acquisition subelement is used for obtaining the position of each vessel cross-sections on the blood vessel bearing of trend;
Described diameter computation subunit is for the diameter that calculates each vessel cross-sections;
Described curve fitting subelement is used for the position of vessel cross-sections as the diameter of X coordinate, vessel cross-sections fitting to curve as the Y coordinate.
Described low spot extraction unit comprises: a flex point is determined subelement, a chooser unit and a judgment sub-unit; Wherein,
Described flex point is determined subelement, is used for first derivative and the second dervative of every bit on the calculated curve, with first derivative be 0 and second dervative be that 0 point is as flex point;
Described chooser unit is used for N adjacent point of flex point left side on the trade-off curve, N adjacent point of flex point right side on the trade-off curve, and wherein, N is positive integer, is preferably more than 3 and less than 10 positive integer;
Whether described judgment sub-unit, the Y coordinate that is used for judging flex point less than selected 2N the Y coordinate of putting each point, if, then with this flex point as low spot.
Described narrow zone determining unit comprises: left margin determines that subelement and a right margin determine subelement; Wherein,
Described left margin is determined subelement, is used for the left margin of the suspected blood vessel narrow area at definite suspected blood vessel narrow point place;
Described right margin is determined subelement, is used for the right margin of the suspected blood vessel narrow area at definite suspected blood vessel narrow point place.
Described left margin determines that subelement comprises: left consecutive points chooser module, left slope meter operator module and one get a submodule left; Wherein,
Described left consecutive points chooser module is used for choosing the adjacent point on the left of the suspected blood vessel narrow point on the described curve;
Described left slope meter operator module, for the tangent slope that calculates selected point, if the absolute value of this slope is less than or equal to a preset value, the corresponding vessel position of X coordinate that then will put is as the left margin of suspected blood vessel narrow area; Get a submodule left otherwise carry out;
A described submodule of getting left is used for choosing next point left, then carries out left slope meter operator module.
Described right margin determines that subelement comprises: right consecutive points chooser module, right slope meter operator module and one get a submodule to the right; Wherein,
Described right consecutive points chooser module is for the adjacent point of choosing suspected blood vessel narrow point right side on the described curve;
Described right slope meter operator module, for the tangent slope that calculates selected point, if the absolute value of this slope is less than or equal to a preset value, the corresponding vessel position of X coordinate that then will put is as the right margin of suspected blood vessel narrow area; Get a submodule to the right otherwise carry out;
A described submodule of getting to the right is used for choosing next point to the right, then carries out right slope meter operator module.
Described left margin determines that subelement comprises: left consecutive points chooser module, left slope meter operator module and one get a submodule left; Wherein,
Described left consecutive points chooser module is used for choosing the adjacent point on the left of the suspected blood vessel narrow point on the described curve;
Described left slope meter operator module, be used for calculating the tangent slope of selected point, if the absolute value of the slope of each point is less than or equal to a preset value in the continuous N point, the corresponding vessel position of X coordinate of the point (preferably first point in this M point) in then this M being put is as the left margin of suspected blood vessel narrow area; Get a submodule left otherwise carry out;
A described submodule of getting left is used for choosing next point left, then carries out left slope meter operator module;
Described right margin determines that subelement comprises: right consecutive points chooser module, right slope meter operator module and one get a submodule to the right; Wherein,
Described right consecutive points chooser module is for the adjacent point of choosing suspected blood vessel narrow point right side on the described curve;
Described right slope meter operator module, be used for calculating the tangent slope of selected point, if the absolute value of the slope of each point is less than or equal to a preset value in the continuous N point, the corresponding vessel position of X coordinate of the point (preferably first point in this M point) in then this M being put is as the right margin of suspected blood vessel narrow area; Get a submodule to the right otherwise carry out;
A described submodule of getting to the right is used for choosing next point to the right, then carries out right slope meter operator module; Wherein, M is positive integer, is preferably more than 1 and less than 5 positive integer.
As seen from the above technical solution, the low spot of the curve of suspected blood vessel narrow area localization method provided by the present invention and the match of extraction institute of system is as the angiostenosis point, and the suspected blood vessel narrow area at definite angiostenosis point place, so just need not the doctor rule of thumb positions suspected blood vessel narrow area, improved the locating speed of suspected blood vessel narrow area, simultaneously also can avoid omitting suspected blood vessel narrow area as far as possible, improve the positioning accuracy of suspected blood vessel narrow area.
Description of drawings
Fig. 1 is the flow chart of localization method of a kind of suspected blood vessel narrow area of the embodiment of the invention.
The flow chart of curve-fitting method when Fig. 2 carries out curve fitting to the vessel cross-sections diameter for the image according to the examine blood vessel.
Curve fitting sketch map when Fig. 3 carries out curve fitting for the diameter according to the position of the vessel cross-sections of obtaining and the vessel cross-sections of calculating.
The flow chart of the method for low spot on the extraction curve that Fig. 4 adopts when being the low spot that extracts on the curve.
The flow chart of the method for definite narrow zone that Fig. 5 adopts when determining the suspected blood vessel narrow area at suspected blood vessel narrow point place.
Fig. 6 is the sketch map of the suspected blood vessel narrow area that indicates.
Fig. 7 is the structure chart of navigation system of a kind of suspected blood vessel narrow area of the embodiment of the invention.
The specific embodiment
For making purpose of the present invention, technical scheme and advantage clearer, referring to the accompanying drawing embodiment that develops simultaneously, the present invention is described in more detail.
Fig. 1 is the flow chart of localization method of a kind of suspected blood vessel narrow area of the embodiment of the invention.As shown in Figure 1, the method may further comprise the steps:
Step 101 scans the examine blood vessel, and obtains the image of examine blood vessel.
In this step, medical imaging device scans the examine blood vessel, and obtains the three-dimensionalreconstruction image of examine blood vessel.
Step 102 carries out curve fitting to the vessel cross-sections diameter according to the image of examine blood vessel.
The flow chart of curve-fitting method when Fig. 2 carries out curve fitting to the vessel cross-sections diameter for the image according to the examine blood vessel.As shown in Figure 2, this approximating method may further comprise the steps:
The position of examine blood vessel initial cross sectional is 0 coordinate, and the position of other vessel cross-sections is the distance with respect to initial cross sectional.In actual applications, the initial cross sectional of examine blood vessel can be any end in the two ends of examine blood vessel.
Calculate the distance (for example distance is what pixels or a plurality of millimeters) of vessel cross-sections and two intersection points of blood vessel wall, get final product to get the diameter of vessel cross-sections.
Curve fitting sketch map when Fig. 3 carries out curve fitting for the diameter according to the position of the vessel cross-sections of obtaining and the vessel cross-sections of calculating, as shown in Figure 3, the X coordinate of curve is the position of vessel cross-sections, the Y coordinate of curve is that the X coordinate of curve is the position of vessel cross-sections, the unit of X coordinate and Y coordinate is a pixel, be that the position of cross section and the diameter of cross section all represent with number of pixels, can certainly represent with long measure (for example millimeter).
Need to prove, in actual applications, also can be every two or more pixels between the adjacent cross section.
Step 103 is carried out smothing filtering to the curve of institute's match.Like this, the noise in the curve can be removed.The particular content of smothing filtering is prior art, and it will not go into details.
The flow chart of the method for low spot on the extraction curve that Fig. 4 adopts when being the low spot that extracts on the curve.As shown in Figure 4, the method may further comprise the steps:
Step 1041, the first derivative of every bit and second dervative on the calculated curve, with first derivative be 0 and second dervative be that 0 point is as flex point.
Step 1042, N adjacent point of flex point left side on the trade-off curve, N adjacent point of flex point right side on the trade-off curve, wherein, N is positive integer, is preferably greater than 3 and less than 10 positive integer.
Whether step 1043, the Y coordinate of judging flex point less than the Y coordinate of each point in the selected 2N of top step the point, if, then with this flex point as low spot.
Step 105 is determined the suspected blood vessel narrow area at suspected blood vessel narrow point place.
The flow chart of the method for definite narrow zone that Fig. 5 adopts when determining the suspected blood vessel narrow area at suspected blood vessel narrow point place.Shown in 5, the method may further comprise the steps:
Choose the adjacent point in suspected blood vessel narrow point left side on the curve, and calculate at this tangent slope, if the absolute value of this slope is less than or equal to a preset value, the corresponding vessel position of X coordinate that then will put is as the left margin of suspected blood vessel narrow area, need to prove, preset value can be 0, also can be to be approximately 0 value for example 0.1; Otherwise, choose next point left, and calculate at this tangent slope, if be less than or equal to a preset value at the absolute value of this tangent slope, the corresponding vessel position of X coordinate that then will put is as the left margin of suspected blood vessel narrow area; Otherwise, by that analogy, continue to choose next point left, until find the absolute value of the slope of curve to be less than or equal to the point of a preset value, with the corresponding vessel position of X coordinate of this some left margin as suspected blood vessel narrow area.
Choose the adjacent point that is positioned at suspected blood vessel narrow point right side on the curve, and calculate at this tangent slope, if the absolute value of this slope is less than or equal to a preset value, the corresponding vessel position of X coordinate that then will put is as the right margin in angiostenosis zone, similarly, preset value can be 0, also can be to be approximately 0 value for example 0.1; Otherwise, choose next point to the right, and calculate at this tangent slope, if be less than or equal to a preset value at the absolute value of this tangent slope, the corresponding vessel position of X coordinate that then will put is as the right margin of suspected blood vessel narrow area; Otherwise, by that analogy, continue to choose next point to the right, until find the absolute value of the slope of curve to be less than or equal to the corresponding point of a preset value, with the corresponding vessel position of X coordinate of this some right margin as suspected blood vessel narrow area.
Need to prove, in actual applications, the diameter of doubtful narrow blood vessel is less than the diameter of normal blood vessels, and the diameter of normal blood vessels is approximately a fixed value, therefore, at suspected blood vessel narrow area, tangent slope might be positive number, also might be negative, and in the normal blood vessels zone, tangent slope is approximately 0, so, the oblique line that need to seek tangent line when determining the left margin of doubtful narrow zone is approximately corresponding point of 0 (absolute value is less than or equal to a preset value), also needs to seek tangent slope and be approximately corresponding point of 0 (absolute value is less than or equal to a preset value) when determining the right margin of doubtful narrow zone, and the zone between left margin and the right margin is suspected blood vessel narrow area.
In step 105, also can adopt other method to determine the border of suspected blood vessel narrow area.For example, with the left margin of P the point in suspected blood vessel narrow point left side as suspected blood vessel narrow area, suspected blood vessel narrow is put Q the point in right side as the right margin of suspected blood vessel narrow area, wherein, P, Q are positive integer, can equate also can not wait.P, Q are preferably 5 to 20.For example, when determining the left margin of doubtful narrow zone, choose the adjacent point in suspected blood vessel narrow point left side, if the absolute value of the slope of each point is less than or equal to a preset value in the continuous N point, wherein, M is positive integer, be preferably greater than 1 and less than 5 positive integer, the corresponding vessel position of X coordinate of any one point (being preferably from right to left first point) is as the left margin of suspected blood vessel narrow area in then this M being put, correspondingly, right margin determine also take similar method.
Step 106 according to the border of suspected blood vessel narrow area, indicates suspected blood vessel narrow area at the image of examine blood vessel.
In the middle of can suspected blood vessel narrow area being enclosed in shapes such as rectangle frame, circle, ellipses, thereby indicate suspected blood vessel narrow area, wherein, the border of suspected blood vessel narrow area can limit the left and right sides end points of rectangle frame, circle, ellipse.Fig. 6 is the sketch map of the suspected blood vessel narrow area that indicates, and wherein, rectangle frame R1, R2, R3 are the suspected blood vessel narrow area that indicates.
Based on above-mentioned suspected blood vessel narrow area localization method, Fig. 7 is the structure chart of navigation system of a kind of suspected blood vessel narrow area of the embodiment of the invention.As shown in Figure 7, this system comprises: an image acquisition unit 701, a match unit 702, a low spot extraction unit 704 and a narrow zone determining unit 705.
Wherein, image acquisition unit 701 obtains the image of examine blood vessel; Match unit 702 carries out curve fitting to the vessel cross-sections diameter according to the image of examine blood vessel; The low spot of the curve of low spot extraction unit 704 extraction institutes match is as the suspected blood vessel narrow point; Narrow zone determining unit 705 is determined the suspected blood vessel narrow area at described suspected blood vessel narrow point place.
This system further comprises: one indicates unit 706.Indicate unit 706 and indicate suspected blood vessel narrow area at the image of examine blood vessel.
This system further comprises a smothing filtering unit 703.The curve of the 703 pairs of institutes in smothing filtering unit match carries out smothing filtering.
Low spot extraction unit 704 comprises: a flex point is determined subelement 7041, a chooser unit 7042 and a judgment sub-unit 7043.Wherein, flex point is determined first derivative and the second dervative of every bit on subelement 7041 calculated curves, with first derivative be 0 and second dervative be that 0 point is as flex point; N adjacent point of flex point left side on chooser unit 7042 trade-off curves, N adjacent point of flex point right side on the trade-off curve, wherein, N is positive integer, is preferably greater than 3 and less than 10 positive integer; Whether the Y coordinate that judgment sub-unit 7043 is judged flex point less than the Y coordinate of each point in selected 2N the point, if, then with this flex point as low spot.
Narrow zone determining unit 705 comprises: left margin determines that subelement 7051 and a right margin determine subelement 7052.Wherein, left margin is determined the left margin of the suspected blood vessel narrow area at subelement 7051 definite suspected blood vessel narrow point places; Right margin is determined the right margin of the suspected blood vessel narrow area at subelement 7052 definite suspected blood vessel narrow point places.
Left margin determines that subelement 7051 comprises: left consecutive points chooser module 70511, left slope meter operator module 70512 and one get a submodule 70513 left.Wherein, left consecutive points chooser module 70511 is chosen the adjacent point in suspected blood vessel narrow point left side on the described curve; Left slope meter operator module 70512 is calculated the tangent slope of selected point, if the absolute value of this slope is less than or equal to a preset value, the corresponding vessel position of X coordinate that then will put is as the left margin of suspected blood vessel narrow area; Get a submodule 70513 left otherwise carry out; Get a submodule 70513 left and continue to choose next point left, then carry out left slope meter operator module 70512.
Right margin determines that subelement 7052 comprises: a right consecutive points chooser module 70521, a right slope meter operator module 70522) and one get a submodule 70523 to the right.Wherein, right consecutive points chooser module 70521 is chosen the adjacent point on suspected blood vessel narrow point right side on the described curve; Right slope meter operator module 70522 is calculated the tangent slope of selected point, if the absolute value of this slope is less than or equal to a preset value, the corresponding vessel position of X coordinate that then will put is as the right margin of suspected blood vessel narrow area; Get a submodule 70523 to the right otherwise carry out; Get a submodule 70523 to the right and continue to choose next point to the right, then carry out right slope meter operator module 70522.
Left margin determines that subelement 7051 comprises: left consecutive points chooser module 70511, left slope meter operator module 70512 and one get a submodule 70513 left.Wherein, left consecutive points chooser module 70511 is chosen the adjacent point in suspected blood vessel narrow point left side on the described curve; Left slope meter operator module 70512 is calculated the tangent slope of selected point, if the absolute value of the slope of each point is less than or equal to a preset value in the continuous N point, the corresponding vessel position of X coordinate of the point (being preferably from right to left first point) in then this M being put is as the left margin of suspected blood vessel narrow area; Get a submodule 70513 left otherwise carry out; Get a submodule 70513 left and choose next point left, then carry out left slope meter operator module 70512.
Right margin determines that subelement 7052 comprises: right consecutive points chooser module 70521, right slope meter operator module 70522 and one get a submodule 70523 to the right.Wherein, right consecutive points chooser module 70521 is chosen the adjacent point on suspected blood vessel narrow point right side on the described curve; Right slope meter operator module 70522 is calculated the tangent slope of selected point, if the absolute value of the slope of each point is less than or equal to a preset value in the continuous N point, the corresponding vessel position of X coordinate of the point (being preferably from left to right first point) in then this M being put is as the right margin of suspected blood vessel narrow area; Get a submodule 70523 to the right otherwise carry out; Get a submodule 70523 to the right and choose next point to the right, then carry out right slope meter operator module 70522; Wherein, M is positive integer, is preferably more than 1 and less than 5 positive integer.
The embodiment of the navigation system of a kind of suspected blood vessel narrow area provided by the present invention specifies the respective description that please refer in the embodiment of the method shown in Figure 1, repeats no more herein.
As seen, based on above-mentioned suspected blood vessel narrow area localization method and system, at first obtain the image of examine blood vessel, then the image according to the examine blood vessel carries out curve fitting to the vessel cross-sections diameter, next extracts the low spot of curve of institute's match as the suspected blood vessel narrow point, and the suspected blood vessel narrow area at definite suspected blood vessel narrow point place, so just can not omit arbitrary suspected blood vessel narrow point, can not omit arbitrary suspected blood vessel narrow area yet, improved the positioning accuracy of suspected blood vessel narrow area, simultaneously, do not need the doctor rule of thumb suspected blood vessel narrow area to be positioned yet, improved the locating speed of suspected blood vessel narrow area.
The above is preferred embodiment of the present invention only, is not for limiting protection scope of the present invention.Within the spirit and principles in the present invention all, any modification of doing, be equal to and replace and improvement etc., all should be included within protection scope of the present invention.
Claims (14)
1. the localization method of a suspected blood vessel narrow area, the method comprises:
Obtain the image of examine blood vessel;
Image according to the examine blood vessel carries out curve fitting to the vessel cross-sections diameter;
The low spot of the curve of extraction institute match is as the suspected blood vessel narrow point;
Determine left margin and the right margin of the suspected blood vessel narrow area at suspected blood vessel narrow point place; Wherein:
The left margin of determining the suspected blood vessel narrow area at suspected blood vessel narrow point place comprises:
Choose the adjacent point in suspected blood vessel narrow point left side on the described curve; And
Calculate the tangent slope of selected point, if the absolute value of this slope is less than or equal to a preset value, the corresponding vessel position of X coordinate that then will put is as the left margin of suspected blood vessel narrow area, otherwise, choose next point left, then repeat step and the subsequent step thereof of above-mentioned slope calculations; Perhaps, calculate the tangent slope of selected point, if the absolute value of the slope of a continuous N point is less than or equal to a preset value, the corresponding vessel position of X coordinate of a point in then this M being put is as the left margin of suspected blood vessel narrow area, otherwise, choose again next point left, then repeat step and the subsequent step thereof of above-mentioned slope calculations;
The right margin of determining the suspected blood vessel narrow area at suspected blood vessel narrow point place comprises:
Choose the adjacent point on suspected blood vessel narrow point right side on the described curve; And
Calculate the tangent slope of selected point, if this slope is less than or equal to preset value, the corresponding vessel position of X coordinate that then will put is as the right margin of suspected blood vessel narrow area, otherwise, choose next point to the right, then repeat step and the subsequent step thereof of above-mentioned slope calculations; Perhaps, calculate the tangent slope of selected point, if the absolute value of the slope of a continuous N point is less than or equal to a preset value, the corresponding vessel position of X coordinate of a point in then this M being put is as the right margin of suspected blood vessel narrow area, otherwise, choose again next point to the right, then repeat step and the subsequent step thereof of above-mentioned slope calculations;
Wherein, M is positive integer.
2. method according to claim 1 is characterized in that, the method further comprises: the image at the examine blood vessel indicates suspected blood vessel narrow area.
3. method according to claim 1 is characterized in that, the method further comprises: the curve to institute's match carries out smothing filtering.
4. method according to claim 1 is characterized in that, described curve fitting comprises:
For each vessel cross-sections in the image of examine blood vessel, obtain the position of this vessel cross-sections on the blood vessel bearing of trend, and calculate the diameter of this vessel cross-sections;
With the position of vessel cross-sections as the diameter of X coordinate, vessel cross-sections as the Y coordinate, fit to curve.
5. method according to claim 1 is characterized in that, the low spot of the curve of described extraction institute match comprises:
The first derivative of every bit and second dervative on the calculated curve, with first derivative be 0 and second dervative be that 0 point is as flex point;
N adjacent point of flex point left side on the trade-off curve, N adjacent point of flex point right side on the trade-off curve, wherein, N is positive integer;
Whether the Y coordinate of judging flex point less than the Y coordinate of each point in selected 2N the point, if, then with this flex point as low spot.
6. method according to claim 5 is characterized in that, described N is greater than 3 and less than 10 positive integer.
7. method according to claim 1 is characterized in that, described M is greater than 1 and less than 5 positive integer.
8. the navigation system of a suspected blood vessel narrow area, this system comprises: an image acquisition unit (701), a match unit (702), a low spot extraction unit (704) and a narrow zone determining unit (705); Wherein:
Described image acquisition unit (701) is for the image that obtains the examine blood vessel;
Described match unit (702) is used for according to the image of examine blood vessel the vessel cross-sections diameter being carried out curve fitting;
Described low spot extraction unit (704) be used for to extract the low spot of curve of institute's match as the suspected blood vessel narrow point;
Described narrow zone determining unit (705) is for left margin and the right margin of the suspected blood vessel narrow area of determining suspected blood vessel narrow point place; Described narrow zone determining unit (705) comprises that a left margin determines that subelement (7051) and a right margin determine subelement (7052), wherein:
Described left margin is determined subelement (7051), is used for the left margin of the suspected blood vessel narrow area at definite suspected blood vessel narrow point place; Described left margin determines that subelement (7051) comprising: a left consecutive points chooser module (70511), a left slope meter operator module (70512) and one get a submodule (70513) left; Wherein, described left consecutive points chooser module (70511) is used for choosing the adjacent point on the left of the suspected blood vessel narrow point on the described curve; Described left slope meter operator module (70512), be used for calculating the tangent slope of selected point, if the absolute value of this slope is less than or equal to a preset value, the corresponding vessel position of X coordinate that then will put is as the left margin of suspected blood vessel narrow area, get a submodule (70513) left otherwise carry out, perhaps be used for calculating the tangent slope of selected point, if the absolute value of the slope of a continuous N point is less than or equal to a preset value, the corresponding vessel position of X coordinate of a point in then this M being put is got a submodule (70513) as the left margin of suspected blood vessel narrow area left otherwise carry out; The described submodule (70513) of getting left is used for choosing next point left, then carries out left slope meter operator module (70512);
Described right margin is determined subelement (7052), is used for the right margin of the suspected blood vessel narrow area at definite suspected blood vessel narrow point place; Described right margin determines that subelement (7052) comprising: a right consecutive points chooser module (70521), a right slope meter operator module (70522) and one get a submodule (70523) to the right; Wherein, described right consecutive points chooser module (70521) is for the adjacent point of choosing suspected blood vessel narrow point right side on the described curve; Described right slope meter operator module (70522), be used for calculating the tangent slope of selected point, if the absolute value of this slope is less than or equal to a preset value, the corresponding vessel position of X coordinate that then will put is as the right margin of suspected blood vessel narrow area, get a submodule (70523) to the right otherwise carry out, perhaps be used for calculating the tangent slope of selected point, if the absolute value of the slope of a continuous N point is less than or equal to a preset value, the corresponding vessel position of X coordinate of a point in then this M being put is got a submodule (70523) as the right margin of suspected blood vessel narrow area to the right otherwise carry out; The described submodule (70523) of getting to the right is used for choosing next point to the right, then carries out right slope meter operator module (70522);
Wherein, M is positive integer.
9. system according to claim 8 is characterized in that, this system comprises that further indicates a unit (706), is used for indicating suspected blood vessel narrow area at the image of examine blood vessel.
10. system according to claim 8 is characterized in that, this system further comprises a smothing filtering unit (703), is used for the curve of institute's match is carried out smothing filtering.
11. system according to claim 8 is characterized in that, described match unit (702) comprising:
A position acquisition subelement (7021) is used for obtaining the position of each vessel cross-sections on the blood vessel bearing of trend;
A diameter computation subunit (7022) is for the diameter that calculates each vessel cross-sections;
A curve fitting subelement (7023) is used for the position of vessel cross-sections as the diameter of X coordinate, vessel cross-sections fitting to curve as the Y coordinate.
12. system according to claim 8 is characterized in that, described low spot extraction unit (704) comprising:
A flex point is determined subelement (7041), is used for first derivative and the second dervative of every bit on the calculated curve, with first derivative be 0 and second dervative be that 0 point is as flex point;
A chooser unit (7042) is used for N adjacent point of flex point left side on the trade-off curve, N adjacent point of flex point right side on the trade-off curve, and wherein, N is positive integer;
Whether a judgment sub-unit (7043), the Y coordinate that is used for judging flex point less than selected 2N the Y coordinate of putting each point, if, then with this flex point as low spot.
13. system according to claim 12 is characterized in that, described N is greater than 3 and less than 10 positive integer.
14. system according to claim 8 is characterized in that, described M is greater than 1 and less than 5 positive integer.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 200910169381 CN101991420B (en) | 2009-08-27 | 2009-08-27 | Method and system for positioning suspected blood vessel narrow area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 200910169381 CN101991420B (en) | 2009-08-27 | 2009-08-27 | Method and system for positioning suspected blood vessel narrow area |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101991420A CN101991420A (en) | 2011-03-30 |
| CN101991420B true CN101991420B (en) | 2013-01-16 |
Family
ID=43782530
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 200910169381 Expired - Fee Related CN101991420B (en) | 2009-08-27 | 2009-08-27 | Method and system for positioning suspected blood vessel narrow area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101991420B (en) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3277183B1 (en) * | 2015-03-31 | 2020-04-22 | Agency For Science, Technology And Research | Method and apparatus for assessing blood vessel stenosis |
| CN108451540B (en) * | 2017-02-17 | 2021-08-31 | 深圳先进技术研究院 | Fractional flow reserve measurement method and device |
| EP3607884A1 (en) * | 2018-08-07 | 2020-02-12 | Koninklijke Philips N.V. | Examination of a blood vessel based on nuclear resonant absorption |
| CN109979593B (en) * | 2018-09-24 | 2021-04-13 | 北京科亚方舟医疗科技股份有限公司 | Method for predicting healthy radius of blood vessel route, method for predicting stenosis candidate of blood vessel route, and device for predicting degree of blood vessel stenosis |
| CN109523515A (en) * | 2018-10-18 | 2019-03-26 | 深圳市孙逸仙心血管医院（深圳市心血管病研究所） | The calculation method of blood flow reserve score |
| CN109758132A (en) * | 2019-02-26 | 2019-05-17 | 数坤（北京）网络科技有限公司 | A kind of detection method and equipment of hemadostewnosis |
| CN109872321A (en) * | 2019-02-26 | 2019-06-11 | 数坤（北京）网络科技有限公司 | A kind of hemadostewnosis detection method and equipment |
| CN112472112B (en) * | 2020-11-25 | 2024-02-27 | 苏州润迈德医疗科技有限公司 | Method, system and storage medium for regulating vascular stenosis |
| CN112487342A (en) * | 2020-11-25 | 2021-03-12 | 苏州润迈德医疗科技有限公司 | Method, system and storage medium for accurately acquiring narrow lesion interval |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1551033A (en) * | 2003-03-20 | 2004-12-01 | ��ʽ���綫֥ | Processor for analyzing tubelike structure such as blood vessel |
-
2009
- 2009-08-27 CN CN 200910169381 patent/CN101991420B/en not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1551033A (en) * | 2003-03-20 | 2004-12-01 | ��ʽ���綫֥ | Processor for analyzing tubelike structure such as blood vessel |
Non-Patent Citations (2)
| Title |
|---|
| JP特开2004-283373A 2004.10.14 |
| JP特开平7-65154A 1995.03.10 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101991420A (en) | 2011-03-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101991420B (en) | Method and system for positioning suspected blood vessel narrow area | |
| CN109754387B (en) | Intelligent detection and positioning method for whole-body bone imaging radioactive concentration focus | |
| CN107833248B (en) | Medical image scanning method and medical imaging equipment | |
| CN103829965B (en) | The method and apparatus of CT scan is guided using marked body | |
| US9123096B2 (en) | Information processing apparatus and control method thereof | |
| CN108109170B (en) | Medical image scanning method and medical imaging equipment | |
| CN100463655C (en) | Image measuring device, method and image instrumentation system of glomerular filtration rate | |
| US8050475B2 (en) | Detection and localization of vascular occlusion from angiography data | |
| CN112022201A (en) | Machine guided imaging techniques | |
| US11455720B2 (en) | Apparatus for ultrasound diagnosis of liver steatosis using feature points of ultrasound image and remote medical-diagnosis method using the same | |
| EP2463822A1 (en) | Image processing apparatus, image processing method and image processing program | |
| KR20130109838A (en) | Apparatus and method for supporting lesion diagnosis | |
| CN113674257A (en) | Method, device and equipment for measuring scoliosis angle and storage medium | |
| EP2271256A1 (en) | Method and system for measuring an object of interest | |
| CN111583212B (en) | Method and device for determining brain midline shift | |
| CN110910335A (en) | Image processing method, image processing device and computer readable storage medium | |
| CN107622491B (en) | Optical fiber bundle image analysis method and device | |
| JP2007159934A (en) | Comparative diagnostic reading supporting apparatus | |
| KR102389067B1 (en) | Method for evaluating scoliosis and device evaluating scoliosis using the same | |
| CN111127550B (en) | Image slice positioning method, device and system based on bed board and storage medium | |
| CN108335732A (en) | A kind of the case recommendation method and its system of OCT images | |
| EP2950699A1 (en) | Navigation using a pre-acquired image | |
| CN110060254A (en) | Interactive vessel measurement method, device, computer readable storage medium and system | |
| CN109350062B (en) | Medical information acquisition method, medical information acquisition device and non-volatile computer storage medium | |
| CN111000580B (en) | Intervertebral disc scanning method, device, console equipment and CT system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130116 Termination date: 20210827 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |