2. 程式人生 > >VTK讀取DICOM醫學圖片進行體繪製的幾種方法

VTK讀取DICOM醫學圖片進行體繪製的幾種方法

阿新 • • 發佈：2019-01-02

注意的是：VTK中不同的vtkVolumeMapper支援不同的資料型別。比如vtkVolumeRayCastMapper和vtkVolumeTextureMapper2D只能支援單分組VTK_UNSIGNED_CHAR和VTK_UNSIGNED_SHORT型別資料，因此當讀入其他型別的影象資料時，需要對資料進行轉換，比如採用vtkImageCast或者vtkImageScale:而vtkVolumeTextMapper3D則支援任意資料型別，但是必須是單分組資料或者多元獨立資料。vtkFixedPointVolumeRayCastMapper靈活性最高，可以支援所有型別資料，最高四元資料。

使用vtkFixedPointVolumeRayCastMapper類

#include "vtkRenderer.h"  
#include "vtkRenderWindow.h"  
#include "vtkRenderWindowInteractor.h"  
#include "vtkActor.h"  
#include "vtkSmartPointer.h"  
#include "vtkProperty.h"  
#include "vtkCamera.h"  
#include "vtkDICOMImageReader.h"  
#include "vtkImageCast.h"  
#include "vtkPiecewiseFunction.h"
 
  
#include "vtkColorTransferFunction.h"  
#include "vtkVolumeProperty.h"  
#include "vtkVolumeRayCastCompositeFunction.h"  
#include "vtkVolumeRayCastMapper.h"  
#include "vtkVolume.h" 
#include "vtkAutoInit.h"
#include "vtkFixedPointVolumeRayCastMapper.h"

VTK_MODULE_INIT(vtkRenderingOpenGL);
VTK_MODULE_INIT(vtkInteractionStyle);
VTK_MODULE_INIT(vtkRenderingFreeType);
VTK_MODULE_INIT(vtkRenderingVolumeOpenGL);

int main()
{
    std::string
 
 str = "F:\\CT3";
    vtkSmartPointer<vtkDICOMImageReader> reader = vtkSmartPointer<vtkDICOMImageReader>::New();
    reader->SetDirectoryName(str.c_str());
    reader->Update();




    vtkSmartPointer<vtkPiecewiseFunction> opacityFun = vtkSmartPointer<vtkPiecewiseFunction>::New();
    opacityFun->AddPoint(120,0.0);
    opacityFun->AddPoint(250,1.0);
    opacityFun->AddPoint(520,1.0);
    opacityFun->AddPoint(650,0.0);


    vtkSmartPointer<vtkColorTransferFunction> TransferFun = vtkSmartPointer<vtkColorTransferFunction>::New();
    TransferFun->AddRGBPoint(120, 255 / 255.0, 98 / 255.0, 98 / 255.0);
    TransferFun->AddRGBPoint(250, 255 / 255.0, 255 / 255.0, 180 / 255.0);
    TransferFun->AddRGBPoint(520, 1.0, 1.0, 1.0);
    TransferFun->AddRGBPoint(650, 1.0, 1.0, 1.0);

    vtkSmartPointer<vtkPiecewiseFunction> grideFun = vtkSmartPointer<vtkPiecewiseFunction>::New();
    grideFun->AddPoint(120, 2.0);
    grideFun->AddPoint(250, 2.0);
    grideFun->AddPoint(520, 0.1);
    grideFun->AddPoint(650, 0.1);

    vtkSmartPointer<vtkVolumeProperty> volumeProperty = vtkSmartPointer<vtkVolumeProperty>::New();
    volumeProperty->SetColor(TransferFun);
    volumeProperty->SetScalarOpacity(opacityFun);
    volumeProperty->SetGradientOpacity(grideFun);
    volumeProperty->ShadeOn();

    volumeProperty->SetAmbient(0.2);
    volumeProperty->SetDiffuse(0.9);
    volumeProperty->SetSpecular(0.2);
    volumeProperty->SetSpecularPower(10);

    vtkSmartPointer<vtkVolumeRayCastCompositeFunction> rayCastFun = vtkSmartPointer<vtkVolumeRayCastCompositeFunction>::New();

    vtkSmartPointer<vtkFixedPointVolumeRayCastMapper> mapper1 = vtkSmartPointer<vtkFixedPointVolumeRayCastMapper>::New();
    mapper1->SetInputConnection(reader->GetOutputPort());
    mapper1->SetAutoAdjustSampleDistances(0);





    vtkSmartPointer<vtkVolume> volume = vtkSmartPointer<vtkVolume>::New();
    volume->SetMapper(mapper1);
    volume->SetProperty(volumeProperty);
    vtkSmartPointer<vtkRenderer> ren = vtkSmartPointer<vtkRenderer>::New();
    ren->AddVolume(volume);
    ren->SetBackground(1.0,1.0,1.0);


    vtkSmartPointer<vtkRenderWindow> renwin = vtkSmartPointer<vtkRenderWindow>::New();
    renwin->AddRenderer(ren);
    renwin->SetSize(740, 480);

    vtkSmartPointer<vtkRenderWindowInteractor> iren = vtkSmartPointer<vtkRenderWindowInteractor>::New();
    iren->SetRenderWindow(renwin);

    iren->Initialize();
    renwin->Render();
    iren->Start();
    return 0;
}

使用tkVolumeRayCastCompositeFunction類

#include "vtkRenderer.h"  
#include "vtkRenderWindow.h"  
#include "vtkRenderWindowInteractor.h"  
#include "vtkActor.h"  
#include "vtkSmartPointer.h"  
#include "vtkProperty.h"  
#include "vtkCamera.h"  
#include "vtkDICOMImageReader.h"  
#include "vtkImageCast.h"  
#include "vtkPiecewiseFunction.h"  
#include "vtkColorTransferFunction.h"  
#include "vtkVolumeProperty.h"  
#include "vtkVolumeRayCastCompositeFunction.h"  
#include "vtkVolumeRayCastMapper.h"  
#include "vtkVolume.h" 
#include "vtkAutoInit.h"

VTK_MODULE_INIT(vtkRenderingOpenGL);
VTK_MODULE_INIT(vtkInteractionStyle);
VTK_MODULE_INIT(vtkRenderingFreeType);
VTK_MODULE_INIT(vtkRenderingVolumeOpenGL);

int main()
{
    std::string str = "F:\\CT3";
    vtkSmartPointer<vtkDICOMImageReader> reader = vtkSmartPointer<vtkDICOMImageReader>::New();
    reader->SetDirectoryName(str.c_str());
    reader->SetDataByteOrderToLittleEndian();// 小端位元組序
    reader->Update();

    vtkSmartPointer<vtkImageCast> imageCast = vtkSmartPointer<vtkImageCast>::New();
    imageCast->SetInputConnection(reader->GetOutputPort());
    imageCast->SetOutputScalarTypeToUnsignedShort();
    imageCast->Update();

    vtkSmartPointer<vtkPiecewiseFunction> opacityFun = vtkSmartPointer<vtkPiecewiseFunction>::New();
    opacityFun->AddPoint(120,0.0);
    opacityFun->AddPoint(250,1.0);
    opacityFun->AddPoint(520,1.0);
    opacityFun->AddPoint(650,0.0);


    vtkSmartPointer<vtkColorTransferFunction> TransferFun = vtkSmartPointer<vtkColorTransferFunction>::New();
    TransferFun->AddRGBPoint(120, 255 / 255.0, 98 / 255.0, 98 / 255.0);
    TransferFun->AddRGBPoint(250, 255 / 255.0, 255 / 255.0, 180 / 255.0);
    TransferFun->AddRGBPoint(520, 1.0, 1.0, 1.0);
    TransferFun->AddRGBPoint(650, 1.0, 1.0, 1.0);

    vtkSmartPointer<vtkPiecewiseFunction> grideFun = vtkSmartPointer<vtkPiecewiseFunction>::New();
    grideFun->AddPoint(120, 2.0);
    grideFun->AddPoint(250, 2.0);
    grideFun->AddPoint(520, 0.1);
    grideFun->AddPoint(650, 0.1);

    vtkSmartPointer<vtkVolumeProperty> volumeProperty = vtkSmartPointer<vtkVolumeProperty>::New();
    volumeProperty->SetColor(TransferFun);
    volumeProperty->SetScalarOpacity(opacityFun);
    volumeProperty->SetGradientOpacity(grideFun);
    volumeProperty->ShadeOn();

    volumeProperty->SetAmbient(0.2);
    volumeProperty->SetDiffuse(0.9);
    volumeProperty->SetSpecular(0.2);
    volumeProperty->SetSpecularPower(10);

    vtkSmartPointer<vtkVolumeRayCastCompositeFunction> rayCastFun = vtkSmartPointer<vtkVolumeRayCastCompositeFunction>::New();

    vtkSmartPointer<vtkVolumeRayCastMapper> mapper1 = vtkSmartPointer<vtkVolumeRayCastMapper>::New();
    mapper1->SetVolumeRayCastFunction(rayCastFun);
    mapper1->SetInputConnection(imageCast->GetOutputPort());
    mapper1->SetAutoAdjustSampleDistances(0);





    vtkSmartPointer<vtkVolume> volume = vtkSmartPointer<vtkVolume>::New();
    volume->SetMapper(mapper1);
    volume->SetProperty(volumeProperty);
    vtkSmartPointer<vtkRenderer> ren = vtkSmartPointer<vtkRenderer>::New();
    ren->AddVolume(volume);
    ren->SetBackground(1.0,1.0,1.0);


    vtkSmartPointer<vtkRenderWindow> renwin = vtkSmartPointer<vtkRenderWindow>::New();
    renwin->AddRenderer(ren);
    renwin->SetSize(740, 480);

    vtkSmartPointer<vtkRenderWindowInteractor> iren = vtkSmartPointer<vtkRenderWindowInteractor>::New();
    iren->SetRenderWindow(renwin);

    iren->Initialize();
    renwin->Render();
    iren->Start();
    return 0;
}

使用vtkPiecewiseFunction 設定不透明度和梯度不透明vtkTransferFunctiony設定顏色傳輸。

這裡寫圖片描述

這裡寫圖片描述

相關推薦

VTK讀取DICOM醫學圖片進行體繪製的幾種方法

注意的是：VTK中不同的vtkVolumeMapper支援不同的資料型別。比如vtkVolumeRayCastMapper和vtkVolumeTextureMapper2D只能支援單分組VTK_UNSIGNED_CHAR和VTK_UNSIGNED_SHORT型別

關於圖片拉伸的幾種方法總結

系統至ios6之後，關於圖片拉伸的方法已經擴充套件至3個函式： 1.ios4提供的方法： - (UIImage *)stretchableImageWithLeftCapWidth:(NSInteger)leftCapWidth topCapHeight:(NSInt

讀取Spring的配置檔案applicationContext.xml幾種方法

1).利用ClassPathXmlApplicationContext，這種方式配置檔案應該放在類路徑下，否則Spring將找不到該檔案。 Java程式碼： ApplicationContext ctx=new ClassPathXmlApplicationContext("

圖片模糊處理的幾種方法

Blurring for Beginners  code download: http://download.csdn.net/detail/u012910985/7243241 Introduction This is a short tutorial on blurring techniques fo

如何加快VTK進行體繪製的速度

通過cmake編譯VTK,注意:configure的時候一定要選上BUILD_TESTING此選項(後邊有用),另外需要注意的是Directx SDK對應的選項時候已經更新(見下圖2)，然後不斷的configure，generate直到沒有錯誤產生。然後通過你的Visual Studio來build

【原始碼】以GUI讀取Dicom, V3D, ISI, NifTI, Analyze等16種醫學資料

函式ReadData3D允許使用者開啟醫學3D檔案。 This function ReadData3D allows the user to open medical 3D files. 本函式支援以下格式： Dicom Files ( .dcm , .dicom )

Android載入html中的svg格式圖片進行顯示的兩種方式

最近做的一個專案是把assets目錄中的html顯示出來,但是因為html裡面有一些工程圖片,雖然我用ViewPager和PhotoView,進行顯示放大了,但是因為工程圖片的線條較多還是比較模糊.所以後來就想用svg圖片來進行顯示,至於svg是什麼,我這裡就不

vtk讀取檔案並顯示的幾種方法

1.用vtkDICOMImageReader  #include "vtkRenderer.h" #include "vtkRenderWindow.h" #include "vtkRenderWindowInteractor.h" #include "vtkImageA

windows系統字體安裝使用的幾種方法

windows ps 字體簡易安裝方法：1、打開系統盤的 windows 下的 fonts 文件夾；2、用另一個窗口打開字體存放的文件夾；3、在存放字體的文件夾中選擇要安裝的字體（部分選或全選均可）；4、用鼠標右鍵將所選字體拖到windows 下的 fonts 文件夾上，松開鼠標，在彈出的提示框上選擇

Python 下載圖片的幾種方法

load color quest content 利用 pytho get .com 圖片 總結下： url = ‘http://img.pconline.com.cn/images/upload/upc/tx/photoblog/1707/31/c14/54293429_

QT中給各控件增加背景圖片(可縮放可旋轉)的幾種方法

.net detail eight iou rotate art board 按鈕 previous 1. 給QPushButton 增加背景圖片：背景圖片可根據Button大小自由縮放。 [cpp] view plain copy vo

HTML5 Canvas中繪製橢圓的幾種方法

1.canvas自帶的繪製橢圓的方法 ellipse(x, y, radiusX, radiusY, rotation, startAngle, endAngle, anticlockwise)是後來新增的， 引數的意思：(起點x.起點y,半徑x,半徑y,旋轉的角度，起始角，結果角，順時針還是逆時針)

python os 設定讀取環境變數的幾種方法

專案中的跳板機是用的使用者公私鑰鑑權的，登入跳板機需要使用的每個人的name，因此如果將封裝的MySQL或者Redis方法中傳入USERNAME作為引數的話，感覺有點傻。因此，考慮將名字寫入環境變數。 兩種方法： 1、使用os.environ.setdefault方法寫入(或者直接

如何放大圖片，教您三種方法！

在日常生活中，我們經常會對一些尺寸不符合的照片進行放大處理，在這裡我們就需要一些特殊的方法，接下來我們介紹三種方法。   最懶的方法：直接拉伸小圖擷取   簡易程度：★★★★★   圖片質量：★   適用人群：只需要圖片尺寸不需要精度的懶人  

圖片上傳前端的幾種方法

      1.路徑 <img src="balabala/123.jpg"> 2.img標籤指向servlet方法 <img src="/ImgServlet?imgId=00001" alt="picture">

css圖片垂直居中 讓html img圖片垂直居中的三種方法

本教程為thinkcss教大家三種讓img元素圖片在盒子內垂直居中的方法教程，根據程式碼與文章教程理解掌握並加以使用。 一、使用flex實現垂直居中 利用css flex實現垂直居中。flex可能不是實現垂直居中最好的選擇，因為IE8,9並不支援它。 現在，為了用flex實現垂

R語言入門到放棄 R語言讀取不同檔案型別中資料的4種方法

R語言入門到放棄 R語言讀取檔案中資料的4中方法 方法一：直接讀取 > getwd() # 獲取當前的工作路徑 將要讀取的檔案剪下到當前的工作路徑，下面讀取檔案 x<-read.table("text.txt") Warning message

給定兩個變數將兩個變數的值內容進行交換的三種方法

建立一個void swap(int a,int b){}函式，賦值作為被調函式 void swap(int *p1,int *p2) { int tmp = *p1; int *p1 = *p2; int *p2 = tmp; } int main() { in

Java讀取檔案的幾種方法

public class ReadFromFile { /** * 以位元組為單位讀取檔案，常用於讀二進位制檔案，如圖片、聲音、影像等檔案。 */ public static void readFileByBytes(String fileName) { Fil

QT 中設定按鈕圖片和文字的兩種方法

1.使用QpushButton自帶的API實現： void setIcon(const QIcon &icon) void setText(const QString &text) 該方法使用了自帶的函式設定，缺點是無法自由調整圖示和文字之間的間距 2.使用Qlabel和佈局