How can i get the divergence of vector field of a gray scale image??

How can i get the divergence of vector field of a gray scale image??

Answer 1

[Thorsten]

In this example the vector field is computed as the gradient in x and y direction:

 function D = imdiv(I)
Ix = diff([I I(:, end)]')';
Iy = diff([I; I(end, :)]);
[X Y] = meshgrid(1:size(I, 2), 1:size(I, 1));
D = divergence(X, Y, Ix, Iy);

Answer 2

[Sean de Wolski]

doc divergence
doc ndgrid

Build a grid of every row/col page and then call divergence with these and your displacement vector values.