Design of an endoscopic OCT probe based on piezoelectric tube with quartered outside electrodes

1 Introduction

Optical Coherence Tomography (OCT) is a non-invasive high-resolution tomography method that utilizes the low-coherence interference principle of light to image inside bio-logical tissue to be detected (Huang et al., 1991). The combination of the endoscopy head and OCT has a wide range of applications, such as fundus retinal imaging (Joos and Shen, 2013), vascular endoscopy imaging (Yu et al., 2017), gastrointestinal endoscopy imaging (Masci et al., 2009), and tooth root crack detection (Wei et al., 2019) and other clinical applications, which can realize non-invasive living tissue high-resolution tomography imaging of human internal organs.

Ophthalmic endoscopy OCT is a special kind of OCT that uses the endoscopic OCT probe mechanism to bypass some opaque medium and directly perform structural imaging. It can display the cross-sectional outline of the large transverse field of view on the retina while reducing the influence of cloudy fluid on the fundus tissue structure imaging so that the surgeon can obtain enough information during surgery to complete intraoperative decision-making. For example, the observation of separation during retinal detachment surgery. At the same time, with the rapid development of robots today, ophthalmic surgical robots have been able to complete clinical surgical trials such as retinal stem cell injection. Endoscopic OCT will effectively solve the problem that the field of vision is not enough during robot surgery, which can ensure the safety of clinical surgical trials.

At present, the fundus OCT imaging system usually uses external irradiation. The light beam is used to image the macula, optic nerve, retinal nerve fiber, and choroid of the fundus (posterior pole of the eye) through the optical system of the cornea, pupil, and crystalline lens (Joos and Shen, 2013), and the field of view is severely limited. The low resolution of images obtained by external OCT may also occur due to optical aberration of media in the eye (Kang et al., 2010). In addition, distortion or turbidity of subjects’ eye media, such as dense cataracts and other factors, will reduce externally transmitted OCT signals and reduce the signal-to-noise ratio of the fundus image (Kang et al., 2010). In addition, endoscopy OCT can display the cross-sectional outline of the upper retinal membrane, avoiding the forced use of potentially toxic dyes to visualize the membrane due to insufficient technology (Almony et al., 2012). The use of ophthalmic endoscopy OCT imaging can well solve the three problems mentioned above in the fundus OCT imaging system of external irradiation. At the same time, optical coherence tomography (OCT), a non-invasive high-resolution non-destructive tomography technology, is introduced into ophthalmic surgery to improve the success rate of ophthalmic diseases and avoid using contrast media. Reduce surgical complications.

In some robots capable of performing ophthalmic microsurgery, OCT integrates surgical instruments as an important auxiliary part of the surgical robot, which realizes cooperation with the operator to accurately achieve synchronous and relative movement. The probe enters the eye through the scleral opening with a maximum acceptable diameter of 3 mm (Hayashi et al., 2014) for endoscopy imaging, which is more conducive to real-time observation of surgery and surgical conditions, to guide tissues to avoid Angle, lens and other important ocular tissue structures, avoid secondary damage (Yu et al., 2015), and improve the safety and controllability of surgical robot operation.

Using the probe to enter the eyeball can properly solve the problem that the opaque tissue and turbidity media in the eye limit the imaging of the fundus retina (Bouma et al., 2000) and other structures. The external handheld probe has a very compact structure, which can be entered into the eyeball to perform OCT imaging on the fundus tissue structure after minimally invasive treatment of the eyeball, and the problem of mutual interference with the volume and position of optical devices such as optical microscopes is greatly optimized.

In 2005, N. V. Iftimia et al. (Iftimia et al., 2005) developed a handheld OCT A-scan imaging system with a central wavelength of 1310 nm to measure depth-resolved information in solid tissue. The designed probe diameter was 250 μm, but this design could not provide two-dimensional information unless the entire probe was moved to produce a scan.

In 2013, Joos et al. (Joos and Shen, 2013) proposed a smaller, independent 25G handheld forward imaging B-scanning OCT probe. The scanning probe uses a coil magnetic oscillator to drive the fiber for scanning, and the mechanism to generate the fiber tip scanning is formed by sliding a 28G straight thin-walled stainless-steel tube along the curved part of the 34G probe tube. When the 28G probe slides along the 34G probe, it forces the tip of the 34G probe to drive the 125 μm single-mode fiber to vibrate in the aerated space for lateral scanning to achieve endoscopic OCT imaging.

In 2016, Asami et al. (2016) proposed a 23G intraoperative handheld OCT probe with optical fiber scanning driven by motors. With this handheld OCT probe, researchers successfully achieved in vitro OCT imaging of pig eye and rabbit eye retinal tissues. Subsequently, it was successfully applied to human clinical cases to distinguish the retina, optic disc, choroid, and other fine intraocular structures through OCT imaging.

In 2016, Mura et al. (2016) performed an OCT endoscopy with a 23G side-scanning SD-OCT probe, confirming that the use of intraocular SD-OCT can expand visual cues during surgery, help the decision-making process, and improve the outcome of surgery. Intraoperative OCT images of 7 patients with vitreoretinal disease were obtained by 20G probes.

Endoscopic OCT probes can be divided into two types according to different scanning modes: side-imaging (Mura et al., 2016) and forward-imaging (Joos and Shen, 2013). Side-imaging means that the scan is emitted and received by the side of the probe, and forward-imaging means that the scan light is emitted and received by the front end of the probe. From the perspective of safety, applicability, image quality, and entrance incision, the forward mode is more suitable for ophthalmic endoscopic OCT imaging. At present, the driving methods of endoscopic OCT probes mainly include Micro-Electro-Mechanical Systems (MEMS) (Zhang et al., 2015), motors (Asami et al., 2016), electromagnetic driving (Joos and Shen, 2013), piezoelectric material driving, etc. MEMS and motors are generally suitable for lateral scanning, and their size limits the size reduction of the endoscope needle.

Due to the selection of driving mode, the current endoscopic OCT probes used for ophthalmic imaging are constrained by the reduction of probe size, and the scanning mode is constrained by lateral scanning, which is not conducive to the safety and improvement of imaging quality of ophthalmic endoscopy. In this paper, the driving mode of the four-part piezoelectric ceramic tube is selected, and the input signal control is more convenient. The dimensions of the PZT tube used in this study were 35 mm in length, 5.7 mm in outer diameter, and 0.67 mm in wall thickness. The size of the PZT tube does not limit the size reduction of the probe, which greatly reduces the fabrication cost of the probe. The biggest difference with our endoscopic OCT probe is that we can change the cannula to a minimum size of 25G (outer diameter 0.51 mm). In theory, it is possible to make a smaller probe with a suitable cannula. In fact, due to the characteristics of the structural design, our endoscopic OCT probe can be replaced with different sleeves to achieve any size from 13G (external diameter 2.41 mm) to 25G (external diameter 0.51 mm) to meet different visual field imaging requirements and different safety requirements.

In this study, a micro intraocular OCT probe was fabricated using a piezoelectric tube with quartered outside electrodes as a driver, which could realize Lissajous scanning of the movement track at the end of a single mode fiber and enter the eyeball through the scleral opening to perform real-time tomography of the fundus retina. Due to the particularity of the structural design, the outer diameter of the endoscopic OCT probe in this study can be replaced by 13G (outer diameter 2.41 mm)- 25G (outer diameter 0.51 mm), to achieve different imaging field sizes and ensure safe access to the ocular scleral incision. In this paper, the principle, theoretical derivation, and feasibility analysis of the application of a piezoelectric tube with quartered outside electrodes on the snooping probe in OCT are introduced, as well as the design characteristics, practical application and imaging situation of the endoscopic OCT probe are introduced, and the future direction of further research has prospected.

2 Principle derivation and feasibility analysis

2.1 Introduction to the principle of piezoelectric ceramics

Piezoelectric ceramics are a class of crystalline materials (Li et al., 1999). The compression or elongation caused by mechanical pressure will cause the two ends of the piezoelectric ceramics to produce different charges, resulting in voltage differences (Ravez and Simon, 2000). As a reversible, adding different voltages to the two ends of the piezoelectric ceramic produces a voltage difference between the two ends, which will also result in mechanical displacement or stress (Fatikow, 1996).

