CN102254097A - Method for identifying fissure on lung CT (computed tomography) image - Google Patents
Method for identifying fissure on lung CT (computed tomography) image Download PDFInfo
- Publication number
- CN102254097A CN102254097A CN2011101907972A CN201110190797A CN102254097A CN 102254097 A CN102254097 A CN 102254097A CN 2011101907972 A CN2011101907972 A CN 2011101907972A CN 201110190797 A CN201110190797 A CN 201110190797A CN 102254097 A CN102254097 A CN 102254097A
- Authority
- CN
- China
- Prior art keywords
- lung
- fissure
- segmentation
- plane
- splits
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30061—Lung
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The invention provides a method for automatically detecting and cutting a fissure on a lung CT (computed tomography) image efficiently and accurately based on three-dimensional plane fitting. In the method, a CT image is taken as a group of three-dimensional space point clouds; a lung area is divided into small spherical subdivision bodies, and a fissure plane is discovered in the small subdivision bodies by utilizing the three-dimensional plane fitting method, so that the fissure detection with free-form surface characteristics is converted into detection on planes, and the complexity of the problem is remarkably reduced. The method has the advantage of being not sensitive to noise or abnormal values. Furthermore, in a fissure detection process, fissures (oblique fissure in left and right lungs and horizontal fissure) in different types can be automatically identified by a simple clustering method. Compared with other methods, the method provided by the invention has the characteristics of good accuracy, stability, high efficiency and the like.
Description
Technical field
The present invention relates to based on segmentation three-dimensional planar approximating method, from brand-new angle, the lung of discerning efficiently on the chest CT image splits.
Background technology
Lung interlobar fissure (being called for short lung splits) is the important sign that embodies human lung's structure.The integrity degree that lung is split and the full appreciation of structure, pulmonary disease detects in early days, and classification on PD and the treatment of diseases, all has very big value for clinical application.Therefore, the accurate identification split of lung is very important concerning clinical.It's a pity, the high resolution CT examination comprises a large amount of images usually, allow the expert manually one by one picture to go to indicate that lung splits the locus, place very consuming time, and lung splits and shows to such an extent that be not very clear on the CT image usually, often can't guarantee result's accuracy and consistance.
At present, proposed a lot of lungs in the world wide and split recognition methods, wherein most of method comprises two stages, i.e. initial identification and further optimization improvement.In the specific implementation process, initial detecting/identification is generally by determining that region of interest (ROI:region-of-interest) realizes.Such as, general people such as grade [1] handles by a simple thresholding and discerns area-of-interest; People such as Wiemker [2] utilize Hessian matrix and structure tensor filtering to determine the prime area.In addition in document [5-8], the prime area be by lung split and blood vessel, bronchus between relation on dissecting determine.After determining the prime area, next step is exactly to determine that specifically the locus that lung splits, the method for this respect are by determining straight line or realize at three dimensions searching curved surface or plane in two-dimensional space substantially.Two dimension straight line (2D) method has comprised growth curve method [9-10], and the linear feature of Vanderbrug detects [11], rim detection [8] and Gaussian distribution and/or be mean curvature analysis [12].Consider in identification aspect lung splits, three-dimension curved surface or plane (3D) than two-dimentional straight line (2D) more comprehensively, accurately, people such as Ukil [7] have proposed a kind of ridge degree mensuration [13] of the 3D of having graph search function and sought the plane in region of interest; Generally wait people [1] to propose the computational geometry method that a kind of usefulness comprises the level and smooth and expansion Gaussian image (EGI:extended Gaussian image) of Laplce; People such as Rikxoort [3-4] propose to utilize second-order logic to calculate point are combined into the plane; Kuhnigk[6] etc. the people a kind of interactive three-dimensional watershed algorithm has been proposed.Different with these methods, open to wait a people [14-16] to propose a kind of masterplate matching process and seek lung and split based on picture library.
Summary of the invention
The lung interlobar fissure has two characteristics on the CT image, and the firstth, split lung tissue on every side with lung and compare, lung splits relative higher brightness, and this also is why human eye can be seen the reason that lung splits; The secondth, lung splits higher density in the image space based on voxel.In order fully to excavate these two characteristics, this invention has proposed a kind of based on plane fitting, discerns the method that lung splits on the CT image efficiently and accurately.
Consider that three-dimensional free surface can accurately approach with a lot of little planes, we propose three-dimensional lung areas is subdivided into a lot of little spheries segmentation bodies (overlapping mutually between these segmentation bodies), and the identification problem split of lung just is converted to a series of plane problems of searching in each spheroid like this.This method has comprised five basic steps (as Fig. 1 to Fig. 5): (1) lung areas is cut apart (Fig. 1), (2) lung areas space segmentation (Fig. 2), (3) medium filtering (Fig. 3), (4) lung of discerning in each segmentation body with planar fit method splits (Fig. 4), and (5) lung splits classification (Fig. 5).Filtration is to have utilized higher relatively lung to split contrast, and plane fitting has utilized lung to split the characteristics of density higher relatively in the voxel space.
In this method, the CT image is counted as one group of three dimensional point cloud, and lung areas is subdivided into many little spheroids, is called the segmentation body, and the identification that lung splits is by seek corresponding plane by plane fitting in each segmentation spheroid.
The present invention has a lot of characteristics: at first, this method is very general, both can discern the straight line in the 2D space, also can discern the plane in the 3d space.
The second, aspect stable, this method is for noise or exceptional value and insensitive.
The 3rd, on counting yield, whole lung splits identifying only needs about 8 minutes, and the method [1] of the general invention of mentioning before such as the people of grade probably needs 20 minutes, and needs about 90 minutes by the method [17] of people such as van Rikxoort invention.
The 4th, on accuracy, neoteric method has higher accuracy.
The 5th, the implementation procedure of this method is very succinct, as long as planar fit method is directly applied to segmentation body in lung areas.In addition, dissimilar lungs split can finish simultaneously by a simple clustering method when being identified in lung and splitting testing process.As far as we know, now also there is not any method can carry out identification and the classification that lung splits simultaneously.
Description of drawings
Concrete embodiment will be described in detail by following accompanying drawing.
Fig. 1 to Fig. 5 is the basic step that lung splits dividing method:
Fig. 1 lung areas is cut apart;
The segmentation of Fig. 2 lung areas;
Fig. 3 medium filtering;
Fig. 4 splits with plane fitting analytic approach identification lung;
The dissimilar lungs of Fig. 5 split classification.
Fig. 6 to Figure 11 is by using the performance of two and three dimensions example illustrated planar approximating method:
Fig. 6 and Fig. 7 have shown the match of the straight line that is assembled by point in the two-dimension picture.The result that the green line representative is drawn by least squares approach, and red line has been represented the result who is drawn by the density approximating method of describing among the present invention;
Fig. 8 shows the three-dimensional point cloud atlas of the spheroid that a pair is segmented out from chest CT;
The point cloud of Fig. 9 among visual Fig. 8 of inner product (formula (6) just) of vector sum normal vector;
Figure 10 splits testing result (representing with redness) by using the resulting lung of the present invention;
Figure 11 shows the function of three-dimensional point cloud in the corresponding diagram 8
Figure 12 to Figure 14 is clear, and the lung of having represented splits identification and sorting result:
Figure 12 one secondary CT image;
Figure 13 is through identification, and the lung of classification splits (with different color showings);
Figure 14 is that the final lung that obtains after having removed among Figure 13 the non-lung split plot represented by arrow splits that image is cut apart and classification chart.
Figure 15 to Figure 17 shows that by example institute's inventive method is in the performance of identification lung aspect splitting.
Figure 18 to Figure 20 is the example of splitting at the lung that is recognized in the case that serious bronchiectasis (cyst cystic fibrosis) arranged.
Figure 21 to Figure 23 is that this inventive method identification lung splits example, and the lung in this example splits on image very fuzzy.
Embodiment
Planar fit method is this neoteric key component, below we with clear description.
The first step is the selection of objective function.
In Euclidean three dimensions, set N point
, our target is to find a plane
, come to minimize on the energy following objective function:
Wherein,
Be a little
p i To the plane
FDistance,
(equaling 1 in this research) is specified point
p i Weight (equal 1, just will somewhat same weighting in this research).From formula 1, can release distance function
In minimizing, play a crucial role.It should be noted that in this research, with the point that disperses
p i Be positioned at (in the segmentation body) in the predefined spheroid.The restriction of this employing encirclement ball can not influence the versatility of algorithm, because for one group of known point in the space, has one ' encirclement ball ' forever.In Euclidean three dimensions, the degree of freedom on plane is 3, and this plane can be decided by three parameters.The plane
FMay be defined as
Wherein,
It is the definition plane
FThree parameters.Though just can directly solve formula (1) with least squares approach, wherein
Although, the efficient height, the result who is solved by least squares approach but is easy to be subjected to the interference of noise or exceptional value (as Fig. 6-Fig. 7).For this reason, we select to use the distance function by formula (3) definition
D, it has relatively accurate solution, and to exceptional value insensitive [18-19]:
Wherein,
uIt is scalar.In this research,
uSplit thickness with lung and be closely connected, will
uBe made as described lung and split the only about half of value of thickness (for example from 1.0 to 1.5 millimeters).Here, we can be formulated as formula (1) again
Bring formula (1) and formula (3) into formula (4), by finding one group of parameters optimization (just
) after, the plane fitting process is from energy minimization
(formula (1) just) is converted into the energy maximization
(formula (5) just).
(5)。
From formula (5), as can be seen, has only the point of working as to the plane
FDistance less than
uThe time, these points can have influence on objective function
E'Since in the space density of point be used for discerning lung in the segmentation body split the plane (if any), doing, we advise using the function of formula (5)
E' Divided by plane F
Cross sectional area S.Like this, the problem that will solve based on the plane fitting of density is to find parameters optimization by maximization formula (6)
,
Wherein ". " represents inner product,
Be with respect to spheroid midplane F
Cross sectional area,
oWith
rCenter and the radius of representing spheroid respectively.With
Divided by S is in order to ensure objective function
Be independent of outside the S, so objective function
Can only be to determine by near the dot density the plane.
Because target function type (6) is not a convex function, therefore can not come parameters optimization with the gradient descent method.Though exhaustive search also is feasible, computation complexity is quite high.Such as for grid system
, temporal computation complexity is
, in practical operation, be difficult to be accepted.Therefore, need an effective optimization algorithm.
Next be to come parameters optimization with the method for integration
ρ
(7)。
If variable
mWith a vector of unit length
n, when
,
, it is in the spheroid of radius that point exists with r
(9)。
Wherein,
=2r/
,
It is the compartment distance of discrete ρ.Here, we will
Be made as a little value (for example, 0.1 millimeter), A, B, C are a in the formula (10), b, the integration of c
Wherein,
,
Here
ρ – mWith
uBe
Integral multiple; Therefore,
lWith
kIt is integer.With the method for integration, the parameter that we can be optimized
ρ,And can with computation complexity from
Reduce to
Though objective function
Be not with
Relevant convex function, but right at optimized parameter
Near zone (
) functional value with (describing among Figure 11) zone functional value big equally.Utilize this characteristics, it is right that we propose to search for optimized parameter effectively with refinement strategy progressively
We at first use a relatively large error
(search spacing just,
Arrive
)
Near-optimization parameter in the search type in the space (12) is right
Then, use the less error of predefine again
(for example
) by
It is right to search for the final optimization pass parameter in the segmentation space of definition
Consider the characteristics that lung splits, for fear of limited in local minimum, for initial ranging
=0.1 is just enough.Such as, at Fig. 8-11, the sphere segmentation body that the some cloud that takes out from chest CT is checked constitutes,
Figure 11 shows the three-dimensional point cloud objective function in the corresponding diagram 8
, wherein,
Utilize progressively refinement strategy, search
Time from
Reduce to
Such as,
With
,
The actual search time by original general 10
5Drop to 10
3
In order to raise the efficiency and guarantee desired accuracy, homogenize on the half unit spherome surface (just unified sampling) correspondence
Candidate's normal vector
Otherwise, when
θValue hour (for example,
θ=0, even establish
Be different values,
Also be identical.), will there be many similarity direction search to be repeated.In the present invention, for each self-generating
To with normal vector n, specify θ=0:
: π and α=0:
/ sin (θ): π, the search parameters optimization is right
Time by
Drop to
Fig. 6-Figure 11 has verified the performance of fitting algorithm by the example of two and three dimensions.Although there are many exceptional values to disturb, line in the picture and surface (plane) come out by accurate recognition.Fig. 6 and Fig. 7 have shown the fitting a straight line of two-dimension picture mid point set.The result that the green line representative is drawn by least squares approach, and red line has been represented the result who is drawn by the density approximating method of being invented.Fig. 8 shows the three-dimensional point cloud atlas of the spheroid that a pair extracts from chest CT.The point cloud of Fig. 9 among visual Fig. 8 of inner product (formula (6) just) of vector sum normal vector.Figure 10 has shown fit Plane (representing with redness) by the method for using the present invention.Figure 11 shows the function of three-dimensional point cloud in the corresponding diagram 8
Lung splits the basic step of cutting apart: extremely shown in Figure 5 as Fig. 1.
Lung segmentation: it is to dwindle the search volume that lung splits identification that lung areas is cut apart (as Fig. 2) fundamental purpose, thereby improves counting yield, has eliminated and produce the possibility that lung splits identification error (false positive) beyond lung areas.
Lung areas segmentation: the lung areas of cutting apart is further segmented.At first, be that lung areas among each CT is created one and segmented grid system with axle alignment bounding box AABB, each cell size is 5
*5
*5 cubic millimeters.The center of segmentation spheroid is located at the summit of case system, is radius (as Fig. 3) with r.Because lung splits and generally has relatively low curvature, we are located at 10 millimeters with the r value.Why we be because the rotational invariance of spheroid and can efficient calculation spheroid cross sectional area in any direction with spheroid rather than cube.If compartment is divided into 5mm, in Euclidean three dimensions, any point
GTo nearest summit
ODistance
(d)Less than
The millimeter or 4.4 millimeters.
Point filters: in order to dwindle the search volume that lung splits, we are with two filter process lung areas.At first, because the brightness value that splits of lung (I (
p)) generally be from-900HU(Hausfeld unit) to-300HU, like this, the voxel outside this range of luminance values is disallowable as non-crack voxel.The second, because territory, most of lung split plot has higher brightness than circle zone (voxel) on every side, therefore some can be proposed in the regional area less than the voxel of intermediate intensity value.Regional area size definition interested is the spheroid size in the lung areas segmentation.
Lung splits identification: split by plane fitting identification lung in each segmentation body.Suppose one with
o(
x 0 , y 0 , z 0) be center, radius
rBe 10 millimeters spheroid, through the point set to be selected of remainder in spheroid after the filter operation
The objective function of application formula (6), come match by
Optimal planar after the parametrization
FSince represent voxel intensity value that lung splits in different detections with zones of different in identical detection in all very inequality, the weighted value of points all in the formula (6) is made as identical value (just
), just only consider the plane fitting Density Distribution.For the planar structure (being that lung splits) in accurate each segmentation body of identification, only consider spheroid (
o) distance at center is less than 5 millimeters plane
, just
Mm has wherein defined in formula (2)
f(o).Guaranteed spheroid like this
In
FCross sectional area greater than
Mm
2Any point
GThe distance that arrives nearest lattice summit so just has from the ball centre to the point-of-interest less than 5 millimeters
(G)Less than 5 millimeters segmentation bodies (center is located at place, nearest summit).Therefore, considered in the fit procedure that lung splits have a few.Suppose that the optimal planar in the spheroid that identifies is
, simultaneously
With
Show normal vector and (circle) transversal section center relevant respectively, if energy function with spherical volume
Greater than predetermined threshold
T, the circular cross sections that identifies (plane) is just regarded the part that (supposition is) lung splits as.In this invention, the application self-adapting threshold value
T(for example,
).Density basis fitting algorithm is used for all segmentation bodies will be identified lung and split (as Figure 12-shown in Figure 14).Figure 12-Figure 14 has represented that lung splits identification and classification results: Figure 12 one secondary CT image; Figure 13 splits (covering with different color) through the lung of identification, classification, and Figure 14 is that the final lung that obtains after having removed among Figure 13 the non-lung split plot represented by arrow splits that image is cut apart and classification chart.
Lung splits classification: except the identification lung splits, this new invention can be split lung simultaneously and classify.Suppose and in two adjacent spheroids, identify two lung broken face sheets
F i With
F j , be met as one of following condition, just suppose it two be belong to identical gathering (just
),
Suppose in adjacent spheroid and identify
F i With
F j Flat blocks, parameter
Be used for controlling
F i With
F j Between curvature or allow the angle, parameter
Be used for controlling
F i With
F j Between the continuity of distance.Here, by rule of thumb will
Be made as 0.98 millimeter, similarly will
Be made as 2 millimeters.Split the series arrangement of successively decreasing of quantity between accumulation area by the lung that is comprised in each gathering.In the lung of a left side, first (maximum) gathering is used for representing left oblique segmentation.In the right lung, first and second (maximum) assemble be used for respectively representing right oblique segmentation and laterally lung split (as Figure 12-Figure 14).Other all gathering is removed as non-lung split plot.
Test and result: the performance of inventive method in order to assess, we have collected 200 routine chronic obstructive pulmonary disease (COPD) chest CT data, and the lung that detects in these data with this method splits, and the result is compared on same group of data with a kind of method before.Recognition performance is measured with Hausdorff distance and tired knowledge error distance distribution (CEDD).The lung that identifies with these two kinds of methods splits nearly 95% within 2 millimeters distances, between average error probably at 0.6 millimeter, average crack maximum error is 16 millimeters.Computing time the aspect, this new method (being configured to the Duo i3 of Intel, 3.2 hertz of central processing units, the desktop computer of 6.0G internal memory) on Daepori energising brain was split identification general 8 minutes consuming time to the lung on the CT image.
In order to show that visually the new invention method splits performance in the identification at lung, we provide the example as Figure 15-Figure 17.In the example, the crack of identifying black line sign.
The example that Figure 18 to Figure 20 provides has shown that lung splits identification in serious bronchiectasis the is arranged abnormality detection of (cyst cystic fibrosis).
Figure 21 to Figure 23 is example that detection difficulty is bigger in splitting with new invention method identification lung.It is very difficult wanting the clear definite position of seeing that clearing lung-heat is split.And method [1] before is applied in this detection, be complete failure.
Here, described a full automatic lung and split identification and sorting technique, by utilize that lung splits in the image space the characteristics of Density Distribution.This method is the Free Surface of complexity to be discerned be converted into serial plane fitting, has characteristics such as versatility, high efficiency in method.
Though by being described with reference to its concrete enforcement, those of ordinary skill in the art will be appreciated that, can realize a lot of modifications, reinforcement and/or variation under the situation that does not break away from the spirit or scope of the present invention in the present invention.The present invention is intended to be limited by the scope of its claims and equivalent thereof.
List of references
[1] J. Pu, J. K. Leader, B. Zheng, F. Knollmann, C. Fuhrman, F. C. Sciurba, and D. Gur, “
A computational geometry approach to automated pulmonary fissure segmentation in ct examinations,” IEEE Transactions on Medical Imaging, vol. 28, no. 5, pp. 710–719, 2009.
[2] R. Wiemker, T. Blow, and T. Blaffert,
“Unsupervised extraction of the pulmonary interlobar fissures from high resolution thoracic ct data,” in Proc. 19th Int. Congress ExhibitionComputer Assist. Radiol. Surg. (CARS), Berlin, Germany, 2005, pp. 1121–1126.
[3] E. M. van Rikxoort, B. van Ginneken, M. Klik, and M. Prokop,
“Supervised enhancement filters: Application to fissure detection in chest ct scans,” IEEE Transactions on Medical Imaging, vol. 27, no. 1, pp. 1–10, Jan. 2008.
[4] E. M. van Rikxoort, B. de Hoop, S. van de Vorst, M. Prokop, and B. van Ginneken,
“Automatic segmentation of pulmonary segments from volumetric chest ct scans,” IEEE Transactions on Medical Imaging, vol. 28, no. 4, pp. 621–630, 2009.
[5] S. Saita, M. Kubo, Y. Kawata, N. Niki, H. Ohmatsu, and N. Moriyama,
“An algorithm for the extraction of pulmonary fissures from low-dose multislice ct image,” System and Computers in Japan, vol. 37, no. 9, pp. 63–76, 2006.
[6] J.-M. Kuhnigk, H. K. Hahn, M. Hindennach, V. Dicken, S. Krass, and H.-O. Peitgen,
“Lung lobe segmentation by anatomy-guided 3d watershed transform,” in SPIE, 2003, pp. 1482–1490.
[7] S. Ukil and J. M. Reinhardt,
“Anatomy-guided lung lobe segmentation in x-ray ct images,” IEEE Transactions on Medical Imaging, vol. 28, no. 2, pp. 202–214, 2009.
[8] X. Zhou, T. Hayashi, T. Hara, H. Fujita, and R. Yokoyama,
“Automatic recognition of lung lobes and fissures from multi-slice ct images,” in SPIE, 2004, pp. 1629–1633.
[9] J. Wang, M. Betke1, and J. P. Ko,
“Shape-based curve growing model and adaptive regularization for pulmonary fissure segmentation in ct,” in MICCAI, LNCS 3216, 2004, pp. 541–548.
[10] J. Wang, M. Betke, and J. Ko,
“Pulmonary fissure segmentation on CT,” Medical Image Analysis, vol. 10, pp. 530–547, 2006.
[11] M. Kubo, N. Niki, S. Nakagawa, K. Eguchi, M. Kaneko, N. Moriyama, H. Omatsu, R. K. R, and N. Yamaguchi,
“Extraction algorithm of pulmonary fissures from thin-section ct images based on linear feature detector method,” IEEE Transactions on Nuclear Science, vol. 46, no. 6, pp. 2128–2133, 1999.
[12] M. Kubo, Y. Kawata, N. Niki, K. Eguchi, H. Omatsu, R. K. R, M. Kaneko, M. Kusumoto, N. Moriyama, K. Mori, and H. Nishiyama,
“Automatic extraction of pulmonary fissures from multidetector-row ct images,” in Proceedings of the IEEE International Conference on Image Processing (ICIP01), 2001, pp. 1091–1094.
[13] K. Li, X. Wu, D. Chen, and M. Sonka,
“Optimal surface segmentation in volumetric images-a graph-theoretic approach,” IEEE transactions on Pattern Analysis and Machine Intelligence, Vol. 28, no. 1, pp. 119–134, 2006.
[14] L. Zhang and J. Reinhardt,
“Detection of pulmonary lobar fissures using fuzzy logic,” in Proc. of SPIE 1999, pp. 188–199.
[15] L. Zhang, E. Hoffman, and J. Reinhardt,
“Pulmonary lobe segmentation by graph search with 3-d shape constraints,” in SPIE, 2001, pp. 204–215.
[16] L. Zhang, E. A. Hoffman, and J. M. Reinhardt,
“Atlas-driven lung lobe segmentation in volumetric x-ray ct images,” IEEE Transactions on Medical Imaging, vol. 25, no. 1, pp. 1–16, 2006.
[17] E. van Rikxoort, M. Prokop, B. de Hoop, M. Viergever, J. Pluim, and B. van Ginneken,
“Automatic segmentation of pulmonary lobes robust against incomplete fissures,” IEEE Transactions on Medical Imaging, vol. 29, no. 6, pp. 1286–1296, 2010.
[18] P. Torr and D. Murray,
“The development and comparison of robust methods for estimating the fundamental matrix,” International Journal of Computer Vision, vol. 24, no. 3, pp. 271–300, 1997.
[19] S. Choi, T. Kim, and W. Yu,
“Performance evaluation of RANSAC family,” in Proc. of the 20th British Machine Vision Conference, 2009; pp. 1-12.
[20] J. Pu, D. S. Paik, X. Meng, J. E. Roos, and G. D. Rubin,
“Shape break-and-repair strategy and its application to computerized medical image segmentation,”IEEE Transactions on Visualization & Computer Graphics, vol. 17, no. 1, pp. 115-124, 2011.
Claims (9)
2. in accordance with the method for claim 1, wherein the selection of objective function comprises:
3. parameters optimization wherein in accordance with the method for claim 1,
ρ,Also comprise: the method with integration is come parameters optimization
ρ
5. also comprise in accordance with the method for claim 4:
With a relatively large error
(search spacing
Arrive
)
Search near-optimal parameter is right in the space
6. the divided method of lung areas comprises and uses axle alignment bounding box AABB to build a segmentation system for each CT detects.
7. with the method for two filter process through the lung areas of over-segmentation.
8. in segmentation body to be selected, discern the planar fit method that lung splits.
9. split in the testing process at lung and split clustering method in order to discern different types of lung simultaneously.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011101907972A CN102254097A (en) | 2011-07-08 | 2011-07-08 | Method for identifying fissure on lung CT (computed tomography) image |
PCT/CN2011/079226 WO2013007063A1 (en) | 2011-07-08 | 2011-09-01 | Method for identifying lung fissure on ct image of lungs |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011101907972A CN102254097A (en) | 2011-07-08 | 2011-07-08 | Method for identifying fissure on lung CT (computed tomography) image |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102254097A true CN102254097A (en) | 2011-11-23 |
Family
ID=44981358
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011101907972A Pending CN102254097A (en) | 2011-07-08 | 2011-07-08 | Method for identifying fissure on lung CT (computed tomography) image |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN102254097A (en) |
WO (1) | WO2013007063A1 (en) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013078782A1 (en) * | 2011-12-02 | 2013-06-06 | Pu Jiantao | The airway wall recognition method |
CN106373118A (en) * | 2016-08-30 | 2017-02-01 | 华中科技大学 | A complex curved surface part point cloud reduction method capable of effectively keeping boundary and local features |
CN106875379A (en) * | 2017-01-10 | 2017-06-20 | 陕西渭南神州德信医学成像技术有限公司 | Lung splits integrity degree appraisal procedure, device and system |
CN104574319B (en) * | 2015-01-22 | 2017-10-13 | 深圳大学 | The blood vessel Enhancement Method and system of a kind of lung CT image |
CN107392910A (en) * | 2017-07-06 | 2017-11-24 | 沈阳东软医疗系统有限公司 | A kind of lobe of the lung dividing method and device based on CT images |
CN107909581A (en) * | 2017-11-03 | 2018-04-13 | 杭州依图医疗技术有限公司 | Lobe of the lung section dividing method, device, system, storage medium and the equipment of CT images |
CN109300122A (en) * | 2018-09-14 | 2019-02-01 | 沈阳东软医疗系统有限公司 | Image procossing and Threshold, device and equipment |
CN109658425A (en) * | 2018-12-12 | 2019-04-19 | 上海联影医疗科技有限公司 | A kind of lobe of the lung dividing method, device, computer equipment and storage medium |
CN111449669A (en) * | 2014-07-02 | 2020-07-28 | 柯惠有限合伙公司 | System and method for detecting trachea |
CN111986299A (en) * | 2019-05-24 | 2020-11-24 | 北京京东尚科信息技术有限公司 | Point cloud data processing method, device, equipment and storage medium |
US10949977B2 (en) | 2017-06-30 | 2021-03-16 | Shanghai United Imaging Healthcare Co., Ltd. | Methods and systems for image segmentation |
CN113160186A (en) * | 2021-04-27 | 2021-07-23 | 青岛海信医疗设备股份有限公司 | Lung lobe segmentation method and related device |
CN114092470A (en) * | 2021-12-08 | 2022-02-25 | 之江实验室 | Deep learning-based automatic detection method and device for pulmonary fissure |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6634186B2 (en) | 2016-08-02 | 2020-01-22 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | Image processing system and image processing method |
EP3968271A1 (en) | 2020-09-11 | 2022-03-16 | Bayer AG | Analysis of intrapulmonary branches |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1391005A (en) * | 2001-12-30 | 2003-01-15 | 吴国雄 | Automatic design simulating system for plane figure of highway in mountain area |
CN1598868A (en) * | 2004-09-06 | 2005-03-23 | 南京大学 | Changing method abstrated characteristic in mode recognition |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008142481A (en) * | 2006-12-13 | 2008-06-26 | Med Solution Kk | Apparatus and program for carrying out segmentation of lungs to units of pulmonary segment automatically |
CN102460471B (en) * | 2009-06-19 | 2015-03-18 | 美国医软科技公司 | Systems for computer aided lung nodule detection in chest tomosynthesis imaging |
-
2011
- 2011-07-08 CN CN2011101907972A patent/CN102254097A/en active Pending
- 2011-09-01 WO PCT/CN2011/079226 patent/WO2013007063A1/en active Application Filing
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1391005A (en) * | 2001-12-30 | 2003-01-15 | 吴国雄 | Automatic design simulating system for plane figure of highway in mountain area |
CN1598868A (en) * | 2004-09-06 | 2005-03-23 | 南京大学 | Changing method abstrated characteristic in mode recognition |
Non-Patent Citations (2)
Title |
---|
胡金水等: "高分辨率CT图像肺部裂纹的自动检测", 《四川大学计算机学报》 * |
范立南: "基于CT影像的肺裂计算机辅助检测方法", 《仪器仪表用户》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013078782A1 (en) * | 2011-12-02 | 2013-06-06 | Pu Jiantao | The airway wall recognition method |
CN111449669B (en) * | 2014-07-02 | 2023-04-18 | 柯惠有限合伙公司 | System and method for detecting trachea |
CN111449669A (en) * | 2014-07-02 | 2020-07-28 | 柯惠有限合伙公司 | System and method for detecting trachea |
CN104574319B (en) * | 2015-01-22 | 2017-10-13 | 深圳大学 | The blood vessel Enhancement Method and system of a kind of lung CT image |
CN106373118A (en) * | 2016-08-30 | 2017-02-01 | 华中科技大学 | A complex curved surface part point cloud reduction method capable of effectively keeping boundary and local features |
CN106875379A (en) * | 2017-01-10 | 2017-06-20 | 陕西渭南神州德信医学成像技术有限公司 | Lung splits integrity degree appraisal procedure, device and system |
US11710242B2 (en) | 2017-06-30 | 2023-07-25 | Shanghai United Imaging Healthcare Co., Ltd. | Methods and systems for image segmentation |
US10949977B2 (en) | 2017-06-30 | 2021-03-16 | Shanghai United Imaging Healthcare Co., Ltd. | Methods and systems for image segmentation |
CN107392910B (en) * | 2017-07-06 | 2020-01-07 | 东软医疗系统股份有限公司 | Lung lobe segmentation method and device based on CT image |
CN107392910A (en) * | 2017-07-06 | 2017-11-24 | 沈阳东软医疗系统有限公司 | A kind of lobe of the lung dividing method and device based on CT images |
CN107909581A (en) * | 2017-11-03 | 2018-04-13 | 杭州依图医疗技术有限公司 | Lobe of the lung section dividing method, device, system, storage medium and the equipment of CT images |
CN109636811B (en) * | 2017-11-03 | 2020-06-12 | 杭州依图医疗技术有限公司 | Integration method and device for lung lobe segment segmentation of CT (computed tomography) image |
CN109636811A (en) * | 2017-11-03 | 2019-04-16 | 杭州依图医疗技术有限公司 | Integration method, the device of the lobe of the lung section segmentation of CT images |
CN107909581B (en) * | 2017-11-03 | 2019-01-29 | 杭州依图医疗技术有限公司 | Lobe of the lung section dividing method, device, system, storage medium and the equipment of CT images |
CN109300122B (en) * | 2018-09-14 | 2023-04-07 | 东软医疗系统股份有限公司 | Image processing and threshold determining method, device and equipment |
CN109300122A (en) * | 2018-09-14 | 2019-02-01 | 沈阳东软医疗系统有限公司 | Image procossing and Threshold, device and equipment |
CN109658425A (en) * | 2018-12-12 | 2019-04-19 | 上海联影医疗科技有限公司 | A kind of lobe of the lung dividing method, device, computer equipment and storage medium |
CN109658425B (en) * | 2018-12-12 | 2021-12-28 | 上海联影医疗科技股份有限公司 | Lung lobe segmentation method and device, computer equipment and storage medium |
CN111986299A (en) * | 2019-05-24 | 2020-11-24 | 北京京东尚科信息技术有限公司 | Point cloud data processing method, device, equipment and storage medium |
CN111986299B (en) * | 2019-05-24 | 2024-03-01 | 北京京东乾石科技有限公司 | Point cloud data processing method, device, equipment and storage medium |
CN113160186A (en) * | 2021-04-27 | 2021-07-23 | 青岛海信医疗设备股份有限公司 | Lung lobe segmentation method and related device |
CN114092470A (en) * | 2021-12-08 | 2022-02-25 | 之江实验室 | Deep learning-based automatic detection method and device for pulmonary fissure |
Also Published As
Publication number | Publication date |
---|---|
WO2013007063A1 (en) | 2013-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102254097A (en) | Method for identifying fissure on lung CT (computed tomography) image | |
El-Regaily et al. | Survey of computer aided detection systems for lung cancer in computed tomography | |
El-Baz et al. | Three-dimensional shape analysis using spherical harmonics for early assessment of detected lung nodules | |
Basavanhally et al. | Incorporating domain knowledge for tubule detection in breast histopathology using O'Callaghan neighborhoods | |
Ge et al. | Computer‐aided detection of lung nodules: false positive reduction using a 3D gradient field method and 3D ellipsoid fitting | |
Seff et al. | 2D view aggregation for lymph node detection using a shallow hierarchy of linear classifiers | |
US8126229B2 (en) | Reduction of lymph tissue false positives in pulmonary embolism detection | |
US20120070074A1 (en) | Method and System for Training a Landmark Detector using Multiple Instance Learning | |
EP1588316A4 (en) | A method, system, and computer program product for computer-aided detection of nodules with three dimensional shape enhancement filters | |
CN101449291A (en) | Process and system for automatically recognising preneoplastic abnormalities in anatomical structures, and corresponding computer program | |
Gonçalves et al. | Learning lung nodule malignancy likelihood from radiologist annotations or diagnosis data | |
Moreno et al. | Study of medical image processing techniques applied to lung cancer | |
Bermejo-Peláez et al. | Emphysema classification using a multi-view convolutional network | |
Duggan et al. | A technique for lung nodule candidate detection in CT using global minimization methods | |
CN113902978B (en) | Depth learning-based interpretable SAR image target detection method and system | |
Saha Tchinda et al. | Automated extraction of the intestinal parasite in the microscopic images using active contours and the Hough transform | |
Devaki et al. | Study of computed tomography images of the lungs: A survey | |
Gu et al. | Identification of pulmonary fissures using a piecewise plane fitting algorithm | |
Krishnaveni et al. | Image Based Group Classifier for Brain Tumor Detection Using Machine Learning Technique. | |
Kayalvizhi et al. | Analysis of ventricle regions in Alzheimer’s brain MR images using level set based methods | |
Zhao et al. | Leveraging tree statistics for extracting anatomical trees from 3D medical images | |
Shaziya et al. | Comprehensive review of automatic lung segmentation techniques on pulmonary CT images | |
Hibbard et al. | Automated recognition and mapping of immunolabelled neurons in the developing brain | |
Ge et al. | Computer-aided detection of lung nodules: false positive reduction using a 3D gradient field method | |
Li et al. | A novel surface-based geometric approach for 3d dendritic spine detection from multi-photon excitation microscopy images |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C53 | Correction of patent of invention or patent application | ||
CB03 | Change of inventor or designer information |
Inventor after: Pu Jiantao Inventor after: Meng Xin Inventor before: Pu Jiantao |
|
COR | Change of bibliographic data |
Free format text: CORRECT: INVENTOR; FROM: PU JIANTAO TO: PU JIANTAO MENG XIN |
|
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20111123 |