CN103914845B - The method obtaining initial profile in Ultrasound Image Segmentation based on active contour model - Google Patents
The method obtaining initial profile in Ultrasound Image Segmentation based on active contour model Download PDFInfo
- Publication number
- CN103914845B CN103914845B CN201410141568.5A CN201410141568A CN103914845B CN 103914845 B CN103914845 B CN 103914845B CN 201410141568 A CN201410141568 A CN 201410141568A CN 103914845 B CN103914845 B CN 103914845B
- Authority
- CN
- China
- Prior art keywords
- ultrasonoscopy
- subimage
- tumor
- vector
- standard
- 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
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
Abstract
The invention discloses a kind of method obtaining initial profile in Ultrasound Image Segmentation based on active contour model, comprise the steps of the textural characteristics of training ultrasonoscopy tumor region, and constitute a standard vector together with tumor priori size in standard ellipse empirical value, ultrasonoscopy to be measured;Ultrasonoscopy to be measured is carried out pretreatment;Image is carried out dynamic threshold segmentation;Extract all closed contours produced in dynamic threshold segmentation result, form corresponding subimage;Calculate the textural characteristics 24 of each subimage, ellipse fitting result parameter 1 and subimage closed contour interior pixels number, and constitute a vector by these 26 data;Calculate the distance of gained vector and the standard vector of tumor region training gained, determine that the closed contour in the corresponding subimage of vector that distance is minimum is the profile of tumor in split ultrasonoscopy.The present invention overcomes the many noises of ultrasonoscopy, the problem of obscurity boundary, initial profile accuracy rate is high.
Description
Technical field
A kind of method that the present invention relates to obtain initial profile in Ultrasound Image Segmentation based on active profile, belongs to ultrasonic
Image processing field.
Background technology
Along with the appearance of high intensity focused ultrasound, non-invasively treating oncologic application more and more extensive, this treatment system
The biggest advantage is had if do not operated in clinic, preventing from scar, noinvasive or Micro trauma, can detect treatment, in real time not by swollen
Tumor size limits, and total cost is low.These advantages all determine that high-strength focus supersonic therapeutic system has the biggest development and application
Prospect.And in high-strength focus supersonic therapeutic system most critical exactly patient's tumor is carried out real-time navigation.Navigation procedure needs
Tumor in real-time ultrasonic image is positioned.Traditional manual positioning mode is changed into and is automatically positioned, this will be controlled
The therapeutic efficiency for the treatment of system and treatment accuracy bring the biggest help.
Ultrasound Image Segmentation be lesion detection, diagnose and treat in be widely used, its segmentation accuracy directly close
It is tied to the location of tumor.Through long-term further investigation, it has been suggested that ultrasonic image division method have based on region, based on
Border and based on statistical information, different types of dividing method has different segmentation features.Active contour model
(Active Contour Model), be otherwise known as Snake, is that a kind of objective contour proposed by Andrew professor Blake is retouched
State method, be mainly used in Target Segmentation based on shape.This model will be appreciated that it for far-ranging one is
Row visual problem gives unified solution, and in nearest more than ten years, it is by the success of increasing researcher
Be applied to many fields of computer vision, such as edge extracting, image segmentation and classification, motion tracking, three-dimensional reconstruction, three-dimensional
Vision matching etc..
But this model needs initial profile, restraining again if manually obtaining initial profile, i.e. wasting time and energy, the most not
Enough accurate.Therefore the method for design acquisition initial profile automatically is very important.
Summary of the invention
The technical problem to be solved is to provide in a kind of Ultrasound Image Segmentation based on active contour model
The method obtaining initial profile.Operation is simple and reliable for the method, speed fast, the method the ultrasonoscopy tumor obtained initially is taken turns
Wide accuracy is the highest.
The method automatically extracting tumor initial profile in a kind of Ultrasound Image Segmentation of the present invention, comprises the steps of
(1) choosing more than 60 width existing ultrasonoscopy composition training set, described existing ultrasonoscopy contains and treats point
Cut the tumor that ultrasonoscopy is identical, from training set, every width ultrasonoscopy extracts the 24 of tumor region with gray level co-occurrence matrix
Individual textural characteristics, in utilizing mathematical linear to return minimum distance method obtain the standard value of 24 textural characteristics, and and standard
In oval empirical value, ultrasonoscopy to be measured, tumor priori size constitutes a standard vector A together0;
(2) ultrasonoscopy to be split is carried out the pretreatment of anisotropy parameter and histogram equalization;
(3) (2) gained image is carried out dynamic threshold segmentation, and extract all closing of producing in dynamic threshold segmentation result
Closing profile, each closed contour forms a corresponding subimage;
(4) utilize gray level co-occurrence matrix to extract 24 textural characteristics of each subimage, calculate ellipse fitting result ginseng
Number and subimage closed contour interior pixels number, and constitute a vectorial A by these 26 data1;
(5) each subimage gained vector A is calculated1Standard vector A with tumor region training gained0Distance D, make
Obtain the minimum vectorial A of D1Closed contour in corresponding subimage is the initial profile of tumor in ultrasonoscopy to be split.
In step (1), standard ellipse empirical value is defined as 0.7.
Above-mentioned 24 textural characteristics include image 0 °, 45 °, 90 °, contrast on 135 ° of four directions, energy, relevant
Property, unfavourable balance square, entropy and each six features of non-similarity.
The acquisition of initial profile to be made is more accurate, is determining standard vector A0During, need to use reliably
Mathematical measure.Such as during training, supersonic tumor image has 60 width, every piece image to be desirable that to obtain a texture feature vector
Ai, then how to process these 60 vectors to obtain a standard vector A0Need to think better of.A is utilized after in conjunction with0Come
Method used when determining tumor region: calculate characteristic vector A extracted from each subimage1With A0Distance, then
The most also this reference of available range obtains standard vector A0.Concrete execution method is the linear regression in mathematics, will certain
60 eigenvalues that individual textural characteristics is had correspond in coordinate system, then obtain data in coordinate system so that these are 60 years old
Individual value is the shortest to the distance sum of this data place straight line, then these data are exactly the standard value of this textural characteristics, successively may be used
In the standard value the most corresponding to obtain these 24 textural characteristics, these standard values and standard ellipse data, ultrasonoscopy to be measured
Tumor priori size constitutes a standard vector A together0。
When in ultrasonoscopy the most to be measured, tumor priori size is by ultrasound detection, doctor or expert's Direct Mark go out tumor
Region major and minor axis, unit is generally millimeter, then according to physical length in ultrasonoscopy and pixel ratio, such as every millimeter 4
Pixel, can be calculated the priori size of tumor region.
Above-mentioned each subimage can only comprise a closed contour.
Method used by above-mentioned computed range D is χ2Range formula is tried to achieve:
A0(k) and A1K () represents vector A respectively0And A1Middle kth parameter, k takes the integer between 0~25.
The invention have benefit that:
(1) in place of changing the clumsiness that traditional manual obtains initial profile.
(2) initial profile acquisition methods is simple, quickly, it is adaptable to movable contour model
(3) initial profile of gained is clearly accurate, has the most excellent as active contour afterwards for convergence
Point, substantially increases efficiency and the accuracy of Ultrasound Image Segmentation.
Accompanying drawing explanation
Fig. 1 is the flow chart of automatic profile acquisition methods of the present invention.
Fig. 2 is a width liver tumor image initial profile results figure of present example 1.Fig. 2 (a) is the ultrasonic figure of liver neoplasm
As original graph;Fig. 2 (b) is original image pre-processed results after anisotropy parameter and histogram equalization;Fig. 2 (c) is
Pretreatment image dynamic threshold segmentation result;Fig. 2 (d) is the initial profile result obtained.
Fig. 3 is two width hysteromyoma ultrasonoscopy initial profile result figures of present example 2, and result is as actively
The initial profile of skeleton pattern, obtains the result figure of final tumor boundaries after the convergence under the GGVF field of force drives.
Detailed description of the invention
Carry out the more detailed description present invention below by way of some detailed description of the invention, but the present invention is not limited to these in fact
Execute example.
Embodiment 1:
Through liver neoplasm ultrasonoscopy existing to 60 width be trained the standard feature obtaining liver tumor ultrasonoscopy to
Amount A0.The most each width imagery exploitation gray level co-occurrence matrixes is calculated its 24 textural characteristics values, then 60 can be obtained
These all characteristic vectors are carried out mathematical linear homing method process, obtain a mark by the group characteristic vector containing 24 elements
Quasi-characteristic vector, i.e. by the vectorial A that searching one to the average distance of this stack features vector is the shortest, correspondence in this vector A
Each textural characteristics value is exactly the standard value that we need to obtain.Then this feature vector A is added standard ellipse fitting data
0.7 and ultrasonoscopy to be measured in the calculated liver of tumor region major and minor axis that marked by doctor during ultrasound detection swell
Two elements of tumor priori size, then just obtained a standard feature vector A containing 26 elements0。
Fig. 2 (a) is liver tumor ultrasonoscopy I to be split, and Fig. 2 (b) is for carry out anisotropy parameter and Nogata by original image
The pre-processed results of gained after figure equilibrium treatment;Fig. 2 (c) is for carry out the knot of dynamic threshold segmentation to gained image in Fig. 2 (b)
Really, and extract all closed contours produced in dynamic threshold segmentation result, form corresponding subimage;Utilize gray level symbiosis square
Battle array calculates the textural characteristics 24 of each subimage, calculates each subimage ellipse fitting result parameter and subimage Guan Bi wheel
Wide interior pixels number, and constitute characteristic vector A by these 26 data1;Utilize χ2Range formula calculates each subimage
Gained vector A1Standard vector A with tumor region training gained0Distance D, determine the vectorial A minimum so that D1Corresponding subgraph
Closed contour in Xiang is the initial profile of tumor in ultrasonoscopy to be split, is the initial profile of acquisition shown in Fig. 2 (d)
Result.
Embodiment 2:
It is trained obtaining the standard feature of hysteromyoma ultrasonoscopy through hysteromyoma ultrasonoscopy existing to 60 width
Vector A0.The most each width imagery exploitation gray level co-occurrence matrixes is calculated its 24 textural characteristics values, then can obtain
These all characteristic vectors are carried out mathematical linear homing method process, obtain one by 60 groups of characteristic vectors containing 24 elements
Standard feature vector is i.e. by the vectorial A that searching one to the average distance of this stack features vector is the shortest, corresponding in this vector A
Each textural characteristics value be exactly we need obtain standard value.Then this feature vector A is added standard ellipse matching number
Calculated according to the tumor region major and minor axis marked by doctor during ultrasound detection in 0.7 and ultrasonoscopy to be measured
Two elements of hysteromyoma priori size, then just obtained a standard feature vector A containing 26 elements0。
Fig. 3 (a1) is hysteromyoma ultrasonoscopy I to be split1, Fig. 3 (b1) for original image carried out anisotropy parameter and
Histogram equalization processes the pre-processed results of gained later;Fig. 3 (c1) divides for gained image in Fig. 3 (b1) is carried out dynamic threshold
The result cut, and extract all closed contours produced in dynamic threshold segmentation result, form corresponding subimage;Utilize gray level
Co-occurrence matrix calculates the textural characteristics 24 of each subimage, calculates each subimage ellipse fitting result parameter and subimage
Closed contour interior pixels number, and constitute characteristic vector A by these 26 data1;Utilize χ2Range formula calculates each
Subimage gained vector A1Standard vector A with tumor region training gained0Distance D, determine the vectorial A minimum so that D1Right
Answering the closed contour in subimage is the initial profile of tumor in ultrasonoscopy to be split, is acquisition shown in Fig. 3 (d1)
Initial profile result.Fig. 3 (e1) is then as the initial profile of active contour model using the initial profile of gained in Fig. 3 (d1),
It is then passed through GGVF external force field and carries out restraining the final hysteromyoma border result figure of gained.
Fig. 3 (a2) is hysteromyoma ultrasonoscopy I to be split2, Fig. 3 (b2) for original image carried out anisotropy parameter and
Histogram equalization processes the pre-processed results of gained later;Fig. 3 (c2) divides for gained image in Fig. 3 (b2) is carried out dynamic threshold
The result cut, and extract all closed contours produced in dynamic threshold segmentation result, form corresponding subimage;Utilize gray level
Co-occurrence matrix calculates the textural characteristics 24 of each subimage, calculates each subimage ellipse fitting result parameter and subimage
Closed contour interior pixels number, and constitute characteristic vector A by these 26 data1;Utilize χ2Range formula calculates each
Subimage gained vector A1Standard vector A with tumor region training gained0Distance D, determine the vectorial A minimum so that D1Right
Answering the closed contour in subimage is the initial profile of tumor in ultrasonoscopy to be split, is acquisition shown in Fig. 3 (d2)
Initial profile result.Fig. 3 (e2) is then as the initial profile of active contour model using the initial profile of gained in Fig. 3 (d2),
It is then passed through GGVF external force field and carries out restraining the final hysteromyoma border result figure of gained.
In embodiment, the extraction of initial profile does not need implementer manually to provide initial profile, and this clearly alleviates
The workload of implementer, improves the efficiency of Ultrasound Image Segmentation.
Claims (2)
1. the method obtaining initial profile in a Ultrasound Image Segmentation based on active contour model, it is characterised in that comprise
Following steps:
(1) choosing more than 60 width existing ultrasonoscopy composition training set, described existing ultrasonoscopy contains super with to be split
The tumor that acoustic image is identical, extracts 24 stricture of vaginas of tumor region from training set with gray level co-occurrence matrix in every width ultrasonoscopy
Reason feature, in utilizing mathematical linear to return minimum distance method obtain the standard value of 24 textural characteristics, and and standard ellipse
In empirical value, ultrasonoscopy to be measured, tumor priori size constitutes a standard vector A together0;
(2) ultrasonoscopy to be split is carried out the pretreatment of anisotropy parameter and histogram equalization;
(3) (2) gained image is carried out dynamic threshold segmentation, and extract all Guan Bi wheels produced in dynamic threshold segmentation result
Exterior feature, each closed contour forms a corresponding subimage;
(4) utilize gray level co-occurrence matrix to extract 24 textural characteristics of each subimage, calculate ellipse fitting result parameter with
And subimage closed contour interior pixels number, and constitute a vectorial A by these 26 data1;
(5) each subimage gained vector A is calculated1Standard vector A with tumor region training gained0Distance D so that D is
Little vectorial A1Closed contour in corresponding subimage is the initial profile of tumor in ultrasonoscopy to be split;
Described 24 textural characteristics include image 0 °, 45 °, 90 °, contrast on 135 ° of four directions, energy, dependency, inverse
Difference square, entropy and six features of non-similarity.
In method the most according to claim 1, it is characterised in that in step (1), standard ellipse empirical value is defined as 0.7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410141568.5A CN103914845B (en) | 2014-04-09 | 2014-04-09 | The method obtaining initial profile in Ultrasound Image Segmentation based on active contour model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410141568.5A CN103914845B (en) | 2014-04-09 | 2014-04-09 | The method obtaining initial profile in Ultrasound Image Segmentation based on active contour model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103914845A CN103914845A (en) | 2014-07-09 |
CN103914845B true CN103914845B (en) | 2016-08-17 |
Family
ID=51040501
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410141568.5A Expired - Fee Related CN103914845B (en) | 2014-04-09 | 2014-04-09 | The method obtaining initial profile in Ultrasound Image Segmentation based on active contour model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103914845B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104504689B (en) * | 2014-12-11 | 2017-06-06 | 深圳先进技术研究院 | Adventitia dividing method and device in ivus image |
CN105719278B (en) * | 2016-01-13 | 2018-11-16 | 西北大学 | A kind of medical image cutting method based on statistics deformation model |
CN105931226A (en) * | 2016-04-14 | 2016-09-07 | 南京信息工程大学 | Automatic cell detection and segmentation method based on deep learning and using adaptive ellipse fitting |
CN112166474A (en) * | 2018-05-16 | 2021-01-01 | 皇家飞利浦有限公司 | Automated tumor identification during surgery using machine learning |
CN111242965B (en) * | 2020-01-10 | 2023-04-07 | 西安电子科技大学 | Genetic algorithm-based breast tumor contour dynamic extraction method |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1793350A1 (en) * | 2005-12-01 | 2007-06-06 | Medison Co., Ltd. | Ultrasound imaging system and method for forming a 3D ultrasound image of a target object |
CN101599174A (en) * | 2009-08-13 | 2009-12-09 | 哈尔滨工业大学 | Method for outline extraction of level set medical ultrasonic image area based on edge and statistical nature |
CN102509286A (en) * | 2011-09-28 | 2012-06-20 | 清华大学深圳研究生院 | Target region sketching method for medical image |
CN102855483A (en) * | 2011-06-30 | 2013-01-02 | 北京三星通信技术研究有限公司 | Method and device for processing ultrasonic images and breast cancer diagnosis equipment |
CN103413307A (en) * | 2013-08-02 | 2013-11-27 | 北京理工大学 | Method for image co-segmentation based on hypergraph |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2954986B1 (en) * | 2010-01-05 | 2012-02-10 | St Microelectronics Grenoble 2 | METHOD FOR DETECTION OF CONTOUR ORIENTATION. |
-
2014
- 2014-04-09 CN CN201410141568.5A patent/CN103914845B/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1793350A1 (en) * | 2005-12-01 | 2007-06-06 | Medison Co., Ltd. | Ultrasound imaging system and method for forming a 3D ultrasound image of a target object |
CN101599174A (en) * | 2009-08-13 | 2009-12-09 | 哈尔滨工业大学 | Method for outline extraction of level set medical ultrasonic image area based on edge and statistical nature |
CN102855483A (en) * | 2011-06-30 | 2013-01-02 | 北京三星通信技术研究有限公司 | Method and device for processing ultrasonic images and breast cancer diagnosis equipment |
CN102509286A (en) * | 2011-09-28 | 2012-06-20 | 清华大学深圳研究生院 | Target region sketching method for medical image |
CN103413307A (en) * | 2013-08-02 | 2013-11-27 | 北京理工大学 | Method for image co-segmentation based on hypergraph |
Also Published As
Publication number | Publication date |
---|---|
CN103914845A (en) | 2014-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103914845B (en) | The method obtaining initial profile in Ultrasound Image Segmentation based on active contour model | |
Hou et al. | Unsupervised histopathology image synthesis | |
CN112270660B (en) | Nasopharyngeal carcinoma radiotherapy target area automatic segmentation method based on deep neural network | |
CN106204587A (en) | Multiple organ dividing method based on degree of depth convolutional neural networks and region-competitive model | |
CN110310287A (en) | It is neural network based to jeopardize the automatic delineation method of organ, equipment and storage medium | |
CN105913432A (en) | Aorta extracting method and aorta extracting device based on CT sequence image | |
CN102663416B (en) | Segmentation method of viscera and internal blood vessels thereof in surgical planning system | |
CN104809723A (en) | Three-dimensional liver CT (computed tomography) image automatically segmenting method based on hyper voxels and graph cut algorithm | |
CN106097298A (en) | The coronary artery automatic segmentation divided based on spherical space and anatomic landmarks method | |
CN104933711A (en) | Automatic fast segmenting method of tumor pathological image | |
CN109508787A (en) | Neural network model training method and system for ultrasound displacement estimation | |
CN102682449A (en) | Automatic segmentation method based on self-adaptive external force level set for magnetic resonance images (MRIs) of soft tissue and realization method thereof | |
CN102496156A (en) | Medical image segmentation method based on quantum-behaved particle swarm cooperative optimization | |
CN107507189A (en) | Mouse CT image kidney dividing methods based on random forest and statistical model | |
Wu et al. | Ultrasound image segmentation method for thyroid nodules using ASPP fusion features | |
CN106127794A (en) | Based on probability FCM algorithm MRI tumor image dividing method and system | |
WO2017062453A1 (en) | Image segmentation of organs depicted in computed tomography images | |
CN108154176A (en) | A kind of 3D human body attitude algorithm for estimating for single depth image | |
CN103793916B (en) | Fibroid ultrasonic image division method in a kind of HIFU Treatment | |
Furuzuki et al. | A detection method for liver cancer region based on faster R-CNN | |
CN102881021A (en) | Aortic valve ultrasonic image segmentation method based on probability distribution and continuous maximum flow | |
CN105225234A (en) | Based on the lung tumor identification method of support vector machine MRI Iamge Segmentation | |
Cruz-Aceves et al. | Multiple active contours guided by differential evolution for medical image segmentation | |
Li et al. | A deeply supervised convolutional neural network for brain tumor segmentation | |
CN104268878B (en) | Probe in detecting and localization method are moved in a kind of ultrasonoscopy |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160817 Termination date: 20200409 |