数字图像处理图象重建.doc

数字图像处理上机实习报告 （第五专题） 学生姓名： 班 级： 学 号： 指导老师： 专题五、图象重建 题目要求 根据傅立叶反变换实施图象重建。 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 1 1 1 0 0 0 0 0 0 1 1 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 1 1 1 0 1 1 1 1 0 0 0 0 0 0 0 1 1 1 0 0 1 1 1 1 0 0 0 0 0 0 1 1 1 0 0 0 1 1 1 1 0 0 0 0 0 1 1 1 0 0 0 0 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 算法设计 原理简介 图像投影，就是说将图像在某一方向上做线性积分（或理解为累加求和）。如果将图像看成二维函数f(x, y)，则其投影就是在特定方向上的线性积分，比如f(x, y)在垂直方向上的线性积分就是其在x轴上的投影；f(x, y)在水平方向上的线积分就是其在y轴上的投影。通过这些投影，可以获取图像在指定方向上的突出特性，这在图像模式识别等处理中可能会用到。 Radon变换（拉东变换），就是将数字图像矩阵在某一指定角度射线方向上做投影变换。这就是说可以沿着任意角度theta来做Radon变换。 具体实现 首先对图像作0到180度的投影，并且角度间隔固定。然后进行Radon反变换。下图为图像重建的流程框图。 load D:\matlab2011\work\mat R = radon(mat,0:179); P = iradon(R,0:179); f = im2bw(P); subplot(1,3,1);imshow(mat);title(原图像) subplot(1,3,2);imshow(P);title(复原后的图像) subplot(1,3,3);imshow(f);title(复原后二值化的图像) figure; [R, xt] = radon(mat, [0 45 90]); % 在0、45、90度方向做radon变换 subplot(131); plot(xt, R(:, 1)); title(水平方向的radon变换曲线); subplot(132); plot(xt, R(:, 2)); title(45度方向的radon变换曲线); subplot(133); plot(xt, R(:, 3)); title(垂直方向的radon变换曲线); 4.结果分析 （1）0到180度的投影 （2）0到90度的投影 （3）0到45度的投影 （4）水平，45度和垂直方向的Randon变换曲线 由于复原后的图像的灰度值已经离散化，为了原图像进行对比，故将复原后的图像进行二值化。通过以上图像对比可知，在进行Radon变换时，采样的角度范围过小，会使图像重建时产生失真。