Cartilage diagnostic device using OCT

Abstract

To put cartilage diagnosis using OCT to a practical use more easily.SOLUTION: A cartilage diagnosis device 1 includes: an optical unit 2 including an optical system using OCT; a probe 4 connected to the optical unit 2, whose tip is configured so as to be inserted into a joint cavity, and which includes a light permeable contact surface that comes in contact with the cartilage and can transmit deformation energy and an optical mechanism for guiding the light from the optical unit 2 to the cartilage; a control and calculation part 6 for controlling the drive of the optical mechanism, calculating a Doppler velocity by processing an optical interference signal output from the optical unit 2 according to imparting of the deformation energy to the cartilage, and calculating a diagnostic evaluation value on a cartilage tissue on the basis of the Doppler velocity; and a display device 8 for displaying the diagnostic evaluation value on the cartilage tissue.SELECTED DRAWING: Figure 1

Description

[Technical field]

[0001] The present invention relates to a device for diagnosing the degree of cartilage degeneration. [Background technology]

[0002] Cartilage plays an important role in absorbing load shock and improving joint gliding, but it is a tissue that has difficulty in self-healing due to the lack of blood circulation. Many elderly people develop osteoarthritis (OA) due to wear of the cartilage, and there is a need to establish diagnostic and treatment methods for this condition, especially in countries with an aging society. Cartilage tissue is composed of 80% water and 20% matrix (extracellular matrix), which contains collagen and proteoglycan. In particular, proteoglycan bound to collagen fibers is thought to be greatly involved in the excellent viscoelastic properties of cartilage, such as determining the flow properties inside the cartilage. OA develops due to the loss of the viscoelastic properties of the cartilage.

[0003] With the recent advancement of medical diagnostic technology, optical coherence tomography (hereinafter referred to as "OCT") has been developed. OCT allows for micro-tomographic visualization of the inside of biological tissues in a non-invasive and non-contact manner. In addition, the acquisition rate of two-dimensional OCT tomographic images is equal to or higher than the video rate, and has high time resolution. A method for tomographically visualizing the mechanical properties of cartilage using OCT has also been proposed (Patent Document 1). According to this method, a predetermined stress is applied to the cartilage to perform tomographic visualization of the degree of deformation of the cartilage tissue, making it possible to distinguish between normal cartilage and degenerated cartilage. [Prior art documents] [Patent documents]

[0004] [Patent Document 1] Patent No. 6623163 Summary of the Invention [Problem to be solved by the invention]

[0005] The technology described in Patent Document 1 is considered to be one method for diagnosing early-stage cartilage diseases using OCT, but the diagnostic images are basically displayed by post-processing after acquiring a large number of OCT tomographic images continuously over time. Therefore, there is room for improvement in terms of practical application from the viewpoint of realizing rapid diagnosis in the medical field. In addition, the technology described in Patent Document 1 requires that the cross-section for acquiring the OCT tomographic image be fixed for signal processing reasons, so that it is difficult to put it into practical use in the actual clinical field in terms of body movement and doctor's technique.

[0006] The present invention has been made in view of the above problems, and one of its objects is to make cartilage diagnosis using OCT more practical. [Means for solving the problem]

[0007] One aspect of the present invention is a cartilage diagnostic device for diagnosing articular cartilage. This device includes an optical unit including an optical system using optical coherence tomography, a probe connected to the optical unit and configured such that a tip portion can be inserted into an articular cavity and has a light-transmitting contact surface capable of contacting cartilage to transmit deformation energy and an optical mechanism for guiding light from the optical unit to the cartilage, a control and calculation unit that controls the drive of the optical mechanism and processes an optical interference signal output from the optical unit in response to the application of deformation energy to the cartilage to calculate a Doppler velocity and calculates a diagnostic evaluation value of the cartilage tissue based on the Doppler velocity, and a display device that displays the diagnostic evaluation value of the cartilage tissue. Effect of the Invention

[0008] According to the present invention, cartilage diagnosis using OCT can be made more practical. [Brief description of the drawings]

[0009] [Figure 1] 1 is a diagram illustrating a schematic configuration of a cartilage diagnostic device according to an embodiment. [Diagram 2] FIG. 2 is a diagram illustrating a schematic configuration of a probe. [Diagram 3] FIG. 1 is a diagram showing a schematic configuration of an experimental apparatus. [Figure 4] FIG. 13 is a diagram showing experimental results. [Diagram 5] FIG. 13 is a diagram showing the results of an experiment verifying the relationship between Doppler velocity and tissue fluid flow velocity and the degree of cartilage degeneration. [Figure 6] 13 is a flowchart showing the flow of cartilage diagnosis processing. [Figure 7] 13 is a diagram illustrating a schematic configuration of a probe according to a modified example. FIG. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0010] One embodiment of the present invention is a cartilage diagnostic device that uses OCT. This cartilage diagnostic device captures OCT tomographic images while applying deformation energy to cartilage, and provides images for cartilage diagnosis based on the displacement speed of cartilage tissue calculated from the OCT tomographic images. A user such as a doctor can diagnose the degree of cartilage degeneration by viewing the images.

