超小自聚焦光纤探头的研究进展

王驰;旷滨;孙建美;朱俊;毕书博;蔡楹;于瀛洁

【摘 要】梯度折射率(Gradient-index,GRIN)光纤探头是一种全光纤型超小光学镜头,在心血管等狭小空间组织内窥影像检测中具有广阔的应用前景.但其发展一直缺少系统的理论体系.本文讨论探头设计、制作和性能测试等方面的关键问题.基于GRIN光纤探头聚焦性能的特征参数,对解析设计方法与数值仿真设计方法进行比较分析.针对超小GRIN光纤探头的制作难题,研究一种光纤熔接和切割的高精度一体化集成装置,描述GRIN光纤探头的制作方法.此外,分析了超小GRIN光纤探头聚焦性能检测的方法及装置.本文为超小GRIN光纤探头的设计、制作及性能测试提供了一个方法体系.

【期刊名称】《中国光学》 【年(卷),期】2018(011)006 【总页数】14页(P875-888)

【关键词】光学相干层析;微小光学探头;自聚焦光纤镜头;GRIN光纤探头;光纤熔接与切割

【作 者】王驰;旷滨;孙建美;朱俊;毕书博;蔡楹;于瀛洁

【作者单位】上海大学精密机械工程系,上海200072;近地面探测技术重点实验室,江苏无锡214035;上海大学精密机械工程系,上海200072;上海大学精密机械工程系,上海200072;近地面探测技术重点实验室,江苏无锡214035;上海大学精密机械工程系,上海200072;近地面探测技术重点实验室,江苏无锡214035;上海大学精密机械工程系,上海200072

【正文语种】中 文

【中图分类】TP394.1;TH691.9

1 Introduction

Optical coherence tomography(OCT[1-3]) is becoming increasingly promising in the fields such as biomedical and surgical operations due to its advantages including rapid imaging speed, high resolutions, and non-contact detection, etc.. Following the X-ray computed tomography(XCT) and magnetic resonance imaging(MRI), OCT is another tomography technology which can be used to obtain the internal structure information of the samples by detecting interference signal. In OCT system, optical probe is a key component that transmits and focuses light beam(e.g. Gaussian beam) inside a sample of interest, and then collects the reflected or scattered light carrying information about the sample and sends the information to the signal processing system. Its focusing performance plays an important role in determining the imaging quality of the OCT system. For example, the focal waist location and spot size of the probe respectively determine the penetration depth and lateral resolution to a large extent. However, most biological tissues are optically nontransparent, so the detection depth of OCT is limited, generally in the range of 1-3 mm[4]. The development of small endoscopic optical probes has become an important growing branch of OCT technology.

Gradient-index(GRIN) fiber probe is defined to be as an all-fiber-type

ultra-small optical probe, which typically consists of a single-mode fiber(SMF), a no-core fiber(NCF), and a GRIN fiber lens. It could be suitable for the imaging detection in deep, narrow tissues or organs such as cardiovascular system, and has been investigated by some researchers in recent years. Lin[5] presented a novel application approach in the study of two-segment lensed fiber collimator. Swanson[6] was granted an invention patent for ultra-small GRIN fiber probe. Then, Reed and Jafri et al. developed in-vivo OCT imaging system [7-8]. Since 2007, Dr. Mao et al. studied the fabrication method and performance testing method of such GRIN fiber probe[9-10]. Jung[11]analyzed a miniaturized OCT probe model, comprising of a “SMF+NCF+GRIN lens”, with the ABCD matrix algorithm for a Gaussian beam. Lorenser[12]modeled the ultra-small GRIN fiber probe with the Beam Propagation Method(BPM). In recent years, the group headed by Prof. Sampson D. D. at the University of Western Australia has obtained application solutions of breast cancer and

preliminary testing results of optical image with the OCT system based on GRIN fiber probe[13]. In addition, Schmitt R. studied the testing scheme of micro-deep holes by using endoscopic detection system based on this probe[14-15] .

Ultra-small GRIN fiber probe has many advantages on endoscope detection. However, the design theory of GRIN fiber probe has been largely overlooked or only simplified in the above studies. In addition, this kind of probe is difficult to fabricate due to its ultra-small size. Since 2011, our research group has researched the design theory of the GRIN fiber

probe through analytical method and numerical simulation technology[16-22]. And the fabrication method of the fiber probe based on fiber cutting and fiber welding were preliminary studied[23]. These design methods published previously are described in this paper. Advantages of each method are presented. Based on the existing researches, we have recently developed a high precision optical integration machine with functions of fiber welding and fiber cutting. In addition, detection method for the focusing performance of the probe is overviewed together with the coupling efficiency testing method. 2 Optical model of the GRIN fiber probe

