Wie erstelle ich einen Diamanten, der mit einem in der Mitte einer Matrix gefüllt ist?

Question

Ich kenne den Code unten: mit 1. gefüllt werden Was soll ich tun

A= 
    0 0 1 0 0 
    0 1 0 1 0 
    1 0 0 0 1 
    0 1 0 1 0 
    0 0 1 0 0 

Aber ich möchte, dass die Diamant:

N = 5; 
assert(N>1 && mod(N,2)==1); 
A = zeros(N); 

% diamond mask 
N2 = fix(N/2); 
[I,J] = meshgrid(-N2:N2); 
mask = (abs(I) + abs(J)) == N2; 

% fill with zeros 
A(mask) = 1; 

die Matrix A dazu verwandelt?

Answer 1

Fügen Sie einfach ein for loop und alle Diagonalen füllen:

N = 5; 
assert(N>1 && mod(N,2)==1); 
A = zeros(N); 

% diamond mask 
N2 = fix(N/2); 
[I,J] = meshgrid(-N2:N2); 
for id = 0:N2 
    A((abs(I) + abs(J)) == id) = 1; 
end 

Answer 2

Was Sie brauchen, ist eine 1 oder eine 0 auf dem Manhattan distance von jedem Array Lage zum Zentrum Ihres Diamanten

N = 5; 
assert(N>1 && mod(N,2)==1); 
A = false(N); 

[m, n] = size(A); %dimensions of A 
X = floor([m, n]/2); %floored division gives integer indices of center of array 
x = X(1); y = X(2); 
radius = m/2;  %half the height gives the radius 
for a = 1 : m 
    for b = 1 : n 
     A(a,b) = abs(a-x)+abs(b-y) <= radius; %test if manhatten distance <= radius 
    end 
end  

Basis zurückzukehren

Dies muss natürlich bearbeitet werden, um Ihrem speziellen Fall gerecht zu werden ... Insbesondere kann die Mitte Ihres Diamanten realistisch platziert werden, indem Sie x, y ändern, und der Radius kann entweder kleiner oder größer als die halbe Breite des Arrays sein also wähle.

Answer 3

Hier ist ein vektorisiert Ansatz bsxfun -

Nh = (N+1)/2; 
range_vec = [1:Nh Nh-1:-1:1]; 
out = bsxfun(@plus,range_vec(:),range_vec) > Nh 

Beispielläufe -

1) N = 5:

out = 
    0  0  1  0  0 
    0  1  1  1  0 
    1  1  1  1  1 
    0  1  1  1  0 
    0  0  1  0  0 

2) N = 9:

out = 
    0  0  0  0  1  0  0  0  0 
    0  0  0  1  1  1  0  0  0 
    0  0  1  1  1  1  1  0  0 
    0  1  1  1  1  1  1  1  0 
    1  1  1  1  1  1  1  1  1 
    0  1  1  1  1  1  1  1  0 
    0  0  1  1  1  1  1  0  0 
    0  0  0  1  1  1  0  0  0 
    0  0  0  0  1  0  0  0  0 

Answer 4

Sie können 0 verwendenund flip Funktionen:

mat = tril(ones(N), round((N-1)/2)) - tril(ones(N), round((-N-1)/2)); 
out = mat & flip(mat) 

Odd Werte von N:

% N = 5; 

out = 

    0  0  1  0  0 
    0  1  1  1  0 
    1  1  1  1  1 
    0  1  1  1  0 
    0  0  1  0  0 

Auch Werte von N:

% N = 4; 

out = 

    0  1  1  0 
    1  1  1  1 
    1  1  1  1 
    0  1  1  0