exa8_18.asv

来自「有关matlab r2007的实用教程」· ASV 代码 · 共 14 行

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