The piezoelectric tube with quartered outside electrodes is a hollow cylindrical piezoelectric ceramic tube. A conductive coating is coated on the inner cylinder surface of the piezoelectric ceramic tube to cover the electrode inside the piezoelectric ceramic tube as a common electrode. In the outer layer, the circumference of the piezoelectric ceramic tube is evenly divided into four quadrants, each of which is coated with electrodes, the area is equal, but the electrodes are insulated from each other, and the electrodes are independent of each other. As shown in Figure 1, when the voltage signal with the same amplitude and opposite symbol is added to the -X and +X poles, the piezoelectric ceramic can be deflected in the X direction. If the sinusoidal voltage signal of a certain frequency is applied, the voltage ceramic tube can be repeatedly deflected by a certain frequency to achieve vibration in the X direction. Similarly, a voltage signal of a certain frequency with the same amplitude and opposite symbol is applied to the poles of -Y and +Y in the Y direction, and the voltage ceramic tube can achieve repeated deflection according to a certain frequency to achieve vibration in the Y direction. When two voltage signals of a certain frequency are applied in the X and Y directions at the same time, the two-dimensional vibration deflection in the X and Y directions of the piezoelectric tube is realized. To use the piezoelectric tube with quartered outside electrodes correctly and rationally, the dynamics analysis is carried out below (Qin et al., 2007).

Figure 1. Piezoelectric tube with quartered outside electrodes.

2.2 Kinetic analysis and derivation of piezoelectric ceramic tube

The following is the theoretical analysis and calculation of the deflection of the piezoelectric tube with quartered outside electrodes. As shown in Figure 2A, the extension of the z-axis can be expressed as Eq. 1 (Chen, 1992)

Figure 2. (A) Geometric relationship before and after X and Y-axis deformation; (B) Force relation after deformation of the Y-axis.

Where $d_{31}$ is the axial strain coefficient of the piezoelectric ceramic tube; $L$ is the length of the piezoelectric ceramic tube; $U$ is the voltage between electrodes outside the piezoelectric ceramic tube; $t$ is the wall thickness of the ceramic pipe.

When the piezoelectric ceramics are extended axially under the action of the electric field, the wall thickness will also change, and the wall thickness change can be expressed as Eq. 2

Where $d_{33}$ is the strain coefficient of the piezoelectric ceramic tube in the polarization direction.

The z-axis strain can be expressed as Eq. 3

For the X and Y directions, when a voltage with the same amplitude and opposite sign is applied to the opposite two electrodes, the part applying a positive electric field will extend $\Delta S_1$, while the part responsible for applying a negative electric field will shorten $\Delta S_2$, which is expressed as Eq. 4

when the positive and negative electric fields are equal.

As shown in Figure 2A, the displacement $△y$ in the Y-axis direction is represented as Eq. 5

When the Y-axis is moving, the coupling caused in the Z-axis is expressed as Eq. 6

Let $\Delta S^0$ be the displacement generated when the electric field in the same electric field strength is applied, which can be expressed as

As shown in Figure 2B, the differential element $dl$ was taken for analysis, and the deformation produced by the differential element $dl$ was $\pm dS$. The external forces generated by the surface of the PZT tube were $\pm F$ when the PZT tube was operated under driving voltages of the same magnitude and opposite sign. This pair of forces of equal magnitude and opposite directions form a couple whose equivalent bending moment is let be $\Delta M$. The PZT tube undergoes bending deformation under the action of this couple, and $d\theta$ is the end Angle. According to the mechanics of materials, Eqs 8, 9 is established:

Where D is the outer diameter of the PZT tube; d is the inner diameter of the PZT tube; $D_0$ is the distance between the centroid of the two deformed members in the Y-axis direction, $D_0=\left(D+d\right)/2$;

According to the relationship between stress and strain in material mechanics, Eq. 10 can be obtained:

Where $E$ is the elastic modulus of the PZT tube; $I$ is the Second moment of area of the PZT tube, $I=\left[\pi \left(D^4-d^4\right)\right]/64$; $\mathrm{A}$ is the area under force, $A=\left[\pi \left(D^2-d^2\right)\right]/16$; $L$ is the length of the PZT tube; $\sigma$ represents the stress of the material in pascals (Pa) and represents the ratio of the force applied inside the material to its unit area; $\epsilon$ represents the strain of the material, no unit, and represents the ratio of the change in the length of the material to the original length.

