Oppgave 3
Mandag 01.10.01
Tema: Clutterfilter ved ultralyd fargedoppler
Oppgaven tar utgangspunkt i artikkelen "Clutter Rejection Filters in Color Flow Imaging "
Endel av dette stoffet er også beskrevet i "Waves...", kap. 10.3 E; Bakgrunnsstoff
Eksempel på ultralyd signal fra fargedoppler ligger på oppg3Data.mat
Fila inneholder kompleks demodulert signal iq fra et bilde tatt opp med lineær array probe.
Senterfrekvens 7.5 MHz. For hver av de 66 strålene er det skutt 7 pulser, med prf=750 Hz
iq(k,r,b) gir sample fra puls nr. k, range r, stråle nr. b
Oppgaven går ut på å designe et filter som ekstraherer signalet fra blod i bevegelse, og vise det som et gråtone-bilde, ved siden av det ufilterte signalet.
Merknader:
iq er lagret som 16 bit heltall for å spare plass. Utfør iq=double(iq); for å konvertere til flyttall presisjon.
Amplitude-bilde av vevssignalet kan man få ved å beregne effekt i signalet fra hver "packet"", f.eks
amp=10*log10(mean(abs(iq));
Se på real og imaginærdel til signalet fra et punkt i bildet:
r=100;b=33;t=(0:6)/prf;
iq1=squeeze(iq(:,r,b));
plot(t,real(iq1),t,imag(iq1));
Filtrer iq1 med høypass filter, og plot resultatet
Ufør filtering for alle punkter i bildet, og vis effekt av utgangs-signalet som (log) gråtonebilde.