2025-11-21T00:28:16.024918

UltraScatter: Ray-Based Simulation of Ultrasound Scattering

Duelmer, Azampour, Navab
Traditional ultrasound simulation methods solve wave equations numerically, achieving high accuracy but at substantial computational cost. Faster alternatives based on convolution with precomputed impulse responses remain relatively slow, often requiring several minutes to generate a full B-mode image. We introduce UltraScatter, a probabilistic ray tracing framework that models ultrasound scattering efficiently and realistically. Tissue is represented as a volumetric field of scattering probability and scattering amplitude, and ray interactions are simulated via free-flight delta tracking. Scattered rays are traced to the transducer, with phase information incorporated through a linear time-of-flight model. Integrated with plane-wave imaging and beamforming, our parallelized ray tracing architecture produces B-mode images within seconds. Validation with phantom data shows realistic speckle and inclusion patterns, positioning UltraScatter as a scalable alternative to wave-based methods.
academic

UltraScatter: অতিশব্দ বিক্ষেপণের রশ্মি-ভিত্তিক সিমুলেশন

মৌলিক তথ্য

  • পেপার আইডি: 2510.10612
  • শিরোনাম: UltraScatter: Ray-Based Simulation of Ultrasound Scattering
  • লেখক: Felix Duelmer, Mohammad Farid Azampour, Nassir Navab
  • শ্রেণীবিভাগ: physics.med-ph cs.CV
  • প্রকাশনার সময়: ২০২৫ সালের অক্টোবর ১২ (arXiv প্রি-প্রিন্ট)
  • পেপার লিংক: https://arxiv.org/abs/2510.10612

সারসংক্ষেপ

ঐতিহ্যবাহী অতিশব্দ সিমুলেশন পদ্ধতি তরঙ্গ সমীকরণ সংখ্যাগতভাবে সমাধান করে উচ্চ নির্ভুলতা অর্জন করে, কিন্তু গণনার খরচ অত্যন্ত বেশি। পূর্বনির্ধারিত আবেগ প্রতিক্রিয়া কনভোলিউশনের উপর ভিত্তি করে দ্রুত বিকল্প পদ্ধতিগুলি এখনও তুলনামূলকভাবে ধীর, সাধারণত সম্পূর্ণ B-মোড ইমেজ তৈরি করতে কয়েক মিনিট প্রয়োজন। এই পত্রটি UltraScatter উপস্থাপন করে, একটি সম্ভাব্যতামূলক রশ্মি ট্রেসিং ফ্রেমওয়ার্ক যা দক্ষতার সাথে এবং বাস্তবসম্মতভাবে অতিশব্দ বিক্ষেপণ মডেল করতে পারে। টিস্যুকে বিক্ষেপণ সম্ভাবনা এবং বিক্ষেপণ প্রশস্ততার ভলিউম ক্ষেত্র হিসাবে প্রতিনিধিত্ব করা হয়, রশ্মি মিথস্ক্রিয়া মুক্ত-ফ্লাইট ডেল্টা ট্রেসিং এর মাধ্যমে সিমুলেট করা হয়। বিক্ষিপ্ত রশ্মিগুলি সেন্সরে ট্রেস করা হয়, পর্যায় তথ্য রৈখিক ফ্লাইট সময় মডেলের মাধ্যমে একীভূত করা হয়। সমতল তরঙ্গ ইমেজিং এবং বিম গঠনের সাথে মিলিত, সমান্তরালকৃত রশ্মি ট্রেসিং আর্কিটেকচার কয়েক সেকেন্ডের মধ্যে B-মোড ইমেজ তৈরি করে। ফ্যান্টম ডেটা যাচাইকরণ বাস্তবসম্মত দাগ এবং অন্তর্ভুক্তি প্যাটার্ন প্রদর্শন করে, UltraScatter কে তরঙ্গ সমীকরণ পদ্ধতির একটি স্কেলেবল বিকল্প হিসাবে অবস্থান করে।

গবেষণা পটভূমি এবং প্রেরণা

সমস্যা সংজ্ঞা

