optimized_multi_coil/example.m at master · hbbdyzf/optimized_multi_coil · 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

%
% This code shows examples of coils optimization
% Please see 'READ_ME' for more details
%
% Author: Ali Aghaeifar <[email protected]>
% Created: April 2019
%

addpath(genpath('.\dependencies'));

%% Generate a Profile
% a profile contains some of the required settings for optimization. Rest of the settings can be changed by modyfing 'opt_profile_manager.m'
profile  = opt_profile_manager('shim_mode', {'global'},...              % Skope of optimization, global or slice-wise
                               'shim_algo_inner', {'pinv'},...          % The algorithm used for calculation of shim currents; 'pinv' for unconstrained shimming and 'lsqlin' for constrained shimming
                               'shim_algo_outer', {'sqp'},...		    % The solver used for optimization of coils' position and size
                               'coilShape', {'square'}, ...             % Shape of the coils, 'square' or 'circle'
                               'nonlconCoef', [0],...				    % Nonlinear constraint weight. Here is the coils overlapping.
                               'shimDCL', [true],...                    % Adjust coil's current based on the coil's size.
                               'coil_no_row', [32,4,30], ...            % Initial coils arrangment, [number of coils, number of rows, size of the coils]
                               'prefix', 'Unconstrained_CoilSize30');   % A prefix for the name of destination folder (results will be saved there)

opt_main(profile); % execute optimization based on the profile

%%
% Several profiles can be generated together
profiles = opt_profile_manager('shim_mode', {'global'},...
                               'shim_algo_inner', {'lsqlin', 'pinv'},...
                               'shim_algo_outer', {'sqp'},...
                               'coilShape', {'square'}, ...
                               'nonlconCoef', [0],...
                               'shimDCL', [true],...
                               'coil_no_row', [96,6,25; 32,4,30], ...
                               'prefix', 'Multiple_profiles');

% Eexecute profiles one by one (internally). Parallelization will be on individual coils (recommended).
opt_main(profiles);

% or
parfor i=1:numel(profiles)
   opt_main(profiles{i})
end