METHOD AND SYSTEM OF QUANTITATIVE INTRAVASCULAR OPTICAL COHERENCE TOMOGRAPHY

Abstract

The present disclosure relates to a method and a system of quantitative intravascular optical coherence tomography. The method includes: acquiring and graying an intravascular optical coherence tomography cross-sectional image, and converting it into a polar coordinate view; calculating a confocal function value and a sensitivity function value according to parameters of the imaging system; establishing a model of an intensity of a backscattered signal being attenuated depending on the detection depth, according to the confocal function value and the sensitivity function value to obtain a theoretical value of the backscattered signal; determining a measured value of the backscattered signal; constructing an error function according to the theoretical value and the corresponding measured value; minimizing the error function to determine a light attenuation coefficient corresponding to each polar coordinate; and determining a spatial distribution diagram of the light attenuation coefficient according to the light attenuation coefficient corresponding to each polar coordinate.

Claims

1. A method of quantitative intravascular optical coherence tomography, comprising: acquiring an intravascular optical coherence tomography cross-sectional image output by an imaging system; wherein, the intravascular optical coherence tomography cross-sectional image is a transverse view; a center of the intravascular optical coherence tomography cross-sectional image corresponds to a center of an imaging catheter, and a vascular lumen, a vascular wall intima/media and a vascular wall adventitia are around the center of the image in sequence; graying the intravascular optical coherence tomography cross-sectional image, and converting the grayed intravascular optical coherence tomography cross-sectional image into a polar coordinate view; respectively calculating a confocal function value and a sensitivity function value of each polar coordinate in the polar coordinate view according to parameters of the imaging system; wherein, the parameters of the imaging system comprise: a center wavelength of an incident light, a resolution of a spectrometer, and a wavelength interval; establishing a model of an intensity of a backscattered signal being attenuated depending on a detection depth, according to the confocal function value and the sensitivity function value of each polar coordinate to obtain a theoretical value of the backscattered signal; determining a measured value of the backscattered signal according to a gray value of each pixel in the grayed intravascular optical coherence tomography cross-sectional image; constructing an error function according to the theoretical value of the backscattered signal and the measured value of the backscattered signal corresponding to each polar coordinate; minimizing the error function to determine a light attenuation coefficient corresponding to each polar coordinate in the polar coordinate view; and determining a spatial distribution diagram of the light attenuation coefficient according to the light attenuation coefficient corresponding to each polar coordinate in the polar coordinate view.

2. The method according to claim 1, wherein graying the intravascular optical coherence tomography cross-sectional image, and converting the grayed intravascular optical coherence tomography cross-sectional image into the polar coordinate view comprises: graying the intravascular optical coherence tomography cross-sectional image to determine the grayed intravascular optical coherence tomography cross-sectional image; establishing an XOY rectangular coordinate system with a center of the imaging catheter as a coordinate origin O, the horizontal right direction as a positive direction of an X-axis, and an upward vertical direction as a positive direction of a Y-axis, according to the grayed intravascular optical coherence tomography cross-sectional image; establishing a θ-l polar coordinate system with O as a coordinate origin according to the XOY rectangular coordinate system, and converting the grayed intravascular optical coherence tomography cross-sectional image into a polar coordinate view; wherein θ is a polar angle, l is a polar diameter, l-axis is a radial thickness of a multilayered vascular wall tissue, and a surface of the vascular wall is parallel to θ-axis and perpendicular to the l-axis.

