Discussion :

[jerome828]

Bonjour à tous,

Je bute sur un problème à priori simple mais j'ai vraiment du mal à m'en sortir.
J'aimerai générer une image avec une fonction en cos^4(teta) qui part du centre de l'image. En imagerie on apelle ça un dome radiométrique et répresente la variation d'éclairement sur un capteur CCD en fonction de l'angle theta entre un pixel et le centre optique.

ça suit une loi en E(teta) = Ecentre*cos(teta)^4

avec Ecentre l'éclairement au centre du capteur
En gros plus l'on séloigne du centre et moins le capteur reçoit d'éclairement.

J'utilise ce code là mais ce n'est pas ce que je souhaite obtenir :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
nx = 9;
ny = 9;
cx = ceil(nx/2);
cy = ceil(ny/2);
x = [1:nx];
y = [1:ny];
Ecentre = 1;
 
[X,Y]=meshgrid(x,y);
rho = sqrt(X.^2+Y.^2);
Y-cy
X-cx
theta = atan2(Y-cy,X-cx);
E = Ecentre.*(cos(theta).^4);
figure(),imagesc(E)
figure(),mesh(E)

Quelqu'un aurait une idée pour générer ce dome qui part de centre de l'image?

Merci d'avance!

[Invité]

Bonjour,

Petite question mathématique: pour faire un dôme, dans quel intervalle crois-tu qu'il faille se trouver?

[jerome828]

Merci, et la lumière fut...

C'est invariant en theta mais bien en rho!
Des rappels de collège ne me feraient pas de mal je crois :$

Code corrigé :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
nx = 99;
ny = 99;
cx = ceil(nx/2);
cy = ceil(ny/2);
x = [1:nx];
y = [1:ny];
Ecentre = 1;
 
[X,Y]=meshgrid(x,y);
rho = sqrt((X-cx).^2+(Y-cy).^2);
E = Ecentre.*(cos(2*pi.*rho./360).^4);
figure(),imagesc(E)
figure(),mesh(E)

[Jerome Briot]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
x = [1:nx];
y = [1:ny];

Les crochet sont inutiles

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
E = Ecentre.*(cos(2*pi.*rho./360).^4);

Les points en rouge sont inutiles