Reconstruction of the Magnetic Particle Imaging System Matrix Using Symmetries and Compressed Sensing

<div>Several system-function components of the <svg height="6.1673pt" id="M38" style="vertical-align:-0.2063904pt" version="1.1" viewbox="-0.0498162 -5.96091 7.39387 6.1673" width="7.39387pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M536 404C536 423 520 448 491 448C445 448 398 404 308 283L286 338C255 416 242 448 217 448C182 448 138 402 92 341L111 321C149 368 169 378 178 378C188 378 198 363 214 320L254 212C185 117 138 65 107 65C98 65 85 69 82 75C79 82 73 86 65 86C44 86 23 60 23 39C23 7 44 -12 71 -12C119 -12 168 33 265 177L306 61C321 17 347 -12 373 -12C413 -12 465 33 507 96L491 119C459 84 432 60 413 60C395 60 378 92 358 148L321 250C341 279 369 310 389 332C417 363 439 382 456 382C466 382 475 376 481 368C486 361 492 358 496 358C513 358 536 381 536 404Z" id="g113-121"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="559" vert-adv-y="559"></glyph.data></g></svg>-channel at frequencies <svg height="12.7178pt" id="M39" style="vertical-align:-3.42947pt" version="1.1" viewbox="-0.0498162 -9.28833 171.021 12.7178" width="171.021pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M619 670C619 686 593 712 555 712S459 686 410 634S335 504 320 430H250L219 400L222 388H312L258 73C223 -133 201 -166 187 -180C175 -191 158 -199 140 -199C123 -199 88 -188 74 -172C68 -166 63 -164 54 -171C38 -185 23 -201 23 -215C23 -236 60 -261 93 -261C122 -261 161 -247 207 -200C268 -138 300 -71 337 94C365 220 376 277 394 387L501 399L521 430H401C432 623 464 665 501 665C524 665 544 651 567 627C577 617 583 618 592 625C601 631 619 651 619 670Z" id="g113-103"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="642" vert-adv-y="642"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,11.984,0)"><path d="M448 1V51H364C248 51 153 129 140 230H448V280H140C153 381 248 459 364 459H448V509H365C208 509 80 395 80 255S208 1 365 1H448Z" id="g117-173"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="528" vert-adv-y="528"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,22.48,0)"><path d="M293 -169V-141C218 -131 209 -90 209 -44C209 19 222 85 222 152C222 207 198 255 139 269V273C199 286 222 334 222 388C222 454 209 523 209 588C209 632 218 671 293 681V709C234 709 148 695 148 577C147 513 155 438 155 372C155 337 152 291 64 285V256C152 250 155 204 155 169C156 105 148 31 148 -41C148 -157 234 -169 293 -169Z" id="g113-124"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="347" vert-adv-y="347"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,26.991,0)"><path d="M384 0V27C293 34 287 42 287 114V635C232 613 172 594 109 583V559L157 557C201 555 205 550 205 499V114C205 42 199 34 109 27V0H384Z" id="g113-50"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,33.231,0)"><path d="M249 635C141 635 70 555 70 471C70 401 114 353 179 316C143 294 106 267 90 252C68 231 45 202 45 157C45 50 130 -12 237 -12C322 -12 435 52 435 169C435 256 372 304 303 343C349 374 375 398 383 407C401 429 411 458 411 487C411 569 344 635 249 635ZM238 603C285 603 337 567 337 482C337 422 310 385 276 358C205 393 145 426 145 500C145 552 179 603 238 603ZM248 20C183 20 125 70 125 163C125 218 158 268 206 300C284 261 355 217 355 143C355 66 308 20 248 20Z" id="g113-57"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,39.471,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z" id="g113-49"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,47.96,0)"><path d="M513 0V26C475 28 448 43 413 81C379 116 277 234 258 264C301 314 319 334 344 359C383 396 404 407 473 418V443L269 434V408C317 403 320 395 299 364C280 337 250 298 228 276C213 260 195 247 165 241V712C126 699 69 682 21 677V653C82 648 87 642 87 576V104C87 40 81 32 14 26V0H237V26C173 32 165 38 165 104V215C171 217 176 218 181 218C189 218 198 212 212 196C267 129 330 49 376 -2L513 0Z" id="g121-105"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="496" vert-adv-y="496"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,54.63,0)"><path d="M729 650H468V622C546 615 554 610 554 527V365H213V527C213 610 221 616 297 622V650H38V622C120 616 128 611 128 527V123C128 39 120 34 40 28V0H303V28C221 34 213 40 213 123V322H554V123C554 39 546 34 460 28V0H728V28C647 34 639 39 639 123V527C639 610 647 615 729 622V650Z" id="g121-70"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="766" vert-adv-y="766"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,64.458,0)"><path d="M387 138L374 114C336 54 316 33 240 33H141C223 167 303 295 385 423L379 437H149C112 437 100 443 94 458H73C65 415 50 352 39 314L64 305L78 332C106 388 119 403 191 403C214 403 246 403 280 403C198 270 117 142 31 15L38 0H367C382 40 399 86 413 125L387 138Z" id="g121-120"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="420" vert-adv-y="420"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,69.918,0)"><path d="M95 130C70 130 46 113 46 88C46 72 54 64 59 64C93 55 121 33 121 -3C121 -41 93 -68 44 -88L55 -117C117 -98 186 -56 186 22C186 91 131 130 95 130Z" id="g113-45"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="228" vert-adv-y="228"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,75.061,0)"><path d="M285 378C315 398 338 416 353 432C373 451 384 474 384 503C384 579 325 635 236 635H235C182 635 136 610 108 579L65 516L85 496C110 533 150 575 205 575C258 575 300 543 300 481C300 407 232 369 141 339L147 310C163 315 188 321 211 321C268 321 338 284 338 192C338 94 288 40 217 40C160 40 119 68 93 91C85 98 77 97 69 91C60 84 47 71 46 58C44 46 48 35 62 22C75 10 116 -12 162 -12C234 -12 424 62 424 224C424 297 373 359 285 376V378Z" id="g113-52"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,81.301,0)"><path d="M153 550H386L412 615L406 623H120L82 318C104 327 142 338 184 338C294 338 347 275 347 187C347 112 305 39 221 39C160 39 119 71 97 89C88 97 80 96 71 90C59 80 50 67 49 57C48 45 52 36 66 23C80 9 123 -12 169 -12C221 -11 288 15 342 59C403 109 431 165 431 225C431 308 366 395 238 395C212 395 165 379 127 364L153 550Z" id="g113-54"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,87.541,0)"><path d="M244 635C114 635 38 519 38 422C38 317 111 240 217 240C236 240 255 244 277 256L345 292C311 140 203 39 59 15L64 -15C89 -15 150 -5 204 17C339 72 440 202 440 386C440 521 368 635 244 635ZM228 602C326 602 352 479 352 390C352 370 351 347 348 324C327 308 293 296 258 296C174 296 124 369 124 458C124 517 152 602 228 602Z" id="g113-58"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,96.03,0)"><path d="M513 0V26C475 28 448 43 413 81C379 116 277 234 258 264C301 314 319 334 344 359C383 396 404 407 473 418V443L269 434V408C317 403 320 395 299 364C280 337 250 298 228 276C213 260 195 247 165 241V712C126 699 69 682 21 677V653C82 648 87 642 87 576V104C87 40 81 32 14 26V0H237V26C173 32 165 38 165 104V215C171 217 176 218 181 218C189 218 198 212 212 196C267 129 330 49 376 -2L513 0Z" id="g121-105"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="496" vert-adv-y="496"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,102.7,0)"><path d="M729 650H468V622C546 615 554 610 554 527V365H213V527C213 610 221 616 297 622V650H38V622C120 616 128 611 128 527V123C128 39 120 34 40 28V0H303V28C221 34 213 40 213 123V322H554V123C554 39 546 34 460 28V0H728V28C647 34 639 39 639 123V527C639 610 647 615 729 622V650Z" id="g121-70"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="766" vert-adv-y="766"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,112.527,0)"><path d="M387 138L374 114C336 54 316 33 240 33H141C223 167 303 295 385 423L379 437H149C112 437 100 443 94 458H73C65 415 50 352 39 314L64 305L78 332C106 388 119 403 191 403C214 403 246 403 280 403C198 270 117 142 31 15L38 0H367C382 40 399 86 413 125L387 138Z" id="g121-120"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="420" vert-adv-y="420"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,117.987,0)"><path d="M95 130C70 130 46 113 46 88C46 72 54 64 59 64C93 55 121 33 121 -3C121 -41 93 -68 44 -88L55 -117C117 -98 186 -56 186 22C186 91 131 130 95 130Z" id="g113-45"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="228" vert-adv-y="228"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,123.13,0)"><path d="M456 178V225H360V632H320C217 496 115 347 20 206V178H280V106C280 40 276 34 189 27V0H445V27C364 34 360 39 360 106V178H456ZM280 225H82C149 335 214 431 278 520H280V225Z" id="g113-53"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,129.371,0)"><path d="M384 0V27C293 34 287 42 287 114V635C232 613 172 594 109 583V559L157 557C201 555 205 550 205 499V114C205 42 199 34 109 27V0H384Z" id="g113-50"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,135.611,0)"><path d="M412 140C382 77 369 73 315 73H129L270 222C362 320 402 379 402 466C402 571 322 635 234 635C177 635 130 609 99 576L42 495L64 475C90 514 133 568 201 568C274 568 318 519 318 435C318 349 255 267 193 193C144 135 87 78 32 23V0H405C417 45 427 89 440 131L412 140Z" id="g113-51"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,144.1,0)"><path d="M513 0V26C475 28 448 43 413 81C379 116 277 234 258 264C301 314 319 334 344 359C383 396 404 407 473 418V443L269 434V408C317 403 320 395 299 364C280 337 250 298 228 276C213 260 195 247 165 241V712C126 699 69 682 21 677V653C82 648 87 642 87 576V104C87 40 81 32 14 26V0H237V26C173 32 165 38 165 104V215C171 217 176 218 181 218C189 218 198 212 212 196C267 129 330 49 376 -2L513 0Z" id="g121-105"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="496" vert-adv-y="496"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,150.769,0)"><path d="M729 650H468V622C546 615 554 610 554 527V365H213V527C213 610 221 616 297 622V650H38V622C120 616 128 611 128 527V123C128 39 120 34 40 28V0H303V28C221 34 213 40 213 123V322H554V123C554 39 546 34 460 28V0H728V28C647 34 639 39 639 123V527C639 610 647 615 729 622V650Z" id="g121-70"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="766" vert-adv-y="766"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,160.597,0)"><path d="M387 138L374 114C336 54 316 33 240 33H141C223 167 303 295 385 423L379 437H149C112 437 100 443 94 458H73C65 415 50 352 39 314L64 305L78 332C106 388 119 403 191 403C214 403 246 403 280 403C198 270 117 142 31 15L38 0H367C382 40 399 86 413 125L387 138Z" id="g121-120"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="420" vert-adv-y="420"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,166.057,0)"><path d="M283 255V284C195 290 192 337 192 375C192 436 200 508 200 580C200 696 116 709 54 709V681C127 671 139 634 139 588C138 524 125 452 125 387C125 333 149 286 209 272V268C148 253 125 206 125 151C126 87 139 16 139 -47C139 -92 127 -131 54 -141V-169C115 -169 200 -152 200 -41C200 32 192 104 192 164C192 202 195 249 283 255Z" id="g113-126"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="347" vert-adv-y="347"></glyph.data></g></svg>, which are obtained with compressed sensing at different sampling rates. Compressed sensing is performed with the DCT in combination with the FISTA algorithm.</div>

Advances in Mathematical Physics

fig2

Figure 2

Figure 2: Reconstruction of the Magnetic Particle Imaging System Matrix Using Symmetries and Compressed Sensing 

数学物理进步

重构磁粒映射系统矩阵使用对称和压缩

图2