Fig.1 displays a typical model of GRIN fiber probes, which is composed of a single-mode fiber(SMF), a no-core fiber(NCF), and a GRIN fiber lens. The SMF is connected with the detection arm of the OCT system, and guides the light beam into the NCF. The NCF is a kind of special fiber with homogeneous refractive index. Adding the NCF is able to improve the focusing performance of the probe by expanding the beam and thus to overcome the problem of limited mode field diameter of the SMF.

However, the length of NCF should be chosen appropriately. On one hand, when the NCF is too long, overflowing of some light energy may exist off the probe side which will finally reduce beam coupling efficiency. On the other hand, when the NCF is too short, it may result in the failure of expanding beam and improving the focusing performance of the probe. GRIN fiber lens is the most critical part of the probe that performs a self-focusing function due to the continuous change of refractive index. When

the length of the GRIN fiber lens goes close to 1/4 pitch(or the integer times), the lens will get a really short focal length for strong focusing performance; When the length is close to 1/2 pitch(or the integer times), a longer focal length will be got with a bigger spot size. In the study of the imaging performance of an OCT system, a longer focal length is needed in order to obtain a deeper detection depth. In addition, focusing spot size is as small as possible to get a higher lateral resolution. Therefore, a tradeoff is adopted to determine the length of GRIN fiber lens. An ideal design for ultra-small GRIN fiber probe can be obtained by using a NCF spacer with the characteristics of beam-expanding to improve the focusing

performance. The length of NCF and the length of GRIN fiber lens are in the submillimeter range; therefore it is a challenging problem to achieve high precision cutting and welding of such short fiber components. It is noted that, the refractive indices of the NCF and the core of GRIN fiber at the axis should match that of the SMF at the core, in order to minimize the influence of the complicated reflection of input beam by different optical interfaces on OCT imaging quality. Fig.1 Model of GRIN fiber probe 3 Design methods of GRIN fiber probe 3.1 Analytical design method

The optical characteristic parameters is the theoretical basis for design of GRIN fiber probe, which are defined as in Ref.[16]. Working distance is the distance between the output end of the probe and the focal plane, characterizing the detection depth of the OCT system. Focused spot size is

the waist diameter of Gauss beam focused by the probe, characterizing the lateral resolution of the OCT system. Depth of field is twice of the Rayleigh range of the Gauss beam, characterizing the range of the detection depth of the OCT system. A brief overview is described below. As shown in Fig.1, an approximate Gaussian beam is output from the SMF with a wavelength of λ and a radius of ω0, then transmits into the NCF with an index of n1, and then transmits into the GRIN fiber lens with an index profile described by the formula (1). The light beam is eventually focused into a spot with a waist radius of W and a waist location of Zw in air. The refractive index of air is n2. The refractive index at the center of the GRIN fiber lens is n0. L0 and L represent the length of the NCF and the length of the GRIN fiber lens, respectively. Planes 1 to 6 denote the input plane, two interfaces between the NCF and the GRIN fiber lens, two interfaces between the GRIN fiber lens and the air, and the focal plane, respectively. The refractive index of a GRIN fiber lens can be described as the following quadratic equation. (1)

Where r is the radius, and g is the gradient constant. Assigning a=λ/n1πw2, according to Ref.[16], the expression of the working distance can be described as formula (2) by adopting the complex beam parameter matrix transformation method. (2)

Where

S1=-2n1n2L0a2

S4=2n0n1gL0a2

The expressions of the spot size can be described as formula (3). (3) Where

The expressions of the Rayleigh range can be described as (4): (4)

According to formulas above, we can comprehensively analyze the relationship between the focusing performance and the structural parameters of the GRIN fiber probe, the results can be used to optimize ultra-small GRIN fiber probes with specific optical performance requirements. For example, MATLAB is adopted to further solve these expressions of the characteristic parameters by means of drawing the function plots between the characteristics and the length of the fiber spacer L0 and the length of GRIN fiber lens L. As shown in Fig.2, Fig.2(a)

illustrates the contour relationship plots of working distance and spot size. Fig.2(b) shows the relationships between working distance and spot size and the length of NCF L0, given a constant length of the GRIN fiber lens L=0.1 mm. Fig.2(c) shows the relationships between working distance and spot size and the length of a GRIN fiber lens L, given a constant NCF length L0=0.36 mm. According to Fig.2, it is obvious that NCF can improve the focusing performance of the GRIN fiber probe and increase the working distance while still keeping the focusing spot small. In addition, the working distance shows periodicity with the increasing of the length of the GRIN fiber lens, which depends on the pitch length of the lens. Fig.2 Relationship between characteristics and the lengths of probe components

Generally speaking, one tends to design an OCT system with working distance and depth of field as long as possible and lateral resolution as high as possible, respectively. In addition, the lateral resolution decreases with the increase of spot size. Therefore, the length of the NCF and the length of the GRIN fiber lens should be chosen to achieve a trade-off between the working distance and the lateral resolution of an OCT imaging system in order to optimize the design of GRIN fiber probe. As a special and useful case of interest, according to Fig.2(b), when the NCF length ranges from 0.32 to 0.4 mm, given the length of GRIN fiber lens L=0.1 mm, the working distance is located into 0.5 to 0.76 mm, and correspondingly the spot size is located into 26 to 40 μm. Therefore, in order to meet such different requirements of the OCT imaging probe, the

use of NCF length, 0.32 to 0.4 mm, may be a good tradeoff, which means that the working distance is greater than 0.5 mm and the spot size less than 40 μm.

From Fig.2, we find another important conclusion that the characteristic parameter can change sharply within some length ranges of the NCF or GRIN fiber lens. Since the length of the GRIN fiber probe is less than 1 mm, any small error in either the length of the NCF or that of the GRIN fiber lens during fabrication can have a great impact on its optical focusing performance. In Ref.[21], the partial derivatives of the mathematical expressions and 2-D curvature graphs of characteristics are adopted to analyze this mutation in detail when designing the lengths of the NCF and the GRIN fiber lens. The selected lengths of the two components of the GRIN fiber probe should meet the specific performance and avoid obviously mutational range.

3.2 Numerical design method by using GLAD

As presented above, analytical design method can help analyze and design such ultra-small GRIN fiber probe with special demands of optical properties by means of the mathematical expressions of characteristic parameters. The use of NCF can improve the focusing performance of GRIN fiber probe by means of its effect of expanding beam. However, it is difficult to obtain the beam profile at locations of interest, e.g. input plane, two interfaces between the NCF and the GRIN fiber lens, two interfaces between the GRIN fiber lens and the air, and the focal plane, etc.. The beam parameters and the beam profile within the ultra-small GRIN fiber

probe path are very difficult or even impossible to measurement by experimental method due to its very limited size. In order to solve such issues, the numerical design method was proposed to model GRIN fiber probes in Ref.[18] by using the optical software of GLAD which has the capability of modeling almost any type of physical optical systems. GLAD treats optical beams as complex amplitude distribution, and gives a much more powerful analysis capability. It models the GRIN fiber probe by assigning relative parameters of light source and optical components. Take a special case of the GRIN fiber probe as an example, Fig.3 shows the full profile of all points along the optical path through the probe with the length of NCF L0=0.36 mm and the length of GRIN L=1.24 mm. Fig.3(a) displays beam profile as a function of distance from the surface of the SMF to a total distance of 3.52 mm. Fig.3(b) demonstrates beam width vs. axial position. From Fig.3, we can acquire the best focus position at about 2.35 mm, and thus the working distance 0.75 mm, and the spot size 36 μm. It is obvious that the beam is expanded through the NCF, then focused and thus narrowed down through the GRIN fiber lens.

Fig.3 Beam profile of all points along the optical path going through the probe

3.3 Numerical design method by using VirtualLab

GLAD can help obtain the beam profile at all locations along the propagation direction in addition to the internal probe. However, it is a kind of programing language and is very difficult for a general user to get a quick grasp. And the design process of GRIN fiber probe is very

complicated. Therefore, the authors proposed another numerical design method by means of the commercial optical software VirtualLab as described in Refs.[19-20]. VirtualLab is a unified optical modeling platform and numerical analysis software for physical optics, designed by German LightTrans Company. Based on the theory of electromagnetic field, it can be used to model all light sources and arbitrary transmission by the usage of field tracking. Model is built in the form of Light Path Diagram(LPD). The light source, components and detectors used in the system, can be called from the optical element library and their optical characteristic can be edited flexibly.

Fig.4 Light path diagram of GRIN fiber probe

Fig.4 is the LPD to model GRIN fiber probe. The element “Gaussian Wave” represents the input Gaussian beam. The output of SMF can be considered as the incident position of the Gaussian beam. Optical element “NCF” represents the coreless fiber. Optical element “GRIN fiber lens” indicates the self-focus fiber lens. Optical element “Detector” indicates the beam parameter detector for detecting the waist size of Gaussian beam. Optical element “Virtual Screen” is to observe the intensity distribution of Gaussian beam. By double-clicking the corresponding symbol, we can flexibly change parameters to complete the configuration of the light source and optical elements. Fig.5(a) shows the intensity distribution of Gaussian beam at the output end of a GRIN fiber probe by setting the length of NCF as 0.36 mm and the GRIN fiber length as 0.11 mm. Fig.5(b) shows the plot of beam radius versus the axial position. It can be

monitored that the Gaussian beam is focused and then gradually diverged at the exit end of the probe. The position with minimal beam radius along the Z-axis is the waist position of the output Gaussian beam. The radius of the waist position is spot size.

Fig.5 Beam focal performance at different locations along the propagation path

4 Fabrication of GRIN fiber probe

Fabricaion of GRIN fiber probe is a very difficult and important task due to its limited size, it needs high-precision manufacturing devices and methods. In order to solve this problem, we have recently developed an all-in-one integration fabrication device of ultra-small GRIN fiber probe as shown in Fig.6, which is made of fiber welding unit and fiber cutting unit. The characteristic of this encapsulation device is to achieve the integration of fabrication device of ultra-small GRIN fiber probe. The steps of fiber cutting and fiber welding are operated on the same device during fabricating the ultra-small GRIN fiber probe, which further simplifies the fabrication process and improves the fabrication efficiency of ultra-small fiber probe.

Fig.6 All-in-one integration device for fiber welding and cutting This machine can cut and weld bare fibers with high cutting precision about 5 μm. The length of NCF and GRIN fiber can also be measured by this device. And the measurement accuracy is about 1 μm. The fabrication process of the GRIN fiber probe is shown in Fig.7, and the main steps are as follows:(1)Fuse the NCF to the single-mode optical fiber; (2)Cut the NCF

to certain length as fiber spacer by taking the fusion point A as the origin between SMF and NCF; (3)Fuse the GRIN fiber lens to the fiber spacer; (4)Cut GRIN fiber to pre-calculated length as the focusing lens by taking the welding joint B as the origin between the fiber spacer and the GRIN fiber lens. Through the above methods, we have fabricated six groups with different sizes of the probe and measured the lengths of probe

components with high magnification rate microscopic, as shown in Tab.1. Fig.7 Process of fabricating a GRIN fiber probeTab.1 Preset length and the measured length of probe components

GroupNCF length/mmGRIN fiber lens length/mmPreset

length100.360200.41030.1600.20040.2400.14050.3600.11060.3000.150Measured

length100.356200.40730.1640.20040.2420.14150.3560.10860.3000.146 5 Properties detection methods of GRIN fiber probe

The miniaturization of dimension and focusing spot is the main advantage of GRIN fiber probe, however, it is difficult to detect the focusing performance of the probe with high precision and speed by traditional beam detection methods, i.e. the mechanical scanning and array detection method, owing to the low efficiency and poor real-time performance in mechanical scanning method[24-25] and the resolution limited by pixel sizes in array detection method[26-27]. Based on the characters of focusing performances, this paper describes a non-contact detection method with infinity optical transformation technology as presented in Ref.[28]. The accurate position of the probe is realized through the precise

adjustment mechanism. Finally, the focusing performances are obtained by using curving fitting method with measured data captured from the measurement system.

The working distance of GRIN fiber probe generally does not exceed 1 mm, and the focus spot size is about 40 μm. The above characters indicate that the focusing position is too close to the output plane of lens with a micro focus spot. Hence, during the measurement of focusing performances, the output plane of the lens should avoid direct contact with the detector, due to the limitation of dimension and beam parameters of the probe. Fig.8 presents theproperty testing method. The distribution of light intensity at various distances along the direction of propagation after the lens is captured through changing the distance between the lens and the detector along the optical axis. Then, the working distance and the spot size are calculated from the measured intensity distribution.

Fig.8 Schematic diagram of the property testing method for the GRIN fiber probe

In order to achieve the reliable mind about detecting the focusing performance of the micro fiber lens as shown in Fig.8, an infinity optical transformation system is presented in Fig.9 based on precise adjustment mechanism. The center of the light from a source through the lens is aligned with the center of the detector by accurately adjusting the lens position. Then, the light is captured with magnification by the infinity optical transformation system in a non-contact mode. The CCD camera converts the optical signals to electrical signals which are then transmited

into a computer. Finally, the computer software analysis the beam information and provides illuminating displays for beam parameters, such as beam width, peak power, and intensity distribution, etc.

Fig.9 Schematic diagram of the focusing performance measurement of the GRIN fiber probe

The measurement system corresponding to Fig.9, can be composed of a Beam Analyzer USB(Duma Optronics Ltd.) with 1.0 μm resolution, a superluminous diode source(SLD-1310-18, Fiberlabs), a microscopic objective lens(M Plan Apo NIR 10X, Mitutoyo), and a precise adjusting platform, etc. In the process of detection, the small fiber probe is fixed on the adjusting platform. The light source is connected to the lens by single-mode fiber, and the CCD detector will capture the magnified beam under the effect of the optical transformation system for analyzing the beam parameters. Before measuring the beam profile from the lens, we use the precise adjustment mechanism(attached a fiber jig) to adjust the position of the probe and make exit beam vertical to the CCD receiving plane. Finally, in order to collect the beam parameters at various distances along the direction of propagation, the probe moving to change the distance between the probe and the detector by adjusting the platform.

Tab.2 shows the properties comparison between the testing data and the calculation data using different design methods. It is obvious that, the calculation results of working distance and spot size are almost the same as those experimental data. Therefore, the presented design methods for investigating GRIN fiber probe are feasible and effective. In terms of the

differences between the experimental and calculation results, there could be the following reasons. Firstly, the cutting lengths of NCF and GRIN fiber lens cannot be completely precisely consistent with the pre-calculated lengths, results in errors of beam expanding and focusing between the experimental and calculation data. Secondly, the actual indices of the center of GRIN fiber lens and the NCF as well as the core of SMF are not same, which leads to beam back reflection between different interfaces. However, this phenomenon is not taken into account in the calculation programe. Thirdly, the measurement error of the experimental system is another reason. Moreover, there may be some other factors to be investigated.

Tab.2 Properties comparison between the testing values and the simulating dataTypesSamplesNCF/mmGRIN fiber lens/mmWorking distance/mmSpot size/μmExperimental 10.360.100.7531results[9-10]20.360.110.602930.360.120.5026Analytical

calculation10.360.100.7532.4results20.360.110.6323.730.360.120.5118.4Numerical calculation10.360.100.7333results using GLAD20.360.110.642830.360.120.5225Numerical calculation10.360.100.7532results using Virtual Lab20.360.110.632430.360.120.5119

Fig.10 Detection scheme of the coupling efficiency Fig.11 Detection system of the coupling efficiency

In addition, the coupling efficiency of GRIN fiber probes, which reflects the comprehensive focusing performance of probes, is an important

characteristic parameter for determining the detection performance of OCT system and it affects the sensitivity of OCT detection system. In Ref.[29-35], a theoretical equation for the coupling efficiency of probes is derived using an analytical approach based on the optical model of ultra-small GRIN fiber probe and transmission characteristics of the Gaussian beam. Variation and influencing factors were analyzed and verified by establishing the corresponding experiment system for testing. The detection scheme of coupling efficiency was investigated by the theoretical calculation method of coupling efficiency of the probe, and Fig.11 shows the actual detection system corresponding to Fig.10. The output and receiving ends of the probe with fiber clamps are placed on two five-dimensional adjustment platforms and fixed to the optical isolation platform. The output and receiving ends of the probe are connected to the output laser and power receiver respectively. The microscope is used to adjust the relative position of the two fiber probes to ensure that the central axes of the two probes overlap with each other. During the experiment, the fixed five-dimensional adjustment platform is used for placing the output end of the probe, and that for placing the receiving end of the probe is used to change the axial distance between the two probes. During the experiment, the adjustment platform for placing the receiving fiber probe is gradually adjusted to only change the axial distance between fibers. The corresponding power value received by the power meter at every position is recorded until the power value drops near zero. After recording all the power values, the coupling efficiency is

calculated to analyze the relationship between the coupling efficiency and the axial distance of the probe. It is noted that, the actual coupling efficiency of ultra-small GRIN fiber probe was investigated by introducing the energy loss factor and position error factor and a method was proposed to estimate the energy loss factor based on the simulated annealing algorithm. The result indicates that the ultra-small GRIN fiber probe has a superior focusing performance, and the mechanism of energy loss of the probe itself should be further explored in the future. 6 Conclusions

Study of the gradient-index lens based small optical probes is an

important topic for the miniaturization of OCT systems[36-40]. GRIN fiber probe is an ultra-small optical probe that may be used for the imaging of small lumen, narrow space in the deep tissues and organs(e.g.

Cardiovascular) of human beings and small animals. In this paper, some different design methods are described. The calculating data are well agreement with the previously published experimental results and thus demonstrate the effectiveness of the design methods. Therefore, the authors argue that these design methods can be used as a design theory of GRIN fiber probes[41-46].

From this paper, we give the other conclusions. Different analytical and numerical methods can solve the problem of the design of the GRIN fiber probe, and each one enjoys its own advantages. The analytical design method, by using the expressions of optical characteristics parameters of GRIN fiber probe, can help analyze the function relationships between the

