VTK 结构化点集 vtkStructuredPoints

目录

引言：

Reader

---

引言：

今天朋友给了一份数据是OCT 采集到的数据， VTK 格式，但是我打不开；

用文本打开看了一下；

vtk DataFile Version 2.0

BINARY
DATASET STRUCTURED_POINTS
DIMENSIONS 512 512 512
SPACING 0.0035 0.00791947 0.010472
ORIGIN 0 0 0
POINT_DATA 134217728
SCALARS volume_scalars float 1
LOOKUP_TABLE default

搜索了一下，是结构化点集数据；

vtkStructuredPoints

vtkStructuredPoints类继承自vtkImageData，是图像数据的一种，要求数据范围与更新范围完全匹配。矢量图数据允许数据范围大于更新范围。StructuredPoints类对原点的定义也与vtkImageData不同。对于结构化点（structured points），原点是第一个点的位置。而图像将原点定义为点(0,0,0)的位置。图像原点存储在ivar中，结构化点具有Set/GetOrigin/Extents的特殊方法。
 

可以拿 它当 vtkImageData用；

Reader

vtkStructuredPointsReader

读取，转polydata

			 auto reader = vtkSmartPointer<vtkStructuredPointsReader>::New();
			 reader->SetFileName(pathTwo.toStdString().c_str());
			 reader->Update();

显示：

#include <vtkStructuredPointsReader.h>
#include <vtkStructuredPoints.h>
    
 // init imageview
	mImageView = vtkSmartPointer<vtkImageViewer2>::New();
	auto interactor = vtkSmartPointer<vtkRenderWindowInteractor>::New();
	mImageView->SetRenderWindow(mVTKWidget->GetRenderWindow());
	mImageView->SetupInteractor(interactor);

 mImageView->SetInputData(reader->GetOutput());
			 mImageView->SetColorLevel(100);
			 mImageView->SetColorWindow(200);
	 
			 mImageView->SetSliceOrientationToYZ();
			 mImageView->GetRenderer()->SetBackground(0, 0, 0);
			 mImageView->SetSlice((mImageView->GetSliceMin() + mImageView->GetSliceMax()) / 2);
			 mImageView->Render();

转vtkPolydata 显示；

		 vtkSmartPointer<vtkImageDataGeometryFilter> geometryFilter =
				 vtkSmartPointer<vtkImageDataGeometryFilter>::New();
			 geometryFilter->SetInputConnection(reader->GetOutputPort());
			 geometryFilter->Update();

or

	vtkSmartPointer< vtkMarchingCubes > boneExtractor =
		vtkSmartPointer< vtkMarchingCubes >::New();
	boneExtractor->SetInputData(mImageOne);
 
	boneExtractor->SetValue(0, 30);
	//boneExtractor->SetNumberOfContours(1);
	boneExtractor->Update();
 

	//mVtkPolyData->DeepCopy(connectivityFilter->GetOutput());
	//剔除旧的或废除的数据单元，提高绘制速度(可略去这一步)  
	vtkSmartPointer< vtkStripper > boneStripper =
		vtkSmartPointer< vtkStripper >::New(); //三角带连接  
	boneStripper->SetInputConnection(boneExtractor->GetOutputPort());
	boneStripper->Update();
	mVtkPolyData->DeepCopy(boneStripper->GetOutput());


   		 mPolyMapper->SetInputData(mVtkPolyData);
		 mPolyMapper->ScalarVisibilityOff();
		 mPolyMapper->Update();
		 m3DRenderer->ResetCamera();
		 m3DRenderer->ResetCameraClippingRange();
		 m3DRenderer->ResetCamera();
		 m3DRenderer->Render();