裁剪Dicom系列

Question

我需要裁剪Dicom系列，并将结果再次保存为Dicom系列，以维护有关位置的信息。我的问题与DICOM File compression类似，但不幸的是，所有的答案都是关于压缩的，而不是裁剪。

另一个问题是，原始系列非常大，我无法将其全部加载到计算机的内存中(但如果需要，我可以使用Unix服务器进行计算)。

我尝试过使用ITK:在编写之前，我使用cropimageFilter复制了示例DicomSeriesReadSeriesWrite.cxx (在ITK软件指南中描述)。问题是，使用指令

seriesWriter->SetMetaDataDictionaryArray(reader->GetMetaDataDictionaryArray() );

在结果中，图像位置的第三个维度(0020,0037)始终设置为0。如果我注释该指令，图像位置是正确的，但图像编号(0020,0013)和SliceLocation (0020,1041)不正确，且结果系列不会加载到某些软件中。

我也尝试过使用Dicom2 (http://www.barre.nom.fr/medical/dicom2/ )。同样，图像被裁剪得很好，但元数据并不一致(它们与原始图像相同)。

我试过使用Amira，但是，除了内存问题，将裁剪的数据保存为dicom的过程忽略了原始的dicom标签(例如，图像位置是-1.#ND/-1.#ND/-1.#ND，其他标签也是“系列日期”更改)。

Answer 1

由于您知道标签的原始值以及图像的裁剪方式，因此可以在保存dicom文件之前将标签修改为正确的值。