focusing properties of the probe and its different influencing factors comprehensively. For example, the three-dimensional graphs of characteristics parameters can be adopted to intuitively judge the influence of different lengths of NCF and GRIN fiber lens on the probe focusing performance. The application of GLAD makes it easy to analyze the propagation performance of the beam at all locations going through the probe directly. The utilization of VirtualLab is convenient for probe modeling and parameter setting.

This paper proposed a high precision fiber welding-cutting device with the welding unit and cutting unit. The device is utilized to fabricate several groups of the GRIN fiber probes and measure the lengths of their components with high magnification rate microscope. The result shows that the measured lengths are similar to the preset lengths, which further verifies the designed device requirements about fabrication of the probe, and it can be used in the researches of miniaturized optical probe and OCT system.

In addition, it is very difficult to detect the focusing performance of the probe with high precision and speed by traditional beam detection methods due to the miniaturization of dimension and focusing spot. Based on the characters of focusing performances, this paper describes a non-contact detection method with infinity optical transformation technology. And the detection scheme of coupling efficiency of GRIN fiber probes was described, which reflects the comprehensive focusing performance of probes. As a result, a system of design, fabrication and testing methods is

provided for the miniaturization research of optical probes. It is noted that, GRIN fiber probe can be used to many other fields in addition to endoscopic optical coherence tomography. By combining it with other techniques or systems, such as those described in Refs.[30-46], relevant issues can be explored in different applications in the future. 参考文献：

【相关文献】

[1] HUANG D,SWANSON E A,LIN C P. Optical coherence tomography[J]. Science,1991,254(5035):1178-1181.

[2] FEICHERL A F,DREXLERL W,HITZENBERGERL C K,et al.. Optical coherence tomography-principles and applications[J]. Reports on Progress in Physics,2003,66(2003):239-303. [3] JIANG Y,TOMOV I,WANG Y M,et al.. Second-harmonic optical coherence tomography[J]. Optics Letters,2004,29(10):1090-1092.

[4] ROLLINS A M,KULKARNI M D,YAZDANFAR S,et al.. In vivo video rate optical coherence tomography[J]. Opt. Express,1998,3(6):219-229.

[5] LIN Y,CHANG CH. Characteristics of two-segment lensed fiber collimator[J]. Microwave & Optical Technology Letters,2010,52(8):1846-1848.

[6] SWANSON E,PETERSEN C L,MCNAMARA E,et al.. Ultra-small optical probes, imaging optics, and methods for using same:US,US6445939[P]. 2002.

[7] REED W A,YAN M F ,SCHNITZER M J. Gradient-index fiber-optic microprobes for minimally invasive in vivo low-coherence interferometry[J]. Opt. Lett.,2002,27:1794-1796 [8] JAFRI M S,FARHANG S,TANG R S,et al.. Optical coherence tomography in the diagnosis and treatment of neurological disorders[J]. Biomed. Opt.,2005,10(5):051603.

[9] MAO Y X,CHANG S D,SHERIF S,et al.. Graded-index fiber lens proposed for ultrasmall probes used in biomedical imaging[J]. Applied Optics,2007,46(23):5887-5894

[10] MAO Y X,CHANG S D,FLUERARU C. Fiber lens for ultra-small probes used in optical coherent tomography[J]. Biomedical Science and Engineering,2010,3:27-34.

[11] JUNG W,BENALCAZAR W,AHNAD A,et al.. Numerical analysis of gradient index lens based optical coherence tomography imaging probes[J]. Journal of Biomedical

Optics,2010,15(6):066027.

[12] LORENSER D,YANG X,SAMPSON D D. Accurate modeling and design of gradient-index fiber probes for optical coherence tomography using the beam propagation method[J]. IEEE Photonics Journal,2013,5(2):3900015.

[13] MCLAUGHLIN R A,QUIRK B C,CURATOLO A,et al.. Imaging of breast cancer with optical coherence tomography needle probes:feasibility and initial results[J]. IEEE J. Sel. Top. Quantum Electron,2012,18(3):1184-1191.

[14] PFEIER T,SCHMITTR R,KONIG N,et al. Interferometric measurement of injection nozzles using ultra-small fiber-optical probes[J]. Chin. Opt. Lett.,2011,9(7):071202. [15] SCHMITT R,PFEIER T,DEPIEREUX F,et al.. Novel fiber-optical interferometer with miniaturized probe for in-hole measurements[J]. Optoelectron. Lett.,2008,4:140-142. [16] WANG C,MAO Y X,FANG C,et al.. Analytical method for designing gradient-index fiber probes[J]. Optical Engineering,2011,50(9):094202-1-9.

