![](https://img-blog.csdnimg.cn/20201014180756913.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
ct重建
文章平均质量分 53
猪猪虾的业余生活
这个作者很懒,什么都没留下…
展开
-
FFT相关学习记录
https://zhuanlan.zhihu.com/p/415814805原创 2023-04-10 19:47:04 · 457 阅读 · 0 评论 -
CT投影矩阵计算问题汇总
1.获得体膜二维投影图像之后,需要对图像进行分割,提取出小球的二维坐标,可以使用函数matlab 函数regionprops(img,‘指令’);对已经分割的小球进行处理,指令 依次为’centroid’, ‘Area’, 'Perimeter’可以分别获得小球的依次编码的质心,小球面积,以及各个小球的周长。但是需要注意的是,使用这些指令之前,需要将分割好的图像转换成二值图,并调用[label_img,num] = bwlabel(img); 图像对二值图进行编码,之后再调用上述的那些指令求相应的参数,否原创 2021-10-02 16:01:55 · 621 阅读 · 0 评论 -
FBP 滤波反投影重建
一.为什么光有反投影不足以得到质量较好的图像,下面这个图像能很好的说明问题,图像的边缘模糊二…中心切片定理二维图像的中心切片定理指出:二维函数 f(x, y) 的投影 p(s) 之傅里叶变换 P(ω) 等于函数 f(x, y) 的傅里叶变换 F(ωx, ωy) 沿与探测器平行的方向过原点的片段1.一维连续傅里叶变换对:二维连续傅里叶变换对:2.根据中心切片的定义,如果探测器绕物体旋转至少180°,物体的二维傅里叶变换 ) , F (ωx, ωy) 所对应于探测器方向的中心片段就能覆盖整原创 2021-08-19 09:18:43 · 3692 阅读 · 1 评论 -
matlab2019 C混编经验总结
常见报错:1.数组超过预设的最大数解决方案 :加-compatibleArrayDims’ ,和 ’没有空格mex filename.c ‘-compatibleArrayDims’原创 2021-07-29 19:36:02 · 681 阅读 · 2 评论 -
CT二维平行光束重建
clc,clear;%% 创建图像及投影M=512; %模型图像的长宽 %建议先和接收传感器的个数一样,这里不一样P = phantom(M); % 创建一个 shepp-logan 模型,原图像theta=0:0.5:179.999; %角度% theta=1:1:180; %角度[R,xp] = radon(P,(theta*180)/pi);%利用radon变换获得不同方向上的投影xp_offset = abs(min(xp) + 1); %探测器的中心% plot(xp,R(:,5原创 2021-03-10 11:46:13 · 815 阅读 · 2 评论