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;}