-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathrecon_all.m
70 lines (48 loc) · 1.39 KB
/
recon_all.m
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
% recon_all( dataDir )
%% Get list data directories to process
if ~exist( 'dataDir', 'var' ),
dataDir = fullfile( fileparts( mfilename ), '..', 'data' );
if ~exist( dataDir, 'dir' )
error( 'data directory doesn''t exist, %s\n', dataDir )
end
end
D = dir( dataDir );
isInc = [ D.isdir ]; % include all directories
% loop to identify hidden directories
ind = find(isInc);
for iD = ind,
% on OSX, hidden directories start with a dot
isInc(iD) = ~strcmp(D(iD).name(1),'.');
if isInc(iD) && ispc
% check for hidden Windows directories - only works on Windows
[~,stats] = fileattrib(fullfile(dataDir,D(iD).name));
if stats.hidden
isInc(iD) = false;
end
end
if isInc(iD),
isInc(iD) = strcmp( D(iD).name(1:4), 'fcmr' );
end
end
D = D(isInc);
%% Reconstruct
isFailed = false( 1, numel(D) );
for iD = 1:numel(D),
fcmrNo = 0;
seriesNo = 0;
try
fcmrNo = str2double( D(iD).name(5:7) );
seriesNo = str2double( D(iD).name(9:10) );
main_fcmr_cine( fcmrNo, seriesNo );
catch ME
isFailed(iD) = true;
warning( 'Caught error reconstructing %s:', D(iD).name )
disp( ME ),
end
end
fprintf( '\n\n\nFailed to reconstruct:\n' )
indFailed = find ( isFailed );
for iF = indFailed,
fprintf( ' %s\n', D(iF).name )
end
% end