%ricostruzione con 8 telecamere mediante minimi quadrati pesati con
%eliminazione dei raggi

%Costruisco le telecamere

T1P=[a{1}' b{1}' c{1}' num_pixel_centri{1}' x_f_as(1)*ones(length(a{1}),1) y_f_as(1)*ones(length(a{1}),1)];
T2P=[a{2}' b{2}' c{2}' num_pixel_centri{2}' x_f_as(2)*ones(length(a{2}),1) y_f_as(2)*ones(length(a{2}),1)];
T3P=[a{3}' b{3}' c{3}' num_pixel_centri{3}' x_f_as(3)*ones(length(a{3}),1) y_f_as(3)*ones(length(a{3}),1)];
T4P=[a{4}' b{4}' c{4}' num_pixel_centri{4}' x_f_as(4)*ones(length(a{4}),1) y_f_as(4)*ones(length(a{4}),1)];
T5P=[a{5}' b{5}' c{5}' num_pixel_centri{5}' x_f_as(5)*ones(length(a{5}),1) y_f_as(5)*ones(length(a{5}),1)];
T6P=[a{6}' b{6}' c{6}' num_pixel_centri{6}' x_f_as(6)*ones(length(a{6}),1) y_f_as(6)*ones(length(a{6}),1)];
T7P=[a{7}' b{7}' c{7}' num_pixel_centri{7}' x_f_as(7)*ones(length(a{7}),1) y_f_as(7)*ones(length(a{7}),1)];
T8P=[a{8}' b{8}' c{8}' num_pixel_centri{8}' x_f_as(8)*ones(length(a{8}),1) y_f_as(8)*ones(length(a{8}),1)];

%Primo passo dell'albero unisco solo 2 telecamere

R11Pa2=cell(1,2);
R21Pa2=cell(1,2);
R31Pa2=cell(1,2);
R41Pa2=cell(1,2);

%Associo le telecamere ricordando che al primo passo nn vi sono
%ricostruzioni

R11Pa2{1}=T1P;
R11Pa2{2}=T2P;
R21Pa2{1}=T3P;
R21Pa2{2}=T4P;
R31Pa2{1}=T5P;
R31Pa2{2}=T6P;
R41Pa2{1}=T7P;
R41Pa2{2}=T8P;

P11Pa2=[-1 -10 -1];

%dimensione stanza
dim_stanza_x=10000;
dim_stanza_y=10000;

%secondo passo dell'albero
set(S.tP1R1,'string','-----','visible','on')
set(S.tP1R2,'string','-----','visible','on')
set(S.tP1R3,'string','-----','visible','on')
set(S.tP1R4,'string','-----','visible','on')

tic
[R12Pa2 P12Pa2]=associa2_pes(R11Pa2, R21Pa2, P11Pa2, P11Pa2,epsilon,dim_stanza_x,dim_stanza_y,d_max);
t1=toc;
set(S.tP2R1,'string',t1,'visible','on')
tic
[R22Pa2 P22Pa2]=associa2_pes(R31Pa2, R41Pa2, P11Pa2, P11Pa2,epsilon,dim_stanza_x,dim_stanza_y,d_max);
t2=toc;
set(S.tP2R1,'string',t2,'visible','on')
%terzo passo dell'albero
tic
[R14Pa2 P14Pa2]=associa2_pes(R12Pa2, R22Pa2, P12Pa2, P22Pa2,epsilon,dim_stanza_x,dim_stanza_y,d_max);
t3=toc;
set(S.tP3R1,'string',t3,'visible','on')
tmpo=t3+max([t1 t2]);
set(S.Ttot,'string',tmpo,'visible','on')



