matlab验证中心切片定理

刚看到别人用matlab的程序验证中心切片定理,但是有两处不明,然后自己做了修改

clear all 
clc
I=phantom(256);%新建一个shapp_logan模型I
figure,imshow(I);%见于图1
H=imrotate(I,270);%将图像旋转270°,相当于θ=90°,
                  %使得X线以与X轴平行的方向进行投影
figure,imshow(H);%见于图2
p=sum(H,1);%进行投影
figure,plot(1:256,p);%见1
f=fft(p);
fab=f.*conj(f);
%f3=fftshift(fab);%**修改处**
figure,plot(1:256,fab);%见2
f2d=fft2(I);%此处是直接对原始图像进行二维傅立叶变换并进行验证
fab2d=f2d.*conj(f2d);
%f4=fftshift(fab2d);%**修改处**
figure,plot(fab2d(:,1));%见3

傅里叶变换后不使用fftshift,直接输出傅里叶变换后的结果,进行比较

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值