<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<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="doi">10.3389/fphy.2019.00192</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>Calibration of Spectral Imaging Devices With Oxygenation-Controlled Phantoms: Introducing a Simple Gel-Based Hemoglobin Model</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Le</surname> <given-names>V. N. Du</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/679066/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Manser</surname> <given-names>Myla</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/783577/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Gurm</surname> <given-names>Sunny</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Wagner</surname> <given-names>Ben</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Hayward</surname> <given-names>Joseph E.</given-names></name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Fang</surname> <given-names>Qiyin</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/303010/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Department of Biomedical Engineering, University of California, Davis</institution>, <addr-line>Davis, CA</addr-line>, <country>United States</country></aff>
<aff id="aff2"><sup>2</sup><institution>Engineering Physics, McMaster University</institution>, <addr-line>Hamilton, ON</addr-line>, <country>Canada</country></aff>
<aff id="aff3"><sup>3</sup><institution>Christie Medical Holdings</institution>, <addr-line>Kitchener, ON</addr-line>, <country>Canada</country></aff>
<aff id="aff4"><sup>4</sup><institution>School of Biomedical Engineering, McMaster University</institution>, <addr-line>Hamilton, ON</addr-line>, <country>Canada</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Zhen Cheng, Stanford University, United States</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Roberta Frass-Kriegl, Medical University of Vienna, Austria; Ashwin B. Parthasarathy, University of South Florida, United States</p></fn>
<corresp id="c001">&#x0002A;Correspondence: V. N. Du Le <email>vdule&#x00040;ucdavis.edu</email></corresp>
<fn fn-type="other" id="fn001"><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>21</day>
<month>11</month>
<year>2019</year>
</pub-date>
<pub-date pub-type="collection">
<year>2019</year>
</pub-date>
<volume>7</volume>
<elocation-id>192</elocation-id>
<history>
<date date-type="received">
<day>31</day>
<month>07</month>
<year>2019</year>
</date>
<date date-type="accepted">
<day>04</day>
<month>11</month>
<year>2019</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2019 Le, Manser, Gurm, Wagner, Hayward and Fang.</copyright-statement>
<copyright-year>2019</copyright-year>
<copyright-holder>Le, Manser, Gurm, Wagner, Hayward and Fang</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>An optical phantom model with stable oxygen saturation (StO2) and optical properties (OPs) is the key for quick and accurate calibration of spectral imaging devices. Toward this goal, we have developed gel-based phantoms using Agarose, polystyrene microsphere, purified hemoglobin, and yeast. Theoretical values of OPs were calculated by applying Beer Lambert law and Mie theory whereas a robust fitting algorithm was used to extract StO2. Stability of OPs and StO2 in the phantoms was validated experimentally over a period of 2 weeks and the effect of pre-heating at 65&#x000B0;C degrees on StO2 was observed. Multi-spectral imaging (750, 780, 850, and 945 nm) of the phantoms were also performed on the FDA-approved commercial system to demonstrate the utility of the phantoms.</p></abstract>
<kwd-group>
<kwd>multi-spectral imaging (MSI)</kwd>
<kwd>hemoglobin</kwd>
<kwd>optical properties</kwd>
<kwd>oxygen saturation</kwd>
<kwd>optical phantom</kwd>
<kwd>scattering</kwd>
<kwd>polystyrene microsphere</kwd>
</kwd-group>
<counts>
<fig-count count="5"/>
<table-count count="1"/>
<equation-count count="1"/>
<ref-count count="22"/>
<page-count count="7"/>
<word-count count="3396"/>
</counts>
</article-meta> 
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<p>Optical phantoms play an important role in calibrating and evaluating the performance optical spectroscopy and spectral imaging systems. Many different models of liquid-based, gelatin-based, and solid optical phantoms have been developed for this purpose [<xref ref-type="bibr" rid="B1">1</xref>]. In our previous works, different materials were selected to fabricate liquid-based phantoms for study of diffuse reflectance spectroscopy and fluorescence spectroscopy [<xref ref-type="bibr" rid="B2">2</xref>, <xref ref-type="bibr" rid="B3">3</xref>]. In general, desirable characteristics of optical phantoms are well-controlled optical properties, and the ability to mimic the spectral shape of tissue scattering and absorption. While hemoglobin (Hb) absorption represents a key source of crosstalk in optical properties extraction and distortion in collected signal [<xref ref-type="bibr" rid="B3">3</xref>, <xref ref-type="bibr" rid="B4">4</xref>], its optical signatures over a wide spectral range are the keys to non-invasive assessment of tissue oxygenation with multi-spectral imaging (MSI) systems [<xref ref-type="bibr" rid="B5">5</xref>, <xref ref-type="bibr" rid="B6">6</xref>]. Although powdered Hb (i.e., H0267, Sigma-Aldrich) has been widely used to fabricate phantoms for calibration of fiber-optics-based techniques [<xref ref-type="bibr" rid="B2">2</xref>&#x02013;<xref ref-type="bibr" rid="B4">4</xref>, <xref ref-type="bibr" rid="B7">7</xref>], it is not an ideal material for varying and accurately controlling of oxygen saturation. This is due to the autoxidation of the present heme group leading to the formation of the ferric derivative (methemoglobin) over time in addition to the already existed small amount of methemoglobin in the powder. Furthermore, powdered Hb is expensive whereas wide field near-infrared imaging technique requires large-volume phantoms to meet semi-infinite condition for photon propagation. Thus, the cost of phantom production significantly increases in such cases.</p>
<p>Meanwhile, although real human blood can be a good alternative to study oxygenation, it requires complicated sample preparation and storage. For example, whole blood needs to be in a pH and temperature controlled medium, and requires oxygenators and blood analyzers to change and test oxygen levels [<xref ref-type="bibr" rid="B8">8</xref>, <xref ref-type="bibr" rid="B9">9</xref>]. Another alternative is to use non-biological dye to mimic the absorption of whole blood in its deoxygenated state [<xref ref-type="bibr" rid="B10">10</xref>]. However, such approach is limited to a few near infrared wavelengths.</p>
<p>The purpose of this study is to develop a simple gelatin-based phantom model for quick calibration of spectral imaging systems at a broad-band spectrum. These phantoms should exhibit the spectral absorption characteristics of Hb and can be used to assess oxygen information. Beer Lambert law and Mie theory were used to guide the design of the phantom and the resulting optical properties and StO2 were validated experimentally. Measurements of the phantoms are also performed on a FDA-approved multi-spectral imaging (MSI) system.</p></sec>
<sec sec-type="methods" id="s2">
<title>Methods</title>
<sec>
<title>Optical Phantoms</title>
<p>To simulate tissue scattering, polystyrene microspheres with diameter of 1 &#x003BC;m (07310-15, Polysciences Inc., PA) were used. To simulate tissue absorption, purified human hemoglobin (Multi-4 CO-Oximeter Control Level 2, Instrumentation Laboratory Co., Bedford, MA) was used. The stock concentration of the Hb solution is 150 mg/ml with an average oxygen saturation (StO2) value of 96.9% [<xref ref-type="bibr" rid="B11">11</xref>, <xref ref-type="bibr" rid="B12">12</xref>]. Polystyrene microspheres were selected as scatterers because of their stability over time, and of their well-controlled size and index of refraction [<xref ref-type="bibr" rid="B1">1</xref>]. To vary StO2 value, different yeast concentrations (Y0375, Sigma-Aldrich, Missouri) were added to the Hb solution [<xref ref-type="bibr" rid="B13">13</xref>]. As shown in <xref ref-type="table" rid="T1">Table 1</xref>, two groups of phantoms were created: control (Ctrl.) group without microspheres for transmission measurement and an experimental (Expt.) group with microspheres for MSI measurement.</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>List of optical phantoms and their compositions.</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th valign="top" align="left"><bold>Phantom label</bold></th>
<th valign="top" align="center"><bold>Microsphere</bold></th>
<th valign="top" align="center"><bold>Agarose</bold></th>
<th valign="top" align="center"><bold>Hemoglobin</bold></th>
<th valign="top" align="center"><bold>Yeast</bold></th>
</tr>
<tr>
<th/>
<th valign="top" align="center"><bold>(w/v)</bold></th>
<th valign="top" align="center"><bold>(w/v)</bold></th>
<th valign="top" align="center"><bold>(mg/ml)</bold></th>
<th valign="top" align="center"><bold>(w/v)</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Ctrl.HbY0</td>
<td valign="top" align="center">0</td>
<td valign="top" align="center">0.75%</td>
<td valign="top" align="center">7</td>
<td valign="top" align="center">0</td>
</tr>
<tr>
<td valign="top" align="left">Ctrl.HbY5</td>
<td valign="top" align="center">0</td>
<td valign="top" align="center">0.75%</td>
<td valign="top" align="center">7</td>
<td valign="top" align="center">5%</td>
</tr>
<tr>
<td valign="top" align="left">Ctrl.HbY10</td>
<td valign="top" align="center">0</td>
<td valign="top" align="center">0.75%</td>
<td valign="top" align="center">7</td>
<td valign="top" align="center">10%</td>
</tr>
<tr>
<td valign="top" align="left">Ctrl.HbY15</td>
<td valign="top" align="center">0</td>
<td valign="top" align="center">0.75%</td>
<td valign="top" align="center">7</td>
<td valign="top" align="center">15%</td>
</tr>
<tr>
<td valign="top" align="left">Expt.HbY5</td>
<td valign="top" align="center">0.73%</td>
<td valign="top" align="center">0.75%</td>
<td valign="top" align="center">7</td>
<td valign="top" align="center">5%</td>
</tr>
<tr>
<td valign="top" align="left">Expt.HbY10</td>
<td valign="top" align="center">0.73%</td>
<td valign="top" align="center">0.75%</td>
<td valign="top" align="center">7</td>
<td valign="top" align="center">10%</td>
</tr>
<tr>
<td valign="top" align="left">Expt.HbY15</td>
<td valign="top" align="center">0.73%</td>
<td valign="top" align="center">0.75%</td>
<td valign="top" align="center">7</td>
<td valign="top" align="center">15%</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Concentration of Hb was kept constant at 7 mg/ml in all phantoms, and the concentration of microspheres was 0.73% <italic>w/v</italic> in the experimental group. Hemoglobin concentration of 7 mg/ml was chosen to obtain enough reflection with the multi-spectral imaging system. This concentration is similar to previous studies [<xref ref-type="bibr" rid="B4">4</xref>, <xref ref-type="bibr" rid="B13">13</xref>]. Furthermore, our goal is to measure StO2, and total Hb concentration in the phantom doesn&#x00027;t affect StO2 level. The concentrations of yeast in the final volume were 5, 10, and 15% <italic>w</italic>/<italic>v</italic>, and were used to create phantoms with approximate StO2 values of 90, 85, and 50%, respectively. Polystyrene microsphere has higher density than water, causing it to sink to the bottom of a liquid-based phantom after a period of time. This problem can be avoided using gel-based phantoms. The Agarose PCR grade (R9012HR, Cedarlane Laboratories, Burlington, ON) was selected to fabricate gel-based phantom because of its high gel strength and 100% transmittance at gel concentration of 0.75% <italic>w</italic>/<italic>v</italic>.</p>
<p>The use of yeast to reduce StO2 was previously described by Jang et al. [<xref ref-type="bibr" rid="B13">13</xref>]. The preparation of yeast in the current study is also similar so that yeast was completely dissolved in water at 65&#x000B0;C before mixing with hemoglobin solution. The main differences are: Agar gel is used in this study and the curing process is not required, and microsphere is used in this study as scatterers instead of Intralipid. The concentrated 4% <italic>w/v</italic> gel solution was previously prepared by mixing 0.78 g of Agar in distilled water at 65&#x000B0;C, at which Agar is completely dissolved in 2 h. The concentrated Agar solution was then added to the freshly made mixture of hemoglobin, microsphere and yeast which was preheated at the same temperature in 30 s. The final concentration of Agar in the phantom was 0.75% w/v (<xref ref-type="table" rid="T1">Table 1</xref>) which produces solid phantoms at room temperature. To prevent the formation of air bubbles, a 100 &#x003BC;l-pipette was used to slowly mix the solutions.</p>
<p>In all phantoms, the total absorption coefficient (&#x003BC;<sub>a</sub>) is calculated by applying Beer-Lambert law to the measured transmission from a spectrophotometer (Ultraspec 3,000, Pharmacia Biotech Inc., NJ) whereas &#x003BC;<sub>s</sub>&#x00027; is calculated by applying Mie theory to 1-&#x003BC;m-diameter polystyrene microsphere [<xref ref-type="bibr" rid="B2">2</xref>, <xref ref-type="bibr" rid="B3">3</xref>, <xref ref-type="bibr" rid="B7">7</xref>]. The final volume of each phantom was 15 ml, and all phantoms were stored at 4&#x000B0;C degrees. In addition, a fitting approach was applied to fit the total absorption coefficient, &#x003BC;<sub>a</sub> and to obtain the oxygen saturation value [<xref ref-type="bibr" rid="B4">4</xref>, <xref ref-type="bibr" rid="B7">7</xref>], such that:</p>
<disp-formula id="E1"><label>(1)</label><mml:math id="M1"><mml:mtable class="eqnarray" columnalign="right center left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>&#x003BC;</mml:mi></mml:mrow><mml:mrow><mml:mtext>a</mml:mtext></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mtext>ln</mml:mtext><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x000B7;</mml:mo><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>&#x000B7;</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>&#x000B7;</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B5;</mml:mi></mml:mrow><mml:mrow><mml:mtext>HbO</mml:mtext><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x0002B;</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x000B7;</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003B5;</mml:mi></mml:mrow><mml:mrow><mml:mtext>Hb</mml:mtext></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo>&#x003BB;</mml:mo></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>&#x003BC;</mml:mi></mml:mrow><mml:mrow><mml:mtext>a</mml:mtext><mml:mo>,</mml:mo><mml:mtext>yeast</mml:mtext></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>In Equation (1), &#x003BC;<sub>a, yeast</sub> (cm<sup>&#x02212;1</sup>) is the absorption coefficient of yeast in Agar gel solution, &#x003BB; is wavelength, <italic>f</italic> <sub>1</sub> (mol.liter<sup>&#x02212;1</sup>) is the total concentration of hemoglobin, <italic>f</italic> <sub>2</sub> (dimensionless) is oxygen saturation (StO2), and &#x003B5;<sub>HbO2</sub> (cm<sup>&#x02212;1</sup>.M<sup>&#x02212;1</sup>) and &#x003B5;<sub>Hb</sub> (cm<sup>&#x02212;1</sup>.M<sup>&#x02212;1</sup>) are molar extinction coefficients of oxygenated and deoxygenated hemoglobin, respectively. To obtain <italic>f</italic> <sub>2</sub>, the least square fitting routine &#x0201C;fminsearch&#x0201D; in MATLAB&#x000AE; was used to fit the broadband measurement of &#x003BC;<sub>a</sub> and &#x003BC;<sub>a,yeast</sub> into Equation (1). This optimization method is based on the Nelder-Mead simplex algorithm and has been widely used for spectral analysis [<xref ref-type="bibr" rid="B14">14</xref>&#x02013;<xref ref-type="bibr" rid="B16">16</xref>]. The volume of the phantoms satisfied semi-infinite condition for photon migration. This can be estimated by comparing the phantom thickness to the optical penetration depth obtained from diffusion theory.</p></sec>
<sec>
<title>Multi-Spectral Imaging (MSI) System</title>
<p>A commercial multi-spectral imaging system (Presygen&#x02122;/si-1, Christie Digital Systems, Kitchener, ON) was used to estimate StO2 in the phantoms. The system includes four near Infrared (NIR) LEDs with illumination wavelengths at 740, 780, 850, and 945 nm (&#x000B1;10 nm), and a NIR camera mounted on an imaging head (<xref ref-type="fig" rid="F1">Figure 1</xref>). Each LED is turned on sequentially, and an image is captured using custom controller software. A background image is captured where none of the LEDs are turned on, and is subtracted from NIR images to account for ambient light. The MSI system utilizes a previously developed scatter model to calculate StO2 values of the experimental phantoms [<xref ref-type="bibr" rid="B17">17</xref>]. The calculated StO2 values were compared to those of control phantoms measured with spectrophotometer.</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Image of the commercial multispectral imaging system (Presygen&#x02122;/si-1).</p></caption>
<graphic xlink:href="fphy-07-00192-g0001.tif"/>
</fig></sec></sec>
<sec id="s3">
<title>Results and Discussion</title>
<p><xref ref-type="fig" rid="F2">Figure 2</xref> shows &#x003BC;<sub>a</sub> spectra of four control phantoms over period of 2 weeks whereas (<xref ref-type="fig" rid="F3">Figure 3</xref>) compares the normalized &#x003BC;<sub>a</sub> spectra of four phantoms in the wavelength region of 450&#x02013;1,000 nm. In general, the characteristic absorption peaks at 540 and 580 nm were observed in all phantoms (<xref ref-type="fig" rid="F2">Figures 2</xref>, <xref ref-type="fig" rid="F3">3A</xref>), and no significant change in &#x003BC;<sub>a</sub> spectra was observed in the period of 2 weeks in all phantoms (<xref ref-type="fig" rid="F2">Figure 2</xref>). <xref ref-type="fig" rid="F2">Figure 2E</xref> shows the reduced scattering. In addition, a noticeable &#x003BC;<sub>a</sub> peak at the 750 nm region and a shift upward in &#x003BC;<sub>a</sub> spectrum in the 650 nm region were observed in phantom with high yeast concentrations of 15% (<xref ref-type="fig" rid="F3">Figure 3A</xref>). These are well-known features when oxy-Hb changes into deoxy- Hb [<xref ref-type="bibr" rid="B18">18</xref>, <xref ref-type="bibr" rid="B19">19</xref>]. <xref ref-type="fig" rid="F3">Figure 3B</xref> shows StO2 values in four control phantoms over 2 weeks. In each phantom, the percentage difference between StO2 value in day 1 and week 2 was within 6%. On average, the calculated StO2 value was 94.6% &#x000B1; 2.7%, 92.2% &#x000B1; 3.2%, 81.93 &#x000B1; 5.2%, 48.2 &#x000B1; 5.9% for phantom Ctrl.HbY0, CTrl.HbY5, Ctrl.HbY10, and Ctrl.HbY15, respectively. Error bars in <xref ref-type="fig" rid="F3">Figure 3B</xref> represent the standard deviation for measurements of three similar phantoms that were fabricated separately but with the same yeast concentration (and everything else). Such averaging accounts for errors in fabrication procedure. Comparing week 1 and week 2, StO2 increases 5% in Ctrl.HbY10, and decreases 3% in Ctrl. HbY5, and 7% in Ctrl. HbY15.</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>&#x003BC;<sub>a</sub> spectra over time of four control phantoms <bold>(A)</bold> hemoglobin without yeast (Ctr.HbY0), <bold>(B)</bold> with yeast 5% (Ctr.HbY5), <bold>(C)</bold> with yeast 10% (Ctr.HbY10), <bold>(D)</bold> with yeast 15% (Ctr.HbY15). Hemoglobin concentration in all phantoms was kept constant at 7 mg/ml and &#x003BC;<sub>s</sub>&#x00027; spectrum of baseline microsphere 0.73% w/v <bold>(E)</bold>.</p></caption>
<graphic xlink:href="fphy-07-00192-g0002.tif"/>
</fig>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Oxygen saturation affects hemoglobin absorption spectrum: <bold>(A)</bold> Normalized &#x003BC;<sub>a</sub> spectra with respect to &#x003BC;<sub>a</sub> at 450 nm and <bold>(B)</bold> oxygen saturation (StO2) over time in four control phantoms: hemoglobin without yeast (Ctr.HbY0), with yeast 5% (Ctr.HbY5), 10% (Ctr.HbY10), and 15% (Ctr.HbY15). Hemoglobin concentration in all phantoms was kept constant at 7 mg/ml.</p></caption>
<graphic xlink:href="fphy-07-00192-g0003.tif"/>
</fig>
<p>All the gel phantoms were prepared at 65&#x000B0;C in 30 s. Therefore, to determine the effect of heating duration on StO2, a set of Ctrl.HbY0 phantoms were made, then heated at 65&#x000B0;C in different time periods up to 1 min. In general, prolonging heating time will increase absorption due to the loss of water through evaporation (<xref ref-type="fig" rid="F4">Figure 4A</xref>). More specifically, &#x003BC;<sub>a</sub> increases by 2.5 and 25.4% compared to the reference &#x003BC;<sub>a</sub> if heating period is 30 s and 1 min, respectively. However, &#x003BC;<sub>a</sub> spectral shape was not affected by heating (<xref ref-type="fig" rid="F4">Figure 4B</xref>). Because StO2 is determined by of the relative &#x003BC;<sub>a</sub> values on the spectrum [<xref ref-type="bibr" rid="B4">4</xref>, <xref ref-type="bibr" rid="B20">20</xref>], the stability of &#x003BC;<sub>a</sub> spectral shape ensures stability of the StO2 value. In these phantoms, the average StO2 value was 96 &#x000B1; 3%.</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>Examples of heating can increase hemoglobin concentration but do not affect the shape of absorption spectrum: Effect of prolonged heating time up to 1 min on <bold>(A)</bold> &#x003BC;<sub>a</sub> values at four wavelengths, and <bold>(B)</bold> &#x003BC;<sub>a</sub> spectral shape. Data from phantom Ctrl.HbY0 was presented. Heating temperature is kept at 65&#x000B0;C degrees. The &#x0201C;reference&#x0201D; curve represents a liquid phantom that was not pre-heated.</p></caption>
<graphic xlink:href="fphy-07-00192-g0004.tif"/>
</fig>
<p>The MSI system was used to capture the experimental phantoms so that 2D spatial images of StO2 were generated for each phantom over course of 1 week to account for day to day variations and aging. Each phantom was imaged at the date of manufacture, day 3, and day 7. For each phantom, averaging over three regions of interests (ROIs) was performed to account for spatial variation. (<xref ref-type="fig" rid="F5">Figures 5A&#x02013;C</xref>). As shown in <xref ref-type="fig" rid="F5">Figure 5D</xref>, a good agreement between spectrophotometer measurements and MSI measurement is obtained with an averaged percentage difference of 8%.</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>Measurement with the multi-spectral imaging system (MSI): <bold>(A&#x02013;C)</bold> example multi-spectral reflectance images of three selected phantoms with ROIs, <bold>(D)</bold> the corresponding StO2 measured with MSI measurement vs. transmission method. Red dotted line in <bold>(D)</bold> represents 95% confidence interval.</p></caption>
<graphic xlink:href="fphy-07-00192-g0005.tif"/>
</fig>
<p>In this paper, the phantoms are selected so that the reflectance signal is within the detection capability of the commercial multi-spectral imaging system. Currently, this system can extract useful information from highly turbid media with minimal absorption (where diffusion theory is valid and where signal to noise ratio is high enough). A further increase in yeast concentration and/or hemoglobin concentration can increase the total absorption coefficient, reducing the accuracy of the Presygen system, and hence is avoided. Meanwhile, vaporization is significant with prolonged aging, causing an increase in absorber concentration, thus increasing total absorption of the phantoms. Therefore, longer aging is not considered when using the Presygen system. Nevertheless, the selected ranges of optical properties and oxygen saturation in the current phantom model are within literature for tissues and optical phantoms [<xref ref-type="bibr" rid="B7">7</xref>, <xref ref-type="bibr" rid="B21">21</xref>, <xref ref-type="bibr" rid="B22">22</xref>].</p></sec>
<sec sec-type="conclusions" id="s4">
<title>Conclusions</title>
<p>In this study, we have demonstrated a gel-based optical phantom model with controlled oxygenation that can be used for quick calibration of spectral imaging devices without the use of whole blood or powdered hemoglobin, and that yeast can be used to decrease the oxygen saturation of a commercial blood product that is fully oxygenated and is more cost effective than powdered hemoglobin. In addition, although increasing yeast concentration decreases oxygen saturation of the phantom, the phenomenon does not occur linearly. Also, a short period of heating can vary the concentration of hemoglobin but does not have major effect on the phantom oxygen saturation. In other words, heating changes the magnitude of absorption but does not change its spectral shape. Specifically, &#x003BC;<sub>a</sub> increases about 2.5% due to heating duration of 30 s at 65&#x000B0;C whereas the absorption spectral shape of the selected hemoglobin was not affected. In addition, the current gel model is able to produce the consistent results in the period of 2 weeks with a relative uncertainty of &#x0003C;6%. The phantom is low cost and simple to make with excellent repeatability, which is critical for its practical application. Therefore, the current gel phantoms have desirable characteristics for calibrating multi-spectral imaging systems, including cost effectiveness, well-controlled optical properties and oxygen saturation, and the ability to retain stability over an adequate period of time.</p></sec>
<sec sec-type="data-availability-statement" id="s5">
<title>Data Availability Statement</title>
<p>The datasets generated for this study are available on request to the corresponding author.</p></sec>
<sec id="s6">
<title>Author Contributions</title>
<p>VL supervised the production of all optical phantoms, applied theories and developed algorithms to control the optical properties, and oxygen saturation in the phantom. MM assisted with data measurement and phantom fabrication. SG and BW measured the phantoms with the spectral imaging system. JH and QF advised the progress of the study.</p>
<sec>
<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></sec>
</body>
<back>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pogue</surname> <given-names>BW</given-names></name> <name><surname>Patterson</surname> <given-names>MS</given-names></name></person-group>. <article-title>Review of tissue simulating phantoms for optical spectroscopy, imaging and dosimetry</article-title>. <source>J Biomed Opt.</source> (<year>2006</year>) <volume>11</volume>:<fpage>041102</fpage>. <pub-id pub-id-type="doi">10.1117/1.2335429</pub-id><pub-id pub-id-type="pmid">16965130</pub-id></citation></ref>
<ref id="B2">
<label>2.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Du Le</surname> <given-names>VN</given-names></name> <name><surname>Nie</surname> <given-names>Z</given-names></name> <name><surname>Hayward</surname> <given-names>JE</given-names></name> <name><surname>Farrell</surname> <given-names>TJ</given-names></name> <name><surname>Fang</surname> <given-names>Q</given-names></name></person-group>. <article-title>Measurements of extrinsic fluorescence in Intralipid and polystyrene microspheres</article-title>. <source>Biomed Opt Exp</source>. (<year>2014</year>) <volume>5</volume>:<fpage>2726</fpage>&#x02013;<lpage>35</lpage>. <pub-id pub-id-type="doi">10.1364/BOE.5.002726</pub-id><pub-id pub-id-type="pmid">25136497</pub-id></citation></ref>
<ref id="B3">
<label>3.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Du Le</surname> <given-names>VN</given-names></name> <name><surname>Patterson</surname> <given-names>MS</given-names></name> <name><surname>Farrell</surname> <given-names>TJ</given-names></name> <name><surname>Hayward</surname> <given-names>JE</given-names></name> <name><surname>Fang</surname> <given-names>Q</given-names></name></person-group>. <article-title>Experimental recovery of intrinsic fluorescence and fluorophore concentration in the presence of hemoglobin: spectral effect of scattering and absorption on fluorescence</article-title>. <source>J Biomed Opt.</source> (<year>2015</year>) <volume>20</volume>:<fpage>127003</fpage>. <pub-id pub-id-type="doi">10.1117/1.JBO.20.12.127003</pub-id><pub-id pub-id-type="pmid">26720881</pub-id></citation></ref>
<ref id="B4">
<label>4.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>Q</given-names></name> <name><surname>Le</surname> <given-names>D</given-names></name> <name><surname>Ramella-Roman</surname> <given-names>J</given-names></name> <name><surname>Pfefer</surname> <given-names>J</given-names></name></person-group>. <article-title>Broadband ultraviolet-visible optical property measurement in layered turbid media</article-title>. <source>Biomed Opt Exp.</source> (<year>2012</year>) <volume>3</volume>:<fpage>1226</fpage>&#x02013;<lpage>40</lpage>. <pub-id pub-id-type="doi">10.1364/BOE.3.001226</pub-id><pub-id pub-id-type="pmid">22741070</pub-id></citation></ref>
<ref id="B5">
<label>5.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>Q</given-names></name> <name><surname>He</surname> <given-names>X</given-names></name> <name><surname>Wang</surname> <given-names>Y</given-names></name> <name><surname>Liu</surname> <given-names>H</given-names></name> <name><surname>Xu</surname> <given-names>D</given-names></name> <name><surname>Guo</surname> <given-names>F</given-names></name></person-group>. <article-title>Review of spectral imaging technology in biomedical engineering: achievements and challenges</article-title>. <source>J Biomed Opt.</source> (<year>2013</year>) <volume>18</volume>:<fpage>100901</fpage>. <pub-id pub-id-type="doi">10.1117/1.JBO.18.10.100901</pub-id><pub-id pub-id-type="pmid">24114019</pub-id></citation></ref>
<ref id="B6">
<label>6.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xu</surname> <given-names>RX</given-names></name> <name><surname>Allen</surname> <given-names>DW</given-names></name> <name><surname>Huang</surname> <given-names>J</given-names></name> <name><surname>Gnyawali</surname> <given-names>S</given-names></name> <name><surname>Melvin</surname> <given-names>J</given-names></name> <name><surname>Elgharably</surname> <given-names>H</given-names></name> <etal/></person-group>. <article-title>Developing digital tissue phantoms for hyperspectral imaging of ischemic wounds</article-title>. <source>Biomed Opt Exp.</source> (<year>2012</year>) <volume>3</volume>:<fpage>1433</fpage>&#x02013;<lpage>45</lpage>. <pub-id pub-id-type="doi">10.1364/BOE.3.001433</pub-id><pub-id pub-id-type="pmid">22741088</pub-id></citation></ref>
<ref id="B7">
<label>7.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Du Le</surname> <given-names>VN</given-names></name> <name><surname>Provias</surname> <given-names>J</given-names></name> <name><surname>Murty</surname> <given-names>N</given-names></name> <name><surname>Patterson</surname> <given-names>MS</given-names></name> <name><surname>Nie</surname> <given-names>Z</given-names></name> <name><surname>Hayward</surname> <given-names>JE</given-names></name> <etal/></person-group>. <article-title>Dual-modality optical biopsy of glioblastomas multiforme with diffuse reflectance and fluorescence: <italic>ex vivo</italic> retrieval of optical properties</article-title>. <source>J Biomed Opt.</source> (<year>2017</year>) <volume>22</volume>:<fpage>027002</fpage>. <pub-id pub-id-type="doi">10.1117/1.JBO.22.2.027002</pub-id><pub-id pub-id-type="pmid">28157245</pub-id></citation></ref>
<ref id="B8">
<label>8.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Faber</surname> <given-names>DJ</given-names></name> <name><surname>Aalders</surname> <given-names>MC</given-names></name> <name><surname>Mik</surname> <given-names>EG</given-names></name> <name><surname>Hooper</surname> <given-names>BA</given-names></name> <name><surname>van Gemert</surname> <given-names>MJ</given-names></name> <name><surname>van Leeuwen</surname> <given-names>TG</given-names></name></person-group>. <article-title>Oxygen saturation-dependent absorption and scattering of blood</article-title>. <source>Phys Rev Lett.</source> (<year>2004</year>) <volume>93</volume>:<fpage>028102</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevLett.93.028102</pub-id><pub-id pub-id-type="pmid">15323954</pub-id></citation></ref>
<ref id="B9">
<label>9.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>PJ</given-names></name> <name><surname>Li</surname> <given-names>WC</given-names></name> <name><surname>Xi</surname> <given-names>GM</given-names></name> <name><surname>Wang</surname> <given-names>HQ</given-names></name> <name><surname>Zhang</surname> <given-names>ZH</given-names></name> <name><surname>Yao</surname> <given-names>BC</given-names></name> <etal/></person-group>. <article-title>Biomechanical study of hepatic portal vein in humans and pigs and its value in liver transplantation</article-title>. <source>Transplant Proc.</source> (<year>2009</year>) <volume>41</volume>:<fpage>1906</fpage>&#x02013;<lpage>10</lpage>. <pub-id pub-id-type="doi">10.1016/j.transproceed.2008.10.097</pub-id><pub-id pub-id-type="pmid">19545754</pub-id></citation></ref>
<ref id="B10">
<label>10.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Akl</surname> <given-names>TJ</given-names></name> <name><surname>Long</surname> <given-names>R</given-names></name> <name><surname>McShane</surname> <given-names>MJ</given-names></name> <name><surname>Ericson</surname> <given-names>MN</given-names></name> <name><surname>Wilson</surname> <given-names>MA</given-names></name> <name><surname>Cot&#x000E9;</surname> <given-names>GL</given-names></name></person-group>. <article-title>Optimizing probe design for an implantable perfusion and oxygenation sensor</article-title>. <source>Biomed Opt Exp.</source> (<year>2011</year>) <volume>2</volume>:<fpage>2096</fpage>&#x02013;<lpage>109</lpage>. <pub-id pub-id-type="doi">10.1364/BOE.2.002096</pub-id><pub-id pub-id-type="pmid">21833350</pub-id></citation></ref>
<ref id="B11">
<label>11.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Steinke</surname> <given-names>JM</given-names></name> <name><surname>Shepherd</surname> <given-names>AP</given-names></name></person-group>. <article-title>Effects of temperature on optical absorbance spectra of oxy-, carboxy-, and deoxyhemoglobin</article-title>. <source>Clin Chem</source>. (<year>1992</year>) <volume>38</volume>:<fpage>1360</fpage>&#x02013;<lpage>4</lpage>. <pub-id pub-id-type="pmid">1623605</pub-id></citation></ref>
<ref id="B12">
<label>12.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bosschaart</surname> <given-names>N</given-names></name> <name><surname>Edelman</surname> <given-names>GJ</given-names></name> <name><surname>Aalders</surname> <given-names>MC</given-names></name> <name><surname>van Leeuwen</surname> <given-names>TG</given-names></name> <name><surname>Faber</surname> <given-names>DJ</given-names></name></person-group>. <article-title>A literature review and novel theoretical approach on the optical properties of whole blood</article-title>. <source>Lasers Med Sci.</source> (<year>2014</year>) <volume>29</volume>:<fpage>453</fpage>&#x02013;<lpage>79</lpage>. <pub-id pub-id-type="doi">10.1007/s10103-013-1446-7</pub-id><pub-id pub-id-type="pmid">24122065</pub-id></citation></ref>
<ref id="B13">
<label>13.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jang</surname> <given-names>H</given-names></name> <name><surname>Pfefer</surname> <given-names>TJ</given-names></name> <name><surname>Chen</surname> <given-names>Y</given-names></name></person-group>. <article-title>Solid hemoglobin-polymer phantoms for evaluation of biophotonic systems</article-title>. <source>Opt Lett.</source> (<year>2015</year>) <volume>40</volume>:<fpage>4321</fpage>&#x02013;<lpage>4</lpage>. <pub-id pub-id-type="doi">10.1364/OL.40.004321</pub-id><pub-id pub-id-type="pmid">26371926</pub-id></citation></ref>
<ref id="B14">
<label>14.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jacques</surname> <given-names>SL</given-names></name> <name><surname>Samatham</surname> <given-names>R</given-names></name> <name><surname>Choudhury</surname> <given-names>N</given-names></name></person-group>. <article-title>Rapid spectral analysis for spectral imaging</article-title>. <source>Biomed Opt Exp.</source> (<year>2010</year>) <volume>1</volume>:<fpage>157</fpage>&#x02013;<lpage>64</lpage>. <pub-id pub-id-type="doi">10.1364/BOE.1.000157</pub-id><pub-id pub-id-type="pmid">21258455</pub-id></citation></ref>
<ref id="B15">
<label>15.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Twardowski</surname> <given-names>MS</given-names></name> <name><surname>Boss</surname> <given-names>ES</given-names></name> <name><surname>Sullivan</surname> <given-names>JM</given-names></name> <name><surname>Donaghay</surname> <given-names>PL</given-names></name></person-group>. <article-title>Modeling the spectral shape of absorption by chromophoric dissolved organic matter</article-title>. <source>Mar Chem</source>. (<year>2004</year>) <volume>89</volume>:<fpage>69</fpage>&#x02013;<lpage>88</lpage>. <pub-id pub-id-type="doi">10.1016/j.marchem.2004.02.008</pub-id></citation></ref>
<ref id="B16">
<label>16.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhong</surname> <given-names>X</given-names></name> <name><surname>Wen</surname> <given-names>X</given-names></name> <name><surname>Zhu</surname> <given-names>D</given-names></name></person-group>. <article-title>Lookup-table&#x02013;based inverse model for human skin reflectance spectroscopy: two-layered Monte Carlo simulations and experiments</article-title>. <source>Opt Exp.</source> (<year>2014</year>) <volume>22</volume>:<fpage>1852</fpage>&#x02013;<lpage>64</lpage>. <pub-id pub-id-type="doi">10.1364/OE.22.001852</pub-id><pub-id pub-id-type="pmid">24515194</pub-id></citation></ref>
<ref id="B17">
<label>17.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Savelieva</surname> <given-names>TA</given-names></name> <name><surname>Stratonnikov</surname> <given-names>AA</given-names></name> <name><surname>Loschenov</surname> <given-names>V</given-names></name></person-group>. <article-title>Multi-spectral imaging of oxygen saturation</article-title>, in <source>Proceedings of SPIE - The International Society for Optical Engineering</source>, <fpage>7022</fpage>. <pub-id pub-id-type="doi">10.1117/12.803905</pub-id></citation></ref>
<ref id="B18">
<label>18.</label>
<citation citation-type="web"><person-group person-group-type="author"><name><surname>Prahl</surname> <given-names>SA</given-names></name></person-group>. <article-title>Tabulated molar extinction coefficient for hemoglobin in water</article-title>. In: <person-group person-group-type="editor"><name><surname>Gratzer</surname> <given-names>WB</given-names></name> <name><surname>Kollias</surname> <given-names>N</given-names></name></person-group> editors. Available online at <ext-link ext-link-type="uri" xlink:href="http://omlc.org/spectra/hemoglobin/summary.html">http://omlc.org/spectra/hemoglobin/summary.html</ext-link> (accessed March 4, 1998).</citation></ref>
<ref id="B19">
<label>19.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Prahl</surname> <given-names>SA</given-names></name> <name><surname>van Gemert</surname> <given-names>M</given-names></name> <name><surname>Welch</surname> <given-names>AJ</given-names></name></person-group>. <article-title>Determining the optical properties of turbid media by using the adding&#x02013;doubling method</article-title>. <source>Appl Opt.</source> (<year>1993</year>) <volume>32</volume>:<fpage>559</fpage>&#x02013;<lpage>68</lpage>. <pub-id pub-id-type="doi">10.1364/AO.32.000559</pub-id></citation></ref>
<ref id="B20">
<label>20.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Seo</surname> <given-names>I</given-names></name> <name><surname>Bargo</surname> <given-names>PR</given-names></name> <name><surname>Kollias</surname> <given-names>N</given-names></name></person-group>. <article-title>Simultaneous assessment of pulsating and total blood in inflammatory skin lesions using functional diffuse reflectance spectroscopy in the visible range</article-title>. <source>J Biomed Opt.</source> (<year>2010</year>) <volume>15</volume>:<fpage>060507</fpage>. <pub-id pub-id-type="doi">10.1117/1.3524191</pub-id><pub-id pub-id-type="pmid">21198148</pub-id></citation></ref>
<ref id="B21">
<label>21.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bargo</surname> <given-names>PR</given-names></name> <name><surname>Prahl</surname> <given-names>SA</given-names></name> <name><surname>Goodell</surname> <given-names>TT</given-names></name> <name><surname>Sleven</surname> <given-names>RA</given-names></name> <name><surname>Koval</surname> <given-names>G</given-names></name> <name><surname>Blair</surname> <given-names>G</given-names></name> <etal/></person-group>. <article-title><italic>In vivo</italic> determination of optical properties of normal and tumor tissue with white light reflectance and an empirical light transport model during endoscopy</article-title>. <source>J Biomed Opt.</source> (<year>2005</year>) <volume>10</volume>:<fpage>034018</fpage>&#x02013;<lpage>1815</lpage>. <pub-id pub-id-type="doi">10.1117/1.1921907</pub-id><pub-id pub-id-type="pmid">16229662</pub-id></citation></ref>
<ref id="B22">
<label>22.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jacques</surname> <given-names>SL</given-names></name></person-group>. <article-title>Optical properties of biological tissues: a review</article-title>. <source>Phys Med Biol.</source> (<year>2013</year>) <volume>58</volume>:<fpage>R37</fpage>. <pub-id pub-id-type="doi">10.1088/0031-9155/58/11/R37</pub-id><pub-id pub-id-type="pmid">23666068</pub-id></citation></ref>
</ref-list> 
</back>
</article>