%This script will process all of the ultrasound data for the flow phantoms.

%For full details see "Ultrasound Imaging Velocimetry with interleaved
%images for improved pulsatile arterial flow measurements: a new correction
%method, experimental and in vivo validation", J Royal Soc Interface, 2017

%K H Fraser, C Poelma 2016

%--------phantom a--------------
expectedHR=60;
depth=0.02;
flow_region_y=500:720;

delta=15;
frames=160;
filename = ['../data/flow_phantom/a/UIV/2cm/' int2str(delta) '/phantom_a_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=31;
frames=155;
filename = ['../data/flow_phantom/a/UIV/2cm/' int2str(delta) '/phantom_a_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=47;
frames=150;
filename = ['../data/flow_phantom/a/UIV/2cm/' int2str(delta) '/phantom_a_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=63;
frames=145;
filename = ['../data/flow_phantom/a/UIV/2cm/' int2str(delta) '/phantom_a_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=79;
frames=140;
filename = ['../data/flow_phantom/a/UIV/2cm/' int2str(delta) '/phantom_a_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=95;
frames=135;
filename = ['../data/flow_phantom/a/UIV/2cm/' int2str(delta) '/phantom_a_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

frames=330;
filename = ['../data/flow_phantom/a/UIV/2cm/conventional/phantom_a_128lines_conventional_20000depth7000000freq'];
pulsepiv1_phantom_f(depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_phantom_f(depth,frames,filename)

%--------phantom b--------------
expectedHR=200;
depth=0.02;
flow_region_y=400:650;

delta=7;
frames=165;
filename = ['../data/flow_phantom/b/UIV/2cm/' int2str(delta) '/phantom_b_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=15;
frames=160;
filename = ['../data/flow_phantom/b/UIV/2cm/' int2str(delta) '/phantom_b_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=31;
frames=155;
filename = ['../data/flow_phantom/b/UIV/2cm/' int2str(delta) '/phantom_b_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=47;
frames=150;
filename = ['../data/flow_phantom/b/UIV/2cm/' int2str(delta) '/phantom_b_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=63;
frames=145;
filename = ['../data/flow_phantom/b/UIV/2cm/' int2str(delta) '/phantom_b_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

frames=330;
filename = ['../data/flow_phantom/b/UIV/2cm/conventional/phantom_b_128lines_conventional_20000depth7000000freq'];
pulsepiv1_phantom_f(depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_phantom_f(depth,frames,filename)

%--------phantom c--------------
expectedHR=200;
depth=0.02;
flow_region_y=450:680;

delta=7;
frames=165;
filename = ['../data/flow_phantom/c/UIV/2cm/' int2str(delta) '/phantom_c_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=15;
frames=160;
filename = ['../data/flow_phantom/c/UIV/2cm/' int2str(delta) '/phantom_c_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=31;
frames=155;
filename = ['../data/flow_phantom/c/UIV/2cm/' int2str(delta) '/phantom_c_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=47;
frames=150;
filename = ['../data/flow_phantom/c/UIV/2cm/' int2str(delta) '/phantom_c_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=63;
frames=145;
filename = ['../data/flow_phantom/c/UIV/2cm/' int2str(delta) '/phantom_c_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

frames=330;
filename = ['../data/flow_phantom/c/UIV/2cm/conventional/phantom_c_128lines_conventional_20000depth7000000freq'];
pulsepiv1_phantom_f(depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_phantom_f(depth,frames,filename)

%--------phantom d--------------
expectedHR=60;
depth=0.02;
flow_region_y=500:750;

delta=7;
frames=165;
filename = ['../data/flow_phantom/d/UIV/2cm/' int2str(delta) '/phantom_d_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=15;
frames=160;
filename = ['../data/flow_phantom/d/UIV/2cm/' int2str(delta) '/phantom_d_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=31;
frames=155;
filename = ['../data/flow_phantom/d/UIV/2cm/' int2str(delta) '/phantom_d_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=47;
frames=150;
filename = ['../data/flow_phantom/d/UIV/2cm/' int2str(delta) '/phantom_d_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=63;
frames=145;
filename = ['../data/flow_phantom/d/UIV/2cm/' int2str(delta) '/phantom_d_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

frames=330;
filename = ['../data/flow_phantom/d/UIV/2cm/conventional/phantom_d_128lines_conventional_20000depth7000000freq'];
pulsepiv1_phantom_f(depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_phantom_f(depth,frames,filename)

%--------phantom e--------------
expectedHR=60;
depth=0.02;
flow_region_y=600:900;

delta=7;
frames=165;
filename = ['../data/flow_phantom/e/UIV/2cm/' int2str(delta) '/phantom_e_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=15;
frames=160;
filename = ['../data/flow_phantom/e/UIV/2cm/' int2str(delta) '/phantom_e_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=31;
frames=155;
filename = ['../data/flow_phantom/e/UIV/2cm/' int2str(delta) '/phantom_e_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=47;
frames=150;
filename = ['../data/flow_phantom/e/UIV/2cm/' int2str(delta) '/phantom_e_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=63;
frames=145;
filename = ['../data/flow_phantom/e/UIV/2cm/' int2str(delta) '/phantom_e_128lines_' int2str(delta) 'delta_20000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

frames=330;
filename = ['../data/flow_phantom/e/UIV/2cm/conventional/phantom_e_128lines_conventional_20000depth7000000freq'];
pulsepiv1_phantom_f(depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_phantom_f(depth,frames,filename)

depth=0.06;
flow_region_y=600:900;

delta=7;
frames=80;
filename = ['../data/flow_phantom/e/UIV/6cm/' int2str(delta) '/phantom_e_128lines_' int2str(delta) 'delta_60000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=15;
frames=75;
filename = ['../data/flow_phantom/e/UIV/6cm/' int2str(delta) '/phantom_e_128lines_' int2str(delta) 'delta_60000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=31;
frames=73;
filename = ['../data/flow_phantom/e/UIV/6cm/' int2str(delta) '/phantom_e_128lines_' int2str(delta) 'delta_60000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=47;
frames=68;
filename = ['../data/flow_phantom/e/UIV/6cm/' int2str(delta) '/phantom_e_128lines_' int2str(delta) 'delta_60000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

delta=63;
frames=64;
filename = ['../data/flow_phantom/e/UIV/6cm/' int2str(delta) '/phantom_e_128lines_' int2str(delta) 'delta_60000depth7000000freq'];
pulsepiv1_interlaced_phantom_f(delta,depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_interlaced_phantom_f(delta,depth,frames,filename)

frames=126;
filename = ['../data/flow_phantom/e/UIV/6cm/conventional/phantom_e_128lines_conventional_60000depth7000000freq'];
pulsepiv1_phantom_f(depth, frames,filename, flow_region_y,expectedHR)
pulsepiv_phantom_f(depth,frames,filename)