GPU Monte Carlo scatter calculations for Cone Beam Computed Tomography

• Main Author: Adler, Jonas
• Format: Others
• Language: English
• Published: KTH, Numerisk analys, NA 2014
• Online Access: http://urn.kb.se/resolve?urn=urn:nbn:se:kth:diva-140221

A GPU Monte Carlo code for x-ray photon transport has been implemented and extensively tested. The code is intended for scatter compensation of cone beam computed tomography images. The code was tested to agree with other well known codes within 5% for a set of simple scenarios. The scatter compensation was also tested using an artificial head phantom. The errors in the reconstructed Hounsfield values were reduced by approximately 70%. Several variance reduction methods have been tested, although most were found infeasible on GPUs. The code is nonetheless fast, and can simulate approximately 3 ·109 photons per minute on a NVIDIA Quadro 4000 graphics card. With the use of appropriate filtering methods, the code can be used to calculate patient specific scatter distributions for a full CBCT scan in approximately one minute, allowing scatter reduction in clinical applications. === En GPU Monte Carlo kod för transport av röntgenfotoner har implementerats och utförligt testats. Koden är avsed för spridningskorrektion av CBCT-bilder. Koden har testats mot PENELOPE och resultaten överenstämmer inom 5% för ett antal enklare geometrier. Koden testades också i en verklig uppställning med ett artificiellt huvud. De resulterande felen i de beräknade Hounsfieldvärdena minbskade med ca 70%. Ett antal variansreduktionstekniker har också testats, men de flesta gav ingen förbättring på GPU. Koden är trots detta avsevärt snabb och kan simulera ca 3 · 109 photoner per minut med ett Quadro 4000 grafik-kort. Med hjälp av väl valda filtreringsmetoder kan koden användas för att beräkna patientspecifika spridningsfördelningar för ett fullständigt CBCT-scan på under en minut. Detta är tillräkligt för spridningskorrektion i kliniska tillämpningar.