Radon Transform in x-y

Radon Transform in x-y, Matlab function

function [rx,xa,ry,ya]=Radonxy(original_image)

% Dr. Miguel Mora

% Radonxy ver 1.0

% Función que realiza la transformada de Radón manualmente en direcciones

% ortogonales x-y, esto con la finalidad de entregar dos vectores con

% tamaño de las dimensiones de la imagen rx y ry, así como otros dos

% vectores que representan sus repectivos ejes xa e ya, o x-axis e y-axis.

%

% Inicio: 14-jul-22

% Fin:    14-jul-22

% Validando matriz de entrada 

[Nx,Ny,color]=size(original_image);

if (color==1)

    iNG=double(original_image);

    fprintf('\nImagen de entrada en niveles de gris: Nx=%d y Ny=%d\n',Nx,Ny);

else

    fprintf('\nImagen de entrada a color: Nx=%d, Ny=%d y color=%d\n',Nx,Ny,color);

    iNG=rgb2gray(double(original_image));

end 

% Transformada de Radón

rx(Nx)=zeros;

ry(Ny)=zeros;

for i=1:Nx

    xa(i)=i;

    for j=1:Ny

        rx(i)=rx(i)+iNG(i,j);

    end

end

for j=1:Ny

    ya(j)=j;

    for i=1:Nx

        ry(j)=ry(j)+iNG(i,j);

    end

end 

end