Extended Finite Element Method with Simplified Spherical Harmonics Approximation for the Forward Model of Optical Molecular Imaging

<table class="table-group" id="tab4"><tr><td><table class="table"><tr><td class="thead-hr" colspan="4"><hr/></td></tr><tr class="thead"><td align="left">Method</td><td align="center">FEM</td><td align="center">XFEM</td><td align="center">MC method</td></tr><tr><td class="thead-hr" colspan="4"><hr/></td></tr><tr><td align="left">Number of nodes</td><td align="center">24906</td><td align="center">6541</td><td align="center">/</td></tr><tr><td align="left">Number of elements</td><td align="center">132202</td><td align="center">32398</td><td align="center">/</td></tr><tr><td align="left"><span style="width: 143.89999px;"><svg height="20.137501" id="M179" style="vertical-align:-4.8335pt" version="1.1" viewbox="0 0 143.89999 20.137501" width="143.89999" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,14.05)"><path d="M614 175l29 -10q-33 -109 -57 -154q-121 -26 -184 -26q-90 0 -160.5 29t-112.5 77t-63.5 105.5t-21.5 119.5q0 157 108 253t277 96q36 0 71.5 -5t69 -13.5t36.5 -8.5q15 -102 20 -150l-29 -8q-20 79 -66.5 114t-128.5 35q-119 0 -187.5 -86t-68.5 -207
q0 -140 73.5 -227.5t188.5 -87.5q73 0 119.5 37.5t86.5 116.5z" id="x43"></path></g><g transform="matrix(.017,-0,0,-.017,11.469,14.05)"><path d="M381 665q131 0 226.5 -93t95.5 -239q0 -158 -96 -253t-238 -95q-137 0 -231 95t-94 238q0 142 92.5 244.5t244.5 102.5zM359 629q-89 0 -151 -74.5t-62 -208.5q0 -141 69 -233t175 -92q90 0 150.5 74t60.5 211q0 152 -69 237.5t-173 85.5z" id="x4F"></path></g><g transform="matrix(.017,-0,0,-.017,24.031,14.05)"><path d="M631 18l-4 -26q-24 0 -35 1q-57 3 -92.5 23t-69.5 68q-21 30 -101 160q-13 22 -30.5 31t-54.5 9h-32v-159q0 -63 14 -77.5t74 -19.5v-28h-262v28q62 5 76 19.5t14 77.5v402q0 62 -13 76t-73 19v28h248q109 0 159 -34q67 -42 67 -131q0 -115 -127 -172q27 -52 89 -144
q49 -71 81 -104q31 -36 72 -47zM212 316h48q73 0 108 29q55 42 55 124q0 76 -42 112t-104 36q-41 0 -54 -10q-11 -7 -11 -44v-247z" id="x52"></path></g><g transform="matrix(.017,-0,0,-.017,34.587,14.05)"><path d="M631 18l-4 -26q-24 0 -35 1q-57 3 -92.5 23t-69.5 68q-21 30 -101 160q-13 22 -30.5 31t-54.5 9h-32v-159q0 -63 14 -77.5t74 -19.5v-28h-262v28q62 5 76 19.5t14 77.5v402q0 62 -13 76t-73 19v28h248q109 0 159 -34q67 -42 67 -131q0 -115 -127 -172q27 -52 89 -144
q49 -71 81 -104q31 -36 72 -47zM212 316h48q73 0 108 29q55 42 55 124q0 76 -42 112t-104 36q-41 0 -54 -10q-11 -7 -11 -44v-247z" id="x52"></path></g><g transform="matrix(.017,-0,0,-.017,45.143,14.05)"><path d="M300 -147l-18 -23q-106 71 -159 185.5t-53 254.5v1q0 139 53 252.5t159 186.5l18 -24q-74 -62 -115.5 -173.5t-41.5 -242.5q0 -130 41.5 -242.5t115.5 -174.5z" id="x28"></path></g><g transform="matrix(.017,-0,0,-.017,51.025,14.05)"><path d="M434 650l-5 -28q-60 -8 -75 -20.5t-25 -69.5l-71 -395q-20 -107 -57.5 -168t-106.5 -94q-40 -18 -61 -20l-10 32q77 23 113 108q21 50 41 160l68 377q10 57 -1.5 69.5t-76.5 20.5l6 28h261z" id="x1D43D"></path></g>
<g transform="matrix(.012,-0,0,-.012,58.788,5.887)"><path d="M535 230h-212v-233h-58v233h-213v50h213v210h58v-210h212v-50z" id="x2B"></path></g><g transform="matrix(.012,-0,0,-.012,58.788,18.138)"><path d="M775 650l-6 -28q-60 -6 -81.5 -16t-61.5 -54l-175 -191l125 -243q30 -58 48.5 -71t82.5 -19l-5 -28h-275l7 28l35 4q31 4 37 12t-6 34l-108 216q-140 -165 -177 -219q-16 -22 -10.5 -30.5t41.5 -13.5l22 -3l-7 -28h-244l8 28q52 4 75 15.5t67 52.5q48 46 206 231
l-110 215q-26 51 -44 63t-72 17l6 28h250l-6 -28l-27 -4q-30 -5 -35 -10t3 -27q17 -43 95 -190q70 78 154 185q15 21 10 29.5t-33 12.5l-30 4l5 28h236z" id="x1D44B"></path></g><g transform="matrix(.012,-0,0,-.012,68.211,18.138)"><path d="M368 703l-264 -866h-60l265 866h59z" id="x2F"></path></g><g transform="matrix(.012,-0,0,-.012,73.114,18.138)"><path d="M599 650l-18 -160l-30 -2q-2 70 -11 94q-5 16 -26 24t-75 8h-94q-27 0 -33.5 -6.5t-12.5 -32.5l-43 -226h108q52 0 72.5 5t31 18t28.5 52h28l-40 -200h-28q-1 42 -6 56.5t-24.5 21t-70.5 6.5h-108l-32 -177q-13 -65 1 -80.5t88 -22.5l-8 -28h-279l6 28q61 5 77.5 20.5
t29.5 81.5l72 387l7.5 40.5t1 26.5t-4 18t-15.5 10t-24 6t-37 4l8 28h461z" id="x1D439"></path></g><g transform="matrix(.012,-0,0,-.012,80.515,18.138)"><path d="M368 703l-264 -866h-60l265 866h59z" id="x2F"></path></g>
<g transform="matrix(.012,-0,0,-.012,85.425,18.137)"><path d="M861 0h-263v28q61 5 73 22t12 97l-4 408h-4l-244 -548h-20l-224 534h-2l-15 -285q-5 -98 -3 -166q1 -34 19.5 -45t72.5 -17v-28h-233v28q45 5 62.5 17t22.5 43q11 65 19 183l14 208q3 36 3.5 55t-2.5 37.5t-6.5 25t-16.5 13t-25.5 8t-39.5 4.5v28h164l228 -501l232 501
h170v-28q-68 -6 -81 -21.5t-11 -91.5l9 -362q2 -80 15 -97t78 -22v-28z" id="x4D"></path></g><g transform="matrix(.012,-0,0,-.012,95.92,18.137)"><path d="M614 175l29 -10q-33 -109 -57 -154q-121 -26 -184 -26q-90 0 -160.5 29t-112.5 77t-63.5 105.5t-21.5 119.5q0 157 108 253t277 96q36 0 71.5 -5t69 -13.5t36.5 -8.5q15 -102 20 -150l-29 -8q-20 79 -66.5 114t-128.5 35q-119 0 -187.5 -86t-68.5 -207
q0 -140 73.5 -227.5t188.5 -87.5q73 0 119.5 37.5t86.5 116.5z" id="x43"></path></g>
<g transform="matrix(.017,-0,0,-.017,104.45,14.05)"><path d="M95 130q31 0 61 -30t30 -78q0 -53 -38 -87.5t-93 -51.5l-11 29q77 31 77 85q0 26 -17.5 43t-44.5 24q-4 0 -8.5 6.5t-4.5 17.5q0 18 15 30t34 12z" id="x2C"></path></g><g transform="matrix(.017,-0,0,-.017,111.148,14.05)"><path d="M434 650l-5 -28q-60 -8 -75 -20.5t-25 -69.5l-71 -395q-20 -107 -57.5 -168t-106.5 -94q-40 -18 -61 -20l-10 32q77 23 113 108q21 50 41 160l68 377q10 57 -1.5 69.5t-76.5 20.5l6 28h261z" id="x1D43D"></path></g>
<g transform="matrix(.012,-0,0,-.012,118.925,5.887)"><path d="M535 230h-212v-233h-58v233h-213v50h213v210h58v-210h212v-50z" id="x2B"></path></g>
<g transform="matrix(.012,-0,0,-.012,118.925,18.137)"><path d="M861 0h-263v28q61 5 73 22t12 97l-4 408h-4l-244 -548h-20l-224 534h-2l-15 -285q-5 -98 -3 -166q1 -34 19.5 -45t72.5 -17v-28h-233v28q45 5 62.5 17t22.5 43q11 65 19 183l14 208q3 36 3.5 55t-2.5 37.5t-6.5 25t-16.5 13t-25.5 8t-39.5 4.5v28h164l228 -501l232 501
h170v-28q-68 -6 -81 -21.5t-11 -91.5l9 -362q2 -80 15 -97t78 -22v-28z" id="x4D"></path></g><g transform="matrix(.012,-0,0,-.012,129.432,18.137)"><path d="M614 175l29 -10q-33 -109 -57 -154q-121 -26 -184 -26q-90 0 -160.5 29t-112.5 77t-63.5 105.5t-21.5 119.5q0 157 108 253t277 96q36 0 71.5 -5t69 -13.5t36.5 -8.5q15 -102 20 -150l-29 -8q-20 79 -66.5 114t-128.5 35q-119 0 -187.5 -86t-68.5 -207
q0 -140 73.5 -227.5t188.5 -87.5q73 0 119.5 37.5t86.5 116.5z" id="x43"></path></g>
<g transform="matrix(.017,-0,0,-.017,137.962,14.05)"><path d="M275 270q0 -296 -211 -440l-19 23q75 62 116.5 174t41.5 243t-42 243t-116 173l19 24q211 -144 211 -440z" id="x29"></path></g>
</svg></span></td><td align="center">0.86</td><td align="center">0.92</td><td align="center">1</td></tr><tr><td align="left"><span style="width: 67.137497px;"><svg height="15.4375" id="M180" style="vertical-align:-3.41936pt" version="1.1" viewbox="0 0 67.137497 15.4375" width="67.137497" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,11.112)"><path d="M861 0h-263v28q61 5 73 22t12 97l-4 408h-4l-244 -548h-20l-224 534h-2l-15 -285q-5 -98 -3 -166q1 -34 19.5 -45t72.5 -17v-28h-233v28q45 5 62.5 17t22.5 43q11 65 19 183l14 208q3 36 3.5 55t-2.5 37.5t-6.5 25t-16.5 13t-25.5 8t-39.5 4.5v28h164l228 -501l232 501
h170v-28q-68 -6 -81 -21.5t-11 -91.5l9 -362q2 -80 15 -97t78 -22v-28z" id="x4D"></path></g><g transform="matrix(.017,-0,0,-.017,15.208,11.112)"><path d="M631 18l-4 -26q-24 0 -35 1q-57 3 -92.5 23t-69.5 68q-21 30 -101 160q-13 22 -30.5 31t-54.5 9h-32v-159q0 -63 14 -77.5t74 -19.5v-28h-262v28q62 5 76 19.5t14 77.5v402q0 62 -13 76t-73 19v28h248q109 0 159 -34q67 -42 67 -131q0 -115 -127 -172q27 -52 89 -144
q49 -71 81 -104q31 -36 72 -47zM212 316h48q73 0 108 29q55 42 55 124q0 76 -42 112t-104 36q-41 0 -54 -10q-11 -7 -11 -44v-247z" id="x52"></path></g><g transform="matrix(.017,-0,0,-.017,25.765,11.112)"><path d="M719 650v-28q-43 -2 -62 -15t-22 -44q-6 -47 -6 -169v-403h-31l-426 524h-2v-251q0 -111 6 -169q4 -37 24 -50.5t72 -16.5v-28h-237v28q45 2 64.5 16t23.5 49q6 62 6 171v220q0 54 -3 68.5t-17 32.5q-16 19 -34.5 26.5t-54.5 10.5v28h147l418 -502h3v246q0 117 -7 166
q-4 34 -24.5 47t-73.5 15v28h236z" id="x4E"></path></g><g transform="matrix(.017,-0,0,-.017,38.395,11.112)"><path d="M517 162l29 -4q-18 -100 -35 -158h-482v28q68 5 83.5 20t15.5 77v398q0 66 -13 80t-73 19v28h446q0 -8 8 -144l-29 -4q-11 51 -32 85q-11 17 -35 23t-78 6h-76q-23 0 -28 -5t-5 -27v-228h94q67 0 83 12t25 63h29v-191h-29q-10 54 -26 65.5t-81 11.5h-95v-189
q0 -65 24 -81q19 -13 104 -13q57 0 83 7t43 27q28 33 50 94z" id="x45"></path></g>
<g transform="matrix(.012,-0,0,-.012,48.05,15.187)"><path d="M861 0h-263v28q61 5 73 22t12 97l-4 408h-4l-244 -548h-20l-224 534h-2l-15 -285q-5 -98 -3 -166q1 -34 19.5 -45t72.5 -17v-28h-233v28q45 5 62.5 17t22.5 43q11 65 19 183l14 208q3 36 3.5 55t-2.5 37.5t-6.5 25t-16.5 13t-25.5 8t-39.5 4.5v28h164l228 -501l232 501
h170v-28q-68 -6 -81 -21.5t-11 -91.5l9 -362q2 -80 15 -97t78 -22v-28z" id="x4D"></path></g><g transform="matrix(.012,-0,0,-.012,58.545,15.187)"><path d="M614 175l29 -10q-33 -109 -57 -154q-121 -26 -184 -26q-90 0 -160.5 29t-112.5 77t-63.5 105.5t-21.5 119.5q0 157 108 253t277 96q36 0 71.5 -5t69 -13.5t36.5 -8.5q15 -102 20 -150l-29 -8q-20 79 -66.5 114t-128.5 35q-119 0 -187.5 -86t-68.5 -207
q0 -140 73.5 -227.5t188.5 -87.5q73 0 119.5 37.5t86.5 116.5z" id="x43"></path></g>
</svg></span></td><td align="center">44%</td><td align="center">45%</td><td align="center">0</td></tr><tr><td align="left">Time cost [<i>s</i>]</td><td align="center">2675</td><td align="center">367</td><td align="center">7292</td></tr><tr class="table-tr"><td colspan="4"><hr class="tbody-hr"/></td></tr></table></td></tr></table>

Comparison among the FEM, XFEM and MC method.

Computational and Mathematical Methods in Medicine

tab4

Table 4

Table 4: Extended Finite Element Method with Simplified Spherical Harmonics Approximation for the Forward Model of Optical Molecular Imaging 