DBP:Harvard:Software:Wish:3T 8
From NAMIC Wiki
Home < DBP:Harvard:Software:Wish:3T 8
9/28/05
- Updates, KangUk?
(9/21/05)
- 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.