অতিশব্দ সিমুলেশন পুনর্নির্মাণ অ্যালগরিদম উন্নত করতে, সেন্সর ডিজাইন অপ্টিমাইজ করতে এবং নিয়ন্ত্রিত শর্তে মেশিন লার্নিং মডেল প্রশিক্ষণ দিতে গুরুত্বপূর্ণ। তবে বিদ্যমান পদ্ধতিগুলি গণনার দক্ষতা এবং শারীরিক বাস্তবতার মধ্যে ট্রেড-অফের সম্মুখীন হয়।

বিদ্যমান পদ্ধতির সীমাবদ্ধতা

  1. তরঙ্গ সমীকরণ সমাধানকারী (যেমন FDTD, k-space সিউডোস্পেক্ট্রাল পদ্ধতি): উচ্চ নির্ভুলতা কিন্তু বিশাল গণনার খরচ, স্কেলিং কঠিন
  2. Field II এর মতো দ্রুত পদ্ধতি: স্থানিক আবেগ প্রতিক্রিয়া পূর্বনির্ধারণ করে কনভোলিউশনের মাধ্যমে, একক B-মোড ইমেজ তৈরি করতে এখনও কয়েক মিনিট প্রয়োজন
  3. SIMUS: জ্যামিতিক রাউন্ড-ট্রিপ দূরত্বের উপর ভিত্তি করে বিলম্ব সমষ্টি অপারেশন, সীমিত দক্ষতা

গবেষণা প্রেরণা

শারীরিক বাস্তবতা এবং দ্রুত গণনা একত্রিত করে এমন একটি পদ্ধতির প্রয়োজন, যা রিয়েল-টাইম বা ইন্টারেক্টিভ অতিশব্দ সিমুলেশন সম্ভব করে। কম্পিউটার গ্রাফিক্সে শারীরিক ভিত্তিক রেন্ডারিং (PBR) প্রযুক্তি এর জন্য অনুপ্রেরণা প্রদান করে, বিশেষত অত্যন্ত অপ্টিমাইজড CUDA ত্বরিত রশ্মি ট্রেসিং অ্যালগরিদম।

মূল অবদান

  1. মডুলার উচ্চ-কর্মক্ষমতা ফ্রেমওয়ার্ক: অংশগ্রহণকারী মাধ্যমে ক্ষয়, শোষণ এবং একাধিক বিক্ষেপণ মডেলিং
  2. ট্রান্সমিটার স্যাম্পলিং কৌশল: প্রতিটি দৃশ্য মিথস্ক্রিয়া সমস্ত সেন্সর উপাদানের সাথে সংযুক্ত করে
  3. সম্পূর্ণ ট্রান্সমিট-রিসিভ বিম গঠন শৃঙ্খল: সরাসরি সিমুলেটেড ইকো B-মোড ইমেজে রূপান্তরিত করে
  4. উল্লেখযোগ্য কর্মক্ষমতা বৃদ্ধি: ঐতিহ্যবাহী ফ্রিকোয়েন্সি-ডোমেইন কোড বাস্তবায়নের তুলনায় প্রায় ৭০ গুণ গতি বৃদ্ধি

পদ্ধতি বিস্তারিত

কাজের সংজ্ঞা

লেবেল ম্যাপ থেকে শুরু করে, প্রতিটি টিস্যু বিভাগে বিক্ষেপণ বৈশিষ্ট্য নির্ধারণ করে, চাপ তরঙ্গ নির্গমন, বিক্ষেপণ এবং ক্ষয় মডেল করতে মন্টে কার্লো রশ্মি ট্রেসিং স্কিম ব্যবহার করে, সেন্সরে ফিরে আসা ইকো নির্দিষ্ট উপাদানের RF বাফারে লিখে, তারপর প্রচলিত ডিজিটাল বিম ফর্মার দ্বারা চূড়ান্ত B-মোড ইমেজ তৈরি করে।

মডেল আর্কিটেকচার

মূল গাণিতিক মডেল

সেন্সর উপাদান e তে সময় t তে পৌঁছানো চাপ সংকেত P সংজ্ঞায়িত করা হয়:

P(e,t)=ΩAPi(x,t,ωi)fd(ωi)dωdaP(e, t) = \int_\Omega \int_A P_i(x, t, \omega_i) f_d(\omega_i) d\omega da

যেখানে Pi(x,t,ωi)P_i(x, t, \omega_i) দৃশ্য অবস্থান x, দিক ωi\omega_i থেকে আসা আপতিত চাপ, fdf_d দিকত্ব ফাংশন।

