Method of bioimage data processing for revealing more meaningful anatomic features of diseased tissues

• US 7,668,342 B2
• Filed: 09/09/2005
• Issued: 02/23/2010
• Est. Priority Date: 09/09/2005
• Status: Active Grant

• Current Assignee: Carl Zeiss Meditec Incorporated (Carl-Zeiss-Stiftung)
• Original Assignee: Carl Zeiss Meditec Incorporated (Carl-Zeiss-Stiftung)
• Inventors: Stetson, Paul F., Zhou, Yan, Meyer, Scott A., Everett, Matthew J.
• Primary Examiner(s): Tabatabai; Abolfazl
• Application Number:
  US11/223,549
• Publication Number:
  US 20070103693A1
• Time in Patent Office: 1,628 Days
• Field of Search: 382/106, 382/128, 382/129, 382/130, 382/131, 382/132, 382/133, 382/134, 600/318, 600/356, 600/383, 600400-402, 600/407, 600/587, 600/592, 600/595, 600/442, 600/445, 600/459, 351/205, 351/211, 351/212, 351/204, 356/479, 356/482, 128/916, 128/920, 128/922
• US Class Current: 382/106
• CPC Class Codes: A61B 5/0066 Optical coherence imaging
• CPC Class Codes: G06T 19/00 Manipulating 3D models or images for computer graphics
• CPC Class Codes: G06T 2200/24 involving graphical user interfaces [GUIs]
• CPC Class Codes: G06T 2207/10101 Optical tomography; Optical coherence tomography [OCT]
• CPC Class Codes: G06T 2207/30041 Eye; Retina; Ophthalmic
• CPC Class Codes: G06T 2210/41 Medical
• CPC Class Codes: G06T 7/0014 using an image reference approach

Abstract

The present invention discloses a method for generating elevation maps or images of a tissue layer/boundary with respect to a fitted reference surface, comprising the steps of finding and segmenting a desired tissue layer/boundary; fitting a smooth reference surface to the segmented tissue layer/boundary; calculating elevations of the same or other tissue layer/boundary relative to the fitted reference surface; and generating maps of elevation relative to the fitted surface. The elevation can be displayed in various ways including three-dimensional surface renderings, topographical contour maps, contour maps, en-face color maps, and en-face grayscale maps. The elevation can also be combined and simultaneously displayed with another tissue layer/boundary dependent set of image data to provide additional information for diagnostics.

Claims

1. A method for generating image maps from a 3D image data set, said data set obtained from imaging a sample using an imaging device, said sample having at least one boundary surface associated therewith, said method comprising the steps of:

identifying a subset of the image data corresponding to the boundary surface;

fitting a substantially smooth reference surface to the subset of image data;

calculating the distance between points on the reference surface to points of interest within the sample identified from the 3D image data set;

generating an image map of the calculated distances; and

displaying the image map.

2. A method as recited in claim 1, wherein the points of interest correspond to points on said boundary surface.

3. A method as recited in claim 1, wherein the points of interest correspond to points on a different boundary surface.

4. A method as recited in claim 1, wherein the sample is skin.

5. A method as recited in claim 1, wherein the sample is a blood vessel.

6. A method as recited in claim 1, wherein the sample is an eye.

7. A method as recited in claim 6, wherein the boundary surface is the cornea.

8. A method as recited in claim 6, wherein the boundary surface is the retina.

9. A method as recited in claim 1, wherein the 3D image data set is obtained using an OCT system.

10. A method as recited in claim 1, wherein the 3D image data set is obtained from a confocal microscope.

11. A method as recited in claim 1, wherein the 3D image data set is obtained from an ultrasound scanner.

12. A method for generating information from a 3D image data set, said 3D data set obtained from imaging a sample using an imaging device, said sample having at least one boundary surface associated therewith, said method comprising the steps of:

identifying a subset of the image data corresponding to the boundary surface;

fitting a substantially smooth reference surface to the subset of image data;

calculating the distance between points on the reference surface to points of interest within the sample identified from the 3D image data set; and

generating a 2D data set of the calculated distances, wherein the identifying, fitting, calculating and generating steps are performed in a processor.

13. A method as recited in claim 12, further including the step of generating an image map from the 2D data set.

14. A method as recited in claim 12, wherein the 3D image data set is obtained using an OCT system.

15. A method as recited in claim 14, wherein the sample is an eye.

16. A method as recited in claim 15, wherein the points of interest correspond to points on said boundary surface.

17. A method as recited in claim 15, wherein the points of interest correspond to points on a different boundary surface.

18. A method as recited in claim 15, wherein the boundary surface is the cornea.

19. A method as recited in claim 15, wherein the boundary surface is the retina.

20. A method as recited in claim 12, wherein regions of the boundary surface that are significantly impacted by disease are omitted from the fitted surface to minimize perturbations of the fitted surface.

21. A method for generating image maps from a 3D image data set, said data set obtained from imagine a sample using an imaging device, said sample having at least one boundary surface associated therewith, said method comprising the steps of:

identifying a subset of the image data corresponding to a first boundary surface;

fitting a substantially smooth reference surface to the subset of image data

identifying a second boundary surface;

calculating the distance between points on the reference surface to corresponding points on the second boundary surface;

generating an image map of the calculated distances; and

displaying the image map.

22. A method as recited in claim 21, wherein the first and second boundary surfaces are the same.

