WO2015193191A1 - Multi-core alignment in optical shape sensing - Google Patents

Multi-core alignment in optical shape sensing Download PDF

Info

Publication number
WO2015193191A1
WO2015193191A1 PCT/EP2015/063176 EP2015063176W WO2015193191A1 WO 2015193191 A1 WO2015193191 A1 WO 2015193191A1 EP 2015063176 W EP2015063176 W EP 2015063176W WO 2015193191 A1 WO2015193191 A1 WO 2015193191A1
Authority
WO
WIPO (PCT)
Prior art keywords
optical
cores
fiber cores
optical fiber
length
Prior art date
Application number
PCT/EP2015/063176
Other languages
French (fr)
Inventor
Jinfeng Huang
Gert Wim 't Hooft
Merel Daniëlle LEISTIKOW
Milan Jan Henri MARELL
Jeroen Jan Lambertus Horikx
Original Assignee
Koninklijke Philips N.V.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips N.V. filed Critical Koninklijke Philips N.V.
Publication of WO2015193191A1 publication Critical patent/WO2015193191A1/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/16Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
    • G01B11/161Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge by interferometric means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02001Interferometers characterised by controlling or generating intrinsic radiation properties
    • G01B9/02002Interferometers characterised by controlling or generating intrinsic radiation properties using two or more frequencies
    • G01B9/02004Interferometers characterised by controlling or generating intrinsic radiation properties using two or more frequencies using frequency scans
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02015Interferometers characterised by the beam path configuration
    • G01B9/02027Two or more interferometric channels or interferometers
    • G01B9/02028Two or more reference or object arms in one interferometer
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/0209Low-coherence interferometers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L1/00Measuring force or stress, in general
    • G01L1/24Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet
    • G01L1/242Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet the material being an optical fibre
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B6/00Light guides; Structural details of arrangements comprising light guides and other optical elements, e.g. couplings
    • G02B6/02Optical fibres with cladding with or without a coating
    • G02B6/02042Multicore optical fibres