3. The method according to claim 1, wherein respectively calculating the confocal function value and the sensitivity function value of each polar coordinate in the polar coordinate view according to the parameters of the imaging system comprises: determining the confocal function value at a polar coordinate (θ.sub.i, m) by using a formula: $h\left({\theta }_i,m\right)={\int }_{m\Delta }^{\left(m+1\right)\Delta }{\left\{{\left[\frac{{\lambda }_0\left(z-z_0\right)}{{\alpha }_{R}n\pi w_0^2}\right]}^2+1\right\}}^{-1}\mathrm{dz};$ determining the sensitivity function value at the polar coordinate (θ.sub.i, m) by using a formula: $f\left({\theta }_i,m\right)={\int }_{m\Delta }^{\left(m+1\right)\Delta }\left[{\left(\frac{\sin \zeta }{\zeta }\right)}^2\exp \left(-\frac{\delta \lambda /\Delta \lambda }{2\ln 2}{\zeta }^2\right)\right]\mathrm{dz};$ wherein, h(θ.sub.i, m) is the confocal function value at the polar coordinate (θ.sub.i, m), f(θ.sub.i, m) is the sensitivity function value at the polar coordinate (θ.sub.i, m) θ.sub.i∈[0, 2π], m=0, 1, . . . , M, M is a radius (unit: pixel) of a frame of the intravascular optical coherence tomography cross-sectional image, Δ is a size (unit: mm) of the pixel in the image, α.sub.R is a derivative parameter of a secondary laser beam formed by an interaction of an incident beam with scatterers in a medium (α.sub.R=1 for a specular reflection; and α.sub.R=2 for a diffuse reflection), n is a refractive index of a tissue to light, z.sub.0 is a focal plane position (z.sub.0 is negative when the focal plane is above the tissue, z.sub.0 is positive when the focal plane is inside the tissue), w.sub.0 is a light spot radius after the incident light is focused, λ.sub.0 is the center wavelength of the incident light, δλ is the resolution of the spectrometer, Δλ is the wavelength interval, λ.sub.0, δλ and Δλ are the parameters of the imaging system, ζ is a normalized value of a depth z with respect to a maximum detection depth, wherein, $\zeta =2\pi z\frac{\Delta \lambda }{{\lambda }_0^2}.$

4. The method according to claim 3, wherein establishing the model of the intensity of the backscattered signal being attenuated depending on the detection depth, according to the confocal function value and the sensitivity function value of each polar coordinate to obtain the theoretical value of the backscattered signal comprises: determining an irradiance of the incident light at the polar coordinate (θ.sub.i, m) by using a formula: $L_{\mathrm{fh}}\left({\theta }_i,m\right)=L_{0}f\left({\theta }_i,m\right)h\left({\theta }_i,m\right)\exp \left[-\Delta \underset{i=0}{\overset{m}{.Math.}}{\mu }_i\left({\theta }_i,l\right)\right];$ determining an irradiance of the backscattered light at the polar coordinate (θ.sub.i, m) by using a formula: $\begin{array}{c}B\left({\theta }_i,m\right)={\mathrm{\alpha \mu }}_i\left({\theta }_i,m\right)L_{\mathrm{fh}}\left({\theta }_i,m\right)\\ ={\mathrm{\alpha \mu }}_i\left({\theta }_i,m\right)L_{0}f\left({\theta }_i,m\right)h\left({\theta }_i,m\right)\exp \left[-\Delta \underset{i=0}{\overset{m}{.Math.}}{\mu }_i\left({\theta }_i,l\right)\right];\end{array}$ determining an irradiance of the backscattered light reaching a detector surface at the polar coordinate (θ.sub.i, m) by using a formula: $\begin{array}{c}A\left({\theta }_i,m\right)=B\left({\theta }_{i}m\right)\exp \left[-\Delta \underset{l=0}{\overset{m}{.Math.}}{\mu }_i\left({\theta }_i,l\right)\right]\\ ={\mathrm{\alpha \mu }}_i\left({\theta }_i,m\right)L_{0}f\left({\theta }_i,m\right)h\left({\theta }_i,m\right)\exp \left[-2\Delta \underset{i=0}{\overset{m}{.Math.}}{\mu }_i\left({\theta }_i,l\right)\right];\end{array}$ and determining the theoretical value of the backscattered signal at the polar coordinate (θ.sub.i, m) by using a formula: $\begin{array}{c}{I}_b\left({\theta }_i,m\right)=n\left({\theta }_i,m\right)+\beta A\left({\theta }_i,m\right)\\ =n\left({\theta }_i,m\right)+\beta {\mathrm{\alpha \mu }}_i\left({\theta }_i,m\right)L_{0}f\left({\theta }_i,m\right)\\ h\left({\theta }_i,m\right)\exp \left[-2\Delta \underset{i=0}{\overset{m}{.Math.}}{\mu }_i\left({\theta }_i,l\right)\right];\end{array}$ wherein, L.sub.fh(θ.sub.i, m) is the irradiance of the incident light at the polar coordinate (θ.sub.i, m) in an imaging plane, L.sub.0 is an incident light power (unit: mW or dBm), μ.sub.t(θ.sub.i, l) is a tissue light attenuation coefficient at the polar coordinate (θ.sub.i, l), B(θ.sub.i, m) is the irradiance of the backscattered light at the polar coordinate (θ.sub.i, m), α is a proportion of the backscattered light in the incident light, A(θ.sub.i, m) is the irradiance of the backscattered light reaching the detector surface at the polar coordinate (θ.sub.i, m), I.sub.b(θ.sub.i, m) is the theoretical value of the backscattered signal at the polar coordinate (θ.sub.i, m), n(θ.sub.i, m) is a noise power value introduced by measurement at the polar coordinate (θ.sub.i, m), and β is an A/D conversion factor of the imaging system.

