zoukankan      html  css  js  c++  java
  • VTK 图像基本操作_图像类型转换

    1.vtkImageCast

    图像数据类型转换在数字图像处理中会被频繁地用到。一些常用到的图像算子(例如梯度算子)在计算时出于精度的考虑,会将结果存储为float或者double类型。但是在图像显示时,一般需要图像为unsigned char类型。

    VTK中最常用到的是VTKImageCast类型。其使用方法如下:

    1 vtkSmartPointer<vtkImageCast> imgCast =
    2 vtkSmartPointer<vtkImageCast>::New();
    3 imgCast->SetInput( (vtkDataObject*)reader->GetOutput() );
    4 imgCast->SetOutputScalarTypeToFloat();
    5 imgCast->Update();

    当我们需要把图像转换成不同的类型进行计算时,只需要把SetOutputScalarTypeToXXX()设置成相应的输出类型即可。

    需要注意的是,这个类中还有一个变量是ClampQverflow,用来标识是否组要截断数据。默认情况下,该变量值为0;当设置为1时,输出的像素值不能超过输出类型的最大值。超过时会自动截断到最大值。

    此外还需要知道的就是,这个类进行类型转换时,只是将类型进行强制转换,而没有进行比例的缩放,因此使用起来会比较受限,因此VTK中时不推荐使用这个类的,例如一副double类型的图像,图像值的范围[-1,1]。如果需要将图像转化为unsigned char类型,这种方法是不可行的!

    2.vtkImageShiftScale

    这个类时用起来要比上一类更加的灵活。他可以指定偏移和比例参数来对输入图像数据进行操作,例如一副double类型的图像,其数值范围为[-1,1],如果将其转换为unsigned char类型,需要设置shift=+1,Scale=127.5;那么输入图像的数据-1可以被映射为(-1+1)*127.5=0;+1可以被映射为(1+1)*127.5=255。具体使用方法如下所示:

    1 vtkSmartPointer<vtkImageShiftScale> ShiftScale =
    2 vtkSmartPointer<vtkImageShiftScale>::New();
    3 ShiftScale->SetInputConnection(img->GetOProducerPort());
    4 ShiftScale->SetOutputScalarTypeToUnsignedChar();
    5 ShiftScale->SetShift(1);
    6 ShiftScale->SetScale(127.5);
    7 ShiftScale->Update();

    3.实验源码

     1 #include <vtkAutoInit.h>
     2 VTK_MODULE_INIT(vtkRenderingOpenGL);
     3  
     4 #include <vtkSmartPointer.h>
     5 #include <vtkImageCast.h>
     6 #include <vtkImageShiftScale.h>
     7 #include <vtkBMPReader.h>
     8 #include <vtkMetaImageReader.h>
     9 #include <vtkImageViewer2.h>
    10 #include <vtkRenderer.h>
    11 #include <vtkRenderWindow.h>
    12 #include <vtkRenderWindowInteractor.h>
    13  
    14 int main()
    15 {
    16     vtkSmartPointer<vtkMetaImageReader> reader =
    17         vtkSmartPointer<vtkMetaImageReader>::New(); //读*.mhd文件
    18     reader->SetFileName("brain.mhd");
    19     reader->Update();
    20     //第一种类型转换方案
    21     vtkSmartPointer<vtkImageCast> imgCast =
    22         vtkSmartPointer<vtkImageCast>::New();
    23     imgCast->SetInputData( (vtkDataObject *) reader->GetOutput() );
    24     imgCast->SetOutputScalarTypeToFloat();
    25     //第二种类型转换方案
    26     vtkSmartPointer<vtkImageShiftScale> imgShiftScale =
    27         vtkSmartPointer<vtkImageShiftScale>::New();
    28     imgShiftScale->SetInputData((vtkDataObject*)reader->GetOutput());
    29     imgShiftScale->SetOutputScalarTypeToFloat();
    30  
    31     vtkSmartPointer<vtkImageViewer2> imgViewer =
    32         vtkSmartPointer<vtkImageViewer2>::New();
    33     imgViewer->SetInputConnection(imgCast->GetOutputPort());
    34  
    35     vtkSmartPointer<vtkRenderWindowInteractor> rwi =
    36         vtkSmartPointer<vtkRenderWindowInteractor>::New();
    37     imgViewer->SetupInteractor(rwi);
    38  
    39     imgViewer->SetColorLevel(500);
    40     imgViewer->SetColorWindow(2000);
    41     imgViewer->SetSlice(20);
    42     imgViewer->SetSliceOrientationToXY();
    43     imgViewer->Render();
    44  
    45     /*imgViewer->GetRenderer()->ResetCamera();
    46     imgViewer->Render();*/
    47     imgViewer->GetRenderer()->SetBackground(1.0, 1.0, 1.0);
    48     imgViewer->SetSize(640, 480);
    49     imgViewer->GetRenderWindow()->SetWindowName("ImageTypeTrans");
    50  
    51     rwi->Start();
    52     return 0;
    53 }
  • 相关阅读:
    [转]Native进程的运行过程
    android ARM 汇编学习 —— hello world
    android ARM 汇编学习—— 在 android 设备上编译c/cpp代码并用objdump/readelf等工具分析
    Linux buffer/cache异同
    分布式消息队列RocketMQ与Kafka架构上的巨大差异之1 -- 为什么RocketMQ要去除ZK依赖?
    kafka对比RocketMQ(转)
    Jmeter测试webocket协议
    Linux atop监控工具部署
    MySql计算两个日期的时间差函数
    Python selenium —— 一定要会用selenium的等待,三种等待方式解读(转)
  • 原文地址:https://www.cnblogs.com/ybqjymy/p/14241242.html
Copyright © 2011-2022 走看看