Matlab - Sum of surrounding elements

Question

I want to calculate the sum of the elements surrounding a given element in a matrix. So far, I have written these lines of code:

for i=1:m,
        rij(1:n)=0
        for j=1:n,
            alive = tijdelijk(i-1,j)+tijdelijk(i+1,j)+tijdelijk(i-1,j-1)+tijdelijk(i+1,j-1)+tijdelijk(i,j+1)+tijdelijk(i,j-1)+tijdelijk(i-1,j+1)+tijdelijk(i+1,j+1)

This results in an error because, for example, i-1 becomes zero for i=1. Anyone got an idea how to do this without getting this error?

petrichor · Answer 1 · 2013-03-12T10:21:48.380

13

You can sum the elements via filtering. conv2 can be used for this manner.

Let me give an example. I create a sample matrix

>> A = reshape(1:20, 4, 5)

A =

 1     5     9    13    17
 2     6    10    14    18
 3     7    11    15    19
 4     8    12    16    20

Then, I create a filter. The filter is like a mask where you put the center on the current cell and the locations corresponding to the 1's on the filter are summed. For eight-connected neighbor case, the filter should be as follows:

>> B = [1 1 1; 1 0 1; 1 1 1]

B =

 1     1     1
 1     0     1
 1     1     1

Then, you simply convolve the matrix with this small matrix.

>> conv2(A, B, 'same')

ans =

13    28    48    68    45
22    48    80   112    78
27    56    88   120    83
18    37    57    77    50

If you want four-connected neighbors, you can make the corners of your filter 0. Similarly, you can design any filter for your purpose, such as for averaging all neighbors instead of summing them.

For details, please see the convolution article in Wikipedia.

edited Mar 12 '13 at 10:21

answered Mar 12 '13 at 10:14

petrichor

6,459
4
36
48

1

Hah, was just writing the same thing! – tdc Mar 12 '13 at 10:18
1

:) If it worked for you, you can accept the answer as a solution for your question. – petrichor Mar 12 '13 at 10:39
Was very helpful, so did that :) – user2089012 Mar 12 '13 at 10:54
I think you are not familiar with Stackoverflow. Accepting an answer is done by clicking to the tick mark on the left corner of the answer. – petrichor Mar 12 '13 at 21:08

score 1 · Answer 2 · answered Mar 12 '13 at 10:20

1

Two possibilities : change the limits of the loops to i=k:(m-k) and j=k:(n-k) or use blkproc

ex :

compute the 2-D DCT of each 8-by-8 block

I = imread('cameraman.tif');
fun = @dct2;
J = blkproc(I,[8 8],fun);
imagesc(J), colormap(hot)

answered Mar 12 '13 at 10:20

lucasg

10,734
4
35
57

*or use, if you have the Image Processing Toolbox, blkproc* methinks – High Performance Mark Mar 12 '13 at 10:29

score 0 · Answer 3 · answered Mar 12 '13 at 15:56

There are lots of things you can do at the edges. Which you do depends very specifically on your problem and is different from usage case to usage case. Typical things to do:

If (i-1) or (i+1) is out of range, then just ignore that element. This is equivalent to zero padding the matrix with zeros around the outside and adjusting the loop limits accordingly
Wrap around the edges. In other words, for an MxN matrix, if (i-1) takes you to 0 then instead of taking element (i-1, j) = (0, j) you take element (M, j).

Since your code mentions "your teacher" I'd guess that you can ask what should happen at the edges (or working it out in a sensible manner may well be part of the task!!).

So how about wrapping around the edges? Can in be done in a simple manner using either **conv2** or **blkproc**? — Pranasas, May 08 '13 at 17:03
Honestly, I'm not sure. I think conv2 zero-pads around the edges, but I couldn't say for sure. — FakeDIY, May 09 '13 at 09:37

Matlab - Sum of surrounding elements

3 Answers3

Linked