-
Notifications
You must be signed in to change notification settings - Fork 13
/
bmode_carotid_experiment.m
61 lines (46 loc) · 1.85 KB
/
bmode_carotid_experiment.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
% Reproduce the carotid experiments of Section V.C of the paper
% "Non-stationary Blur Evaluation for Ultrasound Image Restoration",
% submitted to IEEE Transactions on Computational imaging
% Author: Adrien Besson
% Signal Processing Laboratory 5 (LTS5), EPFL
% email address: adrien.besson@epfl.ch
% July 2018
clear all;
close all;
clc
%-- Add path
addpath(genpath('utils'));
%-- Output filename
list_filename_output = {'results/trf_varying_invivo_fr_1.5.mat','results/trf_est_invivo_fr_1.5.mat', 'results/trf_constant_invivo_fr_1.5.mat'; 'results/trf_varying_invivo_1.3.mat','results/trf_est_invivo_1.3.mat', 'results/trf_constant_invivo_1.3.mat'};
%-- Parameters
flag_carotid_list = [2 1];
flag_psf_meth = [3, 2, 1];
list_p = [1.5, 1.3];
maximum_iterations = 100;
flag_display = 0;
regularization_parameters = [5e-3, 700, 700; 3e-2, 700, 1500];
%-- Loop to generate the B-mode images of the carotids
for kk = 1:numel(flag_carotid_list)
%-- Considered carotid experiment
flag_carotid = flag_carotid_list(kk);
%-- Considered regularization parameters
reg_param_caro = regularization_parameters(kk, :);
%-- Considered value of p
p = list_p(kk);
%-- Considered output filenames
output_filenames = list_filename_output(kk, :);
for pp = 1:numel(flag_psf_meth)
%-- Current PSF method
psf_meth = flag_psf_meth(pp);
%-- Current output filename
filename_out = output_filenames{pp};
%-- Current regularization parameter
reg_parameter = reg_param_caro(pp);
%-- Experiment
deconvolution_carotid(flag_carotid, psf_meth, 0, filename_out, p, reg_parameter, maximum_iterations);
end
%-- Display B-mode images
display_bmode(output_filenames)
%-- Compute metrics
compute_metrics_carotid(output_filenames);
end