CN104665872B - A kind of carotid intimal medial thickness measuring method and device based on ultrasonoscopy - Google Patents
A kind of carotid intimal medial thickness measuring method and device based on ultrasonoscopy Download PDFInfo
- Publication number
- CN104665872B CN104665872B CN201410841210.3A CN201410841210A CN104665872B CN 104665872 B CN104665872 B CN 104665872B CN 201410841210 A CN201410841210 A CN 201410841210A CN 104665872 B CN104665872 B CN 104665872B
- Authority
- CN
- China
- Prior art keywords
- mai
- lii
- subgraph
- ultrasonic
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0858—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
Abstract
The invention provides a kind of carotid intimal medial thickness measuring method and device based on ultrasonoscopy, for middle film in Accurate Segmentation arteria carotis, so as to effectively measure IMT.Technical solution of the present invention includes：Obtain the ultrasonic subgraph comprising interior middle film；The tube chamber intima boundary LII of the interior middle film is obtained from the ultrasonic subgraph using digital image processing techniques；Film epicardial border MAI in being obtained initially according to the LII and characteristics of image；By developing, the initial MAI obtains target MAI；Internal-media thickness IMT is calculated according to the LII and target MAI.By implementing the present invention program, it is possible to increase the precision of LII and MAI identifications, so as to Accurate Segmentation goes out interior middle film so that the accuracy of measurement of IMT is greatly improved.
Description
Technical field
The present invention relates to ultrasonic diagnostic technique field, more particularly to a kind of carotid intimal medial thickness based on ultrasonoscopy
Measuring method and device.
Background technology
Carotid intimal medial thickness (Intima-Media Thickness, IMT) is the important of prediction cardiovascular and cerebrovascular disease
Index, it is generally recognized that IMT is abnormal plump more than 1.1mm, and abnormal plump IMT then typically indicates that patient has cardiovascular and cerebrovascular disease
Disease.Fig. 1 is referred to, in arteria carotis rear wall ultrasonoscopy, includes tube chamber, inner membrance, middle film and adventitia successively.Wherein, tube chamber and inner membrance
Between border be tube chamber intima boundary (Lumen-Intima Interface, LII), during the border between middle film and adventitia is
Film epicardial border (Media-Adventitia Interface, MAI), and the distance between LII and MAI is aforesaid IMT.
Therefore, it is critical only that for IMT measurements accurately recognize LII and MAI, so as to calculate distance between the two.
In patent document (Unexamined Patent 11-318896 publication), it is proposed that in being split using the half-tone information of image
The method of middle film and measurement IMT thickness.
But, the method the inner membrance gray scale of vascular wall is relatively low or tube chamber in noise it is larger in the case of, it is impossible to correctly carry
LII borders and MAI borders are taken, thus causes IMT measurement results incorrect.In addition, the method is carried out for the IMT of Healthy People
Measurement, lack the situation for considering that blood vessel generation pathology produces patch.Therefore, above-mentioned technical proposal is divided merely with half-tone information
In cutting in arteria carotis, film there will be larger error, thus reduce the accuracy of measurement of IMT.
The content of the invention
In order to solve the above problems, the present invention provides a kind of carotid intimal medial thickness measuring method based on ultrasonoscopy
And device, for middle film in Accurate Segmentation arteria carotis, so as to effectively measure IMT.By implementing technical solution of the present invention, Neng Gouti
The precision of high LII and MAI identifications, so as to Accurate Segmentation goes out interior middle film so that the accuracy of measurement of IMT is greatly improved.
A kind of carotid intimal medial thickness measuring method based on ultrasonoscopy, including：
Obtain the ultrasonic subgraph comprising interior middle film；
The tube chamber intima boundary of the interior middle film is obtained from the ultrasonic subgraph using digital image processing techniques
LII；
Film epicardial border MAI in being obtained initially according to the LII and characteristics of image；
By developing, the initial MAI obtains target MAI；
Internal-media thickness IMT is calculated according to the LII and target MAI.
A kind of carotid intimal medial thickness measurement apparatus based on ultrasonoscopy, including：
First acquisition unit, for obtaining the ultrasonic subgraph comprising interior middle film；
Second acquisition unit, for the interior middle film is obtained from the ultrasonic subgraph using digital image processing techniques
Tube chamber intima boundary LII；
3rd acquiring unit, for film epicardial border MAI in being obtained initially according to the LII and characteristics of image；
4th acquiring unit, for obtaining target MAI by the initial MAI that develops；
Computing unit, for calculating Internal-media thickness IMT according to the LII and target MAI.
The invention has the beneficial effects as follows, with reference to features of ultrasound pattern, initial MAI is obtained according to LII first, then by developing
Initial MAI obtains target MAI, calculates IMT finally according to LII and target MAI.Technical solution of the present invention is to LII and target MAI
Identification accuracy is higher than the scheme being only identified using half-tone information such that it is able to be more accurately partitioned into interior middle film so that
The accuracy of measurement of IMT is improved.
Description of the drawings
Fig. 1 is arteria carotis rear wall ultrasonoscopy schematic diagram；
Fig. 2 is a kind of carotid intimal medial thickness measuring method flow chart based on ultrasonoscopy of the present invention；
Fig. 3 is detail flowchart of the present invention based on the carotid intimal medial thickness measuring method step 202 of ultrasonoscopy；
Fig. 4 is detail flowchart of the present invention based on the carotid intimal medial thickness measuring method step 204 of ultrasonoscopy；
Fig. 5 is detail flowchart of the present invention based on the carotid intimal medial thickness measuring method step 203 of ultrasonoscopy；
Fig. 6 is a kind of carotid intimal medial thickness measurement apparatus structure chart based on ultrasonoscopy of the present invention；
Fig. 7 is second acquisition unit structure of the present invention based on the carotid intimal medial thickness measurement apparatus of ultrasonoscopy
Figure；
Fig. 8 is fourth acquiring unit structure of the present invention based on the carotid intimal medial thickness measurement apparatus of ultrasonoscopy
Figure；
Fig. 9 is threeth acquiring unit structure of the present invention based on the carotid intimal medial thickness measurement apparatus of ultrasonoscopy
Figure.
Specific embodiment
Below in conjunction with the Figure of description in the present invention, the technical scheme in invention is clearly and completely described,
Obviously, described embodiment is only a part of embodiment of the invention, rather than the embodiment of whole.Based in the present invention
Embodiment, the every other embodiment obtained under the premise of creative work is not made by those of ordinary skill in the art, all
Belong to the scope of protection of the invention.
The embodiment of the present invention provides a kind of carotid intimal medial thickness measuring method based on ultrasonoscopy, for accurate point
Film in cutting in arteria carotis, so as to effectively measure IMT.By implementing technical solution of the present invention, it is possible to increase LII and MAI identifications
Precision, so as to Accurate Segmentation goes out interior middle film so that the accuracy of measurement of IMT is greatly improved.The embodiment of the present invention is also provided and this
A kind of related carotid intimal medial thickness measurement apparatus based on ultrasonoscopy of method, are carried out specifically to which below respectively
It is bright.
First embodiment of the invention will be carried out to a kind of carotid intimal medial thickness measuring method based on ultrasonoscopy in detail
Describe in detail bright, Fig. 2 referred to based on the carotid intimal medial thickness measuring method idiographic flow of ultrasonoscopy described in the present embodiment,
Including step：
201st, obtain the ultrasonic subgraph comprising interior middle film.
The ultrasonic subgraph that main body is interior middle film is obtained from overall ultrasonoscopy, to reduce what unrelated images information was caused
Interference, improves treatment effeciency.In this step, the mode that one of which obtains ultrasonic subgraph can be the manual sectional drawing of user,
Here it is not especially limited.
Under general default situations in arteria carotis in film segmentation be both for rear wall and carry out, but ultrasonic subgraph under individual cases
It is antetheca as be possible to presentation, therefore, also need to after obtaining the ultrasonic subgraph comprising interior middle film：Sentenced according to characteristics of image
In in disconnected, whether film is antetheca, if the determination result is YES, then carries out 180 degree upset to ultrasonic subgraph, if judged result is no,
Then do not processed.
Judge that principle is：Ultrasonic subgraph is divided into into 3 sections in a lateral direction, respectively count first paragraph, second segment and
3rd section of gray value is 0 number of pixels, if first paragraph or number of pixels that second segment gray value is 0 are more, illustrate be
Rear wall, if the 3rd section of gray value be 0 number of pixels it is more, illustrate it is antetheca.
Ultrasonic subgraph after for upset, records overturning mark, the middle film segmentation in the arteria carotis for complete the present embodiment
After need to reduce ultrasonic subgraph according to overturning mark.
202nd, the tube chamber intima boundary LII of interior middle film is obtained using digital image processing techniques from ultrasonic subgraph.
The digital image processing techniques utilized by this step are included but is not limited to：Threshold method, morphological method and searching
One kind or several combination in largest connected region method.
203rd, film epicardial border MAI in being obtained initially according to LII and characteristics of image.
MAI is located at LII rear ends, therefore, using the LII acquired in previous step, and the characteristics of image of ultrasonic subgraph
Initial MAI can be extracted in LII rear ends substantially.In the present embodiment, characteristics of image refers mainly to the gray feature of ultrasonic subgraph,
The variation characteristic that local gray-value is combined in the preset range of LII rear ends recognizes initial MAI.
204th, target MAI is obtained by the initial MAI that develops.
Specifically, initial MAI deformations are controlled with the marginal information of ultrasonoscopy, obtains the evolution curve for closing, due to MAI
Must be the curve of a non-closed, therefore also need to post-processing operation be carried out to the evolution curve of the closure, including closure is drilled
Change deconsolidation process, fitting and interpolation of curve etc., to obtain target MAI.
205th, Internal-media thickness IMT is calculated according to LII and target MAI.
Specifically, calculate the distance between LII and the multigroup correspondence profile point of target MAI set to characterize IMT, it is preferable that
The distance between same group of profile point adopts Euclidean distance.Wherein, with the maximum in distance set as IMT maximum gauges, with
Minimum of a value in distance set be IMT minimum thickness, the average thickness with the average distance of distance set as IMT.Distance set
Can be also used for calculating IMT variances etc..
Second embodiment of the invention by first embodiment step " 202, using digital image processing techniques from ultrasound
The tube chamber intima boundary LII of film in obtaining in subgraph " is described in detail, and the step particular flow sheet refers to Fig. 3, bag
Include step：
301st, obtain the binary image of ultrasonic subgraph.
As vessel segments are generally in high echo, and blood is in low echo, therefore can use upper and lower threshold method to ultrasonic subgraph
As carrying out binary conversion treatment, image is divided into into blood vessel and blood generally.
302nd, Morphological scale-space is carried out to remove flash removed to binary image, the Morphological scale-space is opening operation.
LII is general relatively flat, smooth, but due to the impact of noise, may near intimal surface in binary image
Some little burrs, it is therefore desirable to carry out morphologic opening operation to binary image, i.e., first corrode and expand afterwards, to reach elimination
Little burr, the purpose for substantially not changing its area at the very thin point while separating objects, smooth larger object border.
303rd, the largest connected region in binary image is obtained, the gray value of pixel is in the largest connected region
255。
After morphology post processing, binary image includes inner membrance, middle film, adventitia and noise.Due to noise variance
The white portion of binary image in the blood vessel, is divided into different connected regions hence with the connectedness of image by distribution.Root
Understand according to anatomical knowledge, inside and outside diaphragm area accounts for the principal status of public economy in white portion, therefore includes maximum area using finding
Connected region obtaining the overall region of interior adventitia, to generate new binary image.In new binary image, most Dalian
In logical region, the gray value of pixel is 255, and the gray value of rest of pixels is 0.
304th, LII is obtained by scanning the largest connected region.
In new binary image, black region is angiosomes, and white portion is interior adventitia overall region, they it
Between line of demarcation be LII borders.Vertical line scanning is carried out to the new bianry image, that is, divides the image into row, from the of each column
One pixel sets out, and downward search sweep records the pixel till the pixel that gray value is 255 is run into, and will
Which finally obtains complete LII as the profile point on LII.
Third embodiment of the invention by first embodiment step " 204, by developing, initial MAI obtains target
MAI " is described in detail, and the step particular flow sheet refers to Fig. 4, including step：
401st, develop the initial MAI to obtain the evolution curve of closure using evolutionary model.
Wherein, evolutionary model is GAC-GVF models or GAC models.
In the present embodiment, evolutionary model is GAC-GVF models, it is characterized in that GVF can both increase the capture model to border
Enclose, while reducing the requirement to initial evolution curve.GAC-GVF model mathematic(al) representations are：
Wherein, β is parameter of the smooth degree with border capture for controlling the profile of equilibrium, and τ is constant,For level set letter
Number,For normalized GVF,For gradient operator, g is Edge-stopping function.
Due to the top edge of adventitia be the present invention segmentation object, and the gradient direction of top edge be it is downward, therefore
When calculating g, only consider downward gradient direction, and do not consider gradient direction upwards.
By taking pixel (x, y) as an example, the definition of the g in formula (1) is：
Wherein,Represent the edge graph of ultrasonic subgraph, FyFor F partial derivatives in vertical direction.
402nd, the evolution curve is post-processed.
As the evolution curve obtained by GAC-GVF model evolutions is closure, but real MAI borders are that closure is drilled
Change the lower boundary of curve, it is therefore desirable to which the evolution curve of closure is carried out into deconsolidation process.First, closure evolution curve head and the tail are removed
4 end points, to the marginal point on remaining evolution curve, if abscissa is identical, take the larger point of ordinate value as mesh
The candidate contours point of mark MAI.Then, these candidate contours points are carried out with conic fitting, the larger point of deviation is removed, protects
The accuracy of card final curves, while using interpolation method keeping the continuity of curve.
403rd, using the evolution curve after post processing as target MAI.
Fourth embodiment of the invention by first embodiment step " 205, obtain initial according to LII and characteristics of image
Middle film epicardial border MAI " is described in detail, and the step particular flow sheet refers to Fig. 5, including step：
501st, several sampled points are equally spaced obtained on LII.
In this example, the sampling interval is 10, i.e., take a sampled point every 10 points, and record the coordinate of these points.
502nd, offset target point corresponding with sampled point is obtained according to preset rules.
Two kinds of film during film is divided to healthy interior middle film and occurs in pathology in usual, the IMT of normal person is 10 pixels or so,
Pathology be may have occurred more than 10 meanings.Generally, in normal, film is the relatively low region of gray value, and occur pathology it is interior in
Film is the higher region of gray value.
According to above-mentioned principle, preset rules are set as：Obtain sampled point average in 10 pixel regions of outer radial
Gray value.If the average gray value is less than first threshold, start to search offset target point, the offset target successively from sampled point
Point meets the local gray-value of continuous 4 pixels thereafter and is all higher than Second Threshold.If average gray value is more than first threshold,
Start to search offset target point successively from the 11st pixel, the offset target point meets its local gray-value more than continuous thereafter
The local gray-value of 3 pixels.
503rd, the curve for linking all offset target points successively is used as initial MAI.
Fifth embodiment of the invention will be carried out to a kind of carotid intimal medial thickness measurement apparatus based on ultrasonoscopy in detail
Describe the bright carotid intimal medial thickness measurement side based on ultrasonoscopy in device and previous embodiment described in the present embodiment in detail
Method is related, therefore the part or all of description in previous embodiment is applied to the present embodiment and subsequent embodiment, the present embodiment institute
The device concrete structure stated refers to Fig. 6, including：
First acquisition unit 601, second acquisition unit 602, the 3rd acquiring unit 603, the 4th acquiring unit 604 and calculating
Unit 605.Wherein, first acquisition unit 601, second acquisition unit 602, the 3rd acquiring unit 603, the 4th acquiring unit 604,
Computing unit 605 is communicated to connect successively.
First acquisition unit 601, for obtaining the ultrasonic subgraph comprising interior middle film.
In embodiments of the present invention, according to characteristics of image judge in film whether be antetheca, it is if the determination result is YES, then right
Ultrasonic subgraph carries out 180 degree upset, if judged result is no, is not processed.
Therefore, the device described in the present embodiment can also include judging unit 606 and roll-over unit 607.Wherein, first obtain
Take unit 601, judging unit 606, roll-over unit 607 to communicate to connect successively.
Judging unit 606, for judging whether the interior middle film is antetheca according to characteristics of image.
Roll-over unit 607, if the judged result for the judging unit 606 is yes, enters to the ultrasonic subgraph
Row 180 degree overturns, and is not otherwise processed.
Second acquisition unit 602, for being obtained from the ultrasonic subgraph in described using digital image processing techniques
The tube chamber intima boundary LII of middle film.
The digital image processing techniques utilized by second acquisition unit 602 are included but is not limited to：Threshold method, morphology side
One kind or several combination in method and the largest connected region method of searching.
3rd acquiring unit 603, for film epicardial border MAI in being obtained initially according to the LII and characteristics of image.
MAI is located at LII rear ends, therefore, using the LII acquired in second acquisition unit 602, and the figure of ultrasonic subgraph
As feature, the 3rd acquiring unit 603 substantially can extract initial MAI in LII rear ends.In the present embodiment, characteristics of image is main
Refer to the gray feature of ultrasonic subgraph, the variation characteristic identification that local gray-value is combined in the preset range of LII rear ends is initial
MAI。
4th acquiring unit 604, for obtaining target MAI by the initial MAI that develops.
Specifically, initial MAI deformations are controlled with the marginal information of ultrasonoscopy, obtains the evolution curve for closing, due to MAI
Must be the curve of a non-closed, therefore also need to post-processing operation be carried out to the evolution curve of the closure, including closure is drilled
Change deconsolidation process, fitting and interpolation of curve etc., to obtain target MAI.
Computing unit 605, for calculating Internal-media thickness IMT according to the LII and target MAI.
Specifically, calculate the distance between LII and the multigroup correspondence profile point of target MAI set to characterize IMT, it is preferable that
The distance between same group of profile point adopts Euclidean distance.Wherein, with the maximum in distance set as IMT maximum gauges, with
Minimum of a value in distance set be IMT minimum thickness, the average thickness with the average distance of distance set as IMT.Distance set
Can be also used for calculating IMT variances etc..
Sixth embodiment of the invention will be described in detail to the second acquisition unit 602 in the 5th embodiment.Second obtains
The concrete structure for taking unit 602 refers to Fig. 7, including：
Binaryzation subelement 6021, form sub-units 6022, connected domain subelement 6023 and LII subelements 6024.Its
In, binaryzation subelement 6021, form sub-units 6022, connected domain subelement 6023, LII subelements 6024 communication link successively
Connect.
Binaryzation subelement 6021, for obtaining the binary image of ultrasonic subgraph.
As vessel segments are generally in high echo, and blood is in low echo, therefore can use upper and lower threshold method to ultrasonic subgraph
As carrying out binary conversion treatment, image is divided into into blood vessel and blood generally.
Form sub-units 6022, for carrying out Morphological scale-space to go flash removed, the shape to the binary image
State is processed as opening operation.
LII is general relatively flat, smooth, but due to the impact of noise, may near intimal surface in binary image
Some little burrs, it is therefore desirable to carry out morphologic opening operation to binary image, i.e., first corrode and expand afterwards, to reach elimination
Little burr, the purpose for substantially not changing its area at the very thin point while separating objects, smooth larger object border.
Connected domain subelement 6023, it is for obtaining the largest connected region in the binary image, described largest connected
In region, the gray value of pixel is 255.
After morphology post processing, binary image includes inner membrance, middle film, adventitia and noise.Due to noise variance
The white portion of binary image in the blood vessel, is divided into different connected regions hence with the connectedness of image by distribution.Root
Understand according to anatomical knowledge, inside and outside diaphragm area accounts for the principal status of public economy in white portion, therefore includes maximum area using finding
Connected region obtaining the overall region of interior adventitia, to generate new binary image.In new binary image, most Dalian
In logical region, the gray value of pixel is 255, and the gray value of rest of pixels is 0.
LII subelements 6024, for obtaining LII by scanning the largest connected region.
In new binary image, black region is angiosomes, and white portion is interior adventitia overall region, they it
Between line of demarcation be LII borders.Vertical line scanning is carried out to the new bianry image, that is, divides the image into row, from the of each column
One pixel sets out, and downward search sweep records the pixel till the pixel that gray value is 255 is run into, and will
Which finally obtains complete LII as the profile point on LII.
Seventh embodiment of the invention will be described in detail to the 4th acquiring unit 604 in the 5th embodiment.4th obtains
The concrete structure for taking unit 604 refers to Fig. 8, including：
Evolution subelement 6041 and post processing subelement 6042.Wherein, evolution subelement 6041, post processing subelement 6042
Communication connection.
Evolution subelement 6041, for developing the initial MAI to obtain the evolution curve of closure using evolutionary model.
Wherein, evolutionary model is GAC-GVF models or GAC models.
In the present embodiment, evolutionary model is GAC-GVF models, it is characterized in that GVF can both increase the capture model to border
Enclose, while reducing the requirement to initial evolution curve.GAC-GVF model mathematic(al) representations are：
Wherein, β is parameter of the smooth degree with border capture for controlling the profile of equilibrium, and τ is constant,For level set letter
Number,For normalized GVF,For gradient operator, g is Edge-stopping function.
Due to the top edge of adventitia be the present invention segmentation object, and the gradient direction of top edge be it is downward, therefore
When calculating g, only consider downward gradient direction, and do not consider gradient direction upwards.
By taking pixel (x, y) as an example, the definition of the g in formula (1) is：
Wherein,Represent the edge graph of ultrasonic subgraph, FyFor F partial derivatives in vertical direction.
Post processing subelement 6042, for post-processing to the evolution curve, the 4th acquiring unit is by rear place
Evolution curve after reason is used as target MAI.
As the evolution curve obtained by GAC-GVF model evolutions is closure, but real MAI borders are that closure is drilled
Change the lower boundary of curve, it is therefore desirable to which the evolution curve of closure is carried out into deconsolidation process.First, closure evolution curve head and the tail are removed
4 end points, to the marginal point on remaining evolution curve, if abscissa is identical, take the larger point of ordinate value as mesh
The candidate contours point of mark MAI.Then, these candidate contours points are carried out with conic fitting, the larger point of deviation is removed, protects
The accuracy of card final curves, while using interpolation method keeping the continuity of curve.
Finally, the 4th acquiring unit 604 will post-process the evolution curve after subelement 6042 is post-processed as target
MAI。
Eighth embodiment of the invention will be described in detail to the 3rd acquiring unit 603 in the 5th embodiment.3rd obtains
The concrete structure for taking unit 603 refers to Fig. 9, including：
Sampled point subelement 6031, offset point subelement 6032 and point link subelement 6033.Wherein, sampled point subelement
6031st, offset point subelement 6032, point link subelement 6033 and communicate to connect successively.
Sampled point subelement 6031, for several sampled points are equally spaced obtained on the LII.
In this example, the sampling interval is 10, i.e., take a sampled point every 10 points, and record the coordinate of these points.
Offset point subelement 6032, for obtaining offset target point corresponding with sampled point according to preset rules.
Two kinds of film during film is divided to healthy interior middle film and occurs in pathology in usual, the IMT of normal person is 10 pixels or so,
Pathology be may have occurred more than 10 meanings.Generally, in normal, film is the relatively low region of gray value, and occur pathology it is interior in
Film is the higher region of gray value.
According to above-mentioned principle, the preset rules of 6032 foundation of offset point subelement are：Sampled point is obtained along outer radial 10
Average gray value in pixel region.If the average gray value is less than first threshold, start to search successively partially from sampled point
Impact point is moved, the offset target point meets the local gray-value of continuous 4 pixels thereafter and is all higher than Second Threshold.If average ash
Angle value is more than first threshold, then start to search offset target point successively from the 11st pixel, and the offset target point meets its office
Local gray-value of portion's gray value more than continuous 3 pixels thereafter.
Point links subelement 6033, for the curve that links all offset target points successively as initial MAI.
The invention has the beneficial effects as follows, with reference to features of ultrasound pattern, initial MAI is obtained according to LII first, then by developing
Initial MAI obtains target MAI, calculates IMT finally according to LII and target MAI.Technical solution of the present invention is to LII and target MAI
Identification accuracy is higher than the scheme being only identified using half-tone information such that it is able to be more accurately partitioned into interior middle film so that
The accuracy of measurement of IMT is improved.
A kind of carotid intimal medial thickness measuring method based on ultrasonoscopy that above embodiment of the present invention is provided
It is described in detail with device, but the explanation of above example is only intended to help the structure and its core that understand the present invention to think
Think, should not be construed as limiting the invention.Those skilled in the art the invention discloses technical scope in, can be light
The change or replacement being readily conceivable that, should all be included within the scope of the present invention.
Claims (12)
1. a kind of carotid intimal medial thickness measuring method based on ultrasonoscopy, it is characterised in that include：
Obtain the ultrasonic subgraph comprising interior middle film；
The tube chamber intima boundary LII of the interior middle film is obtained from the ultrasonic subgraph using digital image processing techniques；
Film epicardial border MAI in being obtained initially according to the LII and characteristics of image；Described image is characterized as the ash of ultrasonic subgraph
Degree feature；
By developing, the initial MAI obtains target MAI；
Internal-media thickness IMT is calculated according to the LII and target MAI.
2. method according to claim 1, it is characterised in that the utilization digital image processing techniques are from ultrasound
The tube chamber intima boundary LII that the interior middle film is obtained in image includes：
Obtain the binary image of the ultrasonic subgraph；
Morphological scale-space is carried out to the binary image to remove flash removed, the Morphological scale-space is opening operation；
The largest connected region in the binary image is obtained, the gray value of pixel is 255 in the largest connected region；
LII is obtained by scanning the largest connected region.
3. method according to claim 1 and 2, it is characterised in that described that target is obtained by the initial MAI that develops
MAI includes：
The initial MAI develop using evolutionary model to obtain the evolution curve of closure, the evolutionary model is GAC-GVF models
Or GAC models；
The evolution curve is post-processed；
Using the evolution curve after post processing as target MAI.
4. method according to claim 1 and 2, it is characterised in that described to obtain initial according to the LII and characteristics of image
Middle film epicardial border MAI includes：
Several sampled points are obtained on the LII equally spaced；
Offset target point corresponding with sampled point is obtained according to preset rules；
The curve that all offset target points are linked successively is used as initial MAI.
5. method according to claim 4, it is characterised in that the preset rules are：
Obtain average gray value of the sampled point in 10 pixel regions of outer radial；
If the average gray value is less than first threshold, start to search offset target point successively from sampled point, the skew mesh
Punctuate meets the local gray-value of continuous 4 pixels thereafter and is all higher than Second Threshold；
If the average gray value is more than first threshold, start to search offset target point successively from the 11st pixel, it is described
Offset target point meets local gray-value of its local gray-value more than continuous 3 pixels thereafter.
6. method according to claim 1, it is characterised in that after ultrasonic subgraph of the acquisition comprising interior middle film and also
Including：
Judge whether the interior middle film is antetheca according to characteristics of image；
If the determination result is YES, then 180 degree upset is carried out to the ultrasonic subgraph, if judged result is no, is not located
Reason.
7. a kind of carotid intimal medial thickness measurement apparatus based on ultrasonoscopy, it is characterised in that include：
First acquisition unit, for obtaining the ultrasonic subgraph comprising interior middle film；
Second acquisition unit, for obtaining the pipe of the interior middle film using digital image processing techniques from the ultrasonic subgraph
Endometrial cavity border LII；
3rd acquiring unit, for film epicardial border MAI in being obtained initially according to the LII and characteristics of image；Described image is special
Levy the gray feature for ultrasonic subgraph；
4th acquiring unit, for obtaining target MAI by the initial MAI that develops；
Computing unit, for calculating Internal-media thickness IMT according to the LII and target MAI.
8. device according to claim 7, it is characterised in that the second acquisition unit includes：
Binaryzation subelement, for obtaining the binary image of the ultrasonic subgraph；
Form sub-units, for carrying out Morphological scale-space to remove flash removed, the Morphological scale-space to the binary image
For opening operation；
Connected domain subelement, for obtaining the largest connected region in the binary image, picture in the largest connected region
The gray value of vegetarian refreshments is 255；
LII subelements, for obtaining LII by scanning the largest connected region.
9. the device according to claim 7 or 8, it is characterised in that the 4th acquiring unit includes：
Evolution subelement, for developing the initial MAI to obtain the evolution curve of closure, the evolution mould using evolutionary model
Type is GAC-GVF models or GAC models；
Post processing subelement, for post-processing to the evolution curve, the 4th acquiring unit by post processing after drill
Change curve as target MAI.
10. the device according to claim 7 or 8, it is characterised in that the 3rd acquiring unit includes：
Sampled point subelement, for several sampled points are equally spaced obtained on the LII；
Offset point subelement, for obtaining offset target point corresponding with sampled point according to preset rules；
Point links subelement, for the curve that links all offset target points successively as initial MAI.
11. devices according to claim 10, it is characterised in that the preset rules of the offset point subelement foundation are：
Obtain average gray value of the sampled point in 10 pixel regions of outer radial；
If the average gray value is less than first threshold, start to search offset target point successively from sampled point, the skew mesh
Punctuate meets the local gray-value of continuous 4 pixels thereafter and is all higher than Second Threshold；
If the average gray value is more than first threshold, start to search offset target point successively from the 11st pixel, it is described
Offset target point meets local gray-value of its local gray-value more than continuous 3 pixels thereafter.
12. devices according to claim 7, it is characterised in that described device also includes：
Judging unit, for judging whether the interior middle film is antetheca according to characteristics of image；
Roll-over unit, if the judged result for the judging unit is yes, carries out 180 degree to the ultrasonic subgraph and turns over
Turn, otherwise do not processed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410841210.3A CN104665872B (en) | 2014-12-29 | 2014-12-29 | A kind of carotid intimal medial thickness measuring method and device based on ultrasonoscopy |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410841210.3A CN104665872B (en) | 2014-12-29 | 2014-12-29 | A kind of carotid intimal medial thickness measuring method and device based on ultrasonoscopy |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104665872A CN104665872A (en) | 2015-06-03 |
| CN104665872B true CN104665872B (en) | 2017-04-05 |
Family
ID=53301841
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410841210.3A Active CN104665872B (en) | 2014-12-29 | 2014-12-29 | A kind of carotid intimal medial thickness measuring method and device based on ultrasonoscopy |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104665872B (en) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105380682B (en) * | 2015-12-24 | 2018-06-08 | 飞依诺科技(苏州)有限公司 | Bilateral intima-media membrane measuring method and system based on ultrasonic images |
| CN105976335B (en) * | 2016-05-09 | 2018-11-09 | 苏州佳世达电通有限公司 | Image treatment method |
| CN106570856A (en) * | 2016-08-31 | 2017-04-19 | 天津大学 | Common carotid artery intima-media thickness measuring device and method combining level set segmentation and dynamic programming |
| CN109688939A (en) * | 2016-09-12 | 2019-04-26 | 富士胶片株式会社 | The control method of ultrasonic diagnostic system and ultrasonic diagnostic system |
| CN106570871B (en) * | 2016-11-02 | 2019-10-18 | 天津大学 | Fuzzy C-mean algorithm carotid ultrasound image Internal-media thickness measurement method and system |
| US11341633B2 (en) | 2017-05-31 | 2022-05-24 | Edan Instruments, Inc. | Systems and methods for adaptive enhancement of vascular imaging |
| CN109493383B (en) * | 2018-11-23 | 2022-02-11 | 深圳市威尔德医疗电子有限公司 | Method for measuring intima-media thickness in ultrasonic image, server and storage medium |
| CN110517263B (en) * | 2019-09-02 | 2022-05-20 | 青岛海信医疗设备股份有限公司 | Method and device for determining inner middle film thickness and storage medium |
| CN112545566A (en) * | 2020-12-07 | 2021-03-26 | 武汉中旗生物医疗电子有限公司 | Fetal neck diaphragmatic layer measuring method and device |
| CN112932542B (en) * | 2021-01-27 | 2022-11-08 | 深圳中科乐普医疗技术有限公司 | Method and system for measuring thickness of intravascular medium membrane and ultrasonic imaging equipment |
| CN113171125A (en) * | 2021-03-31 | 2021-07-27 | 上海深至信息科技有限公司 | Carotid intima-media thickness measurement system and method |
| CN113160265B (en) * | 2021-05-13 | 2022-07-19 | 四川大学华西医院 | Construction method of prediction image for brain corpus callosum segmentation for corpus callosum state evaluation |
| CN113643319B (en) * | 2021-10-18 | 2021-12-10 | 成都汇声科技有限公司 | Fat thickness starting line calibrated on image and fat thickness measuring method |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6835177B2 (en) * | 2002-11-06 | 2004-12-28 | Sonosite, Inc. | Ultrasonic blood vessel measurement apparatus and method |
| WO2008023618A1 (en) * | 2006-08-21 | 2008-02-28 | Tohoku University | Ultrasonograph |
| CN102163326B (en) * | 2010-12-22 | 2013-03-06 | 武汉沃生科学技术研究中心有限公司 | Method for automatic computerized segmentation and analysis on thickness uniformity of intima media of carotid artery blood wall in sonographic image |
| CN103996194B (en) * | 2014-05-23 | 2016-08-31 | 华中科技大学 | A kind of based on film automatic division method middle in ultrasonic carotid images |
-
2014
- 2014-12-29 CN CN201410841210.3A patent/CN104665872B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN104665872A (en) | 2015-06-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104665872B (en) | A kind of carotid intimal medial thickness measuring method and device based on ultrasonoscopy | |
| CN103514605B (en) | Based on the choroid layer automatic division method of HD-OCT retina image | |
| CN102890824B (en) | Motion object outline tracking and device | |
| CN105261015B (en) | Eye fundus image blood vessel automatic division method based on Gabor filter | |
| CN105488796A (en) | Lung segmentation method | |
| CN107766874B (en) | Measuring method and measuring system for ultrasonic volume biological parameters | |
| US20150310638A1 (en) | System and method of medical imaging | |
| CN103810709A (en) | Vessel-based registration method for eye fundus image and SD-OCT projection image | |
| CN108257126B (en) | Blood vessel detection and registration method, equipment and application of three-dimensional retina OCT image | |
| JP2019192215A (en) | 3d quantitative analysis of retinal layers with deep learning | |
| CN109214397A (en) | The dividing method of Lung neoplasm in a kind of lung CT image | |
| CN105913442A (en) | Method for automatically matching pulmonary nodules | |
| CN110782428B (en) | Method and system for constructing clinical brain CT image ROI template | |
| US20160051214A1 (en) | Three-dimensional cardiac image segmentation method and apparatus | |
| CN103186901A (en) | Full-automatic image segmentation method | |
| CN106570856A (en) | Common carotid artery intima-media thickness measuring device and method combining level set segmentation and dynamic programming | |
| WO2018113282A1 (en) | Method and system for processing three-dimensional pelvic floor ultrasound image | |
| CN104361554B (en) | A kind of externa automatic testing method based on ivus image | |
| CN109087310B (en) | Meibomian gland texture region segmentation method and system, storage medium and intelligent terminal | |
| CN101425140A (en) | Encephalic angioma image recognizing and detecting method based on framework characteristic | |
| CN100383811C (en) | Method for extracting tongue body from tongue images | |
| CN104490475B (en) | Probe emission lines angle automatic adjusting method, apparatus and system in one kind puncture | |
| JPWO2019073962A1 (en) | Image processing apparatus and program | |
| CN110517264B (en) | Nidus extraction method and device based on blood vessel segmentation | |
| CN109948622A (en) | A kind of neck body aneurysm detection method, apparatus and computer readable storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP02 | Change in the address of a patent holder |
Address after: 518057 Shenzhen Nanshan District, Guangdong Province Patentee after: Biologic medical Science and Technology Co., Ltd. is opened in Shenzhen Address before: 518051 Fourth Floor, Yizhe Building, Yuquan Road, Nanshan District, Shenzhen City, Guangdong Province Patentee before: Biologic medical Science and Technology Co., Ltd. is opened in Shenzhen |
|
| CP02 | Change in the address of a patent holder |