বিক্ষিপ্ত চাপ মডেল করা হয়:

Pscattered(x,t,ωo)=Ωa(x)p(ωi,ωo)Pi(x,t,ωi)dωiP_{scattered}(x, t, \omega_o) = \int_\Omega a(x) p(\omega_i, \omega_o)P_i(x, t, \omega_i) d\omega_i

যেখানে a(x)a(x) বিক্ষেপণ প্রশস্ততা প্রতিনিধিত্ব করে, p(ωi,ωo)p(\omega_i, \omega_o) কোণীয় পুনর্বিতরণ নিয়ন্ত্রণ করে এমন পর্যায় ফাংশন।

রশ্মি ট্রেসিং অ্যালগরিদম

১. মুক্ত-ফ্লাইট দূরত্ব স্যাম্পলিং: s=smin1μln(1ξ)s = s_{min} - \frac{1}{\mu} \ln(1-\xi)

যেখানে μ\mu সম্পূর্ণ মাধ্যমে সর্বোচ্চ বিলুপ্তি মূল্য, ξU(0,1)\xi \sim U(0,1) র্যান্ডম সংখ্যা।

২. মিথস্ক্রিয়া সিদ্ধান্ত: ξ<σt(x)μপ্রকৃত মিথস্ক্রিয়া, অন্যথায় শূন্য মিথস্ক্রিয়া\xi' < \frac{\sigma_t(x)}{\mu} \Rightarrow \text{প্রকৃত মিথস্ক্রিয়া, অন্যথায় শূন্য মিথস্ক্রিয়া}

३. ট্রান্সমিটার স্যাম্পলিং কৌশল: প্রতিটি বিক্ষেপণ ঘটনায় সমস্ত সেন্সর উপাদানের জন্য স্যাম্পলিং, গৌণ রশ্মির একটি সেট নির্গত করে, প্রতিটি একটি সেন্সর উপাদানের দিকে।

४. দূরক্ষেত্র আনুমানিক প্রক্রিয়াকরণ: প্রতিটি সেন্সর উপাদানকে ν\nu টি অভিন্ন উপ-উপাদানে বিভক্ত করে: ν=2bλmin\nu = \left\lceil \frac{2b}{\lambda_{min}} \right\rceil

যেখানে bb একক সেন্সর উপাদানের অ্যাজিমুথাল প্রস্থ, λmin\lambda_{min} ট্রান্সমিশন পালসে ন্যূনতম তরঙ্গদৈর্ঘ্য।

প্রযুক্তিগত উদ্ভাবন পয়েন্ট

  1. Woodcock ট্রেসিং: শূন্য মিথস্ক্রিয়া সহ মুক্ত-ফ্লাইট ডেল্টা ট্রেসিং গ্রহণ, গণনা দক্ষতা উন্নত করে
  2. সম্পূর্ণ উপাদান ট্রান্সমিটার স্যাম্পলিং: নির্বাচনী পদ্ধতির তুলনায়, স্পষ্টতর ইমেজ এবং তীক্ষ্ণতর কাঠামো সীমানা উৎপাদন করে
  3. GPU সমান্তরালকরণ: Mitsuba 3 ফ্রেমওয়ার্কের উপর ভিত্তি করে, NVIDIA OptiX ব্যবহার করে দক্ষ সমান্তরাল গণনা
  4. পর্যায় তথ্য একীকরণ: রৈখিক ফ্লাইট সময় মডেলের মাধ্যমে পর্যায় সামঞ্জস্য বজায় রাখে

পরীক্ষামূলক সেটআপ

ডেটাসেট

CIRS সার্বজনীন ফ্যান্টম (মডেল 054GS) যাচাইকরণের জন্য ব্যবহার করা হয়, নির্মাতার স্পেসিফিকেশন অনুযায়ী এর অভ্যন্তরীণ কাঠামো পুনর্নির্মাণ করে, এবং অনুপস্থিত পরামিতি সামঞ্জস্য করে যতক্ষণ না সিমুলেশন রেফারেন্স ক্লিনিকাল B-মোডের সাথে মেলে।

