Create The Image Laplacian Matrix Effectively

Question

Royi Avital 2014-6-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/132236-create-the-image-laplacian-matrix-effectively

评论： Royi Avital 2014-6-8

采纳的回答： Matt J

Hello,

I want to build the Spatial Laplacian of a given operation on an image.

The Matrix is given by:

The matrix Dx / Dy is the forward difference operator -> Hence its transpose is the backward difference operator.

The matrix Ax / Ay is diagonal matrix with weights which are function of the gradient of the image.

It is defined by:

Where Ix(i) is the horizontal gradient of the input image at the i-th pixel.

As said above Ax(i, j) = 0, i ~= j.

It is the same for Ay with the direction modification.

Assuming input Image G -> g = vec(G) = G(:).

I want to find and image U -> u = vec(U) = U(:) s.t.:

How can I solve it most efficiently in MATLAB?

How should I build the sparse Matrices?

Thank You.

2 个评论
显示无隐藏无

Matt J 2014-6-4

It looks like deconvreg in the Image Processing Toolbox does the above (or something similar), but without linear algebraic methods, probably.

Royi Avital 2014-6-8

I don't think deconvreg is related.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Matt J 2014-6-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/132236-create-the-image-laplacian-matrix-effectively#answer_139300

在 MATLAB Online 中打开

 [M,N]=size(inputImage);
 g=inputImage(:);
 Dx=diff(speye(N),1,1);
  Dx=kron(Dx,speye(M));
 Dy=diff(speye(M),1,1);
  Dy=kron(speye(N),Dy);
    sp=@(V) spdiags(V(:),0,numel(V),numel(V));
 Ax=sp(Dx*g);
 Ay=sp(Dy*g);
 Lg=Dx.'*Ax*Dx + Dy.'*Ay*Dy;
 u=(speye(size(Lg))+lambda*Lg)\g;

2 个评论
显示无隐藏无

Royi Avital 2014-6-3

编辑：Royi Avital 2014-6-3

Hi Mat,

Thank you for your answer.

Few remarks:

Wouldn't be faster to use 'diff' and then put it into the sparse matrix instead of doing it by matrix multiplication?
Let's say the weights are given by exponent weight (See my update to the question). How would you do it?