CT图像重建系统

CT图像重建系统（源码） 图像重建是图像处理的一个重要的研究分支，其重要意义是在于获取被检测物体内部结构的图像而不是对物体造成任何物理上的损伤。

#include “vtkRenderer.h“
#include “vtkRenderWindow.h“
#include “vtkRenderWindowInteractor.h“
#include “vtkVolume16Reader.h“
#include “vtkPolyDataMapper.h“
#include “vtkActor.h“
#include “vtkProperty.h“
#include “vtkOutlineFilter.h“
#include “vtkCamera.h“
#include “vtkPolyDataMapper.h“
#include “vtkStripper.h“
#include “vtkPolyDataNormals.h“
#include “vtkMarchingCubes.h“

//-------------------------------------------------------
// 主函数
//-------------------------------------------------------
int main （int argc char **argv）
{

    // 参数个数检查
    if（argc < 2）
    {
        cout << “用法： “ << argv[0] << “ DATA/headsq/quarter“ << endl;
        return 1;
    }

    // ---------读取数据---------

    // 创建2D数据读者对象
    vtkVolume16Reader *v16 = vtkVolume16Reader::New（）;
    // 设置数据的各维大小
    v16->SetDataDimensions（6464）;
    // 设置数据的字节排列方式
    v16->SetDataByteOrderToLittleEndian（）;
    // 设置文件前辍，它将结合下一参数，以“FilePrefix.%d”方式读取随后文件
    v16->SetFilePrefix （argv[1]）;
    // 设置文件读取范围
    v16->SetImageRange（1 93）;
    // 设置CT数据的切片间距和间隔
    v16->SetDataSpacing （3.2 3.2 1.5）;  
    
    // ---------重建皮肤---------

    // 创建一个Marching Cubes算法的对象
    vtkMarchingCubes *skinExtractor = vtkMarchingCubes::New（）;
    // 获得所读取的CT数据
    skinExtractor->SetInputConnection（v16->GetOutputPort（））;
    // 提取出CT值为500的皮肤数据
    skinExtractor->SetValue（0 500）;
    // 重新计算法向量
    vtkPolyDataNormals *skinNormals = vtkPolyDataNormals::New（）;
    skinNormals->SetInputConnection（skinExtractor->GetOutputPort（））;
    skinNormals->SetFeatureAngle（60.0）;
    // 创建三角带对象
    vtkStripper *skinStripper = vtkStripper::New（）;
    // 将生成的三角片连接成三角带
    skinStripper->SetInputConnection（skinNormals->GetOutputPort（））;
    // 创建一个数据映射对象
    vtkPolyDataMapper *skinMapper = vtkPolyDataMapper::New（）;
    // 将三角带映射为几何数据
    skinMapper->SetInputConnection（skinStripper->GetOutputPort（））; 
    skinMapper->ScalarVisibilityOff（）;
    // 创建一个代表皮肤的演员对象
    vtkActor *skin = vtkActor::New（）;
    // 获得皮肤几何数据的属性
    skin->SetMapper（skinMapper）;
    // 设置皮肤颜色的属性
    skin->GetProperty（）->SetDiffuseColor（1 .49 .25）;
    // 设置反射率
    skin->GetProperty（）->SetSpecular（.3）;
    // 设置反射光强度 
    skin->GetProperty（）->SetSpecularPower（20）;
    // 设置不透明度
    skin->GetProperty（）->SetOpacity（1.0）;

    // ---------重建骨骼---------

    // 创建一个Marching Cubes算法的对象
    vtkMarchingCubes *boneExtractor = vtkMarchingCubes::New（）;
    // 获得所读取的CT数据
    boneExtractor->SetInputConnection（v16->GetOutputPort（））;
    // 提取出CT值为1150的骨骼数据
    boneExtractor->SetValue（0 1150）;
    // 重新计算法向量
    vtkPolyDataNormals *boneNormals = vtkPolyDataNormals::New（）;
    boneNormals->SetInputConnection（boneExtractor->GetOutputPort（））;
    boneNormals->SetFeatureAngle（60.0）;
    // 创建三角带对象
    vtkStripper *boneStripper = vtkStripper::New（）;
    // 将生成的三角片连接成三角带
    boneStripper->SetInputConnection（boneNormals->GetOutputPort（））;
    // 创建一个数据映射对象
    vtkPolyDataMapper *b

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-08-10 09:27  CT图像重建系统\
     目录           0  2017-08-10 09:27  CT图像重建系统\Data\
     目录           0  2017-08-10 09:27  CT图像重建系统\Data\headsq\
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.1
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.10
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.11
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.12
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.13
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.14
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.15
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.16
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.17
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.18
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.19
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.2
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.20
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.21
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.22
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.23
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.24
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.25
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.26
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.27
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.28
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.29
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.3
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.30
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.31
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.32
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.33
     文件        8192  2011-08-24 09:49  CT图像重建系统\Data\headsq\quarter.34
............此处省略117个文件信息