[17] WANG C,MAO Y X,TANG Z,et al.. Numerical analysis of GRIN lens based miniature probes for optical coherence tomography[J]. Opt. Precision Eng.,2011,19(9):2300-2307. [18] WANG C,MAO Y X,TANG Z,et al.. Numerical simulation of gradient-index fibre probe and its properties of light propagation[J]. Chinese Physics B,2011,20(11):114218-1-7. [19] 王驰,毕书博,王利,等.超小自聚焦光纤探头研究用场追迹数值模拟技术[J].物理学报,2013,62(2):024217.

WANG C,BI S B,WANG L,et al.. Field-tracing based numerical simulation technique for the investigation of ultra-small self-focusing optical fiber probe[J]. Acta Physica Sinica,2013,62(2):024217.(in Chinese)

[20] WANG C,BI S B,XIA X Q,et al.. Field-tracing modeling of the ultra-small gradient-index fiber probe[J]. Optik,2013,124:6437-6443.

[21] WANG C,BI S B,XIA X Q,et al.. Further analysis of focusing performance of an ultra-small gradient-index fiber probe[J]. Optical Engineering,2014,53(1):013106. [22] 王驰,许婷婷,毕书博,等.测量自聚焦光纤透镜聚焦常数的曲线拟合算法[J].光学 精密工程,2015,23(12):3309-3315.

WANG C,XU T T,BI S B,et al.. Curve-fitting algorithm of measuring focusing constant of gradient-index fiber lens[J]. Opt. Precision Eng.,2015,23(12):3309-3315.(in Chinese) [23] WANG C,ZHANG F,BI S B,et al.. Fabrication method of ultra-small gradient-index fiber probe[J]. Advances in Manufacturing,2014,2(4):327-332.

[24] KOKODII N G. Knife gauge for measurement of the intensity distribution in a beam of optical radiation[J]. Measurement Techniques,2003,46(3):240-244.

[25] BALLEGAARD H P,ANDERSEN B,BOEGH N S. Multibeam internal drum scanning system:WO,WO 2005012978 A1[P]. 2005.

[26] CUI X Q,XIN H,WU J G,et al.. Slanted Hole Array Beam Profiler(SHArP)-a high-

resolution portable beam profiler based on a slanted linear aperture array[J]. Optics Letters,2006,31(21):3161-3163.

[27] PANG M,RONG J,YUAN X W,et al.. Research on the measurement method for a large laser beam profile based on CCD diffuse transmission imaging[J]. Measurement Science & Technology,2013,24(12):125202.

[28] BI S B,WANG C,ZHU J,et al.. Detection method for the focusing performance of small fiber lens[C]. International Conference on Optical and Photonics Engineering,2017:104491E.

[29] WANG C,SUN J M,SUN F,et al.. Coupling efficiency of ultra-small gradient-index fiber probe[J]. Optics Communications,2017,389:265-269.

[30] 刘震,王雪梅,倪文波.基于MEMS传感器的高精度姿态角测量研究[J].中国测试,2017,43(2):6-12.

LIU Z,WANG X,NI W. Research on attitude angle measurement with high precision based on MEMS sensors[J]. China Measurement & Test,2017,43(2):6-12.(in Chinese)

[31] 程琳,刘影,田彦.基于最小均方误差的全光纤电流互感器调制相位扰动补偿方法研究[J].中国测试,2016,42(12):116-120.

CHENG L,LIU Y,TIAN Y. Research on modulation phase disturbance compensation method for FOCT based on minimum mean square error[J]. China Measurement & Test,2016,42(12):116-120.(in Chinese)

[32] 蔡苹杨,黄承韫,安金玲.基于纳米氧化铟表面的二硫化碳传感器研究[J].中国测试,2013,39(5):42-45.

CAI P Y,HUANG C Y,AN J L,et al.. Research of gas sensor based on cataluminescence emission on the surface of nano-In2O3[J]. China Measurement & Test,2013,39(5):42-45.(in Chinese)

[33] 叶廷东,程韬波,刘桂雄,等.MEMS气敏传感信息的动态预测补偿方法[J].中国测试,2014,40(4):1-5.

YE T D,CHENG T B,LIU G X,et al.. Dynamic predictive compensation method of MEMS gas sensing information[J]. China Measurement & Test,2014,40(4):1-5.(in Chinese) [34] 胡林亭,李佩军,姚志军.提高外场重频激光光斑测量距离的研究[J].液晶与显示,2016,31(12):1137-1142.