তুলনা পদ্ধতি

  1. ক্লিনিকাল স্ক্যানার: Siemens Acuson Juniper 12L3 রৈখিক অ্যারে সহ (192 উপাদান, 2.9-11.5 MHz)
  2. SIMUS: ফ্রিকোয়েন্সি-ডোমেইন সময়-সুসংগত সিমুলেটর, দূরক্ষেত্র এবং প্যারাক্সিয়াল অনুমানের উপর ভিত্তি করে

বাস্তবায়ন বিবরণ

  • কেন্দ্র ফ্রিকোয়েন্সি: 6.2 MHz
  • ইমেজিং গভীরতা: 10 cm
  • গতিশীল পরিসীমা: 60 dB
  • স্যাম্পলিং হার: 25 MHz
  • F সংখ্যা: 1.0
  • প্রতি উপ-উপাদান 100,000 রশ্মি নির্গমন
  • হার্ডওয়্যার: Intel Core i7-12700 CPU + NVIDIA RTX 4070 Ti GPU

পরীক্ষামূলক ফলাফল

প্রধান ফলাফল

গণনা কর্মক্ষমতা তুলনা:

  • SIMUS: 634 ± 3 সেকেন্ড
  • UltraScatter: 9.3 ± 0.8 সেকেন্ড
  • গতি বৃদ্ধি: প্রায় 70 গুণ

ইমেজ গুণমান মূল্যায়ন:

  1. সামগ্রিক শারীরিক কাঠামো: তিনটি পদ্ধতি (ক্লিনিকাল, UltraScatter, SIMUS) অনুরূপ সামগ্রিক শারীরিক কাঠামো ভাগ করে
  2. পটভূমি দাগ পরিসংখ্যান: UltraScatter প্রকৃত স্ক্যানের অনুরূপ দাগ পরিসংখ্যান উৎপাদন করে
  3. দূরান্ত ছায়া: প্রকৃত স্ক্যানের ছায়া প্রভাব সঠিকভাবে প্রতিফলিত করে

কেস বিশ্লেষণ

UltraScatter বৈশিষ্ট্য:

  • প্রোব সরাসরি নীচের অঞ্চলে অত্যধিক আলোকসজ্জা বিদ্যমান, যা নিকটক্ষেত্র লাভ অবশিষ্ট অতিমূল্যায়নের ফলাফল
  • আরও তীক্ষ্ণ অন্তর্ভুক্তি সীমানা উৎপাদন করে, সম্পূর্ণ উপাদান ট্রান্সমিটার স্যাম্পলিং কৌশলের জন্য দায়ী
  • গভীর সিলিন্ডারের জন্য SIMUS এর চেয়ে স্পষ্টতর রেজোলিউশন

SIMUS বৈশিষ্ট্য:

  • ক্লিনিকাল সিস্টেমের অক্ষীয় তীব্রতা ক্ষয় আরও বিশ্বস্তভাবে পুনরুৎপাদন করে
  • অন্তর্ভুক্তি প্রান্ত সংজ্ঞা সামান্য অস্পষ্ট

সম্পর্কিত কাজ

ঐতিহ্যবাহী অতিশব্দ সিমুলেশন পদ্ধতি

  1. সম্পূর্ণ-তরঙ্গ সমাধানকারী: FDTD, k-space সিউডোস্পেক্ট্রাল পদ্ধতি - নির্ভুল কিন্তু গণনা ব্যয়বহুল
  2. কনভোলিউশন পদ্ধতি: Field II - পূর্বনির্ধারিত আবেগ প্রতিক্রিয়ার উপর ভিত্তি করে
  3. জ্যামিতিক পদ্ধতি: SIMUS - বিলম্ব সমষ্টি অপারেশনের উপর ভিত্তি করে

কনভোলিউশন রশ্মি ট্রেসিং (CRT)

সম্প্রতি উদ্ভূত দ্রুত বিকল্প পদ্ধতি, বৃহৎ রশ্মি সেট ট্রেসিং করে শব্দ প্রচার আনুমানিক করে, ম্যাক্রোস্কোপিক সীমানায় প্রতিফলন এবং প্রতিসরণ, তারপর বিভাজনযোগ্য পয়েন্ট ছড়িয়ে ফাংশন কনভোলিউশনের মাধ্যমে সূক্ষ্ম-স্কেল বিক্ষেপণ মডেলিং।

শারীরিক রেন্ডারিং প্রযুক্তি

কম্পিউটার গ্রাফিক্সে PBR প্রযুক্তি, রশ্মি মার্চিং, ফটন ম্যাপিং এবং ডেল্টা ট্রেসিং সহ, বিষমজাত মাধ্যমে আলো প্রচারের জন্য অত্যন্ত অপ্টিমাইজড বাস্তবায়ন প্রদান করে।

উপসংহার এবং আলোচনা

প্রধান উপসংহার

UltraScatter সফলভাবে রশ্মি ট্রেসিং অতিশব্দ সিমুলেশন বাস্তবায়ন করে, ফ্রিকোয়েন্সি-ডোমেইন সমাধানকারীর জায়গায় মন্টে কার্লো পথ স্যাম্পলিং ব্যবহার করে, একই হার্ডওয়্যারে সাধারণ CIRS ফ্যান্টম ভিউ রেন্ডার করতে প্রায় 9 সেকেন্ড প্রয়োজন, যখন SIMUS 10 মিনিটের বেশি প্রয়োজন, উল্লেখযোগ্য কর্মক্ষমতা বৃদ্ধি অর্জন করে।

সীমাবদ্ধতা

  1. দাগ পরিসংখ্যান: র্যান্ডম স্যাম্পলিং দ্বারা চালিত, দুটি স্বাধীন রান বিভিন্ন প্যাটার্ন উৎপাদন করে, প্রোব পার্শ্বীয় গতির সময় টেক্সচার ডিকোহেরেন্স
  2. রশ্মি কনফিগারেশন: বর্তমান প্রোটোটাইপ একক অক্ষীয় লাইন থেকে ঘনক গ্রিডে রশ্মি নির্গত করে, আরও বাস্তবসম্মত কনফিগারেশন সেন্সর জুড়ে সম্পূর্ণ elevational অ্যাপারচার প্রয়োজন
  3. সময় সামঞ্জস্য: ছোট প্রোব গতিতে দাগ স্থিতিশীলতা বজায় রাখতে সম্পর্কিত র্যান্ডম ক্রম বা ভাগ করা বীজ ম্যাপ প্রবর্তন প্রয়োজন

ভবিষ্যত দিকনির্দেশনা

  1. Elevational ফোকাসিং: গতিশীল গ্রহণ বিলম্ব বা অতিরিক্ত রশ্মি প্রজেকশনের মাধ্যমে বাস্তবায়ন
  2. উন্নত শারীরিক মডেলিং: প্রতিটি ভক্সেলে স্বাধীন ঘনত্ব, শব্দ গতি এবং ক্ষয় নির্ধারণ
  3. অরৈখিক ঘটনা: ম্যাক্রোস্কোপিক প্রতিসরণ, পর্যায় বিকৃতি এমনকি সুরেলা প্রজন্ম মডেলিং
  4. গতিশীল দৃশ্য সমর্থন: দাগ সামঞ্জস্য এবং elevational রেজোলিউশন উন্নত করা

গভীর মূল্যায়ন

সুবিধা

  1. উল্লেখযোগ্য কর্মক্ষমতা বৃদ্ধি: 70 গুণ গতি বৃদ্ধি রিয়েল-টাইম সিমুলেশন সম্ভব করে
  2. শারীরিক বাস্তবতা: ক্লিনিকাল স্ক্যানের অনুরূপ ইমেজ গুণমান এবং দাগ প্যাটার্ন উৎপাদন করে
  3. মডুলার আর্কিটেকচার: Mitsuba 3 ভিত্তিক ডিজাইন সম্প্রসারণ এবং উন্নতি সহজ করে
  4. GPU সমান্তরালকরণ: আধুনিক হার্ডওয়্যারের সমান্তরাল গণনা ক্ষমতা সম্পূর্ণভাবে ব্যবহার করে
  5. সম্পূর্ণ পাইপলাইন: রশ্মি ট্রেসিং থেকে B-মোড ইমেজ প্রজন্মের সম্পূর্ণ সমাধান

অপূর্ণতা

  1. নিকটক্ষেত্র প্রভাব: নিকটক্ষেত্র লাভ অতিমূল্যায়ন সমস্যা বিদ্যমান
  2. দাগ সামঞ্জস্য: সময় সামঞ্জস্যের অভাব, গতিশীল দৃশ্যের জন্য অনুপযুক্ত
  3. মাত্রা সীমাবদ্ধতা: বর্তমানে শুধুমাত্র 2D ইমেজিং সমর্থন করে, সম্পূর্ণ 3D elevational প্রক্রিয়াকরণের অভাব
  4. যাচাইকরণ পরিসীমা: একক ফ্যান্টমে শুধুমাত্র যাচাইকৃত, বিস্তৃত যাচাইকরণ ডেটার অভাব

প্রভাব

একাডেমিক অবদান:

  • কম্পিউটার গ্রাফিক্সের রশ্মি ট্রেসিং প্রযুক্তি অতিশব্দ সিমুলেশনে সফলভাবে প্রয়োগ করা প্রথম
  • অতিশব্দ সিমুলেশন ক্ষেত্রে নতুন প্রযুক্তিগত পথ প্রদান করে

ব্যবহারিক মূল্য:

  • রিয়েল-টাইম অতিশব্দ সিমুলেশন সম্ভব করে, চিকিৎসা শিক্ষা এবং প্রশিক্ষণে সহায়তা করে
  • মেশিন লার্নিং মডেলের জন্য দ্রুত ডেটা প্রজন্ম সরঞ্জাম প্রদান করে
  • সেন্সর ডিজাইন অপ্টিমাইজেশনের দ্রুত পুনরাবৃত্তি সমর্থন করে

পুনরুৎপাদনযোগ্যতা:

  • খোলা উৎস Mitsuba 3 ফ্রেমওয়ার্কের উপর ভিত্তি করে
  • বিস্তারিত বাস্তবায়ন বিবরণ এবং পরামিতি সেটিংস প্রদান করে
  • কোড এবং ডেটার প্রাপ্যতা আরও নিশ্চিতকরণ প্রয়োজন

প্রযোজ্য পরিস্থিতি

  1. চিকিৎসা শিক্ষা: রিয়েল-টাইম অতিশব্দ সিমুলেশন প্রশিক্ষণ সিস্টেম
  2. অ্যালগরিদম উন্নয়ন: গভীর শিক্ষা মডেলের জন্য প্রশিক্ষণ ডেটা দ্রুত প্রজন্ম
  3. ডিভাইস ডিজাইন: সেন্সর পরামিতি অপ্টিমাইজেশনের দ্রুত প্রোটোটাইপ যাচাইকরণ
  4. গবেষণা প্রয়োগ: প্রচুর সিমুলেশন ডেটা প্রয়োজনীয় বৈজ্ঞানিক গবেষণা প্রকল্প

সংদর্ভ

পত্রটি 19টি গুরুত্বপূর্ণ সংদর্ভ উদ্ধৃত করে, যা অন্তর্ভুক্ত করে:

  • ঐতিহ্যবাহী অতিশব্দ সিমুলেশন পদ্ধতি (FDTD, k-Wave, Field II, SIMUS)
  • রশ্মি ট্রেসিং অতিশব্দ সিমুলেশন সম্পর্কিত কাজ
  • কম্পিউটার গ্রাফিক্স রেন্ডারিং প্রযুক্তি (Mitsuba 3, OptiX)
  • মন্টে কার্লো পদ্ধতি এবং ডেল্টা ট্রেসিং অ্যালগরিদম

সারসংক্ষেপ: UltraScatter অতিশব্দ সিমুলেশন ক্ষেত্রে একটি গুরুত্বপূর্ণ প্রযুক্তিগত অগ্রগতি প্রতিনিধিত্ব করে, কম্পিউটার গ্রাফিক্সের উন্নত রশ্মি ট্রেসিং প্রযুক্তি চিকিৎসা অতিশব্দ সিমুলেশনে সফলভাবে প্রবর্তন করে, গণনা দক্ষতা এবং শারীরিক বাস্তবতার মধ্যে ভাল ভারসাম্য অর্জন করে। কিছু সীমাবদ্ধতা থাকলেও, এর উল্লেখযোগ্য কর্মক্ষমতা বৃদ্ধি এবং মডুলার ডিজাইন ভবিষ্যত গবেষণার জন্য দৃঢ় ভিত্তি স্থাপন করে।