WO2022113409A1 - 画像処理方法、画像処理装置、及びプログラム - Google Patents

画像処理方法、画像処理装置、及びプログラム Download PDF

Info

Publication number
WO2022113409A1
WO2022113409A1 PCT/JP2021/023441 JP2021023441W WO2022113409A1 WO 2022113409 A1 WO2022113409 A1 WO 2022113409A1 JP 2021023441 W JP2021023441 W JP 2021023441W WO 2022113409 A1 WO2022113409 A1 WO 2022113409A1
Authority
WO
WIPO (PCT)
Prior art keywords
image
image processing
light
volume data
oct
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.)
Ceased
Application number
PCT/JP2021/023441
Other languages
English (en)
French (fr)
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.)
Nikon Corp
Original Assignee
Nikon Corp
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 Nikon Corp filed Critical Nikon Corp
Publication of WO2022113409A1 publication Critical patent/WO2022113409A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions

Definitions

  • the techniques of the present disclosure relate to image processing methods, image processing devices, and programs.
  • U.S. Patent No. 10238281 relates to a technique for generating volume data of an eye to be inspected using an optical coherence tomography. Conventionally, it has been desired to analyze the volume data of the eye to be inspected.
  • a first aspect of the technique of the present disclosure is an image processing method performed by a processor, in which a step of acquiring OCT volume data of a region including a vortex vein and a plurality of different depths based on the OCT volume data. Includes a step of generating a plurality of en-face images corresponding to a surface.
  • the image processing apparatus of the second aspect of the technique of the present disclosure includes a memory and a processor connected to the memory, wherein the processor acquires OCT volume data of a region including a vortex vein, and the OCT volume.
  • a step of generating a plurality of en-face images corresponding to a plurality of surfaces having different depths based on the data is executed.
  • the program of the third aspect of the technique of the present disclosure includes a step of acquiring OCT volume data of a region including a vortex vein in a computer, and a plurality of surfaces corresponding to a plurality of surfaces having different depths based on the OCT volume data. And to execute the step of generating the en-face image of.
  • FIG. 1 shows a schematic configuration of an ophthalmic system 100.
  • the ophthalmology system 100 includes an ophthalmology device 110, a server device (hereinafter referred to as “server”) 140, and a display device (hereinafter referred to as “viewer”) 150.
  • the ophthalmic apparatus 110 acquires a fundus image.
  • the server 140 corresponds to a patient ID with a plurality of fundus images obtained by photographing the fundus of a plurality of patients by the ophthalmologic apparatus 110 and an axial length measured by an axial length measuring apparatus (not shown). And remember.
  • the viewer 150 displays the fundus image and the analysis result acquired by the server 140.
  • the server 140 is an example of the "image processing device" of the technique of the present disclosure.
  • the ophthalmic apparatus 110, the server 140, and the viewer 150 are connected to each other via the network 130.
  • the network 130 is an arbitrary network such as LAN, WAN, the Internet, and a wide area Ethernet network.
  • a LAN can be adopted for the network 130.
  • the viewer 150 is a client in a client-server system, and a plurality of viewers are connected via a network. Further, a plurality of servers 140 may be connected via a network in order to ensure system redundancy.
  • the ophthalmology device 110 has an image processing function and an image viewing function of the viewer 150, the ophthalmology device 110 can acquire an image of the fundus of the eye, process an image, and view an image in a stand-alone state. In this case, the viewer 150 may be omitted.
  • the server 140 has an image viewing function of the viewer 150, the configuration of the ophthalmic apparatus 110 and the server 140 enables acquisition, image processing, and image viewing of the fundus image. In this case, the viewer 150 may be omitted.
  • a diagnostic support device that performs image analysis using other ophthalmic devices (inspection devices such as visual field measurement and intraocular pressure measurement) and AI (Artificial Integrity) is connected to the ophthalmic device 110, the server 140, and the viewer via the network 130. It may be connected to 150.
  • ophthalmic devices such as visual field measurement and intraocular pressure measurement
  • AI Artificial Integrity
  • the scanning laser ophthalmoscope is referred to as "SLO”.
  • the optical coherence tomography is referred to as "OCT”.
  • the horizontal direction is the "X direction” and the direction perpendicular to the horizontal plane is the "Y direction", connecting the center of the pupil of the anterior segment of the eye 12 to the center of the eyeball.
  • the direction is "Z direction”. Therefore, the X, Y, and Z directions are perpendicular to each other.
  • the ophthalmology device 110 includes a photographing device 14 and a control device 16.
  • the imaging device 14 includes an SLO unit 18 and an OCT unit 20 to acquire a fundus image of the fundus of the eye to be inspected 12.
  • the two-dimensional fundus image acquired by the SLO unit 18 is referred to as an SLO image.
  • a tomographic image of the retina or a frontal image (en-face image) created based on the OCT data acquired by the OCT unit 20 is referred to as an OCT image.
  • the control device 16 includes a computer having a CPU (Central Processing Unit) 16A, a RAM (Random Access Memory) 16B, a ROM (Read-Only memory) 16C, and an input / output (I / O) port 16D. ing.
  • CPU Central Processing Unit
  • RAM Random Access Memory
  • ROM Read-Only memory
  • I / O input / output
  • the control device 16 includes an input / display device 16E connected to the CPU 16A via the I / O port 16D.
  • the input / display device 16E has a graphic user interface for displaying an image of the eye to be inspected 12 and receiving various instructions from the user.
  • the graphic user interface includes a touch panel display.
  • control device 16 includes an image processor 17 connected to the I / O port 16D.
  • the image processor 17 generates an image of the eye 12 to be inspected based on the data obtained by the photographing apparatus 14.
  • the image processor 17 may be omitted, and the CPU 16A may generate an image of the eye 12 to be inspected based on the data obtained by the photographing apparatus 14.
  • the control device 16 is connected to the network 130 via the communication interface 16F.
  • the control device 16 of the ophthalmic device 110 includes the input / display device 16E, but the technique of the present disclosure is not limited to this.
  • the control device 16 of the ophthalmic apparatus 110 may not include the input / display device 16E, but may include an input / display device that is physically independent of the ophthalmic apparatus 110.
  • the display device includes an image processing processor unit that operates under the control of the display control unit 204 of the CPU 16A of the control device 16.
  • the image processor unit may display an SLO image or the like based on the image signal output instructed by the display control unit 204.
  • the photographing device 14 operates under the control of the CPU 16A of the control device 16.
  • the imaging device 14 includes an SLO unit 18, an imaging optical system 19, and an OCT unit 20.
  • the photographing optical system 19 includes an optical scanner 22 and a wide-angle optical system 30.
  • the optical scanner 22 two-dimensionally scans the light emitted from the SLO unit 18 in the X direction and the Y direction.
  • the optical scanner 22 may be any optical element capable of deflecting the luminous flux, and for example, a polygon mirror, a galvano mirror, or the like can be used. Moreover, it may be a combination thereof.
  • the wide-angle optical system 30 synthesizes the light from the SLO unit 18 and the light from the OCT unit 20.
  • the wide-angle optical system 30 may be a catadioptric system using a concave mirror such as an elliptical mirror, a catadioptric system using a wide-angle lens or the like, or a catadioptric system in which a concave mirror or a lens is combined.
  • a wide-angle optical system using an elliptical mirror, a wide-angle lens, or the like it is possible to photograph the retina not only in the central part of the fundus but also in the peripheral part of the fundus.
  • the wide-angle optical system 30 realizes observation in a wide field of view (FOV: Field of View) 12A at the fundus.
  • the FOV 12A indicates a range that can be photographed by the photographing apparatus 14.
  • the FOV12A can be expressed as a viewing angle.
  • the viewing angle may be defined by an internal irradiation angle and an external irradiation angle in the present embodiment.
  • the external irradiation angle is an irradiation angle in which the irradiation angle of the luminous flux emitted from the ophthalmic apparatus 110 to the eye 12 to be inspected is defined with reference to the pupil 27.
  • the internal irradiation angle is an irradiation angle in which the irradiation angle of the luminous flux irradiated to the fundus F is defined with reference to the center O of the eyeball.
  • the external irradiation angle and the internal irradiation angle have a corresponding relationship. For example, when the external irradiation angle is 120 degrees, the internal irradiation angle corresponds to about 160 degrees. In this embodiment, the internal irradiation angle is set to 200 degrees.
  • UWF-SLO fundus image the SLO fundus image obtained by taking a picture with an internal irradiation angle of 160 degrees or more is referred to as a UWF-SLO fundus image.
  • UWF is an abbreviation for UltraWide Field (ultra-wide-angle).
  • FOV viewing angle
  • the ophthalmic apparatus 110 can image a region 12A having an internal irradiation angle of 200 ° with the eyeball center O of the eye 12 to be inspected as a reference position.
  • the internal irradiation angle of 200 ° is 110 ° at the external irradiation angle based on the pupil of the eyeball of the eye 12 to be inspected. That is, the wide-angle optical system 30 irradiates the laser beam from the pupil with an angle of view of 110 ° as an external irradiation angle, and photographs the fundus region of 200 ° with an internal irradiation angle.
  • the SLO system is realized by the control device 16, the SLO unit 18, and the photographing optical system 19 shown in FIG. Since the SLO system includes a wide-angle optical system 30, it enables fundus photography with a wide FOV12A.
  • the SLO unit 18 includes a B light (blue light) light source 40, a G light (green light) light source 42, an R light (red light) light source 44, and an IR light (infrared light (for example, near infrared light)). It includes a light source 46 and optical systems 48, 50, 52, 54, 56 that reflect or transmit light from the light sources 40, 42, 44, 46 and guide them into one optical path.
  • the optical systems 48 and 56 are mirrors, and the optical systems 50, 52 and 54 are beam splitters.
  • B light is reflected by the optical system 48, is transmitted through the optical system 50, is reflected by the optical system 54, G light is reflected by the optical systems 50 and 54, and R light is transmitted through the optical systems 52 and 54.
  • IR light is reflected by the optical systems 52 and 56 and guided to one optical path, respectively.
  • the SLO unit 18 is configured to be able to switch a combination of a light source that emits laser light having a different wavelength or a light source that emits light, such as a mode that emits R light and G light and a mode that emits infrared light.
  • a light source 40 for B light includes four light sources, a light source 40 for B light, a light source 42 for G light, a light source 44 for R light, and a light source 46 for IR light, but the technique of the present disclosure is not limited thereto.
  • the SLO unit 18 is further provided with a light source of white light, and may emit light in various modes such as a mode of emitting G light, R light, and B light, and a mode of emitting only white light. good.
  • the light incident on the photographing optical system 19 from the SLO unit 18 is scanned in the X direction and the Y direction by the optical scanner 22.
  • the scanning light is applied to the fundus through the wide-angle optical system 30 and the pupil 27.
  • the reflected light reflected by the fundus is incident on the SLO unit 18 via the wide-angle optical system 30 and the optical scanner 22.
  • the SLO unit 18 transmits G light among the light from the rear eye portion (fundus) of the eye 12 to be examined, the beam splitter 64 that reflects the B light and transmits other than the B light, and the light transmitted through the beam splitter 64.
  • a beam splitter 58 that reflects and transmits other than G light is provided.
  • the SLO unit 18 includes a beam splitter 60 that reflects R light and transmits other than R light among the light transmitted through the beam splitter 58.
  • the SLO unit 18 includes a beam splitter 62 that reflects IR light among the light transmitted through the beam splitter 60.
  • the SLO unit 18 detects the B light detection element 70 that detects the B light reflected by the beam splitter 64, the G light detection element 72 that detects the G light reflected by the beam splitter 58, and the R light reflected by the beam splitter 60. It includes an R light detection element 74 and an IR light detection element 76 that detects the IR light reflected by the beam splitter 62.
  • the light incident on the SLO unit 18 via the wide-angle optical system 30 and the optical scanner 22 is reflected by the beam splitter 64 and received by the B light detection element 70.
  • the beam splitter 58 In the case of G light, it is reflected by the beam splitter 58 and received by the G light detection element 72.
  • the incident light passes through the beam splitter 58, is reflected by the beam splitter 60, and is received by the R photodetection element 74.
  • the incident light passes through the beam splitters 58 and 60, is reflected by the beam splitter 62, and is received by the IR photodetection element 76.
  • the image processor 17 operating under the control of the CPU 16A displays a UWF-SLO image using the signals detected by the B photodetection element 70, the G photodetection element 72, the R photodetection element 74, and the IR photodetection element 76. Generate.
  • control device 16 controls the light sources 40, 42, and 44 so as to emit light at the same time.
  • a G color fundus image, an R color fundus image, and a B color fundus image in which the positions correspond to each other can be obtained.
  • An RGB color fundus image can be obtained from a G color fundus image, an R color fundus image, and a B color fundus image.
  • the control device 16 controls the light sources 42 and 44 so as to emit light at the same time, and the fundus of the eye to be inspected 12 is simultaneously photographed by the G light and the R light, so that the G color fundus image and the R color corresponding to each other at each position are photographed at the same time.
  • a fundus image is obtained.
  • An RG color fundus image can be obtained from the G color fundus image and the R color fundus image.
  • the viewing angle (FOV: Field of View) of the fundus is set to an ultra-wide angle, and the region beyond the equator from the posterior pole of the fundus of the eye 12 to be inspected can be photographed.
  • the OCT system is realized by the control device 16, the OCT unit 20, and the photographing optical system 19 shown in FIG. Since the OCT system includes a wide-angle optical system 30, it enables OCT imaging of the peripheral portion of the fundus in the same manner as the acquisition of the SLO fundus image described above. That is, the wide-angle optical system 30 in which the viewing angle (FOV) of the fundus is an ultra-wide angle enables OCT imaging of a region beyond the equator portion 178 from the posterior pole portion of the fundus of the eye 12 to be inspected. OCT data of structures existing around the fundus such as vortex veins can be acquired, and tomographic images of vortex veins and 3D structures of vortex veins can be obtained by image processing of OCT data.
  • FOV viewing angle
  • the OCT unit 20 includes a light source 20A, a sensor (detection element) 20B, a first optical coupler 20C, a reference optical system 20D, a collimating lens 20E, and a second optical coupler 20F.
  • the light emitted from the light source 20A is branched by the first optical coupler 20C.
  • One of the branched lights is converted into parallel light by the collimated lens 20E as measurement light, and then incident on the photographing optical system 19.
  • the measurement light is applied to the fundus through the wide-angle optical system 30 and the pupil 27.
  • the measurement light reflected by the fundus is incident on the OCT unit 20 via the wide-angle optical system 30 and incident on the second optical coupler 20F via the collimating lens 20E and the first optical coupler 20C.
  • the other light emitted from the light source 20A and branched by the first optical coupler 20C is incident on the reference optical system 20D as reference light, and is incident on the second optical coupler 20F via the reference optical system 20D. do.
  • the image processor 17 that operates under the control of the image processing unit 206 performs signal processing such as Fourier transform on the detection signal detected by the sensor 20B, and generates OCT data. It is also possible to generate an OCT image such as a tomographic image or an en-face image by the image processor 17 based on the OCT data.
  • the OCT unit 20 can scan a predetermined range (for example, a rectangular range of 6 mm ⁇ 6 mm) in one OCT imaging.
  • the predetermined range is not limited to 6 mm ⁇ 6 mm, but may be a square range of 12 mm ⁇ 12 mm or 23 mm ⁇ 23 mm, or a rectangular range such as 14 mm ⁇ 9 mm or 6 mm ⁇ 3.5 mm, and may be an arbitrary rectangular range. can.
  • the diameter may be in the range of a circular diameter such as 6 mm, 12 mm, and 23 mm.
  • the ophthalmologic apparatus 110 can scan a region 12A having an internal irradiation angle of 200 °.
  • the ophthalmic apparatus 110 can generate OCT data by the OCT imaging. Therefore, the ophthalmic apparatus 110 is an OCT image, a tomographic image of the fundus including a vortex vein (B-scan image), an OCT volume data including the vortex vein, and an en-face image (OCT) which is a cross section of the OCT volume data. It is possible to generate a front image) generated from volume data. Needless to say, the OCT image includes an OCT image of the central part of the fundus (the posterior pole of the eyeball where the macula and the optic disc are present).
  • the OCT data (or the image data of the OCT image) is sent from the ophthalmologic device 110 to the server 140 via the communication interface 16F, and is stored in the storage device 254.
  • the light source 20A exemplifies a wavelength sweep type SS-OCT (Swept-Source OCT), but there are various types such as SD-OCT (Spectral-Domain OCT) and TD-OCT (Time-Domain OCT). It may be an OCT system of various types.
  • the server 140 includes a computer main body 252.
  • the computer body 252 has a CPU 262, a RAM 266, a ROM 264, and an input / output (I / O) port 268.
  • a storage device 254, a display 256, a mouse 255M, a keyboard 255K, and a communication interface (I / F) 258 are connected to the input / output (I / O) port 268.
  • the storage device 254 is composed of, for example, a non-volatile memory.
  • the input / output (I / O) port 268 is connected to the network 130 via the communication interface (I / F) 258.
  • the server 140 can communicate with the ophthalmic apparatus 110 and the viewer 150.
  • the image processing program shown in FIG. 6 is stored in the ROM 264 or the storage device 254.
  • the ROM 264 or the storage device 254 is an example of a "memory" of the technique of the present disclosure.
  • the CPU 262 is an example of a "processor” of the technique of the present disclosure.
  • the image processing program is an example of a "program" of the technique of the present disclosure.
  • the server 140 stores each data received from the ophthalmic device 110 in the storage device 254.
  • the image processing program includes a display control function, an image processing function, and a processing function.
  • the CPU 262 executes an image processing program having each of these functions, the CPU 262 functions as a display control unit 204, an image processing unit 206, and a processing unit 208.
  • the image processing unit 206 acquires the OCT volume data 400, which is the OCT data, from the storage device 254.
  • the OCT volume data 400 is OCT volume data having a predetermined area including the vortex vein VV.
  • the OCT volume data 400 is obtained by OCT imaging one of a plurality of vortex vein VVs with an ophthalmologic apparatus 110.
  • the predetermined area is, for example, a rectangular area of 6 mm ⁇ 6 mm.
  • N faces having different depths from the first plane f401 to the Nth plane f40N are set for the OCT volume data 400.
  • the OCT volume data 400 may be obtained by OCT imaging of each of the plurality of vortex vein VVs existing in the eye to be inspected by the ophthalmologic apparatus 110.
  • the position of each of the plurality of vortex veins is estimated based on the position where the choroidal blood vessels are gathered, for example, by extracting the choroidal blood vessels from the choroidal blood vessel image, estimating the movement direction (blood vessel running direction) of each choroidal blood vessel, and estimating the position. Will be done.
  • the choroidal vessel image is an R-UWF-SLO image (R color fundus image) taken with red light (laser light with a wavelength of 630 to 660 nm) and a G- taken with green light (laser light with a wavelength of 500 to 550 nm). It is generated by image processing the image data of the UWF-SLO image (G color fundus image).
  • a choroidal blood vessel image is generated by extracting retinal blood vessels from the G-colored fundus image, removing the retinal blood vessels from the R-colored fundus image, and performing image processing for enhancing the choroidal blood vessels.
  • step 602 the image processing unit 206 sets the parameter n to 1.
  • the parameter n is a parameter indicating the number of en-face images.
  • the image processing unit 206 analyzes the OCT volume data 400 and is below a predetermined number of pixels from the retinal pigment epithelial cell layer (hereinafter referred to as RPE layer) in the OCT volume data 400, for example, 10 pixels.
  • the lower surface is set as the first surface.
  • the image processing unit 206 first specifies the RPE layer as a reference plane.
  • the RPE layer is specified by performing a predetermined segmentation process on the OCT volume data 400. Further, the image processing unit 206 may specify the RPE layer as the most bright layer in the OCT volume data 400 as the RPE layer.
  • the image processing unit 206 sets the surface 10 pixels below the RPE layer as the first surface. This is because the region deeper than the RPE layer (the region farther than the RPE layer when viewed from the center of the eyeball) is the choroidal region, and an en-face image of the region where the choroidal blood vessels are present is generated.
  • the surface 10 pixels below the RPE layer is not limited to being set as the first surface, and for example, the surface 10 pixels below the Bruch film existing immediately below the RPE layer may be set as the first surface.
  • Bruch's membrane is also identified by subjecting the OCT volume data 400 to another predetermined segmentation process different from that for the RPE layer.
  • it may be 10 pixels below in the direction of the A-scan when the OCT volume data is generated.
  • the surface 10 pixels below the RPE layer or Bruch film is not limited to being specified as the first surface, and may be set with an arbitrary number of pixels. Further, it may be defined by a length such as millimail or nanometer instead of the definition by the number of pixels. Further, as the reference plane, a spherical surface maintained at a certain distance from the pupil or the center of the eyeball may be defined as the reference plane.
  • the first aspect is an example of "other aspects" of the technique of the present disclosure.
  • the RPE layer or Bruch film is an example of the "reference plane" of the technique of the present disclosure.
  • the image processing unit 206 generates a first en-face image 7a (see FIG. 7) corresponding to the set first surface.
  • the en-face image 7a may be generated from the pixel values of the pixels existing on the first surface, or the pixel group in the shallow direction and the pixel group in the deep direction including the first surface are extracted from the OCT volume data 400 and these are extracted.
  • the pixel value may be obtained as the average value or the median value of the brightness values of the pixel group of. Image processing such as noise removal may be used when obtaining the pixel value.
  • the first en-face image corresponding to the generated first surface is stored in the RAM 266 by the processing unit 208. It should be noted that the en-face image corresponding to the RPE layer or Bruch film which is the reference plane may be generated.
  • the image processing unit 206 determines whether or not the parameter n is the preset maximum number N.
  • the image processing unit 206 generates an en-face image corresponding to the above.
  • the maximum number N can be arbitrarily set by the user depending on the thickness of the retina near the vortex vein and the length of the interval between the en-face images.
  • step 608 If it is determined that the parameter n is the maximum number N (step 608: YES), the image processing proceeds to step 610.
  • the processing unit 208 stores the generated N-face images in the RAM 266 or the storage device 254, and ends the processing.
  • step 608 the image processing proceeds to step 612.
  • the image processing unit 206 sets the parameter n to n + 1.
  • step 614 the image processing returns to step 606.
  • the image processing unit 206 generates the nth (eg, 2) en-face image 7b (see FIG. 7) corresponding to the set nth (eg, 2) plane.
  • the generated second en-face image 7b is stored in the RAM 266.
  • the image processing unit 206 repeats the loop of step 606, step 608, step 612, and step 614 until the parameter n reaches the maximum number N.
  • the Nth en-face image 7d from the first en-face image 7a shown in FIG. 7 is generated.
  • These plurality of en-face images are a group of en-face images generated for each predetermined number of pixels in the direction in which the choroid is deepened.
  • the vortex vein For example, comparing the first en-face image 7a, the second en-face image 7b, ... the N-1 en-face image 7c, and the Nth en-face image 7d in order, the vortex vein. It is possible to specify the traveling direction of the vehicle, particularly the traveling direction in the depth direction. In the example shown in FIG. 7, it can be seen that the vortex vein runs in the upper left direction of the paper.
  • the vortex vein is obtained by image processing the en-face image group as shown in FIG. 7, or by artificial intelligence using a trained model in which the traveling direction of the en-face image group and the vortex vein is machine-learned. You may try to identify the vascular running of.
  • the blood vessel running of the vortex vein may be specified by artificial intelligence using the OCT volume data and the learned model in which the running direction of the vortex vein is machine-learned. Based on the specified blood vessel running, information indicating the direction of blood vessel running may be generated, and the information indicating the direction of blood vessel running may be superimposed and displayed on the first to Nth en-face images. Further, it may be superimposed and displayed on the OCT volume data 400.
  • the information indicating the direction of blood vessel travel includes vector information, mark information such as an arrow indicating the direction, and information for blinking so as to flow in the traveling direction. Further, the region of the pixel indicating the blood vessel in each en-face image of the first to Nth en-face images is obtained, and the position of the center of gravity of the blood vessel is calculated. Then, a vector connecting the positions of the centers of gravity calculated from the first en-face image to the Nth en-face image (from the upper en-face image to the lower en-face image) is obtained from the first en-face image. It may be displayed in all the en-face images of the third Nen-face image. Further, the position of the center of gravity in each en-face image may be displayed in the OCT volume data.
  • a first surface below a predetermined number of pixels (10 pixels in the above example) is set from the RPE layer, and then surfaces are set from the second surface to the eighth surface in order, for a total of eight images.
  • An en-face image corresponding to each surface is generated.
  • a surface 40 pixels below the RPE layer may be set as a reference surface, and other surfaces may be set above and below every 10 pixels from the reference surface to generate an en-face image of each surface. good.
  • the display screen is generated by the display control unit 204 of the server 140 based on the user's instruction, and is output as an image signal to the viewer 150 by the processing unit 208.
  • the viewer 150 displays the display screen on a monitor or the like based on the image signal.
  • FIG. 8 shows the first display screen 500A. As shown in FIG. 8, the first display screen 500A has an information area 502 and an image display area 504A.
  • the information area 502 has a patient ID display field 512, a patient name display field 514, an age display field 516, a visual acuity display field 515, a right eye / left eye display field 520, and an axial length display field 522.
  • the viewer 150 displays each information based on the information received from the server 140.
  • the image display area 504A is an area for displaying an image to be inspected or the like.
  • the image display area 504A is provided with the following display fields, specifically, UWF fundus image display field 542, choroidal blood vessel image display field 544, B-scan image display field 546, and en-face image display.
  • the comment field 530 is a remarks column in which the result observed by the ophthalmologist who is the user or the diagnosis result can be arbitrarily input.
  • the UWF-SLO fundus image SL taken by SLO with the ophthalmologic apparatus 110 is displayed on the fundus of the eye to be inspected.
  • the UWF-SLO fundus image SL is superposed with a line SLV showing a vertical watershed and a line SLH showing a horizontal watershed.
  • the rectangular region OCV1 indicating the acquisition position of the OCT volume data 400 around the vortex vein is superimposed and displayed.
  • the line OCB1 indicating the cross-sectional position of the B-scan image OCB displayed in the B-scan image display field 526 is superimposed and displayed.
  • the choroidal blood vessel image CV obtained by image processing the UWF-SLO fundus image SL1 is displayed.
  • the choroidal blood vessel image CV is superposed with a line CVV showing a vertical watershed and a line CVH showing a horizontal watershed.
  • the rectangular region OCV2 indicating the acquisition position of the OCT volume data 400 around the vortex vein is superimposed and displayed.
  • the line OCB2 indicating the cross-sectional position of the B-scan image OCB displayed in the B-scan image display field 526 is superimposed and displayed.
  • the vortex vein to be displayed in the B-scan image display field 546 and the en-face image display field 548 may be selected.
  • the rectangular region OCV2 is superimposed and displayed in the range of the OCT volume data of the selected vortex vein.
  • the B-scan image OCB generated from the OCT volume data 400 is displayed in the B-scan image display field 546.
  • the curve OCL superimposed on the B-scan image OCB is a curve indicating the position of the acquisition surface of the en-face image in the depth direction.
  • the curve OCL is the superimposed position of the curve OCL. Also changes.
  • the moving image generated from the en-face image group is displayed in the en-face image display field 548, the superimposed position of the curve OCL changes with the reproduction of the moving image.
  • en-face image display field 548 a plurality of en-face images from the first en-face image 7a to the Nth en-face image 7d are displayed in order like a slide show, or are displayed as a moving image. A specific en-face image can be displayed or displayed as a still image.
  • the superimposed position of the curve OCL on the B-scan image OCB displayed in the B-scan image display field 546 changes according to the en-face image displayed in the en-face image display field 548.
  • the line OCB3 indicating the cross-sectional position of the B-scan image OCB may be superimposed and displayed on the en-face image displayed in the n-face image display field 548.
  • the line OCB3 may be selected to be superimposed / not superimposed according to the user's instruction. By horizontally moving this line OCB3, the B-scan image OCB to be displayed in the B-scan image display field 546 can be changed.
  • the line OCB3 can be set not only in the horizontal direction but also in any direction such as a vertical direction. Further, the above-mentioned information indicating the direction of blood vessel traveling may be superimposed and displayed on the first to Nth en-face images.
  • FIG. 9 shows a second display screen 500B. Since the second display screen 500B has the same fields as the fields of the first display screen 500A, the same fields are designated by the same reference numerals, the description thereof will be omitted, and different parts will be described.
  • the display screen 500B has an en-face image display field 560 for displaying N sheets of en-face images side by side.
  • the rectangular frame 560f shown by the dotted line is superimposed and displayed on the second en-face image 7b.
  • the position of the surface of the selected en-face image is displayed as a curve OCL in the B-scan image OCB of the B-scan image display field 546.
  • a three-dimensional image obtained by image-processing the OCT volume data 400 may be displayed on the first display screen 500A or the second display screen 500B.
  • the OCT volume data 400 the blood vessel area is black, and the space (interstitium) other than the blood vessel is binarized to extract the blood vessel area, and the blood vessel area is displayed by executing the 3D volume rendering process. May be.
  • it may be displayed after being converted into a file format having a polygonal mesh structure such as an STL (Standard Triangulated Language) file format.
  • STL Standard Triangulated Language
  • each of the plurality of en-face images generated from the OCT volume data 400 is binarized to generate a binarized en-face image. Then, the binarized OCT volume data 400 may be obtained by resynthesizing the binarized en-face image.
  • image processing for visualizing the vortex vein and the choroidal blood vessels connected to the vortex vein may be performed on the OCT volume data 400 to display a three-dimensional image emphasizing the blood vessels around the vortex vein.
  • the retina and choroid around the fundus may be thinner than the central part, and it is difficult to observe the three-dimensional structure of the vortex vein when the actual size is displayed. Therefore, the OCT volume data 400 may be subjected to image processing in which the OCT volume data 400 is stretched only in the depth direction (the display unit is changed only in the depth direction) to display a three-dimensional image.
  • the state of the vortex vein traveling from the choroid to the sclera is visualized by generating an en-face image from the OCT volume data in which the region including the vortex vein is photographed.
  • the vortex vein changes significantly in the depth direction compared to the lateral change, but it was difficult to understand just by visualizing the choroidal vascular structure in the OCT volume data by using the en-face image generated from the OCT volume data.
  • the morphological change in the depth direction of the vortex vein can be visualized and quantified with high resolution.
  • the vortex vein is visualized without using the angiography process by the Argolim of OCT angiography.
  • the image processing (FIG. 5) is executed by the server 140, but the technique of the present disclosure is not limited to this, and additional additions provided in the ophthalmic apparatus 110, the viewer 150, or the network 130.
  • the image processing device may execute.
  • each component may be present alone or in combination of two or more as long as there is no contradiction.
  • image processing is realized by a software configuration using a computer
  • the technique of the present disclosure is not limited to this.
  • the image processing may be executed only by the hardware configuration such as FPGA (Field-Programmable Gate Array) or ASIC (Application Specific Integrated Circuit). Some of the image processing may be performed by the software configuration and the rest may be performed by the hardware configuration.
  • the technique of the present disclosure includes the following techniques because it includes cases where image processing is realized and cases where image processing is not realized by a software configuration using a computer.
  • the acquisition unit that acquires OCT volume data for the area including the vortex vein, A generation unit that generates a plurality of en-face images corresponding to a plurality of surfaces having different depths based on the OCT volume data.
  • An image processing device comprising.
  • Image processing methods including.
  • the image processing unit 206 is an example of the "acquisition unit” and the “generation unit” of the technique of the present disclosure.
  • a computer program product for image processing comprises a computer-readable storage medium that is not itself a temporary signal.
  • a program is stored in the computer-readable storage medium.
  • the program On the computer Steps to acquire OCT volume data for areas containing vortex veins, Based on the OCT volume data, a step of generating a plurality of en-face images corresponding to a plurality of surfaces having different depths, and To execute, Computer program product.
  • the server 140 is an example of the "computer program product" of the technique of the present disclosure.

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Ophthalmology & Optometry (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Eye Examination Apparatus (AREA)

Abstract

被検眼のボリュームデータを解析する。　プロセッサが行う画像処理方法は、渦静脈を含む領域のＯＣＴボリュームデータを取得するステップと、前記ＯＣＴボリュームデータに基づいて、深さが異なる複数の面に対応した複数のｅｎ－ｆａｃｅ画像を生成するステップと、を含む。

Description

画像処理方法、画像処理装置、及びプログラム
　本開示の技術は、画像処理方法、画像処理装置、及びプログラムに関する。
　米国特許第１０２３８２８１号公報には、光干渉断層計を用いて被検眼のボリュームデータを生成する技術に関する。従来、被検眼のボリュームデータを解析することが望まれている。
　本開示の技術の第１の態様は、プロセッサが行う画像処理方法であって、渦静脈を含む領域のＯＣＴボリュームデータを取得するステップと、前記ＯＣＴボリュームデータに基づいて、深さが異なる複数の面に対応した複数のｅｎ－ｆａｃｅ画像を生成するステップと、を含む。
　本開示の技術の第２の態様の画像処理装置は、メモリと、前記メモリに接続するプロセッサとを備え、前記プロセッサは、渦静脈を含む領域のＯＣＴボリュームデータを取得するステップと、前記ＯＣＴボリュームデータに基づいて、深さが異なる複数の面に対応した複数のｅｎ－ｆａｃｅ画像を生成するステップと、を実行する。
　本開示の技術の第３の態様のプログラムは、コンピュータに、渦静脈を含む領域のＯＣＴボリュームデータを取得するステップと、前記ＯＣＴボリュームデータに基づいて、深さが異なる複数の面に対応した複数のｅｎ－ｆａｃｅ画像を生成するステップと、を実行させる。
本実施形態の眼科システムの概略構成図である。 本実施形態の眼科装置の概略構成図である。 サーバの概略構成図である。 サーバのＣＰＵにおいて、画像処理プログラムによって実現される機能の説明図である。 サーバによる画像処理を示したフローチャートである。 ＯＣＴボリュームデータとｅｎ－ｆａｃｅ画像の関係を示す図である。 ＯＣＴボリュームデータから得られる複数枚のｅｎ－ｆａｃｅ画像を示す図である。 ｅｎ－ｆａｃｅ画像を用いた表示画面の第１の例である。 ｅｎ－ｆａｃｅ画像を用いた表示画面の第２の例である。 眼球と渦静脈の位置との関係を示す模式図である。
　以下、本発明の実施形態に係る眼科システム１００について図面を参照して説明する。図１には、眼科システム１００の概略構成が示されている。図１に示すように、眼科システム１００は、眼科装置１１０と、サーバ装置（以下、「サーバ」という）１４０と、表示装置（以下、「ビューワ」という）１５０と、を備えている。眼科装置１１０は、眼底画像を取得する。サーバ１４０は、眼科装置１１０によって複数の患者の眼底が撮影されることにより得られた複数の眼底画像と、図示しない眼軸長測定装置により測定された眼軸長とを、患者ＩＤに対応して記憶する。ビューワ１５０は、サーバ１４０により取得した眼底画像や解析結果を表示する。
　サーバ１４０は、本開示の技術の「画像処理装置」の一例である。
　眼科装置１１０、サーバ１４０、ビューワ１５０は、ネットワーク１３０を介して、相互に接続されている。ネットワーク１３０は、ＬＡＮ、ＷＡＮ、インターネットや広域イーサ網等の任意のネットワークである。例えば、眼科システム１００が１つの病院に構築される場合には、ネットワーク１３０にＬＡＮを採用することができる。
　ビューワ１５０は、クライアントサーバシステムにおけるクライアントであり、ネットワークを介して複数台が接続される。また、サーバ１４０も、システムの冗長性を担保するために、ネットワークを介して複数台が接続されていてもよい。又は、眼科装置１１０が画像処理機能及びビューワ１５０の画像閲覧機能を備えるのであれば、眼科装置１１０がスタンドアロン状態で、眼底画像の取得、画像処理及び画像閲覧が可能となる。この場合、ビューワ１５０は省略してもよい。また、サーバ１４０がビューワ１５０の画像閲覧機能を備えるのであれば、眼科装置１１０とサーバ１４０との構成で、眼底画像の取得、画像処理及び画像閲覧が可能となる。この場合、ビューワ１５０は省略してもよい。
　なお、他の眼科装置（視野測定、眼圧測定などの検査機器）やＡＩ（Ａｒｔｉｆｉｃｉａｌ　Ｉｎｔｅｌｌｉｇｅｎｃｅ）を用いた画像解析を行う診断支援装置がネットワーク１３０を介して、眼科装置１１０、サーバ１４０、及びビューワ１５０に接続されていてもよい。
　次に、図２を参照して、眼科装置１１０の構成を説明する。
　説明の便宜上、走査型レーザ検眼鏡（Ｓｃａｎｎｉｎｇ　Ｌａｓｅｒ　Ｏｐｈｔｈａｌｍｏｓｃｏｐｅ）を「ＳＬＯ」と称する。また、光干渉断層計（Ｏｐｔｉｃａｌ　Ｃｏｈｅｒｅｎｃｅ　Ｔｏｍｏｇｒａｐｈｙ）を「ＯＣＴ」と称する。
　なお、眼科装置１１０が水平面に設置された場合の水平方向を「Ｘ方向」、水平面に対する垂直方向を「Ｙ方向」とし、被検眼１２の前眼部の瞳孔の中心と眼球の中心とを結ぶ方向を「Ｚ方向」とする。従って、Ｘ方向、Ｙ方向、およびＺ方向は互いに垂直である。
　眼科装置１１０は、撮影装置１４および制御装置１６を含む。撮影装置１４は、ＳＬＯユニット１８およびＯＣＴユニット２０を備えており、被検眼１２の眼底の眼底画像を取得する。以下、ＳＬＯユニット１８により取得された二次元眼底画像をＳＬＯ画像と称する。また、ＯＣＴユニット２０により取得されたＯＣＴデータに基づいて作成された網膜の断層画像や正面画像（ｅｎ－ｆａｃｅ画像）などをＯＣＴ画像と称する。
　制御装置１６は、ＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ（中央処理装置））１６Ａ、ＲＡＭ（Ｒａｎｄｏｍ　Ａｃｃｅｓｓ　Ｍｅｍｏｒｙ）１６Ｂ、ＲＯＭ（Ｒｅａｄ－Ｏｎｌｙ　ｍｅｍｏｒｙ）１６Ｃ、および入出力（Ｉ／Ｏ）ポート１６Ｄを有するコンピュータを備えている。
　制御装置１６は、Ｉ／Ｏポート１６Ｄを介してＣＰＵ１６Ａに接続された入力／表示装置１６Ｅを備えている。入力／表示装置１６Ｅは、被検眼１２の画像を表示したり、ユーザから各種指示を受け付けたりするグラフィックユーザインターフェースを有する。グラフィックユーザインターフェースとしては、タッチパネル・ディスプレイが挙げられる。
　また、制御装置１６は、Ｉ／Ｏポート１６Ｄに接続された画像処理器１７を備えている。画像処理器１７は、撮影装置１４によって得られたデータに基づき被検眼１２の画像を生成する。画像処理器１７を省略し、ＣＰＵ１６Ａが、撮影装置１４によって得られたデータに基づき被検眼１２の画像を生成してもよい。なお、制御装置１６は、通信インターフェース１６Ｆを介してネットワーク１３０に接続される。
　上記のように、図２では、眼科装置１１０の制御装置１６が入力／表示装置１６Ｅを備えているが、本開示の技術はこれに限定されない。例えば、眼科装置１１０の制御装置１６は入力／表示装置１６Ｅを備えず、眼科装置１１０とは物理的に独立した別個の入力／表示装置を備えるようにしてもよい。この場合、当該表示装置は、制御装置１６のＣＰＵ１６Ａの表示制御部２０４の制御下で動作する画像処理プロセッサユニットを備える。画像処理プロセッサユニットが、表示制御部２０４が出力指示した画像信号に基づいて、ＳＬＯ画像等を表示するようにしてもよい。
　撮影装置１４は、制御装置１６のＣＰＵ１６Ａの制御下で作動する。撮影装置１４は、ＳＬＯユニット１８、撮影光学系１９、およびＯＣＴユニット２０を含む。撮影光学系１９は、光学スキャナ２２、および広角光学系３０を含む。
　光学スキャナ２２は、ＳＬＯユニット１８から射出された光をＸ方向、およびＹ方向に２次元走査する。光学スキャナ２２は、光束を偏向できる光学素子であればよく、例えば、ポリゴンミラーや、ガルバノミラー等を用いることができる。また、それらの組み合わせであってもよい。
　広角光学系３０は、ＳＬＯユニット１８からの光とＯＣＴユニット２０からの光とを合成する。
　なお、広角光学系３０は、楕円鏡などの凹面ミラーを用いた反射光学系や、広角レンズなどを用いた屈折光学系、あるいは、凹面ミラーやレンズを組み合わせた反射屈折光学系でもよい。楕円鏡や広角レンズなどを用いた広角光学系を用いることにより、眼底中心部だけでなく眼底周辺部の網膜を撮影することが可能となる。
　楕円鏡を含むシステムを用いる場合には、国際公開ＷＯ２０１６／１０３４８４あるいは国際公開ＷＯ２０１６／１０３４８９に記載された楕円鏡を用いたシステムを用いる構成でもよい。国際公開ＷＯ２０１６／１０３４８４の開示および国際公開ＷＯ２０１６／１０３４８９の開示の各々は、その全体が参照により本明細書に取り込まれる。
　広角光学系３０によって、眼底において広い視野（ＦＯＶ：Ｆｉｅｌｄ　ｏｆ　Ｖｉｅｗ）１２Ａでの観察が実現される。ＦＯＶ１２Ａは、撮影装置１４によって撮影可能な範囲を示している。ＦＯＶ１２Ａは、視野角として表現され得る。視野角は、本実施形態において、内部照射角と外部照射角とで規定され得る。外部照射角とは、眼科装置１１０から被検眼１２へ照射される光束の照射角を、瞳孔２７を基準として規定した照射角である。また、内部照射角とは、眼底Ｆへ照射される光束の照射角を、眼球中心Ｏを基準として規定した照射角である。外部照射角と内部照射角とは、対応関係にある。例えば、外部照射角が１２０度の場合、内部照射角は約１６０度に相当する。本実施形態では、内部照射角は２００度としている。
　ここで、内部照射角で１６０度以上の撮影画角で撮影されて得られたＳＬＯ眼底画像をＵＷＦ－ＳＬＯ眼底画像と称する。なお、ＵＷＦとは、ＵｌｔｒａＷｉｄｅ　Ｆｉｅｌｄ（超広角）の略称を指す。眼底の視野角（ＦＯＶ）を超広角な角度とした広角光学系３０により、被検眼１２の眼底の後極部から赤道部を超える領域を撮影することができ、渦静脈などの眼底周辺部に存在する構造物を撮影できる。
　眼科装置１１０は、被検眼１２の眼球中心Ｏを基準位置として内部照射角が２００°の領域１２Ａを撮影することができる。なお、２００°の内部照射角は、被検眼１２の眼球の瞳孔を基準とした外部照射角では１１０°である。つまり、広角光学系３０は外部照射角１１０°の画角で瞳からレーザ光を照射させ、内部照射角で２００°の眼底領域を撮影する。
　ＳＬＯシステムは、図２に示す制御装置１６、ＳＬＯユニット１８、および撮影光学系１９によって実現される。ＳＬＯシステムは、広角光学系３０を備えるため、広いＦＯＶ１２Ａでの眼底撮影を可能とする。
　ＳＬＯユニット１８は、Ｂ光（青色光）の光源４０、Ｇ光（緑色光）の光源４２、Ｒ光（赤色光）の光源４４、およびＩＲ光（赤外線（例えば、近赤外光））の光源４６と、光源４０、４２、４４、４６からの光を、反射又は透過して１つの光路に導く光学系４８、５０、５２、５４、５６とを備えている。光学系４８、５６は、ミラーであり、光学系５０、５２、５４は、ビームスプリッタ―である。Ｂ光は、光学系４８で反射し、光学系５０を透過し、光学系５４で反射し、Ｇ光は、光学系５０、５４で反射し、Ｒ光は、光学系５２、５４を透過し、ＩＲ光は、光学系５２、５６で反射して、それぞれ１つの光路に導かれる。
　ＳＬＯユニット１８は、Ｒ光およびＧ光を発するモードと、赤外線を発するモードなど、波長の異なるレーザ光を発する光源あるいは発光させる光源の組合せを切り替え可能に構成されている。図２に示す例では、Ｂ光の光源４０、Ｇ光の光源４２、Ｒ光の光源４４、およびＩＲ光の光源４６の４つの光源を備えるが、本開示の技術は、これに限定されない。例えば、ＳＬＯユニット１８は、更に、白色光の光源を更に備え、Ｇ光、Ｒ光、およびＢ光を発するモードや、白色光のみを発するモード等の種々のモードで光を発するようにしてもよい。
　ＳＬＯユニット１８から撮影光学系１９に入射された光は、光学スキャナ２２によってＸ方向およびＹ方向に走査される。走査光は広角光学系３０および瞳孔２７を経由して、眼底に照射される。眼底により反射された反射光は、広角光学系３０および光学スキャナ２２を経由してＳＬＯユニット１８へ入射される。
　ＳＬＯユニット１８は、被検眼１２の後眼部（眼底）からの光の内、Ｂ光を反射し且つＢ光以外を透過するビームスプリッタ６４、ビームスプリッタ６４を透過した光の内、Ｇ光を反射し且つＧ光以外を透過するビームスプリッタ５８を備えている。ＳＬＯユニット１８は、ビームスプリッタ５８を透過した光の内、Ｒ光を反射し且つＲ光以外を透過するビームスプリッタ６０を備えている。ＳＬＯユニット１８は、ビームスプリッタ６０を透過した光の内、ＩＲ光を反射するビームスプリッタ６２を備えている。ＳＬＯユニット１８は、ビームスプリッタ６４により反射したＢ光を検出するＢ光検出素子７０、ビームスプリッタ５８により反射したＧ光を検出するＧ光検出素子７２、ビームスプリッタ６０により反射したＲ光を検出するＲ光検出素子７４、およびビームスプリッタ６２により反射したＩＲ光を検出するＩＲ光検出素子７６を備えている。
　広角光学系３０および光学スキャナ２２を経由してＳＬＯユニット１８へ入射された光（眼底により反射された反射光）は、Ｂ光の場合、ビームスプリッタ６４で反射してＢ光検出素子７０により受光され、Ｇ光の場合、ビームスプリッタ５８で反射してＧ光検出素子７２により受光される。上記入射された光は、Ｒ光の場合、ビームスプリッタ５８を透過し、ビームスプリッタ６０で反射してＲ光検出素子７４により受光される。上記入射された光は、ＩＲ光の場合、ビームスプリッタ５８、６０を透過し、ビームスプリッタ６２で反射してＩＲ光検出素子７６により受光される。ＣＰＵ１６Ａの制御下で動作する画像処理器１７は、Ｂ光検出素子７０、Ｇ光検出素子７２、Ｒ光検出素子７４、およびＩＲ光検出素子７６で検出された信号を用いてＵＷＦ－ＳＬＯ画像を生成する。
　また、制御装置１６が、同時に発光するように光源４０、４２、４４を制御する。Ｂ光、Ｇ光およびＲ光で同時に被検眼１２の眼底が撮影されることにより、各位置が互いに対応するＧ色眼底画像、Ｒ色眼底画像、およびＢ色眼底画像が得られる。Ｇ色眼底画像、Ｒ色眼底画像、およびＢ色眼底画像からＲＧＢカラー眼底画像が得られる。制御装置１６が、同時に発光するように光源４２、４４を制御し、Ｇ光およびＲ光で同時に被検眼１２の眼底が撮影されることにより、各位置が互いに対応するＧ色眼底画像およびＲ色眼底画像が得られる。Ｇ色眼底画像およびＲ色眼底画像からＲＧカラー眼底画像が得られる。
　広角光学系３０により、眼底の視野角（ＦＯＶ：Ｆｉｅｌｄ　ｏｆ　Ｖｉｅｗ）を超広角な角度とし、被検眼１２の眼底の後極部から赤道部を超える領域を撮影することができる。
　ＯＣＴシステムは、図２に示す制御装置１６、ＯＣＴユニット２０、および撮影光学系１９によって実現される。ＯＣＴシステムは、広角光学系３０を備えるため、上述したＳＬＯ眼底画像の撮影と同様に、眼底周辺部のＯＣＴ撮影を可能とする。つまり、眼底の視野角（ＦＯＶ）を超広角な角度とした広角光学系３０により、被検眼１２の眼底の後極部から赤道部１７８を超える領域のＯＣＴ撮影を行うことができる。渦静脈などの眼底周辺部に存在する構造物のＯＣＴデータを取得でき、渦静脈の断層像や、ＯＣＴデータを画像処理することにより渦静脈の３Ｄ構造を得ることができる。
　ＯＣＴユニット２０は、光源２０Ａ、センサ（検出素子）２０Ｂ、第一の光カプラ２０Ｃ、参照光学系２０Ｄ、コリメートレンズ２０Ｅ、および第２の光カプラ２０Ｆを含む。
　光源２０Ａから射出された光は、第一の光カプラ２０Ｃで分岐される。分岐された一方の光は、測定光として、コリメートレンズ２０Ｅで平行光にされた後、撮影光学系１９に入射される。測定光は広角光学系３０および瞳孔２７を経由して、眼底に照射される。眼底により反射された測定光は、および広角光学系３０を経由してＯＣＴユニット２０へ入射され、コリメートレンズ２０Ｅおよび第一の光カプラ２０Ｃを介して、第２の光カプラ２０Ｆに入射する。
　光源２０Ａから射出され、第一の光カプラ２０Ｃで分岐された他方の光は、参照光として、参照光学系２０Ｄへ入射され、参照光学系２０Ｄを経由して、第２の光カプラ２０Ｆに入射する。
　第２の光カプラ２０Ｆに入射されたこれらの光、即ち、眼底で反射された測定光と、参照光とは、第２の光カプラ２０Ｆで干渉されて干渉光を生成する。干渉光はセンサ２０Ｂで受光される。画像処理部２０６の制御下で動作する画像処理器１７は、センサ２０Ｂで検出された検出信号に対してフーリエ変換などの信号処理を行い、ＯＣＴデータを生成する。当該ＯＣＴデータに基づいて断層画像やｅｎ－ｆａｃｅ画像などのＯＣＴ画像を画像処理器１７で生成することも可能である。
　ここで、ＯＣＴユニット２０は、所定範囲（例えば６ｍｍ×６ｍｍの矩形範囲）を一回のＯＣＴ撮影で走査することができる。当該所定範囲は６ｍｍ×６ｍｍに限らず、１２ｍｍ×１２ｍｍや２３ｍｍ×２３ｍｍの正方形の範囲でもよいし、１４ｍｍ×９ｍｍ、６ｍｍ×３．５ｍｍなど長方形の範囲でもよく、任意の矩形範囲とすることができる。また、直径６ｍｍ、１２ｍｍ、２３ｍｍなどの円径の範囲であってもよい。広角光学系３０を用いることにより、眼科装置１１０は、内部照射角が２００°の領域１２Ａが走査対象とすることができる。つまり、光学スキャナ２２を制御することにより、渦静脈を含む所定範囲のＯＣＴ撮影を行う。眼科装置１１０は、当該ＯＣＴ撮影によってＯＣＴデータを生成することが可能となる。よって、眼科装置１１０は、ＯＣＴ画像である、渦静脈を含む眼底の断層画像（Ｂ－スキャン画像）、渦静脈を含むＯＣＴボリュームデータや、当該ＯＣＴボリュームデータの断面であるｅｎ－ｆａｃｅ画像（ＯＣＴボリュームデータから生成された正面画像）を生成することができる。なお、ＯＣＴ画像には、眼底中心部（黄斑や視神経乳頭などが存在する眼球の後極部）のＯＣＴ画像を含まれることは言うまでもない。
　ＯＣＴデータ（あるいはＯＣＴ画像の画像データ）は、通信インターフェース１６Ｆを介して眼科装置１１０からサーバ１４０へ送付され、記憶装置２５４に記憶される。
　なお、本実施形態では、光源２０Ａが波長掃引タイプのＳＳ－ＯＣＴ（Ｓｗｅｐｔ－Ｓｏｕｒｃｅ　ＯＣＴ）を例示するが、ＳＤ－ＯＣＴ（Ｓｐｅｃｔｒａｌ－Ｄｏｍａｉｎ　ＯＣＴ）、ＴＤ－ＯＣＴ（Ｔｉｍｅ－Ｄｏｍａｉｎ　ＯＣＴ）など、様々な方式のＯＣＴシステムであってもよい。
　次に、図３を参照して、サーバ１４０の電気系の構成を説明する。図３に示すように、サーバ１４０は、コンピュータ本体２５２を備えている。コンピュータ本体２５２は、ＣＰＵ２６２、ＲＡＭ２６６、ＲＯＭ２６４、入出力（Ｉ／Ｏ）ポート２６８を有する。入出力（Ｉ／Ｏ）ポート２６８には、記憶装置２５４、ディスプレイ２５６、マウス２５５Ｍ、キーボード２５５Ｋ、および通信インターフェース（Ｉ／Ｆ）２５８が接続されている。記憶装置２５４は、例えば、不揮発メモリで構成される。入出力（Ｉ／Ｏ）ポート２６８は、通信インターフェース（Ｉ／Ｆ）２５８を介して、ネットワーク１３０に接続されている。従って、サーバ１４０は、眼科装置１１０、およびビューワ１５０と通信することができる。
　ＲＯＭ２６４又は記憶装置２５４には、図６に示す画像処理プログラムが記憶されている。
　ＲＯＭ２６４又は記憶装置２５４は、本開示の技術の「メモリ」の一例である。ＣＰＵ２６２は、本開示の技術の「プロセッサ」の一例である。画像処理プログラムは、本開示の技術の「プログラム」の一例である。
　サーバ１４０は、眼科装置１１０から受信した各データを、記憶装置２５４に記憶する。
　サーバ１４０のＣＰＵ２６２が画像処理プログラムを実行することで実現される各種機能について説明する。図４に示すように、画像処理プログラムは、表示制御機能、画像処理機能、および処理機能を備えている。ＣＰＵ２６２がこの各機能を有する画像処理プログラムを実行することで、ＣＰＵ２６２は、表示制御部２０４、画像処理部２０６、および処理部２０８として機能する。
　次に、図５を用いて、サーバ１４０による画像処理を詳細に説明する。サーバ１４０のＣＰＵ２６２が画像処理プログラムを実行することで、図５のフローチャートに示された画像処理（画像処理方法）が実現される。
　ステップ６００で、画像処理部２０６は、ＯＣＴデータであるＯＣＴボリュームデータ４００を記憶装置２５４から取得する。ＯＣＴボリュームデータ４００は、図６に示すように、渦静脈ＶＶを含む所定面積のＯＣＴボリュームデータである。被検眼には、複数の渦静脈ＶＶが存在する。ＯＣＴボリュームデータ４００は、複数の渦静脈ＶＶの１つを眼科装置１１０によりＯＣＴ撮影をして得られる。所定面積は、例えば、６ｍｍ×６ｍｍの矩形領域である。後述するように、ＯＣＴボリュームデータ４００に対して第１面ｆ４０１から第Ｎ面ｆ４０Ｎまでの深さの異なるＮ枚の面が設定される。ＯＣＴボリュームデータ４００は、被検眼に複数存在する渦静脈ＶＶの各々が眼科装置１１０によりＯＣＴ撮影して得ておくようにしてもよい。
　なお、複数の渦静脈の各々の位置は、例えば、脈絡膜血管画像から脈絡膜血管を抽出し、各脈絡膜血管の移動方向（血管走行方向）を推定し、脈楽膜血管が集まる位置に基づいて推定される。
　脈絡膜血管画像は、赤色光（波長６３０～６６０ｎｍのレーザ光）で撮影されたＲ－ＵＷＦ－ＳＬＯ画像（Ｒ色眼底画像）と緑色光（波長５００～５５０ｎｍのレーザ光）で撮影されたＧ－ＵＷＦ－ＳＬＯ画像（Ｇ色眼底画像）の画像データを画像処理することにより生成される。具体的には、Ｇ色眼底画像から網膜血管を抽出し、Ｒ色眼底画像から網膜血管を除去し、脈絡膜血管を強調処理する画像処理を行うことにより脈絡膜血管画像が生成される。国際公開ＷＯ２０１９／１８１９８１の開示は、その全体が参照により本明細書に取り込まれる。
　ステップ６０２で、画像処理部２０６は、パラメータｎを１に設定する。当該パラメータｎはｅｎ－ｆａｃｅ画像の枚数を示すパラメータである。
　ステップ６０４では、画像処理部２０６は、ＯＣＴボリュームデータ４００を解析しＯＣＴボリュームデータ４００中の網膜色素上皮細胞層（Ｒｅｔｉｎａｌ　Ｐｉｇｍｅｎｔ　Ｅｐｉｔｈｅｌｉｕｍ、以下、ＲＰＥ層と称する）から所定画素数下、例えば、１０画素下の面を第１面として設定する。具体的には、画像処理部２０６は、まず、ＲＰＥ層を基準面として特定する。ＲＰＥ層は、ＯＣＴボリュームデータ４００に対し所定のセグメンテーション処理を行うことにより特定される。また、画像処理部２０６は、ＲＰＥ層を、ＯＣＴボリュームデータ４００の中で最も高輝層をＲＰＥ層として、特定するようにしてもよい。
　そして、画像処理部２０６は、ＲＰＥ層から１０画素下の面を第１面として設定する。これは、ＲＰＥ層より深い領域（眼球の中心から見てＲＰＥ層より遠い領域）が脈絡膜の領域であり、脈絡膜血管が存在する領域のｅｎ－ｆａｃｅ画像を生成するためである。
　ＲＰＥ層から１０画素下の面を第１面として設定することに限定されず、例えば、ＲＰＥ層のすぐ下に存在するブルッフ膜から１０画素下の面を第１面としてもよい。ブルッフ膜も、ＯＣＴボリュームデータ４００に対しＲＰＥ層のためのものとは異なる他の所定のセグメンテーション処理を行うことにより特定される。なお、１０画素下の位置を特定するために、ＯＣＴボリュームデータを生成したときのＡ－スキャンの方向の１０画素下としてもよい。
　ＲＰＥ層又はブルッフ膜から１０画素下の面を、第１面として特定することに限定されず、任意の画素数で設定するようにしてもよい。また画素数での定義ではなく、ミリメールやナノメートルなどの長さで定義してもよい。
　また、基準面として、瞳孔あるいは眼球中心から、一定距離を保った球面を基準面として定義してもよい。
　第１面は、本開示の技術の「他の面」の一例である。
　ＲＰＥ層又はブルッフ膜は、本開示の技術の「基準面」の一例である。
　ステップ６０６で、画像処理部２０６は、設定された第１面に対応する第１のｅｎ－ｆａｃｅ画像７ａ（図７参照）を生成する。ｅｎ－ｆａｃｅ画像７ａは第１面に存在する画素の画素値から生成されてもよいし、第１面を含む浅い方向の画素群および深い方向の画素群をＯＣＴボリュームデータ４００から抽出し、これらの画素群の輝度値の平均値あるいは中央値として画素値を求めてもよい。画素値を求める場合にノイズ除去などの画像処理を用いてもよい。
　生成された第１面に対応する１枚目のｅｎ－ｆａｃｅ画像は、処理部２０８によりＲＡＭ２６６に保存される。
　なお、基準面であるＲＰＥ層又はブルッフ膜に対応するｅｎ－ｆａｃｅ画像を生成するようにしてもよい。
　ステップ６０８で、画像処理部２０６は、パラメータｎがあらかじめ設定された最大数Ｎであるか否かを判断する。最大数Ｎはｅｎ－ｆａｃｅ画像の生成枚数を示し、８、１０、１５など任意の整数である。本実施の形態では、Ｎ＝８としている。つまり、ＲＰＥ層から１０画素下の第１面、さらにＲＰＥ層から２０画素下の第２面と続き、ＲＰＥ層から８０画素下の第８面まで、面が設定され、合計８枚の各面に対応するｅｎ－ｆａｃｅ画像を画像処理部２０６が生成する。なお、最大数Ｎは、渦静脈付近の網膜の厚さやｅｎ－ｆａｃｅ画像の間隔の長さによってユーザにより任意に設定可能である。
　パラメータｎが最大数Ｎであると判断された場合（ステップ６０８：ＹＥＳ）は、画像処理はステップ６１０に進む。ステップ６１０で、処理部２０８は、生成されたＮ枚のｅｎ－ｆａｃｅ画像をＲＡＭ２６６あるいは記憶装置２５４に保存し、処理を終了する。
　しかし、パラメータｎが最大数Ｎではないと判断された場合（ステップ６０８：ＮＯ）は、画像処理はステップ６１２に進む。ステップ６１２で、画像処理部２０６は、パラメータｎをｎ＋１と設定する。
　次のステップ６１４で、画像処理部２０６は第ｎ面を設定する。例えば、ｎ＝２の場合、すでに設定された第１面から深さ方向に１０画素シフトした面を第２面として設定する。あるいは、ｎ＝２の場合、ＲＰＥ層から２０画素（１０画素×２）シフトした面を第２面とするようにしてもよい。
　ステップ６１４の処理が終了すると、画像処理は、ステップ６０６に戻る。ステップ６０６で、画像処理部２０６は、設定された第ｎ（例えば、２）面に対応する第ｎ（例えば、２）のｅｎ－ｆａｃｅ画像７ｂ（図７参照）を生成する。生成された第２のｅｎ－ｆａｃｅ画像７ｂは、ＲＡＭ２６６に保存される。
　このようにして、画像処理部２０６は、ステップ６０６、ステップ６０８、ステップ６１２、ステップ６１４のループをパラメータｎが最大数Ｎになるまで繰り返す。
　図５のフローチャートに示された画像処理を、画像処理部２０６が実行することにより、図７に示す第１のｅｎ－ｆａｃｅ画像７ａからの第Ｎのｅｎ－ｆａｃｅ画像７ｄが生成される。これらの複数のｅｎ－ｆａｃｅ画像は、脈絡膜を、深さが深くなる方向に所定画素数ごとに生成されたｅｎ－ｆａｃｅ画像群である。このｅｎ－ｆａｃｅ画像群を並べて表示したり、１つのｅｎ－ｆａｃｅ画像を１フレームとして、Ｎフレームからなる動画像として表示したりすることにより、脈絡膜の中における渦静脈の走行方向や、渦静脈の太さなどを可視化することができる。
　例えば、第１のｅｎ－ｆａｃｅ画像７ａ、第２のｅｎ－ｆａｃｅ画像７ｂ、・・・第Ｎ-１のｅｎ－ｆａｃｅ画像７ｃ、及び第Ｎのｅｎ－ｆａｃｅ画像７ｄを順に比較すると、渦静脈の走行方向、特に、深さ方向の走行方向を特定することができる。図７に示す例では紙面の左上の方向に渦静脈が走行していることがわかる。
　なお、図７に示したようなｅｎ－ｆａｃｅ画像群を画像処理することにより、あるいは、ｅｎ－ｆａｃｅ画像群と渦静脈の走行方向を機械学習した学習済みモデルを用いた人工知能により、渦静脈の血管走行を特定するようにしてもよい。あるいは、ＯＣＴボリュームデータと渦静脈の走行方向を機械学習した学習済みモデルを用いた人工知能により、渦静脈の血管走行を特定するようにしてもよい。
　特定された血管走行に基づいて、血管走行の方向を示す情報を生成し、血管走行の方向を示す情報を第１から第Ｎのｅｎ－ｆａｃｅ画像に重畳表示してもよい。また、ＯＣＴボリュームデータ４００に重畳表示するようにしてもよい。
　血管走行の方向を示す情報としては、ベクトル情報、方向を示す矢印などのマーク情報や、走行方向に流れていくような点滅を行うための情報などが含まれる。
　また、第１から第Ｎのｅｎ－ｆａｃｅ画像の各ｅｎ－ｆａｃｅ画像における血管を示す画素の領域を求め、血管の重心位置を計算する。そして、第１のｅｎ－ｆａｃｅ画像から第Ｎのｅｎ－ｆａｃｅ画像にかけて（上のｅｎ－ｆａｃｅ画像から下のｅｎ－ｆａｃｅ画像にかけて）計算された重心位置をつなぐベクトルを、第１ｅｎ－ｆａｃｅ画像から第Ｎｅｎ－ｆａｃｅ画像の全てのｅｎ－ｆａｃｅ画像に表示するようにしてもよい。また、各ｅｎ－ｆａｃｅ画像における当該重心位置をＯＣＴボリュームデータに表示するようにしてもよい。
　本実施の形態では、ＲＰＥ層から所定画素数（上記例では１０画素）下の第１面を設定し、さらに順に第２面・・・第８面まで、面が設定され、合計８枚の各面に対応するｅｎ－ｆａｃｅ画像が生成される。しかし、本開示の技術はこれに限定されない。例えば、ＲＰＥ層から４０画素数下の面を基準面として設定し、基準面から１０画素数毎に、上下に他の面を設定し、各面のｅｎ－ｆａｃｅ画像を生成するようにしてもよい。
　以下、生成されたｅｎ－ｆａｃｅ画像群を表示するためのディスプレイ・スクリーンについて説明する。当該ディスプレイ・スクリーンは、ユーザの指示に基づきサーバ１４０の表示制御部２０４により生成され、処理部２０８により、ビューワ１５０に画像信号として出力される。ビューワ１５０は、当該画像信号に基づき、ディスプレイ・スクリーンをモニターなどに表示する。
　図８には、第１のディスプレイ・スクリーン５００Ａが示されている。図８に示すように、第１のディスプレイ・スクリーン５００Ａは、インフォメーションエリア５０２と、イメージディスプレイエリア５０４Ａとを有する。
　インフォメーションエリア５０２には、患者ＩＤディスプレイフィールド５１２、患者名ディスプレイフィールド５１４、年齢ディスプレイフィールド５１６、視力ディスプレイフィールド５１５、右眼／左眼ディスプレイフィールド５２０、及び眼軸長ディスプレイフィールド５２２を有する。患者ＩＤディスプレイフィールド５１２から眼軸長ディスプレイフィールド５２２の各表示領域には、ビューワ１５０が、サーバ１４０から受信した情報に基づいて、各々の情報を表示する。
　イメージディスプレイエリア５０４Ａは、被検眼像等を表示する領域である。イメージディスプレイエリア５０４Ａには、以下の各表示フィールドが設けられている、具体的には、ＵＷＦ眼底画像表示フィールド５４２、脈絡膜血管画像表示フィールド５４４、Ｂ－スキャン画像表示フィールド５４６、ｅｎ－ｆａｃｅ画像表示フィールド５４８、及びコメントフィールド５５０がある。
　コメントフィールド５３０は、ユーザである眼科医が観察した結果、又は診断結果を任意に入力できる備考欄である。
　ＵＷＦ眼底画像表示フィールド５４２には、被検眼の眼底を眼科装置１１０でＳＬＯ撮影ＵＷＦ－ＳＬＯ眼底画像ＳＬが表示されている。ＵＷＦ－ＳＬＯ眼底画像ＳＬには縦方向の分水嶺を示すラインＳＬＶ、横方向の分水嶺を示すラインＳＬＨが重畳されている。さらに、渦静脈周辺のＯＣＴボリュームデータ４００の取得位置を示す矩形領域ＯＣＶ１が重畳表示されている。さらに、矩形領域ＯＣＶ１の中に、Ｂ－スキャン画像表示フィールド５２６に表示されているＢ－スキャン画像ＯＣＢの断面位置を示すラインＯＣＢ１が重畳表示されている。
　脈絡膜血管画像表示フィールド５４４には、ＵＷＦ－ＳＬＯ眼底画像ＳＬ１を画像処理することにより得られた脈絡膜血管画像ＣＶが表示されている。脈絡膜血管画像ＣＶには縦方向の分水嶺を示すラインＣＶＶ、及び横方向の分水嶺を示すラインＣＶＨが重畳されている。さらに、渦静脈周辺のＯＣＴボリュームデータ４００の取得位置を示す矩形領域ＯＣＶ２が重畳表示されている。さらに、矩形領域ＯＣＶ２の中に、Ｂ－スキャン画像表示フィールド５２６に表示されているＢ－スキャン画像ＯＣＢの断面位置を示すラインＯＣＢ２が重畳表示されている。ユーザの指示により、Ｂ－スキャン画像表示フィールド５４６およびｅｎ－ｆａｃｅ画像表示フィールド５４８に表示する渦静脈を選択するようにしてもよい。この場合、選択された渦静脈のＯＣＴボリュームデータの範囲に矩形領域ＯＣＶ２が重畳表示される。
　Ｂ－スキャン画像表示フィールド５４６には、ＯＣＴボリュームデータ４００から生成されたＢ－スキャン画像ＯＣＢが表示されている。Ｂ－スキャン画像ＯＣＢに重畳されている曲線ＯＣＬは、ｅｎ－ｆａｃｅ画像の取得面の深さ方向の位置を示す曲線である。この曲線ＯＣＬは、ｅｎ－ｆａｃｅ画像表示フィールド５４８に表示されるｅｎ－ｆａｃｅ画像が第１のｅｎ－ｆａｃｅ画像７ａから第Ｎのｅｎ－ｆａｃｅ画像７ｄのいずれかに変化すると、曲線ＯＣＬの重畳位置も変化する。また、ｅｎ－ｆａｃｅ画像表示フィールド５４８にｅｎ－ｆａｃｅ画像群から生成された動画像が表示される場合は、動画像の再生とともに曲線ＯＣＬの重畳位置も変化する。
　ｅｎ－ｆａｃｅ画像表示フィールド５４８には、第１のｅｎ－ｆａｃｅ画像７ａから第Ｎのｅｎ－ｆａｃｅ画像７ｄまでの複数のｅｎ－ｆａｃｅ画像を順にスライドショーのように表示させたり、動画表示させたり、特定のｅｎ－ｆａｃｅ画像を静止画として表示させたり、することができる。ｅｎ－ｆａｃｅ画像表示フィールド５４８に表示されるｅｎ－ｆａｃｅ画像に応じて、Ｂ－スキャン画像表示フィールド５４６に表示されているＢ－スキャン画像ＯＣＢ上の、曲線ＯＣＬの重畳位置が変化する。また、ｎ－ｆａｃｅ画像表示フィールド５４８に表示されるｅｎ－ｆａｃｅ画像上に、Ｂ－スキャン画像ＯＣＢの断面位置を示すラインＯＣＢ３を重畳表示するようにしてもよい。ラインＯＣＢ３はユーザの指示により重畳する／重畳しないを選択するようにしてもよい。このラインＯＣＢ３を水平移動させることにより、Ｂ－スキャン画像表示フィールド５４６に表示させるＢ-スキャン画像ＯＣＢを変更することができる。なお、ラインＯＣＢ３は水平方向だけでなく、垂直方向など任意の方向に設定することもできる。
　また、前述した、血管走行の方向を示す情報を第１から第Ｎのｅｎ－ｆａｃｅ画像に重畳表示してもよい。
　図９には、第２のディスプレイ・スクリーン５００Ｂが示されている。第２のディスプレイ・スクリーン５００Ｂは、第１のディスプレイ・スクリーン５００Ａのフィールドと同じフィールドを有するので、同一のフィールドには同一の符号を付してその説明を省略し、異なる部分を説明する。
　ディスプレイ・スクリーン５００Ｂは、Ｎ枚のｅｎ－ｆａｃｅ画像を並べて表示させるｅｎ－ｆａｃｅ画像表示フィールド５６０を有する。ｅｎ－ｆａｃｅ画像表示フィールド５６０には、Ｎ＝８の場合における、第１のｅｎ－ｆａｃｅ画像から第８のｅｎ－ｆａｃｅ画像が表示されている。ユーザの入力指示に応じて、例えば、第２のｅｎ－ｆａｃｅ画像７ｂが選択されると、点線で示す矩形枠５６０ｆが第２のｅｎ－ｆａｃｅ画像７ｂに重畳表示される。それとともに、選択されたｅｎ－ｆａｃｅ画像の面の位置が、Ｂ－スキャン画像表示フィールド５４６のＢ－スキャン画像ＯＣＢに曲線ＯＣＬとして表示される。
　ところで、第１のディスプレイ・スクリーン５００Ａ、あるいは、第２のディスプレイ・スクリーン５００Ｂには、ＯＣＴボリュームデータ４００を画像処理した３次元画像を表示するようにしてもよい。
　例えば、ＯＣＴボリュームデータ４００、血管領域を黒、血管以外の空間（間質）を白に二値化処理により、血管領域を抽出し、その血管領域を３Ｄボリュームレンダリング処理を実行し表示するようにしてもよい。表示の際にＳＴＬ（Standard Triangulated Language）ファイル形式などのポリゴン型メッシュ構造のファイル形式に変換してから表示させてもよい。
　また、ＯＣＴボリュームデータ４００から生成された複数のｅｎ－ｆａｃｅ画像の各々を二値化処理して二値化したｅｎ－ｆａｃｅ画像を生成する。そして二値化したｅｎ－ｆａｃｅ画像を再合成することにより二値化したＯＣＴボリュームデータ４００を得るようにしてもよい。
　さらには、渦静脈や渦静脈に接続する脈絡膜血管を可視化する画像処理をＯＣＴボリュームデータ４００に対して実行し、渦静脈周辺の血管を強調した３次元画像を表示するようにしてもよい。この場合、眼底周辺部の網膜や脈絡膜は中心部と比較して薄い場合があり、実寸表示すると渦静脈の立体構造が観察しづらい。そのため、ＯＣＴボリュームデータ４００を深さ方向だけ伸ばす（表示単位を深さ方向だけ変える）画像処理を行い、３次元画像を表示するようにしてもよい。
　以上説明したように本実施の形態では、渦静脈含む領域が撮影されたＯＣＴボリュームデータからｅｎ－ｆａｃｅ画像を生成することにより、脈絡膜から強膜に向けて走行する渦静脈の状況を可視化することが可能となる。渦静脈は横方向の変化に比べ深さ方向の変化が著しいが、ＯＣＴボリュームデータから生成されたｅｎ－ｆａｃｅ画像を用いることにより、ＯＣＴボリュームデータ中の脈絡膜血管構造を可視化しただけではわかりにくかった、渦静脈の深さ方向の形態変化を高解像度で可視化及び定量化することができるようになる。
　また、本実施の形態では、ＯＣＴアンジオグラフィーのアルゴリムによる血管造影処理を用いることなく、渦静脈を可視化する。そのため、ＯＣＴアンジオグラフィーのように異なる時刻のＯＣＴボリュームデータを取得することなく、１つのＯＣＴボリュームデータから渦静脈を可視化することができる。
　さらに、図１０に示すように眼球１２Ｖの外へ延びる渦静脈１２Ｖ１や１２Ｖ２について、脈絡膜内の渦静脈の状況を観察することが可能となる。脈絡膜内の渦静脈の状況を観察することにより、眼球の変形などによる渦静脈の変形に伴う血流の変化を推定することが可能となる。　
　上記実施の形態では、画像処理（図５）は、サーバ１４０が実行しているが、本開示の技術はこれに限定されず、眼科装置１１０、ビューワ１５０、又は、ネットワーク１３０に更に設けた追加画像処理装置が実行してもよい。
　本開示において、各構成要素（装置等）は、矛盾が生じない限りは、１つのみ存在しても２つ以上存在してもよい。
　以上説明した各例では、コンピュータを利用したソフトウェア構成により画像処理が実現される場合を例示したが、本開示の技術はこれに限定されるものではない。例えば、コンピュータを利用したソフトウェア構成に代えて、ＦＰＧＡ（Ｆｉｅｌｄ－Ｐｒｏｇｒａｍｍａｂｌｅ　Ｇａｔｅ　Ａｒｒａｙ）またはＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ　Ｓｐｅｃｉｆｉｃ　Ｉｎｔｅｇｒａｔｅｄ　Ｃｉｒｃｕｉｔ）等のハードウェア構成のみによって、画像処理が実行されるようにしてもよい。画像処理のうちの一部の処理がソフトウェア構成により実行され、残りの処理がハードウェア構成によって実行されるようにしてもよい。
　このように本開示の技術は、コンピュータを利用したソフトウェア構成により画像処理が実現される場合とされない場合とを含むので、以下の技術を含む。
（第１の技術）
　渦静脈を含む領域のＯＣＴボリュームデータを取得する取得部と、
　前記ＯＣＴボリュームデータに基づいて、深さが異なる複数の面に対応した複数のｅｎ－ｆａｃｅ画像を生成する生成部と、
　を備える画像処理装置。
（第２の技術）
　取得部が、渦静脈を含む領域のＯＣＴボリュームデータを取得するステップと、
　生成部が、前記ＯＣＴボリュームデータに基づいて、深さが異なる複数の面に対応した複数のｅｎ－ｆａｃｅ画像を生成するステップと、
　を含む画像処理方法。
　画像処理部２０６は、本開示の技術の「取得部」及び「生成部」の一例である。
　以上の開示内容から以下の技術が提案される。
（第３の技術）
　画像処理するためのコンピュータープログラム製品であって、
　前記コンピュータープログラム製品は、それ自体が一時的な信号ではないコンピュータ可読記憶媒体を備え、
　前記コンピュータ可読記憶媒体には、プログラムが格納されており、
　前記プログラムは、
　コンピュータに、
　渦静脈を含む領域のＯＣＴボリュームデータを取得するステップと、
　前記ＯＣＴボリュームデータに基づいて、深さが異なる複数の面に対応した複数のｅｎ－ｆａｃｅ画像を生成するステップと、
　を実行させる、
　コンピュータープログラム製品。
　サーバ１４０は、本開示の技術の「コンピュータープログラム製品」の一例である。
　以上説明した各画像処理はあくまでも一例である。従って、主旨を逸脱しない範囲内において不要なステップを削除したり、新たなステップを追加したり、処理順序を入れ替えたりしてもよいことは言うまでもない。
　本明細書に記載された全ての文献、特許出願、及び技術規格は、個々の文献、特許出願、及び技術規格が参照により取り込まれることが具体的にかつ個々に記載された場合と同様に、本明細書中に参照により取り込まれる。

Claims (10)

  1. 　プロセッサが行う画像処理方法であって、
    　渦静脈を含む領域のＯＣＴボリュームデータを取得するステップと、
    　前記ＯＣＴボリュームデータに基づいて、深さが異なる複数の面に対応した複数のｅｎ－ｆａｃｅ画像を生成するステップと、
    　を含む画像処理方法。
  2. 　前記複数の面は、脈絡膜内に設定されることを特徴とする請求項１に記載の画像処理方法。
  3. 　前記複数の面は、所定間隔を置いて設定されることを特徴とする請求項１または２に記載の画像処理方法。
  4. 　前記複数の面は、基準面より深い位置に設定されることを特徴とする請求項１ないし３の何れか１項に記載の画像処理方法。
  5. 　前記基準面は、網膜色素上皮細胞層、あるいは、ブルッフ膜であることを特徴とする請求項４に記載の画像処理方法。
  6. 　前記複数のｅｎ－ｆａｃｅ画像を表示する表示画面を生成するステップを更に含む、請求項１から請求項５の何れか１項に記載の画像処理方法。
  7. 　前記複数のｅｎ－ｆａｃｅ画像を用いて動画像を生成するステップを更に含む請求項１から請求項６の何れか１項に記載の画像処理方法。
  8. 　前記表示画面に、前記複数のｅｎ－ｆａｃｅ画像が、前記動画像により表示される、請求項６に従属する請求項７に記載の画像処理方法。
  9. 　メモリと、前記メモリに接続するプロセッサとを備え、
    　前記プロセッサは、
    　渦静脈を含む領域のＯＣＴボリュームデータを取得するステップと、
    　前記ＯＣＴボリュームデータに基づいて、深さが異なる複数の面に対応した複数のｅｎ－ｆａｃｅ画像を生成するステップと、
    　を実行する画像処理装置。
  10. 　コンピュータに、
    　渦静脈を含む領域のＯＣＴボリュームデータを取得するステップと、
    　前記ＯＣＴボリュームデータに基づいて、深さが異なる複数の面に対応した複数のｅｎ－ｆａｃｅ画像を生成するステップと、
    　を実行させるプログラム。
PCT/JP2021/023441 2020-11-26 2021-06-21 画像処理方法、画像処理装置、及びプログラム Ceased WO2022113409A1 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2020-196330 2020-11-26
JP2020196330 2020-11-26

Publications (1)

Publication Number Publication Date
WO2022113409A1 true WO2022113409A1 (ja) 2022-06-02

Family

ID=81754436

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2021/023441 Ceased WO2022113409A1 (ja) 2020-11-26 2021-06-21 画像処理方法、画像処理装置、及びプログラム

Country Status (1)

Country Link
WO (1) WO2022113409A1 (ja)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014200093A1 (ja) * 2013-06-14 2014-12-18 国立大学法人名古屋大学 光断層画像撮影装置
JP2016026521A (ja) * 2014-06-30 2016-02-18 株式会社ニデック 光コヒーレンストモグラフィ装置、及びデータ処理プログラム
JP2017047111A (ja) * 2015-09-04 2017-03-09 キヤノン株式会社 眼科装置、表示制御方法およびプログラム
WO2019203311A1 (ja) * 2018-04-18 2019-10-24 株式会社ニコン 画像処理方法、プログラム、及び画像処理装置
JP2020513983A (ja) * 2017-03-23 2020-05-21 ドヘニー アイ インスティテュート 光干渉断層撮影の多重ｅｎ−ｆａｃｅ血管造影平均化システム、方法および装置
JP2020166814A (ja) * 2019-03-11 2020-10-08 キヤノン株式会社 医用画像処理装置、医用画像処理方法及びプログラム
US20200394789A1 (en) * 2019-06-12 2020-12-17 Carl Zeiss Meditec Inc Oct-based retinal artery/vein classification

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014200093A1 (ja) * 2013-06-14 2014-12-18 国立大学法人名古屋大学 光断層画像撮影装置
JP2016026521A (ja) * 2014-06-30 2016-02-18 株式会社ニデック 光コヒーレンストモグラフィ装置、及びデータ処理プログラム
JP2017047111A (ja) * 2015-09-04 2017-03-09 キヤノン株式会社 眼科装置、表示制御方法およびプログラム
JP2020513983A (ja) * 2017-03-23 2020-05-21 ドヘニー アイ インスティテュート 光干渉断層撮影の多重ｅｎ−ｆａｃｅ血管造影平均化システム、方法および装置
WO2019203311A1 (ja) * 2018-04-18 2019-10-24 株式会社ニコン 画像処理方法、プログラム、及び画像処理装置
JP2020166814A (ja) * 2019-03-11 2020-10-08 キヤノン株式会社 医用画像処理装置、医用画像処理方法及びプログラム
US20200394789A1 (en) * 2019-06-12 2020-12-17 Carl Zeiss Meditec Inc Oct-based retinal artery/vein classification

Similar Documents

Publication Publication Date Title
JP2023009530A (ja) 画像処理方法、画像処理装置、及びプログラム
JP2026053471A (ja) 画像処理方法、プログラム、及び画像処理装置
JP7441783B2 (ja) 画像処理方法、プログラム、眼科装置、及び脈絡膜血管画像生成方法
JP2019177032A (ja) 眼科画像処理装置、および眼科画像処理プログラム
JP7585800B2 (ja) 画像処理方法、画像処理装置、及びプログラム
JP2018019771A (ja) 光コヒーレンストモグラフィ装置、および光コヒーレンストモグラフィ制御プログラム
JP2024105697A (ja) 画像処理方法、画像処理装置、及び画像処理プログラム
JP2023030184A (ja) 画像処理方法、プログラム、画像処理装置、及び眼科システム
CN111989030A (zh) 图像处理方法、程序及图像处理装置
JP2025122054A (ja) 画像処理方法、画像処理装置、及びプログラム
JP2022089086A (ja) 画像処理方法、画像処理装置、及び画像処理プログラム
JP7715186B2 (ja) 画像処理方法、画像処理装置、及びプログラム
CN115052514A (zh) 图像处理方法、图像处理装置、以及程序
WO2023282339A1 (ja) 画像処理方法、画像処理プログラム、画像処理装置及び眼科装置
JP7419946B2 (ja) 画像処理方法、画像処理装置、及び画像処理プログラム
JP2023066198A (ja) 情報出力装置、眼底画像撮影装置、情報出力方法、及び情報出力プログラム
JP7272453B2 (ja) 画像処理方法、画像処理装置、およびプログラム
JP7800663B2 (ja) 画像処理方法、画像処理装置、及びプログラム
JP7758171B2 (ja) 画像処理方法、画像処理装置、及びプログラム
JP2021058285A (ja) 画像処理装置、画像処理方法、及びプログラム
JP7677411B2 (ja) 画像処理方法、画像処理装置、及びプログラム
WO2022181729A1 (ja) 画像処理方法、画像処理装置、及び画像処理プログラム
WO2021210281A1 (ja) 画像処理方法、画像処理装置、及び画像処理プログラム

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 21897389

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 21897389

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: JP