Discussion :

[naomie]

Bonjour,

Je cherche a faire des fit gaussiens en 3D, est ce possible en matlab ? Je dois fitter des images comme celle en attachement et je cherche la position (x0, y0) du centre de ma gaussienne.
Je n'ai pas les coordonnees (x,y,z) de ma courbe, je pars simplement d'une matrice 2D, resultat de la commande imread('mon-image.tif')

Merci d'avance,


Naomie.

[mr_samurai]

Salut,

Essayes ceci :

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
% IM est la matrice de ton image ...
% ... un lissage peux être 
IM = rand(10, 30);
 
% Gaussienne N(M,sigma) , M(xm, ym)
% fonction objective: x [a, sigma, xm, ym]
 
[X,Y] = meshgrid(1:size(IM, 2), 1:size(IM, 1));
 
myFun = @(x)norm( IM - x(1) * ...
    exp(-x(2)*((X-x(3)).^2 + (Y-x(4)).^2)));
% L'expression matricielle avec sigma matrice peut servir ...
 
x0 =  [1, 1, size(IM,1)/2 , size(IM,2)/2]
x = fminsearch(myFun,x0)

++ bonne chance