23. A method as recited in claim 21, wherein the first and second boundary surfaces are different.

24. A method as recited in claim 21, wherein the 3D image data set is obtained using an OCT system.

25. A method as recited in claim 24, wherein the sample is an eye.

26. A method as recited in claim 25, wherein the boundary surface is the cornea.

27. A method as recited in claim 25, wherein the boundary surface is the retina.

28. A method as recited in claim 21, wherein regions of the first boundary surface that are significantly impacted by disease are omitted from the fitted surface to minimize perturbations of the fitted surface.

29. A method for generating an elevational image map from a 3D image data set, said data set derived from an OCT scan of an eye, said eye having at least one boundary surface associated therewith, said method comprising the steps of:

identifying a subset of the image data corresponding to the boundary surface;

fitting a substantially smooth reference surface to the subset of image data;

calculating the elevation between points on the reference surface to points of interest within the sample identified from the 3D data set;

generating an elevational image map from the calculated elevational points; and

displaying the image map.

30. A method as recited in claim 29, wherein the points of interest correspond to points on said boundary surface.

31. A method as recited in claim 29, wherein the points of interest correspond to points on a different boundary surface.

32. A method as recited in claim 29, wherein regions of the boundary surface that are significantly impacted by disease are omitted from the fitted surface to minimize perturbations of the fitted surface.

33. A method as recited in claim 29, wherein the fitting step is performed by one of a second-order polynomial fit, Zernike polynomials, Chebyshev polynomials, Bessel functions.

34. A method as recited in claim 29, wherein the subset of image data is fit to portion of a sphere.

35. A method as recited in claim 29, wherein the fitting step is performed using a low-order polynomial fit in x and y.

37. A method as recited in claim 29, wherein the elevational points are calculated between the reference surface and the retinal pigment epithelium.

38. A method as recited in claim 29, wherein the elevational points are calculated between the reference surface and the inner limiting membrane.

39. A method as recited in claim 29, wherein the elevational information is presented as a contour map.

40. A method as recited in claim 29, wherein the elevational information is presented as a color map.

41. A method as recited in claim 29, wherein the elevational information is presented as 3D surface rendering.

42. A method as recited in claim 29, further including the step of generating secondary image information related to a structural feature of the eye and wherein said map comprises both the elevational points and the secondary image information.

43. A method as recited in claim 42, wherein the elevational points are presented in a form selected from the group consisting of contours, color variations, grey scale, 3D surface rendering and wherein the secondary image information is presented in a form selected from the same group as the elevational points but using a different form than used for the elevational points.

44. A method as recited in claim 42, wherein the elevational points are presented as a color map and the secondary image information is shown as a 3D surface rendering.

45. A method as recited in claim 42, wherein the elevational points are presented as 3D surface rendering and the secondary image information is presented as a color map.

46. A method as recited in claim 42, wherein the secondary image information corresponds to one of the elevation between actual boundaries, actual thickness of a tissue layer, image signal strength, birefringence and data from a fit to a different reference surface.

47. A method as recited in claim 42, wherein the elevational points are presented as a color map and said secondary image information is image signal strength and is presented as variations in brightness in the color image.

48. A method as recited in claim 42, wherein the image map comprises a contour map having variations in color superimposed thereon with the elevation points corresponding to either the contour or color variations and with the secondary image information corresponding to remaining image variable.

49. A method as recited in claim 42, wherein the boundary surface is the RPE and wherein the secondary information corresponds to the distance from the ILM to the RPE and is displayed in the form of a contour map and wherein the elevational points correspond to the distance from ILM to the reference surface and are displayed in the form of a color map superimposed on the contour map.

36. A method as recited in 29, claim wherein the reference surface is created by fitting to 3D data corresponding to the retinal pigment epithelium.

50. A method for generating an elevational map from a 3D image data set, said data set derived from an OCT scan of an eye, said eye having at least one boundary surface associated therewith, said method comprising the steps of:

identifying a subset of the image data corresponding to the boundary surface;

fitting a substantially smooth reference surface to the subset of image data;

calculating the elevation between points on the reference surface to points of interest within the sample identified from the 3D data set; and

generating a 2D set of the calculated elevational points, wherein the identifying, fitting, calculating and generating steps are performed in a processor.

51. A method for generating image maps from a 2D image data set, said data set obtained from imaging a sample using an imaging device, said sample having at least one boundary surface associated therewith, said method comprising the steps of:

identifying a subset of the image data corresponding to the boundary surface;

fitting a substantially curved reference line to the subset of image data;

calculating the distance between points on the reference line to points of interest within the sample identified from the 2D image data set;

generating an image map of the calculated distances; and

displaying the image map.

52. A method as recited in claim 51, wherein the image data set includes 3D information.

53. A method for generating images from a 2D image data set, said data set obtained from imaging a sample using an imaging device, said sample having at least one boundary surface associated therewith, said method comprising the steps of:

identifying a subset of the image data corresponding to the boundary surface;

fitting a substantially curved reference line to the subject of image data;

calculating the distance between points on the reference line to points of interest within the sample identified from the 2D image data set; and

generating a 2D data set of the calculated distances, wherein the identifying, fitting, calculating and generating steps are performed in a processor.

54. A method as recited in claim 53, wherein the image data set includes 3D information.

55. A method as recited in claim 54, further including the step of generating an image map from the 2D data set.