APPARATUS, DEVICES AND METHODS FOR OBTAINING OMNIDIRECTIONAL VIEWING BY A CATHETER
CROSS-REFERENCE TO RELATED APPLICATIONS) |0 OI ] This application relates to and claims priority from U.S. Patent Application Serial os. 61 /836,716 filed on June 19, 2013, 61/934,454 filed on January 31, 2014, and
61 /905,893 filed on November 19, 2013, the entire disclosures of which are incorporated herein by reference,
FIELD OF THE DISCLOSURE
{'0002] The present disclosure relates generally to viewin by a catheter, and more specifically, to exemplary embodiments of exemplary devices, apparati and methods for omnidirectional (e.g., laser speckle} imaging, and viewing by a catheter.
BACKGROUND INFORMATION
{'0603] Despite major ad vances in coronary interventions and pharmacotherapies, acute myocardial- infarction ("AMP') remains the leading cause of death, annually claiming over 10 million lives worldwide. AMI can be caused by coronary thrombosis that can frequently result from the rupture of vulnerable plaque. If vulnerable plaques could be identified and treated prior to rupture, the incidence of AMI could be substantially reduced, and tens of thousands of lives could be saved. A key challenge in realizi ng this preventative paradigm can be that plaques wit vulnerable morphology occur at multiple sites in the coronary tree, and therefore, additional knowledge of piaque mechanical stability can be imperative in order to accurately identify plaques with the highest risk of rupture.
{'0004] Laser speckle patterns (see, e.g., References 96 and 97) can be granular intensity patterns that can arise from the interference of coherent light scattered from randomly distributed light scattering particles. The scattered photons can experience different path lengths. The phase difference between partial waves can cause constructive or destructive interference, and can produce randomly distributed high or low intensity spots called speckles. The moving scatterers can introduce different phase shifts for different partial waves, and can change the interference between partial waves, which, can lead to temporally
varying speckle patterns. The temporal evolution of speckles can provide information of seatterers' .movement (see, e.g. , Reference 98), and can further the information of the media properties which can influence the scattering particles motion (e.g., viseoeiastscity). (See, e.g.. References 99-103). Laser speckle imaging ("LSI") techniques have been applied in medical diagnosis to retrieve information about tissue perfusion (see, e.g., References 98 and 104), and mechanical properties (see, e.g.. References 99- 03) of tissues from dynamic speckle patterns. To perform LSI in vivo, coherent light can be delivered via an optical fiber, and the reflected laser speckle patterns ca be collected and transmitted via. optical fiber bundles ("OFB") incorporated within small diameter endoscopes, (see, e.g.. References 99- 103 and 105-107).
[0005] Due to their small transverse dimensions and flexibility, optical fiber bundles have been widely used i medical endoscopy (see, e.g., References 108- 1 12), and. in other minimally-invasive approaches, to enable the capability of being guided through coronaries or other conduits of human body. The large numerical aperture ("NA") can compare to the common optical fiber, and high cores density give fiber bundles can have high light collection efficiency. The high compact density of cores of fiber bundles can also provide high resolution imaging. However in LSI, the light can be highly coherent unlike the white light endoscopy (see, e.g.. Reference 1 13 in which the interference effect between cores can be neglected. The high density of cores can introduce strong coupling between adjacent fibers, which can severely affect the image quality transmitted through fiber bundles. Each fiber in the fiber bundles can support multiple guided modes, and the field of these modes can extend into the cladding, and can overlap with the mode fields of surrounding fibers. (See, e.g., Reference 1 14), Such overlapping can lead to the coupling, between modes, of individual fibers, and ioterfiher power exchange between adjacent fibers known as the optical crosstalk between fibers. Consequently, the transmitted images, or laser speckles, can be modulated by the inter-fiber crosstalk in fiber bundles due to mode coupling. During the in vivo LSI, the movement of fiber bundles due to the bulk motion of surrounding tissue can be hard to prevent. The movement of a fiber bundle can cause the core coupling changing with time, and the modulation to the transmitted speckles can be varying with time. As a result, the time-varying coupling between cores can cause erroneous speckle temporal statistics, and ca reduce the accuracy of an LSI analysis. (See, e.g.. Reference 100).
■7
}0006| The coupling between fibers modes of di fferent fibers m fiber bundles has been extensively studied based on the fiber core size, core spacing, NA and non-uniformity of cores. (See, e.g., References ί ! 5-1 17). However these studies mainly focused on the coupling between fundamental modes of neighboring fibers. Only a few numerical
simulations (see, e.g., Reference 118) and experiments (see, e.g.. References 100. 1 16 and 119) have been conducted to show the fiber crosstalk in multimode fiber bundles and its influence on image transmission. These numerical simulations (see, e.g., Reference 118) only simulated fields propagating a few rnro along the length of bundles due to the intensive computing required. Previous experiments have demonstrated leached fiber bundles can effectively reduce the cross talk between cores, and can obtain a relatively stable temporal decollation function of transmitted speckles during bundles motion because of large core- to-core separation due to manufacturing processes of leached fiber bundle, (See, e.g..
Reference 100). However, the effect of mode coupling between neighboring optical, fibers on the transmission of laser speckles may not be well understood.
[0007.1 In order to conduct laser speckle imaging via a catheter, light can be guided through an optical fiber and distal optical components to illuminate a single spot on the cylindrical lumen to collect reflected speckle patterns via a single fiber or collection of optical fibers (e.g., a fiber bundle). To conduct circumferential mapping, the catheter can be rotated during pull-back. However, this can introduce motion artifacts during catheter rotation that can confound the ability to accurately analyze laser speckle patterns from tissue. |ββ08} Thus, it may be beneficial to provide an exemplary device, apparatus and method for viewing by a catheter, which can overcome at least some of the deficiencies described herein above. SUMMARY OF EXEMPLARY EMBODIMENTS
J0009] These and other objects of the present disclosure can be achieved by provision of an apparatus for obtaining information regarding a biological structure^), which can include, for example a. Sight guiding arrangement which can include a. fiber through which an
electromagnetic radiaiion(s) can be propagated, where the electromagnetic radiation can be provided to or from the structure. An at least partially reflective .arrangement can have multiple surfaces, where the reflecting arrangement can be situated with respect to the optical
arrangement such, thai the surfaces thereof each can recei ve a(s) beam of the electromagnetic radiations instantaneously, and a receiving anangeraent(s) which can be configured to receive the reflected radiation from the surfaces which include speckle patterns,
f OO!Oj In some exemplary embodiments, a polarizing arrangements ) can be incl ded which can receive the electromagnetic radiation, and prevent receipt of a same polarization from returning to the receiving arrangements). The reflective arrangement can have a portion(s) with a shape of a cone, a poly gon or a pyramid. An optical arrangement can be included which can be configured, to receive the electromagnetic radiation(s). The optical arrangement can include a GRIM lens, a ball lens, or an imaging lens. The number surfaces of the reflective arrangement can be 2 or more, 4 or more, or 6 or more. The light guiding arrangement can include a configuratio which can split the electromagnetic radiation to further radiations having different wavelengths where the multiple surfaces can reflect the further radiations, and where the receiving arrangement(s) can be further configured to receive the reflected further radiations provided at the different wavelengths,
[0011] In another exemplary embodiment of the present disclosure can be an apparatus for obtaining information regarding a biological stxucture(s), which can include, for example, catheter arrangement which can include a fiber(s) through which an electromagnetic radiation(s) can be propagated, where the electromagnetic radiation can be provided to or from the structure. A pu!!back arrangement can be configured to facilitate a puliback of the catheter arrangement, and a detector arransenient can includes a plurality of sensors, the sensors being coupled to a surface of a portion(s) of the catheter arrangement, and configured to move together with the pu!Jback arrangement, and receive optical information associated with the electromagnetic radiation (s) provided from the structure so as to generate the information.
[0012j In some exemplary embodiments of the present disclosure, the sensors can be directly attached to the surface of the portion(s) of the catheter arrangement. The detector arrangement can include a CMOS sensor, a CCD sensor, a photodetector or a photodetector array. The fiber(s) can include a plurality of fibers, or a fiber bundle. For example, at an illumination wavelength between about 630~720nm, the fiber bundle can have (i) a core diameter of 3.0 μ,οι ± 0.3 μ«ι with a fluctuation in core diameter of ± 0.03 gra to ± 0.3 μηι, (ii) a numerical aperture of at least 0.3 S, and (iii) a core spacing of 8.0 urn ± 0.5 μχη, The
puiiback arrangement cars, he controlled by a motor(s). The motor(s) can control the puliback arrangement such that, the puliback arrangement can move die catheter and detector
arrangements In a stepped manner. The motor(s) can control the puliback arrangement to rotate a drive shaft or distal optics. The motor(s) can be configured to keep the catheter stationary.
f 0013] in some exemplary embodiments of the present: disclosure, adjacent times for the movements of the catheter and detector arrangement can be between 5 msec and 100 msec. The sensors can receive the optical information that can be associated with the
electromagnetic radiation provided at different wavelengths, A filter arrangement can be configured to filter tlie optical information based on the electromagnetic radiation provided at different wavelengths. The catheter arrangement can include a drive shaft arrangement which can hold the fiberfs) and can be directly connected to the puliback arrangement. The driv shaft arrangement can further hold distai optics. The motor(s) can control the puliback arrangement such that the puiiback arrangement can move the catheter and detector arrangements continuously at a predetermined speed or a variable speed,
10014] In another exemplary embodiment of the present disclosure is an apparatus for imaging portion(s) of a biological structure, which can include, for example a radiation providing arrangement which can be configured to forward a first electromagnetic radiation(s) to the structure at multiple illumination locations. A detector arrangement can be is configured to receive a second electromagnetic radiation(s) from the multiple locations of the structure. A puliback arrangement which, during the forwarding of the first
electromagnetic radiation, can be configured to pull back tire radiation arrangements) of the detector arrangement. The detector arrangement can be further configured to image the portion(s) of the structure based on the second electromagnetic radiation(s), without a rotation of the radiation providing arrangement,
('0015] In some exemplary embodiments of the present disclosure, the first
electromagnetic radiation(s) can be forwarded to the structure at multiple illumination locations substantially simultaneously. The second electromagnetic radiation(s) can be received from the multiple locations of the structure substantially simultaneously. The second electromagnetic radiation(s) ca provides informatio regarding a speckle pattern reflected from the portion(s) of the structure. The speckle pattern can have an intensity that
can vary in time. The variation of the ntensity of the speckle paitera can provide information regarding mechanical properties of the portion(s) of the structure, which can be determined by the detector arrangement. The piiilback arrangement can be controlled, by a motor(s). The rnotor(s) can control the puilback arrangement to rotate a drive shaft or distal optics.
| 0Θ16| In another exemplary embodiment of the present disclosure is a method for imaging portion(s) of a biological structure, which ca include, for example, using a radiation providing arrangement, forwarding a first electromagnetic radiation(s) to the structure at multiple illumination locations, using a detector arrangement, receiving a second
electromagnetic radiations s) from the multiple locations of the structure, and pulling back the radiation arrangements) of the detector arrangement. The detector arrangement c an be configured to image substantially an entire surface of the portion(s) of the structure based on the second electromagnetic radiation(s), without a rotation of the radiation providing
arrangement.
f 0017] In some exemplary embodiments of the present disclosure, the first
electromagnetic radiation(s) can be forwarded to the structure at multiple illumination locations substantially simultaneously. The second electromagnetic^) radiation can be received from the multiple locations of the structure substantially simultaneously.
iO lSJ In another embodiment of the present disclosure, a system, method and computer- accessible medium can be provided for obtaining information regarding a biological structure(s), which can include, for example, receiving information related to a radiation(s) reflected from the biological structure(s including a speckle paitern(s), and generating an image of the biological structure^) based on the information. PixUation artifacts can be removed from the speckle pattem(s), A speckle intensity fluctuation of the speckle pattern(s) can be determined by, for example measuring a change of multiple mirror facets over rime. A background fluctuation or a source fluctuation from can be removed from the speckie pattern(s). Non-fluctuating speckles can be filtered from the speckle pattern(s). A phase fluctuation of the reflected radiaiion(s) can be determined to, for example, characterize the tissue.
[0019] in another embodiment of the present disclosure, a method of reducing inter-fiber crosstalk in a fiber optic bundle can be provided, which can include, for example providing a fiber optic bundle comprising a plurality of core fibers each having a core diameter of 3.0 μπι
±■ 0.3 μι« with a flacft.iatio.ft in core diameter of £ 0.03 μνα to 0.3 μηι, a numerical aperture of at least 0.35, and the fiber optic bundle having a core spacing of 8.0 urn ± 0.5 μηι.
Receiving a light into the fiber optic bundle, where the fiber optic bundle has a reduced inter- fiber crosstalk. The diameter of the core fibers can be 3.0 pm 0.2 μτη. The core fibers can be 3.0 μιη ± 0.1 μηι. The fluctuation of the core diameter can be ± 0.06 pra to ± 0.2 pm. The fluctuation of the core diameter can be approximately ± 0.1 μιη. The numerical aperture can be between 0.38 and 0.41. The core spacing can be 8.0 μνα. ± 0.3 μιη. The core spacing can be 8.0 pm ± 0.2 pro. The fiber optic bundle can have an inter-fiber crosstalk, that can be at least .10% less than the inter-fiber crosstalk, within a leached fiber optic image bundle defined as SCHOTT North America Type i at a propagation distance of 0.5 m using 690 ran radiation. The fiber optic bundle can have an inter-fiber crosstalk that can be negligible.
[0020] In another exemplary embodiment of the present disclosure is, an apparatus can be provided for laser speckle imaging that has low inter-fiber crosstalk, which can include, for example a coherent radiation source, a fiber optic bundle configured to receive radiation from the coherent radiation source including a plurality of core fibers having a core diameter, a fluctuation in core diameter, a numerical aperture, and a core spacing, where each of the core diameter, the fluctuation in core diameter, the numerical aperture, and the core spacing can be determined using coupled mode theory ("C T "). One or more optical elements can be configured to direct coherent radiation from the fiber optic bundle to a tissue and collect radiation from the tissue. A detector can be configured to receive a speckle pattern from the one or more optical elements. The diameter of the core fibers can be 3.0 pm ± 0.3 urn The diameter of the core fibers can be 3.0 μηι 0.2 pm. The diameter of the core fibers can be 3.0 pm ± 0, 1 μπι The fluctuation of the core diameter can be ± 0,03 prn to ± 0,3 pm. The fluctuation of the core diameter ca be * 0.05 pro to ± 0.2 μηι. The fluctuation of the core diameter can be approximately ± 0.3 μηχ The numerical aperture can be at least 0.35. The numerical aperture can be between 0.38 and 0. 1. The core spacing can be 8.0 pm ·±· 0,5 pra. The core spacing can be 8.0 μτη ± 0.3 μηι. The core spacing can be 8.0 pro ·.!.· 0.2 pm. The fiber optic bundle can include a core diameter of 3.0 pra ± 0.3 pm with a fluctuation in core diameter of '·±· 0.1. pm to ± 0.3 pm, and a numerical aperture of at least 0.35, and the fiber optic bundle having a core spacing of 8.0 pm ± 0,5 μη . The numerical aperture and the core
spacing can depend on a wave length of the coherent radiation source. The core diameter can depend on the fiber size, the core spacing or the numerical aperture.
(0021 ] In another exemplary embodiment of the present disclosure, a method can be provided for tissue analysis, which can include, for example illuminating a first cylindrical section(s) of a lumen wall with coherent or partially coherent light by passing the light through a facet(s) of a multiple-faceted pyramidal mirror, receiving light, reflected from the first cylindrical section of a lumen wall at the mirror, illuminating a second cylindrical section(s) of a lumen wall with coherent or partially coherent light, at a time different .from the first illuminating step by passing the light through a second faeef(s) of the multiple-faceted pyramidal mirror, receiving light reflected from the second cy lindrical sec tion of a lumen wall at the mirror, receiving light reflected from the mirror at a detector and forming series of speckle patterns, and analyzing changes in the speckle patterns at time intervals sufficient to measure changes caused by microscopic motion of objects within the tissue,
f 0022] According to a still further exemplary embodiments of the present disclosure, the illumination can occ ur by first illumina ting cylindrical section of a lumen wall through either a single facet of the pyramidal mirror at a tune, or multiple facets of the pyramidal mirror at one time, where the facets are not adjacent to each other. The multiple faceted pyramidal mirror can be a four-sided mirror and the cylindrical section of a lumen wall can be illuminated through two non-adjacent facets simultaneously and then the cylindrical section of a lumen wall can be illuminated through the other two non-adjacent facets simultaneously. The multiple faceted pyramidal mirror can be a six-sided mirror and the cylindrical section of a lumen wall i illuminated through two or three non-adjacent facets simultaneously.
j'0923] These and other objects, features and advantages of the exemplary embodiments of the present disclosure will become apparent upon reading the following detailed description of the exemplary embodiments of the present disclosure, when taken in conjunction with the appended claims.
BRIEF DESCRIPTION OF THE DRAWINGS
{092 ] Further objects, features and advantages of the present disclosure will become apparent, from the following detailed description taken in conjunction with the accompanying Figures showing illustrative embodiments of the present disclosure, in which:
[0025] FIG . 1 is a set of exemplary images of a plaque and a color map of the plaque according to an exemplary embodiment of the present disclosure;
(0026] FIG. 2 is an exemplary speckle image according to an exemplary embodiment of the present disclosure;
f 027] FIG. 3 is an exemplary graph illustrating g:2(t) corves according to a exemplary embodiment of the presen t di sclosure;
(0028] FIG. 4 is an exemplary graph illustrating mean T for different plaque groups according to an exemplary embodiment of the present disclosure;
('8029] FIG. 5 is an exemplary graph illustrating r alues according to an exemplary embodiment of the present disclosure;
(0030] FIG, 6 is an exemplary graph illustrating an estimation of the overall bulk modulus of a necrotic core fibroatheroma as a function of fibrous cap thickness according to an exemplary embodiment of the present disclosure;
(0031 ] FIGS. 7A-7C are exemplary graphs illustrating the evaluation of spatial heterogeneity by beam scanning according to an exemplary embodiment of the present disclosure;
[0032] FIG. 8 is an exemplary colormap illustrating depth imaging in a thin cap fibroath.ero.ma according to an exemplary embodiment of the present disclosure;
(0033] FIG. 9 is an exemplary graph illustrating average plaque X measured via an exemplary ieached fiber bundle according to an exemplary embodiment of the present disclosure;
{0934] FIG. 10 is an exemplary schematic of an exemplary ILSO catheter according to an exemplary embodiment of the present disclosure;
(0035] FIG. 1 1 is an exemplary image of an exemplary LSI catheter sheath according to an exemplary embodiment of the present disclosure;
{0936] FIG. 12 is an exemplary schematic of an exemplary ILSO procedure in a swine xenograft model according to an exemplary embodiment of the present disclosure;
(0037] FIG. 1 is an exemplary graph illustrating average T calculated for 3 plaque groups according to an exemplary embodiment of the present disclosure;
(0038] FIG. 14 is an exemplary graph illustrating X calculated in a swine using exemplary PBO procedures according to an exemplary embodiment of the present disclosure;
(6639) FIG. 15 is an exemplary schematic of an exemplary motor drive assembly for helical scanning according to an exemplary embodiment of the present disclosure;
{0040] FIGS 16 A and I6B are exemplary images and color maps of T over two NC plaques according to an exemplary embodiment of the present disclosure;
J 0041J FIG. 17 A is an exemplary graph illustrating a 3 D distribution of mean penetration depths collected over a catheter according to an exemplary embodiment of the present disclosure;
|M42] FIGS. 17B and 17C are exemplary co!ormaps illustrating cross-sectional distributions along x and y of FIG. 17A according to an exemplary embodiment of the present disclosure;
{6643 j FIGS. 18A and 1 SB are exemplary graphs illustrating spatial resolution estimated using Monte-Carlo Ray Tracing according to an exemplary embodiment of the present disclosure;
{'0644] FIG. 19 is an exemplary OFDi image obtained during visipaque flushing according to an exemplary embodiment of the present disclosure;
{6045] FIG , 20 A. is an image of cross sec tion of a leached fiber bundle accordi ng to an exemplary embodiment of the present disclosure;
{'0046] FIG. 2GB is a schematic of fiber bundle in numerical calculations according to an exemplaty embodiment of the present disclosure;
{©047] FIG. 21 A is an exemplary image illustrating the amplitude of coupling coefficient x between all the 19x7 modes according to an exemplary embodiment of the present disclosure;
{'0048] FIGS. 21 B-2 IF are exemplary graphs illustrating the intensity of different order modes of central fiber coupled to the corresponding modes of surround fibers with propagaticm distance z for I st, 2nd, 6th, 9th and 10th mode respecti ely according to an exemplary embodiment of the present disclosure;
{0049] FIG. 22A-22C are exemplary graphs illustrating core spacing according to a exemplary embodiment of the present disclosure;
[0050] FIGS. 22D-22I are exemplary graphs illustrating thai the coupling strength can increase as core size increase, core spacing decrease and A decrease according to an exemplary embodiment of the present disclosure;
0051] FIGS. 23Α-23Ϊ are exemplary graphs and exemplary images illustrating illustrate how speckle pattern change with propagation distance due to crosstalk between neighboring cores according to an exemplary embodiment of the present disclosure;
{'0052] FIGS. 24A-24C are exemplary graphs illustrating the reduced change of intensity in each fiber of optical fiber bundles according to an exemplary embodiment of the present disclosure;
[0053} FIGS. 25A-25C are exemplary graphs illustrating that the intensity in each core of 7 core structure can with propagation according to an exemplar embodiment of the present disclosure;
[0054] FIG. 26A is an exemplary image of a small region of an optical fiber bundle cross section according to an exemplary embodiment of the presen t disc l osure;
[0055] FIG. 26B and 26C are exemplary recorded raw speckle images and its Fourier transform according to an exemplary embodiment of the pr esent disc los ure ;
[0056] FIG, 26D is an exemplary image of a Fourier transformed speckle pattern superposed by a Butterworth filter according to an exemplary embodiment of the present disclosure;
|005?|
[0058] FIGS. 27A and 27B are an exemplary image and its corresponding exemplary graph illustrating the notch filter according to an exemplar}' embodiment of the present disclosure;
[0059] FIGS, 28A and 28B are exemplary graphs illustrating the temporal response of the total intensity of speckle patterns according to an exemplary embodiment of the present disclosure;
[0060} FIG, 29 is an exemplary coJormap illustrating the spatially smoothed speckle pattern average over time according to an exemplary embodiment of the present disclosure;
[0061] FIG. 3 OA is an exemplary image of an exemplary speckle pattern with pi elation a rtifact removed according to an exemplary embodiment of the presen t disc l osure;
[0062] FIG. 3 OB is an exemplary graph illustrating the autocovariance curves of speckles within small windows according to an exemplary embodiment of the present disclosure;
[0063] FIG. 3 1 A is an exemplary image of an acrylaraide gel phantom in a 3D printed mold according to an exemplary embodiment of the present disclosure;
[0064] FIG. 3 IB is a set of 8 T maps of gels A, B, B, and C from FIG. 31 A according to an exemplar embodiment of the present disclosure;
[0065] FIG. 3.1C is an exemplary image of swine aorta with butter injected in between die aorta layers according to an exemplary embodiment of the present disclosure; {0066} FIG. 3 1 D is a set of the two longitudinal stitched τ maps of a tube according to an exemplary embodiment of the present disclosure;
[0067] FIGS. 32A~A2C shows illustrations of an example of wrapping 2D time constan maps onto a cylinder to form a cylindrical view of the time constant maps according to an exemplary embodiment of the present disclosure;
|0068] FIGS. 32D and 32E are exemplary coiormaps of a speckle intensity pattern and the retrieved phase pattern using the exemplary 2D Hilbert transform according to an exemplary embodiment of the present, disclosure;
[0069] FIG. 32F is an exemplary image illustrating locations of the optical vortices according to an exemplary embodiment of the present disclosure;
10070} FIG. 32G is an exemplary graph illustrating speckle intensity autocorrelations for two different speckle sequence according to an exemplary embodiment of the present disclosure;
[0071] FIG. 3211 is an exemplary graph illustrating the locations of the opiical vortex at different speckle frames for the fast varying speckle sequence according to an exemplary embodiment of the present disclosure;
j'0072] FIG. 321 is an exemplary graph illustrating the locations of the optical vortex at different speckle frames for the slow varying speckle sequence according to an exemplary embodiment of the presen t disclosure;
[0073} FIGS. 33A-33G are exemplar images of exemplar patterns according to an exemplary embodiment of the present disclosure:
[0074] FIG. 34 is an image of an exemplary speckle pattern according to an exemplary embodiment of the present disclosure;
[0075] FIG . 35 is an exemplary image of an exemplar)? colon-nap according to an exemplary embodiment of the present disclosure:
{0076{ FIG. 36 is an exemplary schematic illustrating an exemplary mechanism for the displacement of blood during imaging according to an exemplary embodiment of the present disclosure;
{'0077] FIG. 37 is an exemplary image of exemplary M-mode OFDI according to an exemplary embodiment of the present disclosure;
{0078 FIGS. 3SA-3SL are exemplary schematic diagrams of exemplary catheters according to an exemplary embodiment of the present disclosure;.
{'0079] FIGS. 39A~39H are exemplary images of exemplary laser spots according to an exemplary embodiment of the present disclosure; and
{0080] FIG. 40 is an illustration of an exemplary block diagram of an exemplary system in accordance with certain exemplary embodiments of the present disclosure.
{0081] Throughout the drawings, the same reference numerals and characters, unless otherwise stated, are used to denote like features, elements, components or portions of the illustrated embodiments. Moreover, while the present disclosure will now be described in detail with reference to the figures, it is done so in connection with the illustrative
embodiments and is not limited by the particular embodiments illustrated in the figures and appended claims.
DETAILED DESCRIPTION OF EXEMPLARY EMBODIMENTS
{0082] One exemplary object of the presen t disclosure can be to provide, for patient use, an optical system and method that can be termed Intracoronary Laser Speckle Imaging ("ILSI"), which can evaluate plaque viscoelastic properties, known to be intimately linked with the risk of coronary plaque rapture. It has been determined that plaque rupture can occur when the atheroma, with severely compromised viscoelastic properties, can fail to withstand stresses exerted upon it. Therefore, an important ability of an I LSI exemplary system and method, according to an exemplary embodiment of the present disclosure, ca be to evaluate plaque viscoe la sticky, to facilitate an improved understanding of plaque stability, and advance clinical capability for the detection of vulnerable plaques with the highest risk of rupture in patients.
{0083] The exemplary ILSI technology, according to an exemplary embodiment of the present disclosure, can be based on an exemplary laser speckle approach. For example, laser
speckle, a grainy pattern formed by the interference of laser light scattered from tissue, can be dynamically modulated by endogenous particular Brownian motion governed by the
mechanicai susceptibility of tissue. It has been previously demonstrated that the time scale of speckle modulations, defined by the speckle decoi relation time constant, can provide a highly sensitive metric of viscoelasticity that can be closely related with plaque composition and mechanical moduli. Given the potential impact of ILSi in measuring a key mechanicai metric of plaque stability, it can be possible to utilize the exemplary ILSI systems and methods for evaluating coronary plaques in patients. The large size (e.g., approximately 1.5mm), and limited point sampling capabi lity of existing ILSI devices, however., can render it less than optimal for human use. Therefore, it can be possible to provide a miniaturized ILSI catheter for human use which can facilitate scanning of the entire circumference and length of the coronar artery to evaluate maps of arterial viscoelasticity distribution.
Exemplary Object 1: Provide ILSI Catheter And Console Saitable For Human Use
[0084] The exemplary ILSI technology can utilize a miniaturized intracoronary catheter (e.g., i .0 mm) to acquire speckle images from the arterial wall. The catheter can be interfaced with a high-speed console to facilitate helical scanning of the coronary artery. Speckle analysis and visualization procedures can be implemented to reconstruct cylindrical maps of arterial time constants over the circumference and length of coronary segments. Performance benchmarks including catheter size and mechanical characteristics, imaging time and spatial resolution can be optimized and verified in human cadaveric hearts. The exemplary ILSI catheter performance can be evaluated in living s wine using a human to swine coronary graft model to facilitate imaging of human eoronaries under physiologic conditions.
Exemplary Object 2: Indicate Safety And Feasibility Of ILSI Technology In Patients
[0085] Regulatory procedures involving catheter evaluation, hiocompaiibilhy testing., laser exposure evaluation and safety studies in living swine can be performed in preparation for hospital IRB and FDA applications to conduct clinical studies. ISO 10555 requirements can be followed to fabricate sterile, single-use catheters in a controlled. Class 10000 GLP clean room facility at the Welhnan Center. Following .regulatory approval, ILSI can be conducted in 20 patients undergoing percutaneous coronary intervention. Cylindrical maps of arterial
mechanical properties measured by ILSI can be registered, and correlated with microstructural information obtained using intracoronary OFDI,
10086] The exemplary ILSI can provide a too! to significantly advance current scientific understanding of vulnerable plaque instability in patients. It can also provide a powerful diagnostic role within a comprehensive clinical paradigm of AMI management to fecilitate an identification of plaques with the highest risk of rupture for treatment prior to adverse events in patients,
EXEMPLARY STRATEGY
[U087J Exemplary Vulnerable plaque detection: AMI, frequently caused by the rapture of vulnerable coronary plaque, claims more lives worldwide than cancer, accidents and AIDS combined. Autopsy studies reveal a type of plaque, the thin cap fiutopsy stud ("TCFA") implicated at the site of culprit thrombi in >70% of patients who have succumbed to AMI. (see, e.g., References 1 and 2), TCFA's can be most frequently found within the proximal approximately 5 cm of the major coronary arteries and can be histologically hallmarked by the presence of a thin fibrous cap (e.g., < 65pm), rich in macrophages, overlying a large necrotic lipid pool. (See, e.g.. References 1 -5). A number of technologies such as optical coherence tomography ("OCT'), virtual histology intravascular ultrasound ("VB-IVUS"), computed tomography ("CT"), angioscopy and near infrared spectroscopy ("MRS") have been investigated in patients to evaluate key morphologic features such as fibrous cap thickness, plaque burden, calcific nodules and lipid content. (See, e.g., References 6-26) An iirm rtant cha
patients can, be that plaques, with similar vulnerable .morphologic features d not all possess an equal likelihood of rupture. For example, in 70% of patients dying from AMI, multiple TCP A 's can be found without rupture at sites remote from the culprit plaque and in non- culprit arteries,(see, e.g.. Reference 2) and can appear with similar frequency i stable patients with asymptomatic coronary artery disease ("CAD"). (See, e.g.. References 2, 27- 29). Moreover, in approximately 20% of cases, plaque rapture can be observed in necrotic core ("NC") lesions with thicker fibrous caps (e.g., >I00pm), intra-plaque hemorrhage or calcific nodules, (See, e.g., References 2. 28, 30-31). These findings call into question the current detection paradigm that relies entirely on morphologic criteria, and highlights, the
need to augment morphologic findings with important, surrogate metrics, such as mechanical metrics, in order to accurately evaluate the risk of plaque rupture. (See, e.g., References 1 and 2),
|00S8| Knowledge of mechanical metrics can be important to accurately determine the risk of p laque rupture. The atheroma can be viscoelastk in nature, exhibiting both liquid (e.g., viscous) and solid (e.g., elastic) behavior. During the pathogenesis of atherosclerosis, from lesion initiation to rupture, the vtscoelastie properties of the plaque can be altered by a complex milieu of hemodynamic and biochemical processes. The ultimate event of plaque rupture can be a biomechanieal failure that can occur when a plaque with severely
compromised mechanical properties can be unable to withstand loads exerted on it. (see, e.g., References 32-41 ). Therefore, in order to identify plaques with the highest risk of .rupture, it can be important to complement morphologic information provided by current technologies with knowledge of viscoelasiic properties.
[0089] Current knowledge of plaque viscoelasiicity however can be l imited as it can largely be derived from ex vivo mechanic l testing of cadaveric and animal arteries. These measurements can provide only a retrospective snapshot of bulk properties, limiting the understanding of how mechanical metrics can be altered during the plaque remodeling in vivo. Therefore;, important estimates of plaque viscoelasiic properties predisposed to the final event of rupture can be currently uaknowa. Crucial questions remain on how current knowledge of plaque mechanical stability translates in vivo, restricting the opportunity for accurate detection of high-risk vulnerable plaques in patients. Toge ther, these factors can highlight a important barrier in the field: the ability to detect plaques with the highest risk of rupture can be significantly hindered by the absence of tools for the mechanical
characterization of coronary plaques in patients.
[8090] An exemplary embodiment of an 1 LSI system and method according to the present disclosure can be provided for a clinical use that can be used to evaluate the viscoelasiic characteristics of coronary plaques in patien ts. The exemplary ILSI systems and methods can measures plaque viscoelasiicity by utilizing an exemplary laser speckle approach developed in a laboratory, which can interrogate the ensemble Brownian. motion dynamics of light scattering particles intimately linked with the micromechanical behavior of the atheroma.
The exemplary ILSI systems and methods can measure an index of viscoelasiicity defined by
the speckle decorrelation time constant (* ) that can be highly sensitive to minute alterations in the viscoelasiic properties of the atheroma (e.g. , Section C). (See, e.g., References 42-46). The exemplary ILSI systems and methods can provide an improved understanding of human CAD and advance clinical capability to detect plaques with the highest risk of rupture in patients as discussed below.
10091] Exemplary Understanding of CAD bv ILSi: The exemplary ILSi technology can provide important mechanical metrics implicated in plaque instability in animals and patients. The miniaturized ILSi catheter (e.g., 1mm) can facilitate evaluation of smal l coronary vessels and flo w limiting lesions. The reconstruction of 2D maps (e.g., FIG. I) can provide knowledge of viscoelasiic ity distributions over the circumference and length of the coronary vessel. The capability to evaluate depth-resolved 3D information at high spatial resolutions (e.g., approximately 100 pin) can be provided to facilitate an important understanding of the mechanical properties of the lipid pool and fibrous cap in NC plaques of highest clinical relevance. The superior sensitivity of the exemplary ILSI systems and methods described herein to minute alterations in viscoelastictt can be utilized for plaque remodeling during the natural history of coronary atherosclerosis leading to rupture, it can be known that in early lesions, inflammatory processes can influence the accumulation of low viscosity lipid. (See, e.g., References 47 and 48). in advanced plaques, apoptosis of foam ceils and mtrap!aque hemorrhage can result in large necrotic lipid pools of further reduced viscosity. (See, e.g., References 49 and 50). Furthermore, lipid pool viscosity can also be influenced by cholesterol, phospholipids and triglyceride content. (See, e.g.. Reference 50). ILSI measurements of lipid pool viscosity can provide insights on the load bearing properties of the atheroma, and can offer a likely explanation for why TCFAs do not all possess the equal likelihood of rupture. The mechanical properties and morphology of the fibrous ca can be radically altered by a net reduction in collagen content that can occur due to an imbalance in collagen proteolysis by matrix metalloproteinases Γ'ΜΜΡ") and synthesis due to apoptosis of smooth muscle cells. (See, e.g., References 51-53). ILSi can provide knowledge of important estimates of fibrous cap viscoelasticity related with the final event of plaque rupture.
f 0092] Finite element ('ΤΙ ') studies of coronary cross- sections derived from histology sections, or 1VUS and OCT images can show that peak stresses associated with plaque
rupture can. he dependent on the geometry and viscoe.1as.t-o propert ies of the fibrous cap and lipid pool, and plaque rapture can become imminent when the peak stress in the plaque surpasses an important amplitude, (See, e.g.. References 32-41, 54 and 55). Precise measurement of peak stress amplitudes predisposed to rupture needs accurate estimates of the viscoelastic properties of plaque components in situ. ILSI can help address this challenge; combining FEA approaches with ILSi maps of viscoelasticity distributions can prov de a powerful new method for accurate evaluation of peak stress in situ.
| O93] The spontaneous rupture of coronary plaques leading to. AMI can be unique in human CAD. Because there can be no realistic animal models available that can mimic this event under physiologic conditions, many key hvpotheses that relate mechanical metrics with the final event of plaque rupture can only be best studied in human patients. Exemplary embodiments of the present disclosure address this challenge by providing translating ILSI for use in patients.
ΙΌ994] The exemplary ILSI systems and methods can be used for the detection of vulnerable plaques in patients at risk for AMI, Recent clinical studies show that 10% of patients undergoing PCI and statin therapy following the first acute event develop a second adverse event due to plaque rupture within 3 years. (See, e.g.. References 56 and 57). The exemplar}'' I LSI systems and methods can be used by interventional cardiologists to detect potential plaques such that a second major adverse event can be prevented. Thus, over 100,000 people annually in the USA alone can benefit by ILSI screening.
{0Θ95] In order to reduce mortality due to AMI in the general population, new
preventative paradigms for AMI management can be realized. These exemplary paradigms can use a comprehensive screening strategy to identify at risk patients and detect high-risk vulnerable plaques in these patients such that they can he treated, prior to AMI. Non-invasive imaging of the coronary tree using computed tomography ("CP) and magnetic resonance imaging ("MRI") approaches can be important, in identifying asymptomatic patients at highest risk of AMI. (See, e.g.., References 58 and 59). These approaches, however, lack sufficient sensitivity and resolution to evaluate mechanical and morphological characteristics to detect vulnerable coronary plaques. A second level, of intracoronary screening using ILSI can be used to evaluate plaques with compromised mechanical stability likely to cause AMI in asymptomatic patients at risk.
[0096] Furthermore, the exemplary ILSI systems and methods elegantly can be used for an integration with other intracoronary technologies such as optical coherence tomography
(' CX"} and optical frequency domain imaging ("OFDi") or intravascular ultrasound ('IVUS"), (see, e.g., References 60-62'X'see, e.g., Reference 57) to render powerful approaches that can place mechanical findings within a morphologic context for a composite evaluation of plaque stability.
(0097] Further, treatments for stabilization, including low force self-expanding and bio- absorbable stents, vascular tissue implants, stem cell and. photodynamic therapy, can be de veloped by a number of companies and groups. These therapeutic interventions can utilize diagnost ic tools for the accurate diagnosis or determination of rupture-prone coronary plaques prior to treatment.
Brief Description of Exemplary ILSI Systems and Mettods
|0O98] It can be important to provide a tool to evaluate the viscoelasiie properties of coronary plaques in patients, iVUS-based cystography has been developed to measure plaque strains in response to intraluminal stress. However, evaluation of plaque viscoelasttc properties can be intractable using this approach. (See, e.g.. References 63 and 64), While recent studies have utilized inverse methods and deformabie curves to reconstruct Young's moduli from e'lastography strain maps, the approximation of linear elastic behavior b these methods can restrict accurate evaluation of the load bearing properties of viscoelastic tissue components and Sow viscosity lipid pools. It can be possible to apply elastography approaches to OCT to provide higher resolution for strain estimation relative to IVUS. Loss of OCT signal in lipid rich tissue, however, can prec lude strain assessment in NC plaques, thus significantly limiting clinical utility.65-67 There are no known techniques that can evaluate the viscoelasiic properties of coronary plaques in patients.
[0099] The following benefits can be provided by the exemplary ILSI systems and methods: I ) a measurement of plaque viscoe!asticity that cannot be accomplished by any other previously known technique. 2) facilitation of a clinical grade ILSI device for use in patients, 3) The use of the exemplary ILSI device in human translation.
[00100f The exemplary ILSI device can facilitate a comprehensive screening of the arterial circumference over long coronary segments to evaluate plaque vtscoelasticity maps at spatial resolution approximately Ι ΟΟμηι. It can be possible to provide an exemplary miniaturized
ILSI catheter (e.g., 2.4-3.0 F) including one or more low cross-talk: fiber bundles with sufficient motion tolerance to evaluate the coronary wail in vivo. To achieve capability for helical scanning, it can be possible to provide an exemplary optical rotary junction and motor drive assembly that can couple and recei ve light from multiples cores of the fiber bundle while simultaneously rotating and translating the catheter during imaging. In order to facilitate a collection of arterial speckle decorrelation information, a programmable stepper motor can be utilized to encode and transmit torque to the catheter in discrete increments, which can facilitate sufficien t sampling of the coronary circumference at a rotational rate of approximately I Hz. The exemplary ILSI device can utilize a high-speed complementary metal oxide semiconductor ("CMOS") camera (e.g., 2kHz frame rate) to obtain 7
measurements over very short time scales (e.g., 25 ms) over which the influence of low frequency arterial deformations induced by cardiac (approximately 1Hz) or respiratory (approximately 0.2 Hz) motion can be largely mitigated. Using this exemplary approach, ILSI measurements can be accomplished without the need for electrocardiogram ('ΈΚ(Γ') gating in vivo.
[001011 The exemplary ILSI device does not need apriori approximations on plaque geometrv or loading conditions to measure viscoeiasticitv, therefore, automated ILSI analvsis can be rapidly accomplished rendering ease of use in the catheterization suite. Cylindrical 2D maps of plaque viscoeiasticity can be provided from ILSI data to measure the influence of spatial heterogeneities. Exemplary methods can be provided that can utilize a spatio- temporal speckle analysis in conjunction with Monte Carlo models of light propagation to provide a ne technique for depth-resolved ILSI in NC plaques in vivo. By combining circumferential scanning with depth-resolved ILSf the complete 3D determination of plaque viscoeiasticity distributions can be achieved at selected sites to furnish information on the load bearing properties of the lipid-pool and fibrous cap. Because ILSI measurements can be based on phase shifts of multiply scattered light caused by minute scatterer displacements, this exemplary technique can be highly sensitive to small changes in plaque viscoelasttc properties, and can render high precision for the evaluation of lipid pools.
[00102] It can also be possible to provide a clinical translation of the ILSI technology according to an exemplary embodiment of the present disclosure.
Ex mplary Approach
[ΘΘΪΘ3) i xetnplary Overview of ILSi can be based on an exemplary laser speckle approach thai has been deveioped to evaluate the viseoelastic properties of tissue. (See, e.g.. References 42-45, 65 and 68-70). For example, laser speckle (e.g., FIG. 2), (see, e.g.. Reference 71) a grainy intensity pattern that occurs by the interference of coherent light scattered from tissue, can be modulated by the Brownian motions of endogenous particles within tissue. In can be well known that the extent of particular Brownian motion can be intimately related with the micromechanical susceptibility of the medium, and particles can exhibit larger motions when their local environment can be less viscous. (See, e.g.,
References 72-74). Consequently, in an atherom due to the low viscosity of lipid, scatterers can exhibit rapid Brow nian motions, eliciting rapid speckle intensity fluctuations compared to stiffer fibrous regions. The extent of speckle fluctuations can be quantified from the speckle decorrelatioii curve, g2(t), which can be obtained by calculating the normalized cross- correlation coefficient over a time series of laser speckle patterns (e.g., FIG. 3). The rate of speckle modulation given by die speckle decorrelation time constant T can provide a highly precise index of plaque viscoeiasiiciry that ca be closely related to plaque composition and viseoelastic moduli. (See, e.g.. References 42, 6)
( 104) (i) Plaque Chamcierimtiofi: Studies have been conducted to demonstrate the capability of LSI for evaluating the index of viseoelasticity, r, in cadaveric plaques.'1'' Time- varying speckle patterns were obtained from approximately 100 arterial samples using a Helium Neon source (e. g., 632 rsm) and a CMOS camera to evaluate girt). The time constant, Γ was measured by exponential fitting of g2(() for each plaque (e.g., FIG. 4).
Exemplary results show that t cait provide highly sensitive discrimination of plaque type (e.g., p<O.00I). In particular, LSI can demonstrate exquisite sensitivity (e.g., 100%) and specificity (e.g., 92%) for discriminating the viseoelastic properties ofTCFAs (e.g., r ;;; 45ms) due to rapid particular Brownian motion within low viscosity lipid pool (e.g., p<0,Q0Qi ). Similaiiy, stiffer fibrous and fibrocalcific lesicms can elicit significantly larger T values (e.g., p : 0 001 >
{'00105] (ii) Exemplary Relationship between τ and, plaque composition; Since viscoelastic properties can be highly dependent on composition, t showed high correlation with plaque collagen content (e.g., ::::0.73; pO.OOGl) and consequently with cap thickness (e.g., R:::0.87; p<O.O01). (See, e.g.. Reference 42). Given the low viscosity of lipid, a strong negative correlation (e.g., R--0.8 ! ; pO.OOOl) between t and lipid content was observed. These exemplary results demonstrate that LSI can measure a index of viscoe!astieity, τ, closely related with compositional metrics associated with plaque stability.
[00106! (H- T . a d, yjscoe last c mo dn ius : In order to demonstrate the potential of LSI in estimating the viscoelastic properties of samples, the relationship has been eval uated between the modulus of viscoelasilcity, G, measured by mechanical testing and LSI time constant, t , usin (a) homogenous gels, and (b) atherosclerotic plaques.. In the present disclosure, the term, 'bulk' modulus, G, can be used to define the overall modulus which integrates over the sample volume.
[00107] (a) Homogenous feels: LSI was performed on collagen, PDMS, PEG and Matrigel substrates of varying concentrations. Corresponding mechanic al testing measurements were performed on all samples using a strain-controlled rheorneter (e.g., ARG2, TA Instruments Inc., MA) to measure modulus G. The samples were loaded between the parallel plates of the rheorneter and an oscillatory strain (e.g., 1 %) was applied over a frequency range of about 0.1 -5Hz. High correlation between T and G (e.g., R:::0. 2, p<0.0Ol) was observed over the linear frequency range in ail samples. These results confirm that t can provide a highly accurate estimate of sample viscoelastic properties (e.g., FIG. 5). To evaluate the
measurement sensitivity of .LSI, time lapse measurements of X were compared with G values measured during slow coring of PDMS gels over 24 hours. High correspondence between X and G was observed (e.g., R ~ 0.95, p<0.01 ), confirming the high sensitivity of the LSI approach to changes in viscoelastic properties of the sample. "3
[
'00108]
LSI was conducted by averaging r alues over 3mm disks of aortic sites, histologically confirmed as calcific, fibrous and lipid~.nch. Mechanical testing was performed as above, which revealed distinct G values between plaque groups: 2.27x1
(f Pa (calcific), 3.65x10 Pa (fibrous) and. 2.23x10 Pa (NCFA). Analysis of variance ("ANOVA") tests showed statistically significant differences in G for the plaque types (e.g., pO.OOi .). These values also correspond with previously published reports. (See, e.g.,
Reference 76), For all plaques, t correlated well with G (11=0.97, p<0.G01), establishing the close relationship between f and plaques viscoelastie properties, and suggesting that r
cail pro vide a key metric related to the mechanical strength of the plaque.
Exemplary Influence of spatial heterogeneities:
5 [00109] (i) Modeling Studies: To evaluate the influence of structural parameters on the bulk modulus, a plaque was modeled as a tnultilayered cylinder of thickness, L and viscoelastie modulus, G, For the purpose of this model, it can he assumed that viscoelastie modulus, O* G" (elastic modulus), supported by previous reports; (See, e.g.., References 76 and 77).
Given its clinical significance, it can be possible to consider a NC plaque with a fibrous cap 10 and NC of thicknesses Li and L2_ and moduli Gl and G2, loaded between the parallel plates of a rheometer. The twisting moment M applied by the rheometer can be determined by the distribution of shear stresses integrated across the plaque. (See, e.g.. Reference 50). By equating M with the polar moment of inertia and the angular displacement of the sample, i can be possible to deduce the expression, for example; i t) G *—···-· ( I )
L7G + L^Gf
} 00110] Eqn. (1) shows that the overall bulk modulus of the plaque can be related to the thickness and viscoelastie modulus of each layer. Eqn. (2) below can be applied to evaluate the relationship between the bulk modulus G and fibrous cap thickness in a NC plaque, using previously reported values (see, e.g.. Reference 76) of G l - 496 kPa, and G2 - 222 kPa, for 0 fibrous and lipid rich tissue and can evaluate the influence of varying fibr s cap thickness (e.g., 0-500 μιη) on bulk G (e.g., FIG.3), This exemplary model can be also extended to include multiple layers of varying depth-dependent viscoelasticity by using the following exemplary generalized equatio
5 [001.11 J These studies indicate that the fibrous cap thickness can greatly influences the overall bulk viscoelasticity of the plaque (e.g., FIG. 6), and also indicate that the
measurement of bulk viscoelastie properties can provide a key metric closely related with plaque stability.
f Ml ί 2} (it) Lateral scanning in LSI: To evaluate the capabilit of LSI in measuring heterogeneities, laser speckle images of plaques were obtained by scanning the He Ne spot at 300μηι increments and the spatial distribution of - was measured. Figure 7 demonstrates the lateral variation of fas a function of beam location. As the beam was scanned across each lesion, T varied significantly depending on tissue type: X was low (e.g., 20 - 50 m$) in the low viscosity NC regions (e.g., FIG J A) and higher in the stiffer calcific (e.g., approximately 2200 ms in FIG.7B) and fibrous (e.g., approximately 800ms in FIG.7C) regions. Similarly 2D maps of the spatial τ distributions were obtained by beam scanning over the region of interest ("RQI") (e.g., FIG. 1), to facilitate a detection of heterogeneities such as calci fic nodules and lipid pools to facilitate comprehensive coronary screening.
{'06.1.13] (Hi) Depth-dependent heterogeneities: Due to the diffusive properties of light propagation in tissue, photons returning from deeper regions have a higher probability of remittance farther away from the illumination location.''**"*1 In such publications, f was computed over the entire speckle pattern. Therefore, Brownian motion was integrated over all optical depths and information about tissue heterogeneity was lost. By combining LSI with an exemplary Monte Carlo analysis of light propagation, depth information can be obtained. (See, e.g.. Reference 43). In this study, the capabilit to measure fibrous cap thickness was demonstrated by analyzing variation in T as a function of radial distance, p, from the illumination location in each speckle image. Fibrous cap thickness estimates obtained using this method, were highly correlated with Histologic measurements (e.g., FIG. 8). These findings indicate the potential of obtaining depth information using LSI, which can be further explored for in vivo intracoronary use in the current proposal.
Exemplary Intracoronary 1LSI
Exemplary 1LS1 catheter construction and testing
{00114| (i) Fiber bundle selection: Optical fiber bundles form an important part of the exemplary ILSI catheter to transmit speckle patterns. One challenge can be that speckle modulation can be influenced by inter-fiber light leakage (e.g., cross-talk) which can likely be exacerbated during motion. A study*4 was performed to investigate the influence of motion on the diagnostic efficacy of fiber bundle based LSI in 75 at tenal plaques, whi le cyclically modulating the flexible length of the bundle to mimic cardiac motion and tortuosity. A variety of fiber undles were tested. The bundle with the highest motion tolerance was
selected as having ihe (a) highest correlation, (b) lowest error, and (c) minimal statistically significant difference in measuring plaque τ values under stationary and moving conditions.
Low cross-talk leached fiber bundles provided the best motion stability (e.g., SCHOTT, inc.), likely due to the manufacturing (e.g., leaching) process which can result in large separations between fiber cores and reduced cross-talk. (See, e.g.. Reference 44). In particular, the leached bundle with the smallest partial core size of approximately 0.36 (e.g., core area ÷ fiber area) provided the best results for the above three criteria (e.g., FIG. 9), Based on these findings, miniaturized leached fiber bundles with low partial core sizes (e.g., <0.4) can be incorporated in the clinical-grade ILSl catheter proposed in this grant.
J00115J (u) Exemplary !LSJ catheter: An exemplary ILSl catheter (e.g.. dia::::1.57mm; can be provided that can include an inner optical core and custom-designed external sheath.46 The optical core (e.g., FIG. 10) can consist of an optical fiber to illuminate the arterial wall and a leached optical fiber bundle to collect arterial speck le patterns. The exemplary desig of the catheter distal optics for light delivery and speckle image transmission was optimized using ZEMAX (e.g., ZEMAX Development Corporation) for an approximate 500 um field of view ("FOV"). The optical elements (e.g., GRIN lens, polarizer and mirror) were assembled at the distal face of the fiber bundle within a clear tube (e.g., FIG. 10) and the proximal bundle face was imaged via an objective lens and CMOS camera. To house the optical core, it can be possible to use a double-lumen catheter sheath (e.g., FIG. 1 1 ). Since blood presents an impediment to intracoronary optical approaches, the sheath can include an occiiision balloon which can facilitate the comparison of the effectiveness of proximal balloon occlusion
("PBO") with flushing techniques during the exemplary ITS procedure. The sheath can also have radio-opaque marker at the distal end for fluoroscopic guidance and a rapid exchange guidewire port The catheter performance in evaluating cadaveric plaques was compared with free-space LSI: high correlation (e.g., R = 0.79, p<0.0l) was attained between ILSl and free space Ί measurements. For in vivo testing, the catheter was interfaced with a portable console for intravascular evaluation in the aorta of a living rabbit. Distinct differences in arterial f measured at normal aortic and stented sites confirmed in vivo feasibility. (See e.g., Reference 46).
Exemplary Fea sibi j t y of intra co ronary LSI in living swine
|0β 16) The feasibility of ILSI has been reviewed for coronary evaluation i vivo, and to determine the influence of cardiac motion, and blood displacement approaches.
[00117| (i) Human to swine coronary xenograft model: Exemplary choice of animal model can be motivated by two key requirements: (1) feasibility of ILSI can be best tested on human coronary disease, and (2) testing must be performed under conditions that mimic human cardiac physiology. This model has been previously described to test intraeoronary optical technologies. (See. e.g.. References 82 and 83). Cadaveric hearts (e.g., N-3) from patients who died of AMI were obtained (e.g., DRI), L AD and RCA coronaries (e.g., proximal 5 cm) were prosecuted and side branches Ugated. Coronar grafts were marked with India ink on the adventitial side to identify discrete sites for co-registration with Histopathology. hi anesthetized swine (e.g., N-3), the chest was opened, the grafts were sutured on the heating swine heart, and blood flow was redirected through the graft via an aorto-atrtal conduit. A total of 24 discrete sites in 6 grafts were evaluated using [LSI in 3 living swine.
{00118] (ii) Exemplary ILSI procedure: A portable console was developed, which incorporated a Helium-Neon source (e.g., 632nm, 30mW) and a CMOS camera to capture speckle images at frame rate approximately I kHz (e.g., 512x512 pixels). The ILSI catheter was manually advanced under fluoroscopic guidance over a guide wire via the left carotid and to each discrete lesion by co-registering the illumination spot with the visible India ink mark on the artery. Prior to imaging, the proximal occlusion balloon was engaged while flushing with Laciated Ringers (";LR") to ensure that blood did not re-enter the FOV. To evaluate the influence of cardiac motion, acquisition of the first speckle image was triggered on the R-wave of the swine EKG signal, followed by asynchronous acquisition of subsequent, frames over approximately 5 cardiac cycles. Following the exemplary ILSI procedure, the swine were sacrificed, and the grafts explanted and processed for Histopathological evaluation. Plaques (e.g., N::::24 were diagnosed as lipid pool (e.g., n~3), pathological intima! thickening ("PIT") (e.g., n~7) and fibrous (e.g., n~14) plaques. (See, e.g., Reference 1 ). I LSI analysis was performed as detailed below.
{00.119] iiii) Influence of c rdiac moti on : In order to achieve clinical viability in patients, the ILSI technology can. facilitate rapid coronary screening while retaining adequate motion stability over the cardiac cycle. While EKG gating can be implemented to mitigate the influence of cardiac motion, this approach can add significant time to the imaging procedure.
Instead, a non-gated approach can permit rapid imaging of long coronary segments facilitating the use of the ILSJ device in patients. The studies below were performed to investigate the influence of cardiac motion and compare EKG-gated versus non-gated ILSI measurements.
(001201 To evaluate the EKG gating approach to conduct ILSI, the t value for each plaque was calculated at the mid-diastole phase of the cardiac cycle (e.g., approximately 600 ms after onset of R-wave). To evaluate the non- gated approach, the T value for eac plaque was computed at a time point during the cardiac cycle that was randomly selected by software. For both cases, t was calculated by exponential .fitting of 50 ms of the initial decorrelation of the g2(t) curve. FIG. .13 A shows an exemplary illustration of the average r computed for the plaque groups using the EKG-gated and non-gated approaches, and the results of the pairwise comparisons between plaque groups are shown in FIG. 13B. Using both approaches, differences in f between the three plaque groups were highly significant. Demonstrating that plaque viscoelasticity could be well distinguished even under conditions of cardiac motion. This can be because sufficient motion stability can be achieved by employing rapid image acquisition rates (e.g., l kBz) using a high speed CMOS detector to measure laser speckle fluctuations over very short time scales. In vivo plaque time constants were about <25ms (e.g., FIG. 1.3), indicating that imaging durations of about. 25 ms can sufficiently enable plaque discrimination. Given the low frequency of cardiac motion (e.g.,
approximately 1 Hz) relative to the high rate of speckle decorrelation over short time scales, ILSI can be conducted during the cardiac cycle without: the need for EKG gating. A key result can be that differences in t measured within the same plaque group using the two exemplary approaches were not significantly different (e.g., FIG. 13 A). This can
demonstrate that non-gated ILSI works just as well as EKG-gated ILSI in vivo. From the results of this study, it can be possible to infer that: (a) an imaging duration <25ms can be sufficient to measure speckle decorrelation for plaque evaluation in vivo, and (b) ILSI can be conducted in vivo without E G-gating.
[00121 j (iv) Intracoronary flushing: Similar to other intravascular optical techniques, in ILSI the presence of blood can hinder imaging. Proximal balloon occlusion ("PBO") and purging with flushing media can be two exemplary methods routinely used in conjunction with angioscopy and OCT to displace blood during the imaging procedure. (See, e.g.,
"?7
References 19 and 61 ). While PBO ca routinely be employed in Japan, the risk of ST~ segment elevation can limit the widespread adoption of this method in the USA. Instead, flushing with contrast agent (e.g., Visipaque) or Lactated Ringers ("LR") solution can routinely be used as a safe alternative during imaging."4 Therefore, in order to assess feasibility of ILSI for patient use, studies were conducted in native coronaries of living swine to compare PBO and flushing approaches as detailed below.
[00122| (a) Balloon occlusion versus flushing in living swine: A 3mm coronary stent was deployed into the nati ve LAD of anesthetized swine, and LSI was conducted at normal arterial sites, and within the stent, while the proximal occlusion balloon was engaged. The balloon was then disengaged, and the sites were evaluated, in conjunction with a 30 cc Visipaque flush. Using both PBO and flushing approaches, differences in τ betwee the normal unstented and stented sites were highly significant (e.g., p<0.01), demonstrating that ILSI can be conducted using either of the two exemplary approaches to displace blood during imaging. In addition, differences in τ measured within the same location with both PBO and flushing were not significantly different (e.g., FIG. 1.4). This can demonstrate that LSI can be conducted in conjunction with flushing to sufficiently displace blood during imaging.
[00123] (b) jrrfinence of residual blood: To test the influence of residual blood cells on r values, tt can be possible to perform LSI on tour aortic plaques within a flow ceil through a 3 mm intervening layer of whole blood (e.g. , HCT-30%), seriall di luted using PBS. For example, t values at HCT < 0, 1 % were similar to those values measured without any intervening medium. Subsequently OCT imaging was performed and it was determined that at HCT >0.03%. hackscattering from blood cells was clearl evident in OCT images, in clinical studies using intracoronary OCT and recent swine studies (e.g., FIG. 1 } no hackscattering from blood cells can be observed during flushing. Since blood does not affect LSI at a HCT <0.!% and purging in patients can apparently reduce the intracoronary HCT to <O.03%, levels of residual blood cells during flushing can be sufficiently low to conduct ILSI.
[00124] Summary of Exemplary Studies: Through certain studies, the exemplary LSI systems and methods have been developed and validated as a powerful tool to evaluate plaque viscoelastic properties. These exemplary studies have demonstrated, 1) The LSI time constant, ,. can provide a metric that can intimately be linked with plaque viscoelastic
properties, 2) LSI can enable highly precise differentiation of plaque type, and can have exquisite sensitivity for the evaluation of TCFAs. 3} LSI can facilitate the measurement of spatial and depth-dependent heterogeneities, 4} Intracoronary LSI can be conducted in vivo at high imaging rates in conjunction with flushing. Given the high clinical impact of measuring coronary plaque viscoeiasticity and supported by the success of exemplary results in the current disclosure, it can be possible to extend LSI for intracoronary evaluation in patients. It can also be possible to provide, according to an exemplary embodiment of the present disclosure, clinical grade ILSI technology, and conduct the first I human feasibility studies as detailed below,
Exemplary Design and Methods
[Q0125| Overview of Exemplar design: Efforts have been directed towards developing clinical-grade ILSi catheters suitable for human use and a console to enable helical scanning over long coronary segments. Preclinical validation of the new ILSi device can be conducted to evaluate coronary plaque viscoeiasticity in living swine. Further, for human clinical studies can be conducted, for example, in 20 patients to assess the safety and utility of 1LSL It can also he possible to obtain an exemplary tool that can improve an understanding of human CAD,
Exemplary Methods:
[00.t26| The exemplars' ILSI catheter described in exemplary studies above enabled the demonstration of hi vivo feasibility for intracoronary evaluation. Its functionality for patient use, however, can be restricted given its large size (e.g., approximately 4.5F/ L57mm), in addition, the existing ILSI devices may only be permit limited point, sampling of discrete sites, therefore precluding the capability for comprehensive intracoronary screening to evaluate arterial viscoeiasticity distributions. Furthermore, because the exemplary device can utilize illumination over an extended beam (e.g., approximately 250 pm), and the index of viscoeiasticity, t , evaluated over the entire speckle pattern, depth-dependent information can be lost or degraded. These issues can be solved according to certain exemplary embodiments described herein below.
[001271 In order to achieve clinical utility, for example, a miniaturized exemplary ILSI catheter (e.g. , approximately 2.4F- 3.OF / 0.8- 1.0mm) can be provided that can access small flow-limiting coronary arteries of patients, and can conduct rapid helical scanning of
coronary segments. Speckle analysis, and visualization methods can be implemented to reconstruct arterial viscoelasticity distributions. This can facilitate comprehensive circumferential screening of about 3,0-5.0 cm of the major coronary arteries with a longitudinal image spacing (e.g., pitch) of about 0.25 .1.Omia, while administering a safe total amount (e.g., < 100 cc) of flushing media.
(60128) Exemplary /LSI device: Exemplary modifications of the exemplary device can. be focused on certain components thereof, for example: (i) catheter, (ii) motor drive assembly for helical scanning, and (iii) console. The catheter can. include an inner cable that can. boose the optical core. During imaging, the motor drive assembly can rotate and simultaneously pul!back the inner cable within an outer stationary sheath to accomplish helical scanning (e.g., FIG.1.5).
[001291 Exemplary // ,S7 catheter: It can be possible to provide a miniaturized leached fiber bundles (e.g., diameter approximately 250 p.m, length :::: 1m). Utilizing a fiber size of approximately u.ri, with a partial core area of approximately 0.4, approximately 2000 collection fibers can be incorporated to obtain a fiber bundle with sufficiently low cross-talk to transmit speckle patterns. A central light delivery fiber can be included for illumination. Micro-optical components including a focusing lens, custom polarizer and rod mirror can be optimized, tested and affixed to the distal bundle face. A variety of different lenses can be investigated, including GRIN lenses and custom- fabricated ball lenses, and optimized to provide a focused illumination spot size of approximately 20μηι and imaging FOV of approximately 500μηι. Miniaturization and fabrication of optical components can be conducted to achieve a target optical core size of approximately 300μηι. The optical core can be affixed within a driveshaft cable {e.g., Asahi Intec, CA) to convey torque from a motor to enable helical scanning. A transparent rapid-exchange sheath with a guide wire port can house the catheter cable assembly, and can be tested for optical clarity.
[001301 ( i t) Exemplary motor drive assembly can include an optical rotary junction.
C RJ'") that can couple light with the rotating optical core (e.g., FIG. 15). Excellent, rotational uniformity (e.g., < 10% modulation) and low transmission loss (e.g.. < 1 dB) in can be provided with ORJs provided in the exemplary OCT/OFDI systems. (See, e.g.,
References 60, 61 and 85). For example, the ORJ was designed to couple with single optical fiber within the OCT/OFDI catheter while continuously spinning at speeds of
approximately 6000 rpm. The ORJ can be provided for the use with the exemplary ILSI device such that, (a) it can facilitate coupling of light with a rotating optical fiber bundle consisting of multiple optical fibers, and (b) the exemplary ILSI catheter would, not spin continuously, instead in order to permit acquisition of the speckle image time series over about 25ms at each circumferential location (e.g... based on studies shown in FiCl 13), a stepper motor can be incorporated to rotate the optical core at discrete steps with a residence time of about 25 ms per step. The exemplary ORJ can include a collimating lens (e.g., L2) affixed at the proximal end of the optical core and a motor coupled with the drives haft to enable rotation. A CMOS sensor (e.g., Mikrotron 1310} can be housed directly within the ORJ, and transmitted speckle patterns can be imaged, via a stationary lens (e.g., LI). The rotational rate of the catheter can be 1 Hz. A linear puliback stage can facilitate a
translation puUback during imaging over speeds of about 0.25-1.0mm s. Rotational distortion (e.g., < 10%) can be measured by comparing T values of aortic plaques with a stationary catheter (e.g., Table.1)
' k J ®mn tmt emh & i fw t mce in Aim I
[00131! (ii-i) ¾e p or table . console can be modified to facilitate helical imaging and data visualization. Engineering tasks can include: a) interface to control the motor drive assembly and automated flush devices, and b) software interlace design. Similar to the exemplary preliminary studies, a He Ne light source (e.g.. 632 nm, 30 mW) can be used for illumination. Time-varying laser speckle images can be collected at an approximately 2kHz frame rate (e.g., 512x512 pixels).
(001321 Reconstruction of arterial visco lastirity maps: To obtain sufficient spatial sampling during catheter rotation, a lateral spacing of about <250 pm can be utilized between
rotational steps. Considering the typical coronary circumference of about J.0mra, and the catheter FOV of about 500 μ.∞, 40 discrete steps can facilitate adequate spatial overlap for sufficient circumferential sampling at about a lHz rotational rate. The longitudinal scan pitch and total imaging time can be determined by the pull-back speed (e.g.. Table 1 ).
jJOO.133} Exemplary 2D reconstruction: To evaluate 2D arterial viscoelasticity maps, at each site, f can be computed over each speckle image by exponential fitting of the g2(t) curve using previously reported techniques. (See e.g., References 42 and 46). The resulting 2D army of discrete r values can be processed using spatial filtering and bilinear image interpolation approaches to reconstruct maps corresponding to arterial viscoelasticity distributions.^' NC plaques of high clinical relevance identified by low t values (e.g., approximately 5- 10ms) can be selected to explore depth-resolved analysis.
[00134} Exemplary Depth analysis: The capability of ILS1 to provide 3D depth-resolved distribution of t values in NC plaques in vivo is described below. For example, at each location (x,y) over the FOV, windowed cross-correlation can be performed over the speckle time series to obtain g2(t). To ensure sufficient ensemble averaging, g2(t) can be measured by averaging several cross-correlation functions that evolve in time over about a 25 ms imaging duration and over neighboring pixels, which can influence the measured spatial resolution for mapping. The resulting 2D distribution of (x,y) can be obtained (e.g., FIG.16) by exponential fitting of g2(t) curves. Due to light transport properties, f (x,y) farther from the beam location can be influenced by longer optical paths. Using a Monte- Carlo Ray Tracing ("MCRT") algorithm, a look up table of the 3D distribution of mean penetration depths (z) over the FOV remittance plane can be created {e.g., F G.17), and the corresponding depths for each t (x.y) can be determined to provide the depth-resol ved distribution of T . The process can be repeated at each circumferential beam location to reconstruct, the full 3D viscoelasticity distribution of NC plaques.
[00.135} Estimated Resolution: Axial resolution can be estimated by the MS width-half maximum ("FWH " ) of the penetration depth distribution and the lateral resolu tion can be determined from the FWHM of the radial scattering PDF. Estimated values using MCRT can be plotted (e.g., FIG, 18). Spatial resolution can degrade with depth (e.g., Table I ).
However, over superficial depths, the estimated spatial resolution about <1 ΟΟμιη can be sufficient to evaluate thin caps that can be most clinically relevant. At. deeper depths (abou
>1 0μπι). resolution approximately about !00-200(um can be sufficient to evaluate large necrotic cores of highest, significance. Exemplary methods described herein can be tested on human arteries and phantoms of spatial and depth-varying properties. Axial resolution can be measured by scanning a sample of known G within scattering media using a motorized stage, Lateral resolution can be verified using a patterned PDMS resolution target*'""89 Utilizing exemplary beam scanning in conjunction with depth-resolved LSI can provide an important understanding of the viscoelastic properties of the fibrous cap and NC layers to estimate the load bearing capabilities of clinically significant NC plaques.
} 0θί 36| n f Testing and miidatioft in swine: The human to swim coronary xenograft model (e.g., preliminary studies) can be used to validate the ILSI device for coronary screening. Human coronary grafts (e.g., 2 per heart x 10 hearts) can be grafted in
anaesthetized swine (e.g., -iO) for ILSI validation. The distal start and. end of scan locations can be .marked by India ink corresponding with the visible ILSI beam for co- registration with Histology, Scanning can be performed over an approximately 5cm pull- back in conjunction with a Visipaque flush. Following ILSI, the grafts can be evaluated using intracoronary OFDI in vivo. Histology sections can be obtained at 2mm increments and co-registered with the corresponding ILSI cross-section. For example, a total of 500 ILSI-OFDI-Histology correlated cross-sections can be analyzed (e.g., 25 sections/artery x 2 arteries x 10 hearts). Plaque type can be diagnosed at approximately 250um spacing using both Histology and OFDI as, for example, TCFA, HFA, PIT, Fibrous or fibrocalcific, and compared with t at each site. In NC plaques, fibrous cap thickness can be measured by depth-resolved ILSI and can be compared with Histology. Success can be determined by ANOVA tests to evaluate X difference between groups, based on OFDI and Histology diagnosis, p<(U)5 can be considered statistically significan
Exemplary Alternate Embodiments:
[0β137[ Exemplary optical rotary junction: in the unlikely event that about > 10% deviation in τ can be observed during catheter rotation, an alternative approach (e.g., recently
demonstrated in OCT) (see e.g.. Reference 90) can be implemented in which the optical core can be maintained stationary, and torque can be conveyed to the distal mirror via a dri es!iaft. It can also be possible to the use of cone mirrors to conduct LSI. Assuming that cone mirrors
ccaann bbee s suuffifficciieennttllyy mmiinniiaattuuririzzeedd.,, tthheeyy ccaann l liikkeellyy p prroovviiddee a a. vviiaabbllee o oppttiioonn t too eennaabbllee
oommnniiddiirreeccttiioonnaall v viieewwiinngg i inn tthhee IILLSSII ccaatthheetteerr..
[ {0000113388[| EExxeemmppllaarryy AA UUeerrnnaaiiee eexxeemmppllaarryy ccaatthheetteerr ddeessiiggnnss;; EExxeemmppllaarryy IILLSSII p prroocceedduurreess ccaann bbee ccoonndduucctteedd iinn tthhee ccoonnjjuunnccttiioonn w wiithh ssaalliinnee fflluusshhiinngg.. IInn t thhee uunnlliikkeellyy eevveenntt tthhaatt ssaalliinnee 55 f flluusshhiinngg ddooeess nnoott ssuufffificciieennttllyy ddiissppllaaccee bblloooodd,, aa m muullttii--pprraonngg ccoonnttaacctt bbaasseedd ddeessiiggnn ccaann bbee
eemmppllooyyeedd tthhaatt ccaann . mmaaiinnttaaiinn eeooddoolluummiinnaall ssuurrffaaccee ccoonnttaacctt dduurriinngg iimmaaggiinngg.. SSiimmiillaarr ccoonnttaacctt bbaasseedd ccaatthheetteerrss c caann b bee uuttiilliizzeedd i inn tthheerrmmooggrraapphhyy s sttuuddiieess aanndd ccaann hhee aapppprroovveedd foforr uussee i inn ppaattiieennttss.. ((SSeeee,, ee..gg.. RReeffeerreennccee 2244))..
{ {0000.113399]] EExxeemmppllaarryy DDeepptthh--aannaallyyssiiss:: TThhee iinn vviivvoo ffeeaassiibbiilliittyy ooff 33DD aannaallyyssiiss ccaann bbee
1100 ppeerrffoorrmmeedd,, aanndd tthhee ppeerrffoorrmmaannccee mmeettrriiccss ccaann bbee bbaasseedd oonn 22DD mmaappss ooff bbuullkk rr mmeeaassuurreemmeennttss,, bbaasseedd oonn tthhee rreessuullttss ooff pprerevviioouuss e exxeemmppllaarryy ssttuuddiieess tthhaatt eessttaabblliisshh tthhee ssiiggnniifificcaannccee ooff bbuullkk rr fo forr aasssseessssiinngg hhiigghh--rriisskk ppllaaqquueess..
J 014 1 Optimal flushing parameters to conduct ILSI in patients; ILSI can be
conducted in vivo while flushing with Visipaque to displace blood. To calculate the total imaging duration over which clear viewing of the arterial wall can be achieved, further studies have been conducted in living swine. Flushing with Visipaque was performed at flow rates of about 2-4cc s, commonly used in patients, and OFDI was simultaneously conducted to evaluate blood scattering within the lumen. For a single 10 cc flush at about 3cc/s, optimal blood clearance and unobstructed viewing of the arterial wall was achieved over
approximately 6 s (e.g., FIG.19). From these results, it was inferred that to conduct ILSI in patients, 8 intermittent flushes (e.g., lOcc/flush) of Visipaque can facilitate sufficient blood displacement to scan an approximately 5cm long coronary segment in less than a minute (e.g., at a scan pitch ::: lmm). Thus a low total volume of approximately 80cc of Visipaque can he administered. The average volume safely administered in patients can be reported to be about 265^130ml. (See, e.g.. References 92 and 93).
{'0014.1] Mmn&tt ILSI study. F ollowing regulatory approval, ft can he possible to evaluate coronary plaque viscoeiasticity using ILSI in a cohort, of 20 patients with native CAD who present at the GH cardiac catheterization laboratory for percutaneous coronary intervention
("'PCF"). In order to test the feasibility of the ILSI approach in patients, intraeoronary OFDI can be used to provide a nncrostruetura! context for ILSI results. Briefly, the culprit lesion can be determined from the patient's angiogram. The OFDI catheter can be advanced over a guide wire just distal to the culprit lesion. The maximum coronary length scanned can be about 5.0cm (e.g., range; 2. OS.0cm, imaging/flush parameters calculated below are based on maximum length). During a 3s. 3cc/s flush, the OFDI catheter can be withdrawn at a pullback speed of about 20 mm s to scan a 5cm segment. Following the OFDI procedure, ILSI can be conducted. The ILSI catheter can be similarly advanced distal to the culprit lesion under fluoroscopic guidance. Safety can be evaluated by monitoring hemodynamic parameters, E G and development of symptoms during the exemplary ILSI procedure. The ILSI catheter's rotational rate can be about 1.0Hz and imaging can be conducted in conjunction with 8 intermittent flushes (e.g., lOcc) at about 3cc/s as detailed above to image a matching 5.0 cm length in <50s. The total amount of Visipaque administered for the entire imaging procedure can be <100ce. It can be expected that the exemplary procedure can add 15 - 20 minutes to the routine PCI procedure (e.g., typical duration of 120 minutes),
[001421 i «ta co-registration and analysis: To determine the feasibility of ILSI in patients, ILSI 2D viscoelasticity maps can be compared with plaque type and microstructural information obtained from OFDI, in order to accomplish accurate data comparisons, digital coronary angiography can be conducted at the start and end of both OFDI and ILSI procedures to permit data co-registration. Additional landmarks, including the guiding catheter, stent edges and side-branch vessels can be used to improve registration accuracy, i 1 Co-registration in the circumferential direction can be done by reading the motor encoder positions on the OFDI and ILSI rotary junctions. OFDI images can be interpreted using previously established methods to characterize coronary plaques as: TCFA, THFA, PIT, Fibrous or fihroeakific. (See e.g.. References 12, 62, 94 and 95). ILSl-OFDI correlations can be evaluated using ANOVA tests to assess the feasibility of ILSI in measuring distinct τ values based on plaque type. The feasibility of measuring depth-resoived viscoelasticity can be evaluated in NC plaques by co-registering ILSI 2D cross-sectional maps of τ distributions with corresponding OFDI cross-sections.
Exemplary Potential Problems and Alternative Strategies■;
[00143) Blood: Blood in the FOV can cause rapid blurring of speckle due to moving blood cells. Real-time speckle analysis can be implemented and scan repeated if τ < 1 ms. An alternative solution to detect blood can be to incorporate simultaneous coronary viewing via the same catheter with a white light source and color camera.
[00144) Cardiac motion: For example, ILSI can be conducted without ERG gating. In the unlikely event that cardiac motion can be problematic, EKG gating can be utilized, and the feasibi lity of ILSI can be tested by evaluating discrete arterial sites predetermined by OFDI,
[00145) Nephrotoxicity: In. patients with renal impairment, .Lactated Ringers can. be used which has provided good ILSI results in exemplary studies. In these patients imaging can be restricted to a <3.0cm segment.
|00146| Culprit lesion; In the event that the culprit site can be inaccessible, OFDI and S LSI can be performed post-PCl.
[00147) OFDI: OFDI- ILSI comparisons can be verified. Since, no intracoronary technology exists to measure plaque viscoelasticity metrics In patients, in vivo ILSI feasibility can be tested using OFDI findings that have been well established for plaque evaluation. (See, e.g.. References 12, 62, 94 and 95),
[00148) The developed by Snyder (see, e.g.. Reference 120) can be applied to determine the various parameters of the fiber optic bundle as described herein, CM T can be an approximate analytical approach to study optical crosstalk between neighboring waveguides in terms of the coupling between guided modes of neighboring waveguides, to fully investigate coupling between all modes of adjacent fibers. The influence of multiple fiber bundle parameters on inter-fiber crosstalk and the modulation of transmitted laser speckles can be quantified.
Furthermore, fiber bimdle parameters can be defined to considerably reduce the modulation of transmitted speckle patterns caused by mode coupling between and within multi-mode cores.
Table 2. Specifications of two commercially available OFB
Ex mplary Devices And Methods- For Achieving Omul-Directional Viewing }0§149| I» an exemplary embodiment (e.g., FIG. 1.5) described and shown herein, the motor drive assembly can be used to conduct helical scanning of the vessel. The motor drive assembly can be modified to ach ie ve a 360-degree roiaiiori of the catheter, or it can be rota ted over a limited, or partial angle, to illuminate and image a section or sector of the lumen circumference at one time. The exemplary design can include a optical rotary junction ('ORJ") that can couple light with the rotating optical core. In OC T catheters, the ORJ can be designed to couple light with a single optical fiber while continuously spinning at speeds of approximately 6000 rpm_ The ORJ provided for the exemplary ITS! device can have two exemplary features: (i) it can facilitate coupling of light with a rotating catheter, and can include a fiber bundle with multiple optical fibers, and (ii) the ILSi catheter can be prevented from spinning continuously. In order to permit measurement of speckle decorrelation over about 25ms at each circumferential location, a motor drive can be incorporated to rotate the optical core at discrete steps with a residence time of about 25 ms per step. The ORJ can include a coMirnating lens at the proximal end of the optical core to couple light into a central illumination fiber, and a motor coupled with the driveshaft to enable rotation. An exemplary CMOS sensor can be housed directly within or connected to the ORJ, and can transmit speckle patterns imaged via a stationary lens. The exemplary rotational rate of the catheter can be about 1 Hz. A linear puilback stage can facilitate translation/pullback during imaging over speeds between, but not limited to, about 0.25-1.Onim s. In this exemplary design, the inner optical core can be affixed within a driveshaft cable to convey torque from a motor, to facilitate helical scanning. Some or all of the inner cable (e.g., the optical fiber bundle and distal optics) can rotate.
}00150| In another exemplary embodiment of the present disclosure, the inner optical core can remain stationary, and mechanical torque can be conveyed only to die distal mirror that ca be affixed to the driveshaft cable, in still another exemplary embodiment of the present disclosure, a ring of illumination fibers surrounding the collection bundle can be used to illuminate the tissue, and the distal mirror can be rotated. Via a ring of illumination fibers, the tissue can be illuminated using lisht with a sinele wavelength, o with, multiple .fibers .illuminating the tissue using different wavelengths of light. This can facilitate a better separation and a more robust analysis of speckle patterns. There can also be no are no
moving parts. Instead a multi-faceted .mirror (e.g., figures described below) can be
incorporated at the distal end for onraidirectional. viewing of the entire circumference of the lumen (e.g., 360 degree omnidirectional viewing). The multi-faceted mirror can be a cone mirror. Alternatively, or in addition, a cone-polygon / pyramidal shaped mirror can be used in which one or more of the reflecting surfaces can include one or more flattened reflective facets. Multiple illumination fibers can illuminate di fferent facets of surfaces of the multi- faceted mi ror, and speckle images can he collected simultaneously from 2 or more facets. During image processing, images obtained from multiple facets can be unwrapped and reconstructed to visualize the entire circumference of the luminal tissue of interest as shown below,
[0015JJ In another exemplary embodiment of the present disclosure, the optical core can remain stationary, and a rotating galvo-mirror can be incorporated at the distal end. The mi ror can be provided to fit within a less than about a I mm catheter sheath.
{'00.152) In some or all of the exemplary embodiments, an optional circular polarizer can be included to reduce the influence of back-reflections or specular reflections emanating f om surfaces of the catheter sheath, or from the surface of the tissue of interest. Specular reflections can be removed using software during post-processing of speckle images. This can be achieved by, for example, thresholding the image based on the temporal staiistics of speckle fluctuations where pixels with negligible speckle fluctuation can be masked out during analysis. This can ensure that only light, or other electromagnetic radiation, that has undergone multiple scattering can be analyzed to measure an index of tissue viseoelastkity. {'00.153') Preventing a receipt of the same polarization from .returning in the radiation (e.g. , light) can be beneficial in reducing back-reflected light of the similar polarization state that has scattered only once, or a few times, from the catheter surfaces and/or surface of tissue, which can otherwise increase the strong background intensity and confound the sensitivity of the device in measuring laser speckle intensity fluctuations scattered fro tissue. The polarizer can b replaced by computer software, or other methods, which can include spatial and temporal filtering that can similarly prevent back-reflections of light of the same polarization state, Filtering (e.g., to replace the polarizer) can be achieved by removing pixels in the image in which the intensity fluctuation can be zero, or negligible, over time caused by reflected light that has maintained its polarization state following a single or few
scattering events. Thus, fluctuating speckles causes by depolarized light, which has undergone multiple scattering through tissue, can be analyzed to measure the mechanical properties of tissue. xemp ary mage Processing And Visualisation
fufl154| Exemplary image processing procedures according to an exemplary embodiment of the present disclosure can include image unwrapping (e.g., FIG. 33B and 33C) removal of pixilation artifact (e.g., FIG. 33F), spatio-temporal analysis of speckle fluctuations and visualization using a time constant color map and display. Various procedures can be used depending on the complexity of measurement that is required.
[ 0I55| An exemplary procedure can include measuring measure the speckle decorreladon curve. g2(t }, by cross-correlation of multiple speckle frames obtained over the time series, conducting spatial and temporal averaging over multiple g2(t) curves and determining the time constant by exponential fitting over short time scales. The speckle time constant can be reported as an index of tissue viscoeiasticity. To extract 3D depth information, it can be possible to use a hybrid approach that combines Monte-Carlo ray tracing ("MCR'P) with spatio-temporal windowed analysis of speckle patterns. This exemplary procedure has been previously been validated in necrotic core coronary plaque, (See e.g., Reference 43). it can also be possible to modify this exemplary procedure to account for changes in the number of sc attering events as a function of depth.
{00156| Additional exemplary procedures can be provided to measure the elastic and viscous moduli of plaques directly from laser speckle patterns. g2(t) can be related to mean square displacement ("MSD") of !ight scattering particles within the plaque, and the MSD can be related to elastic and viscous moduli via the Stokes Einstein's formalisms, it can also be possible to display 2D maps of time constant by using spatial averaging, spatial filtering along with bilinear image interpolation techniques. It can be possible to further apply the above apparatus and methods for use in an exemplary helical/cylindrical display and for use co-registration for intra-coronary mapping. (See e.g., FIG. 1).
{00157} in order to validate the use of the exemplary LSI to measure tissue mechanical properties, LSI results of time constant (e.g., on test phantoms and tissue) can be compared with mechanical testing using a rheometer, which has been previously shown to exhibit excellent correlation in these studies (e.g., R
::::0.79, p<0.05). (See e.g., Reference 42).
f OISSI In order to validate capability of LSI to discriminate between plaque mechanical properties, LSI rime constants compared with Histopathological diagnosis of tissue type can be performed by a Pathologist. Differences between time constant measurements for different tissue types can be evaluated using ANOVA tests. Both ex vivo and in vivo studies show distinction can be good between NC plaques and other plaque types (e.g., including normal, fibrous, calcific and pathological intimal thickening). (See e.g.. References 46, 70 and 97). Since plaque mechanical properties can be dependent on collagen and lipid, correlation between time constant and collagen and lipid content within the measurement area of interest can be performed. Collagen content can be measured using Picrosinus staining, polarized Sight microscopy measurements and lipid using oil-red O, as well as
immunohistochemical staining to detect Apolipoprotehi B complex on LDL ciiolesterol, (See e.g., References 42, 45 and 95).
{'001591 Sensitivity and Specificity of the exemplary LSI has been measured previously in ex vivo validation studies. (See e.g., Reference 42). This can be done b receiver operating characteristic ("ROC") analysis. The exemplary test can evaluate the capability of LSI to distinguish mechanical properties of thin cap fibroatheroma ("TCFA") plaques as these can be considered more unstable plaques of clinical significanc e. The presence of TCFA can be considered - ve diagnosis, and all other tissue types can be considered -vediagnosis. Both sensitivity (e.g., 100%) and spec (e.g., 92%) can be maximized, which can be used with a diagnostic threshold of time constant of about 76ms. These exemplary studies can be similarly performed for in vivo studies.
{'001601 Sensitivity - True Positive/(True positive + False Negative) and specificity - True Negative/(False Positive +True Negative). Flushing using contrast agents, lactated ringers or dextran solution can be routinely used in the catheterization laboratory for OCT and angioseopy imaging. Additionally, a contrast agent can be routinely used for angiography for a number of years. No major technical challenges can be expected in the flushing process, and this mechanism can be manual or automated. It can be possible to provide procedures to
optimize flushing for ILSi (e.g., parameters: type of flushing agent, rate of flush, volume of flush, etc.) similar to OCT/angioscopy.
(001611 A practical challenge can potentially be inadequate flushing. Usually, the presence of blood can be easily detected as it can cause very rapid speckle decollation, and can provide a distinct time constant signature. To detect problems with inadequate flushing, it can be possible to include white light source to conduct color atigioscopy in tandem through the same catheter. Alternately, various other exemplary methods can be used (e.g., a dual wavelength illumination to measure absorption due to presence of blood).
(001 *21 i f flushing still poses a challenge., proximal balloon occlusion, can be used for a short period of time. Flushing for clearing blood from the field of view during optical imaging can be routinely employed in atigioscopy as well as and OCT/OFDL Over 1000 studies have been published, and this exemplary method is well accepted by clinicians.
Furthermore, flushing the coronary tree wi h contrast agent has been routinely used for many decades in conventional angiography procedures,
[00163] ILSI can be conducted, in vivo, while flushing with contrast agent or lactated ringers can be used to displace blood. The exemplary flushing mechanism is described in FIG. 36, Using calculations based on exemplary OCT studies in swine (e.g., FIG. 37) to conduct ILSI in patients, the use 8-10 intermittent flushes (e.g., !0ec/flush) of diluted contrast agent or lactated ringers Visipaque can permit sufficient blood displacement to scan an approximately 5cm long coronary segment in less than a minute (e.g., at a scan pitch ::: 1 mm). Thus, it can be expected that a low total volume about 80 - 100 cc of flushing agent can be administered during ILSI, which can be below the average volume that is safely administered in patients. (See e.g., References 92 and 93).
[00164] According to a further exemplary embodiment of the present disclosure, it. is possible to provide a miniaturized (e.g., <lmm) ILSi catheter that can be safel guided through the coronary artery to conduct intracoronary mapping. It can be beneficial to keep the exemplary device as similar to a commercially available (e.g., regulatory approved) 1VUS catheter and system as possible. It can also be possible to confirm ILSI catheter
characteristics (e.g., damage to endothelium, trackability, pushabiiity and ease of use) are similar to an exemplary IV US catheter.
Ex mplary Analysis Of Omni-Djrectional Mirror Configurations:
ίίθΙ 65| Exemplary embodiments of exemplary omni-directional catheters can include reflective arrangements or at least partially-reflective arrangement thai can include multiple facets at the distal tip of the catheter to direct electromagnetic radiation to die cylindrical. lumen, and to collect reflected speckle patterns from multiple sites of the lumen
circumference without rotating the catheter.
[00166| FIGS. 38A-38C illustrate an exemplary cone-poiygon/pyraniidal mirror for omnidirectional (e.g., laser speckle, etc.) imaging. The image is at the bottom of the image plane for the object that is at the top of the mirror. At the image plane, the centra! part can have more aberrations and a larger spot radius, while the edge can have less aberrations and smaller spot radius. The spot size at the edge can be smaller than a fiber's cross-section surface. Additionally; the off-axis object can cause overlap of the images if the off-axis object has an enough large distance.
{00.167} FIGS, 38E-38H illustrate an exemplar}* cone mirror-side view for vertical focal plane. The image is at bottom of the image plane for the object at the top of the mirror. At the image plane, the central part can have more aberrations more aberration and a larger spot radius, while the edge can have less aberrations and smaller spot radius. The spot size at the edge can be smaller than a fiber's cross-section surface. The off-axis object can cause overlap of the images if the off-axis object has an enough large distance. The horizontal aberration can be very strong due to curvature of the cone mirror.
[00168[ FIGS. 38I-38L illustrate an exemplary cone mirror top view tor horizontal focal plane. The vertical focal plane and horizontal focal plane can be at different location, (approximately 1 mm difference. Strong horizontal image aberrations can be seen, and can cause severe image overlap horizontally. Also present, is a big spot size, and an inadequate horizontal resolution.
[00169) FIGS. 39A-39H illustrate exemplary images obtained using various exemplary omni -directional mirror configurations. Exemplary selections of fiber bundle parameters can be used to reduce inter-fiber cross-talk during laser speckle imaging.
[00170J Optical fiber bundles can typically incorporate thousands- hexagonai!y arranged individual optical fiber cores as shown in FIG. 20A. The analysis of mode coupling between all of the fiber cores can be far too complicated and numerically intensive to be calculated.
However a simplified system of 7 parallel fibers can be used to model the coupling between the modes of these fibers (see, e.g.. References 1 15-117) and the result can be easily extended to an entire fiber bundle. Here, a multi-core optical fiber system of 7 hexagonally arranged cores embedded in a uniform cladding materia! as shown in FIG. 20B, can be used. The fiber bundle specifications can be based on two commercially available leached fiber bundles (e.g., SCHOTT North America) and are listed i n Table 2 above. These two types of fiber bundles were chosen because their specifications can be typical for the fiber bundles used in LSI (See, e.g., Reference 100).
}0θί71| Coupled mode theory (see, e.g.. References 1 14 and 120-122) can be a common theoretical model used to obtain approximate solutions to the coupling between waveguides of multiple waveguides systems. Compared to the normal mode expansion method (see, e.g.. Reference 1 15), in which the field can be expanded in terras of normal modes solved from Maxwell's equations with the boundary conditions of the entire complicated structure, in CMT the field can be decomposed into the modes of each individual waveguides (see, e.g., Reference 1 14):
Ε(.χ, ·.2} « Γ ο, ( ?)€,.{>·, *) ex.p "#.2)
Η<Α·, .>>, <-) « (2)hs. {.«,/) ex ' )
(3)
where av can be a complex amplitude of vth mode; ev and hv can be electric and magnetic components of normalized mode field of eac individual fiber, respectively; βμ can be the mode propagation constant of mode μ; z can be the propagation distance along the fiber bundle and the summation over v runs through all modes of all individual fibers. The effective refractive index of mode Er can be defined as neff
can be the wave number. For complex structure, the complete set of the normal modes can be difficult to solve out (see, e.g., Reference 1 .15) while in CMT, modes of each core of fiber bundle can be solved independently. The complex amplitude of modes can be obtained by solving the coupled mode equatio (see, e.g.. Reference 1 14) which can describe how the amplitude can vary with propagation distance z along with the length of the coupled waveguides, where, for example:
where ~~ *A* can be the mode coupling coefficient between mode v and μ. The coupling coefficient can be directly related to the degree of overlapping of mode field. The coupling coefficient kvp along with the difference in mode propagation constant
Αβνμ~~ βμ~
"βν d cide the coupling strength of μίΐτ, mode to vt mode. The mode coupling coefficient c
f[ can be determined by the overlap coefficient of mode fields (ev, h ) and (ep, hp) and the perturbation
1/1 of mode μ to the mode v. Here the element of matrix
and c
YU for the normalized mode field by definition.
The element of matrix x can be iven by,
where > cart be the angular frequency of the laser light and
can be the difference between the dielectric constant of the whole multi-core structure and the dielectric constant of the structure with only the individual fiber supporting the mode μ.
(00172| To evaluate the modulation of laser speckle patterns during transmission through the optical bundles,, laser speckle fields can first be numerically generated (see, e.g.,
Reference 97) by Fourier transform the field with random phase. The polarization of speckles can be chosen along with the linear polarization of fundamental modes of fibers. The generated speckle fields can then be decomposed into HE, EH, TE and TM fiber modes of individual fibers. The complex amplitude of each guided fiber mode at
~:iK av 0), ca be given by for example:
(5)
where Eo can be the generated speckle electric field. By solving the Eq. (4) for each propagating mode with the initial value of ay(O), the complex amplitude at propagation distance z can be obtained. The transmitted speckle patterns can then be reconstructed by linearly combining the fields of ail fiber modes with its amplitude. The modulation of the transmitted speckles can then be evaluated by the correlation coefficien t of the intensity
patterns between reconstructed speckle patterns at different positions along the length of bundles and the reconstructed speckle patterns at z::: 0 (see, e.g.. Reference 27), where, for example:
^ , ( /<*, y, ::} - / (?))(/(*, = - J(z - 0})
C r) * -^ — L
(6) where l(x,y , ) can be the intensity of speckle electric field
cm te Ae average and standard deviation of the intensity patterns at different z, respectively. Here x, y can be the transverse coordinates of the points within the 7 core areas. C - 1 can indicate that two speckle patterns can have same spatial fluctuations and so totally correlated when C
:::: 0, the speckles patterns can have no correlation. Thus the average of C over 20 speckle realizations can he used to measure the speckle modulation.
Exemplary Speckle Image Processing
Exemplary Elimination fPtxelation Artifact
[00173 To conduct ILSl, a small-diameter, flexible optical fiber bundle can be used to transmit the laser speckle patterns reflected from the coronary wall to the high speed CMOS camera at the proximal end of the imaging catheter. However, the hexagonai!y assembled optical fibers can create a honeycomb-like pixel ation artifact, as shown in FIG. 2 A. Each white round area 2605 is a fiber core. The dark gaps 261 between cores are the fiber cladding. Due to these gaps, the speckle images may not be continuous. These gaps can also reduce the Dumber of pixels covered by each speckle, and can therefore reduce the efficiency of spatial average in calculating the temporal statistics of speckles patterns, such as the speckle autocorrelation, within a certain spatial area. Thus, the area can be enlarged to include more pixels to obtain an adequate spatial average. Consequently tlie spatial
resolution o the exemplary maps of the arterial viscoeiasiicity distribution, constructed from, the speckle fluctuations can be degraded. This degradation can limit the ability to distinguish morphological features of tissues, such as the size and shape of plaques. To eliminate the pixelation anifact, two exemplary numerical methods can be applied for two distinct speckle
size regimes (e.g., speckle sizes larger man core spacnigs and speckle sizes smaller than core spacings).
|0Oi74J According to the Nyqttist~Sha.n.non sampling theorem, if the speckle size can be larger than the core spacing between two neighboring cores, the spatial frequencies of the speckle patterns can be lower than that of the hexagonal patten.'! of fibers. Therefore, the hexagonal pattern of fiber cores in Fourier domai can be removed by applying a low pass filter whose cut-off frequency can be no less than the highest spatial frequency of the speckle pattern ,
[ 0I75| The recorded raw images can be transformed (e.g., using a Fourier transform) to spatial frequency domain and then multiplied by a low pass filter HB(u,v) (e.g., a Butterworth tow pass filter), which can rovide, for example:
D(u, v) - [(« - w0 ):' + (v ~v0>? ,
where u and v can be the coordinates in the Fourier domain, uo and vt) can be the center of the filter. Da can be the cut-off frequency and n can be a positive integer. A Butterworth filter can be used because it is a low pass filter with minimal ringing artifacts induced by the shape of the cutting edge owing to the Gibbs phenomenon. Then the product of the Fourier transform of the speckle pattern and the Butterworth filter can be Fourier transformed, back to spatial domain to reconstruct the speckle patterns.
[00176| FIG. 26B illustrates a raw speckle images obtained by an exemplary ILSi catheter from a coronary phantom.. Areas 2615 and 2620 are the speckle patterns reflected from the two opposite area in the phantom. The honeycomb-like pixelation artifact can be easily seen in the Fig. 26 B. Fig. 26C shows the Fourier transform of the raw image. The hexagonal pattern 2625 of the local maximums in Fig. 26C can be due to the hexagonal assembled optical fiber cores. In Fig.26D the Fourier transform is superposed by a Butterworth filter. The filter cutoff frequency can be equal to the spatial frequenc of the fiber cores. Area 2630 gray area is the rejected high frequency area by the low pass filter. The 6 first order hexagonal arranged dots 2635 are at the cutoff region of the low pass filter. If the filter cuiof
frequency is even smaller, the entire periodic pattern can be filtered out such that the
pixelation artifact cart be removed.
j00i7?J For the speckle patterns with speckle size smaller than core spacing, the spatial frequencies of the speckle pattern can be higher than that of the hexagonal pattern of fibers. Thus, simply applying the low pass filter can also remove the high frequenc components of the original speckle patterns. The reconstructed speckle image can also be heavily blurred due to loss of high frequency information. A notch band-rejected filter can be applied for selectively eliminating hexagonal pattern in the Fourier domain. (See e.g.. Reference 126). A notch reject filter can be formed as the product of .multiple Butterworth band-reject filters whose centers are the centers of hexagonal bright spots in the Fourier doma in. The notch filter l½p can be designed as, for example:
where ¾ and v
fe can be the center of the kth bright spot in the Fourier domain and Π can be the multiplication symbol An example of the notch filter is shown in FIGS . 27A and 27B. As shown in FIG. 2? A, hexagonal -arranged maximums of the Fourier transform of the raw speckle image can be covered by dots 2705. The periodic dots 2705 in FIG. 27A are the rejected areas of the notch filter. A 3D view of the exemplary notch filter is shown in FIG. 278. After the notch filter ca be applied, the pixelation artifact can be removed. However, the reconstructed speckle patterns can contain the components whose spatial frequencies ca be higher than the spatial frequencies of the original speckle patterns. To remove the unnecessary high frequency components, an additional Butterworth low pass filter can be applied to the speckle patterns retrieved by using the notch filters. The cutoff frequency of the low pass filter can be set to be larger than the spatial frequencies of the original speckle patterns. The reconstructed speckle pattern is shown in the FIG. 30A, Area 3005 of FIG.
30A can be the area where the pixel intensity can be zero. Outlined regions 3 10 and 3015 are the speckle patterns that can have enough intensity to calculate their temporal statistics.
Exemjpla^
Exemplary Temporal And Spatial Normalisation Of Speckle Patterns,
[00178} in addition to the Brownian motion of light scattering particles, various other effects, such as the fluctuations of output power of laser source, can also cause the
fluctuations of speckle intensity. In order to precisely measure the rate of speckle intensity, temporal fluctuations due to the motion of light scattering, the intensit of each pixel can be divided by the spatially averaged intensity of the corresponding frame. The averaged intensity for each frame can be calculated by averaging the intensity over all pixels. The averaged intensity ca also be temporally smoothed to remove the random noise. FIG. 2SA shows the variation of the total intensity of the speckles with time. Line 2805 represents the smoothed total intensity. FIG. 28B shows the same total intensity over the imaging time afte the pixel intensity is divided by the smoothed average intensity.
[001791 To construct the 2D maps of the viscoelasticity of vessel walls, the spatial variation of speckle intensity lue to the spatial profile of the illumination light can also affect the precision of the measurement of the speckle .fluctuation rate. This can be because the statistics of speckle fluctuations can be dominated by the pixels with high intensity. Thus, the pixel with strong intensity can have more weight than the pixel with low intensity in calculating the statistics of speckle fluctuations. To remove this effect, the averaged speckle patterns over all frames can be calculated. Then the averaged speckle pattern ca be spatially smoothed to remove the residual granular patterns of speckles, A spatially smoothed speckle pattern average over frame sequence is shown in FIG. 29. The intensity of each pixel can.be divided by the corresponding pixel intensity of the spatially smoothed speckle pattern average over imaging time. Therefore, all the pixels can equally contribute to the calculation of the temporal statistics of speckle fluctuations.
Exemplary Speckle Intensity Autocorrelation
f Ml 80| In order to characterize the rate of speckle temporal fluctuations and further ma the viscoelastic properties of vessel wails, the temporal autocorrelation of the speckle intensities gj(At) c
where I(t) and I(t+At) can be the pixel intensities at times t and ί+Δί, and < >pixels and < >
t can indicate spatial and temporal averaging over all the pixels and over the imaging time respectively. However the direct light reflection from the outer sheath and/or other stray light in the 1LS1 catheter can lead to the constant background which can introduce erroneous speckle intensity correlation and the high plateau level of g2(&t) curve. To resolve this issue, the autocovaiiaii.ee (see e.g., Reference 127) of the speckle patterns g2(Al) can be calculated, where, for e
Ci At ) can determine the correlation between the fluctuations around average of the intensity. Cite) can calculate the correlation between the intensity fluctuations around its ensemble average instead of between the intensity itself in g2(At). Because the intensity can include both the speckle intensity and the intensity of the background, if the background light cannot he neglected, the constant background between the intensity can lead to imprecise g2(At). Since the fluctuations of the intensity can come from the time-varying speckle, the correlation between the intensity fluctuations can more precisely measure the rate of the speckle
temporal fluctuations. At the end, the g2(At) or C(At) can be fitted to an exponential function l Ai)::::a*exp(-tA/r)+c where t can be the time, the fitting parameter t can be the decay rate of the speckle correlation functions, a and c are the other fitting parameters, t can also be termed as time constant. This exemplary process can be repeated to calculate spatial and temporal speckle fluctuations from all facet of the onmi-directional mirror incorporated in the exemplary 1LS1 catheter.
Exemplary lime Constant apping And Visualisation
[00181 J To construct 2D maps of the viscoelasticity of tissues, whole imaging area can be divided, as shown in. FIG. 27A, into .multiple small windows (e.g., 40 by 40 pixel windows). The autocorrelation, or the autoco variance of the speckles within each window can be calculated similar to the above. Each windo can have an approximately 50% area overlapped with its 4 neighbors (e.g., top, bottom, left and right neighbors). The different C(At) curves for different small windows in the regio outlined by area 2720 in FIG. 27A are shown in FIG. 30B. Each curves 3020 is a C(.A(a curves s IG. 27A r Each curve 3025 is the exponential fit to the corresponding blue€(Λ( the ex.. Then, the time constants for all windows can be retrieved from the exponential fit. The spatially discrete time constants can then be bi-lmearly interpolated to construct a smooth map of the time constants,
[Q0182| in order to test the exemplary image processing, an Acrylamide gel phantom in a 3D printed mold with 5 slots can be prepared. Each slot can be filled with different gel with different viscoelasticity. The exemplary mold and the exemplary gel filled in are shown, in FIG. 31 A. The gel A contains 4% Acrylamide and 0.025% of bisacrylamide. Gel B contains 5% Acrylamide and 0.025% of bisacrylamide. Gel C contains 5% Acrylamide and 0.055% of bisacrylamide. Gel A has low viscosity while gel C has high viscosity. 3 different time constant maps at 3 different positions in each slot are shown in the FIG. 31 B, As shown in FIG. 31 B a big difference between the maps of the gel A and€ can be observed, as well as between and between the maps of gel A and B. The differences between the 3 maps be obtained at different positions of the same gel are relati vely small.
{00183{ To test the above exemplary methods, a phantom can be prepared using a small piece of swine aorta. A small amount of fat. emulsion can be injected with low viscosity between layers of the aorta to mimic the lipid pool of the coronary plaques. Then, the piece of aorta, can be wrapped into a small tube (e.g., appro imately 3-4 mm in diameter). The swine aorta with injected tat is shown in FIG. 31 C. The exemplary ILSf catheter can be inserted into the tube of aorta and the time varying speckle patterns reflected from the areas of the tube illuminated by the illumination fibers of the catheter can be recorded. At each longitudinal position along the coronary, four r maps can be constructed. Then the catheter can be pullback. a short ncrement to a new position and the imaging can be performed again, The all the τ maps at different positions along the coronary can be longitudinal stitched
together to form 4 long τ maps. AH the τ maps can he stitched together and wrapped on the surface of a cylinder to create 2D cylindrical maps of the viscoeiasticity of the coronary. (00184| An example of wrapping a 2D r map to form a cylindrical view of the maps is shown in FIG. 32A. It can he wrapped onto the surface of a cylinder to form a cylindrical view of the arterial viscoeiasticity map (e.g., FIG. 328). At each longitudinal position, the circumferential distribution of the τ values can be displayed by cross-sectional ring (e.g., FIG. 32C).
(00185) At each position along the coronary, four τ maps can be constructed. All the τ maps can he stitched together and wrapped on the surface of a cy linder to create 2D cylindrical maps of viscoeiasticity of the coronary. An example of wrapping 2D maps to form a cylindrical view of the maps is shown in F G. 32.
Exemplary Other Measure Of The Speckle Flactaation Rate
(00186) Time-varying speckle fields can arise from the interference of laser light scattered by the moving particles in a complex media such as tissoe contain locations of zero intensity. Since both the in- and out-of-phase components of the field can vanish at the position where the intensity can be null, the phase can be undefined there. The locations with zero intensity and undefined phase can be called phase singularities, also called an optical vortex. In addition to the temporal intensify fluctuations of the speckle patterns, the Brownian motion of light scattering particles in tissue can also cause the phase of the speckle field. Therefore, the locations of the optical vortices can also change with time. Thus, the speckle fluctuation rate and the displacement of the optical vortices between speckle frames can be strongly
correlated. The spatial locations of the phase singularities can be tracked over all frames of the speckle sequence. The averaged mean squared displacement, of the speckle vortices can serve as another measure of the speckle fluctuation rate, can measure the viscoeiasticity of tissues. AS the phase of speckle patterns may not be measured using the current ILSl catheter, an exemplary Hilbert transform can he used to generate the pseudo-field U('x,y) (see e.g., Reference 128), where, for example
where I(x,y) ca be the speckle intensity pattern and M {{( ,γ)} can be the Hilheit transform of l(x,y) . Then the phase of the U(x,y) can be called the pseudo-phase <p(xty), which can be, for example: φ(.ν. V ) ~ tan — - -
I(x,,y)
|001 S7| The temporal-spatial behavior of the optical vortices of the pseudo-phase ca be similar to the behavior of the optical vortex of the real phases. (See e.g., Reference 128), The locations of the phase singularity can be obtained by calculating the phase change in a complete counterclockwise circuit around the phase singularity, if there can be a singularity within the closed circuit, the phase change can be ±2π rad. This phase singularity can be described in terms of a topological charge of.*. I . Phase singularities of opposite signs can be created or annihilated in pairs with the evo'lvement of the speckle field.
[00188] FIGS. 32D and 32E show exemplary intensity pattern and the pseudo-phase of this speckle Intensity patterns, respectively. The locations of phase singularities with positive and negative charge are indicated by element 3205 red ***** and element 3210 in the FIG. 32F. The underground area 32.15 is the pseudo-phase of the speckle pattern.
[00189] Two speckle pattern sequences with 50 and 100 frames can be selected. Their temporal autocorrelation g2(t) of the intensity patterns are shown i FIG. 32G. From FIG. 32G, it can be seen that the g2(t) curve of the speckle sequence with 50 frames can decay much faster than the g2(t) curve of the speckle sequence with 1 0 frames. For both sequences, their pseudo-phase can be generated, and the locations of the vortices of all frames can be determined. These locations are then plotted in FIG. 32H. The positively charged vortices are plotted as stars 3220, and the negatively charge vortices are plotted as circles 3225, The locations of each individual vortex over several frames can trace a path called a vortex trail. One example of a trail of an optical vortex of each speckle sequence is outlined in FIGS. 32H and 321. By comparing FIGS. 32H and 321, the trails of the vortices can be seen, and are quite long and straight in a slowly varying sequence (e.g., FIG. 321). In the rapidly decorrelating speckle sequence (e.g.. FIG. 32H). the vortices trails are shorter and tortuous. The straight and long trail can mean thai the vortices stay at the same position for long time and the displacement of the vortex between two consecutive frames can be small. Therefore, the mean squared displacement of the optical vortices can be inversely related to
the time constant of the autocorrelation of the speckle intensity patterns, and can serve as an additional measure of the viscoelasticity of the tissue.
An advantage of utilizing the temporal-spatial behavior of the optical vortices can be that it may only a need few frames to obtain the adequate statistics of the mean squared
displacement of the phase singularities. Therefore, it can greatly shorten the imaging time, while calculating the decorrelation of the speckle frames can require long imaging time that has to be few times longer than the decorrelation time of the speckles.
Exem lary Results And Disc ssion OF iber B» ndle Sel eett n
(ΏΘ1 0) The coupling can be integrated between ail guided modes in different individual cores of a 7-core structure with specifications of type 1 fiber bundle listed in. Table 2 above. Since eac core in this structure can support 1 guided modes at a wavelength of about 690 iim, the total number of guided fiber modes can be 7x19-133, and therefore, the dimensions of the matrices of coupling coefficient between ail modes can be 133 by 133. The amplitudes of the mode coupling coefficients 1 * are shown in FIG. 2 i A. Here the mode inde v can run through all 133 guided modes, it is noted that the first 1 modes, which can propagate in the central fiber, can have large coupling coefficients with the higher order modes of all surrounding fibers while the coupling coefficients between the modes of each of the 6 surrounding fiber and the modes of its 3 nearest neighbors can he much larger than the coupling coefficients to the modes of further cores. Thus, the coupling coeffici ent of modes of each core can be dominant by the coupling to the modes of its nearest neighbors. FIGS.
21B-21F show the intensity in each core, which can be the summation of squared mode amplitudes over all guided modes in the core, which can oscillate between the central fiber and surrounding fiber with propagation distance z. As shown in FIGS, 21B-21 F, the mode 1, 2, 6, 9, 10 of central fiber can initially be excited at 2^0. The mode amplitudes changing with z up to 1 m which can roughly be a typical length of fiber bundles used in medic al endoscopy can then be calculated.
100191 J As the order of excited mode of central fiber can increase from 1 to 10 as shown in FIGS. 21B-21 F, the coupling distance defined as the oscillation period of intensity along with propagation distance z becomes shorter. The intensity in central core represented by line 2105 can't couple to the surrounding cores whose intensity represented by lines 21 10 when
only the fun amental mode of central core can be excited as shown in FIG. 21 8. While FIG . 21F shows thai there can. be multiple coupling distances within 1 m which can indicate strong core-to-core coupling when only mode 10 of central fiber can be initially excited. For the modes with same order of different identical cores, the difference between the propagation constant Δβ of these modes can be 0, and the coupling strength may only depend on the mode coupling coefficient between these modes with same order in each fiber. Since the fields of a higher order modes can extend more into the cladding, the overlapping of higher order mode field can be stronger, and coupling between higher order modes of identical cores can be stronger. For the modes with different order of adjacent cores., due to the di fference in propagation constant, the cross order mode coupling can be neglected which can be observed in FIG . 21 F. if there can be cross order mode coupling, the intensity oscillation between central and surround cores can be more complex than the siraple one period oscillation shown in FIG. 2 IF. Thus, if the number of guided modes in each fiber can be reduced to less than 10, the coupling between cores can be suppressed.
[00192} In order to better understand the effect of different fiber bundles specifications on the coupling efficiency of ail modes, the total intensity coupled from central core to
surrounding cores along with propagation distance for the fiber bundles with different specifications, including core sizes, core spacings and NA, as shown in FIG. 22, can be investigated, fn the exemplary calculation, the initial value of amplitudes of all guided modes of central fiber can be set to be equal, al(0) :::: a2(0) aN(G) ::: (l/N)i/2, and the initial value of the mode amplitudes of surrounding fibers can be all set to zeros. The coupling between fibers for fiber bundles with 3 different core sizes (e.g., 2 μχα, 3 μιη and 4 pro), 3 different core spacings {e.g., 6 μιη, 7 pm and 8 pm), and 3 different NA, (e.g., 0.22, 0.32 and 0.40) can be calculated. As shown in In FIGS. 22A-22C, coupling strength can be stronger as core sizes increase from 2 pm, 3 μτη to 4 μηι because the fibers can support more higher order modes whose coupling can be strong and the overlap of lower order mode ca also be stronger since they can be closer when core size increases. FIGS. 22D-22F show that the large core-to-cote spacing can lead to reduced coupling due to the large separation between mode fie lds because of the less mode field overlaps when the cores can be closer. FIGS. 22G-22I show that the larger NA indicating larger refractive index contrast between
core and cladding material can lead to stonger confinement of mode fields and ca reduce overlapping of modal fields of neighboring fibers.
(001931 The modulation to the transmitted speckle patterns due to the core coupling is shown in FIG. 23A-23C. FIGS. 23A-23C show the numerically generated speckle patterns, the speckle patterns coupled in the fiber bundle at z=0. and speckle pattern at z - 1 m for fiber bundle wi th 3 prn core size, 6 μιη core spacing and 0.40 NA . Strong modulation of the speckle patterns can be observed from the difference of speckle patterns in FIGS. 23B and 23C. The ensemble a verage over 20 speckle realization of correlation functio for fiber bundles with different core sizes, core spacings and NAs are shown in FIGS. 23D-23F, respectively. The large core-to-core separation, small core size and large refractive index contrast between core and cladding materia! can be essential to reliably transmitted speckle patterns. Based on the results shown in FIGS. 22A-22F, and 23A-23I, fiber bundles with 3 μιη core size, S tra core spacing and 0.40 NA can have moderate crosstalk between fibers, and its specifications can be close to those commercially available, such that it can be relatively easy to manufacture. The fiber with 3 tm core size can support 9 modes to avoid strong coupling of higher order diodes, The transmitted speckle patterns at z-O, 1 and 100 cm axe shown in FIGS. 23G-23I, respectively. The modulation of speckle patterns along with z can be less than the modulation of the speckle patterns shown in FIGS. 23 B and 23C. it shows again that the relative large separation can help to suppress core coupling and modulation to speckle patterns. The higher N A ca confine mode field in the core better, but higher A can also increase the number of guided modes and 0,40 A can be the highest currently available contrast of refractive index between core and cladding material of fiber bundles,
(00194) An additional parameter that can influence mode coupling can be the non- uniformity of fibers such as fluctuations of core size and irregular core shape. This non- uniformity can introduce the mismatch in propagation constant R between cores and a small amount of mismatch can extensively reduce mode coupling between fibers. (See e.g..
References 1 15 and 1 17). This great reduction can be observed in FIGS. 24B and 24C, in which the total intensity transferred from central fiber to surrounding fibers for fiber bundles with same core size and 1% and 2% randomness in core size are shown. However even though the fluctuation of core size can introduce mismatch of propagation constants of same
order modes, it can introduce the possibility that different order modes in adjacent fibers have nearly equal and can cause strong coupling between cross order modes of neighboring cores. FIG. 24B shows one example that 5th mode of central fiber and 6th mode of one neighboring fiber is almost the same, such that there can be strong coupling between these two modes. Thus it can be important to utilize the non-uniformity to reduce the cross talk between fibers since if can introduce the possibility of cross order mode coupling. To reduce this possibility, reducing the guided modes of each fiber supported can be performed.
(00195} By combining all above exemplary results, the parameters of fiber bundles can include a core diameter of 3.0 μηι ± 0.3 μηι, or 3.0 urn ± 0.2 μηι, or 3.0 pm ± 0.1 μιη, or a core diameter of 3.0 μιη within measurable error.
[Q0196| The exemplary diameter of the core can have a fluctuation of* 0.02 urn to .·.·£·.· 0.4 μηι; ± 0.02 μχη to ± 0.3 μηι, ± 0.03 urn to ± 0.3 μηα; 0.05 μηι to ± 0.2 μη , or approximately ± 0.1 pm. In some embodiments, the core fluctuation can be approximately 0.06 tm (e.g., 2.0%). An even larger mismatch {e.g., larger than ± 0.4 μτη) could also be used to introduce an even larger mismatch between modes of the cores. However, such a large mismatch in core fluctuation can preferably be used with smaller core diameters (e.g., a core diameter of 2.7 μιΐί, 2.8 μτη, 2.9 μιιι or 3.0 μηι) instead of larger core diameters. Tire exemplary bundle specifications can be used at, and can be based on, a wavelengt of between about 630- 720nrn. The bundle specifications can also be dependent on the illumination wavelength, and can be selected to reduce crosstalk between optical fibers in the exemplary fiber bundle. In some embodiments, the manufacture of the core can provide for such a fluctuation in the core diameter as inherent in the formation process. Thus, it can be an. aspect of the present disclosure that the fluctuation in the core diameter can be defined by the formation of the fiber bundle. In other embodiments, an increased fluctuation as compared to the minimal fluctuation that can be formed can be preferred.
('00197} The fiber bundle can also include a core spacing of 8.0 um ± 0.7 μηι, 8.0 μηι ± 0.5 μαχ, 8.0 pm .·).·. 0.4 μιη, 8.0 μχη ± 0.3 μιιι, 8i) p *. 0.2 μτη, or 8.0 pm-.t- 0.1 μηι, or 8.0 pm within measurable error. The fiber bundle can also include a numerical aperture of at least 0.35, at least 0.36, at least 0.37, at least 0.38, at least 0.39, or at least 0.40. In one
embodiment, the numerical aperture ca be between 0.37 and 0.41 or between 0.38 and 0.41. While the highest NA for current commercially available fiber optic bundles can be
approximately 0.40, higher NA. can be preferred for reducing crosstalk, and the higher MA can be used should they become available, such as a NA of about 0.42, 0,43, 0,44, or 0.45. jjn one embodiment, the fiber bundle has a core diameter of 3.0 gm ± 0.1 um with fluctuations in the core size of*- 0.1 μτη to ± 0.2 p , a core spacing of 8.0 pm ± 0.5 μηι, and a numerical aperture of between 0.38 and 0.41.
[00198} While there ca be variability in each of the parameters as described herein, it can be understood that each of the parameters can be interrelated,, and if it is desirable to change one parameter in the formation of the optical fiber, it can also be advisable to change one or more other parameter to compensate for the initial change.
{00199) The fiber bundle as described herein can reliably transmit speckle patterns at wavelength 690 nrn.
{00200] The fiber optic bundles of the present di sclosure can have reduced inter- fiber crosstalk. In some embodiments, the reduction in inter-fiber crosstalk at a propagation distance of 0.5 m using 690 nm radiation can be at least 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, or more compared to an available fiber optic bundle, such either of the
SCHOTT Type 1 or SCHOTT Type 2 leached image bundles described in Table 2 above, In some embodiments, the inter-fiber crosstalk of the fiber optic bundle can be at a negligible level. While the variance in the fiber size can facilitate that, two adjacent fibers have identical core diameters and thus a relative increase in inter-fiber crosstalk, in some embodiments, the average inter-fiber crosstalk for the fiber optic bundle can be insignificant and provides a near- zero negative contribution to the image quality.
{'00201') The exemplary coupled intensity in cores of one configuration of type I and type 11 fiber bundles with 0Λ pm randomness in core size and exemplary fiber bundles according to the present disclosure are shown in FIGS. 25A-25C, respecti ely. The strong coupling can be seen in both type 1 and type II. fiber bundles while the coupling in fiber bundles of the present disclosure may not be obvious. When fiber bundles are moving, the coupling between cores could change with time because the motion can change the mode overlapping and introduce modulation to the extra phase difference between cores due to bending and twisting of fiber bundles. (See, e.g.. References 124 and 125). However for totally decoupled cores, the effect of fiber bundles motion can be weak, and can be neglected. So, to conduct in vivo LSI, a fiber bundle with fully decoupled cores can be preferred to eliminate
the influence of bundle .motion. The .fiber bundle, as described herein, has shown the small coupling between, cores so that it should not be sensitive to the bundle motion.
(00202| Optical fiber bundles have been demonstrated to be a key component to conduct endoscopic LSI The transmitted laser speckles can be modulated, by inter-fiber coupling reduci ng the accuracy of speckle temporal statistics. As described herein, coupled mode theory can be applied, and the influence of fiber core size, core spacing, numerical aperture and variations in core size on mode coupling and speckle modulation has been analyzed. The analysis of the speckle intensity autocorrelation, of time-resolved speckle frames illustrated that a fiber bundle with about 3±0.1 pm core size, about Slim core spacing and about 0.40 A, can facilitate reliable speckle transmission to conduct endoscopic LSI at about 690 nm. The exemplary results can provide solutions and recommendations for the design, selection and optimization of fiber bundles to conduct endoscopic LSI.
(00203] Figure 40 shows a block diagram of an exemplary embodiment of a system according to the present disclosure. For example, exemplary procedures in accordance with the present disclosure described herein can be performed by a processing arrangement and/or a computing arrangement 4002, Such processing/computing arrangement 4002 can be, for example entirely or a part of, or include, but not limited to, a computer/processor 4004 that can include, for example one or more microprocessors, and use instructions stored on a computer-accessible medium (e.g., RAM, ROM, hard drive, or other storage device).
(00204! As shown in Figure 40. for example a computer-accessible medium 4006 (e.g., as described herein above, a storage device such as a hard disk, floppy disk, memory stick, CD- ROM, RAM, ROM, etc. , or a. collection thereof) can be provided (e.g., in communication with the processing arrangement 4002), The computer-accessible medium 4006 can contain executable instructions 4008 thereon. In addition or alternatively, storage arrangement 4010 can be provided separately from the computer-accessible medium 4006, which can provide the instructions to the processing arrangement 4002 so as to configure the processing arrangement to execute certain exemplary procedures, processes and methods, as described herein above, for exampie.
[00205! Further, the exemplary processing arrangement.4002 can be provided with or include an input/output arrangement 4014. which can include, for example a wired network, a wireless network, the internet, an intranet, a data collection probe, a sensor, etc. As shown in
Figure 40, the exemplary processing arrangement 4002 can be in communication with an exemplary display arrangement 40.12, which, according to certain exemplary embodiments of the present disclosure, can be a touch-screen configured tor inputting information to the processing arrangement in addition to outputting information from the processing
arrangement, for example. Further, the exemplary display 4012 and/or a storage arrangement 4010 can be used to display and/or store data in a user-accessible format and/or user-readable format.
{ 00206] The foregoing merely illustrates the principles of the disclosure. Various modifications and alterations to the described embodiments will be apparent to those skilled in the art in view of the teachings herein. It will thus be appreciated that those skilled in the art will be able to devise numerous systems, arrangements, and procedures which, although not explicitly shown or described herein, embody the principles of the disclosure and can be thus within the spirit and scope of the disclosure. Various different exemplary embodiments can be used together with one another, as wel l as interchangeably therewith, as should be understood by those having ordinary skill in the art. in addition, certain terms used in the present disclosure, including the specification, drawings and claims thereof, can be used synonymously in certain instances, including, but not limited to, for example, data and information, it should be understood that, while these words, and/or other words that can be synonymous to one another, can be used synonymously herein, that there can be instances when such words can be intended to not be used synonymously. Further, to the extent thai the prior art knowledge has not been explicitly incorporated by reference herein above, it is explicitly incorporated herein in its entirety. All publications .referenced are incorporated herein by reference in their entireties.
EXEMPLARY REFERENCES
The following references are hereby incorporated by reference in their entirety.
1. Virmani R, Koiodgie FD, Burke AP, Farb A} Schwartz SM. Lessons from sudden coronary death: A comprehensive morphological classification scheme for atherosclerotic lesions, Arterkateter hramh Vase Biol. 2000;20: 1262-1275
2. Finn AV, Nakano M, Narula J, Koiodgie FD, Virmani R. Concept of
vulnerable/unstable plaque. Afheroscler Throm Vase Biol. 2010;30: 1282-1292
3. Schroeder AP, Faik E. Vulnerable and dangerous coronary plaques.
Atherosclerosis, 1995; 1 18 Suppl:S 141 - 149 4. Koiodgie FD, Burke AP, Farb A, Gold H , Yuan 3. The thin-cap fibroatheronia: A type of vulnerable plaque: The major precursor to acute coronary syndromes. C rr Opi'n Cardiol. 2001 ;16:285-292
5. Virmani R, Burke AP, Koiodgie FD, Farb A. Vulnerable plaque: The pathology of unstable coronary lesions. J btterv Cardial, 2002; 15:439-446 6. Garcia-'Garcia Ft, Mintz GS, Leram A, Vinee D, argolis M, van ED G, Morel
M, ir A, Virmani R, Burke A, Stone G, Serruys PW. Tissue characterization using intravascular radiofrequency data analysis: Recommendations for acquisition, analysis, interpretation and reporting. Euroimervention. 2009;5: 177- 189
7. Murray S, Stables R, Palmer N. Virtual histology imaging in acute coronary syndromes: Useful or j list a research too I , J In vasim Cardial. 2009 ;22 ; 84-91
8. Sangiori G, Clementi F, Cola C, Biondt-Zoccai G. Plaque vulnerability and related coronary event prediction by intravascular ultrasound with virtual histology: its a long way to tipperary. Catheter Cardi c Interv. 2007:70:203-210
9. Jang IK, Tearney G, Bouma B. Visualization of tissue prolapse between coronary stent struts by optical coherence tomography Comparison with intravascular ultrasound.
Circulation. 2001 ;! 4:2754
1.0. Grube E, Gercfcens U, Buellesfel L, Fitzgerald PJ. images in cardiovascular medicine. Intracoronary imaging with optical coherence tomography: A new ' high-resolution technology providing striking visualization in the coronary artery. Circulation.
2002;106:2409-2410 1 1. Jang IK, Bouma BE, Kang DH, Park SJ, Park SW. Seung KB, Choi KB, S shkov M, SchlendorfK, Pomerantsev E, Houser SL, Aretz HT, Tearaey GJ. Visualization of coronary atherosclerotic plaques in patients using optical coherence tomography: Comparison with intra vascular ultrasound, J Am Coll Cardiol, 2002:39:604-609
12. Yahushita H, Bouma BE, Houser SL, Aretz HT, Jang IK, Schlendorf KH,
Kauffman CR, Shishkov M, Kang DH, Halpem EF, Teamey GJ. Characterization of human atherosclerosis by optical coherence tomography. Circulation. 2002; 106: 1640- 1 45
1 . Liebso PR, Klein LW. Intravascular ultrasound in coronary atherosclerosis: A. new approach to clinical assessment. Am Heart . 1992; 123; 1643 - 1660
14. Rogers WJ, Prichard JW, Hu YL, Olson PR, Benckart DH, Kramer CM, Vido DA, Reichek N. Characterization of signal properties in atherosclerotic plaque components by intravascular mri. Jrt rioseier Thromb Vase Biol. 2000;20: 1824-1830
15. Brezinski ME, Teamey GJ, Bouma BE, Izatt J A, Hee MR, Swanson. EA, Southern JF. Fujimoto JG. Optical coherence tomography for optical, biopsy. Properties and demonstration of vascular pathology.. Circulation. 1996;93:1206-1213 16. Teamey G, , Bouma BE. Atherosclerotic plaque characterization by spatial and temporal speckle pattern analysis. Optics Letters. 2002;27:533-535
17, Schmermitnd A, Rodermann J, Erbel R. Intracoronary thermography. Hen.
2003;28:505-512
18. Stefanadis C, Toutouzas K, Tsiamis E, Pitsavos C, Papadimitriou L, Toutouzas P. Identification and stabilization of vulnerable atherosclerotic plaques: The role of coronary thermography and external heat delivery. Indian Heart J. 2001;53:104-109
1.9. Uchicia Y, Fujimori Y, 10 rose J, Oshima. T. Percutaneous coronary angioscopy. Jpn Heart J. 1 92;33:271-294
20. Casscells W, Haihorn B5 David M; Krabach T, Vaughn W McAllister HA, Bearraaa G, Willerson JT. Thermal detection of cellular infiltrates in living atherosclerotic plaques: Possible implications for plaque rupture and thrombosis. Lancet. 1996 347: 1447- 1451
23. Moreno PR, Lodder RA, Pumshotharnan KR} Charash WE, O'Connor WNS Mullet JE, Detection of lipid pool, thin fibrous cap, and inflammatory cells in human aortic atherosclerotic plaques by near-infrared spectroscopy. Circulation. 2002; 105:923-927 22. Waxman. S, Dixon S, P LA, J W M , JL P, D C, I T, R.W , al e. In vi o validation of a catheter-based near- infrared spectroscopy system for detection of lipid core coronary plaques: Initial results and exploratory analysis of the spectroscopic assessment of coronary lipid (spectacl) raultkenter study. . J Am Coil Cardiol Jmg, 2009;2:858-868.
23. Ishibashi F, Yokoyarna S, iyahara , Dabreo A, Weiss ER, lafrati M, Takano M5 Qkamatsu K, Miztsno K, Waxman S. Quanti tati ve colorimetry of atherosclerotic plaque using the l*a*b* color space during angioscopy for the detection of lipid cores underneath thin fibrous caps, hit J Card/ovasc Imaging, 2007;23:679-691
24. Madjid M, Willerson JT, Casscells SW. Intracorohary thermography for detection of high-risk vulnerable plaques. J Am Call Cardial 2006;47;€80~85 25. Arbab-Zadeh A, Hoe J. Quantification of coronary arterial stenoses by
miiltidetecior ct angiography in comparison with coiiventiooai angiography methods, caveats, and implications. JACC Cardiovasc imaging. 201 1 ;4; 1 1-202
26, Chang H, Chung N, Clinical perspective of coronary computed tomographic angiography in diagnosis of coronary artery disease. re J. 201 1;75:246-252 27. Afbustini E, Grasso M, Diegoli M, Pucci A„ Bramerio M. Coronary atherosclerotic plaques with and without thrombus in ischemic heart syndromes: A morphologic,
immunohistochemical and biochemical study. Am J Cardiol. 199 ! ;68.36B-50B
28. Cheruvu PK, Finn AV, Gardner C, Caplan J, Goldstein I. Frequency and distribution of thin-cap fibroatherotna and ruptured plaques in human coronary arteries. J Am
Coll Cardiol 2007; 50; 40-949
29. Serruys PW, Garcia-Garcia HM, Regar E. From postmortem characterization to the in vivo detection of thin capped tlhroatheroroas: The missing link towards percutaneous treatment. J Am Coll Cardiol. 2007;50:950-952
30. Mizttkoshi M5 Imanishi T, Tanaka A, ubo 'Γ, Liu Y, Takarada S. Kitabata H, Tanimoto T, I omukai K, Ishibashi , Akasaka t Clinical classification and plaque •morphology determined by opiical coherence tomography in unstable angina pectoris. Am J Cardiol 20 i 0; i 06:323-328
3 . Tanaka A, Imanishi T, Kitabata H, kubo T, Takarada S, Tanimoto T, Kuroi A, Tsujioka H5 a! e. Lipid rich plaque and myocardial perfusion after successful stenting in patients with non-st segment elevation acure coronary syndrome; An optical coherence tomography study. Eur Heart J. 2009; 30 32, Lee R. Atherosclerotic lesion mechanics versus biology, Z Kardiol. 2000;89:80-84
33. Richardson PD. Da vies ML Born GV. Influence of plaque configuration and stress distribution on Assuring of coronary atherosclerotic plaques. Lancet. - 1 89;2:941-944
34. Williamson SD, Lam. Y, Younis HP, Huang H. On the sensitivity of wall stresses in diseased arteries to variable material properties. J Biomech Eng. 2003:125:147- 155 35. Shah PK. Mechanisms of plaque vulnerability and. rupture. J Am Coll Cardial 20O3;41 ;15S-22S
36. Arroyo L, Lee R. Mechanisms of plaque rupture: Mechanical and biologic interactions. . Cardiovasc Res. 1 99;41 :369-375
37. Richardson PD. Biomechanics of plaque rupture: Progress, problems, and new frontiers, Ann Biomed Eng. 2002:30:524-536
38. f met G, Ohay'on J, ioufol G, L floc S, Tracqui P. Ditbreuli O, Tabib A.
Morphological and biomechanical aspects of vulnerable coronary plaque. Arch Mai Coeur Votes. 2007; 100:547-553
39. Ohayon J, Dubreuil O, Tracqui P, Le Floc'h S, Rioufol G, Chalabreysse L, Thtvolet. F, Pettigrew RI, Fine* G. Influence of residual stress/strain on the biomechanical siabiiiiy of vulnerable coronary plaques: Potential impact for evaluating the risk of plaque rupture. .Am J Physiol Heart Ore Physio!, 2007 ;293:H 1987-1996
40, Sadat U; Teiig Z, Gillard JH. Biomechanical structural stresses of the
atherosclerotic plaque. Expert Rev Cardiovase Ther. 2010;8:1469-1481 41 , Finei G, Ohayon J, Rioufol G. Biomechanical interaction between cap thickness, lipid core composition and blood pressure in vitinerable coronary plaque: Impact on stability or instability. Coron Artery Dis. 2004; 15: 13-20
42. Nadkarni SK, Bouma. BE, Helg T, Chan R, Halpern E, Chau A, Minsky MS, Motz JT, Houser SL, Tearney GJ. Characterization of atherosclerotic plaques by iaser speckle imaging. Circulation. 2005; 1 12:885-892
43. Nadkarni SK, Bilenca A, Bouma BE, Teamey GJ. Measurement of fibrous cap thickness in atherosclerotic plaques by spatiotemporal analysis of laser speckle images. J iomed Opt. 2006:1 :2 006
44. Nadkarni SK, Boum BE, D Y, Tearney G, J. Laser speckle imaging of
atherosclerotic plaques through optica! fiber bundles. J Biomed Opt. 2008; 13 ;054G 16
45. Nadkami SK, Bouma BE, de Boer J, Tearney GJ. Evaluation of collagen in atherosclerotic plaques: The use of two coherent laser-based imaging methods. Lasers Med Sc 2009:24:439-445
46. Hajjarian Z, Xi JQ. Jaffer FA. Teamey GJ, Nadkarni SK. Intravascular laser speckle imaging for the mechanical evaluation of the arterial wall. J Biomed Opt. 20H. ;[epub ahead of print]
47r Koskinas , Feldman C, Chatzizisis Y, Coskun A, Jonas M5 Mayaard C, Baker A, Papafaklis M, Edelman E, Stone P . Natural history of experimental coronary atherosclerosis and vascular remodeling in relation to eiuiothelial shear stress; A serial, in vivo intravascular ultrasound study. Gradation, 2010; 121 :2092-2101 48. Papafaklis , Koskinas , Chateizisis Y, Stone P, Feldman C. In-vivo assessment of the natural history of coronary atherosclerosis: Yasciilar remodeling and endothelial shear stress determine the complexity of atherosclerotic disease progression. Curr Opin Cardiol
2010;25:627-638
49. Koiodgie F, Burke AN, G, Cheng Q, Xu X, Virmani R, Free cholesterol in atherosclerotic plaques: Where does it come from? Curr Opin LipkioL 2007; i 8:500-507
50. Loree H , Tobias BJ, Gibson LX Kamm D, Small DM, Lee RT. Mechanical properties of model atherosclerotic lesion lipid pools. Arierioscier Thromh, 1994;14:230-234
51. Baoriedel G, Eutter R, Welsch U, Bach R, Sievert H, Luderitz B. Role of smooth muscle cell death in advanced coronary primary lesions: Implications for plaque instability. Cardiomsc Res, 99; 1 : 80-488
52. Newby AC, Zaltsman AB. Fibrous cap formation or destruction— the important importance of vascular smooth muscle cel l proliferation, migration and matrix formation.
Cardi vasc Res. 1999;41 :345-360
53. Rekhter MD, Hicks GW; Brammer DWS Hallak H, Kindt E? Chert I, Rosebury WS, Anderson MK, Kuipers PJ, Ryan MJ. Eypercholesterolemia causes mechanical weakening of rabbit atheroma : Local collagen loss as a prerequisite of plaque rupture. Ore Res.
2000;86: 101 -108
54. Chan AH, Chan C„ Shishkov M, acNeill B. Mechanical analysis of
atherosclerotic plaques based on optical coherence tomography. Annals of Biomedical Engineering, 2004;32: 1494-1503
55. Tang D, Yang C, Kobsyasht S, u DN. Effect of a lipid poo! on stress/strain distributions in stenotic arteries: 3~d flwid-strueiure interactions (fsi) models. J Bi mech Eng.
2004; 126:363-370
56. Stone G, Lansky A, Carlier S, al e. A prospective, natural history study of
multimodalily invasive imaging to characterize vulnerable plaque: First report of the baseline ■findings from the prospect trial. Journal of the American College of Cardiology 2007;49: 1 B
57. Stone GW, Serryus PW, de Bruyne B. The prospect trial: An imaging study in patients with unstable atherosclerotic lesions. wwwxlmica rk s.gav. 2009
58. Motoyama S S , Harigaya H, Anno H, Inooe K, Hara T, N'arose H, Ishii J,
Htshida H, Wong ND, Virmani R, Kondo T, Ozaki Y, Narula J. Computed tomographic angiography characteristics of atherosclerotic plaques subsequently resulting in acute coronary syndrome.. J Am Coll Cardiol 200 ;54 -.49-57
59. Manka R JC, ozerke S, Vitanis V, Ctelier G, Gebker R„ Schnackenbarg B,
Boesiger P, -Fleck E» Paetseh I, Dynamic 3-din ertsiona! stress cardiac magnetic resonance perfusion imaging: Detection of coronary artery disease and volumetry of myocardial hypoenhancement before and after coronary stenting.. J Am Coll Cardiol. 201 1 ;57:437-444
60. Tearney GJ, Waxman S, Shishkov M, Vakoc BJ, Suter MS. Three-dimensional coronary artery microscopy by intracoronary optical frequency domai imaging. JACC Cardiovac .Imaging. 2008;1 :752-761 6.1. Jang IK, Tearney GJ, MacNeiH B, Takano iVf , oselewski F. In vivo
characterization of coronary atherosclerotic plaque by use of optical coherence tomography. Circulation. 2005; 11 1 : 1551 - 1556
62, aeNeill BD, Jang IK, Bouma BE, Iftimia N, Takano , Y abushita H, Shishkov M, Kauffman CR, Houser SL, Aretz HT, DeJoseph D, Haipern EF, Tearney GJ. Focal and multi-focal plaque macrophage distributions in patients with acute and stable presentations of coronary artery disease. J Am Coll Cardiol. 2004;44:972-979
63. Baldewsing A, Schaar JA, de Korte CL, Mastik F, Serriiys PW, van der Steen AF. Intra vascular ultrasound eiastography: A clinician's tool for assessing vulnerability and materia! composition of plaques. Stud Health Tecfwol Inform, 2005: 1 13:75-96
64. Schaar JA, De Korte CL, Mastik F, Strijder C, Fasterkamp G, Boersraa E, Semiys PW, Van Der Steen AF. Characterizing vulnerable plaque features with intravascular eiastography . Circulation . 2003 ; 108 :2636-2641
65. Chan C, Chan AH. Karl WG, Nadkarai SK, Khal.il AS. Oci-based arterial eiastography; Rocust estimation exploiting tissue biomechanics. Optics Express.
2004; 12:4558-4572 66. Baldewsing R, Danilonchkine M, Mastik F, Schaar J, Semiys P, van der Steen .
An inverse method for imaging the local elasticity of atherosclerotic coronary plaques. IEEE Trans InfTe hnol Biomed. 2008;12:277-289
67. Baldewsing , Schaar JA, Mastik F, Van der Steen AF. Local elasticity imaging of vulnerable atherosclerotic coronary plaques. Adv Cardiol 2007;44:35- 1 68. Ong W. Laser speckle imaging for the mechanical evaluation of cartilage.
Biomedical Science. 2010;MSc
69. Hajjarian Z, and Nadkarai SK. Are preparing a manuscript titled 'measurement of substrate viscoelastic properties from laser speckle fluctuations,'. 2010
70. Hajjarian Z, Jaffer FA, Tearney GJ, Nadkarai SK. hitracoronary laser speckle imaging for the mechanical evaluation of the arterial wall, submitted to J Biomed Opt, 2010
71. Goodman JW. Statistical optics. Wiley lnterscien.ee; 2000:347-356.
72. Lewis JC. On the einstetn-stoke diffusion coefficient for browman motion in two dimensions. Ph s Letts, 1973;44:245-247
73. Mason TO, Weitz DA. Optical measurements of frequency-dependent linear viscoelastic sty moduli of complex fluids. Physical Review letters. 1995;74: 1250-1253
74. Dasgupta BR, Weitz DA. Microrheology of cross-linked polyacrylamide networks. P ys Rev £ Stat Nonlin Sofi Matter Phys. 2005;? 1 :021504
75. Eajjarian Z, Nadkarni S .. Are preparing a manuscript titled 'measurement, of substrate mechanical properties from laser speckle fluctuations.'. 201 1 76. de Korte CL, Schaar J A, astik F„ Serrays PW, van der Steen AF. Intravascular elastography: From bench to bedside. J lnterv Cardiol, 2003;16:253-259
77. Lee R.T, Richardson SG, Loree HM, Grodzinsky AJ, Gharib SA, Schoen Fj, Pandtan N. Prediction of mechanical properties of human, atherosclerotic tissue by high- frequency intravascular ultrasound imaging. An in vitro study. Arlerioscler Thromb,
1992;12: 1-5
78. Gonik MM, Mishin AB, Zimnyakov DA. Visualization of blood microcirculation parameters in human tissues by time-integrated dynamic speckles analysis. Ann N YAcad Sc 2002;972:325-330
79. Wang L, Jacques SL, Zheng L. Mcml— monte carlo modeling of light transport in multi-layered tissues. Comput Methods Programs Blamed. 1 95;47: 31-146
80. Sadhwani A, Schomacker , T,, Tearney G, J., 'Nishioka NS. Determination of teflon thickness with laser speckle. 1. Potential for burn depth diagnosis. Applied Optics. 1996;35:5727-5735
81. Boas DA, Nishimura G, Yodh AG Diffusing temporal light correlation for burn diagnosis. SPIE. 1999;2979:468-474
82. Waxman S, habbaz K, Connolly R, Tan J, Dabreo A, Egerhei L, ishibashi F> [Vluikr JE. Tearney GM. Intravascular imaging of atherosclerotic human coronaries in a porcine model: A feasibility study. The International Journal of Cardiovascular imaging. 2008;24:37-44
83. Chan AH, Motz JT, Gardecki J A, Waxraan S, Bouma BE, Tearney GJ. Fingerprint and high-wavenumber raman spectroscopy in a human-swine coronary xenograft in vivo. Journal of Biomedical Optics, 2008; 13:04050 i
84. Tearney GJ, Waxman S, Shishkov M, Vakoc BJ, Suter J, Freilich Ml, Desjardins AE, Oil WY, Bart!eit LA, Rosenberg M, Bouma BE. Three-dimensional, coronary artery microscopy by intraeoronary optical frequency domain imaging. JACC Cardiovasc Imaging.
2008;1 :752-761
85. Tearney G, J., Boppart SA, Bouma BE, Brezmski M. Weissman NJ. Scanning single mode fiber optic catheter-endoscope for optical coherence tomograpphy. Optics Letters. 1996;21 :543-545
86. Pratt W. Digital image processing. John Wiley and Sons, Inc.
87. Desai A, Yang MT. Sniadecki NJ, Legant W , Chen CS. Microfabricated post array detectors (rapads): An approach to isolate mechanical forces. Jtmrnal of Visualized
Experiments. 2007 88. "Ian JL, Tien J, Pirone DM, Gray DS, Bhadrirahu , Chen CS. Ceils lying on a bed of microneedles: An approach to isolate mechanical force. Proc Nail Acad Set USA.
2003; 100: 1484- 1489
89. Irimia D. Mierofluidic technologies for temporal perturbations of cheinotaxis. Amtu Rev Biomed Eng. 2010; 12:259-284 90. Bonnema GT, O. Cardinal K, Williams SK, Barton J . A concentric three element radial scanning optical coherence tomography endoscope. J of iaph tonics. 2009;2:353-356
91. Smith-Vaughan Hs Byun R,. Balpm S, Nadkami MA, Jacques A, Hunter N, Morris PS, Leach AJ. Interventions for prevention of otitis media ay be most effective if implemented in the first weeks of life. IntJ P diatr OiarhmolaryngoL 2008;72:57-61 92. Marenzi G, Assanelli E, Campodomco J, Lauri 0, Marana Ϊ, De Metrio M,
Moitrasio M, Grazi. M, ubin o , Vegiia F, Fabbioeebi F, Bartorelli AL. Contrast volume
daring primary percutaneous coronary intervention and subsequent c-on.trast-indu.ced nephropathy and mortality. Ann Intern Med. 2009; 150:170-177
93. Barlis P, Gonzal© "N, Di Mario C, Prati F, Bue!lesfeld L, R.ieber J, Dalhy M'C, Ferrante G . Cera , Grube E, Serruys PW, Regar E. A multicentre ev aluation of the safety of intrac-oronary optical coherence tomography. £urolntervention. 2009;5:90-95
94. Tearney G J, Yabushita H, Houser SL, Areiz HT, Jang IK, Sehlendorf H, Kauffman CR, Sliishkov M, Halpem EF, Bouma BE. Quantification of macrophage content in atherosclerotic plaques by optical coherence tomography; Circulation, 2003;107:113-1 1
95. Nadkami S , Pierce MC, Park BH, de Boer JF, Whittaker P, Bouma BE, Bressner .I E, Halpem E, Houser SL, Tearney GJ. Measurement of collagen and smooth, muscle cell content in atherosclerotic plaques using polarization-sensitive optical coherence tomography. Journal of ike American College of Cardiology. 2007;49:1474-1481
96. J. C, Dainty, Laser Speckle and. Related Phenomena, 2nd enl. (Springer- Verfag, 1984), Vol. v. 9, p. xv, 342 p. 97. J. W. Goodman, Speckle Phenomena in Optics: Theory and Applications (Roberts & Company Publishers, 2007), p. xvi, 387 p.
98. D. Ά Boas and A. K. Dunn, "Laser speckle contrast imaging in biomedical optics.," J. Biomed. Opt 1 5. 01 1 109 (2010).
99. S. . Nadkami, B. E. Bouma, T. flelg, R. Chan, E, Halpem, A. Chan, M. S.
Minsky, J, T, Mntz, S. L. Houser, and G, J. Tearney, "Characterization of atherosclerotic plaques by laser speckle imaging.," Circulation 1 1.2, 885-92 (2005).
100. S. Nadkami, "Laser speckle imaging of aterosclerotic plaques through optical fiber bundles," J. Biomed. Opt 13, 1 -14 (2008).
102, Z. Hajjarian, X XL F, a Jaffer, G. J, Tearney, and S. K. Nadkami, "Intravascular laser speckle imaging catheter tor the mechanical evaluation of the arterial wall.," J. Biomed, Opt. 16, 026005 (201 1 ).
103. . Hajjarian and S. K, Nadkami, "Evaluating the viscoelastie properties of tissue from laser speckle fluctuations.," Sci. Rep. 2, 316 (201.2).
1.04. . Dunn, Ή. Boiay, M. a Moskowit/,, and D. a Boas, "Dynamic imaging of cerebral blood flow using laser speckle.," J. Cereb. Blood Flow Metab. 21, 195-201 (2001).
105. L. Song and D. Elson, "Endoscopic laser speckle contrast imaging system using a fibre image guide," Proc. SPIE 7907, 7 070E-7 070F-9 (201 1).
306, D. A. Zimnyakov, A. P. Svirklov, L, V, Kuzuetsova, S. A. Baranov, and N, Y. Ignatieva, "Monitoring of tissue thermal modification with a bundle-based full- field speckle analyzer;' Appl. Opt. 45, 4480 (2C 107. H. Zhang, P. Li, N. Feng, J. Qiu, B. Li, W. Luo, and Q. Luo, "Correcting the detrimental effects of nonuniform intensity distribution on fiber-transmitting laser speckle imaging of blood flow.," Opt. Express 20, 508-17 (2012).
108. A. F. Gmitro and D, Aziz, "Confocai microscopy through a fiber-optic imaging bundle," Opt. Lett. 18, 565-567 (1993). 109. T. Xie, D. Mukai, S. Guo, M. Brenner, and Z. Chen, "Fiber-optic-biindle-based optical coherence tomography," Opt. Lett. 30, 1803-1805 (2005).
1 10, V. Dubaj, A, MazzoHni, A, Wood, and M. Harris, "Optic fibre bundle contact imaging probe employing a laser scanning confocai microscope," j . Microsc, 207, 108-1 17
11 1. W. O. be I, J. N. D. Kerr, A. Nimmerjahn, and F. Helmchen, "Miniaturized two- photon microscope based on a flex ible coherent fiberhundle and a gradient- index lens objective," Opt. Lett. 29, 2521-2523 (2004).
.1 12. R . Juskattts, T. Wilson, and T. F. Watson, "Real-time white light reflection confocai microscopy using a fibre-optic bundle," Scanning 19, 15-19 (1 97).
1 13, R. Juskaitis, T. Wilson, and T. Watson, "Real-time white light reflection confocai. microscopy using a fibre-optic bundle," Scanning 19, 15-19 (1 97),
i 14. Jia-miiig Liu, Photonic Devices (Cambridge University Press, 2005).
1 15. K. L. Reichenbach and C, Xu, "Numerical analysis of light propagation in image fibers or coherent fiber handles.," Opt. Express 15, 2151-65 (2007).
1.16. X. Chen, . L. Reichenbach, and C. Xu, "Experimental and theoretical analysis of core-to-core coupling on fiber bundle imaging,," Opt. Express 16, 21598-607 (2008).
117. N. Ortega-Quijano, F. FanjiU-Vetez, and J. L. Arce- Diego, "Opttcai crosstalk influence n fiber imaging endoscopes design," Opt. Common. 283, 633-638 (2010).
1 18. A. SagUmi, H. Ford, and R, Tatam, "Numerical modell ing of imaging fibre bundles and their application in. Optical Coherence Tomography," Proc. SPiE 7753, 775350 (2011 ).
1 1 . J.-H. Han and J. U. Kang, "Effect of multimodal coupling in imaging micro- endoscopic fiber bundle on optical coherence tomography.." Appi. Phys. B. 106, 635-643 (2012).
120. A. W. SNYDER, "Coupled-Mode Theory for Optical Fibers," J. Opt. Soc. Am. 62, 1267-1277 (1972).
121. A. Snyder and J. Love, Optical Waveguide Theory (1983).
122. A. Hardy and W. Stretfer, "Coupled mode theory of parallel waveguides," J. Light. Techno!. LT-3, 1 135 (1985).
123. J. Goodman, Statistical Optics (New York, WUey-Interscience,, 2000). 124. M. Skorobogatiy, K. Saitoh, and M. Koshiba, "Fiilt-vectorial coupled mode theory for the evaluation of macro-bending loss in niultimode fibers, application to the hollow-core photonic bandgap fibers.," Opt, Express 16, 14945-53 (2008).
125. M, Koshiba, K. Sai toh, K. Takeftaga, and S, Matsuo, "Mulii-core fiber design and analysts: coupied-inode theory and coupled-power theory.," Opt. Express 19, BJ02~'i J (201 1).
126. C.-Y. Lee and J .-H. Ban, "integrated spatio-spectrai method for efficiently suppressing honeycomb pattern artifact in imaging fiber bundle microscopy," Opt. Common. 306, 67-73 (2013).
127. D, a Boas and A. . Dunn, "Laser speckle contrast imaging in b oniedicaS optics.," i. Biomed, Opt. 15, 01 1 109 (2010),
128, S. J. irkpatrick, , Khaksari, D. Thomas, and D. D. Duncan, "Optical vorte behavior in dynamic speckle fields.," X Biomed. Opt, 17, 050504 (2012).