A MATLAB-Based Boundary Data Simulator for Studying the Resistivity Reconstruction Using Neighbouring Current Pattern

<table>Image reconstruction of multiple inhomogeneities (<svg height="16.299999" id="M339" style="vertical-align:-3.3907pt" version="1.1" viewbox="0 0 95.525002 16.299999" width="95.525002" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,12.013)"><path d="M548 455q-32 -77 -74 -77q-37 0 -103 10l-3 -2q36 -30 48.5 -62t12.5 -80q0 -103 -75 -179.5t-175 -76.5q-70 0 -113 45t-43 128q0 115 83.5 200.5t209.5 85.5q31 0 77.5 -4t67.5 -4q36 0 62 30zM350 274q0 54 -17 86q-19 35 -57 35q-50 0 -90 -37.5t-59 -91.5t-19 -109
q0 -61 24.5 -96.5t65.5 -35.5q51 0 87.5 46.5t50.5 100.5t14 102z" id="x1D70E"></path></g>
<g transform="matrix(.012,-0,0,-.012,9.663,16.088)"><path d="M452 333q0 -82 -46 -164t-116 -128q-81 -53 -158 -53q-26 0 -51 13t-37 32q-26 41 -6 134l91 431q7 40 2.5 47t-34.5 7h-33l2 25q36 4 72.5 13t58.5 15.5t28 6.5q10 0 4 -29l-91 -391h2q65 77 130 116.5t105 39.5q36 0 56.5 -32t20.5 -83zM365 316q0 76 -35 76
q-31 0 -93.5 -49t-107.5 -110q-12 -37 -17 -70q-9 -65 12 -96.5t59 -31.5q31 0 56 14q53 29 89.5 105t36.5 162z" id="x1D44F"></path></g>
<g transform="matrix(.017,-0,0,-.017,20.588,12.013)"><path d="M535 323h-483v50h483v-50zM535 138h-483v50h483v-50z" id="x3D"></path></g><g transform="matrix(.017,-0,0,-.017,35.291,12.013)"><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><g transform="matrix(.017,-0,0,-.017,43.451,12.013)"><path d="M113 -12q-24 0 -39.5 16t-15.5 42q0 24 16 40.5t40 16.5t40 -16.5t16 -40.5q0 -26 -16 -42t-41 -16z" id="x2E"></path></g><g transform="matrix(.017,-0,0,-.017,47.327,12.013)"><path d="M412 140l28 -9q0 -2 -35 -131h-373v23q112 112 161 170q59 70 92 127t33 115q0 63 -31 98t-86 35q-75 0 -137 -93l-22 20l57 81q55 59 135 59q69 0 118.5 -46.5t49.5 -122.5q0 -62 -29.5 -114t-102.5 -130l-141 -149h186q42 0 58.5 10.5t38.5 56.5z" id="x32"></path></g><g transform="matrix(.017,-0,0,-.017,55.486,12.013)"><path d="M384 0h-275v27q67 5 81.5 18.5t14.5 68.5v385q0 38 -7.5 47.5t-40.5 10.5l-48 2v24q85 15 178 52v-521q0 -55 14.5 -68.5t82.5 -18.5v-27z" id="x31"></path></g><g transform="matrix(.017,-0,0,-.017,66.467,12.013)"><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,74.525,12.013)"><path d="M368 703l-264 -866h-60l265 866h59z" id="x2F"></path></g><g transform="matrix(.017,-0,0,-.017,81.528,12.013)"><path d="M797 0h-214v26q46 5 58 17.5t12 63.5v167q0 116 -88 116q-48 0 -112 -52q3 -18 3 -42v-195q0 -46 11 -58t53 -17v-26h-214v26q48 5 59.5 17t11.5 60v172q0 115 -88 115q-47 0 -109 -52v-235q0 -47 10 -59.5t48 -17.5v-26h-211v26q50 5 62.5 17t12.5 60v207q0 48 -9 59.5
t-52 18.5v23q79 14 139 41v-81q3 2 29 20.5t44 29.5q45 28 81 28t66 -22t46 -61q25 20 83 57q39 26 76 26q59 0 93 -40t34 -113v-193q0 -49 10.5 -60.5t54.5 -16.5v-26z" id="x6D"></path></g>
</svg> and <svg height="16.299999" id="M340" style="vertical-align:-3.3907pt" version="1.1" viewbox="0 0 101.35 16.299999" width="101.35" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,12.013)"><path d="M548 455q-32 -77 -74 -77q-37 0 -103 10l-3 -2q36 -30 48.5 -62t12.5 -80q0 -103 -75 -179.5t-175 -76.5q-70 0 -113 45t-43 128q0 115 83.5 200.5t209.5 85.5q31 0 77.5 -4t67.5 -4q36 0 62 30zM350 274q0 54 -17 86q-19 35 -57 35q-50 0 -90 -37.5t-59 -91.5t-19 -109
q0 -61 24.5 -96.5t65.5 -35.5q51 0 87.5 46.5t50.5 100.5t14 102z" id="x1D70E"></path></g>
<g transform="matrix(.012,-0,0,-.012,9.663,16.088)"><path d="M244 607q0 -25 -15.5 -43t-37.5 -18q-19 0 -32 13t-13 35q0 21 15 41t39 20q20 0 32 -14t12 -34zM222 91q-29 -33 -79 -68t-75 -35q-13 0 -19 7.5t-6 31t10 65.5l62 253q5 26 -1 26q-21 0 -72 -43l-13 24q43 40 91 68t71 28q30 0 10 -78l-71 -274q-8 -30 3 -30
q16 0 76 48z" id="x1D456"></path></g>
<g transform="matrix(.017,-0,0,-.017,18.262,12.013)"><path d="M535 323h-483v50h483v-50zM535 138h-483v50h483v-50z" id="x3D"></path></g><g transform="matrix(.017,-0,0,-.017,32.949,12.013)"><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><g transform="matrix(.017,-0,0,-.017,41.109,12.013)"><path d="M113 -12q-24 0 -39.5 16t-15.5 42q0 24 16 40.5t40 16.5t40 -16.5t16 -40.5q0 -26 -16 -42t-41 -16z" id="x2E"></path></g><g transform="matrix(.017,-0,0,-.017,44.985,12.013)"><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><g transform="matrix(.017,-0,0,-.017,53.144,12.013)"><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><g transform="matrix(.017,-0,0,-.017,61.303,12.013)"><path d="M153 550l-26 -186q79 31 111 31q90 0 141.5 -51t51.5 -119q0 -93 -89 -166q-85 -69 -173 -71q-32 0 -61.5 11.5t-41.5 23.5q-18 17 -17 34q2 16 22 33q14 9 26 -1q61 -50 124 -50q60 0 93 43.5t33 104.5q0 69 -41.5 110t-121.5 41q-53 0 -102 -20l38 305h286l6 -8
l-26 -65h-233z" id="x35"></path></g><g transform="matrix(.017,-0,0,-.017,72.302,12.013)"><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,80.359,12.013)"><path d="M368 703l-264 -866h-60l265 866h59z" id="x2F"></path></g><g transform="matrix(.017,-0,0,-.017,87.362,12.013)"><path d="M797 0h-214v26q46 5 58 17.5t12 63.5v167q0 116 -88 116q-48 0 -112 -52q3 -18 3 -42v-195q0 -46 11 -58t53 -17v-26h-214v26q48 5 59.5 17t11.5 60v172q0 115 -88 115q-47 0 -109 -52v-235q0 -47 10 -59.5t48 -17.5v-26h-211v26q50 5 62.5 17t12.5 60v207q0 48 -9 59.5
t-52 18.5v23q79 14 139 41v-81q3 2 29 20.5t44 29.5q45 28 81 28t66 -22t46 -61q25 20 83 57q39 26 76 26q59 0 93 -40t34 -113v-193q0 -49 10.5 -60.5t54.5 -16.5v-26z" id="x6D"></path></g>
</svg>): (a) simulated domain with two circular objects (<svg height="15.1375" id="M341" style="vertical-align:-3.3907pt" version="1.1" viewbox="0 0 77.587502 15.1375" width="77.587502" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,10.862)"><path d="M393 379q-9 -16 -28 -29q-15 -10 -23 -2q-19 19 -36 19q-21 0 -52 -38q-57 -72 -82 -126l-40 -197q-23 -3 -75 -18l-7 7q49 196 74 335q7 43 -2 43q-7 0 -30 -14.5t-47 -37.5l-16 23q37 42 82 73t67 31q41 0 15 -113l-11 -50h4q41 71 85 117t77 46q29 0 45 -26
q13 -21 0 -43z" id="x1D45F"></path></g>
<g transform="matrix(.012,-0,0,-.012,7.238,14.938)"><path d="M244 607q0 -25 -15.5 -43t-37.5 -18q-19 0 -32 13t-13 35q0 21 15 41t39 20q20 0 32 -14t12 -34zM222 91q-29 -33 -79 -68t-75 -35q-13 0 -19 7.5t-6 31t10 65.5l62 253q5 26 -1 26q-21 0 -72 -43l-13 24q43 40 91 68t71 28q30 0 10 -78l-71 -274q-8 -30 3 -30
q16 0 76 48z" id="x1D456"></path></g>
<g transform="matrix(.017,-0,0,-.017,15.825,10.862)"><path d="M535 323h-483v50h483v-50zM535 138h-483v50h483v-50z" id="x3D"></path></g><g transform="matrix(.017,-0,0,-.017,30.512,10.862)"><path d="M412 140l28 -9q0 -2 -35 -131h-373v23q112 112 161 170q59 70 92 127t33 115q0 63 -31 98t-86 35q-75 0 -137 -93l-22 20l57 81q55 59 135 59q69 0 118.5 -46.5t49.5 -122.5q0 -62 -29.5 -114t-102.5 -130l-141 -149h186q42 0 58.5 10.5t38.5 56.5z" id="x32"></path></g><g transform="matrix(.017,-0,0,-.017,38.671,10.862)"><path d="M153 550l-26 -186q79 31 111 31q90 0 141.5 -51t51.5 -119q0 -93 -89 -166q-85 -69 -173 -71q-32 0 -61.5 11.5t-41.5 23.5q-18 17 -17 34q2 16 22 33q14 9 26 -1q61 -50 124 -50q60 0 93 43.5t33 104.5q0 69 -41.5 110t-121.5 41q-53 0 -102 -20l38 305h286l6 -8
l-26 -65h-233z" id="x35"></path></g><g transform="matrix(.017,-0,0,-.017,49.67,10.862)"><path d="M797 0h-214v26q46 5 58 17.5t12 63.5v167q0 116 -88 116q-48 0 -112 -52q3 -18 3 -42v-195q0 -46 11 -58t53 -17v-26h-214v26q48 5 59.5 17t11.5 60v172q0 115 -88 115q-47 0 -109 -52v-235q0 -47 10 -59.5t48 -17.5v-26h-211v26q50 5 62.5 17t12.5 60v207q0 48 -9 59.5
t-52 18.5v23q79 14 139 41v-81q3 2 29 20.5t44 29.5q45 28 81 28t66 -22t46 -61q25 20 83 57q39 26 76 26q59 0 93 -40t34 -113v-193q0 -49 10.5 -60.5t54.5 -16.5v-26z" id="x6D"></path></g><g transform="matrix(.017,-0,0,-.017,63.591,10.862)"><path d="M797 0h-214v26q46 5 58 17.5t12 63.5v167q0 116 -88 116q-48 0 -112 -52q3 -18 3 -42v-195q0 -46 11 -58t53 -17v-26h-214v26q48 5 59.5 17t11.5 60v172q0 115 -88 115q-47 0 -109 -52v-235q0 -47 10 -59.5t48 -17.5v-26h-211v26q50 5 62.5 17t12.5 60v207q0 48 -9 59.5
t-52 18.5v23q79 14 139 41v-81q3 2 29 20.5t44 29.5q45 28 81 28t66 -22t46 -61q25 20 83 57q39 26 76 26q59 0 93 -40t34 -113v-193q0 -49 10.5 -60.5t54.5 -16.5v-26z" id="x6D"></path></g>
</svg>, 180° apart), (b) reconstructed image of the domain shown in (a), (c) simulated domain with three circular objects (<svg height="15.1375" id="M342" style="vertical-align:-3.3907pt" version="1.1" viewbox="0 0 77.587502 15.1375" width="77.587502" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,10.862)"><path d="M393 379q-9 -16 -28 -29q-15 -10 -23 -2q-19 19 -36 19q-21 0 -52 -38q-57 -72 -82 -126l-40 -197q-23 -3 -75 -18l-7 7q49 196 74 335q7 43 -2 43q-7 0 -30 -14.5t-47 -37.5l-16 23q37 42 82 73t67 31q41 0 15 -113l-11 -50h4q41 71 85 117t77 46q29 0 45 -26
q13 -21 0 -43z" id="x1D45F"></path></g>
<g transform="matrix(.012,-0,0,-.012,7.238,14.938)"><path d="M244 607q0 -25 -15.5 -43t-37.5 -18q-19 0 -32 13t-13 35q0 21 15 41t39 20q20 0 32 -14t12 -34zM222 91q-29 -33 -79 -68t-75 -35q-13 0 -19 7.5t-6 31t10 65.5l62 253q5 26 -1 26q-21 0 -72 -43l-13 24q43 40 91 68t71 28q30 0 10 -78l-71 -274q-8 -30 3 -30
q16 0 76 48z" id="x1D456"></path></g>
<g transform="matrix(.017,-0,0,-.017,15.825,10.862)"><path d="M535 323h-483v50h483v-50zM535 138h-483v50h483v-50z" id="x3D"></path></g><g transform="matrix(.017,-0,0,-.017,30.512,10.862)"><path d="M412 140l28 -9q0 -2 -35 -131h-373v23q112 112 161 170q59 70 92 127t33 115q0 63 -31 98t-86 35q-75 0 -137 -93l-22 20l57 81q55 59 135 59q69 0 118.5 -46.5t49.5 -122.5q0 -62 -29.5 -114t-102.5 -130l-141 -149h186q42 0 58.5 10.5t38.5 56.5z" id="x32"></path></g><g transform="matrix(.017,-0,0,-.017,38.671,10.862)"><path d="M153 550l-26 -186q79 31 111 31q90 0 141.5 -51t51.5 -119q0 -93 -89 -166q-85 -69 -173 -71q-32 0 -61.5 11.5t-41.5 23.5q-18 17 -17 34q2 16 22 33q14 9 26 -1q61 -50 124 -50q60 0 93 43.5t33 104.5q0 69 -41.5 110t-121.5 41q-53 0 -102 -20l38 305h286l6 -8
l-26 -65h-233z" id="x35"></path></g><g transform="matrix(.017,-0,0,-.017,49.67,10.862)"><path d="M797 0h-214v26q46 5 58 17.5t12 63.5v167q0 116 -88 116q-48 0 -112 -52q3 -18 3 -42v-195q0 -46 11 -58t53 -17v-26h-214v26q48 5 59.5 17t11.5 60v172q0 115 -88 115q-47 0 -109 -52v-235q0 -47 10 -59.5t48 -17.5v-26h-211v26q50 5 62.5 17t12.5 60v207q0 48 -9 59.5
t-52 18.5v23q79 14 139 41v-81q3 2 29 20.5t44 29.5q45 28 81 28t66 -22t46 -61q25 20 83 57q39 26 76 26q59 0 93 -40t34 -113v-193q0 -49 10.5 -60.5t54.5 -16.5v-26z" id="x6D"></path></g><g transform="matrix(.017,-0,0,-.017,63.591,10.862)"><path d="M797 0h-214v26q46 5 58 17.5t12 63.5v167q0 116 -88 116q-48 0 -112 -52q3 -18 3 -42v-195q0 -46 11 -58t53 -17v-26h-214v26q48 5 59.5 17t11.5 60v172q0 115 -88 115q-47 0 -109 -52v-235q0 -47 10 -59.5t48 -17.5v-26h-211v26q50 5 62.5 17t12.5 60v207q0 48 -9 59.5
t-52 18.5v23q79 14 139 41v-81q3 2 29 20.5t44 29.5q45 28 81 28t66 -22t46 -61q25 20 83 57q39 26 76 26q59 0 93 -40t34 -113v-193q0 -49 10.5 -60.5t54.5 -16.5v-26z" id="x6D"></path></g>
</svg>, 120° apart), and (d) reconstructed image of the domain shown in (c).</table>

Journal of Medical Engineering

fig8

Figure 8

Figure 8: A MATLAB-Based Boundary Data Simulator for Studying the Resistivity Reconstruction Using Neighbouring Current Pattern 