资源描述：
《mr原理、技术与设备（冯衍秋）《图像重建实验》实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、《图像重建实验》实验报告姓名黄仪强学号3128031076专业年级12级医工影像实验题目图像重建实验实验目的熟练图像重建，加深对课程学习的了解熟悉傅里叶变换的过程和强化matlab的编程能力实验内容图像重建作业headraw.mat1.傅里叶逆变换重建图像，2.将数据第1到64行数据置0,然后重建3.将数据1到64行与193到256行置0,…4.选择数据奇数行，…5•选择数据偶数行，…6.将数据放到512*512大小的矩阵中心，其他位置处为0,…实验步骤实验一：(傅里叶逆变换重建图像)>>load(：，F:医学影橡设-§
2、-Ji!Rheadraw.mat?)I=headraw;subplot(b2,1)jimshow(I);title(.^图经过二维傅里叶变化后的图像');C=fftshift(ifft2(I));subplot(1,292)5imshow(C);title(?重建后的图像');实验二：(将数据第1到64行数据置0,然后重建)>>loadCF:医学影像设备MRheadraw.mat')I=headraw;1(1:64,:)=0;subplot(1,2,1),imshow(I);title('将数据第1到64行数据
3、置0的二维FFT);B=fftshift(ifft2(1));subplot(I,2,2),imshow(B);title('重建后');实验三：(将数据1到64行与193到256行置0,…)>>load('F:医学影像设备MRheadraw.mat')I=headraw;1(1:64,:)=0;1(193:256,:)=0;subplot(1,2,1),imshow(I);titleC将数据1到64行与193到256行置0的二维FFT');B=fftshift(ifft2(D);Subplot(1,292)3ims
4、how(B);titleC重建后');nr•iii•・实验四：(选择数据奇数行，…)»loadCF:医学戢像设备MRheadraw.mat')I=headraw;1(1:2:endj:)=0;subplot(132,1),imshow(I);titled选择数据奇数行的二维FFT');B=fftshift(ifft2(I));subplot(1,2,2)jimshow(B);title('重建后');■■■••■■■■实验五：(选择数据偶数行，…)>>load('F:医学影像设备MRheadraw.mat')I
5、=headraw;I(2:2:endj:)=0;subplot(b2,1),imshow(I);title选择数据偶数行的二维FFT');B=fftshift(ifft2(I));subplot(1,2,2),imshow(B);title('重建后');实验六：(将数据放到512*512大小的矩阵中心，其他位置处为0,…)>>load('F:医学影像设备ffiheadraw.mat');I=headraw;D=zeros(512$512);D(129:384,129:384)=1;subplot(lj2j1),im
6、show(D);title('将数据放到512*512大小的矩阵中心,其他位畫处为0的二目B=fftshift(ifft2(D));subplot(1,2,2)imshow(B);title('重建后');附录原图经过二维傅里叶变化后的图像重建后的图像••3验实将数据第1到64行数据置0的二维FFT-…：乡>-/将数据倒64行与193到256行置0的二维FFT重建后实验4:选择数据奇数行的二维FFT重建后实验5:选择数据偶数行的二维FFT重建启实验6:将数据放到512*512t小的矩阵中心，其他位置处为0的二雉FFT重建后
7、实验分析实验的重点是要搞淸楚傅里叶变换和反变换的原理，还有也要搞清楚重建的步骤，另外熟悉matIab的相关函数也很重要。在二维傅里叶变换和频谱转移再傅里叶变换的过程中要细心，步骤不能乱。另外由于matlab的知识忘了很多，很多函数都不记得，要翻书本，花费不少时间，做得头都大了，发现真是书到用时方恨少。重点掌握好matIab的知识也是很重要的。第二次作业1.模拟Bloch方程的受激(Excitation)过程。wO=l;wl=3;t=0:0•001:30;Mz0=100;RMz=MzO・*cos(wl*t/180*pi);M
8、z=RMz;RMxy=MzO.*sin(wl*t/180*pi);Mxy=RMxy•*exp(-i*w0*t);Mx=real(Mxy);My=imag(Mxy);plot3(Mx,My,Mz)1002模拟Bloch方程的弛豫(Relaxation)过程。t=l:0.01:2000;Tl=15;T2=1
