Matlab - Patron de moteado

X = [4,3]';

Y = [2,1]';

centers = [X,Y];

resta = (X-Y).^2;

d = sqrt(resta);

data = ones(300,300);

center = [130, 140; 220, 230; 30, 220];

radio = 15;

for i = 1:300

    for j = 1:300

        if any(sqrt(sum((center - [i, j]).^2, 2)) < radio)

            data(i, j) = 0;

        end

    end

end

imshow(data);

sqrt(sum((puntos - centro).^2, 2))

% Patron de moteado

%La idea es sacar una posicion aleatoria para el centro del circulo (la primera la pones directamente)

%Luego sacas nuevas posiciones aleatorias y compruebas la distancia entre

%centros, de la nueva posicion y la que tienes puesta, si la distancia es

%mayor que el diametro el nuevo circulo vale, si no, lo descartas y creas

%otro nuevo, asi hasta tener suficientes

%Resolucion de la camara 4608x3072 pixels

P = 4608;

Q = 3072;

%Tamaño folio A4 210x297 mm

b = 210;

h = 297;

%Tamaño del radio sabiendo que cada circulo debe ocupar 2 pixels

r = (b/P)*(sqrt(2/pi));

%Como el radio es menor a 1mm multiplico radio por 100 para poder trabajar

%con numeros sin decimales, la matriz de 1 se amplia tambien por 10.

radio = round(100*r);

%Diametro

D = radio*2;

%Numero optimo de motas para garantizar un 50:50 de ratio blanco negro

motas = round((b*h)/(2*D));

%Creacion de matriz tamaño folio de 1 que serán el fondo blanco, cuando se

%cambien 1 por 0 se pondra de negro

data = ones(100*b,100*h);

%Crear buble hasta motas para crear el Nº de centros requerido, haciendo

%que la distancia entre centros sea mayor al diametro

for k = 1:motas

        X(k) = 100*b*rand;

        Y(k) = 100*h*rand;

        center(k,k) = [X(k),Y(k)];

    for i = 1:100*b

     for j = 1:100*h

        %La distancia entre dos puntos se calcula como el calculo de la

        %hipotenusa de un triangulo formado entre 2 puntos.

        distancia = sqrt(((X(k)-X(k+1))^2)+((Y(k)-Y(k+1))^2),2);

        %para pintar un circulo negro se cambian 1 por 0 en la funcion

        %imshow

        pintar = (sqrt(sum((center - [i, j]).^2, 2));

        %Cambiar solo los 1 a 0 que sean menores a una distancia del radio

        %y aquellos cuya distancia a oreos centros sea mayor al diametro

        if any(pintar < radio) and (distancia > D)

            data(i,j) = 0;

        end

     end

    end

end

imshow(data);