CN107865642A - A kind of glaucoma filtering operation avascular filtering bleb evaluation method based on OCT Angiographies - Google Patents
A kind of glaucoma filtering operation avascular filtering bleb evaluation method based on OCT Angiographies Download PDFInfo
- Publication number
- CN107865642A CN107865642A CN201710902054.0A CN201710902054A CN107865642A CN 107865642 A CN107865642 A CN 107865642A CN 201710902054 A CN201710902054 A CN 201710902054A CN 107865642 A CN107865642 A CN 107865642A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- dimensional
- oct
- avascular
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0062—Arrangements for scanning
- A61B5/0066—Optical coherence imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Eye Examination Apparatus (AREA)
Abstract
A kind of glaucoma filtering operation avascular filtering bleb evaluation method based on OCT Angiographies, realize that anterior ocular segment is contactless, long depth, high-speed image sampling, pass through the method for reconstructing of anterior ocular segment SS OCTA three-dimensional capillary networks, establish suitable for filtering function after analysis glaucoma filtration surgery related avascular filtering bleb three-dimensional capillary network structure and objective quantification detection method, two-dimentional capilary distribution to anterior ocular segment tissue carries out non-contact, high-resolution, perspective plane imaging containing depth information, and avoid after subject injects contrast agent and allergy occur, the risk of the complication such as shock, obtain three-dimensional capilary overall structure information, express the objective quantification parameters such as three-dimensional capillary network complex degree of structure.It is most important to the three-dimensional capillary network follow-up investigation of glaucoma filtering operation its avascular filtering bleb function correlation by the foundation of three-dimensional anterior ocular segment capillary network reconstruction technique and corresponding objective quantification index.
Description
Technical field
Present invention relates particularly to ophthalmology medical science and photoelectron technical field, and in particular to one kind is based on OCT angiogram skills
The glaucoma filtering operation avascular filtering bleb evaluation method of art.
Background technology
Glaucoma turns into the blinding illness in eye of the world's second largest, and due to its pathologic Bulbi hypertonia, optic nerve can be caused can not
Inverse property is damaged, and causes visual disorder.Glaucoma filtration surgery is the chief surgical mode for controlling intraocular pressure treatment glaucoma.Operation
Whether successful key be it is postoperative can be formed functional outflow filtration road, including conjunctiva under avascular filtering bleb and shallow-layer scleral flap
Lower drainage channel.It is reported that the mortality of 2 years reaches 15%-30%, the cicatrization and filter of avascular filtering bleb after glaucoma filtration surgery
It is the main reason for causing operative failure that passageway, which is blocked,.One overall excessively vascularization of the avascular filtering bleb formed for a long time means it
Filtering function is lost, operation control intraocular pressure failure.However, the most avascular filtering bleb evaluation of classification marks used at present in clinical position
Standard, the mark of evaluation filtering function is only used as by avascular filtering bleb surface vascularisation degree.However, follow-up observation and objective evaluation green grass or young crops
The overall capillary network structure of avascular filtering bleb after light eye filtration surgery, it can be postoperative judgement avascular filtering bleb function, take processing to arrange early
Offer more secured clinical is applied to prompt.
Blood supply in anti-glaucoma surgery filtration road repair process is essentially from tarsus artery retrotarsal fold branch and ciliary
Prerolandic artery Rolando episclera branch double route.Sought in every way to from other Post operations and promote wound healing different, anti-glaucoma surgery
The maintenance in filtration road is the process of a suppression operative incision reparation, that is, reduces avascular filtering bleb vascularization and filtration dao.
During the follow-up of glaucoma filtering operation, clinical the most frequently used method is by slit-lamp microscope (slit lamp
Microscope) take a picture, avascular filtering bleb profile and surface blood supply observed and recorded, and establish a variety of parting categorizing systems,
Including：Kronfeld parting standards, Migdal parting standards, McCulloch parting standards, IBAGS (Indiana
Bleb grading appearance scale) categorizing system and MBGS (Moorfields bleb grading system)
Categorizing system etc..However, slit-lamp microscope photograph is only capable of shooting avascular filtering bleb vascular surface, depth letter can not be provided
Breath, and related parting categorizing system is only qualitative criteria, relies on and checks diagosis person's experience, as a result subjective, nothing
The objective quantitative of legally constituted authority one records.Further study for contacting between clear and definite avascular filtering bleb internal structure and its function, usually through
Anterior ocular segment ultrasound biomicroscopy UBM (ultrasound biomicroscopy) observes avascular filtering bleb two dimensional slice structure.However, by
Limit, avascular filtering bleb capillary network can not be precisely imaged, and it uses contact inspection in the resolution ratio of ultrasound biomicroscopy UBM
Mode, avascular filtering bleb inflammatory reaction after easy stimulation, increase infection risk.In addition, Laser Scanning Confocal Microscope (confocal
Microscopy) though micro-cell level observation can be carried out to avascular filtering bleb conjunctival epithelium and upper veins beneath the skin, due to its observation
Scope and depth are limited, can not provide macroscopical Global Information of avascular filtering bleb capillary network structure, and due to its contact reviewing party
Formula, high is required to patient compliance's degree, clinical manipulation relative difficulty.It would therefore be highly desirable to develop it is a kind of it is objective, comprising depth information, height
Resolution ratio, contactless test mode, with aid in glaucoma filtering operation avascular filtering bleb integrally three-dimensional capillary network it is clinical with
Visit observation.
The progress of anterior ocular segment means of optical coherence tomography (optical coherence tomography, OCT), promote
It is developed the high-resolution imaging of anterior ocular segment tissue depth information.In recent years, it is emerging to be based on OCT three-dimensional registration weights
Technology OCT angiographic imagings (OCT angiography, OCTA) are built, can be right while being imaged to anterior ocular segment institutional framework
Vessel information is extracted, and the two-dimentional capillary network perspective plane (en face) of parallel different depth is rebuild.It is glimmering compared to anterior ocular segment
Light element blood light radiography (fluorescein angiographic, FA) and Indocyanine-Green (indocyanine
Greenangiographic, ICGC), OCTA can carry out non-contact, high score to the two-dimentional capilary distribution of anterior ocular segment tissue
Resolution, the perspective plane imaging containing depth information, and avoid after subject injects contrast agent and the complication such as allergy, shock occur
Risk.However, clinical OCT technology can only obtain two-dimensional projection's face OCTA images at present, lack overall three-dimensional capillary network and rebuild
Technology, three-dimensional capilary overall structure information can not be obtained, and lack the visitors such as the three-dimensional capillary network complex degree of structure of expression
See quantization parameter.However, the foundation of three-dimensional anterior ocular segment capillary network reconstruction technique and corresponding objective quantification index, is filtered to glaucoma
It is most important to cross the related three-dimensional capillary network follow-up investigation of its postoperative avascular filtering bleb function.
The content of the invention
The defects of in order to overcome above-mentioned prior art to exist, the present invention provide a kind of green grass or young crops based on OCT Angiographies
Light eye filtering surgery avascular filtering bleb evaluation method.
The technical solution that the present invention uses is：A kind of glaucoma filtering operation filter based on OCT Angiographies
Bubble evaluation method is crossed, described evaluation method comprises the following steps：
(1) the avascular filtering bleb image that glaucoma filtering operation is obtained using anterior ocular segment swept light source OCT (SS-OCT) system is adopted
Collection：SS-OCT IMAQs use BM scan patterns, human eye avascular filtering bleb B-scan images are obtained, using phase related algorithm pair
The B-scan images of adjacent twice sweep carry out registration；
(2) using acquisition of the Doppler variance algorithm realization to avascular filtering bleb capilary image information based on intensity, and
Avascular filtering bleb three-dimensional microvessel structure is rebuild using the dynamic algorithm of three-dimensional extended on the basis of two-dimensional projection；
(3) the avascular filtering bleb microvessel network based on reconstruction, retouched using the box counting dimension algorithm based on three-dimensional fractal dimensional analysis
The spatial complex degree of SS-OCTA three-dimensional capillary networks is stated, establishes the quantitative indices for representing avascular filtering bleb function.
Anterior ocular segment swept light source OCT (SS-OCT) system in described step (1) uses centre wavelength as 1310nm, band
Wide 100nm swept light source, sweep speed 200KHz, the OCT B-scan images size for gathering acquisition are 1024 × 512 pictures
Element.
Realized in described step (2) using the Doppler variance algorithm based on intensity to avascular filtering bleb capilary image information
Acquisition comprise the following steps：Using the time difference T of the increase method acquisition B-Scan scannings of measurement interval twice, and pixel
Dot image method for registering corrects the transverse movement x and lengthwise movement z of eyeball between different B-Scan scanning, thus, acquisition it is same
Difference (σ 2) between one position difference B-Scan is expressed as：
Wherein, M represents the pixel of mean scan depth, and N represents the number of same position scanning；AN, mRepresent OCT numbers
According to range value.
Two-dimensional projection's face image of the avascular filtering bleb capilary of described acquisition passes through gaussian filtering and median filter method pair
Image carries out denoising, carries out registration to different scanning location drawing picture using subpixel registration algorithm, is obtained using maximum value projection method
The corneoscleral junction capillary network for obtaining human eye SS-OCTA two-dimensional projections face is rebuild.
Avascular filtering bleb three-dimensional microvessel structure is rebuild including following using the dynamic algorithm of three-dimensional extended in described step (2)
Step：The OCTA data of three-dimensional are split using the dynamic programming algorithm of three-dimensional extended, obtain the three-dimensional micro- blood of SS-OCTA
Pipe network structure image, the border of destination organization is detected using the dynamic programming algorithm based on three-dimensional extended, by asking for frontier probe
The problem of topic is converted into finding shortest path expense, optimal path can efficiently be found using dynamic programming algorithm, and keep preferable
Robustness, for the OCT blood vessel datas v (x, y, z) of three-dimensional, used three-dimensional extended state transition equation can represent
For：
Wherein, C1And C2The accumulative expense in the path being illustrated respectively on x-z-plane and y-z plane, parameter d1, d2And α1,
α2For the smoothness for the target surface for controlling detection.W (x+i, x) then represents that point (x+i, y-i, z) arrives the road of point (x, y, z)
The weight in footpath, weight are obtained by following formula：
W (x+i, x)=2-r (x+i, y-1, z)-r (x, y, z)；
Wherein, r (x, y, z) is the Grad at point (x, y, z) place, and normalizes 0 to 1 section.
It is micro- using the box counting dimension arthmetic statement SS- OCTA three-dimensionals based on three-dimensional fractal dimensional analysis in described step (3)
The spatial complex degree of rete vasculosum, the quantitative indices that foundation represents avascular filtering bleb function comprise the following steps：Using based on three-dimensional
The spatial complex degree of the box counting dimension arthmetic statement SS-OCTA three-dimensional capillary networks of Fractal Dimension analysis, box counting dimension algorithm pass through
Divide the mesh generation that blood vessel segmentation is a length of r of the several bodies of N to shape object, become more meticulous as grid division is infinite, expressed point
Shape dimensional relationships are as follows：
The beneficial effects of the invention are as follows：The invention provides a kind of glaucoma filtering surgery based on OCT Angiographies
Avascular filtering bleb evaluation method afterwards, contactless anterior ocular segment, long depth, high-speed image sampling are realized, it is three-dimensional to pass through anterior ocular segment SS-OCTA
The method for reconstructing of capillary network, establish suitable for the three-dimensional micro- blood of the related avascular filtering bleb of filtering function after analysis glaucoma filtration surgery
Pipe network structure and objective quantification detection method, the two-dimentional capilary distribution to anterior ocular segment tissue carry out non-contact, high-resolution, contained
The perspective plane imaging of depth information, and avoid after subject injects contrast agent and the risk of the complication such as allergy, shock occur, obtain
Three-dimensional capilary overall structure information is taken, expresses the objective quantification parameters such as three-dimensional capillary network complex degree of structure.Pass through three-dimensional
The foundation of anterior ocular segment capillary network reconstruction technique and corresponding objective quantification index, to its avascular filtering bleb function phase of glaucoma filtering operation
The three-dimensional capillary network follow-up investigation of pass is most important.
Brief description of the drawings
Fig. 1 is the swept light source OCT system schematics that Function detection is steeped for glaucoma filtration.
Fig. 2 is the three-dimensional reconstruction figure based on SS-OCT glaucoma filtrations bubble；It is that avascular filtering bleb is scanned in SS-OCT sections wherein to scheme A
B-scan images；Scheme B avascular filtering bleb three-dimensional reconstruction images.
Fig. 3 is avascular filtering bleb SS-OCTA three-dimensional fractal dimensional analysis results；It is avascular filtering bleb conjunctival surface vessel graph wherein to scheme A
Picture；Figure B is SS-OCTA vessel extraction results；Figure C is SS-OCTA blood vessel Skeleton box counting dimension algorithm Dbox quantized results.
Embodiment
With reference to embodiment, the present invention will be described in detail, and embodiment is only the preferred embodiment of the present invention,
It is not limitation of the invention.
SS-OCTA Angiographies involved in the present invention are based on anterior ocular segment SS-OCT systems, use centre wavelength for
1310nm, bandwidth 100nm swept light source (swept source), Systems Theory axial resolution is up to 7.4 μm.SS-OCT systems
System sweep speed is 200KHz, and the OCT B-scan images size for gathering acquisition is 1024 × 512 pixels, and highest image taking speed will
Reach 200 frames/second,
SS-OCT IMAQs use BM scan patterns, obtain human eye avascular filtering bleb B-scan images.B-Scan scanning gaps
The dynamic extraction still on image information of axial eye cause necessarily to influence, move caused influence to reduce eye in scanning process, adopt
Registration is carried out to the B-scan images of adjacent twice sweep with phase related (phase correlation) algorithm.By based on
Doppler variance (intensity-based Doppler variance) algorithm of intensity is realized to avascular filtering bleb capilary image
Acquisition of information.Matched somebody with somebody using the time difference (T) of the increase method acquisition B-Scan scannings of measurement interval twice, and pixel dot image
Quasi- method corrects the transverse movement (x) and lengthwise movement (z) of eyeball between different B-Scan scannings.Thus, the same position of acquisition
The difference (σ 2) put between different B-Scan is expressed as：
Wherein, M represents the pixel of mean scan depth, and N represents the number of same position scanning.Work as M in calculating process
During with N corresponding numerical value increase, signal to noise ratio (the signal to noise ratio) increase of OCT image.Wherein M is arranged to 2, N
It is arranged to 8.AN, mRepresent the range value of OCT data.Image is eliminated using the threshold value analyzed based on histogram (histogram) to make an uproar
Sound.
The acquisition of the capillary network image in two-dimensional projection face, pass through gaussian filtering (Gauss filter) and medium filtering
(median filter) method carries out denoising to image, using subpixel registration (subpixel registration) algorithm pair
Different scanning location drawing picture carries out registration.Obtained using maximum value projection (maximum intensity projection) method
The corneoscleral junction capillary network in human eye SS-OCTA two-dimensional projections face is rebuild.
On the basis of two-dimensional projection's face capillary network is rebuild, this research will further use the Dynamic Programming of three-dimensional extended
(dynamic programming) algorithm is split to the OCTA data of three-dimensional, obtains SS- OCTA three-dimensional capillary network knots
Composition picture.Using the border of the dynamic programming algorithm detection destination organization based on three-dimensional extended, its general principle is to visit border
The problem of the problem of survey, is converted into finding shortest path (shortest path) expense.Dynamic programming algorithm can efficiently be found most
Shortest path, and keep preferable robustness.For OCT blood vessel datas/(x, the y, z) of three-dimensional, used three-dimensional extended state
Equation of transfer can be expressed as：
Wherein, C1And C2The accumulative expense in the path being illustrated respectively on x-z-plane and y-z plane, parameter d, d2And α1,
α2For the smoothness for the target surface for controlling detection.W (x+i, x) then represents that point (x+i, y-i, z) arrives the path of point (x, y, z)
Weight, can be obtained by following formula：
W (x+i, x)=2-r (x+i, y-1, z)-r (x, y, z)
Formula (4)；
Wherein, r (x, y, z) is the Grad at point (x, y, z) place, and normalizes 0 to 1 section.
Because SS-OCTA three-dimensionals capillary network can express the 3 D stereo capillary network structural information of more horn of plenty.Using
Box counting dimension algorithm (box counting, Dbox) description SS- based on three-dimensional fractal dimensional analysis (fractal analysis)
The spatial complex degree of OCTA three-dimensional capillary networks.Box counting dimension algorithm is by the way that the grid that blood vessel segmentation is a length of r of the several bodies of N is drawn
Divide shape object, become more meticulous as grid division is infinite, expressed Fractal Dimension relation is as follows：
Described above is only the preferred embodiment of the present invention, and protection scope of the present invention is not limited merely to above-mentioned implementation
Example, all technical schemes belonged under thinking of the present invention belong to protection scope of the present invention.It should be pointed out that for the art
Those of ordinary skill for, some improvements and modifications without departing from the principles of the present invention, these improvements and modifications
It should be regarded as protection scope of the present invention.
Claims (6)
1. a kind of glaucoma filtering operation avascular filtering bleb evaluation method based on OCT Angiographies, it is characterised in that described
Evaluation method comprises the following steps：
(1) the avascular filtering bleb IMAQ of glaucoma filtering operation is obtained using anterior ocular segment swept light source OCT (SS-OCT) system：
SS-OCT IMAQs use BM scan patterns, human eye avascular filtering bleb B-scan images are obtained, using phase related algorithm to adjacent
The B-scan images of twice sweep carry out registration；
(2) acquisition to avascular filtering bleb capilary image information is realized using the Doppler variance algorithm based on intensity, and in two dimension
Avascular filtering bleb three-dimensional microvessel structure is rebuild using the dynamic algorithm of three-dimensional extended on the basis of projection；
(3) the avascular filtering bleb microvessel network based on reconstruction, using the box counting dimension arthmetic statement SS- based on three-dimensional fractal dimensional analysis
The spatial complex degree of OCTA three-dimensional capillary networks, establish the quantitative indices for representing avascular filtering bleb function.
A kind of 2. glaucoma filtering operation avascular filtering bleb evaluation side based on OCT Angiographies according to claim 1
Method, it is characterised in that anterior ocular segment swept light source OCT (SS-OCT) system in described step (1) use centre wavelength for
1310nm, bandwidth 100nm swept light source, sweep speed 200KHz, the OCT B-scan image sizes for gathering acquisition are
1024 × 512 pixels.
A kind of 3. glaucoma filtering operation avascular filtering bleb evaluation side based on OCT Angiographies according to claim 1
Method, it is characterised in that realized in described step (2) using the Doppler variance algorithm based on intensity to avascular filtering bleb capilary figure
As the acquisition of information comprises the following steps：The time difference T of B-Scan scannings is obtained using the method for increasing measurement interval twice, with
And pixel method for registering images corrects the transverse movement x and lengthwise movement z of eyeball between different B-Scan scannings, thus, obtains
Difference (σ 2) between the same position difference B-Scan obtained is expressed as：
<mrow>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msup>
<mi>T</mi>
<mn>2</mn>
</msup>
</mfrac>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>m</mi>
</mrow>
</msub>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>m</mi>
</mrow>
</msub>
<mo>|</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mo>|</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>n</mi>
<mo>,</mo>
<mi>m</mi>
</mrow>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>+</mo>
<mo>|</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>m</mi>
</mrow>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mo>;</mo>
</mrow>
Wherein, M represents the pixel of mean scan depth, and N represents the number of same position scanning；AN, mRepresent the width of OCT data
Angle value.
A kind of 4. glaucoma filtering operation avascular filtering bleb evaluation side based on OCT Angiographies according to claim 3
Method, it is characterised in that two-dimensional projection's face image of the avascular filtering bleb capilary of described acquisition passes through gaussian filtering and medium filtering
Method carries out denoising to image, carries out registration to different scanning location drawing picture using subpixel registration algorithm, is thrown using maximum
The corneoscleral junction capillary network that shadow method obtains human eye SS-OCTA two-dimensional projections face is rebuild.
A kind of 5. glaucoma filtering operation avascular filtering bleb evaluation side based on OCT Angiographies according to claim 1
Method, it is characterised in that avascular filtering bleb three-dimensional microvessel structure bag is rebuild using the dynamic algorithm of three-dimensional extended in described step (2)
Include following steps：The OCTA data of three-dimensional are split using the dynamic programming algorithm of three-dimensional extended, it is three-dimensional to obtain SS-OCTA
Capillary network structural images, the border of destination organization is detected using the dynamic programming algorithm based on three-dimensional extended, by frontier probe
The problem of be converted into find shortest path expense the problem of, optimal path can efficiently be found using dynamic programming algorithm, and keep
Preferable robustness, for the OCT blood vessel datas v (x, y, z) of three-dimensional, used three-dimensional extended state transition equation can be with table
It is shown as：
<mrow>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mrow>
<mo>-</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mo>&le;</mo>
<mi>i</mi>
<mo>&le;</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
</mrow>
</munder>
<mo>{</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mi>i</mi>
<mo>,</mo>
<mi>y</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
<mi>i</mi>
<mo>|</mo>
<mo>+</mo>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mi>i</mi>
<mo>,</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mrow>
<mo>-</mo>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
<mo>&le;</mo>
<mi>i</mi>
<mo>&le;</mo>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
</mrow>
</munder>
<mo>{</mo>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>y</mi>
<mo>+</mo>
<mi>j</mi>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<mo>|</mo>
<mi>j</mi>
<mo>|</mo>
<mo>}</mo>
<mo>+</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Wherein, C1And C2The accumulative expense in the path being illustrated respectively on x-z-plane and y-z plane, parameter d1, d2And α1, α2With
In the smoothness of the target surface of control detection.W (x+i, x) then represents point (x+i, y-i, z) to the power in the path of point (x, y, z)
Weight, weight are obtained by following formula：
W (x+i, x)=2-r (x+i, y-1, z)-r (x, y, z)；
Wherein, r (x, y, z) is the Grad at point (x, y, z) place, and normalizes 0 to 1 section.
A kind of 6. glaucoma filtering operation avascular filtering bleb evaluation side based on OCT Angiographies according to claim 1
Method, it is characterised in that the box counting dimension arthmetic statement SS-OCTA based on three-dimensional fractal dimensional analysis is used in described step (3)
The spatial complex degree of three-dimensional capillary network, the quantitative indices that foundation represents avascular filtering bleb function comprise the following steps：Using base
In the spatial complex degree of the box counting dimension arthmetic statement SS-OCTA three-dimensional capillary networks of three-dimensional fractal dimensional analysis, box counting dimension is calculated
Method by by blood vessel segmentation be a length of r of the several bodies of N mesh generation divide shape object, become more meticulous as grid division is infinite, institute's table
The Fractal Dimension relation reached is as follows：
<mrow>
<msub>
<mi>dim</mi>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mi>lim</mi>
<mrow>
<mi>r</mi>
<mo>&RightArrow;</mo>
<mn>0</mn>
</mrow>
</munder>
<mfrac>
<mrow>
<mi>log</mi>
<mi> </mi>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>log</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mi>r</mi>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>.</mo>
</mrow>
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710902054.0A CN107865642A (en) | 2017-09-28 | 2017-09-28 | A kind of glaucoma filtering operation avascular filtering bleb evaluation method based on OCT Angiographies |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710902054.0A CN107865642A (en) | 2017-09-28 | 2017-09-28 | A kind of glaucoma filtering operation avascular filtering bleb evaluation method based on OCT Angiographies |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107865642A true CN107865642A (en) | 2018-04-03 |
Family
ID=61756742
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710902054.0A Pending CN107865642A (en) | 2017-09-28 | 2017-09-28 | A kind of glaucoma filtering operation avascular filtering bleb evaluation method based on OCT Angiographies |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107865642A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108902013A (en) * | 2018-05-31 | 2018-11-30 | 山东大学齐鲁医院 | A kind of method for building up of macular edema animal model |
| CN111243087A (en) * | 2020-01-20 | 2020-06-05 | 上海市第一人民医院 | Three-dimensional reconstruction method and device for fundus blood vessels and electronic equipment |
| CN112085830A (en) * | 2019-06-14 | 2020-12-15 | 北京大学 | Optical coherent angiography imaging method based on machine learning |
-
2017
- 2017-09-28 CN CN201710902054.0A patent/CN107865642A/en active Pending
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108902013A (en) * | 2018-05-31 | 2018-11-30 | 山东大学齐鲁医院 | A kind of method for building up of macular edema animal model |
| CN108902013B (en) * | 2018-05-31 | 2022-02-15 | 山东大学齐鲁医院 | Method for establishing retinal edema animal model |
| CN112085830A (en) * | 2019-06-14 | 2020-12-15 | 北京大学 | Optical coherent angiography imaging method based on machine learning |
| CN112085830B (en) * | 2019-06-14 | 2024-02-27 | 北京大学 | Optical coherence angiography imaging method based on machine learning |
| CN111243087A (en) * | 2020-01-20 | 2020-06-05 | 上海市第一人民医院 | Three-dimensional reconstruction method and device for fundus blood vessels and electronic equipment |
| CN111243087B (en) * | 2020-01-20 | 2023-11-21 | 上海市第一人民医院 | Three-dimensional reconstruction method and device for fundus blood vessel and electronic equipment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP7250653B2 (en) | Image processing device, image processing method and program | |
| Petroll et al. | In vivo confocal microscopy of the cornea: new developments in image acquisition, reconstruction, and analysis using the HRT-Rostock corneal module | |
| CN104271031B (en) | The analysis and visualization of OCT angiographic datas | |
| CN112601487A (en) | Medical image processing apparatus, medical image processing method, and program | |
| JP2014527434A (en) | Feature motion correction and normalization in optical coherence tomography | |
| JP6758826B2 (en) | Image processing device and image processing method | |
| JP6843521B2 (en) | Image processing device and image processing method | |
| Migacz et al. | Imaging of vitreous cortex hyalocyte dynamics using non-confocal quadrant-detection adaptive optics scanning light ophthalmoscopy in human subjects | |
| JP7374615B2 (en) | Information processing device, information processing method and program | |
| JP6940349B2 (en) | Ophthalmic equipment | |
| CN107865642A (en) | A kind of glaucoma filtering operation avascular filtering bleb evaluation method based on OCT Angiographies | |
| JP7199236B2 (en) | ophthalmic equipment | |
| JP7348374B2 (en) | Ophthalmology information processing device, ophthalmology imaging device, ophthalmology information processing method, and program | |
| JP7090438B2 (en) | Ophthalmologic imaging equipment, its control method, programs, and recording media | |
| JP2021122559A (en) | Image processing device, image processing method, and program | |
| JP2020163100A (en) | Image processing apparatus and image processing method | |
| Khalil et al. | An overview of automated glaucoma detection | |
| Bhandari et al. | Quantitative analysis of the lamina cribrosa in vivo using a scanning laser opthalmoscope | |
| WO2019230643A1 (en) | Information processing device, information processing method, and program | |
| WO2020075719A1 (en) | Image processing device, image processing method, and program | |
| JP2023014190A (en) | Ophthalmology imaging apparatus | |
| WO2019157113A1 (en) | Segmentation-based corneal mapping | |
| Dhommati et al. | Automated 2D-3D quantitative analysis of corneal graft detachment post DSAEK based on AS-OCT images | |
| WO2022196583A1 (en) | Grade assessment device, opthalmic imaging device, program, recording medium and grade assessment method | |
| Pchelkina et al. | An Algorithm for Automatic Image Segmentation Using the Sobel Method for an Optical Coherence Tomography |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180403 |
|
| RJ01 | Rejection of invention patent application after publication |