Görüntü Rekonstrüksiyon ve İşleme

8.1 Filtrelenmiş Geri Projeksiyon (FBP)

8.1.1 Temel Kavramlar: Projeksiyon ve Sinogram1 soru

Nükleer tıp görüntüleme sistemlerinde (SPECT ve PET), hastadan farklı açılardan toplanan projeksiyon verileri kullanılarak kesitsel (tomografik) görüntüler oluşturulur. Bu süreç görüntü rekonstrüksiyonu olarak adlandırılır.

Projeksiyon verisi:

Sinogram:

Şekil 8.1: Projeksiyon ve sinogram: merkezdeki nokta kaynak düz dikey çizgi, merkez dışı kaynak sinüzoidal eğri oluşturur.

8.1.2 Basit Geri Projeksiyon (Simple Back Projection)

Basit geri projeksiyon, rekonstrüksiyonun en temel ve sezgisel yöntemidir:

Prensip:

1. Her açıdan elde edilen projeksiyon verisi, o açı boyunca görüntü matrisine geri yansıtılır (back project)

2. Tüm açılardan geri yansıtılan veriler toplanır 3. Yeterli sayıda açıdan veri toplandığında, aktivite dağılımının bir yaklaşımı elde edilir

Temel sorun – Yıldız artefaktı (Star artifact / 1/r blurring):

Şekil 8.2: Basit geri projeksiyonda açı sayısı arttıkça yıldız artefaktı 1/r bulanıklığına dönüşür; rampa filtreli FBP bu bulanıklığı düzeltir.

8.1.3 Filtrelenmiş Geri Projeksiyon (FBP)1 soru

FBP, basit geri projeksiyondaki 1/r bulanıklığını ortadan kaldırmak için projeksiyon verilerine geri projeksiyon öncesinde bir rampa filtresi (ramp filter) uygular.

Matematiksel temel – Fourier Kesit Teoremi (Central Slice Theorem):

Rampa filtresi (Ramp filter):

FBP algoritması adımları:

1. Her açıdan projeksiyon verisi toplanır 2. Her projeksiyon profiline 1D Fourier dönüşümü uygulanır 3. Frekans uzayında rampa filtresi (+ pencere fonksiyonu) uygulanır 4. Ters Fourier dönüşümü ile filtrelenmiş projeksiyon elde edilir 5. Filtrelenmiş projeksiyonlar geri yansıtılır 6. Tüm açılardan geri yansıtılan veriler toplanarak kesitsel görüntü oluşturulur

Şekil 8.3: Fourier kesit teoremi ve FBP algoritmasının adımları (şematik).

FBP'nin avantajları:

FBP'nin dezavantajları:

8.2 İteratif Rekonstrüksiyon (OSEM, MLEM)

8.2.1 İteratif Rekonstrüksiyonun Temel Prensibi2 soru

İteratif rekonstrüksiyon, görüntüyü adım adım (iterasyon) iyileştirerek gerçek aktivite dağılımına yaklaşan bir yöntemdir. FBP'den farklı olarak analitik bir çözüm yerine ardışık tahmin ve düzeltme döngüsü kullanır.

Genel algoritma:

1. Başlangıç tahmini: Homojen (uniform) bir görüntü ile başlanır 2. İleri projeksiyon (forward projection): Tahmin edilen görüntüden yapay projeksiyon verileri hesaplanır

3. Karşılaştırma: Hesaplanan projeksiyon verileri, gerçek (ölçülen) projeksiyon verileriyle karşılaştırılır

4. Güncelleme: Fark veya oran kullanılarak görüntü tahmini güncellenir 5. Tekrar: Adımlar 2-4, belirli bir yakınsama kriteri sağlanana veya belirlenen iterasyon sayısına ulaşılana kadar tekrarlanır

