<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Phys.</journal-id>
<journal-title>Frontiers in Physics</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Phys.</abbrev-journal-title>
<issn pub-type="epub">2296-424X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">613342</article-id>
<article-id pub-id-type="doi">10.3389/fphy.2020.613342</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Physics</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>
<italic>In Vivo</italic> Pulse Wave Measurement Through a Multimode Fiber Diffuse Speckle Analysis System</article-title>
<alt-title alt-title-type="left-running-head">Teng et al.</alt-title>
<alt-title alt-title-type="right-running-head">New Way Measuring Pulse Wave</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Teng</surname>
<given-names>Zhongshuai</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="http://loop.frontiersin.org/people/1101640/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Gao</surname>
<given-names>Feng</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Xia</surname>
<given-names>Hua</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Chen</surname>
<given-names>Wenliang</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Li</surname>
<given-names>Chenxi</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="http://loop.frontiersin.org/people/1103786/overview"/>
</contrib>
</contrib-group>
<aff id="aff1">
<label>
<sup>1</sup>
</label>School of Precision Instruments and Opto-electronics Engineering, Tianjin University, <addr-line>Tianjin</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<label>
<sup>2</sup>
</label>Tianjin Key Laboratory of Biomedical Detecting Techniques and Instruments, <addr-line>Tianjin</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/776030">Chao Tian</ext-link>, University of Science and Technology of China, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1130775">Peng Li</ext-link>, Zhejiang University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/939603">Cheng Wang</ext-link>, University of Shanghai for Science and Technology, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Chenxi Li, <email>lichenxi@tju.edu.cn</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Medical Physics and Imaging, a section of the journal Frontiers in Physics</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>19</day>
<month>01</month>
<year>2021</year>
</pub-date>
<pub-date pub-type="collection">
<year>2020</year>
</pub-date>
<volume>8</volume>
<elocation-id>613342</elocation-id>
<history>
<date date-type="received">
<day>02</day>
<month>10</month>
<year>2020</year>
</date>
<date date-type="accepted">
<day>09</day>
<month>12</month>
<year>2020</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2021 Teng, Gao, Xia, Chen and Li.</copyright-statement>
<copyright-year>2021</copyright-year>
<copyright-holder>Teng, Gao, Xia, Chen and Li</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Continuous monitoring of <italic>in vivo</italic> pulsatile blood flow and pulse wave velocity (PWV) is important for clinical applications. These parameters are correlated with physiological parameters, such as blood pressure and elasticity of blood vessels. A multimode fiber diffuse speckle contrast analysis (MMF-DSCA) system was developed for fast measurement of <italic>in vivo</italic> pulsatile blood flow and pulse wave velocity. With MMF and CCD sensor, the diffuse speckle could be captured and processed with higher temporal resolution of 3&#xa0;ms. We also induced for the first time an MMF-DSCA for evaluation of PWV, which allows estimation of the blood pressure continuously. To validate its performance, both phantom and <italic>in vivo</italic> experiments were conducted. The results demonstrate that MMF-DSCA could achieve fast pulsatile blood flow measurement with detailed information of the pulse wave profile and velocity. Taking the advantages of being simple and cost-effective, the flexible system can be easily adapted for continuous monitoring of vital biosigns, such as heart rate, pulse wave, and blood pressure.</p>
</abstract>
<kwd-group>
<kwd>pulse wave velocity</kwd>
<kwd>blood flow</kwd>
<kwd>multimode optical fibers</kwd>
<kwd>blood pressure</kwd>
<kwd>diffuse speckle</kwd>
</kwd-group>
<contract-num rid="cn001">81871396 81971657 81871393 81671727</contract-num>
<contract-num rid="cn002">19JCYBJC29100 19JCTPJC42200</contract-num>
<contract-sponsor id="cn001">National Natural Science Foundation of China<named-content content-type="fundref-id">10.13039/501100001809</named-content>
</contract-sponsor>
<contract-sponsor id="cn002">Natural Science Foundation of Tianjin City<named-content content-type="fundref-id">10.13039/501100006606</named-content>
</contract-sponsor>
</article-meta>
</front>
<body>
<sec id="s1">
<title>Introduction</title>
<p>Cardiopulmonary parameters, such as pulse wave velocity (PWV) and heart rate, are crucial for clinical diagnosing and daily healthcare monitoring [<xref ref-type="bibr" rid="B1">1</xref>]. PWV is the speed of the pulse wave generated by the heart and transfers along the arterials. It is considerably high (5&#x2013;15&#xa0;m/s) and carries information of cardiovascular function and vessel viability [<xref ref-type="bibr" rid="B2">2</xref>]. The monitoring and analysis of PWV provide a good vital biomarker to assess the status of cardiovascular system and microcirculation [<xref ref-type="bibr" rid="B3">3</xref>
<xref ref-type="bibr" rid="B4 B5">&#x2013;</xref>
<xref ref-type="bibr" rid="B6">6</xref>]. In clinic applications, PWV in the aorta has been estimated by measuring the delay in the foot of the wave between ascending aorta and femoral artery. But the intelligent instruments and skillful operators are needed to obtain reliable results. On the other hand, assessment of PWV is also important for daily healthcare monitoring, especially for infants and elderly [<xref ref-type="bibr" rid="B7">7</xref>].</p>
<p>Taking the advantages of noninvasive and real-time measurement, optical approach is an attractive way to measure <italic>in vivo</italic> blood perfusion, heart rate, and pulse wave. Most of the optical modalities require a coherent light source and follow the working principle of dynamic light scattering [<xref ref-type="bibr" rid="B8">8</xref>]. Laser speckle contrast imaging (LSCI) is a powerful tool for wide-field blood flow imaging of superficial tissue [<xref ref-type="bibr" rid="B9">9</xref>
<xref ref-type="bibr" rid="B10">&#x2013;</xref>
<xref ref-type="bibr" rid="B11">11</xref>]. But the penetration depth limits its applications in deep tissue blood flow and PWV measurement. Considering the highly scattering properties of biological tissue, diffuse optical methods have been developed for blood flow index (BFI) measurement and obtaining the fast pulsatile blood flow in deep tissue. Among these methods, diffuse correlation spectroscopy (DCS) system [<xref ref-type="bibr" rid="B12">12</xref>
<xref ref-type="bibr" rid="B13">&#x2013;</xref>
<xref ref-type="bibr" rid="B14">14</xref>] adopts high-sensitivity single-photon counting and correlator, which increases the hardware cost significantly with multiple channels. However, the fiber-based diffuse speckle contrast analysis (DSCA<bold>)</bold> system [<xref ref-type="bibr" rid="B15">15</xref>, <xref ref-type="bibr" rid="B16">16</xref>] can be extended into multiple channels without significant additional cost and is still able to obtain deep tissue blood perfusion information at satisfied accuracy. Similar to LSCI, DSCA could perform with spatial processing algorithm, which will provide better temporal resolution.</p>
<p>Taking the advantages of flexibility and robustness, fiber-based diffuse optical methods are very popular in clinical applications of deep tissue blood flow measurement. Many systems use multimode (MM) source fiber [<xref ref-type="bibr" rid="B7">7</xref>, <xref ref-type="bibr" rid="B17">17</xref>] to obtain the contrast information from the diffuse laser speckles and extract blood flow information. Therefore, it is also possible to use MM detection fiber combined with area array camera for diffuse pattern detection [<xref ref-type="bibr" rid="B18">18</xref>]. That makes the measurement rate of deep tissue blood flow the same as the frames per second (fps) of camera. Although sensing of pulsatile blood flow and heart rate has been previously demonstrated, other parameters such as PWV have not been further investigated. The previous works also indicated that PWV is highly related to the blood pressure (BP) and age [<xref ref-type="bibr" rid="B19">19</xref>]. Studying the properties and velocity of pulse wave of macro- and microcirculations may lead to an early diagnosis of many disorders [<xref ref-type="bibr" rid="B20">20</xref>].</p>
<p>In this paper, we presented a method for fast pulsatile blood flow and pulse wave velocity measurement in deep tissue. With the MM fiber delivering speckles pattern into the CCD sensor, the diffuse speckle contrast could be calculated spatially at each frame. The MMF-DSCA system achieves 300&#xa0;Hz simultaneous measurement of pulsatile blood flow, which is further used to determine the pulse shapes and temporal delays propagation through the arterial tree. Both phantom validation and <italic>in vivo</italic> blood flow measurement are demonstrated. Thus, the linear regression model gives a good approximation between the BP and PWV. It is demonstrated that MMF-DSCA is one of the fastest noninvasive methods for deep tissue blood flow and pulse wave measurement. Taking the advantages of being simple and cost-effective, MMF-DSCA system can be easily adapted for clinical applications, such as continuous monitoring of heart rate, pulse wave, and blood pressure.</p>
</sec>
<sec sec-type="materials|methods" id="s2">
<title>Materials and Methods</title>
<sec id="s2-1">
<title>Theoretical Background</title>
<p>From theoretical analysis, DCS and LSCI probe different aspects of the field autocorrelation curve, which shows how fast the optical signal loses its self-similarity [<xref ref-type="bibr" rid="B21">21</xref>]. The decay rate of this curve is a good indicator of flow speed. From DCS measurements, it is necessary to extract the speed by calculating the electric field temporal autocorrelation function <inline-formula id="inf1">
<mml:math>
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>&#x3c4;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> [<xref ref-type="bibr" rid="B22">22</xref>], as follows:<disp-formula id="e1">
<mml:math>
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>&#x3c4;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>3</mml:mn>
<mml:msub>
<mml:msup>
<mml:mi>&#x3bc;</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>s</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3c4;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3c4;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
<disp-formula id="e2">
<mml:math>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3c4;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mn>3</mml:mn>
<mml:msub>
<mml:mi>&#x3bc;</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:msub>
<mml:msup>
<mml:mi>&#x3bc;</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
<mml:msub>
<mml:msup>
<mml:mi>&#x3bc;</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mi>k</mml:mi>
<mml:mn>0</mml:mn>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mrow>
<mml:mo>&#x2329;</mml:mo>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3c4;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x232a;</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
<disp-formula id="e3">
<mml:math>
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:msup>
<mml:mi>&#x3c1;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>&#x3bc;</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>s</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
<disp-formula id="e4">
<mml:math>
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:msup>
<mml:mi>&#x3c1;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>&#x3bc;</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>s</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mn>4</mml:mn>
<mml:mrow>
<mml:mn>3</mml:mn>
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>&#x3bc;</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>s</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>where <inline-formula id="inf2">
<mml:math>
<mml:mi>&#x3c1;</mml:mi>
</mml:math>
</inline-formula> is the distance between the source fiber and the detector fiber, <inline-formula id="inf3">
<mml:math>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bc;</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the absorption coefficient of the tissue, and <inline-formula id="inf4">
<mml:math>
<mml:mrow>
<mml:msub>
<mml:msup>
<mml:mi>&#x3bc;</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mi>s</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the reduced scattering coefficient of the tissue. <inline-formula id="inf5">
<mml:math>
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1.440</mml:mn>
<mml:msup>
<mml:mi>n</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.710</mml:mn>
<mml:msup>
<mml:mi>n</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.668</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.00636</mml:mn>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> represents the effective reflection coefficient of the medium and the <inline-formula id="inf6">
<mml:math>
<mml:mi>n</mml:mi>
</mml:math>
</inline-formula> is refractive index of tissue relative to air, <inline-formula id="inf7">
<mml:math>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mtext>&#x3c0;</mml:mtext>
<mml:mo>/</mml:mo>
<mml:mtext>&#x3bb;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> is the wavenumber of light in the medium, <inline-formula id="inf8">
<mml:math>
<mml:mrow>
<mml:mtext>&#xa0;&#x3bb;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula> is wavelength of the incident light, <inline-formula id="inf9">
<mml:math>
<mml:mi>&#x3b1;</mml:mi>
</mml:math>
</inline-formula> is the fraction of dynamic photon scattering events in the medium, and <inline-formula id="inf10">
<mml:math>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x2329;</mml:mo>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3c4;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x232a;</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the mean square displacement of the moving scatterers in a delay time of <inline-formula id="inf11">
<mml:math>
<mml:mi>&#x3c4;</mml:mi>
</mml:math>
</inline-formula>. Brown model approximation is used in analysis and research; it defines <inline-formula id="inf12">
<mml:math>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x2329;</mml:mo>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3c4;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x232a;</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>6</mml:mn>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mi>&#x3c4;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, where <inline-formula id="inf13">
<mml:math>
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mi>B</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the effective diffuse coefficient [<xref ref-type="bibr" rid="B23">23</xref>].</p>
<p>When a coherent light illuminates the blood perfused tissue, the speckle pattern is decorrelated. The level of blurring is quantified by the speckle contrast value, which could be calculated by the following equation [<xref ref-type="bibr" rid="B24">24</xref>]:<disp-formula id="e5">
<mml:math>
<mml:mrow>
<mml:msup>
<mml:mi>K</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>T</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
<mml:mi>T</mml:mi>
</mml:mfrac>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x222b;</mml:mo>
</mml:mstyle>
<mml:mn>0</mml:mn>
<mml:mi>T</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mi>&#x3c4;</mml:mi>
<mml:mo>/</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msubsup>
<mml:mi>g</mml:mi>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3c4;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>&#x3c4;</mml:mi>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>where <inline-formula id="inf14">
<mml:math>
<mml:mi>T</mml:mi>
</mml:math>
</inline-formula> refers to the exposure time, <inline-formula id="inf15">
<mml:math>
<mml:mi>&#x3b2;</mml:mi>
</mml:math>
</inline-formula> is the coherence factor determined by the ratio of detector pixel size to speckle size, and <inline-formula id="inf16">
<mml:math>
<mml:mrow>
<mml:msub>
<mml:mi>g</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3c4;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>&#x3c4;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>,</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the normalized electric field temporal autocorrelation function. In practice, <inline-formula id="inf17">
<mml:math>
<mml:mrow>
<mml:mi>K</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mtext>&#x3c3;</mml:mtext>
<mml:mo>/</mml:mo>
<mml:mi>I</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, where <inline-formula id="inf18">
<mml:math>
<mml:mtext>&#x3c3;</mml:mtext>
</mml:math>
</inline-formula> and <inline-formula id="inf19">
<mml:math>
<mml:mi>I</mml:mi>
</mml:math>
</inline-formula> are the standard deviation and the mean value of speckle pattern, respectively.</p>
<p>Previous work established the fact that both speckle contrast and intensity autocorrelation carried information about the blood flow. Under certain conditions, the recovery of blood flow using both speckle contrast and field/intensity autocorrelation is equivalent. The contrast has the maximal sensitivity if the exposure time is of the order of the correlation time. And there is a linear correlation [<xref ref-type="bibr" rid="B25">25</xref>] between <inline-formula id="inf20">
<mml:math>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:msup>
<mml:mi>K</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> and BFI.</p>
<p>When the MMF is used, photons from the source fiber experience multiple scattering. The diffuser photons that reach the detector fiber will carry the information of blood flow in deep tissue. Because the speckle pattern output by an MMF is not ideal, it is necessary to correct the raw speckle pattern. In each measurement, we average 3,000 images and normalize them to the maximum intensity as the background intensity I<sub>M</sub> (x,y). The raw speckle pattern I<sub>o</sub> (x,y) will be divided by I<sub>M</sub> (x,y) to obtain the corrected result I (x,y), which is used for spatial speckle contrast calculation [<xref ref-type="bibr" rid="B18">18</xref>].</p>
</sec>
<sec id="s2-2">
<title>Pulse Wave Analysis</title>
<p>With MMF-DSCA system, BFI can be calculated spatially from each frame captured by the CCD camera. That makes the sampling rate much faster than DCS systems. The pulsatile changes of BFI are related to the pulse wave, which is an indicator of cardiovascular status. The frequency spectrum of pulse wave provides additional information about the speed transverse to the beam axis, and this initial result merits a more detailed investigation. In this study, the frequency characteristic was analyzed and extracted by the Fourier transform, as<disp-formula id="e6">
<mml:math>
<mml:mrow>
<mml:mtext>P</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>N</mml:mi>
</mml:mfrac>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:munderover>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>where f represents the frequency of signal, N is the length of signal, and <inline-formula id="inf21">
<mml:math>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the signal amplitude.</p>
<p>The diastolic time is an important parameter that is highly related to the systolic blood pressure (SBP) and diastolic blood pressure (DBP). To determine the diastolic time, the first-order derivative algorithm was applied to the time-domain pulse waveform to obtain the characteristic points within one cycle. Then, linear regression analysis was performed between the diastolic time and BP.</p>
<p>As the pulse wave propagates through the arterial tree from the heart to the periphery, pulse wave velocity could be estimated by recording the pressure wave transition time (PTT) between two selected areas in the arterial tree. In this study, we choose the forearm and fingertip as the measuring locations. Considering time-varying BFI signal obtained at adjacent locations in forearm (<italic>F1</italic>) and fingertip (<italic>F2</italic>), when the pulse wave travels from <italic>F1</italic> to <italic>F2</italic>, the BFI feature keeps constant in the stream, resulting in the time series of <italic>F1</italic>
<sub>
<italic>B</italic>
</sub>
<italic>(t)</italic> and <italic>F2</italic>
<sub>
<italic>B</italic>
</sub>
<italic>(t)</italic> appearing in nearly identical shapes with a time lag of <italic>&#x3c4;</italic>. Then, the cross-correlation function is used to determine the time delay between two time-varying dynamic speckle signals, as follows:<disp-formula id="e7">
<mml:math>
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>k</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:msub>
<mml:mn>1</mml:mn>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:msub>
<mml:mn>1</mml:mn>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:msub>
<mml:mn>2</mml:mn>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:msub>
<mml:mn>2</mml:mn>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msqrt>
<mml:mrow>
<mml:msubsup>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:msub>
<mml:mn>1</mml:mn>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:msub>
<mml:mn>1</mml:mn>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msubsup>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:msub>
<mml:mn>2</mml:mn>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:msub>
<mml:mn>2</mml:mn>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
<disp-formula id="e8">
<mml:math>
<mml:mrow>
<mml:mtext>PTT</mml:mtext>
<mml:mo>&#x3d;</mml:mo>
<mml:munder>
<mml:mrow>
<mml:mi>max</mml:mi>
</mml:mrow>
<mml:mi>k</mml:mi>
</mml:munder>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>k</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>where <inline-formula id="inf22">
<mml:math>
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the correlation coefficient, k is the time difference of two signals, N is the maximum time difference, and the <inline-formula id="inf23">
<mml:math>
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:msub>
<mml:mn>1</mml:mn>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf24">
<mml:math>
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:msub>
<mml:mn>2</mml:mn>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="true">&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> are mean value of amplitudes of two signals, respectively. The average PWV can be evaluated with the following equation [<xref ref-type="bibr" rid="B26">26</xref>]:<disp-formula id="e9">
<mml:math>
<mml:mrow>
<mml:mtext>PWV</mml:mtext>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mtext>L</mml:mtext>
<mml:mo>/</mml:mo>
<mml:mrow>
<mml:mtext>PTT</mml:mtext>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>where <italic>L</italic> is the arterial length between two selected areas.</p>
</sec>
<sec id="s2-3">
<title>Experimental Setup</title>
<p>The setup of MMF-DSCA system is shown in <xref ref-type="fig" rid="F1">Figure 1</xref>. A long coherence length (&#x3e;10&#xa0;m) laser diode (785&#xa0;nm, 15&#xa0;mW, LP785-SAV50, Thorlabs, United States) was used to illuminate the sample. In order to achieve synchronization, a 50:50 1 &#xd7; 2 fiber coupler was used to split the light into two identical MMFs (d &#x3d; 200&#xa0;&#x3bc;m, NA &#x3d; 0.22). For <italic>in vivo</italic> measurement, the first fiber was attached to the subject&#x2019;s forearm, near the brachial artery, while the second fiber was attached to the subject&#x2019;s fingertip. Two multimode fibers were used as detector fibers with the other ends touching onto the CCD directly. The distance of source-detector fiber can be adjusted in the range of 5&#x2013;20&#xa0;mm. To avoid the need to synchronize two cameras, the diffuse speckle pattern of both fibers was collected by a single CCD camera (Basler aca1920&#x2013;155um, Germany). The magnification was adjusted to make sure the speckle patterns from the two fibers are projected without overlapping and meet the requirement of diffuse speckle sampling.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Structure diagram of the MMF-DSCA system.</p>
</caption>
<graphic xlink:href="fphy-08-613342-g001.tif"/>
</fig>
<p>For speckle analysis, both the spatial resolution and exposure time are important for the sensitivity and dynamic range of flow measurement. In the MMF-DSCA system, the multimode fibers with core diameter of 200 um were used to collect the diffuse speckle of <italic>in vivo</italic> tissue, resulting in the spatial resolution of 0.2&#xa0;mm. The size of a single speckle is about 11&#xa0;&#x3bc;m, which is twice the size of a single pixel size of 5.6&#xa0;&#x3bc;m. That satisfies the Nyquist sampling criterion and maximizes the contrast of the imaged speckle pattern. However, each MMF speckle pattern contains more than 1,000 speckles, which provides a sufficient statistical sample to analyze the spatial speckle contrast. The exposure time is also important for the sensitivity and SNR of blood flow measurement. Dunn&#x2019;s results suggested that any exposure time greater than 2&#xa0;ms will provide optimal sensitivity to blood flow changes [<xref ref-type="bibr" rid="B27">27</xref>]. However, the exposure time should also meet the requirement of monitoring the pulsatile blood flow. Considering temporal resolution for pulse wave monitoring, we set the exposure time to be 3&#xa0;ms.</p>
<p>In this paper, the diffuse speckle contrast was calculated spatially from each frame, resulting in that the temporal resolution for blood flow measurement is approximately 3.3&#xa0;ms. The corresponding cross-correlation is calculated with <xref ref-type="disp-formula" rid="e7">Eq. 7</xref>, which has the same temporal resolution of 3.3&#xa0;ms. According to the distance between two locations and the speed of PWV, the time delay of pulse waves is generally over 30&#xa0;ms. It is demonstrated that the MMF-DSCA system has a sufficiently high temporal resolution for PWV measurement.</p>
<p>Most of the arteries are in the subcutaneous layer with depth varying from 2 to 10&#xa0;mm below the surface. According to the theory of diffuse optics, the effective detection depth is around a half of source-detector separation. Since source-detector separation is 15&#xa0;mm with <italic>in vivo</italic> measurement, the diffuse speckle signal is mainly from the depth range around 7.5&#xa0;mm. The penetration depth is also demonstrated with Monte Carlo simulations and phantom experiments. The cloud map of the light intensity with Monte Carlo simulation is shown in <xref ref-type="fig" rid="F1">Figure 1</xref>. The simulation results indicated that the penetration depth of diffuse light collected by the detection fibers is in the range of 5&#x2013;12&#xa0;mm. Based on the distribution of photon number and penetration depth, the average penetration depth is further calculated to be 7&#xa0;mm, which agrees with the theory of diffuse optics. The depth of MMF-DSCA measurement is also validated by the phantom experiments.</p>
</sec>
<sec id="s2-4">
<title>Phantom Experiments</title>
<p>To verify the performance of the system in deep tissue flow measurement, a phantom experiment was designed with a hollow plastic tube that was embedded inside a solid scattering phantom body (<italic>&#x3bc;</italic>
<sub>a</sub> &#x3d; 0.01&#xa0;mm<sup>&#x2212;1</sup>, <italic>&#x3bc;</italic>
<sub>s</sub> &#x3d; 0.4&#xa0;mm<sup>&#x2212;1</sup>). The tube was buried 5&#xa0;mm underneath the phantom surface. Liquid with the similar scattering coefficient as blood was pumped through the tube by a peristaltic pump. Since the liquid was pumped by pinching the rubber tube through the rotating rollers of the peristaltic pump, a higher flow rate could be achieved by increasing the pinching frequency. During the experiment, the pumping rate was set from 0 to 0.56 ml/s with step size of 0.07&#xa0;ml/s.</p>
</sec>
<sec id="s2-5">
<title>
<italic>In Vivo</italic> Measurements</title>
<p>For <italic>in vivo</italic> measurements, the subjects were requested to sit on a chair and keep still during the experiments. The experiments were carried out on 10 healthy volunteers (including seven males and three females) aged 24&#x2013;26. To evaluate the relationship between PWV and BP, the experiments were performed on five individuals during the recovery period after the stair climbing exercise, and the interval of each measurement was 10&#xa0;min (5 times in total). The blood pressure was also measured by a commercial BP monitor.</p>
</sec>
</sec>
<sec sec-type="results|discussion" id="s3">
<title>Results and Discussion</title>
<sec id="s3-1">
<title>Phantom Result</title>
<p>As shown in <xref ref-type="fig" rid="F2">Figure 2</xref>, the phantom results demonstrate the good correlation between the flow rate and <inline-formula id="inf25">
<mml:math>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:msup>
<mml:mi>K</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula>. Since the system provided a sampling rate at 300&#xa0;Hz, the flow waveforms at each pumping period could be resolved clearly, as shown in <xref ref-type="fig" rid="F2">Figure 2A</xref>. A good linear relationship between the pumping rate and the averaged BFI can be observed in <xref ref-type="fig" rid="F2">Figure 2B</xref>. Because the MMF-DSCA system provides high temporal resolution, the periodical changes of the flow inside the phantom could be resolved clearly. Therefore, the roller pinching frequency can be calculated from the peak-to-peak time interval of the flow waveform. <xref ref-type="fig" rid="F2">Figure 2C</xref> demonstrates a very good linear relationship between the flow rate and the measured pinching frequency.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Experimental results at different flow rates. <bold>(A)</bold> Flow rate waveform at different pumping rates (only three velocities), <bold>(B)</bold> Averaged flow measurement at each flow rate. <bold>(C)</bold> Frequencies of flow rate waveforms at each flow rate.</p>
</caption>
<graphic xlink:href="fphy-08-613342-g002.tif"/>
</fig>
</sec>
<sec id="s3-2">
<title>
<italic>In Vivo</italic> Result</title>
<sec id="s3-2-1">
<title>Pulse Wave Measurement</title>
<p>The <italic>in vivo</italic> BFI measured with MMF-DSCA system are shown in <xref ref-type="fig" rid="F3">Figure 3A</xref>. The temporal blood flow profile is in good agreement with previous research [<xref ref-type="bibr" rid="B28">28</xref>]. The pulse waveform in the forearm has more features within each cycle. Compared with the standard pulse waveform, the two obvious peaks are related to the main wave and the repulse wave. However, the pulse waveform from the wrist has sharper peaks and smoother repulse wave.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Blood flow measurement of a healthy subject&#x2019;s <bold>(A)</bold> forearm and wrist; <bold>(B)</bold> the power spectrum of pulse wave.</p>
</caption>
<graphic xlink:href="fphy-08-613342-g003.tif"/>
</fig>
<p>To quantify this difference, the time-domain pulse wave signals were Fourier transformed as shown in <xref ref-type="fig" rid="F3">Figure 3B</xref>. The power spectra of the pulse waves demonstrate that the energy in higher frequency of pulse wave on forearm is greater than wrist, which is consistent with temporal profile. Meanwhile, the energy of both pulse waves is concentrated in the range of 0.5&#x2013;10&#xa0;Hz. The peaks in 1.4, 2.8, and 4.2&#xa0;Hz are related to the heartbeat rate and its higher harmonic frequencies. It can be concluded that the frequency spectrum distribution of pulse waves in different parts of the human body is relatively consistent.</p>
<p>Compared with forearm, the radial artery at the wrist is shallower and far from heart. That makes its ascending branch steeper. As shown in <xref ref-type="fig" rid="F3">Figure 3</xref>, the position and amplitude of the repulse wave descend, and the sharp corners of the waveform are smoother. The time-domain and frequency-domain characteristics of the pulse waves measured with our system also agree with the pulse theory [<xref ref-type="bibr" rid="B29">29</xref>].</p>
</sec>
<sec id="s3-2-2">
<title>Correlation Between Blood Pressure and Pulse Wave</title>
<p>In order to analyze the features of the pulse wave in more detail, the individual waveforms of one heart beat were averaged during a period of five adjacent pulse cycles. After baseline correction and noise reduction, the division of temporal pulse wave characteristics is shown as <xref ref-type="fig" rid="F4">Figure 4A</xref>. With systolic peak, dicrotic notch, and diastolic runoff, the well-known characteristics for arterial pulsation behavior could be observed clearly. The waveform of the directed movement shows three peaks located at 0.93&#xa0;ms, 1.08&#xa0;ms, and 1.15&#xa0;ms and relative heights of 87.25, 83.30, and 84.31, respectively. Among them, the first peak (P1) represents maximum vascular pressure and there is a shoulder compared to a second peak (P2) and is followed by a shallower third peak (P3). The diastolic time could be calculated with first-order differential algorithm and selected to validate the correlation between the temporal characteristic and BP.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>
<bold>(A)</bold> Recognition of characteristic points of pulse wave and division of time periods. <bold>(B)</bold> Correlation of the diastolic time and blood pressure.</p>
</caption>
<graphic xlink:href="fphy-08-613342-g004.tif"/>
</fig>
<p>The linear relationship between the diastolic time and BP (SBP and DBP) is shown in <xref ref-type="fig" rid="F4">Figure 4B</xref>. Because the diastolic period occupies most of the time in a cardiac cycle, the amount of blood transported by the large artery to the periphery may decrease with the reduction of diastolic time. As a result, the vascular part is filled with blood and induces a much higher pressure. These results were measured during the resting period on seven volunteers. Diastolic times (t) and BP measured with different volunteers all passed the <italic>t</italic> test with significant <italic>p</italic> &#x3c; 0.05. The result demonstrates that the diastolic time could be preliminarily used for BP prediction.</p>
<p>To test the reliability and consistency of this approach, the <italic>in vivo</italic> experiments were conducted with seven healthy subjects. Every subject was tested more than 10 times. The statistical parameters, such as the average value and standard deviation, were calculated with each subject&#x2019;s results. As shown in <xref ref-type="fig" rid="F4">Figure 4</xref>, the diastolic time varied in the range of 5% with stable PWV measurement.</p>
</sec>
<sec id="s3-2-3">
<title>Correlation Between Blood Pressure and Pulse Wave Velocity</title>
<p>To obtain the PWV, the pulse waves of forearm and fingertip were measured simultaneously. As shown in <xref ref-type="fig" rid="F5">Figure 5</xref>, in addition to the different pulse shape, the pulse wave of fingertip is slightly delayed to forearm. The time delay between two locations is caused by the pressure wave propagation along the artery tree. At the cardiac frequency, the mean delay time is about 33.34&#xa0;ms. Accordingly, the length of the arterial tree pass between the measured spots was about 34.50&#xa0;cm. Based on <xref ref-type="disp-formula" rid="e9">Eq. 9</xref>, the PWV is about 10.36&#xa0;m/s. Given the PWV determined in the human vascular system, as well as the theory of pulse propagation, our results seem realistic.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>
<bold>(A)</bold> Simultaneous recording of pulse wave from the forearm and the fingertip and <bold>(B)</bold> PTT between them.</p>
</caption>
<graphic xlink:href="fphy-08-613342-g005.tif"/>
</fig>
<p>To determine the relationship between BP and PWV, the experiments were performed on five subjects. During exercise, the cardiac output and blood velocities are both increasing. That may cause additional flow resistance and consequent raise of blood pressure. The experimental results were approximated and validated by linear regression and presented in <xref ref-type="fig" rid="F6">Figure 6</xref>. The coefficient of determination values between the SBP and PWV is in the range of 0.88&#x2013;0.95. Thus, the linear regression model gives a good approximation between BP and PWV. However, SBP is highly related to the blood vessel elasticity, because it is caused by the heart&#x27;s active and instantaneous pumping of blood. A rise in SBP causes temporary stiffening of the vessels, which results in a higher PWV. Meanwhile, DBP, which is related to the diastolic time and vascular recovery period, is not sensitive to changes in blood vessel elasticity. As a result, linear approximation of the relation between DBP and PWV shows lower level of correlation.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Correlation between blood pressure and PWV in five subjects: <bold>(A)</bold> SBP vs. PWV; <bold>(B)</bold> DBP vs. PWV.</p>
</caption>
<graphic xlink:href="fphy-08-613342-g006.tif"/>
</fig>
<p>These results demonstrated that MMF-DSCA system achieved fast pulsatile blood flow measurement and extraction of the PWV and BP. However, this system and method also have some disadvantages. The fiber probe needs to be fixed on the location, and the measurement may be influenced by the motion artifacts. In clinical applications, the artifacts can be eliminated with improvements of fiber probe design and signal process. The fiber probe should be carefully designed and combined with special cuffs to fix the probe more stably on the arm and fingertip. Because the motion artifacts are generally embodied as DC or low-frequency signals, high-pass filtering can be used to eliminate the artifacts. Therefore, modern signal process methods, such as wavelet and empirical transverse decomposition filtering method, should also be used for further processing. Furthermore, with multiwavelength light sources, this system could be extended to measure oxygen saturation.</p>
</sec>
</sec>
</sec>
<sec sec-type="conclusion" id="s4">
<title>Conclusion</title>
<p>We present an MMF-DSCA system for measuring instantaneous blood flow and pulse wave velocity in deep tissue. With a CCD sensor, the diffuse speckle could be captured and processed at each frame. This approach enabled the observation of PWV in deep tissue with a temporal resolution of 3&#xa0;ms. Optical synchronization of the speckle patterns measured in different body parts was achieved by using 2&#xa0;MM fibers and a single CCD camera. The simultaneous measurement of pulsatile blood flow in different parts could be further analyzed to obtain the pulse wave characteristics. In this work, we also induced for the first time MMF-DSCA system for evaluation of PWV, which allows estimation of the systolic blood pressure. Taking advantages of being simple and cost-effective, the flexible system will be further developed for continuous monitoring of vital biosigns, such as heart rate and blood pressure, as part of individual healthcare.</p>
</sec>
<sec id="s5">
<title>Data Availability Statement</title>
<p>The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s6">
<title>Ethics Statement</title>
<p>Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. The patients/participants provided their written informed consent to participate in this study.</p>
</sec>
<sec id="s7">
<title>Author Contributions</title>
<p>CL was responsible for conceptualization, project administration, and funding acquisition. ZT was responsible for methodology, software, validation, formal analysis, investigation, data curation, and writing&#x2014;original draft preparation. WC was responsible for resources and visualization. HX was responsible for writing&#x2014;review and editing. FG was responsible for supervision. All authors have read and agreed to the published version of the manuscript.</p>
</sec>
<sec id="s8" sec-type="COI-statement">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
</body>
<back>
<ack>
<p>The authors thank every volunteer who offered help with the experiments. The authors acknowledge the financial support provided by the National Natural Science Foundation of China (81871396, 81971657, 81871393, and 81671727) and Tianjin Natural Science Foundation (19JCYBJC29100 and 19JCTPJC42200).</p>
</ack>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1.</label>
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Hoilett</surname>
<given-names>SO</given-names>
</name>
<name>
<surname>Twibell</surname>
<given-names>AM</given-names>
</name>
<name>
<surname>Rohit</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Linnes</surname>
<given-names>CJ</given-names>
</name>
</person-group>. <article-title>Kick ll: a smartwatch for monitoring respiration and heart rate using photoplethysmography</article-title>. In: <conf-name>2018 40th Annual international conference of the IEEE engineering in medicine and biology society</conf-name>; <conf-date>2018 July 18&#x2013;21</conf-date>; <conf-loc>Honolulu, HI</conf-loc>. <publisher-name>Hawaii: IEEE</publisher-name> (<year>2018). p</year>. <fpage>3821</fpage>&#x2013;<lpage>3824</lpage>. <pub-id pub-id-type="doi">10.1109/EMBC.2018.8513356</pub-id>
</citation>
</ref>
<ref id="B2">
<label>2.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Townsend</surname>
<given-names>RR</given-names>
</name>
<name>
<surname>Wimmer</surname>
<given-names>NJ</given-names>
</name>
<name>
<surname>Chirinos</surname>
<given-names>JA</given-names>
</name>
<name>
<surname>Parsa</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Weir</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Perumal</surname>
<given-names>K</given-names>
</name>
<etal/>
</person-group> <article-title>Aortic PWV in chronic kidney disease: a CRIC ancillary study</article-title>. <source>Am J Hypertens</source> (<year>2010</year>) <volume>23</volume>(<issue>3</issue>):<fpage>282</fpage>&#x2013;<lpage>9</lpage>. <pub-id pub-id-type="doi">10.1038/ajh.2009.240</pub-id> </citation>
</ref>
<ref id="B3">
<label>3.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>JoffeLash</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Morrell</surname>
<given-names>CH</given-names>
</name>
<name>
<surname>Fegatelli</surname>
<given-names>DA</given-names>
</name>
<name>
<surname>Fiorillo</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Delitala</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Orru&#x27;</surname>
<given-names>M</given-names>
</name>
<etal/>
</person-group> <article-title>Arterial stiffness and multiple organ damage: a longitudinal study in population</article-title>. <source>Aging Clin Exp Res</source> (<year>2020</year>) <volume>32</volume>(<issue>5</issue>):<fpage>781</fpage>&#x2013;<lpage>8</lpage>. <pub-id pub-id-type="doi">10.1007/s40520-019-01260-0</pub-id> </citation>
</ref>
<ref id="B4">
<label>4.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>CuccaMarongiu</surname>
<given-names>XN</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>HQ</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>BY</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>YB</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>ZM</given-names>
</name>
<etal/>
</person-group> <article-title>Pulse wave velocity as a marker of arteriosclerosis and its comorbidities in Chinese patients</article-title>. <source>Hypertens Res</source> (<year>2007</year>) <volume>30</volume>(<issue>3</issue>):<fpage>237</fpage>&#x2013;<lpage>42</lpage>. <pub-id pub-id-type="doi">10.1291/hypres.30.237</pub-id> </citation>
</ref>
<ref id="B5">
<label>5.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>WangGao</surname>
<given-names>AY</given-names>
</name>
<name>
<surname>Mohamed</surname>
<given-names>MS</given-names>
</name>
<name>
<surname>Ibrahim</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Hun</surname>
<given-names>TM</given-names>
</name>
<name>
<surname>Musa</surname>
<given-names>KI</given-names>
</name>
<name>
<surname>Yusof</surname>
<given-names>Z</given-names>
</name>
</person-group>. <article-title>Pulse wave velocity as a marker of severity of coronary artery disease</article-title>. <source>J Clin Hypertens</source> (<year>2010</year>) <volume>11</volume>(<issue>1</issue>):<fpage>17</fpage>&#x2013;<lpage>21</lpage>. <pub-id pub-id-type="doi">10.1111/j.1751-7176.2008.00061.x</pub-id> </citation>
</ref>
<ref id="B6">
<label>6.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname>
<given-names>HL</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>SH</given-names>
</name>
</person-group>. <article-title>Pulse wave velocity in atherosclerosis</article-title>. <source>Front Cardiovasc Med</source> (<year>2019</year>) <volume>6</volume>:<fpage>41</fpage>. <pub-id pub-id-type="doi">10.3389/fcvm.2019.00041</pub-id> </citation>
</ref>
<ref id="B7">
<label>7.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bennett</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Beiderman</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Agdarov</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Beiderman</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Hendel</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Straussman</surname>
<given-names>B</given-names>
</name>
<etal/>
</person-group> <article-title>Monitoring of vital bio-signs by analysis of speckle patterns in a fabric-integrated multimode optical fiber sensor</article-title>. <source>Optic Express</source> (<year>2020</year>) <volume>28</volume>(<issue>14</issue>):<fpage>20830</fpage>&#x2013;<lpage>44</lpage>. <pub-id pub-id-type="doi">10.1364/OE.384423</pub-id> </citation>
</ref>
<ref id="B8">
<label>8.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zalevsky</surname>
<given-names>B</given-names>
</name>
<name>
<surname>Jing</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Loo</surname>
<given-names>PC</given-names>
</name>
<name>
<surname>Kijoon</surname>
<given-names>L</given-names>
</name>
</person-group>. <article-title>Optical methods for blood perfusion measurement--theoretical comparison among four different modalities</article-title>. <source>J Optic Soc Am A</source> (<year>2015</year>) <volume>32</volume>(<issue>5</issue>):<fpage>860</fpage>&#x2013;<lpage>6</lpage>. <pub-id pub-id-type="doi">10.1364/JOSAA.32.000860</pub-id> </citation>
</ref>
<ref id="B9">
<label>9.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Boas</surname>
<given-names>DA</given-names>
</name>
<name>
<surname>Dunn</surname>
<given-names>AK</given-names>
</name>
</person-group>. <article-title>Laser speckle contrast imaging in biomedical optics</article-title>. <source>J Biomed Optic</source> (<year>2010</year>) <volume>15</volume>(<issue>1</issue>):<fpage>011109</fpage>. <pub-id pub-id-type="doi">10.1117/1.3285504</pub-id> </citation>
</ref>
<ref id="B10">
<label>10.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ring</surname>
<given-names>LL</given-names>
</name>
<name>
<surname>Strandby</surname>
<given-names>RB</given-names>
</name>
<name>
<surname>Henriksen</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Ambrus</surname>
<given-names>R</given-names>
</name>
<name>
<surname>S&#xf8;rensen</surname>
<given-names>H</given-names>
</name>
<name>
<surname>G&#xf8;tze</surname>
<given-names>JP</given-names>
</name>
<etal/>
</person-group> <article-title>Laser speckle contrast imaging for quantitative assessment of facial flushing during mesenteric traction syndrome in upper gastrointestinal surgery</article-title>. <source>J Clin Monit Comput</source> (<year>2019</year>) <volume>33</volume>(<issue>5</issue>):<fpage>903</fpage>&#x2013;<lpage>10</lpage>. <pub-id pub-id-type="doi">10.1007/s10877-018-0226-0</pub-id> </citation>
</ref>
<ref id="B11">
<label>11.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>AchiamSvendsen</surname>
<given-names>JH</given-names>
</name>
<name>
<surname>Nerup</surname>
<given-names>N</given-names>
</name>
<name>
<surname>Strandby</surname>
<given-names>RB</given-names>
</name>
<name>
<surname>Svendsen</surname>
<given-names>MBS</given-names>
</name>
<name>
<surname>Ambrus</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Svendsen</surname>
<given-names>LB</given-names>
</name>
<etal/>
</person-group> <article-title>Laser speckle contrast imaging and quantitative fluorescence angiography for perfusion assessment</article-title>. <source>Langenbeck&#x2019;s Arch Surg</source> (<year>2019</year>) <volume>404</volume>(<issue>4</issue>):<fpage>1</fpage>&#x2013;<lpage>11</lpage>. <pub-id pub-id-type="doi">10.1007/s00423-019-01789-8</pub-id> </citation>
</ref>
<ref id="B12">
<label>12.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Buckley</surname>
<given-names>EM</given-names>
</name>
<name>
<surname>Parthasarathy</surname>
<given-names>AB</given-names>
</name>
<name>
<surname>Grant</surname>
<given-names>PE</given-names>
</name>
<name>
<surname>Yodh</surname>
<given-names>AG</given-names>
</name>
<name>
<surname>Franceschini</surname>
<given-names>MA</given-names>
</name>
</person-group>. <article-title>Diffuse correlation spectroscopy for measurement of cerebral blood flow: future prospects</article-title>. <source>Neurophotonics</source> (<year>2014</year>) <volume>1</volume>(<issue>1</issue>):<fpage>011009</fpage>. <pub-id pub-id-type="doi">10.1117/1.NPh.1.1.011009</pub-id> </citation>
</ref>
<ref id="B13">
<label>13.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Verdecchia</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Diop</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Lawrence</surname>
<given-names>KS</given-names>
</name>
</person-group>. <article-title>Investigation of the best model to characterize diffuse correlation spectroscopy measurements acquired directly on the brain</article-title>. <source>Biomed Appl Light Scatter</source> (<year>2015</year>) <volume>Ix</volume>:<fpage>9333</fpage>. <pub-id pub-id-type="doi">10.1117/12.2079499</pub-id> </citation>
</ref>
<ref id="B14">
<label>14.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kyle</surname>
<given-names>V</given-names>
</name>
<name>
<surname>Mamadou</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Albert</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Morrison</surname>
<given-names>LB</given-names>
</name>
<name>
<surname>Ting-Yim</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Keith</surname>
<given-names>SL</given-names>
</name>
</person-group>. <article-title>Assessment of a multi-layered diffuse correlation spectroscopy method for monitoring cerebral blood flow in adults</article-title>. <source>Biomed Optic Express</source> (<year>2016</year>) <volume>7</volume>(<issue>9</issue>):<fpage>3659</fpage>&#x2013;<lpage>74</lpage>. <pub-id pub-id-type="doi">10.1364/BOE.7.003659</pub-id> </citation>
</ref>
<ref id="B15">
<label>15.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bi</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>K</given-names>
</name>
</person-group>. <article-title>Deep tissue flowmetry based on diffuse speckle contrast analysis</article-title>. <source>Opt Lett</source> (<year>2013</year>) <volume>38</volume>(<issue>9</issue>):<fpage>1401</fpage>&#x2013;<lpage>3</lpage>. <pub-id pub-id-type="doi">10.1364/OL.38.001401</pub-id> </citation>
</ref>
<ref id="B16">
<label>16.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Seong</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Morgan</surname>
<given-names>JP</given-names>
</name>
<name>
<surname>Mazdeyasna</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>JG</given-names>
</name>
<name>
<surname>Hastings</surname>
<given-names>JT</given-names>
</name>
<etal/>
</person-group> <article-title>Low-cost compact diffuse speckle contrast flowmeter using small laser diode and bare charge-coupled-device</article-title>. <source>J Biomed Optic</source> (<year>2016</year>) <volume>21</volume>(<issue>8</issue>):<fpage>80501</fpage>. <pub-id pub-id-type="doi">10.1117/1.JBO.21.8.080501</pub-id> </citation>
</ref>
<ref id="B17">
<label>17.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yu</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Beiderman</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Agdarov</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Beiderman</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Zalevsky</surname>
<given-names>Z</given-names>
</name>
</person-group>. <article-title>Fiber sensor for non-contact estimation of vital bio-signs</article-title>. <source>Optic Commun</source> (<year>2017</year>) <volume>391</volume>:<fpage>63</fpage>&#x2013;<lpage>7</lpage>. <pub-id pub-id-type="doi">10.1016/j.optcom.2017.01.013</pub-id> </citation>
</ref>
<ref id="B18">
<label>18.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bi</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Singh</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Ho</surname>
<given-names>CJ</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Attia</surname>
<given-names>ABE</given-names>
</name>
<etal/>
</person-group> <article-title>Fast pulsatile blood flow measurement in deep tissue through a multimode detection fiber</article-title>. <source>J Biomed Optic</source> (<year>2020</year>) <volume>25</volume>(<issue>5</issue>):<fpage>1</fpage>&#x2013;<lpage>10</lpage>. <pub-id pub-id-type="doi">10.1117/1.JBO.25.5.055003</pub-id> </citation>
</ref>
<ref id="B19">
<label>19.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>OlivoLi</surname>
<given-names>EJ</given-names>
</name>
<name>
<surname>Park</surname>
<given-names>CG</given-names>
</name>
<name>
<surname>Park</surname>
<given-names>JS</given-names>
</name>
<name>
<surname>Suh</surname>
<given-names>SY</given-names>
</name>
<name>
<surname>Choi</surname>
<given-names>CU</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>JW</given-names>
</name>
<etal/>
</person-group> <article-title>Relationship between blood pressure parameters and pulse wave velocity in normotensive and hypertensive subjects: invasive study</article-title>. <source>J Hum Hypertens</source> (<year>2007</year>) <volume>21</volume>(<issue>2</issue>):<fpage>141</fpage>&#x2013;<lpage>8</lpage>. <pub-id pub-id-type="doi">10.1038/sj.jhh.1002120</pub-id> </citation>
</ref>
<ref id="B20">
<label>20.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>OhKim</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>Z</given-names>
</name>
<name>
<surname>Wenming</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Zibin</surname>
<given-names>Y</given-names>
</name>
</person-group>. <article-title>A research on characteristic information of pulse wave</article-title>. <source>J Beijing Polytech Univ</source> (<year>1996</year>) <volume>22</volume>(<issue>1</issue>):<fpage>71</fpage>&#x2013;<lpage>9</lpage>. </citation>
</ref>
<ref id="B21">
<label>21.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>D</surname>
</name>
<name>
<surname>urduran</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Choe</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Baker</surname>
<given-names>BW</given-names>
</name>
<name>
<surname>Yodh</surname>
<given-names>GA</given-names>
</name>
</person-group>. <article-title>Diffuse optics for tissue monitoring and tomography. Reports on progress in physics</article-title>. <source>Phys Soc</source> (<year>2010</year>) <volume>73</volume>(<issue>7</issue>):<fpage>76701</fpage>&#x2013;<lpage>43</lpage>. <pub-id pub-id-type="doi">10.1088/0034-4885/73/7/076701</pub-id> </citation>
</ref>
<ref id="B22">
<label>22.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xie</surname>
<given-names>JB</given-names>
</name>
<name>
<surname>He</surname>
<given-names>XD</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>LM</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>ZP</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>F</given-names>
</name>
</person-group>. <article-title>Diffuse correlation spectroscopy towards dynamic topography of blood flow index in deep tissues: a multi-channel system and experiment validation</article-title>. <source>Infrared Phys Technol</source> (<year>2020</year>) <volume>107</volume>:<fpage>103298</fpage>. <pub-id pub-id-type="doi">10.1016/j.infrared.2020.103298</pub-id> </citation>
</ref>
<ref id="B23">
<label>23.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Boas</surname>
<given-names>DA</given-names>
</name>
<name>
<surname>Yodh</surname>
<given-names>AG</given-names>
</name>
</person-group>. <article-title>Spatially varying dynamical properties of turbid media probed with diffusing temporal light correlation</article-title>. <source>J Opt Soc Am A</source> (<year>1997</year>) <volume>14</volume>(<issue>1</issue>):<fpage>192</fpage>&#x2013;<lpage>215</lpage>. <pub-id pub-id-type="doi">10.1364/JOSAA.14.000192</pub-id> </citation>
</ref>
<ref id="B24">
<label>24.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bandyopadhyay</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Gittings</surname>
<given-names>AS</given-names>
</name>
<name>
<surname>Suh</surname>
<given-names>SS</given-names>
</name>
<name>
<surname>Dixon</surname>
<given-names>PK</given-names>
</name>
<name>
<surname>Durian</surname>
<given-names>DJ</given-names>
</name>
</person-group>. <article-title>Speckle-visibility spectroscopy: a tool to study time-varying dynamics</article-title>. <source>Rev Sci Instrum</source> (<year>2005</year>) <volume>76</volume>(<issue>9</issue>):<fpage>093110</fpage>. <pub-id pub-id-type="doi">10.1063/1.2037987</pub-id> </citation>
</ref>
<ref id="B25">
<label>25.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>JG</given-names>
</name>
</person-group>. <article-title>Development of simple diffuse optical metabolic spectroscopy for tissue metabolism measurement</article-title>. <source>Biomed Optic Express</source> (<year>2019</year>) <volume>10</volume>(<issue>6</issue>):<fpage>2956</fpage>&#x2013;<lpage>66</lpage>. <pub-id pub-id-type="doi">10.1364/BOE.10.002956</pub-id> </citation>
</ref>
<ref id="B26">
<label>26.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ding</surname>
<given-names>X</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>BP</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>YT</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>N</given-names>
</name>
<name>
<surname>Tsang</surname>
<given-names>HK</given-names>
</name>
</person-group>. <article-title>Pulse transit time based continuous cuffless blood pressure estimation: a new extension and A comprehensive evaluation</article-title>. <source>Sci Rep</source> (<year>2017</year>) <volume>7</volume>(<issue>1</issue>):<fpage>11554</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-017-11507-3</pub-id> </citation>
</ref>
<ref id="B27">
<label>27.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yuan</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Devor</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Boas</surname>
<given-names>DA</given-names>
</name>
<name>
<surname>Dunn</surname>
<given-names>AK</given-names>
</name>
</person-group>. <article-title>Determination of optimal exposure time for imaging of blood flow changes with laser speckle contrast imaging</article-title>. <source>Appl Optic</source> (<year>2005</year>) <volume>44</volume>(<issue>10</issue>):<fpage>1823</fpage>&#x2013;<lpage>30</lpage>. <pub-id pub-id-type="doi">10.1364/ao.44.001823</pub-id> </citation>
</ref>
<ref id="B28">
<label>28.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Michael</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Rice</surname>
<given-names>TB</given-names>
</name>
<name>
<surname>Bruce</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>White</surname>
<given-names>SM</given-names>
</name>
<name>
<surname>Tromberg</surname>
<given-names>BJ</given-names>
</name>
</person-group>. <article-title>Wearable speckle plethysmography (SPG) for characterizing microvascular flow and resistance</article-title>. <source>Biomed Optic Express</source> (<year>2018</year>) <volume>9</volume>(<issue>8</issue>):<fpage>3937</fpage>&#x2013;<lpage>52</lpage>. <pub-id pub-id-type="doi">10.1364/BOE.9.003937</pub-id> </citation>
</ref>
<ref id="B29">
<label>29.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jin</surname>
<given-names>W</given-names>
</name>
<name>
<surname>Sang</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Xin</surname>
<given-names>C</given-names>
</name>
</person-group>. <article-title>Study on pressure pulsation and flow pulsation in pulse diagnosis</article-title>. <source>Clin J Chin Med</source> (<year>2014</year>) <volume>6</volume>(<issue>06</issue>):<fpage>33</fpage>&#x2013;<lpage>4</lpage> [<comment>in Chinese, with English summary</comment>]. <pub-id pub-id-type="doi">10.3969/j.issn.1674-7860.</pub-id> </citation>
</ref>
</ref-list>
</back>
</article>