matlab中滤波反投影/ iradon的定性和定量分析

时间:2015-06-08 06:29:01

标签: matlab image-processing tomography-reconstruction quantitative back-projection

我想知道是否有人遇到过这个问题。

我可以从matlab重建类似于原始图像的图像,但实际值总是不同。

例如,原始图像的矩阵值在0到1之间,而我的重建图像的范围从-0.2到0.4。

重建的图像看起来与原始图像类似,只是图像中的数据具有不同的比例。

这是我的意思的示例代码。

           <?php if (count($room['room_images']) > 0 ) { ?>  
        <p><a class="fancybox" rel="gallery<?php echo $counter; ?>" href="<?php echo $image['image']['sizes']['large']; ?>" title="<?php echo $image['caption']; ?>" >Click here</a> to view more room images</p>
        <div class="image">
          <?php // $isFirst = true;
                            foreach ($room['room_images'] as $image ) {
                                // if (!$isFirst) { ?>
          <div class="foto"> <a class="fancybox" rel="gallery<?php echo $counter; ?>" href="<?php echo $image['image']['sizes']['large']; ?>" title="<?php echo $image['caption']; ?>" ><img style="width:150px; height:96px;" src="<?php echo $image['image']['sizes']['thumbnail']; ?>" alt="" /></a> </div>

1 个答案:

答案 0 :(得分:0)

这些结果并不是我发现的结果。

>> min(min(ir))

-0.0583

>> max(max(ir))

0.9658

请记住,反向雷达变换只能逼近原始图像的重建。只有180个视图,肯定会有一些差异。

Radon变换固有地导致一些信息丢失,因为必须将像素投影到新的坐标系上并重新分组 - 在投影和反投影期间。这导致重建的图像稍微降低。 Radon变换不像傅立叶变换那样是可逆的。

为了获得更好的效果,请尝试使用更大的图像尺寸和更多的视角。

p=phantom(256);
theta=0:0.01:179;

并尝试使用不同的滤波器(F.B.P.中的F),例如Shepp-Logan,它可以降低高频并减少过冲。

ir=iradon(r,theta,'linear','Shepp-Logan');