[0011] This cartilage diagnostic device includes an optical unit including an OCT optical system, and a probe connected to the optical unit. The tip of this probe can be inserted into the joint cavity, and has a contact surface capable of contacting the cartilage to transmit deformation energy (load), and an optical mechanism for directing light from the optical unit to the cartilage. The tip of the probe may be formed of a syringe needle. In order not to damage the surface of the cartilage, a light-transmitting resin, glass, or the like may be disposed at the tip of the syringe needle. A flexible sheath may be coaxially disposed inside the syringe needle, and glass, or the like, may be disposed at the tip of the sheath.

[0012] This probe may be one that can be operated manually by a doctor, or one that is automatically driven as part of a device. The former is preferable from the viewpoint of practical use in a wide range of medical settings. Since the contact surface of the probe is optically transparent, the light from the optical unit passes through the contact surface and is irradiated onto the cartilage. The light reflected from the cartilage passes through the contact surface and is captured by the probe.

[0013] When the doctor lightly presses the tip of the probe against the patient's cartilage, deformation energy is transmitted. The degree of cartilage degeneration can be evaluated by measuring the response at this time using OCT. This is because, as will be described later, the displacement speed of the cartilage tissue (including the displacement speed of the matrix and the flow speed of tissue fluid) is expressed in different ways depending on the degree of degeneration.

[0014] The optical unit combines the object light reflected by the cartilage with the reference light, and sends the resulting optical interference signal to the control and calculation unit. The control and calculation unit controls the driving of the optical mechanism, while also processing the optical interference signal to calculate the Doppler velocity (described below), and calculates a diagnostic evaluation value for the cartilage tissue based on the Doppler velocity. The display device displays the diagnostic evaluation value on a screen.

[0015] The "diagnostic evaluation value" may be the Doppler velocity itself. This is because a correlation between the Doppler velocity and the degree of cartilage degeneration has been observed from the experimental results described below. By displaying the cross-sectional distribution of the Doppler velocity, the degree of cartilage degeneration can be intuitively grasped. This Doppler velocity is obtained as the displacement velocity of the cartilage tissue, which is the combination of the displacement velocity of the matrix (extracellular matrix) of the cartilage tissue and the velocity (flow velocity) of the tissue fluid (water).

[0016] On the other hand, as described below, it is believed that the velocity of tissue fluid has a large effect on the early stage of cartilage degeneration. The displacement velocity of the matrix can be calculated from OCT tomographic measurements. Therefore, the velocity of tissue fluid can be calculated based on the difference between the Doppler velocity and the displacement velocity of the matrix, and the tomographic distribution of the calculated velocity can be displayed as the "diagnostic evaluation value."

[0017] The display format of the "diagnostic evaluation value" can be set appropriately taking into consideration the user's visibility. For example, the numerical value itself may be displayed, or a numerical range (which numerical range the value belongs to) may be displayed. Specifically, the visibility of the Doppler velocity and interstitial fluid velocity may be changed, for example, by color-coding each certain velocity range.

[0018] These tomographic distributions may be displayed in a time series manner by displaying images obtained by scanning only the depth direction (z direction) of the cartilage (i.e., one-dimensional images in one axis direction) over time. Alternatively, two-dimensional or three-dimensional images obtained by scanning in the depth direction and perpendicular directions may be displayed. When performing real-time diagnosis in a medical setting, it is preferable to display one-dimensional images.

[0019] More specifically, the control and calculation unit may calculate the diagnostic evaluation value based on two lines of tomographic signals in the depth direction that are continuously acquired over time from the cartilage, and may update the screen of the display device every time the diagnostic evaluation value is calculated based on the two lines of tomographic signals.

[0020] The method of transmitting deformation energy (load) may be based on a compression test method in which a constant strain rate is applied to the measurement object (cartilage), or a stress relaxation method in which a constant strain is applied and then stopped. Alternatively, it may be based on a dynamic viscoelasticity method in which a dynamic strain is applied to the measurement object. Alternatively, it may be based on a creep method in which a constant stress is applied to the measurement object.

[0021] Furthermore, a cartilage diagnostic method may be constructed using the above technology, which may include the steps of transmitting a predetermined deformation energy (load) to the cartilage, calculating a Doppler velocity by processing an optical interference signal of the OCT corresponding to the application of the deformation energy (load application), and outputting a diagnostic evaluation value of the cartilage tissue based on the Doppler velocity.

[0022] DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Hereinafter, preferred embodiments of the present invention will be described in detail with reference to the accompanying drawings. [Example] FIG. 1 is a diagram showing a schematic configuration of a cartilage diagnostic device according to an embodiment. The cartilage diagnostic device of this embodiment displays a diagnostic evaluation value for evaluating the degree of degeneration of cartilage tissue based on an OCT tomographic image of articular cartilage.

