Question

Matlab task：Realizing radon transform without using built-in matlab fuction

Answer 1

`Hey,

Note: If you have any queries related the answer please do comment. I would be very happy to resolve all your queries.

clc%clears screen
clear all%clears history
close all%closes all files
format long
I = zeros(100,100);
I(1:25,1:25) = 1;
figure;imshow(I)
[R] = own_radon(I, 90);
figure;plot(R)

function processed_img = own_radon(img, theta)

[n_row, n_column] = size(img);
Dia = sqrt(n_row + n_column);
row_Pad = ceil(Dia - n_row) + 2;
col_Pad = ceil(Dia - n_column) + 2;
pad = zeros(n_row+row_Pad, n_column+col_Pad);
pad(ceil(abs(row_Pad/2)):(ceil(abs(row_Pad/2))+abs(n_row)-1), ...   
ceil(abs(col_Pad/2)):(ceil(abs(col_Pad/2))+abs(n_column)-1)) = img;
angle = length(theta);
processed_img = zeros(size(pad,2), angle);
for i = 1:angle
tic
tmpimg = imrotate(pad, 90-theta(i), 'bilinear', 'crop');
processed_img(:,i) = (sum(tmpimg))';
theta(i)
toc
end
end

Kindly revert for any queries

Thanks.