Fast and Analytical EAP Approximation from a 4th-Order Tensor

<table>Synthetic dataset experiment. Two fiber bundles intersecting with the DWI signal corrupted by a Rician noise of <svg height="11.6875" id="M187" style="vertical-align:-0.20473pt" version="1.1" viewbox="0 0 66.900002 11.6875" width="66.900002" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,11.363)"><path d="M409 504l-29 -5q-16 60 -44.5 96t-86.5 36q-54 0 -82.5 -32t-28.5 -77q0 -51 30.5 -82.5t96.5 -65.5l35.5 -18t33 -19.5t33.5 -23.5l27 -25.5t24.5 -32t13.5 -36.5t6 -45q0 -80 -62 -134.5t-160 -54.5q-47 0 -98 15q-22 7 -50 21q-8 23 -27 155l30 7q7 -27 18 -52
t30 -51.5t49 -42.5t67 -16q56 0 88 32.5t32 87.5q0 51 -31.5 82t-98.5 67q-80 44 -110 73q-55 53 -55 124q0 75 56 126.5t150 51.5q53 0 126 -23z" id="x53"></path></g><g transform="matrix(.017,-0,0,-.017,7.916,11.363)"><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,20.546,11.363)"><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,35.811,11.363)"><path d="M535 323h-483v50h483v-50zM535 138h-483v50h483v-50z" id="x3D"></path></g><g transform="matrix(.017,-0,0,-.017,50.515,11.363)"><path d="M285 378v-2q65 -13 102 -54.5t37 -97.5q0 -57 -30.5 -104.5t-74 -75t-85.5 -42t-72 -14.5q-31 0 -59.5 11t-40.5 23q-19 18 -16 36q1 16 23 33q13 10 24 0q58 -51 124 -51q55 0 88 40t33 112q0 64 -39 96.5t-88 32.5q-29 0 -64 -11l-6 29q77 25 118 57.5t41 84.5
q0 45 -26.5 69.5t-68.5 24.5q-67 0 -120 -79l-20 20l43 63q51 56 127 56h1q66 0 107 -37t41 -95q0 -42 -31 -71q-22 -23 -68 -54z" id="x33"></path></g><g transform="matrix(.017,-0,0,-.017,58.674,11.363)"><path d="M241 635q53 0 94 -28.5t63.5 -76t33.5 -102.5t11 -116q0 -58 -11 -112.5t-34 -103.5t-63.5 -78.5t-94.5 -29.5t-95 28t-64.5 75t-34.5 102.5t-11 118.5q0 58 11.5 112.5t34.5 103t64.5 78t95.5 29.5zM238 602q-32 0 -55.5 -25t-35.5 -68t-17.5 -91t-5.5 -105
q0 -76 10 -138.5t37 -107.5t69 -45q32 0 55.5 25t35.5 68.5t17.5 91.5t5.5 105t-5.5 105.5t-18 92t-36 68t-56.5 24.5z" id="x30"></path></g>
</svg>. (a) <svg height="16.512501" id="M188" style="vertical-align:-1.32405pt" version="1.1" viewbox="0 0 35.974998 16.512501" width="35.974998" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,1.212,12.312)"><path d="M1974 560l-5 -28l-982 111l-982 -111l-5 28l163.5 24.5l176 26.5l159 24.5t171.5 28l150.5 27t155.5 30.5h22q263 -55 976 -161z" id="xF5C09"></path></g><g transform="matrix(.017,-0,0,-.017,.062,16.2)"><path d="M673 0h-245v28q47 6 58 16t1 39l-55 154h-213q-31 -87 -36 -129q-16 -47 -5 -61t68 -19v-28h-226v28q50 7 68.5 23.5t43.5 80.5l214 524l24 9l193 -536q23 -63 41.5 -79.5t68.5 -21.5v-28zM418 280l-81 255h-2q-74 -197 -101 -255h184z" id="x41"></path></g><g transform="matrix(.017,-0,0,-.017,11.809,16.2)"><path d="M43 650h280q194 0 282 -85q85 -82 85 -215q0 -173 -123 -267q-107 -83 -297 -83h-236v27q65 5 79.5 19.5t14.5 77.5v402q0 45 -6 62.5t-21.5 24t-57.5 9.5v28zM213 547v-416q0 -55 24.5 -75t81.5 -20q133 0 203 86t70 214q0 80 -30 143t-88 97q-66 41 -170 41
q-62 0 -78 -13q-13 -9 -13 -57z" id="x44"></path></g><g transform="matrix(.017,-0,0,-.017,24.592,16.2)"><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> from 4th-order tensors. (b) Analytically approximated EAP, <svg height="22.9125" id="M189" style="vertical-align:-4.9589pt" version="1.1" viewbox="0 0 30.612499 22.9125" width="30.612499" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,16.675)"><path d="M46 650h247q112 0 169 -41q64 -49 64 -141q0 -91 -56.5 -143t-138.5 -62q-23 -3 -45 -1l-74 19v-157q0 -63 15.5 -77t83.5 -19v-28h-271v28q60 5 74 20t14 77v400q0 65 -12.5 79t-69.5 18v28zM212 559v-245q31 -12 73 -12q59 0 102 38t43 124q0 77 -42.5 115.5
t-108.5 38.5q-44 0 -56 -12q-11 -8 -11 -47z" id="x50"></path></g><g transform="matrix(.017,-0,0,-.017,9.48,16.675)"><path d="M181 342h2q63 107 121 107q24 0 41 -16t17 -34q0 -34 -32 -49q-17 -7 -26 2q-21 20 -43 20q-20 0 -42 -22t-38 -63v-183q0 -50 12.5 -62t69.5 -16v-26h-230v26q47 5 58 17t11 61v207q0 49 -9 60.5t-53 16.5v23q74 12 141 40v-109z" id="x72"></path></g>
<g transform="matrix(.012,-0,0,-.012,15.788,8.513)"><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(.012,-0,0,-.012,19.904,8.513)"><path d="M447 623l8 -12l-283 -613l-74 -10l-7 11q174 283 297 551h-216q-48 0 -62.5 -12t-33.5 -63h-29q10 60 18 148h382z" id="x37"></path></g><g transform="matrix(.012,-0,0,-.012,25.616,8.513)"><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><g transform="matrix(.012,-0,0,-.012,15.788,21.45)"><path d="M456 178h-96v-72q0 -51 12.5 -62.5t72.5 -16.5v-27h-256v27q65 5 78 17t13 62v72h-260v28q182 271 300 426h40v-407h96v-47zM280 225v295h-2q-107 -148 -196 -295h198z" id="x34"></path></g><g transform="matrix(.012,-0,0,-.012,21.499,21.45)"><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(.012,-0,0,-.012,24.212,21.45)"><path d="M456 178h-96v-72q0 -51 12.5 -62.5t72.5 -16.5v-27h-256v27q65 5 78 17t13 62v72h-260v28q182 271 300 426h40v-407h96v-47zM280 225v295h-2q-107 -148 -196 -295h198z" id="x34"></path></g>
</svg>.</table>

International Journal of Biomedical Imaging

Fast and Analytical EAP Approximation from a 4th-Order Tensor

Figure 2