0

いくつかのガウス分布を作成して画像に配置しようとしています。ガウスはランダムに作成されたパラメーター (振幅、位置、および標準偏差) を持っています。最初にパラメータをベクトルまたは行列に入れ、次に ngrid() 関数を使用してガウス分布を作成するための 2 次元空間を取得していますが、エラーが発生します (ngrid 値を使用した数学的操作は自明ではない可能性があるため...)。エラーは次のとおりです。

     ??? Error using ==> minus
Integers can only be combined
with integers of the same class,
or scalar doubles.
Error in ==> ss_gauss_fit at 23
  gauss = amp(i)*
  exp(-((x-xc).^2 +
  (y-yc).^2)./(2*std(i)));

コードはここにあります：

clear all;
image = uint8(zeros([300 300]));
imsize=size(image);

noOfGauss=10;
maxAmpGauss=160;
stdMax=15;
stdMin=3;

for i=1:noOfGauss
    posn(:,:,i)=[  uint8(imsize(1)*rand())   uint8(imsize(2)*rand())  ];
    std(i)=stdMin+uint8((stdMax-stdMin)*rand());
    amp(i)= uint8(rand()* maxAmpGauss);
end

% draw the gaussians on blank image
for i=1:noOfGauss
    [x,y] = ndgrid(1:imsize(1), 1:imsize(2));
    xc = posn(1,1,i);
    yc = posn(1,2,i);
    gauss = amp(i)* exp(-((x-xc).^2 + (y-yc).^2)./(2*std(i)));

    image = image + gauss;
end

これを修正する方法を教えてください.2dガウスをパラメータベクトルでプロットしてください...よろしくお願いします

4

1 に答える 1

1

私がよく理解していない「画像に描く」ことについての狂気は別として、グリッド上に多数の個別のガウス分布を追加しようとしていると思います。これが私があなたのコードで行ったことです。二変量ガウスは適切に正規化されておらず、以前は標準偏差ではなく分散を使用していたことに注意してください。後者を修正しました。ただし、とにかくそれぞれに振幅値を掛けているため、正規化は気にしませんでした。

clear all;
xmax = 50;
ymax = 50;

noOfGauss=10;
maxAmpGauss=160;
stdMax=10;
stdMin=3;

posn = zeros(noOfGauss, 2);
std = zeros(noOfGauss, 1);
amp = zeros(noOfGauss, 1);

for i=1:noOfGauss
    posn(i,:)=[ xmax*rand()  ymax*rand() ];
    std(i)=stdMin+(stdMax-stdMin)*rand();
    amp(i)= rand()* maxAmpGauss;
end

% draw the gaussians
[x,y] = ndgrid(1:xmax, 1:ymax);
z = zeros(xmax, ymax);

for i=1:noOfGauss    
    xc = posn(i,1);
    yc = posn(i,2);    
    z = z + amp(i)* exp(-((x-xc).^2 + (y-yc).^2)./(2*std(i)^2));
end 

surf(x, y, z);

ランダム出力:

ここに画像の説明を入力

于 2013-03-06T01:41:22.097 に答える