例9-18 .m

这是一个MATLAB7.0基础与提高例题的所有源码

％首先利用edge函数计算二进制图像的边界
I = fitsread('solarspectra.fts');
I = mat2gray(I);
BW = edge(I);
imshow(I), figure, imshow(BW)
％接着计算边界图像的Radon变换
theta = 0:179;
[R,xp] = radon(BW,theta);
figure, imagesc(theta, xp, R); colormap(hot);
xlabel('\theta (degrees)'); ylabel('x\prime');
title('R_{\theta} (x\prime)');
colorbar