US20150124218A1 - Image processing apparatus, ophthalmologic imaging apparatus, image processing method, and storage medium - Google Patents

Image processing apparatus, ophthalmologic imaging apparatus, image processing method, and storage medium Download PDF

Info

Publication number
US20150124218A1
US20150124218A1 US14/597,127 US201514597127A US2015124218A1 US 20150124218 A1 US20150124218 A1 US 20150124218A1 US 201514597127 A US201514597127 A US 201514597127A US 2015124218 A1 US2015124218 A1 US 2015124218A1
Authority
US
United States
Prior art keywords
image
fundus
eye
information
unit
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.)
Abandoned
Application number
US14/597,127
Inventor
Keiko Yonezawa
Kazuhide Miyata
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.)
Canon Inc
Original Assignee
Canon Inc
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 Canon Inc filed Critical Canon Inc
Priority to US14/597,127 priority Critical patent/US20150124218A1/en
Publication of US20150124218A1 publication Critical patent/US20150124218A1/en
Priority to US16/287,926 priority patent/US11058293B2/en
Abandoned legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/0016Operational features thereof
    • A61B3/0025Operational features thereof characterised by electronic signal processing, e.g. eye models
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/0016Operational features thereof
    • A61B3/0041Operational features thereof characterised by display arrangements
    • A61B3/0058Operational features thereof characterised by display arrangements for multiple images
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/12Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes
    • A61B3/1225Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes using coherent radiation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/14Arrangements specially adapted for eye photography
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10056Microscopic image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20056Discrete and fast Fourier transform, [DFT, FFT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30041Eye; Retina; Ophthalmic

Definitions

  • Fundus cameras and scanning laser ophthalmoscopes are among the apparatuses used for the inspection of the fundus of the eye.
  • the fundus camera captures an image of a fundus of the eye by receiving reflected light of a light beam which has entered the fundus of the eye.
  • the SLO is an ophthalmologic apparatus that uses the principle of confocal laser scanning microscope.
  • fundus cameras and SLOs including an adaptive optical system have been developed and are used for acquiring fundus planar images of high lateral resolution.
  • the adaptive optical system measures an aberration of a subject's eye by a wavefront sensor in real time and corrects the aberration of the measuring beam and the return beam that occurs at the subject's eye by a wavefront correcting device. Further, attempts are being made to capture images of photoreceptor cells of a retina using by these apparatuses and make a diagnosis of a disease or evaluate drug response.
  • an image processing apparatus includes an identification unit configured to identify periodicity of a fundus image obtained by capturing an image of a fundus of an eye, and an information acquisition unit configured to acquire information indicating an imaging state of photoreceptor cells in the fundus image based on the periodicity.
  • FIG. 1 illustrates a functional configuration of an image processing apparatus according to a first exemplary embodiment.
  • FIG. 2 is a flowchart illustrating processing procedures of the image processing apparatus according to the first exemplary embodiment.
  • FIG. 3 is a schematic view of a high-precision planar image of a fundus of the eye including photoreceptor cells acquired by an adaptive optics SLO apparatus.
  • FIG. 4 illustrates an example of a Fourier image obtained by frequency conversion of the planar image of the fundus of the eye.
  • FIGS. 5A and 5B illustrate a method for calculating a structure that reflects an arrangement of the photoreceptor cells from a Fourier image and a graph organizing the result of the calculation.
  • FIGS. 6A and 6B illustrate an example of a Fourier image obtained by frequency conversion of the planar image of the fundus of the eye in a case where the signal is weak and a graph indicating the arrangement of the photoreceptor cells.
  • FIG. 7 illustrates a feature quantity acquired from a Fourier image.
  • FIGS. 8A , 8 B, and 8 C illustrate an example of a Fourier image of low resolution and feature quantities acquired from the image.
  • FIGS. 9A , 9 B, and 9 C illustrate an example of a Fourier image of even lower resolution and feature quantities acquired from the image.
  • FIGS. 10A and 10B illustrate an image quality index acquired from a feature amount extracted from a Fourier image.
  • FIG. 11 is a functional configuration of the image processing apparatus according to a second exemplary embodiment.
  • FIG. 12 is a flowchart illustrating processing procedures of the image processing apparatus according to the second exemplary embodiment.
  • FIG. 13 illustrates an example of dividing the planar image of the fundus of the eye into a plurality of local planar images of the fundus of the eye.
  • FIGS. 14A and 14B illustrate examples of indices acquired from the local planar images of the fundus of the eye.
  • FIG. 15 illustrates a hardware configuration of the image processing apparatus according to another exemplary embodiment.
  • FIG. 16 illustrates a configuration of an ophthalmologic imaging system.
  • FIG. 17 illustrates a configuration of an ophthalmologic imaging apparatus.
  • FIGS. 18A , 18 B, and 18 C illustrate an image acquisition method for optical coherence tomography (OCT).
  • OCT optical coherence tomography
  • FIGS. 19A , 19 B, 19 C, and 19 D illustrate an image acquisition method for SLO.
  • an image quality index value is obtained according to a comparison of a noise level of a region being an evaluation target and a noise level of a different region.
  • the comparison method cannot be used for the evaluation.
  • an index that objectively evaluates an image quality of an image of photoreceptor cells of a retina acquired by an ophthalmologic apparatus such as an adaptive optics SLO is presented.
  • the adjustment of the imaging conditions and the diagnosis of the image can be easily performed.
  • an index that quantitatively indicates the image quality of the acquired image is calculated and presented. More specifically, a spatial frequency image is acquired by discrete Fourier transform using a planar image of a fundus of the eye which has been acquired by the adaptive optics SLO.
  • the acquired spatial frequency image is hereinafter referred to as a Fourier image or a frequency image.
  • a feature quantity of a periodic structure that reflects the regular arrangement of the photoreceptor cells is extracted from the acquired Fourier image, and the index of the image quality is acquired from the extracted feature quantity.
  • the user can determine whether the image needs to be captured again. Further, if a diagnosis is to be performed based on the density of the photoreceptor cells, the user can determine whether it is adequate to use such an image for the diagnosis.
  • FIG. 3 schematically illustrates a planar image of the fundus of the eye captured by the adaptive optics SLO.
  • a small region having relatively high luminance can be extracted in a distinguishable manner as a photoreceptor cell PR.
  • a blood vessel region V having low luminance compared to the luminance of the photoreceptor cell may be extracted.
  • the blood vessel region V is the shadow of the blood vessel above the photoreceptor cells.
  • FIG. 4 illustrates an example of a Fourier image acquired by the discrete Fourier transform of spatial frequency components of the planar image of the fundus of the eye described above. As illustrated in FIG. 4 , a ring that corresponds to the periodicity of the photoreceptor cells is formed according to the periodic arrangement of the photoreceptor cells.
  • FIG. 1 illustrates a functional configuration of an image processing apparatus 10 according to the present embodiment.
  • An image acquisition unit 100 in FIG. 1 acquires a planar image of a fundus of the eye.
  • the acquired planar image is an image of photoreceptor cells of the fundus of the eye.
  • it is an image acquired by a fundus imaging apparatus which includes an aberration measurement unit such as a Hartmann-Shack wavefront sensor and an adaptive optical system which corrects the aberration.
  • the adaptive optical system includes, for example, a reflective display panel having liquid crystals arranged on a silicon substrate (liquid crystal on silicon (LCOS)) or a deformable mirror.
  • LCOS liquid crystal on silicon
  • an input information acquisition unit 110 acquires information of the subject's eye which is being captured.
  • the acquired image is stored in a storage unit 130 via a control unit 120 .
  • An image processing unit 140 includes a frequency conversion unit 141 , a feature extraction unit 142 , and an index calculation unit 143 .
  • the frequency conversion unit 141 performs frequency conversion of the planar image of the fundus of the eye and obtains a frequency image. For example, the discrete Fourier transform is used for converting frequencies, as described above.
  • the feature extraction unit 142 extracts the feature quantity indicating the periodic structure of the planar image of the fundus of the eye from the frequency image.
  • the periodic structure of the planar image of the fundus of the eye is, for example, appears in a ring structure having a point of origin at the center.
  • the ring of the frequency image presents a specified frequency band corresponding to the structure of the photoreceptor cells.
  • the feature quantity of the periodic structure is, for example, the feature quantity concerning this ring structure such as a peak luminance value in the ring structure and a sharpness value that indicates the extension of the ring structure. Further, there is the peak position of the luminance value in the ring structure. Additionally, if the image quality is poor, a disk-like structure may appear in the frequency image rather than the ring structure. In this case, the extension and the position of the disk region, and the peak luminance value in the disk region are extracted as the feature quantities. These are used for acquiring information of the imaging state and the distribution information of the photoreceptor cells in the planar image of the fundus of the eye.
  • the periodicity can be identified using information other than the frequency image. For example, the peak value of the luminance near the center of the photoreceptor cells is detected and an average of the profile of luminance distribution around each detection point is acquired. In this manner, the periodicity of the image of the fundus of the eye can be identified.
  • the peak value can be detected using maximum value detection processing.
  • the profile in this context is in the radial direction having the detection point at the center, and the mean value is taken in the angular direction of the polar coordinates having the detection point at the center. Then, the feature extraction unit 142 evaluates the shape. If the shape is periodic, a peak corresponding to a detection point adjacent to the detection point appears on the average profile.
  • co-occurrence matrix and fractal dimension known as the texture feature quantity can be used.
  • the identification of the periodicity by the Fourier transform is one exemplary embodiment, and the above-described different methods can be used in identifying the periodicity.
  • the frequency conversion unit 141 functions as one exemplary embodiment of an identification unit of the periodicity.
  • the index calculation unit 143 functions as an information acquisition unit configured to acquire information of the imaging state of the photoreceptor cells and distribution information of the photoreceptor cells from the image indicating the periodicity such as a Fourier image. By using the feature quantity extracted by the feature extraction unit 142 , the index calculation unit 143 can obtain accurate calculation results of the imaging state and the distribution of the photoreceptor cells.
  • An output unit 150 outputs the information of the imaging state and the distribution information of the photoreceptor cells acquired by the index calculation unit 143 .
  • the information is output to a display unit 160 as well as an external database and an output apparatus.
  • the information is, for example, printed out by the output apparatus.
  • Information of the imaging state can be displayed as it is or together with the planar image of the fundus of the eye on the output apparatus.
  • the information of the imaging state which is displayed there is the image quality index value of the planar image of the fundus of the eye where the photoreceptor cells have been captured, and information asking the user to perform operation to improve the imaging state.
  • a frequency image including the information of the imaging state can be displayed together with the image of the fundus of the eye.
  • a graph showing a relation between a distance from the center position of the ring structure or from the point of origin of the frequency image and the luminance value is displayed.
  • a control value corresponding to the imaging state of the photoreceptor cells is output to the ophthalmologic imaging apparatus which has captured the planar image of the fundus of the eye.
  • the image processing unit 140 generates a Fourier image from the acquired planar image of the fundus of the eye, calculates an index of the image quality from the feature quantity extracted from the Fourier image, and stores the index in the storage unit 13 .
  • the output unit 150 outputs the calculated index to a monitor. Further, the output unit 150 outputs the result of the processing stored in the storage unit 130 to the database.
  • the image processing apparatus 10 acquires the planar image of the fundus of the eye directly from the adaptive optics SLO in the present embodiment, the image can be acquired via a network.
  • a plurality of planar images of the fundus of the eye captured by the adaptive optics SLO are stored in a database connected to the adaptive optics SLO via the network, and the image processing apparatus 10 acquires the images from the database via the network.
  • step S 210 the image acquisition unit 100 acquires a plurality of planar images of the fundus of the eye to be analyzed from the adaptive optics SLO connected to the image processing apparatus 10 or a database where the planar images of the fundus of the eye captured by the apparatus are stored.
  • the acquired planar images of the fundus of the eye are stored in the storage unit 130 via the control unit 120 .
  • the image acquisition unit 100 acquires imaging parameter information of the planar images of the fundus of the eye when they are captured and stores the information in the storage unit 130 via the control unit 120 .
  • the imaging parameter information is, for example, position information of the fixation lamp when the imaging is performed.
  • the imaging parameter information such as the position information of the fixation lamp may be included in an information file which is added to the planar images of the fundus of the eye, but may also exist as tag information of the images.
  • step S 220 the input information acquisition unit 110 acquires information of the subject's eye from the database or from the input of the operator via an input unit (not illustrated).
  • the information of the subject's eye is information of the patient (e.g., patient ID, name, age, and sex), whether the examination target is right/left eye, and shooting date and time.
  • the acquired information is stored in the storage unit 130 via the control unit 120 .
  • step S 230 the frequency conversion unit 141 acquires a spatial frequency image by the discrete Fourier transform using the planar images of the fundus of the eye acquired by the adaptive optics SLO and stored in the storage unit 130 .
  • a greater part of each planar image of the fundus of the eye is populated with regularly-arranged photoreceptor cells observed as small regions with high luminance.
  • a Fourier image obtained from such a planar image of the fundus of the eye by spatial frequency conversion has a ring structure as illustrated in FIG. 4 .
  • FIG. 5B illustrates a function I(r) acquired from the Fourier image in FIG. 5A .
  • the function I(r) in FIG. 5B includes a lot of information regarding the arrangement of the photoreceptor cells. For example, when the crystalline lens of the subject's eye becomes cloudy due to a disease, the signals of the photoreceptor cells become extremely weak (see FIG. 6A ). In such a case, as illustrated in FIG. 6B , the value of the function I(r) itself becomes small compared to the function I(r) in FIG. 5B . Further, if the photoreceptor cells are partially absent in the imaging region and a periodic structure is not produced, a similar function I(r) is obtained. Thus, as a feature quantity indicating the intensity of the periodic structure of the photoreceptor cells, as illustrated in FIG. 7 , Imax (i.e., a maximum value of I(r)) and Isum (i.e., a sum of I(r)) can be used.
  • Imax i.e., a maximum value of I(r)
  • Isum i.e., a sum
  • I max max r ⁇ ⁇ I ⁇ ( r )
  • I sum ⁇ r ⁇ ⁇ I ⁇ ( r )
  • rmax which is a value of r of Imax, corresponds to the periodicity of the arrangement of the photoreceptor cells.
  • r max arg ⁇ ⁇ max r ⁇ ⁇ I ⁇ ( r )
  • rmax becomes greater. Conversely, if the photoreceptor cells are distant from the macula lutea and the density of the photoreceptor cells is low, rmax becomes smaller.
  • FIG. 8A shows a case where the periodic structure of the photoreceptor cells may not be clearly determined due to low resolution.
  • the shape of the function I(r) is broadened toward the smaller function r.
  • FIG. 8C illustrates the calculation regions.
  • FIG. 9A illustrates the imaging condition in which the imaging condition is poor and the resolution is low.
  • FIG. 9B illustrates the signals of I(r) in such a case.
  • Imax is acquired
  • accurate peak_ratio cannot be calculated.
  • a peak position which is assumed from the acquired visual fixation information (position r_expected_peak) will be used in place of rmax in calculating the peak_ratio. In this manner, the resolution of the image can be evaluated more accurately.
  • step S 250 from the feature quantity acquired in step S 240 , the index calculation unit 143 acquires information of the photoreceptor cells such as the information of the imaging state and the distribution information of the photoreceptor cells of the image captured by the adaptive optics SLO.
  • the acquired image quality index is stored in the storage unit 130 via the control unit 120 .
  • An example of a method for calculating the image quality index from the feature quantity acquired in step S 240 will be described below.
  • the index which is acquired is not limited to the calculation method described below.
  • the indices associated with the intensity of the periodic structure of the photoreceptor cell such as Imax and Isum and those associated with the resolution such as the peak_ratio are important in terms of image quality.
  • the adaptive optics SLO even if the health level of the subject's eye is good and high signal can be obtained, good resolution is not always obtained.
  • an image of good quality may be acquired if, for example, the aberration is corrected again using the wavefront sensor.
  • FIGS. 10A and 10B illustrate examples of the index.
  • the vertical axis represents Imax, which is an index associated with the intensity of the periodic structure and the horizontal axis represents the peak_ratio, which is an index associated with the resolution.
  • the peak_ratio which is an index associated with the resolution.
  • step S 240 a case where rmax is used for calculating the peak_ratio and another case where the position r_expected_peak is used for calculating the peak_ratio have been described.
  • the signal strength is reduced, the noise is increased, and it becomes difficult to acquire an accurate rmax.
  • FIGS. 10A and 10B a case where the peak_ratio is calculated using the r_expected_peak is presented.
  • an image in a region A which is a high index region regarding the intensity of the periodic structure and the resolution
  • An image in a region B which is a low index region regarding the intensity of the periodic structure and the resolution
  • an image in a region C which is a high index region regarding the intensity of the periodic structure but a low index region regarding the resolution
  • the state of the subject's eye is not bad but the resolution may be poor. If such an image is acquired, by performing the aberration correction again, the image quality may be improved.
  • Such an image may also include both a region where the structure of the photoreceptor cells is clear a region where the structure of the photoreceptor cells is not clear.
  • Such a portion may correspond to a lesion portion.
  • an image in a region D which is a low index region regarding the intensity of the periodic structure and a high index region regarding the resolution, has good resolution but weak signals in general. If an image is such, a clearer image may be acquired by increasing the signal strength.
  • step S 260 the output unit 150 acquires the information indicating the imaging state, such as the image quality index, stored in the storage unit in step S 250 and displays it on the display unit 160 .
  • the acquired feature quantity is displayed as it is as the information of the imaging state.
  • a threshold value of the image quality index is set in advance when the image index value is smaller than the threshold value
  • a message such as “low image quality of photoreceptor cell region” can be displayed on the display unit 160 by the output unit 150 as the information of the imaging state.
  • the information is not necessarily character information and an icon or an image corresponding to the character information may be displayed in place of the character information.
  • the inspector can easily notice that the image has a quality problem. If the image index value exceeds the threshold value, a message such as “good image quality” may be displayed. Then, the user can understand that the image has no quality problems.
  • the output unit 150 causes the display unit 160 to display a message indicating that the ophthalmologic imaging apparatus needs correction of aberration. According to this message, the user can understand that the correction of aberration is necessary. Further, if the value of the peak_ratio is high and the value of the Imax is low, a message indicating that an increase in imaging light quantity is necessary is displayed on the display unit 160 . In this manner, the user can understand that alight quantity adjustment is necessary.
  • the density information of the photoreceptor cells may be displayed based on the position of the rmax. Further, a character or a graphic that notifies the user that a blood vessel or a lesion region is extracted instead of the photoreceptor cells on the planar image of the fundus of the eye may be displayed.
  • the output unit 150 also causes the display unit 160 to display a frequency image, such as the one illustrated in FIG. 4 , together with the image of the fundus of the eye. Further, a graph, such as the one illustrated in FIG. 7 , may also be displayed. In this manner, detailed information of the imaging state can be presented to the user. Additionally, information of a distribution state of the photoreceptor cells such as density information may be displayed on the display unit. This information is based on the diameter of the ring structure extracted by the feature extraction unit 142 . Further, a message informing the user that the image includes a region not including photoreceptor cells (i.e., a blood vessel or a lesion portion) can be displayed on the display unit.
  • a region not including photoreceptor cells i.e., a blood vessel or a lesion portion
  • Feature quantity and other information stored in the storage unit 130 in steps S 210 to S 250 are stored in a database.
  • an index useful in objectively evaluating the image quality of a planar image of a fundus of the eye captured by the adaptive optics SLO apparatus can be presented.
  • an objective criterion for determination can be presented.
  • a Fourier image is obtained by a frequency conversion of the entire planar image of the fundus of the eye obtained by the adaptive optics SLO, and an index useful for evaluating the image quality of the entire planar image of the fundus of the eye is calculated by extracting various feature quantities from the obtained Fourier image.
  • an index useful for evaluating the image quality of the entire planar image of the fundus of the eye is calculated by extracting various feature quantities from the obtained Fourier image.
  • the planar image of the fundus of the eye is segmented into a plurality of local regions. Then, a Fourier image is acquired for each of the planar images of the fundus of the eye acquired by the segmentation, and the planar image is analyzed using the feature quantities extracted from the Fourier images.
  • FIG. 11 illustrates a functional configuration of the image processing apparatus 10 according to the present embodiment. Since the functional configurations of units other than the image processing unit 140 are similar to those illustrated in FIG. 1 , their descriptions are not repeated.
  • the image processing unit 140 further includes an image segmentation unit 1140 and an integration unit 1144 as well as the frequency conversion unit 141 , the feature extraction unit 142 , and the index calculation unit 143 . After the planar image of the fundus of the eye is segmented into a plurality of regions, a feature quantity of each region is extracted. The acquired feature quantities are integrated into an integrated index. Then, the image is evaluated according to the integrated index.
  • the image segmentation unit 1140 segments the planar image of the fundus of the eye into a plurality of partial regions.
  • the frequency conversion unit 141 performs the frequency conversion for each partial region and obtains a partial frequency image.
  • the feature extraction unit 142 extracts a feature quantity based on the ring structure that appears on the partial frequency image.
  • the index calculation unit 143 acquires the information of the photoreceptor cells in the plurality of partial regions.
  • the image can be segmented into a partial region where, for example, a blood vessel is extracted but photoreceptor cells are not extracted, and a partial region other than such a region. Extremely fine blood vessels and lesion portions buried in the noise and unable to determine from the feature quantity of the frequency image are not considered in the segmentation.
  • the calculation of the image quality index is performed with respect to the entire planar image of the fundus of the eye acquired by the adaptive optics SLO.
  • the planar image of the fundus of the eye is segmented into a plurality of local regions, and the index of each region is calculated. Then, the indices are combined and used for the evaluation of the entire image.
  • the image processed in steps S 230 , S 240 , and S 250 is a local planar image of the fundus of the eye obtained by the segmentation of the planar image of the fundus of the eye.
  • step S 1230 the image segmentation unit 1140 segments the planar image of the fundus of the eye acquired by the adaptive optics SLO and stored in the storage unit 130 into a plurality of local regions.
  • Various methods can be used for the segmentation. Although the local difference can be more noticeable if the planar image of the fundus of the eye is segmented into a great number of small regions, the accuracy of the information obtained from each local region will be reduced. Further, since the frequency conversion of a plurality of regions requires more time and cost, it is important to use an image size of data of 2 to the n-th power suitable for high speed Fourier transform.
  • a plurality of local planar images of the fundus of the eye with a size of 128 ⁇ 128 are acquired as described in FIG. 13 .
  • Each of the local planar images includes an overlapping portion.
  • the segmentation method is not limited to such an example.
  • Sixteen local planar images of the fundus of the eye generated in this manner are stored in the storage unit 130 via the control unit 120 .
  • the processing in steps S 230 , S 240 , and S 250 is similar to the processing in the first exemplary embodiment, it is different in that the processing is performed for each of sixteen local planar images of the fundus of the eye generated in step 1230 . Accordingly, sixteen indices are obtained. Further, each of the obtained indices and the feature quantities acquired in the course of process is stored in the storage unit 130 in association with the corresponding local planar image of the fundus of the eye.
  • step S 1260 the integration unit 1144 generates an index by integrating the indices acquired from each local planar image of the fundus of the eye. Then, the obtained index is displayed on the monitor via the output unit 150 . Further, the feature quantity and other information stored in the storage unit 130 in steps S 210 to S 1260 are stored in the database.
  • FIGS. 14A and 14B illustrate examples of sixteen indices acquired in step S 250 .
  • all of the sixteen indices show poor image quality. If the state of the subject's eye is not good and, for example, the crystalline lens is very cloudy, very few signals can be obtained. In such a case, all the indices acquired from the local planar images of the fundus of the eye show bad values. Thus, if indices such as those illustrated in FIG. 14A are obtained, the quality of the planar image of the fundus of the eye is considered as low as a whole.
  • the indices of both high values and of low values are determined as illustrated in FIG. 14B , it can be considered that some of the local regions failed in the acquisition of the periodic structure of the photoreceptor cells although the imaging itself has been successful. Such a case indicates that the photoreceptor cells are not extracted due to, for example, the existence of a blood vessel, or the photoreceptor cells do not exist due to, for example, a lesion.
  • planar image can be evaluated by a combination of the first and the second exemplary embodiments, for example, a combination of the overall evaluation and the partial evaluation.
  • the evaluation will be useful since the user can determine the reason for the poor evaluation (e.g., poor imaging condition, lesion portion was included). This is useful as if any part of the partial evaluation indicates high image quality, the possibility of the imaging condition being the cause of the poor image quality is low.
  • an optimum index can be presented by combining the evaluation methods.
  • the planar image of the fundus of the eye acquired by the adaptive optics SLO apparatus is segmented into a plurality of local regions. Further, the image quality indices obtained from the local planar images of the fundus of the eye are integrated to form an integrated index. According to this index, a case where the overall image quality of the fundus of the eye is low and a case where the image quality itself is not low but a region where the periodic structure of the photoreceptor cells is not extracted exists due to the presence of a blood vessel or a lesion region can be presented.
  • the ophthalmologic imaging apparatus can be configured such, from the aspect of hardware, that the index calculation unit 143 automatically determines a lesion region. Then, based on the result of the determination, the output unit 150 outputs information regarding the imaging region to the apparatus. On receiving the information, the apparatus automatically captures the image again by setting the center of the image to that imaging region.
  • the output unit 150 outputs an aberration correction instruction or a control value of an adjustment value of a light quantity according to the imaging state to the ophthalmologic imaging apparatus, then the ophthalmologic imaging apparatus controls the aberration measurement unit and the correction optical system based on the input control value. Accordingly, the user does not need to adjust the apparatus.
  • the information of the imaging state is displayed and the user performs the operation according to the displayed information.
  • the information of the imaging state is used for forming and selecting an image. Since the configurations of the apparatus are similar to those described in the first and the second exemplary embodiments, their descriptions are not repeated and the points different from the above-described exemplary embodiments will be mainly described.
  • the image acquisition unit 100 acquires a plurality of planar images of the fundus of the eye obtained by imaging of different or substantially same positions. Since the fundus of the eye moves due to the involuntary eye movement such as a saccade, even if a tracking function is provided, the tracking may not be successful. Thus, even if the apparatus is set to capture an image of the same position, actually, an image of a different position may be captured. This is why the expression “substantially the same” is used.
  • the image processing unit 140 functions as a determination unit configured to determine the image quality.
  • the image processing unit 140 acquires the information of the imaging state acquired by the index calculation unit 143 and the information of the threshold value of the imaging state determined in advance. Then, the image processing unit 140 determines whether the imaging state obtained for each image exceeds of falls below the threshold value.
  • the information of the imaging state for example, the values of the peak_ratio and Imax described in the first exemplary embodiment or an image quality index value obtained by a combination of such values can be used.
  • the image processing unit 140 functions as a selection unit configured to select an image. If an image is determined by the determination unit as having an image quality below the threshold value, the image processing unit 140 determines that it is inappropriate for the presentation to the user or not suitable for the superimposing processing and does not select the image. In other words, the image processing unit 140 selects only the images that exceed the threshold value. Since the images below the threshold value are not selected, images with aberration, images with insufficient light quantity, and images captured when the eye instantaneously moved due to blinking or the like are not selected by the image processing unit 140 .
  • the image processing unit 140 further aligns the selected images and performs the superimposing processing.
  • the image processing unit 140 also functions as a generation unit configured to generate an image with reduced random noise. Accordingly, the image quality of a plurality of images can be automatically determined and a planar image of a fundus of the eye with good image quality can be obtained.
  • the control unit 120 outputs only the image selected from the acquired plurality of images to the output unit 150 .
  • the control unit 120 outputs such an image together with the generation unit or in place of the generation unit.
  • the output unit 150 causes the display unit 160 to display only the selected image.
  • the unselected images are not transmitted to the output unit 150 but is stored in the storage unit 130 .
  • the unselected images can be deleted as a failed image with the approval of the user.
  • the output unit 150 can separately display the images which have not been selected on the display unit 160 and further display a window that accepts the approval of the user regarding the deletion.
  • the control unit 120 can delete the approved images from the storage unit 130 in one operation. Further, in place of the deletion, the control unit 120 can cause only the approved images not to be transferred to the external database. In this manner, the user can automatically determine unnecessary images and can easily perform the deletion processing or the non-transfer processing.
  • the function of the above-described image processing apparatus 10 can be realized by the hardware illustrated in FIG. 15 and using the software and hardware together.
  • the image processing apparatus illustrated in FIG. 15 includes a central processing unit (CPU) 1501 , a random access memory (RAM) 1502 , a read only memory (ROM) 1503 , an external storage unit 1504 , and a communication interface 1505 . These units are connected to each other via a bus 1509 . Further, a monitor 1506 , a keyboard 1507 , and a mouse 1508 are connected to the image processing apparatus 10 . A program including an instruction that can cause the processing illustrated in the flowchart of FIG. 2 or 12 is stored in the ROM 1503 or the external storage unit 1504 .
  • the function of the image processing apparatus 10 of the above-described exemplary embodiment is realized by the CPU 1501 reading out a stored program, loading it into the RAM 1502 , and executing the instruction included in the program.
  • the image processing apparatus 10 in the imaging system is connected to an ophthalmologic imaging apparatus 30 and a database 50 via a local area network (LAN) 40 .
  • the image processing apparatus 10 can be wirelessly connected to the LAN 40 .
  • the display unit 160 is connected to the image processing apparatus 10 .
  • the configuration of the imaging system is not limited to the above-described example.
  • an ophthalmologic imaging apparatus which includes an ophthalmologic imaging unit, an image processing unit, and the display unit 160 may be used.
  • the ophthalmologic imaging unit includes the function of the ophthalmologic imaging apparatus 30 and the image processing unit includes the function of the image processing apparatus 10 . If the imaging system has such a configuration, the quality of the image of the photoreceptor cells can be determined by the ophthalmologic imaging apparatus 30 alone, and the system can be downsized.
  • a composite apparatus composed of an SLO apparatus and an OCT apparatus according to the embodiments of the present invention is described as an ophthalmologic imaging apparatus.
  • the ophthalmologic imaging apparatus 30 which includes an adaptive optical system and can capture both planar images (SLO images) of high lateral resolution as well as tomographic images (OCT images) for a retina and acquires a planar image of a fundus of the eye will be described.
  • the ophthalmologic imaging apparatus 30 includes an SLO apparatus and an OCT apparatus.
  • the SLO apparatus acquires a planar image of a fundus of the eye by correcting the optical aberration of the subject's eye and using a spatial light modulator.
  • the OCT apparatus employs Fourier domain imaging when it acquires a tomographic image. According to these apparatuses, the ophthalmologic imaging apparatus 30 can obtain a good planar image of the fundus of the eye and a tomographic image of the subject's eye regardless of the visibility and the optical aberration.
  • a light source 201 is divided into a reference beam 205 and a measuring beam 206 by an optical coupler 231 .
  • the measuring beam 206 is guided to a subject's eye 207 as an object to be observed via a single mode fiber 230 - 4 , a spatial light modulator 259 , an XY scanner 219 , an X scanner 221 , and spherical mirrors 260 - 1 to 260 - 9 .
  • the measuring beam 206 is reflected or scattered by the subject's eye 207 being an object to be observed and returned as a return beam 208 .
  • the return beam 208 enters a detector 238 or a line sensor 239 .
  • the detector 238 converts the light intensity of the return beams 208 into a voltage signal.
  • a planar image of the fundus of the eye of the subject's eye 207 is generated.
  • the return beam 208 is combined with the reference beam 205 , and the combined light is caused to enter the line sensor 239 . Accordingly, a tomographic image of the subject's eye 207 is formed. Furthermore, by using a plurality of acquired tomographic images, a three-dimensional course of blood vessels can be extracted.
  • a spatial light modulator is used as a device for correcting the wavefront aberration
  • any device can be used so long as the wavefront aberration can be corrected.
  • a variable shape mirror can be used.
  • the light source 201 uses a super luminescent diode (SLD) being a typical low-coherent light source.
  • SLD super luminescent diode
  • the central wavelength and the bandwidth are 830 nm and 50 nm, respectively.
  • a low-coherent light source is selected.
  • the SLD is selected as a type of the light source, a different light source can be used so long as low coherent light can be emitted.
  • an amplified spontaneous emission (ASE) light source can be used.
  • near-infrared light is suitable as a wavelength from the viewpoint of measurement of eyes. Further, since the wavelength affects the lateral resolution of the acquired planar image of the fundus of the eye, a shorter wavelength is desirable. In the descriptions below, a wavelength of 830 nm is used. Other wavelengths may be selected depending on the measuring portion of the object to be observed.
  • the composite apparatus 30 further includes polarization controllers 253 - 1 and 253 - 2 .
  • the reference beam 205 divided by the optical coupler 231 is adjusted so that it is guided to a lens 235 - 1 via a single mode fiber 230 - 2 and adjusted so as to become parallel light with a beam diameter of 4 mm.
  • the reference beam 205 is guided to a mirror 214 being a reference mirror by mirrors 257 - 1 to 257 - 4 . Since the optical path length of the reference beam 205 is adjusted to be substantially the same as that of the measuring beam 206 , the reference beam 205 and the measuring beam 206 interfere with each other.
  • the reference beam 205 is reflected by the mirror 214 and guided again to the optical coupler 231 .
  • a dispersion compensation glass 215 through which the reference beam 205 passes, compensates the dispersion that occurs when the measuring beam 206 travels to the subject's eye 207 and returns from the subject's eye 207 with respect to the reference beam 205 .
  • a motorized stage 217 - 1 is movable in the directions indicated by an arrow to allow the optical path length of the reference beam 205 to be adjusted and controlled. Further, the motorized stage 217 - 1 is driven by a motorized stage driver 283 in a driver unit 281 under the control of a personal computer 225 .
  • the measuring beam 206 split by the optical coupler 231 is guided to a lens 235 - 4 via the single mode fiber 230 - 4 and adjusted so as to become parallel light with a beam diameter of 4 mm.
  • the polarization controller 253 - 1 or 253 - 2 can adjust the polarization state of the measuring beam 206 .
  • the polarization controller 253 - 1 or 253 - 2 adjusts the polarization state of the measuring beam 206 to be linearly polarized in a direction parallel to the drawing surface of FIG. 17 .
  • the measuring beam 206 passes through a beam splitter 258 and a movable beam splitter 261 (also referred to as a splitting unit) and reaches the spatial light modulator 259 via the spherical mirrors 260 - 1 and 260 - 2 . Then, the measuring beam 206 is modulated at the spatial light modulator 259 .
  • the spatial light modulator 259 utilizes the orientation characteristics of the liquid crystal. More specifically, the spatial optical modulator 259 is arranged in a direction where the spatial optical modulator 259 can modulate the phase of linear polarization parallel to the drawing surface of FIG. 17 (i.e., the P polarization) to coincide with the polarization orientation of the measuring beam 206 .
  • the measuring beam 206 passes through a polarizing plate 273 and reaches a mirror of the X scanner 221 via the spherical mirrors 260 - 3 and 260 - 4 .
  • the polarizing plate 273 has a role to guide only the linear polarization parallel to the drawing surface of FIG. 17 , of the return beams 208 , to the spatial light modulator 259 .
  • the X scanner 221 scans the measuring beam 206 in a direction parallel to the drawing surface of FIG. 17 .
  • the X scanner 221 is a resonant scanner having a drive frequency of approximately 7.9 kHz.
  • the measuring beam 206 is incident on a mirror of the XY scanner 219 via the spherical mirrors 260 - 5 and 260 - 6 .
  • the XY scanner 219 includes two mirrors (X-scanning mirror and Y-scanning mirror) arranged close to each other.
  • the measuring beam 206 is adjusted in such a manner that its center coincides with the center of rotation of the mirror of the XY scanner 219 .
  • the drive frequency of the XY scanner 219 is changeable in the range up to 500 Hz.
  • the spherical mirrors 260 - 7 to 260 - 9 are optical systems that cause the measuring beam 203 to scan a retina 227 . Having a point near a cornea 226 as a support point, the measuring beam 206 scans the retina 227 . Although the measuring beam 206 has a beam diameter of 4 mm, the beam diameter may be increased for acquisition of a high resolution tomographic image.
  • a motorized stage 217 - 2 is configured to move in a direction indicated by an arrow in FIG. 17 and to adjust and control the position of a spherical mirror 260 - 8 attached thereto. Similarly to the motorized stage 217 - 1 , the motorized stage 217 - 2 is controlled by the motorized stage driver 283 .
  • the position of the spherical mirror 260 - 8 By adjusting the position of the spherical mirror 260 - 8 , it is possible to focus the measuring beam 206 on a predetermined layer of the retina 227 of the subject's eye 207 and observe the subject's eye 207 . In the initial state, the position of the spherical mirror 260 - 8 is adjusted so that the measuring beam 206 can enter the cornea 226 as parallel light.
  • the ophthalmologic imaging apparatus 30 according to the present embodiment can cope with the subject's eye 207 having a refractive error.
  • the measuring beam 206 When the measuring beam 206 enters the subject's eye 207 , the measuring beam 206 becomes the return beam 208 due to reflection or scattering from the retina 227 . Then, the return beam 208 reaches the line sensor 239 by being guided again by the optical coupler 231 . A part of the return beam 208 is reflected by the movable beam splitter 261 and guided to the detector 238 via a lens 235 - 5 .
  • a light blocking plate 272 having a pinhole has a function of blocking unnecessary light not focused on the retina 227 among the return beams 208 . Further, the light blocking plate 272 is arranged in a position conjugate with the focusing position of the lens 235 - 5 . The diameter of the pinhole of the light blocking plate 272 is, for example, 50 ⁇ m.
  • the detector 238 for example, an avalanche photodiode (APD), which is a high-speed, highly-sensitive light sensor, is used.
  • a part of the return beam 208 split by the beam splitter 258 enters a wavefront sensor 255 .
  • the wavefront sensor 255 is a Shack-Hartmann wavefront sensor.
  • the spherical mirrors 260 - 1 to 260 - 9 are arranged so that they are optically conjugate with the XY scanner 219 , the X scanner 221 , the cornea 226 , the wavefront sensor 255 , and the spatial light modulator 259 .
  • the wavefront sensor 255 can measure the aberration of the subject's eye 207 .
  • the spatial light modulator 259 can correct the aberration of the subject's eye 207 .
  • by controlling the spatial light modulator 259 in real time based on the obtained aberration it is possible to correct the aberration that occurs in the subject's eye 207 and to acquire tomographic images with higher lateral resolution.
  • the ophthalmologic imaging apparatus 30 can acquire tomographic images (OCT images) as well as planar images of the fundus of the eye (SLO images).
  • the return beams 208 are combined by the optical coupler 231 .
  • the combined light (combined light 242 ) is guided to a transmissive grating 241 via a single mode fiber 230 - 3 and a lens 235 - 2 , and is dispersed for each wavelength. Then, the light enters the line sensor 239 via a lens 235 - 3 .
  • the line sensor 239 converts the light intensity for each position (wavelength) into a voltage signal.
  • the voltage signal is converted into a digital value by a frame grabber 240 so that tomographic images of the subject's eye 207 are formed in the personal computer 225 .
  • the line sensor 239 includes 1024 pixels and can obtain the intensity of the combined light 242 for each wavelength (a segmentation of 1024).
  • a part of the return beam 208 is reflected by the movable beam splitter 261 . After unnecessary light is blocked by the light blocking plate 272 , the reflected light reaches the detector 238 , and the light intensity is converted into an electric signal. Data of the obtained electric signal is processed by the personal computer 225 in synchronization with the scanning signal of the X scanner 221 and the XY scanner 219 . Accordingly, a planar image of the fundus of the eye is formed.
  • An image signal obtained by the wavefront sensor 255 is supplied to the personal computer 225 , and the aberration is calculated.
  • the obtained aberration is expressed by Zernike polynomial, which represents the aberration of the subject's eye 207 .
  • the Zernike polynomial includes a tilt term, a defocus term, an astigmatism term, a coma term, and a trefoil term.
  • the ophthalmologic imaging apparatus 30 acquires a tomographic image of the retina 227 by controlling the XY scanner 219 and acquiring an interference fringe by the line sensor 239 using the X scanner 221 as a fixed mirror.
  • the movable beam splitter 261 is controlled so that the return beams 208 are not guided to the detector 238 .
  • the X scanner 221 and the XY scanner 219 are controlled by an optical scanner driver 282 in the driver unit 281 from the personal computer 225 .
  • the method for acquiring a tomographic image (a plane parallel to the optical axis) of the retina 227 will be described.
  • FIG. 18A is a schematic view of the subject's eye 207 and illustrates a state where the subject's eye 207 is observed by the ophthalmologic imaging apparatus 30 .
  • the measuring beam 206 enters the retina 227 via the cornea 226 , the measuring beam 206 becomes the return beams 208 due to reflection or scattering at various positions. Then, the return beam 208 reaches the line sensor 239 with delays at the respective positions.
  • the line sensor 239 can detect an interference fringe when the optical path length of the reference beam path is approximately equal to the optical path length of the measuring beam path.
  • the interference fringe acquired by the line sensor 239 is an interference fringe in a spectrum region on the wavelength axis.
  • the interference fringe which is information of the wavelength axis, is converted into an interference fringe on an optical frequency axis considering the characteristics of the line sensor 239 and the transmissive grating 241 . Further, by performing inverse Fourier transform of the obtained interference fringe on the optical frequency axis, it is possible to obtain information in the depth direction. As illustrated in FIG. 18B , if the interference fringe is detected while the XY scanner 219 is driven, the interference fringe can be obtained at each position on the X-axis. More specifically, the information in the depth direction can be obtained at each position on the X-axis. As a result, a two-dimensional distribution of the intensity of the return beams 208 on the XZ plane is obtained, which is a tomographic image 232 (see FIG. 18C ).
  • the tomographic image 232 is an image obtained by arranging the intensity of each of the return beams 208 in an array.
  • the tomographic image 232 is a gray scale intensity image.
  • the length of the tomographic image 232 in the X direction is 700 ⁇ m, which is similar to the SLO image described below.
  • FIG. 18C only the boundaries in the obtained tomographic image are expressed in lines.
  • the illustration includes a retinal pigment epithelial layer 246 , a stratum opticum 247 , and a blood vessel 278 . Further, if a plurality of tomographic images is acquired at each position on the Y-axis, it is possible to visualize a three-dimensional course of the blood vessel.
  • the ophthalmologic imaging apparatus 30 can acquire a planar image of the fundus of the eye of the retina 227 by controlling and operating the XY scanner 219 in only the Y-axis direction and the X scanner 221 while fixing the X-axis direction of the XY scanner 219 , and acquiring the intensity of the return beam 208 using the detector 238 .
  • the X scanner 221 and the XY scanner 219 are controlled by the optical scanner driver 282 in the driver unit 281 from the personal computer 225 (see FIG. 17 ). Further, the ophthalmologic imaging apparatus 30 can acquire a planar image of the fundus of the eye while correcting the aberration that occurs at the subject's eye 207 by controlling the spatial light modulator 259 .
  • the spatial light modulator 259 can be controlled by using the aberration of the subject's eye 207 measured by the wavefront sensor 255 . Furthermore, the ophthalmologic imaging apparatus 30 can acquire a planar image of the fundus of the eye while controlling the spatial light modulator 259 in real time.
  • the ophthalmologic imaging apparatus 30 acquires a planar image of a fundus of the eye of the retina 227 by controlling the XY scanner 219 and acquiring the intensity of the return beam 208 by the detector 238 .
  • An acquisition method of a planar image of the fundus of the eye of the retina 227 (planar image perpendicular to the optical axis) will now be described.
  • FIG. 19A is a schematic view of the subject's eye 207 and illustrates a state where the subject's eye 207 is observed by the ophthalmologic imaging apparatus 30 .
  • the measuring beam 206 when the measuring beam 206 enters the retina 227 via the cornea 226 , the measuring beam 206 becomes the return beam 208 due to reflection or scattering at various positions. Then, the return beam 208 reaches the detector 238 . Further, as illustrated in FIG. 19B , if the intensity of the return beam 208 is detected while the XY scanner 209 is moved in the X direction, information at each position on the X-axis can be obtained.
  • the XY scanner 209 is simultaneously moved in both the X-axis and Y-axis directions with respect to an image capturing range 292 where the retina 227 is present. Then, the measuring beam 206 is raster-scanned along a trajectory 293 on the image capturing range 292 . In this state, if the intensity of the return beams 208 is detected, a two-dimensional distribution of the intensity of the return beams 208 can be obtained. Accordingly, a planar image of the fundus of the eye 277 (see FIG. 19D ) is acquired.
  • the measuring beam 206 is scanned from a point S at the upper right to a point E at the bottom left. While the scanning is performed, the intensity of the return beams 208 is used in forming the planar image of the fundus of the eye 277 .
  • the trajectory 293 from the point E to the point S is the movement of the measuring beam 206 as a preparation for the imaging of the planar image of the fundus of the eye 277 which is to be formed next.
  • the time required for the scanning takes 84% for the point S to the point E and 16% for the point E to the point S with respect to the trajectory 293 in FIG. 19C . This ratio is based on the duty ratio of the drive waveform of the Y scanner described above. Further, in FIG. 19C , for the sake of simplicity, the number of times of the scanning in the X direction with respect to the trajectory 293 is smaller than the actual number of times.
  • the planar image of the fundus of the eye 277 has a size of 700 ⁇ 350 ⁇ m.
  • the time required for the acquisition is approximately 15.6 ms. This time is based on the drive frequency of the Y scanner.
  • a photoreceptor cell group 279 where the intensity of the return beam 208 is relatively high, is light colored, whereas the blood vessel 278 , where the intensity is relatively low, is dark colored.
  • blood cells (not illustrated) are light colored in the blood vessel 278 . If the planar image of the fundus of the eye 277 is continuously acquired, the movement of the blood cells through the blood vessel 278 can be visualized. Further, spatio-temporal images may be generated by extracting the blood vessel 278 , from which the blood cells are visualized, from the planar images 277 of the fundus of the eye which are continuously acquired, and superimposing the extracted planar images 277 of the fundus of the eye in the order they have been captured. Movement of the blood cells and the blood speed can be easily obtained.
  • An imaging method for obtaining an image useful for obtaining the blood speed and an imaging method for obtaining an image of the photoreceptor cells are designated by the user via an operation unit of the image processing apparatus 10 .
  • a designation unit (not illustrated) of the image processing apparatus designates the selected imaging method and transmits associated information to the ophthalmologic imaging apparatus 30 .
  • the ophthalmologic imaging apparatus 30 outputs an image of photoreceptor cells which is easy to observe or an image having a bloodstream as an object to be observed.
  • the embodiment can also be realized by a computer of a system or apparatus (or devices such as a CPU or MPU) that reads out and executes a program recorded on a memory device to perform the functions of the above-described embodiment(s), and by a method, the steps of which are performed by a computer of a system or apparatus by, for example, reading out and executing a program recorded on a memory device to perform the functions of the above-described embodiment(s).
  • the program is provided to the computer for example via a network or from a recording medium of various types serving as the memory device (e.g., computer-readable storage medium).

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Ophthalmology & Optometry (AREA)
  • Biomedical Technology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Quality & Reliability (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Signal Processing (AREA)
  • Eye Examination Apparatus (AREA)

Abstract

An image processing apparatus includes an identification unit configured to identify periodicity of a fundus image obtained by capturing an image of a fundus of an eye, and an information acquisition unit configured to acquire information indicating an imaging state of photoreceptor cells in the fundus image based on the periodicity.

Description

    CROSS REFERENCE OF RELATED APPLICATIONS
  • This application is a Continuation of U.S. patent application Ser. No. 13/605,873 filed on Sep. 6, 2012 which claims the benefit of Japanese Application No. 2011-204653 filed Sep. 20, 2011, which is hereby incorporated by reference herein in its entirety.
  • BACKGROUND
  • Examination of a fundus of the eye is widely accepted as important in the early diagnosis of lifestyle-related diseases and diseases likely to cause blindness. Fundus cameras and scanning laser ophthalmoscopes (SLO) are among the apparatuses used for the inspection of the fundus of the eye. The fundus camera captures an image of a fundus of the eye by receiving reflected light of a light beam which has entered the fundus of the eye. The SLO is an ophthalmologic apparatus that uses the principle of confocal laser scanning microscope. In recent years, fundus cameras and SLOs including an adaptive optical system have been developed and are used for acquiring fundus planar images of high lateral resolution. The adaptive optical system measures an aberration of a subject's eye by a wavefront sensor in real time and corrects the aberration of the measuring beam and the return beam that occurs at the subject's eye by a wavefront correcting device. Further, attempts are being made to capture images of photoreceptor cells of a retina using by these apparatuses and make a diagnosis of a disease or evaluate drug response.
  • As an example of visualization of the photoreceptor cells using an adaptive optics SLO, Kaccie Y. Li and Austin Roorda, “Automated identification of cone photoreceptors in adaptive optics retinal images” J. Opt. Soc. Am. A, May 2007, Vol. 24, No. 5, 1358 discusses an ophthalmologic imaging apparatus which is capable of automated extraction of photoreceptor cells by acquiring a planar image of the fundus of the eye regarding the retina. According to this technique, a fundus planar image of a retina with high lateral resolution is acquired by preprocessing the acquired planar image, in other words, removing high frequency components from the planar image using periodicity of the arrangement of the photoreceptor cells visualized in the image.
  • SUMMARY
  • According to some embodiments of the present invention, an image processing apparatus includes an identification unit configured to identify periodicity of a fundus image obtained by capturing an image of a fundus of an eye, and an information acquisition unit configured to acquire information indicating an imaging state of photoreceptor cells in the fundus image based on the periodicity.
  • Further features and aspects of the present disclosure will become apparent from the following detailed description of exemplary embodiments with reference to the attached drawings.
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate exemplary embodiments, features, and aspects of the invention and, together with the description, serve to explain the principles of the embodiments of the present invention.
  • FIG. 1 illustrates a functional configuration of an image processing apparatus according to a first exemplary embodiment.
  • FIG. 2 is a flowchart illustrating processing procedures of the image processing apparatus according to the first exemplary embodiment.
  • FIG. 3 is a schematic view of a high-precision planar image of a fundus of the eye including photoreceptor cells acquired by an adaptive optics SLO apparatus.
  • FIG. 4 illustrates an example of a Fourier image obtained by frequency conversion of the planar image of the fundus of the eye.
  • FIGS. 5A and 5B illustrate a method for calculating a structure that reflects an arrangement of the photoreceptor cells from a Fourier image and a graph organizing the result of the calculation.
  • FIGS. 6A and 6B illustrate an example of a Fourier image obtained by frequency conversion of the planar image of the fundus of the eye in a case where the signal is weak and a graph indicating the arrangement of the photoreceptor cells.
  • FIG. 7 illustrates a feature quantity acquired from a Fourier image.
  • FIGS. 8A, 8B, and 8C illustrate an example of a Fourier image of low resolution and feature quantities acquired from the image.
  • FIGS. 9A, 9B, and 9C illustrate an example of a Fourier image of even lower resolution and feature quantities acquired from the image.
  • FIGS. 10A and 10B illustrate an image quality index acquired from a feature amount extracted from a Fourier image.
  • FIG. 11 is a functional configuration of the image processing apparatus according to a second exemplary embodiment.
  • FIG. 12 is a flowchart illustrating processing procedures of the image processing apparatus according to the second exemplary embodiment.
  • FIG. 13 illustrates an example of dividing the planar image of the fundus of the eye into a plurality of local planar images of the fundus of the eye.
  • FIGS. 14A and 14B illustrate examples of indices acquired from the local planar images of the fundus of the eye.
  • FIG. 15 illustrates a hardware configuration of the image processing apparatus according to another exemplary embodiment.
  • FIG. 16 illustrates a configuration of an ophthalmologic imaging system.
  • FIG. 17 illustrates a configuration of an ophthalmologic imaging apparatus.
  • FIGS. 18A, 18B, and 18C illustrate an image acquisition method for optical coherence tomography (OCT).
  • FIGS. 19A, 19B, 19C, and 19D illustrate an image acquisition method for SLO.
  • DESCRIPTION OF THE EMBODIMENTS
  • In the following detailed description, for purposes of explanation, numerous specific details are set forth in order to provide a thorough understanding of the disclosed embodiments. It will be apparent, however, that one or more embodiments may be practiced without these specific details. In other instances, well-known structures and devices are schematically shown in order to simplify the drawing.
  • When images of photoreceptor cells are captured, it is useful to evaluate the imaging state of the image of the photoreceptor cells since such information helps the adjustment of the imaging conditions and the diagnosis of the obtained image. Generally, in evaluating an imaging state, an image quality index value is obtained according to a comparison of a noise level of a region being an evaluation target and a noise level of a different region. However, since the imaging area of the image of the photoreceptor cells is generally the retina region, the comparison method cannot be used for the evaluation.
  • In one embodiment of the present invention, an index that objectively evaluates an image quality of an image of photoreceptor cells of a retina acquired by an ophthalmologic apparatus such as an adaptive optics SLO is presented.
  • According to an exemplary embodiment of the present invention, since information of the imaging state of the photoreceptor cells can be obtained from processing of an image of the fundus of the eye, the adjustment of the imaging conditions and the diagnosis of the image can be easily performed.
  • According to a first exemplary embodiment, when an image is obtained by imaging of photoreceptor cells of a retina by an adaptive optics SLO, an index that quantitatively indicates the image quality of the acquired image is calculated and presented. More specifically, a spatial frequency image is acquired by discrete Fourier transform using a planar image of a fundus of the eye which has been acquired by the adaptive optics SLO. The acquired spatial frequency image is hereinafter referred to as a Fourier image or a frequency image. Then, a feature quantity of a periodic structure that reflects the regular arrangement of the photoreceptor cells is extracted from the acquired Fourier image, and the index of the image quality is acquired from the extracted feature quantity.
  • If the image quality index acquired in such a manner is presented when a user is capturing an image of the fundus of the subject's eye, the user can determine whether the image needs to be captured again. Further, if a diagnosis is to be performed based on the density of the photoreceptor cells, the user can determine whether it is adequate to use such an image for the diagnosis.
  • FIG. 3 schematically illustrates a planar image of the fundus of the eye captured by the adaptive optics SLO. As illustrated in FIG. 3, a small region having relatively high luminance can be extracted in a distinguishable manner as a photoreceptor cell PR. Further, a blood vessel region V having low luminance compared to the luminance of the photoreceptor cell may be extracted. The blood vessel region V is the shadow of the blood vessel above the photoreceptor cells.
  • FIG. 4 illustrates an example of a Fourier image acquired by the discrete Fourier transform of spatial frequency components of the planar image of the fundus of the eye described above. As illustrated in FIG. 4, a ring that corresponds to the periodicity of the photoreceptor cells is formed according to the periodic arrangement of the photoreceptor cells.
  • FIG. 1 illustrates a functional configuration of an image processing apparatus 10 according to the present embodiment.
  • An image acquisition unit 100 in FIG. 1 acquires a planar image of a fundus of the eye. The acquired planar image is an image of photoreceptor cells of the fundus of the eye. For example, it is an image acquired by a fundus imaging apparatus which includes an aberration measurement unit such as a Hartmann-Shack wavefront sensor and an adaptive optical system which corrects the aberration. The adaptive optical system includes, for example, a reflective display panel having liquid crystals arranged on a silicon substrate (liquid crystal on silicon (LCOS)) or a deformable mirror.
  • When the planar image of the fundus of the eye is captured by the adaptive optics SLO, an input information acquisition unit 110 acquires information of the subject's eye which is being captured. The acquired image is stored in a storage unit 130 via a control unit 120. An image processing unit 140 includes a frequency conversion unit 141, a feature extraction unit 142, and an index calculation unit 143.
  • The frequency conversion unit 141 performs frequency conversion of the planar image of the fundus of the eye and obtains a frequency image. For example, the discrete Fourier transform is used for converting frequencies, as described above.
  • The feature extraction unit 142 extracts the feature quantity indicating the periodic structure of the planar image of the fundus of the eye from the frequency image. In the frequency image, the periodic structure of the planar image of the fundus of the eye is, for example, appears in a ring structure having a point of origin at the center. The ring of the frequency image presents a specified frequency band corresponding to the structure of the photoreceptor cells.
  • The feature quantity of the periodic structure is, for example, the feature quantity concerning this ring structure such as a peak luminance value in the ring structure and a sharpness value that indicates the extension of the ring structure. Further, there is the peak position of the luminance value in the ring structure. Additionally, if the image quality is poor, a disk-like structure may appear in the frequency image rather than the ring structure. In this case, the extension and the position of the disk region, and the peak luminance value in the disk region are extracted as the feature quantities. These are used for acquiring information of the imaging state and the distribution information of the photoreceptor cells in the planar image of the fundus of the eye.
  • The periodicity can be identified using information other than the frequency image. For example, the peak value of the luminance near the center of the photoreceptor cells is detected and an average of the profile of luminance distribution around each detection point is acquired. In this manner, the periodicity of the image of the fundus of the eye can be identified.
  • The peak value can be detected using maximum value detection processing. The profile in this context is in the radial direction having the detection point at the center, and the mean value is taken in the angular direction of the polar coordinates having the detection point at the center. Then, the feature extraction unit 142 evaluates the shape. If the shape is periodic, a peak corresponding to a detection point adjacent to the detection point appears on the average profile.
  • Additionally, as a method for identifying the periodicity, co-occurrence matrix and fractal dimension known as the texture feature quantity can be used.
  • The identification of the periodicity by the Fourier transform is one exemplary embodiment, and the above-described different methods can be used in identifying the periodicity. In this regard, the frequency conversion unit 141 functions as one exemplary embodiment of an identification unit of the periodicity.
  • The index calculation unit 143 functions as an information acquisition unit configured to acquire information of the imaging state of the photoreceptor cells and distribution information of the photoreceptor cells from the image indicating the periodicity such as a Fourier image. By using the feature quantity extracted by the feature extraction unit 142, the index calculation unit 143 can obtain accurate calculation results of the imaging state and the distribution of the photoreceptor cells.
  • An output unit 150 outputs the information of the imaging state and the distribution information of the photoreceptor cells acquired by the index calculation unit 143. The information is output to a display unit 160 as well as an external database and an output apparatus. The information is, for example, printed out by the output apparatus. Information of the imaging state can be displayed as it is or together with the planar image of the fundus of the eye on the output apparatus.
  • As the information of the imaging state which is displayed, there is the image quality index value of the planar image of the fundus of the eye where the photoreceptor cells have been captured, and information asking the user to perform operation to improve the imaging state. Further, a frequency image including the information of the imaging state can be displayed together with the image of the fundus of the eye. Furthermore, as illustrated in FIG. 5B described below, a graph showing a relation between a distance from the center position of the ring structure or from the point of origin of the frequency image and the luminance value is displayed. Additionally, a control value corresponding to the imaging state of the photoreceptor cells is output to the ophthalmologic imaging apparatus which has captured the planar image of the fundus of the eye.
  • The image processing unit 140 generates a Fourier image from the acquired planar image of the fundus of the eye, calculates an index of the image quality from the feature quantity extracted from the Fourier image, and stores the index in the storage unit 13. The output unit 150 outputs the calculated index to a monitor. Further, the output unit 150 outputs the result of the processing stored in the storage unit 130 to the database.
  • Although the image processing apparatus 10 acquires the planar image of the fundus of the eye directly from the adaptive optics SLO in the present embodiment, the image can be acquired via a network. In such a case, a plurality of planar images of the fundus of the eye captured by the adaptive optics SLO are stored in a database connected to the adaptive optics SLO via the network, and the image processing apparatus 10 acquires the images from the database via the network.
  • Next, processing procedures of the image processing apparatus 10 according to the present embodiment will be described with reference to the flowchart of FIG. 2.
  • In step S210, the image acquisition unit 100 acquires a plurality of planar images of the fundus of the eye to be analyzed from the adaptive optics SLO connected to the image processing apparatus 10 or a database where the planar images of the fundus of the eye captured by the apparatus are stored. The acquired planar images of the fundus of the eye are stored in the storage unit 130 via the control unit 120.
  • Further, the image acquisition unit 100 acquires imaging parameter information of the planar images of the fundus of the eye when they are captured and stores the information in the storage unit 130 via the control unit 120. The imaging parameter information is, for example, position information of the fixation lamp when the imaging is performed. The imaging parameter information such as the position information of the fixation lamp may be included in an information file which is added to the planar images of the fundus of the eye, but may also exist as tag information of the images.
  • In step S220, the input information acquisition unit 110 acquires information of the subject's eye from the database or from the input of the operator via an input unit (not illustrated). The information of the subject's eye is information of the patient (e.g., patient ID, name, age, and sex), whether the examination target is right/left eye, and shooting date and time. The acquired information is stored in the storage unit 130 via the control unit 120.
  • In step S230, the frequency conversion unit 141 acquires a spatial frequency image by the discrete Fourier transform using the planar images of the fundus of the eye acquired by the adaptive optics SLO and stored in the storage unit 130. As illustrated in FIG. 3, a greater part of each planar image of the fundus of the eye is populated with regularly-arranged photoreceptor cells observed as small regions with high luminance. Thus, even if the image partially includes a blood vessel or a region of lesion, a Fourier image obtained from such a planar image of the fundus of the eye by spatial frequency conversion has a ring structure as illustrated in FIG. 4.
  • In step S240, the feature extraction unit 142 extracts a feature quantity of the ring structure that shows the periodicity of the arrangement of the photoreceptor cells from the Fourier image obtained in step S230. More specifically, as illustrated in FIG. 5A, if the Fourier image is a square image having a number N of pixels in the vertical and horizontal directions (i.e., a square of N×N pixels), polar coordinates (r, θ) having the center of the Fourier image (coordinates (N/2, N/2)) as the point of origin is considered. Then, a function I(r), which is obtained by calculating the value of each pixel of the Fourier image in the 0 direction, is calculated where r=0, 1, 2 . . . N/2. Since the Fourier image is not a continuous image and each pixel has a value, when the function I(r) is calculated, if r of each pixel is, for example, 4.5 or greater and smaller than 5.5, the value of I(5) is used. Subsequently, the function I(r) is smoothed, for example, by acquiring a mean value of adjacent points. FIG. 5B illustrates a function I(r) acquired from the Fourier image in FIG. 5A.
  • The function I(r) in FIG. 5B includes a lot of information regarding the arrangement of the photoreceptor cells. For example, when the crystalline lens of the subject's eye becomes cloudy due to a disease, the signals of the photoreceptor cells become extremely weak (see FIG. 6A). In such a case, as illustrated in FIG. 6B, the value of the function I(r) itself becomes small compared to the function I(r) in FIG. 5B. Further, if the photoreceptor cells are partially absent in the imaging region and a periodic structure is not produced, a similar function I(r) is obtained. Thus, as a feature quantity indicating the intensity of the periodic structure of the photoreceptor cells, as illustrated in FIG. 7, Imax (i.e., a maximum value of I(r)) and Isum (i.e., a sum of I(r)) can be used.
  • I max = max r I ( r ) I sum = r I ( r )
  • Further, rmax, which is a value of r of Imax, corresponds to the periodicity of the arrangement of the photoreceptor cells.
  • r max = arg max r I ( r )
  • If the distance between adjacent photoreceptor cells is small and the photoreceptor cells are densely arranged, rmax becomes greater. Conversely, if the photoreceptor cells are distant from the macula lutea and the density of the photoreceptor cells is low, rmax becomes smaller.
  • Further, even if the signal strength is strong enough, the periodic structure of the photoreceptor cells may not be clearly determined due to low resolution (see FIG. 8A). In such a case, as illustrated in FIG. 8B, the shape of the function I(r) is broadened toward the smaller function r. As a feature quantity that represents such poor resolution, a ratio of a sum of values around rmax (±n) and a sum of values around r smaller than rmax is calculated. FIG. 8C illustrates the calculation regions.
  • peak_ratio = I sum _ peak I sum _ low = r = r max - n r max + n I ( r ) r = r max - 3 n r max - n I ( r )
  • If the imaging condition is poor and the resolution is low (FIG. 9A), it may be difficult to acquire the peak I(r). FIG. 9B illustrates the signals of I(r) in such a case. In such a case, even if Imax is acquired, accurate peak_ratio cannot be calculated. In such a case, a peak position which is assumed from the acquired visual fixation information (position r_expected_peak) will be used in place of rmax in calculating the peak_ratio. In this manner, the resolution of the image can be evaluated more accurately.
  • In step S250, from the feature quantity acquired in step S240, the index calculation unit 143 acquires information of the photoreceptor cells such as the information of the imaging state and the distribution information of the photoreceptor cells of the image captured by the adaptive optics SLO. The acquired image quality index is stored in the storage unit 130 via the control unit 120. An example of a method for calculating the image quality index from the feature quantity acquired in step S240 will be described below. The index which is acquired, however, is not limited to the calculation method described below.
  • Among a plurality of indices extracted in step S240, the indices associated with the intensity of the periodic structure of the photoreceptor cell such as Imax and Isum and those associated with the resolution such as the peak_ratio are important in terms of image quality. Regarding the adaptive optics SLO, even if the health level of the subject's eye is good and high signal can be obtained, good resolution is not always obtained.
  • In such a case, an image of good quality may be acquired if, for example, the aberration is corrected again using the wavefront sensor. Thus, it is important to differentiate such an image from an image of a subject's eye with a cloudy crystalline lens less likely to produce good image. Thus, by combining the index associated with the intensity of the periodic structure and the index associated with the resolution in a two-dimensional manner, information which can more accurately present the imaging state can be acquired compared to a case where information is acquired from each index.
  • FIGS. 10A and 10B illustrate examples of the index. In FIGS. 10A and 10B, the vertical axis represents Imax, which is an index associated with the intensity of the periodic structure and the horizontal axis represents the peak_ratio, which is an index associated with the resolution. In step S240, a case where rmax is used for calculating the peak_ratio and another case where the position r_expected_peak is used for calculating the peak_ratio have been described. Generally, when the signal strength is reduced, the noise is increased, and it becomes difficult to acquire an accurate rmax. Thus, in FIGS. 10A and 10B, a case where the peak_ratio is calculated using the r_expected_peak is presented.
  • As illustrated in FIG. 10B, an image in a region A, which is a high index region regarding the intensity of the periodic structure and the resolution, is considered as an image of good quality. An image in a region B, which is a low index region regarding the intensity of the periodic structure and the resolution, is considered as an image of poor quality. Regarding an image in a region C, which is a high index region regarding the intensity of the periodic structure but a low index region regarding the resolution, the state of the subject's eye is not bad but the resolution may be poor. If such an image is acquired, by performing the aberration correction again, the image quality may be improved. Such an image may also include both a region where the structure of the photoreceptor cells is clear a region where the structure of the photoreceptor cells is not clear. Such a portion may correspond to a lesion portion. Further, an image in a region D, which is a low index region regarding the intensity of the periodic structure and a high index region regarding the resolution, has good resolution but weak signals in general. If an image is such, a clearer image may be acquired by increasing the signal strength.
  • In step S260, the output unit 150 acquires the information indicating the imaging state, such as the image quality index, stored in the storage unit in step S250 and displays it on the display unit 160.
  • For example, the acquired feature quantity is displayed as it is as the information of the imaging state. Although detailed information can be presented to the user, if a threshold value of the image quality index is set in advance when the image index value is smaller than the threshold value, a message such as “low image quality of photoreceptor cell region” can be displayed on the display unit 160 by the output unit 150 as the information of the imaging state. The information is not necessarily character information and an icon or an image corresponding to the character information may be displayed in place of the character information.
  • In this manner, for example, even if the inspector is not used to the operation, the inspector can easily notice that the image has a quality problem. If the image index value exceeds the threshold value, a message such as “good image quality” may be displayed. Then, the user can understand that the image has no quality problems.
  • Further, if the value of the peak_ratio is low and the value of the Imax is high as described above, the output unit 150 causes the display unit 160 to display a message indicating that the ophthalmologic imaging apparatus needs correction of aberration. According to this message, the user can understand that the correction of aberration is necessary. Further, if the value of the peak_ratio is high and the value of the Imax is low, a message indicating that an increase in imaging light quantity is necessary is displayed on the display unit 160. In this manner, the user can understand that alight quantity adjustment is necessary.
  • Further, as information other than the information of the imaging state, the density information of the photoreceptor cells may be displayed based on the position of the rmax. Further, a character or a graphic that notifies the user that a blood vessel or a lesion region is extracted instead of the photoreceptor cells on the planar image of the fundus of the eye may be displayed.
  • The output unit 150 also causes the display unit 160 to display a frequency image, such as the one illustrated in FIG. 4, together with the image of the fundus of the eye. Further, a graph, such as the one illustrated in FIG. 7, may also be displayed. In this manner, detailed information of the imaging state can be presented to the user. Additionally, information of a distribution state of the photoreceptor cells such as density information may be displayed on the display unit. This information is based on the diameter of the ring structure extracted by the feature extraction unit 142. Further, a message informing the user that the image includes a region not including photoreceptor cells (i.e., a blood vessel or a lesion portion) can be displayed on the display unit.
  • Feature quantity and other information stored in the storage unit 130 in steps S210 to S250 are stored in a database.
  • According to the above-described configuration, an index useful in objectively evaluating the image quality of a planar image of a fundus of the eye captured by the adaptive optics SLO apparatus can be presented. According to the presentation of such an index, for example, when a diagnosis or a determination of an effect of a treatment is performed using the density of the photoreceptor cells, an objective criterion for determination can be presented.
  • According to the first exemplary embodiment, a Fourier image is obtained by a frequency conversion of the entire planar image of the fundus of the eye obtained by the adaptive optics SLO, and an index useful for evaluating the image quality of the entire planar image of the fundus of the eye is calculated by extracting various feature quantities from the obtained Fourier image. However, since the image quality of the planar image of the fundus of the eye is not always consistent and the lesion portions may be unevenly distributed in one planar image, the intensity of the periodicity structure of the photoreceptor cells may be different depending on the portion.
  • According to a second exemplary embodiment, in order to obtain such a local difference, the planar image of the fundus of the eye is segmented into a plurality of local regions. Then, a Fourier image is acquired for each of the planar images of the fundus of the eye acquired by the segmentation, and the planar image is analyzed using the feature quantities extracted from the Fourier images.
  • FIG. 11 illustrates a functional configuration of the image processing apparatus 10 according to the present embodiment. Since the functional configurations of units other than the image processing unit 140 are similar to those illustrated in FIG. 1, their descriptions are not repeated. According to the present embodiment, the image processing unit 140 further includes an image segmentation unit 1140 and an integration unit 1144 as well as the frequency conversion unit 141, the feature extraction unit 142, and the index calculation unit 143. After the planar image of the fundus of the eye is segmented into a plurality of regions, a feature quantity of each region is extracted. The acquired feature quantities are integrated into an integrated index. Then, the image is evaluated according to the integrated index.
  • The image segmentation unit 1140 segments the planar image of the fundus of the eye into a plurality of partial regions. The frequency conversion unit 141 performs the frequency conversion for each partial region and obtains a partial frequency image. The feature extraction unit 142 extracts a feature quantity based on the ring structure that appears on the partial frequency image. The index calculation unit 143 acquires the information of the photoreceptor cells in the plurality of partial regions.
  • By obtaining the frequency image for each partial region, if a region where the photoreceptor cells are not extracted is included in the original image, the image can be segmented into a partial region where, for example, a blood vessel is extracted but photoreceptor cells are not extracted, and a partial region other than such a region. Extremely fine blood vessels and lesion portions buried in the noise and unable to determine from the feature quantity of the frequency image are not considered in the segmentation.
  • Processing procedures of the image processing apparatus 10 according to the present embodiment is described with reference to the flowchart of FIG. 12. Since the processing procedures in steps S210, S220, S230, S240, and S250 are similar to the processing procedures described according to the first exemplary embodiment, their descriptions are not repeated.
  • According to the first exemplary embodiment, the calculation of the image quality index is performed with respect to the entire planar image of the fundus of the eye acquired by the adaptive optics SLO. According to the present embodiment, the planar image of the fundus of the eye is segmented into a plurality of local regions, and the index of each region is calculated. Then, the indices are combined and used for the evaluation of the entire image. Thus, the image processed in steps S230, S240, and S250 is a local planar image of the fundus of the eye obtained by the segmentation of the planar image of the fundus of the eye.
  • Next, each step will be described in detail.
  • In step S1230, the image segmentation unit 1140 segments the planar image of the fundus of the eye acquired by the adaptive optics SLO and stored in the storage unit 130 into a plurality of local regions. Various methods can be used for the segmentation. Although the local difference can be more noticeable if the planar image of the fundus of the eye is segmented into a great number of small regions, the accuracy of the information obtained from each local region will be reduced. Further, since the frequency conversion of a plurality of regions requires more time and cost, it is important to use an image size of data of 2 to the n-th power suitable for high speed Fourier transform. In the present embodiment, for example, with respect to the original planar image of the fundus of the eye having a size of 400×400, a plurality of local planar images of the fundus of the eye with a size of 128×128 are acquired as described in FIG. 13. Each of the local planar images includes an overlapping portion. The segmentation method, however, is not limited to such an example.
  • Sixteen local planar images of the fundus of the eye generated in this manner are stored in the storage unit 130 via the control unit 120. Although the processing in steps S230, S240, and S250 is similar to the processing in the first exemplary embodiment, it is different in that the processing is performed for each of sixteen local planar images of the fundus of the eye generated in step 1230. Accordingly, sixteen indices are obtained. Further, each of the obtained indices and the feature quantities acquired in the course of process is stored in the storage unit 130 in association with the corresponding local planar image of the fundus of the eye.
  • In step S1260, the integration unit 1144 generates an index by integrating the indices acquired from each local planar image of the fundus of the eye. Then, the obtained index is displayed on the monitor via the output unit 150. Further, the feature quantity and other information stored in the storage unit 130 in steps S210 to S1260 are stored in the database.
  • FIGS. 14A and 14B illustrate examples of sixteen indices acquired in step S250. In FIG. 14A, all of the sixteen indices show poor image quality. If the state of the subject's eye is not good and, for example, the crystalline lens is very cloudy, very few signals can be obtained. In such a case, all the indices acquired from the local planar images of the fundus of the eye show bad values. Thus, if indices such as those illustrated in FIG. 14A are obtained, the quality of the planar image of the fundus of the eye is considered as low as a whole.
  • On the other hand, if the indices of both high values and of low values are determined as illustrated in FIG. 14B, it can be considered that some of the local regions failed in the acquisition of the periodic structure of the photoreceptor cells although the imaging itself has been successful. Such a case indicates that the photoreceptor cells are not extracted due to, for example, the existence of a blood vessel, or the photoreceptor cells do not exist due to, for example, a lesion.
  • Regarding the extraction of a blood vessel from a planar image of a fundus of the eye acquired by an adaptive optics SLO, there is a known method such as the one discussed in Tam, J., Martin, J. A., Roorda, A., “Non-invasive visualization and analysis of parafoveal capillaries in humans” Invest. Ophthalmol. Vis. Sci. 51 (3): 1691-1698 (2010). If indices such as those illustrated in FIG. 14B are obtained according to a combination of such methods, whether the possibility of the local region, indicating the low index value, having a blood vessel or a lesion portion is high can be determined.
  • Further, the planar image can be evaluated by a combination of the first and the second exemplary embodiments, for example, a combination of the overall evaluation and the partial evaluation. By presenting a result obtained from the combined evaluation, even if the overall evaluation is not good, the evaluation will be useful since the user can determine the reason for the poor evaluation (e.g., poor imaging condition, lesion portion was included). This is useful as if any part of the partial evaluation indicates high image quality, the possibility of the imaging condition being the cause of the poor image quality is low. Thus, an optimum index can be presented by combining the evaluation methods.
  • According to the present embodiment, the planar image of the fundus of the eye acquired by the adaptive optics SLO apparatus is segmented into a plurality of local regions. Further, the image quality indices obtained from the local planar images of the fundus of the eye are integrated to form an integrated index. According to this index, a case where the overall image quality of the fundus of the eye is low and a case where the image quality itself is not low but a region where the periodic structure of the photoreceptor cells is not extracted exists due to the presence of a blood vessel or a lesion region can be presented.
  • Further, the ophthalmologic imaging apparatus can be configured such, from the aspect of hardware, that the index calculation unit 143 automatically determines a lesion region. Then, based on the result of the determination, the output unit 150 outputs information regarding the imaging region to the apparatus. On receiving the information, the apparatus automatically captures the image again by setting the center of the image to that imaging region.
  • If the output unit 150 outputs an aberration correction instruction or a control value of an adjustment value of a light quantity according to the imaging state to the ophthalmologic imaging apparatus, then the ophthalmologic imaging apparatus controls the aberration measurement unit and the correction optical system based on the input control value. Accordingly, the user does not need to adjust the apparatus.
  • According to the first exemplary embodiment, the information of the imaging state is displayed and the user performs the operation according to the displayed information. According to a third exemplary embodiment, the information of the imaging state is used for forming and selecting an image. Since the configurations of the apparatus are similar to those described in the first and the second exemplary embodiments, their descriptions are not repeated and the points different from the above-described exemplary embodiments will be mainly described.
  • The image acquisition unit 100 acquires a plurality of planar images of the fundus of the eye obtained by imaging of different or substantially same positions. Since the fundus of the eye moves due to the involuntary eye movement such as a saccade, even if a tracking function is provided, the tracking may not be successful. Thus, even if the apparatus is set to capture an image of the same position, actually, an image of a different position may be captured. This is why the expression “substantially the same” is used.
  • The image processing unit 140 functions as a determination unit configured to determine the image quality. The image processing unit 140 acquires the information of the imaging state acquired by the index calculation unit 143 and the information of the threshold value of the imaging state determined in advance. Then, the image processing unit 140 determines whether the imaging state obtained for each image exceeds of falls below the threshold value. As the information of the imaging state, for example, the values of the peak_ratio and Imax described in the first exemplary embodiment or an image quality index value obtained by a combination of such values can be used.
  • Additionally, the image processing unit 140 functions as a selection unit configured to select an image. If an image is determined by the determination unit as having an image quality below the threshold value, the image processing unit 140 determines that it is inappropriate for the presentation to the user or not suitable for the superimposing processing and does not select the image. In other words, the image processing unit 140 selects only the images that exceed the threshold value. Since the images below the threshold value are not selected, images with aberration, images with insufficient light quantity, and images captured when the eye instantaneously moved due to blinking or the like are not selected by the image processing unit 140.
  • The image processing unit 140 further aligns the selected images and performs the superimposing processing. Thus, the image processing unit 140 also functions as a generation unit configured to generate an image with reduced random noise. Accordingly, the image quality of a plurality of images can be automatically determined and a planar image of a fundus of the eye with good image quality can be obtained.
  • The control unit 120 outputs only the image selected from the acquired plurality of images to the output unit 150. The control unit 120 outputs such an image together with the generation unit or in place of the generation unit. The output unit 150 causes the display unit 160 to display only the selected image. The unselected images are not transmitted to the output unit 150 but is stored in the storage unit 130.
  • The unselected images can be deleted as a failed image with the approval of the user. In addition to the display of the image for the diagnosis, the output unit 150 can separately display the images which have not been selected on the display unit 160 and further display a window that accepts the approval of the user regarding the deletion. The control unit 120 can delete the approved images from the storage unit 130 in one operation. Further, in place of the deletion, the control unit 120 can cause only the approved images not to be transferred to the external database. In this manner, the user can automatically determine unnecessary images and can easily perform the deletion processing or the non-transfer processing.
  • The function of the above-described image processing apparatus 10 can be realized by the hardware illustrated in FIG. 15 and using the software and hardware together.
  • The image processing apparatus illustrated in FIG. 15 includes a central processing unit (CPU) 1501, a random access memory (RAM) 1502, a read only memory (ROM) 1503, an external storage unit 1504, and a communication interface 1505. These units are connected to each other via a bus 1509. Further, a monitor 1506, a keyboard 1507, and a mouse 1508 are connected to the image processing apparatus 10. A program including an instruction that can cause the processing illustrated in the flowchart of FIG. 2 or 12 is stored in the ROM 1503 or the external storage unit 1504.
  • The function of the image processing apparatus 10 of the above-described exemplary embodiment is realized by the CPU 1501 reading out a stored program, loading it into the RAM 1502, and executing the instruction included in the program.
  • An example of an imaging system which acquires a planar image of the fundus of the eye described above will be described with reference to FIG. 16. The image processing apparatus 10 in the imaging system is connected to an ophthalmologic imaging apparatus 30 and a database 50 via a local area network (LAN) 40. The image processing apparatus 10 can be wirelessly connected to the LAN 40. Further, the display unit 160 is connected to the image processing apparatus 10.
  • The configuration of the imaging system is not limited to the above-described example. For example, an ophthalmologic imaging apparatus which includes an ophthalmologic imaging unit, an image processing unit, and the display unit 160 may be used. In this case, the ophthalmologic imaging unit includes the function of the ophthalmologic imaging apparatus 30 and the image processing unit includes the function of the image processing apparatus 10. If the imaging system has such a configuration, the quality of the image of the photoreceptor cells can be determined by the ophthalmologic imaging apparatus 30 alone, and the system can be downsized.
  • Configuration of the ophthalmologic imaging apparatus 30 will be described with reference to FIG. 17.
  • A composite apparatus composed of an SLO apparatus and an OCT apparatus according to the embodiments of the present invention is described as an ophthalmologic imaging apparatus. In particular, the ophthalmologic imaging apparatus 30 which includes an adaptive optical system and can capture both planar images (SLO images) of high lateral resolution as well as tomographic images (OCT images) for a retina and acquires a planar image of a fundus of the eye will be described. The ophthalmologic imaging apparatus 30 includes an SLO apparatus and an OCT apparatus. The SLO apparatus acquires a planar image of a fundus of the eye by correcting the optical aberration of the subject's eye and using a spatial light modulator. The OCT apparatus employs Fourier domain imaging when it acquires a tomographic image. According to these apparatuses, the ophthalmologic imaging apparatus 30 can obtain a good planar image of the fundus of the eye and a tomographic image of the subject's eye regardless of the visibility and the optical aberration.
  • First, an overall configuration of the ophthalmologic imaging apparatus 30 according to the present embodiment will be described in detail with reference to FIG. 17. Light emitted from a light source 201 is divided into a reference beam 205 and a measuring beam 206 by an optical coupler 231. The measuring beam 206 is guided to a subject's eye 207 as an object to be observed via a single mode fiber 230-4, a spatial light modulator 259, an XY scanner 219, an X scanner 221, and spherical mirrors 260-1 to 260-9.
  • The measuring beam 206 is reflected or scattered by the subject's eye 207 being an object to be observed and returned as a return beam 208. The return beam 208 enters a detector 238 or a line sensor 239. The detector 238 converts the light intensity of the return beams 208 into a voltage signal. Then, based on the voltage signal, a planar image of the fundus of the eye of the subject's eye 207 is generated. Further, the return beam 208 is combined with the reference beam 205, and the combined light is caused to enter the line sensor 239. Accordingly, a tomographic image of the subject's eye 207 is formed. Furthermore, by using a plurality of acquired tomographic images, a three-dimensional course of blood vessels can be extracted.
  • Although a spatial light modulator is used as a device for correcting the wavefront aberration, in the present embodiment, any device can be used so long as the wavefront aberration can be corrected. Thus, for example, a variable shape mirror can be used.
  • Next, the periphery of the light source 201 will be described. The light source 201 uses a super luminescent diode (SLD) being a typical low-coherent light source. The central wavelength and the bandwidth are 830 nm and 50 nm, respectively. In order to acquire a planar image of the fundus of the eye with small speckle noise, a low-coherent light source is selected. Although the SLD is selected as a type of the light source, a different light source can be used so long as low coherent light can be emitted. For example, an amplified spontaneous emission (ASE) light source can be used.
  • Further, near-infrared light is suitable as a wavelength from the viewpoint of measurement of eyes. Further, since the wavelength affects the lateral resolution of the acquired planar image of the fundus of the eye, a shorter wavelength is desirable. In the descriptions below, a wavelength of 830 nm is used. Other wavelengths may be selected depending on the measuring portion of the object to be observed.
  • Light emitted from the light source 201 is divided into the reference beam 205 and the measuring beam 206 at a ratio of 96:4 via a single mode fiber 230-1 and the optical coupler 231. The composite apparatus 30 further includes polarization controllers 253-1 and 253-2.
  • Next, an optical path of the reference beam 205 will be described.
  • The reference beam 205 divided by the optical coupler 231 is adjusted so that it is guided to a lens 235-1 via a single mode fiber 230-2 and adjusted so as to become parallel light with a beam diameter of 4 mm.
  • Next, the reference beam 205 is guided to a mirror 214 being a reference mirror by mirrors 257-1 to 257-4. Since the optical path length of the reference beam 205 is adjusted to be substantially the same as that of the measuring beam 206, the reference beam 205 and the measuring beam 206 interfere with each other. Next, the reference beam 205 is reflected by the mirror 214 and guided again to the optical coupler 231. A dispersion compensation glass 215, through which the reference beam 205 passes, compensates the dispersion that occurs when the measuring beam 206 travels to the subject's eye 207 and returns from the subject's eye 207 with respect to the reference beam 205. In the following descriptions, for example, an average diameter of an oculus of a Japanese being L1=23 mm is used.
  • A motorized stage 217-1 is movable in the directions indicated by an arrow to allow the optical path length of the reference beam 205 to be adjusted and controlled. Further, the motorized stage 217-1 is driven by a motorized stage driver 283 in a driver unit 281 under the control of a personal computer 225.
  • Next, the optical path of the measuring beam 206 will be described. The measuring beam 206 split by the optical coupler 231 is guided to a lens 235-4 via the single mode fiber 230-4 and adjusted so as to become parallel light with a beam diameter of 4 mm. Further, the polarization controller 253-1 or 253-2 can adjust the polarization state of the measuring beam 206. In the present embodiment, the polarization controller 253-1 or 253-2 adjusts the polarization state of the measuring beam 206 to be linearly polarized in a direction parallel to the drawing surface of FIG. 17.
  • The measuring beam 206 passes through a beam splitter 258 and a movable beam splitter 261 (also referred to as a splitting unit) and reaches the spatial light modulator 259 via the spherical mirrors 260-1 and 260-2. Then, the measuring beam 206 is modulated at the spatial light modulator 259. In the present embodiment, the spatial light modulator 259 utilizes the orientation characteristics of the liquid crystal. More specifically, the spatial optical modulator 259 is arranged in a direction where the spatial optical modulator 259 can modulate the phase of linear polarization parallel to the drawing surface of FIG. 17 (i.e., the P polarization) to coincide with the polarization orientation of the measuring beam 206.
  • Further, the measuring beam 206 passes through a polarizing plate 273 and reaches a mirror of the X scanner 221 via the spherical mirrors 260-3 and 260-4. In the present embodiment, the polarizing plate 273 has a role to guide only the linear polarization parallel to the drawing surface of FIG. 17, of the return beams 208, to the spatial light modulator 259. Further, in the present embodiment, the X scanner 221 scans the measuring beam 206 in a direction parallel to the drawing surface of FIG. 17. For example, the X scanner 221 is a resonant scanner having a drive frequency of approximately 7.9 kHz.
  • Further, the measuring beam 206 is incident on a mirror of the XY scanner 219 via the spherical mirrors 260-5 and 260-6. Although the number of the mirrors in the XY scanner 219 is illustrated as one, actually, the XY scanner 219 includes two mirrors (X-scanning mirror and Y-scanning mirror) arranged close to each other. Further, the measuring beam 206 is adjusted in such a manner that its center coincides with the center of rotation of the mirror of the XY scanner 219. The drive frequency of the XY scanner 219 is changeable in the range up to 500 Hz.
  • The spherical mirrors 260-7 to 260-9 are optical systems that cause the measuring beam 203 to scan a retina 227. Having a point near a cornea 226 as a support point, the measuring beam 206 scans the retina 227. Although the measuring beam 206 has a beam diameter of 4 mm, the beam diameter may be increased for acquisition of a high resolution tomographic image.
  • A motorized stage 217-2 is configured to move in a direction indicated by an arrow in FIG. 17 and to adjust and control the position of a spherical mirror 260-8 attached thereto. Similarly to the motorized stage 217-1, the motorized stage 217-2 is controlled by the motorized stage driver 283.
  • By adjusting the position of the spherical mirror 260-8, it is possible to focus the measuring beam 206 on a predetermined layer of the retina 227 of the subject's eye 207 and observe the subject's eye 207. In the initial state, the position of the spherical mirror 260-8 is adjusted so that the measuring beam 206 can enter the cornea 226 as parallel light. The ophthalmologic imaging apparatus 30 according to the present embodiment can cope with the subject's eye 207 having a refractive error.
  • When the measuring beam 206 enters the subject's eye 207, the measuring beam 206 becomes the return beam 208 due to reflection or scattering from the retina 227. Then, the return beam 208 reaches the line sensor 239 by being guided again by the optical coupler 231. A part of the return beam 208 is reflected by the movable beam splitter 261 and guided to the detector 238 via a lens 235-5.
  • A light blocking plate 272 having a pinhole has a function of blocking unnecessary light not focused on the retina 227 among the return beams 208. Further, the light blocking plate 272 is arranged in a position conjugate with the focusing position of the lens 235-5. The diameter of the pinhole of the light blocking plate 272 is, for example, 50 μm. As the detector 238, for example, an avalanche photodiode (APD), which is a high-speed, highly-sensitive light sensor, is used. A part of the return beam 208 split by the beam splitter 258 enters a wavefront sensor 255. The wavefront sensor 255 is a Shack-Hartmann wavefront sensor.
  • The spherical mirrors 260-1 to 260-9 are arranged so that they are optically conjugate with the XY scanner 219, the X scanner 221, the cornea 226, the wavefront sensor 255, and the spatial light modulator 259. Thus, the wavefront sensor 255 can measure the aberration of the subject's eye 207. Further, the spatial light modulator 259 can correct the aberration of the subject's eye 207. Furthermore, by controlling the spatial light modulator 259 in real time based on the obtained aberration, it is possible to correct the aberration that occurs in the subject's eye 207 and to acquire tomographic images with higher lateral resolution.
  • Next, the configuration of the measurement system will be described. The ophthalmologic imaging apparatus 30 can acquire tomographic images (OCT images) as well as planar images of the fundus of the eye (SLO images).
  • First, the measurement system for the tomographic images will be described. The return beams 208 are combined by the optical coupler 231. The combined light (combined light 242) is guided to a transmissive grating 241 via a single mode fiber 230-3 and a lens 235-2, and is dispersed for each wavelength. Then, the light enters the line sensor 239 via a lens 235-3.
  • The line sensor 239 converts the light intensity for each position (wavelength) into a voltage signal. The voltage signal is converted into a digital value by a frame grabber 240 so that tomographic images of the subject's eye 207 are formed in the personal computer 225. The line sensor 239 includes 1024 pixels and can obtain the intensity of the combined light 242 for each wavelength (a segmentation of 1024).
  • Next, the measurement system for the planar image of the fundus of the eye will be described. A part of the return beam 208 is reflected by the movable beam splitter 261. After unnecessary light is blocked by the light blocking plate 272, the reflected light reaches the detector 238, and the light intensity is converted into an electric signal. Data of the obtained electric signal is processed by the personal computer 225 in synchronization with the scanning signal of the X scanner 221 and the XY scanner 219. Accordingly, a planar image of the fundus of the eye is formed.
  • A part of the return beam 208 split by the beam splitter 258 enters the wavefront sensor 255, and the aberration of the return beam 208 is measured. An image signal obtained by the wavefront sensor 255 is supplied to the personal computer 225, and the aberration is calculated. The obtained aberration is expressed by Zernike polynomial, which represents the aberration of the subject's eye 207. The Zernike polynomial includes a tilt term, a defocus term, an astigmatism term, a coma term, and a trefoil term.
  • Next, an acquisition method for tomographic images (OCT images) using the ophthalmologic imaging apparatus 30 will be described with reference to FIGS. 18A to 18C. The ophthalmologic imaging apparatus 30 acquires a tomographic image of the retina 227 by controlling the XY scanner 219 and acquiring an interference fringe by the line sensor 239 using the X scanner 221 as a fixed mirror. The movable beam splitter 261 is controlled so that the return beams 208 are not guided to the detector 238. Further, the X scanner 221 and the XY scanner 219 are controlled by an optical scanner driver 282 in the driver unit 281 from the personal computer 225. In the present embodiment, the method for acquiring a tomographic image (a plane parallel to the optical axis) of the retina 227 will be described.
  • FIG. 18A is a schematic view of the subject's eye 207 and illustrates a state where the subject's eye 207 is observed by the ophthalmologic imaging apparatus 30. As illustrated in FIG. 10A, when the measuring beam 206 enters the retina 227 via the cornea 226, the measuring beam 206 becomes the return beams 208 due to reflection or scattering at various positions. Then, the return beam 208 reaches the line sensor 239 with delays at the respective positions.
  • Since the light source 201 has a large bandwidth and a short coherence length, the line sensor 239 can detect an interference fringe when the optical path length of the reference beam path is approximately equal to the optical path length of the measuring beam path. As described above, the interference fringe acquired by the line sensor 239 is an interference fringe in a spectrum region on the wavelength axis.
  • Subsequently, the interference fringe, which is information of the wavelength axis, is converted into an interference fringe on an optical frequency axis considering the characteristics of the line sensor 239 and the transmissive grating 241. Further, by performing inverse Fourier transform of the obtained interference fringe on the optical frequency axis, it is possible to obtain information in the depth direction. As illustrated in FIG. 18B, if the interference fringe is detected while the XY scanner 219 is driven, the interference fringe can be obtained at each position on the X-axis. More specifically, the information in the depth direction can be obtained at each position on the X-axis. As a result, a two-dimensional distribution of the intensity of the return beams 208 on the XZ plane is obtained, which is a tomographic image 232 (see FIG. 18C).
  • As described above, the tomographic image 232 is an image obtained by arranging the intensity of each of the return beams 208 in an array. For example, the tomographic image 232 is a gray scale intensity image. The length of the tomographic image 232 in the X direction is 700 μm, which is similar to the SLO image described below.
  • In FIG. 18C, only the boundaries in the obtained tomographic image are expressed in lines. The illustration includes a retinal pigment epithelial layer 246, a stratum opticum 247, and a blood vessel 278. Further, if a plurality of tomographic images is acquired at each position on the Y-axis, it is possible to visualize a three-dimensional course of the blood vessel.
  • Next, an acquisition method of a planar image of the fundus of the eye (SLO image) using the ophthalmologic imaging apparatus 30 will be described. The ophthalmologic imaging apparatus 30 can acquire a planar image of the fundus of the eye of the retina 227 by controlling and operating the XY scanner 219 in only the Y-axis direction and the X scanner 221 while fixing the X-axis direction of the XY scanner 219, and acquiring the intensity of the return beam 208 using the detector 238.
  • The X scanner 221 and the XY scanner 219 are controlled by the optical scanner driver 282 in the driver unit 281 from the personal computer 225 (see FIG. 17). Further, the ophthalmologic imaging apparatus 30 can acquire a planar image of the fundus of the eye while correcting the aberration that occurs at the subject's eye 207 by controlling the spatial light modulator 259. The spatial light modulator 259 can be controlled by using the aberration of the subject's eye 207 measured by the wavefront sensor 255. Furthermore, the ophthalmologic imaging apparatus 30 can acquire a planar image of the fundus of the eye while controlling the spatial light modulator 259 in real time.
  • Next, the acquisition method of the planar image of the fundus of the eye (SLO image) will be described with reference to FIGS. 19A to 19D.
  • The ophthalmologic imaging apparatus 30 acquires a planar image of a fundus of the eye of the retina 227 by controlling the XY scanner 219 and acquiring the intensity of the return beam 208 by the detector 238. An acquisition method of a planar image of the fundus of the eye of the retina 227 (planar image perpendicular to the optical axis) will now be described.
  • FIG. 19A is a schematic view of the subject's eye 207 and illustrates a state where the subject's eye 207 is observed by the ophthalmologic imaging apparatus 30.
  • As illustrated in FIG. 19A, when the measuring beam 206 enters the retina 227 via the cornea 226, the measuring beam 206 becomes the return beam 208 due to reflection or scattering at various positions. Then, the return beam 208 reaches the detector 238. Further, as illustrated in FIG. 19B, if the intensity of the return beam 208 is detected while the XY scanner 209 is moved in the X direction, information at each position on the X-axis can be obtained.
  • As illustrated in FIG. 19C, the XY scanner 209 is simultaneously moved in both the X-axis and Y-axis directions with respect to an image capturing range 292 where the retina 227 is present. Then, the measuring beam 206 is raster-scanned along a trajectory 293 on the image capturing range 292. In this state, if the intensity of the return beams 208 is detected, a two-dimensional distribution of the intensity of the return beams 208 can be obtained. Accordingly, a planar image of the fundus of the eye 277 (see FIG. 19D) is acquired.
  • In FIG. 19D, the measuring beam 206 is scanned from a point S at the upper right to a point E at the bottom left. While the scanning is performed, the intensity of the return beams 208 is used in forming the planar image of the fundus of the eye 277. The trajectory 293 from the point E to the point S is the movement of the measuring beam 206 as a preparation for the imaging of the planar image of the fundus of the eye 277 which is to be formed next. The time required for the scanning takes 84% for the point S to the point E and 16% for the point E to the point S with respect to the trajectory 293 in FIG. 19C. This ratio is based on the duty ratio of the drive waveform of the Y scanner described above. Further, in FIG. 19C, for the sake of simplicity, the number of times of the scanning in the X direction with respect to the trajectory 293 is smaller than the actual number of times.
  • The planar image of the fundus of the eye 277 has a size of 700×350 μm. The time required for the acquisition is approximately 15.6 ms. This time is based on the drive frequency of the Y scanner.
  • In the planar image 277 of the fundus of the eye, a photoreceptor cell group 279, where the intensity of the return beam 208 is relatively high, is light colored, whereas the blood vessel 278, where the intensity is relatively low, is dark colored. Further, blood cells (not illustrated) are light colored in the blood vessel 278. If the planar image of the fundus of the eye 277 is continuously acquired, the movement of the blood cells through the blood vessel 278 can be visualized. Further, spatio-temporal images may be generated by extracting the blood vessel 278, from which the blood cells are visualized, from the planar images 277 of the fundus of the eye which are continuously acquired, and superimposing the extracted planar images 277 of the fundus of the eye in the order they have been captured. Movement of the blood cells and the blood speed can be easily obtained.
  • An imaging method for obtaining an image useful for obtaining the blood speed and an imaging method for obtaining an image of the photoreceptor cells are designated by the user via an operation unit of the image processing apparatus 10. When the user selects the imaging method via the operation unit, a designation unit (not illustrated) of the image processing apparatus designates the selected imaging method and transmits associated information to the ophthalmologic imaging apparatus 30. Then, according to the designated imaging method, the ophthalmologic imaging apparatus 30 outputs an image of photoreceptor cells which is easy to observe or an image having a bloodstream as an object to be observed.
  • The embodiment can also be realized by a computer of a system or apparatus (or devices such as a CPU or MPU) that reads out and executes a program recorded on a memory device to perform the functions of the above-described embodiment(s), and by a method, the steps of which are performed by a computer of a system or apparatus by, for example, reading out and executing a program recorded on a memory device to perform the functions of the above-described embodiment(s). For this purpose, the program is provided to the computer for example via a network or from a recording medium of various types serving as the memory device (e.g., computer-readable storage medium).
  • While the embodiment of the present invention has been described with reference to exemplary embodiments, it is to be understood that the invention is not limited to the disclosed exemplary embodiments. The scope of the following claims is to be accorded the broadest interpretation so as to encompass all modifications, equivalent structures, and functions.

Claims (11)

What is claimed is:
1. An image processing apparatus comprising:
a conversion unit configured to acquire a frequency image by performing conversion of a fundus image obtained by capturing an image of a fundus of an eye;
an extraction unit configured to extract a feature amount from the frequency image; and
an information acquisition unit configured to acquire information indicating a periodicity of an arrangement of photoreceptor cells in the fundus image based of the extracted feature amount.
2. The apparatus according to claim 1, further comprising a segmentation unit configured to segment the fundus image into a plurality of partial images,
wherein the conversion unit obtains partial frequency images by performing conversion on each of the plurality of partial images,
wherein the extraction unit extracts the feature amount from each of the partial frequency images, and
wherein information acquisition unit acquires the information indicating the information of each of the plurality of partial images.
3. The apparatus according to claim 2, wherein the information acquisition unit acquires, based on the feature amount corresponding to each partial image, information about a partial image including a region from which photoreceptor cells are not extracted.
4. The apparatus according to claim 1, wherein the extraction unit extracts a peak value of luminance of a ring structure of the frequency image as the feature amount.
5. The apparatus according to claim 1, wherein the extraction unit extracts a sharpness of a ring structure of the frequency image as the feature amount.
6. The apparatus according to claim 1, further comprising an output unit configured to output the information indicating the periodicity.
7. The apparatus according to claim 6, wherein the output unit causes a display unit to display the fundus image and the information indicating the periodicity.
8. The apparatus according to claim 6, wherein the output unit outputs a control value corresponding to the feature amount to an imaging apparatus configured to capture the fundus image.
9. The apparatus according to claim 1, further comprising an adapted optics scanning laser opthalmoscope configured to capture the fundus image.
10. An image processing method comprising:
obtaining a frequency image by performing conversion of a fundus image obtained by capturing an image of a fundus of an eye;
extracting a feature amount from the frequency image; and
acquiring information indicating a periodicity of an arrangement of photoreceptor cells in the fundus image based of the extracted feature amount.
11. Anon-transitory computer-readable storage medium storing computer-executable instructions for causing a computer to execute the image processing method according to claim 10.
US14/597,127 2011-09-20 2015-01-14 Image processing apparatus, ophthalmologic imaging apparatus, image processing method, and storage medium Abandoned US20150124218A1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US14/597,127 US20150124218A1 (en) 2011-09-20 2015-01-14 Image processing apparatus, ophthalmologic imaging apparatus, image processing method, and storage medium
US16/287,926 US11058293B2 (en) 2011-09-20 2019-02-27 Image processing apparatus, ophthalmologic imaging apparatus, image processing method, and storage medium

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2011204653A JP5843542B2 (en) 2011-09-20 2011-09-20 Image processing apparatus, ophthalmologic photographing apparatus, image processing method, and program
JP2011-204653 2011-09-20
US13/605,873 US8960906B2 (en) 2011-09-20 2012-09-06 Image processing apparatus, ophthalmologic imaging apparatus, image processing method, and storage medium
US14/597,127 US20150124218A1 (en) 2011-09-20 2015-01-14 Image processing apparatus, ophthalmologic imaging apparatus, image processing method, and storage medium

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
US13/605,873 Continuation US8960906B2 (en) 2011-09-20 2012-09-06 Image processing apparatus, ophthalmologic imaging apparatus, image processing method, and storage medium

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US16/287,926 Continuation US11058293B2 (en) 2011-09-20 2019-02-27 Image processing apparatus, ophthalmologic imaging apparatus, image processing method, and storage medium

Publications (1)

Publication Number Publication Date
US20150124218A1 true US20150124218A1 (en) 2015-05-07

Family

ID=47880360

Family Applications (3)

Application Number Title Priority Date Filing Date
US13/605,873 Expired - Fee Related US8960906B2 (en) 2011-09-20 2012-09-06 Image processing apparatus, ophthalmologic imaging apparatus, image processing method, and storage medium
US14/597,127 Abandoned US20150124218A1 (en) 2011-09-20 2015-01-14 Image processing apparatus, ophthalmologic imaging apparatus, image processing method, and storage medium
US16/287,926 Active 2033-01-26 US11058293B2 (en) 2011-09-20 2019-02-27 Image processing apparatus, ophthalmologic imaging apparatus, image processing method, and storage medium

Family Applications Before (1)

Application Number Title Priority Date Filing Date
US13/605,873 Expired - Fee Related US8960906B2 (en) 2011-09-20 2012-09-06 Image processing apparatus, ophthalmologic imaging apparatus, image processing method, and storage medium

Family Applications After (1)

Application Number Title Priority Date Filing Date
US16/287,926 Active 2033-01-26 US11058293B2 (en) 2011-09-20 2019-02-27 Image processing apparatus, ophthalmologic imaging apparatus, image processing method, and storage medium

Country Status (3)

Country Link
US (3) US8960906B2 (en)
JP (1) JP5843542B2 (en)
CN (1) CN103006174B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107292877A (en) * 2017-07-05 2017-10-24 北京至真互联网技术有限公司 A kind of right and left eyes recognition methods based on eye fundus image feature
US11315241B2 (en) * 2018-06-08 2022-04-26 Shanghai Sixth People's Hospital Method, computer device and storage medium of fundus oculi image analysis

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5610884B2 (en) 2010-07-09 2014-10-22 キヤノン株式会社 Optical tomographic imaging apparatus and optical tomographic imaging method
JP5835938B2 (en) * 2011-05-10 2015-12-24 キヤノン株式会社 Aberration correction method, fundus imaging method using the method, and fundus imaging apparatus
JP6146951B2 (en) 2012-01-20 2017-06-14 キヤノン株式会社 Image processing apparatus, image processing method, photographing apparatus, and photographing method
JP5936368B2 (en) 2012-01-20 2016-06-22 キヤノン株式会社 Optical coherence tomography apparatus and method for operating the same
JP6061554B2 (en) 2012-01-20 2017-01-18 キヤノン株式会社 Image processing apparatus and image processing method
JP6039185B2 (en) 2012-01-20 2016-12-07 キヤノン株式会社 Imaging device
JP5988772B2 (en) * 2012-01-20 2016-09-07 キヤノン株式会社 Image processing apparatus and image processing method
JP2013148509A (en) 2012-01-20 2013-08-01 Canon Inc Image processing device and image processing method
WO2013143113A1 (en) * 2012-03-30 2013-10-03 Intel Corporation Techniques for media quality control
US9483111B2 (en) 2013-03-14 2016-11-01 Intel Corporation Techniques to improve viewing comfort for three-dimensional content
JP6241091B2 (en) * 2013-06-27 2017-12-06 株式会社ニデック Image processing apparatus and image processing program
JP6410498B2 (en) * 2014-07-08 2018-10-24 キヤノン株式会社 Ophthalmic apparatus and control method thereof
JP6506518B2 (en) * 2014-08-29 2019-04-24 キヤノン株式会社 Ophthalmic apparatus and control method thereof
JP2016150090A (en) * 2015-02-17 2016-08-22 キヤノン株式会社 Imaging apparatus and control method of the same
KR101825923B1 (en) * 2015-09-03 2018-03-22 주식회사 이오테크닉스 Laser processing apparatus and laser processing method
KR101746763B1 (en) 2016-02-01 2017-06-14 한국과학기술원 Diagnostic technique and imaging method of retinal and choroidal optical coherence tomography angiography
JP2017185057A (en) * 2016-04-06 2017-10-12 キヤノン株式会社 Ophthalmic imaging apparatus and control method therefor, and program
CN107423726A (en) * 2017-08-12 2017-12-01 合肥全飞信息科技有限公司 A kind of illness in eye identifies cell phone application design method
CN108272434B (en) * 2017-12-07 2020-06-19 柯鑫 Method and device for processing fundus images
CN109238131B (en) * 2018-08-09 2020-12-22 江苏度微光学科技有限公司 Transverse ultrahigh-resolution optical coherence tomography method and system
WO2021255779A1 (en) * 2020-06-15 2021-12-23 株式会社ニコン Information display method, information display device, and program
CN111751013B (en) * 2020-07-07 2021-09-03 中国科学院长春光学精密机械与物理研究所 Aberration measuring method for optical imaging and optical imaging method
CN113974964B (en) * 2021-12-28 2022-04-22 广东麦特维逊医学研究发展有限公司 Laser capsulorhexis device and working method thereof

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3030073A (en) * 1958-08-11 1962-04-17 Auto Specialties Mfg Co Positioning apparatus
WO2003030073A1 (en) * 2001-10-03 2003-04-10 Retinalyze Danmark A/S Quality measure
US20100195048A1 (en) * 2009-01-15 2010-08-05 Physical Sciences, Inc. Adaptive Optics Line Scanning Ophthalmoscope

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1427328B1 (en) 2001-08-30 2014-07-02 University Of Rochester Adaptive optics in a scanning lase ophtalmoscope
EP2130486B1 (en) * 2008-06-02 2016-03-23 Nidek Co., Ltd. Ophthalmic Photographing Apparatus
JP5209377B2 (en) * 2008-06-02 2013-06-12 株式会社ニデック Fundus photographing device
CN101833697A (en) 2009-03-13 2010-09-15 新奥特硅谷视频技术有限责任公司 Monitoring system based on iris identification

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3030073A (en) * 1958-08-11 1962-04-17 Auto Specialties Mfg Co Positioning apparatus
WO2003030073A1 (en) * 2001-10-03 2003-04-10 Retinalyze Danmark A/S Quality measure
US20100195048A1 (en) * 2009-01-15 2010-08-05 Physical Sciences, Inc. Adaptive Optics Line Scanning Ophthalmoscope

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107292877A (en) * 2017-07-05 2017-10-24 北京至真互联网技术有限公司 A kind of right and left eyes recognition methods based on eye fundus image feature
US11315241B2 (en) * 2018-06-08 2022-04-26 Shanghai Sixth People's Hospital Method, computer device and storage medium of fundus oculi image analysis

Also Published As

Publication number Publication date
US11058293B2 (en) 2021-07-13
CN103006174A (en) 2013-04-03
US20190191987A1 (en) 2019-06-27
US8960906B2 (en) 2015-02-24
JP2013063215A (en) 2013-04-11
JP5843542B2 (en) 2016-01-13
CN103006174B (en) 2016-07-06
US20130070202A1 (en) 2013-03-21

Similar Documents

Publication Publication Date Title
US11058293B2 (en) Image processing apparatus, ophthalmologic imaging apparatus, image processing method, and storage medium
US8602556B2 (en) Image processing apparatus and image processing system for displaying information about ocular blood flow
US8899749B2 (en) Image processing apparatus, image processing method, image processing system, SLO apparatus, and program
US9044167B2 (en) Image processing device, imaging system, image processing method, and program for causing computer to perform image processing
US9456748B2 (en) Ophthalmological apparatus, alignment method, and non-transitory recording medium
JP5721412B2 (en) Ophthalmic apparatus, blood flow velocity calculation method and program
JP5721411B2 (en) Ophthalmic apparatus, blood flow velocity calculation method and program
US9082010B2 (en) Apparatus and a method for processing an image of photoreceptor cells of a fundus of an eye
US10993615B2 (en) Image processing apparatus and image processing system for displaying information about ocular blood flow
JP6322042B2 (en) Ophthalmic photographing apparatus, control method thereof, and program
US20200375453A1 (en) Tomographic image acquisition apparatus and tomographic image acquisition method
JP2022040372A (en) Ophthalmologic apparatus
JP2022189963A (en) Ophthalmologic apparatus
US20170340203A1 (en) Imaging apparatus and control method for the same
JP6882242B2 (en) Ophthalmic equipment and its control method
JP2020031873A (en) Ophthalmologic apparatus, control method thereof, program, and recording medium
JP7339011B2 (en) Ophthalmic device, ophthalmic information processing device, program, and recording medium
JP2018158153A (en) Ophthalmology imaging device
JP2024099210A (en) Optical coherence tomography device, control method of the same, and program
JP2017176716A (en) Image processing apparatus and image processing method
JP2019176970A (en) Ophthalmologic apparatus and ophthalmologic information processing program

Legal Events

Date Code Title Description
STPP Information on status: patent application and granting procedure in general

Free format text: FINAL REJECTION MAILED

STCB Information on status: application discontinuation

Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION