资源简介
环形伪影是ct图像常见的伪影之一，本文的算法用sobel算子，插值处理去除伪影
代码片段和文件信息
p=phantom（256）;
imtool（p[]）;
[m n]=size（p）;
x0=m/2;
y0=n/2;
radius=sqrt（x0^2+y0^2）;
for x=1:m
for y=1:n
if （x-x0）^2+（y-y0）^2>=2025 && （x-x0）^2+（y-y0）^2<=2500
p（xy）=0.5;
end
end
end
imtool（p[]）; %加环形伪影
BW1=edge（p‘sobel‘0.01）; %sobel图像边缘提取
imtool（BW1[]）;
tt=zeros（360181）;
for theta=1:360
for r=1:181
i=x0+round（r*cos（theta/180*pi））;
j=y0+round（r*sin（theta/180*pi））;
if i>=1 && i<=m && j>=1 && j<=n
tt（thetar）=BW1（ij）;
end
en- 上一篇：直接序列扩频通信系统Matlab代码仿真
- 下一篇：自动洗衣机multisim
川公网安备 51152502000135号
评论
共有 条评论