DBP:Harvard:Software:Wish:3T 8-24-05

From NAMIC Wiki
Jump to: navigation, search
Home < DBP:Harvard:Software:Wish:3T 8-24-05

9/28/05

  • Updates, KangUk?



(9/21/05)

raw_17
Enlarge
raw_17
smoothing_17
Enlarge
smoothing_17
raw_20
Enlarge
raw_20
smoothing_20
Enlarge
smoothing_20
  • Image parameter

1.7mm thickness:
0.86 x 2.29 x 1.7
2.0mm thickness:
1.05 x 3.18 x 2.0

  • Matlab Code

if 0

 fName = '/projects/fmri/pfiles-and-structural-files/dti/For-Bob/recon_17_iso/slicer/as_8_15_05_17_iso_DWI'
 N = 259
 outName = '/projects/fmri/pfiles-and-structural-files/dti/For-Bob/recon_17_iso/slicer/filt_as_8_15_05_17_iso_DWI'

else

 fName = '/projects/fmri/pfiles-and-structural-files/dti/For-Bob/recon_20_iso/slicer/as_8_15_05_20_iso_DWI'
 N = 364
 outName = '/projects/fmri/pfiles-and-structural-files/dti/For-Bob/recon_20_iso/slicer/filt_as_8_15_05_20_iso_DWI'

end

tmp = readSlices(fName,1,N);

gy = fspecial('gauss',[1 7],1); gy = gy./sum(gy(:)); gx = gy'; gz = permute(gx,[3 2 1]); surfl(gx*gy);

s=size(tmp,3)/7; imsize = [256 256 s]; c = ones(imsize);

ca = convn(convn(convn(c,gx,'same'),gy,'same'),gz,'same');

for k=1:7

 k
 im = tmp(:,:,k:7:end);
 size(im)
 res = convn(convn(convn(im,gx,'same'),gy,'same'),gz,'same');
 res = res./ca;
 %put back into tmp
 tmp(:,:,k:7:end) = res;

end

writeSlices(tmp,outName);

(9/14/05)

Gordon's details here


(8-24-05)

KangUk and Marek have been experimenting with higher resolution scans using the 3T Magnet, and KangUk will work on comparing quantitative measures between these new protocols.