《Radon变换自定义函数实现》由会员分享，可在线阅读，更多相关《Radon变换自定义函数实现（3页珍藏版）》请在金锄头文库上搜索。
1、% Share Information % Tag: Matlab, Medical Imaging, DICOM % 标题：Radon 变换自定义函数实现% 内容：自定义函数进行雷登变换% 作者：Destiny% 日期：2012-10-05% 是否转载：否%自定义函数实现 Radon 变换%Use dicominfo() to get the metadatafilename = dir(*.dcm);r c=size(filename);a=78; %选取第 n 层则令 a = nstr = sprintf(filename(a).name);info = dicominfo(str);P
2、 = dicomread(info);figure,imagesc(P);colormap(hot);colorbar;title(原图像 );m,n = size(P);%先计算图像最大尺寸 mmax*nmax thm=45;%45 度时会出现最大尺寸pre = imrotate(P,thm);mmax,nmax = size(pre);s = 1; %定义步长%创建一个 180*nmax 的空白图片，用以存储投影后的线状图片Final = zeros(180/s,nmax);%这里 180 代表 180 角度，每个角度投影成为一条线t = 1;for theta = 1:s:179%对原图
3、旋转一个角度，求和（线积分）Protate = imrotate(P,theta);Pf = sum(Protate,1);%figure,imagesc(Pf);%colormap(hot);%colorbar;mreal,nreal=size(Pf);%计算实际尺寸%确定起始点if (nmax - nreal)/2-floor(nmax - nreal)/2) = 0From = floor(nmax - nreal)/2 + 1);%总点数为偶数时elseFrom = floor(nmax - nreal)/2) + 1;%总点数为奇数时end%确定结束点End = floor(nmax
4、-nreal)/2) + nreal;%将一个角度 Radon 变换后的线状图存入结果图像的某一行Final(180/s-t,From:End) = Pf; %从最底下一行开始存起%上移一行t = t + 1;end%再逆时针旋转以后与系统函数实现的 Radon 变换相对应Final = imrotate(Final,90);figure,imagesc(Final);colormap(hot);colorbar;title(自定义函数实现 Radon 空间信息采集);%对自定义函数实现的 Radon 空间图像进行逆变换测试theta = 0:179;iFinal = iradon(Final
5、,theta);figure,imagesc(iFinal);colormap(hot);colorbar;title(对自定义函数实现的 Radon 空间图像进行逆变换测试);%利用系统函数实现 Radon 和 Radon 逆变换theta = 1:179;pt = radon(P,theta);figure,imagesc(pt);colormap(hot);colorbar;title(系统函数实现 Radon 空间信息采集); ipt,xp = iradon(pt,theta);figure,imagesc(ipt);colormap(hot);colorbar;title(对系统函数实现的 Radon 空间图像进行逆变换测试);% % Contact: % Share and Access - Destiny %
