VTK: 转化DICOM ST 数据为vtkPolyData显示

原创

已于 2022-09-26 11:42:06 修改 · 1.4k 阅读

2 ·

CC 4.0 BY-SA版权

文章标签：

#c++ #visual studio #算法

于 2022-09-26 11:41:17 首次发布

本文介绍如何将DICOM-RT文件中的放疗靶区数据转换为vtkPolyData格式，包括坐标转换、掩码生成及三维渲染等关键步骤，并讨论了存在的问题。

一、简介

放疗的靶区等勾画数据通常会存储为DICOM-RT文件，该文件通常存储勾画点坐标数据。新的需要中DICOM-RT中存储了一个固定模版，需要进行3D打印，因此需要对其格式进行转化，以便进行三维渲染、编辑、保存和打印功能。

处理流程：读取DICOM-RT文件中的坐标信息（设备），将坐标转化为图象坐标系下，根据每层坐标，构建闭合曲线，填充获得一张数据（和原始DICOM数据同大小，数据类型为unsigned char的vtkImageData），依次完成对所有层的操作，最后将三维数据转化为vtkPolyData。

二、坐标转化并生成vtkImageData，验证其正确性

坐标转换：读取DICOM文件中的[ (0020,0032), DS, Image Position(Patient) ]，对从小到大排序；读取DICOM-RT文件中的所有坐标信息，将其通过平移变换到图像空间。坐标变换到图象空间为： X Y Z

将点坐标转化为闭合的掩码，此步骤采用opencv来完。

// pData 为指向某层图像的指针
cv::Mat dstMat(mDim[1], mDim[0], CV_8UC1, pData);
// 绘制轮廓
cv::drawContours(dstMat, smoothContours, 0, cv::Scalar(255, 255, 255), -1);

将绘制的图像保存为jpg,以

最低0.47元/天解锁文章