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 PDF

Info

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
Application number
CN201810922391.0A
Other languages
Chinese (zh)
Other versions
CN109166108B (en
Inventor
张新财
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Rongda Information Technology Co Ltd
Original Assignee
Shanghai Rongda Information Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Rongda Information Technology Co Ltd filed Critical Shanghai Rongda Information Technology Co Ltd
Priority to CN201810922391.0A priority Critical patent/CN109166108B/en
Publication of CN109166108A publication Critical patent/CN109166108A/en
Application granted granted Critical
Publication of CN109166108B publication Critical patent/CN109166108B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • G06T7/0014Biomedical image inspection using an image reference approach
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/194Segmentation; Edge detection involving foreground-background segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74Image or video pattern matching; Proximity measures in feature spaces
    • G06V10/75Organisation 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/751Comparing pixel values or logical combinations thereof, or feature values having positional relevance, e.g. template matching
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30061Lung

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

A kind of automatic identifying method of CT images pulmonary abnormalities tissue
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.
CN201810922391.0A 2018-08-14 2018-08-14 Automatic identification method for abnormal lung tissue of CT (computed tomography) image Active CN109166108B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (6)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
MARINA SERGEJEVA等: ""Anatomical landmarks for registration of experimental image data to volumetric rodent brain atlasing templates"", 《JOURNAL OF NEUROSCIENCE METHODS》 *
米子川: ""基于图形分析方法的函数型数据异常值检验实证研究"", 《统计与信息论坛》 *

Cited By (1)

* Cited by examiner, † Cited by third party
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