雷当变换的实现

数字图像处理实验报告实验名称雷当变换的实现一、实验目的1.理解Radon变换的基本原理；2.掌握显示变换图象的方法；3.了解使用MATLAB进行Radon变换的编程方法。二、实验原理图像投影，就是说将图像在某一方向上做线性积分（或理解为累加求和）。如果将图像看成二维函数f(x,y)，则其投影就是在特定方向上的线性积分，比如f(x,y)在垂直方向上的线性积分就是其在x轴上的投影；f(x,y)在水平方向上的线积分就是其在y轴上的投影。Radon变换，就是将数字图像矩阵在某一指定角度射线方向上做投影变换。这就是说可以沿着任意角度theta来做Radon变换。三、实验代码及结果1.实验代码(1)clc;clearall;closeall;I=zeros(256,256);[r,c]=size(I);I(floor(1/5*r:4/5*r),floor(3/5*c:4/5*c))=1;figure;subplot(2,2,1);imshow(I);title('原图像');[R,xt]=radon(I,[04590]);%在0、45、90度方向做radon变换subplot(2,2,2);plot(xt,R(:,1));title('水平方向的radon变换曲线');subplot(2,2,3);plot(xt,R(:,2));title('45度方向的radon变换曲线');subplot(2,2,4);plot(xt,R(:,3));title('垂直方向的radon变换曲线');(2)I=zeros(100,100);I(25:75,25:75)=1;imshow(I)[R,xp]=radon(I,[045]);figure;plot(xp,R(:,1));title('R_{0^o}(x\prime)')theta=0:180;[R,xp]=radon(I,theta);imagesc(theta,xp,R);title('R_{\theta}(X\prime)');xlabel('\theta(degrees)');ylabel('X\prime');set(gca,'XTick',0:20:180);colormap(hot);colorbar2.实验结果(1)(2)结果原图四、总结将图像变换到按角度投影区域，可以应用于检测直线；通过将图像矩阵在多角度做积分投影，再对得到的数据做统计分析，可以确定出图像的一些基本性质。