例9-18 .m

来自「对初学matlab的有很大帮助 认真学应该能学好」· M 代码 · 共 13 行

％首先利用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