Definitions

  • the present invention relates to the field of optical shape sensing (OSS). Especially, the invention provides an OSS system and a method for alignment of multiple optical fiber cores in an OSS device.
  • OSS optical shape sensing
  • an optical shape sensor In optical shape sensing (OSS), an optical shape sensor often has a set of helix outer cores and a center core. When the optical shape sensor is curved, each of the outer cores will deform differently.
  • an optical console system providing laser interferometric interrogation of the optical fiber cores in the optical shape sensor, it is possible to provide a 3D shape reconstruction of the optical shape sensor, and thus allows tracking of an elongated device which the optical shape sensor is mounted to follow, e.g. to help navigating a medical interventional device.
  • OSS relies on different deforms of the individual optical fiber cores
  • alignment between the different fiber cores is a key aspect to improve the accuracy of the shape reconstruction.
  • misalignments between the fiber cores are caused by their different refractive index and the helix structure.
  • the OSS system has an increased tolerance towards physical length and refractive index differences of cores in an optical fiber in a multi-core OSS device, i.e. an improved alignment of multiple cores embedded in an optical fiber is advantageous.
  • the solution also leads to OSS with reduced jitter.
  • the solution is preferably stable over time, and with respect to temperature changes.
  • the invention provides an optical shape sensing system arranged to determine a measure of a three-dimensional shape of an associated optical shape sensor comprising a plurality of optical fiber cores arranged between a proximal end and a distal end, wherein the plurality of optical fiber cores comprise optical shape sensing properties, wherein a center core is one of the plurality of optical fiber cores arranged centrally to a plurality of outer cores, wherein the system comprises an optical console system arranged for optical interrogation of the optical fiber cores in the associated optical shape sensor, wherein the optical console system is arranged:
  • Such OSS device with calibration steps a)-c) is advantageous, since it is possible to align physical length as well as refractive differences between multiple optical fiber cores of an optical shape sensor. This improves OSS accuracy, and the calibration measurements only introduce little extra calculation time which can easily satisfy the real- time requirements of an OSS system.
  • the optical console system may be arranged to determine peak positions of start and end points for all of the optical fiber cores in response to a discrete spectral representation, e.g. in the form of a Fast Fourier Transform (FFT), of the optical calibration measurements, and to calculate the individual measures of optical lengths as distances between start and end points in response to the determined peak positions.
  • FFT Fast Fourier Transform
  • the individual optical scan wavelength ranges for each of the outer cores may be determined so as to equal a recognized number of locations along the fiber core of each of the plurality of optical fiber cores, since the determined optical lengths for each of the plurality of optical fiber cores may be considered as an expression of a finite number of single locations along the fiber core which can be optically recognized. This can be equalized between the fiber cores, if their individual optical scan wavelength ranges are adjusted so as to align the optical lengths of the fiber cores.
  • the optical console system is arranged to calculate a measure of optical length ratio for each of the outer cores.
  • the measure of optical length ratio for each of the outer cores is calculated as a ratio between the calculated measure of optical length of the outer core and the calculated measure of optical length of the center core, and to determine the individual optical scan wavelength ranges for each of the outer cores in response to the calculated measure of optical length ratios.
  • This optical length ratio is found to provide a robust measure suited for the aligning, i.e. by determining the individual optical scan wavelength ranges.
  • At least some of a), b), and c) may be implemented in hardware and/or in software. Especially, some of a), b), and c) may be implemented in software, while other part of a), b), and c) are implemented in hardware. Especially, all of a), b), and c) may be implemented in software.
  • the optical console system may be arranged to select a limited portion of sampled data points obtained for an outer core, so as to reduce an individual scan wavelength range for the outer core. This is a rather simple software implementable way of aligning the fiber cores.
  • the limited portion of sampled data points may be selected to be in the beginning, in the end, or in the middle of the sampled data points.
  • the optical console system preferably comprises a controllable laser light source system arranged to generate light for the optical interrogation of the plurality of optical fiber cores with a light wavelength sweep within an individually controllable optical scan wavelength range for each of the plurality of optical fiber cores.
  • the optical console system may in a specific embodiment comprise a laser light source system arranged to generate a light wavelength sweep with a wavelength range within the range 1525-1565 nm. It is to be understood that this is by no means limiting, and the skilled person will know how to translate this to other wavelength ranges as well and optical scan wavelength ranges accordingly, and thus the invention is applicable for a wide variety of laser source systems or other light source systems.
  • the system may comprise an optical shape sensor comprising a plurality of optical fiber cores comprising optical shape sensing properties and being arranged between a proximal end and a distal, wherein the plurality of optical fiber cores comprises a center core and a plurality of outer cores, wherein the center core is arranged centrally to the plurality of outer cores being wound around the center core.
  • the optical shape sensor may be arranged to follow a shape of an elongated medical interventional device, and wherein the system is arranged to display an image generated in response to said measure of three- dimensional shape of the optical shape sensor.
  • the associated optical shape sensor preferably has fiber cores with optical shape sensing properties comprising backscattering properties used for optical shape sensing.
  • the optical fibers may comprise at least one of the properties: Rayleigh scattering and fiber Bragg gratings.
  • the multi-core optical fiber may be arranged within one common coating or cover.
  • the optical shape sensor may be mounted to follow the shape of an elongated device, e.g. a medical interventional device, e.g. in the form of a catheter, a guide wire, an endoscope etc.
  • a medical interventional device e.g. in the form of a catheter, a guide wire, an endoscope etc.
  • the invention is applicable also for non-medical use.
  • the number of optical fiber cores in the associated optical shape sensor may be four, i.e. one center core and three outer cores. However, it is to be understood that other number of optical fiber cores can be used, e.g. a total of three: one center core and two outer cores, e.g. for use in robotic applications with a non-twistable robotic arm. In other applications, four or more than four outer cores may be used together with one center core.
  • the amount of data points it is convenient to keep the amount of data points constant over the full length of the optical shape sensor.
  • a power of 2 may be chosen in order to perform fast Fourier transformations (FFT) and make optimal use of parallel computing during the remaining signal processing steps.
  • FFT fast Fourier transformations
  • sensors of different lengths By adjusting the step size in wavelength during data acquisition and consequently adjusting the total wavelength span it is possible to keep the amount of data points the same for different lengths of sensors. For an increase in sensor length, one needs to decrease the wavelength span proportionally. When length variations are large, the variation in wavelength span will also become large. Since strain causes the spectra to shift in wavelength, it means that adjustment of the wavelength span will cause a change in maximum amount of allowable strain, i.e. a change in minimum bending radius.
  • This extension of the invention means that it is possible to exchange length variations of the sensors for variations in minimum bending radius in vice versa.
  • the invention provides a method for calibrating an optical console system arranged for optical interrogation of an optical shape sensor comprising a plurality of optical fiber cores arranged between a proximal end and a distal, wherein the plurality of optical fiber cores comprise optical shape sensing properties, wherein a center core is one of the plurality of optical fiber cores arranged centrally to a plurality of outer cores, the method comprising
  • the invention provides a computer program product having instructions which when executed cause an optical console system to perform the method according to the second aspect.
  • FIG. 1 shows sketch of an OSS system embodiment
  • FIG. 2 illustrates the basic principle of strain measurements used in OSS
  • FIG. 3 illustrates a typical optical shape sensor for OSS with one center core and three outer cores wound around the center core
  • FIG. 4 illustrates the principle behind OSS: optical measurement of compression and tension of various parts of a curved optical shape sensor with multiple optical fiber cores
  • FIG. 5 illustrates an example of measured optical lengths for four optical fiber cores forming part of one optical shape sensor
  • FIG. 6 illustrates an example of sub-index peak detection
  • FIG. 7 shows a graph of measured length of a center core repeated 10000 times
  • FIG. 8 shows a graph of ratio between an outer core and the center core, repeated 10000 times
  • FIG. 9 shows a graph similar to FIG. 8 but carried out 2 weeks earlier
  • FIG. 10 shows a graph of example data points and their corresponding laser wavelength (here an optical scan wavelength range of 20 nm)
  • FIG. 1 1 shows fiber index number versus amplitude for 4 fiber cores with aligned end peaks
  • FIGs. 12a- 12c shows different experimental shapes tested for an optical shape sensor
  • FIGs. 13a shows a photo of an experiment setup
  • FIGs. 13b and 13c show results of optical shape sensor tip position distribution for 100 scans without and with alignment according to the invention
  • FIGs. 14a and 14b show results for end peaks of four cores as fiber index number versus amplitude for the measured data with and without alignment
  • FIG. 15 illustrates steps of a method embodiment.
  • FIG. 1 illustrates basic parts of an OSS system with an optical shape sensor OSF comprising a plurality of optical fiber cores arranged between a proximal end and a distal end where the optical fiber cores form part of the same fiber.
  • the plurality of optical fiber cores comprise optical shape sensing properties, e.g. fiber Bragg gratings, so as to allow OSS of at least a part of their lengths.
  • a center core is one of the plurality of optical fiber cores which is arranged centrally to a plurality of outer cores, e.g. three outer cores, wound around one center core, thereby in itself introducing a physical length difference between the center core and the outer cores.
  • the optical shape sensor OSF incorporated in an elongated device T, e.g. a medical catheter or guide wire.
  • An optical console C has a controllable laser light source capable of providing a light sweep, and an optical detector for detecting an optical response to such light sweep.
  • the optical console C is connected to the optical shape sensor OSF and arranged to optically interrogate the strain sensing optical therein, and to accordingly determine a measure of a three-dimensional shape of at least a part of the optical fiber OSF and thereby the elongated device T in which the optical shape sensor OSF is placed.
  • a processor P controls the optical console C, and a 3D image I of the optical shape sensor OSF can be generated, e.g. displayed as an image on a monitor in real time.
  • the processor P executes a control algorithm CA which serves to perform a calibration procedure with the purpose of aligning differences in physical length and refractive index between the plurality of optical fiber cores that constitute the optical shape sensor OSF, so as to improve OSS accuracy.
  • the calibration procedure comprises performing optical calibration measurements on the plurality of optical fiber cores with one common optical scan wavelength range ⁇ C for all of the optical fiber cores.
  • an initial ⁇ C can be selected to be such as 20 nm, e.g. corresponding to scanning from 1532 nm to 1552 nm.
  • a measure of optical length of each of the optical fiber cores is calculated. E.g. by determining, for each of the fiber cores, fiber index number between detected start and end peaks in the calibration measurement data sample results. Individual optical scan wavelength ranges ⁇ , ⁇ 2, ⁇ 3 for each of the outer cores are then calculated in response to the individual measures of their optical lengths relative to the measure of optical length of the center core.
  • optical scan wavelength ranges ⁇ , ⁇ 2, ⁇ 3 for each of the outer cores can then be applied in optical interrogation of the plurality of optical fiber cores, and to determine a measure of three-dimensional shape of at least a part of the optical shape sensor OSF accordingly.
  • the originally used optical scan wavelength range ⁇ C for the calibration measurements can be used for the center core.
  • an increased OSS accuracy can be obtained, since the introduced individual scan wavelength ranges for the fiber cores serve to align them at a basic signal processing level.
  • Such calibration procedure is easy to implement in existing equipment, e.g. in pure software, and the extra time required for the calibration procedure does not compromise real-time OSS. Experiments show that both the shape accuracy and system jitter can be improved with the proposed calibration procedure.
  • FIG. 2 illustrates the principle behind detecting of position of 4 discrete optical sensors or sensor elements SI, S2, S3, and S4 on an optical fiber OF.
  • a swept laser source LS emits light into a single mode optical fiber OF via a circulator C.
  • the response is detected via the circulator C at an optical detector DT which also receives the light emitted from the laser source LS via a reference arm RF A.
  • the light is reflected at each point SI, S2, S3, and S4 along the fiber.
  • optical interferometry allows measuring the strain-dependent phase shift with high spatial resolution as a function of distance along the fiber OF.
  • the reflected light at each point SI, S2, S3, and S4 interferes with the light from the reference arm RF A.
  • the interfered signal is in the form of a sine wave for each sensor SI, S2, S3, and S4 if the laser scans linearly over frequency.
  • the FFT of the summation signal can clearly identify positions of sensors SI, S2, S3, and S4. This is illustrated to the right where an example of position PS versus amplitude A is shown as four peaks corresponding to the positions of the four sensors S 1 , S2, S3, and S4.
  • FIG. 3 shows a typical example of an optical shape sensor for OSS applications with four single-mode optical fiber cores: one center core CC, and three outer cores CI, C2, C3 helix would around the center core CC.
  • the upper illustration shows a longitudinal section of the optical shape sensor, while the lower illustration shows a cross sectional sketch showing the central position of the center core CC relative to the outer cores CI, C2, C3 within one common jacket JKT.
  • the four cores CC, CI, C2, C3 can end at different locations. However, these length differences are very small. In the worst scenario, an outer core can end a few ⁇ away from that of the center core. E.g., if the distance between the center core CC and outer cores CI, C2, C3 is 50 ⁇ and angle cleave is 8 °, the maximum difference between the outer core CI, C2, C3 and the center core CC is at most 7 ⁇ (8*50* ⁇ /180) in physical space. This small difference can be neglected, so that in the physical space, all four cores CC, CI, C2, C3 can be assumed to start and to end at equal positions.
  • a small angle e.g. 8°
  • FIG. 4 shows an illustration of the optical shape sensor from FIG. 3 in the situation when it is bend, i.e. when the fiber structure is curved, and thus sensors on three outer-cores behave differently. Some will be compressed CMP and some will be stretched due to tension TNS. When the sensors from different cores are properly aligned, the measured phase differences from the aligned sensors can be translated via strains to a full 3- D shape.
  • an optical fiber with optical shape sensing elements can be considered as a series of continuous sensors along its length.
  • the resolution of a sensor ( ⁇ ) is determined by, where, ⁇ is the center wavelength of the scan (e.g. 1542 nm), ⁇ is the optical scan wavelength range (e.g. 20 nm from 1532 nm to 1552 nm), and n is the refractive index of the core. When n is about 1.49, ⁇ is about 40 ⁇ in this example.
  • sensors from each core distribute evenly along its physical length.
  • the sensor resolutions vary from core to core.
  • the physical lengths of the three outer cores are slightly longer than the center one. Therefore, sensors at the same physical location along the sensing fiber from different cores are not naturally matched with each other.
  • FIG. 5 shows an example illustrating this mismatch between the different cores.
  • a spectral representation in the form of FFT signals of the sensors from different cores shows amplitude versus fiber index F_I.
  • F_I fiber index
  • the end of the center core CC and that of an outer core CI have a difference of 13 fiber index numbers, which approximately corresponds 0.5 mm (13*40 ⁇ ) in the physical space.
  • the proposed calibration procedure according to the invention serves to provide this alignment, namely based on an optical calibration measurement of the optical length difference between the cores, thereby resulting in improved OSS accuracy.
  • the refractive index n see e.g. formula (1), where the sensor resolution ⁇ is determined by ⁇ , ⁇ , and n. If each core has a different refractive index, the sensor resolution can vary between cores. Consequently, the number of points can be different for the same physical length.
  • an outer core is slight longer than the center core.
  • an outer core can be about 0.25 mm ⁇ 0.5*2000*(2 ⁇ *0.050/20) ) ⁇ longer than the center core in a distance of 2 m, if each helix period is about 20 mm and the distance between the center core and an outer core is about 50 ⁇ .
  • the mismatch between sensors from different cores can be compensated, if it is possible to precisely measure the refractive index of each core and the helix period.
  • it requires accuracy better than 2* 10 "6 , in case the mismatch between sensors of a 2 m long sensing fiber should always be less than one fiber index. This is an extremely high accuracy for the measurement of the reflective index in practice.
  • the invention proposes to use the distance information between the start and end peaks of the FFT signal, as shown in FIG. 5, to compensate for the impact of the refractive index and the helix structure, in response to measured optical distances for the fiber cores. Based on formula (1), the sensor resolution is then accordingly adjusted by
  • FIG. 6 illustrates an example of a sub-index detection of a peak which is possible to implement in order to improve accuracy of detection of optical length by means of start and end peak detections.
  • the illustration shows a peak at integer index k, and the maximum estimation error ⁇ of the real peak is at most half an index. Consequently, the length measurement of each core can have maximum one index error.
  • the length of each core can be calculated in the units of fiber index.
  • FIG. 7 shows the result of 10000 consecutive measurements on the length of the center core in an optical shape sensor. It shows that the measurement result is rather stable.
  • the noise band is within ⁇ 0.05 fiber index.
  • There is a slope on the measured length which is caused by the change of the environment temperature (-0.04°). The same trend is observed in the length measurement of the other three outer cores. However, it has been found that if the length ratio between an outer core and the center core is calculated, the temperature influence is cancelled.
  • FIG. 8 shows that the length ratio is always centered at one value, here
  • a preferred measure of optical length of an outer core is to calculate the ratio between the measure optical length of the outer core and the measured optical length of the center core.
  • FIG. 9 serves to illustrate robustness of the optical length ratio for variations over time.
  • the graph shows the result compared with an experiment carried out 2 weeks earlier. It shows that the average length ratio between cores has only very little variance
  • the proposed optical length ratio has also proven to be stable with respect to different shapes.
  • the ratio has a change of 1.4* 10 "6 (from 1.0001706 to 1.0001720) when the straight shape is deformed to a curved shape indicated. Furthermore, the ratio is stable with respect to temperature change. When the temperature changes from 42.9° to 27.7°, the ratio variance is also within 10 "6 .
  • the interfering signals from four fiber cores are sampled by hardware, e.g. using a data acquisition card, and further post-processed in software. Therefore, there are various ways to implement the proposed approach in software or hardware. A few examples are illustrated below.
  • a hardware sampling technique can be applied to sample the interfering signals with equal steps of the laser frequency change (linear K-space sampling). In case the ratios of scan wavelength range between cores is known based on the calibration measurement, the ratios of the frequency steps between cores is also known for a fixed laser scan range.
  • FIG. 10 shows, for instance, after the resampling, the center core has 2 million data points SMP corresponding to the optical scan wavelength range from 1532 nm to 1552 nm. If the length ratio between an outer core and the center core is 1.0001, it means that the scan range of the outer-core is 1/1.0001 of that of the center core.
  • the scan range of this outer-core can be reduced by removing part of samples at the beginning or at the end of the scan. In this example, we can remove either 198 points at the beginning or 200 points at the end. In order to keep the center optical frequency the same, it is also possible to remove a part of points from the beginning and a part from the end.
  • FIG. 1 1 shows amplitude A versus fiber index F_I of an example of end peaks of four fiber cores aligned within one fiber index F_I range when the scan range of each core is adjusted according to their calibration measurements. This is a unique feature of the approach of the present invention.
  • solution b) is similar to the hardware implementation, the resampling step of each core can be adjusted accordingly in software implementation.
  • FIGs. 12a-12c shows shapes that have been tested in an existing OSS system where the invention has been implemented: one circle, an 8-shape, and two overlapping circles. All shapes are essentially 2 dimensional.
  • the out-of plane errors show the accuracy of the reconstructed shapes.
  • the errors have been compared with the same shape between the aligned case (i.e. according to the invention) and in the non-aligned case (standard).
  • the results (out-of-plane error in mm) are for the aligned case: 1.0; 2.5; and 4.2 for the respective three shapes of FIG. 12, while for the non-aligned case the results are: 5.5; 2.9; and 5.5.
  • the errors of the tested shapes are always less when the cores are aligned.
  • FIG. 13a shows a photo of an optical shape sensor placed to form one big half circle statically on a table.
  • FIGs. 13b and 13c show results of tested jitter of the system. The jitter is checked by measuring the position variations of the end tip of the optical shape sensor during continuous scanning. The shapes from 100 scans have been recorded and plotted as position distribution of the end tip projected on the XY plane.
  • FIG. 13a shows for the non- aligned case, and FIG. 13b with alignment according to the invention.
  • the calculated standard errors (mm) for the jitter in three directions (X, Y and Z) are: 0.412; 0.609; and
  • the proposed approach according to the invention aligns different fiber cores based on raw signals, just as explained previously.
  • FFT signals Linear interpolation FFT signals Linear interpolation
  • FFT signals Spline interpolation FFT signals Spline interpolation
  • Phase signals Linear interpolation Phase signals Linear interpolation
  • Phase signals Spline interpolation Phase signals Spline interpolation
  • FIGs. 14a and 14b illustrate an example of the variation of refractive index and the helix structure in the optical lengths of cores which can be different when the laser scans in the same optical scan wavelength range. This is visualized by the FFT signals where end peaks of all four fiber cores are seen. In FIG. 14a, the FFT signals locate at different fiber indexes for all cores (the length of 2 samples SMP is shown by the double arrows). In
  • FIG. 14b the end peaks in FFT signals locate almost at the same fiber index when their raw signals have been adjusted according to proper optical scan wavelength ranges according to the invention.
  • FIG. 15 shows steps of an embodiment for calibrating an optical console system arranged for optical interrogation of an optical shape sensor comprising a plurality of optical fiber cores arranged between a proximal end and a distal end.
  • the method comprises performing optical calibration measurements P CM on the plurality of optical fiber cores with one common optical scan wavelength range for all of the optical fiber cores.
  • Next step is calculating C OL, for each of the optical fiber cores, a measure of its optical length in response to the calibration measurements. Preferably, this is done by calculating C OL R optical length ratios for the measured optical lengths of each of the outer cores and the measured length for the center core.
  • step is determining D I SWR an individual optical scan wavelength range for each of the outer cores in response to the individual measures of their optical lengths relative to the measure of optical length of the center core, so as to compensate for optical length differences between the plurality of optical fiber cores.
  • the system is prepared for providing OSS P OSS, i.e. to perform optical interrogation of the plurality of optical fibers to perform 3D shape reconstruction of the optical shape sensor improved accuracy.
  • the invention provides an OSS system with an optical console system P, C arranged for optical interrogation of the optical fiber cores CC, CI, C2, C3 in the associated optical shape sensor OSF.
  • the optical console system P, C is arranged to perform optical calibration measurements on the fiber cores CC, C 1 , C2, C3 with one common optical scan wavelength range ⁇ C.
  • a measure of its optical length is calculated based on the results of the calibration measurements, preferably as ratios between outer core optical lengths and center core optical length.
  • Individual optical scan wavelength ranges ⁇ , ⁇ 2, ⁇ 3 for each of the outer cores CI, C2, C3 are then determined according to their individual optical lengths relative to the optical length of the center core CC, so as to compensate for optical length differences between the plurality of optical fiber cores CC, CI, C2, C3.
  • OSS interrogation can be carried out with the determined individual optical scan wavelength ranges ⁇ , ⁇ 2, ⁇ 3 for each of the outer cores CI , C2, C3, and a high OSS accuracy can be obtained since both physical length differences and refractive index differences have been aligned.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

An optical shape sensing (OSS) system with an optical console system (P, C) arranged for optical interrogation of the optical fiber cores (CC, C1, C2, C3) in the associated optical shape sensor (OSF). The optical console system (P, C) is arranged to perform optical calibration measurements on the fiber cores (CC, C1, C2, C3) with one common optical scan wavelength range (∆λ C). For each of the optical fiber cores (CC, C1, C2, C3), a measure of its optical length in calculated based on the results of the calibration measurements, preferably as ratios between outer core optical lengths and center core optical length. Individual optical scan wavelength ranges (∆λ1, ∆λ2, ∆λ3) for each of the outer cores (C1, C2, C3) are then determined according to their individual optical lengths relative to the optical length of the center core (CC), so as to compensate for optical length differences between the plurality of optical fiber cores (CC, C1, C2, C3). With this way of calibrating the OSS system, OSS interrogation can be carried out with the determined individual optical scan wavelength ranges (∆λ1, ∆λ2, ∆λ3) for each of the outer cores (C1, C2, C3), and a high OSS accuracy can be obtained since both physical length differences and refractive index differences have been aligned.

Description

Multi-core alignment in Optical Shape Sensing
FIELD OF THE INVENTION
The present invention relates to the field of optical shape sensing (OSS). Especially, the invention provides an OSS system and a method for alignment of multiple optical fiber cores in an OSS device.
BACKGROUND OF THE INVENTION
In optical shape sensing (OSS), an optical shape sensor often has a set of helix outer cores and a center core. When the optical shape sensor is curved, each of the outer cores will deform differently. By means of an optical console system providing laser interferometric interrogation of the optical fiber cores in the optical shape sensor, it is possible to provide a 3D shape reconstruction of the optical shape sensor, and thus allows tracking of an elongated device which the optical shape sensor is mounted to follow, e.g. to help navigating a medical interventional device.
Since OSS relies on different deforms of the individual optical fiber cores, alignment between the different fiber cores is a key aspect to improve the accuracy of the shape reconstruction. Typically, misalignments between the fiber cores are caused by their different refractive index and the helix structure.
US 2012/0069347 Al describes this misalignment issue and proposes to realign different cores by resampling. However, how the resampling is actually performed is not described.
SUMMARY OF THE INVENTION
Following the above description, it would be advantageous to provide an OSS system with an improved accuracy. Preferably, the OSS system has an increased tolerance towards physical length and refractive index differences of cores in an optical fiber in a multi-core OSS device, i.e. an improved alignment of multiple cores embedded in an optical fiber is advantageous. Still, it is preferred that the solution also leads to OSS with reduced jitter. Furthermore, the solution is preferably stable over time, and with respect to temperature changes. In a first aspect, the invention provides an optical shape sensing system arranged to determine a measure of a three-dimensional shape of an associated optical shape sensor comprising a plurality of optical fiber cores arranged between a proximal end and a distal end, wherein the plurality of optical fiber cores comprise optical shape sensing properties, wherein a center core is one of the plurality of optical fiber cores arranged centrally to a plurality of outer cores, wherein the system comprises an optical console system arranged for optical interrogation of the optical fiber cores in the associated optical shape sensor, wherein the optical console system is arranged:
- a) to perform optical calibration measurements on the plurality of optical fiber cores with one common optical scan wavelength range for all of the optical fiber cores,
- b) to calculate, for each of the optical fiber cores, a measure of its optical length in response to the calibration measurements,
- c) to determine an individual optical scan wavelength range for each of the outer cores in response to the individual measures of their optical lengths relative to the measure of optical length of the center core, so as to compensate for optical length differences between the plurality of optical fiber cores, and
- d) to apply the individual optical scan wavelength ranges for each of the outer cores in optical interrogation of the plurality of optical fiber cores, and to determine a measure of three-dimensional shape of at least a part of the optical shape sensor accordingly.
Such OSS device with calibration steps a)-c) is advantageous, since it is possible to align physical length as well as refractive differences between multiple optical fiber cores of an optical shape sensor. This improves OSS accuracy, and the calibration measurements only introduce little extra calculation time which can easily satisfy the real- time requirements of an OSS system.
Experiments have shown that the calibration method is reliable for normal operational situations. Secondly, the proposed re-alignment approach performs better than normal resampling approaches. Experiments show that both the shape accuracy and the system jitter are improved with the proposed approach. Furthermore, experiments show that the improved accuracy results obtained are stable over time, and over temperature changes.
In the following, embodiments and/or preferred features are described.
The optical console system may be arranged to determine peak positions of start and end points for all of the optical fiber cores in response to a discrete spectral representation, e.g. in the form of a Fast Fourier Transform (FFT), of the optical calibration measurements, and to calculate the individual measures of optical lengths as distances between start and end points in response to the determined peak positions. This is a rather simple way of determining optical lengths of each of the fiber cores in response to the calibration measurements with one common optical scan wavelength range for all fiber cores. Thereby, the difference in optical length between the different fiber cores is also known. The peak positions may be determined with a sub-index precision by applying an estimation algorithm to the discrete spectral representation (FFT) of the optical calibration
measurements. Such estimation algorithms are known in the art.
The individual optical scan wavelength ranges for each of the outer cores may be determined so as to equal a recognized number of locations along the fiber core of each of the plurality of optical fiber cores, since the determined optical lengths for each of the plurality of optical fiber cores may be considered as an expression of a finite number of single locations along the fiber core which can be optically recognized. This can be equalized between the fiber cores, if their individual optical scan wavelength ranges are adjusted so as to align the optical lengths of the fiber cores.
Preferably, the optical console system is arranged to calculate a measure of optical length ratio for each of the outer cores. Preferably, the measure of optical length ratio for each of the outer cores is calculated as a ratio between the calculated measure of optical length of the outer core and the calculated measure of optical length of the center core, and to determine the individual optical scan wavelength ranges for each of the outer cores in response to the calculated measure of optical length ratios. This optical length ratio is found to provide a robust measure suited for the aligning, i.e. by determining the individual optical scan wavelength ranges.
At least some of a), b), and c) may be implemented in hardware and/or in software. Especially, some of a), b), and c) may be implemented in software, while other part of a), b), and c) are implemented in hardware. Especially, all of a), b), and c) may be implemented in software.
The optical console system may be arranged to select a limited portion of sampled data points obtained for an outer core, so as to reduce an individual scan wavelength range for the outer core. This is a rather simple software implementable way of aligning the fiber cores. Especially, the limited portion of sampled data points may be selected to be in the beginning, in the end, or in the middle of the sampled data points.
The optical console system preferably comprises a controllable laser light source system arranged to generate light for the optical interrogation of the plurality of optical fiber cores with a light wavelength sweep within an individually controllable optical scan wavelength range for each of the plurality of optical fiber cores. The optical console system may in a specific embodiment comprise a laser light source system arranged to generate a light wavelength sweep with a wavelength range within the range 1525-1565 nm. It is to be understood that this is by no means limiting, and the skilled person will know how to translate this to other wavelength ranges as well and optical scan wavelength ranges accordingly, and thus the invention is applicable for a wide variety of laser source systems or other light source systems.
The system may comprise an optical shape sensor comprising a plurality of optical fiber cores comprising optical shape sensing properties and being arranged between a proximal end and a distal, wherein the plurality of optical fiber cores comprises a center core and a plurality of outer cores, wherein the center core is arranged centrally to the plurality of outer cores being wound around the center core. Especially, the optical shape sensor may be arranged to follow a shape of an elongated medical interventional device, and wherein the system is arranged to display an image generated in response to said measure of three- dimensional shape of the optical shape sensor.
The associated optical shape sensor preferably has fiber cores with optical shape sensing properties comprising backscattering properties used for optical shape sensing. Especially, the optical fibers may comprise at least one of the properties: Rayleigh scattering and fiber Bragg gratings. Especially, the multi-core optical fiber may be arranged within one common coating or cover.
Especially, the optical shape sensor may be mounted to follow the shape of an elongated device, e.g. a medical interventional device, e.g. in the form of a catheter, a guide wire, an endoscope etc. However, it is to be understood that the invention is applicable also for non-medical use.
The number of optical fiber cores in the associated optical shape sensor may be four, i.e. one center core and three outer cores. However, it is to be understood that other number of optical fiber cores can be used, e.g. a total of three: one center core and two outer cores, e.g. for use in robotic applications with a non-twistable robotic arm. In other applications, four or more than four outer cores may be used together with one center core.
On the one hand, it is convenient to keep the amount of data points constant over the full length of the optical shape sensor. In particular a power of 2 may be chosen in order to perform fast Fourier transformations (FFT) and make optimal use of parallel computing during the remaining signal processing steps. On the other hand, it may be preferred to use sensors of different lengths. By adjusting the step size in wavelength during data acquisition and consequently adjusting the total wavelength span it is possible to keep the amount of data points the same for different lengths of sensors. For an increase in sensor length, one needs to decrease the wavelength span proportionally. When length variations are large, the variation in wavelength span will also become large. Since strain causes the spectra to shift in wavelength, it means that adjustment of the wavelength span will cause a change in maximum amount of allowable strain, i.e. a change in minimum bending radius. This extension of the invention means that it is possible to exchange length variations of the sensors for variations in minimum bending radius in vice versa.
In a second aspect, the invention provides a method for calibrating an optical console system arranged for optical interrogation of an optical shape sensor comprising a plurality of optical fiber cores arranged between a proximal end and a distal, wherein the plurality of optical fiber cores comprise optical shape sensing properties, wherein a center core is one of the plurality of optical fiber cores arranged centrally to a plurality of outer cores, the method comprising
- a) performing optical calibration measurements on the plurality of optical fiber cores with one common optical scan wavelength range for all of the optical fiber cores,
- b) calculating, for each of the optical fiber cores, a measure of its optical length in response to the calibration measurements, and
- c) determining an individual optical scan wavelength range for each of the outer cores in response to the individual measures of their optical lengths relative to the measure of optical length of the center core, so as to compensate for optical length differences between the plurality of optical fiber cores.
In a third aspect, the invention provides a computer program product having instructions which when executed cause an optical console system to perform the method according to the second aspect.
It is appreciated that the same advantages and embodiments of the first aspect apply as well for the second and third aspects. In general the first, second and third aspects may be combined and coupled in any way possible within the scope of the invention. These and other aspects, features and/or advantages of the invention will be apparent from and elucidated with reference to the embodiments described hereinafter. BRIEF DESCRIPTION OF THE DRAWINGS
Embodiments of the invention will be described, by way of example only, with reference to the drawings, in which
FIG. 1 shows sketch of an OSS system embodiment,
FIG. 2 illustrates the basic principle of strain measurements used in OSS,
FIG. 3 illustrates a typical optical shape sensor for OSS with one center core and three outer cores wound around the center core,
FIG. 4 illustrates the principle behind OSS: optical measurement of compression and tension of various parts of a curved optical shape sensor with multiple optical fiber cores,
FIG. 5 illustrates an example of measured optical lengths for four optical fiber cores forming part of one optical shape sensor,
FIG. 6 illustrates an example of sub-index peak detection,
FIG. 7 shows a graph of measured length of a center core repeated 10000 times,
FIG. 8 shows a graph of ratio between an outer core and the center core, repeated 10000 times,
FIG. 9 shows a graph similar to FIG. 8 but carried out 2 weeks earlier, FIG. 10 shows a graph of example data points and their corresponding laser wavelength (here an optical scan wavelength range of 20 nm),
FIG. 1 1 shows fiber index number versus amplitude for 4 fiber cores with aligned end peaks,
FIGs. 12a- 12c shows different experimental shapes tested for an optical shape sensor,
FIGs. 13a shows a photo of an experiment setup, and FIGs. 13b and 13c show results of optical shape sensor tip position distribution for 100 scans without and with alignment according to the invention,
FIGs. 14a and 14b show results for end peaks of four cores as fiber index number versus amplitude for the measured data with and without alignment, and
FIG. 15 illustrates steps of a method embodiment.
DESCRIPTION OF EMBODIMENTS
FIG. 1 illustrates basic parts of an OSS system with an optical shape sensor OSF comprising a plurality of optical fiber cores arranged between a proximal end and a distal end where the optical fiber cores form part of the same fiber. The plurality of optical fiber cores comprise optical shape sensing properties, e.g. fiber Bragg gratings, so as to allow OSS of at least a part of their lengths. A center core is one of the plurality of optical fiber cores which is arranged centrally to a plurality of outer cores, e.g. three outer cores, wound around one center core, thereby in itself introducing a physical length difference between the center core and the outer cores. The optical shape sensor OSF incorporated in an elongated device T, e.g. a medical catheter or guide wire. An optical console C has a controllable laser light source capable of providing a light sweep, and an optical detector for detecting an optical response to such light sweep. The optical console C is connected to the optical shape sensor OSF and arranged to optically interrogate the strain sensing optical therein, and to accordingly determine a measure of a three-dimensional shape of at least a part of the optical fiber OSF and thereby the elongated device T in which the optical shape sensor OSF is placed.
A processor P controls the optical console C, and a 3D image I of the optical shape sensor OSF can be generated, e.g. displayed as an image on a monitor in real time. The processor P executes a control algorithm CA which serves to perform a calibration procedure with the purpose of aligning differences in physical length and refractive index between the plurality of optical fiber cores that constitute the optical shape sensor OSF, so as to improve OSS accuracy. The calibration procedure comprises performing optical calibration measurements on the plurality of optical fiber cores with one common optical scan wavelength range ΔλC for all of the optical fiber cores. E.g. an initial ΔλC can be selected to be such as 20 nm, e.g. corresponding to scanning from 1532 nm to 1552 nm. Based on these calibration measurements, a measure of optical length of each of the optical fiber cores is calculated. E.g. by determining, for each of the fiber cores, fiber index number between detected start and end peaks in the calibration measurement data sample results. Individual optical scan wavelength ranges Δλΐ, Δλ2, Δλ3 for each of the outer cores are then calculated in response to the individual measures of their optical lengths relative to the measure of optical length of the center core.
Hereby, it is possible to compensate for optical length differences between the plurality of optical fiber cores. These individual optical scan wavelength ranges Δλΐ, Δλ2, Δλ3 for each of the outer cores can then be applied in optical interrogation of the plurality of optical fiber cores, and to determine a measure of three-dimensional shape of at least a part of the optical shape sensor OSF accordingly. For the center core, the originally used optical scan wavelength range ΔλC for the calibration measurements can be used. Hereby, an increased OSS accuracy can be obtained, since the introduced individual scan wavelength ranges for the fiber cores serve to align them at a basic signal processing level. Such calibration procedure is easy to implement in existing equipment, e.g. in pure software, and the extra time required for the calibration procedure does not compromise real-time OSS. Experiments show that both the shape accuracy and system jitter can be improved with the proposed calibration procedure.
FIG. 2 illustrates the principle behind detecting of position of 4 discrete optical sensors or sensor elements SI, S2, S3, and S4 on an optical fiber OF. A swept laser source LS emits light into a single mode optical fiber OF via a circulator C. The response is detected via the circulator C at an optical detector DT which also receives the light emitted from the laser source LS via a reference arm RF A. The light is reflected at each point SI, S2, S3, and S4 along the fiber. If strain is exerted on the fiber OF, optical interferometry allows measuring the strain-dependent phase shift with high spatial resolution as a function of distance along the fiber OF. The reflected light at each point SI, S2, S3, and S4 interferes with the light from the reference arm RF A. The interfered signal is in the form of a sine wave for each sensor SI, S2, S3, and S4 if the laser scans linearly over frequency. Although the summation of the interfered signals has not obvious relations with locations of the sensors SI, S2, S3, and S4, the FFT of the summation signal can clearly identify positions of sensors SI, S2, S3, and S4. This is illustrated to the right where an example of position PS versus amplitude A is shown as four peaks corresponding to the positions of the four sensors S 1 , S2, S3, and S4.
FIG. 3 shows a typical example of an optical shape sensor for OSS applications with four single-mode optical fiber cores: one center core CC, and three outer cores CI, C2, C3 helix would around the center core CC. The upper illustration shows a longitudinal section of the optical shape sensor, while the lower illustration shows a cross sectional sketch showing the central position of the center core CC relative to the outer cores CI, C2, C3 within one common jacket JKT.
In case the end of the fiber is cleaved with a small angle, e.g. 8°, the four cores CC, CI, C2, C3 can end at different locations. However, these length differences are very small. In the worst scenario, an outer core can end a few μηι away from that of the center core. E.g., if the distance between the center core CC and outer cores CI, C2, C3 is 50 μηι and angle cleave is 8 °, the maximum difference between the outer core CI, C2, C3 and the center core CC is at most 7 μηι (8*50*π/180) in physical space. This small difference can be neglected, so that in the physical space, all four cores CC, CI, C2, C3 can be assumed to start and to end at equal positions.
FIG. 4 shows an illustration of the optical shape sensor from FIG. 3 in the situation when it is bend, i.e. when the fiber structure is curved, and thus sensors on three outer-cores behave differently. Some will be compressed CMP and some will be stretched due to tension TNS. When the sensors from different cores are properly aligned, the measured phase differences from the aligned sensors can be translated via strains to a full 3- D shape.
Ideally, an optical fiber with optical shape sensing elements can be considered as a series of continuous sensors along its length. However, in practice, only a finite number of sensors can be recognized by an optical interrogation system. The resolution of a sensor (Δζ) is determined by,
Figure imgf000010_0001
where, λ is the center wavelength of the scan (e.g. 1542 nm), Δλ is the optical scan wavelength range (e.g. 20 nm from 1532 nm to 1552 nm), and n is the refractive index of the core. When n is about 1.49, Δζ is about 40 μηι in this example.
When a sensing fiber lies straight without applying local strains, sensors from each core distribute evenly along its physical length. However, due to different refractive index of cores in the fiber, the sensor resolutions vary from core to core. Furthermore, due to the helix structure, the physical lengths of the three outer cores are slightly longer than the center one. Therefore, sensors at the same physical location along the sensing fiber from different cores are not naturally matched with each other.
FIG. 5 shows an example illustrating this mismatch between the different cores. A spectral representation in the form of FFT signals of the sensors from different cores shows amplitude versus fiber index F_I. Typically, at the start and end points of each core, there are large peaks in the FFT signals due to the significant reflections at those points. Between the start and end peaks, the plot reveals the FFT signal of the sensors on each of the cores. When all cores are aligned at the start peaks, the end peaks from 4 cores locate at different fiber index numbers F_I (the signals are sampled in a discretized manner and hence its FFT will be discretized). In this example, the end of the center core CC and that of an outer core CI have a difference of 13 fiber index numbers, which approximately corresponds 0.5 mm (13*40 μηι) in the physical space. To improve the 3D shape accuracy of the sensing fiber, the sensors from different cores need to be aligned properly. The proposed calibration procedure according to the invention serves to provide this alignment, namely based on an optical calibration measurement of the optical length difference between the cores, thereby resulting in improved OSS accuracy.
Since all cores have the same physical distance between start and end positions without considering the minor difference caused by the angle cleave, the sensors from different cores are better aligned when the number of fiber index points for each core between the start and end points is the same. As mentioned, two factors influence the number of points:
1) The refractive index n, see e.g. formula (1), where the sensor resolution Δζ is determined by λ, Δλ, and n. If each core has a different refractive index, the sensor resolution can vary between cores. Consequently, the number of points can be different for the same physical length.
2) Due to the helix structure, the outer core is slight longer than the center core. For instance, an outer core can be about 0.25 mm {0.5*2000*(2π*0.050/20) )} longer than the center core in a distance of 2 m, if each helix period is about 20 mm and the distance between the center core and an outer core is about 50 μηι.
In principle, the mismatch between sensors from different cores can be compensated, if it is possible to precisely measure the refractive index of each core and the helix period. In the OSS application, it requires accuracy better than 2* 10"6, in case the mismatch between sensors of a 2 m long sensing fiber should always be less than one fiber index. This is an extremely high accuracy for the measurement of the reflective index in practice. Thus, instead the invention proposes to use the distance information between the start and end peaks of the FFT signal, as shown in FIG. 5, to compensate for the impact of the refractive index and the helix structure, in response to measured optical distances for the fiber cores. Based on formula (1), the sensor resolution is then accordingly adjusted by
individually adjusting the optical scan wavelength range of the laser. By introducing such individual optical scan wavelength range for each fiber core, the number of points of each fiber core over the full length can be adjusted to equal.
FIG. 6 illustrates an example of a sub-index detection of a peak which is possible to implement in order to improve accuracy of detection of optical length by means of start and end peak detections. The illustration shows a peak at integer index k, and the maximum estimation error δ of the real peak is at most half an index. Consequently, the length measurement of each core can have maximum one index error. In literature, there are various ways to estimate the peak position of an FFT signal with sub-index accuracy, such as known by the skilled person. Thus improved accuracy can be obtained by applying such method. Based on the peak measurement, the length of each core can be calculated in the units of fiber index.
FIG. 7 shows the result of 10000 consecutive measurements on the length of the center core in an optical shape sensor. It shows that the measurement result is rather stable. The noise band is within ±0.05 fiber index. There is a slope on the measured length, which is caused by the change of the environment temperature (-0.04°). The same trend is observed in the length measurement of the other three outer cores. However, it has been found that if the length ratio between an outer core and the center core is calculated, the temperature influence is cancelled.
FIG. 8 shows that the length ratio is always centered at one value, here
1.0001727, over 10000 repeated measurements. The maximum variation of the ratio is about 7* 10"7, which is about a mismatch of 0.03 fiber index (7* 10"7*45000 fiber index) between the outer and the center core. Thus, a preferred measure of optical length of an outer core is to calculate the ratio between the measure optical length of the outer core and the measured optical length of the center core.
FIG. 9 serves to illustrate robustness of the optical length ratio for variations over time. The graph shows the result compared with an experiment carried out 2 weeks earlier. It shows that the average length ratio between cores has only very little variance
(6* 10" ) between the two batches of measurements.
The proposed optical length ratio has also proven to be stable with respect to different shapes. The ratio has a change of 1.4* 10"6 (from 1.0001706 to 1.0001720) when the straight shape is deformed to a curved shape indicated. Furthermore, the ratio is stable with respect to temperature change. When the temperature changes from 42.9° to 27.7°, the ratio variance is also within 10"6.
Knowing the optical length ratios between the outer cores and the center core, it is possible to adjust the scan range Δλ for each core based on formula (1) such that the same number of locations (points) for each core between the start and the end peaks are sampled.
During a laser scan, the interfering signals from four fiber cores are sampled by hardware, e.g. using a data acquisition card, and further post-processed in software. Therefore, there are various ways to implement the proposed approach in software or hardware. A few examples are illustrated below.
1) Hardware implementation. A hardware sampling technique can be applied to sample the interfering signals with equal steps of the laser frequency change (linear K-space sampling). In case the ratios of scan wavelength range between cores is known based on the calibration measurement, the ratios of the frequency steps between cores is also known for a fixed laser scan range.
2) Software implementation. When a generic data acquisition card is used for the data sampling, the interfering signals are sampled with the same time interval. The data needs to be resampled in software to ensure that they have the same amount of change in laser frequency between any two adjacent samples. The resampling of the data can be done by the software in FPGA, GPU or CPU. In one solution: a) The data of all cores are first sampled with the same frequency step in software. To obtain the various scan ranges, the amount of the resampled data used for FFT conversion can be adjusted for different cores.
FIG. 10 shows, for instance, after the resampling, the center core has 2 million data points SMP corresponding to the optical scan wavelength range from 1532 nm to 1552 nm. If the length ratio between an outer core and the center core is 1.0001, it means that the scan range of the outer-core is 1/1.0001 of that of the center core. The scan range of this outer-core can be reduced by removing part of samples at the beginning or at the end of the scan. In this example, we can remove either 198 points at the beginning or 200 points at the end. In order to keep the center optical frequency the same, it is also possible to remove a part of points from the beginning and a part from the end.
FIG. 1 1 shows amplitude A versus fiber index F_I of an example of end peaks of four fiber cores aligned within one fiber index F_I range when the scan range of each core is adjusted according to their calibration measurements. This is a unique feature of the approach of the present invention.
As an alternative to the above mentioned solution a), solution b) is similar to the hardware implementation, the resampling step of each core can be adjusted accordingly in software implementation.
FIGs. 12a-12c shows shapes that have been tested in an existing OSS system where the invention has been implemented: one circle, an 8-shape, and two overlapping circles. All shapes are essentially 2 dimensional. The out-of plane errors show the accuracy of the reconstructed shapes. The errors have been compared with the same shape between the aligned case (i.e. according to the invention) and in the non-aligned case (standard). The results (out-of-plane error in mm) are for the aligned case: 1.0; 2.5; and 4.2 for the respective three shapes of FIG. 12, while for the non-aligned case the results are: 5.5; 2.9; and 5.5. Thus, the errors of the tested shapes are always less when the cores are aligned.
FIG. 13a shows a photo of an optical shape sensor placed to form one big half circle statically on a table. FIGs. 13b and 13c show results of tested jitter of the system. The jitter is checked by measuring the position variations of the end tip of the optical shape sensor during continuous scanning. The shapes from 100 scans have been recorded and plotted as position distribution of the end tip projected on the XY plane. FIG. 13a shows for the non- aligned case, and FIG. 13b with alignment according to the invention. The calculated standard errors (mm) for the jitter in three directions (X, Y and Z) are: 0.412; 0.609; and
0.1 12 for the non-aligned case, while 0.377; 0.464; and 0.064 for the aligned case. Thus, the three directions are all reduced about 10%, 20% and 50% respectively.
The proposed approach according to the invention aligns different fiber cores based on raw signals, just as explained previously. There are other ways to align the length of cores by resampling their FFT signals or phase signals afterwards. Examples of other methods are: FFT signals Linear interpolation, FFT signals Spline interpolation, Phase signals Linear interpolation, and Phase signals Spline interpolation. The off plane errors of each approach for four different shapes have been tested, and it has been found that the proposed approach according to the invention gives better and more robust improvement on the accuracy of the system.
FIGs. 14a and 14b illustrate an example of the variation of refractive index and the helix structure in the optical lengths of cores which can be different when the laser scans in the same optical scan wavelength range. This is visualized by the FFT signals where end peaks of all four fiber cores are seen. In FIG. 14a, the FFT signals locate at different fiber indexes for all cores (the length of 2 samples SMP is shown by the double arrows). In
FIG. 14b, the end peaks in FFT signals locate almost at the same fiber index when their raw signals have been adjusted according to proper optical scan wavelength ranges according to the invention.
FIG. 15 shows steps of an embodiment for calibrating an optical console system arranged for optical interrogation of an optical shape sensor comprising a plurality of optical fiber cores arranged between a proximal end and a distal end. The method comprises performing optical calibration measurements P CM on the plurality of optical fiber cores with one common optical scan wavelength range for all of the optical fiber cores. Next step is calculating C OL, for each of the optical fiber cores, a measure of its optical length in response to the calibration measurements. Preferably, this is done by calculating C OL R optical length ratios for the measured optical lengths of each of the outer cores and the measured length for the center core. Next, step is determining D I SWR an individual optical scan wavelength range for each of the outer cores in response to the individual measures of their optical lengths relative to the measure of optical length of the center core, so as to compensate for optical length differences between the plurality of optical fiber cores. Finally, with the hereby calibrated OSS system, the system is prepared for providing OSS P OSS, i.e. to perform optical interrogation of the plurality of optical fibers to perform 3D shape reconstruction of the optical shape sensor improved accuracy.
To sum up, the invention provides an OSS system with an optical console system P, C arranged for optical interrogation of the optical fiber cores CC, CI, C2, C3 in the associated optical shape sensor OSF. The optical console system P, C is arranged to perform optical calibration measurements on the fiber cores CC, C 1 , C2, C3 with one common optical scan wavelength range ΔλC. For each of the optical fiber cores CC, CI, C2, C3, a measure of its optical length in calculated based on the results of the calibration measurements, preferably as ratios between outer core optical lengths and center core optical length.
Individual optical scan wavelength ranges Δλΐ, Δλ2, Δλ3 for each of the outer cores CI, C2, C3 are then determined according to their individual optical lengths relative to the optical length of the center core CC, so as to compensate for optical length differences between the plurality of optical fiber cores CC, CI, C2, C3. With this way of calibrating the OSS system, OSS interrogation can be carried out with the determined individual optical scan wavelength ranges Δλΐ, Δλ2, Δλ3 for each of the outer cores CI , C2, C3, and a high OSS accuracy can be obtained since both physical length differences and refractive index differences have been aligned.
While the invention has been illustrated and described in detail in the drawings and foregoing description, such illustration and description are to be considered illustrative or exemplary and not restrictive; the invention is not limited to the disclosed embodiments. Other variations to the disclosed embodiments can be understood and effected by those skilled in the art in practicing the claimed invention, from a study of the drawings, the disclosure, and the appended claims. In the claims, the word "comprising" does not exclude other elements or steps, and the indefinite article "a" or "an" does not exclude a plurality. A single processor or other unit may fulfill the functions of several items recited in the claims. The mere fact that certain measures are recited in mutually different dependent claims does not indicate that a combination of these measured cannot be used to advantage. A computer program may be stored/distributed on a suitable medium, such as an optical storage medium or a solid-state medium supplied together with or as part of other hardware, but may also be distributed in other forms, such as via the Internet or other wired or wireless
telecommunication systems. Any reference signs in the claims should not be construed as limiting the scope.

