CN109166108A - A kind of automatic identifying method of CT images pulmonary abnormalities tissue - Google Patents
A kind of automatic identifying method of CT images pulmonary abnormalities tissue Download PDFInfo
- Publication number
- CN109166108A CN109166108A CN201810922391.0A CN201810922391A CN109166108A CN 109166108 A CN109166108 A CN 109166108A CN 201810922391 A CN201810922391 A CN 201810922391A CN 109166108 A CN109166108 A CN 109166108A
- Authority
- CN
- China
- Prior art keywords
- block
- template
- images
- lung
- tissue
- 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.)
- Granted
Links
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
- G06T7/0014—Biomedical image inspection using an image reference approach
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/194—Segmentation; Edge detection involving foreground-background segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/74—Image or video pattern matching; Proximity measures in feature spaces
- G06V10/75—Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
- G06V10/751—Comparing pixel values or logical combinations thereof, or feature values having positional relevance, e.g. template matching
-
- 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)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Evolutionary Computation (AREA)
- Computing Systems (AREA)
- Artificial Intelligence (AREA)
- Databases & Information Systems (AREA)
- Multimedia (AREA)
- Geometry (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The present invention relates in medically computer implemented automatic diagosis technology, a kind of automatic identifying method for the CT images pulmonary abnormalities tissue realized using template difference comparison technology.The parameter vector mainly extracted with the CT lung image of one group of arm's length standard is compared, the method for being quickly found out abnormal structure and automatic marking as template group with real-time dynamic CT lung image.
Description
Technical field
The present invention relates in medically computer implemented automatic diagosis technology, more specifically, being related to a kind of using mould
Plate difference comparison technology and the automatic identifying method of CT images pulmonary abnormalities tissue realized.
Background technique
The a large amount of extensive use of CT images in medicine, has had resulted in two problems: on the one hand, to each patient one
The secondary information content for checking image generated is increasing, and the image information that doctor's daily requirement is read is too many, and working strength is got over
Come bigger;On the other hand, diagnostic level is highly dependent on the experience and ability of diagnostician, and mistaken diagnosis happens occasionally.By medicine
The development of image format DICOM standard and computer-aided diagnosis technology realizes that the automation of the abnormal structure of CT images is known
Not, the accuracy rate of diagosis diagnosis is provided, it is more and more important.
It realizes that the main application target of the automatic identifying method of CT images pulmonary abnormalities tissue shows two aspects: 1. existing
On the basis of handling a large amount of patient CT lung image data, suspicious lesion is found out as precisely as possible, emphasis mark,
It is supplied to the reference that doctor's behaviours further diagnose, the workload that doctor reads piece is reduced, the accuracy rate of diagnosis is provided.2. big
Size medical image processing method and modeling is realized on the basis of largely analyzing certain special diseases, easy in iconography
The normal tissue and lesion information obscured are differentiated by the selection of model parameter, and real lesion is screened out to come.
Summary of the invention
Defect in view of the prior art, the present invention relates to adopt in medically computer implemented automatic diagosis technology, one kind
The automatic identifying method for the CT images pulmonary abnormalities tissue realized with template difference comparison technology.It is mainly normal with one group
The parameter vector that the CT lung image of standard extracts is compared, quickly as template group with real-time dynamic CT lung image
The method for finding abnormal structure and automatic marking.
We filter out the lung CT image of 50 Healthy Peoples according to gender, age bracket, fat rank respectively first
As standard form, amount to 50, template number, male numbers from M1, M2 ... ... M25, and women is from F1, F2 ... ... F25
Then number generates template data by step.
CT standard video template characteristic the present invention relates to building be by the two-dimensional polar coordinates of multiple groups, node type, curvature,
The data such as distance and contraction-expansion factor describe each minutiae feature between directionality, efficient association point, and these compare
Template saves in the database in strict accordance with packet numbering.
In actual CT lung images comparison process, we will determine age and the gender of patient first, then carry out with
The identical analysis processing step of template image, until finally obtaining its characteristic.Then corresponding grouping feature mould is found again
Plate data carry out the processing of comparison one by one respectively.
During template matching, three main processing steps are needed:
(1) we obtain and the difference block S with 5 lung image template datas of group respectivelyn(may also not have) and area grade
Other βn(dividing 0 ~ 10 grade by block area);
(2) in SnIn, overlapping region s ' all existing for more than 2 template above(n)Respectively by the validity of block closure
Carry out image removing.If image closure is invalid, from s '(n)It is deleted in list;
(3) each s ' is calculated(n)Central point, using maximum Edge Distance as radius r, finally in original CT images
Circle remembers effective abnormal block out.
Detailed description of the invention
Attached drawing 1: being that CT lung of the present invention template characteristic generates schematic diagram;
Attached drawing 2: being the schematic diagram of CT lung image comparison calculation process of the present invention.
Specific embodiment
The present invention will be further described with reference to the accompanying drawings and examples:
We filter out the lung CT image conduct of 50 Healthy Peoples according to gender, age bracket, fat rank respectively first
Standard form amounts to 50, and template number, male numbers from M1, M2 ... ... M25, and women is numbered from F1, F2 ... ... F25,
Then template data is generated by step, as shown in Figure 1.
CT lung template characteristic in attached drawing 1 generates, and is classified according to gender, by young (0 years old -6 years old), teenager
(7 years old-17 years old), young (18 years old-40 years old), middle aged (41-65 years old) and old (after 66 years old) grouping is established.Each template
There is the ID number of oneself, it is gender (1: male, 2: female) that ID, which be 3: first,;Second represents group number (1-5);Third position
It represents with the serial number in group.
CT standard video template characteristic minutiae point is by two-dimensional polar coordinates, node type, curvature, directionality, effectively pass
The data such as distance and contraction-expansion factor describe each minutiae feature between connection point, and these comparison templates in strict accordance with point
Group # saves in the database.
During template matching, firstly for the CT lung images of any input, it is intended to carry out and template processing step
It is identical to pre-process to obtain the characteristic of minutiae point.Then corresponding templates group is selected, is compared one by one.Finally respectively
Out with the difference block S of same group of 5 lung image template datasn(may also not have) and area rank βnIt (is divided by block area
0 ~ 10 grade), it then finds out comprising SnThe minimum rectangle of block, and four coordinate points are saved, for βnValue also to save.Allow
There are multiple difference blocks, Sn∈ Diff_Area, data structure are as follows:
struct Diff_Area {
Short Block_Number;// difference block serial number
Struct BLOCK_INDEX [Block_Number];// difference block coordinate
Short Scale[Block_Number];// difference block area rank
};
Difference block coordinate data structure is as follows:
Struct BLOCK_INDEX {
Int Top_X;X-coordinate on // block
Int Top_Y;Y coordinate on // block
Int Down_X;X-coordinate under // block
Int Down_Y;Y coordinate under // block
Int Left_X;The left X-coordinate of // block
Int Left_Y;The left Y coordinate of // block
Int Right_X;The right X-coordinate of // block
Int Right_Y;The right Y coordinate of // block
};
In this alignment algorithm, the second step of progress is that method is: respectively by five template matchings go out difference block between into
Capable operation of seeking common ground.In SnIn, overlapping region s ' all existing for more than 2 template above(n)It is accordingly to be regarded as effectively.
For effective block, the operation of image closed is carried out respectively.If superthreshold, it is invalid to be considered as closure, then from
s’(n)It is deleted in list.
Final step is to calculate each s '(n)Central point, using this maximum Edge Distance as radius r, finally in original
Circle remembers that effective abnormal block, alignment algorithm process are as shown in Figure 2 out in the CT images of beginning.
The foundation of CT lung image standard feature template is to need expert doctor to screen from existing CT file store to dispatch troops
Complete modeling work.For the CT lung image to be checked, is then automatically read by computer and complete the institutes such as to compare and mark
There is working link, without manually being intervened.Example below is then with CT lung image standard feature template and comparison process
Objective for implementation.
1) data structure of CT lung image standard feature template is established
Template data be by distance between two-dimensional polar coordinates, node type, curvature, directionality, efficient association point and it is flexible because
The data such as son describe each minutiae feature, and these comparison templates are stored in database in strict accordance with packet numbering
In.Therefore its key data structure is defined as follows
struct CT_ BellowsForm {
char ID[4];// template number
int GroupID;// packet numbering
char Sex;// gender: M: male, F: female
date Birthday;// the date of birth
int Block_ToT_Number;// block total quantity
int Block_Serial_No;// block serial number
struct Block_Area *PArea;// block pointer
struct Block_ Direct *PDir;// block directional diagram pointer
struct Block_ Curvature *PCur;// block directional curvature
struct Block_ KeyPoint *PKeyPoint;// key node coordinate
struct Block_KeyPoint_ Distance *PDist;// key point vector table
struct Block_ StretchFactor *PFact;// key point contraction-expansion factor
int *PScale;// block amount of area rank (1 ~ 10)
}
struct Block_Area {
Short Block_Number;// block serial number
Struct BLOCK_INDEX *PNo;// block coordinate
Short Scale[Block_Number];// block area rank
} *PArea;
Difference block coordinate data structure is as follows:
Struct BLOCK_INDEX {
Int Top_X;X-coordinate on // block
Int Top_Y;Y coordinate on // block
Int Down_X;X-coordinate under // block
Int Down_Y;Y coordinate under // block
Int Left_X;The left X-coordinate of // block
Int Left_Y;The left Y coordinate of // block
Int Right_X;The right X-coordinate of // block
Int Right_Y;The right Y coordinate of // block
} *PNo;
2) CT lung image pre-processes
L baseline drift processing
L image
L image segmentation
L block divides
l … ...
3) realization of alignment algorithm
We are it needs to be determined that refuse to recognize rate (FAR), misclassification rate (FRR) and similarity (Similarity) three matched default first
Amount is three groups:
0.1%、 0.02%、0.75%;
0.05%、0.06%、0.85%;
0.03%、0.08%、0.93%;
It is also required to simultaneously to every important indicator in image, block, minutiae point, sets up reasonable comparison threshold by largely calculating
Value, does not illustrate herein.
The ECGWaveVerify (WaveData, WaveTemplate, VerifyParameter) of alignment algorithm;
4) CT images template reconstructs
Template reconstruction of function name are as follows:
ECGWaveTemplateRefactoring(OldWaveTemplate, NewWaveTemplate ,Parameter);
5) CT images are template optimized
Template optimized function name are as follows:
ECGWaveTemplateOptimize (WaveTemplate, Parameter);
Effect of the invention is designed from the angle of CT lung image template, can be improved CT images template to a certain extent
Precision and reliability promote the accuracy rate diagnosed in advance automatically.
Claims (4)
1. a kind of automatic identifying method for the CT images pulmonary abnormalities tissue realized using template difference comparison technology, feature
It is, the CT lung image standard for establishing one group of arm's length standard form extracts its vector as template characteristic data.
2. the method according to claim 1, wherein template group and the comparison of the CT lung image acquired in real time are
Using difference block cutting techniques, the method for finding out abnormal block.
3. the method according to claim 1, wherein the spy that the vector with CT lung images size is adaptive
Sign.
4. the method according to claim 1, wherein the two-dimensional pole of the multiple groups obtained after CT images are pre-processed is sat
Mark, node type, curvature, directional diagram, distance and contraction-expansion factor between efficient association point, as lung's surrounding tissue template characteristic
Parameter vector.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810922391.0A CN109166108B (en) | 2018-08-14 | 2018-08-14 | Automatic identification method for abnormal lung tissue of CT (computed tomography) image |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810922391.0A CN109166108B (en) | 2018-08-14 | 2018-08-14 | Automatic identification method for abnormal lung tissue of CT (computed tomography) image |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109166108A true CN109166108A (en) | 2019-01-08 |
CN109166108B CN109166108B (en) | 2022-04-08 |
Family
ID=64895564
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810922391.0A Active CN109166108B (en) | 2018-08-14 | 2018-08-14 | Automatic identification method for abnormal lung tissue of CT (computed tomography) image |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109166108B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113450899A (en) * | 2021-06-22 | 2021-09-28 | 上海市第一人民医院 | Intelligent diagnosis guiding method based on artificial intelligence cardiopulmonary examination images |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006114003A1 (en) * | 2005-04-27 | 2006-11-02 | The Governors Of The University Of Alberta | A method and system for automatic detection and segmentation of tumors and associated edema (swelling) in magnetic resonance (mri) images |
CN1934589A (en) * | 2004-03-23 | 2007-03-21 | 美国西门子医疗解决公司 | Systems and methods providing automated decision support for medical imaging |
WO2008107809A2 (en) * | 2007-03-06 | 2008-09-12 | Koninklijke Philips Electronics, N.V. | Automated diagnosis and alignment supplemented with pet/mr flow estimation |
CN101727666A (en) * | 2008-11-03 | 2010-06-09 | 深圳迈瑞生物医疗电子股份有限公司 | Image segmentation method and device, and method for judging image inversion and distinguishing front side and back side of sternum |
CN103479375A (en) * | 2012-06-11 | 2014-01-01 | 柯尼卡美能达株式会社 | Thoracic diagnosis assistance information generation method, thoracic diagnosis assistance system, and dynamic state image processing apparatus |
CN106846317A (en) * | 2017-02-27 | 2017-06-13 | 北京连心医疗科技有限公司 | A kind of feature based extracts the method for retrieving medicine image with Similarity matching |
-
2018
- 2018-08-14 CN CN201810922391.0A patent/CN109166108B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1934589A (en) * | 2004-03-23 | 2007-03-21 | 美国西门子医疗解决公司 | Systems and methods providing automated decision support for medical imaging |
WO2006114003A1 (en) * | 2005-04-27 | 2006-11-02 | The Governors Of The University Of Alberta | A method and system for automatic detection and segmentation of tumors and associated edema (swelling) in magnetic resonance (mri) images |
WO2008107809A2 (en) * | 2007-03-06 | 2008-09-12 | Koninklijke Philips Electronics, N.V. | Automated diagnosis and alignment supplemented with pet/mr flow estimation |
CN101727666A (en) * | 2008-11-03 | 2010-06-09 | 深圳迈瑞生物医疗电子股份有限公司 | Image segmentation method and device, and method for judging image inversion and distinguishing front side and back side of sternum |
CN103479375A (en) * | 2012-06-11 | 2014-01-01 | 柯尼卡美能达株式会社 | Thoracic diagnosis assistance information generation method, thoracic diagnosis assistance system, and dynamic state image processing apparatus |
CN106846317A (en) * | 2017-02-27 | 2017-06-13 | 北京连心医疗科技有限公司 | A kind of feature based extracts the method for retrieving medicine image with Similarity matching |
Non-Patent Citations (2)
Title |
---|
MARINA SERGEJEVA等: ""Anatomical landmarks for registration of experimental image data to volumetric rodent brain atlasing templates"", 《JOURNAL OF NEUROSCIENCE METHODS》 * |
米子川: ""基于图形分析方法的函数型数据异常值检验实证研究"", 《统计与信息论坛》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113450899A (en) * | 2021-06-22 | 2021-09-28 | 上海市第一人民医院 | Intelligent diagnosis guiding method based on artificial intelligence cardiopulmonary examination images |
Also Published As
Publication number | Publication date |
---|---|
CN109166108B (en) | 2022-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Vijayarajeswari et al. | Classification of mammogram for early detection of breast cancer using SVM classifier and Hough transform | |
CN106056595B (en) | Based on the pernicious assistant diagnosis system of depth convolutional neural networks automatic identification Benign Thyroid Nodules | |
CN109389584A (en) | Multiple dimensioned rhinopharyngeal neoplasm dividing method based on CNN | |
US7783094B2 (en) | System and method of computer-aided detection | |
CN102395320B (en) | Medical apparatus and method for controlling the medical apparatus | |
CN108257135A (en) | The assistant diagnosis system of medical image features is understood based on deep learning method | |
US11651850B2 (en) | Computer vision technologies for rapid detection | |
CN109886933A (en) | A kind of medical image recognition method, apparatus and storage medium | |
JP2008259682A (en) | Section recognition result correcting device, method and program | |
JP7082993B2 (en) | Medical image processing equipment, methods and programs, diagnostic support equipment, methods and programs, and medical support systems and methods | |
CN109859168A (en) | A kind of X-ray rabat picture quality determines method and device | |
Vezzetti et al. | 3D Landmarking in multiexpression face analysis: a preliminary study on eyebrows and mouth | |
CN101103924A (en) | Galactophore cancer computer auxiliary diagnosis method based on galactophore X-ray radiography and system thereof | |
Sapate et al. | Breast cancer diagnosis using abnormalities on ipsilateral views of digital mammograms | |
CN109919254B (en) | Breast density classification method, system, readable storage medium and computer device | |
JP2019008349A (en) | Learning data generation support apparatus and learning data generation support method and learning data generation support program | |
JP2010000133A (en) | Image display, image display method and program | |
CN113808125A (en) | Medical image processing method, focus type identification method and related product | |
JP5539478B2 (en) | Information processing apparatus and information processing method | |
CN115862819A (en) | Medical image management method based on image processing | |
JP5640280B2 (en) | Osteoporosis diagnosis support device and osteoporosis diagnosis support program | |
CN109166108A (en) | A kind of automatic identifying method of CT images pulmonary abnormalities tissue | |
CN113539476A (en) | Stomach endoscopic biopsy Raman image auxiliary diagnosis method and system based on artificial intelligence | |
Ull et al. | Regression to Classification: Ordinal prediction of calcified vessels using customized ResNet50 | |
KR102480134B1 (en) | Apparatus for quality managment of medical image interpretation usnig machine learning, and method thereof |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |