1.前言

VTK應用程序所需的數(shù)據(jù)可以通過兩種途徑獲取: 第一種是生成模型?;第二種是從外部存儲介質(zhì)里導入相關的數(shù)據(jù)文件,(如vtkBMPReader讀取 BMP圖像) 。VTK 也可以將程序中處理完成的數(shù)據(jù)寫入單個文件中, 或者將所渲染的場景導出。從可視化管線的角度來看,一般以數(shù)據(jù)的讀取 (或由模型創(chuàng)建數(shù)據(jù))開始,而以數(shù)據(jù)的寫盤操作(或 Mapper)結(jié)束。

前面我們已經(jīng)接觸到了 VTK的 Reader類,將數(shù)據(jù)導入可視化管道的流程：

• 實例化 Reader對象；
• 指定所要讀取的文件名；
• 調(diào)用 Update()方法促使管線執(zhí)行。當管線后續(xù)的 Fiter有 Update0請求時,如調(diào)用Render()方法管線就會讀取相應的圖像文件, 所以這一步可省略 。

類似地, 使用 Writer類的主要步驟如下:

• 實例化 Writer對象;
• 輸入要寫盤的數(shù)據(jù)以及指定待寫盤的文件名;
• 調(diào)用 Write()方法促使 Writer類開始寫盤操作。

2.vtkImageData類

圖像數(shù)據(jù)在 VTK 中是用vtkImageData類表示的,對于不同的圖像文件類型, VTK提供相對應的類對圖像文件進行讀寫操作。比如,前面章節(jié)中所提的vtkBMPReader是用于讀取 BMP圖像, vtkJPEGReader用于讀取 JPG圖像。 VTK除了支持 BMP、 JPG圖像格式之外,還支持其他多種圖像格式的讀寫。例如：
注意： 1.值得注意的是vtkImageReader/vtkImageWriter用于讀寫RAW格式的數(shù)據(jù)(即俗稱的“裸數(shù)據(jù)”) ,該類型的圖像沒有文件信息,因此在讀取此類圖像時,需要指定圖像各個維度的大小、字節(jié)順序(是大端字節(jié)序還是小端字節(jié)序)、存儲像素値的類型等信息,只有指定這些信息、,類vtkImageReader才能正確讀取圖像。
2.類vtkDicomImageReader可用于讀取DICOM圖像,但該類的功能很不完善,雖然VTK 最初是因醫(yī)學圖像可視化而誕生。但VTK對DICOM圖像的讀寫操作卻很不支持,比如該類不支持多幀DICOM圖像的讀取，而且VTK 也沒有實現(xiàn)對 DICOM圖像的寫操作, 即沒有提供類vlkDlCOMlmageWriter。 對 DICOM圖像支持較好的函數(shù)庫主要有 GDCM和 DCMTK。著名的醫(yī)學圖像分割與配準工具包 ITK 就是封裝了 GDCM函數(shù)庫進行 DICOM圖像的讀寫。

3.單幅圖像的讀寫

#include <vtkAutoInit.h> VTK_MODULE_INIT(vtkRenderingOpenGL); #include <vtkSmartPointer.h> #include <vtkPNGReader.h> #include <vtkImageViewer2.h> #include <vtkRenderWindowInteractor.h> #include <vtkJPEGWriter.h> #include <vtkRenderer.h> //定義了Cameraint main() {//讀取PNG圖像vtkSmartPointer<vtkPNGReader> pngread = vtkSmartPointer<vtkPNGReader>::New();pngread->SetFileName("logo.png");//顯示該幅圖像vtkSmartPointer<vtkImageViewer2> ImageViewer =vtkSmartPointer<vtkImageViewer2>::New();ImageViewer->SetInputConnection(pngread->GetOutputPort());vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInter = vtkSmartPointer<vtkRenderWindowInteractor>::New();ImageViewer->SetupInteractor(renderWindowInter);ImageViewer->Render();ImageViewer->GetRenderer()->ResetCamera();ImageViewer->Render();//寫圖像vtkSmartPointer<vtkJPEGWriter> writer = vtkSmartPointer<vtkJPEGWriter>::New();writer->SetFileName("VTK-logo.jpg");writer->SetInputConnection(pngread->GetOutputPort());writer->Write();renderWindowInter->Start(); }? ?? ?
該段代碼先使用vtkPNGReader讀入 PNG圖像,然后用VTK 窗口顯示讀取的 PNG圖像,最后使用類vtkJPEGWriter將讀入的文件寫成JPG圖像。程序中使用 SetFileName方法設置要讀寫的圖像名, 在寫文件操作時要調(diào)用方法 Write方法才會將內(nèi)存中的數(shù)據(jù)寫入到存儲介質(zhì)中。?

可能我們會感覺到很奇怪，數(shù)據(jù)管線已經(jīng)鋪設好了，可是我們是怎么樣實現(xiàn)引擎渲染的呢？已到最后完成了與Windows操作系統(tǒng)完成了交互？在顯示圖像時并沒有用到，vtkRenderWindow、 vtkRenderer、 vtkActor等類,而只是使用了vtkImageViewer2以及設置了交互樣式。其實, VTK 可視化管線相關的幾個類都已經(jīng)封裝在vtkImageViewer2 里。vtkImageViewer2主要是針對二維圖像(特別是醫(yī)學圖像)顯示設計的,實現(xiàn)了圖像縮放、轉(zhuǎn)、平移、窗寬窗位調(diào)節(jié)等功能: 除了可以用于單幅二維圖像的顯示之外, 也可以顯示三維圖像的某個切片, 還可以設置不同的顯示方向。

4.讀取序列圖像文件

醫(yī)學圖像應用程序中常常會處理序列的圖像文件,比如計算機斷層成像或者磁共振成像所成圖像一般都是由多個有順序的二維圖像組成, 應用程序需要一次性導入一個序列的二維圖像。VTK沒有提供專門的類讀取序列圖像文件，但是VTK的圖像Reader類都有提供方法SetFileNames()來設置多個圖像文件名，利用該方法可以實現(xiàn)序列圖像的讀取。
#include <vtkAutoInit.h> VTK_MODULE_INIT(vtkRenderingOpenGL); #include <stdio.h> #include <vtkSmartPointer.h> #include <vtkStringArray.h> #include <vtkJPEGReader.h> #include <vtkImageViewer2.h> #include <vtkRenderWindowInteractor.h> #include <vtkRenderer.h>int main() {//生成文件序列組名vtkSmartPointer <vtkStringArray> fileArray =vtkSmartPointer <vtkStringArray>::New();char fileName[128];for(int i=1; i<100; i++){sprintf(fileName,"Head/head%03d.jpg",i);vtkstd::string fileStr(fileName);fileArray->InsertNextValue(fileStr);}//讀取JPG序列圖像vtkSmartPointer <vtkJPEGReader> reader = vtkSmartPointer <vtkJPEGReader>::New();reader->SetFileNames(fileArray);//顯示vtkSmartPointer<vtkImageViewer2> viewer = vtkSmartPointer<vtkImageViewer2>::New();viewer->SetInputConnection(reader->GetOutputPort());vtkSmartPointer<vtkRenderWindowInteractor> interact =vtkSmartPointer<vtkRenderWindowInteractor>::New();//默認選擇第50張切片viewer->SetSlice(50);//viewer->SetSliceOrientationToXY();viewer->SetSliceOrientationToXZ();//viewer->SetSliceOrientationToYZ();viewer->SetupInteractor(interact);viewer->Render();interact->Start();return 0; } 程序中使用 vtkStringArray生成文件名列表, 然后調(diào)用 vtkJPEGReader的方法。SetFileNames()設置符讀取的序列圖像的文件名 。初始顯示該序列圖像的第50個切片,顯示的方向為 SetSliceOrientationToXY()。
還可以使用 Reader類里的 SetFilePrefix()和 SetFilePattem()等方法讀取序列圖像。

5.參考資料

1.《C++ primer》
2.《The VTK User’s Guide – 11thEdition》
3.《The Visualization Toolkit – AnObject-Oriented Approach To 3D Graphics (4th Edition)》
4. ?張曉東, 羅火靈. VTK圖形圖像開發(fā)進階[M]. 機械工業(yè)出版社, 2015.

總結(jié)

以上是生活随笔為你收集整理的VTK修炼之道13：数据读写_图像数据的读写的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。