Erictoo6 2022-02-20 12:41 采纳率: 66.7%
浏览 110
已结题

vtk读取dicom文件提取像素矩阵


使用vtk读取dicom文件后,想将像素值信息提取出来用OpenCV处理,但是卡在了不知道怎么提取像素值信息。尝试了下述代码

#include <vtkSmartPointer.h>
#include <vtkImageViewer2.h>
#include <vtkImageCast.h>
#include <vtkDICOMImageReader.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
#include <vtkRenderer.h>
#include <vtkAutoInit.h>
#include <opencv2/core.hpp>
#include <opencv2/highgui.hpp>
#include <opencv2/imgproc.hpp>
#include <typeinfo>
#include <iostream>
VTK_MODULE_INIT(vtkRenderingOpenGL2)

using namespace cv;
using namespace std;

void dicomread(string inputFilename, Mat& img, vtkSmartPointer<vtkDICOMImageReader>& reader);

int main(int argc, char* argv[])
{
    // Verify input arguments
    //if (argc != 2)
    //{
    //  std::cout << "Usage: " << argv[0]
    //      << " Filename(.img)" << std::endl;
    //  return EXIT_FAILURE;
    //}

    std::string inputFilename = "C:\\Users\\26327\\Desktop\\dll4sw\\data\\SE3\\IM0.dcm";
    Mat img;
    // Read all the DICOM files in the specified directory.
    vtkSmartPointer<vtkDICOMImageReader> reader =
        vtkSmartPointer<vtkDICOMImageReader>::New();
    reader->SetFileName(inputFilename.c_str());
    reader->Update();
    dicomread(inputFilename, img, reader);
    // 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;
}

void dicomread(string inputFilename, Mat& img, vtkSmartPointer<vtkDICOMImageReader>& reader)
{
    img.create(512, 512, CV_32SC1);

    vtkSmartPointer<vtkImageCast> imageCast =
        vtkSmartPointer<vtkImageCast>::New();

    reader->SetFileName(inputFilename.c_str());
    vtkImageData* p;
    reader->Update();
    cout << typeid(reader->GetOutput()).name() << endl;

    imageCast->SetInputConnection(reader->GetOutputPort());
    imageCast->SetOutputScalarTypeToInt();
    imageCast->Update();

    int dims[3];
    reader->GetOutput()->GetDimensions(dims);

    for (int k = 0; k < dims[2]; k++)
    {
        for (int j = 0; j < dims[1]; j++)
        {
            for (int i = 0; i < dims[0]; i++)
            {
                int* pixel =
                    (int*)(imageCast->GetOutput()->GetScalarPointer(i, j, k)); 
                img.at<int>(j, i) = int(*pixel); 
            }
        }
    }
}

但是这两行都报错,报错信息都是“pointer to incomplete class type "vtkImageData" is not allowed"。请问如何解决


    reader->GetOutput()->GetDimensions(dims);

    int* pixel = (int*)(imageCast->GetOutput()->GetScalarPointer(i, j, k)); 

  • 写回答

1条回答 默认 最新

  • 黑山老妖的笔记本 2022-02-20 13:47
    关注

    加上#include<vtkImageData.h>

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

问题事件

  • 系统已结题 2月28日
  • 已采纳回答 2月20日
  • 创建了问题 2月20日

悬赏问题

  • ¥15 如何删除这个虚拟音频
  • ¥50 hyper默认的default switch
  • ¥15 网站打不开,提示502 Bad Gateway
  • ¥20 基于MATLAB的绝热压缩空气储能系统代码咨询
  • ¥15 R语言建立随机森林模型出现的问题
  • ¥15 中级微观经济学,生产可能性边界问题
  • ¥15 TCP传输时不同网卡传输用时差异过大
  • ¥15 请各位看看我写的属于什么算法,或者有更正确的写法?
  • ¥15 html5 qrcode 扫描器
  • ¥15 爬取网页信息并保存需要完整代码