admmOpticalFlow/loadData.m at master · ndwork/admmOpticalFlow · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100

function [x,eta,rho,offset] = loadData( datacase )

  offset = 0;

  switch datacase
    case 1
      load '/Volumes/Seagate2TB/Data/mriCardiacNavs/Okai/esp3dn25148.mat';
      % The max dimensions of x are (128,128,64,508);
      x = x / 500;
      eta = 3d-9;   % Great for datacase 1
      rho = 1d-3;
      offset = 10;
    case 2
      load '/Volumes/Seagate2TB/Data/cardiacRecons/res2mm/img318143.mat';
      % The max dimensions of m are (120,120,80,20);
      x = m;  clearvars m;
    case 3
      % Not very much of a heartbeat
      load '/Volumes/Seagate2TB/Data/cardiacRecons/res1.2mm/img32013.mat';
      % The max dimensions of m are (256,256,112,6);
      x = m;  clearvars m;
    case 4
      x = loadSerenaData(1);
      x = x / max( x(:));
    case 5
      x = loadSerenaData(2);
      x = x / max( x(:));
    case 6
      x = loadSerenaData(3);
      x = x / max( x(:));
    case 7
      x = loadSerenaData(4);
      x = x / max( x(:));
    case 8
      load '/Volumes/Seagate2TB/Data/mriCardiacNavs/Jieying/esp3dn29696_wvt1_all.mat';
      x = x / 500;
    case 9
      load '/Volumes/Seagate2TB/Data/mriCardiacNavs/Jieying/esp3dn29696_wvt4_all.mat';
      x = x / 300;
      offset = 10;
      eta = 3d-10;    % Great for datacase 9
      rho = 1d-3;
    case 10
      load '/Volumes/Seagate2TB/Data/cardiacRecons/okais/img20130320_44544_1057.mat';
      x = m / 0.06;  m=[];
      eta = 3d-1;   % Great for datacase 1
      rho = 1d-3;
    case 11
      load '/Volumes/Seagate2TB/Data/mriCardiacNavs/Jieying/esp3dn02048_11_60_wvt.mat';
      x = combineCoils( x );
      x = x / 380;
      eta = 8d-10;   % Great for datacase 1
      rho = 1d-3;
  end
end


function imgs = combineCoils( data )
  sImgFiles = size( data );
  nCoils = sImgFiles(4);
  nPhases = sImgFiles(5);

  imgs = zeros([ sImgFiles(1:3), nPhases ]);
  for phase = 1:nPhases
    thisImg = abs(data(:,:,:,1,phase)).^2;
    for coil=2:nCoils
      % Combine images from each coil
      thisImg = thisImg + abs(img(:,:,:,coil,phase)).^2;
    end
    imgs(:,:,:,phase) = sqrt( thisImg );
  end
end


function imgs = loadSerenaData( fileIndx )
  imgFiles = cell(4,1);
  imgFiles{1} = '/Volumes/Seagate2TB/Data/cardiacImgs/P02560_img.mat';
  imgFiles{2} = '/Volumes/Seagate2TB/Data/cardiacImgs/P27648_img.mat';
  imgFiles{3} = '/Volumes/Seagate2TB/Data/cardiacImgs/P62976_img.mat';
  imgFiles{4} = '/Volumes/Seagate2TB/Data/cardiacImgs/P712133_img.mat';

  img = 0;
  load( imgFiles{fileIndx} );

  sImgFiles = size( img );
  nCoils = sImgFiles(4);
  nPhases = sImgFiles(5);
  imgs = zeros([ sImgFiles(1:3), nPhases ]);

  for phase = 1:nPhases
    thisImg = abs(img(:,:,:,1,phase)).^2;
    for coil=2:nCoils
      % Combine images from each coil
      thisImg = thisImg + abs(img(:,:,:,coil,phase)).^2;
    end
    imgs(:,:,:,phase) = sqrt( thisImg );
  end

end