CN103699901A - Automatic detection method for IS/OS (intermediate system/operating system) missing area in 3D (three-dimensional) OCT (optical coherence tomography) retina image based on support vector machine - Google Patents
Automatic detection method for IS/OS (intermediate system/operating system) missing area in 3D (three-dimensional) OCT (optical coherence tomography) retina image based on support vector machine Download PDFInfo
- Publication number
- CN103699901A CN103699901A CN201310692851.2A CN201310692851A CN103699901A CN 103699901 A CN103699901 A CN 103699901A CN 201310692851 A CN201310692851 A CN 201310692851A CN 103699901 A CN103699901 A CN 103699901A
- Authority
- CN
- China
- Prior art keywords
- retina
- test result
- support vector
- vector machine
- region
- 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
- 210000001525 retina Anatomy 0.000 title claims abstract description 29
- 238000012706 support-vector machine Methods 0.000 title claims abstract description 13
- 238000012014 optical coherence tomography Methods 0.000 title abstract description 14
- 238000001514 detection method Methods 0.000 title abstract description 3
- 238000012360 testing method Methods 0.000 claims abstract description 42
- 230000004256 retinal image Effects 0.000 claims abstract description 13
- 238000012549 training Methods 0.000 claims abstract description 13
- 238000000034 method Methods 0.000 claims abstract description 11
- 230000008034 disappearance Effects 0.000 claims description 24
- 238000000605 extraction Methods 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 3
- 230000032798 delamination Effects 0.000 claims description 3
- 230000004438 eyesight Effects 0.000 abstract description 3
- 230000000007 visual effect Effects 0.000 abstract description 2
- 230000002207 retinal effect Effects 0.000 abstract 4
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 abstract 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 208000017442 Retinal disease Diseases 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000000574 ganglionic effect Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Landscapes
- Eye Examination Apparatus (AREA)
Abstract
The invention discloses an automatic detection method of an IS/OS (intermediate system/operating system) missing area in a 3D OCT (three-dimensional optical coherence tomography) retinal image based on a support vector machine, which comprises the steps of acquiring the retinal three-dimensional image data of a subject; layering the three-dimensional image data of the retina to find out an IS/OS area of the retina; manually and manually calibrating a missing area of the retinal IS/OS area to serve as a gold standard; selecting features to extract the features of each pixel point of the retinal IS/OS region to construct a feature set; constructing a training set and a test set; training and testing the training set and the testing set by using a support vector machine to obtain a testing result of the testing set; and comparing the test result with the gold standard to obtain the accuracy information of the test result, and feeding the test result back to the retina three-dimensional image data to finish the visual display of the test result. The method is simple and feasible; so that the quantitative expression of the relationship between retinal IS/OS loss and vision IS obtained.
Description
Technical field
The present invention relates to IS/OS disappearance region automatic testing method in a kind of 3D OCT retinal images based on support vector machine, belong to Biologic Medical Image processing technology field.
Background technology
Growing perfect along with retina checkout equipment, 3D OCT(optical coherence tomography, optical coherence photography) used more and more widely.3D OCT helps ophthalmologist, oculist to observe better amphiblestroid form, and the generation development of retinal disease is made to judgement more accurately.
At present, IS/OS(Inner segment/outer segment layer in retinal images, inside/outside ganglionic layer) lack the artificial visual check that method for detecting area mainly depends on doctor, greatly strengthened doctor's workload.。
In addition, many experts and scholars, by research, find that the disappearance of retina IS/OS and the disappearance of human eyesight exist certain relation both at home and abroad, but at present, do not have any documents and materials and show, someone did the research of quantitative relationship between IS/OS disappearance and anopsia.
Summary of the invention
The deficiency existing for prior art, the object of the invention is to provide IS/OS disappearance region automatic testing method in a kind of 3D OCT retinal images based on support vector machine of simple possible.
To achieve these goals, the present invention realizes by the following technical solutions:
The present invention includes following step:
(1) by 3D OCT equipment, obtain at least ten experimenters' retina three-dimensional view data;
(2) utilize retina delamination software, described retina three-dimensional view data is carried out to layering, find retina IS/OS region;
(3) utilize CAVASS software, the disappearance region in described retina IS/OS region is demarcated manually, as goldstandard;
(4) choose at least seven kinds of features on different stage, each pixel in described retina IS/OS region is carried out to feature extraction, construction feature collection;
(5) use leaving-one method, using the feature set of retina three-dimensional view data described in one of them experimenter as test data at every turn, using all the other remaining feature sets as training set, complete the structure of training set and test set;
(6) utilize support vector machine to carry out training and testing to described training set and test set, can obtain the test result of test set;
(7) described test result and goldstandard are compared, obtain the accuracy rate information of described test result, described test result is fed back in the retina three-dimensional view data in step (1) simultaneously, can provide one and show intuitively.
Above-mentioned experimenter is provided with ten, and chooses ten kinds of features on different stage.
Above-mentioned accuracy rate information comprises accuracy rate, true positives and true negative.
The present invention, by the extraction of feature, utilizes support vector machine to train test, can complete disappearance region robotization detection and result and show, the method simple possible; In addition, the present invention is by carrying out feature extraction to the pixel in retina IS/OS region seriatim, make " retina IS/OS disappearance exists certain relation with eyesight " this qualitative results, obtained quantitative expression, indirectly, by this quantitative relationship, assist a physician the ophthalmology disease of retina IS/OS disappearance is judged in advance.
Accompanying drawing explanation
Fig. 1 is the some sectioning images (in figure, indication closed curve scope A is disappearance) in 3D OCT retinal images;
Fig. 2 is workflow diagram of the present invention;
Fig. 3 is 13 displacement vectors that direction is corresponding;
Fig. 4 is that goldstandard result shows (in figure, white portion is shown as pixel);
Fig. 5 is that the test result of respective regions shows (in figure, white portion is shown as pixel).
Embodiment
For technological means, creation characteristic that the present invention is realized, reach object and effect is easy to understand, below in conjunction with embodiment, further set forth the present invention.
Referring to Fig. 1 and Fig. 2, first retina IS/OS disappearance of the present invention detects carries out the collection of 3DOCT view data, by the layering to retinal images, completes determining of retina IS/OS region.The disappearance position of image and non-disappearance position are demarcated, obtained goldstandard.By the feature extraction to each pixel, utilize SVM to train, test, draws last testing result.
The method of the invention is to obtain retinal images on the 3D OCT equipment of analyzing, and IS/OS disappearance detects principle and mainly utilized ten kinds of features on different scale, utilizes this machine learning method of SVM to analyze.
This method by obtaining amphiblestroid 3-D view on 3D OCT equipment, use retina delamination software, find IS/OS region, under experienced expert's help, utilize CAVASS software (a Computer Assisted Visualization and Analysis Software System), mark is manually carried out in disappearance region to retina IS/OS region, as goldstandard.
To each the pixel construction feature collection in IS/OS region, selected altogether ten kinds of features (these ten kinds of features are existing feature) herein, specific as follows:
< feature group 1>
(1) normalized gray-scale value:
Wherein, I
original(i, j, k) is the original gray-scale value of object pixel (i, j, k)
I
minit is the minimum gray-scale value in this retinal image data
I
maxit is the minimum gray-scale value in this retinal image data
I
normalized(i, j, k) is the gray-scale value after object pixel (i, j, k) normalization
(2) average (piece average) in subregion (region of 5*5*5 pixel, object pixel is positioned at center):
M
block(i, j, k) is the mean value of all pixel gray-scale values in the space in the 5*5*5 region centered by object pixel (i, j, k).
(3) standard deviation in subregion (region of 5*5*5 pixel, object pixel is positioned at center):
STD
block(i, j, k) is the standard deviation of the gray-scale value of all pixels in the space in the 5*5*5 region centered by object pixel (i, j, k).
(4) entropy in subregion (region of 5*5*5 pixel, object pixel is positioned at center):
The quantity of information that entropy Description Image has, shows the complicated process of image, and image complexity is high, and entropy is larger, otherwise less
Wherein, r
mit is the gray-scale value of pixel
P(r
m) be the probability that this gray-scale value occurs
ENT
blockit is the information entropy in this 5*5*5 region.
< feature group 2>
Because data are three-dimensional, so 13 directions (referring to Fig. 3) that can be defined as follows.
1. directions X, Y-direction, Z direction, totally 3;
2. X-Y plane diagonal, Y-Z plane diagonal, X-Z plane diagonal, totally 6;
3. body diagonal direction, totally 4.
Be defined as follows shown in Fig. 3: α is the projection of this direction in X-Y direction and the angle of X-axis, β is the angle of this direction and Z axis, and D is step-length.
(1) step-length is 1 o'clock, the gray-scale value antipode of the pixel in 13 directions of object pixel:
(2) step-length is 2 o'clock, the gray-scale value antipode of the pixel in 13 directions of object pixel:
< feature group 3>
By build the region of 5*5*5 centered by object pixel, obtain three-dimensional gray level co-occurrence matrixes, and calculate corresponding specific features, what choose is that direction 1 builds gray level co-occurrence matrixes herein.After completing the structure of gray level co-occurrence matrixes, the gray level co-occurrence matrixes fundamental function that can carry by MATLAB, calculates four category features below:
(7) contrast (contrast) in subregion (region of 5*5*5 pixel, object pixel is positioned at center)
(8) autocorrelation (correlation) in subregion (region of 5*5*5 pixel, object pixel is positioned at center)
(9) energy (energy) in subregion (region of 5*5*5 pixel, object pixel is positioned at center)
(10) homogeneity (homogeneity) in subregion (region of 5*5*5 pixel, object pixel is positioned at center) so far, has completed the structure of feature set.
Because the present embodiment data set used is limited, so locate to use existing method leaving-one method (leave-one-out) to carry out training and testing data.Support vector machines is a kind of basic, novel small-sample learning method that has solid theory, utilize it to carry out training and testing and can obtain test result, add up and can obtain corresponding accuracy rate, true positives, true negative, test result also can feed back in original data and go simultaneously, provides one and shows intuitively (referring to Fig. 4 and Fig. 5).
TP: the pixel count that disappearance correctly detected
FP: flase drop is the pixel count of disappearance
TN: the pixel count that non-disappearance correctly detected
FN: the pixel count that flase drop is non-disappearance
True positives:
True negative:
Accuracy rate:
More than show and described ultimate principle of the present invention and principal character and advantage of the present invention.The technician of the industry should understand; the present invention is not restricted to the described embodiments; that in above-described embodiment and instructions, describes just illustrates principle of the present invention; without departing from the spirit and scope of the present invention; the present invention also has various changes and modifications, and these changes and improvements all fall in the claimed scope of the invention.The claimed scope of the present invention is defined by appending claims and equivalent thereof.
Claims (3)
1. an IS/OS disappearance region automatic testing method in the 3D OCT retinal images based on support vector machine, is characterized in that, comprises following step:
(1) by 3D OCT equipment, obtain at least ten experimenters' retina three-dimensional view data;
(2) utilize retina delamination software, described retina three-dimensional view data is carried out to layering, find retina IS/OS region;
(3) utilize CAVASS software, the disappearance region in described retina IS/OS region is demarcated manually, as goldstandard;
(4) choose at least seven kinds of features on different stage, each pixel in described retina IS/OS region is carried out to feature extraction, construction feature collection;
(5) use leaving-one method, using the feature set of retina three-dimensional view data described in one of them experimenter as test data at every turn, using all the other remaining feature sets as training set, complete the structure of training set and test set;
(6) utilize support vector machine to carry out training and testing to described training set and test set, can obtain the test result of test set;
(7) described test result and goldstandard are compared, obtain the accuracy rate information of described test result, described test result is fed back in the retina three-dimensional view data in step (1) simultaneously, complete the demonstration directly perceived of test result.
2. IS/OS disappearance region automatic testing method in the 3D OCT retinal images based on support vector machine according to claim 1, is characterized in that,
Described experimenter is provided with ten, and chooses ten kinds of features on different stage.
3. IS/OS disappearance region automatic testing method in the 3D OCT retinal images based on support vector machine according to claim 1, is characterized in that,
Described accuracy rate information comprises accuracy rate, true positives and true negative.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310692851.2A CN103699901A (en) | 2013-12-17 | 2013-12-17 | Automatic detection method for IS/OS (intermediate system/operating system) missing area in 3D (three-dimensional) OCT (optical coherence tomography) retina image based on support vector machine |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310692851.2A CN103699901A (en) | 2013-12-17 | 2013-12-17 | Automatic detection method for IS/OS (intermediate system/operating system) missing area in 3D (three-dimensional) OCT (optical coherence tomography) retina image based on support vector machine |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103699901A true CN103699901A (en) | 2014-04-02 |
Family
ID=50361423
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310692851.2A Pending CN103699901A (en) | 2013-12-17 | 2013-12-17 | Automatic detection method for IS/OS (intermediate system/operating system) missing area in 3D (three-dimensional) OCT (optical coherence tomography) retina image based on support vector machine |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103699901A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105551038A (en) * | 2015-12-14 | 2016-05-04 | 苏州大学 | Method for fully automatically classifying and segmenting retinal branch artery obstruction based on three-dimensional OCT image |
CN107495923A (en) * | 2017-08-03 | 2017-12-22 | 苏州大学 | A kind of method for measuring eyeball retina form |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101778593A (en) * | 2007-06-15 | 2010-07-14 | 南加州大学 | Pattern analysis of retinal maps for diagnosis of optic nerve diseases by optical coherence tomography |
US20120148130A1 (en) * | 2010-12-09 | 2012-06-14 | Canon Kabushiki Kaisha | Image processing apparatus for processing tomographic image of subject's eye, imaging system, method for processing image, and recording medium |
US20120165799A1 (en) * | 2010-12-27 | 2012-06-28 | Nidek Co., Ltd. | Ophthalmic laser treatment apparatus |
WO2012170722A2 (en) * | 2011-06-07 | 2012-12-13 | California Institute Of Technology | Enhanced optical angiography using intensity contrast and phase contrast imaging methods |
-
2013
- 2013-12-17 CN CN201310692851.2A patent/CN103699901A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101778593A (en) * | 2007-06-15 | 2010-07-14 | 南加州大学 | Pattern analysis of retinal maps for diagnosis of optic nerve diseases by optical coherence tomography |
US20120148130A1 (en) * | 2010-12-09 | 2012-06-14 | Canon Kabushiki Kaisha | Image processing apparatus for processing tomographic image of subject's eye, imaging system, method for processing image, and recording medium |
US20120165799A1 (en) * | 2010-12-27 | 2012-06-28 | Nidek Co., Ltd. | Ophthalmic laser treatment apparatus |
WO2012170722A2 (en) * | 2011-06-07 | 2012-12-13 | California Institute Of Technology | Enhanced optical angiography using intensity contrast and phase contrast imaging methods |
Non-Patent Citations (3)
Title |
---|
JÓZSEF MOLNÁR .ETC: ""Layer extraction in rodent retinal images acquired by optical coherence tomography"", 《MACHINE VISION AND APPLICATIONS》 * |
ROBERT J. ZAWADZKI .ETC: ""Adaptation of a support vector machine algorithm for segmentation and visualization of retinal structures in volumetric optical coherence tomography data sets"", 《JOURNAL OF BIOMEDICAL OPTICS》 * |
杨丽亚 等: ""黄斑区光感受器内外节层局部微小缺失初步分析"", 《中国实用眼科杂志》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105551038A (en) * | 2015-12-14 | 2016-05-04 | 苏州大学 | Method for fully automatically classifying and segmenting retinal branch artery obstruction based on three-dimensional OCT image |
CN105551038B (en) * | 2015-12-14 | 2018-11-30 | 苏州大学 | Method for fully automatically classifying and segmenting retinal branch artery obstruction based on three-dimensional OCT image |
CN107495923A (en) * | 2017-08-03 | 2017-12-22 | 苏州大学 | A kind of method for measuring eyeball retina form |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109035187B (en) | Medical image labeling method and device | |
CN107480677B (en) | Method and device for identifying interest region in three-dimensional CT image | |
EP3477589B1 (en) | Method of processing medical image, and medical image processing apparatus performing the method | |
CN106909778A (en) | A kind of Multimodal medical image recognition methods and device based on deep learning | |
CN109770932A (en) | The processing method of multi-modal brain neuroblastoma image feature | |
CN104885126B (en) | The Computer assisted identification of tissue of interest | |
CN106415555A (en) | System and method for correlation of pathology reports and radiology reports | |
US9002083B2 (en) | System, method, and software for optical device recognition association | |
CN111862020B (en) | Method and device for predicting physiological age of anterior ocular segment, server and storage medium | |
CN109512464A (en) | A kind of disorder in screening and diagnostic system | |
CN102693353A (en) | Method and computer system for automatically generating a statistical model | |
CN112819818B (en) | Image recognition module training method and device | |
US11049255B2 (en) | Image processing device and method thereof | |
US20230057933A1 (en) | Storage medium, diagnosis support apparatus and diagnosis support method | |
CN109255354A (en) | medical CT-oriented computer image processing method and device | |
CN106296685A (en) | Based on cranium outline feature geometries invariant multi information 3D medical image registration method | |
CN103699901A (en) | Automatic detection method for IS/OS (intermediate system/operating system) missing area in 3D (three-dimensional) OCT (optical coherence tomography) retina image based on support vector machine | |
US11416994B2 (en) | Method and system for detecting chest x-ray thoracic diseases utilizing multi-view multi-scale learning | |
JP2008132320A (en) | Image diagnosis support system | |
CN111951950B (en) | Three-dimensional data medical classification system based on deep learning | |
WO2023036243A1 (en) | Medical devices, methods and systems for monitoring the medical devices | |
CN111904450B (en) | Extraction method, device and system for center of left ventricle and region of interest | |
RU120799U1 (en) | INTEREST AREA SEARCH SYSTEM IN THREE-DIMENSIONAL MEDICAL IMAGES | |
KR20230049938A (en) | Method and apparatus for quantitative analysis of emphysema | |
CN107256544A (en) | A kind of prostate cancer image diagnosing method and system based on VCG16 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140402 |