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!