5. The method according to claim 1, wherein determining the measured value of the backscattered signal according to the gray value of each pixel in the grayed intravascular optical coherence tomography cross-sectional image comprises: determining a normalized measured value of the backscattered signal at a polar coordinate (θ.sub.i, m) by using a formula: $I_m\left({\theta }_i,m\right)=\frac{g\left({\theta }_i,m\right)}{255};$ wherein, g(θ.sub.i, m) is a gray value at the polar coordinate (θ.sub.i, m) in the polar coordinate view of an intravascular optical coherence tomography image, and I.sub.m(θ.sub.i, m) is the normalized measured value of the backscattered signal at the polar coordinate (θ.sub.i, m).

6. The method according to claim 1, wherein constructing the error function according to the theoretical value of the backscattered signal and the measured value of the backscattered signal corresponding to each polar coordinate comprises: determining the error function by using a formula:
f(μ.sub.t(θ.sub.i,m))=Ī.sub.b(θ.sub.i,m)−I.sub.m(θ.sub.i,m); wherein, f(.Math.) is the error function, Ī.sub.b(θ.sub.i, m) is a normalized theoretical value of the backscattered signal at a polar coordinates (θ.sub.i, m).

7. The method according to claim 1, wherein minimizing the error function to determine the light attenuation coefficient corresponding to each polar coordinate in the polar coordinate view comprises: determining an estimated value of the light attenuation coefficient by using a formula: ${\hat{\mu }}_t\left({\theta }_i,m\right)=\underset{\mu \ge 0}{\arg \min }F\left(\mu \right)=\underset{\mu \ge 0}{\arg \min }{.Math.f\left(\mu \right).Math.}^2;$ performing Taylor expansion approximation on an objective function F(μ) near μ by using a formula:
F(μ+Δμ)=∥f(μ+Δμ)∥.sup.2≈∥f(μ)∥.sup.2+f′(μ)Δμ+½f″(μ)(Δμ).sup.2; converting an optimization problem of the an estimated value of the light attenuation coefficient into solving an optimal increment Δμ for each iteration by using a formula as follows, and adding a damping tam to punish Δμ: $\hspace{1em}\{\begin{array}{c}\Delta \hat{\mu }=\underset{\mathrm{\Delta \mu }}{\arg \min }m\left(\Delta \mu \right)\\ m\left(\Delta \mu \right)={.Math.f\left(\mu \right).Math.}^2+f^{\prime }\left(\mu \right)\Delta \mu +\frac{1}{2}{f}^{″}\left(\mu \right){\left(\Delta \mu \right)}^2+\frac{1}{2}{\gamma \left(\mathrm{\Delta \mu }\right)}^2\end{array};$ and solving a formula: $\hspace{1em}\{\begin{array}{c}\Delta \hat{\mu }=\underset{\mathrm{\Delta \mu }}{\arg \min }m\left(\Delta \mu \right)\\ m\left(\Delta \mu \right)={.Math.f\left(\mu \right).Math.}^2+f^{\prime }\left(\mu \right)\Delta \mu +\frac{1}{2}{f}^{″}\left(\mu \right){\left(\Delta \mu \right)}^2+\frac{1}{2}{\gamma \left(\mathrm{\Delta \mu }\right)}^2\end{array},$ to determine the light attenuation coefficient corresponding to each polar coordinate in the polar coordinate: wherein, F(μ) is the objective function, {circumflex over (μ)}.sub.t(θ.sub.i, m) is the estimated value of the light attenuation coefficient, f′(μ) and f″(μ) ae a first derivative and a second derivative of f(μ) respectively.

8. A system of quantitative intravascular optical coherence tomography, comprising: an acquisition module for intravascular optical coherence tomography cross-sectional image, configured to acquire an intravascular optical coherence tomography cross-sectional image output by an imaging system; wherein, the intravascular optical coherence tomography cross-sectional image is a transverse view; a center of the intravascular optical coherence tomography cross-sectional image corresponds to a center of an imaging catheter, and a vascular lumen, a vascular wall intima/media and a vascular wall adventitia are around the canter of the image in sequence; a determination module for a polar coordinate view, configured to gray the intravascular optical coherence tomography cross-sectional image, and convert the grayed intravascular optical coherence tomography cross-sectional image into a polar coordinate view; a determination module for a confocal function value and a sensitivity function value, configured to respectively calculate the confocal function value and the sensitivity function value of each polar coordinate in the polar coordinate view according to parameters of the imaging system wherein, the parameters of the imaging system comprise: a center wavelength of an incident light, a resolution of a spectrometer, and a wavelength interval; a determination module for a theoretical value of the backscattered signal, configured to establish a model of the intensity of the backscattered signal being attenuated depending on a detection depth, according to the confocal function value and the sensitivity function value of each polar coordinate to obtain the theoretical value of the backscattered signal; a determination module for a measured value of the backscattered signal, configured to determine the measured value of the backscattered signal according to a gray value of each pixel in the grayed intravascular optical coherence tomography cross-sectional image; an error function constructing module, configured to construct an error function according to the theoretical value of the backscattered signal and the measured value of the backscattered signal corresponding to each polar coordinate; a determination module for a light attenuation coefficient, configured to minimize the error function to determine the light attenuation coefficient corresponding to each polar coordinate in the polar coordinate view; and a determination module for a spatial distribution diagram of the light attenuation coefficient, configured to determine the spatial distribution diagram of the light attenuation coefficient according to the light attenuation coefficient corresponding to each polar coordinate in the polar coordinate view.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0052] In order to more clearly illustrate embodiments of the present disclosure or technical solutions in the conventional technology, accompanying drawings used in the embodiments will now be described briefly. It is obvious that the drawings in the following description are only some embodiments of the present disclosure, and that those skilled in the art can obtain other drawings from these drawings without involving any inventive effort.

[0053] FIG. 1 is a schematic flow chart of a method of quantitative intravascular optical coherence tomography provided by the present disclosure;

[0054] FIG. 2 is an intravascular optical coherence tomography grayscale image;

[0055] FIG. 3 is a schematic diagram of a cross section of a blood vessel in an intravascular coherence tomography image;

[0056] FIG. 4 is a schematic flow chart of an algorithm minimizing an error function; and

[0057] FIG. 5 is a schematic structural diagram of a system of quantitative intravascular optical coherence tomography provided by the present disclosure.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0058] In the following technical solutions in the embodiments of the present disclosure will be clearly and completely described with reference to the drawings in the embodiments of the present disclosure. Obviously, the described embodiments are only a part of the embodiments of the present disclosure, but not all the embodiments. Based on the embodiments of the present disclosure, all other embodiments obtained by a person of ordinary skill in the art without involving any inventive effort are within the scope of the present disclosure.

