Question

我正在尝试迭代dicom文件中的所有像素。现在我正在使用演示：

int main(int argc, char* argv[])
{
  std::string inputFilename = "MR000000";

//   Read all the DICOM files in the specified directory.
  vtkSmartPointer<vtkDICOMImageReader> reader =
    vtkSmartPointer<vtkDICOMImageReader>::New();
 reader->SetFileName(inputFilename.c_str());
reader->Update();

// Visualize
vtkSmartPointer<vtkImageViewer2> imageViewer =
vtkSmartPointer<vtkImageViewer2>::New();
imageViewer->SetInputConnection(reader->GetOutputPort());
vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor =
vtkSmartPointer<vtkRenderWindowInteractor>::New();
imageViewer->SetupInteractor(renderWindowInteractor);
imageViewer->Render();
imageViewer->GetRenderer()->ResetCamera();
imageViewer->Render();

 renderWindowInteractor->Start();

return EXIT_SUCCESS;
}

这是一种从读者对象中获取像素数组的方法吗？

Answer 1

vtkImageViewer2有GetIntput()方法，可返回vtkImageData个实例。我想这个方法应该在imageViewer->Render()调用之后使用。

通过此ImageData，您可以检索图像的每个组件。您需要通过imageData->GetScalarPointer(x,y,z)调整返回组件的类型;使用imageData->GetScalarType ()和imageData->GetNumberOfScalarComponents()。

看看这个exemple

Answer 2

vtkImageViewer2不是你想要的，vtkDICOMImageReader更适合：

以下是一些应该有帮助的来源（这取决于pcl点，但实际上是相同的）：

vtkSmartPointer<vtkDICOMImageReader> reader = vtkSmartPointer<vtkDICOMImageReader>::New();
vtkSmartPointer<vtkImageData> sliceData = vtkSmartPointer<vtkImageData>::New();
reader->SetFileName(dicomFile.c_str());
reader->Update();
sliceData = reader->GetOutput();

int numberOfDims = sliceData->GetDataDimension();
int * dims = sliceData->GetDimensions();
double* spacingData = reader->GetDataSpacing();

for (int y = 0; y < dims[1]; y++)
        {
            for (int x = 0; x < dims[0]; x++)
            {
                pcl::PointXYZI tempPt = pcl::PointXYZI();
                tempPt.x = x * spacingData[0];
                tempPt.y = y * spacingData[1];
                tempPt.z = 0;
                double tempIntensity = sliceData->GetScalarComponentAsDouble(x, y, z, 0);
                if (!isinf(tempIntensity))
                {
                    tempPt.intensity = tempIntensity;
                }
                else
                {
                    tempPt.intensity = 0;
                }

                cloud->points.push_back(tempPt);
            }
        }

Vtk读取dicom文件

2 个答案: