System and method for low coherence broadband quadrature interferometry

• US 20040239943A1
• Filed: 05/28/2004
• Published: 12/02/2004
• Est. Priority Date: 05/30/2003
• Status: Active Grant

• Current Assignee: Duke University
• Original Assignee: Duke University
• Inventors: Izatt, Joseph A., Yang, Changhuei, Choma, Michael
• Granted Patent: US 7,019,838 B2
• Time in Patent Office: Days
• US Class Current: 356/479
• CPC Class Codes: G01B 2290/45 Multiple detectors for detecting interferometer signals
• CPC Class Codes: G01B 2290/70 Using polarization in the interferometer
• CPC Class Codes: G01B 9/02002 using two or more frequencies
• CPC Class Codes: G01B 9/02004 using frequency scans
• CPC Class Codes: G01B 9/02044 Imaging in the frequency domain, e.g. by using a spectrometer
• CPC Class Codes: G01B 9/02081 simultaneous quadrature detection, e.g. by spatial phase shifting
• CPC Class Codes: G01B 9/02091 Tomographic interferometers, e.g. based on optical coherence
• CPC Class Codes: G01N 21/4795 spatially resolved investigating of object in scattering medium in vivo A61B

Abstract

A quadrature broadband interferometry system and method obtains a complete complex interferometric signal instantaneously in both homodyne and heterodyne systems in a simple, compact, and inexpensive setup. This is accomplished by separating interferometric components from non-interferometric components in each of at least two detector signals of an interferometer having a number of N×N couplers, scaling the interferometric components, and generating real and imaginary parts of a complex interferometric signal from the scaled interferometric components. The detector signals preferably derive from a broadband light source coupled to the interferometer, and the number of N×N couplers may be one or more, with N≧3.

Claims

1. A signal processing method, comprising:

separating interferometric components from non-interferometric components in each of at least two detector signals of an interferometer having a number of N×

N couplers;

scaling the interferometric components; and

generating real and imaginary parts of a complex interferometric signal from the scaled interferometric components.

2. The method of claim 1, wherein the detector signals derive from a broadband light source coupled to the interferometer.

3. The method of claim 1, wherein scaling the interferometric components includes phase shifting the interferometric components of the two detector signals to have a predetermined phase difference.

4. The method of claim 3, wherein the interferometric components are shifted to have orthogonal phases.

5. The method of claim 1, wherein scaling the interferometric components is performed to reduce optical losses.

6. The method of claim 5, wherein the optical losses correspond to at least one of optical path losses and variations in gain between different ones of the detectors.

7. The method of claim 1, wherein separating the interferometric components includes performing high-pass filtering of signals output from the detectors.

8. The method of claim 1, wherein separating the interferometric components includes extracting AC-coupled signals from signals output from the detectors.

9. The method of claim 8, wherein the AC-coupled signals are extracted by:

removing at least one of DC components, source noise components, and auto-correlation terms from the detector outputs prior to re-scaling.

10. The method of claim 1, wherein the imaginary part of the complex interferometric signal is computed based on a method derived from the cosine sum rule.

11. The method of claim 10, wherein the cosine sum rule-derived method is based on the following equation:

12. The method of claim 11, wherein scaling the interferometric components includes adjusting φ

_m−

φ

_nto equal a predetermined angle.

13. The method of claim 12, wherein the predetermined angle is at least substantially 90°

.

14. The method of claim 1, further comprising:

calculating a magnitude and phase from the real and imaginary parts of the complex interferometric signal.

15. The method of claim 14, further comprising:

calculating a Doppler shift value based on a time derivative of the calculated phase.

16. The method of claim 15, wherein the Doppler shift value is calculated based on the following equation:

18. The method of claim 1, wherein the number of N×

N couplers is one.

19. The method of claim 18, wherein N≧

3.

20. The method of claim 19, wherein the detector signals are at least partially derived from outputs of the N×

N coupler.

21. The method of claim 20, wherein scaling the interferometric components includes:

adjusting the interferometric components output from the N×

N coupler to have predetermined phases, wherein the complex interferometric signal is derived from the—

phase-adjusted interferometric components.

22. The method of claim 21, wherein the interferometric components of the N×

N coupler are adjusted to have orthogonal phases.

23. The method of claim 1, further comprising:

generating a plurality of complex interferometric signals from scaled interferometric components derived from different pairs of detectors in the interferometer; and

generating an average complex interferometric signal from said plurality of complex interferometric signals.

24. The method of claim 1, wherein the interferometer includes a first coupler connected to a second coupler by an optical fiber, and wherein the two detector signals respectively correspond to different outputs of the first and second fiber-optic couplers.

25. The method of claim 23, wherein the first coupler is an N×

N coupler and the second coupler is an M×

M coupler, where N≠

M.

26. The method of claim 23, wherein the first coupler is an N×

N coupler and the second coupler is an M×

M coupler, where N=M.

27. The method of claim 1, wherein the interferometer includes an optical circulator connected to one of the fiber-optic couplers, and wherein the two detector signals respectively correspond to outputs of the fiber-optic coupler and optical circulator.

