Improved Ultrafast Power Doppler Imaging by Using Spatiotemporal Non-Local Means Filtering

Lijie Huang,Jingke Zhang,Xingyue Wei,Linkai Jing,Qiong He,Xia Xie,Guihuai Wang,Jianwen Luo

DOI: https://doi.org/10.1109/tuffc.2022.3158611

2022-05-01

Abstract:The change of microvasculature is associated with the occurrence and development of many diseases. Ultrafast power Doppler imaging (uPDI) is an emerging technology for the visualization of microvessels due to the development of ultrafast plane wave (PW) imaging and advanced clutter filters. However, the low signal-to-noise ratio (SNR) caused by unfocused transmit of PW imaging deteriorates the subsequent imaging of microvasculature. Nonlocal means (NLM) filtering has been demonstrated to be effective in the denoising of both natural and medical images, including ultrasound power Doppler images. However, the feasibility and performance of applying an NLM filter on the ultrasound radio frequency (RF) data have not been investigated so far. In this study, we propose to apply an NLM filter on the spatiotemporal domain of clutter filtered blood flow RF data (St-NLM) to improve the quality of uPDI. Experiments were conducted to compare the proposed method with three different methods (under various similarity window sizes), including conventional uPDI without NLM filtering (Non-NLM), NLM filtering on the obtained power Doppler images (PD-NLM), and NLM filtering on the spatial domain of clutter filtered blood flow RF data (S-NLM). Phantom experiments, in vivo contrast-enhanced human spinal cord tumor experiments, and in vivo contrast-free human liver experiments were performed to demonstrate the superiority of the proposed St-NLM method over the other three methods. Qualitative and quantitative results show that the proposed St-NLM method can effectively suppress the background noise, improve the contrast between vessels and background, and preserve the details of small vessels at the same time. In the human liver study, the proposed St-NLM method achieves 31.05-, 24.49-, and 11.15-dB higher contrast-to-noise ratios (CNRs) and 36.86-, 36.86-, and 15.22-dB lower noise powers than Non-NLM, PD-NLM, and S-NLM, respectively. In the human spina- cord tumor, the full-width at half-maximums (FWHMs) of vessel cross Section are 76, 201, and <span class="mjpage"><svg xmlns:xlink="http://www.w3.org/1999/xlink" width="12.99ex" height="2.176ex" style="vertical-align: -0.338ex;" viewBox="0 -791.3 5593 936.9" role="img" focusable="false" xmlns="http://www.w3.org/2000/svg"><g stroke="currentColor" fill="currentColor" stroke-width="0" transform="matrix(1 0 0 -1 0 0)"> <use xlink:href="#MJMAIN-31"></use> <use xlink:href="#MJMAIN-39" x="500" y="0"></use> <use xlink:href="#MJMAIN-37" x="1001" y="0"></use> <use xlink:href="#MJMATHI-6D" x="1501" y="0"></use> <use xlink:href="#MJMATHI-75" x="2380" y="0"></use> <use xlink:href="#MJMATHI-74" x="2952" y="0"></use> <use xlink:href="#MJMATHI-65" x="3314" y="0"></use> <use xlink:href="#MJMATHI-78" x="3780" y="0"></use> <use xlink:href="#MJMATHI-74" x="4353" y="0"></use> <use xlink:href="#MJMATHI-6D" x="4714" y="0"></use></g></svg></span> for St-NLM, Non-NLM, and S-NLM, respectively. The proposed St-NLM method can enhance the microvascular visualization in uPDI and has the potential for the diagnosis of many microvessel-change-related diseases.<svg xmlns="http://www.w3.org/2000/svg" style="display: none;"><defs id="MathJax_SVG_glyphs"><path stroke-width="1" id="MJMAIN-31" d="M213 578L200 573Q186 568 160 563T102 556H83V602H102Q149 604 189 617T245 641T273 663Q275 666 285 666Q294 666 302 660V361L303 61Q310 54 315 52T339 48T401 46H427V0H416Q395 3 257 3Q121 3 100 0H88V46H114Q136 46 152 46T177 47T193 50T201 52T207 57T213 61V578Z"></path><path stroke-width="1" id="MJMAIN-39" d="M352 287Q304 211 232 211Q154 211 104 270T44 396Q42 412 42 436V444Q42 537 111 606Q171 666 243 666Q245 666 249 666T257 665H261Q273 665 286 663T323 651T370 619T413 560Q456 472 456 334Q456 194 396 97Q361 41 312 10T208 -22Q147 -22 108 7T68 93T121 149Q143 149 158 135T173 96Q173 78 164 65T148 49T135 44L131 43Q131 41 138 37T164 27T206 22H212Q272 22 313 86Q352 142 352 280V287ZM244 248Q292 248 321 297T351 430Q351 508 343 542Q341 552 337 562T323 588T293 615T246 625Q208 625 181 598Q160 576 154 546T147 441Q147 358 152 329T172 282Q197 248 244 248Z"></path><path stroke-width="1" id="MJMAIN-37" d="M55 458Q56 460 72 567L88 674Q88 676 108 676H128V672Q128 662 143 655T195 646T364 644H485V605L417 512Q408 500 387 472T360 435T339 403T319 367T305 330T292 284T284 230T278 162T275 80Q275 66 275 52T274 28V19Q270 2 255 -10T221 -22Q210 -22 200 -19T179 0T168 40Q168 198 265 368Q285 400 349 489L395 552H302Q128 552 119 546Q113 543 108 522T98 479L95 458V455H55V458Z"></path><path stroke-width="1" id="MJMATHI-6D" d="M21 287Q22 293 24 303T36 341T56 388T88 425T132 442T175 435T205 417T221 395T229 376L231 369Q231 367 232 367L243 378Q303 442 384 442Q401 442 415 440T441 433T460 423T475 411T485 398T493 385T497 373T500 364T502 357L510 367Q573 442 659 442Q713 442 746 415T780 336Q780 285 742 178T704 50Q705 36 709 31T724 26Q752 26 776 56T815 138Q818 149 821 151T837 153Q857 153 857 145Q857 144 853 130Q845 101 831 73T785 17T716 -10Q669 -10 648 17T627 73Q627 92 663 193T700 345Q700 404 656 404H651Q565 404 506 303L499 291L466 157Q433 26 428 16Q415 -11 385 -11Q372 -11 364 -4T353 8T350 18Q350 29 384 161L420 307Q423 322 423 345Q423 404 379 404H374Q288 404 229 303L222 291L189 157Q156 26 151 16Q138 -11 108 -11Q95 -11 87 -5T76 7T74 17Q74 30 112 181Q151 335 151 342Q154 357 154 369Q154 405 129 405Q107 405 92 377T69 316T57 280Q55 278 41 278H27Q21 284 21 287Z"></path><path stroke-width="1" id="MJMATHI-75" d="M21 287Q21 295 30 318T55 370T99 420T158 442Q204 442 227 417T250 358Q250 340 216 246T182 105Q182 62 196 45T238 27T291 44T328 78L339 95Q341 99 377 247Q407 367 413 387T427 416Q444 431 463 431Q480 431 488 421T496 402L420 84Q419 79 419 68Q419 43 426 35T447 26Q469 29 482 57T512 145Q514 153 532 153Q551 153 551 144Q550 139 549 130T540 98T523 55T498 17T462 -8Q454 -10 438 -10Q372 -10 347 46Q345 45 336 36T318 21T296 6T267 -6T233 -11Q189 -11 155 7Q103 38 103 113Q103 170 138 262T173 379Q173 380 173 381Q173 390 173 393T169 400T158 404H154Q131 404 112 385T82 344T65 302T57 280Q55 278 41 278H27Q21 284 21 287Z"></path><path stroke-width="1" id="MJMATHI-74" d="M26 385Q19 392 19 395Q19 399 22 411T27 425Q29 430 36 430T87 431H140L159 511Q162 522 166 540T173 566T179 586T187 603T197 615T211 624T229 626Q247 625 254 615T261 596Q261 589 252 549T232 470L222 433Q222 431 272 431H323Q330 424 330 420Q330 398 317 385H210L174 240Q135 80 135 68Q135 26 162 26Q197 26 230 60T283 144Q285 150 288 151T303 153H307Q322 153 322 145Q322 142 319 133Q314 117 301 95T267 48T216 6T155 -11Q125 -11 98 4T59 56Q57 64 57 83V101L92 241Q127 382 128 383Q128 385 77 385H26Z"></path><path stroke-width="1" id="MJMATHI-65" d="M39 168Q39 225 58 272T107 350T174 402T244 433T307 442H310Q355 442 388 420T421 355Q421 265 310 237Q261 224 176 223Q139 223 138 221Q138 219 132 186T125 128Q125 81 146 54T209 26T302 45T394 111Q403 121 406 121Q410 121 419 112T429 98T420 82T390 55T344 24T281 -1T205 -11Q126 -11 83 42T39 168ZM373 353Q367 405 305 405Q272 405 244 391T199 357T170 316T154 280T149 261Q149 260 169 260Q282 260 327 284T373 353Z"></path><path stroke-width="1" id="MJMATHI-78" d="M52 289Q59 331 106 386T222 442Q257 442 286 424T329 379Q371 442 430 442Q467 442 494 420T522 361Q522 332 508 314T481 292T458 288Q439 288 427 299T415 328Q415 374 465 391Q454 404 425 404Q412 404 406 402Q368 386 350 336Q290 115 290 78Q290 50 306 38T341 26Q378 26 414 59T463 140Q466 150 469 151T485 153H489Q504 153 504 145Q504 144 502 134Q486 77 440 33T333 -11Q263 -11 227 52Q186 -10 133 -10H127Q78 -10 57 16T35 71Q35 103 54 123T99 143Q142 143 142 101Q142 81 130 66T107 46T94 41L91 40Q91 39 97 36T113 29T132 26Q168 26 194 71Q203 87 217 139T245 247T261 313Q266 340 266 352Q266 380 251 392T217 404Q177 404 142 372T93 290Q91 281 88 280T72 278H58Q52 284 52 289Z"></path></defs></svg>