[0023] The cartilage diagnostic device 1 includes an optical unit 2, a probe 4, a control and calculation unit 6, and a display device 8. The optical unit 2 includes a light source 10, an object arm 12, a reference arm 14, an optical mechanism 16, and a light detection device 18. The probe 4 is connected to the optical unit 2, and constitutes one end of the object arm 12.

[0024] The optical elements of the optical unit 2 are connected to each other by optical fibers. In the illustrated example, an optical system based on a Mach-Zehnder interferometer is shown, but a Michelson interferometer or other optical system can also be used. In this embodiment, a Time Domain OCT (TD-OCT) is used as the OCT, but a Swept Source OCT (SS-OCT), Spectral Domain OCT (SD-OCT), or other OCTs may also be used.

[0025] The light emitted from the light source 10 is split by a coupler 20 (beam splitter), one of which is guided to the object arm 12 to become object light, and the other is guided to the reference arm 14 to become reference light. The object light from the object arm 12 is guided to the probe 4 via a circulator 22 and irradiated onto the cartilage, which is the subject of diagnosis. This object light is reflected as backscattered light on the surface and cross section of the cartilage, taken into the probe 4, and guided to the coupler 24 via the circulator 22.

[0026] An optical fiber 32 extending from the circulator 22 is inserted into the main body 30 of the probe 4. A syringe needle 34 is attached to the tip of the main body 30, and the tip of the probe 4 can be inserted into a living body (the patient's knee K) via the syringe needle 34. The main body 30 is provided with an optical mechanism for directing light from the optical unit 2 to the cartilage of the knee K, and an indenter (contact surface) for contacting the cartilage to transmit a predetermined load (stress). The optical mechanism irradiates light toward the cartilage and directs the reflected light to the object arm 12.

[0027] On the other hand, the reference light from the reference arm 14 is guided to the optical mechanism 16 (functioning as a "second optical mechanism") via a circulator 26. This reference light is reflected by a reference mirror 28 of the optical mechanism 16, returns to the circulator 26, and is guided to the coupler 24. That is, the object light and the reference light are multiplexed (superimposed) by the coupler 24, and the interference light is detected by the photodetector 18. The photodetector 18 detects this as an optical interference signal (a signal indicating the intensity of the interference light). This optical interference signal is input to the control and calculation unit 6 via an A / D converter 36.

[0028] The control and calculation unit 6 controls the entire optical system and performs calculation processing for image output by OCT. Command signals from the control and calculation unit 6 are input to the light source 10, optical mechanism 16, etc. via a D / A converter (not shown). The control and calculation unit 6 processes the optical interference signal input to the photodetector 18 and obtains a tomographic image of the cartilage by OCT. Then, based on the tomographic image data, an evaluation value (referred to as a "diagnostic evaluation value") for diagnosing the degree of degeneration of the cartilage tissue by a method described below is calculated.

[0029] More details are as follows. The light source 10 is a broadband light source consisting of a super luminescent diode (hereinafter referred to as "SLD"). The light emitted from the light source 10 is split into object light and reference light by a coupler 20, and these are guided to an object arm 12 and a reference arm 14, respectively. Single mode fibers are used as optical fibers in the object arm 12 and the reference arm 14.

[0030] The optical mechanism 16 is a rapid scanning optical delay line (RSOD) mechanism, and constitutes the reference arm 14. This RSOD may be a single-pass type in which light is irradiated to a diffraction grating 42 (described later) once in a round trip, or may be a double-pass type in which light is irradiated to a diffraction grating 42 (described later) twice in a round trip. The optical mechanism 16 includes a collimator lens 40, a diffraction grating 42, a lens 44, and a reference mirror 28 (resonant mirror).

[0031] The light that has passed through the circulator 26 is guided to the diffraction grating 42 via the collimator lens 40. This light is split into wavelengths by the diffraction grating 42, and the phases of the wavelengths are adjusted, and then the light is focused at different positions on the reference mirror 28 by the lens 44. By rotating the reference mirror 28 at a small angle, high-speed optical path scanning and frequency modulation are possible in a state in which the phases of the wavelengths are adjusted (a state in which the Doppler modulation frequency can be detected). The reflected light from the reference mirror 28 returns to the circulator 26 as reference light and is guided to the coupler 24. It is then superimposed with the object light and sent to the photodetector 18 as interference light. In a modified example, the light that has passed through the diffraction grating 42 may be focused on the reference mirror 28 by a curved mirror instead of the lens 44.

[0032] The photodetector 18 includes a photodetector 50, a filter 52, and an amplifier 54. The interference light obtained by passing through the coupler 24 is detected as an optical interference signal by the photodetector 50. This optical interference signal has noise removed or reduced by the filter 52, and is input to the control and calculation unit 6 via the amplifier 54 and the A / D converter 36.

[0033] The control and calculation unit 6 includes a CPU, a ROM, a RAM, a hard disk, etc. Using these hardware and software, the control and calculation unit 6 controls the entire optical system and performs calculation processing for image output by OCT. The control and calculation unit 6 controls the drive of the optical mechanism 16, processes the optical interference signal output from the photodetector 18 based on the drive, and obtains a tomographic image of the cartilage by OCT. Then, based on the tomographic image data, the control and calculation unit 6 calculates a diagnostic evaluation value of the cartilage tissue by a method described below.

[0034] The display device 8 is, for example, a liquid crystal display, and displays the diagnostic evaluation value of the cartilage tissue calculated by the control and calculation unit 6 on a screen in a tomographically visualized form.

[0035] FIG. 2 is a diagram illustrating a schematic configuration of the probe 4. As shown in FIG. The probe 4 has a syringe needle 34 attached as a puncture part to the tip of a stepped cylindrical body 30. The lower half of the body 30 is a narrow-diameter part 35 whose outer diameter is slightly larger than that of the syringe needle 34, improving the visibility of the syringe needle 34 when the doctor operates the probe 4.

[0036] An indenter 64 capable of contacting cartilage is coaxially disposed at the tip of the syringe needle 34. The indenter 64 is a light-transmitting member, and its tip surface serves as an abutment surface 65 capable of contacting the surface of the cartilage J. In this embodiment, the indenter 64 is made of glass, but plastic or other materials may be selected as long as the material is strong enough not to break or deform when pressed against the cartilage and is capable of transmitting object light.

[0037] A collimator lens 60 is disposed in the upper half of the main body 30, and an achromatic lens 62 (objective lens) is disposed at the boundary between the small diameter portion 35 and the syringe needle 34. The optical fiber 32 penetrates the wall of the main body 30 and is connected to the collimator lens 60. The collimator lens 60, the achromatic lens 62, and the indenter 64 are arranged on the same straight line, and constitute a "first optical mechanism."

[0038] The light guided to the probe 4 by the optical fiber 32 is emitted via the collimator lens 60, the achromatic lens 62, and the indenter 64, and is irradiated onto the cartilage J. As a result, the light reflected on the surface or inside of the cartilage J is taken in by the probe 4 and guided to the object arm 12 of the optical unit 2 via the optical fiber 32.

[0039] The calculation processing method up to presenting the diagnostic evaluation value of the cartilage tissue will be described below. As described above, the object light (light reflected from the cartilage J) that has passed through the object arm 12 and the reference light that has passed through the reference arm 14 are combined and detected as a low-coherence optical interference signal by the optical detection device 18. The control and calculation unit 6 acquires this optical interference signal as a tomographic image of the cartilage J based on the intensity of the interference light.

[0040] The coherence length l, which is the resolution in the optical axis direction (depth direction) of OCT c is determined by the autocorrelation function of the light source, where the coherence length l c is the half-width at half maximum of the envelope of the autocorrelation function, and can be expressed by the following equation (1).

number

[0041] On the other hand, the resolution in the direction perpendicular to the optical axis (beam scanning direction) is set to 1 / 2 of the beam spot diameter D based on the focusing performance of the focusing lens. The beam spot diameter ΔΩ can be expressed by the following formula (2).

number

[0042] (Doppler velocity detection) Optical interference signal intensity I obtained by OCT det(z) is obtained by the convolution integral shown in the following equation (3).

number

[0043] Optical interference signal intensity I det A fast Fourier transform is applied to (z), and the Hilbert transform of the following equation (4) is applied to obtain a spectrum S^(f) with a phase delay of π / 2 in the frequency domain. An inverse Fourier transform is applied to this spectrum S^(f) to obtain an analytic signal Γ(t)=s(t)+is^(t). The square root of the sum of the squares of the real part s(t) and the imaginary part s^(t) of this analytic signal (instantaneous amplitude) gives the signal intensity (signal envelope) of the OCT tomographic image.

[0044] By using the Hilbert transform, it is possible to obtain not only the instantaneous amplitude but also the instantaneous phase. The instantaneous frequency can also be calculated from the time derivative of the instantaneous phase. For this reason, it is particularly effective when processing high-frequency modulated signals as in this embodiment. This improves the detection resolution of the Doppler modulation frequency, enabling highly accurate and real-time detection to be achieved.

number

[0045] Here, the Doppler modulation frequency f generated by the internal fluctuation of the measurement target (cartilage) is dIn order to obtain in real time, the adjacent autocorrelation method is used. That is, for the optical interference signal (RF interference signal) acquired by the depth z direction scanning (A scan), the adjacent j, j+1 analytic signals in the x direction scanning (or the j, j+1 analytic signals with the x direction fixed) are respectively calculated as Γ j ,Γ j+1 In this case, Γ j+1 is expressed by the following equation (5).

number

[0046] Since the scanning (rotational drive) of the reference mirror 28 by the RSOD is periodic, its Doppler modulation frequency f RSOD is assumed to be constant. Then, the analytic signal Γ j and Γ j+1 The phase difference φ derived by the adjacent autocorrelation with i,j Using the Doppler modulation frequency f d can be calculated from the following formula (6). Here, i represents the number of sampling points of the optical interference signal in the depth z direction, and j represents the number of lines of the detected optical interference signal. i,j * is Gamma i,j The following equation (6) utilizes the fact that instantaneous phase can be detected by the Hilbert transform. Furthermore, the phase difference φ for N lines of A-scan (or N times of A-scan) i,j By performing ensemble averaging on and pixel averaging (averaging per pixel) in the z direction, the Doppler modulation frequency f d The detection accuracy can be improved. In addition, the time t corresponds to the number of sampling points i of the optical interference signal in the depth z direction.

number

[0047] This Doppler modulation frequency f d By calculating the following equation (7) using d The time t corresponds to the number of sampling points i of the optical interference signal in the depth z direction, and therefore represents the Doppler velocity at the depth z position. The Doppler velocity v detected for the cartilage is d is obtained as the combined value of the displacement rate of the cartilage tissue matrix and the flow rate of the tissue fluid.

number

[0048] (Stabilization of Doppler modulation detection ability) In this embodiment, instead of the autocorrelation for each sampling point shown in the above equation (6) (hereinafter also referred to as "point autocorrelation"), an integral-type autocorrelation (autocorrelation based on correlation integral: hereinafter also referred to as "autocorrelation integral") shown in the following equation (8) is calculated.

[0049] Here, when calculating the autocorrelation of two lines of A-scan data that are adjacent in time, the integral interval (-T w / 2≦t≦T w / 2) is set to reduce the fluctuation in the calculation of the phase difference. w represents the window width. Then, in that integration interval, the Doppler modulation frequency f d The degree of noise reduction can be adjusted by setting this integral interval.

number

[0050] Phase difference φ i,jWhen calculating, the ensemble average of N lines of A-scan (N A-scans) is calculated using the following formula (9), thereby making it possible to suppress system noise that fluctuates like shot noise.

number

[0051] (Detection of tissue displacement velocity and interstitial fluid flow velocity) The square root of the sum of the squares of the real part s(t) and the imaginary part s^(t) of the analytic signal Γ(t) is the amplitude A(t) = A i,j Gives.

number

[0052] Therefore, in order to detect the amount of movement (displacement speed) of the cartilage tissue, the correlation coefficient R of the instantaneous amplitude shown in the following formula (11) is used. A i+Δi,j The number of sampling points i, which represents the depth z position of the tissue morphology distribution, is changed by Δi to obtain R A i+Δi,j Calculate.

number

[0053] This correlation coefficient indicates the similarity of the depth pattern in the tissue morphology distribution. Therefore, as shown in the following formula (12), the correlation coefficient R A i+Δi,jThe position i+Δi of the sampling point where is maximum represents the tissue movement in the depth z-axis direction of the tissue morphology distribution.

number

[0054] Therefore, as shown in the following equation (13), the displacement velocity v of the matrix of the cartilage tissue is calculated by taking into account the distance Δz between the sampling points. t Detect.

number

[0055] As described above, the Doppler velocity in the above formula (7) is considered to include the displacement velocity of the matrix in the above formula (13). The displacement velocity of the matrix and the flow velocity of the interstitial fluid are calculated at each time and each position, and are detected as a time-space distribution. Therefore, the velocity v of the liquid component such as the interstitial fluid is calculated as w is given by the difference between the two, as shown in the following equation (14).

number

[0056] The velocity of this liquid component v w is considered to be the tissue fluid flow rate component of tissue permeation. In the living body, it represents the tissue permeation of plasma components leaking from capillaries located just below the epidermis, and in cultured regenerated tissue, it represents the permeation rate of culture fluid, and is considered to be a parameter for evaluating the activity of cell tissue. In the case of cartilage tissue, it represents the flow rate of the tissue fluid and is an index for evaluating the degree of cartilage degeneration.

[0057] In the above method, the movement amount (displacement speed) of the tissue is detected by applying the correlation method to the signal of the amplitude A(t) shown in the above formula (10). det (z)=I i,j Alternatively, the analytic signal Γ of the RF interference signal may be used. i,j For example, the analytic signal Γi,j When using the above equation, the number of sampling points i of the optical interference signal in the z direction is changed by Δi to obtain the correlation coefficient R Γ i+Δi,j Calculate.

number

[0058] This correlation represents the similarity of the depth patterns of low coherence interference signals (RF interference signals). The correlation coefficient R Γ i+Δi,j The position i+Δi of the sampling point where is maximum represents the tissue movement in the z-direction of the j+1th RF interference signal. Therefore, the tissue displacement velocity and the tissue fluid velocity can be detected based on the above equations (12) and (13).

[0059] Thus, the RF interference signal I det (z)=I i,j and the analytic signal Γ i,j When using , the correlation is calculated using a signal that also contains information on tissue fluid flow, so it is possible that noise is being used. However, since the amount of tissue movement is expected to be greater than the amount of tissue fluid flow, the analytical signal Γ is an important factor in evaluating the degree of cartilage degeneration. i,j The use of these methods is also considered to be useful.

[0060] Next, the effectiveness of the above-mentioned technique will be described. As described above, in this embodiment, the Doppler velocity of the cartilage tissue is calculated using OCT, and the degree of degeneration of the cartilage tissue can be evaluated based on the Doppler velocity. Here, the results of an experiment conducted to verify that the cartilage diagnostic device 1 can visually display the difference between normal cartilage and degenerated cartilage will be described.

[0061] In a normal diagnosis, a doctor holds the probe 4 and inserts its tip into the knee joint of a patient. Then, the tip surface (contact surface 65) of the indenter 64 is brought into contact with the surface of the cartilage J and pressed slightly. The amount of pressing is expected to be, for example, about 10 to 20 μm. In this experimental device, this process is mechanically reproduced in place of the doctor.

[0062] FIG. 3 is a diagram showing a schematic configuration of the experimental apparatus. This experimental device is configured by connecting a compression tester 70 to the optical unit 2 shown in Fig. 1. The compression tester 70 is equipped with a mounting table 72 for mounting a test piece W, a support mechanism 74 for supporting a probe 4 in the vertical direction above the mounting table 72, and a load cell 76 capable of applying a downward pressing load (i.e., toward the test piece W) to the rear end of the probe 4. The load cell 76 is driven by the control and calculation unit 6. The pressing force of the probe 4 can also be detected by the load cell 76.

[0063] In this experiment, a cartilage tissue sample from the lateral femoral condyle of a porcine knee joint was collected in a cylindrical shape with a diameter of 7 mm, and this was used as a normal cartilage test piece W. As a comparative example, a simulated denatured cartilage (OA test piece in which the degree of denaturation was artificially controlled to simulate cartilage) was prepared by enzymatically treating this test piece W. For this enzymatic treatment, collagenase was used, which breaks down collagen fibers that maintain the elasticity of articular cartilage. That is, normal cartilage was placed in a solution in which collagenase was dissolved in phosphate buffered saline at pH 7.4, and the solution was left at 37°C for 2 hours to obtain a simulated denatured cartilage.

[0064] Figure 4 shows the experimental results. Figure 4(A) shows the experimental results for normal cartilage, and Figure 4(B) shows the experimental results for simulated degenerated cartilage. The top of each figure shows the time change in OCT cross-sectional images, and the middle shows the time change in the displacement velocity (Doppler velocity) of the cartilage tissue. The bottom of each figure shows the relationship between the pressure acting on the test piece W and the stroke of the probe 4, along with the time change.

[0065] The experimental conditions were as follows: cAn LSD broadband light source with a wavelength of 1317.7 nm and a full width at half maximum Δλ of 102.0 nm was used. The compression amount by the compression tester 70 was 100 μm and the compression speed was 100 μm / sec. The scanning frequency of the reference mirror 28 was about 4.0 kHz, the depth resolution was 7.5 μm, and the depth tomographic scanning range was 1130 μm. The sampling frequency by the OCT was 20 MHz.

[0066] From the bottom of Fig. 4, the stroke of the probe 4 is zero until time t = 0.7 sec, and the compression tester 70 is stopped. Driving of the probe 4 by the compression tester 70 begins at time t = 0.7 sec, and the probe 4 is displaced at a constant speed until approximately time 2.8 sec. During this time, the indenter 64 approaches the test piece W, and comes into contact with the surface of the test piece W at approximately time 1.2 sec. The position of the probe 4 is fixed at approximately time t = 2.8 sec. As a result, stress relaxation occurs in the test piece W from that point on.

[0067] In the upper part of Figure 4, the area approximately 200 μm from the cartilage surface corresponds to the superficial layer, and below that corresponds to the middle layer. In the OCT cross-sectional image (tissue morphology cross-sectional distribution), it can be seen that a high-intensity signal moves diagonally upward from the time when compression begins in the cartilage and a load response can be confirmed until the indenter 64 comes to rest. This indicates that the cartilage tissue is displaced in the direction approaching the indenter 64, and the amount of displacement. From the upper part of Figure 4, it can be qualitatively confirmed that the amount of displacement in the simulated degenerated cartilage (t = 1.3-2.9 sec) is greater than the amount of displacement in normal cartilage (t = 1.2-2.8 sec).

[0068] The time course of the Doppler velocity is displayed in color in the middle of Figure 4. The Doppler velocity during compression of the cartilage shows that the cartilage tissue or tissue fluid moves toward the cartilage surface in both normal cartilage and the simulated degenerated cartilage. This is easily visualized on the screen.

[0069] This figure shows that the Doppler velocity in the middle layer of the cartilage is approximately 20 μm / sec in normal cartilage, while in the simulated degenerated cartilage it is more than twice as high at approximately 50 μm / sec. Since the Doppler velocity has a spatial distribution that increases in the depth z direction, it can also be confirmed that greater strain is occurring in the surface layer of the cartilage compared to the middle layer. This is thought to be due to the fact that the elastic modulus of the surface layer is lower than that of the middle layer, due to the orientation characteristics of the collagen fibers. Furthermore, in the simulated degenerated cartilage, the elastic and viscous properties are further reduced due to the destruction of collagen fibers and the detachment of proteoglycans, which is thought to be why the large displacement velocity was detected.

[0070] As shown in the bottom panel of Figure 4, stress relaxation began after time t = 2.8 sec for normal cartilage and after time t = 2.9 sec for simulated degenerated cartilage. As shown in the middle panel of Figure 4, the Doppler velocity approximately 0.5 sec after the start of stress relaxation was detected to be 5 μm / sec or less for normal cartilage, and movement quickly stopped. In contrast, the Doppler velocity for the simulated degenerated cartilage was as high as about 20 μm / sec. This is also thought to be due to the increased tissue permeability in the simulated degenerated cartilage, i.e., the decreased viscosity characteristics.

[0071] Figure 5 shows the results of an experiment verifying the relationship between Doppler velocity and tissue fluid flow velocity and the degree of cartilage degeneration. Figure 5(A) shows the relationship between Doppler velocity and the degree of cartilage degeneration, and Figure 5(B) shows the relationship between tissue fluid flow velocity and the degree of cartilage degeneration. The vertical axis shows velocity during compression in Figure 5(A) and during stress relaxation in Figure 5(B). On the horizontal axis, "Con" shows normal cartilage, and "2H" shows simulated degenerated cartilage.

[0072] As shown in Figure 5(A), a significant difference was observed between normal cartilage and the simulated degenerated cartilage in terms of Doppler velocity. As shown in Figure 5(B), a significant difference was also observed between normal cartilage and the simulated degenerated cartilage in terms of tissue fluid flow velocity.

[0073] From the above experiments, it was found that there was a correlation between the distribution of Doppler velocity and tissue fluid flow velocity and the mechanical behavior of cartilage, and that they were significantly different between normal cartilage and simulated degenerated cartilage. In other words, it was suggested that the diagnosis of degenerated cartilage using the OCT of this embodiment is effective. For this reason, for example, the distribution of Doppler velocity shown in the middle of Figure 4 can be displayed as a diagnostic evaluation value of cartilage. In addition, although not shown in this experiment, the tissue fluid flow velocity obtained from the difference between the Doppler velocity and the displacement velocity of the matrix can also be displayed as a diagnostic evaluation value. In the early stage of cartilage degeneration, the water retention ability (the function of keeping tissue fluid stable) is lost, which may lead to the discovery of a tendency to degeneration.

[0074] Next, a specific flow of processing executed by the control and calculation unit 6 will be described. 6 is a flowchart showing the flow of the cartilage diagnosis process. Note that this cartilage diagnosis process is performed by a doctor or the like with the tip of the probe 4 inserted into the knee joint of a patient, and is repeatedly executed at a predetermined calculation cycle when the drive of the optical unit 2 is started.

[0075] The control and calculation unit 6 acquires the OCT interference signal (S10) and performs a Fourier transform (FFT) (S12). The control and calculation unit 6 then performs bandpass filtering based on the carrier frequency to improve the signal-to-noise ratio (S14) and performs a Hilbert transform (S16). The analytic signal obtained by the Hilbert transform is used to perform adjacent autocorrelation processing to obtain a phase difference (S18) and a Doppler modulation frequency (S20).

[0076] The control and calculation unit 6 calculates the Doppler velocity v based on the above equation (7) using the Doppler modulation frequency. d On the other hand, (S22) calculates the displacement velocity v of the matrix of the cartilage tissue based on the above formula (12). t Then, the velocity v of the tissue fluid is calculated based on the above formula (13) (S24). w (S26). The control and calculation unit 6 calculates the Doppler velocity v d and the tissue fluid velocity v wEither one or both of the above are displayed on the screen of the display device 8 as a diagnostic evaluation value (S28).

[0077] As described above, according to this embodiment, the diagnostic evaluation value of the cartilage to be diagnosed is visualized in a cross-sectional view based on the finding that there is a correspondence between the Doppler velocity calculated from the OCT cross-sectional image of the cartilage and the degree of degeneration of the cartilage. Therefore, a doctor or the like can easily perform a cartilage diagnosis by checking the cross-sectional visualized image.

[0078] In particular, since it is sufficient to obtain an OCT tomographic image (one-dimensional tomographic distribution) in the depth z direction, diagnostic images can be displayed in a short time, enabling real-time diagnosis in the medical field. Even if the OCT tomographic image (two-dimensional tomographic distribution) in the depth z direction is used, it is also possible to display a two-dimensional tomographic color map of Doppler velocity in real time, since only the correlation of one axis of the A-scan is used. In addition, there is no need to fix the position of the probe 4 relative to the cartilage when taking an image. Therefore, the patient's body movement and the doctor's technique are less likely to affect the diagnosis. Therefore, cartilage diagnosis using OCT can be made more practical in actual clinical settings.

[0079] Although the preferred embodiment of the present invention has been described above, it goes without saying that the present invention is not limited to the specific embodiment, and various modifications are possible within the scope of the technical concept of the present invention.

[0080] [Variations] FIG. 7 is a diagram illustrating a schematic configuration of a probe according to a modified example. The probe 104 of this modified example includes an optical mechanism 80. The optical mechanism 80 constitutes the object arm 12 and includes a mechanism for directing light from the light source 10 to the measurement object (cartilage J) for scanning, and a drive unit for driving the mechanism. The optical mechanism 80 includes a collimator lens 60, a two-axis galvanometer mirror 82 of an optical unit, and an achromatic lens 62. The light that has passed through the circulator 22 is directed to the galvanometer mirror 82 via the collimator lens 60, and is irradiated onto the measurement object while being scanned in the x-axis and y-axis directions. The reflected light from the measurement object is captured by the probe 104.

[0081] Providing such an optical mechanism 80 in the probe 104 enables OCT measurement of two-dimensional tomographic distribution or three-dimensional tomographic distribution, thereby making it possible to display the diagnostic evaluation value of the cartilage in two or three dimensions.

[0082] [Other variations] Although not mentioned in the above embodiment, a sensor such as a load cell may be attached to the probe 4 itself, making it possible to measure the load applied to the cartilage. This makes it possible to display the degree of stress relaxation, for example, as shown in the lower part of Fig. 4, and by reading this together with the displacement velocity image shown in the middle part of Fig. 4, it is possible to display the diagnostic evaluation value of the cartilage tissue from various angles.

[0083] In the above embodiment, as shown in Fig. 2, it is assumed that the pressing force by the probe 4 is applied manually by the user (doctor). In a modified example, a load device may be provided on the probe 4 so that a predetermined deformation energy (load) can be applied to the cartilage. As such a load device, a loading mechanism that applies stress by contact with a piezoelectric element or the like may be used. Alternatively, a load device that applies a load (excitation force) to the cartilage in a non-contact manner by ultrasound (sound pressure), photoacoustic waves, electromagnetic waves, or the like may be used.

[0084] The present invention is not limited to the above-mentioned embodiment and modifications, and the components can be modified without departing from the spirit of the present invention. Various inventions can be formed by appropriately combining multiple components disclosed in the above-mentioned embodiment and modifications. Also, some components can be deleted from all the components shown in the above-mentioned embodiment and modifications. [Explanation of symbols]

[0085] 1 cartilage diagnostic device, 2 optical unit, 4 probe, 6 control and calculation unit, 8 display device, 10 light source, 12 object arm, 14 reference arm, 16 optical mechanism, 18 photodetector, 20 coupler, 22 circulator, 24 coupler, 26 circulator, 28 reference mirror, 30 main body, 32 optical fiber, 34 syringe needle, 35 thin-diameter portion, 40 collimator lens, 42 diffraction grating, 44 lens, 50 photodetector, 52 filter, 54 amplifier, 60 collimator lens, 62 achromatic lens, 64 indenter, 65 contact surface, 70 compression tester, 72 mounting table, 74 support mechanism, 76 load cell, 80 optical mechanism, 82 galvanometer mirror, 104 probe, J cartilage, K knee, W test piece.

Claims

1. A cartilage diagnostic device for diagnosing articular cartilage, comprising: an optical unit including an optical system using optical coherence tomography; a probe connected to the optical unit, the probe having a tip portion configured to be insertable into a joint cavity, the probe having a light-transmitting contact surface capable of contacting cartilage to transmit deformation energy, and an optical mechanism for guiding light from the optical unit to the cartilage; a control and calculation unit that controls the driving of the optical mechanism, processes an optical interference signal output from the optical unit in response to the application of deformation energy to the cartilage, calculates a Doppler velocity, and calculates a diagnostic evaluation value of the cartilage tissue based on the Doppler velocity; A display device that displays the diagnostic evaluation value of the cartilage tissue; Equipped with The control calculation unit is calculating a displacement velocity of the matrix of the cartilage by processing an optical interference signal output from the optical unit in response to the application of deformation energy to the cartilage; A cartilage diagnostic device characterized by calculating a velocity of tissue fluid of the cartilage based on a difference between the Doppler velocity and the displacement velocity of the matrix, and calculating a diagnostic evaluation value of the cartilage tissue based on the velocity of the tissue fluid.

2. A first optical mechanism provided on an object arm passing through the cartilage and functioning as the optical mechanism; a second optical mechanism provided on the reference arm that does not pass through the cartilage and that guides light from a light source to a reference mirror and reflects the light; a light detection device that detects interference light obtained by superimposing an object light reflected by the cartilage and a reference light reflected by the reference mirror; Equipped with The cartilage diagnostic device according to claim 1 , wherein the control and calculation unit performs frequency analysis on the optical interference signal input from the photodetector, and calculates the Doppler velocity using a Doppler modulation frequency.

Images