HU L T,LI P J,YAO ZH J. Improvement of the measuring distance of repetitive-frequency laser spot in field[J]. Chinese Journal of Liquid Crystals and Displays,2016,31(12):1137-1142.(in Chinese)

[35] 王嘉成,孙海江,刘培勋,等.高分辨率多传感器融合图像跟踪系统的设计与实现[J].液晶与显示,2016,31(8):825-830.

WANG J CH,SUN H J,LIU P X,et al.. Design and implementation of high resolution multi-sensors fusion tracking system[J]. Chinese Journal of Liquid Crystals and

Displays,2016,31(8):825-830.(in Chinese)

[36] 丁鹏,张叶,刘让,等.结合形态学和Canny算法的红外弱小目标检测[J].液晶与显示,2016,31(8):793-800.

DING P,ZHANG Y,LIU R,et al.. Infrared small target detection based on adaptive Canny algorithm and morphology[J]. Chinese Journal of Liquid Crystals and Displays,2016,31(8):793-800.(in Chinese)

[37] 王玉坤,贾娜,张锐.激光通信成像光斑处理方法研究[J].液晶与显示,2017,32(9):736-740. WANG Y K,JIA N,ZHANG R. Laser communication spots imaging process method[J]. Chinese Journal of Liquid Crystals and Displays,2017,32(9):736-740.(in Chinese) [38] 吴春婷,李贺,蔡继兴,等.采用马赫-曾德尔干涉法测量单晶硅在线应力损伤[J].光学 精密工程,2017,25(5):1395-1401.

WU CH T,LI H,CAI J X,et al.. Measurement of real-time stress damage of monocrystal silicon by Mach-Zehnder interferometry[J]. Opt. Precision Eng.,2017,25(5):1395-1401.(in Chinese)

[39] 李颖奎,齐冀,张洁,等.基于液晶空间光调制器的变倍率激光扩束技术研究[J].液晶与显示,2018,33(9):764-771.

LI Y K,QI J,ZHANG J,et al.. Laser beam expanding technology based on liquid crystal spatial light modulator[J]. Chinese Journal of Liquid Crystals and Display,2018,33(9):764-771.(in Chinese)

[40] 孔晓波,刘丽娟,刘永刚,等.基于液晶/聚合物光栅的可调谐双波长有机激光器[J].液晶与显示,2018,33(1):49-54.

KONG X B,LIU L J,LIU Y G,et al.. Tunable dual-wavelength organic laser based on

holographic polymer dispersed liquid crystal grating[J]. Chinese Journal of Liquid Crystals and Display,2018,33(1):49-54.(in Chinese)

[41] 乌日娜,王兴,杨帆,等.染料掺杂手性向列相液晶器件中实现随机激光辐射[J].液晶与显示,2018,33(6):464-468.

WU R N,WANG X,YANG F,et al.. Realization of random lasing in dye-doped chiralnematic liquid crystal[J]. Chinese Journal of Liquid Crystals and Display,2018,33(6):464-468.(in Chinese)

[42] 吴迪,葛廷武,秦文斌,等.光纤激光器故障模式分析[J].发光学报,2018,39(7):1002-1007.. WU D,GE T W,QIN W B,et al.. Fault mode analysis on fiber laser[J]. Chinese Journal of Luminesecence,2018,39(7):1002-1007.(in Chinese)

[43] 顾宏灿,黄俊斌,程玲,等.20～1 250 Hz光纤激光加速度传感系统设计[J].中国光学,2017,10(4):469-476.

GU H C,HUANG J B,CHENG L,et al.. 20-1 250 Hz fiber laser acceleration sensing system[J]. Chinese Optics,2017,10(4):469-476.(in Chinese)

[44] 谢检来,郝永芹,张家斌,等.一种具有低折射率的高对比度光栅反射镜的设计[J].发光学

报,2018,39(6):855-861.

XIE J L,HAO Y Q,ZHANG J B,et al.. Design of high contrast grating mirror with low index grating layer[J]. Chinese Journal of Luminesecence,2018,39(6):855-861.(in Chinese) [45] 张磊,刘东,师途,等.光学自由曲面面形检测技术[J].中国光学,2017,10(3):283-299. ZHANG L,LIU D,SHI T,et al.. Optical free-form surfaces testing technologies[J]. Chinese Optics,2017,10(3):283-299.(in Chinese)

[46] 吕强,李文昊,巴音贺希格,等.基于衍射光栅的干涉式精密位移测量系统[J].中国光学,2017,10(1):39-50.

LV Q,LI W H,BAYINHEXIGE,et al.. Interferometric precision displacement measurement system based on diffraction grating[J]. Chinese Optics,2017,10(1):39-50.(in Chinese)