Claims

CLAIMS:
1. An optical shape sensing system arranged to determine a measure of a three- dimensional shape of an associated optical shape sensor (OSF) comprising a plurality of optical fiber cores (CC, CI, C2, C3) arranged between a proximal end and a distal end, wherein the plurality of optical fiber cores (CC, CI, C2, C3) comprise optical shape sensing properties, wherein a center core (CC) is one of the plurality of optical fiber cores (CC, CI, C2, C3) arranged centrally to a plurality of outer cores (CI, C2, C3), wherein the system comprises an optical console system (P, C) arranged for optical interrogation of the optical fiber cores (CC, CI, C2, C3) in the associated optical shape sensor (OSF), wherein the optical console system (P, C) is arranged:
- a) to perform optical calibration measurements on the plurality of optical fiber cores (CC, CI, C2, C3) with one common optical scan wavelength range (Δλ^ for all of the optical fiber cores (CC, C1, C2, C3),
- b) to calculate, for each of the optical fiber cores (CC, CI, C2, C3), a measure of its optical length in response to the calibration measurements,
- c) to determine an individual optical scan wavelength range (Δλΐ, Δλ2, Δλ3) for each of the outer cores (CI, C2, C3) in response to the individual measures of their optical lengths relative to the measure of optical length of the center core (CC), so as to compensate for optical length differences between the plurality of optical fiber cores (CC, CI, C2, C3), and
- d) to apply the individual optical scan wavelength ranges (Δλΐ, Δλ2, Δλ3) for each of the outer cores (CI, C2, C3) in optical interrogation of the plurality of optical fiber cores (CC,
CI, C2, C3), and to determine a measure of three-dimensional shape of at least a part of the optical shape sensor (OSF) accordingly.
2. System according to claim 1, wherein the optical console system (P, C) is arranged to determine peak positions of start and end points for all of the optical fiber cores
(CC, CI, C2, C3) in response to a discrete spectral representation (FFT) of the optical calibration measurements, and to calculate the individual measures of optical lengths as distances between start and end points in response to the determined peak positions.
3. System according to claim 2, wherein said peak positions are determined with a sub-index precision by applying an estimation algorithm to the discrete spectral representation (FFT) of the optical calibration measurements.
4. System according to claim 1, wherein the individual optical scan wavelength ranges (Δλΐ, Δλ2, Δλ3) for each of the outer cores are determined so as to equal a recognized number of index points along the optical fiber core of each of the plurality of optical fiber cores (CC, C1, C2, C3).
5. System according to claim 1, wherein the optical console system (P, C) is arranged to calculate a measure of optical length ratio for each of the outer cores (CI, C2, C3), wherein the measure of optical length ratio for each of the outer cores (CI, C2, C3) is calculated as a ratio between the calculated measure of optical length of the outer core (CI, C2, C3) and the calculated measure of optical length of the center core (CC), and to determine the individual optical scan wavelength ranges (Δλΐ, Δλ2, Δλ3) for each of the outer cores (CI, C2, C3) in response to the calculated measure of optical length ratios.
6. System according to claim 1, wherein at least some of a), b), and c) are implemented in hardware.
7. System according to claim 1, wherein at least some of a), b), and c) are implemented in software.
8. System according to claim 7, wherein some of a), b), and c) are implemented in software, while other ones of a), b), and c) are implemented in hardware.
9. System according to claim 7, wherein all of a), b), and c) are implemented in software.
10. System according to claim 7, where the optical console system is (P, C) arranged to select a limited portion of sampled data points obtained for an outer core, so as to reduce an individual scan wavelength range for the outer core.
1 1. System according to claim 1 , wherein the optical console system (P, C) comprises a controllable laser light source system (LS) arranged to generate light for the optical interrogation of the plurality of optical fiber cores with a light wavelength sweep within an individually controllable optical scan wavelength range (CC, CI, C2, C3) for each of the plurality of optical fiber cores.
12. System according to claim 1, comprising an optical shape sensor (OSF) comprising a plurality of optical fiber cores (CC, CI, C2, C3) comprising optical shape sensing properties and being arranged between a proximal end and a distal, wherein the plurality of optical fiber cores (CC, CI, C2, C3) comprises a center core (CC) and a plurality of outer cores (CI, C2, C3), wherein the center core (CC) is arranged centrally to the plurality of outer cores (CI, C2, C3) being wound around the center core (CC).
13. System according to claim 12, wherein the optical shape sensor (OSF) is arranged to follow a shape of an elongated medical interventional device (T), and wherein the system is arranged to display an image (I) generated in response to said measure of three- dimensional shape of the optical shape sensor (OSF).
14. Method for calibrating an optical console system (P, C) arranged for optical interrogation of an optical shape sensor (OSF) comprising a plurality of optical fiber cores arranged between a proximal end and a distal end, wherein the plurality of optical fiber cores comprise optical shape sensing properties, wherein a center core is one of the plurality of optical fiber cores arranged centrally to a plurality of outer cores, the method comprising:
- performing optical calibration measurements (P CM) on the plurality of optical fiber cores with one common optical scan wavelength range for all of the optical fiber cores,
- calculating (C OL), for each of the optical fiber cores, a measure of its optical length in response to the calibration measurements, and
- determining (D I SWR) an individual optical scan wavelength range for each of the outer cores in response to the individual measures of their optical lengths relative to the measure of optical length of the center core, so as to compensate for optical length differences between the plurality of optical fiber cores.
15. A computer program product having instructions which when executed cause an optical console system (P, C) to perform a calibration comprising: - performing optical calibration measurements (P CM) on a plurality of optical fiber cores with one common optical scan wavelength range for all of the optical fiber cores,
- calculating (C OL), for each of the optical fiber cores, a measure of its optical length in response to the calibration measurements, and
- determining (D I SWR) an individual optical scan wavelength range for each outer cores in response to the individual measures of their optical lengths relative to the measure of optical length of a center core, so as to compensate for optical length differences between the plurality of optical fiber cores.
PCT/EP2015/063176 2014-06-18 2015-06-12 Multi-core alignment in optical shape sensing WO2015193191A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP14172900.4 2014-06-18
EP14172900 2014-06-18

Publications (1)

Publication Number Publication Date
WO2015193191A1 true WO2015193191A1 (en) 2015-12-23

Family

ID=50942607

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/EP2015/063176 WO2015193191A1 (en) 2014-06-18 2015-06-12 Multi-core alignment in optical shape sensing

Country Status (1)

Country Link
WO (1) WO2015193191A1 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017182535A1 (en) 2016-04-20 2017-10-26 Koninklijke Philips N.V. Methods and systems for optically connecting an optical fiber sensor to an optical shape sensing console
CN112703364A (en) * 2018-09-20 2021-04-23 皇家飞利浦有限公司 Optical shape sensing system and method
CN112703365A (en) * 2018-09-20 2021-04-23 皇家飞利浦有限公司 Fiber optic sensors for shape sensing, optical shape sensing devices, systems and methods
CN114279363A (en) * 2016-07-08 2022-04-05 直观外科手术操作公司 Computation of redundant bends in a secure multicore fiber

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120069347A1 (en) * 2010-09-17 2012-03-22 Luna Innovations Incorporated Compensating for non-ideal multi-core optical fiber structure
US20140112615A1 (en) * 2011-06-14 2014-04-24 Luna Innovations Incorporated Co-registration of cores in multicore optical fiber sensing systems
WO2014072845A1 (en) * 2012-11-09 2014-05-15 Koninklijke Philips N.V. Optical frequency domain reflectometry system with multiple fibers per detection chain

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120069347A1 (en) * 2010-09-17 2012-03-22 Luna Innovations Incorporated Compensating for non-ideal multi-core optical fiber structure
US20140112615A1 (en) * 2011-06-14 2014-04-24 Luna Innovations Incorporated Co-registration of cores in multicore optical fiber sensing systems
WO2014072845A1 (en) * 2012-11-09 2014-05-15 Koninklijke Philips N.V. Optical frequency domain reflectometry system with multiple fibers per detection chain

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017182535A1 (en) 2016-04-20 2017-10-26 Koninklijke Philips N.V. Methods and systems for optically connecting an optical fiber sensor to an optical shape sensing console
CN109073820A (en) * 2016-04-20 2018-12-21 皇家飞利浦有限公司 For fibre optical sensor to be optically connected to the method and system of optic shape sensing console
JP2019522228A (en) * 2016-04-20 2019-08-08 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. Method and system for optically connecting an optical fiber sensor to an optical shape sensing console
CN109073820B (en) * 2016-04-20 2020-07-24 皇家飞利浦有限公司 Method and system for optically connecting a fiber optic sensor to an optical shape sensing console
US11079217B2 (en) 2016-04-20 2021-08-03 Koninklijke Philips N.V. Methods and systems for optically connecting an optical fiber sensor to an optical shape sensing console
JP7090554B2 (en) 2016-04-20 2022-06-24 コーニンクレッカ フィリップス エヌ ヴェ Methods and systems for optically connecting fiber optic sensors to optical geometry sensing consoles
CN114279363A (en) * 2016-07-08 2022-04-05 直观外科手术操作公司 Computation of redundant bends in a secure multicore fiber
CN112703364A (en) * 2018-09-20 2021-04-23 皇家飞利浦有限公司 Optical shape sensing system and method
CN112703365A (en) * 2018-09-20 2021-04-23 皇家飞利浦有限公司 Fiber optic sensors for shape sensing, optical shape sensing devices, systems and methods
US11519760B2 (en) 2018-09-20 2022-12-06 Koninklijke Philips N.V. Optical shape sensing system and method
CN112703365B (en) * 2018-09-20 2024-04-09 皇家飞利浦有限公司 Optical fiber sensor for shape sensing, optical shape sensing apparatus, system and method
US12001049B2 (en) 2018-09-20 2024-06-04 Koninklijke Philips N.V. Optical fiber sensor for shape sensing, optical shape sensing device, system and method

Similar Documents

Publication Publication Date Title
US10551173B2 (en) Methods and apparatus to determine a twist parameter and/or a bend angle associated with a multi-core fiber
CN111854811B (en) Registering measured fiber optic interference data using reference fiber optic interference data
JP6270483B2 (en) 3D shape reconstruction for optical tracking of elongated devices
CN104126103B (en) The method and device of motion compensation in interference-type sensor-based system
US8531655B2 (en) Compensating for non-ideal multi-core optical fiber structure
US8400620B2 (en) Registration of an extended reference for parameter measurement in an optical sensing system
US20170160148A1 (en) Stress visualization device and mechanical property value visualization device
JP2018515799A5 (en)
EP2577222B1 (en) Registration of an extended reference for parameter measurement in an optical sensing system
WO2016029110A1 (en) Method and apparatus for multiple localized interferometric measurements
WO2015193191A1 (en) Multi-core alignment in optical shape sensing
WO2014204839A1 (en) Methods and apparatus segmented calibration of a sensing optical fiber
Duncan et al. Characterization of a fiber-optic shape and position sensor
US11815371B2 (en) Overlapping fiber gratings
CN109307482B (en) Temperature insensitive optical fiber shape reconstruction system
JP2019537726A (en) Method and apparatus for optimizing optical performance of an interferometer

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 15727691

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 15727691

Country of ref document: EP

Kind code of ref document: A1