28. The method of claim 1, wherein separating the interferometric components includes:

measuring DC levels of the detector outputs prior to re-scaling; and

subtracting the DC levels from the detector outputs respectively.

17. The method of claim 17, further comprising:

estimating reflector velocity based on

$\frac{<br><br>\mathrm{phase}}{<br><br>t}.$

29. An analyzer, comprising:

an interferometer including at least one fiber-optic coupler;

a number of detectors that detect light from said at least one coupler; and

a processor which acquires interferometric components from signals output from the detectors, scales the interferometric components, and generates real and imaginary parts of a complex interferometric signal from the scaled interferometric components.

30. The analyzer of claim 29, wherein the light from said at least one coupler derives from a broadband light source coupled to the interferometer.

31. The analyzer of claim 29, wherein said at least one coupler is an N×

N coupler, where N≧

2.

32. The analyzer of claim 31, wherein N≧

3.

33. The analyzer of claim 29, wherein the interferometer includes at least a first N×

N coupler and a second M×

M coupler, and wherein the detector signals are generated from light output from the first and second couplers respectively.

34. The analyzer of claim 33, wherein N=M.

35. The analyzer of claim 34, wherein N≠

M.

36. The analyzer of claim 29, wherein, during scaling, the processor phase shifts the interferometric components of the detector signals to have a predetermined phase difference.

37. The analyzer of claim 36, wherein the interferometric components are shifted to have orthogonal phases.

38. The analyzer of claim 29, wherein the processor scales the interferometric components to reduce optical losses.

39. The analyzer of claim 38, wherein the optical losses correspond to at least one of optical path losses and variations in gain between different ones of the detectors.

40. The analyzer of claim 29, wherein the processor separates the interferometric components of each of the detector signals from non-interferometric components.

41. The analyzer of claim 40, wherein the processor separates the interferometric components by performing high-pass or band-pass filtering of the detector signals.

42. The analyzer of claim 40, wherein the processor separates the interferometric components by extracting AC-coupled signals from the detector signals.

43. The analyzer of claim 42, wherein the processor extracts the AC-coupled signals by removing at least one of DC components, source noise components, and auto-correlation terms from the detector outputs prior to re-scaling.

44. The analyzer of claim 29, wherein the processor computes the imaginary part of the complex interferometric signal based on a cosine sum rule.

45. The analyzer of claim 44, wherein the cosine sum rule is based on the following equation:

46. The analyzer of claim 45, wherein the processor scales the interferometric components by adjusting φ

_m−

φ

_n to equal a predetermined angle.

47. The analyzer of claim 46, wherein the predetermined angle is at least substantially 90°

.

48. The analyzer of claim 29, wherein the processor calculates a magnitude and phase from the real and imaginary parts of the complex interferometric signal.

49. The analyzer of claim 48, wherein the processor calculates a Doppler shift value based on a time derivative of the calculated phase.

50. The analyzer of claim 49, wherein the Doppler shift value is calculated based on the following equation:

51. The analyzer of claim 29, wherein the processor generates a plurality of complex interferometric signals from scaled interferometric components derived from different pairs of the detectors, and generates an average complex interferometric signal from said plurality of complex interferometric signals.

52. A method for performing optical coherence tomography, comprising:

generating a complex interferometric signal; and

instantaneously recovering Doppler shift information from the complex interferometric signal.

53. The method of claim 52, wherein the recovering step includes:

recovering, from the complex interferometric signal, a sign and magnitude of a Doppler shift arising from moving scatterers in a sample scanned by light from the interferometer.

54. The method of claim 52, wherein the generating step includes:

acquiring interferometric components from signals output from at least one fiber-optic coupler located along an optical signal path of the interferometer;

scaling the interferometric components; and

generating real and imaginary parts of a complex interferometric signal from the scaled interferometric components.

55. The method of claim 54, wherein the light from said at least one coupler derives from a broadband light source coupled to the interferometer.

56. The method of claim 54, wherein said at least one coupler is an N×

N coupler, where N≧

2.

57. The method of claim 56, wherein N≧

3.

58. The method of claim 52, wherein the interferometer includes at least a first N×

N coupler and a second M×

M coupler, and wherein the detector signals are generated from light output from the first and second couplers respectively.

59. The method of claim 58, wherein N≠

M.

60. The method of claim 54, wherein the generating step includes:

generating the imaginary part of the complex interferometric signal based on a cosine sum rule.

61. The method of claim 60, wherein the cosine sum rule is based on the following equation:

62. The method of claim 61, wherein the first and second detectors detect light along different optical signal paths of the N×

N coupler.

63. The method of claim 61, wherein the processor scales the interferometric components by adjusting φ

_m−

φ

_n to equal a predetermined angle.

64. The method of claim 63, wherein the predetermined angle is at least substantially 90°

.

65. The method of claim 52, further comprising:

calculating a magnitude and phase from the real and imaginary parts of the complex interferometric signal.

66. The method of claim 65, wherein the recovering step includes obtaining a Doppler shift value based on a time derivative of the calculated phase.

67. The method of claim 65, wherein the Doppler shift value is calculated based on the following equation: