%-----------------------------------------------------------
% Carlos Cabral 2003 - Dep.UA

% COEFICIENTE DE CORRELAÇAO usando a linearizaçao no plano complexo
%----------------------------------------------------------

function [cond]=COEFCORR_CIRCLINEAR(c,i6,pontos,cond)
% translaçao dos pontos de modo a que o centro da circunferencia obtida
% pelo metodo de reimann coincida com a origem do sistema de coordenadas 
trans=[1 0 -c(1)
       0 1 -c(2)
       0 0 1];
pontos(3,:)=ones(1,length(pontos(1,:)));   
pontos_trans=trans*pontos;
pontos_trans(3,:)=[];
pontos(3,:)=[];
                
%transformaçao dos pontos em numeros complexos (mudança para o espaço complexo)
pontos_cplx=complex(pontos_trans(1,:),pontos_trans(2,:));
pontos_cplx_log=log(pontos_cplx); %logaritmizaçao dos mesmos pontos
                
%plot(real(pontos_cplx),imag(pontos_cplx),'.'), grid on

%ciclo para rotaçao de todos os pontos para que seja possivel o calculo do coeficiente de correlaçao
pontos_log_rot=[real(pontos_cplx_log);imag(pontos_cplx_log)];
pontos_log_rot(3,:)=ones(1,length(pontos_log_rot));
teta1=-45*pi/180;
rot=[cos(teta1) -sin(teta1) 0
     sin(teta1) cos(teta1)  0
     0          0           1]; 
pontos_rot=rot*pontos_log_rot;
pontos_rot(3,:)=[];
                
%plot(pontos_rot(1,:),pontos_rot(2,:),'.'), axis('equal'),grid on
coef_ideal=corrcoef(pontos_rot');
cond(i6,3)=coef_ideal(1,2);
%------------------------------------------------------------------------       
                