ピクセルごとに非線形最小二乗法を行い、その画像を出力する方法についてご教示ください

7 次查看（过去 30 天）

KT 2015-11-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/254483-

评论： Tohru Kikawada 2017-1-20

各画像のピクセルごとに最小二乗法を行い、最終的にその最適化された値を画像として読み込みたいです。モデル関数は既知で、tとyも既知です。64*64などのmatrixになると関数をどのように設定すればいいか分からず、また、エラーが起こりピクセルごとの計算がうまくいきません。

-------------------------------------------------------------------------------------------------------------------- x1 = double(imread('1.png')); x2 = double(imread('2.png')); x3 = double(imread('3.png'));

t = [zeros(size(x1)) repmat(20,size(x1)) repmat(30,size(x1))];

y = [ones(size(x1)) x2./x1 x3./x1];

f = (exp(-x(1)*t) + exp(-(x(1)+x(2))*t)) - y;

%モデル関数とした場合に、

プログラムで利用できる最大変数サイズを超えました。

エラー lsqncommon (line 139) Jstr = sparse(ones(Jrows,Jcols));

エラー lsqnonlin (line 237) [xCurrent,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...

エラー lm (line 31) [x,resnorm,residual,exitflag,output] = lsqnonlin(@myfun1,x0); % Invoke optimizer

------------------------------------------------------------------------------------------------------------------ 「各ピクセルに対するlsqnonlin関数の使い方」と「各ピクセルごとに連続的に計算を行う」という点が分かっていません。何卒ご教示ください。

3 个评论
显示 1更早的评论隐藏 1更早的评论

Steven Lord 2015-11-12

在 MATLAB Online 中打开

After translation via Google Translate and code formatting:

To perform the least square method for each pixel of each image, you want to read finally its optimized value as an image. Model function is known, t and y are also known. When it comes to the matrix, such as 64 * 64 do not know how should I how to set the function, also, does not work is the calculation of each error occurs pixels.

-------------------------------------------------- -------------------------------------------------- ----------------

 x1 = double (imread ('1.png'));
 x2 = double (imread ('2.png'));
 x3 = double (imread ( '3.png'));
 t = [zeros (size (x1)) repmat (20, size (x1)) repmat (30, size (x1))];
 y = [.. ones (size (x1)) x2 / x1 x3 / x1];
 f = (exp (-x (1) * t) + exp (- (x (1) + x (2)) * t)) - y;

% In the case of the model function,

 You have exceeded the maximum variable size that are available in the program.
 Error lsqncommon (line 139)
 Jstr = sparse (ones (Jrows, Jcols));
 Error lsqnonlin (line 237)
 [xCurrent, Resnorm, FVAL, EXITFLAG, OUTPUT, LAMBDA, JACOB] = ...
 Error lm (line 31)
 [x, resnorm, residual, exitflag, output] = lsqnonlin (@ myfun1, x0);% Invoke optimizer