View source for DBP:Harvard:Software:Wish:3T 8

* Updates, KangUk?

<br />

<br />

=== '''(9/21/05)''' ===

<div class="thumb tleft"><div style="width: 182px">[[Image:Raw_17_fibertracking.png|[[Image:180px-Raw_17_fibertracking.png|raw_17]]]]<div class="thumbcaption"><div class="magnify" style="float: right">[[Image:Raw_17_fibertracking.png|[[Image:magnify-clip.png|Enlarge]]]]</div>raw_17</div></div></div><div class="thumb tleft"><div style="width: 182px">[[Image:Smoothing_17_fibertrac.png|[[Image:180px-Smoothing_17_fibertrac.png|smoothing_17]]]]<div class="thumbcaption"><div class="magnify" style="float: right">[[Image:Smoothing_17_fibertrac.png|[[Image:magnify-clip.png|Enlarge]]]]</div>smoothing_17</div></div></div><div class="thumb tleft"><div style="width: 182px">[[Image:Raw_20_fibertrack.png|[[Image:180px-Raw_20_fibertrack.png|raw_20]]]]<div class="thumbcaption"><div class="magnify" style="float: right">[[Image:Raw_20_fibertrack.png|[[Image:magnify-clip.png|Enlarge]]]]</div>raw_20</div></div></div><div class="thumb tleft"><div style="width: 182px">[[Image:Smooth_20_fibertrack.png|[[Image:180px-Smooth_20_fibertrack.png|smoothing_20]]]]<div class="thumbcaption"><div class="magnify" style="float: right">[[Image:Smooth_20_fibertrack.png|[[Image:magnify-clip.png|Enlarge]]]]</div>smoothing_20</div></div></div>

* '''Image parameter'''

'''1.7mm thickness:'''<br /> 0.86 x 2.29 x 1.7 <br />'''2.0mm thickness:'''<br /> 1.05 x 3.18 x 2.0 <br />

* '''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 [http://lmi.bwh.harvard.edu/personnel/gk/17-20/ here]

<br />

=== '''(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.