VTK学习笔记-2-TIFF图像数据的重切片

玩了半天，终于可以实现对于TIFF图像的三维重切片。首先需要注意的是vtkTIFFReader并不是支持所有的TIFF压缩格式，比如LZ压缩的就不支持。个人觉得，VTK在DICOM格式上具有更强的操作性，貌似其存在便是针对了DICOM格式而进行的，至于其他格式的数据，有点是后期为了迎合兼容性而做出的各种的添加。再者，VTK的资料是有限有限的，除了那些Examples。。。那个用户手册是要钱的，79美金。。。。。。

#include "vtkImageReader2.h"
           
    
#include "vtkImageReslice.h"
#include "vtkWindowLevelLookupTable.h"
#include "vtkImageMapToColors.h"
#include "vtkImageActor.h"
#include "vtkRenderer.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkImageData.h"
#include<vtkTIFFReader.h>

int main ()
{	
  vtkTIFFReader *reader= vtkTIFFReader::New();//用tiff进行读取
   reader->SetFilePrefix("H:\\Crop_test4_4X4\\");//尝试更大的数据
  reader->SetFilePattern("%saaaa_%05d%_bbb.tif");//这个具体的图片的命名，就要看自己的实际操作了
  reader->SetDataExtent(0, 63, 0, 63, 188, 230);
  reader->SetDataSpacing(3.2, 3.2, 1.5);
  reader->SetDataOrigin(0.0, 0.0, 0.0);
  reader->SetDataScalarTypeToUnsignedShort();
  reader->UpdateWholeExtent();

  // 更新并获取图像参数。
  reader->GetOutput()->UpdateInformation();
  int extent[6];
  double spacing[3];
  double origin[3];
  reader->GetOutput()->GetWholeExtent(extent);
  reader->GetOutput()->GetSpacing(spacing);
  reader->GetOutput()->GetOrigin(origin);

  // 计算中心位置。
  double center[3];
  center[0] = origin[0] + spacing[0] * 0.5 * (extent[0] + extent[1]); 
  center[1] = origin[1] + spacing[1] * 0.5 * (extent[2] + extent[3]); 
  center[2] = origin[2] + spacing[2] * 0.5 * (extent[4] + extent[5]); 

   //使用下面三个轴重切后，ZSlice将对应原图像横断面（默认）
   /*static double axialX[3]={ 1, 0, 0};
   static double axialY[3]={ 0, 1, 0};
   static double axialZ[3]={ 0, 0, 1};*/

  // ZSlice对应原图像冠状面（即重切的z轴方向为原y轴方向）
  // static double coronalX[3]={ 1, 0, 0};
  // static double coronalY[3]={ 0, 0,-1};
  // static double coronalZ[3]={ 0, 1, 0};

  // ZSlice对应原图像矢状面（即重切的z轴方向为原x轴方向）
  static double sagittalX[3]={ 0, 1, 0};
  static double sagittalY[3]={ 0, 0,-1};
  static double sagittalZ[3]={-1, 0, 0};

  // ZSlice对应原图像斜面
   static double obliqueX[3]={ 1, 0, 0};
   static double obliqueY[3]={ 0, 0.866025,0.5};
   static double obliqueZ[3]={ 0,-0.5, 0.866025};

  // 按指定的方向抽取片。
  vtkImageReslice *reslice = vtkImageReslice::New();
  reslice->SetInputConnection(reader->GetOutputPort());
  reslice->SetOutputDimensionality(2);
  //reslice->SetResliceAxesDirectionCosines(sagittalX, sagittalY, sagittalZ);
  //reslice->SetResliceAxesDirectionCosines(axialX, axialY, axialZ);
  reslice->SetResliceAxesDirectionCosines(obliqueX, obliqueY, obliqueZ);
  reslice->SetResliceAxesOrigin(center);
  reslice->SetInterpolationModeToLinear();

  vtkWindowLevelLookupTable *table=vtkWindowLevelLookupTable::New();
  /*table->SetWindow(2000);
  table->SetLevel(1000);*/
  table->SetWindow(200);
  table->SetLevel(100);
  vtkImageMapToColors *color = vtkImageMapToColors::New();
  color->SetLookupTable(table);
  color->SetInputConnection(reslice->GetOutputPort());

  vtkImageActor *actor = vtkImageActor::New();
  actor->SetInput(color->GetOutput());

  vtkRenderer *renderer = vtkRenderer::New();
  renderer->AddActor(actor);

  vtkRenderWindow *window = vtkRenderWindow::New();
  window->AddRenderer(renderer);

  vtkRenderWindowInteractor *interactor = vtkRenderWindowInteractor::New();
  window->SetInteractor(interactor);
  window->Render();

  interactor->Start();

  return 0;
}