25 lines
612 B
Matlab
25 lines
612 B
Matlab
clear;
|
|
close all;
|
|
|
|
load spheres;
|
|
|
|
nb_lignes = size(spheres,1);
|
|
nb_colonnes = size(spheres,2);
|
|
|
|
centre = [nb_lignes/2 nb_colonnes/2];
|
|
rayon = min(nb_colonnes,nb_lignes)/2;
|
|
|
|
[X,Y] = meshgrid(1:nb_colonnes,1:nb_lignes);
|
|
marge_rayon = 5;
|
|
exterieur_masque = (X-rayon).^2+(Y-rayon).^2>(rayon-marge_rayon)^2;
|
|
|
|
s = etalonnage(spheres,exterieur_masque,centre,rayon);
|
|
|
|
[theta,phi] = conversion(s);
|
|
plot(theta,phi,'o','Color','r','LineWidth',4,'MarkerSize',10);
|
|
xlabel('$\theta$','Interpreter','Latex','FontSize',20);
|
|
ylabel('$\phi$','Interpreter','Latex','FontSize',20);
|
|
axis([0,pi/2,-pi,pi]);
|
|
|
|
save eclairages theta phi;
|