[0059] The present disclosure intends to provide a method and a system of quantitative intravascular optical coherence tomography, which can quantitatively estimate a spatial distribution of a light attenuation coefficient of a tissue on a cross section of a blood vessel with a standard cross-sectional image output by an imaging system, to improve a capability to identify different tissue components.

[0060] To further clarify the above objects, features and advantages of the present disclosure, a more particular description of the disclosure will be rendered by reference to the accompanying drawings and specific embodiments thereof.

[0061] FIG. 1 is a schematic flow chart of a method of quantitative intravascular optical coherence tomography provided by the present disclosure. As shown in FIG. 1, the method of quantitative intravascular optical coherence tomography provided by the present disclosure may include the following steps S101 to S108.

[0062] S101, an intravascular optical coherence tomography cross-sectional image output by an imaging system is acquired. Wherein, the intravascular optical coherence tomography cross-sectional image is a transverse view; a center of the intravascular optical coherence tomography cross-sectional image corresponds to a center of an imaging catheter, and a vascular lumen, a vascular wall intima/media and a vascular wall adventitia are around the cuter of the image in sequence, as shown in FIG. 2.

[0063] S102, the intravascular optical coherence tomography cross-sectional image is grayed, and the grayed intravascular optical coherence tomography cross-sectional image is converted into a polar coordinate view.

[0064] Step S102 specifically includes:

[0065] graying the intravascular optical coherence tomography cross-sectional image to determine the grayed intravascular optical coherence tomography cross-sectional image.

[0066] As shown in FIG. 3, an XOY rectangular coordinate system, with a center of the imaging catheter as a coordinate origin O, the horizontal right direction as a positive direction of an X-axis, and an upward vertical direction as a positive direction of a Y-axis, is established according to the grayed intravascular optical coherence tomography cross-sectional image.