İteratif rekonstrüksiyonun avantajları (FBP'ye göre):

İteratif rekonstrüksiyonun dezavantajları:

8.2.2 MLEM (Maximum Likelihood Expectation Maximization)1 soru

MLEM, istatistiksel olarak en doğru iteratif rekonstrüksiyon algoritmasıdır.

Prensip:

MLEM güncelleme formülü:

fⱼⁿ⁺¹ = [fⱼⁿ / Σᵢ aᵢⱼ] · Σᵢ [aᵢⱼ · gᵢ / Σₖ aᵢₖ fₖⁿ]

Burada:

MLEM'in özellikleri:

8.2.3 OSEM (Ordered Subsets Expectation Maximization)4 soru

OSEM, MLEM'in hızlandırılmış versiyonudur ve günümüzde hem SPECT hem PET'te standart rekonstrüksiyon algoritmasıdır.

Prensip:

Hızlanma faktörü:

Şekil 8.4: İteratif rekonstrüksiyon döngüsü ve OSEM'de projeksiyonların eşit aralıklı alt kümelere bölünmesi (şematik).

Efektif iterasyon sayısı:

Efektif iterasyon = İterasyon sayısı × Subset sayısı

OSEM parametre seçimi:

8.2.4 PSF Modelleme (Resolution Recovery)2 soru

Modern iteratif rekonstrüksiyon algoritmalarında, dedektörün nokta yayılma fonksiyonu (Point Spread Function – PSF) sistem matrisine dahil edilebilir.

PSF modellemenin etkileri:

Şekil 8.5: İterasyon sayısı arttıkça kontrast ve gürültü birlikte artar; PSF modelleme kenarları keskinleştirir ancak Gibbs aşımı oluşturabilir (simülasyon).

Ticari isimler:

8.2.5 FBP vs İteratif Rekonstrüksiyon Karşılaştırması

8.3 Atenüasyon Düzeltme (CT-bazlı, Chang)

8.3.1 Atenüasyon Sorunu

Hastanın vücudundan yayılan gama fotonları, dedektöre ulaşmadan önce doku tarafından absorbe edilir veya saçılır. Bu süreç atenüasyon olarak adlandırılır ve düzeltilmezse ciddi kantitatif hatalara ve artefaktlara neden olur.

Atenüasyonun etkileri:

Atenüasyon katsayısı:

Şekil 8.6: SPECT'te atenüasyon kaynak derinliğine bağlıdır; PET'te foton çiftinin atenüasyonu yalnızca toplam kalınlığa (D) bağlıdır.

8.3.2 SPECT Atenüasyon Düzeltme Yöntemleri4 soru

1. Chang Yöntemi (Birinci Derece Chang Düzeltmesi)

En basit ve en yaygın kullanılan SPECT atenüasyon düzeltme yöntemidir.

Prensip:

Chang düzeltme formülü:

C(x,y) = 1 / [(1/N) · Σ_i=1^N e^{−µ·d_i(x,y)}]

Burada:

Uygulama:

Chang yönteminin sınırlılıkları:

2. CT-Bazlı Atenüasyon Düzeltme (CTAC) – SPECT/CT

SPECT/CT sistemlerinde CT verileri kullanılarak voksel bazlı atenüasyon haritası oluşturulur.

Prensip:

CTAC'ın Chang'e göre avantajları:

CTAC'ın potansiyel sorunları:

8.3.3 PET Atenüasyon Düzeltme3 soru

PET'te atenüasyonun özelliği:

PET'te iki foton (511 keV) birlikte atenüasyona uğrar. Toplam atenüasyon, foton çiftinin hastanın tüm kalınlığı boyunca geçtiği yol uzunluğuna bağlıdır ve kaynak konumundan bağımsızdır.

ACF = e^+µ·D (atenüasyon faktörü = e^−µ·D)

Burada D, foton çiftinin geçtiği toplam yol uzunluğudur (hasta kalınlığı). Bu özellik PET'e özgüdür ve atenüasyon düzeltmeyi basitleştirir.

PET atenüasyon düzeltme yöntemleri:

1. Transmisyon kaynağı (eski yöntem):

2. CT-bazlı atenüasyon düzeltme (CTAC) – Standart yöntem:

3. MR-bazlı atenüasyon düzeltme (MRAC) – PET/MR:

8.4 Saçılma Düzeltme

8.4.1 Saçılma Sorunu

Compton saçılması, fotonun yön ve enerji değiştirerek dedektöre ulaşmasına neden olur. Saçılmış fotonlar yanlış konumda kaydedilir → görüntüde kontrast kaybı ve kantitatif hata oluşur.

Saçılma fraksiyonu:

8.4.2 SPECT Saçılma Düzeltme Yöntemleri1 soru

1. Enerji Penceresi Yöntemleri:

a) Çift enerji penceresi (Dual Energy Window – DEW / Jaszczak yöntemi):

b) Üçlü enerji penceresi (Triple Energy Window – TEW):

Şekil 8.7: Tc-99m spektrumunda çift (DEW) ve üçlü (TEW) enerji penceresi ile saçılma tahmini (şematik).

2. Model Bazlı Saçılma Düzeltme:

3. Konvolüsyon-Çıkarma Yöntemi:

8.4.3 PET Saçılma Düzeltme Yöntemleri2 soru

1. Single Scatter Simulation (SSS) – Standart yöntem:

2. Enerji penceresi yöntemleri:

8.5 Hareket Düzeltme ve Gating

8.5.1 Hareket Artefaktlarının Kaynakları

Nükleer tıp görüntülemede hareket, görüntü kalitesini ve kantitatif doğruluğu ciddi şekilde etkileyen önemli bir sorundur.

Hareket tipleri:

1. Fizyolojik hareket:

2. İstemli/istemsiz hasta hareketi:

Hareketin görüntüye etkileri:

8.5.2 Gating Teknikleri7 soru

Gating, fizyolojik hareketi "dondurarak" hareket artefaktlarını azaltan bir tekniktir. Akvizisyon sırasında fizyolojik sinyaller (EKG, solunum) kaydedilir ve veriler hareket fazlarına göre ayrılır.

1. Kardiyak Gating (EKG Gating):

Prensip:

Klinik uygulamalar:

LVEF hesaplaması (gated SPECT/PET):

LVEF = (EDV − ESV) / EDV × 100

Gating sorunları:

2. Solunum Gating (Respiratory Gating):

Prensip:

Solunum gating yöntemleri:

Şekil 8.8: EKG gating (8 frame, kısa R-R atımının reddi, hacim eğrisi) ve solunum gatingde faz ve amplitüd bazlı bölme (şematik).

Klinik uygulamalar:

3. Dual Gating (Kardiyak + Solunum):

8.5.3 Hareket Düzeltme Yazılımları1 soru

Gating dışında, yazılım bazlı hareket düzeltme yöntemleri de kullanılmaktadır:

1. Rigid body registration (Katı cisim registrasyonu):

2. Deformable registration (Elastik registrasyon):

3. Data-driven hareket düzeltme:

8.6 Filtreler ve Görüntü İşleme Parametreleri

8.6.1 Frekans Uzayı ve Fourier Dönüşümü

Görüntü işlemede filtreler genellikle frekans uzayında tanımlanır. Fourier dönüşümü, bir görüntüyü uzaysal alandan frekans alanına dönüştürür.

Temel kavramlar:

8.6.2 Filtre Tipleri7 soru

Nükleer tıpta kullanılan filtreler iki ana kategoriye ayrılır:

1. Düzleştirme (Smoothing / Low-pass) Filtreleri:

Yüksek frekansları bastırır → gürültü azalır ancak çözünürlük düşer.

a) Butterworth filtresi:

H(f) = 1 / √(1 + (f/f_c)²ⁿ)

b) Gaussian filtresi:

H(f) = e^{−2π²σ²f²}

c) Hanning (Hann) filtresi:

H(f) = 0.5 + 0.5 cos(π f / f_c), f ≤ f_c

d) Hamming filtresi:

H(f) = 0.54 + 0.46 cos(π f / f_c), f ≤ f_c

2. Keskinleştirme (Sharpening / High-pass) Filtreleri:

3. Rampa filtresi (Ramp filter):

Şekil 8.9: Butterworth filtresinde kesim frekansı ve derecenin etkisi; FBP'de rampa filtresinin pencere fonksiyonuyla birleştirilmesi.

8.6.3 Filtre Seçimi ve Klinik Uygulama

Filtre seçimi, çözünürlük ve gürültü arasındaki dengeye dayanır ve klinik göreve bağlıdır:

Genel ilkeler:

Klinik örnekler:

8.6.4 Görüntü Matris Boyutu ve Piksel Boyutu3 soru

Matris boyutu:

Piksel boyutu ve örnekleme kuralı:

Örnekleme yetersizliği (undersampling):

8.6.5 Projeksiyon Sayısı ve Açısal Örnekleme3 soru

SPECT'te projeksiyon sayısı:

SPECT akvizisyon yörüngesi:

Şekil 8.10: Dairesel ve vücut konturu SPECT yörüngeleri (şematik).

8.7 Veri Toplama Modları ve Nicel Analiz

8.7.1 Veri Toplama Modları1 soru

Temel kavramlar:

Statik veri toplama:

Dinamik veri toplama:

Liste modu (list mode):

Kapılı (gated) veri toplama:

Dijital Görüntü Biçimleri ve Bilgi Sistemleri

8.7.2 Nicel Görüntülemenin Ön Koşulları ve Hata Kaynakları4 soru

ROI/VOI tanımlama ve SUV:

8.7.3 Kinetik Modelleme ve Parametrik Görüntüleme4 soru

Şekil 8.11: İki doku kompartman modeli ile Patlak ve Logan grafiksel analizleri (şematik).

Board İçin Kritik Noktalar