function [r2] = Grille_TPSP(rayon,pts)

%-------------------------------------------------------------------------------------
% Fonction permettant de calculer les positions des capteurs de surface, par découpage
% "équi-réparti" de la surface de la sphère du scalp (découpage selon theta et phi)
%-------------------------------------------------------------------------------------
% entrées :
% rayon -> rayon de la sphère du scalp
% delta -> nombre de valeurs de chaque angle
%
% sortie :
% r     -> matrice (3xN) des coord cartésiennes 3D des capteurs de surface
%-------------------------------------------------------------------------------------

%II. Construction d'une sphère
sphe2 = [rayon*ones(1,pts^2);...
        kron(ones(1,pts),0:2*pi/(pts-1):2*pi);... %theta1 theta2 theta3 theta4
        kron(0:pi/(pts-1):pi,ones(1,pts))];       %phi1   phi2   phi3   phi4


% III. Transformer les coord sphériques en coord cartésiennes 3D :
r2(1,:) = sphe2(1,:) .* sin(sphe2(3,:)) .* cos(sphe2(2,:)) ; % x=r*sin(theta)*cos(phi) 
r2(2,:) = sphe2(1,:) .* sin(sphe2(3,:)) .* sin(sphe2(2,:)) ; % y=r*sin(theta)*sin(phi)
r2(3,:) = sphe2(1,:) .* cos(sphe2(3,:)) ; % z=r*cos(theta)



