CN101599174A - Method for outline extraction of level set medical ultrasonic image area based on edge and statistical nature - Google Patents
Method for outline extraction of level set medical ultrasonic image area based on edge and statistical nature Download PDFInfo
- Publication number
- CN101599174A CN101599174A CNA2009100716415A CN200910071641A CN101599174A CN 101599174 A CN101599174 A CN 101599174A CN A2009100716415 A CNA2009100716415 A CN A2009100716415A CN 200910071641 A CN200910071641 A CN 200910071641A CN 101599174 A CN101599174 A CN 101599174A
- Authority
- CN
- China
- Prior art keywords
- image
- profile
- level set
- edge
- statistical nature
- 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.)
- Pending
Links
Images
Abstract
The invention provides a kind of level set medical ultrasonic image area profile extraction method based on edge and statistical nature, the present invention is directed to the characteristics that the medical ultrasonic image contrast is low, signal to noise ratio (S/N ratio) is little, designed a kind of level set region contour extracting method based on edge and statistical nature.This method is divided three steps: the first step, region contour slightly extracts, and finishes image negate, self-adapting Gaussian function background inhibition, uses otsu image automatic threshold algorithm that image is converted into bianry image, image degrease perturbation operation, closed region profile extraction work; Second step was adopted a kind of selective anisotropic medical ultrasonic image smoothing algorithm, and original image is carried out pre-service; The level set image-region profile that the 3rd step was based on edge and statistical nature accurately extracts.Experimental result shows that the inventive method compares with existing method, can obtain more precise partition result.
Description
(1) technical field
The present invention relates to the Medical Ultrasound Image Processing field, relate in particular to a kind of region contour extracting method of medical ultrasonic image.
(2) background technology
Medical ultrasonic image is being brought into play important effect as one of important medical image form in clinical diagnosis and medical research.In clinical, the morphological feature of area-of-interest and area thereof or volume all are important diagnostic information in the ultrasonoscopy diagnosis, and the calculating of these information depends on the profile in zone.Therefore, extract the key link that region contour is the medical ultrasonic image diagnosis.Yet, because ultrasonic signal decay, speckle, shade, leakage signal and because the defectives such as edge loss that the directivity of Image Acquisition causes cause in the image obscurity boundary between each tissue unclear, thereby make the doctor relatively more difficult to area-of-interest profile identification in the ultrasonoscopy, simultaneously, doctor's manual drawing area-of-interest profile is not only loaded down with trivial details but also have subjectivity, and it is significant that therefore research has high-quality profile extraction method.
The region contour of image extracts inseparable with image segmentation, image segmentation is exactly that image is divided into some treatment technologies that the zone of certain sense is arranged, promptly some characteristic according to image is divided into the different zone of SOME PROPERTIES to image, at each intra-zone identical or close characteristic is arranged, and the characteristic of adjacent area is inequality.According to pixel characteristic in the target area inner character consistance and interregional borderline uncontinuity, medical ultrasonic image commonly used is cut apart employing based on the edge-detected image dividing method with based on image partition method (the J.Alison Noble of provincial characteristics, Djamal Boukerroui.Ultrasound ImageSegmentation:A Survey.IEEE Trans.on Medical Imaging.2006,25:987-1010).Dividing method based on rim detection mainly is the closed boundary of seeking area-of-interest, generally all adopt various differentiating operators to realize, as various gradient operators such as Roberts operator, Sobel operator, Laplace operator and wavelet transformation etc., affected by noise bigger owing to differentiate, so quality is not very high when these class methods are cut apart ultrasonoscopy; Image partition method based on provincial characteristics is that general image is divided into the overlapping zonule of some complementations, certain similarity that makes each regional interior pixel is greater than the similarity between the zone, this type of dividing method comprises threshold method, region-growing method, Watershed algorithm, EM algorithm, Mathematical Morphology Method and statistical nature method etc., these class methods are in carrying out processing procedure, all pixels are all carried out computing, often, therefore be unfavorable for effective extraction at edge with edge fog.
In recent years on the basis that traditional pixel characteristic is cut apart, dividing method based on deformation model has also been proposed, comprise the Snakes model, Optimization Model, how much and statistical nature model, physical model, (M.Alem á n such as random field models, P.Alem á n, L.Alvarez.Semiautomatic Snake-based Segmentation of Solid Breast Nodules onUltrasonography.Lectures Notes in Computer Science (EuroCast 2005), 2005,3634:467-472), and utilize the effective numerical method of level set on solution curve evolvement problem to realize method (the Yu-Len Huang that extracts with profile of cutting apart of image, Yu-Ru Jiang, Dar-Ren Chen.Level Set Contouring forBreast Tumor in Sonography.Journal of Digital Imaging.2007,20 (3): 238-247) etc., but segmentation effect is still good inadequately.
In a word, for the automatic extraction of ultrasonic image area profile, the existing methods quality is not high, can't satisfy the actual needs of present medical diagnosis.
(3) summary of the invention
The object of the present invention is to provide a kind of on the basis of traditional images dividing method, the advantage of performance level set aspect the curve evolvement numerical solution can access the method for outline extraction of level set medical ultrasonic image area based on edge and statistical nature that high-quality ultrasonic image area profile extracts the result.
The object of the present invention is achieved like this: it comprises three parts: slightly extract based on the region contour of original image (1); (2) based on the ultrasonoscopy pre-service of original image; (3) the level set image-region profile based on edge and statistical nature of bound fraction (1) and part (2) extracts.The image-region profile slightly extracts, image pre-service, image-region profile are smart extracts, the image-region profile is thick to be extracted according to the characteristics of selection area interested (as tumour etc.) in ultrasonoscopy, adopt image negate, Gaussian function background inhibition, binaryzation, morphology opening operation to go to disturb, need not the initialization Level Set Method and extract the regional coarse contours of extraction such as profile, as the initial value of the smart extraction algorithm of follow-up profile; A kind of selective anisotropic medical ultrasonic image smoothing algorithm is adopted in the image pre-service; The image-region profile is smart to extract the level set region contour extraction method that adopts edge and statistical nature to combine.
(1) first: region contour slightly extracts.
The thick purpose of extracting of region contour is to do initialization for the meticulous extraction of back facial contour, and it is divided into following five steps:
1) image negate
In ultrasonoscopy, mostly present echoless or low echo property characteristic such as area-of-interests such as ventricle, tumor of breast (being selection area).The purpose of image negate is for these zones present high brightness, and tissue on every side becomes low-light level, for next step is prepared.
The negate formula is as follows:
μ
2(i,j)=255-μ
1(i,j) (1)
I in the formula, j are respectively the horizontal ordinate and the ordinates of certain pixel in the image, μ
1, μ
2It is respectively before the conversion and the image after the conversion.
2) self-adapting Gaussian function background suppresses
The effect of Adaptive Suppression function is to make the area-of-interest that is positioned at picture centre present high brightness, and deep pixel is in low-light level even is zero.So-called self-adaptation is meant the size automatically adjusting of inhibition function window width according to area-of-interest.The present invention has designed a kind of Adaptive Suppression function based on Gaussian function, function expression as shown in the formula:
P is the coordinate of certain some pixel in the image in the formula,
Be the centre coordinate of image,
σ
x 2, σ
y 2Be respectively the axial variance of image x, y, and σ
x=w/4, σ
y=h/4, wherein w, h are respectively the width and the height of image.
After the adaptive Gauss inhibition, area-of-interests such as image tumour will be than surrounding tissue brightness height.
3) utilization otsu image automatic threshold algorithm is converted into bianry image to image
The otsu algorithm is the automatic selection algorithm of a kind of image threshold of classics, and this algorithm is asked for the maximal value (being maximum between-cluster variance) of following formula by exhaustive method, thereby determines the threshold value of image:
ω in the formula
1(t), ω
2(t) being image is divided into the probability of two time-likes, u by threshold value t
1(t), u
2(t) be the average of two classes respectively.
Adopt otsu automatic threshold method to calculate the threshold value of image, thereby image is become bianry image.
4) operation such as image degrease, other material interference
From view picture ultrasonoscopy shearing area-of-interest the time, inevitably can around the area-of-interest profile, attach the image that fat or other tissue are arranged.These tissues belong to the dark hypoecho tissue, brightness is lower in medical ultrasonic image, behind the image inversion operation, the image of these tissues will present high brightness, though suppressing operation its brightness meeting of back through adaptive Gauss decreases, yet behind otsu adaptive thresholding Value Operations, wherein a part still might keep.For the interference that prevents that these from organizing bright spot, adopt the method for morphology opening operation to remove these bright spots, promptly consider ultrasonoscopy resolution, to the bright spot of those areas, all remove as interference less than certain pixel count.
5) the closed region profile extracts
Adopt a kind of initialization Level Set Method that need not to extract profile, can obtain the area-of-interest profile of a sealing, its Evolution Equation is:
Div is a divergence operator in the formula,
Be gradient operator, δ (φ) is a Dirac function, and parameter lambda is positive weight constant,
* be convolution algorithm, standard deviation is the Gaussian function G of σ
σ(x)=(2 π σ)
-1/2Exp (| x|
2/ 4 σ).
This Level Set Method has overcome traditional level set needs continuous initialization level set function in evolution process shortcoming, compares with traditional Level Set Method to have that speed is fast, the simple advantage of Numerical Implementation.
(2) second portion: ultrasonoscopy pre-service.
Because ultrasonoscopy structure more complicated, signal to noise ratio (S/N ratio) are low, the meticulous profile that directly extracts area-of-interests such as tumour is difficulty relatively, so earlier ultrasonoscopy is done pre-service.
At the characteristics of medical ultrasonic image, adopt a kind of selective anisotropic medical ultrasonic image smoothing algorithm as image pre-processing method, its expression formula is:
In the formula
Be gradient operator, * is a convolution algorithm, and div is a divergence operator.Spread function g (| s|) be a monotonic decreasing function, when | s| → ∞, g (| s|) → 0, as the control rate of propagation, to strengthen the image border.Usually spread function is elected as
K wherein is the Grads threshold parameter.G in the formula
σ(x)=(2 π σ)
-1/2Exp (| x|
2/ 4 σ).Because
Essence be image to be done based on the smoothing processing of Gaussian function to remove noise earlier, the Grad of the image after the computing then is according to Grad decision image smoothing degree.
(3) third part: the level set image-region profile based on edge and statistical nature extracts.
This part is that the essence of image outline is extracted.Elder generation is embedded in the coarse contour image that extracts through first finishing the pretreated image of second portion, as the smart level set initialization of extracting of third part, begins the level set region contour extraction of third part based on edge and statistical nature then.
Be based on the curve differentiation method of geometric deformation model based on the level set region contour extracting method of edge and statistical nature.Suppose to develop curve C (s, t)=[X (s, t), Y (s, t)], wherein s is the arbitrary parameter variable, t is time (R.Malladi, J.A.Sethian, B.C.Vemuri.Shape Modeling with Front Propagation:Level Set Approach.IEEE Trans.on Pattern Analysis and Machine Intelligence, 1995,17 (2): 158-175).If developing the interior of curve is N to unit normal vector, curvature is k, and then curve can be represented with following partial differential equation along the evolution process of its unit normal vector direction:
V in the formula (k) is a velocity function, every evolution speed on the decision curve.
The basic thought of geometric deformation model is that speed of deformation (with curvature and/or constant) and view data are combined, and evolving speed is stopped at object edge.How according to view data and geometrical property design proper speed function, make contour curve develop as early as possible and stop at the core content that is based on geometry deformation model image edge extracting on the object boundary exactly.The geometry deformation model is generally realized (R.Malladi with Level Set Method, J.A.Sethian, B.C.Vemuri.Shape Modeling withFront Propagation:Level Set Approach.IEEE Trans.on Pattern Analysis and Machine Intelligence, 1995,17 (2): 158-175).
Suppose a level set function φ (x, y, t), its zero level collection be C (s, t), then have φ (C (and s, t), t)=0, ask its differential, and utilize the chain type rule to obtain t:
Be without loss of generality, can suppose φ in the zero level collection for negative, for just, the interior of level set curve to unit normal vector is so outside the zero level collection
Wushu (6) is brought formula (7) into, so have:
Curvature k at zero level Ji Chu is:
Formula (8) is the fundamental equation of level set, establishes F=V (k), and function F just is based on the curve evolving speed function of curvature so, and it depends on:
1. the forward position local characteristics that develops of curve for example develops the local curvature of curve;
2. with the relevant external parameter of input data, for example gradient;
3. Fu Jia propagation item.
At the characteristics of ultrasonoscopy, synthetic image marginal information and local texture statistics feature, level set equation speed develops function and elects as:
F=g(x,y)[λ(k(x,y)-k)-probn(x,y)] (10)
Wherein g (x, y)=(1-η (x, y))
2, (x y) is point (x, the edge confidence degree of y) locating to η.Based on the rim detection of edge confidence degree than have based on traditional gradient operator edge detection method calculate simple, efficient is high, to advantages such as insensitive for noise (Meer Peter, Georgescu Bogdan.Edge Detection with Embedded Confidence.IEEE Trans.Pattern Analysis andMachine Intelligence.2001,23 (12): 1351-1365).Because (x, y)<1, (x is y) just more near 1, so (x y) can be very little, and the speed of differentiation will approach zero for g when curve develops edge near edge η more for 0<η.K is a mean curvature, and λ is positive parameter.
Probn (x, y) represent current point (x is the possibility size of area-of-interest outside organization y), and it is drawn by following formula:
Wherein p be wait to classify a little coordinate (x, y), mapping function
β is positive adjusting parameter, and f () is the support vector machine decision function (V.Vapnik.Statistical Learning Theory.Wiley, 1998) based on gaussian kernel function, and its expression formula is:
In the formula (11)
Be nine texture feature vectors, they are respectively mean square deviation coefficients (1) and based on the contrast of grey scale difference image, angle second moment, average, entropy (0 °, 90 ° both direction totally 8).Its expression formula is:
1) mean square deviation coefficient
2) contrast
3) the angle second order is apart from (energy)
4) average
5) entropy
Formula (11) finally is:
Wherein K () is a gaussian kernel function, promptly
, λ
i, ω
0Be the parameter that obtains in the training stage in order to the structure decision function, y
iBe the v of i support vector
iThe class label, N
sIt is the number of support vector.
U, σ are respectively the average and the variances of image in the window.
p
Δ(i) distribution probability of presentation video μ difference, it finds the solution as follows: establish (x, y) be in the image a bit, this point and the point that faces mutually (x+dx, grey scale difference y+dy) is defined as:
Δμ(x,y)=μ(x,y)-μ(x+dx,y+dy) (19)
Obtain the difference value of every of entire image, as the brightness of image and do normalized, the grey scale difference image Δ μ that composing images μ is ordered according to the histogram of difference image, obtains p this difference result
Δ(i), 0≤i≤L-1, L are the maximum possible difference values.If i is less, but bigger frequency values is arranged, this situation represents that texture is more coarse; If histogram of difference is more smooth, illustrate that then texture is more careful.Generally when calculating based on the grey scale difference parametric texture, do 0 °, 90 ° both direction grey scale differences, get dx=1, dy=1, the textures windows size elects 15 * 15 as, like this every bit is had 9 proper vectors.
As seen from formula (18), when the probn of certain some value more near 1 the time, just think the tissue of area-of-interest outside.If approach 0, this point belongs to area-of-interest.See (10) formula now again, establish level set initialization curve in area-of-interest the time, curve begins to develop, and at this moment function probn approaches 0, and only ((x y)-k) decides k evolving speed, and speed is bigger by curvature part λ.When function probn approached 1, the curve evolving speed was by curvature that develops curve and the decision of probn difference, and evolving speed reduces.By the method for adding image local textural characteristics, the evolving speed that develops curve is not only determined by the edge at this place, also relevant with the image property that develops wavefront place image-region, the wavefront differentiation in weak edge also can be stopped.
Of low quality at present existing ultrasonic image area profile extraction method, can't satisfy the problem of the actual needs of medical diagnosis, the present invention proposes a kind of region contour extracting method of medical ultrasonic image, this method is considered the characteristics of medical ultrasonic image, organically combining based on rim detection with based on the method for provincial characteristics, on the basis of traditional images dividing method, the advantage of performance level set aspect the curve evolvement numerical solution obtained high-quality ultrasonic image area profile and extracted the result.
The medical ultrasonic image area profile extraction method that the present invention proposes can be divided into three parts: (1) region contour slightly extracts; (2) ultrasonoscopy pre-service; (3) the level set image-region profile based on edge and statistical nature extracts.Wherein first and second portion are based on all that original image carries out, and third part is to realize on first and the synthetic result's of second portion basis.Can carry out first's operation earlier, also can carry out the second portion operation earlier, perhaps two parts carry out simultaneously.
The above concrete steps of the medical ultrasonic image area profile extraction method that proposes for the present invention.The present invention combines rim detection and based on the method for provincial characteristics, and the application level collection carries out the curve evolvement numerical solution, therefore can access higher-quality result images.
(4) description of drawings
Fig. 1 is the method for outline extraction of level set medical ultrasonic image area structured flowchart based on edge and statistical nature;
Fig. 2 is the thick leaching process synoptic diagram of profile;
Fig. 3 is embedded into after the pre-service result schematic diagram in the image for the ultrasonoscopy pre-service and a thick profile that extracts;
Fig. 4 extracts the net result synoptic diagram for the area-of-interest profile;
Fig. 5 is a segmentation result method for measuring similarity synoptic diagram.
(5) embodiment
The present invention is further illustrated with concrete embodiment below in conjunction with accompanying drawing:
Among Fig. 1: 101 slightly extract, 102 for being the meticulous extraction of level set region contour based on edge and statistical nature based on anisotropic ultrasonoscopy pre-service, 103 for the area-of-interest profile; Among Fig. 2: 201 for the tumor of breast original image, 202 for image inversion operation result, 203 for the image result of application self-adapting Gaussian function after handling, 204 for slightly extracting the result for profile for the image, 206 after degrease and the interference of other materials based on the bianry image of OTSU adaptive threshold, 205; Among Fig. 3: 301 are embedded into result in the pretreated image for coarse contour for ultrasonoscopy pre-service result, 302; Among Fig. 4: 401 is that original tumor of breast image, 402 is for the level set medical ultrasonic image area profile based on edge and statistical nature extracts the result, 403 contour extraction methods that propose for the present invention extract result's (redness) and medical expert's manual drawing result (white) contrast; Among Fig. 5: the profile that 501 profiles of cutting apart for the inventive method, 502 are cut apart for the inventive method for the profile of medical expert's manual drawing, 503 and the profile phase of medical expert's manual drawing are relatively, the profile that the area of lost part, 504 is cut apart for the inventive method and the profile phase of medical expert's manual drawing are relatively, the profile that the area of lap, 505 is cut apart for the inventive method and the profile phase of medical expert's manual drawing relatively have more the area of part.
In conjunction with Fig. 1, present embodiment is as follows based on the specific implementation step of the level set medical ultrasonic image area profile extraction method of edge and statistical nature:
(1) the region of interest area image is carried out cutting.The purpose of reducing is to reduce the calculated amount that profile extracts automatically, the method of reducing is to comprise whole image of interest zone, and Edge Distance area-of-interest profile is unsuitable excessive, is generally hundreds of to about several ten thousand pixels, is chosen as 96 * 64 pixels in the present embodiment;
(2) extract the initial profile of area-of-interest, the tumor of breast original image for after reducing shown in 201 among Fig. 2 utilizes formula (1) to carry out the image negate, obtains image 202; Pass through formula (2) Gaussian function background again and suppress, obtain image 203; Adopt otsu automatic threshold method to calculate the threshold value of image, image is become bianry image, shown among Fig. 2 204; Adopt morphology opening operation method to remove the discontinuous bright spot at random of image again, in the present embodiment area is thought that less than the bright spot of 5 pixels bright spot is removed at random, reach the purpose of interference such as degrease, obtain image 205; Utilize formula (4) to realize that the sealing image outline extracts, and can obtain the initial profile of a tumour, shown among Fig. 2 206;
(3) the tumor of breast original image shown in 201 in to Fig. 2 utilizes formula (5) to carry out the anisotropy image smoothing, obtains the pretreated image shown in 301 in the accompanying drawing 3;
(4) initial profile 206 of step (2) gained is embedded in the pretreated image 301, obtains image 302;
(5) utilize formula (19) to calculate the difference image of pretreated image;
(6) utilize formula (13)-(18) to calculate the proper vector of region of interest tract tissue
(7) the improvement iteration cross term verification method [7] that proposes of utilization Carl Staelin is determined the parameter of gaussian kernel function support vector machine (C γ), is got C=512 promptly 2 in the present embodiment
9, γ=0.03125 promptly 2
-5
(8) using formula (10) is cut apart image, extracts the profile of final area-of-interest.
Describe for beneficial effect of the present invention below in conjunction with embodiment.Picked at random three width of cloth malignant galactophore supersonic tumor images and three width of cloth breast ultrasound benign tumour images.The similarity measurement of the profile that the employing aforementioned embodiments is extracted and the tumour profile of medical expert's manual drawing is estimated the profile extraction effect.The medical expert comprises that radiation expert and clinical expert form, so these manual profiles have certain authority, can estimate the effect of automatic extraction.As shown in Figure 5,501 profiles of cutting apart for the inventive method, 502 is the profile of medical expert's manual drawing.The profile phase of 503 profiles of cutting apart for the inventive method and medical expert's manual drawing relatively, the area of lost part is designated as S2; 504 is the area of lap, is designated as S3; 505 for having more the area of part, is designated as S1, and then similarity measurement is defined as follows:
The value of O when O=0, illustrates with two profiles and do not occur simultaneously between 0 and 1, when O=1, illustrates that the profile cut apart and the profile of medical expert's manual drawing overlap fully.The similarity measurement result of aforesaid 6 width of cloth tumor images is respectively: 97.84%, 98.45%, 98.33%, 92.12%, 96.62%, 98.39%, and average similarity tolerance is 96.96%.Adopt aforementioned embodiments that 240 breast ultrasound evaluate image (wherein malignant tumour 128 width of cloth, optimum 112 width of cloth) are carried out the tumour profile and extract research, average similarity tolerance is 91.37%.
Document Yu-Len Huang, Yu-Ru Jiang, Dar-Ren Chen.Level Set Contouring for BreastTumor in Sonography.Journal of Digital Imaging.2007,20 (3): the tumour profile of mammary gland is extracted in the 238-247 employing based on smothing filtering and traditional level set contour extraction method, (3 width of cloth are the benign tumour image to have provided 6 width of cloth images in evaluation of algorithm, 3 width of cloth are the malignant tumour image in addition) profile extraction result, similarity measurement is respectively: 95%, 86%, 95%, 95%, 94%, 91%, this 6 width of cloth image similarity tolerance average is 92.67%.Adopt 118 width of cloth image libraries (pernicious 34 width of cloth, optimum 84 width of cloth) that the algorithm that this article proposes is assessed simultaneously in the document, average similarity tolerance is 87.64%.
From above analysis as can be known, the region contour extracting method of the medical ultrasonic image of the present invention's proposition can obtain higher profile and extract quality.
Claims (5)
1. level set medical ultrasonic image area profile extraction method based on edge and statistical nature, it is characterized in that: this method is divided three steps: the image-region profile slightly extracts, ultrasonoscopy pre-service, image-region profile are smart extracts, the image-region profile slightly extracts according to the characteristics of original image selection area in ultrasonoscopy, adopt image negate, Gaussian function background inhibition, binaryzation, morphology opening operation to go to disturb, need not the initialization Level Set Method and extract the regional coarse contours of extraction such as profile, as the initial value of the smart extraction algorithm of follow-up profile; A kind of selective anisotropic medical ultrasonic image smoothing algorithm is adopted in the image pre-service; The image-region profile is smart to extract the level set region contour extraction method that adopts edge and statistical nature to combine.
2. according to right 1 described level set medical ultrasonic image area profile extraction method based on edge and statistical nature, it is characterized in that: the level set region contour extraction method that described edge and statistical nature combine, calculate the texture statistics feature of every some region on the profile earlier, and be input in the support vector machine that has trained, obtain the symbolic distance between these texture statistics features and classification lineoid, introduce a mapping function simultaneously, these symbolic distances are mapped in the interval, characterize the possibility that this point belongs to normal structure, the curve evolving speed of determining this point in conjunction with the edge confidence degree and the curvature of this point is realized accurately extracting automatically of interesting image regions profile at last.
3. according to right 1 described level set medical ultrasonic image area profile extraction method based on edge and statistical nature, it is characterized in that: the thick extraction step of described image-region profile comprises following five steps:
1) image negate
In ultrasonoscopy, mostly present echoless or low echo property characteristic such as area-of-interests such as ventricle, tumors of breast; The negate formula is as follows:
μ
2(i,j)=255-μ
1(i,j) (1)
I in the formula, j are respectively the horizontal ordinate and the ordinates of certain pixel in the image, μ
1, μ
2It is respectively before the conversion and the image after the conversion;
2) self-adapting Gaussian function background suppresses
The effect of Adaptive Suppression function is to make the area-of-interest that is positioned at picture centre present high brightness, and deep pixel is in low-light level even is zero, sets a kind of Adaptive Suppression function based on Gaussian function, function expression as shown in the formula:
P is the coordinate of certain some pixel in the image in the formula,
Be the centre coordinate of image,
σ
x 2, σ
y 2Be respectively the axial variance of image x, y, and σ
x=w, 4, σ
y=h/4, wherein w, h are respectively the width and the height of image;
3) utilization otsu image automatic threshold algorithm is converted into bianry image to image
The otsu algorithm is the automatic selection algorithm of a kind of image threshold of classics, and the maximal value that this algorithm is asked for following formula by exhaustive method is a maximum between-cluster variance, thereby determines the threshold value of image:
ω in the formula
1(t), ω
2(t) being image is divided into the probability of two time-likes, u by threshold value t
1(t), u
2(t) be the average of two classes respectively; Adopt otsu automatic threshold method to calculate the threshold value of image, thereby image is become bianry image;
4) operation such as image degrease, other material interference
Adopt the method for morphology opening operation to remove subsidiary these bright spots of image that fat or other tissue are arranged of region contour on every side, promptly consider ultrasonoscopy resolution,, all remove as interference to the bright spot of those areas less than certain pixel count;
5) the closed region profile extracts
Adopt a kind of initialization Level Set Method that need not to extract profile, can obtain the area-of-interest profile of a sealing, its Evolution Equation is:
4. the level set medical ultrasonic image area profile extraction method based on edge and statistical nature according to claim 1, it is characterized in that: described ultrasonoscopy pre-treatment step comprises: at the characteristics of medical ultrasonic image, adopt a kind of selective anisotropic medical ultrasonic image smoothing algorithm as image pre-processing method, its expression formula is:
In the formula
Be gradient operator, * is a convolution algorithm, and div is a divergence operator.Spread function g (| s|) be a monotonic decreasing function, when | s| → ∞, g (| s|) → 0, as the control rate of propagation, to strengthen the image border.Usually spread function is elected as
K wherein is the Grads threshold parameter, G in the formula
σ(x)=(2 π σ)
-1/2Exp (| x|
2/ 4 σ), because
Essence be image to be done based on the smoothing processing of Gaussian function to remove noise earlier, the Grad of the image after the computing then is according to Grad decision image smoothing degree.
5. the level set medical ultrasonic image area profile extraction method based on edge and statistical nature according to claim 1, it is characterized in that: described level set image-region profile extraction step based on edge and statistical nature comprises: earlier finishing the pretreated image of second portion ultrasonoscopy, be embedded in the coarse contour image that extracts through first, as the smart level set initialization of extracting of third part, begin the level set region contour extraction of third part then based on edge and statistical nature; Be based on the curve differentiation method of geometric deformation model based on the level set region contour extracting method of edge and statistical nature, suppose to develop curve C (s, t)=[X (and s, t), Y (s, t)], wherein s is the arbitrary parameter variable, t is the time; If developing the interior of curve is N to unit normal vector, curvature is k, and then curve can be represented with following partial differential equation along the evolution process of its unit normal vector direction:
V in the formula (k) is a velocity function, every evolution speed on the decision curve.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2009100716415A CN101599174A (en) | 2009-08-13 | 2009-08-13 | Method for outline extraction of level set medical ultrasonic image area based on edge and statistical nature |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2009100716415A CN101599174A (en) | 2009-08-13 | 2009-08-13 | Method for outline extraction of level set medical ultrasonic image area based on edge and statistical nature |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101599174A true CN101599174A (en) | 2009-12-09 |
Family
ID=41420607
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2009100716415A Pending CN101599174A (en) | 2009-08-13 | 2009-08-13 | Method for outline extraction of level set medical ultrasonic image area based on edge and statistical nature |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101599174A (en) |
Cited By (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101866487A (en) * | 2010-06-11 | 2010-10-20 | 中国科学院深圳先进技术研究院 | Method and system for extracting body part from medical image |
CN102034105A (en) * | 2010-12-16 | 2011-04-27 | 电子科技大学 | Object contour detection method for complex scene |
CN102133109A (en) * | 2010-01-25 | 2011-07-27 | 安克生医股份有限公司 | Method for quantifying and imaging features of a tumor |
CN102306373A (en) * | 2011-08-17 | 2012-01-04 | 深圳市旭东数字医学影像技术有限公司 | Method and system for dividing up three-dimensional medical image of abdominal organ |
CN102499711A (en) * | 2011-09-28 | 2012-06-20 | 无锡祥生医学影像有限责任公司 | Three-dimensional or four-dimensional automatic ultrasound image optimization and adjustment method |
CN102663749A (en) * | 2012-03-28 | 2012-09-12 | 汕头大学 | Ultrasonic image lesion segmentation technology based on phase information |
CN102760236A (en) * | 2012-03-20 | 2012-10-31 | 苏州迪凯尔医疗科技有限公司 | Priori shape modeling method based on combined sparse model |
CN103679685A (en) * | 2012-09-11 | 2014-03-26 | 北京三星通信技术研究有限公司 | Image processing system and image processing method |
CN103914845A (en) * | 2014-04-09 | 2014-07-09 | 武汉大学 | Method for acquiring initial contour in ultrasonic image segmentation based on active contour model |
CN104126194A (en) * | 2012-02-27 | 2014-10-29 | 爱克发医疗保健公司 | Image alignment of breast images |
CN104166976A (en) * | 2013-05-16 | 2014-11-26 | 上海联影医疗科技有限公司 | Segmentation method of foreground in three-dimensional image |
CN104331869A (en) * | 2014-11-25 | 2015-02-04 | 南京信息工程大学 | Image smoothing method with combination of gradient and curvature |
CN104470443A (en) * | 2012-07-18 | 2015-03-25 | 皇家飞利浦有限公司 | Method and system for processing ultrasonic imaging data |
CN104484860A (en) * | 2014-12-18 | 2015-04-01 | 南京信息工程大学 | Image smoothing algorithm based on quasi-normal distribution |
CN104599270A (en) * | 2015-01-18 | 2015-05-06 | 北京工业大学 | Breast neoplasms ultrasonic image segmentation method based on improved level set algorithm |
CN105139410A (en) * | 2015-09-11 | 2015-12-09 | 浙江工商大学 | Brain tumor MRI image segmentation algorithm utilizing space section projection |
CN105261024A (en) * | 2015-10-22 | 2016-01-20 | 四川航天系统工程研究所 | Remote sensing image complex water body boundary extraction method based on improved T-Snake model |
CN105374025A (en) * | 2015-09-30 | 2016-03-02 | 上海联影医疗科技有限公司 | Breast image obtaining method and device and mammography system |
CN105574882A (en) * | 2015-12-30 | 2016-05-11 | 中国科学院深圳先进技术研究院 | Lung segmentation extraction method and system based on CT image of chest cross section |
CN105719275A (en) * | 2015-12-10 | 2016-06-29 | 中色科技股份有限公司 | Parallel combination image defect segmentation method |
CN106548183A (en) * | 2016-10-27 | 2017-03-29 | 南京理工大学 | A kind of method that edge extracting is carried out to fuzzy object quickly |
CN106599861A (en) * | 2016-12-20 | 2017-04-26 | 江苏大学附属医院 | Method and apparatus for automatically reading numerical values of ultrasonic wave report |
CN107169975A (en) * | 2017-03-27 | 2017-09-15 | 中国科学院深圳先进技术研究院 | The analysis method and device of ultrasonoscopy |
CN107610149A (en) * | 2017-09-25 | 2018-01-19 | 北京奇虎科技有限公司 | Image segmentation result edge optimization processing method, device and computing device |
CN107657641A (en) * | 2016-07-26 | 2018-02-02 | 聊城大学 | A kind of intelligent stone cutting method based on machine vision |
CN107767385A (en) * | 2017-08-28 | 2018-03-06 | 江苏理工学院 | A kind of carborundum line grain count method and apparatus based on machine vision |
CN108334879A (en) * | 2018-02-11 | 2018-07-27 | 深圳市朗驰欣创科技股份有限公司 | A kind of method for extracting region, system and terminal device |
CN108471995A (en) * | 2015-09-30 | 2018-08-31 | 上海联影医疗科技有限公司 | The system and method for determining breast area in medical image |
CN109255757A (en) * | 2018-04-25 | 2019-01-22 | 江苏大学 | A kind of machine vision places grape cluster carpopodium region segmentation method naturally |
CN109711392A (en) * | 2019-01-24 | 2019-05-03 | 郑州市现代人才测评与考试研究院 | A kind of talent's assessment method based on recognition of face |
CN109813244A (en) * | 2019-01-22 | 2019-05-28 | 湘潭大学 | A kind of automatic speed transmission electrohydraulic module heating power deformation detecting device and detection method |
CN110378912A (en) * | 2019-07-12 | 2019-10-25 | 深圳码隆科技有限公司 | Wrap up detection method, device, computer readable storage medium and computer equipment |
CN110533682A (en) * | 2019-08-30 | 2019-12-03 | 福建省德腾智能科技有限公司 | A kind of image border real time extracting method based on curvature filtering |
CN110853024A (en) * | 2019-11-14 | 2020-02-28 | 北京推想科技有限公司 | Medical image processing method, medical image processing device, storage medium and electronic equipment |
CN111008987A (en) * | 2019-12-06 | 2020-04-14 | 深圳市碧海扬帆科技有限公司 | Method and device for extracting edge image based on gray background and readable storage medium |
US11020094B2 (en) | 2012-07-18 | 2021-06-01 | Koninklijke Philips N.V. | Method and system for processing ultrasonic imaging data |
CN113140013A (en) * | 2020-01-20 | 2021-07-20 | 纬创资通股份有限公司 | Apparatus and method for post-processing of computed tomography |
CN113948173A (en) * | 2021-10-22 | 2022-01-18 | 昆明理工大学 | Medical auxiliary system based on augmented reality and finite element analysis and use method |
CN114581392A (en) * | 2022-02-28 | 2022-06-03 | 山东省人工智能研究院 | Image segmentation method based on deep learning and anisotropic active contour |
CN115035004A (en) * | 2022-04-15 | 2022-09-09 | 腾讯科技(深圳)有限公司 | Image processing method, apparatus, device, readable storage medium and program product |
-
2009
- 2009-08-13 CN CNA2009100716415A patent/CN101599174A/en active Pending
Cited By (71)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102133109A (en) * | 2010-01-25 | 2011-07-27 | 安克生医股份有限公司 | Method for quantifying and imaging features of a tumor |
CN101866487B (en) * | 2010-06-11 | 2014-01-08 | 中国科学院深圳先进技术研究院 | Method and system for extracting body part from medical image |
CN101866487A (en) * | 2010-06-11 | 2010-10-20 | 中国科学院深圳先进技术研究院 | Method and system for extracting body part from medical image |
CN102034105A (en) * | 2010-12-16 | 2011-04-27 | 电子科技大学 | Object contour detection method for complex scene |
CN102306373A (en) * | 2011-08-17 | 2012-01-04 | 深圳市旭东数字医学影像技术有限公司 | Method and system for dividing up three-dimensional medical image of abdominal organ |
CN102499711A (en) * | 2011-09-28 | 2012-06-20 | 无锡祥生医学影像有限责任公司 | Three-dimensional or four-dimensional automatic ultrasound image optimization and adjustment method |
CN102499711B (en) * | 2011-09-28 | 2013-07-10 | 无锡祥生医学影像有限责任公司 | Three-dimensional or four-dimensional automatic ultrasound image optimization and adjustment method |
CN104126194A (en) * | 2012-02-27 | 2014-10-29 | 爱克发医疗保健公司 | Image alignment of breast images |
CN102760236B (en) * | 2012-03-20 | 2014-10-01 | 苏州迪凯尔医疗科技有限公司 | Priori shape modeling method based on combined sparse model |
CN102760236A (en) * | 2012-03-20 | 2012-10-31 | 苏州迪凯尔医疗科技有限公司 | Priori shape modeling method based on combined sparse model |
CN102663749A (en) * | 2012-03-28 | 2012-09-12 | 汕头大学 | Ultrasonic image lesion segmentation technology based on phase information |
CN104470443A (en) * | 2012-07-18 | 2015-03-25 | 皇家飞利浦有限公司 | Method and system for processing ultrasonic imaging data |
US11020094B2 (en) | 2012-07-18 | 2021-06-01 | Koninklijke Philips N.V. | Method and system for processing ultrasonic imaging data |
CN103679685A (en) * | 2012-09-11 | 2014-03-26 | 北京三星通信技术研究有限公司 | Image processing system and image processing method |
CN103679685B (en) * | 2012-09-11 | 2018-03-27 | 北京三星通信技术研究有限公司 | Image processing system and image processing method |
CN104166976A (en) * | 2013-05-16 | 2014-11-26 | 上海联影医疗科技有限公司 | Segmentation method of foreground in three-dimensional image |
CN103914845A (en) * | 2014-04-09 | 2014-07-09 | 武汉大学 | Method for acquiring initial contour in ultrasonic image segmentation based on active contour model |
CN103914845B (en) * | 2014-04-09 | 2016-08-17 | 武汉大学 | The method obtaining initial profile in Ultrasound Image Segmentation based on active contour model |
CN104331869A (en) * | 2014-11-25 | 2015-02-04 | 南京信息工程大学 | Image smoothing method with combination of gradient and curvature |
CN104331869B (en) * | 2014-11-25 | 2017-10-27 | 南京信息工程大学 | The image smoothing method that gradient is combined with curvature |
CN104484860A (en) * | 2014-12-18 | 2015-04-01 | 南京信息工程大学 | Image smoothing algorithm based on quasi-normal distribution |
CN104484860B (en) * | 2014-12-18 | 2017-05-17 | 南京信息工程大学 | Image smoothing method based on quasi-normal distribution |
CN104599270B (en) * | 2015-01-18 | 2017-10-10 | 北京工业大学 | A kind of Ultrasound Image of Breast Tumor dividing method based on improvement level set algorithm |
CN104599270A (en) * | 2015-01-18 | 2015-05-06 | 北京工业大学 | Breast neoplasms ultrasonic image segmentation method based on improved level set algorithm |
CN105139410A (en) * | 2015-09-11 | 2015-12-09 | 浙江工商大学 | Brain tumor MRI image segmentation algorithm utilizing space section projection |
CN105139410B (en) * | 2015-09-11 | 2017-11-10 | 浙江工商大学 | The brain tumor MRI image dividing method projected using aerial cross sectional |
CN105374025A (en) * | 2015-09-30 | 2016-03-02 | 上海联影医疗科技有限公司 | Breast image obtaining method and device and mammography system |
US10636143B2 (en) | 2015-09-30 | 2020-04-28 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for determining a breast region in a medical image |
US11250567B2 (en) | 2015-09-30 | 2022-02-15 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for determining a breast region in a medical image |
US10297024B2 (en) | 2015-09-30 | 2019-05-21 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for determining a breast region in a medical image |
CN108471995A (en) * | 2015-09-30 | 2018-08-31 | 上海联影医疗科技有限公司 | The system and method for determining breast area in medical image |
CN105374025B (en) * | 2015-09-30 | 2018-05-04 | 上海联影医疗科技有限公司 | Breast image acquisition methods and device, mammography system |
CN105261024B (en) * | 2015-10-22 | 2018-02-16 | 四川航天系统工程研究所 | Based on the remote sensing image Complex water body boundary extraction method for improving T Snake models |
CN105261024A (en) * | 2015-10-22 | 2016-01-20 | 四川航天系统工程研究所 | Remote sensing image complex water body boundary extraction method based on improved T-Snake model |
CN105719275A (en) * | 2015-12-10 | 2016-06-29 | 中色科技股份有限公司 | Parallel combination image defect segmentation method |
CN105574882B (en) * | 2015-12-30 | 2018-08-10 | 中国科学院深圳先进技术研究院 | Lung segmentation extracting method based on chest cross section CT images and system |
CN105574882A (en) * | 2015-12-30 | 2016-05-11 | 中国科学院深圳先进技术研究院 | Lung segmentation extraction method and system based on CT image of chest cross section |
CN107657641B (en) * | 2016-07-26 | 2021-09-10 | 聊城大学 | Intelligent stone cutting method based on machine vision |
CN107657641A (en) * | 2016-07-26 | 2018-02-02 | 聊城大学 | A kind of intelligent stone cutting method based on machine vision |
CN106548183B (en) * | 2016-10-27 | 2019-04-12 | 南京理工大学 | A method of edge extracting quickly being carried out to fuzzy object |
CN106548183A (en) * | 2016-10-27 | 2017-03-29 | 南京理工大学 | A kind of method that edge extracting is carried out to fuzzy object quickly |
CN106599861B (en) * | 2016-12-20 | 2019-07-26 | 江苏大学附属医院 | A kind of automatic read method of ultrasound report numerical value and device |
CN106599861A (en) * | 2016-12-20 | 2017-04-26 | 江苏大学附属医院 | Method and apparatus for automatically reading numerical values of ultrasonic wave report |
CN107169975B (en) * | 2017-03-27 | 2019-07-30 | 中国科学院深圳先进技术研究院 | The analysis method and device of ultrasound image |
WO2018176319A1 (en) * | 2017-03-27 | 2018-10-04 | 中国科学院深圳先进技术研究院 | Ultrasound image analysis method and device |
CN107169975A (en) * | 2017-03-27 | 2017-09-15 | 中国科学院深圳先进技术研究院 | The analysis method and device of ultrasonoscopy |
CN107767385A (en) * | 2017-08-28 | 2018-03-06 | 江苏理工学院 | A kind of carborundum line grain count method and apparatus based on machine vision |
CN107610149A (en) * | 2017-09-25 | 2018-01-19 | 北京奇虎科技有限公司 | Image segmentation result edge optimization processing method, device and computing device |
CN107610149B (en) * | 2017-09-25 | 2020-03-03 | 北京奇虎科技有限公司 | Image segmentation result edge optimization processing method and device and computing equipment |
CN108334879A (en) * | 2018-02-11 | 2018-07-27 | 深圳市朗驰欣创科技股份有限公司 | A kind of method for extracting region, system and terminal device |
CN108334879B (en) * | 2018-02-11 | 2020-08-21 | 深圳市朗驰欣创科技股份有限公司 | Region extraction method, system and terminal equipment |
CN109255757A (en) * | 2018-04-25 | 2019-01-22 | 江苏大学 | A kind of machine vision places grape cluster carpopodium region segmentation method naturally |
CN109255757B (en) * | 2018-04-25 | 2022-01-11 | 江苏大学 | Method for segmenting fruit stem region of grape bunch naturally placed by machine vision |
CN109813244B (en) * | 2019-01-22 | 2020-01-07 | 湘潭大学 | Detection device and detection method for thermal deformation of electro-hydraulic module of automatic transmission |
CN109813244A (en) * | 2019-01-22 | 2019-05-28 | 湘潭大学 | A kind of automatic speed transmission electrohydraulic module heating power deformation detecting device and detection method |
CN109711392A (en) * | 2019-01-24 | 2019-05-03 | 郑州市现代人才测评与考试研究院 | A kind of talent's assessment method based on recognition of face |
CN110378912B (en) * | 2019-07-12 | 2021-11-16 | 深圳码隆科技有限公司 | Package detection method and device, computer readable storage medium and computer equipment |
CN110378912A (en) * | 2019-07-12 | 2019-10-25 | 深圳码隆科技有限公司 | Wrap up detection method, device, computer readable storage medium and computer equipment |
CN110533682A (en) * | 2019-08-30 | 2019-12-03 | 福建省德腾智能科技有限公司 | A kind of image border real time extracting method based on curvature filtering |
CN110533682B (en) * | 2019-08-30 | 2023-02-14 | 福建省德腾智能科技有限公司 | Image edge real-time extraction method based on curvature filtering |
CN110853024A (en) * | 2019-11-14 | 2020-02-28 | 北京推想科技有限公司 | Medical image processing method, medical image processing device, storage medium and electronic equipment |
CN110853024B (en) * | 2019-11-14 | 2020-12-22 | 推想医疗科技股份有限公司 | Medical image processing method, medical image processing device, storage medium and electronic equipment |
CN111008987A (en) * | 2019-12-06 | 2020-04-14 | 深圳市碧海扬帆科技有限公司 | Method and device for extracting edge image based on gray background and readable storage medium |
CN111008987B (en) * | 2019-12-06 | 2023-06-09 | 深圳市碧海扬帆科技有限公司 | Method and device for extracting edge image based on gray background and readable storage medium |
CN113140013A (en) * | 2020-01-20 | 2021-07-20 | 纬创资通股份有限公司 | Apparatus and method for post-processing of computed tomography |
CN113948173A (en) * | 2021-10-22 | 2022-01-18 | 昆明理工大学 | Medical auxiliary system based on augmented reality and finite element analysis and use method |
CN113948173B (en) * | 2021-10-22 | 2024-03-22 | 昆明理工大学 | Medical auxiliary system based on augmented reality and finite element analysis and use method |
CN114581392A (en) * | 2022-02-28 | 2022-06-03 | 山东省人工智能研究院 | Image segmentation method based on deep learning and anisotropic active contour |
CN114581392B (en) * | 2022-02-28 | 2022-10-11 | 山东省人工智能研究院 | Image segmentation method based on deep learning and anisotropic active contour |
CN115035004A (en) * | 2022-04-15 | 2022-09-09 | 腾讯科技(深圳)有限公司 | Image processing method, apparatus, device, readable storage medium and program product |
CN115035004B (en) * | 2022-04-15 | 2023-02-10 | 腾讯科技(深圳)有限公司 | Image processing method, apparatus, device, readable storage medium and program product |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101599174A (en) | Method for outline extraction of level set medical ultrasonic image area based on edge and statistical nature | |
CN109840913B (en) | Method and system for segmenting tumor in mammary X-ray image | |
Zhang et al. | Level set evolution driven by optimized area energy term for image segmentation | |
CN102147861A (en) | Moving target detection method for carrying out Bayes judgment based on color-texture dual characteristic vectors | |
CN105335967A (en) | Back-of-hand vein line extraction method based on local maximum between-class variance and mathematical morphology | |
CN104331683B (en) | A kind of facial expression recognizing method with noise robustness | |
CN107403438A (en) | Improve the ultrasonoscopy focal zone dividing method of fuzzy clustering algorithm | |
CN101901476A (en) | SAR image de-noising method based on NSCT domain edge detection and Bishrink model | |
CN102073994B (en) | Ultrasonic medical image speckle noise inhibition method based on multi-scale anisotropic diffusion | |
Cheng et al. | ACCOMP: Augmented cell competition algorithm for breast lesion demarcation in sonography | |
CN104778670A (en) | Fractal-wavelet self-adaption image denoising method based on multivariate statistical model | |
Chen et al. | A new watershed algorithm for cellular image segmentation based on mathematical morphology | |
CN103914829A (en) | Method for detecting edge of noisy image | |
Ananth | An Exclusive Review of Popular Image Processing Techniques | |
CN110458773A (en) | A kind of anisotropy parameter method for processing noise based on edge enhancement operator | |
CN104331864A (en) | Breast imaging processing based on non-subsampled contourlet and visual salient model | |
Jian et al. | Tumor detection in MRI brain images based on saliency computational modeling | |
CN106651781A (en) | Image noise suppression method for laser active imaging | |
CN102402783A (en) | Method for processing spots of three-dimensional ultrasonic image | |
Gupta et al. | A novel hybrid method for segmentation and analysis of brain MRI for tumor diagnosis | |
CN105825514A (en) | Shear wave coefficient probability density distribution-based image edge detection method | |
Kachouie et al. | An elliptical level set method for automatic TRUS prostate image segmentation | |
Gun et al. | A contour detector with improved corner detection | |
JP3383072B2 (en) | Image processing method and image processing apparatus | |
Srilatha et al. | Automated ultrasound ovarian tumour segmentation and classification based on deep learning techniques |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20091209 |