Please enable JavaScript.
Coggle requires JavaScript to display documents.
IMAGE ANALYSIS AND PROCESSING: METODI TRADIZIONALI - Coggle Diagram
IMAGE ANALYSIS AND PROCESSING: METODI TRADIZIONALI
MEDICAL IMAGING
rappresentazione immagini corpo umano
UTILIZZO RADIAZIONI
IONIZZANTI
alta energia
rischio reazione DNA
raggi gamma
raggi x
ultra violetti
RADIATION DOSE
energia assorbita per unità di massa
$$ D= \frac {E}{M}[mGy]$$
NON IONIZZANTI
bassa energia
rischio zero
ultrasuoni
thermal imaging
MEDICINA NUCLEARE
utilizzo radio isotopi
UTILIZZO RADIOISOTOPI
EMETTO RADIAZIONI
RILEVO EMISSIONE
Immagini a bassa qualità
vediamo solo
distribuzione isotopi
CATEGORIE
Spect
single photon emission computed tomography
detectors riconstruiscono distribuzione spaziale 3d
raggi
gamma
Ultrasuoni
ultrasuoni rimbalzano su superfici
tessuti molli
X ray Imaging
tubo radiogeno
proiezione paziente in una determinata direzione
coefficiente attenuazione dipende da densità e composizione chimica
Interazione
radiazione materia
$$I(Z)= -\int_{z0}^{z1} \mu(z) dz $$
VISUALIZZIAMO IL COEFFICIENTE DI ATTENUAZIONE
Compute Tomography
Tac
movimento complessivo a spirale
algoritmi di back proiezione per ricostruire mappa coeficienti
qualità elevata, rischio radiazioni alto
100 radiation dose
OTTIMIZZARE QUALITA'
ALARA PRINCIPLE
As Low As Reasonable Achievable
Risonanza magnetica
ricezione proprietà diffusione protoni
-campo magnetico 3 Tesla
onde radio 64 Mhz
Pet
Positron emission tomography
tempo di volo brevissimo
generazione 2 fotoni in direzioni opposte
IMAGE QUALITY
Indici quantitativi
Signal To Noise Ratio
$$ SNR= \frac {\sum (x_i - \hat x_{BG})}{\sigma_{BG}} $$
in caso di random noise
SNR proporzionale a radice di N
POSSO DIMINUIRE RUMORE CON ACQUISIZIONI MULTIPLE
dimensioni oggetto influenzano SNR
Contrast to Noise Ratio
$$CNR_{AB} = |SNR_A - SNR_B|$$
Risoluzione Spaziale
la + piccola distanza tra 2 oggetti tali da poter essere risolti
Line Spread Function
approssimazione Gaussiana
$$LSF(y) = \frac {1}{\sqrt {2 \pi \sigma^2}}{exp}(- \frac {(y-y_0)^2}{2 \sigma^2})$$
inversa
FWHM
$$2\sqrt {2ln2\sigma}$$
Point SPread Function
non in una direzione, tante LSF
convoluzione oggetto e PSF per ottenere immagine
kernel applicato a tutti i punti dell'immagine
$$ I(x,y,z) = O(x,y,z) \star PSF(x,y,z)$$
Freqienza Spaziale
n linee per unità di lunghezza
se aumenta frequenza devo aumentare risoluzione
$$ F=1/(2 \Delta) [mm^-1]$$
una qualsiasi immagine può essere decomposta in funzioni sinuisoidali
$$F(x,y)= \int_{-\inf}^\inf \int_{-\inf}^{inf} f(x,y)e^{-j2\pi (ux+vy)}dxdy$$
Modulation Transfer Function
misura response del sistema di imaging a differenti frequenze spaziali
La trasformata di Fourier del PSF
LIMITING SPAZIAL RESOLUTION: circa al 10%
CONVOLUTION THEOREM
$$I=O \star PSF = F^{-1} (F(O) * MTF) $$
più facile fare un prodotto che una convoluzione
Noise
sorgenti
rilevatore
strutturale
radiazione
quantum mottle
caso fortunato
Random Noise
$$ \sigma= \sqrt \frac {\sum_i (I_i - \hat I )^2}{N-1} $$
Segnale
Somma intensità Pixel
Contrasto
differenza segnale fra 2 tessuti del corpo
sparisce la dipendenza dalla dimensione
dipendenza proprietà intrinseche tessuti
$$C_{AB} =|S_{A} -S_{B}|$$
tipologie
Detector
Intrinseco
Agente (corrispondente e specifico)
Tecniche e Settings di Imaging
Blurring
Cause
motion unsharpness
geometric unsharpness
subject unsharpness
receptor unsharpness
IMAGE PROCESSING
IMAGE FILTERING
apllico filtro alla trasformata di Fourier e poi torno indietro
low pass: diminuisco maggior parte del rumore
high pass: bordi più nitidi
SPATIAL FILTERING
convoluzione tra filtro e immagine
es Kernel 1x1, per stretching contrasto o thresholding
negativi: s=L-1-r
logaritmiche: s=c log (1+r)
Power low: s=c r^gamma
HISTOGRAM PROCESSING
conta n° pixel che hanno una determinata intensità
usualmente normalizzato all'image size(distribuzione di probabilità)
Equalizzazione: appiattire il più possibile l'istogramma
Fuzzy: attribuiamo pixel a categorie con una determinata probabilità
FILTRI DENOISE
per diminuire rumore
Averaging: appiatisce immagine, si perde risoluzione spaziale
Median: efficaci in risoluzione sale e pepe
FILTRI SHARPENING
aumento nitidezza bordi utilizzando derivata
gradient: utilizzo derivate prime, edge detection
Lapalissiana: derivate seconde, metto in risalto i bordi