Eqs 11–13 can be obtained by static analysis:

where $D_0$ is the distance between the centroids of the two deformed members in the Y-axis; $D$ is the outer diameter of the piezoelectric ceramic tube; $d$ is the inner diameter of the piezoelectric ceramic tube; $R$ is the bending radius of the ceramic tube; $\theta$ is shown in Figure 2A.

The end of the piezoelectric ceramic tube should be connected to a single-mode fiber cantilever beam. If the length of the single-mode fiber cantilever beam is $l_t$, the transverse displacement is:

Bring Formulas 12–14 into (Eq. 15) to get

In theory, $\Delta \mathrm{x}$ and $\Delta \mathrm{y}$ are not coupled, so it can be considered that they have the same displacement equation. By substituting Eq. 7 into Eq. 16, displacement equations in X and Y directions can be obtained. The displacement equations of the end of the piezoelectric ceramic tube in three directions Eqs 17–19 are obtained by considering the coupling caused by the lateral displacement in the Z direction.

Eq. 19 is approximately the displacement in the Z-axis direction of the PZT tube without deflection, which can be approximately regarded as $\Delta {\mathrm{z}}^{\prime }$ is proportional to $U_Z$.

Through the dynamic analysis of the whole, when the PZT tube drives the end of the single-mode fiber to scan, the certain displacement of the end of the single-mode fiber in the x, y, and z directions is proportional to the $\Delta {\mathrm{x}}^{\prime }$, $\Delta {\mathrm{y}}^{\prime }$, $\Delta {\mathrm{z}}^{\prime }$, respectively. Therefore, the vibration of the piezoelectric tube with quartered outside electrodes and the end of the single-mode fiber can be improved by adjusting the signal voltage.

2.3 Analysis and derivation of vibration model of single-mode fiber cantilever beam

A cantilever beam is a structure with one end fixed and another suspended. Compared with the fixed end that does not generate axial or vertical displacement, the suspended free end can generate both axial and vertical forces. The simplified model is shown in Figure 3. The fiber cantilever beam is a fixed fiber end, and the other end of the suspended fiber makes it a free end. In the case that the single-mode fiber is not affected by gravity and the external environment, the fiber cantilever is in a horizontal state (that is, reconnected with the X-axis). When there is no vibration in the external environment, the fiber is only affected by gravity, currently, the external load is 0, and the load of the cantilever beam is only the gravity of the cantilever beam itself, currently, the cantilever beam is in a free state. As shown in Figure 3, the coordinate system shown above is established. The X-axis overlaps with the cantilever beam when the cantilever beam is stationary and not affected by gravity. O is the origin of the coordinate system, the total length of the cantilever beam is l, the distance from any point of the cantilever beam to the fixed end is x, the uniform load is q, and the arbitrary angle of the fiber optic cantilever beam is θ.

Figure 3. Cantilever beam structure diagram.

The mechanical model of the Euler-Bernoulli beam is used in the theoretical analysis of the natural frequency of a fiber optic cantilever beam. The vibration of the cantilever beam belongs to the vibration of the continuous elastomer in material mechanics, so the vibration of the cantilever beam has infinite degrees of freedom and corresponding natural frequency and principal mode. The final vibration of the cantilever beam can be expressed as the superposition of infinite principal mode. When the cantilever beam is forced to vibrate due to the external environment, only the deformation caused by bending is considered, and the influence of the deformation caused by shear and the moment of inertia is ignored, which conforms to the mechanical model analysis of Euler-Bernoulli beam (Necla and Süleyman, 2016).

As can be seen from the structural diagram of the cantilever beam in Figure 3, one end of the cantilever beam is a fixed end, and the other end is free. The single-mode fiber is recombined with the X-axis without gravity, and the differential equation of the cantilever beam motion (Eq. 20) (Liu, 2003; Ohm et al., 2006; Sun et al., 2009) is:

Where $E$ is the elastic modulus of single-mode fiber, it can be seen from the data that $E=7.4668\times {10}^{10}Pa$; $I$ is the Second moment of area of the optical fiber, $I=\left(\pi D^4\right)/64$; $D$ is the diameter of the optical fiber 125 μm; $\mathrm{A}$ is the optical fiber cross-section area; $\rho$ is the density of conventional communication single-mode fiber, $\rho =2203kg/m^3$; $x$ is the distance of a beam section from the origin O point; $w\left(x,t\right)$ is the displacement of the beam section at time $t$ at the distance $x$ from the origin O point; $t$ there stands for a moment in time.

The differential equation has the fourth partial derivative with respect to $x$ and the second partial derivative with respect to $t$, so it requires four boundary conditions and two initial conditions to solve. The boundary conditions of the cantilever beam (Fang and Zhang, 2011; Yang et al., 2011) are denoted as Eqs 21–24.

The vibration mode of the system is independent of time, so the equation can be solved by the method of separating variables, and the solution of the equation can be separated by variables, and the free vibration solution of the partial differential equation can be obtained as Eq. 25:

By substituting this solution into the differential equation of motion of a cantilever beam, Eq. 26 is obtained:

where ${\beta }^4=\frac{\rho A{w}^2}{\text{EI}}$.

Bring the boundary condition Eqs 21, 22 into the above formula to get Eq. 27:

Further sorting can obtain Eq. 28:

Then the boundary conditions Eqs 23, 24 are brought into the equation:

So, the resulting frequency equation is Eq. 30:

the root ${\beta }_{n}l$ of this equation represents the natural frequency of the vibration system.

The value of ${\beta }_{n}l\left(n=1,2,3\dots \right)$ satisfying the above Eq. 31 is Eq. 32:

If the value of $C_2$ relative to ${\beta }_n$ is expressed as $C_{2n}$, according to $C_{1n}$ in the formula, $C_{2n}$ can be expressed as Eq. 33:

The principal mode function corresponding to the natural frequency of order $n$ (Eq. 34) can be obtained:

Thus, the first five natural frequencies of the cantilever beam (Eq. 35) can be obtained, and n = 1,2,3,4,5 can be brought into the equation:

Where $w=2\pi f$, $f$ is the natural frequency of vibration of the cantilever beam.

As shown in Figure 4, the natural frequencies of fiber cantilever beams with different lengths can be obtained by MATLAB simulation. Since the resonant frequency of the cantilever beam is mainly determined by the first-order natural vibration frequency, the first-order vibration frequency of the cantilever beam with different lengths is shown in Table 1.

Figure 4. Relationship between the length of the single-mode fiber cantilever beam and the first-order natural frequency.

Table 1. First-order vibration frequencies of single-mode fiber cantilever beams of different lengths.

The single-mode fiber cantilever beam is regarded as a stressed object with a uniform load. As shown in Figure 3, the fiber is rejoined with the X-axis without gravity. In the range of $0\le x\le l$, the bending moment equation of the single-mode fiber cantilever beam is Eq. 36:

Where $M\left(x\right)$ is the natural offset of the fiber cantilever of different lengths.

The differential equation and Angle equation of the deflection curve Eqs 37–39 are as follows:

Integrate Formula 39 to obtain:

It can be seen from the boundary conditions that when $x=0$, $y=0$, $\theta =0$. Therefore, the integral constants C = 0 and D = 0 in (Eq. 40) can be determined, that is, the Angle and deflection of the fiber cantilever beam at the fixed point are both 0. Eq. 40 can be simplified as Eq. 41:

The equation that can finally obtain the offset is Eq. 42:

When $x=l$, the offset $y$ (Wang et al., 2013; Zhang et al., 2013; Xu et al., 2014) has a maximum value, which is Eq. 43:

When the external force is not applied to the cantilever beam, the load is only the gravity of the fiber cantilever beam itself. $q$ is the uniform load, that is, the gravity per meter of the cantilever beam, so the above equation can be rewritten as (Beléndez et al., 2008):

Where W is the gravity of the cantilever beam itself, $\mathrm{W}=ql$.

Bring the Formula 45 into (Eq. 44) to simplify:

Where $y_{\max }$ is the maximum offset of the cantilever beam; $g$ is the acceleration of gravity, $g=9.8m/s^2$; $l$ is the length of the fiber optic cantilever; $I$ is the Second moment of area, $I=\left(\pi D^4\right)/64$.

The minus sign in Eq. 46 indicates only the direction. It can be seen from Eq. 46 that as the length of the optical fiber cantilever increases, the offset of the optical fiber cantilever increases. However, when the length of the fiber cantilever beam is very small, the vibration amplitude of the fiber cantilever beam is too small when measuring the vibration signal, the light intensity received by the receiving end will hardly change, and the imaging effect will be greatly reduced, so it is necessary to find the appropriate length of the fiber cantilever beam at the transmitting end. Considering all aspects, the length of single-mode fiber cantilever beam is 22–23 mm.

3 Study and design of endoscopic OCT probe

3.1 Mechanical structure design of endoscope OCT probe

The resonance frequency and the maximum amplitude of the optical fiber cantilever beam are analyzed theoretically. The probe was designed to fix the treated single-mode optical fiber in the center position of the piezoelectric tube with quartered outside electrodes by designing 3D printed units and fixed with glue using precision instruments.

To reduce the cost of the piezoelectric tube with quartered outside electrodes, the probe structure was innovatively optimized. The probe no longer installs the piezoelectric ceramic tube in the probe tube wall, thus removing the limitation of the probe size on the piezoelectric tube with quartered outside electrodes, and there is no need to customize the piezoelectric tube with quartered outside electrodes of specific size, thus greatly reducing the manufacturing cost of the probe. Through the theoretical simulation analysis of the single-mode fiber cantilever beam, the fiber extension length is set at about 22–23 mm, so that the extended single-mode fiber is completely under the protection of the medical stainless-steel probe.

Through Gaussian beam modeling and simulation, the GRIN lens with appropriate parameters is selected to make the divergent beam emitted from the end of the single-mode fiber get a good focus, and at the same time, the optimal distance between the end of the single-mode fiber and the GRIN lens and the optimal imaging distance from the GRIN lens to the imaging object is determined, so as to effectively improve the imaging quality. The GRIN lens is fixed at the front end of the medical stainless-steel probe tube, and the GRIN lens maintains the simulated optimal distance from the end of the single-mode fiber, so that the scanning beam of the single-mode fiber is focused through the GRIN lens and resonates with the single-mode fiber at a certain frequency and a certain trajectory to achieve high-quality imaging, as shown in Figure 5.

Figure 5. (A) Structure diagram of the endoscope OCT probe designed in this study; (B) A schematic depiction of the completed encapsulation process for the endoscopic OCT probe designed in this study; (C) Mechanical structure drawing of the physical object; (D) Some units of the endoscope OCT probe are 3D printed.

3.2 COMSOL modeling was performed for resonance analysis

The equations should be inserted in editable format from the equation editor.

After the preliminary mechanical structure design of the endoscope OCT probe was completed, the modeling analysis was performed using COMSOL, considering that the actual situation was different from the theoretical analysis of the single-mode fiber cantilever beam analyzed previously. Because the single-mode fiber is processed using precision instruments, it is fixed with glue at the center of the piezoelectric tube with quartered outside electrodes, rather than a bare fiber.

The piezoelectric tube with quartered outside electrodes, single-mode optical fibers, some 3D printed units, and adhesions for precision instruments (used in connection sites) were modeled using COMSOL. The density, Young’s modulus, Poisson’s ratio, and other material parameters of each part were respectively input to conduct full-true modeling. The characteristic frequency analysis of the built model was carried out to obtain the resonance frequency at this time.

As shown in Figure 6A, when the frequency is 162.91 Hz and the length of the single-mode fiber cantilever beam is 23.5 mm, COMSOL modeling analysis shows that the fiber end has an obvious vibration amplitude of 2–3 mm in the y-z plane. As shown in Figure 6B, when the frequency is 162.9 Hz and the length of the single-mode fiber cantilever beam is 23.5 mm, the end of the fiber appears obvious vibration in the x-z plane, and the vibration amplitude is 2–3 mm, which meets the predetermined requirements. As shown in Figure 6C, after applying the voltage signal, the obvious vibration at the end of the single-mode fiber can be directly observed. As shown in Figure 6D, the maximum vibration amplitude of single-mode fiber can reach 2360 μm after accurate observation and measurement under the electron microscope by applying the above voltage signal.

Figure 6. (A) In COMSOL modeling, the resonant frequency is 162.91 Hz, and the vibration amplitude of 2–3 mm appears at the end in the y-z plane; (B) The resonant frequency in COMSOL modeling is 162.9 Hz, and the vibration amplitude of 2–3 mm appears at the end in the x-z plane; (C) The vibration of the end of the single-mode fiber is clearly and directly observed with the bare eye; (D) The maximum vibration amplitude of the single-mode fiber can reach 2360 μm after obvious observation and measurement under the electron microscope.

3.3 Study on Lissajous scanning trajectory at the end of single-mode fiber

The Lissajous trajectory was selected for the single-mode optical fiber end scanning trajectory of the OCT probe, which is driven by two orthogonal harmonic vibrations, and determined by the amplitude, frequency, and initial phase of the two harmonic vibrations. The resonant frequency of 162.9 Hz obtained by the previous resonance simulation using COMSOL is taken as a reference, and the simple harmonic vibration frequencies in the x and y directions are set to 163 Hz and 162 Hz respectively. Considering that the output voltage of the signal will be amplified to 60 V in the actual situation, the amplitudes of the two simple harmonic vibrations are set to 60 V. As shown in Figure 7, the input signal in the x direction will be simulated by MATLAB for scanning trajectory, and the adjustment of parameters such as density, filling rate, and lobe number of Lissajous trajectory can be realized by phase adjustment, to achieve the optimal single-mode optical fiber end scanning.

Figure 7. Because of the different phase difference of x and y direction signals of the Lissajous locus, the density of the Lissajous locus and the number of lobes are obviously different. (A) The phase difference between x and y signals is π/2; (B) the phase difference between the x and y signals is π/3.

Considering the actual scanning situation of the endoscopic OCT probe, the end of the single-mode fiber is a focused spot with a certain actual size. By connecting the single-mode fiber to the laser and using the COMS camera for observation, the trajectory of the focused spot at the end of the fiber can be observed. As shown in Figure 8, we show the contrast between the end of a single-mode fiber with a focused spot of a certain size captured by a real CMOS camera and the simulation in MATLAB. The scanning range of the end of the single-mode fiber was 2.1 mm × 2.1 mm at this time of actual measurement.

Figure 8. (A,B) are the movement trajectories of the end of the single-mode fiber taken with the COMS camera, and (C,D) are the simulation plots of the Lissajous trajectory carried out in MATLAB. It can be found that (A) achieves good contrast with (C,B) with (D), respectively.

3.4 Design and construction of optical path for endoscopic OCT probe

The imaging principle of endoscopic OCT probe technology is Optical Coherence Tomography (OCT), which is based on the principle of low coherence interference of light. A reference arm is set as the reference object of vertical and horizontal depth, and the sample arm is moved for scanning. The backward-facing interference light returned from the sample at different depths is processed, and the fault structure of the tissue to be measured is obtained after processing, to achieve A-scanning, that is, a series of one-dimensional longitudinal depth information is obtained (Danielson and Boisrobert, 1991). A series of A-scans are obtained by controlling the sample arm scanning on the transverse plane, and two-dimensional images are obtained through processing and analysis, that is, B-scanning is realized. Also, a series of B-scans obtained by scanning were processed and analyzed to obtain three-dimensional images of the tested samples (Song et al., 2021). As shown in Figure 9, the optical path was designed and built, and the length of each fiber was accurately calculated to ensure the optical path matching to ensure normal imaging.

Figure 9. (A) Design of the optical path structure diagram. (B) The actual optical path structure.

4 Encapsulation technology and imaging experiment of endoscopic OCT probe

4.1 Encapsulation technology experiment of endoscopic OCT probe

For the encapsulation of the endoscopic OCT probe, it is necessary to ensure the accuracy of the encapsulation process. As shown in Figure 10, it is ensured that the scanning optical fiber is in the geometric center of the medical stainless-steel probe, so that the resonant vibration at the end of the fiber is not affected by the medical stainless-steel probe. It is necessary to design some assembly and encapsulation devices to assist the encapsulation work. At the same time, the GRIN lens is installed for focusing, and the distance between the GRIN lens and the end of the single-mode fiber is required to ensure that the best effect of Gaussian light speed simulation is achieved. This adjustment of accuracy requires the installation of a precision displacement platform and the use of an optical power meter to find the strongest optical power to ensure the best imaging quality.

Figure 10. (A) and (B) are single-mode optical fibers located in the geometric center of the medical stainless-steel probe to ensure good imaging quality.

The preliminary encapsulation of the endoscopic OCT probe was obtained, and the finished product was shown in Figure 11.

Figure 11. The encapsulated endoscope OCT probe can be connected to the OCT optical path and apply input signals for OCT imaging.

4.2 Imaging experiment of endoscopic OCT probe

The encapsulated endoscopic OCT probe was first used in the finger for imaging experiments. As shown in Figure 12, OCT tomographic images with good effect and obvious stratification of finger tissue can be obtained. The endoscopic OCT probe has not entered the eyeball in vivo for fundus retinal experiments in this study, which will be carried out in the future.

Figure 12. When observing a finger in vivo using an endoscopic OCT probe, OCT tomographic images with a distinct finger tissue stratification effect can be observed, where the scale bar is 1 mm.

The endoscopic OCT probe in this study can be used to adjust the vibration amplitude of the end of the single-mode fiber by controlling the applied signal voltage, to adapt to different types of 13G (outer diameter 2.41 mm)-25G (outer diameter 0.51 mm) endoscopic probes, which can be seen in Figure 13. The OCT probe size can be changed from 13G to 25G according to the work needs, which will bring different imaging field sizes. It is worth noting that GRIN lenses need to be customized to match the probe size to achieve normal beam focusing, and the encapsulation of endoscope OCT probes needs to ensure a reasonable optical distance.

Figure 13. (A) 13G (outer diameter 2.41 mm) endoscopic OCT probe made in this study; (B) 25G (outer diameter 0.51 mm) endoscopic OCT probe after probe replacement.

5 Prospects for the subsequent development of endoscopic OCT probes

It is worth looking forward to the possibility that the endoscopic OCT probe can be assembled with the RCM robotic arm for subsequent automated imaging. It can be considered to realize two robotic arms for ex vivo ophthalmic surgery. One robotic arm is equipped with a scalpel or injection needle, and the other robotic arm is equipped with the endoscopic OCT probes from this study so that the two robotic arms can be matched for work, where the position, angle, and other parameters of surgical instruments can be displayed in real-time. This can greatly improve the safety and reliability of the robotic arm in ophthalmic surgery.

How to achieve more stable OCT scanning, how to improve the performance parameters such as the resolution of endoscopic OCT, and how to ensure the high-precision cooperation between the endoscopic OCT probe and the surgical robot will become the technical focus and difficulties in the follow-up research.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement

Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

JH: Conceptualization, Investigation, Methodology, Validation, Visualization, Writing–original draft, Writing–review and editing, Formal Analysis. SW: Validation, Writing–review and editing. GS: Project administration, Supervision, Writing–review and editing. JF: Formal Analysis, Project administration, Writing–review and editing. HY: Formal Analysis, Investigation, Writing–review and editing. SC: Formal Analysis, Writing–review and editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This project was funded by: National Key R&D Program of China, Grant Nos.2021YFF0700503. CAS Project for Young Scientists in Basic Research, Grant Nos. YSBR-067. Jiangsu Science and Technology Plan Program, Grant Nos. BK20220263. National Key R&D Program of China, Grant Nos. 2022YFC2404201.

Acknowledgments

The simulations in this study were conducted using COMSOL Multiphysics^® v. [5.6.0.280], COMSOL AB, Stockholm, Sweden.

The data analysis and numerical simulations in this study were performed using MATLAB R2022, MathWorks Inc., Natick, MA, United States.

The probe structure modeling and simulations in this study were performed using SOLIDWORKS 2022, Dassault Systèmes SolidWorks Corporation, Waltham, MA, United States.

The Gaussian beam simulations in this study were conducted using Virtual Lab Fusion ^® v. [7.6], LightTrans International GmbH, Jena, Germany.

We would like to express our sincere thanks to the above software companies.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References