[0067] A θ-l polar coordinate system with O as a coordinate origin is established according to the XOY rectangular coordinate system, and the grayed intravascular optical coherence tomography cross-sectional image is convened into a polar coordinate view. Wherein θ is a polar angle, l is a polar diameter, l-axis is a radial thickness of a multilayered vascular wall tissue, and a surface of the vascular wall is parallel to θ-axis and perpendicular to the l-axis.

[0068] S103, a confocal function value and a sensitivity function value of each polar coordinate in the polar coordinate view are respectively calculated according to parameters of the imaging system. Wherein, the parameters of the imaging system comprise: a center wavelength of an incident light, a resolution of a spectrometer, and a wavelength interval;

[0069] Step S103 specifically includes:

[0070] determining the confocal function value at a polar coordinate (θ.sub.i, m) by using a formula:

[00012]$h\left({\theta }_i,m\right)={\int }_{m\Delta }^{(m+1.Math.\Delta }{\left\{{\left[\frac{{\lambda }_0\left(z-z_0\right)}{{\alpha }_{R}n\pi w_0^2}\right]}^2+1\right\}}^{-1}\mathrm{dz};$

and

[0071] determining the sensitivity function value at the polar coordinate (θ.sub.i, m) by using a formula:

[00013]$f\left({\theta }_i,m\right)={\int }_{m\Delta }^{\left(m+1\right)\Delta }\left[{\left(\frac{\sin \zeta }{\zeta }\right)}^2\exp \left(-\frac{\mathrm{\delta \lambda }/\mathrm{\Delta \lambda }}{2\ln 2}{\zeta }^2\right)\right]\mathrm{dz}.$

[0072] Wherein, h(θ.sub.i, m) is the confocal function value at the polar coordinate (θ.sub.i, m), f(θ.sub.i, m) is the sensitivity function value at the polar coordinate (θ.sub.i, m), θ.sub.i ∈[0, 2π], m=0, 1, . . . , M, M is a radius (unit: pixel) of a frame of the intravascular optical coherence tomography cross-sectional image, Δ is a size (unit: mm) of the pixel in the image, α.sub.R is a derivative parameter of a secondary laser beam formed by an interaction of an incident beam with scatterers in a medium (α.sub.R=1 for a specular reflection; α.sub.R=2 for a diffuse reflection), n is a refractive index of a tissue to light, z.sub.0 is a focal plane position (z.sub.0 is negative when the focal plane is above the tissue, and z.sub.0 is positive when the focal plane is inside the tissue), w.sub.0 is a light spot radius after the incident light is focused, λ0 is the center wavelength of the incident light, δλ is the resolution of the spectrometer, Δλ is the wavelength interval, λ.sub.0, δλ and Δλ are the parameters of the imaging system, and ζ is a normalized value of a depth z with respect to a maximum detection depth,

[00014]$\zeta =2\pi z\frac{\mathrm{\Delta \lambda }}{{\lambda }_0^2}.$

[0073] S104, a model of an intensity of a backscattered signal being attenuated depending on a detection depth is established according to the confocal function value and the sensitivity function value of each polar coordinate to obtain a theoretical value of the backscattered signal.

[0074] Step S104 specifically includes the following steps S1041-1044:

[0075] S1041: an irradiance of the incident light at the polar coordinate (θ.sub.i, M) is determined by using a formula:

[00015]$L_{\mathrm{fh}}\left({\theta }_i,m\right)=L_{0}f\left({\theta }_i,m\right)h\left({\theta }_i,m\right)\exp \left[-\Delta \underset{l=0}{\overset{m}{.Math.}}{\mu }_t\left({\theta }_i,l\right)\right].$

[0076] S1042: an irradiance of the backscattered light at the polar coordinate (θ.sub.i, m) is determined by using a formula:

[00016]$\begin{array}{c}B\left({\theta }_i,m\right)={\mathrm{\alpha \mu }}_t\left({\theta }_i,m\right)L_{\mathrm{fh}}\left({\theta }_i,m\right)\\ ={\mathrm{\alpha \mu }}_t\left({\theta }_i,m\right)L_{0}f\left({\theta }_i,m\right)h\left({\theta }_i,m\right)\exp \left[-\Delta \underset{i=0}{\overset{m}{.Math.}}{\mu }_t\left({\theta }_i,l\right)\right]\end{array}.$

That is, the attenuation of the incident light in the tissue is mainly caused by absorption and backscattering, and the irradiance of the backscattered light is expressed as the above formula assuming that the light attenuation coefficient and the backscattered coefficient of the tissue have a linear relationship.

[0077] S1043: an irradiance of the backscattered light reaching a detector surface at the polar coordinate (θ.sub.i, m) is determined by using a formula:

[00017]$\begin{array}{c}A\left({\theta }_i,m\right)=B\left({\theta }_i,m\right)\exp \left[-\Delta \underset{i=0}{\overset{m}{.Math.}}{\mu }_t\left({\theta }_i,l\right)\right]\\ ={\mathrm{\alpha \mu }}_t\left({\theta }_i,m\right)L_{0}f\left({\theta }_i,m\right)h\left({\theta }_i,m\right)\exp \left[-2\Delta \underset{i=0}{\overset{m}{.Math.}}{\mu }_t\left({\theta }_i,l\right)\right]\end{array}.$

[0078] S1044: the theoretical value of the backscattered signal at the polar coordinate (θ.sub.i, m) is determined by using a formula:

[00018]$\begin{array}{c}{I}_b\left({\theta }_i,m\right)=n\left({\theta }_i,m\right)+\beta A\left({\theta }_i,m\right)\\ =n\left({\theta }_i,m\right)+{\mathrm{\beta \alpha \mu }}_t\left({\theta }_i,m\right)L_{0}f\left({\theta }_i,m\right)h\left({\theta }_i,m\right)\exp \\ \left[-2\Delta \underset{i=0}{\overset{m}{.Math.}}{\mu }_t\left({\theta }_i,l\right)\right]\end{array},$

i.e. taking into account an influence of noise.

[0079] Wherein, L.sub.fh(θ.sub.i, m) is the irradiance of the incident light at the polar coordinate (θ.sub.i, m) in an imaging plane, L.sub.0 is an incident light power (unit: mW or dBm), μt(θi, l) is a tissue light attenuation coefficient at the polar coordinate (θi, l), B(θ.sub.i, m) is the irradiance of the backscattered light at the polar coordinate (θ.sub.i, m), α is a proportion of the backscattered light in the incident light, A(θ.sub.i, m) is the irradiance of the backscattered light reaching the detector surface at the polar coordinate (θ.sub.i, m), I.sub.b(θ.sub.i, m) is the theoretical value of the backscattered signal at the polar coordinate (θ.sub.i, m), n(θ.sub.i, m) is a noise power value introduced by measurement at the polar coordinate (θ.sub.i, m), and β is an A/D conversion factor of the imaging system.

[0080] S105, a measured value of the backscattered signal is determined according to a gray value of each pixel in the grayed intravascular optical coherence tomography cross-sectional image.

[0081] Step S105 specifically includes:

[0082] determining a normalized measured value of the backscattered signal at a polar coordinate (θ.sub.i, m) by using a formula:

[00019]$I_m\left({\theta }_i,m\right)=\frac{g\left({\theta }_i,m\right)}{255}.$

[0083] Wherein, g(θ.sub.i, m) is a gray value at the polar coordinate (θ.sub.i, m) in the polar coordinate view of an intravascular optical coherence tomography image, and I.sub.m(θ.sub.i, m) is the normalized measured value of the backscattered signal at the polar coordinate (θ.sub.i, M).

[0084] S106, an error function is constructed according to the theoretical value of the backscattered signal and the measured value of the backscattered signal corresponding to each polar coordinate.

[0085] Step S106 specifically includes:

[0086] determining the or function by using a formula: f(μ.sub.t(θ.sub.i, m))=Ī.sub.b(θ.sub.i, m)−I.sub.m(θ.sub.i, m).

[0087] Wherein, f(.Math.) is the error function, Ī.sub.b(θ.sub.i, m) is a normalized theoretical value of the backscattered signal at a polar coordinates (θ.sub.i, M).

[0088] 8107, the error function is minimized to determine a light attenuation coefficient corresponding to each polar coordinate in the polar coordinate view.

[0089] Step S107 specifically includes the following steps S1071-1074:

[0090] S1071: an estimated value of the light attenuation coefficient is determined by using a formula:

[00020]${\hat{\mu }}_t\left({\theta }_i,m\right)=\underset{\mu \ge 0}{\arg \min }F\left(\mu \right)=\underset{\mu \ge 0}{\mathrm{argmin}}{.Math.f\left(\mu \right).Math.}^2.$

[0091] S1072: Taylor expansion approximation is performed on an objective function F(μ) near μ by using a formula:

F(μ+Δμ)=∥f(μ+Δμ)∥.sup.2≈∥f(μ)∥.sup.2+f′(μ)Δμ+½f″(μ)(Δμ).sup.2.

[0092] S1073: an optimization of the estimated value of the light attenuation coefficient is converted into solving an optimal increment Δμ for each iteration by using a formula:

[00021]$\{\begin{array}{c}\Delta \hat{\mu }=\underset{\mathrm{\Delta \mu }}{\arg \min m}\left(\mathrm{\Delta \mu }\right)\\ m\left(\mathrm{\Delta \mu }\right)={.Math.f\left(\mu \right).Math.}^2+f^{\prime }\left(\mu \right)\mathrm{\Delta \mu }+\frac{1}{2}{f}^{″}\left(\mu \right){\left(\mathrm{\Delta \mu }\right)}^2+\frac{1}{2}{\gamma \left(\mathrm{\Delta \mu }\right)}^2\end{array},$

and a damping term is added to punish Δμ;

[0093] S1074: a formula:

[00022]$\hspace{1em}\{\begin{array}{c}\Delta \hat{\mu }=\underset{\mathrm{\Delta \mu }}{\arg \min }m\left(\Delta \mu \right)\\ m\left(\Delta \mu \right)={.Math.f\left(\mu \right).Math.}^2+f^{\prime }\left(\mu \right)\Delta \mu +\frac{1}{2}{f}^{″}\left(\mu \right){\left(\Delta \mu \right)}^2+\frac{1}{2}{\gamma \left(\mathrm{\Delta \mu }\right)}^2\end{array}$

is solved, to determine the light attenuation coefficient corresponding to each polar coordinate in the polar coordinate view.

[0094] Wherein, F(μ) is an objective function, {circumflex over (μ)}.sub.t(θ.sub.i, m) is the estimated value of the light attenuation coefficient, f′(μ) and f″(μ) are the first and second derivatives of f(μ) respectively.

[0095] As shown in FIG. 4, the process for solving the formula:

[00023]$\hspace{1em}\{\begin{array}{c}\Delta \hat{\mu }=\underset{\mathrm{\Delta \mu }}{\arg \min }m\left(\Delta \mu \right)\\ m\left(\Delta \mu \right)={.Math.f\left(\mu \right).Math.}^2+f^{\prime }\left(\mu \right)\Delta \mu +\frac{1}{2}{f}^{″}\left(\mu \right){\left(\Delta \mu \right)}^2+\frac{1}{2}{\gamma \left(\mathrm{\Delta \mu }\right)}^2\end{array}$

includes following step 1 to step 4.

[0096] Step 1: the parameters is initialized, i.e. an initial iteration value μ.sub.0 and an initial damping coefficient value γ.sub.0 are set, and the iteration number k is set to be 0.

[0097] Step 2: let m′(Δμ)=0, and the following formula is obtained:

f′(μ.sub.k)+f″(μ.sub.k)Δμ.sub.k+γ.sub.kΔμ.sub.k=0.

[0098] A minimum value of m(Δμ) is solved, where pa is the light attenuation coefficient at the polar coordinate (θ.sub.i, m) after the k.sup.th iteration, Δμ.sub.k is the increment of μ.sub.k, f′(μ.sub.k) and f″(μ.sub.k) are the first and second derivatives of f(μ.sub.k) respectively, and γ.sub.k is the damping coefficient after the k.sup.th iteration.

[0099] The increment is solved according to the following formula:

Δμ.sub.k=−[f″(μ.sub.k)+γ.sub.k].sup.−1f′(μ.sub.k).

[0100] Step 3: an approximation index ρ representing a similarity degree between the linear approximation and the real objective function is calculated:

[00024]$\rho =\frac{F\left({\mu }_k+{\mathrm{\Delta \mu }}_k\right)-F\left({\mu }_k\right)}{f^{\prime }\left({\mu }_k\right){\mathrm{\Delta \mu }}_k+\frac{1}{2}{f}^{″}\left({\mu }_k\right){\left(\Delta {\mu }_k\right)}^2+\frac{1}{2}{{\gamma }_k\left({\mathrm{\Delta \mu }}_k\right)}^2}.$

[0101] Step 4: the damping coefficient is adjusted according to ρ:

[0102] (1) If ρ≤0, i.e., the variation trend of the error being opposite to the optimization objective, which indicates a poor linear approximation, let μ.sub.k+1=μ.sub.k, and decrease the damping coefficient, i.e., let γ.sub.k+1=0.5γ.sub.k, and then go to step 5.

[0103] (2) If ρ>0, i.e. the variation trend of the error being the same as the optimization objective, which indicates that the linear approximation is good at this time and an iterative update should be performed, let μ.sub.k+1=μ.sub.k+Δμ.sub.k. If 0<ρ≤0.25, reduce the damping coefficient, i.e. let γ.sub.k+1=0.5γ.sub.k. If 0.25<ρ≤0.75, keep the damping unchanged, i.e. let γ.sub.k+1=γ.sub.k. If ρ>0.75, increase the damping coefficient, i.e., let γ.sub.k+1=2γ.sub.k.

[0104] Step 5: It is determined whether the iteration is terminated according to the following formula:

[00025]$\underset{1\le i\le k+1}{\max }.Math.f^{\prime }\left({\mu }_i\right)f\left({\mu }_i\right).Math.\ge 0.01.$

[0105] If the above formula is satisfied, turn to step 6; otherwise, the iteration is terminated, aid the estimation result μ.sub.k+1 of the light attenuation coefficient is output.

[0106] Step 6: the iteration times is updated, and let k to be k+1, and turn to step 2.

[0107] S108, a spatial distribution diagram of the light attenuation coefficient is determined according to the light attenuation coefficient corresponding to each polar coordinate in the polar coordinate view.

[0108] S108 specifically includes:

[0109] normalizing and graying the polar coordinate view, and converting it into a rectangular coordinate view, namely a lateral view after the light attenuation coefficient corresponding to each polar coordinate in the polar coordinate view is obtained.

[0110] FIG. 5 is a schematic structural diagram of a system of quantitative intravascular optical coherence tomography provided by the present disclosure. As shown in FIG. 5, the system of quantitative intravascular optical coherence tomography provided by the present disclosure may include: an acquisition module 301 for intravascular optical coherence tomography cross-sectional image, a determination module 502 for a polar coordinate view, a determination module 503 for a confocal function value and a sensitivity function value, a determination module 504 for a theoretical value of the backscattered signal, a determination module 505 for a measured value of the backscattered signal, an error function constructing module 506, a determination module 307 for a light attenuation coefficient and a determination module 508 for a spatial distribution diagram of the light attenuation coefficient.

[0111] The acquisition module 501 for the intravascular optical coherence tomography cross-sectional image is configured to acquire an intravascular optical coherence tomography cross-sectional image output by an imaging system. Wherein, the intravascular optical coherence tomography cross-sectional image is a lateral view; a center of the intravascular optical coherence tomography cross-sectional image corresponds to a center of an imaging catheter, and a vascular lumen, a vascular wall intima/media and a vascular wall adventitia are around the center of the image in sequence.

[0112] The determination module 502 for the polar coordinate view is configured to gray the intravascular optical coherence tomography cross-sectional image, and convert the grayed intravascular optical coherence tomography cross-sectional image into a polar coordinate view.

[0113] The determination module 503 for the confocal function value and the sensitivity function value is configured to respectively calculate the confocal function value and the sensitivity function value of each polar coordinate in the polar coordinate view according to parameters of the imaging system. Wherein, the parameters of the imaging system comprise: a center wavelength of an incident light, a resolution of a spectrometer, and a wavelength interval.

[0114] The determination module 504 for the theoretical value of the backscattered signal is configured to establish a model of the intensity of the backscattered signal being attenuated depending on the detection depth, according to the confocal function value and the sensitivity function value of each polar coordinate to obtain a theoretical value of the backscattered signal.

[0115] The determination module 505 for the measured value of the backscattered signal is configured to determine the measured value of the backscattered signal according to a gray value of each pixel in the grayed intravascular optical coherence tomography cross-sectional image.

[0116] The error function constructing module 506 is configured to construct an error function according to the theoretical value of the backscattered signal and the measured value of the backscattered signal corresponding to each polar coordinate.

[0117] The determination module 507 for the light attenuation coefficient is configured to minimize the error function to determine a light attenuation coefficient corresponding to each polar coordinate in the polar coordinate view.

[0118] The determination module 508 for the spatial distribution diagram of the light attenuation coefficient is configured to determine the spatial distribution diagram of the light attenuation coefficient according to the light attenuation coefficient corresponding to each polar coordinate in the polar coordinate view.

[0119] The embodiments in the present specification am described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. The system disclosed by the embodiment corresponds to the method disclosed by the embodiment and thus is briefly described, and the relevant parts can refer to the portion of the method.

[0120] The principles and implementation of the present disclosure have been described herein with specific examples, and the above embodiments are presented to aid in the understanding of the methods and core concepts of the present disclosure; meanwhile, those skilled in the art may make some changes in both the detailed description and an application scope according to the teachings of this disclosure. In conclusion, the contents of the description should not be construed as limiting the disclosure.