Commit f4111c66 authored by p.jhagema's avatar p.jhagema
Browse files

adapt parameters for mg wire exabple

parent e3916c41
Pipeline #266010 passed with stage
in 52 seconds
......@@ -96,8 +96,8 @@ p.PMset.plotInterval = 10;
p.conSet = PConstraints.defaultSettings();
p.conSet.upperAmplitude = 1.0101;
p.conSet.lowerAmplitude = 0.553;
p.conSet.smoothAmplitudesBy = 2;
p.conSet.smoothPhasesBy = 2;
p.conSet.smoothAmplitudesBy = 0;
p.conSet.smoothPhasesBy = 0;
p.conSet.projectOnNegativity = 0;
p.conSet.unwrapPhases = 0;
......
......@@ -24,16 +24,17 @@ seq_conv_hist{1} = reco.stats.convergenceHistory;
%% we refinie the support from the final reconstruchtion
if(1)
p.supp = (angle(reco.x)) < -0.2;
p.supp = abs(angle(reco.x)) > 0.2;
p.supp = padToSize(cropToCenter(p.supp, p.data_size), p.recSize, 0);
se = strel('disk', 70);
p.supp = p.supp | imdilate(supp, se);
p.supp = bwareafilt(p.supp, 1, 'largest');
se = strel('disk', 10);
p.supp = imerode(p.supp, se);
se = strel('disk', 10);
p.supp = imdilate(p.supp, se);
se = strel('disk', 10);
p.supp = imerode(p.supp, se);
p.supp = imfill(p.supp, 'holes');
p.supp = bwareafilt(p.supp, 1, 'largest');
......@@ -51,11 +52,13 @@ p.algSet.gamma = 0.9;
% set or unset phase unwrapping, be warned it take time...
p.conSet.unwrapPhases = 1;
if(p.conSet.unwrapPhases)
p.conSet.projectOnNegativity = 1;
p.conSet.smoothAmplitudesBy = 0;
p.conSet.smoothPhasesBy = 0;
p.conSet.projectOnNegativity = 0;
end
PS = PConstraints(p.conSet) * PSupp;% * PStat;
p.algSet.iterations = 1000;
p.algSet.iterations = 10000;
projAlg = projectionAlgorithm(PM, PS, guess, p.algSet);
reco = projAlg.execute(p.algSet.iterations);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment